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

Albert P. Bartók

Albert P. Bartók appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.chem-phphysics.comp-phcond-mat.stat-mechcond-mat.soft

Catalog footprint

What is connected

13works
5topics
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

13 published item(s)

preprint2022arXiv

Improved uncertainty quantification for Gaussian process regression based interatomic potentials

The error estimation capability of machine learning interatomic potentials (MLIPs) based on probabilistic learning methods such as Gaussian process regression (GPR) is currently under-exploited, because of the tendancy of the predicted errors to overestimate the true error. We present approaches based on maximising either the marginal likelihood or an alternative likelihood constructed using leave-one-out cross validation to provide improved error estimates for interatomic potentials based on GPR. We benchmarked these approaches on models representing the Ar trimer, showing significant improvements in the robustness of the predicted error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Combining phonon accuracy with high transferability in Gaussian approximation potential models

Machine learning driven interatomic potentials, including Gaussian approximation potential (GAP) models, are emerging tools for atomistic simulations. Here, we address the methodological question of how one can fit GAP models that accurately predict vibrational properties in specific regions of configuration space, whilst retaining flexibility and transferability to others. We use an adaptive regularization of the GAP fit that scales with the absolute force magnitude on any given atom, thereby exploring the Bayesian interpretation of GAP regularization as an "expected error", and its impact on the prediction of physical properties for a material of interest. The approach enables excellent predictions of phonon modes (to within 0.1-0.2 THz) for structurally diverse silicon allotropes, and it can be coupled with existing fitting databases for high transferability. These findings and workflows are expected to be useful for GAP-driven materials modeling more generally.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparing molecules and solids across structural and alchemical space

Evaluating the (dis)similarity of crystalline, disordered and molecular compounds is a critical step in the development of algorithms to navigate automatically the configuration space of complex materials. For instance, a structural similarity metric is crucial for classifying structures, searching chemical space for better compounds and materials, and driving the next generation of machine-learning techniques for predicting the stability and properties of molecules and materials. In the last few years several strategies have been designed to compare atomic coordination environments. In particular, the Smooth Overlap of Atomic Positions (SOAP) has emerged as an elegant framework to obtain translation, rotation and permutation-invariant descriptors of groups of atoms, driven by the design of various classes of machine-learned inter-atomic potentials. Here we discuss how one can combine such local descriptors using a Regularized Entropy Match (REMatch) approach to describe the similarity of both whole molecular and bulk periodic structures, introducing powerful metrics that enable the navigation of alchemical and structural complexity within a unified framework. Furthermore, using this kernel and a ridge regression method we can predict atomization energies for a database of small organic molecules with a mean absolute error below 1kcal/mol, reaching an important milestone in the application of machine-learning techniques to the evaluation of molecular properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient prediction of Nucleus Independent Chemical Shifts for polycyclic aromatic hydrocarbons

Nuclear Magnetic Resonance (NMR) is one of the most powerful experimental techniques to characterize the structure of molecules and confined liquids. Nevertheless, the complexity of the systems under investigation usually requires complementary computational studies to interpret the NMR results. In this work we focus on polycyclic aromatic hydrocarbons (PAHs), an important class of organic molecules which have been commonly used as simple analogues for the spectroscopic properties of more complex systems, such as porous disordered carbons. We use Density Functional Theory (DFT) to calculate 13C chemical shifts and Nucleus Independent Chemical Shifts (NICS) for 34 PAHs. The results show a clear molecular size dependence of the two quantities, as well as the convergence of the 13C NMR shifts towards the values observed for graphene. We then present two computationally cheap models for the prediction of NICS in simple PAHs. We show that while a simple dipolar model fails to produce accurate values, a perturbative tight-binding approach can be successfully applied for the prediction of NICS in this series of molecules, including some non-planar ones containing 5- and 7-membered rings. This model, one to two orders of magnitudes faster than DFT calculations, is very promising and can be further refined in order to study more complex systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gaussian Approximation Potentials: a brief tutorial introduction

We present a swift walk-through of our recent work that uses machine learning to fit interatomic potentials based on quantum mechanical data. We describe our Gaussian Approximation Potentials (GAP) framework, discussing a variety of descriptors, how to train the model on total energies and derivatives and the simultaneous use of multiple models. We also show a small example using QUIP, the software sandbox implementation of GAP that is available for non-commercial use.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ultrasoft pseudopotentials with kinetic energy density support: implementing the modified Becke-Johnson potential

We extend the Vanderbilt ultrasoft pseudopotential scheme by adding kinetic energy density terms, in order to use meta-GGA exchange potentials, such as the Becke-Johnson or Tran-Blaha potentials, in the planewave-pseudopotential implementation of Density Functional Theory. Having implemented kinetic energy augmentation and non-linear core correction terms in the CASTEP density functional package, we evaluate the validity of our approach by comparing the calculated electronic structure of isolated atoms and semiconductor crystals to all-electron benchmark calculations. Based on our results, we provide recommendations for the practical use of the Tran-Blaha exchange in planewave-pseudopotential codes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Determining pressure-temperature phase diagrams of materials

We extend the nested sampling algorithm to simulate materials under periodic boundary and constant pressure conditions, and show how it can be used to determine the complete equilibrium phase diagram, for a given potential energy function, efficiently and in a highly automated fashion. The only inputs required are the composition and the desired pressure and temperature ranges, in particular, solid-solid phase transitions are recovered without any a priori knowledge about the structure of solid phases. We benchmark and showcase the algorithm on the periodic Lennard-Jones system, aluminium and NiTi.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accuracy and transferability of GAP models for tungsten

We introduce interatomic potentials for tungsten in the bcc crystal phase and its defects within the Gaussian Approximation Potential (GAP) framework, fitted to a database of first principles density functional theory (DFT) calculations. We investigate the performance of a sequence of models based on databases of increasing coverage in configuration space and showcase our strategy of choosing representative small unit cells to train models that predict properties only observable using thousands of atoms. The most comprehensive model is then used to calculate properties of the screw dislocation, including its structure, the Peierls barrier and the energetics of the vacancy-dislocation interaction. All software and raw data are available at www.libatoms.org.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nested sampling for materials: the case of hard spheres

The recently introduced nested sampling algorithm allows the direct and efficient calculation of the partition function of atomistic systems. We demonstrate its applicability to condensed phase systems with periodic boundary conditions by studying the three dimensional hard sphere model. Having obtained the partition function, we show how easy it is to calculate the compressibility and the free energy as functions of the packing fraction and local order, verifying that the transition to crystallinity has a very small barrier, and that the entropic contribution of jammed states to the free energy is negligible for packing fractions above the phase transition. We quantify the previously proposed schematic phase diagram and estimate the extent of the region of jammed states. We find that within our samples, the maximally random jammed configuration is surprisingly disordered.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On representing chemical environments

We review some recently published methods to represent atomic neighbourhood environments, and analyse their relative merits in terms of their faithfulness and suitability for fitting potential energy surfaces. The crucial properties that such representations (sometimes called descriptors) must have are differentiability with respect to moving the atoms, and invariance to the basic symmetries of physics: rotation, reflection, translation, and permutation of atoms of the same species. We demonstrate that certain widely used descriptors that initially look quite different are specific cases of a general approach, in which a finite set of basis functions with increasing angular wave numbers are used to expand the atomic neighbourhood density function. Using the example system of small clusters, we quantitatively show that this expansion needs to be carried to higher and higher wave numbers as the number of neighbours increases in order to obtain a faithful representation, and that variants of the descriptors converge at very different rates. We also propose an altogether new approach, called Smooth Overlap of Atomic Positions (SOAP), that sidesteps these difficulties by directly defining the similarity between any two neighbourhood environments, and show that it is still closely connected to the invariant descriptors. We test the performance of the various representations by fitting models to the potential energy surface of small silicon clusters and the bulk crystal.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient sampling of atomic configurational spaces

We describe a method to explore the configurational phase space of chemical systems. It is based on the nested sampling algorithm recently proposed by Skilling [Skilling J. (2004) In AIP Conference Proceedings, vol. 735, p. 395.; Skilling J. (2006) J of Bayesian Analysis 1:833-860.] and allows us to explore the entire potential energy surface (PES) efficiently in an unbiased way. The algorithm has two parameters which directly control the trade-off between the resolution with which the space is explored and the computational cost. We demonstrate the use of nested sampling on Lennard-Jones (LJ) clusters. Nested sampling provides a straightforward approximation for the partition function, thus evaluating expectation values of arbitrary smooth operators at arbitrary temperatures becomes a simple post-processing step. Access to absolute free energies allows us to determine the temperature-density phase diagram for LJ cluster stability. Even for relatively small clusters, the efficiency gain over parallel tempering in calculating the heat capacity is more than an order of magnitude. Furthermore, by analysing the topology of the resulting samples we are able to visualise the PES in a new and illuminating way. We identify a discretely valued order parameter with basins and supra-basins of the PES allowing a straightforward and unambiguous definition of macroscopic states of an atomistic system and the evaluation of the associated free energies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics

Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the computational resources required to solve the quantum mechanical equations limits the use of Quantum Mechanics at most to a few hundreds of atoms and only to a small fraction of the available configurational space. This thesis presents the results of my research on the development of a new interatomic potential generation scheme, which we refer to as Gaussian Approximation Potentials. In our framework, the quantum mechanical potential energy surface is interpolated between a set of predetermined values at different points in atomic configurational space by a non-linear, non-parametric regression method, the Gaussian Process. To perform the fitting, we represent the atomic environments by the bispectrum, which is invariant to permutations of the atoms in the neighbourhood and to global rotations. The result is a general scheme, that allows one to generate interatomic potentials based on arbitrary quantum mechanical data. We built a series of Gaussian Approximation Potentials using data obtained from Density Functional Theory and tested the capabilities of the method. We showed that our models reproduce the quantum mechanical potential energy surface remarkably well for the group IV semiconductors, iron and gallium nitride. Our potentials, while maintaining quantum mechanical accuracy, are several orders of magnitude faster than Quantum Mechanical methods.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Gaussian Approximation Potentials: the accuracy of quantum mechanics, without the electrons

We introduce a class of interatomic potential models that can be automatically generated from data consisting of the energies and forces experienced by atoms, derived from quantum mechanical calculations. The resulting model does not have a fixed functional form and hence is capable of modeling complex potential energy landscapes. It is systematically improvable with more data. We apply the method to bulk carbon, silicon and germanium and test it by calculating properties of the crystals at high temperatures. Using the interatomic potential to generate the long molecular dynamics trajectories required for such calculations saves orders of magnitude in computational cost.

0 citations0 reviewsNo code linkedOpen access