Insights into Uranyl Chemistry from Molecular Dynamics Simulations

An adsorbate biased dynamic 3D porous framework for inverse CO2 sieving over C2H2

Separating carbon dioxide (CO2) from acetylene (C2H2) is one of the most critical and complex industrial separations due to similarities in physicochemical properties and molecular dimensions. Herein, we report a novel Ni-based three-dimensional framework {[Ni4(μ3-OH)2(μ2-OH2)2(1,4-ndc)3](3H2O)}n (1,4-ndc = 1,4-naphthalenedicarboxylate) with a one-dimensional pore channel (3.05 × 3.57 Å2), that perfectly matches with the molecular size of CO2 and C2H2. The dehydrated framework shows structural transformation, decorated with an unsaturated Ni(ii) centre and pendant oxygen atoms. The dynamic nature of the framework is evident by displaying a multistep gate opening type CO2 adsorption at 195, 273, and 298 K, but not for C2H2. The real time breakthrough gas separation experiments reveal a rarely attempted inverse CO2 selectivity over C2H2, attributed to open metal sites with a perfect pore aperture. This is supported by crystallographic analysis, in situ spectroscopic inspection, and selectivity approximations. In situ DRIFTS measurements and DFT-based theoretical calculations confirm CO2 binding sites are coordinatively unsaturated Ni(ii) and carboxylate oxygen atoms, and highlight the influence of multiple adsorption sites.

On the Diffusion of Ionic Liquids in ILs@ZIF-8 Composite Materials: A Density Functional Theory Study

Recently, composite materials consisting of ionic liquids (ILs) and metal-organic frameworks (MOFs) have attracted a great deal of attention due to their fantastic properties. Many theoretical studies have been performed on their special structures and gas separation applications. Yet, the mechanism for the diffusion of ILs inside MOF channels still remains unclear. Here, the DFT calculations (e.g., rigid and relaxed potential energy surface, PES, scan) together with frontier orbital analysis, natural charge analysis, and energy decomposition analysis were performed to investigate the diffusion behavior of a typical IL, [C4mim][PF6], into the ZIF-8 SOD cage. The PES profiles indicate that it is quite difficult for the cation [C4min]+ to diffuse into the cage of ZIF-8 through the pristine pores because of the large imidazole steric hindrance, which results in a large energy barrier of ca. 40 kcal·mol-1 at the least. Interestingly, the PES reveals that a successful diffusion could be obtained by thermal contributions, which enlarge the pore size through swing effects at higher temperatures. For example, both [C4mim]+ and [PF6]- could easily diffuse through the channel of the ZIF-8 SOD cage when the pore size was increased to 6.9 Å. Subsequently, electronic structure analyses reveal that the main interactions between [PF6]- or [C4mim]+ and ZIF-8 are the steric repulsion interactions. Finally, the effects of the amounts of [C4mim][PF6] on the ZIF-8 structures were investigated, and the results show that two pairs of [C4mim][PF6] per SOD cage are the most stable in terms of the interaction between energies and structural changes. With these findings, we propose that the high-temperature technique could be employed during the synthesis of IL@MOF membranes, to enrich their family members and their industrial applications.

Nature and coordination geometry of geologically relevant aqueous Uranium(VI) complexes up to 400 ºC: A review and new data

The structure of the uranyl aqua ion (UO₂²⁺) and a number of its inorganic complexes (specifically, UO₂Cl⁺, UO₂Cl₂⁰, UO₂SO₄⁰, [Formula: see text] , [Formula: see text] and UO₂OH₄^2-) have been characterised using X-Ray absorption spectroscopy/extended X-Ray absorption fine structure (XAS/EXAFS) at temperatures ranging from 25 to 326 ºC. Results of ab initio molecular dynamics (MD) calculations are also reported for uranyl in chloride and sulfate-bearing fluids from 25 to 400 ºC and 600 bar to 20 kilobar (kb). These results are reported alongside a comprehensive review of prior structural characterisation work with particular focus given to EXAFS works to provide a consistent and up-to-date view of the structure of these complexes under conditions relevant to U mobility in ore-forming systems and around high-grade nuclear waste repositories. Regarding reported EXAFS results, average equatorial coordination was found to decrease in uranyl and its sulfate and chloride complexes as temperature rose - the extent of this decrease differed between species and solution compositions but typically resulted in an equatorial coordination number of ∼3-4 at temperatures above 200 ºC. The [Formula: see text] complex was observed at temperatures from 25 to 247 ºC and exhibited no major structural changes over this temperature range. UO₂OH₄^2- exhibited only minor structural changes over a temperature range from 88 to 326 ºC and was suggested to manifest fivefold coordination with four hydroxyl molecules and one water molecule around its equator. Average coordination values derived from fits of the reported EXAFS data were compared to average coordination values calculated using the experimentally derived thermodynamic data for chloride complexes reported by Dargent et al. (2013) and Migdisov et al. (2018b), and for sulfate complexes reported by Alcorn et al. (2019) and Kalintsev et al. (2019). Sulfate EXAFS data were well described by available thermodynamic data, and chloride EXAFS data were described well by the thermodynamic data of Migdisov et al. (2018b), but not by the data of Dargent et al. (2013). The ab initio molecular dynamics calculations confirmed the trends in equatorial coordination observed with EXAFS and were also able to provide an insight into the effect of pressure in equatorial water coordination - for a given temperature, higher pressures appear to lead to a greater number of equatorially bound waters counteracting the temperature effect.

Contribution of Molecular Dynamics in pNMR for the Structural Determination of AnV and AnVI Complexes in Solution

In this study, we propose to use classical molecular dynamics (MD) coupled with ¹H NMR spectroscopy to study the conformations of different actinyl An^VI (An = U, Np, and Pu) and An^V (An = Np) complexes with tetra-ethyl dyglicolamide (TEDGA) ligands in order to have a better representation of such complexes in solution. Molecular dynamics simulations showed its effectiveness in interpreting the experiments by the calculation of geometric factors needed for the determination of magnetic properties of these complexes. We demonstrated that different conformations of the An^V and An^VI complexes with TEDGA exist in solution with different coordination modes, which is experimentally confirmed by ¹H NMR and EXAFS spectroscopies. Furthermore, MD simulations provide additional insights into the structures of complexes in solution since conformations with fast exchanges, which are not accessible from NMR experiments, have been observed by MD simulations.

Hydrolysis of Small Oxo/Hydroxo Molecules Containing High Oxidation State Actinides (Th, Pa, U, Np, Pu): A Computational Study

The energetics of hydrolysis reactions for high oxidation states of oxo/hydroxo monomeric actinide species (Th^IVO₂, Pa^IVO₂, U^IVO₂, Pa^VO₂(OH), U^VO₂(OH), U^VIO₃, Np^VIO₃, Np^VIIO₃(OH), and Pu^VIIO₃(OH)) were calculated at the CCSD(T) level. The first step is the formation of a Lewis acid/base adduct with H₂O (hydration), followed by a proton transfer to form a dihydroxide molecule (hydrolysis); this process is repeated until all oxo groups are hydrolyzed. The physisorption (hydration) for each H₂O addition was predicted to be exothermic, ca. -20 kcal/mol. The hydrolysis products are preferred energetically over the hydration products for the +IV and +V oxidation states. The compounds with An^VI are a turning point in terms of favoring hydration over hydrolysis. For An^VIIO₃(OH), hydration products are preferred, and only two waters can bind; the complete hydrolysis process is now endothermic, and the oxidation state for the An in An(OH)₇ is +VI with two OH groups each having one-half an electron. The natural bond order charges and the reaction energies provide insights into the nature of the hydrolysis/hydration processes. The actinide charges and bond ionicity generally decrease across the period. The ionic character decreases as the oxidation state and coordination number increase so that covalency increases moving to the right in the actinide period.

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.

Molecular Dynamics Simulation Studies on Structure, Dynamics, and Thermodynamics of Uranyl Nitrate Solution at Various Acid Concentrations

The structural and dynamical characteristics of uranyl ions in an aqueous acidic environment are of immense importance in the field of nuclear fuel reprocessing. In view of that, the structural and dynamical behavior of the uranyl ion in water has been investigated by performing molecular dynamics (MD) simulations using different force fields. All the force fields have depicted similar structural and dynamical properties except the free energy of hydration where the Guilbaud-Wipff (GW) model performs well over the others. The calculated density using MD simulations is found to be in excellent agreement with the measured experimental density, which ensures the accuracy of the adopted GW force field. The calculated surface tension and shear viscosity are seen to be increased with uranyl nitrate concentrations. At a higher concentration of about 4.0 mol/L, the supersaturation effect has been captured by an inflection in the plot of surface tension and shear viscosity against concentration because of the solution heterogeneity, which was correlated by an inflection in the scattering intensity observed by performing the dynamic light scattering experiment. The binding mode of nitrate ions with the uranyl ion is found to be concentration-dependent, and at higher concentration, it is predominantly monodentate.

UO22+ structure in solvent extraction phases resolved at molecular and supramolecular scales: a combined molecular dynamics, EXAFS and SWAXS approach

We developed a polarizable force field for unraveling the UO₂²⁺ structure in both aqueous and solvent extraction phases.

Uranyl Arsenate Complexes in Aqueous Solution: Insights from First-Principles Molecular Dynamics Simulations

In this study, the structures and acidity constants (p K_a's) of uranyl arsenate complexes in solutions have been revealed by using the first principle molecular dynamics technique. The results show that uranyl and arsenate form stable complexes with the U/As ratios of 1:1 and 1:2, and the bidentate complexation between U and As is highly favored. Speciation-pH distributions are derived based on free energy and p K_a calculations, which indicate that for the 1:1 species, UO₂(H₂AsO₄)(H₂O)₃⁺ is the major species at pH < 7, while UO₂(HAsO₄)(H₂O)₃⁰ and UO₂(AsO₄)(H₂O)₃^- dominate in acid-to-alkaline and extreme alkaline pH ranges. For the 1:2 species, UO₂(H₂AsO₄)₂(H₂O)⁰ is dominant under acid-to-neutral pH conditions, while UO₂(HAsO₄)(H₂AsO₄)(H₂O)^-, UO₂(HAsO₄)(HAsO₄)(H₂O)^2-, and UO₂(AsO₄)(HAsO₄)(H₂O)^3- become the major forms in the pH range of 7.2-10.7, 10.7-12.1, and >12.1, respectively.

A molecular dynamics investigation of actinyl–ligand speciation in aqueous solution

Actinyl ions (AnO₂ⁿ⁺), the form in which actinides are commonly found in aqueous solution, are important species in the nuclear fuel cycle.

Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics

Dynamic structural effects, NMR ligand chemical shifts, and ¹⁷O NMR quadrupolar relaxation rates are investigated in the series of complexes UO₂²⁺, UO₂(CO₃)₃^4-, and (UO₂)₃(CO₃)₆^6-. Car-Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO₂²⁺ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl ¹⁷O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO₂(CO₃)₃^4- and (UO₂)₃(CO₃)₆^6- complexes show a solvent dependence of the terminal carbonate ¹⁷O and ¹³C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin-orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique ¹³C chemical shifts is correctly reproduced for (UO₂)₃(CO₃)₆^6-. NMR relaxation rate data support large ¹⁷O peak widths, but remain below those noted in the experimental literature. Comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.

X-Ray crystallographic and computational study on uranyl-salophen complexes bearing nitro groups

In the solid state, salophen-UO₂ complexes bearing one, two, or three NO₂ groups lack the pronounced ligand curvature that represents a structural hallmark for this class of compounds. A detailed structural study based on single-crystal X-ray crystallography and computational methods, comprising molecular dynamics, gas-phase Hartree Fock, and DFT calculations, was carried out to investigate the coordination properties of the uranyl cation.

Computational study of Th(4+) and Np(4+) hydration and hydrolysis of Th(4+) from first principles

The aqueous solvation of Th and Np in the IV oxidation state was examined using cluster models generated by Monte Carlo simulations and density functional theory embedded within the COSMO continuum model to approximate the effect of bulk water. Our results suggest that the coordination number (CN) for both Th(IV) and NP(IV) should be 9, in accordance to some experimental and theoretical results from the literature. The structural values for average oxygen-metal distances are within 0.01 Å compared to experimental data, and also within the experimental error. The calculated ΔG _Sol⁰ are in very good agreement with experimental reported values, with deviations at CN = 9 lower than 1% for both Th(IV) and Np(IV). The hydrolysis constants are also in very good agreement with experimental values. Finally, this [corrected] methodology has the advantage of using a GGA functional (BP86) that not only makes the calculations more affordable computationally than hybrid functional or ab initio molecular dynamics simulations (Car-Parrinello) calculations, but also opens the perspective to use resolution of identity (RI) calculations for more extended systems.

Synthesis, structure and bonding of actinide disulphide dications in the gas phase

CASPT2 computations reveal that gas-phase AnS₂²⁺ ions have ground states with triangular geometries and linear thio-actinyl structures are higher in energy, with a difference that increases upon moving from U to Pu.

First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes

Computational chemistry is a valuable aid to complement experimental studies of organometallic systems and their reactivity. It allows probing mechanistic hypotheses and investigating molecular structures, shedding light on the behavior and properties of molecular assemblies at the atomic scale. When approaching a chemical problem, the computational chemist has to decide on the theoretical approach needed to describe electron/nuclear interactions and the composition of the model used to approximate the actual system. Both factors determine the reliability of the modeling study. The community dedicated much effort to developing and improving the performance and accuracy of theoretical approaches for electronic structure calculations, on which the description of (inter)atomic interactions rely. Here, the importance of the model system used in computational studies is highlighted through examples from our recent research focused on organometallic systems and homogeneous catalytic processes. We show how the inclusion of explicit solvent allows the characterization of molecular events that would otherwise not be accessible in reduced model systems (clusters). These include the stabilization of nascent charged fragments via microscopic solvation (notably, hydrogen bonding), transfer of charge (protons) between distant fragments mediated by solvent molecules, and solvent coordination to unsaturated metal centers. Furthermore, when weak interactions are involved, we show how conformational and solvation properties of organometallic complexes are also affected by the explicit inclusion of solvent molecules. Such extended model systems may be treated under periodic boundary conditions, thus removing the cluster/continuum (or vacuum) boundary, and require a statistical mechanics simulation technique to sample the accessible configurational space. First-principles molecular dynamics, in which atomic forces are computed from electronic structure calculations (namely, density functional theory), is certainly the technique of choice to investigate chemical events in solution. This methodology is well established and thanks to advances in both algorithms and computational resources simulation times required for the modeling of chemical events are nowadays accessible, though the computational requirements use to be high. Specific applications reviewed here include mechanistic studies of the Shilov and Wacker processes, speciation in Pd chemistry, hydrogen bonding to metal centers, and the dynamics of agostic interactions.

Theoretical investigation on the mechanism and dynamics of oxo exchange of neptunyl(VI) hydroxide in aqueous solution

Four types of reaction mechanisms for the oxo ligand exchange of monomeric and dimeric neptunyl(VI) hydroxide in aqueous solution were explored computationally using density functional theory (DFT) and ab initio classical molecular dynamics. The obtained results were compared with previous studies on the oxo exchange of uranyl hydroxide, as well as with experiments. It is found that the stable T-shaped [NpO3(OH)3](3-) intermediate is a key species for oxo exchange in the proton transfer in mononuclear Path I and binuclear Path IV, similar to the case of uranyl(VI) hydroxide. Path I is thought to be the preferred oxo exchange mechanism for neptunyl(VI) hydroxide in our calculations, due to the lower activation energy (22.7 and 13.1 kcal mol(-1) for ΔG(‡) and ΔH(‡), respectively) of the overall reaction. Path II via a cis-neptunyl structure assisted by a water molecule might be a competitive channel against Path I with a mononuclear mechanism, owing to a rapid dynamical process occurring in Path II. In Path IV with the binuclear mechanism, oxo exchange is accomplished via the interaction between [NpO2(OH)4](2-) and T-shaped [NpO3(OH)3](3-) with a low activation energy for the rate-determining step, however, the overall energy required to fulfill the reaction is slightly higher than that in mononuclear Path I, suggesting a possible binuclear process in the higher energy region. The chemical bonding evolution along the reaction pathways was discussed by using topological methodologies of the electron localization function (ELF).

Molecular dynamics simulation study of distribution and dynamics of aqueous solutions of uranyl ions: the effect of varying temperature and concentration

Investigating the characteristics of actinyl ions has been of great interest due to their direct relevance in the nuclear fuel cycle. All-atom molecular dynamics simulations have been employed to study the orientational structure and dynamics of aqueous solutions of uranyl ions of various concentrations. The orientational structure of water around a uranyl ion has been thoroughly investigated by calculating different orientational probability distributions corresponding to different molecular axes of water. The orientational distribution of water molecules in the first coordination shell of a uranyl ion is found to be markedly different from that in bulk water. Analysis of counterion distribution around the uranyl ion reveals the presence of nitrate ions along with water molecules in the first solvation shell. From the comparison of the number of coordinated water and nitrate ions at various uranyl nitrate concentrations, it is evident that these two species compete for occupying the first solvation shell of the uranyl ion. Orientational dynamics of water molecules about different molecular axes of water in the vicinity of uranyl ions have also been investigated and decreasing orientational mobility of water with increasing uranyl concentration has been found. However, it is observed that the orientational dynamics remains more or less the same whether we consider all the water molecules in the aqueous solution or only the solvation shell water molecules. The effect of temperature on the translational and orientational characteristics of the aqueous uranyl solutions has also been studied in detail.

Uranyl Solvation by a Three-Dimensional Reference Interaction Site Model

We report an implementation of the three-dimensional reference interaction site model (3D RISM) that in particular addresses the treatment of the long-range Coulomb field of charged species, represented by point charges and/or a distributed charge density. A comparison of 1D and 3D results for atomic ions demonstrates a reasonable accuracy, even for a moderate size of the unit cell and a moderate grid resolution. In an application to uranyl complexes with 4-6 explicit aqua ligands and an implicit bulk solvent modeled by RISM, we show that the 3D technique is not susceptible to the deficiencies of the 1D technique exposed in our previous work [Li, Matveev, Krüger, Rösch, Comp. Theor. Chem. 2015, 1051, 151]. The 3D method eliminates the artificial superposition of explicit aqua ligands and the RISM medium and predicts essentially the same values for uranyl and uranyl-water bond lengths as a state-of-the-art polarizable continuum model. With the first solvation shell treated explicitly, the observables are nearly independent of the order of the closure relationship used when solving the set of integral equations for the various distribution functions. Furthermore, we calculated the activation barrier of water exchange with a hybrid approach that combines the 3D RISM model for the bulk aqueous solvent and a quantum mechanical description (at the level of electronic density functional theory) of uranyl interacting with explicitly represented water molecules. The calculated result agrees very well with experiment and the best theoretical estimates.

[UO2(NH3)5]Br2·NH3: synthesis, crystal structure, and speciation in liquid ammonia solution by first-principles molecular dynamics simulations

X-Ray diffraction and Car–Parrinello molecular dynamics simulations furnish insights into the speciation of uranyl(vi) in liquid ammonia, calling special attention to the effect of solvation on the U–N bond length and bond strength.

Promising density functional theory methods for predicting the structures of uranyl complexes

By examining the overall accuracy of different theoretical methods in predicting the U–X bond distances (of a series uranyl complexes), we found that both the global-hybrid meta-GGA functional of BB1K and the range-seperated LC-BLYP functional are fairly good (even better than the popular B3LYP method).

Coordination and hydrolysis of plutonium ions in aqueous solution using Car-Parrinello molecular dynamics free energy simulations

Car-Parrinello molecular dynamics (CPMD) simulations have been used to examine the hydration structures, coordination energetics, and the first hydrolysis constants of Pu(3+), Pu(4+), PuO2(+), and PuO2(2+) ions in aqueous solution at 300 K. The coordination numbers and structural properties of the first shell of these ions are in good agreement with available experimental estimates. The hexavalent PuO2(2+) species is coordinated to five aquo ligands while the pentavalent PuO2(+) complex is coordinated to four aquo ligands. The Pu(3+) and Pu(4+) ions are both coordinated to eight water molecules. The first hydrolysis constants obtained for Pu(3+) and PuO2(2+) are 6.65 and 5.70, respectively, all within 0.3 pH unit of the experimental values (6.90 and 5.50, respectively). The hydrolysis constant of Pu(4+), 0.17, disagrees with the value of -0.60 in the most recent update of the Nuclear Energy Agency Thermochemical Database (NEA-TDB) but supports recent experimental findings. The hydrolysis constant of PuO2(+), 9.51, supports the experimental results of Bennett et al. [Radiochim. Acta 1992, 56, 15]. A correlation between the pKa of the first hydrolysis reaction and the effective charge of the plutonium center was found.

Direct calculation of 1-octanol-water partition coefficients from adaptive biasing force molecular dynamics simulations

The 1-octanol-water partition coefficient log K(ow) of a solute is a key parameter used in the prediction of a wide variety of complex phenomena such as drug availability and bioaccumulation potential of trace contaminants. In this work, adaptive biasing force molecular dynamics simulations are used to determine absolute free energies of hydration, solvation, and 1-octanol-water partition coefficients for n-alkanes from methane to octane. Two approaches are evaluated; the direct transfer of the solute from 1-octanol to water phase, and separate transfers of the solute from the water or 1-octanol phase to vacuum, with both methods yielding statistically indistinguishable results. Calculations performed with the TIP4P and SPC∕E water models and the TraPPE united-atom force field for n-alkanes show that the choice of water model has a negligible effect on predicted free energies of transfer and partition coefficients for n-alkanes. A comparison of calculations using wet and dry octanol phases shows that the predictions for log K(ow) using wet octanol are 0.2-0.4 log units lower than for dry octanol, although this is within the statistical uncertainty of the calculation.

Force field development for actinyl ions via quantum mechanical calculations: an approach to account for many body solvation effects

Advances in computational algorithms and methodologies make it possible to use highly accurate quantum mechanical calculations to develop force fields (pair-wise additive intermolecular potentials) for condensed phase simulations. Despite these advances, this approach faces numerous hurdles for the case of actinyl ions, AcO2(n+) (high-oxidation-state actinide dioxo cations), mainly due to the complex electronic structure resulting from an interplay of s, p, d, and f valence orbitals. Traditional methods use a pair of molecules (“dimer”) to generate a potential energy surface (PES) for force field parametrization based on the assumption that many body polarization effects are negligible. We show that this is a poor approximation for aqueous phase uranyl ions and present an alternative approach for the development of actinyl ion force fields that includes important many body solvation effects. Force fields are developed for the UO2(2+) ion with the SPC/Fw, TIP3P, TIP4P, and TIP5P water models and are validated by carrying out detailed molecular simulations on the uranyl aqua ion, one of the most characterized actinide systems. It is shown that the force fields faithfully reproduce available experimental structural data and hydration free energies. Failure to account for solvation effects when generating PES leads to overbinding between UO2(2+) and water, resulting in incorrect hydration free energies and coordination numbers. A detailed analysis of arrangement of water molecules in the first and second solvation shell of UO2(2+) is presented. The use of a simple functional form involving the sum of Lennard-Jones + Coulomb potentials makes the new force field compatible with a large number of available molecular simulation engines and common force fields.

Bromide complexation by the Eu(III) lanthanide cation in dry and humid ionic liquids: a molecular dynamics PMF study

We report a molecular dynamics study on the EuBr(n)(3-n) complexes (n=0 to 6) formed upon complexation of Br(-) by Eu(3+) in the [BMI][PF(6)], [BMI][Tf(2)N] and [MeBu(3)N][Tf(2)N] ionic liquids (ILs), to compare the effect of the IL anion (PF(6)(-) versus Tf(2)N(-)), the IL cation (BMI(+) versus MeBu(3)N(+)) and the "IL humidity" on their solvation and stability. In "dry" solutions all complexes remain stable and the first coordination shell of Eu(3+) is purely anionic (Br(-) and IL anions), surrounded by IL cations (BMI(+) or MeBu(3)N(+) ions). Long range "onion type" solvation features (up to 20 Å from Eu(3+)), with alternating cation-rich and anion-rich solvent shells, are observed around the different complexes. The comparison of gas phase-optimized structures of EuBr(n)(3-n) complexes (that are unstable for n=5 and 6) with those observed in solution points to the importance of solvation forces on the nature of the complex, with a higher stabilization by imidazolium- than by ammonium-based dry ILs. Adding water to the IL has different effects, depending on the IL. In the highly hygroscopic [BMI][PF(6)] IL, Br(-) ligands are displaced by water, to finally form Eu(H(2)O)(9)(3+). In the less "humid" [BMI][Tf(2)N], the EuBr(n)(3-n) complexes do not dissociate and coordinate at most 1-2 H(2)O molecules. We also calculated the free-energy profiles (Potential of Mean Force calculations) for the stepwise complexation of Br(-), and found significant solvent effects. EuBr(6)(3-) is predicted to form in both [BMI][PF(6)] and [BMI][Tf(2)N], but not in [MeBu(3)N][Tf(2)N], mainly due to weaker interactions with the cationic solvation shell. First steps are found to be more exergonic in the PF(6)(-)- than in the Tf(2)N(-)-based IL. Molecular dynamics (MD) comparisons between ILs and classical solvents (acetonitrile and water) are also reported, affording good agreement with the experimental observations of Br(-) complexation by trivalent lanthanides in these classical solvents.