No Arabic abstract
We propose a systematic method to generate a complete orthonormal basis set of multipole expansion for magnetic structures in arbitrary crystal structure. The key idea is the introduction of a virtual atomic cluster of a target crystal, on which we can clearly define the magnetic configurations corresponding to symmetry-adapted multipole moments. The magnetic configurations are then mapped onto the crystal so as to preserve the magnetic point group of the multipole moments, leading to the magnetic structures classified according to the irreducible representations of crystallographic point group. We apply the present scheme to pyrhochlore and hexagonal ABO3 crystal structures, and demonstrate that the multipole expansion is useful to investigate the macroscopic responses of antiferromagnets.
It has been a challenge to accurately simulate Li-ion diffusion processes in battery materials at room temperature using {it ab initio} molecular dynamics (AIMD) due to its high computational cost. This situation has changed drastically in recent years due to the advances in machine learning-based interatomic potentials. Here we implement the Deep Potential Generator scheme to textit{automatically} generate interatomic potentials for LiGePS-type solid-state electrolyte materials. This increases our ability to simulate such materials by several orders of magnitude without sacrificing {it ab initio} accuracy. Important technical aspects like the statistical error and size effects are carefully investigated. We further establish a reliable protocol for accurate computation of Li-ion diffusion processes at experimental conditions, by investigating important technical aspects like the statistical error and size effects. Such a protocol and the automated workflow allow us to screen materials for their relevant properties with much-improved efficiency. By using the protocol and automated workflow developed here, we obtain the diffusivity data and activation energies of Li-ion diffusion that agree well with the experiment. Our work paves the way for future investigation of Li-ion diffusion mechanisms and optimization of Li-ion conductivity of solid-state electrolyte materials.
Molecular dynamics is a powerful simulation tool to explore material properties. Most of the realistic material systems are too large to be simulated with first-principles molecular dynamics. Classical molecular dynamics has lower computational cost but requires accurate force fields to achieve chemical accuracy. In this work, we develop a symmetry-adapted graph neural networks framework, named molecular dynamics graph neural networks (MDGNN), to construct force fields automatically for molecular dynamics simulations for both molecules and crystals. This architecture consistently preserves the translation, rotation and permutation invariance in the simulations. We propose a new feature engineering method including higher order contributions and show that MDGNN accurately reproduces the results of both classical and first-principles molecular dynamics. We also demonstrate that force fields constructed by the model has good transferability. Therefore, MDGNN provides an efficient and promising option for molecular dynamics simulations of large scale systems with high accuracy.
The analysis of defects and defect dynamics in crystalline materials is important for fundamental science and for a wide range of applied engineering. With increasing system size the analysis of molecular-dynamics simulation data becomes non-trivial. Here, we present a workflow for semi-automatic identification and classification of defects in crystalline structures, combining a new approach for defect description with several already existing open-source software packages. Our approach addresses the key challenges posed by the often relatively tiny volume fraction of the modified parts of the sample, thermal motion and the presence of potentially unforeseen atomic configurations (defect types) after irradiation. The local environment of any atom is converted into a rotation-invariant descriptive vector (fingerprint), which can be compared to known defect types and also yields a distance metric suited for classification. Vectors which cannot be associated to known structures indicate new types of defects. As proof-of-concept we apply our method on an iron sample to analyze the defects caused by a collision cascade induced by a 10 keV primary-knock-on-atom. The obtained results are in good agreement with reported literature values.
Deep learning has fostered many novel applications in materials informatics. However, the inverse design of inorganic crystals, $textit{i.e.}$ generating new crystal structure with targeted properties, remains a grand challenge. An important ingredient for such generative models is an invertible representation that accesses the full periodic table. This is challenging due to limited data availability and the complexity of 3D periodic crystal structures. In this paper, we present a generalized invertible representation that encodes the crystallographic information into the descriptors in both real space and reciprocal space. Combining with a generative variational autoencoder (VAE), a wide range of crystallographic structures and chemistries with desired properties can be inverse-designed. We show that our VAE model predicts novel crystal structures that do not exist in the training and test database (Materials Project) with targeted formation energies and band gaps. We validate those predicted crystals by first-principles calculations. Finally, to design solids with practical applications, we address the sparse label problem by building a semi-supervised VAE and demonstrate its successful prediction of unique thermoelectric materials
The hyperspherical harmonic basis is used to describe bound states in an $A$--body system. The approach presented here is based on the representation of the potential energy in terms of hyperspherical harmonic functions. Using this representation, the matrix elements between the basis elements are simple, and the potential energy is presented in a compact form, well suited for numerical implementation. The basis is neither symmetrized nor antisymmetrized, as required in the case of identical particles; however, after the diagonalization of the Hamiltonian matrix, the eigenvectors reflect the symmetries present in it, and the identification of the physical states is possible, as it will be shown in specific cases. We have in mind applications to atomic, molecular, and nuclear few-body systems in which symmetry breaking terms are present in the Hamiltonian; their inclusion is straightforward in the present method. As an example we solve the case of three and four particles interacting through a short-range central interaction and Coulomb potential.