Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY OF STATES”

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

Determination of Celestial Body Principal Axes via Gravity Field Estimation

The determination of physical parameters of celestial bodies provides critical navigational and scientific information. Determining the mass, rotation state, and density distribution is an important task of the navigation team on an operational mission. One representation of the internal density distribution, spherical harmonics, may be leveraged to compute the orientation of the principal axis frame of the celestial body, which in turn informs the dynamics of the body’s motion. Using data from NASA’s Dawn mission to the massive asteroids (1) Ceres and (4) Vesta, a method is proposed that refines the knowledge of the principal axis frame by utilizing estimated spherical harmonic coefficients. Applying an iterative rotation scheme to the solved-for body-fixed frame leads to convergence on a frame that is very close to the actual dynamical principal axis frame.

Takahashi, Yu↗
NASA NTRS2021

What it Takes to Compute Highly Accurate Rovibrational Line Lists for Use in Astrochemistry

We review the Best Theory + Reliable High-resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is two-fold: a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH3, or have a high density of states such as SO2; and b) characterization of the atmospheres of exoplanets which will be one of the main areas of research in astronomy in the coming decades. The first motivation arises due to the need to eliminate lines due to common molecules in an astronomical observation in order to identify lines from new molecules, while the second motivation arises due to the need to obtain accurate molecular opacities in order to characterize the atmosphere of an exoplanet. The BTRHE strategy first consists of using high quality ab initio quantum chemical methods to obtain a global potential energy surface (PES) and dipole moment surface (DMS) that contains the proper physics. The global PES is then refined using a subset of the reliable high-resolution experimental data. The refined PES then gives energy level predictions to an accuracy similar to the reproduction accuracy of the experimental data used in the refinement step in the interpolation region (i.e., within the range of the experimental data used in the refinement step). The accuracy of the energy levels will slowly degrade as they are extrapolated to spectral regions beyond the high-resolution experimental data used in the refinement step. However, because the degradation is slow, the predicted energy levels can be used to assign new high-resolution experiments, and the data from these can then be used in a subsequent refinement step. In this way, the global PES eventually can yield highly accurate energy levels for all desired spectral regions including to very high energies and high J values. We show that IR intensities computed with the BTRHE rovibrational wavefunctions and the DMS can be very accurate provided one has minimized the fitting error of the DMS and tested the completeness of the DMS. Some examples of our work on NH3, CO2, and SO2 are given to highlight the usefulness of the BTRHE strategy and to provide ideas on how to further improve its predictive power in the future. In particular, it is shown how successive refinement steps, once new high-resolution data is available, can lead to PESs that yield highly accurate transition energies to larger spectral regions. The importance of including non-adiabatic corrections to reduce the J-dependence of errors for H-containing molecules is shown with work on NH3. Another very important aspect of the BTRHE approach is the consistency across isotopologues, which allows for highly accurate line lists for any isotopologue once one is obtained for the main isotopologue (which has more high-resolution data available for refinement).

Xinchuan Huang↗
NASA NTRSDate not supplied

Globally Optimal Band Structure for Thermoelectrics in Realistic Systems

One of the ultimate questions in the field of thermoelectrics is the optimal band structure for the figure of merit. The observation made here is that a linear band dispersion, of any dimension, under acoustic-phonon-deformation-potential scattering comprises the theoretical regimen for delivering boxcar transport distribution, which has been mathematically proven to deliver optimum figure of merit. This is guaranteed by the interplay of the energy-dependencies of groupvelocity and density of states of linear dispersion as well as the energy-behavior of deformation-potential scattering. Then linear dispersion, capped with optimum bandwidth and combined with provisions of Park et al, npj Comput. Mater. 7, 43 2021, is the globally optimal band structure for thermoelectricity.

Junsoo Park↗
DOE OSTI2020

Higher-order phonon scattering: advancing the quantum theory of phonon linewidth, thermal conductivity and thermal radiative properties

Phonon scattering plays a central role in the quantum theory of phonon linewidth, which in turn governs important properties including infrared spectra, Raman spectra, lattice thermal conductivity, thermal radiative properties, and also signifi- cantly affects other important processes such as hot electron relaxation. Since Maradudin and Fein’s classic work in 1962, three-phonon scattering had been considered as the dominant intrinsic phonon scattering mechanism and has seen tremendous advances. However, the role of the higher-order four-phonon scattering had been persistently unclear and so was ignored. The tremendous complexity of the formalism and computational challenges stood in the way, prohibiting the direct and quantitative treatment of four-phonon scattering. In 2016, a rigorous four-phonon scattering formalism was developed, and the prediction was realized using empirical potentials. In 2017, the method was extended using first-principles calculated force constants, and the thermal conductivities of boron arsenides (BAs), Si and diamond were predicted. The predictions for BAs were later confirmed by several independent experiments. Four-phonon scattering has since been investigated in a range of materials and established as an important intrinsic scattering mechanism for thermal transport and radiative properties. Specifically, four-phonon scattering is important when the fourth-order scattering potential or phase space becomes relatively large. The former scenario includes: (i) nearly all materials when the temperature is high; (ii) strongly anharmonic (low thermal conductivity) materials, including most rocksalt compounds, halides, hydrides, chalcogenides and oxides. The latter scenario includes: (iii) materials with large acoustic–optical phonon band gaps, such as XY compounds with a large atomic mass ratio between X and Y; (iv) two- dimensional materials with reflection symmetry, such as single-layer graphene, single-layer boron nitride and carbon nanotubes; and (v) phonons with a large density of states, such as optical phonons, which are important for Raman, infrared and thermal radiative properties. Four-phonon scattering is expected to gain broad interest in various technologically important materials for thermoelectrics, thermal barrier coatings, thermal energy storage, phase change, nuclear power, ultra-high temperature ceramics, infrared spectra, Raman spectra, radiative transport, hot electron relaxation and radiative cooling. Four-phonon scattering has been, and will continue to be, established as an important intrinsic phonon scattering mechanism beyond three-phonon scattering. The prediction of four-phonon scattering will transition from a breakthrough to a new routine in the next decade.

Feng, Tianli↗
DOE OSTI2022

Comparative critical mass calculations for NNL and ENDF/B-VIII.0 - Zirconium hydride thermal neutron scattering laws

Zirconium hydride (ZrH{sub x}) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high-enriched uranium (HEU) and ZrH{sub x} in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of {sup 235}U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH{sub x}. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Development of Thermal Scattering Law and Cross Sections of Hydrogen in Paraffin Wax

Paraffin wax is frequently used as a neutron moderator and shielding material. The main component of paraffin wax is straight-chain alkanes (n-alkanes). The deposition of paraffin wax is primarily attributed to the crystallization of n-alkanes. It is important to gain a deeper understanding of the mechanisms underlying the behavior of paraffin wax, which would impact the thermal scattering Law (TSL) and cross sections and affect the analysis of neutronic and critical systems. In this work, a classical molecular dynamics (CMD) simulation model was used in LAMMPS to create the TSL and cross sections at room temperature and pressure. To generate the required velocity autocorrelation functions (VACF), previously published data were used to validate the approach and models by comparing them with key model parameters. The phonon density of state (DOS) was calculated using Fourier transformation of the normalized VACF. This DOS was used as the primary input to estimate the TSL (S($a$, $β$)) and cross sections of hydrogen in paraffin wax. The TSL and cross sections of hydrogen were estimated using the Full Law Analysis Scattering System Hub (FLASSH) code. The cross section of hydrogen in paraffin wax is consistent with other hydrocarbon materials such as polyethylene with deviations due to structure in the lowest energy region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Generation of Enrichment-Dependent Thermal Neutron Scattering Data

This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1974

Calculations of the early evolution of Jupiter

The evolution of the protoplanet Jupiter is followed, using a hydrodynamic computer code with radiative energy transport. Jupiter is assumed to have formed as a subcondensation in the primitive solar nebula at a density just high enough for gravitational collapse to occur. The initial state has a density of 0.0015 nanograms per cu cm and a temperature of 43 K; the calculations are carried to an equilibrium state where the central density reaches 0.5 g per cu cm and the central temperature reaches 25,000 K. During the early part of the evolution the object contracts in quasi-hydrostatic equilibrium; later on hydrodynamic collapse occurs, induced by the dissociation of hydrogen molecules. After dissociation is complete, the planet regains hydrostatic equilibrium with a radius of a few times the present value. Further evolution beyond this point is not treated here; however the results are consistent with the existence of a high-luminosity phase shortly after the planet settles into its final quasi-static contraction.-

Bodenheimer, P.↗
DOE OSTI2025

Fundamental Understanding of “Fresh” Lithium Nucleation and Growth in Sulfide-Based Anode-Free Solid-State Batteries: Effects of Substrate, Current Density, and Li + Supply

Understanding the electrochemical extraction and deposition of lithium (Li) from cathode is crucial for advancing anode-free solid-state batteries (AFSSBs). Herein, cryo-transmission electron microscopy (cryo-TEM) and electrochemical studies are employed to investigate how current collector surface properties, current densities, and cathode loadings influence the morphology of fresh electrochemically deposited Li and the electrochemical performance in sulfide-based AFSSBs. Cryo-TEM reveals that Cu current collectors induce irregular, dendritic Li deposits due to their lithiophobic nature and reactivity with Li 5.4 PS 4.4 Cl 1.6 (LPSC), while Ni and Au facilitate more uniform, planar-like Li growth. The morphology of the deposited Li also depends on current density: higher rates produce smaller, porous particles versus larger, denser deposits at lower rates. Importantly, for the first time, we discovered that low cathode loadings lead to poor cycling stability due to insufficient Li + supply for complete anode coverage, contrary to conventional solid-state batteries with preloaded Li metal anodes. This finding establishes a design principle where adequate cathode Li + reservoirs are essential for anode interface stability in AFSSBs. Overall, this work elucidates the interplayed effect of current collectors, current densities, and cathode loadings on Li morphology and cycle stability, offering fundamental insights into cathode-derived Li behavior and practical guidelines for optimizing AFSSBs performance through nucleation control and interface engineering.

25 ENERGY STORAGE↗
DOE OSTI2023

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

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

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1979

Release adiabat measurements on minerals - The effect of viscosity

The paper examines the effect of viscosity in the release adiabat measurements on minerals. The current inversion of pressure-particle velocity data for release from a high-pressure shock state to a pressure-density path usually depends critically upon the assumption that the release process is isentropic. It has been shown that the effective viscosity for geological materials below stresses of 150 GPa must be at least 1000 kg/m/s so that viscous work in the shock state remains small compared to mechanical work recovered upon adiabatic rarefaction. The magnitude of the shear stress in the shock state in minerals and viscosities of engineering materials shocked to pressures below 150 GPa show effective viscosities of about 1000 kg/m/s or less indicating that the conditions for isentropic release of materials from shock states are achieved and that the Riemann integral can be applied to obtain pressure-density states along the release adiabats of minerals.

Jeanloz, R.↗
NASA NTRS1979

Release adiabat measurements on minerals: The effect of viscosity

The current inversion of pressure-particle velocity data for release from a high pressure shock state to a pressure-density path is analyzed. It is assumed that the release process is isentropic. It was shown that for geological materials below stresses of 150 GPa, the effective viscosity must be 1000 kg/m/s in order that the viscous (irreversible) work carried out on the material in the shock state remains small compared to the mechanical work recovered upon adiabatic rarefaction. The available data pertaining to the offset of the Rayleigh line from the Hugoniot for minerals, the magnitude of the shear stress in the high pressure shock state for minerals, and the direct measurements of the viscosities of several engineering materials shocked to pressures below 150 GPa yield effective viscosities of 1000 kg/m/s or less. An inferance that this indicates that the conditions for isentropic release of minerals from shock states are achieved, and a conclusion that the application of the Riemann integral to obtain pressure-density states along the release adiabats of minerals in shock experiments is valid are made.

Jeanloz, R.↗
DOE OSTI2023

Lithiation Gradients and Tortuosity Factors in Thick NMC111-Argyrodite Solid-State Cathodes

Achieving high energy density in all-solid-state lithium batteries will require the design of thick cathodes, and these will need to operate reversibly under normal use conditions. We use high-energy depth-profiling X-ray diffraction to measure the localized lithium content of Li 1-x Ni 1/3 Mn 1/3 Co 1/3 O 2 (NMC111) through the thickness of 110 μm thick composite cathodes. The composite cathodes consisted of NMC111 of varying mass loadings mixed with argyrodite solid electrolyte Li 6 PS 5 Cl (LPSC). During cycling at C/10, substantial lithiation gradients developed, and varying the NMC111 loading altered the nature of these gradients. Microstructural analysis and cathode modeling showed this was due to high tortuosities in the cathodes. This was particularly true in the solid electrolyte phase, which experienced a marked increase in tortuosity factor during the initial charge. Our results demonstrate that current distributions are observed in sulfide-based composites and that these will be an important consideration for practical design of all-solid-state batteries.

25 ENERGY STORAGE↗
NASA NTRS1980

Improved characterization of the Si-SiO2 interface

Refined quasi-static and conductance methods, based on effectively thin composite insulating layers, low-carrier-concentration bulk semiconductors, and low-level illumination, have been applied to an improved characterization of the (100) Si-SiO2 interface. Accurate measurement of both the total density of interface states and its major components as a function of energy in the forbidden gap have been made over four decades (10-billion to 100-trillion states/eV sq cm) on a single sample. The normal U-shaped density of states is resolved into separate valence- and conduction-band-derived contributions as well as impurity-derived contributions corresponding to concentrations on the order of 20 ppm at the interface.

Su, P.↗
DOE OSTI2022

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

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

Spectroscopy and electronic structure of the low-energy states of ThN

The electronic spectrum of ThN over the range 19,600–21,200 cm -1 is remarkably congested, showing at least twenty vibronic bands originating from the ground state zero-point level. Rotationally resolved spectra for eleven of these bands have been examined to probe the underlying reasons for the dense manifold of states. Dispersed fluorescence spectra and fluorescence decay lifetimes were measured to provide addition insights. Our experimental measurements were complemented by electronic structure calculations for the low-energy doublet and quartet states. The ab initio calculations yielded a density of excited states that was consistent with the congested spectrum. These states were derived from the formal configurations Th 3+ (6d)N 3- , Th 2+ (7s 2 )N 2- (2p 5 ) and Th 2+ (7s6d)N 2 -(2p 5 ). The calculations indicated that the rotational and vibrational constants would be characteristic of the electronic configurations. This prediction was not borne out by the experimental data, suggesting that the full range of inter-configurational state interactions were not captured by the computational model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The effect of aspect ratio on the mechanical behavior of Li metal in solid-state cells

The Li metal anode is necessary for solid-state batteries to dramatically exceed the energy density of state-of-the-art conventional Li-ion batteries. To enable the Li metal anode, the mechanical properties of Li metal under compression in a commercially relevant cell configuration must be determined. This work examines the compression creep behavior of Li metal coupled with a relevant model solid-state electrolyte (Li 6.5 La 3 Ta 0.5 Zr 1.5 O 12 ) and current collector (Ni) as a function of applied stress (1–10 MPa) and Li foil thickness (100–700 μm). At 1 MPa applied stress, the strain rate of all Li foil thicknesses under compression was almost four orders of magnitude lower than what is predicted from tensile creep tests. Steady-state creep behavior was observed for Li foil ≥350 μm at 5 MPa and ≥250 μm at 10 MPa. Conversely, sustained creep deformation was not observed for 100 and 150 μm Li foils at any applied stress. It is suggested that friction at the Li/LLZO and Li/Ni interfaces leads to the formation of hydrostatic stress that controls Li deformation at low aspect ratios. Finally, these results have implications for battery pack stack pressure requirements, particularly relating to void formation during stripping and Li metal filament penetration during plating.

25 ENERGY STORAGE↗