A Theoretical Study of the Hydration of Li+ by Monte Carlo Simulations with Refined Ab Initio Based Model Potentials

Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide on the Structure of the Aqueous Solution-Air Interface

We present the results of ab initio molecular dynamics simulations of the solution–air interface of aqueous lithium bromide (LiBr). We find that, in agreement with the experimental data and previous simulation results with empirical polarizable force field models, Br^– anions prefer to accumulate just below the first molecular water layer near the interface, whereas Li⁺ cations remain deeply buried several molecular layers from the interface, even at very high concentration. The separation of ions has a profound effect on the average orientation of water molecules in the vicinity of the interface. We also find that the hydration number of Li⁺ cations in the center of the slab N_{c,Li⁺–H2O} ≈ 4.7 ± 0.3, regardless of the salt concentration. This estimate is consistent with the recent experimental neutron scattering data, confirming that results from nonpolarizable empirical models, which consistently predict tetrahedral coordination of Li⁺ to four solvent molecules, are incorrect. Consequently, disruption of the hydrogen bond network caused by Li⁺ may be overestimated in nonpolarizable empirical models. Overall, our results suggest that empirical models, in particular nonpolarizable models, may not capture all of the properties of the solution–air interface necessary to fully understand the interfacial chemistry.

The hydration structure of the heavy-alkalines Rb+ and Cs+ through molecular dynamics and X-ray absorption spectroscopy: surface clusters and eccentricity

Hydration shells around Rb⁺ and Cs⁺ are not symmetric; the cation and the 1st-shell water mass center are separated by ∼0.4 Å, and this is supported by agreement between the theoretical and experimental EXAFS spectrum.

Specific effects of monovalent counterions on the structural and interfacial properties of dodecyl sulfate monolayers

Shows that NH₄⁺ ions dehydrate the DS⁻ headgroup by displacing hydrogen bonded waters from the interface.

Theoretical study of the aqueous solvation of HgCl2: Monte Carlo simulations using second-order Moller-Plesset-derived flexible polarizable interaction potentials

A study of the solvation of HgCl(2) including ab initio aggregates of up to 24 water molecules and the results of extensive Monte Carlo simulations for the liquid phase using MP2-derived interaction potentials is presented. The interaction potentials are flexible, polarizable, and include non-additive effects. We conclude that a cluster description of the solvation mechanism is limited when compared to the condensed phase. The molecular image derived from the MC simulations is peculiar. It resembles that of a hydrophobic solute, which explains the rather easy passage of this neutral molecule through the cell membrane; however, it also shows an intermittent binding of one, two, or three water molecules to HgCl(2) in the fashion of a hydrophilic solute.

A refined potential for hydroxylamine clusters and the liquid phase

A detailed study including ab initio calculations and classic Monte-Carlo simulations of hydroxylamine in the gas and liquid phases is presented. A classical interaction potential for hydroxylamine, which includes polarizability, many-body effects, and intramolecular relaxation, was constructed. The results of the simulation were compared to the available experimental data in order to validate the model. We conclude that liquid hydroxylamine has a multitude of hydrogen bonds leading to a large density where the existence of cis conformers and clusters of these conformers is possible. This explains the occurrence of the classical [R. Nast and I. Z. Foppl, Z. Anorg. Allg. Chem. 263, 310 (1950)] scheme for the molecule's decomposition at room temperature and its large exothermicity and instability.

The role of Li + , Na + , and K + in the ligand binding inside the human acetylcholinesterase gorge

Alkali cations can affect the catalytic efficiency of enzymes. This is particularly true when dealing with enzymes whose substrate bears a formal positive charge. Computational and biochemical approaches have been combined to shed light on the atomic aspects of the role of Li(+), Na(+), and K(+) on human acetylcholinesterase (hAChE) ligand binding. In this respect, molecular dynamics simulations and our recently developed metadynamics method were applied to study the entrance of the three cations in the gorge of hAChE, and their effect on the dynamical motion of a ligand (tetramethylammonium) from the bulk of the solvent into the deep narrow enzyme gorge. Furthermore, in order to support the theoretical results, K(M) and k(cat) for the acetylcholine hydrolysis in the presence of the three cations were evaluated by using an approach based on the Ellman's method. The combination of computational and biochemical experiments clearly showed that Li(+), Na(+), and K(+) may influence the ligand binding at the hAChE gorge.

Adaptive Partitioning in Combined Quantum Mechanical and Molecular Mechanical Calculations of Potential Energy Functions for Multiscale Simulations

In many applications of multilevel/multiscale methods, an active zone must be modeled by a high-level electronic structure method, while a larger environmental zone can be safely modeled by a lower-level electronic structure method, molecular mechanics, or an analytic potential energy function. In some cases though, the active zone must be redefined as a function of simulation time. Examples include a reactive moiety diffusing through a liquid or solid, a dislocation propagating through a material, or solvent molecules in a second coordination sphere (which is environmental) exchanging with solvent molecules in an active first coordination shell. In this article, we present a procedure for combining the levels smoothly and efficiently in such systems in which atoms or groups of atoms move between high-level and low-level zones. The method dynamically partitions the system into the high-level and low-level zones and, unlike previous algorithms, removes all discontinuities in the potential energy and force whenever atoms or groups of atoms cross boundaries and change zones. The new adaptive partitioning (AP) method is compared to Rode's "hot spot" method and Morokuma's "ONIOM-XS" method that were designed for multilevel molecular dynamics (MD) simulations. MD simulations in the microcanonical ensemble show that the AP method conserves both total energy and momentum, while the ONIOM-XS method fails to conserve total energy and the hot spot method fails to conserve both total energy and momentum. Two versions of the AP method are presented, one scaling as O(2N) and one with linear scaling in N, where N is the number of groups in a buffer zone separating the active high-level zone from the environmental low-level zone. The AP method is also extended to systems with multiple high-level zones to allow, for example, the study of ions and counterions in solution using the multilevel approach.

Ion hydration in nanopores and the molecular basis of selectivity

Using a simple model, it is shown that the cost of constraining a hydrated potassium ion inside a narrow pore is smaller than the cost of constraining hydrated sodium or lithium ions in pores of radius around 1.5 A. The opposite is true for pores of radius around 2.5 A. The reason for the selectivity in the first region is that the potassium ion allows for a greater distortion of its hydration shell and can therefore maintain a better coordination, and the reason for the reverse selectivity in the second region is that the smaller ions retain their hydration shells in these pores. This is relevant to the molecular basis of ion selective channels, and since this mechanism does not depend on the molecular details of the pore, it could also operate in all sorts of nanotubes.