SEARCH · Engineering Papers

Results for “Matrix elements”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Toward the determination of 𝐶⁢𝑃-odd pion-nucleon couplings

The nucleon matrix elements (NMEs) associated with quark chromomagnetic dipole moments (cMDMs) play a crucial role in determining the 𝐶𝑃-odd pion-nucleon couplings induced by quark chromoelectric dipole moments. In recent years, it has been argued that the NMEs of cMDMs can be related to the third moment of the nucleon's higher-twist (specifically, twist-3) parton distribution function (PDF) 𝑒⁡(𝑥), which can, in principle, be measured through dihadron production in semi-inclusive deep inelastic scattering processes. By applying the spin-flavor expansion to the cMDM operators in the large-𝑁 𝑐 limit, where 𝑁 𝑐 is the number of quark colors, we show that the NMEs receive contributions not only from the twist-3 PDF 𝑒⁡(𝑥) but also from an additional, previously neglected nucleon form factor. Incorporating constraints from the spin-flavor expansion, recent experimental data on 𝑒⁡(𝑥), as well as model calculations of 𝑒⁡(𝑥), we estimate the NMEs of the cMDM operators. Our analysis indicates that the NMEs are dominated by the nucleon form factors, and the cMDM contributions to pion-nucleon couplings can be comparable to those from the quark sigma terms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Reexamining circular dichroism in photoemission from a topological insulator

The orbital angular momentum (OAM) of electron states is an essential ingredient for topological and quantum geometric quantities in solids. For example, Dirac surface states with helical spin- and orbital-angular momenta are a hallmark of a 3D topological insulator. Angle-resolved photoemission spectroscopy (ARPES) with variable circular light polarization, known as circular dichroism (CD), has been assumed to be a direct probe of OAM and, by proxy, of the Berry curvature of electronic bands in energy- and momentum-space. Indeed, topological surface states have been shown to exhibit angle-dependent CD (CDAD), and more broadly, CD is often interpreted as evidence of spin-orbit coupling. Meanwhile, it is well-established that CD originates from the photoemission matrix elements, which can have extrinsic contributions related to the experimental geometry and the inherently broken inversion symmetry at the sample surface. Therefore, it is important to broadly examine CD-ARPES to determine the scenarios in which it provides a robust probe of intrinsic material physics. We performed CD-ARPES on the canonical topological insulator Bi 2 ⁢Se 3 over a wide range of incident photon energies. Not only do we observe angle-dependent CD in the surface states, as expected, but we also find CD of a similar magnitude in virtually all bulk bands. Since OAM is forbidden by inversion symmetry in the bulk, we conclude this originates from symmetry-breaking in the photoemission process. Comparison with theoretical calculations supports this view and suggests that “hidden” OAM—localized to atomic sites within each unit cell—contributes significantly. Additional effects, including inter-atomic interference and final-state resonances, are responsible for the rapid variation of the CDAD signal with photon energy.

Sidilkover, Ittai [Tel Aviv University, Tel Aviv (↗

DOE OSTI2026

Skewness-dependent moments of the pion GPD from nonlocal quark-bilinear correlators

We present lattice QCD calculations of the odd Mellin moments of pion valence-quark generalized parton distribution up to fifth order ⟨𝑥 4 ⟩ and for the skewness range [−0.33, 0] using operator product expansion of bilocal quark-bilinear operators. The calculations are performed on an ensemble with lattice spacing 𝑎 = 0.04 fm and valence pion mass 300 MeV, employing boosted pion states with momenta up to 2.428 GeV and momentum transfers reaching 2.748 GeV 2 . We employ ratio-scheme renormalization and next-to-leading logarithmic resummed perturbative matching. At zero skewness, our results are consistent with previous lattice studies. By combining matrix elements at multiple values of skewness and momentum transfer, skewness-dependent moments are obtained through simultaneous polynomiality-constrained fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Efficient many-jet event generation with flow matching

We apply for the first time, to the best of our knowledge, the flow matching method to the problem of phase-space sampling for event generation in high-energy collider physics. By training the model to remap the random numbers used to generate the momenta and helicities of the scattering matrix elements as implemented in the portable partonic event generator pepper, we find substantial efficiency improvements in the studied processes. We focus our study on the highest final-state multiplicities in Drell-Yan and top-antitop pair production used in simulated samples for the Large Hadron Collider, which computationally are the most relevant ones. We find that the unweighting efficiencies improve by factors of 184 and 25, respectively, when compared to the standard approach of using a vegas-based optimization. We also compare continuous normalizing flows trained with flow matching against the previously studied normalizing flows based on coupling layers and find that the former leads to better results, faster training and a better scaling behavior across the studied multiplicity range, while the latter evaluate faster. When combining the advantages of both methods using the regflow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest final-state multiplicities.

Bothmann, E. [CERN; Gottingen U.] (ORCID:000000016↗

DOE OSTI2025

Search for Majorana Neutrinos with the Complete KamLAND-Zen Dataset

We present a search for neutrinoless double-beta (0⁢𝜈⁢𝛽⁢𝛽) decay of 136 Xe using the full KamLAND-Zen 800 dataset with 745 kg of enriched xenon, corresponding to an exposure of 2.1 ton yr of 136 Xe. This updated search benefits from a more than twofold increase in exposure, recovery of photo-sensor gain, and reduced background from muon-induced spallation of xenon. Combining with the search in the previous KamLAND-Zen phase, we obtain a lower limit for the 0⁢𝜈⁢𝛽⁢𝛽 decay halflife of 𝑇$^{0⁢𝜈}_{1/2}$ >3.8 ×10 26 yr at 90% CL, a factor of 1.7 improvement over the previous limit. The corresponding upper limits on the effective Majorana neutrino mass are in the range 28–122 meV using phenomenological nuclear matrix element calculations.

Abe, S. [Tohoku University] (ORCID:000000022110513↗

DOE OSTI2026

Monte Carlo Event Generation with Continuous Normalizing Flows

We apply continuous normalizing flows trained with the flow matching method to the problem of phase-space sampling in Monte Carlo event generation for high-energy collider physics. Focusing on lepton-pair and top-quark pair production with multiple jets, the two computationally most expensive processes at the Large Hadron Collider, we train helicity-conditioned continuous normalizing flows to remap the random numbers used in matrix element evaluation. Compared to standard methods, we achieve unweighting efficiency improvements by factors of up to 184 and 25 for the two processes at their respective highest jet number, at the cost of an increased evaluation time. When combining the advantages of continuous normalizing flows with the fast evaluation times of coupling-layer-based flows, using the RegFlow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest jet numbers. These substantial gains highlight the promise of samplers based on machine learning for next-generation collider experiments.

Bothmann, Enrico [CERN; Gottingen U.] (ORCID:00000↗

DOE OSTI2025

Logical error rates for the surface code under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions

With fault-tolerant quantum computing (FTQC) on the horizon, it is critical to understand sources of logical errors in plausible hardware implementations of quantum error-correcting codes. Detailed error modeling of computational instructions on particular FTQC architectures will enable the better prediction of error propagation in FT-encoded quantum circuits while revealing where greater attention is needed in hardware design. In this work, we consider logical error rates for the surface code implemented on a hypothetical grid-based trapped-ion quantum charge-coupled device architecture. Specifically, we construct logical channels for the idling surface code and examine its diamond error under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions. We include the coherent dephasing noise that is known to accumulate during physical qubit idling and transport in these systems, determining idling and transport durations using the time-resolved output of an open-source trapped-ion surface code compiler. To estimate expectation values of logical Pauli observables following hardware circuits containing non-Clifford sources of noise, we utilize a Monte Carlo technique to sample from an underlying quasiprobability distribution of Clifford circuits that we independently simulate in a phase-sensitive fashion. We verify error suppression up to code distance 𝑑 = 11 at coherent dephasing rates near and below those of current-generation trapped-ion quantum computers and find that logical error rates align with those of analogous fully stochastic simulations in this regime. Exploring higher dephasing rates at 𝑑 = 3−5, we find evidence for growing coherent rotations about all three logical Pauli axes, increased diagonal logical error process matrix elements relative to those of stochastic simulations, and a reduced dephasing rate threshold. Overall, our work paves a way toward realistic hardware emulation of small fault-tolerant quantum processes, e.g., members of an FTQC instruction set.

Quantum benchmarking↗

DOE OSTI2025

Robust ab initio predictions for dimensionless ratios of 𝐸⁢2 and radius observables. II. Estimation of 𝐸⁢2 transition strengths by calibration to the charge radius

Converged results for 𝐸⁢2 observables are notoriously challenging to obtain in ab initio no-core configuration interaction approaches. Matrix elements of the 𝐸⁢2 operator are sensitive to the large-distance tails of the nuclear wave function, which converge slowly in an oscillator basis expansion. Similar convergence challenges beset ab initio prediction of the nuclear charge radius. However, we exploit systematic correlations between the calculated 𝐸⁢2 and radius observables to yield meaningful predictions for relations among these observables. In particular, we examine ab initio predictions for dimensionless ratios of the form 𝐵⁡(𝐸⁢2)/(𝑒 2 ⁢𝑟 4 ) for nuclei throughout the 𝑝 shell. Finally, meaningful predictions for 𝐸⁢2 transition strengths may then be made by calibrating to the ground-state charge radius if experimentally known.

ab initio calculations↗

DOE OSTI2026

Back-to-back dijet production in DIS at arbitrary Bjorken x: TMD gluon distributions to twist-3 accuracy

We derive the gluon transverse-momentum-dependent (TMD) operator structure of back-to-back\\\\r\\\\nquark–antiquark dijet production in deep inelastic scattering at arbitrary Bjorken-x to twist-3 ac\\\\r\\\\ncuracy. Working at leading order in the strong coupling and in the kinematic regime where the\\\\r\\\\ntransverse momentum imbalance of the jets is much smaller than their individual transverse mo\\\\r\\\\nmenta, we perform a systematic gradient expansion of the quark propagator in a background gluon\\\\r\\\\nfield. This expansion organizes multiple interactions with the target in terms of longitudinal Wilson\\\\r\\\\nlines and gauge-invariant field-strength insertions, yielding a TMD description valid beyond the\\\\r\\\\nstrict high-energy eikonal (x → 0) approximation. We obtain explicit cross sections for longitudi\\\\r\\\\nnally and transversely polarized virtual photons, identifying all contributing gluon TMD operators\\\\r\\\\nup to twist-3, including structures involving F+−, Fij, and three-gluon correlators. The full lon\\\\r\\\\ngitudinal phase eixP+z− associated with Bjorken-x is retained throughout. In the small-x limit,\\\\r\\\\nour results reproduce the known sub-eikonal expressions obtained in the Color Glass Condensate\\\\r\\\\nframework, establishing a direct connection between the general-x TMD expansion and high-energy\\\\r\\\\nfactorization. We further reduce the operator basis using equations of motion, minimizing the num\\\\r\\\\nber of independent nonperturbative matrix elements entering the cross section. This work provides\\\\r\\\\na systematic foundation for extending TMD analyses of dijet production beyond leading twist, es\\\\r\\\\ntablishing a unified operator framework valid at arbitrary Bjorken-x that smoothly interpolates\\\\r\\\\nbetween moderate- and small-x descriptions of gluon TMDs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Electron correlation effects in attosecond photoionization of CO 2

A technique for measuring photoionization time delays with attosecond precision is combined with calculations of photoionization matrix elements to demonstrate how multielectron dynamics affect photoionization time delays in carbon dioxide. Electron correlation is observed to affect the time delays through two mechanisms: autoionization of molecular Rydberg states and accelerated escape from a continuum shape resonance.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Orbital-resolved calculations of two-center interferences in linear triatomic molecules

In this report we perform ab initio calculations of high-harmonic spectroscopy (HHS) of two-center interference phenomena in oriented carbon-dichalcogen molecules, using time-dependent density functional theory (TDDFT). We show that by resolving the harmonic response into contributions from individual Kohn-Sham orbitals, we can extract target-specific characteristics for both the spectral amplitude and phase. We also discuss that this extraction is predicated on a careful analysis and normalization of the harmonic spectrum. Finally, we present field-free scattering calculations that mimic the recollision step in high-order-harmonic generation and allow the calculation of recombination dipole matrix elements without explicitly calculating the scattering states of a molecule. We show that the orbital-resolved TDDFT HHS results and results based on our field-free scattering calculations are in very good agreement with each other.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Quantum Krylov subspace algorithms for ground- and excited-state energy estimation

Quantum Krylov subspace diagonalization (QKSD) algorithms provide a low-cost alternative to the conventional quantum phase estimation algorithm for estimating the ground- and excited-state energies of a quantum many-body system. While QKSD algorithms typically rely on using the Hadamard test for estimating Krylov subspace matrix elements of the form $\langle \phi_i|e^{-\widehat{H}τ}|\phi_j\rangle$, the associated quantum circuits require an ancilla qubit with controlled multiqubit gates that can be quite costly for near-term quantum hardware. In this paper, we show that a wide class of Hamiltonians relevant to condensed-matter physics and quantum chemistry contain symmetries that can be exploited to avoid the use of the Hadamard test. We propose a multifidelity estimation protocol that can be used to compute such quantities, showing that our approach, when combined with efficient single-fidelity estimation protocols, provides a substantial reduction in circuit depth. In addition, here we develop a unified theory of quantum Krylov subspace algorithms and present three quantum-classical algorithms for the ground- and excited-state energy estimation problems, where each algorithm provides various advantages and disadvantages in terms of total number of calls to the quantum computer, gate depth, classical complexity, and stability of the generalized eigenvalue problem within the Krylov subspace.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Simultaneous stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Furthermore, solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Selective hybridization between the main band and the superstructure band in the Bi 2 Sr 2 CaCu 2 O 8 + δ superconductor

High-resolution laser-based angle-resolved photoemission measurements have been carried out on Bi 2 Sr 2 CaCu 2 O 8+δ (Bi2212) and Bi 2 Sr 2–x La x CuO 6+δ (Bi2201) superconductors. Unexpected hybridization between the main band and the superstructure band in Bi2212 is clearly revealed. In the momentum space where one main Fermi surface intersects with one superstructure Fermi surface, four bands are observed instead of two. The hybridization exists in both the superconducting state and the normal state, and in Bi2212 samples with different doping levels. Such a hybridization is not observed in Bi2201. This phenomenon can be understood by considering the bilayer splitting in Bi2212, the selective hybridization of two bands with peculiar combinations, and the altered matrix element effects of the hybridized bands. These observations indicate that the origin of the superstructure bands is intrinsic to the CuO 2 planes. Furthermore, understanding physical properties and superconductivity in Bi2212 should consider the complete Fermi surface topology which involves the main bands, the superstructure bands, and their interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Diffusive hydrodynamics from integrability breaking

We describe here the crossover from generalized to conventional hydrodynamics in nearly integrable systems. Integrable systems have infinitely many conserved quantities, which spread ballistically, in general. When integrability is broken, only a few of these conserved quantities survive. The remaining conserved quantities are generically transported diffusively; we further derive a compact and general diffusion equation for these. The diffusion constant depends on the matrix elements of the integrability-breaking perturbation; for a certain class of integrability-breaking perturbations, including long-range interactions, the diffusion constant can be expressed entirely in terms of generalized hydrodynamic data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Quantum dynamics in strongly driven random dipolar magnets

The random dipolar magnet LiHo x Y 1– x F 4 enters a strongly frustrated regime for small Ho 3+ concentrations with x < 0.05. In this regime, the magnetic moments of the Ho 3+ ions experience small quantum corrections to the common Ising approximation of LiHo x Y 1– x F 4 , which lead to a $\mathbb{Z}_2$-symmetry breaking and small, degeneracy breaking energy shifts between different eigenstates. Here we show that destructive interference between two almost degenerate excitation pathways burns spectral holes in the magnetic susceptibility of strongly driven magnetic moments in LiHo x Y 1– x F 4 . Such spectral holes in the susceptibility, microscopically described in terms of Fano resonances, can already occur in setups of only two or three frustrated moments, for which the driven level scheme has the paradigmatic Λ shape. For larger clusters of magnetic moments, the corresponding level schemes separate into almost isolated many-body Λ schemes, in the sense that either the transition matrix elements between them are negligibly small or the energy difference of the transitions is strongly off-resonant to the drive. This enables the observation of Fano resonances, caused by many-body quantum corrections to the common Ising approximation also in the thermodynamic limit. Finally, we discuss its dependence on the driving strength and frequency as well as the crucial role that is played by lattice dissipation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Perturbative instability of nonergodic phases in non-Abelian quantum chains

An important challenge in the field of many-body quantum dynamics is to identify nonergodic states of matter beyond many-body localization (MBL). Strongly disordered spin chains with non-Abelian symmetry and chains of non-Abelian anyons are natural candidates, as they are incompatible with standard MBL. In such chains, real space renormalization group methods predict a partially localized, nonergodic regime known as a quantum critical glass (a critical variant of MBL). Here, this regime features a treelike hierarchy of integrals of motion and symmetric eigenstates with entanglement entropy that scales as a logarithmically enhanced area law. We argue that such tentative nonergodic states are perturbatively unstable using an analytic computation of the scaling of off-diagonal matrix elements and accessible level spacing of local perturbations. Our results indicate that strongly disordered chains with non-Abelian symmetry display either spontaneous symmetry breaking or ergodic thermal behavior at long times. We identify the relevant length and timescales for thermalization: Even if such chains eventually thermalize, they can exhibit nonergodic dynamics up to parametrically long timescales with a nonanalytic dependence on disorder strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗