Engineering Papers⌕ Search

Engineering topics

Bulatov, Vasily

Publications and source records attributed to Bulatov, Vasily.

9 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Score-based denoising for atomic structure identification

We propose an effective method for removing thermal vibrations that complicate the task of analyzing complex dynamics in atomistic simulation of condensed matter. Our method iteratively subtracts thermal noises or perturbations in atomic positions using a denoising score function trained on synthetically noised but otherwise perfect crystal lattices. The resulting denoised structures clearly reveal underlying crystal order while retaining disorder associated with crystal defects. Purely geometric, agnostic to interatomic potentials, and trained without inputs from explicit simulations, our denoiser can be applied to simulation data generated from vastly different interatomic interactions. The denoiser is shown to improve existing classification methods, such as common neighbor analysis and polyhedral template matching, reaching perfect classification accuracy on a recent benchmark dataset of thermally perturbed structures up to the melting point. Demonstrated here in a wide variety of atomistic simulation contexts, the denoiser is general, robust, and readily extendable to delineate order from disorder in structurally and chemically complex materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of an interatomic potential for the W–Ta system

A physics-inspired, data-driven interatomic potential framework for multi-element system is presented. This potential is based on the generalization of the embedded atom method potential and systematically incorporates two-, three- and many-body effects. Different atomic environments are described by atom-centered Gaussian basis sets which provides sufficient flexibility to capture diverse atomic environments and allows for easy optimization of the free parameters. An interatomic potential model for the tungsten–tantalum (W–Ta) system is developed using this framework by training on data from ab initio density functional theory (DFT) calculations. The potential is thoroughly tested at various compositions by comparing bulk and defect properties from experimental and DFT data. It is shown that the potential predicts the elastic constants, defect properties, such as vacancy and interstitial formation energies, core structures of dislocations, and melting points for both pure and tantalum and tungsten with an accuracy comparable to other single element machine learning based potentials. In conclusion, the potential model is used to investigate the formation energies of ordered alloys and vacancy formation energies for chemically random W–Ta alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Parallel Dislocation Simulator

ParaDiS, or Parallel Dislocation Simulator, is a simulation tool that performs direct numerical simulation of dislocation ensembles, the carriers of plasticity, to predict the strength in crystalline materials from the fundamental physics of defect motion, evolution, and interaction. The code has been successfully deployed on high performance computing architectures and used to study the origins of strength and strain hardening for cubic crystals, the strength of micro-pillars, and irradiated materials at LLNL. The ParaDiS code has been successfully deployed on more than one hundred thousand CPU's with over ten million active degrees of freedom.

Bulatov, Vasily↗
DOE OSTI2023

One dislocation at a time

The direct observation of enhanced dislocation mobility in iron by in situ electron microscopy offers key insights and adds to the ongoing debate on the mechanisms of hydrogen embrittlement. The tiniest atom in nature, hydrogen, is pervasive and sneaky. When hydrogen atoms, either from gas or liquid sources, come into contact with a metal, they rush into the solid and squeeze through the narrow spaces between the host atoms. Being chemically active, atomic hydrogen can potentially bind to either host or impurity atoms, forming hydrides and thus changing a material’s properties. Metals, which are well known and widely used in industry for their strength and ductility, can become dangerously brittle when exposed to hydrogen. Furthermore, this phenomenon is known as hydrogen embrittlement and can lead to catastrophic failure of a load-bearing part. For example, the hydrogen produced when water molecules break apart in coolants in a nuclear reactor can cause sudden failure of the pressure vessel. In another example, although blending hydrogen in natural-gas pipelines provides a promising pathway to transition into the hydrogen economy, it could lower the fatigue resistance of the pipeline steel, making it more susceptible to crack growth due to cyclic loading induced by pressure fluctuations in the pipeline. Remarkably, perhaps as a result of a multitude of recognized failure mechanisms in which hydrogen can potentially participate, it still remains unknown exactly why hydrogen makes a metal brittle. Now, writing in Nature Materials, Longchao Huang and colleagues present unambiguous experimental observations of hydrogen-enhanced dislocation mobility in iron, offering key insights into the hydrogen embrittlement mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional

The absence of a reliable formulation of kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions which have been proposed in the literature and the appropriate coefficients are calculated via a linear regression framework. To model the nonlocal contributions, we explore two new nonlocal functionals - a functional that captures fluctuations in electronic density and a functional that incorporates gradient information. Since, the analytical functional forms of the kernels present in these nonlocal terms are not known from theory, we propose a basis function expansion to model these seemingly difficult nonlocal quantities. This allows us to easily reconstruct kernels for any system using only a few structures. The proposed method is able to learn kinetic energy densities and total kinetic energies of molecular and periodic systems, such as H 2 , LiH, LiF and a one-dimensional chain of 8 hydrogens using data from Kohn-Sham density functional theory calculations for only a few structures. For the ease of reproduction, codes used to generate the models are provided in the supporting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Coupling 2D atomistic information to 3D kink-pair enthalpy models of screw dislocations in bcc metals

The kink-pair activation enthalpy is a fundamental parameter in the theory of plasticity of body-centered cubic (bcc) metals. It controls the thermally activated motion of screw dislocation at low and intermediate temperatures. While direct atomistic calculations of kink pairs on screw dislocations have reached a high degree of accuracy, they can only be practically performed using semiempirical interatomic force fields, as electronic structure methods have not yet reached the level of efficiency needed to capture the system sizes required to model kink-pair structures. In this context, an alternative approach based on standard three-dimensional elastic models, which are efficient but lack atomic-level information, coupled to a substrate potential that represents the underlying lattice, has been widely applied over the past few years. This class of methods, known as ‘line-on-substrate’ (LOS) models, uses the substrate potential to calculate the lattice contribution to the kink-pair energies. In this work, we introduce the stress dependence of the substrate potential into LOS models to evaluate its impact on kink-pair energies. In addition, we include asymmetric dislocation core energies as an extra descriptor of the dislocation character. This asymmetry is also elevated to the continuum level by adding core energies to the general LOS formulation and used to explain potential energy differences known to exist between left and right kinks in bcc metals. Finally, more importantly, by matching the total LOS energies to previously calculated atomistic energies of kink-pair configurations, we issue a rule to establish the value of the so-called core width in nonsingular elasticity theories and reduce its arbitrariness as a mathematical construct.

36 MATERIALS SCIENCE↗