Structural and Thermodynamic Properties of the CmIII Ion Solvated by Water and Methanol

Solution Structures of Europium Terpyridyl Complexes with Nitrate and Triflate Counterions in Acetonitrile

Lanthanide-ligand complexes are key components of technological applications, and their properties depend on their structures in the solution phase, which are challenging to resolve experimentally or computationally. The coordination structure of the Eu³⁺ ion in different coordination environments in acetonitrile is examined using ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy. AIMD simulations are conducted for the solvated Eu³⁺ ion in acetonitrile, both with or without a terpyridyl ligand, and in the presence of either triflate or nitrate counterions. EXAFS spectra are calculated directly from AIMD simulations and then compared to experimentally measured EXAFS spectra. In acetonitrile solution, both nitrate and triflate anions are shown to coordinate directly to the Eu³⁺ ion forming either ten- or eight-coordinate solvent complexes where the counterions are binding as bidentate or monodentate structures, respectively. Coordination of a terpyridyl ligand to the Eu³⁺ ion limits the available binding sites for the solvent and anions. In certain cases, the terpyridyl ligand excludes any solvent binding and limits the number of coordinated anions. The solution structure of the Eu-terpyridyl complex with nitrate counterions is shown to have a similar arrangement of Eu³⁺ coordinating molecules as the crystal structure. This study illustrates how a combination of AIMD and EXAFS can be used to determine how ligands, solvent, and counterions coordinate with the lanthanide ions in solution.

Hierarchical phenomena in multicomponent liquids: simulation methods, analysis, chemistry

Complex, multicomponent, solutions have often been studied solely through the lens of specific applications of interest. Yet advances to both simulation methodologies (enhanced sampling, etc.) and analysis techniques (network analysis algorithms and others), are creating a trove of data that reveal transcending characteristics across vast compositional phase space. This perspective discusses technical considerations of the reliable and accurate simulations of complex solutions, followed by the advances to analysis algorithms that elucidate coupling of different length and timescale behavior (hierarchical phenomena). The different manifestations of hierarchical phenomena are presented across an array of solution environments, emphasizing fundamental and ongoing science questions. With a more advanced molecular understanding in hand, a quintessential application (solvent extraction) is discussed, where significant opportunities exist to re-imagine the technical scope of an established technology.

PageRank as a collective variable to study complex chemical transformations and their energy landscapes

A reduced set of reaction coordinates is often employed in chemistry to describe the collective change between reactants and products within the context of rare event theories and the exploration of energy landscapes. Yet selecting the proper collective variable becomes increasingly challenging as the systems under study become more complex. Recent advancement of new descriptions of collective molecular coordinates has included graph-theoretical metrics, including social permutation invariant and PageRank (PR) coordinates, based upon the network of interactions about molecules and atoms within a system. Herein we continue the development of PR by (1) presenting a new formulation that is continuous along a reaction path, (2) illustrating that the fluctuations in PR are demonstrative of the fundamental motions of the atoms/molecules, and (3) providing the analytical derivatives with respect to atomic coordinates. The latter is subsequently combined with a harmonic bias to create the potential of mean force (PMF). As an example, we first consider the transformation of tetrahedral [Al(OH)₄]_(aq) ^- to octahedral [Al(OH)₄(H₂O)₂]_(aq) ^- using the PR PMF. Second, we explore the interchange of contact ion pair and solvent separated ion pairs of aqueous Na⋯OH, where the distance-biased PMF is projected onto PR space. In turn, this reveals where solvent rearrangement has the most impact upon the reaction pathway.

Binding affinity of pyridines with AmIII/CmIII elucidated by density functional theory calculations

In recent decades, N-heterocyclic ligands have been extensively used in the separation of lanthanides/actinides, whereas the selective extraction of amercium or curium has been very challenging. Using density functional theory calculations, this study is devoted to the investigation of the binding affinity of a series of modelling pyridine ligands with Am^III and Cm^III. The structure-property correlations between the amercium and curium systems and the binding affinity were obtained, and promising strategies for efficient separation of Am^III/Cm^III were proposed.

Competitive Interactions Within Cm(III) Solvation in Binary Water/Methanol Solutions

Competitive forces exist in multicomponent solutions, and within electrolytes they consist of both ion-solvent and solvent-solvent interactions. These can influence a myriad of processes, including ligand complexation. In the case of water/alcohol solutions, recent work revealed an interesting dilemma regarding the overall solution dynamics and organization as compared to solute-solvent interactions. This is particularly true for highly charged ions in solution, whose ion-solvent interactions were demonstrated to be highly sensitive to the composition of the immediate solvation environment. Faster solvent exchange should be observed about the ion, considering that second-order Møller-Plesset perturbation theory predicts an average decrease in ion-solvent dissociation energy when methanol enters the first solvation shell of Cm³⁺_(aq). Yet the addition of methanol to water causes the dynamic features of the hydrogen-bond network of the entire solution to slow. The apparent competition between these contrary forces was examined using a combination of electronic structure calculations with both ab initio and classical molecular dynamics simulations, using binary water/methanol solutions and Cm³⁺ as a representative solute. This combination of theoretical methods predicts that, among the competitive effects of the solvent-solvent and ion-solvent interactions, the solution-phase dynamics imparted by the addition of methanol to water kinetically restricts the solvation exchange rates about Cm³⁺ in these binary solutions.

Spectroscopic and Computational Characterization of Diethylenetriaminepentaacetic Acid/Transplutonium Chelates: Evidencing Heterogeneity in the Heavy Actinide(III) Series

The chemistry of trivalent transplutonium ions (Am³⁺ , Cm³⁺ , Bk³⁺ , Cf³⁺ , Es³⁺ …) is usually perceived as monotonic and paralleling that of the trivalent lanthanide series. Herein, we present the first extended X-ray absorption fine structure (EXAFS) study performed on a series of aqueous heavy actinide chelates, extending past Cm. The results obtained on diethylenetriaminepentaacetic acid (DTPA) complexes of trivalent Am, Cm, Bk, and Cf show a break to much shorter metal-oxygen nearest-neighbor bond lengths in the case of Cf³⁺ . Corroborating those results, density functional theory calculations, extended to Es³⁺ , suggest that the shorter Cf-O and Es-O bonds could arise from the departure of the coordinated water molecule and contraction of the ligand around the metal relative to the other [M^III DTPA(H₂ O)]^2- (M=Am, Cm, Bk) complexes. Taken together, these experimental and theoretical results demonstrate inhomogeneity within the trivalent transplutonium series that has been insinuated and debated in recent years, and that may also be leveraged for future nuclear waste reprocessing technologies.

On the Origin of Covalent Bonding in Heavy Actinides

Recent reports have suggested the late actinides participate in more covalent interactions than the earlier actinides, yet the origin of this shift in chemistry is not understood. This report considers the chemistry of actinide dipicolinate complexes to identify why covalent interactions become more prominent for heavy actinides. A modest increase in measured actinide:dipicolinate stability constants is coincident with a significant increase in An 5f energy degeneracy with the dipicolinate molecular orbitals for Bk and Cf relative to Am and Cm. While the interactions in the actinide-dipicolinate complex are largely ionic, the decrease in 5f orbital energy across the series manifests in orbital-mixing and, hence, covalency driven by energy degeneracy. This observation suggests the origin of covalency in heavy actinide interactions stems from the degeneracy of 5f orbitals with ligand molecular orbitals rather than spatial orbital overlap. These findings suggest that the limiting radial extension of the 5f orbitals later in the actinide series could make the heavy actinides ideal elements to probe and tune effects of energy degeneracy driven covalency.

MD simulation study of the diffusion and local structure of n-alkanes in liquid and supercritical methanol at infinite dilution

The diffusion coefficients of 14 n-alkanes (ranging from methane to n-tetradecane) in liquid and supercritical methanol at infinite dilution (at a pressure of 10.5 MPa and at temperatures of 299 K and 515 K) were deduced via molecular dynamics simulations. Values for the radial distribution function, coordination number, and number of hydrogen bonds were then calculated to explore the local structure of each fluid. The flexibility of the n-alkane (as characterized by the computed dihedral distribution, end-to-end distance, and radius of gyration) was found to be a major influence and hydrogen bonding to be a minor influence on the local structure. Hydrogen bonding reduces the flexibility of the n-alkane, whereas increasing the temperature enhances its flexibility, with temperature having a greater effect than hydrogen bonding on flexibility. Graphical abstract The flexibility of the alkane is a major influence and the hydrogen bonding is a minor influence on the first solvation shell; the coordination numbers of long-chain n-alkanes in the first solvation shell are rather low.