Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theoretical investigation of surface chemical reactivity of compositionally disordered multi-metallic alloys

Alloys are widely used in both academic and industrial research as heterogeneous catalysts for accelerating technologically important reactions ranging from hydrocarbon reforming, selective hydrogenation, to electrochemical reactions involving hydrogen and oxygen that are fundamental and indispensable to a sustainable future energy infrastructure. Most basic catalysis research has been limited to combinations of two elements that are usually ordered in some manner at the nanoscale because doing so allows observable results to be easily rationalized. It however leaves a huge material space consisting of multi metallic, compositionally disordered alloys unexplored. This research demonstrates a viable theoretical approach combining different types of first principles calculations (including Green’s function-based Korringa–Kohn–Rostoker coherent potential approximation (KKR–CPA) and planewave-based density functional theory) to enable the surface chemical reactivity of compositionally complex alloy surfaces to be investigated. We demonstrate the approach by investigating hydrogen adsorption and evolution on two types of compositionally disordered alloys: a quaternary CoCrFeNi high-entropy alloy, and binary Pt-3d base metal alloys.

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

An Overview of the Molten Salt Thermal Properties Database–Thermophysical, Version 2.1.1 (MSTDB-TP v.2.1.1)

The current status of the Molten Salt Thermodynamic Database- Thermal Physical (MSTDB-TP) is reported. Building off of MSTDB-TP and through the direction of the Roadmap for thermal property measurement of Molten Salt Reactor systems, MSTDB-TP 2.1.1 has now been released containing a total 448 salt entries (data from 140+independent studies) compared to the original commit of 62, containing thermophysical properties including density, viscosity, thermal conductivity, and heat capacity. Along with the new database release, further advancement of Saline, a C++application programming interface (API), and a graphical user interface (GUI) has facilitated increased user/developer interaction with the database. Furthermore, estimation techniques, first principles calculations (Ab-Initio) and interpolation/extrapolation methods (Redlich-Kister/Muggianu), have shown great promise in the future of thermophysical property determination for filling out compositional spaces and investigating experimentally difficult salts (hazardous/expensive). This report describes the database composition, development and advancement of database tools, and the strategy of advancing and implementing estimation data into future iterations of the database for MSTDB-TP.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Final Report

This project focused on investigating the prospects of two-dimensional (2D) materials, transition metal dichalcogenides (TMDs) and graphene, for developing nanostructures that are potentially useful for quantum science. These materials class exhibit a range of fascinating electronic and optical properties and, due to their sub-nanometer thickness, are particularly well suited for confining electronic quantum states in atomic-size structures. The project involved developing first-principles calculations that can predict coherent lifetimes of defects; theoretical schemes for identifying new platforms suitable for studying the dynamics of driven quantum systems; and optical and electronic spectroscopy characterization of novel defects and heterostructures with the aim of creating the next generation of materials for quantum science. This report provides a brief background into the project, highlights some of the major achievements accomplished within the projects and lists the supported publications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Progress Report on Model Development for Aerosol Transport through Divergent Cracks Paths

This report summarizes the progress in developing a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose of this effort is to introduce to a user community—involving researchers, regulators, and industry—a generic, reliable numerical model for the prediction of aerosol transport while accounting for potential deposition and plugging of the leak paths to model spent nuclear fuel (SNF) release from postulated Stress Corrosion Cracks in canisters under storage or transportation. Current work focuses on expanding the model to predict aerosol (and gas) flow through complex microchannel (nozzle) geometries in addition to the rectangular and cylindrical geometries, as we approach more realistic Stress Corrosion Crack conditions. In this regard, a divergent nozzle geometry, used at Sandia National Laboratories (SNL) for testing aerosol release and retention (Durbin et al. 2021), was added to the model. The model was then validated with blowdown data for the particular microchannel, from experiments conducted at SNL. The report also presents preliminary non-benchmarked aerosol penetration fraction and mass flow rate for a monodisperse 10-micron (AED) particle concentration of 1.7e-08 kg/m 3 being released through the divergent microchannel geometry from the SNL’s aerosol experimental tank setup.Future work will involve validating the aerosol model and updating the model’s Graphical User Interface (GUI) to include the divergent nozzle geometry. This is expected to help stakeholders perform quick and easy first principles calculations without needing to understand the underlying MATLAB script.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Modeling and Optimization of Zeolites for Contaminant Removal from Coal Combustion Impoundment Leachates

Coal Combustion Residual (CCR) landfills in the U.S. can lead to elevated contaminant concentrations in groundwater and runoff, particularly of arsenic, boron, and selenium. Zeolites can be readily synthesized from materials contained within the coal combustion ash impoundments and can be cation-exchanged to enhance adsorption capacity, selectivity, and reactivity. However, the optimization of zeolites becomes a daunting task when accounting for the variety of Si:Al ratios, the species of extra-framework cations present, and the zeolite pore structure. Molecular simulations provide methods to study and guide the design of zeolites for the sorption of contaminants from aqueous solution. In this work, models that predict the sorption of arsenic, boron, and selenium from water using cation-exchanged zeolites. Because of the lack of experimental adsorption data for these species of contaminants, models were parametrized to reproduce the results of first-principles calculations and then used to predict sorption for zeolites for a dataset containing 6000 combinations of zeolites and sorption conditions. Machine learning was used to train a model to predict sorption for materials in this database based on the results of the molecular simulations. Next, a genetic algorithm was used to optimize zeolites for the removal of each contaminant from aqueous solution for individual impoundment sites based on contaminant concentrations reported by the Electric Power Research Institute.

Findley, John↗
DOE OSTI2024

Next-Generation Materials Design: Quantum Mechanics and Data-Driven Modeling

The future of materials design is rapidly advancing through the combination of quantum mechanics and data-driven modeling. These approaches integrate quantum principles with advanced data analysis, enabling precise insights into material behavior. This talk will highlight recent progress in using these methods for computational design, particularly in high-entropy alloy catalysts, emphasizing the role of hierarchical machine-learning architectures for accurate predictions. Additionally, I will discuss our work on developing machine learning interatomic potentials (MLPs) for single-element metals, metal oxides, and alloys under extreme conditions, focusing on melting behavior and phase properties at high temperatures and pressures. We have also refined our MLP models to capture dynamic surface interactions, such as CO2 and CO adsorption on MgO, using both static and molecular dynamics simulations. These models maintain high accuracy while significantly reducing computational costs compared to first-principles calculations. By enabling efficient and accurate simulations, this work supports broader community adoption, optimizes datasets for materials discovery, and extends the accessible time, size, and environmental conditions beyond the limits of experiments and traditional simulations.

machine learning↗
DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effect of Dilute Magnetism in a Topological Insulator

Three-dimensional (3D) topological insulator (TI) has emerged as a unique state of quantum matter and generated enormous interests in condensed matter physics. The surfaces of a 3D TI consist of a massless Dirac cone, which is characterized by the Z 2 topological invariant. Introduction of magnetism on the surface of a TI is essential to realize the quantum anomalous Hall effect and other novel magneto-electric phenomena. Here, by using a combination of first-principles calculations, magneto-transport and angle-resolved photoemission spectroscopy (ARPES), we study the electronic properties of gadolinium (Gd)-doped Sb 2 Te 3 . Our study shows that Gd doped Sb 2 Te 3 is a spin-orbit-induced bulk band-gap material, whose surface is characterized by a single topological surface state. Our results provide a new platform to investigate the interactions between dilute magnetism and topology in magnetic doped topological materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of microstructure and neutron irradiation defects on deuterium retention in SiC

Retention of hydrogen isotopes is a critical concern for operating fusion reactors as retained tritium both activates components and removes scarce fuel from the fuel cycle. Radiation-induced displacement damage in SiC influences the retention of hydrogen isotopes compared to pristine SiC. Deuterium retention in neutron irradiated high purity SiC has been compared to different microstructures of non-irradiated high purity SiC using thermal desorption spectroscopy after gas charging and low energy ion implantation. Experimental results show lower deuterium retention in single crystal SiC than in polycrystal SiC indicating that grain boundaries are key trapping features in unirradiated SiC. Deuterium is released at lower temperatures in neutron irradiated polycrystal SiC compared to pristine polycrystal SiC, suggesting weaker trapping by radiation-induced defects compared to grain boundary trapping sites in the pristine materials. Low energy ion implantation caused a high deuterium release temperature, highlighting the sensitivity of deuterium release behaviour to radiation defect characteristics. First principles calculations have been conducted to identify energetically favourable trapping sites in SiC at the H ABc V Si and H TSi V C complexes, and migration barriers between interstitial sites. This helps interpret experimental results and derive effective diffusivity of hydrogen isotopes in SiC in the presence of vacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase Stability in U-6Nb Alloy Doped with Ti from the First Principles Theory

First-principles calculations within the density-functional-theory (DFT) approach are conducted in order to explore and explain the effect of small amounts of titanium on phase stability in the U-6Nb alloy. During rapid quenching from high to room temperature, metastable phases α' (orthorhombic), α" (monoclinic), and γ 0 (tetragonal) can form, depending on Nb concentration. Important mechanical properties depend on the crystal structure and, therefore, an understanding of the effect of impurities on phase stability is essential. Insights on this issue are obtained from quantum-mechanical DFT calculations. The DFT framework does not rely on any material-specific assumptions and is therefore ideal for an unbiased investigation of the U-Nb system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Topological Dirac Semimetal Phase in Bismuth Based Anode Materials for Sodium-Ion Batteries

Bismuth has recently attracted interest in connection with Na-ion battery anodes due to its high volumetric capacity. It reacts with Na to form Na 3 Bi which is a prototypical Dirac semimetal with a nontrivial electronic structure. Density-functional-theory based first-principles calculations are playing a key role in understanding the fascinating electronic structure of Na 3 Bi and other topological materials. In particular, the strongly-constrained-and-appropriately-normed (SCAN) meta-generalized-gradient-approximation (meta-GGA) has shown significant improvement over the widely used generalized-gradient-approximation (GGA) scheme in capturing energetic, structural, and electronic properties of many classes of materials. Here, we discuss the electronic structure of Na 3 Bi within the SCAN framework and show that the resulting Fermi velocities and s-band shift around the Γ point are in better agreement with experiments than the corresponding GGA predictions. SCAN yields a purely spin-orbit-coupling (SOC) driven Dirac semimetal state in Na 3 Bi in contrast with the earlier GGA results. Our analysis reveals the presence of a topological phase transition from the Dirac semimetal to a trivial band insulator phase in Na 3 Bi x Sb 1-x alloys as the strength of the SOC varies with Sb content, and gives insight into the role of the SOC in modulating conduction properties of Na 3 Bi.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Site Occupancy Preference and Magnetic Properties in Nd 2 (Fe,Co) 14 B

Partial replacement of Fe by Co is an effective method to increase Curie temperature (T C ) , which improves the thermal stability of magnetic properties in Nd 2 Fe 14 B-based permanent magnets. The correlation between Fe substitution and magnetic properties has been studied in Nd 2 (Fe,Co) 14 B via a first-principles calculation. The calculated Fe substitution energies indicate that the Co atoms avoid the 8 j 2 site, which agrees with the experiments. The Co atoms are ferromagnetically coupled with Fe sublattice and show magnetic moments of about 1.2 to 1.7 μ B at different crystallographic sites, less than that of Fe (2.1–2.7 μ B ), resulting in the decrease in total magnetization at ground state (0 K) with increasing Co content. The effective exchange interaction parameter, derived from the energy difference between varied magnetic structures, increases from 7.8 meV to 17.0 meV with increasing Co content from x = 0 to x = 14 in Nd 2 Fe 14–x Co x B. This change in the effective exchange interaction parameter is responsible for the enhancement of T C in Nd 2 (Fe,Co) 14 B. The total magnetization at 300 K, derived from mean-field theory, shows a peak maximum value at x = 1 in Nd 2 Fe 14–x Co x B. The phenomenon results from the interplay between the reduction of the magnetic moment in the Fe(Co) sublattice and the enhancement of T C with increasing Co content.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson’s theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs–Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born–Oppenheimer molecular dynamics performed at a single temperature. In conclusion, this advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Enhanced Phase Stability of Sm 2 (Fe, Al) 17 C x

Aluminum doping can improve the phase stability of metastable compound Sm 2 Fe 17 C x with a high carbon content (x > 1.5). We investigated the preferential site substitution of Al, chemical bonding, and structural stability in Sm 2 (Fe,Al) 17 C 3 using first-principle calculations. Our results reveal a strong correlation between the preferential substitution of Fe by Al and the atomic site chemical environment, which affects the overall phase stability. Specifically, Al preferentially occupies the 9d site in Sm 2 (Fe,Al) 17 C 3 . At the same time, Al prefers the site 6c in its parent phase Sm 2 (Fe,Al) 17 . Partial replacement of Fe with Al leads to a more negative formation energy, indicating enhanced thermodynamic stability. Crystal Orbital Hamilton Population (COHP) and Crystal Orbital Bond Index (COBI) analysis suggest that insertion of carbon weakens the bonding strength of Sm-Fe (18f) and Sm-Fe (18h), resulting in metastability of Sm 2 Fe 17 C x . Doping Al strengthens Al-Fe, Al-Sm, Sm-Fe (18f, 18h) and Fe–C bonding in Sm 2 (Fe,Al) 17 C 3 , as revealed by calculated COHP and COBI. These effects contribute to improved phase stability in the Al-doped 2:17 interstitial compound.

chemical bonding↗
DOE OSTI2023

Vibronic Relaxation Pathways in Molecular Spin Qubit Na 9 [Ho(W 5 O 18 ) 2 ]·35H 2 O under Pressure

In order to explore how spectral sparsity and vibronic decoherence pathways can be controlled in a model qubit system with atomic clock transitions, we combined diamond anvil cell techniques with synchrotron-based far infrared spectroscopy and first-principles calculations to reveal the vibrational response of Na 9 [Ho(W 5 O 18 ) 2 ]·35H 2 O under compression. Because the hole in the phonon density of states acts to reduce the overlap between the phonons and f manifold excitations in this system, we postulated that pressure might move the HoO 4 rocking, bending, and asymmetric stretching modes that couple with the M J = ±5, ±2, and ±7 levels out of resonance, reducing their interactions and minimizing decoherence processes, while a potentially beneficial strategy for some molecular qubits, pressure slightly hardens the phonons in Na 9 [Ho(W 5 O 18 ) 2 ]·35H 2 O and systematically fills in the transparency window in the phonon response. The net result is that the vibrational spectrum becomes less sparse and the overlap with the various M J levels of the Ho 3+ ion actually increases. These findings suggest that negative pressure, achieved using chemical means or elongational strain, could further open the transparency window in this rare earth-containing spin qubit system, thus paving the way for the use of device surfaces and interface elongational/compressive strains to better manage decoherence pathways.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Designing Quaternary and Quinary Refractory-Based High-Entropy Alloys: Statistical Analysis of Their Lattice Distortion, Mechanical, and Thermal Properties

The rapid evolution in materials science has resulted in a significant interest in high-entropy alloys (HEAs) for their unique properties. This study focuses on understanding both quaternary and quinary body-centered cubic (BCC) of 12 refractory-based HEAs, and on analysis of their electronic structures, lattice distortions, mechanical, and thermal properties. A comprehensive assessment is undertaken by means of density functional theory (DFT)-based first principles calculations. It is well known that multiple constituents lead to notable lattice distortions, especially in quinary HEAs. This distortion, in turn, has significant implications on the electronic structure that ultimately affect mechanical and thermal behaviors of these alloys such as ductility, lattice thermal conductivity, and toughness. Our in-depth analysis of their electronic structures revealed the role of valence electron concentration and its correlation with bond order and mechanical properties. Local lattice distortion (LD) was investigated for these 12 HEA models. M1 (WTiVZrHf), M7 (TiZrHfW), and M12 (TiZrHfVNb) have the highest LD whereas the models M3 (MoTaTiV), M5 (WTaCrV), M6 (MoNbTaW), and M9 (NbTaTiV) have the less LD. Furthermore, we investigated the thermal properties focusing on Debye temperature (ΘD), thermal conductivity (κ), Grüneisen parameter (γα), and dominant phonon wavelength (λdom). The NbTaTiV(M9) and TiVNbHf(M10) models have significantly reduced lattice thermal conductivities (κL). This reduction is due to the mass increase and strain fluctuations, which in turn signify lattice distortion. The findings not only provide an understanding of these promising materials but also offer guidance for the design of next-generation HEAs with properties tailored for potential specific applications.

Materials Science↗
DOE OSTI2021

Lithium Diffusion in Silicon Encapsulated with Graphene

The model of a graphene (Gr) sheet putting on a silicon (Si) substrate is used to simulate the structures of Si microparticles wrapped up in a graphene cage, which may be the anode of lithium-ion batteries (LIBS) to improve the high-volume expansion of Si anode materials. The common low-energy defective graphene (d–Gr) structures of DV5–8–5, DV555–777 and SV are studied and compared with perfect graphene (p–Gr). First-principles calculations are performed to confirm the stable structures before and after Li penetrating through the Gr sheet or graphene/Si-substrate (Gr/Si) slab. The climbing image nudged elastic band (CI-NEB) method is performed to evaluate the diffusion barrier and seek the saddle point. The calculation results reveal that the d–Gr greatly reduces the energy barriers for Li diffusion in Gr or Gr/Si. The energy stability, structural configuration, bond length between the atoms and layer distances of these structures are also discussed in detail.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anisotropic Thermal Conductivity of Inkjet-Printed 2D Crystal Films: Role of the Microstructure and Interfaces

Two-dimensional (2D) materials are uniquely suited for highly anisotropic thermal transport, which is important in thermoelectrics, thermal barrier coatings, and heat spreaders. Solution-processed 2D materials are attractive for simple, low-cost, and large-scale fabrication of devices on, virtually, any substrate. However, to date, there are only few reports with contrasting results on the thermal conductivity of graphene films, while thermal transport has been hardly measured for other types of solution-processed 2D material films. In this work, inkjet-printed graphene, h-BN and MoS2 films are demonstrated with thermal conductivities of ∼10 Wm−1K−1 and ∼0.3 Wm−1K−1 along and across the basal plane, respectively, giving rise to an anisotropy of ∼30, hardly dependent on the material type and annealing treatment. First-principles calculations indicate that portion of the phonon spectrum is cut-off by the quality of the thermal contact for transport along the plane, yet the ultra-low conductivity across the plane is associated with high-transmissivity interfaces. These findings can drive the design of highly anisotropic 2D material films for heat management applications.

2D materials↗