Dielectric Constant of Liquid Water Determined with Neural Network Quantum Molecular Dynamics

Self-learning path integral hybrid Monte Carlo with mixed ab initio and machine learning potentials for modeling nuclear quantum effects in water

The introduction of machine learned potentials (MLPs) has greatly expanded the space available for studying Nuclear Quantum Effects computationally with ab initio path integral (PI) accuracy, with the MLPs' promise of an accuracy comparable to that of ab initio at a fraction of the cost. One of the challenges in development of MLPs is the need for a large and diverse training set calculated by ab initio methods. This dataset should ideally cover the entire phase space, while not searching this space using ab initio methods, as this would be counterproductive and generally intractable with respect to computational time. In this paper, we present the self-learning PI hybrid Monte Carlo Method using a mixed ab initio and ML potential (SL-PIHMC-MIX), where the mixed potential allows for the study of larger systems and the extension of the original SL-HMC method [Nagai et al., Phys. Rev. B 102, 041124 (2020)] to PI methods and larger systems. While the MLPs generated by this method can be directly applied to run long-time ML-PIMD simulations, we demonstrate that using PIHMC-MIX with the trained MLPs allows for an exact reproduction of the structure obtained from ab initio PIMD. Specifically, we find that the PIHMC-MIX simulations require only 5000 evaluations of the 32-bead structure, compared to the 100 000 evaluations needed for the ab initio PIMD result.

How to Assess and Predict Electrical Double Layer Properties. Implications for Electrocatalysis

The electrical double layer (EDL) plays a central role in electrochemical energy systems, impacting charge transfer mechanisms and reaction rates. The fundamental importance of the EDL in interfacial electrochemistry has motivated researchers to develop theoretical and experimental approaches to assess EDL properties. In this contribution, we review recent progress in evaluating EDL characteristics such as the double-layer capacitance, highlighting some discrepancies between theory and experiment and discussing strategies for their reconciliation. We further discuss the merits and challenges of various experimental techniques and theoretical approaches having important implications for aqueous electrocatalysis. A strong emphasis is placed on the substantial impact of the electrode composition and structure and the electrolyte chemistry on the double-layer properties. In addition, we review the effects of temperature and pressure and compare solid-liquid interfaces to solid-solid interfaces.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

Nonlinear Effects of Hydrophobic Confinement on the Electronic Structure and Dielectric Response of Water

Fundamental studies of the dielectrics of confined water are critical to understand the ion transport across biological and synthetic nanochannels. The relevance of these fundamental studies, however, surmounts the difficulty of probing water's dielectric constant as a function of a fine variation in confinement. In this work, we explore the computational efficiency of machine learning potentials to derive the confinement effects on the dielectric constant, polarization, and dipole moment of water. Our simulations predict an enhancement of the axial dielectric constant of water under extreme confinement, arising from either the formation of ferroelectric structures of ordered water or larger dipole fluctuations facilitated by the disruption of water's H-bond network. Our study highlights the impact of hydrophobic nanoconfinement on the dielectric constant and on the ionic and electronic structure of water molecules, pointing to the importance of geometric flexibility and electronic polarizability to properly model confinement effects on water.

Learning Electronic Polarizations in Aqueous Systems

The polarization of periodically repeating systems is a discontinuous function of the atomic positions, a fact which seems at first to stymie attempts at their statistical learning. Two approaches to build models for bulk polarizations are compared: one in which a simple point charge model is used to preprocess the raw polarization to give a learning target that is a smooth function of atomic positions and the total polarization is learned as a sum of atom-centered dipoles and one in which instead the average position of Wannier centers around atoms is predicted. For a range of bulk aqueous systems, both of these methods perform perform comparatively well, with the former being slightly better but often requiring an extra effort to find a suitable point charge model. As a challenging test, we also analyze the performance of the models at the air-water interface. In this case, while the Wannier center approach delivers accurate predictions without further modifications, the preprocessing method requires augmentation with information from isolated water molecules to reach similar accuracy. Finally, we present a simple protocol to preprocess the polarizations in a data-driven way using a small number of derivatives calculated at a much lower level of theory, thus overcoming the need to find point charge models without appreciably increasing the computation cost. We believe that the training strategies presented here help the construction of accurate polarization models required for the study of the dielectric properties of realistic complex bulk systems and interfaces with ab initio accuracy.

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.

Tensorial Properties via the Neuroevolution Potential Framework: Fast Simulation of Infrared and Raman Spectra

Infrared and Raman spectroscopy are widely used for the characterization of gases, liquids, and solids, as the spectra contain a wealth of information concerning, in particular, the dynamics of these systems. Atomic scale simulations can be used to predict such spectra but are often severely limited due to high computational cost or the need for strong approximations that limit the application range and reliability. Here, we introduce a machine learning (ML) accelerated approach that addresses these shortcomings and provides a significant performance boost in terms of data and computational efficiency compared with earlier ML schemes. To this end, we generalize the neuroevolution potential approach to enable the prediction of rank one and two tensors to obtain the tensorial neuroevolution potential (TNEP) scheme. We apply the resulting framework to construct models for the dipole moment, polarizability, and susceptibility of molecules, liquids, and solids and show that our approach compares favorably with several ML models from the literature with respect to accuracy and computational efficiency. Finally, we demonstrate the application of the TNEP approach to the prediction of infrared and Raman spectra of liquid water, a molecule (PTAF-), and a prototypical perovskite with strong anharmonicity (BaZrO3). The TNEP approach is implemented in the free and open source software package gpumd, which makes this methodology readily available to the scientific community.

Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+ Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Ruru Ma Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Nitish Baradwaj Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Ken-Ichi Nomura Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Aravind Krishnamoorthy Department of Mechanical Engineering, Texas A&M University, College Station, Texas 77843-3123, USA
Rajiv K Kalia Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Aiichiro Nakano Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Priya Vashishta Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA

Collapse

Chromium removal from chromium gypsum through microwave hydrothermal crystal phase regulation

Chromium gypsum (CG) is a common hazardous waste formed in chromium salt or electroplating industries. The trapped or lattice-doped CrO42- in gypsum crystals are difficult to be reduced or removed by traditional methods, which will be re-oxidized or slowly released during long-term hypaethral storage. In this study, microwave hydrothermal treatment was applied to remove chromium in CG. Under optimal conditions (solid-liquid ratio of 1:5, 0.1 M sulfuric acid as liquid media, and 110 °C), over 99% of the chromium in CG can be removed within 10 min. XRD spectra indicated that 59.8% gypsum was transformed to from dihydrate gypsum to hemihydrate gypsum. The toxicity leaching test shows that chromium in CG is 377.0 mg/L before detoxification and 0.55 mg/L after detoxification, which proves that chromium in CG lattice can be efficiently removed. This work enables to significantly advance the dehydration phase transformation process of gypsum and release the heavy metal impurities within it more quickly and provides new possibilities to treat similar solid waste containing gypsum or minerals with hydration water.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

Wrinkles, Ridges, Miura-Ori, and Moiré Patterns in MoSe2 Using Neural Networks

Effects of lateral compression on out-of-plane deformation of two-dimensional MoSe₂ layers are investigated. A MoSe₂ monolayer develops periodic wrinkles under uniaxial compression and Miura-Ori patterns under biaxial compression. When a flat MoSe₂ monolayer is placed on top of a wrinkled MoSe₂ layer, the van der Waals (vdW) interaction transforms wrinkles into ridges and generates mixed 2H and 1T phases and chain-like defects. Under a biaxial strain, the vdW interaction induces regions of Miura-Ori patterns in bilayers. Strained systems analyzed using a convolutional neural network show that the compressed system consists of semiconducting 2H and metallic 1T phases. The energetics, mechanical response, defect structure, and dynamics are analyzed as bilayers undergo wrinkle-ridge transformations under uniaxial compression and moiré transformations under biaxial compression. Our results indicate that in-plane compression can induce self-assembly of out-of-plane metasurfaces with controllable semiconducting and metallic phases and moiré patterns with unique optoelectronic properties.

Local Molecular Field Theory for Coulomb Interactions in Aqueous Solutions

Coulomb interactions play a crucial role in a wide array of processes in aqueous solutions but present conceptual and computational challenges to both theory and simulations. We review recent developments in an approach addressing these challenges─local molecular field (LMF) theory. LMF theory exploits an exact and physically suggestive separation of intermolecular Coulomb interactions into strong short-range and uniformly slowly varying long-range components. This allows us to accurately determine the averaged effects of the long-range components on the short-range structure using effective single particle fields and analytical corrections, greatly reducing the need for complex lattice summation techniques used in most standard approaches. The simplest use of these ideas in aqueous solutions leads to the short solvent (SS) model, where both solvent-solvent and solute-solvent Coulomb interactions have only short-range components. Here we use the SS model to give a simple description of pairing of nucleobases and biologically relevant ions in water.

Squishing Skyrmions: Symmetry-Guided Dynamic Transformation of Polar Topologies Under Compression

Mechanical controllability of recently discovered topological defects (e.g., skyrmions) in ferroelectric materials is of interest for the development of ultralow-power mechano-electronics that are protected against thermal noise. However, fundamental understanding is hindered by the "multiscale quantum challenge" to describe topological switching encompassing large spatiotemporal scales with quantum mechanical accuracy. Here, we overcome this challenge by developing a machine-learning-based multiscale simulation framework─a hybrid neural network quantum molecular dynamics (NNQMD) and molecular mechanics (MM) method. For nanostructures composed of SrTiO₃ and PbTiO₃, we find how the symmetry of mechanical loading essentially controls polar topological switching. We find under symmetry-breaking uniaxial compression a squishing-to-annihilation pathway versus formation of a topological composite named skyrmionium under symmetry-preserving isotropic compression. The distinct pathways are explained in terms of the underlying materials' elasticity and symmetry, as well as the Landau-Lifshitz-Kittel scaling law. Such rational control of ferroelectric topologies will likely facilitate exploration of the rich ferroelectric "topotronics" design space.

Dielectric dispersion characteristics of the phospholipid bilayer with subnanometer resolution from terahertz to mid-infrared

There is growing interest in whether the myelinated nerve fiber acts as a dielectric waveguide to propagate terahertz to mid-infrared electromagnetic waves, which are presumed stable signal carrier for neurotransmission. The myelin sheath is formed as a multilamellar biomembrane structure, hence insights into the dielectric properties of the phospholipid bilayer is essential for a complete understanding of the myelinated fiber functioning. In this work, by means of atomistic molecular dynamics simulations of the dimyristoylphosphatidylcholine (DMPC) bilayer in water and numerical calculations of carefully layered molecules along with calibration of optical dielectric constants, we for the first time demonstrate the spatially resolved (in sub-nm) dielectric spectrum of the phospholipid bilayer in a remarkably wide range from terahertz to mid-infrared. More specifically, the membrane head regions exhibit both larger real and imaginary permittivities than that of the tail counterparts in the majority of the 1–100 THz band. In addition, the spatial variation of dielectric properties suggests advantageous propagation characteristics of the phospholipid bilayer in a relatively wide band of 55–85 THz, where the electromagnetic waves are well confined within the head regions.

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures (e.g., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO3 nanowires embedded in SrTiO3. We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Exploring far-from-equilibrium ultrafast polarization control in ferroelectric oxides with excited-state neural network quantum molecular dynamics

Ferroelectric materials exhibit a rich range of complex polar topologies, but their study under far-from-equilibrium optical excitation has been largely unexplored because of the difficulty in modeling the multiple spatiotemporal scales involved quantum-mechanically. To study optical excitation at spatiotemporal scales where these topologies emerge, we have performed multiscale excited-state neural network quantum molecular dynamics simulations that integrate quantum-mechanical description of electronic excitation and billion-atom machine learning molecular dynamics to describe ultrafast polarization control in an archetypal ferroelectric oxide, lead titanate. Far-from-equilibrium quantum simulations reveal a marked photo-induced change in the electronic energy landscape and resulting cross-over from ferroelectric to octahedral tilting topological dynamics within picoseconds. The coupling and frustration of these dynamics, in turn, create topological defects in the form of polar strings. The demonstrated nexus of multiscale quantum simulation and machine learning will boost not only the emerging field of ferroelectric topotronics but also broader optoelectronic applications.

Self-consistent determination of long-range electrostatics in neural network potentials

Machine learning has the potential to revolutionize the field of molecular simulation through the development of efficient and accurate models of interatomic interactions. Neural networks can model interactions with the accuracy of quantum mechanics-based calculations, but with a fraction of the cost, enabling simulations of large systems over long timescales. However, implicit in the construction of neural network potentials is an assumption of locality, wherein atomic arrangements on the nanometer-scale are used to learn interatomic interactions. Because of this assumption, the resulting neural network models cannot describe long-range interactions that play critical roles in dielectric screening and chemical reactivity. Here, we address this issue by introducing the self-consistent field neural network - a general approach for learning the long-range response of molecular systems in neural network potentials that relies on a physically meaningful separation of the interatomic interactions - and demonstrate its utility by modeling liquid water with and without applied fields.