Efficient crystal structure prediction based on the symmetry principle

Crystal structure prediction (CSP) is an evolving field aimed at discerning crystal structures with minimal prior information. Despite the success of various CSP algorithms, their practical applicability remains circumscribed, particularly for large and complex systems. Here, to address this challenge, we show an evolutionary structure generator within the MAGUS (Machine Learning and Graph Theory Assisted Universal Structure Searcher) framework, inspired by the symmetry principle. This generator extracts both global and local features of explored crystal structures using group and graph theory. By integrating an on-the-fly space group miner and fragment reorganizer, augmented by symmetry-kept mutation, our approach generates higher-quality initial structures, reducing the computational costs of CSP tasks. Benchmarking tests show up to fourfold performance improvements. The method also proves valid in complex phosphorus allotrope systems. Furthermore, we apply our approach to the diamond-silicon (111)-(7 × 7) surface system, identifying up to 42 metastable structures within an 18 meV Å-2 energy range, demonstrating the efficacy of our approach in navigating challenging search spaces.

Data-efficient fine-tuning of foundational models for first-principles quality sublimation enthalpies

Calculating sublimation enthalpies of molecular crystal polymorphs is relevant to a wide range of technological applications. However, predicting these quantities at first-principles accuracy - even with the aid of machine learning potentials - is a challenge that requires sub-kJ mol-1 accuracy in the potential energy surface and finite-temperature sampling. We present an accurate and data-efficient protocol for training machine learning interatomic potentials by fine-tuning the foundational MACE-MP-0 model and showcase its capabilities on sublimation enthalpies and physical properties of ice polymorphs. Our approach requires only a few tens of training structures to achieve sub-kJ mol-1 accuracy in the sublimation enthalpies and sub-1% error in densities at finite temperature and pressure. Exploiting this data efficiency, we perform preliminary NPT simulations of hexagonal ice at the random phase approximation level and demonstrate a good agreement with experiments. Our results show promise for finite-temperature modelling of molecular crystals with the accuracy of correlated electronic structure theory methods.

Assessing the Structural Diversity of Form II Red Phosphorus via Stepwise Crystal Structure Search

Recently, the long less-known Form II red phosphorus (RP) (viz. Type II RP) was ascertained by the state-of-the-art 3-dimensional electron diffraction technique with a triclinic lattice, completely distinct from other known elemental phosphorus and leaving atomic coordinates not determined. The cell composed of ∼250 atoms might exceed the capacity of current readily available crystal structure search packages, which have been widely applied to systems with several tens of atoms. Besides, mistaking Form II RP for violet phosphorus is still surprisingly common in the studies on allotropic phosphorus due to misinterpretations on JCPDS card #00-044-0906. Herein, by reproducing annealing synthesis and cell relaxation of structures obtained in the literature, we verified two former crystal structures for Form II RP to be wrong and explained how the misinterpretations have occurred. Then, on the basis of experimental lattice data, we provided possible Form II RP models containing atomic positions by a nearly exhaustive high-throughput stepwise crystal structure search approach optimized by molecular mechanics, machine-learned force field, and density functional theory in succession. The energetic stability of Form II RP was found to rank between white phosphorus and black phosphorus, similar to the nanorod modifications.

Machine-Learned Potentials by Active Learning from Organic Crystal Structure Prediction Landscapes

A primary challenge in organic molecular crystal structure prediction (CSP) is accurately ranking the energies of potential structures. While high-level solid-state density functional theory (DFT) methods allow for mostly reliable discrimination of the low-energy structures, their high computational cost is problematic because of the need to evaluate tens to hundreds of thousands of trial crystal structures to fully explore typical crystal energy landscapes. Consequently, lower-cost but less accurate empirical force fields are often used, sometimes as the first stage of a hierarchical scheme involving multiple stages of increasingly accurate energy calculations. Machine-learned interatomic potentials (MLIPs), trained to reproduce the results of ab initio methods with computational costs close to those of force fields, can improve the efficiency of the CSP by reducing or eliminating the need for costly DFT calculations. Here, we investigate active learning methods for training MLIPs with CSP datasets. The combination of active learning with the well-developed sampling methods from CSP yields potentials in a highly automated workflow that are relevant over a wide range of the crystal packing space. To demonstrate these potentials, we illustrate efficiently reranking large, diverse crystal structure landscapes to near-DFT accuracy from force field-based CSP, improving the reliability of the final energy ranking. Furthermore, we demonstrate how these potentials can be extended to more accurately model structures far from lattice energy minima through additional on-the-fly training within Monte Carlo simulations.

Accelerating the prediction of stable materials with machine learning

Despite the rise in computing power, the large space of possible combinations of elements and crystal structure types makes large-scale high-throughput surveys of stable materials prohibitively expensive, especially for complex materials and materials subject to environmental conditions such as finite temperature. When physics-based computational methods and labor-intensive experiments are not feasible, machine learning (ML) methods can be a rapid and powerful alternative. Owing to a wealth of experimental and first-principles data as well as improved ML frameworks designed for materials modeling, ML is shown to be effective in predicting stability parameters and accelerating the discovery of new stable materials. In this Review, we summarize the most recent advancements in applying ML methodologies in predicting materials stability, focusing particularly on predictions of zero- and finite-temperature stability. We also highlight the need for more ML development in predictions of other thermodynamic knobs, such as pressure and surface/interfacial energy, which practically impact materials stability.

Insights into the Effects of Co-Regulated Factors on Li1.3Al0.3Ti1.7(PO4)3 Solid Electrolyte Preparation: Sources, Calcination Temperatures, and Sintering Temperatures

The ionic conductivity, phase components, and microstructures of LATP depend on its synthesis process. However, their relative importance and their interactions with synthesis process parameters (such as source materials, calcination temperature, and sintering temperature) remain unclear. In this work, different source materials were used to prepare LATP via the solid-state reaction method under different calcination and sintering temperatures, and an analysis via orthogonal experiments and machine learning was used to systematically study the effects of the process parameters. Sintering temperatures had the greatest effect on the total ionic conductivity of LATP pellets, followed by the sources and calcination temperatures. Sources, as the foundational factors, directly determine the composition of a major secondary phase of LATP pellets, which influences the whole process. The calcination temperature had limited impact on the ion conductivity of LATP pellets if pellets were sintered under the optimal temperature. The sintering temperature is the most important factor that influences the ion conductivity by eliminating most secondary phases and altering the microstructure of LATP, including the intergranular contact, grain size, relative densities, etc. This work offers a novel perspective to comprehend the synthesis of solid-state electrolytes beyond LATP.

Renaissance of elemental phosphorus materials: properties, synthesis, and applications in sustainable energy and environment

The polymorphism of phosphorus-based materials has garnered much research interest, and the variable chemical bonding structures give rise to a variety of micro and nanostructures. Among the different types of materials containing phosphorus, elemental phosphorus materials (EPMs) constitute the foundation for the synthesis of related compounds. EPMs are experiencing a renaissance in the post-graphene era, thanks to recent advancements in the scaling-down of black phosphorus, amorphous red phosphorus, violet phosphorus, and fibrous phosphorus and consequently, diverse classes of low-dimensional sheets, ribbons, and dots of EPMs with intriguing properties have been produced. The nanostructured EPMs featuring tunable bandgaps, moderate carrier mobility, and excellent optical absorption have shown great potential in energy conversion, energy storage, and environmental remediation. It is thus important to have a good understanding of the differences and interrelationships among diverse EPMs, their intrinsic physical and chemical properties, the synthesis of specific structures, and the selection of suitable nanostructures of EPMs for particular applications. In this comprehensive review, we aim to provide an in-depth analysis and discussion of the fundamental physicochemical properties, synthesis, and applications of EPMs in the areas of energy conversion, energy storage, and environmental remediation. Our evaluations are based on recent literature on well-established phosphorus allotropes and theoretical predictions of new EPMs. The objective of this review is to enhance our comprehension of the characteristics of EPMs, keep abreast of recent advances, and provide guidance for future research of EPMs in the fields of chemistry and materials science.

MAGUS: machine learning and graph theory assisted universal structure searcher

Crystal structure predictions based on first-principles calculations have gained great success in materials science and solid state physics. However, the remaining challenges still limit their applications in systems with a large number of atoms, especially the complexity of conformational space and the cost of local optimizations for big systems. Here, we introduce a crystal structure prediction method, MAGUS, based on the evolutionary algorithm, which addresses the above challenges with machine learning and graph theory. Techniques used in the program are summarized in detail and benchmark tests are provided. With intensive tests, we demonstrate that on-the-fly machine-learning potentials can be used to significantly reduce the number of expensive first-principles calculations, and the crystal decomposition based on graph theory can efficiently decrease the required configurations in order to find the target structures. We also summarized the representative applications of this method on several research topics, including unexpected compounds in the interior of planets and their exotic states at high pressure and high temperature (superionic, plastic, partially diffusive state, etc.); new functional materials (superhard, high-energy-density, superconducting, photoelectric materials), etc. These successful applications demonstrated that MAGUS code can help to accelerate the discovery of interesting materials and phenomena, as well as the significant value of crystal structure predictions in general.

Reversible assembly of nanoparticles: theory, strategies and computational simulations

The significant advances in synthesis and functionalization have enabled the preparation of high-quality nanoparticles that have found a plethora of successful applications. The unique physicochemical properties of nanoparticles can be manipulated through the control of size, shape, composition, and surface chemistry, but their technological application possibilities can be further expanded by exploiting the properties that emerge from their assembly. The ability to control the assembly of nanoparticles not only is required for many real technological applications, but allows the combination of the intrinsic properties of nanoparticles and opens the way to the exploitation of their complex interplay, giving access to collective properties. Significant advances and knowledge gained over the past few decades on nanoparticle assembly have made it possible to implement a growing number of strategies for reversible assembly of nanoparticles. In addition to being of interest for basic studies, such advances further broaden the range of applications and the possibility of developing innovative devices using nanoparticles. This review focuses on the reversible assembly of nanoparticles and includes the theoretical aspects related to the concept of reversibility, an up-to-date assessment of the experimental approaches applied to this field and the advanced computational schemes that offer key insights into the assembly mechanisms. We aim to provide readers with a comprehensive guide to address the challenges in assembling reversible nanoparticles and promote their applications.

Atomistic Global Optimization X: A Python package for optimization of atomistic structures

Modelling and understanding properties of materials from first principles require knowledge of the underlyingatomistic structure. This entails knowing the individual chemical identity and position of all atoms involved.Obtaining such information for macro-molecules, nano-particles, clusters, and for the surface, interface, andbulk phases of amorphous and solid materials represents a difficult high-dimensional global optimizationproblem. The rise of machine learning techniques in materials science has, however, led to many compellingdevelopments that may speed up structure searches. The complexity of such new methods has prompted aneed for an efficient way of assembling them into global optimization algorithms that can be experimentedwith. In this paper, we introduce the Atomistic Global Optimization X (AGOX) framework and code, asa customizable approach that enables efficient building and testing of global optimization algorithms. Amodular way of expressing global optimization algorithms is described and modern programming practicesare used to enable that modularity in the freely available AGOX python package. A number of examplesof global optimization approaches are implemented and analyzed. This ranges from random search andbasin-hopping to machine learning aided approaches with on-the-fly learnt surrogate energy landscapes. Themethods are show-cased on problems ranging from supported clusters over surface reconstructions to largecarbon clusters and metal-nitride clusters incorporated into graphene sheets.

Structure of an Ultrathin Oxide on Pt3 Sn(111) Solved by Machine Learning Enhanced Global Optimization

Determination of the atomic structure of solid surfaces typically depends on comparison of measured properties with simulations based on hypothesized structural models. For simple structures, the models may be guessed, but for more complex structures there is a need for reliable theory-based search algorithms. So far, such methods have been limited by the combinatorial complexity and computational expense of sufficiently accurate energy estimation for surfaces. However, the introduction of machine learning methods has the potential to change this radically. Here, we demonstrate how an evolutionary algorithm, utilizing machine learning for accelerated energy estimation and diverse population generation, can be used to solve an unknown surface structure-the (4×4) surface oxide on Pt₃ Sn(111)-based on limited experimental input. The algorithm is efficient and robust, and should be broadly applicable in surface studies, where it can replace manual, intuition based model generation.

Structure of an Ultrathin Oxide on Pt3Sn(111) Solved by Machine Learning Enhanced Global Optimization

Determination of the atomic structure of solid surfaces typically depends on comparison of measured properties with simulations based on hypothesized structural models. For simple structures, the models may be guessed, but for more complex structures there is a need for reliable theory-based search algorithms. So far, such methods have been limited by the combinatorial complexity and computational expense of sufficiently accurate energy estimation for surfaces. However, the introduction of machine learning methods has the potential to change this radically. Here, we demonstrate how an evolutionary algorithm, utilizing machine learning for accelerated energy estimation and diverse population generation, can be used to solve an unknown surface structure-the (4×4) surface oxide on Pt3Sn(111)-based on limited experimental input. The algorithm is efficient and robust, and should be broadly applicable in surface studies, where it can replace manual, intuition based model generation.

Systematic Comparison of Genetic Algorithm and Basin Hopping Approaches to the Global Optimization of Si(111) Surface Reconstructions

We present a systematic study of two widely used material structure prediction methods, the Genetic Algorithm and Basin Hopping approaches to global optimization, in a search for the 3 × 3, 5 × 5, and 7 × 7 reconstructions of the Si(111) surface. The Si(111) 7 × 7 reconstruction is the largest and most complex surface reconstruction known, and finding it is a very exacting test for global optimization methods. In this paper, we introduce a modification to previous Genetic Algorithm work on structure search for periodic systems, to allow the efficient search for surface reconstructions, and present a rigorous study of the effect of the different parameters of the algorithm. We also perform a detailed comparison with the recently improved Basin Hopping algorithm using Delocalized Internal Coordinates. Both algorithms succeeded in either resolving the 3 × 3, 5 × 5, and 7 × 7 DAS surface reconstructions or getting “sufficiently close”, i.e., identifying structures that only differ for the positions of a few atoms as well as thermally accessible structures within k_BT/unit area of the global minimum, with T = 300 K. Overall, the Genetic Algorithm is more robust with respect to parameter choice and in success rate, while the Basin Hopping method occasionally exhibits some advantages in speed of convergence. In line with previous studies, the results confirm that robustness, success, and speed of convergence of either approach are strongly influenced by how much the trial moves tend to preserve favorable bonding patterns once these appear.

Cluster Fragments in Amorphous Phosphorus and their Evolution under Pressure

Amorphous phosphorus (a-P) has long attracted interest because of its complex atomic structure, and more recently as an anode material for batteries. However, accurately describing and understanding a-P at the atomistic level remains a challenge. Here, it is shown that large-scale molecular-dynamics simulations, enabled by a machine-learning (ML)-based interatomic potential for phosphorus, can give new insights into the atomic structure of a-P and how this structure changes under pressure. The structural model so obtained contains abundant five-membered rings, as well as more complex seven- and eight-atom clusters. Changes in the simulated first sharp diffraction peak during compression and decompression indicate a hysteresis in the recovery of medium-range order. An analysis of cluster fragments, large rings, and voids suggests that moderate pressure (up to about 5 GPa) does not break the connectivity of clusters, but higher pressure does. The work provides a starting point for further computational studies of the structure and properties of a-P, and more generally it exemplifies how ML-driven modeling can accelerate the understanding of disordered functional materials.

Nanohardness from First Principles with Active Learning on Atomic Environments

We propose a methodology for the calculation of nanohardness by atomistic simulations of nanoindentation. The methodology is enabled by machine-learning interatomic potentials fitted on the fly to quantum-mechanical calculations of local fragments of the large nanoindentation simulation. We test our methodology by calculating nanohardness, as a function of load and crystallographic orientation of the surface, of diamond, AlN, SiC, BC₂N, and Si and comparing it to the calibrated values of the macro- and microhardness. The observed agreement between the computational and experimental results from the literature provides evidence that our method has sufficient predictive power to open up the possibility of designing materials with exceptional hardness directly from first principles. It will be especially valuable at the nanoscale where the experimental measurements are difficult, while empirical models fitted to macrohardness are, as a rule, inapplicable.

Gaussian Process Regression for Materials and Molecules

We provide an introduction to Gaussian process regression (GPR) machine-learning methods in computational materials science and chemistry. The focus of the present review is on the regression of atomistic properties: in particular, on the construction of interatomic potentials, or force fields, in the Gaussian Approximation Potential (GAP) framework; beyond this, we also discuss the fitting of arbitrary scalar, vectorial, and tensorial quantities. Methodological aspects of reference data generation, representation, and regression, as well as the question of how a data-driven model may be validated, are reviewed and critically discussed. A survey of applications to a variety of research questions in chemistry and materials science illustrates the rapid growth in the field. A vision is outlined for the development of the methodology in the years to come.

Perspective on integrating machine learning into computational chemistry and materials science

Machine learning (ML) methods are being used in almost every conceivable area of electronic structure theory and molecular simulation. In particular, ML has become firmly established in the construction of high-dimensional interatomic potentials. Not a day goes by without another proof of principle being published on how ML methods can represent and predict quantum mechanical properties-be they observable, such as molecular polarizabilities, or not, such as atomic charges. As ML is becoming pervasive in electronic structure theory and molecular simulation, we provide an overview of how atomistic computational modeling is being transformed by the incorporation of ML approaches. From the perspective of the practitioner in the field, we assess how common workflows to predict structure, dynamics, and spectroscopy are affected by ML. Finally, we discuss how a tighter and lasting integration of ML methods with computational chemistry and materials science can be achieved and what it will mean for research practice, software development, and postgraduate training.

Machine-learned potentials for next-generation matter simulations

The choice of simulation methods in computational materials science is driven by a fundamental trade-off: bridging large time- and length-scales with highly accurate simulations at an affordable computational cost. Venturing the investigation of complex phenomena on large scales requires fast yet accurate computational methods. We review the emerging field of machine-learned potentials, which promises to reach the accuracy of quantum mechanical computations at a substantially reduced computational cost. This Review will summarize the basic principles of the underlying machine learning methods, the data acquisition process and active learning procedures. We highlight multiple recent applications of machine-learned potentials in various fields, ranging from organic chemistry and biomolecules to inorganic crystal structure predictions and surface science. We furthermore discuss the developments required to promote a broader use of ML potentials, and the possibility of using them to help solve open questions in materials science and facilitate fully computational materials design.

Deep potential generation scheme and simulation protocol for the Li10GeP2S12-type superionic conductors

Solid-state electrolyte materials with superior lithium ionic conductivities are vital to the next-generation Li-ion batteries. Molecular dynamics could provide atomic scale information to understand the diffusion process of Li-ion in these superionic conductor materials. Here, we implement the deep potential generator to set up an efficient protocol to automatically generate interatomic potentials for Li₁₀GeP₂S₁₂-type solid-state electrolyte materials (Li₁₀GeP₂S₁₂, Li₁₀SiP₂S₁₂, and Li₁₀SnP₂S₁₂). The reliability and accuracy of the fast interatomic potentials are validated. With the potentials, we extend the simulation of the diffusion process to a wide temperature range (300 K-1000 K) and systems with large size (∼1000 atoms). Important technical aspects such as the statistical error and size effect are carefully investigated, and benchmark tests including the effect of density functional, thermal expansion, and configurational disorder are performed. The computed data that consider these factors agree well with the experimental results, and we find that the three structures show different behaviors with respect to configurational disorder. Our work paves the way for further research on computation screening of solid-state electrolyte materials.

PyXtal_FF: a python library for automated force field generation

We present PyXtal_FF—a package based on Python programming language—for developing machine learning potentials (MLPs). The aim of PyXtal_FF is to promote the application of atomistic simulations through providing several choices of atom-centered descriptors and machine learning regressions in one platform. Based on the given choice of descriptors (including the atom-centered symmetry functions, embedded atom density, SO4 bispectrum, and smooth SO3 power spectrum), PyXtal_FF can train MLPs with either generalized linear regression or neural network models, by simultaneously minimizing the errors of energy/forces/stress tensors in comparison with the data from ab-initio simulations. The trained MLP model from PyXtal_FF is interfaced with the Atomic Simulation Environment (ASE) package, which allows different types of light-weight simulations such as geometry optimization, molecular dynamics simulation, and physical properties prediction. Finally, we will illustrate the performance of PyXtal_FF by applying it to investigate several material systems, including the bulk SiO2, high entropy alloy NbMoTaW, and elemental Pt for general purposes. Full documentation of PyXtal_FF is available at https://pyxtal-ff.readthedocs.io.

Structure prediction of crystals, surfaces and nanoparticles

We review the current techniques used in the prediction of crystal structures and their surfaces and of the structures of nanoparticles. The main classes of search algorithm and energy function are summarized, and we discuss the growing role of methods based on machine learning. We illustrate the current status of the field with examples taken from metallic, inorganic and organic systems. This article is part of a discussion meeting issue 'Dynamic in situ microscopy relating structure and function'.

Can we predict materials that can be synthesised?

The discovery of materials is an important element in the development of new technologies and abilities that can help humanity tackle many challenges. Materials discovery is frustratingly slow, with the large time and resource cost often providing only small gains in property performance. Furthermore, researchers are unwilling to take large risks that they will only know the outcome of months or years later. Computation is playing an increasing role in allowing rapid screening of large numbers of materials from vast search space to identify promising candidates for laboratory synthesis and testing. However, there is a problem, in that many materials computationally predicted to have encouraging properties cannot be readily realised in the lab. This minireview looks at how we can tackle the problem of confirming that hypothetical materials are synthetically realisable, through consideration of all the stages of the materials discovery process, from obtaining the components, reacting them to a material in the correct structure, through to processing into a desired form. In an ideal world, a material prediction would come with an associated 'recipe' for the successful laboratory preparation of the material. We discuss the opportunity to thus prevent wasted effort in experimental discovery programmes, including those using automation, to accelerate the discovery of novel materials.

A general-purpose machine-learning force field for bulk and nanostructured phosphorus

Elemental phosphorus is attracting growing interest across fundamental and applied fields of research. However, atomistic simulations of phosphorus have remained an outstanding challenge. Here, we show that a universally applicable force field for phosphorus can be created by machine learning (ML) from a suitably chosen ensemble of quantum-mechanical results. Our model is fitted to density-functional theory plus many-body dispersion (DFT + MBD) data; its accuracy is demonstrated for the exfoliation of black and violet phosphorus (yielding monolayers of "phosphorene" and "hittorfene"); its transferability is shown for the transition between the molecular and network liquid phases. An application to a phosphorene nanoribbon on an experimentally relevant length scale exemplifies the power of accurate and flexible ML-driven force fields for next-generation materials modelling. The methodology promises new insights into phosphorus as well as other structurally complex, e.g., layered solids that are relevant in diverse areas of chemistry, physics, and materials science.

Combining Machine Learning Potential and Structure Prediction for Accelerated Materials Design and Discovery

The theoretical structure prediction method via quantum mechanical atomistic simulations such as density functional theory (DFT), based solely on chemical composition, has already become a routine tool to determine the structures of physical and chemical systems, e.g., solids and clusters. However, the application of DFT to more realistic simulations, to a large extent, is impeded because of the unfavorable scaling of the computational cost with respect to the system size. During recent years, the machine learning potential (MLP) method has been rapidly rising as an accurate and efficient tool for atomistic simulations. In this Perspective, we provide an introduction to the basic principles and advantages of the combination of structure prediction and MLP, as well as the challenges and opportunities associated with this promising approach.

Hierarchically Structured Allotropes of Phosphorus from Data-Driven Exploration

The discovery of materials is increasingly guided by quantum-mechanical crystal-structure prediction, but the structural complexity in bulk and nanoscale materials remains a bottleneck. Here we demonstrate how data-driven approaches can vastly accelerate the search for complex structures, combining a machine-learning (ML) model for the potential-energy surface with efficient, fragment-based searching. We use the characteristic building units observed in Hittorf's and fibrous phosphorus to seed stochastic ("random") structure searches over hundreds of thousands of runs. Our study identifies a family of hierarchically structured allotropes based on a P8 cage as principal building unit, including one-dimensional (1D) single and double helix structures, nanowires, and two-dimensional (2D) phosphorene allotropes with square-lattice and kagome topologies. These findings yield new insight into the intriguingly diverse structural chemistry of phosphorus, and they provide an example for how ML methods may, in the long run, be expected to accelerate the discovery of hierarchical nanostructures.

Machine Learning Force Fields and Coarse-Grained Variables in Molecular Dynamics: Application to Materials and Biological Systems

Machine learning encompasses tools and algorithms that 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.

Quantum Chemistry in the Age of Machine Learning

As the quantum chemistry (QC) community embraces machine learning (ML), the number of new methods and applications based on the combination of QC and ML is surging. In this Perspective, a view of the current state of affairs in this new and exciting research field is offered, challenges of using machine learning in quantum chemistry applications are described, and potential future developments are outlined. Specifically, examples of how machine learning is used to improve the accuracy and accelerate quantum chemical research are shown. Generalization and classification of existing techniques are provided to ease the navigation in the sea of literature and to guide researchers entering the field. The emphasis of this Perspective is on supervised machine learning.

Efficient Global Structure Optimization with a Machine-Learned Surrogate Model

We propose a scheme for global optimization with first-principles energy expressions of atomistic structure. While unfolding its search, the method actively learns a surrogate model of the potential energy landscape on which it performs a number of local relaxations (exploitation) and further structural searches (exploration). Assuming Gaussian processes, deploying two separate kernel widths to better capture rough features of the energy landscape while retaining a good resolution of local minima, an acquisition function is used to decide on which of the resulting structures is the more promising and should be treated at the first-principles level. The method is demonstrated to outperform by 2 orders of magnitude a well established first-principles based evolutionary algorithm in finding surface reconstructions. Finally, global optimization with first-principles energy expressions is utilized to identify initial stages of the edge oxidation and oxygen intercalation of graphene sheets on the Ir(111) surface.

A machine learning based intramolecular potential for a flexible organic molecule

Here, we employ the kernel regression machine learning technique to construct an analytical potential that reproduces the quantum mechanical potential energy surface of a small, flexible, drug-like molecule, 3-(benzyloxy)pyridin-2-amine.

Machine Learning Interatomic Potentials as Emerging Tools for Materials Science

Atomic-scale modeling and understanding of materials have made remarkable progress, but they are still fundamentally limited by the large computational cost of explicit electronic-structure methods such as density-functional theory. This Progress Report shows how machine learning (ML) is currently enabling a new degree of realism in materials modeling: by "learning" electronic-structure data, ML-based interatomic potentials give access to atomistic simulations that reach similar accuracy levels but are orders of magnitude faster. A brief introduction to the new tools is given, and then, applications to some select problems in materials science are highlighted: phase-change materials for memory devices; nanoparticle catalysts; and carbon-based electrodes for chemical sensing, supercapacitors, and batteries. It is hoped that the present work will inspire the development and wider use of ML-based interatomic potentials in diverse areas of materials research.

"Teamwork Makes the Dream Work": Tribal Competition Evolutionary Search as a Surrogate for Free-Energy-Based Structural Predictions

Crystal structure prediction has been a grand challenge in material science owing to the large configurational space that one must explore. Evolutionary (genetic) algorithms coupled with first principles calculations are commonly used in crystal structure prediction to sample the ground and metastable states of materials based on configurational energies. However, crystal structure predictions at finite temperature ( T), pressure ( P), and composition ( X) require a free-energy-based search that is often computationally expensive and tedious. Here, we introduce a new machine-learning workflow for structure prediction that is based on a concept inspired by the evolution of human tribes in primitive society. Our tribal genetic algorithm (GA) combines configurational sampling with evolutionary optimization to accurately predict entropically stabilized phases at finite ( T, P, X), at a computational cost that is an order of magnitude smaller than that required for a free-energy-based search. In a departure from standard GA techniques, the populations of individuals are divided into multiple tribes based on a bond-order fingerprint, and genetic operations are modified to ensure that cluster configurations are sampled adequately to capture entropic contributions. Team competition introduced into the evolutionary process allows winning teams (representing a better set of individuals) to expand their sizes; this translates into a more expanded search of the phase space allowing us to explore solutions near possible global minimum. Each team explores a specific section of the structural phase space and avoids bias on solutions arising from the use of individual populations in a purely energy-based search. We demonstrate the efficacy of our approach by performing the structural prediction of a representative two-dimensional two-body system as well as Lennard-Jones clusters over a range of temperatures up to its melting point. Our approach outperforms the standard GA approaches and enables structural search under "real nonambient conditions" on both bulk systems and finite-sized clusters.

Machine-Learned Fragment-Based Energies for Crystal Structure Prediction

Crystal structure prediction involves a search of a complex configurational space for local minima corresponding to stable crystal structures, which can be performed efficiently using atom-atom force fields for the assessment of intermolecular interactions. However, for challenging systems, the limitations in the accuracy of force fields prevent a reliable assessment of the relative thermodynamic stability of potential structures, while the cost of fully quantum mechanical approaches can limit applications of the methods. We present a method to rapidly improve force field lattice energies by correcting two-body interactions with a higher level of theory in a fragment-based approach and predicting these corrections with machine learning. Corrected lattice energies with commonly used density functionals and second order perturbation theory (MP2) all significantly improve the ranking of experimentally known polymorphs where the rigid molecule model is applicable. The relative lattice energies of known polymorphs are also found to systematically improve with the fragment corrections. Predicting two-body interactions with atom-centered symmetry functions in a Gaussian process is found to give highly accurate results using as little as 10-20% of the data for training, reducing the cost of the energy correction by up to an order of magnitude. The machine learning approach opens up the possibility of more widespread use of fragment-based methods in crystal structure prediction, whose increased accuracy at a low computational cost will benefit applications in areas such as polymorph screening and computer-guided materials discovery.