Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERATOMIC FORCES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Compressive sensing lattice dynamics. I. General formalism

Ab initio calculations have been successfully used for evaluating lattice dynamical properties of solids within the (quasi)harmonic approximation (i.e., assuming noninteracting phonons with infinite lifetimes), but it remains difficult to treat anharmonicity in all but the simplest compounds. We detail a systematic information-theory-based approach to deriving ab initio anharmonic force constants: Compressive sensing lattice dynamics (CSLD). The non-negligible terms that are necessary to reproduce the first-principles calculated interatomic forces are automatically selected by minimizing the ℓ 1 norm (sum of absolute values) of the scaled force constants. Here, by using efficient sampling of the configuration space using a modest number of atomic configurations with quasirandom displacements, CSLD is well suited for deriving accurate anharmonic potentials for complex multicomponent compounds with large unit cells. We demonstrate the power and generality of CSLD by calculating the phonon lifetimes and thermal transport properties of type-I Si clathrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Nonthermal Bonding Origin of a Novel Photoexcited Lattice Instability in SnSe

Lattice dynamics measurements are often crucial tools for understanding how materials transform between different structures. Here, we report time-resolved x-ray scattering-based measurements of the nonequilibrium lattice dynamics in SnSe, a monochalcogenide reported to host a novel photoinduced lattice instability. By fitting interatomic force models to the fluence dependent excited-state dispersion, we determine the nonthermal origin of the lattice instability to be dominated by changes of interatomic interactions along a bilayer-connecting bond, rather than of an intralayer bonding network that is of primary importance to the lattice instability in thermal equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Crystalline Order Yet Glass-Like Heat Transport Driven by Hidden Local Distortions as the Structural Origin of Ultralow Thermal Conductivity in AgErTe 2

Given their rich chemical diversity and the interplay among the p-, d-, and f-orbitals of chalcogens, transition metals, and lanthanides, respectively, rare-earth transition-metal chalcogenides exhibit a wide variety of structural, magnetic, and transport phenomena. As a result, they form a particularly appealing platform for investigating structure–property relationships, emergent electronic and magnetic behaviors, and thermal transport. Here we investigate AgErTe 2 as a model system to understand phonon-glass behavior in ordered crystalline solids, which establishes the design principles for thermal barrier coatings and next-generation thermoelectrics. The local bonding asymmetry and lattice softness suppress the inherently low lattice thermal conductivity, resembling the characteristics of amorphous materials. This suppression is significantly influenced by local off-centering of Ag atoms, which breaks lattice periodicity while maintaining global crystallinity. The presence of antibonding states just below the Fermi level, arising from Ag 4d and Te 5p orbital interactions, leads to lattice softening and destabilizes ideal tetrahedral coordination, resulting in a pseudo Jahn–Teller distortion. Furthermore, the coexistence of weaker, more polarizable Ag–Te bonds and stronger Er–Te bonds creates a complex vibrational landscape enriched with low-frequency modes and enhanced phonon scattering. A pronounced disparity in interatomic force constants gives rise to highly localized, low-energy optical phonons linked to Ag rattling. These flat vibrational modes exhibit strong coupling with transverse acoustic phonons, resulting in ultrashort phonon lifetimes and mean free paths approaching interatomic distances. These features collectively enhance phonon scattering across a broad range of length and energy scales. This work offers a framework for engineering suppressed thermal conductivity in crystalline systems without the introduction of alloying elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Anomalous random correlations of force constants on the lattice dynamical properties of disordered Au-Fe alloys

Here, gold iron (Au-Fe) alloys are of immense interest due to their biocompatibility, anomalous Hall conductivity, and applications in various medical treatments. However, irrespective of the method of preparation, they often exhibit a high level of disorder with properties sensitive to the thermal or magnetic annealing temperatures. We calculate the lattice dynamical properties of Au 1–x Fe x alloys using density functional theory methods where, being multisite properties, reliable interatomic force constant (IFC) calculations in disordered alloys remain a challenge. We follow a twofold approach: an accurate IFC calculation in an environment with nominally zero chemical pair correlations to mimic the homogeneously disordered alloy and a configurational averaging for the desired phonon properties (e.g., dispersion, density of states, and entropy). We find an anomalous change in the IFC's and phonon dispersion (split bands) near x=0.19, which is attributed to the local stiffening of the Au-Au bonds when Au is in the vicinity of Fe. Other results based on mechanical and thermophysical properties reflect a similar anomaly: Phonon entropy, e.g., becomes negative below x=0.19, suggesting a tendency for chemical unmixing, reflecting the onset of a miscibility gap in the phase diagram. Our results match fairly well with reported data wherever available.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On-the-fly active learning of interpretable Bayesian force fields for atomistic rare events

Abstract Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important rare events, such as chemical reactions and diffusion. We present an adaptive Bayesian inference method for automating the training of interpretable, low-dimensional, and multi-element interatomic force fields using structures drawn on the fly from molecular dynamics simulations. Within an active learning framework, the internal uncertainty of a Gaussian process regression model is used to decide whether to accept the model prediction or to perform a first principles calculation to augment the training set of the model. The method is applied to a range of single- and multi-element systems and shown to achieve a favorable balance of accuracy and computational efficiency, while requiring a minimal amount of ab initio training data. We provide a fully open-source implementation of our method, as well as a procedure to map trained models to computationally efficient tabulated force fields.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Observation of Extraordinary Vibration Scatterings Induced by Strong Anharmonicity in Lead‐Free Halide Double Perovskites

Abstract Lead‐free halide double perovskites provide a promising solution for the long‐standing issues of lead‐containing halide perovskites, i.e., the toxicity of Pb and the low stability under ambient conditions and high‐intensity illumination. Their light‐to‐electricity or thermal‐to‐electricity conversion is strongly determined by the dynamics of the corresponding lattice vibrations. Here, the measurement of lattice dynamics is presented in a prototypical lead‐free halide double perovskite(Cs 2 NaInCl 6 ). The quantitative measurements and first‐principles calculations show that the scatterings among lattice vibrations at room temperature are at the timescale of ≈1 ps, which stems from the extraordinarily strong anharmonicity in Cs 2 NaInCl 6 . Further the degree of anharmonicity of each type of atom is quantitatively characterized in the Cs 2 NaInCl 6 single crystal, which stems from the interatomic forces, and demonstrate that this strong anharmonicity is synergistically contributed by the bond hierarchy, the tilting of the NaCl 6 and InCl 6 octahedral units, and the rattling of Cs + ions. Consequently, the crystalline Cs 2 NaInCl 6 possesses an ultralow thermal conductivity of ≈0.43 W mK −1 at room temperature, and a weak temperature dependence ofT −0.41 . These findings uncovered the underlying mechanisms behind the dynamics of lattice vibrations in double perovskites, which can largely benefit the design of optoelectronics and thermoelectrics based on halide double perovskites.

Chemistry↗
DOE OSTI2018

Thermal Transport by First-Principles Anharmonic Lattice Dynamics

Predictive methodologies for calculating lattice thermal transport and related properties have been recently developed and benchmarked successfully against measurements in a wide range of materials, from bulk to the nanoscale. This chapter presents an overview of the theoretical underpinnings of typical numerical methods for solving the Peierls-Boltzmann equation (PBE) for phonon transport when coupled with density functional theory (DFT). We discuss how these numerical methods may be employed to advance fundamental understanding of thermal transport, particularly as applied to measurements of anharmonic properties (e.g., linewidths, lattice expansion) in materials. In particular, calculations of interatomic force constants, integration methods, and construction of transport properties will be explained in detail with reference to relevant research in the literature. Furthermore, this chapter will provide an overview of applications and methods for calculations of thermal conductivities beyond perfect infinite crystals (including phonon-defect scattering) and of other measured observables, including mean free path accumulation, linewidths, and densities of states from scattering experiments.

Lindsay, Lucas↗
DOE OSTI2024

Thermal conductivity suppression in uranium-doped thorium dioxide due to phonon-spin interactions

In this work, impact of low level of uranium (U) atom substitution on thermal conductivity of thorium dioxide (ThO 2 ) is investigated. ThO 2 is an electronic insulator with a wide optical band-gap and no unpaired electrons whose thermal transport is governed by phonons. U-substitution introduces unpaired f-electrons resulting in paramagnetic behavior of U-ThO 2 at room temperature, which significantly suppresses its thermal conductivity. A single crystal of U-ThO 2 with graded composition of U is grown using a hydrothermal synthesis method, and thermal conductivity measurements are performed in regions with uniform composition of U at levels of 0%, 6%, 9% and 16%. Measured thermal conductivity profiles over 77–300 K temperature range are analyzed using an analytical expression for phonon-mediated thermal transport based on Klemens-Callaway model. Temperature dependent thermal conductivity is found to deviate significantly from the Rayleigh scattering trend expected for a simple substitutional point defect with a small perturbation to mass and interatomic forces. With the resonant scattering term, observed large suppression of thermal conductivity at low temperatures can be closely reproduced. Additionally, the extracted phonon-spin coupling constants imply a nonlinear relation of phonon-spin interaction intensity with respect to U doping percentage. Our study reveals how phonon-spin scattering contributed by unpaired f-electrons in U atoms influences thermal transport in the U-ThO 2 system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Development and assessment of hierarchical multi-reward reinforcement learning based potential for silicene with state-of-the-art models

We develop a new interatomic force field for Silicene, a 2D material with a buckled hexagonal lattice structure with high polymorphism. We introduce new parameterizations of a Tersoff model using a hierarchical multi-reward reinforcement learning (RL) methodology coupled with a continuous Monte Carlo Tree Search optimization. Our model significantly outperforms existing methods by enhancing the accuracy of predictions for the structural and thermodynamic properties of seven silicene polymorphs-including structure, energy, equation of state, elasticity, and phonon dispersion-when compared to established models. We further make a comprehensive comparison of the various models in predicting the mechanical and thermal properties of silicene. We trace the origin of the improved performance to the description of the angular dependence in the bond-order term, suggesting that modifying the angular terms in short-range models is essential to capture the structural diversity in low dimensional systems.

2D materials↗
DOE OSTI2022

Origin of ultralow phonon transport and strong anharmonicity in lead-free halide perovskites

All-inorganic lead-free halide double perovskites offer a promising avenue toward non-toxic, stable optoelectronic materials, properties that are missing in their prominent lead-containing counterparts. Their large thermopowers and high carrier mobilities also make them promising for thermoelectric applications. Here, in this work, we present a first-principles study of the lattice vibrations and thermal transport behaviors of Cs 2 SnI 6 and γ-CsSnI 3 , two prototypical compounds in this materials class. We show that conventional static zero temperature density functional theory (DFT) calculations severely underestimate the lattice thermal conductivities (κ l ) of these compounds, indicating the importance of dynamical effects. By calculating anharmonic renormalized phonon dispersions, we show that some optic phonons significantly harden with increasing temperature (T), which reduces the scattering of heat carrying phonons and enhances calculated κ l values when compared with standard zero temperature DFT. Furthermore, we demonstrate that coherence contributions to κ l , arising from wave like phonon tunneling, are important in both compounds. Overall, calculated κ l with temperature-dependent interatomic force constants, built from particle and coherence contributions, are in good agreement with available measured data, for both magnitude and temperature dependence. Large anharmonicity combined with low phonon group velocities yield ultralow values, with room temperature values of 0.26 W/m-K and 0.72 W/m-K predicted for Cs 2 SnI 6 and γ-CsSnI 3 , respectively. We further show that the lattice dynamics of these compounds are highly anharmonic, largely mediated by rotation of the SnI6 octahedra and localized modes originating from Cs rattling motion. These thermal characteristics combined with their previously computed excellent electronic properties make these perovskites promising candidates for optoelectronic and room temperature thermoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Pd 2 Se 3 Monolayer: A Promising Two-Dimensional Thermoelectric Material with Ultralow Lattice Thermal Conductivity and High Power Factor

A high power factor and low lattice thermal conductivity are two essential ingredients of highly efficient thermoelectric materials. Although monolayers of transition-metal dichalcogenides possess high power factors, high lattice thermal conductivities significantly impede their practical applications. Our first-principles calculations show that these two ingredients are well fulfilled in the recently synthesized Pd 2 Se 3 monolayer, whose crystal structure is composed of [Se 2 ] 2- dimers, Se 2- anions, and Pd 2+ cations coordinated in a square-planar manner. Our detailed analysis of third-order interatomic force constants reveals that the anharmonicity and soft phonon modes associated with covalently bonded [Se 2 ] 2- dimers lead to ultralow lattice thermal conductivities in Pd 2 Se 3 monolayers (1.5 and 2.9 W m -1 K -1 along the a- and b-axes at 300 K, respectively), which are comparable to those of high-performance bulk thermoelectric materials such as PbTe. Moreover, the “pudding-mold” type band structure, caused by Pd 2+ (d 8 ) cations coordinated in a square-planar crystal field, leads to high power factors in Pd 2 Se 3 monolayers. Consequently, both electron- and hole-doped thermoelectric materials with a considerably high zT can be achieved at moderate carrier concentrations, suggesting that Pd 2 Se 3 is a promising two-dimensional thermoelectric material. Our results suggest that hierarchical chemical bonds, that is, coexistence of different types of chemical bonds, combined with a square-planar crystal field is a promising route for designing high-efficiency thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Simulation of Inelastic Neutron Scattering Spectra Directly from Molecular Dynamics Trajectories

Inelastic neutron scattering (INS) is a widely used technique to study atomic and molecular vibrations. With the increasing complexity of materials and thus the INS spectra, being able to simulate the spectra from various atomistic models becomes an essential step and also a major bottleneck for INS data analysis. The conventional approach using density functional theory and lattice dynamics often falls short when the materials of interest are complex (e.g., defective, disordered, heterogeneous, amorphous, large-scale), for which molecular dynamics driven by an interatomic force field is a more common approach. In this paper, we demonstrate a method to directly convert molecular dynamics trajectories into simulated INS spectra, including not only fundamental but also higher order excitations. In this work, the results are compared with data collected on various representative samples from different neutron spectrometers. This development will open great opportunities by providing the key tool to perform in-depth analysis of INS data and to validate and optimize computer models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Numerical Optimization of Density Functional Tight Binding Models: Application to Molecules Containing Carbon, Hydrogen, Nitrogen, and Oxygen

New parametrizations for semiempirical density functional tight binding (DFTB) theory have been developed by the numerical optimization of adjustable parameters to minimize errors in the atomization energy and interatomic forces with respect to ab initio calculated data. Initial guesses for the radial dependences of the Slater- Koster bond integrals and overlap integrals were obtained from minimum basis density functional theory calculations. The radial dependences of the pair potentials and the bond and overlap integrals were represented by simple analytic functions. The adjustable parameters in these functions were optimized by simulated annealing and steepest descent algorithms to minimize the value of an objective function that quantifies the error between the DFTB model and ab initio calculated data. The accuracy and transferability of the resulting DFTB models for the C, H, N, and O system were assessed by comparing the predicted atomization energies and equilibrium molecular geometries of small molecules that were not included in the training data from DFTB to ab initio data. The DFTB models provide accurate predictions of the properties of hydrocarbons and more complex molecules containing C, H, N, and O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗