No Arabic abstract
Path-integral ab initio molecular dynamics (PI-AIMD) calculations have been employed to probe the nature of chloride ion solvation in aqueous solution. Nuclear quantum effects (NQEs) are shown to weaken hydrogen bonding between the chloride anion and the solvation shell of water molecules. As a consequence, the disruptive effect of the anion on the solvent water structure is significantly reduced compared to what is found in the absence of NQEs. The chloride hydration structure obtained from PI-AIMD agrees well with information extracted from neutron scattering data. Inparticular, the observed satellite peak in the hydrogen-chloride-hydrogen triple angular distribution serves as a clear signature of NQEs. The present results suggest that NQEs are likely to play acrucial role in determining the structure of saline solutions.
Photoelectron spectroscopy experiments in ionic solutions reveal important electronic structure information, in which the interaction between hydrated ions and water solvent can be inferred. Based on many-body perturbation theory with GW approximation, we theoretically compute the quasiparticle electronic structure of chloride anion solution, which is modeled by path-integral $ab$ $initio$ molecular dynamics simulation by taking account the nuclear quantum effects (NQEs). The electronic levels of hydrated anion as well as water are determined and compared to the recent experimental photoelectron spectra. It is found that NQEs improve the agreement between theoretical prediction and experiment because NQEs effectively weaken the hybridization of the between the $rm Cl^-$ anion and water. Our study indicates that NQEs plays a small but non-negligible role in predicting the electronic structure of the aqueous solvation of ions of the Hofmeister series.
The resolving power of solid-state nuclear magnetic resonance (NMR) crystallography depends heavily on the accuracy of the computational prediction of NMR chemical shieldings of candidate structures, which are usually taken to be local minima in the potential energy surface. To test the limits of this approximation, we perform a systematic study of the role of finite-temperature and quantum nuclear fluctuations on $^1$H, $^{13}$C, and $^{15}$N chemical shieldings in molecular crystals -- considering the paradigmatic examples of the different polymorphs of benzene, glycine, and succinic acid. We find the effect of quantum fluctuations to be comparable in size to the typical errors of predictions of chemical shieldings for static nuclei with respect to experimental measurements, and to improve the match between experiments and theoretical predictions, translating to more reliable assignment of the NMR spectra to the correct candidate structure. Thanks to the use of integrated machine-learning models trained on both first-principles configurational energies and chemical shieldings, the accurate sampling of thermal and quantum fluctuations of the structures can be achieved at an affordable cost, setting a new standard for the calculations that underlie solid-state structural determination by NMR.
We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQE) in molecular simulations: the adaptive Quantum Thermal Bath (adQTB). In this method, zero point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation-dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and we report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.
The stationary nonempirical simulations of Na+(H2O)n clusters with n in a range of 28 to 51 carried out at the density functional level with a hybrid B3LYP functional and the Born-Oppenheimer molecular dynamics modeling of the size selected clusters reveal the interrelated structural and energetic peculiarities of sodium hydration structures. Surface, bulk, and transient configurations of the clusters are distinguished with the different location of the sodium nucleus (close to either the spatial center of the structure or one of its side faces) and its consistently changing coordination number (which typically equals five or six). The <rNaO> mean Na-O distances for the first-shell water molecules are found to depend both on the spatial character of the structure and the local coordination of sodium. The <rNaO> values are compared to different experimental estimates, and the virtual discrepancy of the latter is explained based on the results of the cluster simulations. Different coordination neighborhoods of sodium are predicted depending on its local fraction in the actual specimens.
Some aspects of direct ion transfer across the water/1,2-dichloroethane are analyzed using a very simple model based on thermodynamic considerations. It was concluded that ion solvation by water molecules may occur in some particular cases in the organic phase, delivering an important contribution to the Gibbs free energy of ion transfer between the aqueous and the organic phase. In general terms, this particular type of transfer should be favored in the case of highly charged small ions at interfaces with a relatively low surface tension and a large difference between the reciprocal of the corresponding dielectric constants.