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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Quantum Computational Determination of the Weak Mixing Angle in the Standard Model

The weak mixing angle $s_W$ is a fundamental constant in the Standard Model (SM) and measured at the $Z$ boson mass to be $\widehat{s}^2_W(m_Z) = 0.23129 \pm 0.00004$ in the $\overline{\rm MS}$ renormalization scheme, where $m_Z=91.2\ \text{GeV}$. On the other hand, non-stabilizerness - the magic - characterizes the computational advantage of a quantum system over classical computers. We consider the production of magic from stabilizer initial states, which carry zero magic, in the 2-to-2 scattering of charged leptons in the SM at the tree level, which is mediated by the photon and the $Z$ boson. Using the second order stabilizer Rényi entropy, and averaging over all 60 initial stabilizer states and the scattering angle, we compute and minimize the magic production as a function of $s^2_W$ in the Møller scattering $e^-e^-\to e^-e^-$, which is free of kinematic thresholds. At the centre-of-mass energy $\sqrt{s}=m_Z$, there is a unique minimum in magic production at $\mathbf{s}^{2}_W(m_Z)=0.2317$, which agrees with the measured $\widehat{s}^2_W(m_Z)$ at the sub-percent level. At higher energies, the magic-minimizing $\mathbf{s}^{2}_W$ continues to agree with the empirical value at the percent level or better, up to 10 TeV. The finding suggests the electroweak sector of the SM tends to generate minimal quantum resources from the computational viewpoint.

Liu, Qiaofeng [Northwestern U.]↗
DOE OSTI2025

A Quantum Computational Determination of the Weak Mixing Angle in the Standard Model

The weak mixing angle $s_W$ is a fundamental constant in the Standard Model (SM) and measured at the $Z$ boson mass to be $\widehat{s}^2_W(m_Z) = 0.23129 \pm 0.00004$ in the $\overline{\rm MS}$ renormalization scheme, where $m_Z=91.2\ \text{GeV}$. On the other hand, non-stabilizerness - the magic - characterizes the computational advantage of a quantum system over classical computers. We consider the production of magic from stabilizer initial states, which carry zero magic, in the 2-to-2 scattering of charged leptons in the SM at the tree level, which is mediated by the photon and the $Z$ boson. Using the second order stabilizer Rényi entropy, and averaging over all 60 initial stabilizer states and the scattering angle, we compute and minimize the magic production as a function of $s^2_W$ in the Møller scattering $e^-e^-\to e^-e^-$, which is free of kinematic thresholds. At the centre-of-mass energy $\sqrt{s}=m_Z$, there is a unique minimum in magic production at $\mathbf{s}^{2}_W(m_Z)=0.2317$, which agrees with the measured $\widehat{s}^2_W(m_Z)$ at the sub-percent level. At higher energies, the magic-minimizing $\mathbf{s}^{2}_W$ continues to agree with the empirical value at the percent level or better, up to 10 TeV. The finding suggests the electroweak sector of the SM tends to generate minimal quantum resources from the computational viewpoint.

Liu, Qiaofeng [Northwestern U.]↗
DOE OSTI2022

Transverse Single Spin Asymmetries of Heavy Flavor Electrons and Charged Pions in 200 GeV $p+p^{\uparrow}$ Collisions at Midrapidity

Transverse single spin asymmetries of particles produced in p+p^{\uparrow} p + p ↑ collisions provide insight on the partonic spin and momentum structure of hadrons; heavy flavor electrons provide access to initial state spin-momentum correlations of gluons in the proton, while charged pions provide access to initial and final state transverse spin effects. Charged particles are measured at midrapidity at PHENIX using the silicon vertex detector and central arm spectrometer, made of an electromagnetic calorimeter, a ring-imaging Cherenkov detector, and drift and pad chambers. Recent results for both electron and charged pion measurements from the 2015 running period will be presented.

Fitzgerald, Dillon↗
DOE OSTI2024

Adaptive variational ground state preparation for spin-1 models on qubit-based architectures

Here, we apply the adaptive variational quantum imaginary time evolution (AVQITE) method to prepare ground states of one-dimensional spin S = 1 models. We compare different spin-to-qubit encodings (standard binary, Gray, unary, and multiplet) with regard to the performance and quantum resource cost of the algorithm. Using state-vector simulations, we study two well-known spin-1 models: the Blume-Capel model of transverse-field Ising spins with single-ion anisotropy, and the XXZ model with single-ion anisotropy. We consider system sizes of up to 20 qubits, which corresponds to spin-1 chains up to length 10. We determine the dependence of the number of CNOT gates in the AVQITE state preparation circuit on the encoding, the initial state, and the choice of operator pool in the adaptive method. Independent of the choice of encoding, we find that the CNOT gate count scales cubically with the number of spins for the Blume-Capel model and quartically for the anisotropic XXZ model. However, the multiplet and Gray encodings present smaller prefactors in the scaling relations. These results provide useful insights for the implementation of AVQITE on quantum hardware.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Power-law inflation satisfies Penrose’s Weyl curvature hypothesis

Based on entropy considerations and the arrow of time Penrose argued that the Universe must have started in a special initial singularity with vanishing Weyl curvature. This is often interpreted to be at odds with inflation. Here we argue just the opposite, that Penrose’s persuasions are in fact consistent with inflation. Using the example of power law inflation, we show that inflation begins with a past null singularity, where Weyl tensor vanishes when the metric is initially exactly conformally flat. This initial state precisely obeys Penrose’s conditions. The initial null singularity breaks T-reversal spontaneously and picks the arrow of time. It can be regulated and interpreted as a creation of a universe from nothing, initially fitting in a bubble of Planckian size when it materializes. Penrose’s initial conditions are favored by the initial O(4) symmetry of the bubble, selected by extremality of the regulated Euclidean action. The predicted observables are marginally in tension with the data, but they can fit if small corrections to power law inflation kick in during the last 60⁢e-folds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

ADoSPSiQS (Automated Detection of Symmetry-Protected Subspaces in Quantum Simulations ) [SWR-23-55]

Symmetry-protected subspaces in quantum systems are subspaces demarcated by a conserved quantity of that quantum system. We introduce two classical algorithms, which efficiently compute and elucidate features of these subspaces. The first algorithm explores the entire symmetry-protected subspace of an initial state in time complexity linear to the size of the subspace by closing local basis state-to-basis state transitions. The second algorithm determines, with bounded error, if a given measurement outcome of a dynamically-generated state is within the symmetry-protected subspace of the state in which the dynamical system is initialized. This software is paired with the paper of the same name at https://arxiv.org/pdf/2302.08586.pdf.

Rotello, Caleb↗
DOE OSTI2023

Time-Energy Uncertainty Relation for Noisy Quantum Metrology

Detection of very weak forces and precise measurement of time are two of the many applications of quantum metrology to science and technology. To sense an unknown physical parameter, one prepares an initial state of a probe system, allows the probe to evolve as governed by a Hamiltonian 𝐻 for some time 𝑡, and then measures the probe. If 𝐻 is known, we can estimate 𝑡 by this method; if 𝑡 is known, we can estimate classical parameters on which 𝐻 depends. The accuracy of a quantum sensor can be limited by either intrinsic quantum noise or by noise arising from the interactions of the probe with its environment. In this work, we introduce and study a fundamental trade-off, which relates the amount by which noise reduces the accuracy of a quantum clock to the amount of information about the energy of the clock that leaks to the environment. Specifically, we consider an idealized scenario in which a party Alice prepares an initial pure state of the clock, allows the clock to evolve for a time that is not precisely known, and then transmits the clock through a noisy channel to a party Bob. Meanwhile, the environment (Eve) receives any information about the clock that is lost during transmission. We prove that Bob’s loss of quantum Fisher information about the elapsed time is equal to Eve’s gain of quantum Fisher information about a complementary energy parameter. We also prove a similar, but more general, trade-off that applies when Bob and Eve wish to estimate the values of parameters associated with two noncommuting observables. We derive the necessary and sufficient conditions for the accuracy of the clock to be unaffected by the noise, which form a subset of the Knill-Laflamme error-correction conditions. A state and its local time-evolution direction, if they satisfy these conditions, are said to form a metrological code. We provide a scheme to construct metrological codes in the stabilizer formalism. We show that there are metrological codes that cannot be written as a quantum error-correcting code with similar distance in which the Hamiltonian acts as a logical operator, potentially offering new schemes for constructing states that do not lose any sensitivity upon application of a noisy channel. We discuss applications of the trade-off relation to sensing using a quantum many-body probe subject to erasure or amplitude-damping noise.

metrology↗
DOE OSTI2021

TMD fragmentation functions at N 3 LO

We compute the unpolarized quark and gluon transverse-momentum dependent fragmentation functions (TMDFFs) at next-to-next-to-next-to-leading order (N 3 LO) in perturbative QCD. The calculation is based on a relation between the TMDFF and the limit of the semi-inclusive deep inelastic scattering cross section where all final-state radiation becomes collinear to the detected hadron. The required cross section is obtained by analytically continuing our recent computation of the Drell-Yan and Higgs boson production cross section at N 3 LO expanded around the limit of all final-state radiation becoming collinear to one of the initial states. Our results agree with a recent independent calculation by Luo et al.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The BEST framework for the search for the QCD critical point and the chiral magnetic effect

The Beam Energy Scan Theory (BEST) Collaboration was formed with the goal of providing a theoretical framework for analyzing data from the Beam Energy Scan (BES) program at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory. The physics goal of the BES program is the search for a conjectured QCD critical point as well as for manifestations of the chiral magnetic effect. We describe progress that has been made over the previous five years. This includes studies of the equation of state and equilibrium susceptibilities, the development of suitable initial state models, progress in constructing a hydrodynamic framework that includes fluctuations and anomalous transport effects, as well as the development of freezeout prescriptions and hadronic transport models. Lastly, we address the challenge of integrating these components into a complete analysis framework. This document describes the collective effort of the BEST Collaboration and its collaborators around the world.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Neutrino many-body flavor evolution: The full Hamiltonian

We study neutrino flavor evolution in the quantum many-body approach using the full neutrino-neutrino Hamiltonian, including the usually neglected terms that mediate nonforward scattering processes. Working in the occupation number representation with plane waves as single-particle states, we explore the time evolution of simple initial states with up to N = 10 neutrinos. We discuss the time evolution of the Loschmidt echo, one body flavor and kinetic observables, and the one-body entanglement entropy. For the small systems considered, we observe “thermalization” of both flavor and momentum degrees of freedom on comparable time scales, with results converging towards expectation values computed within a microcanonical ensemble. We also observe that the inclusion of nonforward processes generates a faster flavor evolution compared to the one induced by the truncated (forward) Hamiltonian. Published by the American Physical Society 2024

Cirigliano, Vincenzo (ORCID:000000029056754X)↗
DOE OSTI2022

Study of Short-Range nuclear Correlations in light nuclei using the BeAGLE event generator

Nuclear dynamics at short distances among nucleons is one of the most outstanding phenomena in nuclear physics. Understanding the role of QCD in generating nuclear forces is important for uncovering the underlying physics of Short-Range Correlations (SRCs). In recent years, SRCs has been observed from light to heavy nuclei using fixed target experiments at Jefferson lab via high energy electron-nucleus scattering. In this talk, I will talk about the opportunity of studying SRCs using light nuclei with collider experiments, e.g., the Electron-Ion Collider (EIC). The experimental technique of studying the light nuclei can be based on exclusive processes with tagging final-state particles, in order to fully control the initial state of the target wavefunction. In particular, incoherent diffractive production of J/\psi J / ψ particle off deuteron will be presented. In addition, the spectral function in light nuclei has been recently modeled in the BeAGLE event generator, where the decay kinematics of the light nuclei and their influence on the very forward detector design at the EIC will be discussed.

Tu, Zhoudunming↗
DOE OSTI2023

Evaluating skill in predicting the Interdecadal Pacific Oscillation in initialized decadal climate prediction hindcasts in E3SMv1 and CESM1 using two different initialization methods and a small set of start years

Abstract It is a daunting challenge to conduct initialized hindcasts with enough ensemble members and associated start years to form a drifted climatology from which to compute the anomalies necessary to quantify the skill of the hindcasts when compared to observations. This limits the ability to experiment with case studies and other applications where only a few initial years are needed. Here we run a set of hindcasts with CESM1 and E3SMv1 using two different initialization methods for a limited set of start years and use the respective uninitialized free-running historical simulations to form the model climatologies. Since the drifts from the observed initial states in the hindcasts toward the uninitialized model state are large and rapid, after a few years the drifted initialized models approach the uninitialized model climatological errors. Therefore, hindcasts from the limited start years can use the uninitialized climatology to represent the drifted model states after about lead year 3, providing a means to compute forecast anomalies in the absence of a large hindcast sample. There is comparable skill for predicting spatial patterns of multi-year Pacific sea surface temperature anomalies in the domain of the Interdecadal Pacific Oscillation using this method compared to the conventional methodology with a large hindcast data set, though there is a model dependence to the drifts in the two initialization methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗
DOE OSTI2026

Radiative Electronic Bound States in the Continuum from Defects in Semiconductors

Continuum-buried defect states in semiconductors are generally expected to be optically inactive because of their strong coupling to continuum bands. Here, we show that such defects can instead host radiative electronic bound states in the continuum (BICs) using the silicon G center as a prototypical example. Hybrid functional first-principles calculations with a Hubbard U correction reveal that a localized defect state, initially buried below the valence band maximum (VBM) in the ground state, undergoes exchange-driven energy-level reordering under optical excitation and shifts above the VBM. This exchange-induced transition suppresses nonradiative decay and enables a robust radiative emission. By computing temperature-dependent nonradiative lifetimes and comparing them to experimental photoluminescence (PL) lifetimes, we quantitatively reproduce the observed temperature dependence of the emission. These results uncover a stabilization mechanism for continuum-embedded defect states and establish electronic BICs as a general paradigm for designing defect-based optical systems, including quantum emitters and qubits.

Hong, SeongYun↗
DOE OSTI2021

Ultrafast x-ray pump x-ray probe transient absorption spectroscopy: A computational study and proposed experiment probing core-valence electronic correlations in solvated complexes

Femtosecond X-ray pump X-ray probe experiments are currently possible at free electron lasers such as the linac coherent light source (LCLS), which opens new opportunities for studying solvated transition metal-ligand complexes. In order to make the most effective use of these kinds of experiments, it is necessary to determine which chemical properties an X-ray pump pulse will measure. In this paper, we have combined electron cascade calculations and excited-state time-dependent density functional theory calculations to predict the initial state prepared by an X-ray pump and the subsequent X-ray probe spectra at the Fe K-edge in the solvated model transition metal complex, K 4 Fe II (CN) 6 . We find several key spectral features that report on the ligand-field splitting, solvent-solute interactions, and the 3p and 3d electron interactions. We then show how these features could be measured in an experiment.

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

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Improving SGP4 Orbit Determination with New State Estimation Algorithm

The Simplified General Perturbations 4 Model (SGP4) is a well-known tool for performing satellite orbit determination. However, uncertainties and inaccuracies in the initial state inputs (required by SGP4) degrade the performance of the propagator. We present a new state estimation algorithm that allows for independent computation of these initial inputs using Unscented Kalman Filtering and GPS data from a satellite. The algorithm is tested on real flight data and demonstrates a notable performance improvement over the standard method of orbit determination using SGP4.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Analysis of an Arctic cold air outbreak during autumn and related air mass transformations forced by surface changes and advection in higher altitudes

During marine cold air outbreaks (MCAOs), cold and dry Arctic air masses are transported from the central Arctic southward across the closed sea ice and much warmer open oceans. They experience significant transformations including a rapid heating and moistening, often leading to cloud formation. While intense wintertime MCAOs have been analyzed widely, the air mass transformations during other seasons have been studied sparsely. We address this gap by investigating an MCAO case observed in September 2020. To study the transformation processes, we combine the fifth generation of atmospheric reanalyses of the global climate (ERA5), trajectory calculations, as well as shipborne and airborne measurements. In the central Arctic, observations acquired from aboard the research vessel (RV) Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition characterized the initial state of the air mass over closed sea ice. Trajectories indicated the pathway the air mass took from RV Polarstern southward to the Fram Strait. For the first 24 h of the southbound drift, the air masses remained quasi-stationary. Then, still 15 h ahead of the marginal sea ice zone, differential advection across the boundary layer flow introduced humidity and clouds at higher altitudes between 1.5 and 2.5 km. ERA5-derived temperature and humidity tendencies indicated complex vertical interactions. Radiative cloud-top cooling, entrainment, and turbulence were significantly reduced in the lower and enhanced in the upper advected cloud layer. Eventually, the lower cloud deck dissipated. After this confluence of 2 different air masses, observations gathered by Polar 5 in Fram Strait as part of the MOSAiC Airborne observations in the Central Arctic campaign revealed cloudy, moist layers throughout the lowest 3.5 km and an increasing boundary layer height. Comparing the initial with the final state 48 h later, the largest net heating of +8 K was found close to the surface, yet the largest net moistening of +2.5 g kg -1 at an altitude of 1 km, as the initial profile was exceptionally dry here. We conclude that the observed air mass transformations were driven by the surface changes from sea ice to open ocean but additionally strongly impacted by the differential advection of clouds and moisture across the near-surface MCAO flow.

54 ENVIRONMENTAL SCIENCES↗