Real-time lattice gauge theory actions: Unitarity, convergence, and path integral contour deformations
Not Available
SEARCH · Search NASA
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.
Not Available
The statistical models used to derive the results of experimental analyses are of incredible scientific value and are essential information for analysis preservation and reuse. In this paper, we make the scientific case for systematically publishing the full statistical models and discuss the technical developments that make this practical. By means of a variety of physics cases - including parton distribution functions, Higgs boson measurements, effective field theory interpretations, direct searches for new physics, heavy flavor physics, direct dark matter detection, world averages, and beyond the Standard Model global fits - we illustrate how detailed information on the statistical modelling can enhance the short- and long-term impact of experimental results.
We study a 2+1D lattice gauge theory with fundamental representation scalar fields which has both Higgs and confining regimes with a spontaneously-broken U(1) 0-form symmetry. We show that the Higgs and confining regimes may be distinguished by a natural gauge invariant observable: the phase Ω of a correlation function of a vortex line operator linking with an electric Wilson line. We employ dualities and strong coupling expansions to analytically explore parameter regimes which were inaccessible in previous continuum calculations, and discuss possible implications for the phase diagram.
This paper presents an application to metal matrix composites (MMCs) of an enhanced elasto-viscoplastic Fast Fourier Transform (EVP-FFT) formulation coupled with a phenomenological continuum Mesoscale Field Dislocation Mechanics (MFDM) theory. Contrary to conventional crystal plasticity, which only accounts for plastic flow and hardening induced by statistically stored dislocations (SSDs), MFDM-EVP-FFT also describes the evolution of polarized geometrically necessary dislocation (GND) density and its effect on both plastic flow and hardening. Numerical results for a Fe–TiB 2 MMC made of a ferrite matrix (α-Fe) and elastic ceramic particles (TiB 2 ) are presented. Full-field simulations are performed using synthetic periodic unit cells representative of the MMC, with single-crystalline and polycrystalline matrix, for different particle interspacing distances. A strong dependence of the predicted equivalent stress, cumulated plastic strain and GND density fields with particle interspacing distance is observed, in contrast with conventional crystal plasticity. Correlations between these mechanical fields and microstructural features, and their influence on local and global mechanical behavior are examined for the different MMC microstructures.
Here we present the Aemulusν simulations: a suite of 150 (1.05 h -1 Gpc) 3 N-body simulations with a mass resolution of 3.51 × 10 10 Ω cb /0.3 h -1 M ⊙ in a wνCDM cosmological parameter space. The simulations have been explicitly designed to span a broad range in σ8 to facilitate investigations of tension between large scale structure and cosmic microwave background cosmological probes. Neutrinos are treated as a second particle species to ensure accuracy to 0.5 eV, the maximum neutrino mass that we have simulated. By employing Zel'dovich control variates, we increase the effective volume of our simulations by factors of 10-10 5 depending on the statistic in question. As a first application of these simulations, we build new hybrid effective field theory and matter power spectrum surrogate models, demonstrating that they achieve ≤ 1% accuracy for k ≤ 1 hMpc -1 and 0 ≤ z ≤ 3, and ≤ 2% accuracy for k ≤ 4 hMpc -1 for the matter power spectrum. We publicly release the trained surrogate models, and estimates of the surrogate model errors in the hope that they will be broadly applicable to a range of cosmological analyses for many years to come.
ABSTRACT The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as ‘wide-angle’ effects in redshift-space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work, we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in autocorrelation and in cross-correlation with the density field. Using both linear theory and the Zel'dovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low-redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration space, Fourier space, and spherical analyses; providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.
Production cross sections of heavy neutron-rich isotopes are calculated by employing quantal transport description in 250 Cf + 232 Th collisions. This quantal transport description is based on the stochastic mean-field approach, and it provides a microscopic approach beyond time-dependent Hartree-Fock theory to include mean-field fluctuations. Deexcitation of primary fragments is determined by employing the statistical GEMINI++ code. Here, calculations provide predictions for production cross sections of neutron rich transfermium isotopes without any adjustable parameters.
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 .
ABSTRACT Perturbation theory (PT) is often used to model statistical observables capturing the translation and rotation-invariant information in cosmological density fields. PT produces higher order corrections by integration over linear statistics of the density fields weighted by kernels resulting from recursive solution of the fluid equations. These integrals quickly become high dimensional and naively require increasing computational resources the higher the order of the corrections. Here, we show how to decouple the integrands that often produce this issue, enabling PT corrections to be computed as a sum of products of independent 1D integrals. Our approach is related to a commonly used method for calculating multiloop Feynman integrals in quantum field theory, the Gegenbauer Polynomial x-Space Technique. We explicitly reduce the three terms entering the 2-loop power spectrum, formally requiring 9D integrations, to sums over successive 1D radial integrals. These 1D integrals can further be performed as convolutions, rendering the scaling of this method Nglog Ng with Ng the number of grid points used for each fast Fourier transform. This method should be highly enabling for upcoming large-scale structure redshift surveys where model predictions at an enormous number of cosmological parameter combinations will be required by Monte Carlo Markov Chain searches for the best-fitting values.
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.
Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.
Model averaging is a useful and robust method for dealing with model uncertainty in statistical analysis. Often, it is useful to consider data subset selection at the same time, in which model selection criteria are used to compare models across different subsets of the data. Two different criteria have been proposed in the literature for how the data subsets should be weighted. We compare the two criteria closely in a unified treatment based on the Kullback-Leibler divergence and conclude that one of them is subtly flawed and will tend to yield larger uncertainties due to loss of information. Here, analytical and numerical examples are provided.
The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.
Excited-state effects lead to hard-to-quantify systematic uncertainties in lattice quantum chromodynamics (LQCD) spectroscopy calculations when computationally accessible imaginary times are smaller than inverse excitation gaps, as often arises for multi-hadron systems with signal-to-noise problems. Lanczos residual bounds address this by providing two-sided constraints on energies that do not require assumptions beyond Hermiticity, but often give very conservative systematic uncertainty estimates. Here, a more-constraining set of gap bounds is introduced for hadron spectroscopy. These bounds provide tighter constraints whose validity requires an explicit assumption about an energy gap. Exactly solvable lattice field theory correlators are used to test the utility of residual and gap bounds at finite and infinite statistics. Two-sided bounds and other analysis methods are then applied to a high-statistics LQCD calculation of nucleon-nucleon scattering at $m_π\sim 800$ MeV. Generalized eigenvalue problem (GEVP) and Lanczos energy estimators are compatible when applied to the same correlator data, but analyses including different interpolating operators show statistically significant inconsistencies. However, two-sided bounds from all operators are consistent. Under the assumption that the number of energy levels below $NΔ$ and $ΔΔ$ thresholds is the same as for non-interacting nucleons, gap bounds are sufficient to constrain nucleon-nucleon scattering amplitudes at phenomenologically relevant precision. Lanczos methods further reveal that energy-eigenstate estimates from previously studied asymmetric correlators have not converged over accessible imaginary times. Nevertheless, data-driven examples demonstrate why assumptions are required to draw conclusions about the natures of two-nucleon ground states at these masses.
We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory (χ EFT) that are provided by bound-state observables in the A = 3 and A = 4 sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the χ EFT expansion at next-to-next-to-leading order. A consistent solution for the 3 H binding energy, the 4 He binding energy and radius, and the 3 H β-decay rate can only be obtained if χ EFT truncation errors are included in the analysis. Here, the β-decay rate is the only one of these that yields a nondegenerate constraint on the 3NF low-energy constants, which makes it crucial for the parameter estimation. We use eigenvector continuation for fast and accurate emulation of no-core shell model calculations of the few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the parameters that quantify the truncation error. We find a χ EFT expansion parameter of Q = 0.33 ± 0.06 for these observables.
We analyze the robustness of the DESI 2024 cosmological inference from the full shape of the galaxy power spectrum to uncertainties in the Halo Occupation Distribution (HOD) model of the galaxy-halo connection and the choice of priors on nuisance parameters. We assess variations in the recovered cosmological parameters across a range of mocks populated with different HOD models and find that shifts are often greater than 20% of the expected statistical uncertainties from the DESI data. We encapsulate the effect of such shifts in terms of a systematic covariance term, C HOD , and an additional diagonal contribution quantifying the impact of our choice of nuisance parameter priors on the ability of the effective field theory (EFT) model to correctly recover the cosmological parameters of the simulations. These two covariance contributions are designed to be added to the usual covariance term, C stat , describing the statistical uncertainty in the power spectrum measurement, in order to fairly represent these sources of systematic uncertainty. This novel approach should be more general and robust to the choice of model or additional external datasets used in cosmological fits than the alternative approach of adding systematic uncertainties to the recovered marginalised parameter posteriors. We compare the approaches within the context of a fixed ΛCDM model and demonstrate that our method gives conservative estimates of the systematic uncertainty that nevertheless have little impact on the final posteriors obtained from DESI data.
Bayesian model averaging is a practical method for dealing with uncertainty due to model specification. Use of this technique requires the estimation of model probability weights. Here, we revisit the derivation of estimators for these model weights. Use of the Kullback-Leibler divergence as a starting point leads naturally to a number of alternative information criteria suitable for Bayesian model weight estimation. We explore three such criteria, known to the statistics literature before, in detail: a Bayesian analog of the Akaike information criterion which we call the BAIC, the Bayesian predictive information criterion, and the posterior predictive information criterion (PPIC). We compare the use of these information criteria in numerical analysis problems common in lattice field theory calculations. We find that the PPIC has the most appealing theoretical properties and can give the best performance in terms of model-averaging uncertainty, particularly in the presence of noisy data, while the BAIC is a simple and reliable alternative.
In this work, we develop a density functional treatment of noninteracting Abelian anyons, which is capable, in principle, of dealing with a system of a large number of anyons in an external potential. Comparison with exact results for a few particles shows that the model captures the behavior qualitatively and semiquantitatively, especially in the vicinity of the fermionic statistics. We then study anyons with statistics parameter 1+1/n, which are thought to condense into a superconducting state. An indication of the superconducting behavior is the mean-field result that, for uniform density systems, the ground-state energy increases under the application of an external magnetic field independent of its direction. Our density functional theory based analysis does not find that to be the case for finite systems of anyons, which can accommodate a weak external magnetic field through density transfer between the bulk and the boundary rather than through transitions across effective Landau levels, but the “Meissner repulsion” of the external magnetic field is recovered in the thermodynamic limit as the effect of the boundary becomes negligible. We also consider the quantum Hall effect of anyons, and show that its topological properties, such as the charge and statistics of the excitations and the quantized Hall conductance, arise in a self-consistent fashion.