Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Volume”

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

Conjugate Multimode Heat Transfer Analysis of Cryogenic EXLO Manipulation

Abstract In this study, a conjugate radiation/conduction multimode heat transfer analysis of cryogenic focused ion beam (FIB) milling steps necessary for producing ex situ lift out specimens under cryogenic conditions (cryo-EXLO) is performed. Using finite volume for transient heat conduction and enclosure theory for radiation heat transfer, the analysis shows that as long as the specimen is attached or touching the FIB side wall trenches, the specimen will remain vitreous indefinitely, while actively cooled at liquid nitrogen (LN2) temperatures. To simulate the time needed to perform a transfer step to move the bulk sample containing the FIB-thinned specimen from the cryo-FIB to the cryo-EXLO cryostat, the LN2 temperature active cooling is turned off after steady-state conditions are reached and the specimen is monitored over time until the critical devitrification temperature is reached. Under these conditions, the sample will remain vitreous for >3 min, which is more than enough time needed to perform the cryo-transfer step from the FIB to the cryostat, which takes only ∼10 s. Cryo-transmission electron microscopy images of a manipulated cryo-EXLO yeast specimen prepared with cryo-FIB corroborates the heat transfer analysis.

Materials Science↗
DOE OSTI2022

Modelling the Lyman-α forest with Eulerian and SPH hydrodynamical methods

ABSTRACT We compare two state-of-the-art numerical codes to study the overall accuracy in modelling the intergalactic medium and reproducing Lyman-α forest observables for DESI and high-resolution data sets. The codes employ different approaches to solving both gravity and modelling the gas hydrodynamics. The first code, Nyx, solves the Poisson equation using the Particle-Mesh (PM) method and the Euler equations using a finite-volume method. The second code, CRK-HACC , uses a Tree-PM method to solve for gravity, and an improved Lagrangian smoothed particle hydrodynamics (SPH) technique, where fluid elements are modelled with particles, to treat the intergalactic gas. We compare the convergence behaviour of the codes in flux statistics as well as the degree to which the codes agree in the converged limit. We find good agreement overall with differences being less than observational uncertainties, and a particularly notable ≲1 per cent agreement in the 1D flux power spectrum. This agreement was achieved by applying a tessellation methodology for reconstructing the density in CRK-HACC instead of using an SPH kernel as is standard practice. We show that use of the SPH kernel can lead to significant and unnecessary biases in flux statistics; this is especially prominent at high redshifts, z ∼ 5, as the Lyman-α forest mostly comes from lower-density regions that are intrinsically poorly sampled by SPH particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Scattering Processes from Quantum Simulation Algorithms for Scalar Field Theories

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Hardy, Andrew [Toronto U.] (ORCID:0000000235817382↗
DOE OSTI2026

Rotating synchrotron radiation: Photon emission from magnetized and rotating quark-gluon plasma

This paper investigates the production of nonprompt photons originating from rotating synchrotron radiation (RoSyRa), specifically the emission of photons by a rigidly rotating quark-gluon plasma in thermal equilibrium, in the presence of an external magnetic field. We compute the nonprompt photon spectrum and its elliptic flow (𝑣 2 ) at midrapidity. In particular, we investigate the finite-volume effects. We find that at low transverse momentum, the magnetic field induces a significant 𝑣 2 , while the plasma rotation boosts the synchrotron radiation of negatively charged quarks. These findings make RoSyRa a viable candidate mechanism to resolve the “direct photon puzzle.”

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Two- and three-meson scattering amplitudes with physical quark masses from lattice QCD

We study systems of two and three mesons composed of pions and kaons at maximal isospin using four CLS ensembles with 𝑎 ≈ 0.063 fm, including one with approximately physical quark masses. Using the stochastic Laplacian-Heaviside method, we determine the energy spectrum of these systems including many levels in different momentum frames and irreducible representations. Using the relativistic two- and three-body finite-volume formalism, we constrain the two- and three-meson K matrices, including not only the leading 𝑠 wave, but also 𝑝 and 𝑑 waves. By solving the three-body integral equations, we determine, for the first time, the physical-point scattering amplitudes for 3⁢𝜋 + , 3⁢𝐾 + , 𝜋 + ⁢𝜋 + ⁢𝐾 + , and 𝐾 + ⁢𝐾 + ⁢𝜋 + systems. These are determined for total angular momentum 𝐽 𝑃 = 0 − , 1 + , and 2 − . We also obtain accurate results for 2⁢𝜋 + , 𝜋 + ⁢𝐾 + , and 2⁢𝐾 + phase shifts. We compare our results to chiral perturbation theory and to phenomenological fits.

FOS: Physical sciences↗
DOE OSTI2026

Di-nucleons do not form bound states at heavy pion mass

We perform a high-statistics lattice QCD calculation of the low-energy two-nucleon scattering amplitudes. To address discrepancies in the literature, the calculation is performed at a heavy pion mass in the limit that the light quark masses are equal to the physical strange quark mass, 𝑚 𝜋 = 𝑚 𝐾 ≃ 714 MeV. Using a state-of-the-art momentum space method, we rule out the presence of a bound di-nucleon in both the isospin 0 (deuteron) and 1 (di-neutron) channels, in contrast with many previous results that made use of compact hexaquark creation operators. To diagnose the discrepancy, we add such hexaquark interpolating operators to our basis and find that they do not affect the determination of the two-nucleon finite-volume spectrum, and thus they do not couple to deeply bound di-nucleons that are missed by the momentum-space operators. Furthermore, we perform a high-statistics calculation of the HAL QCD potential on the same gauge ensembles and find qualitative agreement with our main results. We conclude that di-nucleons do not form bound states at heavy pion masses and that previous identification of deeply bound di-nucleons must have arisen from a misidentification of the spectrum from off-diagonal elements of a correlation function.

Physics - Physics of elementary particles and fiel↗
DOE OSTI2022

Radiative properties of an artificial atom coupled to a Josephson-junction array

We study the radiative properties—the Lamb shift, Purcell decay rate, and spontaneous emission dynamics—of an artificial atom coupled to a long, multimode cavity formed by an array of Josephson junctions. Introducing a tunable coupling element between the atom and the array, we demonstrate that such a system can exhibit a crossover from a perturbative to a nonperturbative regime of light-matter interaction as one strengthens the coupling between the atom and the Josephson-junction array. As a consequence, the concept of spontaneous emission as the occupation of the local atomic site being governed by a single complex-valued exponent breaks down. This breakdown, we show, can be interpreted in terms of formation of hybrid atom-resonator modes with radiative losses that are nontrivially related to the effective coupling between individual modes. We develop a singular function expansion approach for the description of the open quantum system dynamics in such a multimode nonperturbative regime. This modal framework generalizes the normal-mode description of quantum fields in a finite volume, incorporating exact radiative losses and incident quantum noise at the delimiting surface. Our results are pertinent to recent experiments with Josephson atoms coupled to high-impedance Josephson-junction arrays.

42 ENGINEERING↗
DOE OSTI2022

Three-body resonances in pionless effective field theory

We investigate the appearance of resonances in three-body systems using pionless effective field theory at leading order with two complementary methods. The Faddeev equation is analytically continued to the unphysical sheet adjacent to the positive real energy axis using a contour rotation. We consider both the three-boson system and the three-neutron system. For the former, we calculate the trajectory of Borromean three-body Efimov states turning into resonances as they cross the three-body threshold. For the latter, we find no sign of three-body resonances or virtual states at leading order. Furthermore, this result is validated by exploring the level structure of three-body states in a finite volume approach.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Lattice model of heavy-light three-body system

We present a study of a $1+1$ dimensional heavy-light three-body system in finite volume. The heavy-light system is simulated by a coupled-channel $\phi^4$ type lattice model, and both ground state and excited states of multiparticle energy spectra are measured on various lattices. The lattice simulation data analysis is performed based on variational approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Dalitz plots and lineshape of a 1 ( 1260 ) from a relativistic three-body unitary approach

We formulate the final-state interaction of the a 1 ( 1260 ) resonance decay in a manifestly three-body unitary parametrization and fit it to the a 1 ( 1260 ) line shape measured by the ALEPH experiment. Dalitz plots calculated from this fit are presented. The work demonstrates the feasibility to numerically solve a previously derived amplitude and its generalization to isobars with spin and coupled channels. The model can also be applied to other meson decays and modified for the finite-volume problem as it arises in lattice QCD due to its manifest unitarity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Moments of nucleon isovector structure functions in 2+1+1-flavor QCD

We present results on the isovector momentum fraction, $\langle x \rangle_{u–d}$ , helicity moment, $\langle x \rangle_{Δ u–Δd}$ , and the transversity moment, $\langle x \rangle_{δu–δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC Collaboration using 2 + 1 + 1 -flavors of dynamical highly improved staggered quarks. The correlation functions are calculated using the Wilson-Clover action, and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI'–MOM scheme. The data have been collected at lattice spacings a ≈ 0.15 , 0.12, 0.09, and 0.06 fm and $M_π$ ≈ 310 , 220, and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{\text{MS}}$ scheme at 2 GeV, are $\langle$ x $\rangle$ u – d = 0.173 ( 14 ) ( 07 ) , $\langle x \rangle_{Δu–Δd}$ = 0.213 ( 15 ) ( 22 ) , and $\langle x \rangle_{δu–δd}$ = 0.208 ( 19 ) ( 24 ) , where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Charged multihadron systems in lattice QCD+QED

Systems with the quantum numbers of up to 12 charged and neutral pseudoscalar mesons, as well as one-, two-, and three-nucleon systems, are studied using dynamical lattice quantum chromodynamics and quantum electrodynamics (QCD+QED) calculations and effective field theory. QED effects on hadronic interactions are determined by comparing systems of charged and neutral hadrons after tuning the quark masses to remove strong isospin breaking effects. A nonrelativistic effective field theory, which perturbatively includes finite-volume Coulomb effects, is analyzed for systems of multiple charged hadrons and found to accurately reproduce the lattice QCD+QED results. QED effects on charged multihadron systems beyond Coulomb photon exchange are determined by comparing the two- and three-body interaction parameters extracted from the lattice QCD+QED results for charged and neutral multihadron systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Real-time quantum calculations of phase shifts using wave packet time delays

We present a method to extract the phase shift of a scattering process using the real-time evolution in the early and intermediate stages of the collision in order to estimate the time delay of a wave packet. This procedure is convenient when using noisy quantum computers for which the asymptotic out-state behavior is unreachable. We demonstrate that the challenging Fourier transforms involved in the state preparation and measurements can be implemented in 1+1 dimensions with current trapped ion devices and IBM quantum computers. We compare quantum computations of the time delays obtained in the one-particle quantum mechanics limit and the scalable quantum field theory formulation with accurate numerical results. We discuss the finite volume effects in the Wigner formula connecting time delays to phase shifts. The results reported involve two- and four-qubit calculations, and we discuss the possibility of larger scale computations in the near future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Extraction of low-energy constants of single- and double-β decays from lattice QCD: A sensitivity analysis

Lattice quantum chromodynamics (LQCD) has the promise of constraining low-energy constants (LECs) of nuclear effective field theories (EFTs) from first-principles calculations that incorporate the dynamics of quarks and gluons. Given the Euclidean and finite-volume nature of LQCD outputs, complex mappings are developed in recent years to obtain the Minkowski and infinite-volume counterparts of LQCD observables. In particular, as LQCD is moving toward computing a set of important few-nucleon matrix elements at the physical values of the quark masses, it is important to investigate whether the anticipated precision of LQCD spectra and matrix elements will be sufficient to guarantee tighter constraints on the relevant LECs than those already obtained from phenomenology, considering the nontrivial mappings involved. With a focus on the leading-order LECs of the pionless EFT, L 1,A and $g$$^{NN}_{v}$, which parametrize, respectively, the strength of the isovector axial two-body current in a single-β decay (and other related processes such pp fusion), and of the isotensor contact two-body operator in the neutrinoless double-β decay within the light neutrino exchange scenario, the expected uncertainty on future extractions of L 1,A and $g$$^{NN}_{v}$ are examined using synthetic data at the physical values of the quark masses. It is observed that achieving small uncertainties in L 1,A will be challenging, and (sub)percent-level precision in the two-nucleon spectra and matrix elements is essential in reducing the uncertainty on this LEC compared to the existing constraints. On the other hand, the short-distance coupling of the neutrinoless double-β decay, $g$$^{NN}_{v}$, is shown to be less sensitive to uncertainties on both LQCD energies and the matrix element, and can likely be constrained with percent-level precision in the upcoming LQCD calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Light-quark connected intermediate-window contributions to the muon 𝑔 − 2 hadronic vacuum polarization from lattice QCD

We present a lattice-QCD calculation of the light-quark connected contribution to window observables associated with the leading-order hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon, $a^{HVP,LO}_μ$. We employ the MILC Collaboration’s isospin-symmetric QCD gauge-field ensembles, which contain four flavors of dynamical highly improved staggered quarks with four lattice spacings between $a ≈ 0.06-0.15$ fm and close-to-physical quark masses. We consider several effective-field-theory-based schemes for finite volume and other lattice corrections and combine the results via Bayesian model averaging to obtain robust estimates of the associated systematic uncertainties. After unblinding, our final results for the intermediate and "W2" windows are $a^{U,W}_μ$(conn) = 206.6(1.0) x 10 -10 and $a^{U,W2}_μ$ (conn) = 100.7(3.2) x 10 -10 , respectively

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Lattice QCD study of π Σ − K ¯ N scattering and the Λ ( 1405 ) resonance

A lattice QCD computation of the coupled channel π Σ – K ¯ N scattering amplitudes in the Λ ( 1405 ) region is detailed. Results are obtained using a single ensemble of gauge field configurations with N f = 2 + 1 dynamical quark flavors and m π ≈ 200 MeV and m K ≈ 487 MeV . Hermitian correlation matrices using both single baryon and meson-baryon interpolating operators for a variety of different total momenta and irreducible representations are used. Several parametrizations of the two-channel scattering K -matrix are utilized to obtain the scattering amplitudes from the finite-volume spectrum. The amplitudes, continued to the complex energy plane, exhibit a virtual bound state below the π Σ threshold and a resonance pole just below the K ¯ N threshold. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dynamics of the O ( 4 ) critical point in QCD: Critical pions and diffusion in model G

We present a detailed study of the finite momentum dynamics of the O ( 4 ) critical point of QCD, which lies in the dynamic universality class of “model G.” The critical scaling of the model is analyzed in multiple dynamical channels. For instance, the finite momentum analysis allows us to precisely extract the pion dispersion curve below the critical point. The pion velocity is in striking agreement with the predictions relation and static universality. The pion damping rate and velocity are both consistent with the dynamical critical exponent ζ = 3 / 2 of model G. Similarly, although the critical amplitude for the diffusion coefficient of the conserved O ( 4 ) charges is small, it is clearly visible both in the restored phase and with finite explicit symmetry breaking, and its dynamical scaling is again consistent with ζ = 3 / 2 . We determine a new set of universal dynamical critical amplitude ratios relating the diffusion coefficient to a suitably defined order parameter relaxation time. We also show that in a finite volume simulation, the chiral condensate diffuses on the coset manifold in a manner consistent with dynamical scaling, and with a diffusion coefficient that is determined by the transport coefficients of hydrodynamic pions. Finally, the amplitude ratios (together with other nonuniversal amplitudes also reported here) compile all relevant information for further studies of model G both in and out of equilibrium. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Trace anomaly form factors from lattice QCD

The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the pion up to Q 2 ∼ 4.3 GeV 2 and the nucleon up to Q 2 ∼ 1 GeV 2 . The calculations are performed on a domain wall fermion ensemble with overlap valence quarks at seven valence pion masses varying from ∼ 250 to ∼ 540 MeV , including the unitary point ∼ 340 MeV . We calculate the radius of the glue trace anomaly for the pion and the nucleon from the z expansion. By performing a two-dimensional Fourier transform on the glue trace anomaly form factors in the infinite momentum frame with no energy transfer, we also obtain their spatial distributions for several valence quark masses. The results are qualitatively extrapolated to the physical valence pion mass with systematic errors from the unphysical sea quark mass, discretization effects in the renormalization sum rule, and finite-volume effects to be addressed in the future. We find the pion’s form factor changes sign, as does its spatial distribution, for light quark masses. This explains how the trace anomaly contribution to the pion mass approaches zero toward the chiral limit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗