Showing 1–50 of 79 results for author: Car, R

Deuteration removes quantum dipolar defects from KDP crystals

Authors: Bingjia Yang, Pinchen Xie, Roberto Car

Abstract: The structural, dielectric, and thermodynamic properties of the hydrogen-bonded ferroelectric crystal potassium dihydrogen phosphate ($\mathbf{KH_2PO_4}$), KDP for short, differ significantly from those of DKDP ($\mathbf{KD_2PO_4}$). It is well established that deuteration affects the interplay of hydrogen-bond switches and heavy ion displacements that underlie the emergence of macroscopic polariz… ▽ More The structural, dielectric, and thermodynamic properties of the hydrogen-bonded ferroelectric crystal potassium dihydrogen phosphate ($\mathbf{KH_2PO_4}$), KDP for short, differ significantly from those of DKDP ($\mathbf{KD_2PO_4}$). It is well established that deuteration affects the interplay of hydrogen-bond switches and heavy ion displacements that underlie the emergence of macroscopic polarization, but a detailed microscopic model is missing. Here we show that all-atom path integral molecular dynamics simulations can predict the isotope effects, revealing the microscopic mechanism that differentiates KDP and DKDP. Proton tunneling in the hydrogen bonds generates phosphate configurations that do not contribute to the polarization. These dipolar defects are always present in KDP, but disappear at low temperatures in DKDP, which behaves more classically. Quantum disorder confers residual entropy to ferroelectric KDP, explaining its lower spontaneous polarization and transition entropy relative to DKDP. Tunneling should also contribute to the anomalous heat capacity observed in KDP near absolute zero. The prominent role of quantum fluctuations in KDP is related to the unusual strength of the hydrogen bonds in this system and should be equally important in the other crystals of the KDP family, which exhibit similar isotope effects. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Quantum effects in the H-bond symmetrization and in the thermodynamic properties of high pressure ice

Authors: Marco Cherubini, Lorenzo Monacelli, Bingjia Yang, Roberto Car, Michele Casula, Francesco Mauri

Abstract: We investigate the structural and thermodynamic properties of high-pressure ice by incorporating quantum anharmonicity at a non-perturbative level. Quantum fluctuations reduce the critical pressure of the phase transition between phase VIII (with asymmetric H-bonds) and phase X (with symmetric H-bonds) by 65 GPa from its classical value of 116 GPa at 0K. Moreover, quantum effects make it temperatu… ▽ More We investigate the structural and thermodynamic properties of high-pressure ice by incorporating quantum anharmonicity at a non-perturbative level. Quantum fluctuations reduce the critical pressure of the phase transition between phase VIII (with asymmetric H-bonds) and phase X (with symmetric H-bonds) by 65 GPa from its classical value of 116 GPa at 0K. Moreover, quantum effects make it temperature-independent over a wide temperature range (0K-300K), in agreement with experimental estimates obtained through vibrational spectroscopy and in striking contrast to the strong temperature dependence found in the classical approximation. The equation of state shows fingerprints of the transition in accordance with experimental evidence. Additionally, we demonstrate that, within our approach, proton disorder in phase VII has a negligible impact on the occurrence of phase X. Finally, we reproduce with high accuracy the 10 GPa isotope shift due to the hydrogen-to-deuterium substitution. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Critical behavior in a chiral molecular model

Authors: Pablo M. Piaggi, Roberto Car, Frank H. Stillinger, Pablo G. Debenedetti

Abstract: 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 an… ▽ More 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 $T\approx0.985 T_c$ and our findings indicate that the free energy barrier for this transformation increases with system size as $N^{2/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. △ Less

Submitted 11 June, 2023; originally announced June 2023.

Comments: 9 pages, 7 figures

Journal ref: J. Chem. Phys. 159, 114502 (2023)

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar

Authors: Pablo M. Piaggi, Annabella Selloni, Athanassios Z. Panagiotopoulos, Roberto Car, Pablo G. Debenedetti

Abstract: 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 at 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… ▽ More 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 at 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. We apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in 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 1 nm 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 on ice nucleation. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: 12 pages, 6 figures

Journal ref: Faraday Discuss., 2024, Advance Article

Thermal Conductivity of Water at Extreme Conditions

Authors: Cunzhi Zhang, Marcello Puligheddu, Linfeng Zhang, Roberto Car, Giulia Galli

Abstract: 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… ▽ More 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. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Melting curves of ice polymorphs in the vicinity of the liquid-liquid critical point

Authors: Pablo M. Piaggi, Thomas E. Gartner III, Roberto Car, Pablo G. Debenedetti

Abstract: The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate in part due to the challenges associated with providing definitive experimental evidence. Pioneering work by Mishima and Stanley [Nature 392, 164 (1998) and Phys. Rev. Lett. 85, 334 (2000)] sought to shed light on this problem by studying the melting curves of different ice polymorphs… ▽ More The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate in part due to the challenges associated with providing definitive experimental evidence. Pioneering work by Mishima and Stanley [Nature 392, 164 (1998) and Phys. Rev. Lett. 85, 334 (2000)] 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 location of the liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in slope of the melting curves, Mishima suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. Here, we explore this conjecture using molecular dynamics simulations with a purely-predictive machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII using this model and find that the melting lines of all the studied ice polymorphs are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in slope of the melting lines upon crossing of the locus of maximum compressibility of the liquid. Finally, we analyze critically the literature in light of our findings, and conclude that the scenario in which melting curves are supercritical is favored by the most recent computational and experimental evidence. Thus, although the preponderance of experimental and computational evidence is consistent with the existence of a second critical point in water, the behavior of the melting lines of ice polymorphs does not provide strong evidence in support of this viewpoint, according to our calculations. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: 21 pages, 4 figures, supplementary information

Journal ref: J. Chem. Phys. 159, 054502 (2023)

Ab Initio Generalized Langevin Equation

Authors: Pinchen Xie, Roberto Car, Weinan E

Abstract: We introduce a machine learning-based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables 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 form… ▽ More We introduce a machine learning-based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables 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 multi-scale 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 collective variables 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. △ Less

Submitted 15 February, 2024; v1 submitted 11 November, 2022; originally announced November 2022.

Liquid-liquid transition in water from first principles

Authors: Thomas E. Gartner III, Pablo M. Piaggi, Roberto Car, Athanassios Z. Panagiotopoulos, Pablo G. Debenedetti

Abstract: A longstanding 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 cor… ▽ More A longstanding 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. △ Less

Submitted 2 November, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 14 pages, 3 figures. Supplemental material contains 11 pages, 8 figures

Ab initio multi-scale modeling of ferroelectrics: The case of PbTiO3

Authors: Pinchen Xie, Yixiao Chen, Weinan E, Roberto Car

Abstract: We report an ab initio multi-scale study of lead titanate using the Deep Potential (DP) models, a family of machine learning-based atomistic models, trained on first-principles density functional theory data, to represent potential and polarization surfaces. Our approach includes anharmonic effects beyond the limitations of reduced models and of the linear approximation for the polarization. The c… ▽ More We report an ab initio multi-scale study of lead titanate using the Deep Potential (DP) models, a family of machine learning-based atomistic models, trained on first-principles density functional theory data, to represent potential and polarization surfaces. Our approach includes anharmonic effects beyond the limitations of reduced models and of the linear approximation for the polarization. The calculated enthalpy, spontaneous polarization, specific heat and dielectric susceptibility agree well with experiments on single crystals. In addition, we study how the free energy depends on the polarization with enhanced sampling methods, further supporting the first-order and order-disorder character of the transition. The latter is evidenced by persistence of local dipoles above the transition temperature. The simulated free energy surface as a function of the global polarization leads to a Landau-Devonshire theory of the single domain crystal. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Homogeneous ice nucleation in an ab initio machine learning model of water

Authors: Pablo M. Piaggi, Jack Weis, Athanassios Z. Panagiotopoulos, Pablo G. Debenedetti, Roberto Car

Abstract: 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-funct… ▽ More 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 (DFT) 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. △ Less

Submitted 24 August, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 20 pages, 5 figures

Journal ref: Proc. Natl. Acad. Sci. 2021, 119 (33) e2207294119

Viscosity in water from first-principles and deep-neural-network simulations

Authors: Cesare Malosso, Linfeng Zhang, Roberto Car, Stefano Baroni, Davide Tisi

Abstract: We report on an extensive study of the viscosity of liquid water at near-ambient conditions, performed within the Green-Kubo theory of linear response and equilibrium ab initio molecular dynamics (AIMD), based on density-functional theory (DFT). In order to cope with the long simulation times necessary to achieve an acceptable statistical accuracy, our ab initio approach is enhanced with deep-neur… ▽ More We report on an extensive study of the viscosity of liquid water at near-ambient conditions, performed within the Green-Kubo theory of linear response and equilibrium ab initio molecular dynamics (AIMD), based on density-functional theory (DFT). In order to cope with the long simulation times necessary to achieve an acceptable statistical accuracy, our ab initio approach is enhanced with deep-neural-network potentials (NNP). This approach is first validated against AIMD results, obtained by using the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional and paying careful attention to crucial, yet often overlooked, aspects of the statistical data analysis. Then, we train a second NNP to a dataset generated from the Strongly Constrained and Appropriately Normed (SCAN) functional. Once the error resulting from the imperfect prediction of the melting line is offset by referring the simulated temperature to the theoretical melting one, our SCAN predictions of the shear viscosity of water are in very good agreement with experiments. △ Less

Submitted 2 July, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Journal ref: npj Comput Mater 8, 139 (2022)

Many-Body Effects in the X-ray Absorption Spectra of Liquid Water

Authors: Fujie Tang, Zhenglu Li, Chunyi Zhang, Steven G. Louie, Roberto Car, Diana Y. Qiu, Xifan Wu

Abstract: X-ray absorption spectroscopy (XAS) is a powerful experimental technique to probe the local order in materials with core electron excitations. Experimental interpretation requires supporting theoretical calculations. For water, these calculations are very demanding and, to date, could only be done with major approximations that limited the accuracy of the calculated spectra. This prompted an inten… ▽ More X-ray absorption spectroscopy (XAS) is a powerful experimental technique to probe the local order in materials with core electron excitations. Experimental interpretation requires supporting theoretical calculations. For water, these calculations are very demanding and, to date, could only be done with major approximations that limited the accuracy of the calculated spectra. This prompted an intense debate on whether a substantial revision of the standard picture of tetrahedrally bonded water was necessary to improve the agreement of theory and experiment. Here, we report a new first-principles calculation of the XAS of water that avoids the approximations of prior work thanks to recent advances in electron excitation theory. The calculated XAS spectra, and their variation with changes of temperature and/or with isotope substitution, are in excellent quantitative agreement with experiments. The approach requires accurate quasi-particle wavefunctions beyond density functional theory approximations, accounts for the dynamics of quasi-particles and includes dynamic screening as well as renormalization effects due to the continuum of valence-level excitations. The three features observed in the experimental spectra are unambiguously attributed to excitonic effects. The pre-edge feature is associated to a bound intramolecular exciton, the main-edge feature is associated to an exciton localized within the coordination shell of the excited molecule, while the post-edge one is delocalized over more distant neighbors, as expected for a resonant state. The three features probe the local order at short, intermediate, and longer range relative to the excited molecule. The calculated spectra are fully consistent with a standard tetrahedral picture of water. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 11 pages, 3 figures

Journal ref: Proc. Natl. Acad. Sci. U.S.A., 2022, 119, e2201258119

Heat transport in liquid water from first-principles and deep-neural-network simulations

Authors: Davide Tisi, Linfeng Zhang, Riccardo Bertossa, Han Wang, Roberto Car, Stefano Baroni

Abstract: We compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represe… ▽ More We compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way, that a DNN model, trained on meta-GGA (SCAN) data, reduce the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two. △ Less

Submitted 23 December, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Journal ref: Phys. Rev. B 104 (2021), 224202

Manifestations of metastable criticality in the long-range structure of model water glasses

Authors: Thomas E. Gartner III, Salvatore Torquato, Roberto Car, Pablo G. Debenedetti

Abstract: 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… ▽ More 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. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: 9 pages, 5 figures. Published in Nature Communications (2021)

Journal ref: Nat Commun 12, 3398 (2021)

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Authors: Pablo M. Piaggi, Roberto Car

Abstract: Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, even below this temperature the formation of ice XI is kinetically hindered and experimentally it is obtained by doping ice with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton… ▽ More Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, even below this temperature the formation of ice XI is kinetically hindered and experimentally it is obtained by doping ice with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping in molecular dynamics simulations using a bias potential that enhances the formation of ionic defects. The recombination of the ions thus formed proceeds through fast migration of the hydroxide and results in the jump of protons along a hydrogen bond loop. This provides a physical and expedite way to change the proton configuration, and to accelerate diffusion in proton configuration space. A key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modeling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter to distinguish the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of our calculation. △ Less

Submitted 14 May, 2024; v1 submitted 22 January, 2021; originally announced January 2021.

Comments: 17 pages, 10 figures

Journal ref: Mol. Phys. 2021, e1916634

Phase equilibrium of water with hexagonal and cubic ice using the SCAN functional

Authors: Pablo M. Piaggi, Athanassios Z. Panagiotopoulos, Pablo G. Debenedetti, Roberto Car

Abstract: 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 towards studying ice nucleation. The machine learning model is based on deep neura… ▽ More 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 towards 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 ice 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. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: 20 pages, 9 figures

Journal ref: J. Chem. Theory Comput. 2021, 17, 5, 3065-3077

Band Engineering of Dirac Semimetals using Charge Density Waves

Authors: Shiming Lei, Samuel M. L. Teicher, Andreas Topp, Kehan Cai, Jingjing Lin, Fanny Rodolakis, Jessica L. McChesney, Maxim Krivenkov, Dmitry Marchenko, Andrei Varykhalov, Christian R. Ast, Roberto Car, Jennifer Cano, Maia G. Vergniory, N. Phuan Ong, Leslie M. Schoop

Abstract: 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 performed 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… ▽ More 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 performed 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 we introduce a new mechanism that is based on charge density waves and non-symmorphic symmetry to design an idealized Dirac semimetal. We then show 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. △ Less

Submitted 21 June, 2021; v1 submitted 1 September, 2020; originally announced September 2020.

Continuous-time Monte Carlo Renormalization Group

Authors: Yantao Wu, Roberto Car

Abstract: We implement Monte Carlo Renormalization Group (MCRG) in the continuous-time Monte Carlo simulation of a quantum system. We demonstrate numerically the emergent isotropy between space and time at large distances for the systems that exhibit Lorentz invariance at quantum criticality. This allows us to estimate accurately the sound velocity for these quantum systems. $Q$-state Potts models in one an… ▽ More We implement Monte Carlo Renormalization Group (MCRG) in the continuous-time Monte Carlo simulation of a quantum system. We demonstrate numerically the emergent isotropy between space and time at large distances for the systems that exhibit Lorentz invariance at quantum criticality. This allows us to estimate accurately the sound velocity for these quantum systems. $Q$-state Potts models in one and two space dimensions are used to illustrate the method. △ Less

Submitted 5 October, 2020; v1 submitted 19 April, 2020; originally announced April 2020.

Journal ref: Phys. Rev. B 102, 014456 (2020)

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

Authors: Pablo M. Piaggi, Roberto Car

Abstract: 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 a… ▽ More 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. △ Less

Submitted 11 May, 2020; v1 submitted 17 April, 2020; originally announced April 2020.

Comments: 9 pages, 6 figures

Raman Spectrum and Polarizability of Liquid Water from Deep Neural Networks

Authors: Grace M. Sommers, Marcos F. Calegari Andrade, Linfeng Zhang, Han Wang, Roberto Car

Abstract: 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… ▽ More 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. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance

Authors: Hsin-Yu Ko, Junteng Jia, Biswajit Santra, Xifan Wu, Roberto Car, Robert A. DiStasio Jr

Abstract: By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semi-local density functional theory (DFT), and thereby furnish a more accurate and reliable description of the electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX… ▽ More By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semi-local density functional theory (DFT), and thereby furnish a more accurate and reliable description of the 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 have devised a linear-scaling yet formally exact 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 at the hybrid DFT level. We provide a comprehensive description of the exx algorithm, which is currently implemented in the Quantum ESPRESSO program and employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize high-performance computing (HPC) resources. This is followed by a critical assessment of the accuracy and parallel performance of this approach when performing AIMD simulations of liquid water in the canonical ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT based AIMD simulations of condensed-phase systems containing 500-1000 atoms with a walltime cost that is comparable to semi-local DFT. In doing so, exx takes us closer to routinely performing AIMD simulations of large-scale condensed-phase systems for sufficiently long timescales at the hybrid DFT level of theory. △ Less

Submitted 6 May, 2021; v1 submitted 24 November, 2019; originally announced November 2019.

Comments: 34 pages and 10 figures

Influence of point defects on the electronic and topological properties of monolayer WTe$_2$

Authors: Lukas Muechler, Wei Hu, Chao Yang, Lin Lin, Roberto Car

Abstract: In some topological insulators, such as graphene and WTe$_2$, band inversion originates from chemical bonding and space group symmetry, in contrast to materials such as Bi$_2$Se$_3$, where the band inversion derives from relativistic effects in the atoms. In the former, band inversion is susceptible to changes of the chemical environment, e.g. by defects, while the latter are less affected by defe… ▽ More In some topological insulators, such as graphene and WTe$_2$, band inversion originates from chemical bonding and space group symmetry, in contrast to materials such as Bi$_2$Se$_3$, where the band inversion derives from relativistic effects in the atoms. In the former, band inversion is susceptible to changes of the chemical environment, e.g. by defects, while the latter are less affected by defects due to the larger energy scale associated with atomic relativistic effects. Motivated by recent experiments, we study the effect of Te-vacancies and Te-adatoms on the electronic properties of WTe$_2$. We find that the Te-vacancies have a formation energy of $2.21$ eV, while the formation energy of the Te-adatoms is much lower with $0.72$ eV. The vacancies strongly influence the band structure and we present evidence that band inversion is already reversed at the nominal composition of WTe$_{1.97}$. In contrast, we show that the adatoms do not change the electronic structure in the vicinity of the Fermi level and thus the topological properties. Our findings indicate that Te-adatoms should be present in thin films that are grown in a Te-rich environment, and we suggest that they have been observed in scanning tunneling microscopy experiments. △ Less

Submitted 8 November, 2019; v1 submitted 5 November, 2019; originally announced November 2019.

Journal ref: Phys. Rev. B 102, 041103 (2020)

Deep neural network for the dielectric response of insulators

Authors: Linfeng Zhang, Mohan Chen, Xifan Wu, Han Wang, Weinan E, Roberto Car

Abstract: We introduce a deep neural network to model in a symmetry preserving way the environmental dependence of the centers of the electronic charge. The model learns from ab-initio density functional theory, wherein the electronic centers are uniquely assigned by the maximally localized Wannier functions. When combined with the Deep Potential model of the atomic potential energy surface, the scheme pred… ▽ More We introduce a deep neural network to model in a symmetry preserving way the environmental dependence of the centers of the electronic charge. The model learns from ab-initio density functional theory, wherein the electronic centers are uniquely assigned by the maximally localized Wannier functions. When combined with the Deep Potential model of the atomic potential energy surface, the scheme predicts the dielectric response of insulators for trajectories inaccessible to direct ab-initio simulation. The scheme is non-perturbative and can capture the response of a mutating chemical environment. We demonstrate the approach by calculating the infrared spectra of liquid water at standard conditions, and of ice under extreme pressure, when it transforms from a molecular to an ionic crystal. △ Less

Submitted 9 June, 2020; v1 submitted 27 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. B 102, 041121 (2020)

Determination of the Critical Manifold Tangent Space and Curvature with Monte Carlo Renormalization Group

Authors: Yantao Wu, Roberto Car

Abstract: 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 c… ▽ More 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. △ Less

Submitted 5 October, 2020; v1 submitted 19 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. E 100, 022138 (2019)

Active Learning of Uniformly Accurate Inter-atomic Potentials for Materials Simulation

Authors: Linfeng Zhang, De-Ye Lin, Han Wang, Roberto Car, Weinan E

Abstract: An active learning procedure called Deep Potential Generator (DP-GEN) is proposed for the construction of accurate and transferable machine learning-based models of the potential energy surface (PES) for the molecular modeling of materials. This procedure consists of three main components: exploration, generation of accurate reference data, and training. Application to the sample systems of Al, Mg… ▽ More An active learning procedure called Deep Potential Generator (DP-GEN) is proposed for the construction of accurate and transferable machine learning-based models of the potential energy surface (PES) for the molecular modeling of materials. This procedure consists of three main components: exploration, generation of accurate reference data, and training. Application to the sample systems of Al, Mg and Al-Mg alloys demonstrates that DP-GEN can produce uniformly accurate PES models with a minimal number of reference data. △ Less

Submitted 6 February, 2019; v1 submitted 28 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. Materials 3, 023804 (2019)

Monte Carlo Renormalization Group for Systems with Quenched Disorder

Authors: Yantao Wu, Roberto Car

Abstract: 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… ▽ More 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. △ Less

Submitted 5 October, 2020; v1 submitted 22 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. Lett. 125, 190601 (2020)

End-to-end Symmetry Preserving Inter-atomic Potential Energy Model for Finite and Extended Systems

Authors: Linfeng Zhang, Jiequn Han, Han Wang, Wissam A. Saidi, Roberto Car, Weinan E

Abstract: Machine learning models are changing the paradigm of molecular modeling, which is a fundamental tool for material science, chemistry, and computational biology. Of particular interest is the inter-atomic potential energy surface (PES). Here we develop Deep Potential - Smooth Edition (DeepPot-SE), an end-to-end machine learning-based PES model, which is able to efficiently represent the PES for a w… ▽ More Machine learning models are changing the paradigm of molecular modeling, which is a fundamental tool for material science, chemistry, and computational biology. Of particular interest is the inter-atomic potential energy surface (PES). Here we develop Deep Potential - Smooth Edition (DeepPot-SE), an end-to-end machine learning-based PES model, which is able to efficiently represent the PES for a wide variety of systems with the accuracy of ab initio quantum mechanics models. By construction, DeepPot-SE is extensive and continuously differentiable, scales linearly with system size, and preserves all the natural symmetries of the system. Further, we show that DeepPot-SE describes finite and extended systems including organic molecules, metals, semiconductors, and insulators with high fidelity. △ Less

Submitted 20 December, 2018; v1 submitted 23 May, 2018; originally announced May 2018.

Journal ref: Conference on Neural Information Processing Systems (NeurIPS), 2018

Root-Growth of Boron Nitride Nanotubes: Experiments and \textit{Ab Initio} Simulations

Authors: Biswajit Santra, Hsin-Yu Ko, Yao-Wen Yeh, Fausto Martelli, Igor Kaganovich, Yevgeny Raitses, Roberto Car

Abstract: We have synthesized boron nitride nanotubes (BNNTs) in an arc in presence of boron and nitrogen species only, without transition metals. We find that BNNTs are often attached to pure boron nanoparticles, suggesting that root-growth is a likely mechanism for their formation. To gain further insight into this process we have studied key mechanisms for root growth of BNNTs on the surface of a liquid… ▽ More We have synthesized boron nitride nanotubes (BNNTs) in an arc in presence of boron and nitrogen species only, without transition metals. We find that BNNTs are often attached to pure boron nanoparticles, suggesting that root-growth is a likely mechanism for their formation. 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. △ Less

Submitted 6 September, 2018; v1 submitted 30 March, 2018; originally announced March 2018.

Comments: supporting inofrmation inlcuded

Thermal Expansion in Dispersion-Bound Molecular Crystals

Authors: Hsin-Yu Ko, Robert A. DiStasio Jr., Biswajit Santra, Roberto Car

Abstract: We explore how anharmonicity, nuclear quantum effects (NQE), many-body dispersion interactions, and Pauli repulsion influence thermal properties of dispersion-bound molecular crystals. Accounting for anharmonicity with $ab$ $initio$ molecular dynamics yields cell parameters accurate to within 2% of experiment for a set of pyridine-like molecular crystals at finite temperatures and pressures. From… ▽ More We explore how anharmonicity, nuclear quantum effects (NQE), many-body dispersion interactions, and Pauli repulsion influence thermal properties of dispersion-bound molecular crystals. Accounting for anharmonicity with $ab$ $initio$ molecular dynamics yields cell parameters accurate to within 2% of experiment for a set of pyridine-like molecular crystals at finite temperatures and pressures. From the experimental thermal expansion curve, we find that pyridine-I has a Debye temperature just above its melting point, indicating sizable NQE across the entire crystalline range of stability. We find that NQE lead to a substantial volume increase in pyridine-I ($\approx 40$% more than classical thermal expansion at $153$ K) and attribute this to intermolecular Pauli repulsion promoted by intramolecular quantum fluctuations. When predicting delicate properties such as the thermal expansivity, we show that many-body dispersion interactions and sophisticated treatments of Pauli repulsion are needed in dispersion-bound molecular crystals. △ Less

Submitted 1 March, 2018; originally announced March 2018.

Journal ref: Phys. Rev. Materials 2, 055603 (2018)

Why does hydronium diffuse faster than hydroxide in liquid water?

Authors: Mohan Chen, Lixin Zheng, Biswajit Santra, Hsin-Yu Ko, Robert A. DiStasio Jr., Michael L. Klein, Roberto Car, Xifan Wu

Abstract: Proton transfer via hydronium and hydroxide ions in water is ubiquitous. It underlies acid-base chemistry, certain enzyme reactions, and even infection by the flu. Despite two-centuries of investigation, the mechanism underlying why hydronium diffuses faster than hydroxide in water is still not well understood. Herein, we employ state of the art Density Functional Theory based molecular dynamics,… ▽ More Proton transfer via hydronium and hydroxide ions in water is ubiquitous. It underlies acid-base chemistry, certain enzyme reactions, and even infection by the flu. Despite two-centuries of investigation, the mechanism underlying why hydronium diffuses faster than hydroxide in water is still not well understood. Herein, we employ state of the art Density Functional Theory based molecular dynamics, with corrections for nonlocal van der Waals interactions, and self-interaction in the electronic ground state, to model water and the hydrated water ions. At this level of theory, structural diffusion of hydronium preserves the previously recognized concerted behavior. However, by contrast, proton transfer via hydroxide is dominated by stepwise events, arising from a stabilized hyper-coordination solvation structure that discourages proton transfer. Specifically, the latter exhibits non-planar geometry, which agrees with neutron scattering results. Asymmetry in the temporal correlation of proton transfer enables hydronium to diffuse faster than hydroxide. △ Less

Submitted 23 January, 2018; originally announced January 2018.

Comments: accepted by Nature Chemistry, more information will be updated once published

Journal ref: Nature Chemistry 10, 413-419 (2018)

Prediction of a magnetic Weyl semimetal without spin-orbit coupling and strong anomalous Hall effect in the Heusler compensated ferrimagnet Ti2MnAl

Authors: Wujun Shi, Lukas Muechler, Kaustuv Manna, Yang Zhang, Klaus Koepernik, Roberto Car, Jeroen van den Brink, Claudia Felser, Yan Sun

Abstract: We predict a magnetic Weyl semimetal in the inverse Heusler Ti2MnAl, a compensated ferrimagnet with a vanishing net magnetic moment and a Curie temperature of over 650 K. Despite the vanishing net magnetic moment, we calculate a large intrinsic anomalous Hall effect (AHE) of about 300 S/cm. It derives from the Berry curvature distribution of the Weyl points, which are only 14 meV away from the Fer… ▽ More We predict a magnetic Weyl semimetal in the inverse Heusler Ti2MnAl, a compensated ferrimagnet with a vanishing net magnetic moment and a Curie temperature of over 650 K. Despite the vanishing net magnetic moment, we calculate a large intrinsic anomalous Hall effect (AHE) of about 300 S/cm. It derives from the Berry curvature distribution of the Weyl points, which are only 14 meV away from the Fermi level and isolated from trivial bands. Different from antiferromagnets Mn3X (X= Ge, Sn, Ga, Ir, Rh, and Pt), where the AHE originates from the non-collinear magnetic structure, the AHE in Ti2MnAl stems directly from the Weyl points and is topologically protected. The large anomalous Hall conductivity (AHC) together with a low charge carrier concentration should give rise to a large anomalous Hall angle. In contrast to the Co-based ferromagnetic Heusler compounds, the Weyl nodes in Ti2MnAl do not derive from nodal lines due to the lack of mirror symmetries in the inverse Heusler structure. Since the magnetic structure breaks spin-rotation symmetry, the Weyl nodes are stable without SOC. Moreover, because of the large separation between Weyl points of opposite topological charge, the Fermi arcs extent up to 75% of the reciprocal lattice vectors in length. This makes Ti2MnAl an excellent candidate for the comprehensive study of magnetic Weyl semimetals. It is the first example of a material with Weyl points, large anomalous Hall effect and angle despite a vanishing net magnetic moment. △ Less

Submitted 24 September, 2018; v1 submitted 10 January, 2018; originally announced January 2018.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 97, 060406 (2018)

From colossal to zero: Controlling the Anomalous Hall Effect in Magnetic Heusler Compounds via Berry Curvature Design

Authors: Kaustuv Manna, Lukas Muechler, Ting-Hui Kao, Rolf Stinshoff, Yang Zhang, Johannes Gooth, Nitesh Kumar, Guido Kreiner, Klaus Koepernik, Roberto Car, Jürgen Kübler, Gerhard H. Fecher, Chandra Shekhar, Yan Sun, Claudia Felser

Abstract: Since the discovery of the anomalous Hall effect (AHE), the anomalous Hall conductivity (AHC) has been thought to be zero when there is no net magnetization. However, the recently found relation between the intrinsic AHE and the Berry curvature predicts other possibilities, such as a large AHC in non-colinear antiferromagnets with no net magnetization but net Berry curvature. Vice versa, the AHE i… ▽ More Since the discovery of the anomalous Hall effect (AHE), the anomalous Hall conductivity (AHC) has been thought to be zero when there is no net magnetization. However, the recently found relation between the intrinsic AHE and the Berry curvature predicts other possibilities, such as a large AHC in non-colinear antiferromagnets with no net magnetization but net Berry curvature. Vice versa, the AHE in principle could be tuned to zero, irrespective of a finite magnetization. Here, we experimentally investigate this possibility and demonstrate that, the symmetry elements of Heusler magnets can be changed such that the Berry curvature and all the associated properties are switched while leaving the magnetization unaffected. This enables us to tune the AHC from 0 Ω-1cm-1 up to 1600 Ω-1cm-1 with an exceptionally high anomalous Hall angle up to 12 %, while keeping the magnetization same. Our study shows that the AHC can be controlled by selectively changing the Berry curvature distribution, independent of the magnetization. △ Less

Submitted 13 December, 2018; v1 submitted 29 December, 2017; originally announced December 2017.

Comments: Published in Physical Review X. 16 pages, 5 figures

Journal ref: Phys. Rev. X 8, 041045 (2018)

Notes on the Hybrid Monte Carlo Method

Authors: Jeremy C. Palmer, Amir Haji-Akbari, Rakesh S. Singh, Fausto Martelli, Roberto Car, Athanassios Z. Panagiotopoulos, Pablo G. Debenedetti

Abstract: We discuss the detailed balance condition for hybrid Monte Carlo method We discuss the detailed balance condition for hybrid Monte Carlo method △ Less

Submitted 21 December, 2017; originally announced December 2017.

Ab initio theory and modeling of water

Authors: Mohan Chen, Hsin-Yu Ko, Richard C. Remsing, Marcos F. Calegari Andrade, Biswajit Santra, Zhaoru Sun, Annabella Selloni, Roberto Car, Michael L. Klein, John P. Perdew, Xifan Wu

Abstract: Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der… ▽ More Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice I{\it h} at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner. △ Less

Submitted 29 September, 2017; originally announced September 2017.

Journal ref: Proc. Natl. Acad. Sci., 114, 10846 (2017)

Advanced capabilities for materials modelling with Quantum ESPRESSO

Authors: P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. Buongiorno Nardelli, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, M. Cococcioni, N. Colonna, I. Carnimeo, A. Dal Corso, S. de Gironcoli, P. Delugas, R. A. DiStasio Jr., A. Ferretti, A. Floris, G. Fratesi, G. Fugallo, R. Gebauer, U. Gerstmann, F. Giustino, T. Gorni, J. Jia , et al. (25 additional authors not shown)

Abstract: Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudo-potential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wi… ▽ More Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudo-potential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement theirs ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: Psi-k highlight September 2017: psi-k.net/dowlnload/highlights/Highlight_137.pdf; J. Phys.: Condens. Matter, accepted

Journal ref: J. Phys.: Condens. Matter 29, 465901 (2017)

A Variational Approach to Monte Carlo Renormalization Group

Authors: Yantao Wu, Roberto Car

Abstract: We present a Monte Carlo method for computing the renormalized coupling constants and the critical exponents within renormalization theory. The scheme, which derives from a variational principle, overcomes critical slowing down, by means of a bias potential that renders the coarse grained variables uncorrelated. The 2D Ising model is used to illustrate the method. We present a Monte Carlo method for computing the renormalized coupling constants and the critical exponents within renormalization theory. The scheme, which derives from a variational principle, overcomes critical slowing down, by means of a bias potential that renders the coarse grained variables uncorrelated. The 2D Ising model is used to illustrate the method. △ Less

Submitted 26 July, 2017; originally announced July 2017.

Comments: 4 pages, 3 figures, 1 table

Journal ref: Phys. Rev. Lett. 119, 220602 (2017)

Large-scale structure and hyperuniformity of amorphous ices

Authors: Fausto Martelli, Salvatore Torquato, Nicolas Giovambattista, Roberto Car

Abstract: We investigate the large-scale structure of amorphous ices and transitions between their different forms by quantifying their large-scale density fluctuations. Specifically, we simulate the isothermal compression of low-density amorphous ice (LDA) and hexagonal ice (Ih) to produce high-density amorphous ice (HDA). Remarkably, both HDA and LDA are nearly hyperuniform, meaning that they are characte… ▽ More We investigate the large-scale structure of amorphous ices and transitions between their different forms by quantifying their large-scale density fluctuations. Specifically, we simulate the isothermal compression of low-density amorphous ice (LDA) and hexagonal ice (Ih) to produce high-density amorphous ice (HDA). Remarkably, both HDA and LDA are nearly hyperuniform, meaning that they are characterized by an anomalous suppression of large-scale density fluctuations. By contrast, in correspondence with both non-equilibrium phase transitions to HDA, the presence of structural heterogeneities strongly suppresses the hyperuniformity and, remarkably, the system becomes hyposurficial (devoid of "surface-area" fluctuations). Our investigation challenges the largely accepted "frozen-liquid" picture, which views glasses as structurally arrested liquids. Beyond implications for water, our findings enrich our understanding of the structural transformations that occur in glasses with varying pressures. △ Less

Submitted 28 May, 2017; originally announced May 2017.

Journal ref: Phys. Rev. Lett. 119, 136002 (2017)

A well-scaling natural orbital theory

Authors: Ralph Gebauer, Morrel H. Cohen, Roberto Car

Abstract: We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latt… ▽ More We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the one-particle density matrix. △ Less

Submitted 26 November, 2016; originally announced November 2016.

Comments: http://www.pnas.org/cgi/doi/10.1073/pnas.1615729113. arXiv admin note: text overlap with arXiv:1309.3929

Journal ref: Proceedings of the National Academy of Sciences of the United States of America, vol. 113, no. 46, pp. 12913-12918 (2016)

A local order metric for condensed phase environments

Authors: Fausto Martelli, Hsin-Yu Ko, Erdal C. Oguz, Roberto Car

Abstract: We introduce a local order metric (LOM) that measures the degree of order in the neighborhood of an atomic or molecular site in a condensed medium. The LOM maximizes the overlap between the spatial distribution of sites belonging to that neighborhood and the corresponding distribution in a suitable reference system. The LOM takes a value tending to zero for completely disordered environments and t… ▽ More We introduce a local order metric (LOM) that measures the degree of order in the neighborhood of an atomic or molecular site in a condensed medium. The LOM maximizes the overlap between the spatial distribution of sites belonging to that neighborhood and the corresponding distribution in a suitable reference system. The LOM takes a value tending to zero for completely disordered environments and tending to one for environments that match perfectly the reference. The site averaged LOM and its standard deviation define two scalar order parameters, $S$ and $δS$, that characterize with excellent resolution crystals, liquids, and amorphous materials. We show with molecular dynamics simulations that $S$, $δS$ and the LOM provide very insightful information in the study of structural transformations, such as those occurring when ice spontaneously nucleates from supercooled water or when a supercooled water sample becomes amorphous upon progressive cooling. △ Less

Submitted 11 September, 2016; originally announced September 2016.

Journal ref: Phys. Rev. B 97, 064105 (2018)

Assessment of a nonempirical semilocal density functional on solids and surfaces

Authors: Yuxiang Mo, Roberto Car, Viktor N. Staroverov, Gustavo E. Scuseria, Jianmin Tao

Abstract: Recently, Tao and Mo developed a new nonempirical semilocal exchange-correlation density functional. The exchange part of this functional is derived from a density matrix expansion corrected to reproduce the fourth-order gradient expansion in the slowly varying limit, while the correlation part is based on the TPSS correlation model with a modification for the low-density limit. In the present wor… ▽ More Recently, Tao and Mo developed a new nonempirical semilocal exchange-correlation density functional. The exchange part of this functional is derived from a density matrix expansion corrected to reproduce the fourth-order gradient expansion in the slowly varying limit, while the correlation part is based on the TPSS correlation model with a modification for the low-density limit. In the present work, the Tao-Mo functional is assessed by calculations on a variety of solids and jellium surfaces. This includes 22 lattice constants and bulk moduli, 7 cohesive energies, and jellium surface exchange and correlation energies for the density parameter rs in the range from 2 to 3 bohrs. Our calculations show that this meta-generalized gradient approximation can yield consistently remarkable accuracy for the properties considered here, with mean absolute errors of 0.017 Å for lattice constants, 7.0 GPa for bulk moduli, 0.08 eV for cohesive energies, and 35 erg/cm2 for surface exchange-correlation energies, substantially improving upon existing nonempirical semilocal density functionals. △ Less

Submitted 11 August, 2016; v1 submitted 18 July, 2016; originally announced July 2016.

Journal ref: Phys. Rev. B 95, 035118 (2017)

Tilted Dirac Fermions

Authors: Lukas Muechler, Aris Alexandradinata, Titus Neupert, Roberto Car

Abstract: We introduce the notion of a band-inverted, topological semimetal in two-dimensional nonsymmorphic crystals. This notion is materialized in the monolayers of MTe$_2$ (M $=$ W, Mo) if spin-orbit coupling is neglected. We characterize the Dirac band touching topologically by the Wilson loop of the non-Abelian Berry gauge field. An additional feature of the Dirac cone in monolayer MTe$_2$ is that it… ▽ More We introduce the notion of a band-inverted, topological semimetal in two-dimensional nonsymmorphic crystals. This notion is materialized in the monolayers of MTe$_2$ (M $=$ W, Mo) if spin-orbit coupling is neglected. We characterize the Dirac band touching topologically by the Wilson loop of the non-Abelian Berry gauge field. An additional feature of the Dirac cone in monolayer MTe$_2$ is that it tilts over in a Lifshitz transition to produce electron and hole pockets, a type-II Dirac cone. These pockets, together with the pseudospin structure of the Dirac electrons, suggest a unified, topological explanation for the recently-reported, non-saturating magnetoresistance in WTe$_2$, as well as its circular dichroism in photoemission. We complement our analysis and first-principle bandstructure calculations with an $\textit{ab-initio}$-derived-derived tight-binding model for the WTe$_2$ monolayer. △ Less

Submitted 16 May, 2016; v1 submitted 5 April, 2016; originally announced April 2016.

Comments: Changed Title, improved introduction and added references

Journal ref: Phys. Rev. X 6, 041069 (2016)

Fermionic Symmetry-Protected Topological Phase in a Two-Dimensional Hubbard Model

Authors: Cheng-Chien Chen, Lukas Muechler, Roberto Car, Titus Neupert, Joseph Maciejko

Abstract: We study the two-dimensional (2D) Hubbard model using exact diagonalization for spin-1/2 fermions on the triangular and honeycomb lattices decorated with a single hexagon per site. In certain parameter ranges, the Hubbard model maps to a quantum compass model on those lattices. On the triangular lattice, the compass model exhibits collinear stripe antiferromagnetism, implying $d$-density wave char… ▽ More We study the two-dimensional (2D) Hubbard model using exact diagonalization for spin-1/2 fermions on the triangular and honeycomb lattices decorated with a single hexagon per site. In certain parameter ranges, the Hubbard model maps to a quantum compass model on those lattices. On the triangular lattice, the compass model exhibits collinear stripe antiferromagnetism, implying $d$-density wave charge order in the original Hubbard model. On the honeycomb lattice, the compass model has a unique, quantum disordered ground state that transforms nontrivially under lattice reflection. The ground state of the Hubbard model on the decorated honeycomb lattice is thus a 2D fermionic symmetry-protected topological phase. This state -- protected by time-reversal and reflection symmetries -- cannot be connected adiabatically to a free-fermion topological phase. △ Less

Submitted 25 August, 2016; v1 submitted 10 March, 2016; originally announced March 2016.

Comments: Original: 5+3 pages, 4+2 figures, including the supplemental material. Revised: updated title, content, and references in version 2

Journal ref: Phys. Rev. Lett. 117, 096405 (2016)

Inverse design of disordered stealthy hyperuniform spin chains

Authors: Eli Chertkov, Robert A. DiStasio Jr., Ge Zhang, Roberto Car, Salvatore Torquato

Abstract: Positioned between crystalline solids and liquids, disordered many-particle systems which are stealthy and hyperuniform represent new states of matter that are endowed with novel physical and thermodynamic properties. Such stealthy and hyperuniform states are unique in that they are transparent to radiation for a range of wavenumbers around the origin. In this work, we employ recently developed in… ▽ More Positioned between crystalline solids and liquids, disordered many-particle systems which are stealthy and hyperuniform represent new states of matter that are endowed with novel physical and thermodynamic properties. Such stealthy and hyperuniform states are unique in that they are transparent to radiation for a range of wavenumbers around the origin. In this work, we employ recently developed inverse statistical-mechanical methods, which seek to obtain the optimal set of interactions that will spontaneously produce a targeted structure or configuration as a unique ground state, to investigate the spin-spin interaction potentials required to stabilize disordered stealthy hyperuniform one-dimensional (1D) Ising-like spin chains. By performing an exhaustive search over the spin configurations that can be enumerated on periodic 1D integer lattices containing $N=2,3,\ldots,36$ sites, we were able to identify and structurally characterize \textit{all} stealthy hyperuniform spin chains in this range of system sizes. Within this pool of stealthy hyperuniform spin configurations, we then utilized such inverse optimization techniques to demonstrate that stealthy hyperuniform spin chains can be realized as either unique or degenerate disordered ground states of radial long-ranged (relative to the spin chain length) spin-spin interactions. Such exotic ground states are distinctly different from spin glasses in both their inherent structural properties and the nature of the spin-spin interactions required to stabilize them. As such, the implications and significance of the existence of such disordered stealthy hyperuniform ground state spin systems warrants further study, including whether their bulk physical properties and excited states, like their many-particle system counterparts, are singularly remarkable, and can be experimentally realized. △ Less

Submitted 30 August, 2015; originally announced August 2015.

Comments: 11 pages, 9 figures

Journal ref: Phys. Rev. B 93, 064201 (2016)

Local Structure Analysis in $Ab$ $Initio$ Liquid Water

Authors: Biswajit Santra, Robert A. DiStasio Jr., Fausto Martelli, Roberto Car

Abstract: Within the framework of density functional theory, the inclusion of exact exchange and non-local van der Waals/dispersion (vdW) interactions is crucial for predicting a microscopic structure of ambient liquid water that quantitatively agrees with experiment. In this work, we have used the local structure index (LSI) order parameter to analyze the local structure in such highly accurate $ab$… ▽ More Within the framework of density functional theory, the inclusion of exact exchange and non-local van der Waals/dispersion (vdW) interactions is crucial for predicting a microscopic structure of ambient liquid water that quantitatively agrees with experiment. In this work, we have used the local structure index (LSI) order parameter to analyze the local structure in such highly accurate $ab$ $initio$ liquid water. At ambient conditions, the LSI probability distribution, P($I$), was unimodal with most water molecules characterized by more disordered high-density-like local environments. With thermal excitations removed, the resultant bimodal P($I$) in the inherent potential energy surface (IPES) exhibited a 3:1 ratio between high- and low-density-like molecules, with the latter forming small connected clusters amid the predominant population. By considering the spatial correlations and hydrogen bond network topologies $among$ water molecules with the same LSI identities, we demonstrate that the signatures of the experimentally observed low- (LDA) and high-density (HDA) amorphous phases of ice are present in the IPES of ambient liquid water. Analysis of the LSI autocorrelation function uncovered a persistence time of $\sim$ 4 ps---a finding consistent with the fact that natural thermal fluctuations are responsible for transitions between these distinct yet transient local aqueous environments in ambient liquid water. △ Less

Submitted 27 February, 2015; originally announced March 2015.

Comments: 12 pages, 6 figures

Journal ref: Molecular Physics 113, 829 (2015)

3D Dirac semimetals: current materials, design principles and predictions of new materials

Authors: Quinn D. Gibson, Leslie M. Schoop, Lukas Muechler, Lilia S. Xie, Maximillian Hirschberger, Nai Phuan Ong, Roberto Car, Robert J. Cava

Abstract: Design principles and novel predictions of new 3D Dirac semimetals are presented, along with the context of currently known materials. Current materials include those based on a topological to trivial phase transition, such as in TlBiSe$_{2-x}$S$_x$ and Hg$_{1-x}$Cd$_x$Te, Bi$_{1-x}$Sb$_x$, Bi$_{2-x}$In$_x$Se$_3$, and Pb$_{1-x}$Sn$_x$Se. Some more recently revealed materials, Na$_3$Bi and Cd$_3$As… ▽ More Design principles and novel predictions of new 3D Dirac semimetals are presented, along with the context of currently known materials. Current materials include those based on a topological to trivial phase transition, such as in TlBiSe$_{2-x}$S$_x$ and Hg$_{1-x}$Cd$_x$Te, Bi$_{1-x}$Sb$_x$, Bi$_{2-x}$In$_x$Se$_3$, and Pb$_{1-x}$Sn$_x$Se. Some more recently revealed materials, Na$_3$Bi and Cd$_3$As$_2$, are 3D Dirac semimetals in their native composition. The different design principles presented each yield novel predictions for new candidates. For Case I, 3D Dirac semimetals based on charge balanced compounds, BaAgBi, SrAgBi, YbAuSb, PtBi$_2$ and SrSn$_2$As$_2$ are identified as candidates. For Case II, 3D Dirac semi-metals in analogy to graphene, BaGa$_2$ is identified as a candidate, and BaPt and Li$_2$Pt are discussed. For Case III, 3D Dirac semi-metals based on glide planes and screw axes, TlMo$_3$Te$_3$ and the AMo$_3$X$_3$ family in general (A=K, Na, In, Tl, X=Se,Te) as well as the Group IVb trihalides such as HfI$_3$ are identified as candidates. Finally we discuss conventional intermetallic compounds with Dirac cones, and identify Cr$_2$B as a potentially interesting material. △ Less

Submitted 31 October, 2014; originally announced November 2014.

Comments: submitted 10/30

Möbius molecules and fragile Mott insulators

Authors: Lukas Muechler, Joseph Maciejko, Titus Neupert, Roberto Car

Abstract: Motivated by the concept of Möbius aromatics in organic chemistry, we extend the recently introduced concept of fragile Mott insulators (FMI) to ring-shaped molecules with repulsive Hubbard interactions threaded by a half-quantum of magnetic flux ($hc/2e$). In this context, a FMI is the insulating ground state of a finite-size molecule that cannot be adiabatically connected to a single Slater dete… ▽ More Motivated by the concept of Möbius aromatics in organic chemistry, we extend the recently introduced concept of fragile Mott insulators (FMI) to ring-shaped molecules with repulsive Hubbard interactions threaded by a half-quantum of magnetic flux ($hc/2e$). In this context, a FMI is the insulating ground state of a finite-size molecule that cannot be adiabatically connected to a single Slater determinant, i.e., to a band insulator, provided that time-reversal and lattice translation symmetries are preserved. Based on exact numerical diagonalization for finite Hubbard interaction strength $U$ and existing Bethe-ansatz studies of the one-dimensional Hubbard model in the large-$U$ limit, we establish a duality between Hubbard molecules with $4n$ and $4n+2$ sites, with $n$ integer. A molecule with $4n$ sites is an FMI in the absence of flux but becomes a band insulator in the presence of a half-quantum of flux, while a molecule with $4n+2$ sites is a band insulator in the absence of flux but becomes an FMI in the presence of a half-quantum of flux. Including next-nearest-neighbor-hoppings gives rise to new FMI states that belong to multidimensional irreducible representations of the molecular point group, giving rise to a rich phase diagram. △ Less

Submitted 23 September, 2014; originally announced September 2014.

Journal ref: Phys. Rev. B 90, 245142 (2014)

Response to Comment [arXiv:1407.6854] on Palmer et al., Nature, 510, 385, 2014

Authors: Jeremy C. Palmer, Pablo G. Debenedetti, Roberto Car, Athanassios Z. Panagiotopoulos

Abstract: We respond to a Comment [arXiv:1407.6854 (2014)] on our recent Nature paper [Nature, 510, 385 (2014)]. We categorically disprove the arguments provided in arXiv:1407.6854 (2014) and thereby further substantiate the evidence we presented in our recent study, demonstrating the existence of a metastable liquid-liquid transition in a molecular model of water. We will make our code publicly available s… ▽ More We respond to a Comment [arXiv:1407.6854 (2014)] on our recent Nature paper [Nature, 510, 385 (2014)]. We categorically disprove the arguments provided in arXiv:1407.6854 (2014) and thereby further substantiate the evidence we presented in our recent study, demonstrating the existence of a metastable liquid-liquid transition in a molecular model of water. We will make our code publicly available shortly along with proper user documentation that is currently under development. △ Less

Submitted 29 July, 2014; originally announced July 2014.

Comments: 10 pages, 2 figures

The Individual and Collective Effects of Exact Exchange and Dispersion Interactions on the Ab Initio Structure of Liquid Water

Authors: Robert A. DiStasio Jr., Biswajit Santra, Zhaofeng Li, Xifan Wu, Roberto Car

Abstract: In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local vdW/dispersion interactions, via a fully self-consistent den… ▽ More In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local vdW/dispersion interactions, via a fully self-consistent density-dependent dispersion correction, and approximate nuclear quantum effects (aNQE), via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx, vdW, and aNQE as resulting from a large-scale AIMD simulation of (H$_2$O)$_{128}$ at the PBE0+vdW level of theory, significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, $S_{\rm OO}(Q)$, and corresponding oxygen-oxygen radial distribution function, $g_{\rm OO}(r)$, that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment as demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of the PBE0 hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. △ Less

Submitted 20 May, 2014; originally announced May 2014.

Journal ref: J. Chem. Phys. 141, 084502 (2014)

A favorably-scaling natural-orbital functional theory based on higher-order occupation probabilities

Authors: Ralph Gebauer, Morrel H. Cohen, Roberto Car

Abstract: We introduce a novel energy functional for ground-state electronic-structure calculations. Its fundamental variables are the natural spin-orbitals of the implied singlet many-body wave function and their joint occupation probabilities. The functional derives from a sequence of controlled approximations to the two-particle density matrix. Algebraic scaling of computational cost with electron number… ▽ More We introduce a novel energy functional for ground-state electronic-structure calculations. Its fundamental variables are the natural spin-orbitals of the implied singlet many-body wave function and their joint occupation probabilities. The functional derives from a sequence of controlled approximations to the two-particle density matrix. Algebraic scaling of computational cost with electron number is obtainable in general, and Hartree-Fock scaling in the seniority-zero version of the theory. Results obtained with the latter version for saturated small molecular systems are compared with those of highly-accurate quantum-chemical computations. The numerical results are variational, capturing most of the correlation energy from equilibrium to dissociation. Their accuracy is considerably greater than that obtainable with current density-functional theory approximations and with current functionals of the one-particle density matrix only. △ Less

Submitted 20 June, 2015; v1 submitted 16 September, 2013; originally announced September 2013.

On the Accuracy of van der Waals Inclusive Density-Functional Theory Exchange-Correlation Functionals for Ice at Ambient and High Pressures

Authors: Biswajit Santra, Jiří Klimeš, Alexandre Tkatchenko, Dario Alfè, Ben Slater, Angelos Michaelides, Roberto Car, Matthias Scheffler

Abstract: Density-functional theory (DFT) has been widely used to study water and ice for at least 20 years. However, the reliability of different DFT exchange-correlation (xc) functionals for water remains a matter of considerable debate. This is particularly true in light of the recent development of DFT based methods that account for van der Waals (vdW) dispersion forces. Here, we report a detailed study… ▽ More Density-functional theory (DFT) has been widely used to study water and ice for at least 20 years. However, the reliability of different DFT exchange-correlation (xc) functionals for water remains a matter of considerable debate. This is particularly true in light of the recent development of DFT based methods that account for van der Waals (vdW) dispersion forces. Here, we report a detailed study with several xc functionals (semi-local, hybrid, and vdW inclusive approaches) on ice Ih and six proton ordered phases of ice. Consistent with our previous study [Phys. Rev. Lett. 107, 185701 (2011)] which showed that vdW forces become increasingly important at high pressures, we find here that all vdW inclusive methods considered improve the relative energies and transition pressures of the high-pressure ice phases compared to those obtained with semi-local or hybrid xc functionals. However, we also find that significant discrepancies between experiment and the vdW inclusive approaches remain in the cohesive properties of the various phases, causing certain phases to be absent from the phase diagram. Therefore, room for improvement in the description of water at ambient and high pressures remains and we suggest that because of the stern test the high pressure ice phases pose they should be used in future benchmark studies of simulation methods for water. △ Less

Submitted 13 August, 2014; v1 submitted 12 July, 2013; originally announced July 2013.

Comments: 13 pages, 5 figures, 4 tables

Journal ref: J. Chem. Phys. 139, 154702 (2013)