-
Robust linear-scaling optimization of compact localized orbitals in density functional theory
Authors:
Yifei Shi,
Jessica Karaguesian,
Rustam Z. Khaliullin
Abstract:
Locality of compact one-electron orbitals expanded strictly in terms of local subsets of basis functions can be exploited in density functional theory (DFT) to achieve linear growth of computation time with systems size, crucial in large-scale simulations. However, despite advantages of compact orbitals the development of practical orbital-based linear-scaling DFT methods has long been hindered be…
▽ More
Locality of compact one-electron orbitals expanded strictly in terms of local subsets of basis functions can be exploited in density functional theory (DFT) to achieve linear growth of computation time with systems size, crucial in large-scale simulations. However, despite advantages of compact orbitals the development of practical orbital-based linear-scaling DFT methods has long been hindered because a compact representation of the electronic ground state is difficult to find in a variational optimization procedure. In this work, we show that the slow and unstable optimization of compact orbitals originates from the nearly-invariant mixing of compact orbitals that are mostly but not completely localized within the same subsets of basis functions. We also construct an approximate Hessian that can be used to identify the problematic nearly-invariant modes and obviate the variational optimization along them without introducing significant errors into the computed energies. This enables us to create a linear-scaling DFT method with a low computational overhead that is demonstrated to be efficient and accurate in fixed-nuclei calculations and molecular dynamics simulations of semiconductors and insulators.
△ Less
Submitted 30 September, 2021; v1 submitted 9 April, 2020;
originally announced April 2020.
-
CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
Authors:
Thomas D. Kühne,
Marcella Iannuzzi,
Mauro Del Ben,
Vladimir V. Rybkin,
Patrick Seewald,
Frederick Stein,
Teodoro Laino,
Rustam Z. Khaliullin,
Ole Schütt,
Florian Schiffmann,
Dorothea Golze,
Jan Wilhelm,
Sergey Chulkov,
Mohammad Hossein Bani-Hashemian,
Valéry Weber,
Urban Borstnik,
Mathieu Taillefumier,
Alice Shoshana Jakobovits,
Alfio Lazzaro,
Hans Pabst,
Tiziano Müller,
Robert Schade,
Manuel Guidon,
Samuel Andermatt,
Nico Holmberg
, et al. (14 additional authors not shown)
Abstract:
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achiev…
▽ More
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2k to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
△ Less
Submitted 11 March, 2020; v1 submitted 8 March, 2020;
originally announced March 2020.
-
On the role of interfacial hydrogen bonds in "on-water" catalysis
Authors:
Kristof Karhan,
Rustam Z. Khaliullin,
Thomas D. Kühne
Abstract:
Numerous experiments have demonstrated that many classes of organic reactions exhibit increased reaction rates when performed in heterogeneous water emulsions. Despite enormous practical importance of the observed "on-water" catalytic effect and several mechanistic studies, its microscopic origins remains unclear. In this work, the second generation Car-Parrinello molecular dynamics method is exte…
▽ More
Numerous experiments have demonstrated that many classes of organic reactions exhibit increased reaction rates when performed in heterogeneous water emulsions. Despite enormous practical importance of the observed "on-water" catalytic effect and several mechanistic studies, its microscopic origins remains unclear. In this work, the second generation Car-Parrinello molecular dynamics method is extended to self-consistent charge density-functional based tight-binding in order to study "on-water" catalysis of the Diels-Alder reaction between dimethyl azodicarboxylate and quadricyclane. We find that the stabilization of the transition state by dangling hydrogen bonds exposed at the aqueous interfaces plays a significantly smaller role in "on-water" catalysis than has been suggested previously.
△ Less
Submitted 21 August, 2014;
originally announced August 2014.
-
Recent achievements in ab initio modelling of liquid water
Authors:
Rustam Z. Khaliullin,
Thomas D. Kühne
Abstract:
The application of newly developed first-principle modeling techniques to liquid water deepens our understanding of the microscopic origins of its unusual macroscopic properties and behaviour. Here, we review two novel ab initio computational methods: second-generation Car-Parrinello molecular dynamics and decomposition analysis based on absolutely localized molecular orbitals. We show that these…
▽ More
The application of newly developed first-principle modeling techniques to liquid water deepens our understanding of the microscopic origins of its unusual macroscopic properties and behaviour. Here, we review two novel ab initio computational methods: second-generation Car-Parrinello molecular dynamics and decomposition analysis based on absolutely localized molecular orbitals. We show that these two methods in combination not only enable ab initio molecular dynamics simulations on previously inaccessible time and length scales, but also provide unprecedented insights into the nature of hydrogen bonding between water molecules. We discuss recent applications of these methods to water clusters and bulk water.
△ Less
Submitted 8 March, 2013;
originally announced March 2013.
-
Electronic signature of the instantaneous asymmetry in the first coordination shell of liquid water
Authors:
Thomas D. Kühne,
Rustam Z. Khaliullin
Abstract:
Interpretation of the X-ray spectra of water as evidence for its asymmetric structure has challenged the conventional symmetric nearly-tetrahedral model and initiated an intense debate about the order and symmetry of the hydrogen bond network in water. Here, we present new insights into the nature of local interactions in water obtained using a novel energy decomposition method. Our simulations re…
▽ More
Interpretation of the X-ray spectra of water as evidence for its asymmetric structure has challenged the conventional symmetric nearly-tetrahedral model and initiated an intense debate about the order and symmetry of the hydrogen bond network in water. Here, we present new insights into the nature of local interactions in water obtained using a novel energy decomposition method. Our simulations reveal that while a water molecule forms, on average, two strong donor and two strong acceptor bonds, there is a significant asymmetry in the energy of these contacts. We demonstrate that this asymmetry is a result of small instantaneous distortions of hydrogen bonds, which appear as fluctuations on a timescale of hundreds of femtoseconds around the average symmetric structure. Furthermore, we show that the distinct features of the X-ray absorption spectra originate from molecules with high instantaneous asymmetry. Our findings have important implications as they help reconcile the symmetric and asymmetric views on the structure of water.
△ Less
Submitted 19 January, 2013;
originally announced January 2013.
-
Microscopic origins of the anomalous melting behaviour of high-pressure sodium
Authors:
Hagai Eshet,
Rustam Z. Khaliullin,
Thomas D. Kühne,
Jörg Behler,
Michele Parrinello
Abstract:
Recent experiments have shown that sodium, a prototype simple metal at ambient conditions, exhibits unexpected complexity under high pressure. One of the most puzzling phenomena in the behaviour of dense sodium is the pressure-induced drop in its melting temperature, which extends from 1000 K at ~30GPa to as low as room temperature at ~120GPa. Despite significant theoretical effort to understand t…
▽ More
Recent experiments have shown that sodium, a prototype simple metal at ambient conditions, exhibits unexpected complexity under high pressure. One of the most puzzling phenomena in the behaviour of dense sodium is the pressure-induced drop in its melting temperature, which extends from 1000 K at ~30GPa to as low as room temperature at ~120GPa. Despite significant theoretical effort to understand the anomalous melting its origins have remained unclear. In this work, we reconstruct the sodium phase diagram using an ab-initio-quality neural-network potential. We demonstrate that the reentrant behaviour results from the screening of interionic interactions by conduction electrons, which at high pressure induces a softening in the short-range repulsion. It is expected that such an effect plays an important role in governing the behaviour of a wide range of metals and alloys.
△ Less
Submitted 10 October, 2011;
originally announced October 2011.
-
Nucleation mechanism for the direct graphite-to-diamond phase transition
Authors:
Rustam Z. Khaliullin,
Hagai Eshet,
Thomas D. Kuhne,
Jorg Behler,
Michele Parrinello
Abstract:
Graphite and diamond have comparable free energies, yet forming diamond from graphite is far from easy. In the absence of a catalyst, pressures that are significantly higher than the equilibrium coexistence pressures are required to induce the graphite-to-diamond transition. Furthermore, the formation of the metastable hexagonal polymorph of diamond instead of the more stable cubic diamond is favo…
▽ More
Graphite and diamond have comparable free energies, yet forming diamond from graphite is far from easy. In the absence of a catalyst, pressures that are significantly higher than the equilibrium coexistence pressures are required to induce the graphite-to-diamond transition. Furthermore, the formation of the metastable hexagonal polymorph of diamond instead of the more stable cubic diamond is favored at lower temperatures. The concerted mechanism suggested in previous theoretical studies cannot explain these phenomena. Using an ab initio quality neural-network potential we performed a large-scale study of the graphite-to-diamond transition assuming that it occurs via nucleation. The nucleation mechanism accounts for the observed phenomenology and reveals its microscopic origins. We demonstrated that the large lattice distortions that accompany the formation of the diamond nuclei inhibit the phase transition at low pressure and direct it towards the hexagonal diamond phase at higher pressure. The nucleation mechanism proposed in this work is an important step towards a better understanding of structural transformations in a wide range of complex systems such as amorphous carbon and carbon nanomaterials.
△ Less
Submitted 7 January, 2011;
originally announced January 2011.
-
Ab initio quality study of the graphite-diamond phase coexistence
Authors:
Rustam Z. Khaliullin,
Hagai Eshet,
Thomas D. Kühne,
Jörg Behler,
Michele Parrinello
Abstract:
An interatomic potential for the diamond and graphite phases of carbon has been created using a neural-network (NN) representation of the ab initio potential energy surface. The NN potential combines the accuracy of a first-principle description of both phases with the efficiency of empirical force fields and allows one to perform, for the first time, a molecular dynamics study, of ab initio qua…
▽ More
An interatomic potential for the diamond and graphite phases of carbon has been created using a neural-network (NN) representation of the ab initio potential energy surface. The NN potential combines the accuracy of a first-principle description of both phases with the efficiency of empirical force fields and allows one to perform, for the first time, a molecular dynamics study, of ab initio quality, of the thermodynamics of graphite-diamond coexistence. Good agreement between the experimental and calculated coexistence curves is achieved if nuclear quantum effects are included in the simulation.
△ Less
Submitted 20 April, 2010; v1 submitted 15 February, 2010;
originally announced February 2010.
-
Ab initio quality neural-network potential for sodium
Authors:
Hagai Eshet,
Rustam Z. Khaliullin,
Thomas D. Kuhne,
Jorg Behler,
Michele Parrinello
Abstract:
An interatomic potential for high-pressure high-temperature (HPHT) crystalline and liquid phases of sodium is created using a neural-network (NN) representation of the ab initio potential energy surface. It is demonstrated that the NN potential provides an ab initio quality description of multiple properties of liquid sodium and bcc, fcc, cI16 crystal phases in the P-T region up to 120 GPa and 1…
▽ More
An interatomic potential for high-pressure high-temperature (HPHT) crystalline and liquid phases of sodium is created using a neural-network (NN) representation of the ab initio potential energy surface. It is demonstrated that the NN potential provides an ab initio quality description of multiple properties of liquid sodium and bcc, fcc, cI16 crystal phases in the P-T region up to 120 GPa and 1200 K. The unique combination of computational efficiency of the NN potential and its ability to reproduce quantitatively experimental properties of sodium in the wide P-T range enables molecular dynamics simulations of physicochemical processes in HPHT sodium of unprecedented quality.
△ Less
Submitted 19 April, 2010; v1 submitted 15 February, 2010;
originally announced February 2010.