Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Noniterative finite amplitude methods for giant resonances and the application to the neutron radiative capture cross sections

We calculate the electric dipole (E1) and the magnetic dipole (M1) giant resonances with noniterative finite amplitude methods and demonstrate how the fully microscopic density functional theory predicts the giant resonances without any phenomenological parameters. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the result of E1 and M1 transitions and find that the capture cross sections for deformed nuclei are enhanced due to the contribution from the low energy M1 scissors mode.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Correlation-assisted quantized charge pumping

Here, we investigate charge pumping in the vicinity of order-obstructed topological phases, i.e., symmetry-protected topological phases masked by spontaneous symmetry breaking in the presence of strong correlations. To explore this, we study a prototypical Su-Schrieffer-Heeger model with finite-range interaction that gives rise to orbital charge density wave order and characterize the impact of this order on the model's topological properties. In the ordered phase, where the many-body topological invariant loses quantization, we find that not only is quantized charge pumping still possible, but it is even assisted by the collective nature of the orbital charge density wave order. Remarkably, we show that the Thouless pump scenario may be used to uncover the underlying topology of order-obstructed phases.

1-dimensional systems↗
DOE OSTI2024

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S = 1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule. Here, we anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Probing the Neutron Skin with Ultrarelativistic Isobaric Collisions

Particle production in ultrarelativistic heavy ion collisions depends on the details of the nucleon density distributions in the colliding nuclei. We demonstrate that the charged hadron multiplicity distributions in isobaric collisions at ultrarelativistic energies provide a novel approach to determine the poorly known neutron density distributions and thus the neutron skin thickness in finite nuclei, which can in turn put stringent constraints on the nuclear symmetry energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Stability, metallicity, and magnetism in niobium silicide nanofilms

Modern superconducting qubits based on two-dimensional (2D) transmons typically involve the growth of Nb thin films on high-resistivity Si substrates. Since imperfections at the Nb-Si heterointerface have been implicated as a source of two-level systems that limit quantum coherence times, detailed characterization and understanding of niobium silicide interfacial layers are critical to improving superconducting qubit technology. While bulk binary intermetallic niobium silicide phases are well understood, the thermodynamic phase stability and properties of ultrathin niobium silicides, such as those found at the Nb-Si heterointerface in 2D transmons, have not yet been explored. Here, we report finite-sized effects for ultrathin niobium silicide films using density functional theory calculations and predict nanoscale stabilization of Nb 6 Si 5 over the bulk α-Nb 5 Si 3 phase. This result is consistent with our experimental observations of a niobium silicide interfacial layer between a sputtered Nb thin film and the underlying Si substrate. Furthermore, our calculations show that Nb 6 Si 5 nanofilms are nonmagnetic, making them superior to nanofilms of α-Nb 5 Si 3 that exhibit antiferromagnetic correlations detrimental to long coherence times in superconducting qubits. Furthermore, by providing atomic-scale insight into niobium silicide nanofilms, this paper can help guide ongoing efforts to optimize Nb-Si heterointerfaces for long coherence times in superconducting qubits.

36 MATERIALS SCIENCE↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Comprehensive Studies on Hot Compaction and Vibration-Assisted Compaction Tests of Aluminum Powder

Abstract Aluminum powder compaction was studied using both test and simulation. Cold compaction, hot compaction, and vibration-assisted (cold) compaction tests were conducted to achieve different density ratios. First, the hot compaction test (at 300 °C, compression pressure 140 MPa) improved about 6% compared with cold compaction under the same compression pressure. Second, although the relative density ratio does not obviously improve at a vibration-assisted (cold) compaction, the strength of the specimens made under vibration loading is much better than those of cold compaction. Additionally, finite element models with well-calibrated Drucker–Prager Cap (DPC) material constitutive model were built in abaqus/standard to simulate the powder compaction process. The results of the finite element model have very good correlations with test results up to the tested range, and this finite element model further predicts the loading conditions needed to achieve the higher density ratios. Two exponential equations of the predicted density ratio were obtained by combining the test data and the simulation results. A new analytical solution was developed to predict the axial pressure versus the density ratio for the powder compaction according to DPC material model. The results between the analytical solution and the simulation model have a very good match.

Engineering↗
DOE OSTI2022

Nonequilibrium spintronic transport through Kondo impurities

In this work we analyze the nonequilibrium transport through a quantum impurity (quantum dot or molecule) attached to ferromagnetic leads by using a hybrid numerical renormalization group–time-dependent density matrix renormalization group thermofield quench approach. For this, we study the bias dependence of the differential conductance through the system, which shows a finite zero-bias peak, characteristic of the Kondo resonance and reminiscent of the equilibrium local density of states. In the nonequilibrium settings, the resonance in the differential conductance is also found to decrease with increasing the lead spin polarization. The latter induces an effective exchange field that lifts the spin degeneracy of the dot level. Therefore, as we demonstrate, the Kondo resonance can be restored by counteracting the exchange field with a finite external magnetic field applied to the system. Lastly, we investigate the influence of temperature on the nonequilibrium conductance, focusing on the split Kondo resonance. Our work thus provides an accurate quantitative description of the spin-resolved transport properties relevant for quantum dots and molecules embedded in magnetic tunnel junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Kinetically constrained freezing transition in a dipole-conserving system

Here, we study a stochastic lattice gas of particles in one dimension with strictly finite-range interactions that respect the fractonlike conservation laws of total charge and dipole moment. As the charge density is varied, the connectivity of the system's charge configurations under the dynamics changes qualitatively. We find two distinct phases: Near half filling the system thermalizes subdiffusively, with almost all configurations belonging to a single dynamically connected sector. As the charge density is tuned away from half filling there is a phase transition to a frozen phase, where locally active finite bubbles cannot exchange particles and the system fails to thermalize. The two phases exemplify what has recently been referred to as weak and strong Hilbert space fragmentation, respectively. We study the static and dynamic scaling properties of this weak-to-strong fragmentation phase transition in a kinetically constrained classical Markov circuit model, obtaining some conjectured exact critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Microscopic Theory of Pair Density Waves in Spin-Orbit Coupled Kondo Lattice

We demonstrate that the discommensuration between the Fermi surfaces of a conduction sea and an underlying spin liquid provides a natural mechanism for the spontaneous formation of pair density waves. Using a recent formulation of the Kondo lattice model that incorporates a Yao Lee spin liquid proposed by the authors, we demonstrate that doping away from half filling induces finite-momentum electron-Majorana pair condensation, resulting in amplitude-modulated pair density waves (PDWs). Our approach provides a precise, analytically tractable pathway for understanding the spontaneous formation of PDWs in higher dimensions and offers a natural mechanism for PDW formation in the absence of Zeeman splitting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Global LIS/OTD Climatology of Lightning Flash Extent Density

Previous lightning climatologies derived from Lightning Imaging Sensor (LIS) and Optical Transient Detector (OTD) total lightning measurements have quantified lightning frequency as a Flash Rate Density (FRD). This approach assumes that lightning flashes can be represented as points and quantifies the frequency of lightning centered in each grid cell. However, lightning has a finite extent that can reach hundreds of kilometers. A new climatology based on Flash Extent Density (FED) is constructed for LIS (including ISS-LIS) and OTD that accounts for the horizontal dimension of lightning. The FED climatology documents the frequency that an observer can expect lightning to be visible overhead—regardless of where the flash began or ended. This new FED climatology confirms and elaborates on the previous global LIS/OTD FRD and Americas-only Geostationary Lightning Mapper (GLM) findings. The FED climatology maintains Lake Maracaibo as the global lightning hotspot with an average of 389 flashes/day but designates Karabre in the Democratic Republic of the Congo as the global thunderstorm duty (percent of the total viewtime where lightning is observed) hotspot at 7.29%. Meanwhile, Kuala Lumpur is the national capital city with the most lightning, and its airport (KUL) is the top major airport affected by lightning (in terms of duty). The FED seasonal cycle and month-to-month changes in the “center of lightning” for the three continental chimney regions are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Role of the chiral anomaly in polarized deeply inelastic scattering. III. Wess-Zumino-Witten contributions and chiral Ward identities for finite quark mass

We extend our prior results on the worldline computation of the axial vector-vector-vector (AVV) triangle anomaly in polarized deeply inelastic scattering (DIS) to the finite mass case by computing in addition the pseudoscalar-vector-vector (PVV) triangle graph. For the well-studied QED case, we show explicitly how the off-forward AVV pole exactly cancels an identical PVV pole. We then demonstrate the dramatic difference in QCD due to the chiral condensate, which qualitatively modifies anomalous Ward identities. As in the massless case, the anomaly pole in QCD is canceled by the dynamics of a primordial isosinglet pseudoscalar η ¯ -meson, whose Wess-Zumino-Witten coupling to the topological charge density shifts the pole to the physical η ′ mass, with the finite quark mass contribution differing by O ( 10 % ) from the Witten-Veneziano formula. We obtain a compact analytic expression for the finite mass corrections to Shore and Veneziano’s result that the proton’s net quark helicity Δ Σ ∝ χ QCD ′ | m = 0 ( 0 ) , the forward slope of the topological susceptibility in the chiral limit, and show they are of the order of a few percentages. Our prior prediction that the polarized DIS structure function g 1 is quenched by sphaleronlike topological transitions at small x is unaffected by quark mass effects. Our results illustrate how worldline computations of anomalous processes, in synergy with lattice computations and nonet chiral perturbation theory, can uncover novel nonperturbative features of QCD at the Electron-Ion Collider. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Topological band gap in intercalated epitaxial graphene

The study of functional manipulation of graphene is a critical subject, both for fundamental research and practical applications. In this study, we present that the intercalation of 5d transition metals into epitaxial graphene on SiC is a promising strategy for the realization of topologically nontrivial phases with a finite band gap in graphene. Employing first-principles calculations, grounded in density functional theory, we demonstrate that Re- and Ta-intercalated graphene evolve into two-dimensional topological insulators. These exhibit linear Dirac cones and quadratic bands with topological band gaps, respectively. The emergence of these topological states is attributed to the strong spin–orbit coupling strength of the intercalants. We show that the corresponding topological edge states persist within the finite bulk band gap, aligning with the bulk-boundary correspondence. Furthermore, we explore the spin splitting of the band structure, brought about by the inversion symmetry breaking and the spin–orbit coupling. Finally, our study underscores that the intercalation of graphene is an effective and a feasible approach for manipulating the band gap and the topological nature of graphene. Such intercalated graphene systems hold potential utility for spintronics and low-dimensional quantum device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topological Doping and Superconductivity in Cuprates: An Experimental Perspective

Hole doping into a correlated antiferromagnet leads to topological stripe correlations, involving charge stripes that separate antiferromagnetic spin stripes of opposite phases. The topological spin stripe order causes the spin degrees of freedom within the charge stripes to feel a geometric frustration with their environment. In the case of cuprates, where the charge stripes have the character of a hole-doped two-leg spin ladder, with corresponding pairing correlations, anti-phase Josephson coupling across the spin stripes can lead to a pair-density-wave order in which the broken translation symmetry of the superconducting wave function is accommodated by pairs with finite momentum. This scenario is now experimentally verified by recently reported measurements on La2−xBaxCuO4 with x=1/8. While pair-density-wave order is not common as a cuprate ground state, it provides a basis for understanding the uniform d-wave order that is more typical in superconducting cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Charges on a suspended silicon nitride membrane under a high-energy electron beam

Thin silicon nitride (SiN x ) membranes are widely used in gas and liquid phase transmission electron microscopy (TEM) and as phase plates to enhance imaging contrast. SiN x contains trap sites for both positive and negative charges, which can be manipulated by high-energy electron irradiation, external potential biasing, or light exposure. Charge accumulation on the membrane can significantly affect in situ TEM processes, including chemical and electrochemical reactions, nanoparticle dynamics, and catalytic activity, or introduce unwanted phase shifts when used as a phase plate. Here, in this study, charge accumulation on suspended SiN x membranes was investigated using off-axis electron holography combined with model-free charge analysis, supported by custom finite element analysis (FEA) simulations. An average residual positive charge density of approximately 2.8 × 10 −4 C m −2 was measured. Localized and stable regions of both positive and negative charges were identified on the membrane. The global positive and localized positive/negative charges give rise to strong electric fields and electroosmotic slip velocities at the membrane surface, which are sufficient to induce non-Brownian particle behavior and directional fluid flow, offering a physical explanation for previously observed anomalies in particle dynamics, nucleation, and growth during gas and liquid phase TEM experiments. These results provide a benchmark for understanding charge behavior at SiNx interfaces in gas and liquid phase TEM. Furthermore, the FEA simulations establish a framework for future investigations into charge distribution, electrostatic potentials, and electrical double layers at solid–liquid interfaces, particularly in complex geometries and chemically dynamic environments.

Suspended silicon↗