Walking near a conformal fixed point

We investigate the 2-d O(3) model with a q-term as a toy model for slowly walking 4-d non-Abelian gauge theories. Using the very efficient meron-cluster algorithm, an accurate investigation of the scale dependence of the renormalized coupling is carried out for different values of the vacuum angle q. Approaching q = p, the infrared dynamics of the 2-d O(3) model is determined by a non-trivial conformal fixed point. We provide evidence for a slowly walking behavior near the fixed point and we perform a finite-size scaling analysis of the mass gap.

Hot QCD vs cosmology

It is a long-standing dream to “simulate” cosmology in laboratory through heavy ion collision experiments. Although the QCD epoch itself may not leave major cosmological signatures, theoretical methods developed and tested in the context of heavy ion collision experiments could indeed find applications at other energy scales. Here recent progress in this spirit is reviewed.

Factorization of leading chiral logarithms in the pion form factors

In the recently proposed framework of hard pion chiral perturbation theory, the leading chiral logarithms are predicted to factorize with respect to the energy dependence in the chiral limit. We have scrutinized this assumption in the case of vector and scalar pion form factors FV;S(s) by means of standard chiral perturbation theory and dispersion relations. We show that this factorization property is valid for the elastic contribution to the dispersion integrals for FV;S(s) but it is violated starting at three loops when the inelastic four-pion contributions arise.

Loop formulation of the supersymmetric nonlinear O(N) sigma model

We derive the fermion loop formulation for the supersymmetric nonlinear O(N) sigma model by performing a hopping expansion using Wilson fermions. In this formulation the fermionic contribution to the partition function becomes a sum over all possible closed non-oriented fermion loop configurations. The interaction between the bosonic and fermionic degrees of freedom is encoded in the constraints arising from the supersymmetry and induces flavour changing fermion loops. For N ≥ 3 this leads to fermion loops which are no longer self-avoiding and hence to a potential sign problem. Since we use Wilson fermions the bare mass needs to be tuned to the chiral point. For N = 2 we determine the critical point and present boson and fermion masses in the critical regime.

A first look at maximally twisted mass lattice QCD calculations at the physical point

In this contribution, a first look at simulations using maximally twisted mass Wilson fermions at the physical point is presented. A lattice action including clover and twisted mass terms is presented and the Monte Carlo histories of one run with two mass-degenerate flavours at a single lattice spacing are shown. Measurements from the light and heavy-light pseudoscalar sectors are compared to previous Nf = 2 results and their phenomenological values. Finally, the strategy for extending simulations to Nf = 2+1+1 is outlined.

Experiences with OpenMP in tmLQCD

An overview is given of the lessons learned from the introduction of multi-threading using OpenMP in tmLQCD. In particular, programming style, performance measurements, cache misses, scaling, thread distribution for hybrid codes, race conditions, the overlapping of communication and computation and the measurement and reduction of certain overheads are discussed. Performance measurements and sampling profiles are given for different implementations of the hopping matrix computational kernel.

Recent developments in the tmLQCD software suite

We present an overview of recent developments in the tmLQCD software suite. We summarise the features of the code, including actions and operators implemented. In particular, we discuss the optimisation efforts for modern architectures using the Blue Gene/Q system as an example.

First order thermal phase transition with 126 GeV Higgs mass

We study the strength of the electroweak phase transition in models with two light Higgs doublets and a light SU(3)c triplet by means of lattice simulations in a dimensionally reduced effective theory. In the parameter region considered the transition on the lattice is significantly stronger than indicated by a 2-loop perturbative analysis. Within some ultraviolet uncertainties, the finding applies to MSSM with a Higgs mass mh ≈ 126 GeV and shows that the parameter region useful for electroweak baryogenesis is enlarged. In particular (even though only dedicated analyses can quantify the issue), the tension between LHC constraints after the 7 TeV and 8 TeV runs and frameworks where the electroweak phase transition is driven by light stops, seems to be relaxed.

Towards the continuum limit in transport coefficient computations

The analytic continuation needed for the extraction of transport coefficients necessitates in principle a continuous function of the Euclidean time variable. We report on progress towards achieving the continuum limit for 2-point correlator measurements in thermal SU(3) gauge theory, with specific attention paid to scale setting. In particular, we improve upon the determination of the critical lattice coupling and the critical temperature of pure SU(3) gauge theory, estimating r0Tc ≃ 0.7470(7) after a continuum extrapolation. As an application the determination of the heavy quark momentum diffusion coefficient from a correlator of colour-electric fields attached to a Polyakov loop is discussed.

Temporal mesonic correlators at NLO for any quark mass

We present NLO results for thermal imaginary-time correlators in the vector and scalar channels as a function of the quark mass. The range of quark masses for which a non-relativistic approximation works in the temperature range considered is estimated, and charm quarks turn out to be a borderline case. Comparing with simulation data from fine lattices, we find good agreement in the vector channel but a substantial discrepancy in the scalar one. An explanation for the discrepancy is suggested in terms of physics of the quark-antiquark threshold region. Perturbative predictions for the bottom scalar spectral function around the threshold are also briefly reviewed.

Towards understanding thermal jet quenching via lattice simulations

After reviewing how simulations employing classical lattice gauge theory permit to test a conjectured Euclideanization property of a light-cone Wilson loop in a thermal non-Abelian plasma, we show how Euclidean data can in turn be used to estimate the transverse collision kernel, C(k⊥), characterizing the broadening of a high-energy jet. First results, based on data produced recently by Panero et al, suggest that C(k⊥) is enhanced over the known NLO result in a soft regime k⊥ < a few T. The shape of k3⊥ C(k⊥) is consistent with a Gaussian at small k⊥.

A new Bayesian approach to the reconstruction of spectral functions

We present a novel approach for the reconstruction of spectra from Euclidean correlator data that makes close contact to modern Bayesian concepts. It is based upon an axiomatically justified dimensionless prior distribution, which in the case of constant prior function m(ω) only imprints smoothness on the reconstructed spectrum. In addition we are able to analytically integrate out the only relevant overall hyper-parameter α in the prior, removing the necessity for Gaussian approximations found e.g. in the Maximum Entropy Method. Using a quasi-Newton minimizer and high-precision arithmetic, we are then able to find the unique global extremum of P[ρ|D] in the full Nω » Nτ dimensional search space. The method actually yields gradually improving reconstruction results if the quality of the supplied input data increases, without introducing artificial peak structures, often encountered in the MEM. To support these statements we present mock data analyses for the case of zero width delta peaks and more realistic scenarios, based on the perturbative Euclidean Wilson Loop as well as the Wilson Line correlator in Coulomb gauge.

Benchmarking the Bayesian reconstruction of the non-perturbative heavy ǪǬ potential

The extraction of the finite temperature heavy quark potential from lattice QCD relies on a spectral analysis of the real-time Wilson loop. Through its position and shape, the lowest lying spectral peak encodes the real and imaginary part of this complex potential. We benchmark this extraction strategy using leading order hard-thermal loop (HTL) calculations. I.e. we analytically calculate the Wilson loop and determine the corresponding spectrum. By fitting its lowest lying peak we obtain the real- and imaginary part and confirm that the knowledge of the lowest peak alone is sufficient for obtaining the potential. We deploy a novel Bayesian approach to the reconstruction of spectral functions to HTL correlators in Euclidean time and observe how well the known spectral function and values for the real and imaginary part are reproduced. Finally we apply the method to quenched lattice QCD data and perform an improved estimate of both real and imaginary part of the non-perturbative heavy ǪǬ potential.

Lattice NRQCD study of in-medium bottomonium states using Nf = 2+1, 48³ × 12 HotQCD configurations

The behavior of bottomonium state correlators at non-zero temperature, 140.4(β = 6.664) ≤ T ≤ 221(β = 7.280) (MeV), where the transition temperature is 154(9) (MeV), is studied, using lattice NRQCD on 48³ ×12 HotQCD HiSQ action configurations with light dynamical Nf = 2+1 (mu,s/ms = 0.05) staggered quarks. In order to understand finite temperature effects on quarkonium states, zero temperature behavior of bottomonium correlators is compared based on 32⁴ (β = 6.664,6.800 and 6.950) and 48³ ×64 (β = 7.280) lattices. We find that temperature effects on S-wave bottomoniumstates are small but P-wave bottomoniumstates show a noticeable temperature dependence above the transition temperature.

Crystalline Confinement

We show that exotic phases arise in generalized lattice gauge theories known as quantum link models in which classical gauge fields are replaced by quantum operators. While these quantum models with discrete variables have a finite-dimensional Hilbert space per link, the continuous gauge symmetry is still exact. An efficient cluster algorithm is used to study these exotic phases. The (2+1)-d system is confining at zero temperature with a spontaneously broken translation symmetry. A crystalline phase exhibits confinement via multi stranded strings between chargeanti-charge pairs. A phase transition between two distinct confined phases is weakly first order and has an emergent spontaneously broken approximate SO(2) global symmetry. The low-energy physics is described by a (2 + 1)-d RP(1) effective field theory, perturbed by a dangerously irrelevant SO(2) breaking operator, which prevents the interpretation of the emergent pseudo-Goldstone boson as a dual photon. This model is an ideal candidate to be implemented in quantum simulators to study phenomena that are not accessible using Monte Carlo simulations such as the real-time evolution of the confining string and the real-time dynamics of the pseudo-Goldstone boson.