Aqueous Solvation of SmI2: A Born–Oppenheimer Molecular Dynamics Density Functional Theory Cluster Approach

Exploring the Dynamic Coordination Sphere of Lanthanide Aqua Ions: Insights from r2SCAN-3c Composite-DFT Born-Oppenheimer Molecular Dynamics Studies

Born-Oppenheimer molecular dynamics (BOMD) simulations were performed to investigate the structure and dynamics of the first hydration shells of five trivalent lanthanide ions (Ln3+) at room temperature. These ions are relevant in various environments, including the bulk aqueous solution. Despite numerous studies, accurately classifying the molecular geometry of the first hydration sphere remains a challenge. To addres this, a cluster microsolvation approach was employed to study the interaction of Ln3+ ions (La, Nd, Gd, Er, and Lu) with up to 27 explicit water molecules. Electronic structure calculations were performed with the composite r2SCAN-3c method. The results demonstrate that this method offers an optimal balance between precision and computational efficiency. Specifically, it accurately predicts average Ln-O distances (MAE = 0.02 Å) of the first hydration sphere and preferred coordination numbers (CN) for the different lanthanide cations as compared to reported data in bulk. Highly dynamic first hydration shells for the examined Ln3+ ions were found, with noticeable and rapid rearrangements in their coordination geometries, some of which can be recognized as the tricapped trigonal prism (TTP) and the capped square antiprism (CSAP) for CN = 9, and as the square antiprism (SAP), the bicapped trigonal prism (BTP), and the trigonal dodecahedron (DDH) for CN = 8. However, ca. 70% of the nonacoordinated configurations did not meet the criteria of TTP or CSAP structures. For CN = 8, the percentage of configurations that could not be assigned to SAP, BTP, or DDH was lower, around 30%. The theoretical EXAFS spectra obtained from the BOMD simulations are in good agreement with the experimental data and confirm that model microsolvated environments accurately represent the near-solvation structure of these trivalent rare-earth ions. Moreover, this demonstrates that the faster dynamics of the first hydration shell can be studied separately from the dynamics of water exchange in the bulk aqueous solution.

Accuracy of Discrete-Continuum Solvation Model for Cations: A Benchmark Study

Metal ions play important roles in chemistry, biochemistry, and material sciences. Accurately modeling ion solvation is crucial for simulating ion-containing systems. There are different models for ion solvation in computational chemistry, such as the explicit model, continuum model, and discrete-continuum model. Compared to the explicit model and continuum model, the discrete-continuum model of solvation is a hybrid solvation model in which the first solvation shell is described explicitly, and the remainder of the bulk liquid is characterized by a continuum model, which provides an excellent balance between accuracy and computational costs. This work serves as a systematic benchmark of the discrete-continuum model for the solvation of cations with +2, +3, and +4 charges. The calculated hydration free energies (HFEs) of ions were compared to those obtained by the SMD continuum model alone and the available experimental data. The discrete-continuum model showed improved performance over the continuum model alone via a smaller overall error and more consistent performance. Experimentally observed trends, such as the Irving-Williams series, are generally reproduced. In contrast, greater overall error was obtained for Ln3+ ions, and the HFE trend along the Ln3+ series was more difficult to reproduce, indicating these ions are challenging to model by the discrete-continuum model and continuum model. Overall, the discrete-continuum model is recommended to calculate the HFEs of cations when experimental data are not available.

Diastereoselective Dearomatizing Cyclizations of 5-Arylpentan-2-ones by Samarium Diiodide - A Computational Analysis

This study analyzes the samarium diiodide-promoted cyclizations of 5-arylpentan-2-ones to dearomatized bicyclic products utilizing density functional theory. The reaction involves a single electron transfer to the carbonyl group, which occurs synchronously with the rate determining cyclization event, and a second subsequent proton-coupled electron transfer. These redox reactions are accurately computed employing small core pseudo potentials explicitly involving all f-electrons of samarium. Comparison of the energies of the possible final products rules out thermodynamic control of the observed regio- and diastereoselectivities. Kinetic control via appropriate transition states is correctly predicted, but to obtain reasonable energy levels the influence of the co-solvent hexamethylphosphortriamide has to be estimated by using a correction term. The steric effect of the bulky samarium ligands is decisive for the observed stereoselectivity. Carbonyl groups in para-position of the aryl group change the regioselectivity of the cyclization and lead to spiro compounds. The computations suggest again kinetic control of this deviating outcome. However, the standard mechanism has to be modified and the involvement of a complex activated by two SmI2 moieties is proposed in which two electrons are transferred simultaneously to form the new C-C bond. Computation of model intermediates show the feasibility of this alternative+ mechanism.

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.

Coordination Structure of Samarium Diiodide in a Tetrahydrofuran-Water Mixture

Chemoselective reductive conversion of organic and inorganic compounds has been developed by the combination of samarium(II) diiodide (SmI₂) and water. Despite the extensive previous studies to elucidate the role of water in the reactivity of SmI₂, the direct structural data of the reactive Sm²⁺-water complexes, SmI₂(H₂O)_n, in an organic solvent-water mixture have not been reported experimentally so far. Herein, we performed the structure analysis of the Sm²⁺-water complex in tetrahydrofuran (THF) in the presence of water by in situ X-ray absorption spectroscopy using high-energy X-rays (Sm K-edge, 46.8 keV). The analysis revealed the dissociation of the Sm²⁺-I bonds in the presence of ≥ eight equivalents of water in the THF-water mixture. The origin of the peak shift in the UV/visible absorption spectra after the addition of water into SmI₂/THF solution was proposed based on electron transitions simulated with time-dependent density-functional-theory calculations using optimized structures in THF or water. The obtained structural information provides the fundamental insights for elucidating the reactivity and chemoselectivity in the Sm²⁺-water complex system.

Cooperative Intrinsic Basicity and Hydrogen Bonding Render SmI2 More Azaphilic than Oxophilic

It has been recently shown that SmI₂ is more azaphilic than oxophilic. Density functional theory calculations reveal that coordination of 1-3 molecules of ethylenediamine is more exothermic by up to 10 kcal/mol than coordination of the corresponding number of ethylene glycol molecules. Taking into account also hydrogen bonds between ligands and tetrahydrofuran doubles this preference. The intrinsic affinity parallels the order of basicity. The cooperativity with the hydrogen bonding makes SmI₂ more azaphilic than oxophilic.

Coordination-Induced Bond Weakening

Coordination-induced bond weakening is a phenomenon wherein ligand X-H bond homolysis occurs in concert with the energetically favorable oxidation of a coordinating metal complex. The coupling of these two processes enables thermodynamically favorable proton-coupled electron transfer reductions to form weak bonds upon formal hydrogen atom transfer to substrates. Moreover, systems utilizing coordination-induced bond weakening have been shown to facilitate the dehydrogenation of feedstock molecules including water, ammonia, and primary alcohols under mild conditions. The formation of exceptionally weak substrate X-H bonds via small molecule homolysis is a powerful strategy in synthesis and has been shown to enable nitrogen fixation under mild conditions. Coordination-induced bond weakening has also been identified as an integral process in biophotosynthesis and has promising applications in renewable chemical fuel storage systems. This review presents a discussion of the advances made in the study of coordination-induced bond weakening to date. Because of the broad range of metal and ligand species implicated in coordination-induced bond weakening, each literature report is discussed individually and ordered by the identity of the low-valent metal. We then offer mechanistic insights into the basis of coordination-induced bond weakening and conclude with a discussion of opportunities for further research into the development and applications of coordination-induced bond weakening systems.

Hydration of divalent lanthanides, Sm2+ and Eu2+ : A molecular dynamics study with polarizable AMOEBA force field

The chemistry of divalent lanthanides, Ln2+ , is a growing sub-field of heavy element chemistry owing to new synthetic approaches. However, some theoretical aspects of these unusual cations are currently underdeveloped, especially as they relate to their dynamic properties in solution. In this work, we address the hydration of two of the classical Ln2+ cations, Sm2+ and Eu2+ , using atomic multipole optimized energetic for biomolecular applications (AMOEBA) force fields. These cations have not been parameterized to date with AMOEBA, and few studies are available because of their instability with respect to oxidation in aqueous media. Coordination numbers (CN's) of 8.2 and 8.1 respectively for Sm2+ and Eu2+ , and 8.8 for both Sm3+ and Eu3+ have been obtained and are in good agreement with the few available AIMD and X-ray absorption fine structures studies. The decreased CN of Ln2+ compared with Ln3+ arises from progressive water exchange events that indicates the gradual stabilization of 8-coordinate structures with respect to 9-coordinate geometries. Moreover, the effects of the chloride counter anions on the coordination of Ln2+ cations have been studied at different chloride concentrations in this work. Lastly, water exchange times of Ln2+ cations have been calculated to provide a comprehensive understanding of the behavior of Eu2+ and Sm2+ in aqueous chloride media.

Ammonia Solvation vs Aqueous Solvation of Samarium Diiodide. A Theoretical and Experimental Approach to Understanding Bond Activation Upon Coordination to Sm(II)

Coordination-induced desolvation or ligand displacement by cosolvents and additives is a key feature responsible for the reactivity of Sm(II)-based reagent systems. High-affinity proton donor cosolvents such as water and glycols also demonstrate coordination-induced bond weakening of the O-H bond, facilitating reduction of a broad range of substrates. In the present work, the coordination of ammonia to SmI₂ was examined using Born-Oppenheimer molecular dynamics simulations and mechanistic studies, and the SmI₂-ammonia system is compared to the SmI₂-water system. The coordination number and reactivity of the SmI₂-ammonia solvent system were found to be similar to those of SmI₂-water but exhibited an order of magnitude greater rate of arene reduction by SmI₂-ammonia than by SmI₂-water at the same concentrations of cosolvent. In addition, upon coordination of ammonia to SmI₂, the Sm(II)-ammonia solvate demonstrates one of the largest degrees of N-H bond weakening reported in the literature compared to known low-valent transition metal ammonia complexes.

Dynamics of Water in the Solvation Shell of an Iodate Ion: A Born-Oppenheimer Molecular Dynamics Study

The iodate ion has an anisotropic structure and charge distribution. It has a pyramidal shape with the iodine atom located at the peak of the pyramid. The water molecules interact differently with the positively charged iodine and the negatively charged oxygen atoms of this anion, giving rise to two distinct solvation shells. In the present study, we have performed ab initio Born-Oppenheimer molecular dynamics simulations to investigate the dynamics of water molecules in the iodine and oxygen solvation shells of the iodate ion and compared the behavior with those of the bulk. The dynamics of water is calculated for both the BLYP and the dispersion-corrected BLYP-D3 functionals at room temperature. The dynamics of water in the solvation shells at higher temperatures of 353 and 330 K has also been investigated for the BLYP and BLYP-D3 functionals, respectively. The hydrogen bond dynamics, vibrational spectral diffusion, orientational and translational diffusion, and residence dynamics of water molecules in the two solvation shells are looked at in the current study. The ion-water hydrogen bond dynamics is found to be somewhat faster than that for water-water hydrogen bonds in the bulk, which can be attributed to a ring-like electron distribution on the iodate oxygens. The dynamical trends are connected to the water structure making/breaking properties of the positively charged iodine and negatively charged oxygen sites of the anion. Furthermore, orientational jumps of the iodate ion and also those of surrounding water molecules which are hydrogen bonded to the oxygen atoms of the iodate ion are also investigated. It is found that the nature of these orientational jumps can be different from those reported earlier for planar polyoxyanions such as the nitrate ion.

Coordination of tridentate ligands to SmI2: cooperativity and incremental effect on reduction potential and on reactivity

The effect of coordination of a series of tridentate ligands (TDLs) on various features of SmI2 was determined. The TDLs used in this study were diethylene glycol (OOO), diethanolamine (ONO), 2-(2-Aminoethoxy) ethanol (OON), N-(2-Hydroxyethyl) ethylene diamine (ONN) and glycerol (GLY). Of special interest is the effect of these additives on the reduction potential of SmI2. The cyclic voltammograms of the TDLs with nitrogen at the binding sites display simultaneously several peaks, each corresponding to a different coordination level of SmI2, enabling determination of three equilibrium constants. The results are in concert with electronic spectra of SmI2 complexes with these ligands. The second and third equilibrium constants were found to be larger than the first, demonstrating the cooperativity effect. Moreover, the incremental effect of each moiety on the reduction potential of SmI2 was determined. Regarding reactivity of SmI2, excessive coordination of some ligands is shown to have an adverse effect.

Proton donor effects on the reactivity of SmI2. Experimental and theoretical studies on methanol solvation vs. aqueous solvation

Using both computational and experimental data the SmI₂–MeOH system is directly compared to the SmI₂–H₂O system to uncover the basis for their drastic differences in reactivity.

Ab initio molecular dynamics studies of hydroxide coordination of alkaline earth metals and uranyl

Ab initio molecular dynamics (AIMD) simulations of the Mg2+, Ca2+, Sr2+ and UO22+ ions in either a pure aqueous environment or an environment containing two hydroxide ions have been carried out at the density functional level of theory, employing the generalised gradient approximation via the PBE exchange-correlation functional. Calculated mean M-O bond lengths in the first solvation shell of the aquo systems compared very well to existing experimental and computational literature, with bond lengths well within values measured previously and coordination numbers in line with previously calculated values. When applied to systems containing additional hydroxide ions, the methodology revealed increased bond lengths in all systems. Proton transfer events (PTEs) were recorded and were found to be most prevalent in the strontium hydroxide systems, likely due to the low charge density of the ion and the consequent lack of hydroxide coordination. For all alkaline earths, intrashell PTEs which occurred outside of the first solvation shell were most prevalent. Only three PTEs were identified in the entire simulation data of the uranium dihydroxide system, indicating the clear impact of the increased charge density of the hexavalent uranium ion on the strength of metal-oxygen bonds in aqueous solution. Broadly, systems containing more charge dense ions were found to exhibit fewer PTEs than those containing ions of lower charge density.

Coordination-induced O–H bond weakening in Sm(ii)-water complexes

Coordination of water to low-valent Sm leads to O–H bond-weakening that enables PCET to substrates.

Aqueous solvation of Mg(ii) and Ca(ii): A Born-Oppenheimer molecular dynamics study of microhydrated gas phase clusters

The hydration features of [Mg(H₂O)_n]²⁺ and [Ca(H₂O)_n]²⁺ clusters with n = 3-6, 8, 18, and 27 were studied by means of Born-Oppenheimer molecular dynamics simulations at the B3LYP/6-31+G** level of theory. For both ions, it is energetically more favorable to have all water molecules in the first hydration shell when n ≤ 6, but stable lower coordination average structures with one water molecule not directly interacting with the ion were found for Mg²⁺ at room temperature, showing signatures of proton transfer events for the smaller cation but not for the larger one. A more rigid octahedral-type structure for Mg²⁺ than for Ca²⁺ was observed in all simulations, with no exchange of water molecules to the second hydration shell. Significant thermal effects on the average structure of clusters were found: while static optimizations lead to compact, spherically symmetric hydration geometries, the effects introduced by finite-temperature dynamics yield more prolate configurations. The calculated vibrational spectra are in agreement with infrared spectroscopy results. Previous studies proposed an increase in the coordination number (CN) from six to eight water molecules for [Ca(H₂O)_n]²⁺ clusters when n ≥ 12; however, in agreement with recent measurements of binding energies, no transition to a larger CN was found when n > 8. Moreover, the excellent agreement found between the calculated extended X-ray absorption fine structure spectroscopy spectra for the larger cluster and the experimental data of the aqueous solution supports a CN of six for Ca²⁺.

Aqueous Solvation of SmI3: A Born-Oppenheimer Molecular Dynamics Density Functional Theory Cluster Approach

We report the results of Born-Oppenheimer molecular dynamics (BOMD) simulations on the aqueous solvation of the SmI₃ molecule and of the bare Sm³⁺ cation at room temperature using the cluster microsolvation approach including 37 and 29 water molecules, respectively. The electronic structure calculations were done using the M062X hybrid exchange-correlation functional in conjunction with the 6-31G** basis sets for oxygen and hydrogen. For the iodine and samarium atoms, the Stuttgart-Köln relativistic effective-core potentials were utilized with their associated valence basis sets. When SmI₃ is embedded in the microsolvation environment, we find that substitution of the iodine ions by water molecules around Sm(III) cannot be achieved due to an insufficient number of explicit water molecules to fully solvate the four separate metal and halogen ions. Therefore, we studied the solvation dynamics of the bare Sm³⁺ cation with a 29-water molecule model cluster. Through the Sm-O radial distribution function and the evolution of the Sm-O distances, the present study yields a very tightly bound first rigid Sm(III) solvation shell from 2.3 to 2.9 Å whose integration leads to a coordination number of 9 water molecules and a second softer solvation sphere from 3.9 to 5 Å with 12 water molecules. No water exchange processes were found. The theoretical EXAFS spectrum is in excellent agreement with the experimental spectrum for Sm(III) in liquid water. The strong differences between the solvation patterns of Sm(III) vs Sm(II) are discussed in detail.

On the aqueous solvation of AsO(OH)3vs. As(OH)3. Born–Oppenheimer molecular dynamics density functional theory cluster studies

BOMD simulations were used to reveal the hydration features of As(OH)₃ and (for the first time) AsO(OH)₃ in aqueous solution.

A combined NMR and DFT study of conformational dynamics in lanthanide complexes of macrocyclic DOTA-like ligands

A NMR study of the conformational dynamics of paramagnetic Eu³⁺ complexes provides parameters for isomer interconversion with unprecedented accuracy, while DFT reveals mechanistic details at the molecular level.