BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Bernd A. Berg contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
It is shown that whenever the multiplicative normalization of a fitting function is not known, least square fitting by $χ^2$ minimization can be performed with one parameter less than usual by converting the normalization parameter into a function of the remaining parameters and the data.
Recently the Yang-Mills gradient flow of pure SU(3) lattice gauge theory has been calculated in the range from $β=6/g_0^2=6.3$ to~7.5 (Asakawa et al.), where $g_0^2$ is the bare coupling constant of the SU(3) Wilson action. Estimates of the deconfining phase transition are available from $β=5.7$ to~6.8 (Francis et al.). Here it is shown that the entire range from 5.7 to 7.5 is well described by a power series of the lattice spacing $a$ times the lambda lattice mass scale $Λ_L$, using asymptotic scaling in the 2-loop and 3-loop approximations for $aΛ_L$. In both cases identical ratios for gradient flows versus deconfinement observables are obtained. Differences in the normalization constants with respect to $Λ_L$ give a handle on their systematic errors.
With the emergence of the Yang-Mills gradient flow technique there is renewed interest in the issue of scale setting in lattice gauge theory. Here I compare for the SU(3) Wilson gauge action non-perturbative scale functions of Edwards, Heller and Klassen (EHK), Necco and Sommer (NS), both relying on Sommer's method using the quark potential, and the scale function derived by Bazavov, Berg and Velytsky (BBV) from a deconfining phase transition investigation by the Bielefeld group. It turns out that the scale functions are based on mutually inconsistent data, though the BBV scale function is consistent with the EHK data when their low $β$ ($β=5.6$) data point is removed. Besides, only the BBV scale function is consistent with three data points calculated from the gradient flow by Lüscher. In the range for which data exist the discrepancies between the scale functions are only up to $\pm 2$\% of their values, but clearly visible within the statistical accuracy.
Statistical methods for the extraction of a small shift in broad data distributions are examined by means of Monte Carlo simulations. This work was originally motivated by the CERN neutrino beam to Gran Sasso (CNGS) experiment for which the OPERA detector collaboration reported a time shift in a broad distribution with an accuracy of $\pm 7.8\,$ns, while the fluctuation of the average time turns with $\pm 23.8\,$ns out to be much larger. Although the physical result of a big shift has been withdrawn, statistical methods that make an identification in a broad distribution with such a small error possible remain of interest.
Using the geometry of a double-layered torus we investigate the deconfining phase transition of pure SU(3) lattice gauge theory by Markov chain Monte Carlo simulations. In one layer, called "outside", the temperature is set below the deconfining temperature and in the other, called "inside", it is iterated to a pseudo-transition temperature. Lattice sizes are chosen in a range suggested by the physical volumes achieved in relativistic heavy ion collisions and both temperatures are kept close enough to stay in the SU(3) scaling region. Properties of the transition are studied as function of the volume for three outside temperatures. When compared with infinite volume extrapolations, small volume corrections of the deconfining temperature and width become competitive with those found by including quarks. Effective finite size scaling exponents of the specific and Polyakov loop susceptibilities are also calculated.
Fisher zeros are the zeros of the partition function in the complex beta=2N_c/g^2 plane. When they pinch the real axis, finite size scaling allows one to distinguish between first and second order transition and to estimate exponents. On the other hand, a gap signals confinement and the method can be used to explore the boundary of the conformal window. We present recent numerical results for 2D O(N) sigma models, 4D U(1) and SU(2) pure gauge and SU(3) gauge theory with N_f=4 and 12 flavors. We discuss attempts to understand some of these results using analytical methods. We discuss the 2-lattice matching and qualitative aspects of the renormalization group (RG) flows in the Migdal-Kadanoff approximation, in particular how RG flows starting at large beta seem to move around regions where bulk transitions occur. We consider the effects of the boundary conditions on the nonperturbative part of the average energy and on the Fisher zeros for the 1D O(2) model.
It is noted that, in contrast to widespread believes, free fields do not only allow for global, but also for local gauge invariance.
For a free scalar boson field and for U(1) gauge theory finite volume (infrared) and other corrections to the energy-momentum dispersion in the lattice regularization are investigated calculating energy eigenstates from the fall off behavior of two-point correlation functions. For small lattices the squared dispersion energy defined by $E_{\rm dis}^2=E_{\vec{k}}^2-E_0^2-4\sum_{i=1}^{d-1}\sin(k_i/2)^2$ is in both cases negative ($d$ is the Euclidean space-time dimension and $E_{\vec{k}}$ the energy of momentum $\vec{k}$ eigenstates). Observation of $E_{\rm dis}^2=0$ has been an accepted method to demonstrate the existence of a massless photon ($E_0=0$) in 4D lattice gauge theory, which we supplement here by a study of its finite size corrections. A surprise from the lattice regularization of the free field is that infrared corrections do {\it not} eliminate a difference between the groundstate energy $E_0$ and the mass parameter $M$ of the free scalar lattice action. Instead, the relation $E_0=\cosh^{-1} (1+M^2/2)$ is derived independently of the spatial lattice size.
Previously a model of only vector fields with a local U(2) symmetry was introduced for which one finds a massless U(1) photon and a massive SU(2) vector boson in the lattice regularization. Here it is shown that quantization of its classical continuum action leads to perturbative renormalization difficulties. But, non-perturbative Monte Carlo calculations favor the existence of a quantum continuum limit.
Based on cumulative distribution functions, Fourier series expansion and Kolmogorov tests, we present a simple method to display probability densities for data drawn from a continuous distribution. It is often more efficient than using histograms.
Deconfined regions in relativistic heavy ion collisions are limited to small volumes surrounded by a confined exterior. Using the geometry of a double layered torus, we keep an outside temperature slightly lower than the inside temperature, so that both regions are in the SU(3) scaling region. Deconfined volume sizes are chosen to be in a range typical for such volumes created at the BNL RHIC. Even with small temperature differences a dependence of the (pseudo) deconfining temperature on a colder surrounding temperature is clearly visible. For temporal lattice sizes Ntau=4, 6 and 8 we find consistency with SU(3) scaling behavior for the measured transition temperature signals.
Two simple exercises are solved, which educators can use to awake interest of their students in subtleties of the CERN Neutrino beam to Grand Sasso (CNGS) experiment. The first one is about the statistical error of the average departure time of neutrinos from CERN. The second one about a hypothetical bias in the departure times.
Markov chain Monte Carlo simulations of pure SU(2)xU(1) lattice gauge theory show a (zero temperature) deconfining phase transition in the SU(2) gluon sector when a term is added to the SU(2) and U(1) Wilson actions, which requires joint U(2) gauge transformations of the SU(2) and U(1) vector fields. Investigations of this deconfined phase are of interest as it could provide an alternative to the Higgs mechanism.
Deconfined regions in relativistic heavy ion collisions are limited to small volumes surrounded by a confined exterior. Here the geometry of a double layered torus is discussed, which allows for different temperatures in its two layers. This geometry enables one to approach the QCD continuum limit for small deconfined volumes with confined exteriors in a more realistic fashion than by using periodic boundary conditions. Preliminary data from a study for pure SU(3) lattice gauge theory support a substantial increase in a pseudo transition temperature.
Gauge invariance is extended to allow to allow for a U(1)-SU(2) mixing term, which can cause a SU(2) deconfining transition.
We document Fortran MPI checkerboard code for Markov Chain Monte Carlo simulations of pure SU(3) lattice gauge theory with the Wilson action on a D-dimensional double-layered torus. This includes the usual torus with periodic boundary conditions as an optional case. We use Cabibbo-Marinari heatbath checkerboard updating. Parallelization on sublattices is implemented in all D directions and can be restricted to less than D directions. The parallelization techniques of this paper can be used for any model with interactions of link variables defined on plaquettes.
We study the performance of MPI checkerboard code for SU(3) lattice gauge theory as function of the number of MPI processes, which run in parallel on an identical number of CPU cores. Computing platforms explored are a small PC cluster at FSU and the Cray at NERSC.
We study Kolmogorov-Johnson-Mehl-Avrami (KJMA) theory of phase conversion in finite volumes. For the conversion time we find the relationship $τ_{\rm con} = τ_{\rm nu} [1+f_d(q)]$. Here $d$ is the space dimension, $τ_{\rm nu}$ the nucleation time in the volume $V$, and $f_d(q)$ a scaling function. Its dimensionless argument is $q=τ_{\rm ex}/ τ_{\rm nu}$, where $τ_{\rm ex}$ is an expansion time, defined to be proportional to the diameter of the volume divided by expansion speed. We calculate $f_d(q)$ in one, two and three dimensions. The often considered limits of phase conversion via either nucleation or spinodal decomposition are found to be volume-size dependent concepts, governed by simple power laws for $f_d(q)$.
We introduce two simple models with nearest neighbor interactions on 3D hexagonal lattices. Each model allows one to calculate the residual entropy of ice I (ordinary ice) by means of multicanonical simulations. This gives the correction to the residual entropy derived by Linus Pauling in 1935. Our estimate is found to be within less than 0.1% of an analytical approximation by Nagle, which is an improvement of Pauling's result. We pose it as a challenge to experimentalists to improve on the accuracy of a 1936 measurement by Giauque and Stout by about one order of magnitude, which would allow one to identify corrections to Pauling's value unambiguously. It is straightforward to transfer our methods to other crystal systems.
This article is a tutorial on Markov chain Monte Carlo simulations and their statistical analysis. The theoretical concepts are illustrated through many numerical assignments from the author's book on the subject. Computer code (in Fortran) is available for all subjects covered and can be downloaded from the web.
The most efficient MC weights for the calculation of physical, canonical expectation values are not necessarily those of the canonical ensemble. The use of suitably generalized ensembles can lead to a much faster convergence of the simulation. Although not realized by nature, these ensembles can be implemented on computers. In recent years generalized ensembles have in particular been studied for the simulation of complex systems. For these systems it is typical that conflicting constraints lead to free energy barriers, which fragment the configuration space. Examples of major interest are spin glasses and proteins. In my overview I first comment on the strengths and weaknesses of a few major approaches, multicanonical simulations, transition variable methods, and parallel tempering. Subsequently, two applications are presented: a new analysis of the Parisi overlap distribution for the 3d Edwards-Anderson Ising spin glass and the helix-coil transition of amino-acid homo-oligomers.
Monte Carlo simulation with {\it a-priori} unknown weights have attracted recent attention and progress has been made in understanding (i) the technical feasibility of such simulations and (ii) classes of systems for which such simulations lead to major improvements over conventional Monte Carlo simulations. After briefly sketching the history of multicanonical calculations and their range of application, a general introduction in the context of the statistical physics of the d-dimensional generalized Potts models is given. Multicanonical simulations yield canonical expectation values for a range of temperatures or any other parameter(s) for which appropriate weights can be constructed. We shall address in some details the question how the multicanonical weights are actually obtained. Subsequently miscellaneous topics related to the considered algorithms are reviewed. Then multicanonical studies of first order phase transitions are discussed and finally applications to complex systems such as spin glasses and proteins.
We investigate the spectrum of the staggered Dirac operator in 4d quenched U(1) lattice gauge theory and its relationship to random matrix theory. In the confined as well as in the Coulomb phase the nearest-neighbor spacing distribution of the unfolded eigenvalues is well described by the chiral unitary ensemble. The same is true for the distribution of the smallest eigenvalue and the microscopic spectral density in the confined phase. The physical origin of the chiral condensate in this phase deserves further study.
The problem of calculating multicanonical parameters recursively is discussed. I describe in detail a computational implementation which has worked reasonably well in practice.