No Arabic abstract
Lattice QCD simulations of multi-baryon correlation functions can predict the structure and reactions of nuclei without encountering the baryon chemical potential sign problem. However, they suffer from a signal-to-noise problem where Monte Carlo estimates of observables have quantum fluctuations that are exponentially larger than their average values. Recent lattice QCD results demonstrate that the complex phase of baryon correlations functions relates the baryon signal-to-noise problem to a sign problem and exhibits unexpected statistical behavior resembling a heavy-tailed random walk on the unit circle. Estimators based on differences of correlation function phases evaluated at different Euclidean times are discussed that avoid the usual signal-to-noise problem, instead facing a signal-to-noise problem as the time interval associated with the phase difference is increased, and allow hadronic observables to be determined from arbitrarily large-time correlation functions.
The theory of quantum chromodynamics (QCD) encodes the strong interactions that bind quarks and gluons into nucleons and that bind nucleons into nuclei. Predictive control of QCD would allow nuclear structure and reactions as well as properties of supernovae and neutron stars to be theoretically studied from first principles. Lattice QCD can represent generic QCD predictions in terms of well-defined path integrals, but the sign and signal-to-noise problems have obstructed lattice QCD calculations of large nuclei and nuclear matter in practice. This thesis presents a statistical study of lattice QCD correlation functions, with a particular focus on characterizing the structure of the noise associated with quantum fluctuations. The signal-to-noise problem in baryon correlation functions is demonstrated to arise from a sign problem associated with Monte Carlo sampling of complex correlation functions. The phases of complex correlation functions are analyzed in the framework of circular statistics, and the time evolution of the phase is shown to resemble a heavy-tailed random walk on the unit circle. Building on these observations, a new technique called phase reweighting is investigated that involves calculations of phase differences over fixed-length time intervals. Phase reweighting allows results for hadronic observables to be extracted from large-time correlation functions with constant signal-to-noise ratios. The signal-to-noise problem re-emerges as the length of the phase-difference interval is increased. Applications of phase reweighting to meson, baryon, and two-baryon systems are discussed.
Path integrals describing quantum many-body systems can be calculated with Monte Carlo sampling techniques, but average quantities are often subject to signal-to-noise ratios that degrade exponentially with time. A phase-reweighting technique inspired by recent observations of random walk statistics in correlation functions is proposed that allows energy levels to be extracted from late-time correlation functions with time-independent signal-to-noise ratios. Phase reweighting effectively includes dynamical refinement of source magnitudes but introduces a bias associated with the phase. This bias can be removed by performing an extrapolation, but at the expense of re-introducing a signal-to-noise problem. Lattice Quantum Chromodynamics calculations of the $rho$ and nucleon masses and of the $XiXi$ binding energy show consistency between standard results obtained using earlier-time correlation functions and phase-reweighted results using late-time correlation functions inaccessible to standard statistical analysis methods.
We calculate spectral functions of the relativistic $O(4)$ model from real-time lattice simulations in classical-statistical field theory. While in the low and high temperature phase of the model, the spectral functions of longitudinal $(sigma)$ and transverse $(pi)$ modes are well described by relativistic quasi-particle peaks, we find a highly non-trivial behavior of the spectral functions in the cross over region, where additional structures appear. Similarly, we observe a significant broadening of the quasi-particle peaks, when the amount explicit $O(4)$ symmetry breaking is reduced. We further demonstrate that in the vicinity of the $O(4)$ critical point, the spectral functions develop an infrared power law associated with the critical dynamics, and comment on the extraction of the dynamical critical exponent $z$ from our simulations.
Using the generating functional analysis an exact recursion relation is derived for the time evolution of the effective local field of the fully connected Little-Hopfield model. It is shown that, by leaving out the feedback correlations arising from earlier times in this effective dynamics, one precisely finds the recursion relations usually employed in the signal-to-noise approach. The consequences of this approximation as well as the physics behind it are discussed. In particular, it is pointed out why it is hard to notice the effects, especially for model parameters corresponding to retrieval. Numerical simulations confirm these findings. The signal-to-noise analysis is then extended to include all correlations, making it a full theory for dynamics at the level of the generating functional analysis. The results are applied to the frequently employed extremely diluted (a)symmetric architectures and to sequence processing networks.
High-temperature expansions are presently the only viable approach to the numerical calculation of the higher susceptibilities for the spin and the scalar-field models on high-dimensional lattices. The critical amplitudes of these quantities enter into a sequence of universal amplitude-ratios which determine the critical equation of state. We have obtained a substantial extension through order 24, of the high-temperature expansions of the free energy (in presence of a magnetic field) for the Ising models with spin s >= 1/2 and for the lattice scalar field theory with quartic self-interaction, on the simple-cubic and the body-centered-cubic lattices in four, five and six spatial dimensions. A numerical analysis of the higher susceptibilities obtained from these expansions, yields results consistent with the widely accepted ideas, based on the renormalization group and the constructive approach to Euclidean quantum field theory, concerning the no-interaction (triviality) property of the continuum (scaling) limit of spin-s Ising and lattice scalar-field models at and above the upper critical dimensionality.