Comparative CHARMM and AMOEBA Simulations of Lanthanide Hydration Energetics and Experimental Aqueous-Solution Structures

The accurate understanding of metal ion hydration in solutions is a prerequisite for predicting stability, reactivity, and solubility. Herein, additive CHARMM force field parameters were developed to enable molecular dynamics simulations of lanthanide (Ln) speciation in water. Quantitatively similar to the much more resource-intensive polarizable AMOEBA potential, the CHARMM simulations reproduce the experimental hydration free energies and correlations in the first shell (Ln-oxygen distance and hydration number). Comparisons of difference pair-distribution functions obtained from the two simulation approaches with those from high-energy X-ray scattering experiments reveal good agreement of first-coordination sphere correlations for the Lu³⁺ ion (CHARMM only), but further improvement to both approaches is required to reproduce the broad, non-Gaussian distribution seen from the La³⁺ experiment. Second-coordination sphere comparisons demonstrate the importance of explicitly including an anion in the simulation. This work describes the usefulness of less resource-intensive additive potentials in some complex chemical systems such as solution environments where multiple interactions have similar energetics. In addition, 3-dimensional descriptions of the La³⁺ and Lu³⁺ coordination geometries are extracted from the CHARMM simulations and generally discussed in terms of potential improvements to solute-structure modeling within solution environments

Driving Force for Water Permeation Across Lipid Membranes

The permeation of water across lipid membranes is of paramount importance in biological and technological processes. The driving force for such energetically unfavorable processes is explored here. To determine the effect of the lipid membrane conformation, water transport in both liquid-crystalline and ordered gel phases is studied in zwitterionic dipalmitoyl phosphatidylcholine (DPPC) bilayers and anionic 1,2-dilauroyl-<i>sn</i>-glycero-3-phosphol-l-serine (DLPS) bilayers via atomistic molecular dynamics simulations. These phases are accessed by changing the temperature in DPPC membranes and by additionally changing the valency of counterions (i.e., Na⁺ and Zn²⁺) in DLPS membranes. The membrane conformation is found to play a critical function in water permeation, regardless of the type of lipid. The fluctuations in the potential energy are found to have a significant, if not the exclusive, role in the transportation of water across lipid membranes

Comparative CHARMM and AMOEBA Simulations of Lanthanide Hydration Energetics and Experimental Aqueous-Solution Structures

The accurate understanding of metal ion hydration in solutions is a prerequisite for predicting stability, reactivity, and solubility. Herein, additive CHARMM force field parameters were developed to enable molecular dynamics simulations of lanthanide (Ln) speciation in water. Quantitatively similar to the much more resource-intensive polarizable AMOEBA potential, the CHARMM simulations reproduce the experimental hydration free energies and correlations in the first shell (Ln-oxygen distance and hydration number). Comparisons of difference pair-distribution functions obtained from the two simulation approaches with those from high-energy X-ray scattering experiments reveal good agreement of first-coordination sphere correlations for the Lu³⁺ ion (CHARMM only), but further improvement to both approaches is required to reproduce the broad, non-Gaussian distribution seen from the La³⁺ experiment. Second-coordination sphere comparisons demonstrate the importance of explicitly including an anion in the simulation. This work describes the usefulness of less resource-intensive additive potentials in some complex chemical systems such as solution environments where multiple interactions have similar energetics. In addition, 3-dimensional descriptions of the La³⁺ and Lu³⁺ coordination geometries are extracted from the CHARMM simulations and generally discussed in terms of potential improvements to solute-structure modeling within solution environments

Non-governmental organizations structures

Anion exchange at positively charged interfaces plays an important role in a variety of physical and chemical processes. However, the molecular-scale details of these processes, especially with heavy and large anionic complexes, are not well-understood. We studied the adsorption of PtCl₆^2– anionic complexes to floating DPTAP monolayers in the presence of excess Cl^– as a function of the bulk chlorometalate concentration. This system aims to simulate the industrial conditions for heavy metal separations with solvent extraction. In situ X-ray scattering and fluorescence measurements, which are element and depth sensitive, show that the chlorometalate ions only adsorb in the diffuse layer at lower concentrations, while they adsorb predominantly in the Stern layer at higher concentrations. The response of DPTAP molecules to the adsorbed ions is determined independently by grazing incidence X-ray diffraction and supports this picture. Molecular dynamics simulations further elucidate the nanoscale structure of the interfacial complexes. The results suggest that ion hydration and ion–ion correlations play a key role in the competitive adsorption process

Ion Transport Mechanisms in Liquid–Liquid Interface

Interfacial liquid–liquid ion transport is of crucial importance to biotechnology and industrial separation processes including nuclear elements and rare earths. A water-in-oil microemulsion is formulated here with density and dimensions amenable to atomistic molecular dynamics simulation, facilitating convergent theoretical and experimental approaches to elucidate interfacial ion transport mechanisms. Lutetium­(III) cations are transported from the 5 nm diameter water pools into the surrounding oil using an extractant (a lipophilic ligand). Changes in ion coordination sphere and interactions between the interfacial components are studied using a combination of synchrotron X-ray scattering, spectroscopy, and atomistic molecular dynamics simulations. Contrary to existing hypotheses, our model system shows no evidence of interfacial extractant monolayers, but rather ions are exchanged through water channels that penetrate the surfactant monolayer and connect to the oil-based extractant. Our results highlight the dynamic nature of the oil–water interface and show that lipophilic ion shuttles need not form flat monolayer structures to facilitate ion transport across the liquid–liquid interface

Molecular Origins of Mesoscale Ordering in a Metalloamphiphile Phase

Controlling the assembly of soft and deformable molecular aggregates into mesoscale structures is essential for understanding and developing a broad range of processes including rare earth extraction and cleaning of water, as well as for developing materials with unique properties. By combined synchrotron small- and wide-angle X-ray scattering with large-scale atomistic molecular dynamics simulations we analyze here a metalloamphiphile–oil solution that organizes on multiple length scales. The molecules associate into aggregates, and aggregates flocculate into meso-ordered phases. Our study demonstrates that dipolar interactions, centered on the amphiphile headgroup, bridge ionic aggregate cores and drive aggregate flocculation. By identifying specific intermolecular interactions that drive mesoscale ordering in solution, we bridge two different length scales that are classically addressed separately. Our results highlight the importance of individual intermolecular interactions in driving mesoscale ordering

How Hydrogen Bonds Affect the Growth of Reverse Micelles around Coordinating Metal Ions

Extensive research on hydrogen bonds (H-bonds) have illustrated their critical role in various biological, chemical and physical processes. Given that existing studies are predominantly performed in aqueous conditions, how H-bonds affect both the structure and function of aggregates in organic phase is poorly understood. Herein, we investigate the role of H-bonds on the hierarchical structure of an aggregating amphiphile-oil solution containing a coordinating metal complex by means of atomistic molecular dynamics simulations and X-ray techniques. For the first time, we show that H-bonds not only stabilize the metal complex in the hydrophobic environment by coordinating between the Eu­(NO₃)₃ outer-sphere and aggregating amphiphiles, but also affect the growth of such reverse micellar aggregates. The formation of swollen, elongated reverse micelles elevates the extraction of metal ions with increased H-bonds under acidic condition. These new insights into H-bonds are of broad interest to nanosynthesis and biological applications, in addition to metal ion separations

How Hydrogen Bonds Affect the Growth of Reverse Micelles around Coordinating Metal Ions

Extensive research on hydrogen bonds (H-bonds) have illustrated their critical role in various biological, chemical and physical processes. Given that existing studies are predominantly performed in aqueous conditions, how H-bonds affect both the structure and function of aggregates in organic phase is poorly understood. Herein, we investigate the role of H-bonds on the hierarchical structure of an aggregating amphiphile-oil solution containing a coordinating metal complex by means of atomistic molecular dynamics simulations and X-ray techniques. For the first time, we show that H-bonds not only stabilize the metal complex in the hydrophobic environment by coordinating between the Eu­(NO₃)₃ outer-sphere and aggregating amphiphiles, but also affect the growth of such reverse micellar aggregates. The formation of swollen, elongated reverse micelles elevates the extraction of metal ions with increased H-bonds under acidic condition. These new insights into H-bonds are of broad interest to nanosynthesis and biological applications, in addition to metal ion separations

Subtle Effects of Aliphatic Alcohol Structure on Water Extraction and Solute Aggregation in Biphasic Water/<i>n</i>‑Dodecane

Organic phase aggregation behavior of 1-octanol and its structural isomer, 2-ethylhexanol, in a biphasic <i>n</i>-dodecane–water system is studied with a combination of physical measurement, small-angle X-ray scattering (SAXS), and atomistic molecular dynamic simulations. Physical properties of the organic phases are probed following their mixing and equilibration with immiscible water phases. Studies reveal that the interfacial tension decreases as a function of increasing alcohol concentration over the solubility range of the alcohol with no evidence for a critical aggregate concentration (cac). An uptake of water into the organic phases is quantified, as a function of alcohol content, by Karl Fischer titrations. The extraction of water into dodecane was further assessed as a function of alcohol concentration via the slope-analysis method sometimes employed in chemical separations. This method provides a qualitative understanding of solute (water/alcohol) aggregation in the organic phase. The physical results are supported by analyses of SAXS data that reveals an emergence of aggregates in <i>n</i>-dodecane at elevated alcohol concentrations. The observed aggregate structure is dependent on the alcohol tail group geometry, consistent with surfactant packing parameter. The formation of these aggregates is discussed at a molecular level, where alcohol–alcohol and alcohol–water H-bonding interactions likely dominate the occurrence and morphology of the aggregates

Aggregation of Heteropolyanions in Aqueous Solutions Exhibiting Short-Range Attractions and Long-Range Repulsions

Charged colloids and proteins in aqueous solutions interact via short-range attractions and long-range repulsions (SALR) and exhibit complex structural phases. These include homogeneously dispersed monomers, percolated monomers, clusters, and percolated clusters. We report the structural architectures of simple charged systems in the form of spherical, Keggin-type heteropolyanions (HPAs) by small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations. Structure factors obtained from the SAXS measurements show that the HPAs interact via SALR. Concentration and temperature dependences of the structure factors for HPAs with −3<i>e</i> (<i>e</i> is the charge of an electron) charge are consistent with a mixture of nonassociated monomers and associated randomly percolated monomers, whereas those for HPAs with −4<i>e</i> and −5<i>e</i> charges exhibit only nonassociated monomers in aqueous solutions. Our experiments show that the increase in magnitude of the charge of the HPAs increases their repulsive interactions and inhibits their aggregation in aqueous solutions. MD simulations were done to reveal the atomistic scale origins of SALR between HPAs. The short-range attractions result from water or proton-mediated hydrogen bonds between neighboring HPAs, whereas the long-range repulsions are due to the distributions of ions surrounding the HPAs