The microscopic origin of the anomalous isotopic properties of ice relies on the strong quantum anharmonic regime of atomic vibration

Neutron scattering and neural-network quantum molecular dynamics investigation of the vibrations of ammonia along the solid-to-liquid transition

Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs.

Anharmonic Coupling of Stretching Vibrations in Ice: A Periodic VSCF and VCI Description

The anharmonicity of O-H stretching vibrations of water ice is characterized by use of a periodic implementation of the vibrational self-consistent field (VSCF) and vibrational configuration interaction (VCI) methods, which take phonon-phonon couplings explicitly into account through numerical evaluation of high-order terms of the nuclear potential. The low-temperature, proton-ordered phase of water ice (namely, ice XI) is investigated. The net effect of a coupled anharmonic treatment of stretching modes is not just a rigid blue-shift of the respective harmonic spectral frequencies but rather a complex change of their relative spectral positions, which cannot be captured by simple scaling strategies based on harmonic calculations. The adopted techniques allow for a hierarchical treatment of anharmonic terms of the nuclear potential, which is key to an effective identification of leading factors. We show that the anharmonic independent-mode approximation─only describing the "intrinsic anharmonicity" of the O-H stretches─is unable to capture the correct physics, and that couplings among O-H stretches must be described. Inspection of harmonic normal coordinates allows identification of specific features of the O-H stretching motions which most likely enable strong mode-mode couplings. Finally, by coupling O-H stretches to all other possible modes of ice XI (THz collective vibrations, molecular librations, bendings), we identify specific types of motion which significantly affect O-H stretching states: in particular, molecular librations are found to affect the stretching states more than molecular bendings.

Understanding anharmonic effects on hydrogen desorption characteristics of MgnH2n nanoclusters by ab initio trained deep neural network

Magnesium hydride (MgH₂) has been widely studied for effective hydrogen storage. However, its bulk desorption temperature (553 K) is deemed too high for practical applications. Besides doping, a strategy to decrease such reaction energy for releasing hydrogen is the use of MgH₂-based nanoparticles (NPs). Here, we investigate first the thermodynamic properties of Mg_nH_2n NPs (n < 10) from first-principles, in particular by assessing the anharmonic effects on the enthalpy, entropy and thermal expansion by means of the stochastic self consistent harmonic approximation (SSCHA). This method goes beyond previous approaches, typically based on molecular mechanics and the quasi-harmonic approximation, allowing the ab initio calculation of the fully-anharmonic free energy. We find an almost linear dependence on temperature of the interatomic bond lengths - with a relative variation of few percent over 300 K - alongside with a bond distance decrease of the Mg-H bonds. In order to increase the size of Mg_nH_2n NPs toward experiments of hydrogen desorption we devise a computationally effective machine learning model trained to accurately determine the forces and total energies (i.e. the potential energy surfaces), integrating the latter with the SSCHA model to fully include the anharmonic effects. We find a significative decrease of the H-desorption temperature for sub-nanometric clusters Mg_nH_2n with n ≤ 10, with a non-negligible, although little effect due to anharmonicities (up to 10%).