Search

In molecular simulations, neural network force fields aim at achieving ab initio accuracy with reduced computational cost. This work introduces enhancements to the Deep Potential network architecture, integrating a message-passing framework and a new lightweight implementation with various improvements. Our model achieves accuracy on par with leading machine learning force fields and offers significant speed advantages, making it well-suited for large-scale, accuracy-sensitive systems. We also introduce a new iterative model for Wannier center prediction, allowing us to keep track of electron positions in simulations of general insulating systems. We apply our model to study the solvated electron in bulk water, an ostensibly simple system that is actually quite challenging to represent with neural networks. Our trained model is not only accurate, but can also transfer to larger systems. Our simulation confirms the cavity model, where the electron's localized state is observed to be stable. Through an extensive run, we accurately determine various structural and dynamical properties of the solvated electron.

We introduce a machine learning-based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori-Zwanzig formalism, under the constraint of the fluctuation-dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency., Competing Interests: Competing interests statement:The authors declare no competing interest.

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water under deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to a vacuum, liquid water, and ice. We apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to a vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in a vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30 000 atoms. Future work will be directed towards the calculation of nucleation free-energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces in ice nucleation.

West Nile virus (WNV), a flavivirus transmitted by mosquito bites, causes primarily mild symptoms but can also be fatal. Therefore, predicting and controlling the spread of West Nile virus is essential for public health in endemic areas. We hypothesized that socioeconomic factors may influence human risk from WNV. We analyzed a list of weather, land use, mosquito surveillance, and socioeconomic variables for predicting WNV cases in 1-km hexagonal grids across the Chicago metropolitan area. We used a two-stage lightGBM approach to perform the analysis and found that hexagons with incomes above and below the median are influenced by the same top characteristics. We found that weather factors and mosquito infection rates were the strongest common factors. Land use and socioeconomic variables had relatively small contributions in predicting WNV cases. The Light GBM handles unbalanced data sets well and provides meaningful predictions of the risk of epidemic disease outbreaks., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Wan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The chemical equilibrium between self-ionized and molecular water dictates the acid-base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)-based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid-base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance [Formula: see text] is found to converge for simulation cells containing more than 500 H 2 O, and a distance of [Formula: see text] 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

The MB-pol many-body potential accurately predicts many properties of water, including cluster, liquid phase, and vapor-liquid equilibrium properties, but its high computational cost can make applying it in large-scale simulations quite challenging. In order to address this limitation, we developed a "deep potential" neural network (DPMD) model based on the MB-pol potential for water. We find that a DPMD model trained on mostly liquid configurations yields a good description of the bulk liquid phase but severely underpredicts vapor-liquid coexistence densities. By contrast, adding cluster configurations to the neural network training set leads to a good agreement for the vapor coexistence densities. Liquid phase densities under supercooled conditions are also represented well, even though they were not included in the training set. These results confirm that neural network models can combine accuracy and transferability if sufficient attention is given to the construction of a representative training set for the target system.

Understanding the condensed-phase behavior of chiral molecules is important in biology as well as in a range of technological applications, such as the manufacture of pharmaceuticals. Here, we use molecular dynamics simulations to study a chiral four-site molecular model that exhibits a second-order symmetry-breaking phase transition from a supercritical racemic liquid into subcritical D-rich and L-rich liquids. We determine the infinite-size critical temperature using the fourth-order Binder cumulant, and we show that the finite-size scaling behavior of the order parameter is compatible with the 3D Ising universality class. We also study the spontaneous D-rich to L-rich transition at a slightly subcritical temperature of T = 0.985Tc, and our findings indicate that the free energy barrier for this transformation increases with system size as N2/3, where N is the number of molecules, consistent with a surface-dominated phenomenon. The critical behavior observed herein suggests a mechanism for chirality selection in which a liquid of chiral molecules spontaneously forms a phase enriched in one of the two enantiomers as the temperature is lowered below the critical point. Furthermore, the increasing free energy barrier with system size indicates that fluctuations between the L-rich and D-rich phases are suppressed as the size of the system increases, trapping it in one of the two enantiomerically enriched phases. Such a process could provide the basis for an alternative explanation for the origin of biological homochirality. We also conjecture the possibility of observing nucleation at subcritical temperatures under the action of a suitable chiral external field., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164-168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. We explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Measuring the thermal conductivity (κ) of water at extreme conditions is a challenging task, and few experimental data are available. We predict κ for temperatures and pressures relevant to the conditions of the Earth mantle, between 1,000 and 2,000 K and up to 22 GPa. We employ close to equilibrium molecular dynamics simulations and a deep neural network potential fitted to density functional theory data. We then interpret our results by computing the equation of state of water on a fine grid of points and using a simple model for κ. We find that the thermal conductivity is weakly dependent on temperature and monotonically increases with pressure with an approximate square-root behavior. In addition, we show how the increase of κ at high pressure, relative to ambient conditions, is related to the corresponding increase in the sound velocity. Although the relationships between the thermal conductivity, pressure and sound velocity established here are not rigorous, they are sufficiently accurate to allow for a robust estimate of the thermal conductivity of water in a broad range of temperatures and pressures, where experiments are still difficult to perform.

We propose a quantum Monte Carlo approach to solve the many-body Schrödinger equation for the electronic ground state. The method combines optimization from variational Monte Carlo and propagation from auxiliary field quantum Monte Carlo in a way that significantly alleviates the sign problem. In application to molecular systems, we obtain highly accurate results for configurations dominated by either dynamic or static electronic correlation.

In this work, we construct distinct first-principles-based machine-learning models of CO 2 , reproducing the potential energy surface of the PBE-D3, BLYP-D3, SCAN, and SCAN-rvv10 approximations of density functional theory. We employ the Deep Potential methodology to develop the models and consequently achieve a significant computational efficiency over ab initio molecular dynamics (AIMD) that allows for larger system sizes and time scales to be explored. Although our models are trained only with liquid-phase configurations, they are able to simulate a stable interfacial system and predict vapor-liquid equilibrium properties, in good agreement with results from the literature. Because of the computational efficiency of the models, we are also able to obtain transport properties, such as viscosity and diffusion coefficients. We find that the SCAN-based model presents a temperature shift in the position of the critical point, while the SCAN-rvv10-based model shows improvement but still exhibits a temperature shift that remains approximately constant for all properties investigated in this work. We find that the BLYP-D3-based model generally performs better for the liquid phase and vapor-liquid equilibrium properties, but the PBE-D3-based model is better suited for predicting transport properties.

We model, via large-scale molecular dynamics simulations, the isothermal compression of low-density amorphous ice (LDA) to generate high-density amorphous ice (HDA) and the corresponding decompression extending to negative pressures to recover the low-density amorphous phase (LDA HDA ). Both LDA and HDA are nearly hyperuniform and are characterized by a dynamical HBN, showing that amorphous ices are nonstatic materials and implying that nearly hyperuniformity can be accommodated in dynamical networks. In correspondence with both the LDA-to-HDA and the HDA-to-LDA HDA phase transitions, the (partial) activation of rotational degrees of freedom activates a cascade effect that induces a drastic change in the connectivity and a pervasive reorganization of the HBN topology which, ultimately, break the samples' hyperuniform character. Key to this effect is the rapid rate at which changes occur, and not their magnitude. The inspection of structural properties from the short- to the long-range shows that signatures of metastability are present at all length-scales, hence providing further solid evidence in support of the liquid-liquid critical point scenario. LDA and LDA HDA differ in terms of HBN and structural properties, implying that they are distinct low-density glasses. Our work unveils the role of molecular rotations in the phase transitions between amorphous ices and shows how the unfreezing of rotational degrees of freedom generates a cascade effect that propagates over multiple length-scales. Our findings greatly improve our basic understanding of water and amorphous ices and can potentially impact the field of molecular network-forming materials at large.

Onosma hookeri Clarke. var. longiforum Duthie (OHC-LD), one of the traditional Tibetan medicine, has been found many functions, including removing heat to cool blood, nourishing lung and inhibiting bacteria. In order to study the polysaccharides in OHC-LD water extract, the optimal extraction progress of polysaccharides of the roots of OHC-LD by response surface method designed with three-factor three-level Box-Behnken method and the antioxidant capacity and immune activity of the crude polysaccharide were studied in this investigation. Under the best conditions, the extraction yield of polysaccharide was 3.19±0.09% (n = 3). After purification, the crude polysaccharide was obtained with polysaccharide contents of 42.57%, which demonstrated stronger DPPH scavenging activity than BHT at low concentrations (<625 µg/mL), and comparable ABTS radical scavenging activity as BHT at high concentrations (≥1250 µg/mL). Additionally, it also exhibited a certain cell proliferation activity and an enhancement of the phagocytic ability of RAW264.7 cells. This study revealed that the crude polysaccharide from the roots of OHC-LD might be exploited as a natural antioxidant and immune enhance agent in the future in both medical and food industry.

A long-standing question in water research is the possibility that supercooled liquid water can undergo a liquid-liquid phase transition (LLT) into high- and low-density liquids. We used several complementary molecular simulation techniques to evaluate the possibility of an LLT in an ab initio neural network model of water trained on density functional theory calculations with the SCAN exchange correlation functional. We conclusively show the existence of a first-order LLT and an associated critical point in the SCAN description of water, representing the first definitive computational evidence for an LLT in water from first principles.

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

We studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice-liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius-Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

SignificanceIn X-ray absorption spectroscopy, an electron-hole excitation probes the local atomic environment. The interpretation of the spectra requires challenging theoretical calculations, particularly in a system like liquid water, where quantum many-body effects and molecular disorder play an important role. Recent advances in theory and simulation make possible new calculations that are in good agreement with experiment, without recourse to commonly adopted approximations. Based on these calculations, the three features observed in the experimental spectra are unambiguously attributed to excitonic effects with different characteristic correlation lengths, which are distinctively affected by perturbations of the underlying H-bond structure induced by temperature changes and/or by isotopic substitution. The emerging picture of the water structure is fully consistent with the conventional tetrahedral model.

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

New developments in the field of topological matter are often driven by materials discovery, including novel topological insulators, Dirac semimetals, and Weyl semimetals. In the last few years, large efforts have been made to classify all known inorganic materials with respect to their topology. Unfortunately, a large number of topological materials suffer from non-ideal band structures. For example, topological bands are frequently convoluted with trivial ones, and band structure features of interest can appear far below the Fermi level. This leaves just a handful of materials that are intensively studied. Finding strategies to design new topological materials is a solution. Here, a new mechanism is introduced, which is based on charge density waves and non-symmorphic symmetry, to design an idealized Dirac semimetal. It is then shown experimentally that the antiferromagnetic compound GdSb 0.46 Te 1.48 is a nearly ideal Dirac semimetal based on the proposed mechanism, meaning that most interfering bands at the Fermi level are suppressed. Its highly unusual transport behavior points to a thus far unknown regime, in which Dirac carriers with Fermi energy very close to the node seem to gradually localize in the presence of lattice and magnetic disorder., (© 2021 Wiley-VCH GmbH.)

Image registration is a fundamental task in medical image analysis. Recently, many deep learning based image registration methods have been extensively investigated due to their comparable performance with the state-of-the-art classical approaches despite the ultra-fast computational time. However, the existing deep learning methods still have limitations in the preservation of original topology during the deformation with registration vector fields. To address this issues, here we present a cycle-consistent deformable image registration, dubbed CycleMorph. The cycle consistency enhances image registration performance by providing an implicit regularization to preserve topology during the deformation. The proposed method is so flexible that it can be applied for both 2D and 3D registration problems for various applications, and can be easily extended to multi-scale implementation to deal with the memory issues in large volume registration. Experimental results on various datasets from medical and non-medical applications demonstrate that the proposed method provides effective and accurate registration on diverse image pairs within a few seconds. Qualitative and quantitative evaluations on deformation fields also verify the effectiveness of the cycle consistency of the proposed method., Competing Interests: Declaration of Competing Interest All authors have participated in (a) conception and design, or analysis and interpretation of the data; (b) drafting the article or revising it critically for important intellectual content; and (c) approval of the final version. This manuscript has not been submitted to, nor is under review at, another journal or other publishing venue. The authors have no affiliation with any organization with a direct or indirect financial interest in the subject matter discussed in the manuscript., (Copyright © 2021. Published by Elsevier B.V.)

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. The fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII^{''} and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

Much attention has been devoted to water's metastable phase behavior, including polyamorphism (multiple amorphous solid phases), and the hypothesized liquid-liquid transition and associated critical point. However, the possible relationship between these phenomena remains incompletely understood. Using molecular dynamics simulations of the realistic TIP4P/2005 model, we found a striking signature of the liquid-liquid critical point in the structure of water glasses, manifested as a pronounced increase in long-range density fluctuations at pressures proximate to the critical pressure. By contrast, these signatures were absent in glasses of two model systems that lack a critical point. We also characterized the departure from equilibrium upon vitrification via the non-equilibrium index; water-like systems exhibited a strong pressure dependence in this metric, whereas simple liquids did not. These results reflect a surprising relationship between the metastable equilibrium phenomenon of liquid-liquid criticality and the non-equilibrium structure of glassy water, with implications for our understanding of water phase behavior and glass physics. Our calculations suggest a possible experimental route to probing the existence of the liquid-liquid transition in water and other fluids.

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

We explore the role of long-range interactions in atomistic machine-learning models by analyzing the effects on fitting accuracy, isolated cluster properties, and bulk thermodynamic properties. Such models have become increasingly popular in molecular simulations given their ability to learn highly complex and multi-dimensional interactions within a local environment; however, many of them fundamentally lack a description of explicit long-range interactions. In order to provide a well-defined benchmark system with precisely known pairwise interactions, we chose as the reference model a flexible version of the Extended Simple Point Charge (SPC/E) water model. Our analysis shows that while local representations are sufficient for predictions of the condensed liquid phase, the short-range nature of machine-learning models falls short in representing cluster and vapor phase properties. These findings provide an improved understanding of the role of long-range interactions in machine learning models and the regimes where they are necessary.

We extend to quenched-disordered systems the variational scheme for real-space renormalization group calculations that we recently introduced for homogeneous spin Hamiltonians. When disorder is present our approach gives access to the flow of the renormalized Hamiltonian distribution, from which one can compute the critical exponents if the correlations of the renormalized couplings retain finite range. Key to the variational approach is the bias potential found by minimizing a convex functional in statistical mechanics. This potential reduces dramatically the Monte Carlo relaxation time in large disordered systems. We demonstrate the method with applications to the two-dimensional dilute Ising model, the random transverse field quantum Ising chain, and the random field Ising in two- and three-dimensional lattices.

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553-823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

The possible existence of a metastable liquid-liquid transition (LLT) and a corresponding liquid-liquid critical point (LLCP) in supercooled liquid water remains a topic of much debate. An LLT has been rigorously proved in three empirically parametrized molecular models of water, and evidence consistent with an LLT has been reported for several other such models. In contrast, experimental proof of this phenomenon has been elusive due to rapid ice nucleation under deeply supercooled conditions. In this work, we combined density functional theory (DFT), machine learning, and molecular simulations to shed additional light on the possible existence of an LLT in water. We trained a deep neural network (DNN) model to represent the ab initio potential energy surface of water from DFT calculations using the Strongly Constrained and Appropriately Normed (SCAN) functional. We then used advanced sampling simulations in the multithermal-multibaric ensemble to efficiently explore the thermophysical properties of the DNN model. The simulation results are consistent with the existence of an LLCP, although they do not constitute a rigorous proof thereof. We fit the simulation data to a two-state equation of state to provide an estimate of the LLCP's location. These combined results-obtained from a purely first-principles approach with no empirical parameters-are strongly suggestive of the existence of an LLT, bolstering the hypothesis that water can separate into two distinct liquid forms., Competing Interests: The authors declare no competing interest.

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical ( NVT ) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H 2 O) 256 ) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. We analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

We introduce a scheme based on machine learning and deep neural networks to model the environmental dependence of the electronic polarizability in insulating materials. Application to liquid water shows that training the network with a relatively small number of molecular configurations is sufficient to predict the polarizability of arbitrary liquid configurations in close agreement with ab initio density functional theory calculations. In combination with a neural network representation of the interatomic potential energy surface, the scheme allows us to calculate the Raman spectra along 2-nanosecond classical trajectories at different temperatures for H 2 O and D 2 O. The vast gains in efficiency provided by the machine learning approach enable longer trajectories and larger system sizes relative to ab initio methods, reducing the statistical error and improving the resolution of the low-frequency Raman spectra. Decomposing the spectra into intramolecular and intermolecular contributions elucidates the mechanisms behind the temperature dependence of the low-frequency and stretch modes.

Quantum ESPRESSO is an open-source distribution of computer codes for quantum-mechanical materials modeling, based on density-functional theory, pseudopotentials, and plane waves, and renowned for its performance on a wide range of hardware architectures, from laptops to massively parallel computers, as well as for the breadth of its applications. In this paper, we present a motivation and brief review of the ongoing effort to port Quantum ESPRESSO onto heterogeneous architectures based on hardware accelerators, which will overcome the energy constraints that are currently hindering the way toward exascale computing.

TiO 2 is a widely used photocatalyst in science and technology and its interface with water is important in fields ranging from geochemistry to biomedicine. Yet, it is still unclear whether water adsorbs in molecular or dissociated form on TiO 2 even for the case of well-defined crystalline surfaces. To address this issue, we simulated the TiO 2 -water interface using molecular dynamics with an ab initio -based deep neural network potential. Our simulations show a dynamical equilibrium of molecular and dissociative adsorption of water on TiO 2 . Water dissociates through a solvent-assisted concerted proton transfer to form a pair of short-lived hydroxyl groups on the TiO 2 surface. Molecular adsorption of water is Δ F = 8.0 ± 0.9 kJ mol -1 lower in free energy than the dissociative adsorption, giving rise to a 5.6 ± 0.5% equilibrium water dissociation fraction at room temperature. Due to the relevance of surface hydroxyl groups to the surface chemistry of TiO 2 , our model might be key to understanding phenomena ranging from surface functionalization to photocatalytic mechanisms., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

In this paper, we consider the momentum operator of a quantum particle directed along the displacement of two of its neighbors. A modified open-path path integral molecular dynamics is presented to sample the distribution of this directional momentum distribution, where we derive and use a new estimator for this distribution. Variationally enhanced sampling is used to obtain this distribution for an example molecule, malonaldehyde, in the very low temperature regime where deep tunneling happens. We find no secondary feature in the directional momentum distribution and that its absence is due to quantum entanglement through a further study of the reduced density matrix.

Multiple scattering and absorption limit the depth at which biological tissues can be imaged with light. In thick unlabeled specimens, multiple scattering randomizes the phase of the field and absorption attenuates light that travels long optical paths. These obstacles limit the performance of transmission imaging. To mitigate these challenges, we developed an epi-illumination gradient light interference microscope (epi-GLIM) as a label-free phase imaging modality applicable to bulk or opaque samples. Epi-GLIM enables studying turbid structures that are hundreds of microns thick and otherwise opaque to transmitted light. We demonstrate this approach with a variety of man-made and biological samples that are incompatible with imaging in a transmission geometry: semiconductors wafers, specimens on opaque and birefringent substrates, cells in microplates, and bulk tissues. We demonstrate that the epi-GLIM data can be used to solve the inverse scattering problem and reconstruct the tomography of single cells and model organisms.

We show that the critical manifold of a statistical mechanical system in the vicinity of a critical point is locally accessible through correlation functions at that point. A practical numerical method is presented to determine the tangent space and the curvature to the critical manifold with variational Monte Carlo renormalization group. Because of the use of a variational bias potential of the coarse-grained variables, critical slowing down is greatly alleviated in the Monte Carlo simulation. In addition, this method is free of truncation error. We study the isotropic Ising model on square and cubic lattices, the anisotropic Ising model, and the tricritical Ising model on square lattices to illustrate the method.

Reliable prediction of the polymorphic energy landscape of a molecular crystal would yield profound insight into drug development in terms of the existence and likelihood of late-appearing polymorphs. However, the computational prediction of molecular crystal polymorphs is highly challenging due to the high dimensionality of conformational and crystallographic space accompanied by the need for relative free energies to within 1 kJ/mol per molecule. In this study, we combine the most successful crystal structure sampling strategy with the most successful first-principles energy ranking strategy of the latest blind test of organic crystal structure prediction methods. Specifically, we present a hierarchical energy ranking approach intended for the refinement of relative stabilities in the final stage of a crystal structure prediction procedure. Such a combined approach provides excellent stability rankings for all studied systems and can be applied to molecular crystals of pharmaceutical importance.

We have synthesized boron nitride nanotubes (BNNTs) in an arc in the presence of boron and nitrogen species. We find that BNNTs are often attached to large nanoparticles, suggesting that root-growth is a likely mechanism for their formation. Moreover, the tube-end nanoparticles are composed of boron, without transition metals, indicating that transition metals are not necessary for the arc synthesis of BNNTs. To gain further insight into this process we have studied key mechanisms for root growth of BNNTs on the surface of a liquid boron droplet by ab initio molecular dynamics simulations. We find that nitrogen atoms reside predominantly on the droplet surface where they organize to form boron nitride islands below 2400 K. To minimize contact with the liquid particle underneath, the islands assume non-planar configurations that are likely precursors for the thermal nucleation of cap structures. Once formed, the caps are stable and can easily incorporate nitrogen and boron atoms at their base, resulting in further growth. Our simulations support the root-growth mechanism of BNNTs and provide comprehensive evidence of the active role played by liquid boron.

The photocatalytic activity of TiO 2 for water splitting has been known for decades, yet the adsorption structure and hydrogen bonding of water at the interface with TiO 2 have remained controversial. We investigate the prototypical aqueous interface with anatase TiO 2 (101) using ab initio molecular dynamics (AIMD) with the strongly constrained and appropriately normed (SCAN) density functional, recently shown to provide an excellent description of the properties of bulk liquid water. We find that water forms a stable bilayer of intact molecules with ice-like dynamics and enhanced dipole moment and polarizability on the anatase surface. The orientational order and H-bond environment of interfacial water are reflected in the computed sum frequency generation (SFG) spectrum, which agrees well with recent measurements in the OH stretching frequency range (3000-3600 cm -1 ). Additional AIMD simulations for a model interface with 66% of dissociated water in the contact layer show that surface hydroxyls disrupt the order in the bilayer and lead to a much faster orientational dynamics of interfacial water. Nonetheless, the computed SFG spectrum for the hydroxylated surface also agrees with experiment, suggesting that SFG measurements in a wider frequency range would be necessary to unambiguously identify the character of interfacial water on anatase.

We introduce a general framework for constructing coarse-grained potential models without ad hoc approximations such as limiting the potential to two- and/or three-body contributions. The scheme, called the Deep Coarse-Grained Potential (abbreviated DeePCG), exploits a carefully crafted neural network to construct a many-body coarse-grained potential. The network is trained with full atomistic data in a way that preserves the natural symmetries of the system. The resulting model is very accurate and can be used to sample the configurations of the coarse-grained variables in a much faster way than with the original atomistic model. As an application, we consider liquid water and use the oxygen coordinates as the coarse-grained variables, starting from a full atomistic simulation of this system at the ab initio molecular dynamics level. We find that the two-body, three-body, and higher-order oxygen correlation functions produced by the coarse-grained and full atomistic models agree very well with each other, illustrating the effectiveness of the DeePCG model on a rather challenging task.