Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initial state”

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

Final state interactions for high energy scattering off atomic electrons

We consider the scattering of high energy leptons off bound atomic electrons focusing primarily on final state interactions between the outgoing energetic electron, and the heavy residual charged “debris” in the final state. These effects are inherently absent from calculations for a free electron at rest. Coulomb exchanges are enhanced by the large number of electrons in the atomic debris, and are unsuppressed by nonrelativistic velocities in the debris. We find that these exchanges can be resummed using operator methods, and cancel at the level of the cross section until at least 𝑂⁡(𝛼 3 ). Furthermore, we argue that both final and initial state Coulomb exchanges (enhanced by the number of electrons in the atom) do not affect the cross section until at least 𝑂⁡(𝛼 3 ). Transverse photon couplings to nonrelativistic electrons are proportional to their small velocities, and rotational invariance suppresses their contribution to 𝑂⁡(𝛼 3 ). Our results are relevant for precision experiments involving neutrinos, electrons, positrons, and muons scattering off of atomic electrons in a fixed target.

Atoms↗
DOE OSTI2024

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A General Decentralized Dynamic State Estimation With Synchronous Generator Magnetic Saturation

We report the saturation of the nonlinear magnetic circuit of synchronous generators is often neglected when performing Kalman filter-based dynamic state estimation (DSE), yielding significant estimation bias. This letter addresses this problem and proposes a generalized DSE framework to handle magnetic saturation. Moreover, this letter derives a state initialization procedure that improves the Kalman filter tracking speed. The framework is flexible in dealing with different saturation functions and generator models. Numerical results on the Texas 2000-bus system verify the effectiveness of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Solvent Polarity Independent Symmetry-Breaking Charge Separation in a Slip-Stacked Covalent Terrylene Monoimide Dimer

The design of organic materials capable of efficient photoinduced charge separation in low-polarity environments is a critical challenge for advancing organic photovoltaics. Symmetry-breaking charge separation (SB-CS) offers a promising route to efficient charge separation with minimal energy loss; however, in conventional organic materials this process typically relies on polar solvents to stabilize the charge-separated state. Here, we investigate a slip-stacked terrylene monoimide dimer (TMI 2 ) engineered with intrinsic electronic asymmetry to circumvent this limitation. We demonstrate through a comprehensive suite of ultrafast spectroscopic techniques that TMI 2 undergoes remarkably efficient SB-CS regardless of solvent polarity enabled by the permanent dipole moment of TMI and the accessible intramolecular charge-transfer (ICT) character of its monomer units. Furthermore, analysis of vibronic coherences reveals that the initial state mixing is actively driven by a 193 cm –1 intermolecular mode, while the final SB-CS state is marked by a distinct vibrational fingerprint of the radical ion pair product. Our findings reveal how the intrinsic electronic asymmetry of the TMI monomers that constitute TMI 2 creates an efficient charge separation pathway that precludes the need for external stabilization from a polar solvent. Here, this work establishes a powerful molecular design principle for leveraging intrinsic monomeric asymmetry to achieve efficient charge separation in low-dielectric environments, with significant implications for the development of next-generation organic optoelectronic materials.

Aromatic compounds↗
DOE OSTI2022

Measuring the electron Yukawa coupling via resonant s-channel Higgs production at FCC-ee

The Future Circular Collider (FCC-ee) offers the unique opportunity of studying the Higgs Yukawa coupling to the electron, \(y_\mathrm {e}\) , via resonant s -channel production, \(\mathrm {e^+e^-}\rightarrow \mathrm {H}\) , in a dedicated run at \(\sqrt{s} = m_\mathrm {H}\) . The signature for direct Higgs production is a small rise in the cross sections for particular final states, consistent with Higgs decays, over the expectations for their occurrence due to Standard Model (SM) background processes involving \(\mathrm {Z}^*\) , \(\gamma ^*\) , or t -channel exchanges alone. Performing such a measurement is remarkably challenging for four main reasons. First, the low value of the e \(^\pm \) mass leads to a tiny \(y_\mathrm {e}\) coupling and correspondingly small cross section: \(\sigma _\mathrm {ee\rightarrow H} \propto m_\mathrm {e}^2 = 0.57\) fb accounting for initial-state \(\gamma \) radiation. Second, the \(\mathrm {e^+e^-}\) beams must be monochromatized such that the spread of their centre-of-mass (c.m.) energy is commensurate with the narrow width of the SM Higgs boson, \(\varGamma _\mathrm {H} = 4.1\) MeV, while keeping large beam luminosities. Third, the Higgs mass must also be known beforehand with a few-MeV accuracy in order to operate the collider at the resonance peak, \(\sqrt{s} = m_\mathrm {H}\) . Last but not least, the cross sections of the background processes are many orders-of-magnitude larger than those of the Higgs decay signals. A preliminary generator-level study of 11 Higgs decay channels using a multivariate analysis, which exploits boosted decision trees to discriminate signal and background events, identifies two final states as the most promising ones in terms of statistical significance: \(\mathrm {H}\rightarrow gg\) and \(\mathrm {H}\rightarrow \mathrm {W}\mathrm {W}^*\!\rightarrow \ell \nu \) + 2 jets. For a benchmark monochromatization with 4.1-MeV c.m. energy spread (leading to \(\sigma _\mathrm {ee\rightarrow H} = 0.28\) fb) and 10 ab \(^{-1}\) of integrated luminosity, a \(1.3\sigma \) signal significance can be reached, corresponding to an upper limit on the e \(^\pm \) Yukawa coupling at 1.6 times the SM value: |y e | < 1.6|y e SM | at 95% confidence level, per FCC-ee interaction point per year. Directions for future improvements of the study are outlined.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Chiral effects and Joule heating in hot and dense matter

Initial states of dense matter with nonzero electron chiral imbalance could potentially give rise to strong magnetic fields through chiral plasma instability. Previous work indicated that unless chiral chemical potential is as large as the electron vector chemical potential, the growth of magnetic fields due to the instability is washed out by chirality flipping rate enabled by electron mass. We reexamine this claim in a broader range of parameters and find that at higher temperatures the hierarchy is reversed supporting a growing magnetic field for an initial electron chiral chemical potential much smaller than the electron vector chemical potential. Further, we identify a qualitatively new effect relevant for magnetized hot and dense medium where chiral magnetic effect (CME) sourced by density fluctuation acts as a powerful source of Joule heating. Remarkably, even modest chiral chemical potentials (keV) in such environment can deposit energy densities set by the QCD scale in a relatively short time of the order of a few milliseconds or seconds. We speculate how this mechanism makes CME-driven Joule heating a potentially critical ingredient in the dynamics of turbulent density fluctuation of supernovae and neutron star mergers.

Particle astrophysics↗
DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Efficient state preparation for the Schwinger model with a theta term

We present a comparison of different quantum state preparation algorithms and their overall efficiency for the Schwinger model with a theta term. While adiabatic state preparation is proved to be effective, in practice it leads to large gate counts to prepare the ground state. The quantum approximate optimization algorithm (QAOA) provides excellent results while keeping the counts small by design, at the cost of an expensive classical minimization process. We introduce a “blocked” modification of the Schwinger Hamiltonian to be used in the QAOA that further decreases the length of the algorithms as the size of the problem is increased. The rodeo algorithm (RA) provides a powerful tool to efficiently prepare any eigenstate of the Hamiltonian, as long as its overlap with the initial guess is large enough. We obtain the best results when combining the blocked QAOA ansatz and the RA, as this provides an excellent initial state with a relatively short algorithm without the need to perform any classical steps for large problem sizes. Published by the American Physical Society 2025

Bazavov, Alexei (ORCID:0000000321411901)↗
DOE OSTI2021

Discrete Time-Crystalline Order Enabled by Quantum Many-Body Scars: Entanglement Steering via Periodic Driving

The control of many-body quantum dynamics in complex systems is a key challenge in the quest to reliably produce and manipulate large-scale quantum entangled states. Recently, quench experiments in Rydberg atom arrays [Bluvstein et al. Science 371, 1355 (2021)] demonstrated that coherent revivals associated with quantum many-body scars can be stabilized by periodic driving, generating stable subharmonic responses over a wide parameter regime. We analyze a simple, related model where these phenomena originate from spatiotemporal ordering in an effective Floquet unitary, corresponding to discrete time-crystalline behavior in a prethermal regime. Unlike conventional discrete time crystals, the subharmonic response exists only for Néel-like initial states, associated with quantum scars. We predict robustness to perturbations and identify emergent timescales that could be observed in future experiments. Our results suggest a route to controlling entanglement in interacting quantum systems by combining periodic driving with many-body scars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

O 16 O 16 collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider comparing α clustering versus substructure

Collisions of . light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed 16 O 16 O run at the CERN Large Hadron Collider (LHC) and at the BNL Relativistic Heavy Ion Collider (RHIC) we study the potential for finding α clustering in 16 O. Here we use the state-of-the-art iEBE-VISHNU package with 16 O nucleonic configurations from ab initio nuclear lattice simulations. This setup was tuned using a Bayesian analysis on p Pb and PbPb systems. We find that the 16 O 16 O system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between α-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
DOE OSTI2024

Dark Matter searches with photons at the LHC

We unveil blind spot regions in dark matter (DM) direct detection (DMDD), for weakly interacting massive particles with a mass around a few hundred GeV that may reveal interesting photon signals at the LHC. We explore a scenario where the DM primarily originates from the singlet sector within the Z 3 -symmetric Next-to-Minimal Supersymmetric Standard Model (NMSSM). A novel DMDD spin-independent blind spot condition is revealed for singlino-dominated DM, in cases where the mass parameters of the higgsino and the singlino-dominated lightest supersymmetric particle (LSP) exhibit opposite relative signs (i.e., κ < 0), emphasizing the role of nearby bino and higgsino-like states in tempering the singlino-dominated LSP. Additionally, proximate bino and/or higgsino states can act as co-annihilation partner(s) for singlino-dominated DM, ensuring agreement with the observed relic abundance of DM. Remarkably, in scenarios involving singlino-higgsino co-annihilation, higgsino-like neutralinos can distinctly favor radiative decay modes into the singlino-dominated LSP and a photon, as opposed to decays into leptons/hadrons. In exploring this region of parameter space within the singlino-higgsino compressed scenario, we study the signal associated with at least one relatively soft photon alongside a lepton, accompanied by substantial missing transverse energy (E T ) and a hard initial state radiation jet at the LHC. In the context of singlino-bino co-annihilation, the bino state, as the next-to-LSP, exhibits significant radiative decay into a soft photon and the LSP, enabling the possible exploration at the LHC through the triggering of this soft photon alongside large E T and relatively hard leptons/jets resulting from the decay of heavier higgsino-like states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A Bayesian Approach for In-Situ Stress Prediction and Uncertainty Quantification for Subsurface Engineering

Many subsurface engineering applications require accurate knowledge of the in-situ state of stress for their safe design and operation. Existing methods to meet this need primarily include field measurements for estimating one or more of the principal stresses from a borehole, or optimization methods for constructing a 3D geomechanical model in terms of geophysical measurements. These methods, however, often contain considerable uncertainty in estimating the state of stress. Here, in this paper, we build on a Bayesian approach to quantify uncertainty in stress estimations for subsurface engineering applications. This approach can provide an estimate of the 3D distribution of stress throughout the volume of interest and provide an estimate of the uncertainty arising from the stress measurement, the rheology parameters, and a paucity of measurements. The value of this approach is demonstrated using stress measurements from the In Salah carbon storage site, which was one of the world’s first industrial carbon capture and storage projects. This demonstration shows the application of this Bayesian approach for estimating the initial state of stress for In Salah and quantifying the uncertainty in the estimated stress. Also, an assessment of a maximum injection pressure to prevent geomechanical risks from CO 2 injection pressures is provided in terms of the probability distribution of the minimum principal stress quantified by the approach. With the In Salah case study, this paper demonstrates that using the Bayesian approach can provide additional insights for site explorations and/or project operations to make informed-site decisions for subsurface engineering applications.

58 GEOSCIENCES↗
DOE OSTI2021

One-directional polarization transport in electron/nuclear spin chains with loss and gain

Understanding the joint dynamics of electron and nuclear spins is central to core concepts in solid-state magnetic resonance such as spin-lattice relaxation and dynamic nuclear polarization, but a generalization that capitalizes on competing polarization loss and gain channels is still lacking. Here, we study theoretically the non-Hermitian dynamics of hybrid electron/nuclear spin systems in the simultaneous presence of electron spin pumping and spin-lattice relaxation. Focusing on periodic, one-dimensional chains, we find that, by adjusting the electron spin pumping to a critical level, it is possible to steer the flow of nuclear polarization to create site-dependent distributions where either end of the array polarizes in opposite ways, irrespective of the initial state. By contrast, we show that ring-like patterns—where the limit nuclear polarization is uniform—exhibit a non-decaying, externally driven nuclear spin current. Interestingly, cyclic magnetic field modulation can render these processes largely robust to defects in the chain, a response featuring some interesting similarities—and differences—with recent findings in other non-Hermitian physical platforms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Search for the Z γ decay mode of new high-mass resonances in pp collisions at s = 13 TeV with the ATLAS detector

This letter presents a search for narrow, high-mass resonances in the Zγ final state with the Z boson decaying into a pair of electrons or muons. The $\sqrt{s}$ =13 TeV pp collision data were recorded by the ATLAS detector at the CERN Large Hadron Collider and have an integrated luminosity of 140 fb -1 . The data are found to be in agreement with the Standard Model background expectation. Upper limits are set on the resonance production cross section times the decay branching ratio into Zγ. For spin-0 resonances produced via gluon-gluon fusion, the observed limits at 95% confidence level vary between 65.5 fb and 0.6 fb, while for spin-2 resonances produced via gluon-gluon fusion (or quark-antiquark initial states) limits vary between 77.4 (76.1) fb and 0.6 (0.5) fb, for the mass range from 220 GeV to 3400 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Classical optimization with imaginary-time block encoding on quantum computers: The MaxCut problem

Optimization problems in finance, physics, and computer science are typically very hard to tackle in classical computing; quantum computing could help speed up computations and provide efficient methods for tackling large problems. Typically, to treat a problem with a quantum computer, the optimal solution is cast as the ground state of a diagonal Hamiltonian. Here, we develop a method, called imaginary-time evolution block encoding (ITE-BE), based on a recent imaginary-time algorithm, which requires no variational parameter optimization, as all parameters can be derived analytically from the target Hamiltonian. We also demonstrate that our method can be successfully combined with other quantum algorithms such as the quantum approximate optimization algorithm (QAOA). For illustration, here we study the MaxCut problem. We find that the QAOA ansatz increases the postselection success of ITE-BE, and shallow QAOA circuits, when boosted with ITE-BE, achieve better performance than deeper QAOA circuits. For the special case of the transverse initial state, we adapt our block-encoding scheme to allow for a deterministic application of the first layer of the circuit.

Zhong, Dawei [University of Southern California, L↗
DOE OSTI2022

The Role of Seasonality and the ENSO Mode in Central and East Pacific ENSO Growth and Evolution

Abstract A cyclostationary linear inverse model (CSLIM) is used to investigate the seasonal growth of tropical Pacific Ocean El Niño–Southern Oscillation (ENSO) events with canonical, central Pacific (CP), or eastern Pacific (EP) sea surface temperature (SST) characteristics. Analysis shows that all types of ENSO events experience maximum growth toward final states occurring in November and December. ENSO events with EP characteristics also experience growth into May and June, but CP events do not. A single dominant “ENSO mode,” growing from an equatorial heat content anomaly into a characteristic ENSO-type SST pattern in about 9 months (consistent with the delayed/recharge oscillator model of ENSO), is essential for the predictable development of all ENSO events. Notably, its seasonality is responsible for the late-calendar-year maximum in ENSO amplification. However, this ENSO mode alone does not capture the observed growth and evolution of diverse ENSO events, which additionally involve the seasonal evolution of other nonorthogonal Floquet modes. EP event growth occurs when the ENSO mode is initially “covered up” in combination with other Floquet modes. The ENSO mode’s slow seasonal evolution allows it to emerge while the other modes rapidly evolve and/or decay, leading to strongly amplifying and more predictable EP events. CP events develop when the initial state has a substantial contribution from Floquet modes with meridional mode–like SST structures. Thus, while nearly all ENSO events involve the seasonally varying ENSO-mode dynamics, the diversity and predictability of ENSO events cannot be understood without identifying contributions from the remaining Floquet modes. Significance Statement The purpose of this study is to identify structures that lead to seasonal growth of diverse types of El Niño–Southern Oscillation (ENSO) events. An important contribution from this study is that it uses an observationally constrained, empirically derived seasonal model. We find that processes affecting the evolution of diverse ENSO events are strongly seasonally dependent. ENSO events with eastern equatorial Pacific sea surface temperature (SST) characteristics are closely related to a single “ENSO mode” that resembles theoretical models of ENSO variability. ENSO events that have central equatorial Pacific SST characteristics include contributions from additional “meridional mode” structures that evolve via different physical processes. These findings are an important step in evaluating the seasonal predictability of ENSO diversity.

Vimont, Daniel J.↗