Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fermions”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Form factors for semi-leptonic $B_{(s)} \to D_{(s)}^∗\ell \nu_\ell$ decays

Semileptonic $B_{(s)}$ decays are of great phenomenological interestbecause they allow to extract CKM matrix elements or test lepton flavouruniversality. Taking advantage of existing data, we explore extractingform factors for vector final states using the narrow widthapproximation. Based on RBC/UKQCD's set of 2+1 flavour gauge fieldensembles with Shamir domain-wall fermion and Iwasaki gauge fieldaction, we study semileptonic $B_{(s)}$ decays using domain-wallfermions for light, strange and charm quarks, whereas bottom quarks aresimulated with the relativistic heavy quark (RHQ) action. Exploratoryresults for $B_s \to D_s^* \ell \nu_\ell$ are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Simulating Effective QED on Quantum Computers

In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Eigenstate entanglement scaling for critical interacting spin chains

With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ − 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ = polylog ( ϵ − 1 ) .

Tong, Yu↗
DOE OSTI2023

Quantum Gauge Networks: A New Kind of Tensor Network

Although tensor networks are powerful tools for simulating low-dimensional quantum physics, tensor network algorithms are very computationally costly in higher spatial dimensions. We introduce quantum gauge networks: a different kind of tensor network ansatz for which the computation cost of simulations does not explicitly increase for larger spatial dimensions. We take inspiration from the gauge picture of quantum dynamics, which consists of a local wavefunction for each patch of space, with neighboring patches related by unitary connections. A quantum gauge network (QGN) has a similar structure, except the Hilbert space dimensions of the local wavefunctions and connections are truncated. We describe how a QGN can be obtained from a generic wavefunction or matrix product state (MPS). All 2k-point correlation functions of any wavefunction for M many operators can be encoded exactly by a QGN with bond dimension O(M k ). In comparison, for just k = 1, an exponentially larger bond dimension of 2 M/6 is generically required for an MPS of qubits. We provide a simple QGN algorithm for approximate simulations of quantum dynamics in any spatial dimension. The approximate dynamics can achieve exact energy conservation for time-independent Hamiltonians, and spatial symmetries can also be maintained exactly. We benchmark the algorithm by simulating the quantum quench of fermionic Hamiltonians in up to three spatial dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Scattering wave packets of hadrons in gauge theories: Preparation on a quantum computer

Quantum simulation holds promise of enabling a complete description of high-energy scattering processes rooted in gauge theories of the Standard Model. A first step in such simulations is preparation of interacting hadronic wave packets. To create the wave packets, one typically resorts to adiabatic evolution to bridge between wave packets in the free theory and those in the interacting theory, rendering the simulation resource intensive. In this work, we construct a wave-packet creation operator directly in the interacting theory to circumvent adiabatic evolution, taking advantage of resource-efficient schemes for ground-state preparation, such as variational quantum eigensolvers. By means of an ansatz for bound mesonic excitations in confining gauge theories, which is subsequently optimized using classical or quantum methods, we show that interacting mesonic wave packets can be created efficiently and accurately using digital quantum algorithms that we develop. Specifically, we obtain high-fidelity mesonic wave packets in the Z 2 and U(1) lattice gauge theories coupled to fermionic matter in 1+1 dimensions. Our method is applicable to both perturbative and non-perturbative regimes of couplings. The wave-packet creation circuit for the case of the Z 2 lattice gauge theory is built and implemented on the Quantinuum H1-1 trapped-ion quantum computer using 13 qubits and up to 308 entangling gates. The fidelities agree well with classical benchmark calculations after employing a simple symmetry-based noise-mitigation technique. This work serves as a step toward quantum computing scattering processes in quantum chromodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The 1 S 0 Pairing Gap in Neutron Matter

We report ab initio calculations of the S wave pairing gap in neutron matter calculated using realistic nuclear Hamiltonians that include two- and three-body interactions. We use a trial state, properly optimized to capture the essential pairing correlations, from which we extract ground state properties by means of auxiliary field diffusion Monte Carlo simulations. We extrapolate our results to the thermodynamic limit by studying the finite-size effects in the symmetry-restored projected Bardeen-Cooper-Schrieffer (PBCS) theory and compare our results to other ab initio studies done in the past. Our quantum Monte Carlo results for the pairing gap show a modest suppression with respect to the mean-field BCS values. These results can be connected to cold atom experiments, via the unitarity regime where fermionic superfluidity assumes a unified description, and they are important in the prediction of thermal properties and the cooling of neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Variational Quantum Chemistry Programs in JaqalPaq

We present example quantum chemistry programs written with JaqalPaq, a python meta-programming language used to code in Jaqal (Just Another Quantum Assembly Language). These JaqalPaq algorithms are intended to be run on the Quantum Scientific Computing Open User Testbed (QSCOUT) platform at Sandia National Laboratories. Our exemplars use the variational quantum eigensolver (VQE) quantum algorithm to compute the ground state energies of the H2, HeH+, and LiH molecules. Since the exemplars focus on how to program in JaqalPaq, the calculations of the second-quantized Hamiltonians are performed with the PySCF python package, and the mappings of the fermions to qubits are obtained from the OpenFermion python package. Using the emulator functionality of JaqalPaq, we emulate how these exemplars would be executed on an error-free QSCOUT platform and compare the emulated computation of the bond-dissociation curves for these molecules with their exact forms within the relevant basis.

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

Bilocal Field Theory for Composite Scalar Bosons

We give a bilocal field theory description of a composite scalar with an extended binding potential that reduces to the Nambu–Jona-Lasinio (NJL) model in the pointlike limit. This provides a description of the internal dynamics of the bound state and features a static internal wave function, $\phi(\vec{r})$, in the center-of-mass frame that satisfies a Schrödinger–Klein–Gordon equation with eigenvalues m 2 . We analyze the “coloron” model (single perturbative massive gluon exchange) which yields a UV completion of the NJL model. This has a BCS-like enhancement of its interaction, ∝N c the number of colors, and is classically critical with g critical remarkably close to the NJL quantum critical coupling. Negative eigenvalues for m 2 lead to spontaneous symmetry breaking, and the Yukawa coupling of the bound state to constituent fermions is emergent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Flipped Quartification: Product Group Unification with Leptoquarks

The quartification model is an SU(3) 4 extension with a bi-fundamental fermion sector of the well-known SU(3) 3 bi-fundamentalfication model. An alternative “flipped” version of the quartification model is obtained by rearrangement of the particle assignments. The flipped model has two standard (bi-fundamentalfication) families and one flipped quartification family. In contrast to traditional product group unification models, flipped quartification stands out by featuring leptoquarks and thus allows for new mechanisms to explain the generation of neutrino masses and possible hints of lepton-flavor non-universality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Analysis of Near-Field Magnetic Responses on ZrTe 5 through Cryogenic Magneto-THz Nano-Imaging

One manifestation of light-Weyl fermion interaction is the emergence of chiral magnetic effects under magnetic fields. Probing real space magnetic responses at terahertz (THz) scales is challenging but highly desired, as the local responses are less affected by the topologically trivial inhomogeneity that is ubiquitous in spatially averaged measurements. Here, we implement a cryogenic THz microscopy instrument under a magnetic field environment—a task only recently achieved. We explore the technical approach of this system and characterize the magnetic field’s influence on our AFM operation by statistical noise analysis. We find evidence for local near-field spatial variations in the topological semimetal ZrTe 5 up to a 5-Tesla magnetic field and obtain near-field THz spectra to discuss their implications for future studies on the chiral magnetic effect.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure of Odd-A Ag Isotopes Studied via Algebraic Approaches

The structure of the odd-A silver isotopes Ag 103–115 is discussed within the frame of the interacting boson–fermion model (IBFM). An overview of their key properties is presented, with a particular attention paid to the “ J -1 anomaly”, represented by an abnormal ordering of the lowest 7/2 + and 9/2 + states. By examining previously published data and newly performed calculations, it is demonstrated that the experimentally known level schemes and electromagnetic properties of Ag 103–115 can be reproduced well within IBFM-1 by using a consistent set of model parameters. The contribution of different single-particle orbitals to the structure of the lowest-lying excited nuclear states in Ag 103–115 is discussed. Given that the J-1 anomaly brings down the 7/2 + level from the j −3 multiplet to energies, which can be thermally populated in hot stellar environments, the importance of low-lying excited states in odd-A silver isotopes for astrophysical processes is outlined.

IBFM-1↗
DOE OSTI2025

Illuminating Dark Matter Admixed in Neutron Stars with Simultaneous Mass–Radius Constraints

We investigate how simultaneous mass and radius measurements of massive neutron stars can help constrain the properties of dark matter possibly admixed in them. Within a fermionic dark matter model that interacts only through gravitation, along with a well-constrained nuclear matter equation of state, we show that the simultaneous mass and radius measurement of PSRJ0740+6620 reduces the uncertainty of dark matter central energy density by more than 50% compared to the results obtained from using the two observables independently, while other dark matter parameters remain unconstrained. Additionally, we find that the dark matter fraction fD should be smaller than 2% when constrained by the observed neutron star maximum mass alone, and it could be even smaller than 0.3% with the simultaneous measurement of mass and radius, supporting the conclusion that only a small amount of dark matter exists in dark matter admixed neutron stars (DANSs).

NICER↗
DOE OSTI2021

General Quantum Field Theory of Flavor Mixing and Oscillations

We review the canonical transformation in quantum physics known as the Bogoliubov transformation and present its application to the general theory of quantum field mixing and oscillations with an arbitrary number of mixed particles with either boson or fermion statistics. The mixing relations for quantum states are derived directly from the definition of mixing for quantum fields and the unitary inequivalence of the Fock space of energy and flavor eigenstates is shown by a straightforward algebraic method. The time dynamics of the interacting fields is then explicitly solved and the flavor oscillation formulas are derived in a unified general formulation with emphasis on antiparticle content and effect introduced by nontrivial flavor vacuum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

Abstract The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Shawqi, Shafayat (ORCID:0000000210956183)↗
DOE OSTI2022

Diophantine equations with sum of cubes and cube of sum

We solve Diophantine equations of the type $a(x^3+y^3+z^3)=(x+y+z)^3$, where $x$, $y$, $z$ are integer variables, and the coefficient $a \neq 0$ is rational. We show that there are infinite families of such equations, including those where $a$ is any cube or certain rational fractions, that have nontrivial solutions. There are also infinite families of equations that do not have any nontrivial solution, including those where $1/a=1-24/m$ with restrictions on the integer $m$. The equations can be represented by elliptic curves unless $a=9$ or $1$, and any elliptic curve of nonzero $j$-invariant and torsion group $\mathbb{Z}/3k\mathbb{Z}$ for $k=2,3,4$, or $\mathbb{Z}/2\mathbb{Z} \times \mathbb{Z}/6\mathbb{Z}$ corresponds to a particular $a$. We prove that for any $a$ the number of nontrivial solutions is at most $3$ or is infinite, and for integer $a$ it is either $0$ or $\infty$. For $a=9$, we find the general solution, which depends on two integer parameters. These cubic equations are important in particle physics, because they determine the fermion charges under the $U(1)$ gauge group.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The lattice extraction of the TMD soft function using the auxiliary field representation of the Wilson line

The transverse momentum dependent (TMD) soft function can be obtained by formulating the Wilson line in terms of auxiliary 1-dimensional fermion fields on the lattice. In this formulation, the directional vector of the auxiliary field in Euclidean space has the form ˜ 𝑛 = ¹𝑖𝑛0– ®0?– 𝑛3º, where the time component is purely imaginary. The components of these complex directional vectors in the Euclidean space can be mapped directly to the rapidities of the Minkowski space soft function. We present the results of the one-loop calculation of the Euclidean space analog to the soft function using these complex directional vectors. As a result, we show that the calculation is valid only when the directional vectors obey the relation: j𝑟 j = j𝑛3𝑛0j ¡ 1, and that this result corresponds to a computation in Minkowski space with space-like directed Wilson lines. Finally, we show that a lattice calculable object can be constructed that has the desired properties of the soft function.

Francis, Anthony↗
DOE OSTI2025

Surrogate Constructed Scalable Circuits ADAPT-VQE in the Schwinger model

Inspired by recent advancements of simulating periodic systems on quantum computers, we develop a new approach, (SC)$^2$-ADAPT-VQE, to further advance the simulation of these systems. Our approach extends the scalable circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable ``Surrogate Constructed'' method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC)$^2$ methodology. Our approach allows for a wider set of classical computations, on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute properties of the Schwinger model - quantum electrodynamics for a single, massive fermion in $1+1$ dimensions - and show that our method can be used to accurately extrapolate to the continuum limit.

Gustafson, Erik [RIACS, Mtn. View] (ORCID:00000001↗