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

Thomas D. Kühne

Thomas D. Kühne appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.comp-phphysics.chem-phcond-mat.mtrl-scicond-mat.stat-mechcond-mat.softcond-mat.dis-nnquant-phcond-mat.str-elBiological Physicscond-mat.othercond-mat.supr-conDistributed, Parallel, and Cluster Computingmath.NAastro-ph.EPcond-mat.mes-hallmath.RANumerical AnalysisPerformancephysics.class-phQuantitative Methods

Catalog footprint

What is connected

38works
20topics
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

38 published item(s)

preprint2026arXiv

Hitchhiker's guide to second-generation Car-Parrinello ab-initio molecular dynamics

In a recent letter [T. D. Kühne, M. Krack, F. Mohamed and M. Parrinello, Phys. Rev. Lett. 98, 066401 (2007)], we outlined a new Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics. Here, we provide a guide to performing actual calculations using our method and demonstrate this on liquid water at ambient conditions. We do not go into methodological details beyond those necessary for applying this approach, but focus on practical details pertinent to our particular implementation within the CP2K/Quickstep code [T. D. Kühne et al., J. Chem. Phys. 152, 194103 (2020)].

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

The CP-PAW code package for first-principles calculations from a user's perspective

CP-PAW is a combined electronic structure and ab-initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave method with the Car-Parrinello approach to determine not only the electronic and nuclear ground state of condensed matter, but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Breaking the Exascale Barrier for the Electronic Structure Problem in Ab-Initio Molecular Dynamics

The non-orthogonal local submatrix method applied to electronic-structure based molecular dynamics simulations is shown to exceed 1.1 EFLOP/s in FP16/FP32 mixed floating-point arithmetic when using 4,400 NVIDIA A100 GPUs of the Perlmutter system. This is enabled by a modification of the original method that pushes the sustained fraction of the peak performance to about 80%. Example calculations are performed for SARS-CoV-2 spike proteins with up to 83 million atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CP2K on the road to exascale

The CP2K program package, which can be considered as the swiss army knife of atomistic simulations, is presented with a special emphasis on ab-initio molecular dynamics using the second-generation Car-Parrinello method. After outlining current and near-term development efforts with regards to massively parallel low-scaling post-Hartree-Fock and eigenvalue solvers, novel approaches on how we plan to take full advantage of future low-precision hardware architectures are introduced. Our focus here is on combining our submatrix method with the approximate computing paradigm to address the immanent exascale era.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Parallel Quantum Chemistry on Noisy Intermediate-Scale Quantum Computers

A novel parallel hybrid quantum-classical algorithm for the solution of the quantum-chemical ground-state energy problem on gate-based quantum computers is presented. This approach is based on the reduced density-matrix functional theory (RDMFT) formulation of the electronic structure problem. For that purpose, the density-matrix functional of the full system is decomposed into an indirectly coupled sum of density-matrix functionals for all its subsystems using the adaptive cluster approximation to RDMFT. The approximations involved in the decomposition and the adaptive cluster approximation itself can be systematically converged to the exact result. The solutions for the density-matrix functionals of the effective subsystems involves a constrained minimization over many-particle states that are approximated by parametrized trial states on the quantum computer similarly to the variational quantum eigensolver. The independence of the density-matrix functionals of the effective subsystems introduces a new level of parallelization and allows for the computational treatment of much larger molecules on a quantum computer with a given qubit count. In addition, for the proposed algorithm techniques are presented to reduce the qubit count, the number of quantum programs, as well as its depth. The new approach is demonstrated for Hubbard-like systems on IBM quantum computers based on superconducting transmon qubits.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Towards Electronic Structure-Based Ab-Initio Molecular Dynamics Simulations with Hundreds of Millions of Atoms

We push the boundaries of electronic structure-based \textit{ab-initio} molecular dynamics (AIMD) beyond 100 million atoms. This scale is otherwise barely reachable with classical force-field methods or novel neural network and machine learning potentials. We achieve this breakthrough by combining innovations in linear-scaling AIMD, efficient and approximate sparse linear algebra, low and mixed-precision floating-point computation on GPUs, and a compensation scheme for the errors introduced by numerical approximations. The core of our work is the non-orthogonalized local submatrix method (NOLSM), which scales very favorably to massively parallel computing systems and translates large sparse matrix operations into highly parallel, dense matrix operations that are ideally suited to hardware accelerators. We demonstrate that the NOLSM method, which is at the center point of each AIMD step, is able to achieve a sustained performance of 324 PFLOP/s in mixed FP16/FP32 precision corresponding to an efficiency of 67.7% when running on 1536 NVIDIA A100 GPUs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An automated approach for developing neural network interatomic potentials with FLAME

The performance of machine learning interatomic potentials relies on the quality of the training dataset. In this work, we present an approach for generating diverse and representative training data points which initiates with \it{ab initio} calculations for bulk structures. The data generation and potential construction further proceed side-by-side in a cyclic process of training the neural network and crystal structure prediction based on the developed interatomic potentials. All steps of the data generation and potential development are performed with minimal human intervention. We show the reliability of our approach by assessing the performance of neural network potentials developed for two inorganic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Artificial Neural Networks as Trial Wave Functions for Quantum Monte Carlo

Inspired by the universal approximation theorem and widespread adoption of artificial neural network techniques in a diversity of fields, we propose feed-forward neural networks as a general purpose trial wave function for quantum Monte Carlo simulations of continous many-body systems. Whereas for simple model systems the whole many-body wave function can be represented by a neural network, the antisymmetry condition of non-trivial fermionic systems is incorporated by means of a Slater determinant. To demonstrate the accuracy of our trial wave functions, we have studied an exactly solvable model system of two trapped interacting particles, as well as the hydrogen dimer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accurate Sampling with Noisy Forces from Approximate Computing

In scientific computing, the acceleration of atomistic computer simulations by means of custom hardware is finding ever growing application. A major limitation, however, is that the high efficiency in terms of performance and low power consumption entails the massive usage of low-precision computing units. Here, based on the approximate computing paradigm, we present an algorithmic method to rigorously compensate for numerical inaccuracies due to low-accuracy arithmetic operations, yet still obtaining exact expectation values using a properly modified Langevin-type equation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Disordered Crystals from First Principles II: Transport Coefficients

This is the second part of a project on the foundations of first-principle calculations of the electron transport in crystals at finite temperatures, aiming at a predictive first-principles platform that combines ab-initio molecular dynamics (AIMD) and a finite-temperature Kubo-formula with dissipation for thermally disordered crystalline phases. The latter are encoded in an ergodic dynamical system $(Ω,\mathbb G,{\rm d}\mathbb P)$, where $Ω$ is the configuration space of the atomic degrees of freedom, $\mathbb G$ is the space group acting on $Ω$ and ${\rm d}\mathbb P$ is the ergodic Gibbs measure relative to the $\mathbb G$-action. We first demonstrate how to pass from the continuum Kohn-Sham theory to a discrete atomic-orbitals based formalism without breaking the covariance of the physical observables w.r.t. $(Ω,\mathbb G,{\rm d}\mathbb P)$. Then we show how to implement the Kubo-formula, investigate its self-averaging property and derive an optimal finite-volume approximation for it. We also describe a numerical innovation that made possible AIMD simulations with longer orbits and elaborate on the details of our simulations. Lastly, we present numerical results on the transport coefficients of crystal silicon at different temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient Ab-Initio Molecular Dynamic Simulations by Offloading Fast Fourier Transformations to FPGAs

A large share of today's HPC workloads is used for Ab-Initio Molecular Dynamics (AIMD) simulations, where the interatomic forces are computed on-the-fly by means of accurate electronic structure calculations. They are computationally intensive and thus constitute an interesting application class for energy-efficient hardware accelerators such as FPGAs. In this paper, we investigate the potential of offloading 3D Fast Fourier Transformations (FFTs) as a critical routine of plane-wave-based electronic structure calculations to FPGA and in conjunction demonstrate the tolerance of these simulations to lower precision computations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy Transfer within the Hydrogen Bonding Network of Water Following Resonant Terahertz Excitation

Energy dissipation in water is very fast and more efficient than in many other liquids. This behavior is commonly attributed to the intermolecular interactions associated with hydrogen bonding. Here, we investigate the dynamic energy flow in the hydrogen-bond network of liquid water by a pump-probe experiment. We resonantly excite intermolecular degrees of freedom with ultrashort single-cycle terahertz pulses and monitor its Raman response. By using ultrathin sample-cell windows, a background-free bipolar signal whose tail relaxes mono-exponentially is obtained. The relaxation is attributed to the molecular translational motions, using complementary experiments, force-field and ab initio molecular dynamics simulations. They reveal an initial coupling of the terahertz electric field to the molecular rotational degrees of freedom whose energy is rapidly transferred, within the excitation pulse duration, to the restricted-translational motion of neighboring molecules. This rapid energy transfer may be rationalized by the strong anharmonicity of the intermolecular interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modulation of Nearly Free Electron States in Hydroxyl-Functionalized MXenes: A First-Principles Study

The transition metal carbides (namely MXenes) and their functionalized derivatives exhibit various physical and chemical characteristics and offer many potential applications in electronic devices and sensors. Using density functional theory (DFT), it is revealed that the nearly free electron (NFE) states are near the Fermi levels in hydroxyl (OH) functionalized MXenes. Most of the OH-terminated MXene are metallic, but some of them, e.g. Sc2C(OH)2, are semiconductors and the NFE states are conduction bands. In this paper, to investigate the NFE states in MXenes, an attractive image-potential well model is adopted. Compared the solutions of this model with the DFT calculations, it is found that due to the overlap of spatially extensive wave functions of NFE states and their hybridization between the artificial neighboring layers imposed by the periodical boundary conditions (PBCs), the DFT results represent the properties of multiple layers, intrinsically. Based on the DFT calculations, it is found that the energy gap widths are affected by the interlayer distances. We address that the energetics of the NFE states can be modulated by the external electric fields and it is possible to convert semiconducting MXenes into metals. This band-gap manipulation makes the OH-terminated semiconducting MXenes an excellent candidate for electronic switch applications. Finally, using a set of electron transport calculations, I-V characteristics of Sc2C(OH)2 devices are investigated with the gate voltages.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tumbling with a Limp: Local Asymmetry in Water's Hydrogen Bond Network and its Consequences

Ab initio molecular dynamics simulations of liquid water under equilibrium ambient conditions, together with a novel energy decomposition analysis, have recently shown that a substantial fraction of water molecules exhibit a significant asymmetry between the strengths of the two donor and/or the two acceptor interactions. We refer to this recently unraveled aspect as the "local asymmetry in the hydrogen bond network". We discuss how this novel aspect was first revealed, and provide metrics that can be consistently employed on simulated water trajectories to quantify this local heterogeneity in the hydrogen bond network and its dynamics. We then discuss the static aspects of the asymmetry, pertaining to the frozen geometry of liquid water at any given instant of time and the distribution of hydrogen bond strengths therein, and also its dynamic characteristics pertaining to how fast this asymmetry decays and the kinds of molecular motions responsible for this decay. Following this we discuss the spectroscopic manifestations of this asymmetry, from ultrafast X-ray absorption spectra to infrared spectroscopy and down to the much slower terahertz regime. Finally, we discuss the implications of these findings in a broad context and its relation to the current notions about the structure and dynamics of liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Water structure near the surface of Weyl semimetals as catalysts in photocatalytic proton reduction

In this work, second-generation Car-Parrinello-based QM/MM molecular dynamics simulations of small nanoparticles of NbP, NbAs, TaAs and 1T-TaS$_2$ in water are presented. The first three materials are topological Weyl semimetals, which were recently discovered to be active catalysts in photocatalytic water splitting. The aim of this research was to correlate potential differences in the water structure in the vicinity of the nanoparticle surface with the photocatalytic activity of these materials in light induced proton reduction. The results presented herein allow to explain the catalytic activity of these Weyl semimetals: the most active material, NbP, exhibits a particularly low water coordination near the surface of the nanoparticle, whereas for 1T-TaS$_2$, with the lowest catalytic activity, the water structure at the surface is most ordered. In addition, the photocatalytic activity of several organic and metalorganic photosensitizers in the hydrogen evolution reaction was experimentally investigated with NbP as proton reduction catalyst. Unexpectedly, the charge of the photosensitizer plays a decisive role for the photocatalytic performance.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A General Algorithm to Calculate the Inverse Principal $p$-th Root of Symmetric Positive Definite Matrices

We address the general mathematical problem of computing the inverse $p$-th root of a given matrix in an efficient way. A new method to construct iteration functions that allow calculating arbitrary $p$-th roots and their inverses of symmetric positive definite matrices is presented. We show that the order of convergence is at least quadratic and that adaptively adjusting a parameter $q$ always leads to an even faster convergence. In this way, a better performance than with previously known iteration schemes is achieved. The efficiency of the iterative functions is demonstrated for various matrices with different densities, condition numbers and spectral radii.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Metal-Insulator Transition of Solid Hydrogen by the Antisymmetric Shadow Wave Function

We revisit the pressure-induced metal-insulator-transition of solid hydrogen by means of variational quantum Monte Carlo simulations based on the antisymmetric shadow wave function. In order to facilitate studying the electronic structure of large-scale fermionic systems, the shadow wave function formalism is extended by a series of technical improvements, such as a revised optimization method for the employed shadow wave function and an enhanced treatment of periodic systems with long-range interactions. It is found that the superior accuracy of the antisymmetric shadow wave function results in a significantly increased transition pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient "on-the-fly" calculation of Raman spectra from \textit{ab-initio} molecular dynamics: Application to hydrophobic/hydrophilic solutes in bulk water

We present a computational method to accurately calculate Raman spectra from first principles with an at least one order of magnitude higher efficiency. This scheme thus allows to routinely calculate finite-temperature Raman spectra "on-the-fly" by means of \textit{ab-initio} molecular dynamics simulations. To demonstrate the predictive power of this approach we investigate the effect of hydrophobic and hydrophilic solutes in water solution on the infrared and Raman spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Ring-Polymer Contraction Method: Including nuclear quantum effects at no additional computational cost in comparison to ab-initio molecular dynamics

We present a simple and accurate computational method, which facilitates ab-initio path-integral molecular dynamics simulations, where the quantum mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions. This development permits to routinely include nuclear quantum effects in ab-initio molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resonating Valence Bond Quantum Monte Carlo: Application to the ozone molecule

We study the potential energy surface of the ozone molecule by means of Quantum Monte Carlo simulations based on the resonating valence bond concept. The trial wave function consists of an antisymmetrized geminal power arranged in a single-determinant that is multiplied by a Jastrow correlation factor. Whereas the determinantal part incorporates static correlation effects, the augmented real-space correlation factor accounts for the dynamics electron correlation. The accuracy of this approach is demonstrated by computing the potential energy surface for the ozone molecule in three vibrational states: symmetric, asymmetric and scissoring. We find that the employed wave function provides a detailed description of rather strongly-correlated multi-reference systems, which is in quantitative agreement with experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structure and Dynamics of the Instantaneous Water/Vapor Interface Revisited by Path-Integral and Ab-Initio Molecular Dynamics Simulations

The structure and dynamics of the water/vapor interface is revisited by means of path-integral and second-generation Car-Parrinello ab-initio molecular dynamics simulations in conjunction with an instantaneous surface definition [A. P. Willard and D. Chandler, J. Phys. Chem. B 114, 1954 (2010)]. In agreement with previous studies, we find that one of the OH bonds of the water molecules in the topmost layer is pointing out of the water into the vapor phase, while the orientation of the underlying layer is reversed. Therebetween, an additional water layer is detected, where the molecules are aligned parallel to the instantaneous water surface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inverse simulated annealing: Improvements and application to amorphous InSb

An improved inverse simulated annealing method is presented to determine the structure of complex disordered systems from first principles in agreement with available experimental data or desired predetermined target properties. The effectiveness of this method is demonstrated by revisiting the structure of amorphous InSb. The resulting network is mostly tetrahedral and in excellent agreement with available experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear quantum effects in liquid water from path-integral simulations using an ab initio force matching approach

We have applied path integral simulations, in combination with new ab initio based water potentials, to investigate nuclear quantum effects in liquid water. Because direct ab initio path integral simulations are computationally expensive, a flexible water model is parameterized by force-matching to density functional theory-based molecular dynamics simulations. The resulting effective potentials provide an inexpensive replacement for direct ab inito molecular dynamics simulations and allow efficient simulation of nuclear quantum effects. Static and dynamic properties of liquid water at ambient conditions are presented and the role of nuclear quantum effects, exchange-correlation functionals and dispersion corrections are discussed in regards to reproducing the experimental properties of liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the role of interfacial hydrogen bonds in "on-water" catalysis

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Sign Problem of the Fermionic Shadow Wave Function

We present a whole series of novel methods to alleviate the sign problem of the Fermionic Shadow Wave Function in the context of Variational Monte Carlo. The effectiveness of our new techniques is demonstrated on the example of liquid 3He. We found that although the variance is substantially reduced, the gain in efficiency is restricted by the increased computational cost. Yet, this development not only extends the scope of the Fermionic Shadow Wave Function, but also facilitates highly accurate Quantum Monte Carlo simulations previously thought not feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-consistent field theory based molecular dynamics with linear system-size scaling

We present an improved field-theoretic approach to the grand-canonical potential suitable for linear scaling molecular dynamics simulations using forces from self-consistent electronic structure calculations. It is based on an exact decomposition of the grand canonical potential for independent fermions and does neither rely on the ability to localize the orbitals nor that the Hamilton operator is well-conditioned. Hence, this scheme enables highly accurate all-electron linear scaling calculations even for metallic systems. The inherent energy drift of Born-Oppenheimer molecular dynamics simulations, arising from an incomplete convergence of the self-consistent field cycle, is circumvented by means of a properly modified Langevin equation. The predictive power of the present linear scaling \textit{ab-initio} molecular dynamics approach is illustrated using the example of liquid methane under extreme conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ab-Initio Molecular Dynamics

Computer simulation methods, such as Monte Carlo or Molecular Dynamics, are very powerful computational techniques that provide detailed and essentially exact information on classical many-body problems. With the advent of ab-initio molecular dynamics, where the forces are computed on-the-fly by accurate electronic structure calculations, the scope of either method has been greatly extended. This new approach, which unifies Newton's and Schrödinger's equations, allows for complex simulations without relying on any adjustable parameter. This review is intended to outline the basic principles as well as a survey of the field. Beginning with the derivation of Born-Oppenheimer molecular dynamics, the Car-Parrinello method and the recently devised efficient and accurate Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics, which unifies best of both schemes are discussed. The predictive power of this novel second-generation Car-Parrinello approach is demonstrated by a series of applications ranging from liquid metals, to semiconductors and water. This development allows for ab-initio molecular dynamics simulations on much larger length and time scales than previously thought feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anhamonic finite temperature effects on the Raman and Infrared spectra to determine the crystal structure phase III of solid molecular hydrogen

We present theoretical calculations of the Raman and IR spectra, as well as electronic properties at zero and finite temperature to elucidate the crystal structure of phase III of solid molecular hydrogen. We find that anharmonic finite temperature are particularly important and qualitatively influences the main conclusions. While P6$_3$/m is the most likely candidate for phase III at the nuclear ground state, at finite temperature the C2/c structure appears to be more suitable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electronic signature of the instantaneous asymmetry in the first coordination shell of liquid water

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inverse simulated annealing for the determination of amorphous structures

We present a new and efficient optimization method to determine the structure of disordered systems in agreement with available experimental data. Our approach permits the application of accurate electronic structure calculations within the structure optimization. The new technique is demonstrated within density functional theory by the calculation of a model of amorphous carbon.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Recent achievements in ab initio modelling of liquid water

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Absence of Metallization in Solid Molecular Hydrogen

Being the simplest element with just one electron and proton the electronic structure of the Hydrogen atom is known exactly. However, this does not hold for the complex interplay between them in a solid and in particular not at high pressure that is known to alter the crystal as well as the electronic structure. Back in 1935 Wigner and Huntington predicted that at very high pressure solid molecular hydrogen would dissociate and form an atomic solid that is metallic. In spite of intense research efforts the experimental realization, as well as the theoretical determination of the crystal structure has remained elusive. Here we present a computational study showing that the distorted hexagonal P6$_3$/m structure is the most likely candidate for Phase III of solid hydrogen. We find that the pairing structure is very persistent and insulating over the whole pressure range, which suggests that metallization due to dissociation may precede eventual bandgap closure. Due to the fact that this not only resolve one of major disagreement between theory and experiment, but also excludes the conjectured existence of phonon-driven superconductivity in solid molecular hydrogen, our results involve a complete revision of the zero-temperature phase diagram of Phase III.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Liquid methane at extreme temperature and pressure: Implications for models of Uranus and Neptune

We present large scale electronic structure based molecular dynamics simulations of liquid methane at planetary conditions. In particular, we address the controversy of whether or not the interior of Uranus and Neptune consists of diamond. In our simulations we find no evidence for the formation of diamond, but rather sp2-bonded polymeric carbon. Furthermore, we predict that at high tem- perature hydrogen may exist in its monoatomic and metallic state. The implications of our finding for the planetary models of Uranus and Neptune are in detail discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Microscopic origins of the anomalous melting behaviour of high-pressure sodium

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio quality study of the graphite-diamond phase coexistence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

A hybrid approach to Fermi operator expansion

In a recent paper we have suggested that the finite temperature density matrix can be computed efficiently by a combination of polynomial expansion and iterative inversion techniques. We present here significant improvements over this scheme. The original complex-valued formalism is turned into a purely real one. In addition, we use Chebyshev polynomials expansion and fast summation techniques. This drastically reduces the scaling of the algorithm with the width of the Hamiltonian spectrum, which is now of the order of the cubic root of such parameter. This makes our method very competitive for applications to ab-initio simulations, when high energy resolution is required.

0 citations0 reviewsNo code linkedOpen access