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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Beam Wave Considerations for Optical Link Budget Calculations

The bounded beam wave nature of electromagnetic radiation emanating from a finite size aperture is considered for diffraction-based link power budget calculations for an optical communications system. Unlike at radio frequency wavelengths, diffraction effects are very important at optical wavelengths. In the general case, the situation cannot be modeled by supposing isotropic radiating antennas and employing the concept of effective isotropic radiated power. It is shown here, however, that these considerations are no more difficult to treat than spherical-wave isotopic based calculations. From first principles, a general expression governing the power transfer for a collimated beam wave is derived and from this are defined the three regions of near-field, first Fresnel zone, and far-field behavior. Corresponding equations for the power transfer are given for each region. It is shown that although the well-known linear expressions for power transfer in the far-field hold for all distances between source and receiver in the radio frequency case, nonlinear behavior within the first Fresnel zone must be accounted for in the optical case at 1550 nm with typical aperture sizes at source/receiver separations less that 100 km.

Wave Propagation↗
NASA NTRS2016

A Physics-Based Approach to Hypersonic Nonequilibrium Chemistry

NASA has undertaken a concerted effort to generate a complete model for the chemical reactions and other collisional processes that take place in the bow shock layer and wake of spacecraft entering the atmospheres of Earth and Mars. This new model is based on first-principles physics calculations: interaction potentials between colliding species obtained from accurate quantum mechanics calculations and reaction rate coefficients and cross sections based on classical mechanics simulations of heavy particle collisions.This presentation will give an introduction to the underlying theory and computational methods used in this work and describe the range of chemical and physical processes that take place in the close vicinity of spacecraft during atmospheric entry. I will also present some of our recent results for reactions involving nitrogen, carbon monoxide and CO2.

Jaffe, Richard↗
NASA NTRS2021

Spin and Accretion Rate Dependence of Black Hole X-Ray Spectra

We present a survey of how the spectral features of black hole X-ray binary systems depend on spin, accretion rate, viewing angle, and Fe abundance when predicted on the basis of first-principles physical calculations. The power-law component hardens with increasing spin. The thermal component strengthens with increasing accretion rate. The Compton bump is enhanced by higher accretion rate and lower spin. The Fe Kα equivalent width grows sublinearly with Fe abundance. Strikingly, the Kα profile is more sensitive to accretion rate than to spin because its radial surface brightness profile is relatively flat, and higher accretion rate extends the production region to smaller radii. The overall radiative efficiency is at least 30%–100% greater than as predicted by the Novikov–Thorne model.

Brooks E Kinch↗
DOE OSTI2021

Competing energetic states in γ – Fe 2 W O 6 with strong spin-charge-lattice coupling

We report magnetic and electronic properties of γ – Fe 2 W O 6 via neutron powder-diffraction measurements and first-principles density-function theory calculations. We reveal a magnetic ground state with two phases coexisting in which the minor phase is the same as an earlier report, despite the fact that both materials studied have the same space group but with slightly different atomic positions. Interestingly, both magnetic phases are well captured by first-principles calculations. Furthermore, we show that these two magnetic phases are correlated with electronic properties, with one being insulating and the other being metallic. These features suggest that γ – Fe 2 W O 6 exhibits competing energetic states in which spin, charge, and lattice degrees of freedom are strongly coupled to each other.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Thermodynamics Modeling for Actinide Monocarbides and Mononitrides from First Principles

The high-temperature thermodynamical properties for the actinide monocarbides and mononitrides ThC, ThN, UC, UN, PuC, and PuN are calculated from first-principles electronic-structure theory. The electronic structure is modeled with density-functional theory (DFT) and is fully relativistic, including the spin-orbit interaction. Furthermore, the DFT is extended to account for orbital–orbital interactions, by means of a parameter-free orbital-polarization (OP) technique, that has proven to be essential for the 5f electrons in plutonium. Strong anharmonicity and the temperature dependence of the lattice vibrations are captured with the self-consistent ab initio lattice dynamics (SCAILD) method. The calculated free energies and heat capacities are compared to published results from quasi-harmonic (QH) theory, and experiments, where available. For the uranium and plutonium compounds, we make use of CALPHAD assessments to help evaluate the theory. Generally, our anharmonic relativistic approach compares well with both CALPHAD and experiments. For the thorium compounds, our theory is in good accord with QH modeling of the free energy at lower temperatures but for the heat capacity the comparison is less favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Melting and transport properties of Al 2⁢ O 3 at extreme conditions

The high-pressure, temperature phase diagram and transport properties of materials are of broad interest to planetary sciences and high-energy-density sciences and applications. Alumina (Al 2 O 3 ) or its various forms (e.g., solutions with other oxides or silicates) are important constituents in Earth’s and super-Earths’ mantles, common window materials in dynamic compression experiments, and standard pressure calibrators in diamond-anvil-cell experiments. Its structures and transport properties are of particular importance but have not been well studied at above 100 GPa pressures. Based on extensive first-principles molecular dynamics calculations, we obtain atomic level insights on structural differences and accurate results on electrical conductivity (σ dc ), thermal conductivity (κ), and reflectivity of Al 2 O 3 in several phases at pressures of ~130–1300 GPa and temperatures of 4000–20,000 K. We find the solid-to-liquid state changes of Al 2 O 3 are accompanied by an insulator-to-semimetal transition with enhanced σ dc and κ, which is similar to other silicates and oxides and can facilitate magnetic field generation in deep interiors of rocky planets. We also find the solid-to-solid transition (from Rh 2 O 3 (II) to CaIrO 3 ) is accompanied by negligible changes in σ dc and κ at 160 GPa, contrary to previous expectations, whereas another higher pressure phase transition (CaIrO 3 to U 2 S 3 at 450 GPa) is associated with increases in σ dc and κ. Furthermore, we show the transport properties can be enhanced by oxygen vacancy in Al 2 O 3 , but are not significantly affected by forming solutions with MgSiO 3 , particularly when in the solid states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Lattice anharmonicity effects in fluorite oxide single crystals and anomalous increase in phonon lifetime in ceria at elevated temperature

Here, we investigate the temperature dependence of the frequency and linewidth of the triply degenerate T 2g zone-centered optical phonon in flux-grown ceria and hydrothermally synthesized thoria single crystals from room temperature to 1273 K using Raman spectroscopy. Both crystals exhibit an expected increase in the phonon linewidth with temperature due to enhanced phonon–phonon scattering. However, ceria displays an anomalous linewidth reduction in the temperature range of 1023–1123 K. First-principles phonon linewidth calculations considering cubic and quartic phonon interactions within temperature-independent phonon dispersion fail to describe this anomaly. A parameterization of the temperature-dependent second-order interatomic force constants based on previously reported phonon dispersion measured at room and high temperatures predicts a deviation from the monotonic linewidth increase, albeit at temperatures lower than those observed experimentally for ceria. The qualitative agreement in the trend of temperature-dependent linewidth suggests that lattice anharmonicity-induced phonon renormalization plays a role in phonon lifetime. Specifically, a change in the overlap between softened acoustic and optical branches in the dispersion curve reduces the available phonon scattering phase space of the Raman-active mode at the zone center, leading to an increased phonon lifetime within a narrow temperature interval. These findings provide insights into higher-order anharmonic interactions in ceria and thoria, motivating further investigations into the role of anharmonicity-induced phonon renormalization on phonon lifetimes at high temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗
DOE OSTI2020

First-principles study on structural mechanical and thermodynamic properties of HfMoTaTiZr

This study calculates the phase diagram of HfMoTaTiZr using the TCHEA1 database in Thermo-Calc software. Furthermore, it also investigates the structural and mechanical properties of HfMoTaTiZr by applying first-principles calculations. The elastic and other mechanical properties such as Young’s modulus, pugh ratio, hardness, bulk modulus was derived and discussed. These parameters were found to be in good harmony with previously available experimental reports. Vickers hardness of HfMoTaTiZr predicted from the Tians model agrees with the available experimental result. This study also investigates the thermal coefficient of linear expansion and vibrational heat capacity of the HfMoTaTiZr using the Debye–Grüneisen theory.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Near-room-temperature ferromagnetic ordering in the pressure-induced collapsed-tetragonal phase in SrC⁢o 2 ⁢P 2

We present high-pressure electrical transport, magnetization, and single-crystal x-ray diffraction data on SrC⁢o 2 ⁢P 2 single crystals, revealing a transition from the uncollapsed tetragonal (ucT) to the collapsed tetragonal (cT) phase >∼10 GPa. This structural transition is accompanied by ferromagnetic (FM) ordering with a Curie temperature up to 260 K. First-principles-based density functional theory calculations show that there is a first-order transition between ucT and cT phases, with an onset ∼10 GPa as well as the appearance of FM ordering in the cT phase. Above ∼30 GPa, the experimental signatures of the magnetic ordering vanish in a first-order-like manner, consistent with the theoretical calculation results, indicating that SrC⁢o 2 ⁢P 2 is another example of the avoidance of quantum criticality in FM intermetallic compounds with almost room-temperature critical temperature. SrC⁢o 2 ⁢P 2 provides clear evidence that the structural, electronic, and magnetic properties associated with the cT transition are strongly entangled and are not only qualitatively captured by our first-principles-based calculations but are quantitatively reproduced as well.

density functional theory↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First Principles Evaluation of Phase Stability in the In-Sn Binary System

Abstract The In-Sn binary alloy system exhibits several unusual features that challenge crystallographic and thermodynamic expectations. We combine first principles total energy calculation with simple thermodynamic modeling to address two key points. First, we evaluate energies along the Bain path to interpret the discontinuous transition between the phases α-In (Pearson type tI2) and β-In 3 Sn (also Pearson type tI2) that are identical in symmetry. Second, we demonstrate that the solid solution phases β-In 3 Sn and γ-InSn 4 (Pearson type hP1) exist at high temperatures only, and they exhibit eutectoid decompositions at low temperatures.

Chemistry↗
DOE OSTI2021

Orbital ordering in the layered perovskite material CsVF 4

In strongly correlated electronic systems, several novel physical properties are induced by the orbital degree of freedom. In particular, orbital degeneracy near the Fermi level leads to spontaneous symmetry breaking, such as the nematic state in FeSe and the orbital ordering in several perovskite systems. Here, the novel layered perovskite material CsVF 4 , with a 3 d 2 electronic configuration, was systematically studied using density-functional theory and a multiorbital Hubbard model within the Hatree-Fock approximation. Furthermore, our results show that CsVF 4 should be magnetic, with a G-type antiferromagnetic arrangement in the a b plane and weak antiferromagnetic exchange along the c axis, in agreement with experimental results. Driven by the Jahn-Teller distortion in the VF 6 octahedra that shorten the c axis, the system displays an interesting electron occupancy d x y 1 ( d x z d y z ) 1 corresponding to the lower nondegenerate d x y orbital being half-filled and the other two degenerate d y z and d x z orbitals sharing one electron per site. We show that this degeneracy is broken and a novel d y z / d x z staggered orbital pattern is here predicted by both the first-principles and Hubbard model calculations. This orbital ordering is driven by the electronic instability associated with degeneracy removal to lower the energy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab initio study of tungsten-based alloys under fusion power-plant conditions

Tungsten (W) is considered a leading candidate for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. However, tungsten also presents a very low fracture toughness, primarily associated with intergranular failure and bulk plasticity, limiting its applications. In recent years, several families of tungsten-based alloys have been explored to overcome the aforementioned limitations of pure tungsten. These include tungsten-based high-entropy alloys (W-HEAs) and tungsten-based Self-passivating Metal Alloys with Reduced Thermo-oxidation or “SMART alloys” (W-SAs). Given their proximity to the plasma, it is crucial to understand how the exposure of these candidate plasma-facing materials (PFMs) to the neutron fluxes expected in fusion reactors impacts their material behavior over time. In this work, we present a computational approach that combines inventory codes and first-principles DFT electronic structure calculations to understand the behavior of transmuting tungsten-based PFMs. In particular, we calculate the changes in the chemical composition, production uncertainties, the elastic and ductility properties, and the density of states for five tungsten-based PFMs when exposed to EU-DEMO fusion first wall conditions for ten years.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Structure sensitivity of ammonia electro-oxidation on transition metal surfaces: A first-principles study

Ammonia electro-oxidation is a promising catalytic reaction for application in alkaline fuel cells. In this paper we study the trends in this reaction on the (1 0 0) facet of eight fcc transition metals: Au, Ag, Cu, Pt, Pd, Ni, Ir, and Rh. We calculate from first-principles (DFT-GGA-PW91) the energetics for two mechanisms: (i) the N + N mechanism and (ii) the Gerischer-Mauerer mechanism. The onset potentials for both mechanisms are provided for both (1 0 0) and (1 1 1) facets. We also calculate the activation energies for the N-N bond-making non-Faradaic events on Cu, Pt, Pd, Ni, Ir, and Rh. We discover that N* – which is considered a poison for Pt(1 1 1) – is a reactive intermediate on Pt(1 0 0), allowing for a facile dimerization to N2 on Pt(1 0 0) at high potentials. We draw similar analyses for the structure sensitivity on the other metals, and we conclude with some general principles to guide catalyst design for this reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗