Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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 433 records · Page 24

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermoelectric Figure of Merit of a Superatomic Crystal Re 6 Se 8 I 2 Monolayer

Superatomic materials are newly emerging candidates for high-performance thermoelectric (TE) devices due to their intrinsic ultralow thermal conductivities. However, the low TE power factor becomes a huge obstacle to reaching the required the dimensionless figure of merit (ZT) values for practical applications. Here, motivated by the recently synthesized superatomic Re 6 Se 8 I 2 monolayer [J. Am. Chem. Soc. 144, 74 (2022)], we study its superior thermoelectric properties by using density functional theory combined with phonon Boltzmann transport theory and deep potential molecular dynamics. We show that the large mass and anharmonic Re-I bonds introduce strong phonon scattering and result in low lattice thermal conductivity of 1.20 W m -1 K -1 at 300 K, while the strong and harmonic Re-Se network ensures the high TE power factor of 4344 μW m -1 K -2 in b direction for n-type doping. This is an order of magnitude higher than that of other cluster-based materials. Due to the low thermal conductivity and high TE power factor, Re 6 Se 8 I 2 exhibits high ZT values of 1.20 (500 K) and 1.43 (900 K) with n-type doping along the b direction among all the cluster-based thermoelectric materials reported so far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phonons, Q-dependent Kondo spin fluctuations, and 4f phonon resonance in YbAl 3

4f intermediate valence (IV) compounds are canonical hosts of correlated electron physics and can contribute to our understanding of the larger class of correlated electron materials. In this work, we study the prototype IV compound YbAl 3 which exhibits a nonintegral valence with a moderately heavy fermion ground state and a large Kondo temperature (T K ~500–600K). To better characterize the correlated physics of YbAl 3 , we have measured the phonon and the magnetic excitation spectra on single crystals of this material by time-of-flight inelastic neutron scattering and inelastic x-ray scattering. We have also performed theoretical calculations of the phonon spectra. We present three findings of these measurements. First, we observe that the measured phonon spectra can be described adequately by a calculation based on standard DFT+U density functional theory. The calculated energies, however, are 10% too low compared to the measured energies. This discrepancy may reflect a hardening of the phonons due to dynamic 4f correlations. Second, the low-temperature spin fluctuations on the Kondo energy scale k B T K have a momentum (Q) dependence similar to that seen recently in the IV compound CePd 3 . For that system, the Q dependence has been attributed to particle-hole excitations in a coherent itinerant 4f correlated ground state. We suggest a similar origin for the momentum dependence seen in YbAl 3 . This Q dependence disappears as the temperature is raised towards room temperature and the 4f electron band states become increasingly incoherent. Such a coherent/incoherent crossover is expected to be generic for correlated electron systems. Third, a low-temperature magnetic peak observed in the neutron scattering near 30 meV shows dispersion identical to a particular optic-phonon branch. This 4f/phonon resonance disappears for T≥150K. The phonon spectrum appears to be unaffected by the resonance. We discuss several possibilities for the origin of this unusual excitation, which may be unique to YbAl 3 . We suggest that the excitation may arise from the large amplitude beating of the light Al atoms against the heavy Yb atoms, resulting in a dynamic 4f/3p hybridization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic states dressed by an out-of-plane supermodulation in the quasi-two-dimensional kagome superconductor CsV 3 Sb 5

CsV 3 Sb 5 has attracted much recent attention as the first quasi-two-dimensional (2D) kagome superconductor. While the kagome layers are 2D in nature, increasing evidence has pointed to the importance of out-of-plane correlation in this material. However, it remains unclear whether such correlation can change the fundamental electronic structure of the quasi-2D system. Here, we reveal this missing piece of information, using angle-resolved photoemission spectroscopy, complemented by scanning tunneling microscope measurements. The three-dimensional electronic structures in the high-temperature state are revealed, which agree well with density-functional theory calculations. Electron energy bands are observed in the low-temperature state that exhibit additional periodicities along the out-of-plane momentum. Furthermore, these results reveal a direct response to the out-of-plane electronic supermodulation in the single-particle spectral function of CsV 3 Sb 5 , thus establishing an electronic platform to examine emergent phenomena beyond 2D limit in kagome superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Excitonic Condensate in Flat Valence and Conduction Bands of Opposite Chirality

Excitonic Bose-Einstein condensation (EBEC) has drawn increasing attention recently with the emergence of 2D materials. A general criterion for EBEC, as expected in an excitonic insulator (EI) state, is to have negative exciton formation energies in a semiconductor. Here, using exact diagonalization of a multiexciton Hamiltonian modeled in a diatomic kagome lattice, we demonstrate that the negative exciton formation energies are only a prerequisite but insufficient condition for realizing an EI. By a comparative study between the cases of both conduction and valence flat bands (FBs) versus that of a parabolic conduction band, we further show that the presence and increased FB contribution to exciton formation provide an attractive avenue to stabilize the excitonic condensate, as confirmed by calculations and analyses of multiexciton energies, wave functions, and reduced density matrices. In conclusion, our results warrant a similar many-exciton analysis for other known and/or new candidates of EIs and demonstrate the FBs of opposite parity as a unique platform for studying exciton physics, paving the way to material realization of spinor BEC and spin superfluidity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Initializing BSQ with Open-Source ICCING

While it is well known that there is a significant amount of conserved charges in the initial state of nuclear collisions, the production of these due to gluon splitting has yet to be thoroughly investigated. The ICCING (Initial Conserved Charges in Nuclear Geometry) algorithm reconstructs these quark distributions, providing conserved strange, baryon, and electric charges, by sampling a given model for the g → qq¯ splitting function over the initial energy density, which is valid at top collider energies, even when µB = 0. The ICCING algorithm includes fluctuations in the gluon longitudinal momenta, a structure that supports the implementation of dynamical processes, and the c++ version is now open-source. A full analysis of parameter choices on the model has been done to quantify the effect these have on the underlying physics. We find there is a sustained difference across the different charges that indicates sensitivity to hot spot geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ion beam sputter etching and deposition of fluoropolymers

Fluoropolymer etching and deposition techniques including thermal evaporation, RF sputtering, plasma polymerization, and ion beam sputtering are reviewed. Etching and deposition mechanism and material characteristics are discussed. Ion beam sputter etch rates for polytetrafluoroethylene (PTFE) were determined as a function of ion energy, current density and ion beam power density. Peel strengths were measured for epoxy bonds to various ion beam sputtered fluoropolymers. Coefficients of static and dynamic friction were measured for fluoropolymers deposited from ion bombarded PTFE.

Banks, B. A.↗

Ion beam sputter etching and deposition of fluoropolymers

Fluoropolymer etching and deposition techniques including thermal evaporation, RF sputtering, plasma polymerization, and ion beam sputtering are reviewed. Etching and deposition mechanisms and material characteristics are discussed. Ion beam sputter etch rates for polytetrafluoroethylene (PTFE) were determined as a function of ion energy, current density and ion beam power density. Peel strengths were measured for epoxy bonds to various ion beam sputtered fluoropolymers. Coefficients of static and dynamic friction were measured for fluoropolymers deposited from ion bombarded PTFE.

Banks, B. A.↗

Measurements by the plasma diagnostics package on STS-3

A comprehensive set of measurements about the orbiter environment are provided by the plasma diagnostics package (PDP). Ion and electron particle densities, energies, and spatial distribution functions; ion mass for identification of particular molecular ion species; and magnetic fields, electric fields and electromagnetic waves over a broad frequency range are studied. Shuttle environmental measurements will be made both on the pallet and, by use of the remote manipulator system (RMS), the PDP will be maneuvered in and external to the bay area to continue environmental measurements and to carry on a joint plasma experiment with the Utah State University fast-pulsed electron generator. Results of orbiter environment EMI measurements and S-band field strengths as well as preliminary results from wake search operations indicating wake boundary identifiers are reported.

Shawhan, S. D.↗

A Model for Lower Hybrid Wave Excitation Compared with Observations by Viking

The mechanism of lower hybrid wave (LHW) excitation due to the O+ relative drift in a plasma subjected to low-frequency waves (LFWs) is used for analysis of Viking satellite data for events in the cusp/cleft region. In some cases, such a mechanism leads to LHW energy densities and ion distribution functions close to those observed, suggesting the proposed mechanism is a plausible candidate to explain certain classes of LHW generation events in space plasmas.

Khazanov, G. V.↗

Electrons In The Low Density Solar Wind

The recent occurrence of an interval (May 9th to May 12th, 1999) of abnormally low density solar wind has drawn attention to such events. The SWE instrument on the Wind spacecraft observed nine similar events between launch (November 1994) and August 1999: one in 1997, three in 1998, and five in January-August 1999. No such events were observed in 1996, the year of solar minimum. This already suggests a strong dependence upon solar activity. In this paper we discuss observations of the electron strahl, a strong anisotropy in the solar wind electrons above 60 eV directed along the magnetic field and observed continuously during the periods of low density in 1998 and 1999. When the solar wind density was less than 2/cc, the angular width of the strahl was below 3.5 degrees and the temperature deduced from the slope of the electron strahl phase density (as a function of energy in the energy range 200 to 800 eV) was 100 to 150 eV, equivalent to a typical coronal electron temperature. Three examples of this phenomenon, observed on Feb. 20- 22, April 26-27 and May 9-12, 1999, are discussed to show their similarity to one another. These electron observations are interpreted to show that the strahl occurs as a result of the conservation of the first adiabatic invariant, combined with the lack of coulomb collisions as suggested by Fairfield and Scudder, 1985.

Ogilvie, Keith W.↗

Composite Cathodes for Dual-Rate Li-Ion Batteries

Composite-material cathodes that enable Li-ion electrochemical cells and batteries to function at both high energy densities and high discharge rates are undergoing development. Until now, using commercially available cathode materials, it has been possible to construct cells that have either capability for high-rate discharge or capability to store energy at average or high density, but not both capabilities. However, both capabilities are needed in robotic, standby-power, and other applications that involve duty cycles that include long-duration, low-power portions and short-duration, high-power portions. The electrochemically active ingredients of the present developmental composite cathode materials are: carbon-coated LiFePO4, which has a specific charge capacity of about 160 mA h/g and has been used as a high-discharge-rate cathode material and Li[Li(0.17)Mn(0.58)Ni(0.25)]O2, which has a specific charge capacity of about 240 mA h/g and has been used as a high-energy-density cathode material. In preparation for fabricating the composite material cathode described, these electrochemically active ingredients are incorporated into two sub-composites: a mixture comprising 10 weight percent of poly(vinylidine fluoride); 10 weight percent of carbon and 80 weight percent of carbon coated LiFePO4; and, a mixture comprising 10 weight percent of PVDF, and 80 weight percent of Li[Li(0.17)Mn(0.58)Ni(0.25)]O2. In the fabrication process, these mixtures are spray-deposited onto an aluminum current collector. Electrochemical tests performed thus far have shown that better charge/discharge performance is obtained when either 1) each mixture is sprayed on a separate area of the current collector or (2) the mixtures are deposited sequentially (in contradistinction to simultaneously) on the same current-collector area so that the resulting composite cathode material consists of two different sub-composite layers.

Whitacre, Jay↗

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗

The Distribution of Lightning Channel Lengths in Northern Alabama Thunderstorms

Lightning is well known to be a major source of tropospheric NOx, and in most cases is the dominant natural source (Huntreiser et al 1998, Jourdain and Hauglustaine 2001). Production of NOx by a segment of a lightning channel is a function of channel segment energy density and channel segment altitude. A first estimate of NOx production by a lightning flash can be found by multiplying production per segment [typically 104 J/m; Hill (1979)] by the total length of the flash s channel. The purpose of this study is to determine average channel length for lightning flashes near NALMA in 2008, and to compare average channel length of ground flashes to the average channel length of cloud flashes.

Peterson, H. S.↗

Novel Catalysts and Processing Technologies for Production of Aerospace Fuels from Non-Petroleum Raw Materials

Transportation fuels production (including aerospace propellants) from non-traditional sources (gases, waste materials, and biomass) has been an active area of research and development for decades. Reducing terrestrial waste streams simultaneous with energy conversion, plentiful biomass, new low-cost methane sources, and/or extra-terrestrial resource harvesting and utilization present significant technological and business opportunities being realized by a new generation of visionary entrepreneurs. We examine several new approaches to catalyst fabrication and new processing technologies to enable utilization of these nontraditional raw materials. Two basic processing architectures are considered: a single-stage pyrolysis approach that seeks to basically re-cycle hydrocarbons with minimal net chemistry or a two-step paradigm that involves production of supply or synthesis gas (mainly carbon oxides and H2) followed by production of fuel(s) via Sabatier or methanation reactions and/or Fischer-Tröpsch synthesis. Optimizing the fraction of product stream relevant to targeted aerospace (and other transportation) fuels via modeling, catalyst fabrication and novel reactor design are described. Energy utilization is a concern for production of fuels for either terrestrial or space operations; renewable sources based on solar energy and/or energy efficient processes may be mission enabling. Another important issue is minimizing impurities in the product stream(s), especially those potentially posing risks to personnel or operations through (catalyst) poisoning or (equipment) damage. Technologies being developed to remove (and/or recycle) heteroatom impurities are briefly discussed as well as the development of chemically robust catalysts whose activities are not diminished during operation. The potential impacts on future missions by such new approaches as well as balance of system issues are addressed.

space processing↗

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Silicon-based ceramic components for next-generation jet turbine engines offer potential weight savings, as well as higher operating temperatures, both of which lead to increased efficiency and lower fuel costs. Silicon carbide (SiC), in particular, offers low density, good strength at high temperatures, and good oxidation resistance in dry air. However, reaction of SiC with high-temperature water vapor, as found in the hot section of jet turbine engines in operation, can cause rapid surface recession, which limits the lifetime of such components. Environmental Barrier Coatings (EBCs) are therefore needed if long component lifetime is to be achieved. Rare earth silicates such as Yb2Si2O7 and Yb2SiO5 have been proposed for such applications; in an effort to better understand diffusion in such materials, we have performed kinetic Monte Carlo (kMC) simulations of oxygen diffusion in Ytterbium disilicate, Yb2- Si2O7. The diffusive process is assumed to take place via the thermally activated hopping of oxygen atoms among oxygen vacancy sites or among interstitial sites. Migration barrier energies are computed using density functional theory (DFT).

Computer Simulation↗