Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electromagnetic properties of indium isotopes illuminate the doubly magic character of 100 Sn

Understanding the nuclear properties in the vicinity of 100 Sn, which has been suggested to be the heaviest doubly magic nucleus with proton number Z equal to neutron number N, has been a long-standing challenge for experimental and theoretical nuclear physics. In particular, contradictory experimental evidence exists regarding the role of nuclear collectivity in this region of the nuclear chart. Here, we provide further evidence for the doubly magic character of 100 Sn by measuring the ground-state electromagnetic moments and nuclear charge radii of indium (Z = 49) isotopes as N approaches 50 from above using precision laser spectroscopy. Our results span almost the complete range between the two major closed neutron shells at N = 50 and N = 82 and reveal parabolic trends as a function of the neutron number, with a clear reduction towards these two closed neutron shells. Furthermore, a detailed comparison between our experimental results and numerical results from two complementary nuclear many-body frameworks (density functional theory and ab initio methods) exposes deficiencies in nuclear models and establishes a benchmark for future theoretical developments.

100Sn↗
DOE OSTI2025

Bridging material models across scales: An integrated approach to equation of state and molecular dynamics modeling of copper

New uncertainty-aware equation of state (EOS) and electrical conductivity (EC) models for copper have been developed. The multiphase EOS/EC models are fit to experimental solid/liquid EC isobar measurements as well as density-functional theory molecular dynamics (DFT-MD) EC calculations in both expanded and compressed regimes (0.1–16 g/ cm 3 ⁠). The liquid and solid EOS phases were fit to available experimental data along with additional DFT-MD data over the same range as the EC. Leveraging the DFT-MD data, a corresponding machine-learned interatomic potential (MLIAP) for copper was trained using genetic-algorithm optimization. The copper MLIAP was constrained by EOS shock points at high compressions. The final EOS bounded MLIAP proves to be stable over a large density range (approximately 0.1–20 g/ cm 3 ) with good agreement to an isothermal compression curve, shock Hugoniot, and liquid speed of sound measurements at high pressures (100s of GPa).

Acoustic measurements and instrumentation↗
DOE OSTI2020

Temperature-dependent lattice dynamics in iridium

The characterization of simple elemental systems is key to benchmarking first-principles modeling of electronic and vibrational behaviors of materials. A large body of literature has been built for most elemental systems in the periodic table. However, surprisingly little neutron work has been performed to probe the vibrational properties of iridium, likely due to its large neutron absorption cross section. Nonetheless, iridium is of significant scientific and technological interest due to large relativistic electron effects and electron-phonon coupling, particularly in strongly correlated iridate compounds. In this paper, we report temperature-dependent inelastic neutron scattering measurements of the vibrational properties of iridium, from which we extract key thermodynamic properties. To overcome the challenge of the large neutron absorption of iridium, we developed a simple post-processing correction procedure. The measured densities of phonon states compare well with quasiharmonic density functional theory calculations, although the obtained experimental phonon Grüneisen parameters are much larger than expected, reaching as high as γ=4.5, indicating substantial anharmonicity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ground electronic state description of thiourea coordination in homoleptic Zn 2+ , Ni 2+ and Co 2+ complexes using sulfur K -edge X-ray absorption spectroscopy

Sulfur K-edge X-ray absorption spectroscopy (XAS) was employed to experimentally characterize the coordinative bond between the thiourea (TU) or thiocarbamide ligand and transition metal (TM) ions Zn 2+ , Co 2+ and Ni 2+ in distorted tetrahedral and octahedral homoleptic coordination environments. Comparisons of XAS spectra of the free TU ligand and [Zn(TU) 4 ] 2+ , [Co(TU) 4 ] 2+ and [Ni(TU) 6 ] 2+ complexes clearly identify spectral features unique to TM 2+ –S(TU) bonding. Quantitative analysis of pre-edge intensities describes the covalency of Ni 2+ —S(TU) and Co 2+ —S(TU) bonding to be at most 21% and 9% as expressed by the S 3p contributions per TM 3d electron hole. Using relevant Ni 2+ complexes with dithiocarbamate and thioether ligands, we evaluated the empirical S 1s → 3p transition dipole integrals developed for S-donor ligands and their dependence on heteroatom substitutions. With the aid of density functional theory-based ground electronic state calculations, we found evidence for the need of using a transition dipole that is dependent on the presence of conjugated heteroatom (N) substitution in these S-donor ligands.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Low-energy moiré phonons in twisted bilayer van der Waals heterostructures

We develop a low-energy continuum model for phonons in twisted moiré bilayers, based on a configuration-space approach. In this approach, interatomic force constants are obtained from density functional theory (DFT) calculations of untwisted bilayers with various in-plane shifts. This allows for efficient computation of phonon properties for any small twist angle, while maintaining DFT-level accuracy. Based on this framework, we show how the low-energy phonon modes, including interlayer shearing and layer-breathing modes, vary with the twist angle. As the twist angle decreases, the frequencies of the low-energy modes are reordered and the atomic displacement fields corresponding to phonon eigenmodes break translational symmetry, developing periodicity on the moiré length scale. Here, we demonstrate the capabilities of our model by calculating the phonon properties of three specific structures: bilayer graphene, bilayer molybdenum disulfide (MoS 2 ), and molybdenum diselenide-tungsten diselenide (MoSe 2 -WSe 2 ).

2-dimensional systems↗
DOE OSTI2023

Efficient computational design of two-dimensional van der Waals heterostructures: Band alignment, lattice mismatch, and machine learning

Here, we develop a computational database, website applications (web-apps), and machine-learning (ML) models to accelerate the design and discovery of two-dimensional (2D) heterostructures. Using density functional theory (DFT) based lattice parameters and electronic band energies for 674 nonmetallic exfoliable 2D materials, we generate 226 779 possible bilayer heterostructures. We classify these heterostructures into type-I, -II, and -III systems according to Anderson’s rule, which is based on the relative band alignments of the noninteracting monolayers. We find that type II is the most common and type III the least common heterostructure type. We subsequently analyze the chemical trends for each heterostructure type in terms of the Periodic Table of constituent elements. The band alignment data can also be used for identifying photocatalysts and high-work-function 2D metals for contacts. We validate our results by comparing them to experimental data as well as hybrid-functional predictions. Additionally, we carry out DFT calculations of a few selected systems (MoS 2 /WSe 2 , MoS 2 /h-BN, and MoSe 2 /CrI 3 ), to compare the band-alignment description with the predictions from Anderson’s rule. We develop web-apps to enable users to virtually create combinations of 2D materials and predict their properties. Additionally, we use ML tools to predict band-alignment information for 2D materials. The web-apps, tools, and associated data will be distributed through the JARVIS-HETEROSTRUCTURE website. Our analysis, results, and the developed web-apps can be applied to the screening and design applications, such as finding alternative photocatalysts, photodetectors, and high-work-function (WF) 2D-metal contacts.

2-dimensional systems↗
DOE OSTI2020

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tackling Disorder in γ‐Ga 2 O 3

Abstract Ga 2 O 3 and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga 2 O 3 offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. γ‐Ga 2 O 3 presents a particular challenge across synthesis, characterization, and theory due to its inherent disorder and resulting complex structure–electronic‐structure relationship. Here, density functional theory is used in combination with a machine‐learning approach to screen nearly one million potential structures, thereby developing a robust atomistic model of the γ‐phase. Theoretical results are compared with surface and bulk sensitive soft and hard X‐ray photoelectron spectroscopy, X‐ray absorption spectroscopy, spectroscopic ellipsometry, and photoluminescence excitation spectroscopy experiments representative of the occupied and unoccupied states of γ‐Ga 2 O 3 . The first onset of strong absorption at room temperature is found at 5.1 eV from spectroscopic ellipsometry, which agrees well with the excitation maximum at 5.17 eV obtained by photoluminescence excitation spectroscopy, where the latter shifts to 5.33 eV at 5 K. This work presents a leap forward in the treatment of complex, disordered oxides and is a crucial step toward exploring how their electronic structure can be understood in terms of local coordination and overall structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of elevated temperature and oxygen on the capture of radioactive iodine by silver functionalized silica aerogel

Reprocessing is considered a competent strategy for spent nuclear fuel management, yet radioiodine ( 129 I) is emitted in reprocessing off-gas as a hazardous byproduct. Silver functionalized silica aerogel (Ag 0 -aerogel), a promising iodine capture material, experiences a reduction in its capacity after prolonged exposure to off-gas components at elevated temperatures, a phenomenon termed as aging. To fully understand this process, we isolated the contribution of each aging factor, exposing Ag 0 -aerogel samples to N 2 and dry air gas streams, respectively, at 150 °C for different time periods. Aged samples were loaded with I 2 to examine the capacity change and comprehensively characterized to investigate the evolution of their properties. Results show that temperature alone did not alter Ag 0 -aerogel's capacity but triggered Ag 0 nanoparticles sintering and generated organic sulfur species. The presence of O 2 reduced the capacity by ~20 %, causing (i) formation of silver sulfide (Ag 2 S) crystals and (ii) oxidation of Ag-thiolate (Ag-S-r) to Ag sulfonate (Ag-SO 3 -r). Given that Ag 2 S readily adsorbs I 2 , the formation of Ag-SO 3 -r is the major inhibitor for iodine adsorption. This hypothesis was supported by density functional theory (DFT) simulations. These findings unraveled key mechanisms of Ag 0 -aerogel aging, which are useful in the development of materials that withstand realistic spent-nuclear-fuel-reprocessing off-gas conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

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

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗
DOE OSTI2020

Predicting metal–metal interactions. I. The influence of strain on nanoparticle and metal adlayer stabilities

Strain-engineering of bimetallic nanomaterials is an important design strategy for developing new catalysts. In this paper, we introduce an approach for including strain effects into a recently introduced, density functional theory (DFT)-based alloy stability model. The model predicts adsorption site stabilities in nanoparticles and connects these site stabilities with catalytic reactivity and selectivity. Strain-based dependencies will increase the model’s accuracy for nanoparticles affected by finite-size effects. In addition to the stability of small nanoparticles, strain also influences the heat of adsorption of epitaxially grown metal-on-metal adlayers. In this respect, we successfully benchmark the strain-including alloy stability model with previous experimentally determined trends in the heats of adsorption of Au and Cu adlayers on Pt (111). For these systems, our model predicts stronger bimetallic interactions in the first monolayer than monometallic interactions in the second monolayer. We explicitly quantify the interplay between destabilizing strain effects and the energy gained by forming new metal–metal bonds. While tensile strain in the first Cu monolayer significantly destabilizes the adsorption strength, compressive strain in the first Au monolayer has a minimal impact on the heat of adsorption. Hence, this study introduces and, by comparison with previous experiments, validates an efficient DFT-based approach for strain-engineering the stability, and, in turn, the catalytic performance, of active sites in bimetallic alloys with atomic level resolution.

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

A -type antiferromagnetic order in the Zintl-phase insulator EuZn 2 P 2

Zintl phases, containing strongly covalently bonded frameworks with separate ionically bonded ions, have emerged as a critical materials family in which to couple magnetism and strong spin-orbit coupling to drive diverse topological phases of matter. Here we report the single-crystal synthesis, magnetic, thermodynamic, transport, and theoretical properties of the Zintl compound EuZn 2 P 2 that crystallizes in the anti-La 2 O 3 (CaAl 2 Si 2 ) P-3m1 structure, containing triangular layers of Eu 2+ ions. In-plane resistivity measurements reveal insulating behavior with an estimated activation energy of E g = 0.11eV. Specific heat and magnetization measurements indicate antiferromagnetic ordering at T N = 23K. Curie-Weiss analysis of in-plane and out of plane magnetic susceptibility from T = 150 to 300 K yields p eff = 8.61 for μ 0 H⊥c and p eff = 7.74 for μ 0 H//c, close to the expected values for the 4f 7 J = S = 7/2 Eu 2+ ion and indicative of weak anisotropy. Below T N , a significant anisotropy of χ ⊥ /χ // ≈ 2.3 develops, consistent with A-type magnetic order as observed in isostructural analogs and as predicted by the density functional theory calculations reported herein. The positive Weiss temperatures of θ W =19.2K for μ 0 H⊥c and θ W =41.9K for μ 0 H//c show a similar anisotropy and suggest competing ferromagnetic and antiferromagnetic interactions. Comparing Eu magnetic ordering temperatures across trigonal EuM 2 X 2 (M= divalent metal, X= pnictide) shows that EuZn 2 P 2 exhibits the highest ordering temperature, with variations in T N correlating with changes in expected dipolar interaction strengths within and between layers and independent of the magnitude of electrical conductivity. These results provide experimental validation of the crystochemical intuition that the cation Eu 2+ layers and the anionic (M 2 X 2 ) 2– framework can be treated as electronically distinct subunits, enabling further predictive materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tunable Topological Phonon for Next-generation Quantum Transduction

This project aims to understand the critical factors that determine transduction performance of topological phonons across an oxide perovskite/tungsten diselenide heterojunction. We will investigate mechanisms of their propagation and interfacial coupling using modeling and machine learning approaches. Methods include density functional theory, molecular dynamics, numerical transport simulations, and active learning for building up a training dataset for force field development.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core– electron binding energy (CEBE) analysis. We extend the development of realspace pseudopotential methods to calculating 1s, 2s, and 2p 3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn–Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

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

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗
DOE OSTI2024

Machine-Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

To cope with legacy greenhouse gas emissions and to achieve net-zero emissions by 2050, the U.S. Department of Energy (DOE) is funding efforts to develop direct air capture (DAC), a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been studied as DAC sorbent materials because of their structural and chemical diversity. Thermodynamic calculations using classical force fields are often used to screen MOFs for their performance in separations such as CO2 capture. Machine-learned force fields (MLFFs) can use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). These descriptions of forces and energies can be used to improve the accuracy of adsorption calculations. In this work, MLFF models were developed for MOFs to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption properties.

Findley, John↗
DOE OSTI2024

Multimodal fission from self-consistent calculations

When multiple fission modes coexist in a given nucleus, distinct fragment yield distributions appear. Multimodal fission has been observed in a number of fissioning nuclei spanning the nuclear chart, and this phenomenon is expected to affect the nuclear abundances synthesized during the rapid neutron-capture process (𝑟-process). In this study, we generalize the previously proposed hybrid model for fission-fragment yield distributions to predict competing fission modes and estimate the resulting yield distributions. Here, our framework allows for a comprehensive large-scale calculation of fission-fragment yields suited for 𝑟-process nuclear network studies. Nuclear density functional theory is employed to obtain the potential energy and collective inertia tensor on a multidimensional collective space defined by mass multipole moments. Fission pathways and their relative probabilities are determined using the nudged elastic band method. Based on this information, mass and charge fission yields are predicted using the recently developed hybrid model. Fission properties of fermium isotopes are calculated in the axial quadrupole-octupole collective space for three energy density functionals (EDFs). Disagreement between the EDFs appears when multiple fission modes are present. Within our framework, the UNEDF⁢1 HFB EDF agrees best with experimental data. Calculations in the axial quadrupole-octupole-hexadecapole collective space improve the agreement with the experiment for SkM*. We also discuss the sensitivity of fission predictions on the choice of EDF for several superheavy nuclei. Fission-fragment yield predictions for nuclei with multiple fission modes are sensitive to the underlying EDF. For large-scale calculations in which a minimal number of collective coordinates is considered, UNEDF⁢1 HFB provides the best description of experimental data, though the sensitivity motivates robust quantification of the uncertainties of the theoretical model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗