Solvation structure and dynamics of K+ in aqueous ammonia solution: Insights from an ONIOM-XS MD simulation

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.

Adaptive-Partitioning Multilayer Dynamics Simulations: 1. On-the-Fly Switch between Two Quantum Levels of Theory

We propose to generalize the previously developed two-layer permuted adaptive-partitioning quantum-mechanics/molecular-mechanics (QM/MM), which reclassifies atoms as QM or MM on-the-fly in dynamics simulations, to multilayer adaptive-partitioning algorithms that enable multiple levels of theory. In this work, we formulate two new algorithms that smoothly interpolate the energy between two QM (Q1 and Q2) levels of theory. The first "permuted adaptive-partitioning" scheme is based on the weighted many-body expansion of the potential, as in the adaptive-partitioning QM/MM. Unconventional and potentially more efficient, the second "interpolated adaptive-partitioning" method employs alchemical QM calculations with Q1/Q2-mixed basis sets, Fock matrices, and overlap matrices. To our knowledge, this is the first time that such alchemical calculations are performed in QM, although they are routinely done in MM. Test calculations on water-cluster models show that both new algorithms indeed yield smooth energy curves when water molecules shift between Q1 and Q2.