Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local density approximations”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Galactic Potential and Dark Matter Density from Angular Stellar Accelerations

Here, we present an approach to measure the Milky Way (MW) potential using the angular accelerations of stars in aggregate as measured by astrometric surveys like Gaia. Accelerations directly probe the gradient of the MW potential, as opposed to indirect methods using, e.g., stellar velocities. We show that end-of-mission Gaia stellar acceleration data may be used to measure the potential of the MW disk at approximately 3σ significance and, if recent measurements of the solar acceleration are included, the local dark matter density at ~2σ significance. Since the significance of detection scales steeply as t 5/2 for observing time t, future surveys that include angular accelerations in the astrometric solutions may be combined with Gaia to precisely measure the local dark matter density and shape of the density profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Orientation-dependent structural properties during growth and growth mechanism of CoO films

The orientation-dependent local structural properties of CoO films on sapphire substrates during growth were investigated through linearly-polarized extended X-ray absorption fine structure (EXAFS) measurements. Specifically, CoO(1 1 1) and (100) crystals with a rock-salt structure (Fm3m) were epitaxially grown on alpha-Al2O3(0001) and (1012) substrates, respectively, at 700 degrees C using a radio-frequency sputtering system. The local structural properties of CoO films in the in-plane and out-of-plane orientations were quantitatively determined using linearly-polarized EXAFS at the Co K-edge during growth. The EXAFS analysis revealed that during the initial stages of growth, the local structural properties exhibit significant differences compared to thick films, with short atomic distances and large (small) Debye-Waller factors observed in the out-of-plane (in-plane) orientations. The local structural strain mostly diminished when approximately 20 CoO layers accumulated on the substrate. Density functional theory (DFT) calculations further supported these findings by confirming that cobalt atoms initially form stable bonds with the sapphire surface, leading to the simultaneous growth of oxygen and cobalt layers in a coordinated manner through layer-by-layer growth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Relative importance of nonlinear electron-phonon coupling and vertex corrections in the Holstein model

Determining the range of validity of Migdal’s approximation for electron-phonon (e-ph) coupled systems is a long-standing problem. Many attempts to answer this question employ the Holstein Hamiltonian, where the electron density couples linearly to local lattice displacements. When these displacements are large, however, nonlinear corrections to the interaction must also be included, which can significantly alter the physical picture obtained from this model. Using determinant quantum Monte Carlo and the self-consistent Migdal approximation, we compared superconducting and charge-density-wave correlations in the Holstein model with and without second-order nonlinear interactions. We find a disagreement between the two cases, even for relatively small values of the e-ph coupling strength, and, importantly, that this can occur in the same parameter regions where Migdal’s approximation holds. Our results demonstrate that questions regarding the validity of Migdal’s approximation go hand in hand with questions of the validity of a linear e-ph interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Anderson Localization in the Fractional Quantum Hall Effect

The interplay between interaction and disorder-induced localization is of fundamental interest. This article addresses localization physics in the fractional quantum Hall state, where both interaction and disorder have nonperturbative consequences. We provide compelling theoretical evidence that the localization of a single quasiparticle of the fractional quantum Hall state at filling factor ν=n/(2n+1) has a striking quantitative correspondence to the localization of a single electron in the (n+1)th Landau level. By analogy to the dramatic experimental manifestations of Anderson localization in integer quantum Hall effect, this leads to predictions in the fractional quantum Hall regime regarding the existence of extended states at a critical energy, and the nature of the divergence of the localization length as this energy is approached. Within a mean field approximation, these results can be extended to situations where a finite density of quasiparticles is present.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Application of the locally self-consistent embedding approach to the Anderson model with non-uniform random distributions

Highlights: • Typical Medium Theory (TMT) for the Anderson Localization. • Locally Self-Consistent Multiple Scattering Method (LSMS) for Random Disordered Systems. • Linear Scaling Computational Method for Random Disordered Systems. We apply the recently developed embedding scheme for the locally self-consistent method to random disorder electrons systems. The method is based on the locally self-consistent multiple scattering theory and the typical medium theory. The locally self-consistent multiple scattering theory divides a system into many small designated local interaction zones. The subsystem within each local interaction zone is embedded in a self-consistent field from the typical medium theory. This approximation allows the study of random systems with large numbers of sites. We present results for the three dimensional Anderson model with different random disorder potential distributions. Using the typical density of states as an indicator of Anderson localization, we find that the method can capture the localization for commonly studied disorder potentials. These include the uniform distribution, the Gaussian distribution, and even the unbounded Cauchy distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Static dipole polarizabilities of polyacenes using self-interaction-corrected density functional approximations

Density functional approximations are known to significantly overestimate the polarizabilities of long chain-like molecules. We study the static electric dipole polarizabilities and the vertical ionization potentials of polyacenes from benzene to pentacene using the Fermi-L\"owdin orbital based self-interaction corrected (FLOSIC) density functional method. The orbital-by-orbital self-interaction correction corrects for the overestimation tendency of density functional approximations. The polarizabilities calculated with FLOSIC-DFA are however overly corrected. We also tested the recently developed locally-scaled self-interaction correction (LSIC) method on the polyacenes. The local-scaling method applies full SIC in the one-electron regions and restores the proper behavior of the SIC exchange-correlation functionals in the uniform density limit. The results show that LSIC removes the overcorrection tendency of the FLOSIC-DFA and produces results that are in excellent agreement with reference CCSD values. The vertical ionization potentials with LSIC also show good agreement with available experimental values.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Projections of H-mode access and edge pedestal in the SPARC tokamak

In order to inform core performance projections and divertor design, the baseline SPARC tokamak plasma discharge is evaluated for its expected H-mode access, pedestal pressure and edge-localized mode (ELM) characteristics. A clear window for H-mode access is predicted for full field DT plasmas, with the available 25 MW of design auxiliary power. Additional alpha heating is likely needed for H-mode sustainment. Pressure pedestal predictions in the developed H-mode are surveyed using the EPED model. The projected SPARC pedestal would be limited dominantly by peeling modes and may achieve pressures in excess of 0.3 MPa at a density of approximately 3 × 10 20 m -3 . High pedestal pressure is partially enabled by strong equilibrium shaping, which has been increased as part of recent design iterations. Edge-localized modes (ELMs) with >1 MJ of energy are projected, and approaches for reducing the ELM size, and thus the peak energy fluence to divertor surfaces, are under consideration. The high pedestal predicted for SPARC provides ample margin to satisfy its high fusion gain (Q) mission, so that even if ELM mitigation techniques result in a 2x reduction of the pedestal pressure, Q > 2 is still predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Entanglement and correlations in fast collective neutrino flavor oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales t ∝ μ –1 , with μ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number N of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales τ ≈ µ –1 $\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales τ ≈ μ –1 log(N). Here, we present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Estimating Electron Temperature and Density Using Van Allen Probe Data: Typical Behavior of Energetic Electrons in the Inner Magnetosphere

Abstract The Earth's inner magnetosphere contains multiple electron populations influenced by different factors. The cold electrons of the plasmasphere, warm plasma that contributes to the ring current, and the relativistic plasma of the radiation belts often seem to behave independently. Using omni‐directional flux and energy measurements from the HOPE and Magnetic Electron Ion Spectrometer instruments aboard the Van Allen Probes, we provide a detailed density and temperature description of the inner magnetosphere, offering a comprehensive statistical analysis of the entire Van Allen Probe era. While number density and temperature data at geosynchronous orbit are available, this study focuses on the warm plasma in the inner magnetosphere . Values of density and temperature are extracted by fitting energy and phase space density to obtain the distribution function. The fitted distributions are related to the zeroth and second moments to estimate the number density and temperature. Analysis has indicated that a two Maxwellian fit is sufficient over a wide range of and that there are two independent plasma populations. The more energetic population has a median number density of approximately and a temperature of around 130 keV, with a temperature peak observed between L * = 4 and L * = 4.5. This population is relatively uniform in magnetic local time (MLT). In contrast, the less energetic warm electron population has a median number density of about and a temperature of 7.4 keV. Strong statistical trends in density and temperature across both L * and MLT are presented, along with potential sources driving these variations.

58 GEOSCIENCES↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An adaptive Hessian approximated stochastic gradient MCMC method

Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their ability to handle large datasets and the potential to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be inefficient when the random variables in the target posterior densities have scale differences or are highly correlated. Here, we present an adaptive Hessian approximated stochastic gradient MCMC method to incorporate local geometric information while sampling from the posterior. The idea is to apply stochastic approximation (SA) to sequentially update a preconditioning matrix at each iteration. The preconditioner possesses second-order information and can guide the random walk of a sampler efficiently. Instead of computing and saving the full Hessian of the log posterior, we use limited memory of the samples and their stochastic gradients to approximate the inverse Hessian-vector multiplication in the updating formula. Moreover, by smoothly optimizing the preconditioning matrix via SA, our proposed algorithm can asymptotically converge to the target distribution with a controllable bias under mild conditions. To reduce the training and testing computational burden, we adopt a magnitude-based weight pruning method to enforce the sparsity of the network. Our method is user-friendly and demonstrates better learning results compared to standard SG-MCMC updating rules. The approximation of inverse Hessian alleviates storage and computational complexities for large dimensional models. Numerical experiments are performed on several problems, including sampling from 2D correlated distribution, synthetic regression problems, and learning the numerical solutions of heterogeneous elliptic PDE. The numerical results demonstrate great improvement in both the convergence rate and accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Plasmons in two-dimensional topological insulators

Here, we analyze collective excitations in models of two-dimensional topological insulators using the random phase approximation. In a two-dimensional extension of the Su-Schrieffer-Heeger model, edge plasmonic excitations with induced charge-density distributions localized at the boundaries of the system are found in the topologically nontrivial phase, dispersing similarly as one-dimensional bulk plasmons in the conventional Su-Schrieffer-Heeger chain. For two-dimensional bulk collective modes, we reveal regimes of enhanced interband wave function correlations, leading to characteristic hardening and softening of inter- and intraband bulk plasmonic branches, respectively. In the two-dimensional Haldane Chern insulator model, chiral, unidirectional edge plasmons in nano-ribbon architectures are observed, which can be characterized by an effective Coulomb interaction cross section. Bulk collective excitations in the two-dimensional Haldane model are shown to be originated by single-particle band structure details in different topological phases.

2-dimensional systems↗
DOE OSTI2020

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Development of local hybrid density functionals to treat self-interaction error and many-electron effects

The goal of this project was to improve density functional theory (DFT) for systems where conventional semilocal approximations are limited by self-interaction error (SIE) and by near-degeneracy or strong many-electron effects. While DFT remains the only broadly practical first-principles framework for large-scale materials simulations, its predictive accuracy is often challenged in situations involving stretched bonds, charge transfer, transition-metal chemistry, magnetic couplings, band gaps, and correlated electronic states. This project addressed these limitations by developing physically grounded and numerically robust exchange–correlation functionals that retain the efficiency of modern DFT while extending its predictive scope.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Implications from Monitoring Gopher Tortoises at Two Spatial Scales

A problem that conservation biologists face is how to monitor species, given both resource limitations and the inherent challenges of assessing long-term demographic processes. We assessed gopher tortoise (Gopherus polyphemus) abundance at a landscape scale and at the scale of 3 local populations within the Conecuh National Forest (CNF).USA between 1991 and 2017. Landscape-level data were collected 26 from line transect distance sampling arranged uniformly across the CNF and collected during a single season (2011); data for local populations were generated from long-term mark-recapture of individuals at three sites selected based on prior knowledge of high density at each. At a landscape scale, we estimated 5,242 (95% CI = 3,538-7,768) tortoises occurred across the approximately 34,000-ha forest, yielding a density of 0.14-0.32 tortoises/ ha. These low densities across the landscape suggest that, on average, management activities across the property have not allowed tortoise populations to retain social structure needed for long-term persistence. The three local populations, however, contained 25-60 individuals and densities of 1.9–6.9 tortoises/ha. Over the study period, populations at two sites were stable and the third experienced significant population growth. Mean annual survival of individuals was 0.89 and invariant across size classes. Altogether, line transect distance sampling is important for assessing landscape-scale abundance of tortoises but may fail to detect local clusters of high-density sites important for population persistence. Our mark-recapture efforts at the local scale revealed that small populations on these high-density sites can exhibit long-term stability or growth even though they do not meet current established criteria for viability. Improved models that incorporate immigration and emigration and better reflect dynamics of peripheral populations would assist in determining how such populations best contribute to species recovery and regional conservation targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronic structure of mononuclear Cu-based molecule from density-functional theory with self-interaction correction

In this paper, we investigate the electronic structure of a planar mononuclear Cu-based molecule [Cu(C 6 H 4 S 2 ) 2 ] z in two oxidation states ( z = –2, –1) using density-functional theory (DFT) with Fermi–Löwdin orbital (FLO) self-interaction correction (SIC). The dianionic Cu-based molecule was proposed to be a promising qubit candidate. Self-interaction error within approximate DFT functionals renders severe delocalization of electron and spin densities arising from 3 d orbitals. The FLO-SIC method relies on optimization of Fermi–Löwdin orbital descriptors (FODs) with which localized occupied orbitals are constructed to create SIC potentials. Starting with many initial sets of FODs, we employ a frozen-density loop algorithm within the FLO-SIC method to study the Cu-based molecule. We find that the electronic structure of the molecule remains unchanged despite somewhat different final FOD configurations. In the dianionic state (spin S = 1/2), FLO-SIC spin density originates from the Cu d and S p orbitals with an approximate ratio of 2:1, in quantitative agreement with multireference calculations, while in the case of SIC-free DFT, the orbital ratio is reversed. Overall, FLO-SIC lowers the energies of the occupied orbitals and, in particular, the 3 d orbitals unhybridized with the ligands significantly, which substantially increases the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) compared to SIC-free DFT results. The FLO-SIC HOMO–LUMO gap of the dianionic state is larger than that of the monoanionic state, which is consistent with experiment. Our results suggest a positive outlook of the FLO-SIC method in the description of magnetic exchange coupling within 3 d -element-based systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗