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

Alexandre Tkatchenko

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

ResearcherUnclaimed source record
physics.chem-phcond-mat.mtrl-sciphysics.comp-phcond-mat.mes-hallMachine Learningquant-phBiological PhysicsBiomoleculescond-mat.othercond-mat.softcond-mat.dis-nnecon.EMphysics.atm-clusphysics.optics

Catalog footprint

What is connected

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

Stability and Vibrations of Proteins in Vacuum and Water: Bridging Quantum Accuracy and Force-Field Efficiency

Predicting biomolecular thermodynamics and spectroscopy requires accurate relative energies of metastable states and local curvatures on the potential-energy surface. We show that the general-purpose SO3LR machine-learned force field (MLFF) reproduces PBE0+MBD density-functional theory with unprecedented fidelity across bio-relevant molecules spanning sizes and complexities far beyond its training dataset. For 23 small molecules, SO3LR captures harmonic and anharmonic vibrational features, including frequencies, displacement patterns, and IR spectra. We perform detailed dynamical studies of the amino acid oF-Phe+, folding of the alanine-15 peptide, and assembly of monomeric p53 transactivation domains into tetramers, in vacuum and water. SO3LR consistently reproduces DFT-level potential-energy surfaces, vibrational densities of states, and mode eigenvectors, capturing anharmonicity, polarization, and medium-range environment-driven interactions crucial for proteins. Our results establish that MLFF-driven dynamics provide a quantum-accurate picture of metastable minima and vibrational properties, achieving DFT-level accuracy at force-field computational cost and opening new possibilities for the computational study of biomolecules.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Analyzing Atomic Interactions in Molecules as Learned by Neural Networks

While machine learning (ML) models have been able to achieve unprecedented accuracies across various prediction tasks in quantum chemistry, it is now apparent that accuracy on a test set alone is not a guarantee for robust chemical modeling such as stable molecular dynamics (MD). To go beyond accuracy, we use explainable artificial intelligence (XAI) techniques to develop a general analysis framework for atomic interactions and apply it to the SchNet and PaiNN neural network models. We compare these interactions with a set of fundamental chemical principles to understand how well the models have learned the underlying physicochemical concepts from the data. We focus on the strength of the interactions for different atomic species, how predictions for intensive and extensive quantum molecular properties are made, and analyze the decay and many-body nature of the interactions with interatomic distance. Models that deviate too far from known physical principles produce unstable MD trajectories, even when they have very high energy and force prediction accuracy. We also suggest further improvements to the ML architectures to better account for the polynomial decay of atomic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accurate Machine Learned Quantum-Mechanical Force Fields for Biomolecular Simulations

Molecular dynamics (MD) simulations allow atomistic insights into chemical and biological processes. Accurate MD simulations require computationally demanding quantum-mechanical calculations, being practically limited to short timescales and few atoms. For larger systems, efficient, but much less reliable empirical force fields are used. Recently, machine learned force fields (MLFFs) emerged as an alternative means to execute MD simulations, offering similar accuracy as ab initio methods at orders-of-magnitude speedup. Until now, MLFFs mainly capture short-range interactions in small molecules or periodic materials, due to the increased complexity of constructing models and obtaining reliable reference data for large molecules, where long-ranged many-body effects become important. This work proposes a general approach to constructing accurate MLFFs for large-scale molecular simulations (GEMS) by training on "bottom-up" and "top-down" molecular fragments of varying size, from which the relevant physicochemical interactions can be learned. GEMS is applied to study the dynamics of alanine-based peptides and the 46-residue protein crambin in aqueous solution, allowing nanosecond-scale MD simulations of >25k atoms at essentially ab initio quality. Our findings suggest that structural motifs in peptides and proteins are more flexible than previously thought, indicating that simulations at ab initio accuracy might be necessary to understand dynamic biomolecular processes such as protein (mis)folding, drug-protein binding, or allosteric regulation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Anisotropic van der Waals Dispersion Forces in Polymers: Structural Symmetry Breaking Leads to Enhanced Conformational Search

The modeling of conformations and dynamics of (bio)polymers is of primary importance for understanding physicochemical properties of soft matter. Although short-range interactions such as covalent and hydrogen bonding control the local arrangement of polymers, non-covalent interactions play a dominant role in determining the global conformations. Here we focus on how the inclusion of many-body effects in van der Waals dispersion affects the outcome of geometry optimizations and molecular dynamics simulations of model polymers. We find that delocalized force contributions are key to explore the conformational landscape, as they induce an anisotropic polarization response which efficiently guides the conformation towards globally optimized structures. This is in contrast with the commonly used pair-wise approach, where the atomic polarizabilities lack information on the global geometry. We show that such local approximation causes the conformational search to be obstructed by conformations with unphysically limited spatial symmetry, while the many-body formalism strongly reduces the roughness of the potential energy landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Four-Dimensional Scaling of Dipole Polarizability in Quantum Systems

Polarizability is a key response property of physical and chemical systems, which has an impact on intermolecular interactions, spectroscopic observables, and vacuum polarization. The calculation of polarizability for quantum systems involves an infinite sum over all excited (bound and continuum) states, concealing the physical interpretation of polarization mechanisms and complicating the derivation of efficient response models. Approximate expressions for the dipole polarizability, $α$, rely on different scaling laws $α\propto$ $R^3$, $R^4$, or $R^7$, for various definitions of the system radius $R$. Here, we consider a range of single-particle quantum systems of varying spatial dimensionality and having qualitatively different spectra, demonstrating that their polarizability follows a universal four-dimensional scaling law $α= C (4 μq^2/\hbar^2)L^4$, where $μ$ and $q$ are the (effective) particle mass and charge, $C$ is a dimensionless excitation-energy ratio, and the characteristic length $L$ is defined via the $\mathcal{L}^2$-norm of the position operator. %The applicability of this unified formula is demonstrated by accurately predicting the dipole polarizability of 36 atoms and 1641 small organic~molecules. This unified formula is also applicable to many-particle systems, as shown by} accurately predicting the dipole polarizability of 36 atoms, 1641 small organic \rrr{molecules, and Bloch electrons in periodic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Human Wellbeing and Machine Learning

There is a vast literature on the determinants of subjective wellbeing. International organisations and statistical offices are now collecting such survey data at scale. However, standard regression models explain surprisingly little of the variation in wellbeing, limiting our ability to predict it. In response, we here assess the potential of Machine Learning (ML) to help us better understand wellbeing. We analyse wellbeing data on over a million respondents from Germany, the UK, and the United States. In terms of predictive power, our ML approaches do perform better than traditional models. Although the size of the improvement is small in absolute terms, it turns out to be substantial when compared to that of key variables like health. We moreover find that drastically expanding the set of explanatory variables doubles the predictive power of both OLS and the ML approaches on unseen data. The variables identified as important by our ML algorithms - $i.e.$ material conditions, health, and meaningful social relations - are similar to those that have already been identified in the literature. In that sense, our data-driven ML results validate the findings from conventional approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecular Interactions Induced by a Static Electric Field in Quantum Mechanics and Quantum Electrodynamics

By means of quantum mechanics and quantum electrodynamics applied to coupled harmonic Drude oscillators, we study the interaction between two neutral atoms or molecules subject to a uniform static electric field. Our focus is to understand the interplay between leading contributions to field-induced electrostatics/polarization and dispersion interactions, as considered within the employed Drude model for both non-retarded and retarded regimes. For the first case, we present an exact solution for two coupled oscillators obtained by diagonalizing the corresponding quantum-mechanical Hamiltonian and demonstrate that the external field can control the strength of different intermolecular interactions and relative orientations of the molecules. In the retarded regime described by quantum electrodynamics, our analysis shows that field-induced electrostatic and polarization energies remain unchanged (in isotropic and homogeneous vacuum) compared to the nonretarded case. For interacting species modeled by quantum Drude oscillators, the developed framework based on quantum mechanics and quantum electrodynamics yields the leading contributions to molecular interactions under the combined action of external and vacuum fields.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Machine Learning Corrects Classical Force Fields: Stretching DNA Base Pairs in Explicit Solvent

In order to improve the accuracy of molecular dynamics simulations, classical force fields are supplemented with a kernel-based machine learning method trained on quantum-mechanical fragment energies. As an example application, a potential-energy surface is generalised for a small DNA duplex, taking into account explicit solvation and long-range electron exchange--correlation effects. Study of the corrected potential energy versus extension shows that leading classical DNA models have excessive stiffness with respect to stretching. This discrepancy is found to be common across multiple forcefields. The quantum correction is in qualitative agreement to the experimental thermodynamics for larger DNA double helices, providing a candidate explanation for the general and long-standing discrepancy between single molecule stretching experiments and classical calculations of DNA stretching. The new dataset of quantum calculations and the associated Kernel Modified Molecular Dynamics (KMMD) method should be of general utility in biomolecular simulations. KMMD is made available as part of the AMBER22 simulation software.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The three-center two-positron bond

Computational studies have shown that one or more positrons can stabilize two repelling atomic anions through the formation of two-center positronic bonds. In the present work, we study the energetic stability of a system containing two positrons and three hydride anions, namely 2\ce{e^+[H^{3-}_3]}. To this aim, we performed a preliminary scan of the potential energy surface of the system with both electrons and positron in a spin singlet state, with a multi-component MP2 method, that was further refined with variational and diffusion Monte Carlo calculations, and confirmed an equilibrium geometry with \ce{D_${3h}$} symmetry. The local stability of 2\ce{e^+[H^{3-}_3]} is demonstrated by analyzing the vertical detachment and adiabatic energy dissociation channels. Bonding properties of the positronic compound, such as the equilibrium interatomic distances, force constants, dissociation energies, and bonding densities are compared with those of the purely electronic \ce{H^+_3} and \ce{Li^+_3} systems. Through this analysis, we find compelling similarities between the 2\ce{e^+[H^{3-}_3]} compound and the trilithium cation. Our results strongly point out the formation of a non-electronic three-center two-positron bond, analogous to the well known three-center two-electron counterparts, which is fundamentally distinct from the two-center two-positron bond [D. Bressanini, \textit{J. Chem. Phys.} \textbf{155}, 054306 (2021)], thus extending the concept of positron bonded molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

BIGDML: Towards Exact Machine Learning Force Fields for Materials

Machine-learning force fields (MLFF) should be accurate, computationally and data efficient, and applicable to molecules, materials, and interfaces thereof. Currently, MLFFs often introduce tradeoffs that restrict their practical applicability to small subsets of chemical space or require exhaustive datasets for training. Here, we introduce the Bravais-Inspired Gradient-Domain Machine Learning (BIGDML) approach and demonstrate its ability to construct reliable force fields using a training set with just 10-200 geometries for materials including pristine and defect-containing 2D and 3D semiconductors and metals, as well as chemisorbed and physisorbed atomic and molecular adsorbates on surfaces. The BIGDML model employs the full relevant symmetry group for a given material, does not assume artificial atom types or localization of atomic interactions and exhibits high data efficiency and state-of-the-art energy accuracies (errors substantially below 1 meV per atom) for an extended set of materials. Extensive path-integral molecular dynamics carried out with BIGDML models demonstrate the counterintuitive localization of benzene--graphene dynamics induced by nuclear quantum effects and allow to rationalize the Arrhenius behavior of hydrogen diffusion coefficient in a Pd crystal for a wide range of temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Challenges for Machine Learning Force Fields in Reproducing Potential Energy Surfaces of Flexible Molecules

Dynamics of flexible molecules are often determined by an interplay between local chemical bond fluctuations and conformational changes driven by long-range electrostatics and van der Waals interactions. This interplay between interactions yields complex potential-energy surfaces (PES) with multiple minima and transition paths between them. In this work, we assess the performance of state-of-the-art Machine Learning (ML) models, namely sGDML, SchNet, GAP/SOAP, and BPNN for reproducing such PES, while using limited amounts of reference data. As a benchmark, we use the cis to trans thermal relaxation in an azobenzene molecule, where at least three different transition mechanisms should be considered. Although GAP/SOAP, SchNet, and sGDML models can globally achieve chemical accuracy of 1 kcal mol-1 with fewer than 1000 training points, predictions greatly depend on the ML method used as well as the local region of the PES being sampled. Within a given ML method, large differences can be found between predictions of close-to-equilibrium and transition regions, as well as for different transition mechanisms. We identify key challenges that the ML models face in learning long-range interactions and the intrinsic limitations of commonly used atom-based descriptors. All in all, our results suggest switching from learning the entire PES within a single model to using multiple local models with optimized descriptors, training sets, and architectures for different parts of complex PES.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Colossal Atomic Force Response in van der Waals Materials Arising From Electronic Correlations

Understanding static and dynamic phenomena in complex materials at different length scales requires reliably accounting for van der Waals (vdW) interactions, which stem from long-range electronic correlations. While the important role of many-body vdW interactions has been extensively documented when it comes to the stability of materials, much less is known about the coupling between vdW interactions and atomic forces. Here we analyze the Hessian force response matrix for a single and two vdW-coupled atomic chains to show that a many-body description of vdW interactions yields atomic force response magnitudes that exceed the expected pairwise decay by 3-5 orders of magnitude for a wide range of separations between the perturbed and the observed atom. Similar findings are confirmed for graphene and carbon nanotubes. This colossal force enhancement suggests implications for phonon spectra, free energies, interfacial adhesion, and collective dynamics in materials with many interacting atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Comprehensive Quantum Framework for Describing Retarded and Non-Retarded Molecular Interactions in External Electric Fields

We employ various quantum-mechanical approaches for studying the impact of electric fields on both nonretarded and retarded noncovalent interactions between atoms or molecules. To this end, we apply perturbative and non-perturbative methods within the frameworks of quantum mechanics (QM) as well as quantum electrodynamics (QED). In addition, to provide a transparent physical picture of the different types of resulting interactions, we employ a stochastic electrodynamic approach based on the zero-point fluctuating field. Atomic response properties are described via harmonic Drude oscillators - an efficient model system that permits an analytical solution and has been convincingly shown to yield accurate results when modeling non-retarded intermolecular interactions. The obtained intermolecular energy contributions are classified as field-induced (FI) electrostatics, FI polarization, and dispersion interactions. The interplay between these three types of interactions enables the manipulation of molecular dimer conformations by applying transversal or longitudinal electric fields along the intermolecular axis. Our framework combining four complementary theoretical approaches paves the way toward a systematic description and improved understanding of molecular interactions when molecules are subject to both external and vacuum fields.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum-Mechanical Force Balance Between Multipolar Dispersion and Pauli Repulsion in Atomic van der Waals Dimers

The structure and stability of atomic and molecular systems with van der Waals (vdW) bonding are often determined by the interplay between attractive dispersion interactions and repulsive interactions caused by electron confinement. Arising due to different mechanisms -- electron correlation for dispersion and the Pauli exclusion principle for exchange-repulsion -- these interactions do not appear to have a straightforward connection. In this paper, we use a coarse-grained approach for evaluating the exchange energy for two coupled quantum Drude oscillators and investigate the mutual compensation of the attractive and repulsive forces at the equilibrium distance within the multipole expansion of the Coulomb potential. This compensation yields a compact formula relating the vdW radius of an atom to its multipole polarizabilities, $R_{\rm vdW} = A_l^{\,}\, α_l^{{2}/{7(l+1)}}$, where $l$ is the multipole rank and $A_l$ is a conversion factor. Such a relation is compelling because it connects an electronic property of an isolated atom (atomic polarizability) with an equilibrium distance in a dimer composed of two closed-shell atoms. We assess the accuracy of the revealed formula for noble-gas, alkaline-earth, and alkali atoms and show that the $A_l$ can be assumed to be universal constants. Besides a seamless definition of vdW radii, the proposed relation can also be used for the efficient determination of atomic multipole polarizabilities solely based on the corresponding dipole polarizability and the vdW radius. Finally, our work provides a basis for the construction of efficient and minimally-empirical interatomic potentials by combining multipolar interatomic exchange and dispersion forces on an equal footing.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accurate Many-Body Repulsive Potentials for Density-Functional Tight-Binding from Deep Tensor Neural Networks

We combine density-functional tight-binding (DFTB) with deep tensor neural networks (DTNN) to maximize the strengths of both approaches in predicting structural, energetic, and vibrational molecular properties. The DTNN is used to learn a non-linear model for the localized many-body interatomic repulsive energy, which so far has been treated in an atom-pairwise manner in DFTB. Substantially improving upon standard DFTB and DTNN, the resulting DFTB-NN$_{\sf{rep}}$ model yields accurate predictions of atomization and isomerization energies, equilibrium geometries, vibrational frequencies and dihedral rotation profiles for a large variety of organic molecules compared to the hybrid DFT-PBE0 functional. Our results highlight the high potential of combining semi-empirical electronic-structure methods with physically-motivated machine learning approaches for predicting localized many-body interactions. We conclude by discussing future advancements of the DFTB-NN$_{\sf{rep}}$ approach that could enable chemically accurate electronic-structure calculations for systems with tens of thousands of atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Density-functional model for van der Waals interactions: Unifying many-body atomic approaches with nonlocal functionals

Noncovalent van der Waals (vdW) interactions are responsible for a wide range of phenomena in matter. Popular density-functional methods that treat vdW interactions use disparate physical models for these intricate forces, and as a result the applicability of these methods is often restricted to a subset of relevant molecules and materials. Aiming towards a general-purpose density functional model of vdW interactions, here we unify two complementary approaches: nonlocal vdW functionals for polarization and interatomic methods for many-body interactions. The developed nonlocal many-body dispersion method (MBD-NL) increases the accuracy and efficiency of existing vdW functionals and is shown to be broadly applicable to molecules, soft and hard materials including ionic and metallic compounds, as well as organic/inorganic interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fine-Structure Constant Connects the Polarizability of Atoms and Vacuum

We examine the recently derived quantum-mechanical relation between atomic polarizabilities and equilibrium internuclear distances in van der Waals (vdW) bonded diatomic systems [Phys. Rev. Lett. {\bf 121}, 183401 (2018)]. For homonuclear dimers, this relation is described by the compact formula $α_{\rm m}^{\rm q} = ΦR_{\rm vdW}^7$, where the constant factor in front of the vdW radius was determined empirically. Here, we derive $Φ= (4πε_0/a_0^4) \times α^{4/3}$ expressed in terms of the vacuum electric permittivity $ε_0$, the Bohr radius $a_0$, and the fine-structure constant $α$. The validity of the obtained formula is confirmed by estimating the value of the fine-structure constant from non-relativistic quantum-mechanical calculations of atomic polarizabilities and equilibrium internuclear vdW distances. The presented derivation allows to interpret the fine-structure constant as the ratio between the polarizability densities of vacuum and matter, whereas the vdW radius becomes a geometrical length scale of atoms endowed by the vacuum field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fluctuational Electrodynamics in Atomic and Macroscopic Systems: van der Waals Interactions and Radiative Heat Transfer

We present an approach to describing fluctuational electrodynamic (FED) interactions, particularly van der Waals (vdW) interactions as well as radiative heat transfer (RHT), between material bodies of vastly different length scales, allowing for going between atomistic and continuum treatments of the response of each of these bodies as desired. Any local continuum description of electromagnetic (EM) response is compatible with our approach, while atomistic descriptions in our approach are based on effective electronic and nuclear oscillator degrees of freedom, encapsulating dissipation, short-range electronic correlations, and collective nuclear vibrations (phonons). While our previous works using this approach have focused on presenting novel results, this work focuses on the derivations underlying these methods. First, we show how the distinction between "atomic" and "macroscopic" bodies is ultimately somewhat arbitrary, as formulas for vdW free energies and RHT look very similar regardless of how the distinction is drawn. Next, we demonstrate that the atomistic description of material response in our approach yields EM interaction matrix elements which are expressed in terms of analytical formulas for compact bodies or semianalytical formulas based on Ewald summation for periodic media; we use this to compute vdW interaction free energies as well as RHT powers among small biological molecules in the presence of a metallic plate as well as between parallel graphene sheets in vacuum, showing strong deviations from conventional macroscopic theories due to the confluence of geometry, phonons, and EM retardation effects. Finally, we propose formulas for efficient computation of FED interactions among material bodies in which those that are treated atomistically as well as those treated through continuum methods may have arbitrary shapes, extending previous surface-integral techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Machine learning force fields and coarse-grained variables in molecular dynamics: application to materials and biological systems

Machine learning encompasses a set of tools and algorithms which are now becoming popular in almost all scientific and technological fields. This is true for molecular dynamics as well, where machine learning offers promises of extracting valuable information from the enormous amounts of data generated by simulation of complex systems. We provide here a review of our current understanding of goals, benefits, and limitations of machine learning techniques for computational studies on atomistic systems, focusing on the construction of empirical force fields from ab-initio databases and the determination of reaction coordinates for free energy computation and enhanced sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

QM7-X: A comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules

We introduce QM7-X, a comprehensive dataset of 42 physicochemical properties for $\approx$ 4.2 M equilibrium and non-equilibrium structures of small organic molecules with up to seven non-hydrogen (C, N, O, S, Cl) atoms. To span this fundamentally important region of chemical compound space (CCS), QM7-X includes an exhaustive sampling of (meta-)stable equilibrium structures - comprised of constitutional/structural isomers and stereoisomers, e.g., enantiomers and diastereomers (including cis-/trans- and conformational isomers) - as well as 100 non-equilibrium structural variations thereof to reach a total of $\approx$ 4.2 M molecular structures. Computed at the tightly converged quantum-mechanical PBE0+MBD level of theory, QM7-X contains global (molecular) and local (atom-in-a-molecule) properties ranging from ground state quantities (such as atomization energies and dipole moments) to response quantities (such as polarizability tensors and dispersion coefficients). By providing a systematic, extensive, and tightly-converged dataset of quantum-mechanically computed physicochemical properties, we expect that QM7-X will play a critical role in the development of next-generation machine-learning based models for exploring greater swaths of CCS and performing in silico design of molecules with targeted properties.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Mechanical and Tribological Properties of Layered Materials Under High Pressure: Assessing the Importance of Many-Body Dispersion Effects

The importance of many-body dispersion effects in layered materials subjected to high external loads is evaluated. State-of-the-art many-body dispersion density functional theory calculations performed for graphite, hexagonal boron nitride, and their hetero-structures were used to fit the parameters of a classical registry-dependent interlayer potential. Using the latter, we performed extensive equilibrium molecular dynamics simulations and studied the mechanical response of homogeneous and heterogeneous bulk models under hydrostatic pressures up to 30 GPa. Comparison with experimental data demonstrates that the reliability of the many-body dispersion model extends deep into the sub-equilibrium regime. Friction simulations demonstrate the importance of many-body dispersion effects for the accurate description of the tribological properties of layered materials interfaces under high pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Impact of nuclear vibrations on van der Waals and Casimir interactions at zero and finite temperature

Van der Waals (vdW) and Casimir interactions depend crucially on material properties and geometry, especially at molecular scales, and temperature can produce noticeable relative shifts in interaction characteristics. Despite this, common treatments of these interactions ignore electromagnetic retardation, atomism, or contributions of collective mechanical vibrations (phonons) to the infrared response, which can interplay with temperature in nontrivial ways. We present a theoretical framework for computing electromagnetic interactions among molecular structures, accounting for their geometry, electronic delocalization, short-range interatomic correlations, dissipation, and phonons at atomic scales, along with long-range electromagnetic interactions among themselves or in the vicinity of continuous macroscopic bodies. We find that in carbon allotropes, particularly fullerenes, carbyne wires, and graphene sheets, phonons can couple strongly with long-range electromagnetic fields, especially at mesoscopic scales (nanometers), to create delocalized phonon polaritons that significantly modify the infrared molecular response. These polaritons especially depend on the molecular dimensionality and dissipation, and in turn affect the vdW interaction free energies of these bodies above a macroscopic gold surface, producing nonmonotonic power laws and nontrivial temperature variations at nanometer separations that are within the reach of current Casimir force experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adsorption structures and energetics of molecules on metal surfaces: Bridging experiment and theory

Adsorption geometry and stability of organic molecules on surfaces are key parameters that determine the observable properties and functions of hybrid inorganic/organic systems (HIOSs). Despite many recent advances in precise experimental characterization and improvements in first-principles electronic structure methods, reliable databases of structures and energetics for large adsorbed molecules are largely amiss. In this review, we present such a database for a range of molecules adsorbed on metal single-crystal surfaces. The systems we analyze include noble-gas atoms, conjugated aromatic molecules, carbon nanostructures, and heteroaromatic compounds adsorbed on five different metal surfaces. The overall objective is to establish a diverse benchmark dataset that enables an assessment of current and future electronic structure methods, and motivates further experimental studies that provide ever more reliable data. Specifically, the benchmark structures and energetics from experiment are here compared with the recently developed van der Waals (vdW) inclusive density-functional theory (DFT) method, DFT+vdW$^{\mathrm{surf}}$. In comparison to 23 adsorption heights and 17 adsorption energies from experiment we find a mean average deviation of 0.06 Å and 0.16 eV, respectively. This confirms the DFT+vdW$^{\mathrm{surf}}$ method as an accurate and efficient approach to treat HIOSs. A detailed discussion identifies remaining challenges to be addressed in future development of electronic structure methods, for which the here presented benchmark database may serve as an important reference.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermal and electronic fluctuations of flexible adsorbed molecules: Azobenzene on Ag(111)

We investigate the thermal and electronic collective fluctuations that contribute to the finite-temperature adsorption properties of flexible adsorbates on surfaces on the example of the molecular switch azobenzene C$_{12}$H$_{10}$N$_{2}$ on the Ag(111) surface. Using first-principles molecular dynamics simulations we obtain the free energy of adsorption that accurately accounts for entropic contributions, whereas the inclusion of many-body dispersion interactions accounts for the electronic correlations that govern the adsorbate binding. We find the adsorbate properties to be strongly entropy-driven, as can be judged by a kinetic molecular desorption prefactor of 10$^{24}$ s$^{-1}$ that largely exceeds previously reported estimates. We relate this effect to sizable fluctuations across structural and electronic observables. Comparison of our calculations to temperature-programmed desorption measurements demonstrates that finite-temperature effects play a dominant role for flexible molecules in contact with polarizable surfaces, and that recently developed first-principles methods offer an optimal tool to reveal novel collective behavior in such complex systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Many-body dispersion effects in the binding of adsorbates on metal surfaces

A correct description of electronic exchange and correlation effects for molecules in contact with extended (metal) surfaces is a challenging task for first-principles modeling. In this work we demonstrate the importance of collective van der Waals dispersion effects beyond the pairwise approximation for organic--inorganic systems on the example of atoms, molecules, and nanostructures adsorbed on metals. We use the recently developed many-body dispersion (MBD) approach in the context of density-functional theory [Phys. Rev. Lett. 108, 236402 (2012); J. Chem. Phys. 140, 18A508 (2014)] and assess its ability to correctly describe the binding of adsorbates on metal surfaces. We briefly review the MBD method and highlight its similarities to quantum-chemical approaches to electron correlation in a quasiparticle picture. In particular, we study the binding properties of xenon, 3,4,9,10-perylene-tetracarboxylic acid (PTCDA), and a graphene sheet adsorbed on the Ag(111) surface. Accounting for MBD effects we are able to describe changes in the anisotropic polarizability tensor, improve the description of adsorbate vibrations, and correctly capture the adsorbate--surface interaction screening. Comparison to other methods and experiment reveals that inclusion of MBD effects improves adsorption energies and geometries, by reducing the overbinding typically found in pairwise additive dispersion-correction approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Insight into the description of van der Waals forces for benzene adsorption on transition metal (111) surfaces

Exploring the role of van der Waals (vdW) forces on the adsorption of molecules on extended metal surfaces has become possible in recent years thanks to exciting developments in density functional theory (DFT). Among these newly developed vdW-inclusive methods, interatomic vdW approaches that account for the nonlocal screening within the bulk [V. G. Ruiz, W. Liu, E. Zojer, M. Scheffler, and A. Tkatchenko, Phys. Rev. Lett. 108, 146103 (2012)] and improved nonlocal functionals [J. Klimes, D. R. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 022201(2010)] have emerged as promising candidates to account efficiently and accurately for the lack of long-range vdW forces in most popular DFT exchange-correlation functionals. Here we have used these two approaches to compute benzene adsorption on a range of close-packed (111) surfaces upon which it either physisorbs (Cu, Ag, and Au) or chemisorbs (Rh, Pd, Ir, and Pt). We have thoroughly compared the performance between the two classes of vdW-inclusive methods and when available compared the results obtained with experimental data. By examining the computed adsorption energies, equilibrium distances, and binding curves we conclude that both methods allow for an accurate treatment of adsorption at equilibrium adsorbate-substrate distances. To this end, explicit inclusion of electrodynamic screening in the interatomic vdW scheme and optimized exchange functionals in the case of nonlocal vdW density functionals is mandatory. Nevertheless, some discrepancies are found between these two classes of methods at large adsorbate-substrate separations.

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
preprint2013arXiv

Long-range correlation energy calculated from coupled atomic response functions

An accurate determination of the electron correlation energy is essential for describing the structure, stability, and function in a wide variety of systems, ranging from gas-phase molecular assemblies to condensed matter and organic/inorganic interfaces. Even small errors in the correlation energy can have a large impact on the description of chemical and physical properties in the systems of interest. In this context, the development of efficient approaches for the accurate calculation of the long-range correlation energy (and hence dispersion) is the main challenge. In the last years a number of methods have been developed to augment density functional approximations via dispersion energy corrections, but most of these approaches ignore the intrinsic many-body nature of correlation effects, leading to inconsistent and sometimes even qualitatively incorrect predictions. Here we build upon the recent many-body dispersion (MBD) framework, which is intimately linked to the random-phase approximation for the correlation energy. We separate the correlation energy into short-range contributions that are modeled by semi-local functionals and long-range contributions that are calculated by mapping the complex all-electron problem onto a set of atomic response functions coupled in the dipole approximation. We propose an effective range-separation of the coupling between the atomic response functions that extends the already broad applicability of the MBD method to non-metallic materials with highly anisotropic responses, such as layered nanostructures. Application to a variety of high-quality benchmark datasets illustrates the accuracy and applicability of the improved MBD approach, which offers the prospect of first-principles modeling of large structurally complex systems with an accurate description of the long-range correlation energy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Machine Learning of Molecular Electronic Properties in Chemical Compound Space

The combination of modern scientific computing with electronic structure theory can lead to an unprecedented amount of data amenable to intelligent data analysis for the identification of meaningful, novel, and predictive structure-property relationships. Such relationships enable high-throughput screening for relevant properties in an exponentially growing pool of virtual compounds that are synthetically accessible. Here, we present a machine learning (ML) model, trained on a data base of \textit{ab initio} calculation results for thousands of organic molecules, that simultaneously predicts multiple electronic ground- and excited-state properties. The properties include atomization energy, polarizability, frontier orbital eigenvalues, ionization potential, electron affinity, and excitation energies. The ML model is based on a deep multi-task artificial neural network, exploiting underlying correlations between various molecular properties. The input is identical to \emph{ab initio} methods, \emph{i.e.} nuclear charges and Cartesian coordinates of all atoms. For small organic molecules the accuracy of such a "Quantum Machine" is similar, and sometimes superior, to modern quantum-chemical methods---at negligible computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Benzene Adsorbed on Metals: Concerted Effect of Covalency and van der Waals Bonding

The adsorption of aromatic molecules on metal surfaces plays a key role in condensed matter physics and functional materials. Depending on the strength of the interaction between the molecule and the surface, the binding is typically classified as either physisorption or chemisorption. Van der Waals (vdW) interactions contribute significantly to the binding in physisorbed systems, but the role of the vdW energy in chemisorbed systems remains unclear. Here we study the interaction of benzene with the (111) surface of transition metals, ranging from weak adsorption (Ag and Au) to strong adsorption (Pt, Pd, Ir, and Rh). When vdW interactions are accurately accounted for, the barrier to adsorption predicted by standard density functional theory (DFT) calculations essentially vanishes, producing a metastable precursor state on Pt and Ir surfaces. Notably, vdW forces contribute more to the binding of covalently bonded benzene than they do when benzene is physisorbed. Comparison to experimental data demonstrates that some of the recently developed methods for including vdW interactions in DFT allow quantitative treatment of both weakly and strongly adsorbed aromatic molecules on metal surfaces, extending the already excellent performance found for gas-phase molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electrodynamic Response and Stability of Molecular Crystals

We show that electrodynamic dipolar interactions, responsible for long-range fluctuations in matter, play a significant role in the stability of molecular crystals. Density functional theory calculations with van der Waals interactions determined from a semilocal "atom-in-a-molecule" model result in a large overestimation of the dielectric constants and sublimation enthalpies for polyacene crystals from naphthalene to pentacene, whereas an accurate treatment of non-local electrodynamic response leads to an agreement with the measured values for both quantities. Our findings suggest that collective response effects play a substantial role not only for optical excitations, but also for cohesive properties of non-covalently bound molecular crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Interatomic Methods for the Dispersion Energy Derived from the Adiabatic Connection Fluctuation-Dissipation Theorem

Interatomic pairwise methods are currently among the most popular and accurate ways to include dispersion energy in density functional theory (DFT) calculations. However, when applied to more than two atoms, these methods are still frequently perceived to be based on \textit{ad hoc} assumptions, rather than a rigorous derivation from quantum mechanics. Starting from the adiabatic connection fluctuation-dissipation (ACFD) theorem, an exact expression for the electronic exchange-correlation energy, we demonstrate that the pairwise interatomic dispersion energy for an arbitrary collection of isotropic polarizable dipoles emerges from the second-order expansion of the ACFD formula. Moreover, for a system of quantum harmonic oscillators coupled through a dipole--dipole potential, we prove the equivalence between the full interaction energy obtained from the Hamiltonian diagonalization and the ACFD correlation energy in the random-phase approximation. This property makes the Hamiltonian diagonalization an efficient method for the calculation of the many-body dispersion energy. In addition, we show that the switching function used to damp the dispersion interaction at short distances arises from a short-range screened Coulomb potential, whose role is to account for the spatial spread of the individual atomic dipole moments. By using the ACFD formula we gain a deeper understanding of the approximations made in the interatomic pairwise approaches, providing a powerful formalism for further development of accurate and efficient methods for the calculation of the dispersion energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Many-Body Dispersion Interactions in Molecular Crystal Polymorphism

Polymorphs in molecular crystals are often very close in energy, yet they may possess markedly different physical and chemical properties. The understanding and prediction of polymorphism is of paramount importance for a variety of applications, including pharmaceuticals, non-linear optics, and hydrogen storage. Here, we show that the non-additive many-body dispersion (MBD) energy beyond the standard pairwise approximation is crucial for the correct qualitative and quantitative description of polymorphism in molecular crystals. This is rationalized by the sensitive dependence of the MBD energy on the polymorph geometry and the ensuing dynamic electric fields inside molecular crystals. We use the glycine crystal as a fundamental and stringent benchmark case to demonstrate the accuracy of the DFT+MBD method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions

Efficient implementations of electronic structure methods are essential for first-principles modeling of molecules and solids. We here present a particularly efficient common framework for methods beyond semilocal density-functional theory, including Hartree-Fock (HF), hybrid density functionals, random-phase approximation (RPA), second-order Møller-Plesset perturbation theory (MP2), and the $GW$ method. This computational framework allows us to use compact and accurate numeric atom-centered orbitals (popular in many implementations of semilocal density-functional theory) as basis functions. The essence of our framework is to employ the "resolution of identity (RI)" technique to facilitate the treatment of both the two-electron Coulomb repulsion integrals (required in all these approaches) as well as the linear density-response function (required for RPA and $GW$). This is possible because these quantities can be expressed in terms of products of single-particle basis functions, which can in turn be expanded in a set of auxiliary basis functions (ABFs). The construction of ABFs lies at the heart of the RI technique, and here we propose a simple prescription for constructing the ABFs which can be applied regardless of whether the underlying radial functions have a specific analytical shape (e.g., Gaussian) or are numerically tabulated. We demonstrate the accuracy of our RI implementation for Gaussian and NAO basis functions, as well as the convergence behavior of our NAO basis sets for the above-mentioned methods. Benchmark results are presented for the ionization energies of 50 selected atoms and molecules from the G2 ion test set as obtained with $GW$ and MP2 self-energy methods, and the G2-I atomization energies as well as the S22 molecular interaction energies as obtained with the RPA method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Beyond the Random Phase Approximation for the Electron Correlation Energy: The Importance of Single Excitations

The random phase approximation (RPA) for the electron correlation energy, combined with the exact-exchange energy, represents the state-of-the-art exchange-correlation functional within density-functional theory (DFT). However, the standard RPA practice -- evaluating both the exact-exchange and the RPA correlation energy using local or semilocal Kohn-Sham (KS) orbitals -- leads to a systematic underbinding of molecules and solids. Here we demonstrate that this behavior is largely corrected by adding a "single excitation" (SE) contribution, so far not included in the standard RPA scheme. A similar improvement can also be achieved by replacing the non-self-consistent exact-exchange total energy by the corresponding self-consistent Hartree-Fock total energy, while retaining the RPA correlation energy evaluated using Kohn-Sham orbitals. Both schemes achieve chemical accuracy for a standard benchmark set of non-covalent intermolecular interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast and Accurate Modeling of Molecular Atomization Energies with Machine Learning

We introduce a machine learning model to predict atomization energies of a diverse set of organic molecules, based on nuclear charges and atomic positions only. The problem of solving the molecular Schrödinger equation is mapped onto a non-linear statistical regression problem of reduced complexity. Regression models are trained on and compared to atomization energies computed with hybrid density-functional theory. Cross-validation over more than seven thousand small organic molecules yields a mean absolute error of ~10 kcal/mol. Applicability is demonstrated for the prediction of molecular atomization potential energy curves.

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
preprint2010arXiv

Stacking and Registry Effects in Layered Materials: The Case of Hexagonal Boron Nitride

The interlayer sliding energy landscape of hexagonal boron nitride (h-BN) is investigated via a van der Waals corrected density functional theory approach. It is found that the main role of the van der Waals forces is to "anchor" the layers at a fixed distance, whereas the electrostatic forces dictate the optimal stacking mode and the interlayer sliding energy. A nearly free-sliding path is identified, along which bandgap modulations of ~0.6 eV are obtained. We propose a simple geometrical model that quantifies the registry matching between the layers and captures the essence of the corrugated h-BN interlayer energy landscape. The simplicity of this phenomenological model opens the way to the modeling of complex layered structures, such as carbon and boron nitride nanotubes.

0 citations0 reviewsNo code linkedOpen access