Researcher profile

Jörg Behler

Jörg Behler contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2022arXiv

Insights into lithium manganese oxide-water interfaces using machine learning potentials

Unraveling the atomistic and the electronic structure of solid-liquid interfaces is the key to the design of new materials for many important applications, from heterogeneous catalysis to battery technology. Density functional theory (DFT) calculations can in principle provide a reliable description of such interfaces, but the high computational costs severely restrict the accessible time and length scales. Here, we report machine learning-driven simulations of various interfaces between water and lithium manganese oxide (Li$_x$Mn$_2$O$_4$), an important electrode material in lithium ion batteries and a catalyst for the oxygen evolution reaction. We employ a high-dimensional neural network potential (HDNNP) to compute the energies and forces several orders of magnitude faster than DFT without loss in accuracy. In addition, a high-dimensional neural network for spin prediction (HDNNS) is utilized to analyze the electronic structure of the manganese ions. Combining these methods, a series of interfaces is investigated by large-scale molecular dynamics. The simulations allow us to gain insights into a variety of properties like the dissociation of water molecules, proton transfer processes, and hydrogen bonds, as well as the geometric and electronic structure of the solid surfaces including the manganese oxidation state distribution, Jahn-Teller distortions, and electron hopping.

preprint2021arXiv

A Hessian-Based Assessment of Atomic Forces for Training Machine Learning Interatomic Potentials

In recent years, many types of machine learning potentials (MLPs) have been introduced, which are able to represent high-dimensional potential-energy surfaces (PES) with close to first-principles accuracy. Most current MLPs rely on atomic energy contributions given as a function of the local chemical environments. Frequently, in addition to total energies also atomic forces are used to construct the potentials, as they provide detailed local information about the PES. Since many systems are too large for electronic structure calculations, obtaining reliable reference forces from smaller subsystems like molecular fragments or clusters can substantially simplify the construction of the training sets. Here, we propose a method to determine structurally converged molecular fragments providing reliable atomic forces based on an analysis of the Hessian. The method, which serves as a locality test and allows to estimate the importance of long-range interactions, is illustrated for a series of molecular model systems and the metal-organic framework MOF-5 as an example for a complex organic-inorganic hybrid material.

preprint2021arXiv

High-Dimensional Neural Network Potentials for Magnetic Systems Using Spin-Dependent Atom-Centered Symmetry Functions

Machine learning potentials have emerged as a powerful tool to extend the time and length scales of first principles-quality simulations. Still, most machine learning potentials cannot distinguish different electronic spin orientations and thus are not applicable to materials in different magnetic states. Here, we propose spin-dependent atom-centered symmetry functions as a new type of descriptor taking the atomic spin degrees of freedom into account. When used as input for a high-dimensional neural network potential (HDNNP), accurate potential energy surfaces of multicomponent systems describing multiple magnetic states can be constructed. We demonstrate the performance of these magnetic HDNNPs for the case of manganese oxide, MnO. We show that the method predicts the magnetically distorted rhombohedral structure in excellent agreement with density functional theory and experiment. Its efficiency allows to determine the Néel temperature considering structural fluctuations, entropic effects, and defects. The method is general and is expected to be useful also for other types of systems like oligonuclear transition metal complexes.

preprint2020arXiv

A Bin and Hash Method for Analyzing Reference Data and Descriptors in Machine Learning Potentials

In recent years the development of machine learning (ML) potentials (MLP) has become a very active field of research. Numerous approaches have been proposed, which allow to perform extended simulations of large systems at a small fraction of the computational costs of electronic structure calculations. The key to the success of modern ML potentials is the close-to first principles quality description of the atomic interactions. This accuracy is reached by using very flexible functional forms in combination with high-level reference data from electronic structure calculations. These data sets can include up to hundreds of thousands of structures covering millions of atomic environments to ensure that all relevant features of the potential energy surface are well represented. The handling of such large data sets is nowadays becoming one of the main challenges in the construction of ML potentials. In this paper we present a method, the bin-and-hash (BAH) algorithm, to overcome this problem by enabling the efficient identification and comparison of large numbers of multidimensional vectors. Such vectors emerge in multiple contexts in the construction of ML potentials. Examples are the comparison of local atomic environments to identify and avoid unnecessary redundant information in the reference data sets that is costly in terms of both the electronic structure calculations as well as the training process, the assessment of the quality of the descriptors used as structural fingerprints in many types of ML potentials, and the detection of possibly unreliable data points. The BAH algorithm is illustrated for the example of high-dimensional neural network potentials using atom-centered symmetry functions for the geometrical description of the atomic environments, but the method is general and can be combined with any current type of ML potential.

preprint2020arXiv

A criticial view on e$_g$ occupancy as a descriptor for oxygen evolution catalytic activity in LiMn$_2$O$_4$ nanoparticles

We investigate the effect of the surface electronic structure and composition of LiMn$_2$O$_4$ nanoparticles on the electrocatalytic oxygen evolution reaction (OER). Scanning transmission electron microscopy (STEM) electron energy loss spectroscopy (EELS) studies combined with density functional theory (DFT) based simulations of the EEL spectra reveal in pristine nanoparticles a 4 nm thick surface layer with reduced average Mn oxidation state and increased Mn concentration. This is attributed to Mn$^{2+}$ partially replacing Li$^+$ at the tetrahedral sites of the spinel lattice accompanied by Mn 3d-state filling of octahedrally coordinated Mn. During electrocatalytic OER cycling, this near-surface tetrahedral Mn is leached out, thereby increasing the oxidation state of octahedrally coordinated Mn. Using rotating ring-disc electrode (RRDE) based detection of O and Mn during the OER, we show that the oxygen evolution remains constant while the Mn$^{2+}$ is removed, revealing that near-surface tetrahedrally coordinated Mn has no effect on the OER activity of LiMn$_2$O$_4$. This is surprising since the e$_g$ occupancy of Mn in octahedral sites, which is widely used as a descriptor of OER activity, changes significantly in the surface layer during cycling. The fact that e$_g$ emptying fails to correlate with OER activity here indicates that octahedral cation valence is not a fundamental measure of activity, either because the active surface state is not affected by tetrahedral Mn or because other details of the band structure or metal-oxygen bonding character, more strongly regulate the rate-limiting steps for OER.

preprint2020arXiv

A flexible and adaptive grid algorithm for global optimization utilizing basin hopping Monte Carlo

Global optimization is an active area of research in atomistic simulations, and many algorithms have been proposed to date. A prominent example is basin hopping Monte Carlo, which performs a modified Metropolis Monte Carlo search to explore the potential energy surface of the system of interest. These simulations can be very demanding due to the high-dimensional configurational search space. The effective search space can be reduced by utilizing grids for the atomic positions, but at the cost of possibly biasing the results if fixed grids are employed. In this paper, we present a flexible grid algorithm for global optimization that allows to exploit the efficiency of grids without biasing the simulation outcome. The method is general and applicable to very heterogeneous systems, such as interfaces between two materials of different crystal structure or large clusters supported at surfaces. As a benchmark case, we demonstrate its performance for the well-known global optimization problem of Lennard-Jones clusters containing up to 100 particles. In spite of the simplicity of this model potential, Lennard-Jones clusters represent a challenging test case, since the global minima for some "magic" numbers of particles exhibit geometries that are very different from those of clusters with only a slightly different size.

preprint2020arXiv

A Hybrid Density Functional Theory Benchmark Study on Lithium Manganese Oxides

The lithium manganese oxide spinel Li$_x$Mn$_2$O$_4$, with $0\leq x\leq 2$, is an important example for cathode materials in lithium ion batteries. However, an accurate description of Li$_x$Mn$_2$O$_4$ by first-principles methods like density functional theory is far from trivial due to its complex electronic structure, with a variety of energetically close electronic and magnetic states. It was found that the local density approximation as well as the generalized gradient approximation (GGA) are unable to describe Li$_x$Mn$_2$O$_4$ correctly. Here, we report an extensive benchmark for different Li$_x$Mn$_y$O$_z$ systems using the hybrid functionals PBE0 and HSE06, as well as the recently introduced local hybrid functional PBE0r. We find that all of these functionals yield energetic, structural, electronic, and magnetic properties in good agreement with experimental data. The notable benefit of the PBE0r functional, which relies on on-site Hartree-Fock exchange only, is a much reduced computational effort that is comparable to GGA functionals. Furthermore, the Hartree-Fock mixing factors in PBE0r are smaller than in PBE0, which improves the results for (lithium) manganese oxides. The investigation of Li$_x$Mn$_2$O$_4$ shows that two Mn oxidation states, +III and +IV, coexist. The Mn$^\text{III}$ ions are in the high-spin state and the corresponding MnO$_6$ octahedra are Jahn-Teller distorted. The ratio between Mn$^\text{III}$ and Mn$^\text{IV}$ and thus the electronic structure changes with the Li content while no major structural changes occur in the range from $x=0$ to $1$. This work demonstrates that the PBE0r functional provides an equally accurate and efficient description of the investigated Li$_x$Mn$_y$O$_z$ systems.

preprint2020arXiv

An experimentally validated neural-network potential energy surface for H atoms on free-standing graphene in full dimensionality

We present a first principles-quality potential energy surface (PES) describing the inter-atomic forces for hydrogen atoms interacting with free-standing graphene. The PES is a high-dimensional neural network potential that has been parameterized to 75945 data points computed with density-functional theory employing the PBE-D2 functional. Improving over a previously published PES (Jiang et al., Science, 2019, 364, 379), this neural network exhibits a realistic physisorption well and achieves a 10-fold reduction in the RMS fitting error, which is 0.6 meV/atom. We used this PES to calculate about 1.5 million classical trajectories with carefully selected initial conditions to allow for direct comparison to results of H- and D-atom scattering experiments performed at incidence translational energy of 1.9 eV and a surface temperature of 300 K. The theoretically predicted scattering angular and energy loss distributions are in good agreement with experiment, despite the fact that the experiments employed graphene grown on Pt(111). The remaining discrepancies between experiment and theory are likely due to the influence of the Pt substrate only present in the experiment.

preprint2020arXiv

Global Optimization of Copper Clusters at the ZnO(10-10) Surface Using a DFT-based Neural Network Potential and Genetic Algorithms

The determination of the most stable structures of metal clusters supported at solid surfaces by computer simulations represents a formidable challenge due to the complexity of the potential-energy surface. Here we combine a high-dimensional neural network potential, which allows to predict the energies and forces of a large number of structures with first-principles accuracy, with a global optimization scheme employing genetic algorithms. This very efficient setup is used to identify the global minima and low-energy local minima for a series of copper clusters containing between four and ten atoms adsorbed at the ZnO(10$\bar{1}$0) surface. A series of structures with common structural features resembling the Cu(111) and Cu(110) surfaces at the metal-oxide interface has been identified, and the geometries of the emerging clusters are characterized in detail. We demonstrate that the frequently employed approximation of a frozen substrate surface in global optimization can result in missing the most relevant structures.

preprint2020arXiv

Predicting Oxidation and Spin States by High-Dimensional Neural Networks: Applications to Lithium Manganese Oxide Spinels

Lithium ion batteries often contain transition metal oxides like Li$_{x}$Mn$_2$O$_4$ ($0\leq x\leq2$). Depending on the Li content different ratios of Mn$^\text{III}$ to Mn$^\text{IV}$ ions are present. In combination with electron hopping the Jahn-Teller distortions of the Mn$^\text{III}$O$_6$ octahedra can give rise to complex phenomena like structural transitions and conductance. While for small model systems oxidation and spin states can be determined using density functional theory (DFT), the investigation of dynamical phenomena by DFT is too demanding. Previously, we have shown that a high-dimensional neural network potential can extend molecular dynamics (MD) simulations of Li$_{x}$Mn$_2$O$_4$ to nanosecond time scales, but these simulations did not provide information about the electronic structure. Here we extend the use of neural networks to the prediction of atomic oxidation and spin states. The resulting high-dimensional neural network is able to predict the spins of the Mn ions with an error of only 0.03 $\hbar$. We find that the Mn e$_\text{g}$ electrons are correctly conserved and that the number of Jahn-Teller distorted Mn$^\text{III}$O$_6$ octahedra is predicted precisely for different Li loadings. A charge ordering transition is observed between 280 and 300 K, which matches resistivity measurements. Moreover, the activation energy of the electron hopping conduction above the phase transition is predicted to be 0.18 eV deviating only 0.02 eV from experiment. This work demonstrates that machine learning is able to provide an accurate representation of both, the geometric and the electronic structure dynamics of Li$_x$Mn$_2$O$_4$, on time and length scales that are not accessible by ab initio MD.

preprint2020arXiv

Properties of $α$-Brass Nanoparticles I: Neural Network Potential Energy Surface

Binary metal clusters are of high interest for applications in heterogeneous catalysis and have received much attention in recent years. To gain insights into their structure and composition at the atomic scale, computer simulations can provide valuable information if reliable interatomic potentials are available. In this paper we describe the construction of a high-dimensional neural network potential (HDNNP) intended for simulations of large brass nanoparticles with thousands of atoms, which is also applicable to bulk $α$-brass and its surfaces. The HDNNP, which is based on reference data obtained from density-functional theory calculations, is very accurate with a root mean square error of 1.7 meV/atom for total energies and 39 meV/Å for the forces of structures not included in the training set. The potential has been thoroughly validated for a wide range of energetic and structural properties of bulk $α$-brass, its surfaces as well as clusters of different size and composition demonstrating its suitability for large-scale molecular dynamics and Monte Carlo simulations with first principles accuracy.

preprint2019arXiv

A Performance and Cost Assessment of Machine Learning Interatomic Potentials

Machine learning of the quantitative relationship between local environment descriptors and the potential energy surface of a system of atoms has emerged as a new frontier in the development of interatomic potentials (IAPs). Here, we present a comprehensive evaluation of ML-IAPs based on four local environment descriptors --- Behler-Parrinello symmetry functions, smooth overlap of atomic positions (SOAP), the Spectral Neighbor Analysis Potential (SNAP) bispectrum components, and moment tensors --- using a diverse data set generated using high-throughput density functional theory (DFT) calculations. The data set comprising bcc (Li, Mo) and fcc (Cu, Ni) metals and diamond group IV semiconductors (Si, Ge) is chosen to span a range of crystal structures and bonding. All descriptors studied show excellent performance in predicting energies and forces far surpassing that of classical IAPs, as well as predicting properties such as elastic constants and phonon dispersion curves. We observe a general trade-off between accuracy and the degrees of freedom of each model, and consequently computational cost. We will discuss these trade-offs in the context of model selection for molecular dynamics and other applications.

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.