BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Roberto Car

Roberto Car appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phphysics.comp-phcond-mat.softcond-mat.stat-mechcond-mat.str-elcond-mat.dis-nncond-mat.otherquant-phcond-mat.mes-hallmath.NAmtrl-th

Catalog footprint

What is connected

37works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2025arXiv

Ab Initio Melting Properties of Water and Ice from Machine Learning Potentials

Liquid water exhibits several important anomalous properties in the vicinity of the melting temperature ($T_{\mathrm{m}}$) of ice Ih, including a higher density than ice and a density maximum at 4~$^{\circ}$C. Experimentally, an isotope effect on $T_{\mathrm{m}}$ is observed: the melting temperature of H$_2$O is approximately 4~K lower than that of D$_2$O. This difference can only be explained by nuclear quantum effects (NQEs), which can be accurately captured using path integral molecular dynamics (PIMD). Here we run PIMD simulations driven by Deep Potential (DP) models trained on data from density functional theory (DFT) based on SCAN, revPBE0-D3, SCAN0, and revPBE-D3 and a DP model trained on the MB-pol potential. We calculate the \tm of ice, the density discontinuity at melting, and the temperature of density maximum ($T_{\mathrm{dm}}$) of the liquid. We find that the model based on MB-pol agrees well with experiment. The models based on DFT incorrectly predict that NQEs lower $T_{\mathrm{m}}$. For the density discontinuity, SCAN and SCAN0 predict values close to the experimental result, while revPBE-D3 and revPBE0-D3 significantly underestimate it. Additionally, the models based on SCAN and SCAN0 correctly predict that the $T_{\mathrm{dm}}$ is higher than $T_{\mathrm{m}}$, while those based on revPBE-D3 and revPBE0-D3 predict the opposite. We attribute the deviations of the DFT-based models from experiment to the overestimation of hydrogen bond strength. Our results set the stage for more accurate simulations of aqueous systems grounded on DFT.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Assessment of First-Principles Methods in Modeling the Melting Properties of Water

First-principles simulations have played a crucial role in deepening our understanding of the thermodynamic properties of water, and machine learning potentials (MLPs) trained on these first-principles data widen the range of accessible properties. However, the capabilities of different first-principles methods are not yet fully understood due to the lack of systematic benchmarks, the underestimation of the uncertainties introduced by MLPs, and the neglect of nuclear quantum effects (NQEs). Here, we systematically assess first-principles methods by calculating key melting properties using path integral molecular dynamics (PIMD) driven by Deep Potential (DP) models trained on data from density functional theory (DFT) with SCAN, revPBE0-D3, SCAN0 and revPBE-D3 functionals, as well as from the MB-pol potential. We find that MB-pol is in qualitatively good agreement with the experiment in all properties tested, whereas the four DFT functionals incorrectly predict that NQEs increase the melting temperature. SCAN and SCAN0 slightly underestimate the density change between water and ice upon melting, but revPBE-D3 and revPBE0-D3 severely underestimate it. Moreover, SCAN and SCAN0 correctly predict that the maximum liquid density occurs at a temperature higher than the melting point, while revPBE-D3 and revPBE0-D3 predict the opposite behavior. Our results highlight limitations in widely used first-principles methods and call for a reassessment of their predictive power in aqueous systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A deep potential model with long-range electrostatic interactions

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deep neural network for the dielectric response of insulators

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pushing the limit of molecular dynamics with ab initio accuracy to 100 million atoms with machine learning

For 35 years, {\it ab initio} molecular dynamics (AIMD) has been the method of choice for modeling complex atomistic phenomena from first principles. However, most AIMD applications are limited by computational cost to systems with thousands of atoms at most. We report that a machine learning-based simulation protocol (Deep Potential Molecular Dynamics), while retaining {\it ab initio} accuracy, can simulate more than 1 nanosecond-long trajectory of over 100 million atoms per day, using a highly optimized code (GPU DeePMD-kit) on the Summit supercomputer. Our code can efficiently scale up to the entire Summit supercomputer, attaining $91$ PFLOPS in double precision ($45.5\%$ of the peak) and {$162$/$275$ PFLOPS in mixed-single/half precision}. The great accomplishment of this work is that it opens the door to simulating unprecedented size and time scales with {\it ab initio} accuracy. It also poses new challenges to the next-generation supercomputer for a better integration of machine learning and physical modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Raman Spectrum and Polarizability of Liquid Water from Deep Neural Networks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Deep Potential: a general representation of a many-body potential energy surface

We present a simple, yet general, end-to-end deep neural network representation of the potential energy surface for atomic and molecular systems. This methodology, which we call Deep Potential, is "first-principle" based, in the sense that no ad hoc approximations or empirical fitting functions are required. The neural network structure naturally respects the underlying symmetries of the systems. When tested on a wide variety of examples, Deep Potential is able to reproduce the original model, whether empirical or quantum mechanics based, within chemical accuracy. The computational cost of this new model is not substantially larger than that of empirical force fields. In addition, the method has promising scalability properties. This brings us one step closer to being able to carry out molecular simulations with accuracy comparable to that of quantum mechanics models and computational cost comparable to that of empirical potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A well-scaling natural orbital theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Performance of a Nonempirical Density Functional on Molecules and Hydrogen-Bonded Complexes

Recently, Tao and Mo (TM) derived a meta-generalized gradient approximation functional based on a model exchange-correlation hole. In this work, the performance of this functional is assessed on standard test sets, using the 6-311++G(3df,3pd) basis set. These test sets include 223 G3/99 enthalpies of formation, 99 atomization energies, 76 barrier heights, 58 electron affinities, 8 proton affinities, 96 bond lengths, 82 harmonic vibrational frequencies, 10 hydrogen-bonded molecular complexes, and 22 atomic excitation energies. Our calculations show that the TM functional can achieve high accuracy for most properties considered, relative to the LSDA, PBE, and TPSS functionals. In particular, it yields the best accuracy for proton affinities, harmonic vibrational frequencies, hydrogen-bonded dissociation energies and bond lengths, and atomic excitation energies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tilted Dirac Fermions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

Accurate treatment of the long-range electron correlation energy, including van der Waals (vdW) or dispersion interactions, is essential for describing the structure, dynamics, and function of a wide variety of systems. Among the most accurate models for including dispersion into density functional theory (DFT) is the range-separated many-body dispersion (MBD) method [A. Ambrossetti et al., J. Chem. Phys. 140, 18A508 (2014)], in which the correlation energy is modeled at short-range by a semi-local density functional and at long-range by a model system of coupled quantum harmonic oscillators. In this work, we develop analytical gradients of the MBD energy with respect to nuclear coordinates, including all implicit coordinate dependencies arising from the partitioning of the charge density into Hirshfeld effective volumes. To demonstrate the efficiency and accuracy of these MBD gradients for geometry optimizations of systems with intermolecular and intramolecular interactions, we optimized conformers of the benzene dimer and isolated small peptides with aromatic side-chains. We find excellent agreement with the wavefunction theory reference geometries of these systems (at a fraction of the computational cost) and find that MBD consistently outperforms the popular TS and D3(BJ) dispersion corrections. To demonstrate the performance of the MBD model on a larger system with supramolecular interactions, we optimized the C$_{60}$@C$_{60}$H$_{28}$ buckyball catcher host-guest complex. Finally, we find that neglecting the implicit nuclear coordinate dependence arising from the charge density partitioning, as has been done in prior numerical treatments, leads to an unacceptable error in the MBD forces, with relative errors of ~20% (on average) that can extend well beyond 100%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverse design of disordered stealthy hyperuniform spin chains

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Möbius molecules and fragile Mott insulators

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Role of quantum nuclei and local fields in the x-ray absorption spectra of water and ice

We calculate the x-ray absorption spectra of liquid water at ambient conditions and of hexagonal ice close to melting, using a static GW approach that includes approximately local field effects. Quantum dynamics of the nuclei is taken into account by averaging the absorption cross section over molecular configurations generated by path integral simulations. We find that inclusion of quantum disorder is essential to bring the calculated spectra in close agreement with experiment. In particular, the intensity of the pre-edge feature, a spectral signature of broken and distorted hydrogen bonds, is accurately reproduced, in water and ice, only when quantum nuclei are considered. The effect of the local fields is less important but non negligible, particularly in ice.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrogen bonds and van der Waals forces in ice at ambient and high pressures

The first principles approaches, density functional theory (DFT) and quantum Monte Carlo, have been used to examine the balance between van der Waals (vdW) forces and hydrogen (H) bonding in ambient and high pressure phases of ice. At higher pressure, the contribution to the lattice energy from vdW increases and that from H bonding decreases, leading vdW to have a substantial effect on the transition pressures between the crystalline ice phases. An important consequence, likely to be of relevance to molecular crystals in general, is that transition pressures obtained from DFT functionals which neglect vdW forces are greatly overestimated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Momentum distribution, vibrational dynamics and the potential of mean force in ice

By analyzing the momentum distribution obtained from path integral and phonon calculations we find that the protons in hexagonal ice experience an anisotropic quasi-harmonic effective potential with three distinct principal frequencies that reflect molecular orientation. Due to the importance of anisotropy, anharmonic features of the environment cannot be extracted from existing experimental distributions that involve the spherical average. The full directional distribution is required, and we give a theoretical prediction for this quantity that could be verified in future experiments. Within the quasi-harmonic context, anharmonicity in the ground state dynamics of the proton is substantial and has quantal origin, a finding that impacts the interpretation of several spectroscopies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Displaced path integral formulation for the momentum distribution of quantum particles

The proton momentum distribution, accessible by deep inelastic neutron scattering, is a very sensitive probe of the potential of mean force experienced by the protons in hydrogen-bonded systems. In this work we introduce a novel estimator for the end to end distribution of the Feynman paths, i.e. the Fourier transform of the momentum distribution. In this formulation, free particle and environmental contributions factorize. Moreover, the environmental contribution has a natural analogy to a free energy surface in statistical mechanics, facilitating the interpretation of experiments. The new formulation is not only conceptually but also computationally advantageous. We illustrate the method with applications to an empirical water model, ab-initio ice, and one dimensional model systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The isotope-effect in the phase transition of KDP: New insights from ab initio path-integral simulations

We investigate the quantum-mechanical localization of protonated and deterated isotopes in the symmetric low-barrier hydrogen-bonds of potassium dihydrogen phosphate (KDP) crystals in the paraelectric phase. The spatial density distributions of these hydrogen atoms are suspected to be responsible for the surprisingly large isotope effect observed for the ferroelectric phase transition in KDP. We employ ab initio path integral molecular dynamics simulations to obtain the nuclear real-space and momentum-space densities n(R) and n(k) of protons and deuterons, which are compared to experimental Neutron Compton Scattering data. Our results suggest a qualitative difference in the nature of the paraelectic phase in KDP between the two isotopes. We are able to discriminate between real quantum delocalization and vibration-assisted hopping and thus provide evidence for two distinct mechanisms of the ferroelectric phase transition in this class of materials.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hybrid density functional calculations of the band gap of Ga$_x$In$_{1-x}$N

Recent theoretical work has provided evidence that hybrid functionals, which include a fraction of exact (Hartree Fock) exchange in the density functional theory (DFT) exchange and correlation terms, significantly improve the description of band gaps of semiconductors compared with local and semilocal approximations. Based on a recently developed order-$N$ method for calculating the exact exchange in extended insulating systems, we have implemented an efficient scheme to determine the hybrid functional band gap. We use this scheme to study the band gap and other electronic properties of the ternary compound In$_{1-x}$Ga$_{x}$N using a 64-atom supercell model.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Tunneling and delocalization in hydrogen bonded systems: a study in position and momentum space

Novel experimental and computational studies have uncovered the proton momentum distribution in hydrogen bonded systems. In this work, we utilize recently developed open path integral Car-Parrinello molecular dynamics methodology in order to study the momentum distribution in phases of high pressure ice. Some of these phases exhibit symmetric hydrogen bonds and quantum tunneling. We find that the symmetric hydrogen bonded phase possesses a narrowed momentum distribution as compared with a covalently bonded phase, in agreement with recent experimental findings. The signatures of tunneling that we observe are a narrowed distribution in the low-to-intermediate momentum region, with a tail that extends to match the result of the covalently bonded state. The transition to tunneling behavior shows similarity to features observed in recent experiments performed on confined water. We corroborate our ice simulations with a study of a particle in a model one-dimensional double well potential that mimics some of the effects observed in bulk simulations. The temperature dependence of the momentum distribution in the one-dimensional model allows for the differentiation between ground state and mixed state tunneling effects.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

X-ray absorption signatures of the molecular environment in water and ice

The x-ray absorption spectra of water and ice are calculated with a many-body approach for electron-hole excitations. The experimental features, including the small effects of temperature change in the liquid, are quantitatively reproduced from molecular configurations generated by ab-initio molecular dynamics. The spectral difference between the solid and the liquid is due to two major short range order effects. One, due to breaking of hydrogen bonds, enhances the pre-edge intensity in the liquid. The other, due to a non-bonded molecular fraction in the first coordination shell, affects the main spectral edge in the conversion of ice to water. This effect may not involve hydrogen bond breaking as shown by experiment in high-density amorphous ice.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Charge Transfer in Partition Theory

The recently proposed Partition Theory (PT) [J.Phys.Chem.A 111, 2229 (2007)] is illustrated on a simple one-dimensional model of a heteronuclear diatomic molecule. It is shown that a sharp definition for the charge of molecular fragments emerges from PT, and that the ensuing population analysis can be used to study how charge redistributes during dissociation and the implications of that redistribution for the dipole moment. Interpreting small differences between the isolated parts' ionization potentials as due to environmental inhomogeneities, we gain insight into how electron localization takes place in H2+ as the molecule dissociates. Furthermore, by studying the preservation of the shapes of the parts as different parameters of the model are varied, we address the issue of transferability of the parts. We find good transferability within the chemically meaningful parameter regime, raising hopes that PT will prove useful in chemical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Multipole Representation of the Fermi Operator with Application to the Electronic Structure Analysis of Metallic Systems

We propose a multipole representation of the Fermi-Dirac function and the Fermi operator, and use this representation to develop algorithms for electronic structure analysis of metallic systems. The new algorithm is quite simple and efficient. Its computational cost scales logarithmically with $βΔ\eps$ where $β$ is the inverse temperature, and $Δ\eps$ is the width of the spectrum of the discretized Hamiltonian matrix.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

First-principle study of excitonic self-trapping in diamond

We present a first-principles study of excitonic self-trapping in diamond. Our calculation provides evidence for self-trapping of the 1s core exciton and gives a coherent interpretation of recent experimental X-ray absorption and emission data. Self-trapping does not occur in the case of a single valence exciton. We predict, however, that self-trapping should occur in the case of a valence biexciton. This process is accompanied by a large local relaxation of the lattice which could be observed experimentally.

0 citations0 reviewsNo code linkedOpen access