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

Two-point sum-rules in three-dimensional Yang-Mills theory

We compute the stress-tensor two-point function in three-dimensional Yang-Mills theory to three-loops in perturbation theory. Using its calculable shape at high momenta, we test the notion that its Borel transform is saturated at low energies by the lowest glueball state(s). This assumption provides relatively stable estimates for the mass of the lightest glueball that we compare with lattice simulations. We also provide estimates for the coupling of the lightest glueball to the stress tensor. Along the way, we comment on the extent that such estimates are non-rigorous. Lastly, we discuss the possibility of applying the sum-rule analysis to two-point functions of higher-spin operators and obtain a crude approximation for the glueball couplings to these operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A lattice study of ππ scattering at large N c

We present the first lattice study of pion-pion scattering with varying number of colors, N c . We use lattice simulations with four degenerate quark flavors, N f = 4, and N c = 3 – 6. We focus on two scattering channels that do not involve vacuum diagrams. These correspond to two irreducible representations of the SU(4) flavor group: the fully symmetric one, SS, and the fully antisymmetric one, AA. The former is a repulsive channel equivalent to the isospin-2 channel of SU(2). By contrast, the latter is attractive and only exists for N f ≥ 4. A representative state is (|D$^{+}_{s}$π + > – | D + K + >)/√2. Using Lüscher’s formalism, we extract the near-threshold scattering amplitude and we match our results to Chiral Perturbation Theory (ChPT) at large N c . For this, we compute the analytical U(N f ) ChPT prediction for two-pion scattering, and use the lattice results to constrain the N c scaling of the relevant low-energy couplings.

1/N Expansion↗
DOE OSTI2023

Mapping out the thermodynamic stability of a QCD equation of state with a critical point using active learning

The Beam Energy Scan Theory (BEST) collaboration's equation of state (EoS) incorporates a three-dimensional Ising model critical point into the quantum chromodynamics (QCD) equation of state from lattice simulations. Furthermore, it contains four free parameters related to the size and location of the critical region in the QCD phase diagram. Certain combinations of the free parameters lead to acausal or unstable realizations of the EoS that should not be considered. In this work, we use an active learning framework to rule out pathological EoS efficiently. We find that checking stability and causality for a small portion of the parameters' range is sufficient to construct algorithms that perform with > 96 % accuracy across the entire parameter space. Though in this work we focus on a specific case, our approach can be generalized to any EoS containing a parameter space-class correspondence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Symmetry breaking in an extended O(2) model

Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we study an extended-O(2) model in two dimensions that differs from the ordinary O(2) model by the addition of an explicit symmetry breaking term − h q cos ( q φ ) . Its coupling h q allows to smoothly interpolate between the O(2) model ( h q = 0 ) and a q -state clock model ( h q → ∞ ). In the latter case, a q -state clock model can also be defined for noninteger values of q . Thus, such a limit can also be considered as an analytic continuation of an ordinary q -state clock model to noninteger q . In previous work, we established the phase diagram for noninteger q in the infinite coupling limit ( h q → ∞ ). We showed that there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we seek to establish the phase diagram at finite values of the coupling using Monte Carlo and tensor methods. We show that for noninteger q , the second-order phase transition at low temperature and crossover at high temperature persist to finite coupling. For integer q = 2 , 3, 4, we know there is a second-order phase transition at infinite coupling (i.e. the well-known clock models). At finite coupling, we find that the critical exponents for q = 3 , 4 vary with the coupling, and for q = 4 the transition may turn into a Berezinskii-Kosterlitz-Thouless transition at small coupling. We comment on the similarities and differences of the phase diagrams with those of quantum simulators of the Abelian-Higgs model based on ladder-shaped arrays of Rydberg atoms. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2023

Infrared fixed point of the SU(3) gauge theory with 𝑁 𝑓 = 10 flavors

We use lattice simulations and the continuous renormalization group method, based on the gradient flow, to calculate the 𝛽 function and anomalous dimensions of the SU(3) gauge theory with 𝑁 𝑓 = 10 flavors of fermions in the fundamental representation. We employ several improvements to extend the range of available renormalized couplings, including the addition of heavy Pauli-Villars bosons to reduce cutoff effects and the combination of a range of gradient flow transformations. While in the weak coupling regime our result is consistent with those of earlier studies, our techniques allow us to study the system at much stronger couplings than previously possible. We find that the renormalization group 𝛽 function develops a zero, corresponding to an infrared-stable fixed point, at gradient-flow coupling 𝑔 2 = 15.0⁢(5). We also determine the mass and tensor anomalous dimensions: At the fixed point we find 𝛾 𝑚 ≃0.6, suggesting that this system might be deep inside the conformal window.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Stochastic renormalization group and gradient flow

A non-perturbative and continuous definition of RG transformations as stochastic processes is proposed, inspired by the observation that the functional RG equations for effective Boltzmann factors may be interpreted as Fokker-Planck equations. The result implies a new approach to Monte Carlo RG that is amenable to lattice simulation. Long-distance correlations of the effective theory are shown to approach gradient-flowed correlations, which are simpler to measure. The Markov property of the stochastic RG transformation implies an RG scaling formula which allows for the measurement of anomalous dimensions when transcribed into gradient flow expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Stability of superconducting strings

We investigate the stability of superconducting strings as bound states of strings and fermion zero modes at both the classical and quantum levels. The dynamics of these superconducting strings can result in a stable configuration, known as a vorton . We mainly focus on global strings, but the majority of the discussion can be applied to local strings. Using lattice simulations, we study the classical dynamics of superconducting strings and confirm that they relax to the vorton configuration through Nambu-Goldstone boson radiation, with no evidence of over-shooting that would destabilize the vorton. We explore the tunneling of fermion zero modes out of the strings. Both our classical analysis and quantum calculations yield consistent results: the maximum energy of the zero mode significantly exceeds the fermion mass, in contrast to previous literature. Additionally, we introduce a world-sheet formalism to evaluate the decay rate of zero modes into other particles, which constitute the dominant decay channel. We also identify additional processes that trigger zero-mode decay due to non-adiabatic changes of the string configuration. In these decay processes, the rates are suppressed by the curvature of string loops, with exponential suppression for large masses of the final states. We further study the scattering with light charged particles surrounding the string core produced by the zero-mode current and find that a wide zero-mode wavefunction can enhance vorton stability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The electroweak phase transition: a collider target

Determining the thermal history of electroweak symmetry breaking (EWSB) is an important challenge for particle physics and cosmology. Lattice simulations indicate that EWSB in the Standard Model (SM) occurs through a crossover transition, while the presence of new physics beyond the SM could alter this thermal history. The occurrence of a first order EWSB transition would be particularly interesting, providing the needed pre-conditions for generation of the cosmic matter-antimatter asymmetry and sources for potentially observable gravitational radiation. I provide simple, generic arguments that if such an alternate thermal history exists, the new particles involved cannot be too heavy with respect to the SM electroweak temperature, nor can they interact too feebly with the SM Higgs boson. These arguments do not rely on the decoupling limit. I derive corresponding quantitative expectations for masses and interaction strengths which imply that their effects could in principle be observed (or ruled out) by prospective next generation high energy colliders. The simple, generic arguments provide a quantitative, parametric understanding of results obtained in a wide range of explicit model studies; relate them explicitly to the electroweak temperature; and delineate broad contours of collider phenomenology pertaining to a non-standard history of EWSB.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Chirality transfer and chiral turbulence in gauge theories

Chirality transfer between fermions and gauge fields plays a crucial role for understanding the dynamics of anomalous transport phenomena such as the Chiral Magnetic Effect. In this proceeding we present a first principles study of these processes based on classical-statistical real-time lattice simulations of strongly coupled QED (e 2 N f = 64). Furthermore, our simulations demonstrate that a chirality imbalance in the fermion sector triggers chiral plasma instabilities in the gaugefield sector, which ultimately lead to the generation of long range helical magnetic fields via an self-similar turbulent cascade of the magnetic helicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Non-cancellation of the parity anomaly in the strong-field regime of QED 2+1

Quantum fluctuations lead to an anomalous violation of parity symmetry in quantum electrodynamics for an even number of spatial dimensions. While the leading parity-odd electric current vanishes in vacuum, we uncover a non-cancellation of the anomaly for strong electric fields with distinct macroscopic signatures. We perform real-time lattice simulations with fully dynamical gauge fields and Wilson fermions in 2 + 1 space-time dimensions. In the static field limit, relevant at early times, we solve the problem analytically. Our results point out the fundamental role of quantum anomalies for strong-field phenomena, relevant for a wide range of condensed matter and high-energy applications, but also for the next generation of gauge theory quantum simulators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Perturbatively including inhomogeneities in axion inflation

Axion inflation, i.e. an axion-like inflaton coupled to an Abelian gauge field through a Chern-Simons interaction, comes with a rich and testable phenomenology. This is particularly true in the strong backreaction regime, where the gauge field production heavily impacts the axion dynamics. Lattice simulations have recently demonstrated the importance of accounting for inhomogeneities of the axion field in this regime. We propose a perturbative scheme to account for these inhomogeneities while maintaining high computational efficiency. Our goal is to accurately capture deviations from the homogeneous axion field approximation within the perturbative regime as well as self-consistently determine the onset of the nonperturbative regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

PT symmetry, pattern formation, and finite-density QCD

A longstanding issue in the study of quantum chromodynamics (QCD) is its behavior at nonzero baryon density, which has implications for many areas of physics. The path integral has a complex integrand when the quark chemical potential is nonzero and therefore has a sign problem, but it also has a generalized ${\mathscr{PT}}$ symmetry. We review some new approaches to ${\mathscr{PT}}$-symmetric field theories, including both analytical techniques and methods for lattice simulation. We show that ${\mathscr{PT}}$-symmetric field theories with more than one field generally have a much richer phase structure than their Hermitian counterparts, including stable phases with patterning behavior. The case of a ${\mathscr{PT}}$-symmetric extension of a φ 4 model is explained in detail. The relevance of these results to finite density QCD is explained, and we show that a simple model of finite density QCD exhibits a patterned phase in its critical region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nonperturbative and perturbative dynamics of a light QCD axion: Dark matter and the strong 𝐶⁢𝑃 problem

Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶⁢𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶⁢𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.

Axions↗
DOE OSTI2025

Bayesian Inferring Nucleon Gravitational Form Factors via Near-Threshold 𝐽/𝜓 Photoproduction

With Bayesian inference, we investigate the impact of recent near-threshold 𝐽/𝜓 production measurements by the 𝐽/𝜓⁢007 experiment and GlueX Collaboration on the extraction of proton gravitational form factors. We apply the generalized parton distribution framework at the next to leading order and demonstrate a stable expansion for the near-threshold kinematics. We find that the experimental constraints are in good agreement with the state-of-the-art lattice simulations, where negative 𝐶 𝑞 ⁡(𝑡) and 𝐶 𝑔 ⁡(𝑡) are strongly preferred. This highlights a great potential to extract them from future high-precision experiments.

QCD phenomenology↗
DOE OSTI2025

Hidden-strangeness tetraquarks in the dynamical diquark model

The dynamical diquark model describes multiquark exotic hadrons in terms of diquark components nucleated by heavy quarks and successfully explains multiple features of hidden-charm and hidden-bottom exotics. Here, we apply the model to the marginally heavy case of hidden-strange states to probe whether mesons near 2 GeV with peculiar properties, such as 𝜙⁡(2170), 𝑓 2 ⁡(2340), and 𝑋⁡(2370), are possible tetraquark candidates. We calculate spin-multiplet average masses using potentials obtained through lattice simulations and quark models, and we also describe the detailed spectra of the expected multiplets as a diagnostic to discern the nature of future hadrons likely to be discovered in this mass region by experiments at facilities such as BESIII, JLab, and the EIC.

Bound states↗
DOE OSTI2021

O 16 O 16 collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider comparing α clustering versus substructure

Collisions of . light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed 16 O 16 O run at the CERN Large Hadron Collider (LHC) and at the BNL Relativistic Heavy Ion Collider (RHIC) we study the potential for finding α clustering in 16 O. Here we use the state-of-the-art iEBE-VISHNU package with 16 O nucleonic configurations from ab initio nuclear lattice simulations. This setup was tuned using a Bayesian analysis on p Pb and PbPb systems. We find that the 16 O 16 O system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between α-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Minimal entanglement and emergent symmetries in low-energy QCD

Here, we study low-energy scattering of spin-$\frac{1}{2}$ baryons from the perspective of quantum information science, focusing on the correlation between entanglement minimization and the appearance of accidental symmetries. The baryon transforms as an octet under the SU⁡(3) flavor symmetry and its interactions below the pion threshold are described by contact operators in an effective field theory (EFT) of QCD. Despite there being 64 channels in the 2-to-2 scattering, only six independent operators in the EFT are predicted by SU⁡(3). We show that successive entanglement minimization in SU⁡(3)-symmetric channels are correlated with increasingly large emergent symmetries in the EFT. In particular, we identify scattering channels whose entanglement suppression are indicative of emergent SU⁡(6), SO⁡(8), SU⁡(8), and SU⁡(16) symmetries. We also observe the appearance of non-relativistic conformal invariance in channels with unnaturally large scattering lengths. Improved precision from lattice simulations could help determine the degree of entanglement suppression, and consequently the amount of accidental symmetry, in low-energy QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Time scales for nonlinear processes in preheating after multifield inflation with nonminimal couplings

We have conducted extensive lattice simulations to study the postinflation dynamics of multifield models involving nonminimal couplings. In this work, we explore the parameter dependence of preheating in these models and describe the various time scales that control such nonlinear processes as energy transfer, rescattering, and the approach to radiation domination and thermalization. In the limit of large nonminimal couplings ($ξ_I$ ~ 100), we find that efficient transfer of energy from the inflaton condensate to radiative degrees of freedom, emergence of a radiation-dominated equation of state, and the onset of thermalization each consistently occur within $N_{\text{reh}} ≲ 3 e$-folds after the end of inflation, largely independent of the values of the other couplings in the models. The exception is the case of negative ellipticity, in which there is a misalignment between the dominant direction in field space along which the system evolves and the larger of the nonminimal couplings $ξ_I$. In those cases, the field-space-driven parametric resonance is effectively shut off. More generally, the competition between the scalar fields’ potential and the field-space manifold structure can yield interesting phenomena such as two-stage resonances. Across many regions of parameter space, we find efficient re-scattering between the distinct fields, leading to a partial memory loss of the shape of the initial fluctuation spectrum. Despite the explosive particle production, which can lead to a quick depletion of the background energy density, the nonlinear processes do not induce any superhorizon correlations after the end of inflation in these models, which keeps predictions for cosmic microwave background observables unaffected by the late-time amplification of isocurvature fluctuations. Hence the excellent agreement between primordial observables and recent observations is preserved for this class of models, even when we consider postinflation dynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗