Your browser doesn't support javascript.
loading
Mostrar: 20 | 50 | 100
Resultados 1 - 20 de 113
Filtrar
Más filtros

Banco de datos
País/Región como asunto
Tipo del documento
Intervalo de año de publicación
1.
Proc Natl Acad Sci U S A ; 121(14): e2308668121, 2024 Apr 02.
Artículo en Inglés | MEDLINE | ID: mdl-38551836

RESUMEN

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.

2.
Proc Natl Acad Sci U S A ; 120(46): e2302468120, 2023 Nov 14.
Artículo en Inglés | MEDLINE | ID: mdl-37931100

RESUMEN

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 H2O) 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 H2O, 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.

3.
Proc Natl Acad Sci U S A ; 119(33): e2207294119, 2022 Aug 16.
Artículo en Inglés | MEDLINE | ID: mdl-35939708

RESUMEN

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.

4.
Proc Natl Acad Sci U S A ; 119(20): e2201258119, 2022 05 17.
Artículo en Inglés | MEDLINE | ID: mdl-35561212

RESUMEN

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.

5.
Faraday Discuss ; 249(0): 98-113, 2024 Feb 06.
Artículo en Inglés | MEDLINE | ID: mdl-37791889

RESUMEN

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.

6.
Phys Chem Chem Phys ; 26(35): 23080-23088, 2024 Sep 11.
Artículo en Inglés | MEDLINE | ID: mdl-39177036

RESUMEN

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.

7.
J Chem Phys ; 159(5)2023 Aug 07.
Artículo en Inglés | MEDLINE | ID: mdl-37531247

RESUMEN

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.

8.
J Chem Phys ; 159(11)2023 Sep 21.
Artículo en Inglés | MEDLINE | ID: mdl-37712791

RESUMEN

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.

9.
J Chem Phys ; 159(5)2023 Aug 07.
Artículo en Inglés | MEDLINE | ID: mdl-37526163

RESUMEN

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.

10.
Proc Natl Acad Sci U S A ; 117(42): 26040-26046, 2020 10 20.
Artículo en Inglés | MEDLINE | ID: mdl-33008883

RESUMEN

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.

11.
Phys Rev Lett ; 129(25): 255702, 2022 Dec 16.
Artículo en Inglés | MEDLINE | ID: mdl-36608224

RESUMEN

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.

12.
J Chem Phys ; 157(5): 054504, 2022 Aug 07.
Artículo en Inglés | MEDLINE | ID: mdl-35933204

RESUMEN

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.

13.
J Chem Phys ; 156(12): 124107, 2022 Mar 28.
Artículo en Inglés | MEDLINE | ID: mdl-35364869

RESUMEN

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.

14.
Phys Rev Lett ; 126(23): 236001, 2021 Jun 11.
Artículo en Inglés | MEDLINE | ID: mdl-34170175

RESUMEN

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.

15.
J Chem Phys ; 154(3): 034111, 2021 Jan 21.
Artículo en Inglés | MEDLINE | ID: mdl-33499637

RESUMEN

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.

16.
Phys Rev Lett ; 125(19): 190601, 2020 Nov 06.
Artículo en Inglés | MEDLINE | ID: mdl-33216587

RESUMEN

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.

17.
Phys Chem Chem Phys ; 22(19): 10592-10602, 2020 May 21.
Artículo en Inglés | MEDLINE | ID: mdl-32377657

RESUMEN

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 H2O and D2O. 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.

18.
Nature ; 510(7505): 385-8, 2014 Jun 19.
Artículo en Inglés | MEDLINE | ID: mdl-24943954

RESUMEN

Liquid water's isothermal compressibility and isobaric heat capacity, and the magnitude of its thermal expansion coefficient, increase sharply on cooling below the equilibrium freezing point. Many experimental, theoretical and computational studies have sought to understand the molecular origin and implications of this anomalous behaviour. Of the different theoretical scenarios put forward, one posits the existence of a first-order phase transition that involves two forms of liquid water and terminates at a critical point located at deeply supercooled conditions. Some experimental evidence is consistent with this hypothesis, but no definitive proof of a liquid-liquid transition in water has been obtained to date: rapid ice crystallization has so far prevented decisive measurements on deeply supercooled water, although this challenge has been overcome recently. Computer simulations are therefore crucial for exploring water's structure and behaviour in this regime, and have shown that some water models exhibit liquid-liquid transitions and others do not. However, recent work has argued that the liquid-liquid transition has been mistakenly interpreted, and is in fact a liquid-crystal transition in all atomistic models of water. Here we show, by studying the liquid-liquid transition in the ST2 model of water with the use of six advanced sampling methods to compute the free-energy surface, that two metastable liquid phases and a stable crystal phase exist at the same deeply supercooled thermodynamic condition, and that the transition between the two liquids satisfies the thermodynamic criteria of a first-order transition. We follow the rearrangement of water's coordination shell and topological ring structure along a thermodynamically reversible path from the low-density liquid to cubic ice. We also show that the system fluctuates freely between the two liquid phases rather than crystallizing. These findings provide unambiguous evidence for a liquid-liquid transition in the ST2 model of water, and point to the separation of time scales between crystallization and relaxation as being crucial for enabling it.


Asunto(s)
Modelos Moleculares , Agua/química , Temperatura , Termodinámica
19.
J Chem Phys ; 152(2): 024106, 2020 Jan 14.
Artículo en Inglés | MEDLINE | ID: mdl-31941303

RESUMEN

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.

20.
J Chem Phys ; 152(20): 204116, 2020 May 29.
Artículo en Inglés | MEDLINE | ID: mdl-32486691

RESUMEN

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.

SELECCIÓN DE REFERENCIAS
DETALLE DE LA BÚSQUEDA