Hydration structure and dynamics of K+ and Ca2+ in aqueous solution: Comparison of conventional QM/MM and ONIOM-XS MD simulations

Hydration Structure of Na+ and K+ Ions in Solution Predicted by Data-Driven Many-Body Potentials

The hydration structure of Na⁺ and K⁺ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na⁺ and K⁺ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.

Study on the Structure of a Mixed KCl and K2SO4 Aqueous Solution Using a Modified X-ray Scattering Device, Raman Spectroscopy, and Molecular Dynamics Simulation

The microstructure of a mixed KCl and K2SO4 aqueous solution was studied using X-ray scattering (XRS), Raman spectroscopy, and molecular dynamics simulation (MD). Reduced structure functions [F(Q)], reduced pair distribution functions [G(r)], Raman spectrum, and pair distribution functions (PDF) were obtained. The XRS results show that the main peak (r = 2.81 Å) of G(r) shifted to the right of the axis (r = 3.15 Å) with increased KCl and decreased K2SO4. The main peak was at r = 3.15 Å when the KCl concentration was 26.00% and the K2SO4 concentration was 0.00%. It is speculated that this phenomenon was caused by the main interaction changing, from K-OW (r = 2.80 Å) and OW-OW (r = 2.80 Å), to Cl−-OW (r = 3.14 Å) and K+-Cl− (r = 3.15 Å). According to the trend of the hydrogen bond structure in the Raman spectrum, when the concentration of KCl was high and K2SO4 was low, the destruction of the tetrahedral hydrogen bond network in the solution was more serious. This shows that the destruction strength of the anion to the hydrogen bond network structure in solution was Cl− > SO42−. In the MD simulations, the coordination number of OW-OW decreased with increasing KCl concentration, indicating that the tetrahedral hydrogen bond network was severely disrupted, which confirmed the results of the Raman spectroscopy. The hydration radius and coordination number of SO42− in the mixed solution were larger than Cl−, thus revealing the reason why the solubility of KCl in water was greater than that of K2SO4 at room temperature.

Structure of Aqueous KNO3 Solutions by Wide-Angle X-ray Scattering and Density Functional Theory

The structure of aqueous KNO₃ solutions was studied by wide-angle X-ray scattering (WAXS) and density functional theory (DFT). The interference functions were subjected to empirical potential structure refinement (EPSR) modeling to extract the detailed hydration structure information. In aqueous KNO₃ solutions with a concentration from 0.50 to 2.72 mol·dm^-3, water molecules coordinate K⁺ with a mean coordination number (CN) from 6.6 ± 0.9 to 5.8 ± 1.2, respectively, and a K-O_W(H₂O) distance of 2.82 Å. To further describe the hydration properties of K⁺, a hydration factor (f_h) was defined based on the orientational angle between the water O-H vector and the ion-oxygen vector. The f_h value obtained for K⁺ is 0.792, 0.787, 0.766, and 0.765, and the corresponding average hydration numbers (HN) are 5.2, 5.1, 4.8, and 4.5. The reduced HN is compensated by the direct binding of oxygen atoms of NO₃^- with the average CN from 0.3 ± 0.7 to 2.6 ± 1.7, respectively, and the K-O_N(NO₃^-) distance of 2.82 Å. The average number of water molecules around NO₃^- slightly reduces from 10.5 ± 1.5 to 9.6 ± 1.7 with r_N-OW = 3.63 Å. K⁺-NO₃^- ion association is characterized by K-O_N and K-N pair correlation functions (PCFs). A K-N peak is observed at 3.81 Å as the main peak with a shoulder at 3.42 Å in g_K-N(r). This finding indicates that two occupancy sites are available for K⁺, i.e., the edge-shared bidentate (BCIP) and the corner-shared monodentate (MCIP) contact ion pairs. The structure and stability of the KNO₃(H₂O)₁₀ cluster were investigated by a DFT method and cross-checked with the results from WAXS.

Structures of ions accommodated in salty ice Ih crystals

Frozen aqueous electrolytes are ubiquitous and involved in various phenomena occurring in the natural environment. Although salts are expelled from ice during freezing of aqueous solutions, minor amounts of the constituent ions are accommodated in the crystal lattice of ice. This phenomenon was associated with the generation of the Workman-Reynolds freezing potential. Molecular simulations also confirmed the ion incorporation in the crystal lattice of ice Ih upon freezing of aqueous electrolytes and identified possible local structures of the ions. However, no experimental information is available on the structure of ions accommodated in the crystal lattice of ice Ih. In this work, we use X-ray absorption fine structure (XAFS) to study the local structures of K⁺ and Cl^- accommodated in ice Ih single crystals. Previous molecular simulations predicted that ions are trapped in the hexagonal cavities of the ice structure or replace two water molecules in the crystal lattice. Four possible configurations are considered and optimized by the calculations using ONIOM (QM/QM/QM). The results are evaluated in terms of the agreement between the experimental XAFS spectra and those simulated from the optimized structures. The spectra are most reasonably interpreted by assuming that K⁺ replaces one water molecule in the ice crystal lattice and is accommodated in a tetrahedral coordination cage. Similarly, Cl^- probably adopts the same configuration, because it explains the coordination number better than other structures, such as that assuming the replacement of two water molecules belonging to the same hexagonal planes.

Ground-State Structures of Hydrated Calcium Ion Clusters From Comprehensive Genetic Algorithm Search

We searched the lowest-energy structures of hydrated calcium ion clusters Ca2+(H2O)n (n = 10-18) in the whole potential energy surface by the comprehensive genetic algorithm (CGA). The lowest-energy structures of Ca2+(H2O)10-12 clusters show that Ca2+ is always surrounded by six H2O molecules in the first shell. The number of first-shell water molecules changes from six to eight at n = 12. In the range of n = 12-18, the number of first-shell water molecules fluctuates between seven and eight, meaning that the cluster could pack the water molecules in the outer shell even though the inner shell is not full. Meanwhile, the number of water molecules in the second shell and the total hydrogen bonds increase with an increase in the cluster size. The distance between Ca2+ and the adjacent water molecules increases, while the average adjacent O-O distance decreases as the cluster size increases, indicating that the interaction between Ca2+ and the adjacent water molecules becomes weaker and the interaction between water molecules becomes stronger. The interaction energy and natural bond orbital results show that the interaction between Ca2+ and the water molecules is mainly derived from the interaction between Ca2+ and the adjacent water molecules. The charge transfer from the lone pair electron orbital of adjacent oxygen atoms to the empty orbital of Ca2+ plays a leading role in the interaction between Ca2+ and water molecules.

Transport behavior of a single Ca(2+), K(+), and Na(+) in a water-filled transmembrane cyclic peptide nanotube

Molecular dynamics simulations have been performed to investigate the transport properties of a single Ca(2+), K(+), and Na(+) in a water-filled transmembrane cyclic peptide nanotube (CPNT). Two transmembrane CPNTs, i.e., 8×(WL)n=4,5/POPE (with uniform lengths but various radii), were applied to clarify the dependence of ionic transport properties on the channel radius. A huge energy barrier keeps Ca(2+) out of the octa-CPNT, while Na(+) and K(+) can be trapped in two CPNTs. The dominant electrostatic interaction of a cation with water molecules leads to a high distribution of channel water around the cation and D-defects in the first and last gaps, and significantly reduces the axial diffusion of channel water. Water-bridged interactions were mostly found between the artificially introduced Ca(2+) and the framework of the octa-CPNT, and direct coordinations with the tube wall mostly occur for K(+) in the octa-CPNT. A cation may drift rapidly or behave lazily in a CPNT. K(+) behaves most actively and can visit the whole deca-CPNT quickly. The first solvation shells of Ca(2+) and Na(+) are basically saturated in two CPNTs, while the hydration of K(+) is incomplete in the octa-CPNT. The solvation structure of Ca(2+) in the octa-CPNT is most stable, while that of K(+) in the deca-CPNT is most labile. Increasing the channel radius induces numerous interchange attempts between the first-shell water molecules of a cation and the ones in the outer region, especially for the K(+) system.