Search NASA⌕ Search

SEARCH · Search NASA

Results for “lattice simulations”

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 OSTI2024

Shedding light on charmonium

We investigate E1 radiative transitions within charmonium in a relativistic approach based on light-front QCD. In quantum field theory, two sets of processes are pure E1: χ c 0 → J / ψ γ ( ψ → χ c 0 γ ) and h c → η c γ ( η c ′ → h c γ ), both involving the P -wave charmonia. We compute the E1 radiative decay widths as well as the corresponding transition form factors of various processes including those involving 2 P states. These observables provide an access to the microscopic structures of the P -wave charmonium. We show that our parameter-free predictions are in excellent agreement with the experimental measurements as well as lattice simulations whenever available. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Stress inside the pion in holographic light-front QCD

In this work, we propose a method to compute the gravitational form factor D ( Q 2 ) in holographic QCD by exploiting the remarkable correspondence between semiclassical light-front QCD and semiclassical field theories in wrapped spacetime in five dimensions. The use of light-front holography bridges physics at large Q 2 as attained in light-front QCD and physics at small Q 2 where the coupling to scalar and tensor fields, e.g. glueballs, are dominant. As an application, we compute the D -term for the pion and compare the results with recent lattice simulations. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Mixing effects on spectroscopy and partonic observables of heavy mesons with logarithmic confining potential in a light-front quark model

Using the variational principle, we systematically investigate the mass spectra and wave functions of both 1⁢𝑆 and 2⁢𝑆 state heavy pseudoscalar (𝑃) and vector (𝑉) mesons within the light-front quark model. This approach incorporates a Coulomb plus logarithmic confinement potential to accurately describe the constituent quark and antiquark dynamics. Additionally, spin hyperfine interactions are introduced perturbatively to compute the masses of pseudoscalar and vector mesons. The present analyses of the 1⁢𝑆 and 2⁢𝑆 states require the consideration of mixing between them to account for empirical constraints. These constraints include the mass gap Δ⁢𝑀 𝑃 >Δ⁢𝑀 𝑉 , where Δ⁢𝑀 𝑃⁡(𝑉) =𝑀$^{2⁢𝑆}_{𝑃⁡(𝑉)}$−𝑀$^{1⁢𝑆}_{𝑃⁡(𝑉)}$ and the hierarchy of the decay constants 𝑓 1⁢𝑆 >𝑓 2⁢𝑆 . We find the optimal value of the mixing angle to be 𝜃 =1⁢8°, significantly enhancing the consistency between our spectroscopic predictions and the experimental data compiled by the Particle Data Group. Furthermore, based on the predicted mass, the newly observed resonance 𝐵 𝐽⁡ (5840) could be assigned as a 2 1⁢ 𝑆 0 state in the 𝐵 meson family. The study also reports various pertinent observables, including twist-two distribution amplitudes, electromagnetic form factors, charge radii, 𝜉 moments, and transition form factors that are found to be consistent with both available lattice simulations and experimental data. In addition, our predicted branching ratios for the channels of 𝐵 + →𝜏 + ⁢𝜈 𝜏 as well as rare decays of 𝐵 0 and 𝐵$^0_𝑠$ appear in accordance with experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Chiral Instabilities and the Onset of Chiral Turbulence in QED Plasmas

We present a first principles study of chiral plasma instabilities and the onset of chiral turbulence in QED plasmas with strong gauge matter interaction ( e 2 N f = 64 ), far from equilibrium. By performing classical-statistical lattice simulations of the microscopic theory, we show that the generation of strong helical magnetic fields from a helicity imbalance in the fermion sector proceeds via three distinct phases. During the initial linear instability regime the helicity imbalance of the fermion sector causes an exponential growth (damping) of magnetic field modes with right- (left-) handed polarization, for which we extract the characteristic growth (damping) rates. Secondary growth of unstable modes accelerates the helicity transfer from fermions to gauge fields and ultimately leads to the emergence of a self-similar scaling regime characteristic of a decaying turbulence, where magnetic helicity is efficiently transferred to macroscopic length scales. Within this turbulent regime, the evolution of magnetic helicity spectrum can be described by an infrared power spectrum with spectral exponent κ = 10.2 ± 0.5 and dynamical scaling exponents α = 1.14 ± 0.50 and β = 0.37 ± 0.13 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Conjectures about the Chiral Phase Transition in QCD from Anomalous Multi-Instanton Interactions

To date numerical simulations of lattice QCD have not found a chiral phase transition of first order that is expected to occur for sufficiently light pions. We show how the restoration of an exact global chiral symmetry can strongly decrease the breaking of the approximate, anomalous U A ( 1 ) symmetry. This is testable on the lattice through simulations for one through four flavors. In QCD a small breaking of the U A ( 1 ) symmetry in the chirally symmetric phase generates novel experimental signals. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Gravitational Wave Symphony from Oscillating Spectator Scalar Fields

We investigate a generic source of stochastic gravitational wave background due to the parametric resonance of oscillating scalar fields in the early Universe. By systematically analyzing benchmark models through lattice simulations and considering a wide range of parameters, we demonstrate that such a scenario can lead to detectable signals in gravitational wave detectors over a broad frequency range and potentially address the recent findings by pulsar timing array experiments. Furthermore, these models naturally yield ultralight dark matter candidates or dark radiation detectable by cosmic microwave background observatories. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Nuclear effective field theory: Status and perspectives

The nuclear physics landscape has been redesigned as a sequence of effective field theories (EFTs) connected to the Standard Model through symmetries and lattice simulations of Quantum Chromodynamics (QCD). EFTs in this sequence are expansions around different low-energy limits of QCD, each with its own characteristics, scales, and ranges of applicability regarding energy and number of nucleons. We review each of the three main nuclear EFTs---Chiral, Pionless, Halo/Cluster---highlighting their similarities, differences, and connections. Here, in doing so, we survey the structural properties and reactions of nuclei that have been derived from the \abinitio solution of the few- and many-body problem built upon EFT input.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Wires, Solitons and the Big Bang. Final report

The overarching aim of the project was to improve our understanding of inflation and its end. The following were the main goals of the project: 1. Develop a novel statistical framework to understand particle production during inflation and reheating in scenarios with many poorly constrained components and calculate potential observational signatures. The goal was to be achieved in the following steps (i) Develop the theoretical framework (ii) Apply the framework to inflation and reheating 2. Develop analytical and numerical tools to explore the end of inflation and the energy transfer from inflation to daughter fields during reheating. The goal was to be achieved in two steps: (i) Characterize the eq. of state at the end of inflation using non-perturbative techniques and lattice simulations. (ii) Develop a novel numerical algorithm for gauge fields in a cosmological setting.

42 ENGINEERING↗
DOE OSTI2023

ESR Dipole Power Supply Current Ripple and Noise Specifications

This note presents key findings for the ESR main magnet dipole power supplies (PS), where we find the current ripple specification to be close to or beyond the state-of-the-art. These specifications originate from beam-beam considerations, with the requirement to limit the ripple-induced hadron emittance growth to below 10%/hour. Beam dynamics that drive this PS ripple specification arise from the beam motions at the Interaction Point (IP). The frequency of the motions can be separated into "low", compared to the betatron frequency, and "high", i.e. around the betatron frequency and harmonics. In terms of the driving frequency, "low" implies ƒ<<ƒ 0 v x,y and "fast" means ƒ≈{ƒ 0 v x,y , ƒ 0 (1- v x,y , etc.}, where ƒ 0 =1/ T 0 =78.2 kHz is the revolution frequency, and v x,y are the fractional parts of the betatron tunes. Frequencies higher than ƒ 0 /2 will be folded back due to the particles sampling the field once per turn. To provide flexibility for future lattice adjustments and working point variations, we do not consider the tunes as fixed. Instead, we assume a certain margin and allow them to potentially fall within the range of 0.1< v x,y <0.5. In other words, the high-frequency region spans approximately from 8kHz to 40 kHz. Consequently,, we definte the dipole PS ripple in two distinct frequency ranges: the low-frequency range of [1-8000] Hz and the high-frequency range of [8-40] kHz. For the physics effects we analyzed in this note, there is no distinction between the ripple (which can be approximately reproduced in the frequency domain) and the random noise if both have some power within the frequency bandwidth of interest. Therefore, while we will use the term "ripple" for short, it should always be understood that we are referring to "ripple plus noise." Except for the lower end of the low-frequency range, the impact of the rippling PS current on the beam will be considerably reduced due to the eddy currents induced in the walls of the vacuum chamber. We will account for this effect in the PS ripple specifications to follow. The remaining sections of this note are structure as follows: Section II outlines the beam-beam physics requirements for the positional stability of the beam at the IP. Section III describes the anticipated shielding effect of the eddy currents induced in the vacuum chamber. Section IV derives the ripple requirement for the low-frequency range by propagating the closed orbit ripple resulting from the rippling dipoles to the IP (relevant lattice simulation results are summarized in the Appendix). In Section V, we present the analytical criterion for the ripple in the high-frequency range by considering resonant oscillations of the electron beam around a stable closed orbit near the betatron frequency. Finally, Section VI provides a summary of our findings and discusses related work.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Variations to the z-Expansion of the Form Factor Describing the Decay of B Mesons

We examine the decay rate of the particle decay B 0 → D -ℓ +νℓ using data collected from the Belle Collaboration. We studied three parameterizations of the form factor which describe the differential decay rate, the Caprini, Lellouch, and Neubert (CLN) parametrization, the Boyd, Grinstein, and Lebed (BGL) parametrization, and the Bourrely, Caprini, and Lellouch (BCL) parameterization. The form factor is a function of the hadronic recoil variable w, and each parameterization contains unique free parameters which are the focus of this work. We test the extrapolations of the form factor by fitting many different subsets of the low w data and then compare the prediction of the fit to the high w data using a χ 2 -metric. By only fitting the low w data we are able to examine the stability of extrapolations which will be informative for lattice simulations.

Simons, Daniel↗
DOE OSTI2022

Symmetry Breaking in an Extended-O(2) Model

Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we consider an extended-O(2) model that differs from the ordinary O(2) model by an explicit symmetry breaking term. Its coupling allows to smoothly interpolate between the O(2) model (zero coupling) and a $q$-state clock model (infinite coupling). In the latter case, a $q$-state clock model can also be defined for non-integer values of $q$. Thus, such a limit can also be considered as an analytic continuation of an ordinary $q$-state clock model to non-integer $q$. The phase diagram of the extended-O(2) model in the infinite coupling limit was established in our previous work, where it was shown that for non-integer $q$, there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we investigate the model at finite values of the coupling using Monte Carlo and tensor methods. The results may be relevant for configurable Rydberg-atom arrays.

Hostetler, Leon↗
DOE OSTI2022

Qubit Regularization and Qubit Embedding Algebras

Qubit regularization is a procedure to regularize the infinite dimensional local Hilbert space of bosonic fields to a finite dimensional one, which is a crucial step when trying to simulate lattice quantum field theories on a quantum computer. When the qubit-regularized lattice quantum fields preserve important symmetries of the original theory, qubit regularization naturally enforces certain algebraic structures on these quantum fields. We introduce the concept of qubit embedding algebras (QEAs) to characterize this algebraic structure associated with a qubit regularization scheme. We show a systematic procedure to derive QEAs for the O(N) lattice spin models and the SU(N) lattice gauge theories. While some of the QEAs we find were discovered earlier in the context of the D-theory approach, our method shows that QEAs are far richer. A more complete understanding of the QEAs could be helpful in recovering the fixed points of the desired quantum field theories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Stout Smearing on a Quantum Computer

Smearing of gauge-field configurations in lattice field theory improves the results of lattice simulations by suppressing high energy modes from correlation functions. In quantum simulations, high kinetic energy eigenstates are introduced when the time evolution operator is approximated such as Trotterization. While improved Trotter product formulae exist to reduce the errors, they have diminishing accuracy returns with respect to resource costs. Therefore having an algorithm that has fewer resources than an improved Trotter formula is desirable. In this work I develop a representation agnostic method for quantum smearing and show that it reduces the coupling to high energy modes in the discrete nonabelian gauge theory $D_4$

Gustafson, Erik J.↗
DOE OSTI2024

Investigating SU(3) with Nf=8 fundamental fermions at strong renormalized coupling

Lattice simulations have observed a novel strong coupling symmetric mass generation (SMG) phase for the SU(3) gauge system with $N_f=8$ fundamental fermions (represented by two sets of staggered fields) at very large renormalized coupling ($g^2_{GF} \gtrsim 25$). The results of Phys.Rev.D 106 (2022) 014513 suggest that the SMG phase is separated from the weak coupling, conformal phase by a continuous phase transition, implying that the SMG phase exists in the continuum limit. To scrutinize these findings, we are generating a set of large volume zero temperature ensembles using nHYP improved staggered fermions with additional Pauli-Villars fields to tame gauge field fluctuations. We consider the low-lying meson spectrum and verify the existence of the SMG phase. Based on a finite size scaling analysis we predict that the phase transition between the strong and weak coupling phases is likely governed by a merged fixed point that is ultraviolet in the strong coupling but infrared in the weak coupling side. This finding suggests that the SU(3) 8-flavor system sits at the opening of the conformal window

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Lattice Monte Carlo simulations with two impurity worldlines

Abstract We develop the impurity lattice Monte Carlo formalism for the case of two distinguishable impurities in a bath of polarized fermions. The majority particles are treated as explicit degrees of freedom, while the impurities are described by worldlines. The latter serve as localized auxiliary fields, which affect the majority particles. We apply the method to non-relativistic three-dimensional systems of two impurities and a number of majority particles where both the impurity–impurity interaction and the impurity–majority interaction have zero range. We consider the case of an attractive impurity–majority interaction, and we study the formation and disintegration of bound states as a function of the impurity–impurity interaction strength. We also discuss the potential applications of this formalism to other quantum many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Off-Lattice Markov Chain Monte Carlo Simulations of Mechanically Driven Polymers

Here, we develop off-lattice simulations of semiflexible polymer chains subjected to applied mechanical forces by using Markov Chain Monte Carlo. Our approach models the polymer as a chain of fixed length bonds, with configurations updated through adaptive nonlocal Monte Carlo moves. This proposed method enables precise calculation of a polymer’s response to a wide range of mechanical forces, which traditional on-lattice models cannot achieve. Our approach has shown excellent agreement with theoretical predictions of persistence length and end-to-end distance in quiescent states as well as stretching distances under tension. Moreover, our model eliminates the orientational bias present in on-lattice models, which significantly impacts calculations such as the scattering function, a crucial technique for revealing the polymer conformation.

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

Quantum error thresholds for gauge-redundant digitizations of lattice field theories

In the quantum simulation of lattice gauge theories, gauge symmetry can be either fixed or encoded as a redundancy of the Hilbert space. While gauge-fixing reduces the number of qubits, keeping the gauge redundancy can provide code space to mitigate and correct quantum errors by checking and restoring Gauss’s law. In this work, we consider the correctable errors for generic finite gauge groups and design the quantum circuits to detect and correct them. We calculate the error thresholds below which the gauge-redundant digitization with Gauss’s law error correction has better fidelity than the gauge-fixed digitization involving only gauge-invariant states. Our results provide guidance for fault-tolerant quantum simulations of lattice gauge theories. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗