Optimized symmetry functions for machine-learning interatomic potentials of multicomponent systems

OBIWAN: An Element-Wise Scalable Feed-Forward Neural Network Potential

Estimating the potential energy of a molecular system at a quantum level of theory is a task of paramount importance in computational chemistry. The often employed density functional theory approach allows one to accomplish this task, yet most often at significant computational costs. This prompted the community to develop so-called machine learning potentials to achieve near-quantum accuracy at molecular mechanics computational cost. In this paper, we introduce OBIWAN, a feed-forward neural network that bears some relevant structural properties that also led to the definition of a new kind of general-purpose neural network layer. Its featurization process scales efficiently with newly added atomic species. This allows one to seamlessly add new atom types without requiring to change the topology of the network. Also, this allows one to train on new data sets leveraging a previously trained OBIWAN, hence converging very quickly. This avoids training from scratch and renders the approach more compliant with a green computing perspective.

Amorphous MoS2 from a machine learning inter-atomic potential

Amorphous molybdenum disulfide has shown potential as a hydrogen evolution catalyst, but the origin of its high activity is unclear, as is its atomic structure. Here, we have developed a classical inter-atomic potential using the charge equilibration neural network method, and we have employed it to generate atomic models of amorphous MoS2 by melting and quenching processes. The amorphous phase contains an abundance of molybdenum and sulfur atoms in low coordination. Besides the 6-coordinated molybdenum typical of the crystalline phases, a substantial fraction displays coordinations 4 and 5. The amorphous phase is also characterized by the appearance of direct S-S bonds. Density functional theory shows that the amorphous phase is metallic, with a considerable contribution of the 4-coordinated molybdenum to the density of states at the Fermi level. S-S bonds are related to the reduction of sulfur, with the excess electrons spread over several molybdenum atoms. Moreover, S-S bond formation is associated with a distinctive broadening of the 3s states, which could be exploited for experimental characterization of the amorphous phases. The large variety of local environments and the high density of electronic states at the Fermi level may play a positive role in increasing the electrocatalytic activity of this compound.

Physics-Inspired Structural Representations for Molecules and Materials

The first step in the construction of a regression model or a data-driven analysis, aiming to predict or elucidate the relationship between the atomic-scale structure of matter and its properties, involves transforming the Cartesian coordinates of the atoms into a suitable representation. The development of atomic-scale representations has played, and continues to play, a central role in the success of machine-learning methods for chemistry and materials science. This review summarizes the current understanding of the nature and characteristics of the most commonly used structural and chemical descriptions of atomistic structures, highlighting the deep underlying connections between different frameworks and the ideas that lead to computationally efficient and universally applicable models. It emphasizes the link between properties, structures, their physical chemistry, and their mathematical description, provides examples of recent applications to a diverse set of chemical and materials science problems, and outlines the open questions and the most promising research directions in the field.

New strontium titanate polymorphs under high pressure

We report six new dynamically stable structures of SrTiO₃ at various pressures ranging from 0 to 200 GPa. These structures were found by exploring the enthalpy surface with the Minima Hopping structure prediction method. The potential energy surface was generated by a machine learned potential, the charge equilibration via neural network technique (CENT), based on an extensive training data set of highly diverse SrTiO₃ periodic and cluster structures. All our CENT structures were validated at the level of density functional theory. For our new structures, we performed phonon calculations and NVT molecular dynamics calculations to investigate their dynamical stability. Finally, X-ray diffraction patterns were simulated to help to identify our predicted structures in experiments.

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.

Atomic partial charge predictions for furanoses by random forest regression with atom type symmetry function

Furanoses that are components for many important biomolecules have complicated conformational spaces due to the flexible ring and exo-cyclic moieties. Machine learning algorithms, which require descriptors as structural inputs, can be used to efficiently compute conformational adaptive (CA) charges to capture the electrostatic potential variations caused by the conformational changes in the molecular mechanics (MM) calculations. In the present study, we introduced atom type symmetry function (ATSF) developed based on atom centered symmetry function (ACSF) for describing conformations for furanoses, in which atoms were categorized by atom types defined by their properties or connectivity in classic molecular mechanics (MM) force field parameters to generate a suitable coordinate size. Random forest regression (RFR) models with ATSF showed improvements for predicting CA charges and dipole moments for furanoses compared to those with ACSF and atom name symmetry functions where atoms were categorized by their unique atom names. The CA charges predicted by RFR models with ATSF showed more comparable reproductions of the carbohydrate-water and carbohydrate-protein interactions computed with RESP charges individually derived from QM calculations than the ensemble-averaged atomic charge sets commonly employed in molecular mechanics force fields, suggesting that the predicted CA charges were capable of including electrostatic variations in their dynamic charge values. Improvements by ATSF showed that categorizing atoms by atom types introduced chemical structural perceptions to descriptors and produced a suitable coordinate size in ATSF to capture key structural features for furanoses. This categorizing scheme also allows ATSF to be readily adopted by other biomolecules thanks to the broad implementations of MM force fields.