Source author record

M. Laine

M. Laine appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

81works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

81 published item(s)

preprint2026arXiv

Classical equipartition dynamics between axions and non-Abelian gauge fields

Motivated by axion-like inflation and its warm embedding within the Standard Model, we study the early stages of the energy transfer between an axion condensate and an SU(2) gauge ensemble, by employing non-linear classical real-time lattice simulations. The discretized equations of motion are worked out, elaborating on Gauss constraints. A numerical solution is implemented on the CosmoLattice platform. Adopting a quadratic potential, and omitting universe expansion for the moment, we establish initial exponential growth of the low-momentum gauge modes; damping of axion oscillations after some delay; and subsequent energy equipartition between axion and gauge ensembles. A clear difference between the SU(2) and U(1) dynamics is observed, likely associated with non-Abelian self-interactions. We elaborate on what this implies for the possible thermalization of the SU(2) ensemble.

preprint2026arXiv

Evolution of coupled scalar perturbations through smooth reheating. II. Thermal fluctuation regime

Curvature perturbations with short wavelengths exit the Hubble horizon when the universe may contain a thermal plasma in addition to an inflaton field that drives its expansion. We solve the corresponding fluctuation-dissipation dynamics at linear order, building upon a previously established set of gauge-invariant evolution equations. The properties of the noise autocorrelator are constrained via a matching of equilibrium correlators to quantum-statistical physics deep inside the Hubble horizon. The curvature power spectrum is determined numerically, without slow-roll approximations or assumptions about the equilibration of the inflaton field. As applications, we scrutinize two issues from recent literature: the model dependence of the thermally modified power spectrum as a function of freeze-out parameters, and the viability of embedding warm inflation within the Standard Model (we offer support for the latter proposal). The role of pre-horizon-exit acoustic oscillations is illustrated.

preprint2022arXiv

Gravitational wave background from non-Abelian reheating after axion-like inflation

A pseudoscalar inflaton $φ$, coupled to the topological charge density $F\tilde{F}$ of a non-Abelian sector, can decay to gauge bosons ($φ\to g g $), which may thermalize rapidly. The friction felt by $φ$ is then increased by non-Abelian "strong sphalerons", leading to a self-amplifying process that can efficiently heat up the medium. We determine a lower bound for the gravitational wave production rate from such a process, originating via hydrodynamic fluctuations and particle collisions, in terms of a minimal number of parameters. Only a moderate fraction of energy density is converted to gravitational waves, suggesting that non-Abelian models may avoid the overproduction observed in some Abelian cases.

preprint2022arXiv

Lattice study of a magnetic contribution to heavy quark momentum diffusion

Heavy quarks placed within a hot QCD medium undergo Brownian motion, characterized by specific transport coefficients. Their determination can be simplified by expanding them in $T/M$, where $T$ is the temperature and $M$ is a heavy quark mass. The leading term in the expansion originates from the colour-electric part of a Lorentz force, whereas the next-to-leading order involves the colour-magnetic part. We measure a colour-magnetic 2-point correlator in quenched QCD at $T \sim 1.2 ... 2.0 T_{\rm c}^{ }$. Employing multilevel techniques and non-perturbative renormalization, a good signal is obtained, and its continuum extrapolation can be estimated. Modelling the shape of the corresponding spectral function, we subsequently extract the momentum diffusion coefficient, $κ$. For charm (bottom) quarks, the magnetic contribution adds $\sim 30\%$ ($10\%$) to the electric one. The same increases apply also to the drag coefficient, $η$. As an aside, the colour-magnetic spectral function is computed at NLO.

preprint2022arXiv

Mass-suppressed effects in heavy quark diffusion

Many lattice studies of heavy quark diffusion originate from a colour-electric correlator, obtained as a leading term after an expansion in the inverse of the heavy-quark mass. In view of the fact that the charm quark is not particularly heavy, we consider subleading terms in the expansion. Working out correlators up to $O(1/M^2)$, we argue that the leading corrections are suppressed by $O(T/M)$, and one of them can be extracted from a colour-magnetic correlator. The corresponding transport coefficient is non-perturbative already at leading order in the weak-coupling expansion, and therefore requires a non-perturbative determination.

preprint2022arXiv

Sterile neutrino rates for general $M$, $T$, $μ$, $k$: review of a theoretical framework

The temperature of sphaleron freeze-out, $T\approx 130$ GeV, sets a scale for mechanisms for generating lepton asymmetries and converting them to baryon asymmetry (leptogenesis). For Majorana masses $M \gg πT$, leptogenesis takes place in the non-relativistic regime, while for $M \ll πT$, the dynamics is ultrarelativistic. The intermediate case $M \sim πT$ is the most cumbersome, as no effective theory methods are available. We review the definitions of and provide integral representations for all rate coefficients and mass corrections of $\mathcal{O}(h_ν^2)$ in this regime, such that the expressions extrapolate to known limits and are gauge independent to the order computed, both in the symmetric and the Higgs phase, for any helicity, momentum, and set of chemical potentials. Road signs towards a numerical evaluation are offered.

preprint2021arXiv

Smooth interpolation between thermal Born and LPM rates

In a weakly coupled ultrarelativistic plasma, $1 + n \leftrightarrow 2 + n$ scatterings, with $n \ge 0$, need sometimes to be summed to all orders, in order to determine a leading-order interaction rate. To implement this "LPM resummation", kinematic approximations are invoked. However, in cosmological settings, where the temperature changes by many orders of magnitude and both small and large momenta may play a role, such approximations are not always justified. We suggest a procedure to smoothly interpolate between LPM-resummed $1 + n \leftrightarrow 2 + n$ and Born-level $1 \leftrightarrow 2$ results, rendering the outcome applicable to a broader range of masses and momenta. The procedure is illustrated for right-handed neutrino production from a Standard Model plasma, and dilepton production from a QCD plasma.

preprint2020arXiv

A QCD Debye mass in a broad temperature range

The Debye mass sets a scale for the screening of static charges and the scattering of fast charges within a gauge plasma. Inspired by its potential cosmological applications, we determine a QCD Debye mass at 2-loop order in a broad temperature range (1 GeV ... 10 TeV), demonstrating how quark mass thresholds get smoothly crossed. Along the way, integration-by-parts identities pertinent to massive loops at finite temperature are illuminated.

preprint2020arXiv

Gravitational wave background from Standard Model physics: Complete leading order

We compute the production rate of the energy density carried by gravitational waves emitted by a Standard Model plasma in thermal equilibrium, consistently to leading order in coupling constants for momenta $k\sim πT$. Summing up the contributions from the full history of the universe, the highest temperature of the radiation epoch can be constrained by the so-called $N_{\rm eff}$ parameter. The current theoretical uncertainty $ΔN_{\rm eff} \le 10^{-3}$ corresponds to $T_{\rm max} \le 2\times 10^{17}$ GeV. In the course of the computation, we show how a subpart of the production rate can be determined with the help of standard packages, even if subsequently an IR subtraction and thermal resummation need to be implemented.

preprint2020arXiv

Sterile neutrino dark matter via coinciding resonances

It has been proposed that two resonances could coincide in the early universe at temperatures $T \sim 0.2 ... 0.5$ GeV: one between two nearly degenerate GeV-scale sterile neutrinos, producing a large lepton asymmetry through freeze-out and decays; another between medium-modified active neutrinos and keV-scale sterile neutrinos, converting the lepton asymmetry into dark matter. Making use of a framework which tracks three sterile neutrinos of both helicities as well as three separate lepton asymmetries, and scanning the parameter space of the GeV-scale species, we establish the degree of fine-tuning that is needed for realizing this scenario.

preprint2019arXiv

A thermal neutrino interaction rate at NLO

The interaction rate of an ultrarelativistic active neutrino at a temperature below the electroweak crossover plays a role in leptogenesis scenarios based on oscillations between active neutrinos and GeV-scale sterile neutrinos. By making use of a Euclideanization property of a thermal light-cone correlator, we determine the $O(g)$ correction to such an interaction rate in the high-temperature limit $πT \gg m_W$, finding a $\sim 15 ... 40\%$ reduction. For a benchmark point, this NLO correction decreases the lepton asymmetries produced by $\sim 1\%$.

preprint2016arXiv

Heavy ions at the Future Circular Collider

The Future Circular Collider (FCC) Study is aimed at assessing the physics potential and the technical feasibility of a new collider with centre-of-mass energies, in the hadron-hadron collision mode, seven times larger than the nominal LHC energies. Operating such machine with heavy ions is an option that is being considered in the accelerator design studies. It would provide, for example, Pb-Pb and p-Pb collisions at sqrt{s_NN} = 39 and 63 TeV, respectively, per nucleon-nucleon collision, with integrated luminosities above 30 nb^-1 per month for Pb-Pb. This is a report by the working group on heavy-ion physics of the FCC Study. First ideas on the physics opportunities with heavy ions at the FCC are presented, covering the physics of the Quark-Gluon Plasma, of gluon saturation, of photon-induced collisions, as well as connections with other fields of high-energy physics.

preprint2016arXiv

Lattice constraints on the thermal photon rate

We estimate the photon production rate from an SU(3) plasma at temperatures of about 1.1Tc and 1.3Tc. Lattice results for the vector current correlator at spatial momenta k ~ (2-6)T are extrapolated to the continuum limit and analyzed with the help of a polynomial interpolation for the corresponding spectral function, which vanishes at zero frequency and matches to high-precision perturbative results at large invariant masses. For small invariant masses the interpolation is compared with the NLO weak-coupling result, hydrodynamics, and a holographic model. At vanishing invariant mass we extract the photon rate which for k \gsim 3T is found to be close to the NLO weak-coupling prediction. For k \lsim 2T uncertainties remain large but the photon rate is likely to fall below the NLO prediction, in accordance with the onset of a strongly interacting behaviour characteristic of the hydrodynamic regime.

preprint2016arXiv

Neutrino dynamics below the electroweak crossover

We estimate the thermal masses and damping rates of active (m < eV) and sterile (M ~ GeV) neutrinos with thermal momenta k ~ 3T at temperatures below the electroweak crossover (5 GeV < T < 160 GeV). These quantities fix the equilibration or "washout" rates of Standard Model lepton number densities. Sterile neutrinos interact via direct scatterings mediated by Yukawa couplings, and via their overlap with active neutrinos. Including all leading-order reactions we find that the washout rate generally exceeds the Hubble rate for 5 GeV < T < 30 GeV. Therefore it is challenging to generate a large lepton asymmetry facilitating dark matter computations operating at T < 5 GeV, whereas the generation of a baryon asymmetry at T > 130 GeV remains an option. Our differential rates are tabulated in a form suitable for studies of specific scenarios with given neutrino Yukawa matrices.

preprint2016arXiv

Perturbative renormalization of the electric field correlator

The momentum diffusion coefficient of a heavy quark in a hot QCD plasma can be extracted as a transport coefficient related to the correlator of two colour-electric fields dressing a Polyakov loop. We determine the perturbative renormalization factor for a particular lattice discretization of this correlator within Wilson's SU(3) gauge theory, finding a ~12% NLO correction for values of the bare coupling used in the current generation of simulations. The impact of this result on existing lattice determinations is commented upon, and a possibility for non-perturbative renormalization through the gradient flow is pointed out.

preprint2016arXiv

Rapid thermal co-annihilation through bound states in QCD

The co-annihilation rate of heavy particles close to thermal equilibrium, which plays a role in many classic dark matter scenarios, can be "simulated" in QCD by considering the pair annihilation rate of a heavy quark and antiquark at a temperature of a few hundred MeV. We show that the so-called Sommerfeld factors, parameterizing the rate, can be defined and measured non-perturbatively within the NRQCD framework. Lattice measurements indicate a modest suppression in the octet channel, in reasonable agreement with perturbation theory, and a large enhancement in the singlet channel, much above the perturbative prediction. The additional enhancement is suggested to originate from bound state formation and subsequent decay. Making use of a Green's function based method to incorporate thermal corrections in perturbative co-annihilation rate computations, we show that qualitative agreement with lattice data can be found once thermally broadened bound states are accounted for. We suggest that our formalism may also be applicable to specific dark matter models which have complicated bound state structures.

preprint2015arXiv

A non-perturbative estimate of the heavy quark momentum diffusion coefficient

We estimate the momentum diffusion coefficient of a heavy quark within a pure SU(3) plasma at a temperature of about 1.5Tc. Large-scale Monte Carlo simulations on a series of lattices extending up to 192^3*48 permit us to carry out a continuum extrapolation of the so-called colour-electric imaginary-time correlator. The extrapolated correlator is analyzed with the help of theoretically motivated models for the corresponding spectral function. Evidence for a non-zero transport coefficient is found and, incorporating systematic uncertainties reflecting model assumptions, we obtain kappa = (1.8 - 3.4)T^3. This implies that the "drag coefficient", characterizing the time scale at which heavy quarks adjust to hydrodynamic flow, is (1.8 - 3.4) (Tc/T)^2 (M/1.5GeV) fm/c, where M is the heavy quark kinetic mass. The results apply to bottom and, with somewhat larger systematic uncertainties, to charm quarks.

preprint2015arXiv

Charm contribution to bulk viscosity

In the range of temperatures reached in future heavy ion collision experiments, hadronic pair annihilations and creations of charm quarks may take place within the lifetime of the plasma. As a result, charm quarks may increase the bulk viscosity affecting the early stages of hydrodynamic expansion. Assuming thermalization, we estimate the charm contribution to bulk viscosity within the same effective kinetic theory framework in which the light parton contribution has been computed previously. The time scale at which this physics becomes relevant is related to the width of the transport peak associated with the trace anomaly correlator, and is found to be < 20 fm/c for T > 600 MeV.

preprint2015arXiv

Critical point and scale setting in SU(3) plasma: An update

We explore a method developed in statistical physics which has been argued to have exponentially small finite-volume effects, in order to determine the critical temperature Tc of pure SU(3) gauge theory close to the continuum limit. The method allows us to estimate the critical coupling betac of the Wilson action for temporal extents up to Nt ~ 20 with < 0.1% uncertainties. Making use of the scale setting parameters r0 and sqrt{t0} in the same range of beta-values, these results lead to the independent continuum extrapolations Tc r0 = 0.7457(45) and Tc sqrt{t0} = 0.2489(14), with the latter originating from a more convincing fit. Inserting a conversion of r0 from literature (unfortunately with much larger errors) yields Tc / LambdaMSbar = 1.24(10).

preprint2015arXiv

Gravitational wave background from Standard Model physics: Qualitative features

Because of physical processes ranging from microscopic particle collisions to macroscopic hydrodynamic fluctuations, any plasma in thermal equilibrium emits gravitational waves. For the largest wavelengths the emission rate is proportional to the shear viscosity of the plasma. In the Standard Model at T > 160 GeV, the shear viscosity is dominated by the most weakly interacting particles, right-handed leptons, and is relatively large. We estimate the order of magnitude of the corresponding spectrum of gravitational waves. Even though at small frequencies (corresponding to the sub-Hz range relevant for planned observatories such as eLISA) this background is tiny compared with that from non-equilibrium sources, the total energy carried by the high-frequency part of the spectrum is non-negligible if the production continues for a long time. We suggest that this may constrain (weakly) the highest temperature of the radiation epoch. Observing the high-frequency part directly sets a very ambitious goal for future generations of GHz-range detectors.

preprint2015arXiv

Improved determination of sterile neutrino dark matter spectrum

The putative recent indication of an unidentified 3.55 keV X-ray line in certain astrophysical sources is taken as a motivation for an improved theoretical computation of the cosmological abundance of 7.1 keV sterile neutrinos. If the line is interpreted as resulting from the decay of Warm Dark Matter, the mass and mixing angle of the sterile neutrino are known. Our computation then permits for a determination of the lepton asymmetry that is needed for producing the correct abundance via the Shi-Fuller mechanism, as well as for an estimate of the non-equilibrium spectrum of the sterile neutrinos. The latter plays a role in structure formation simulations. Results are presented for different flavour structures of the neutrino Yukawa couplings and for different types of pre-existing lepton asymmetries, accounting properly for the charge neutrality of the plasma and incorporating approximately hadronic contributions.

preprint2015arXiv

Standard Model thermodynamics across the electroweak crossover

Even though the Standard Model with a Higgs mass mH = 125 GeV possesses no bulk phase transition, its thermodynamics still experiences a "soft point" at temperatures around T = 160 GeV, with a deviation from ideal gas thermodynamics. Such a deviation may have an effect on precision computations of weakly interacting dark matter relic abundances if their mass is in the few TeV range, or on leptogenesis scenarios operating in this temperature range. By making use of results from lattice simulations based on a dimensionally reduced effective field theory, we estimate the relevant thermodynamic functions across the crossover. The results are tabulated in a numerical form permitting for their insertion as a background equation of state into cosmological particle production/decoupling codes. We find that Higgs dynamics induces a non-trivial "structure" visible e.g. in the heat capacity, but that in general the largest radiative corrections originate from QCD effects, reducing the energy density by a couple of percent from the free value even at T > 160 GeV.

preprint2014arXiv

A relation between screening masses and real-time rates

Thermal screening masses related to the conserved vector current are determined for the case that the current carries a non-zero Matsubara frequency, both in a weak-coupling approach and through lattice QCD. We point out that such screening masses are sensitive to the same infrared physics as light-cone real-time rates. In particular, on the perturbative side, the inhomogeneous Schrodinger equation determining screening correlators is shown to have the same general form as the equation implementing LPM resummation for the soft-dilepton and photon production rates from a hot QCD plasma. The static potential appearing in the equation is identical to that whose soft part has been determined up to NLO and on the lattice in the context of jet quenching. Numerical results based on this potential suggest that screening masses overshoot the free results (multiples of 2piT) more strongly than at zero Matsubara frequency. Four-dimensional lattice simulations in two-flavour QCD at temperatures of 250 and 340 MeV confirm the non-static screening masses at the 10% level. Overall our results lend support to studies of jet quenching based on the same potential at T > 250 MeV.

preprint2014arXiv

Interpolation of hard and soft dilepton rates

Strict next-to-leading order (NLO) results for the dilepton production rate from a QCD plasma at temperatures above a few hundred MeV suffer from a breakdown of the loop expansion in the regime of soft invariant masses M^2 << (pi T)^2. In this regime an LPM resummation is needed for obtaining the correct leading-order result. We show how to construct an interpolation between the hard NLO and the leading-order LPM expression. The final numerical results are presented in a tabulated form, suitable for insertion into hydrodynamical codes.

preprint2014arXiv

Kubo relations and radiative corrections for lepton number washout

The rates for lepton number washout in extensions of the Standard Model containing right-handed neutrinos are key ingredients in scenarios for baryogenesis through leptogenesis. We relate these rates to real-time correlation functions at finite temperature, without making use of any particle approximations. The relations are valid to quadratic order in neutrino Yukawa couplings and to all orders in Standard Model couplings. They take into account all spectator processes, and apply both in the symmetric and in the Higgs phase of the electroweak theory. We use the relations to compute washout rates at next-to-leading order in g, where g denotes a Standard Model gauge or Yukawa coupling, both in the non-relativistic and in the relativistic regime. Even in the non-relativistic regime the parametrically dominant radiative corrections are only suppressed by a single power of g. In the non-relativistic regime radiative corrections increase the washout rate by a few percent at high temperatures, but they are of order unity around the weak scale and in the relativistic regime.

preprint2014arXiv

Right-handed neutrino production rate at T > 160 GeV

The production rate of right-handed neutrinos from a Standard Model plasma at a temperature above a hundred GeV has previously been evaluated up to NLO in Standard Model couplings (g ~ 2/3) in relativistic (M ~ pi T) and non-relativistic regimes (M >> pi T), and up to LO in an ultrarelativistic regime (M < gT). The last result necessitates an all-orders resummation of the loop expansion, accounting for multiple soft scatterings of the nearly light-like particles participating in 1 <-> 2 reactions. In this paper we suggest how the regimes can be interpolated into a result applicable for any right-handed neutrino mass and at all temperatures above 160 GeV. The results can also be used for determining the lepton number washout rate in models containing right-handed neutrinos. Numerical results are given in a tabulated form permitting for their incorporation into leptogenesis codes. We note that due to effects from soft Higgs bosons there is a narrow intermediate regime around M ~ g^{1/2} T in which our interpolation is phenomenological and a more precise study would be welcome.

preprint2014arXiv

Vector screening masses in the quark-gluon plasma and their physical significance

Static and non-static thermal screening states that couple to the conserved vector current are investigated in the high-temperature phase of QCD. Their masses and couplings to the current are determined at weak coupling, as well as using two-flavor lattice QCD simulations. A consistent picture emerges from the comparison, providing evidence that non-static Matsubara modes can indeed be treated perturbatively. We elaborate on the physical significance of the screening masses.

preprint2013arXiv

Charm mass effects in bulk channel correlations

The bulk viscosity of thermalized QCD matter at temperatures above a few hundred MeV could be significantly influenced by charm quarks because their contribution arises four perturbative orders before purely gluonic effects. In an attempt to clarify the challenges of a lattice study, we determine the relevant imaginary-time correlator (of massive scalar densities) up to NLO in perturbation theory, and compare with existing data. We find discrepancies much larger than in the vector channel; this may hint, apart from the importance of taking a continuum limit, to larger non-perturbative effects in the scalar channel. We also recall how a transport peak related to the scalar density spectral function encodes non-perturbative information concerning the charm quark chemical equilibration rate close to equilibrium.

preprint2013arXiv

Dimensionally regularized Polyakov loop correlators in hot QCD

A popular observable in finite-temperature lattice QCD is the so-called singlet quark-antiquark free energy, conventionally defined in Coulomb gauge. In an effort to interpret the existing numerical data on this observable, we compute it at order O(alpha_s^2) in continuum, and analyze the result at various distance scales. At short distances (r << 1/pi T) the behaviour matches that of the gauge-independent zero-temperature potential; on the other hand at large distances (r >> 1/pi T) the singlet free energy appears to have a gauge-fixing related power-law tail. At infinite distance the result again becomes physical in the sense that it goes over to a gauge-independent disconnected contribution, the square of the expectation value of the trace of the Polyakov loop; we recompute this quantity at O(alpha_s^2), finding for pure SU(Nc) a different non-logarithmic term than in previous literature, and adding for full QCD the quark contribution. We also discuss the value of the singlet free energy in a general covariant gauge, as well as the behaviour of the cyclic Wilson loop that is obtained if the singlet free energy is made gauge-independent by inserting straight spacelike Wilson lines into the observable. Comparisons with lattice data are carried out where possible.

preprint2013arXiv

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 m_h approximately 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.

preprint2013arXiv

Light-cone Wilson loop in classical lattice gauge theory

The transverse broadening of an energetic jet passing through a non-Abelian plasma is believed to be described by the thermal expectation value of a light-cone Wilson loop. In this exploratory study, we measure the light-cone Wilson loop with classical lattice gauge theory simulations. We observe, as suggested by previous studies, that there are strong interactions already at short transverse distances, which may lead to more efficient jet quenching than in leading-order perturbation theory. We also verify that the asymptotics of the Wilson loop do not change qualitatively when crossing the light cone, which supports arguments in the literature that infrared contributions to jet quenching can be studied with dimensionally reduced simulations in the space-like domain. Finally we speculate on possibilities for full four-dimensional lattice studies of the same observable, perhaps by employing shifted boundary conditions in order to simulate ensembles boosted by an imaginary velocity.

preprint2013arXiv

NLO thermal dilepton rate at non-zero momentum

The vector channel spectral function and the dilepton production rate from a QCD plasma at a temperature above a few hundred MeV are evaluated up to next-to-leading order (NLO) including their dependence on a non-zero momentum with respect to the heat bath. The invariant mass of the virtual photon is taken to be in the range K^2 ~ (pi T)^2 ~ (1 GeV)^2, generalizing previous NLO results valid for K^2 >> (pi T)^2. In the opposite regime 0 < K^2 << (pi T)^2 the loop expansion breaks down, but agrees nevertheless in order of magnitude with a previous result obtained through resummations. Ways to test the vector spectral function through comparisons with imaginary-time correlators measured on the lattice are discussed.

preprint2013arXiv

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.

preprint2013arXiv

Thermal 2-loop master spectral function at finite momentum

When considering NLO corrections to thermal particle production in the "relativistic" regime, in which the invariant mass squared of the produced particle is K^2 ~ (pi T)^2, then the production rate can be expressed as a sum of a few universal "master" spectral functions. Taking the most complicated 2-loop master as an example, a general strategy for obtaining a convergent 2-dimensional integral representation is suggested. The analysis applies both to bosonic and fermionic statistics, and shows that for this master the non-relativistic approximation is only accurate for K^2 > (8 pi T)^2, whereas the zero-momentum approximation works surprisingly well. Once the simpler masters have been similarly resolved, NLO results for quantities such as the right-handed neutrino production rate from a Standard Model plasma or the dilepton production rate from a QCD plasma can be assembled for K^2 ~ (pi T)^2.

preprint2013arXiv

Thermal right-handed neutrino production rate in the relativistic regime

The production rate of right-handed neutrinos from a Standard Model plasma at a temperature above a hundred GeV is evaluated up to NLO in Standard Model couplings. The results apply in the so-called relativistic regime, referring parametrically to a mass M ~ pi T, generalizing thereby previous NLO results which only apply in the non-relativistic regime M >> pi T. The non-relativistic expansion is observed to converge for M > 15 T, but the smallness of any loop corrections allows it to be used in practice already for M > 4 T. In the latter regime any non-covariant dependence of the differential rate on the spatial momentum is shown to be mild. The loop expansion breaks down in the ultrarelativistic regime M << pi T, but after a simple mass resummation it nevertheless extrapolates reasonably well towards a result obtained previously through complete LPM resummation, apparently confirming a strong enhancement of the rate at high temperatures (which facilitates chemical equilibration). When combined with other ingredients the results may help to improve upon the accuracy of leptogenesis computations operating above the electroweak scale.

preprint2013arXiv

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 r0*Tc ~ 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.

preprint2013arXiv

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_perp), characterizing the broadening of a high-energy jet. First results, based on data produced recently by Panero et al, suggest that C(k_perp) is enhanced over the known NLO result in a soft regime k_perp < a few T. The shape of k_perp^3 C(k_perp) is consistent with a Gaussian at small k_perp.

preprint2012arXiv

A non-perturbative contribution to jet quenching

It has been argued by Caron-Huot that infrared contributions to the jet quenching parameter in hot QCD, denoted by qhat, can be extracted from an analysis of a certain static-potential related observable within the dimensionally reduced effective field theory. Following this philosophy, the order of magnitude of a non-perturbative contribution to qhat from the colour-magnetic scale, g^2T/pi, is estimated. The result is small; it is probably below the parametrically perturbative but in practice slowly convergent contributions from the colour-electric scale, whose all-orders resummation therefore remains an important challenge.

preprint2012arXiv

Heavy quark chemical equilibration rate as a transport coefficient

Motivated by indications that heavy (charm and bottom) quarks interact strongly at temperatures generated in heavy ion collision experiments, we suggest a non-perturbative definition of a heavy quark chemical equilibration rate as a transport coefficient. Within leading-order perturbation theory (corresponding to 3-loop level), the definition is argued to reduce to an expression obtained from the Boltzmann equation. Around T ~ 400 MeV, an order-of-magnitude estimate for charm yields a rate Gamma^{-1}_{chem} > 60 fm/c which remains too slow to play a practical role in current experiments. However, the rate increases rapidly with T and, due to non-linear effects, also if the initial state contains an overabundance of heavy quarks.

preprint2012arXiv

Massive vector current correlator in thermal QCD

We present an NLO analysis of the massive vector current correlator at temperatures above a few hundred MeV. The physics of this correlator originates from a transport peak, related to heavy quark diffusion, and from the quark-antiquark threshold, related to quarkonium physics. In the bottom case both can be studied with separate effective theories, but for charm these may not be accurate, so a study within the full theory is needed. Working in imaginary time, the NLO correlator can be computed in unresummed perturbation theory; comparing with lattice data, we find good agreement. Subsequently we inspect how non-perturbative modifications of the transport peak would affect the imaginary-time correlator. The massive NLO quark-number susceptibility is also contrasted with numerical measurement.

preprint2012arXiv

Sommerfeld effect in heavy quark chemical equilibration

The chemical equilibration of heavy quarks in a quark-gluon plasma proceeds via annihilation or pair creation. For temperatures T much below the heavy quark mass M, when kinetically equilibrated heavy quarks move very slowly, the annihilation in the colour singlet channel is enhanced because the quark and antiquark attract each other which increases their probability to meet, whereas the octet contribution is suppressed. This is the so-called Sommerfeld effect. It has not been taken into account in previous calculations of the chemical equilibration rate, which are therefore incomplete for T < alphas^2 M. We compute the leading-order equilibration rate in this regime; there is a large enhancement in the singlet channel, but the rate is dominated by the octet channel, and therefore the total effect is small. In the course of the computation we demonstrate how operators that represent the annihilation of heavy quarks in non-relativistic QCD can be incorporated into the imaginary-time formalism.

preprint2012arXiv

Thermal right-handed neutrino production rate in the non-relativistic regime

We consider the next-to-leading order thermal production rate of heavy right-handed neutrinos in the non-relativistic regime m_top < pi*T << M, where m_top refers to the electroweak scale. Rephrasing the problem in an OPE language and making use of different techniques than a previous analysis by Salvio et al, we confirm the general structure of their result and many of the coefficients. We also extend the analysis to the next order in the non-relativistic expansion, thereby revealing the leading non-trivial momentum dependence, as well as to NNLO in couplings, revealing the leading sensitivity to thermal resummations. Our results are expressed as a sum of simple "master" structures, which renders them a suitable starting point for determining the next-to-leading order rate also in the relativistic regime pi*T ~ M.

preprint2012arXiv

Thermal right-handed neutrino self-energy in the non-relativistic regime

Recently the issue of radiative corrections to leptogenesis has been raised. Considering the "strong washout" regime, in which OPE-techniques permit to streamline the setup, we report the thermal self-energy matrix of heavy right-handed neutrinos at NLO (resummed 2-loop level) in Standard Model couplings. The renormalized expression describes flavour transitions and "inclusive" decays of chemically decoupled right-handed neutrinos. Although CP-violation is not addressed, the result may find use in existing leptogenesis frameworks.

preprint2012arXiv

Towards flavour diffusion coefficient and electrical conductivity without ultraviolet contamination

By subtracting from a recent lattice measurement of the thermal vector-current correlator the known 5-loop vacuum contribution, we demonstrate that the remainder is small and shows no visible short-distance divergence. It can therefore in principle be subjected to model-independent analytic continuation. Testing a particular implementation, we obtain estimates for the flavour-diffusion coefficient (2 pi T D \gsim 0.8) and electrical conductivity which are significantly smaller than previous results. Although systematic errors remain beyond control at present, some aspects of our approach could be of a wider applicability.

preprint2011arXiv

A test on analytic continuation of thermal imaginary-time data

Some time ago, Cuniberti et al have proposed a novel method for analytically continuing thermal imaginary-time correlators to real time, which requires no model input and should be applicable with finite-precision data as well. Given that these assertions go against common wisdom, we report on a naive test of the method with an idealized example. We do encounter two problems, which we spell out in detail; this implies that systematic errors are difficult to quantify. On a more positive note, the method is simple to implement and allows for an empirical recipe by which a reasonable qualitative estimate for some transport coefficient may be obtained, if statistical errors of an ultraviolet-subtracted imaginary-time measurement can be reduced to roughly below the per mille level.

preprint2011arXiv

Heavy flavour kinetic equilibration in the confined phase

By making use of a non-perturbative definition of a momentum diffusion coefficient as well as Heavy Meson Chiral Perturbation Theory, we investigate the Brownian motion and kinetic equilibration of heavy quark flavours deep in the confined phase. It appears that the momentum diffusion coefficient can be expressed in terms of known low-energy constants; it increases rapidly at temperatures above 50 MeV, behaving as ~ T^7/F_pi^4 for m_pi / pi << T << F_pi, where m_pi and F_pi are the pion mass and decay constant, respectively. The early increase may suggest a broad peak in kappa/T^3 around the QCD crossover. For a more detailed understanding the computation could be generalized in a number of ways.

preprint2011arXiv

Next-to-leading order thermal spectral functions in the perturbative domain

Motivated by applications in thermal QCD and cosmology, we elaborate on a general method for computing next-to-leading order spectral functions for composite operators at vanishing spatial momentum, accounting for real, virtual as well as thermal corrections. As an example, we compute these functions (together with the corresponding imaginary-time correlators which can be compared with lattice simulations) for scalar and pseudoscalar densities in pure Yang-Mills theory. Our results may turn out to be helpful in non-perturbative estimates of the corresponding transport coefficients, which are the bulk viscosity in the scalar channel and the rate of anomalous chirality violation in the pseudoscalar channel. We also mention links to cosmology, although the most useful results in that context may come from a future generalization of our methods to other correlators.

preprint2011arXiv

Towards a non-perturbative measurement of the heavy quark momentum diffusion coefficient

We report on a lattice investigation of heavy quark diffusion within pure SU(3) plasma above the deconfinement transition, with the quarks treated to leading order in the heavy mass expansion. Using a multilevel algorithm, several volumes and lattice spacings, as well as tree-level improvement and perturbative renormalization, we measure the relevant "colour-electric" Euclidean correlator, finding that it clearly exceeds its perturbative counterpart. Even without analytic continuation, this suggests that at temperatures just above the critical one, non-perturbative interactions felt by the heavy quarks are stronger than within the weak-coupling expansion. After introducing rough modelling of the spectral shape, diffusion coefficients down to D ~ 0.5/T appear possible.

preprint2010arXiv

Colour-electric spectral function at next-to-leading order

The spectral function related to the correlator of two colour-electric fields along a Polyakov loop determines the momentum diffusion coefficient of a heavy quark near rest with respect to a heat bath. We compute this spectral function at next-to-leading order, O(alpha_s^2), in the weak-coupling expansion. The high-frequency part of our result (omega >> T), which is shown to be temperature-independent, is accurately determined thanks to asymptotic freedom; the low-frequency part of our result (omega << T), in which Hard Thermal Loop resummation is needed in order to cure infrared divergences, agrees with a previously determined expression. Our result may help to calibrate the overall normalization of a lattice-extracted spectral function in a perturbative frequency domain T << omega << 1/a, paving the way for a non-perturbative estimate of the momentum diffusion coefficient at omega -> 0. We also evaluate the colour-electric Euclidean correlator, which could be directly compared with lattice simulations. As an aside we determine the Euclidean correlator in the lattice strong-coupling expansion, showing that through a limiting procedure it can in principle be defined also in the confined phase of pure Yang-Mills theory, even if a practical measurement could be very noisy there.

preprint2010arXiv

Intermediate distance correlators in hot Yang-Mills theory

Lattice measurements of spatial correlation functions of the operators FF and FF-dual in thermal SU(3) gauge theory have revealed a clear difference between the two channels at "intermediate" distances, x ~ 1/(pi T). This is at odds with the AdS/CFT limit which predicts the results to coincide. On the other hand, an OPE analysis at short distances (x << 1/(pi T)) as well as effective theory methods at long distances (x >> 1/(pi T)) suggest differences. Here we study the situation at intermediate distances by determining the time-averaged spatial correlators through a 2-loop computation. We do find unequal results, however the numerical disparity is small. Apart from theoretical issues, a future comparison of our results with time-averaged lattice measurements might also be of phenomenological interest in that understanding the convergence of the weak-coupling series at intermediate distances may bear on studies of the thermal broadening of heavy quarkonium resonances.

preprint2010arXiv

On bulk viscosity and moduli decay

This pedagogically intended lecture, one of four under the header "Basics of thermal QCD", reviews an interesting relationship, originally pointed out by Bodeker, that exists between the bulk viscosity of Yang-Mills theory (of possible relevance to the hydrodynamics of heavy ion collision experiments) and the decay rate of scalar fields coupled very weakly to a heat bath (appearing in some particle physics inspired cosmological scenarios). This topic serves, furthermore, as a platform on which a number of generic thermal field theory concepts are illustrated. The other three lectures (on the QCD equation of state and the rates of elastic as well as inelastic processes experienced by heavy quarks) are recapitulated in brief encyclopedic form.

preprint2010arXiv

Ultraviolet asymptotics of scalar and pseudoscalar correlators in hot Yang-Mills theory

Inspired by recent lattice measurements, we determine the short-distance (a << r << 1/pi T) as well as large-frequency (1/a >> omega >> pi T) asymptotics of scalar (trace anomaly) and pseudoscalar (topological charge density) correlators at 2-loop order in hot Yang-Mills theory. The results are expressed in the form of an Operator Product Expansion. We confirm and refine the determination of a number of Wilson coefficients; however some discrepancies with recent literature are detected as well, and employing the correct values might help, on the qualitative level, to understand some of the features observed in the lattice measurements. On the other hand, the Wilson coefficients show slow convergence and it appears uncertain whether this approach can lead to quantitative comparisons with lattice data. Nevertheless, as we outline, our general results might serve as theoretical starting points for a number of perhaps phenomenologically more successful lines of investigation.

preprint2009arXiv

Quarkonium dissociation in the presence of a small momentum space anisotropy

We consider the dissociation of heavy quarkonium in a medium close to thermal equilibrium but with a small momentum space anisotropy. Dissociation is defined to take place when the width of the ground state equals its binding energy. We show that if the anisotropic medium is obtained isentropically from the equilibrium one, then to first order in the anisotropy parameter the dissociation temperature remains unchanged. If, in contrast, the non-equilibrium system has a smaller entropy density than the equilibrium one, then the dissociation temperature increases with respect to the isotropic case, by up to 10% for modest anisotropies.

preprint2007arXiv

Four-loop pressure of massless O(N) scalar field theory

Inspired by the corresponding problem in QCD, we determine the pressure of massless O(N) scalar field theory up to order g^6 in the weak-coupling expansion, where g^2 denotes the quartic coupling constant. This necessitates the computation of all 4-loop vacuum graphs at a finite temperature: by making use of methods developed by Arnold and Zhai at 3-loop level, we demonstrate that this task is manageable at least if one restricts to computing the logarithmic terms analytically, while handling the ``constant'' 4-loop contributions numerically. We also inspect the numerical convergence of the weak-coupling expansion after the inclusion of the new terms. Finally, we point out that while the present computation introduces strategies that should be helpful for the full 4-loop computation on the QCD-side, it also highlights the need to develop novel computational techniques, in order to be able to complete this formidable task in a systematic fashion.

preprint2007arXiv

Real-time static potential in hot QCD

We derive a static potential for a heavy quark-antiquark pair propagating in Minkowski time at finite temperature, by defining a suitable gauge-invariant Green's function and computing it to first non-trivial order in Hard Thermal Loop resummed perturbation theory. The resulting Debye-screened potential could be used in models that attempt to describe the ``melting'' of heavy quarkonium at high temperatures. We show, in particular, that the potential develops an imaginary part, implying that thermal effects generate a finite width for the quarkonium peak in the dilepton production rate. For quarkonium with a very heavy constituent mass M, the width can be ignored for T \lsim g^2 M/12π, where g^2 is the strong gauge coupling; for a physical case like bottomonium, it could become important at temperatures as low as 250 MeV. Finally, we point out that the physics related to the finite width originates from the Landau-damping of low-frequency gauge fields, and could be studied non-perturbatively by making use of the classical approximation.

preprint2006arXiv

Probing the chiral weak Hamiltonian at finite volumes

Non-leptonic kaon decays are often described through an effective chiral weak Hamiltonian, whose couplings ("low-energy constants") encode all non-perturbative QCD physics. It has recently been suggested that these low-energy constants could be determined at finite volumes by matching the non-perturbatively measured three-point correlation functions between the weak Hamiltonian and two left-handed flavour currents, to analytic predictions following from chiral perturbation theory. Here we complete the analytic side in two respects: by inspecting how small ("epsilon-regime") and intermediate or large ("p-regime") quark masses connect to each other, and by including in the discussion the two leading Delta I = 1/2 operators. We show that the epsilon-regime offers a straightforward strategy for disentangling the coefficients of the Delta I = 1/2 operators, and that in the p-regime finite-volume effects are significant in these observables once the pseudoscalar mass M and the box length L are in the regime ML \lsim 5.0.

preprint2006arXiv

Quark mass thresholds in QCD thermodynamics

We discuss radiative corrections to how quark mass thresholds are crossed, as a function of the temperature, in basic thermodynamic observables such as the pressure, the energy and entropy densities, and the heat capacity of high temperature QCD. The indication from leading order that the charm quark plays a visible role at surprisingly low temperatures, is confirmed. We also sketch a way to obtain phenomenological estimates relevant for generic expansion rate computations at temperatures between the QCD and electroweak scales, pointing out where improvements over the current knowledge are particularly welcome.

preprint2005arXiv

On the determination of low-energy constants for $ΔS=1$ transitions

We present our preliminary results for three-point correlation functions involving the operators entering the $Δ{S}=1$ effective Hamiltonian with an active charm quark, obtained using overlap fermions in the quenched approximation. This is the first computation carried out for valence quark masses small enough so as to permit a matching to Quenched Chiral Perturbation Theory in the $ε$-regime. The commonly observed large statistical fluctuations are tamed by means of low-mode averaging techniques, combined with restrictions to individual topological sectors. We also discuss the matching of the resulting hadronic matrix elements to the effective low-energy constants for $Δ{S}=1$ transitions. This involves (a) finite-volume corrections which can be evaluated at NLO in Quenched Chiral Perturbation Theory, and (b) the short-distance renormalization of the relevant four-quark operators in discretizations based on the overlap operator. We discuss perturbative estimates for the renormalization factors and possible strategies for their non-perturbative evaluation. Our results can be used to isolate the long-distance contributions to the $ΔI=1/2$ rule, coming from physics effects around the intrinsic QCD scale.

preprint2005arXiv

Real-time Chern-Simons term for hypermagnetic fields

If non-vanishing chemical potentials are assigned to chiral fermions, then a Chern-Simons term is induced for the corresponding gauge fields. In thermal equilibrium anomalous processes adjust the chemical potentials such that the coefficient of the Chern-Simons term vanishes, but it has been argued that there are non-equilibrium epochs in cosmology where this is not the case and that, consequently, certain fermionic number densities and large-scale (hypermagnetic) field strengths get coupled to each other. We generalise the Chern-Simons term to a real-time situation relevant for dynamical considerations, by deriving the anomalous Hard Thermal Loop effective action for the hypermagnetic fields, write down the corresponding equations of motion, and discuss some exponentially growing solutions thereof.

preprint2004arXiv

Low-energy couplings of QCD from current correlators near the chiral limit

We investigate a new numerical procedure to compute fermionic correlation functions at very small quark masses. Large statistical fluctuations, due to the presence of local ``bumps'' in the wave functions associated with the low-lying eigenmodes of the Dirac operator, are reduced by an exact low-mode averaging. To demonstrate the feasibility of the technique, we compute the two-point correlator of the left-handed vector current with Neuberger fermions in the quenched approximation, for lattices with a linear extent of L~1.5 fm, a lattice spacing a~0.09 fm, and quark masses down to the epsilon-regime. By matching the results with the corresponding (quenched) chiral perturbation theory expressions, an estimate of (quenched) low-energy constants can be obtained. We find agreement between the quenched values of F extrapolated from the p-regime and extracted in the epsilon-regime.

preprint2004arXiv

Plaquette expectation value and gluon condensate in three dimensions

In three dimensions, the gluon condensate of pure SU(3) gauge theory has ultraviolet divergences up to 4-loop level only. By subtracting the corresponding terms from lattice measurements of the plaquette expectation value and extrapolating to the continuum limit, we extract the finite part of the gluon condensate in lattice regularization. Through a change of regularization scheme to MSbar and (inverse) dimensional reduction, this result would determine the first non-perturbative coefficient in the weak-coupling expansion of hot QCD pressure.

preprint2003arXiv

Finite-Size Scaling of Vector and Axial Current Correlators

Using quenched chiral perturbation theory, we compute the long-distance behaviour of two-point functions of flavour non-singlet axial and vector currents in a finite volume, for small quark masses, and at a fixed gauge-field topology. We also present the corresponding predictions for the unquenched theory at fixed topology. These results can in principle be used to measure the low-energy constants of the chiral Lagrangian, from lattice simulations in volumes much smaller than one pion Compton wavelength. We show that quenching has a dramatic effect on the vector correlator, which is argued to vanish to all orders, while the axial correlator appears to be a robust observable only moderately sensitive to quenching.

preprint2001arXiv

Critical behaviour of the Ginzburg-Landau model in the type II region

We study the critical behaviour of the three-dimensional U(1) gauge+Higgs theory (Ginzburg-Landau model) at large scalar self-coupling λ(``type II region'') by measuring various correlation lengths as well as the Abrikosov-Nielsen-Olesen vortex tension. We identify different scaling regions as the transition is approached from below, and carry out detailed comparisons with the criticality of the 3d O(2) symmetric scalar theory.

preprint2001arXiv

Finite baryon density effects on gauge field dynamics

We discuss the effective action for QCD gauge fields at finite temperatures and densities, obtained after integrating out the hardest momentum scales from the system. We show that a non-vanishing baryon density induces a charge conjugation (C) odd operator to the gauge field action, proportional to the chemical potential. Even though it is parametrically smaller than the leading C even operator, it could have an important effect on C odd observables. The same operator appears to be produced by classical kinetic theory, allowing in principle for a non-perturbative study of such processes.

preprint2001arXiv

Measuring infrared contributions to the QCD pressure

For the pressure (or free energy) of QCD, four-dimensional (4d) lattice data is available at zero baryon density up to a few times the critical temperature $T_c$. Perturbation theory, on the other hand, has serious convergence problems even at very high temperatures. In a combined analytical and three-dimensional (3d) lattice method, we show that it is possible to compute the QCD pressure from about $2 T_c$ to infinity. The numerical accuracy is good enough to resolve in principle, e.g., logarithmic contributions related to 4-loop perturbation theory.

preprint2000arXiv

Vortex phases in condensed matter and cosmology

Placing a high-Tc superconductor in an increasing external magnetic field, the flux first penetrates the sample through an Abrikosov vortex lattice, and then a first order transition is observed by which the system goes to the normal phase. We discuss the cosmological motivation for considering the electroweak phase transition in the presence of an external magnetic field, the analogies this system might have with the superconductor behaviour described above, and in particular whether at large physical Higgs masses, corresponding to the high-Tc regime, an analogue of the vortex phase and an associated first order phase transition could be generated.

preprint1999arXiv

CP Violating Bubble Wall Profiles

We solve the equations of motion for a CP violating phase between the two Higgs doublets at the bubble wall of the MSSM electroweak phase transition. Contrary to earlier suggestions, we do not find indications of spontaneous ``transitional'' CP violation in the MSSM. On the other hand, in case there is explicit CP violation in the stop and chargino/neutralino sectors, the relative phase between the Higgses does become space dependent, but only mildly even in the maximal case. We also demonstrate that spontaneous CP violation within the bubble wall could occur, e.g., if the Higgs sector of the MSSM were supplemented by a singlet. Finally we point out some implications for baryogenesis computations.

preprint1999arXiv

Effective Theories of MSSM at High Temperature

We construct effective 3d field theories for the Minimal Supersymmetric Standard Model, relevant for the thermodynamics of the cosmological electroweak phase transition. The effective theories include a 3d theory for the bosonic sector of the original 4d theory; a 3d two Higgs doublet model; and a 3d SU(2)+Higgs model. The integrations are made at 1-loop level. In integrals related to vacuum renormalization we take into account only quarks and squarks of the third generation. Using existing non-perturbative lattice results for the 3d SU(2)+Higgs model, we then derive infrared safe upper bounds for the lightest Higgs boson mass required for successful baryogenesis at the electroweak scale. The Higgs mass bounds turn out to be close to those previously found with the effective potential, allowing baryogenesis if the right-handed stop mass parameter $m_U^2$ is small. Finally we discuss the effective theory relevant for $m_U^2$ very small, the most favourable case for baryogenesis.

preprint1999arXiv

Spontaneous CP violation on the lattice

At finite temperatures around the electroweak phase transition, the thermodynamics of the MSSM can be described by a three-dimensional two Higgs doublet effective theory. This effective theory has a phase where CP is spontaneously violated. We study spontaneous CP violation with non-perturbative lattice simulations, and analyse whether one could end up in this phase for any physical MSSM parameter values.

preprint1998arXiv

A Strong Electroweak Phase Transition up to m_H ~ 105 GeV

Non-perturbative lattice simulations have shown that there is no electroweak phase transition in the Standard Model for the allowed Higgs masses, m_H \gsim 75 GeV. In the Minimal Supersymmetric Standard Model, in contrast, it has been proposed that the transition should exist and even be strong enough for baryogenesis up to m_H ~ 105 GeV, provided that the lightest stop mass is in the range 100...160 GeV. However, this prediction is based on perturbation theory, and suffers from a noticeable gauge parameter and renormalization scale dependence. We have performed large-scale lattice Monte Carlo simulations of the MSSM electroweak phase transition. Extrapolating the results to the infinite volume and continuum limits, we find that the transition is in fact stronger than indicated by 2-loop perturbation theory. This guarantees that the perturbative Higgs mass bound m_H ~ 105 GeV is a conservative one, allows slightly larger stop masses (up to ~ 165 GeV), and provides a strong motivation for further studies of MSSM electroweak baryogenesis.

preprint1998arXiv

The Phase Diagram of Three-Dimensional SU(3) + Adjoint Higgs Theory

We study the phase diagram of the three-dimensional SU(3)+adjoint Higgs theory with lattice Monte Carlo simulations. A critical line consisting of a first order line, a tricritical point and a second order line, divides the phase diagram into two parts distinguished by <Tr A0^3>=0 and /=0. The location and the type of the critical line are determined by measuring the condensates <Tr A0^2> and <Tr A0^3>, and the masses of scalar and vector excitations. Although in principle there can be different types of broken phases, corresponding perturbatively to unbroken SU(2)xU(1) or U(1)xU(1) symmetries, we find that dynamically only the broken phase with SU(2)xU(1)-like properties is realized. The relation of the phase diagram to 4d finite temperature QCD is discussed.

preprint1998arXiv

What's new with the electroweak phase transition?

We review the status of non-perturbative lattice studies of the electroweak phase transition. In the Standard Model, the complete phase diagram has been reliably determined, and the conclusion is that there is no phase transition at all for the experimentally allowed Higgs masses. In the Minimal Supersymmetric Standard Model (MSSM), in contrast, there can be a strong first order transition allowing for baryogenesis. Finally, we point out possibilities for future simulations, such as the problem of CP-violation at the MSSM electroweak phase boundary.

preprint1997arXiv

Non-perturbative Debye mass in finite T QCD

Employing a non-perturbative gauge invariant definition of the Debye screening mass m_D in the effective field theory approach to finite T QCD, we use 3d lattice simulations to determine the leading O(g^2) and to estimate the next-to-leading O(g^3) corrections to m_D in the high temperature region. The O(g^2) correction is large and modifies qualitatively the standard power-counting hierarchy picture of correlation lengths in high temperature QCD.

preprint1997arXiv

The finite temperature real time \hbar^2 corrections in quantum mechanics

We study non-perturbative real time correlation functions at finite temperature. In order to see whether the classical term gives a good approximation in the high temperature limit T >> \hbarω, we consider the first \hbar^2 quantum corrections. We find that for the simplest non-trivial case, the quantum mechanical anharmonic oscillator, the classical result is reliable only for moderately large times: after some time t_* the classical approximation breaks down even at high temperatures. Moreover, the result for the first quantum corrections cannot, in general, be reproduced by modifying the parameters of the classical theory.

preprint1996arXiv

Is There a Hot Electroweak Phase Transition at $m_H\gsim m_W$?

We provide non-perturbative evidence for the fact that there is no hot electroweak phase transition at large Higgs masses, $m_H = 95$, 120 and 180 GeV. This means that the line of first order phase transitions separating the symmetric and broken phases at small $m_H$ has an end point $m_{H,c}$. In the minimal standard electroweak theory 70 GeV $<m_{H,c}<$ 95 GeV and most likely $m_{H,c} \approx 80$ GeV. If the electroweak theory is weakly coupled and the Higgs boson is found to be heavier than the critical value (which depends on the theory in question), cosmological remnants from the electroweak epoch are improbable.

preprint1995arXiv

Supersonic Deflagrations in Cosmological Phase Transitions

The classification of the hydrodynamical growth mechanisms for the spherical bubbles of the low-temperature phase in cosmological phase transitions is completed by showing that the bubbles can grow as supersonic deflagrations. Such deflagrations consist of a Jouguet deflagration, followed by a rarefaction wave. Depending on the amount of supercooling, the maximal velocity of supersonic deflagrations varies between the sound and the light velocities. The solutions faster than supersonic deflagrations are weak detonations.

preprint1994arXiv

Gauge Dependence of the High-Temperature 2-Loop Effective Potential for the Higgs Field

The high-temperature limit of the 2-loop effective potential for the Higgs field is calculated from an effective 3d theory, in a general covariant gauge. It is shown explicitly that a gauge-independent result can be extracted for the equation of state from the gauge-dependent effective potential. The convergence of perturbation theory is estimated in the broken phase, utilizing the gauge dependence of the effective potential.