Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interatomic potentials”

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.

550 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Crystal Prediction and Design of Tunable Light Emission in BTB–Based Metal–Organic Frameworks

Metal-organic frameworks (MOFs) have recently been shown to exhibit unique mechanisms of luminescence based on charge transfer between structural units in the framework. These MOFs have the potential to be structural tuned for targeted emission with little or no metal participation. Here a computationally led, material design and synthesis methodology is presented that elucidates the mechanisms of light emission in interpenetrated structures comprised of metal centers (M = In, Ga, InGa, InEu) and BTB (1,3,5-Tris(4-carboxyphenyl)benzene) linkers, forming unique luminescent M-BTB MOF frameworks. Gas phase and periodic electronic structure calculations indicate that the intensity of the emission and the wavelength are overwhelmingly controlled by a combination of the number of interacting stacked linkers and their interatomic spacings, respectively. In the MOF, the ionic radii of the metal centers primarily control the expansion or shrinkage of the linker stacking distances. Experimentally, multiple M-BTB-based MOFs are synthesized and their photoluminescence was tested. Experiments validated the modeling by confirming that shifts in the crystal structure result in variations in light emission. Through this material design method, the mechanisms of tuning luminescence properties in interpenetrated M-BTB MOFs have been identified and applied to the design of MOFs with specific wavelength emission based on their structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations

Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations M n+ (M n+ = Li + , Na + , K + , Cu + , Cu 2+ , Co 2+ , Zn 2+ , Ni 2+ , Mg 2+ , Ca 2+ , Ti 2+ , Fe 2+ , Fe 3+ , Co 3+ , Cr 3+ , Ti 4+ and Cr 6+ ) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + M n+ ) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different M n+ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions ( R 2 values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.

Gong, Kai↗
DOE OSTI2025

Hund’s flat band in a frustrated spinel oxide

Electronic flat bands associated with quenched kinetic energy and heavy electron mass have attracted great interest for promoting strong electronic correlations and emergent phenomena such as high-temperature charge fractionalization and superconductivity. Intense experimental and theoretical research has been devoted to establishing the rich nontrivial metallic and heavy fermion phases intertwined with such localized electronic states. Here, in this study, we investigate the transition metal oxide spinel LiV 2 O 4 , an enigmatic heavy fermion compound lacking localized f orbital states. We use angle-resolved photoemission spectroscopy and dynamical mean-field theory to reveal a kind of correlation-induced flat band with suppressed interatomic electron hopping arising from intra-atomic Hund’s coupling. The appearance of heavy quasiparticles is ascribed to a proximate orbital-selective Mott state characterized by fluctuating local moments as evidenced by complementary magnetotransport measurements. The spectroscopic fingerprints of long-lived quasiparticles and their disappearance with increasing temperature further support the emergence of a high-temperature “bad” metal state observed in transport data. This work resolves a long-standing puzzle on the origin of heavy fermion behavior and unconventional transport in LiV 2 O 4 . Simultaneously, it opens a path to achieving flat bands through electronic interactions in d-orbital systems with geometrical frustration, potentially enabling the realization of exotic phases of matter such as the fractionalized Fermi liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-Principles Calculations of the Structural, Electronic, Optical, and Mechanical Properties of 21 Pyrophosphate Crystals

Pyrophosphate crystals have a wide array of applications in industrial and biomedical fields. However, fundamental understanding of their electronic structure, optical, and mechanical properties is still scattered and incomplete. In the present research, we report a comprehensive theoretical investigation of 21 pyrophosphates A 2 M (H 2 P 2 O 7 ) 2 •2H 2 O with either triclinic or orthorhombic crystal structure. The molecule H 2 P 2 O 7 is the dominant molecular unit, whereas A = (K, Rb, NH 4 , Tl), M = (Zn, Cu, Mg, Ni, Co, Mn), and H 2 O stand for the cation elements, transition metals, and the water molecules, respectively. The electronic structure, interatomic bonding, partial charge distribution, optical properties, and mechanical properties are investigated by first-principles calculations based on density functional theory (DFT). Most of these 21 crystals are theoretically investigated for the first time. The calculated results show a complex interplay between A, M, H 2 P 2 O 7 , and H 2 O, resulting in either metallic, half-metallic, or semi-conducting characteristics. The novel concept of total bond order density (TBOD) is used as a single quantum mechanical metric to characterize the internal cohesion of these crystals to correlate with the calculated properties, especially the mechanical properties. This work provides a large database for pyrophosphate crystals and a road map for potential applications of a wider variety of phosphates.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2025

Influence of isoelectronic mass modulation on phonon decay and transport in layered pnictogen-carbon systems

This study systematically investigates the structural, electronic, and thermal transport properties of layered pnictogen–carbon (Pn 2 C 2 ; Pn = P, As, Sb, Bi) monolayers, revealing a profound influence of isoelectronic mass modulation on phonon dynamics and thermal conductivity. Through first-principles calculations and lattice dynamics analysis, we demonstrate that increasing the atomic mass of pnictogens leads to elongated Pn-C bonds, weakened interatomic interactions, and enhanced phonon-phonon scattering, resulting in a dramatic reduction in lattice thermal conductivity (κ). Specifically, κ decreases exponentially from 65.6 W/m-K in P 2 C 2 to an ultralow 0.37 W/m-K in Bi 2 C 2 , driven by suppressed acoustic-optical phonon gaps Δ$^{Γ}_{A-0}$, increased anharmonicity, and reduced phonon lifetimes and group velocities. The transition from semiconducting to metallic behavior in heavier pnictogen-based systems further enhances phonon-electron scattering, contributing to thermal conductivity suppression. Structural analysis highlights the strengthening of out-of-plane C—C bonds and the contraction of C—Pn—C bond angles, which disrupt in-plane phonon transport. Further, these findings establish a design framework for engineering ultralow-κ materials through chemical substitution and structural tuning, offering significant potential for thermoelectric and thermal management applications. This work provides fundamental insights into phonon decay mechanisms and thermal transport in 2D materials, paving the way for advanced phononic and energy-efficient devices.

2D materials↗
DOE OSTI2021

Multiphonon diffuse scattering in solids from first principles: Application to layered crystals and two-dimensional materials

Time-resolved diffuse scattering experiments have gained increasing attention due to their potential to reveal nonequilibrium dynamics of crystal lattice vibrations with full momentum resolution. Although progress has been made in interpreting experimental data on the basis of one-phonon scattering, understanding the role of individual phonons can be sometimes hindered by multiphonon excitations. In Ref. [M. Zacharias, H. Seiler, F. Caruso, D. Zahn, F. Giustino, P. C. Kelires, and R. Ernstorfer, Phys. Rev. Lett. 127, 207401 (2021)], we have introduced a rigorous approach for the calculation of the all-phonon inelastic scattering intensity of solids from first-principles. In the present work, we describe our implementation in detail and show that multiphonon interactions are captured efficiently by exploiting translational and time-reversal symmetries of the crystal. We demonstrate its predictive power by calculating the scattering patterns of monolayer molybdenum disulfide (MoS 2 ), bulk MoS 2 , and black phosphorus (bP), and we obtain excellent agreement with our measurements of thermal electron diffuse scattering. Remarkably, our results show that multiphonon excitations dominate in bP across multiple Brillouin zones, while in MoS 2 they play a less pronounced role. We expand our analysis for each system and examine the effect of individual atomic and interatomic vibrational motion on the diffuse scattering signals. We further demonstrate the high-throughput capability of our approach by reporting all-phonon scattering maps of two-dimensional MoSe 2 , WSe 2 , WS 2 , graphene, and CdI 2 , rationalizing in each case the effect of multiphonon processes. As a side point, we show that the special displacement method reproduces the thermally distorted configuration that generates precisely the all-phonon diffuse pattern. In conclusion, the present methodology opens the way for systematic calculations of the scattering intensity in crystals and the accurate interpretation of static and time-resolved inelastic scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Structural and physical properties of 99 complex bulk chalcogenides crystals using first-principles calculations

Chalcogenide semiconductors and glasses have many applications in the civil and military fields, especially in relation to their electronic, optical and mechanical properties for energy conversion and in enviormental materials. However, they are much less systemically studied and their fundamental physical properties for a large class chalcogenide semiconductors are rather scattered and incomplete. Here, we present a detailed study using well defined first-principles calculations on the electronic structure, interatomic bonding, optical, and mechanical properties for 99 bulk chalcogenides including thirteen of these crytals which have never been calculated. Due to their unique composition and structures, these 99 bulk chalcogenides are divided into two main groups. The first group contains 54 quaternary crystals with the structure composition (A2BCQ4) (A = Ag, Cu; B = Zn, Cd, Hg, Mg, Sr, Ba; C = Si, Ge, Sn; Q = S, Se, Te), while the second group contains scattered ternary and quaternary chalcogenide crystals with a more diverse composition (AxByCzQn) (A = Ag, Cu, Ba, Cs, Li, Tl, K, Lu, Sr; B = Zn, Cd, Hg, Al, Ga, In, P, As, La, Lu, Pb, Cu, Ag; C = Si, Ge, Sn, As, Sb, Bi, Zr, Hf, Ga, In; Q = S, Se, Te; x=1, 2, 3; y=0, 1, 2, 5; z=0, 1, 2 and n=3, 4, 5, 6, 9). Moreover, the total bond order density (TBOD) is used as a single quantum mechanical metric to characterize the internal cohesion of these crystals enabling us to correlate them with the calculated properties, especially their mechanical properties. This work provides a very large database for bulk chalcogenides crucial for the future theoretical and experimental studies, opening opportunities for study the properties and potential application of a wide variety of chalcogenides.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Optimization of Thermal Conductance at Interfaces Using Machine Learning Algorithms

We report optimization of thermal transport across the interface of two different materials is critical to micro-/nanoscale electronic, photonic, and phononic devices. Although several examples of compositional intermixing at the interfaces having a positive effect on interfacial thermal conductance (ITC) have been reported, an optimum arrangement has not yet been determined because of the large number of potential atomic configurations and the significant computational cost of evaluation. On the other hand, computation-driven materials design efforts are rising in popularity and importance. Yet, the scalability and transferability of machine learning models remain as challenges in creating a complete pipeline for the simulation and analysis of large molecular systems. In this work we present a scalable Bayesian optimization framework, which leverages dynamic spawning of jobs through the Message Passing Interface (MPI) to run multiple parallel molecular dynamics simulations within a parent MPI job to optimize heat transfer at the silicon and aluminum (Si/Al) interface. We found a maximum of 50% increase in the ITC when introducing a two-layer intermixed region that consists of a higher percentage of Si. Because of the random nature of the intermixing, the magnitude of increase in the ITC varies. We observed that both homogeneity/heterogeneity of the intermixing and the intrinsic stochastic nature of molecular dynamics simulations account for the variance in ITC.

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

Top Loading Helium Cryostat Integrated with High-Pressure Cell with Fast Remote Pressure Control (Final Report for SBIR Phase 1)

A large number of physical phenomena, such as superconductivity and quantum critical phenomena, often appear only at very low temperatures below 5 K. There is an immense interest to investigate these phenomena at high pressure as a means of tuning interatomic distances, and thus the interaction parameters controlling these phenomena, in a continuous and controlled fashion. The current P-T condition for neutron scattering experiments are limited to either relatively low pressures of about 2 GPa at temperatures below 5 K, or to relatively high temperatures at pressures or tens of GPa. The purpose of this research is to develop neutron sample environment instrumentation for reaching 10-20 GPa at 1-2K with rapid and reliable online pressure and temperature control (i.e. without having to interrupt the experiment). The ultimate goal of the project is to design an integrated fast-cooling low-temperature sample environment cryogenic system compatible with state of the art neutron diamond anvil cells suitable for single-crystal neutron scattering experiments for temperatures down to 2K and pressures of several tens of GPa. The integrated system will consist of top-loading Helium flow cryostat with in-situ sample alignment mechanisms, large-volume diamond anvil cells (DAC) made from novel superalloy Pascalloy and optimized for fast cooling and heating, and a compact remote pressure control mechanism for the DAC based on a novel concept of inflatable bellows integrated with a lever-arm based force amplifier. In Phase I, we have designed, manufactured, and tested prototypes of the novel compact force-amplified pneumatic pressure control mechanism for Neutron Diamond Anvil Cell (nDAC), which allows to use pneumatic bellows system for smooth remote pressure control in the nDAC inside a top-loading cryostat with bore size of 70 mm or larger. This allows significant minimum temperature decrease in remotely controlled nDAC from 5-10 K down to 2 K. We also studied mechanical properties of a novel non-magnetic superalloy Pascalloy with different heat-treatment conditions and preliminary results indicate that up to date this is probably the strongest and most suitable material for making more compact and lightweight cryogenic nDACs for neuron scattering experiments at extreme conditions, allowing much faster cooling and heating. Potential applications: The new developments will allow to create a wide range of compact Diamond Anvil Cells for neutron diffraction which can reach several tens of GPA pressure at 2-4 K while preserving accurate remote pressure control capabilities. The new concepts can be used for developing sample environment instrumentation outside the neutron scattering field. The new development have very strong potential for expanding experimental capabilities in materials sciences and have high commercialization potential. The proposed new cryogenic high-pressure system or any of its individual components will be in high demand not only in neutron scattering facilities, but also at synchrotron beamlines and other high-pressure research facilities around the world.

36 MATERIALS SCIENCE↗