No Arabic abstract
We study how to numerically simulate quantum fermions out of thermal equilibrium, in the context of electroweak baryogenesis. We find that by combining the lattice implementation of Aarts and Smit [1] with the low cost fermions of Borsanyi and Hindmarsh [2], we are able to describe the dynamics of a classical bosonic system coupled to quantum fermions, that correctly reproduces anomalous baryon number violation. To demonstrate the method, we apply it to the 1+1 dimensional axial U(1) model, and perform simulations of a fast symmetry breaking transition. Compared to solving all the quantum mode equations as in [1], we find that this statistical approach may lead to a significant gain in computational time, when applied to 3+1 dimensional physics.
We study the spectral properties of a highly occupied non-Abelian non-equilibrium plasma appearing ubiquitously in weak coupling descriptions of QCD matter. The spectral function of this far-from-equilibrium plasma is measured by employing linear response theory in classical-statistical real-time lattice Yang-Mills simulations. We establish the existence of transversely and longitudinally polarized quasiparticles and obtain their dispersion relations, effective mass, plasmon frequency, damping rate and further structures in the spectral and statistical functions. Our new method can be interpreted as a non-perturbative generalization of hard thermal loop (HTL) effective theory. We see indications that our results approach leading order HTL in the appropriate limit. The method can also be employed beyond the range of validity of HTL.
We present a new method to obtain spectral properties of a non-Abelian gauge theory in the region where occupation numbers are high. The method to measure the (single-particle) spectral function is based on linear response theory and classical-statistical lattice simulations. Although we apply it to a system far from equilibrium in a self-similar regime, the extracted spectral function can be understood within the hard thermal loop (HTL) formalism and can thus be connected to thermal equilibrium at high temperatures. This allows us to obtain quantities like the lifetime of quasiparticles that are beyond the leading order and difficult to compute within HTL. The approach has the potential to measure transport coefficients, to study the earliest stages of heavy-ion collisions in a controlled way and it can be employed beyond the range of validity of HTL.
We compute the baryon asymmetry created in a tachyonic electroweak symmetry breaking transition, focusing on the dependence on the source of effective CP-violation. Earlier simulations of Cold Electroweak Baryogenesis have almost exclusively considered a very specific CP-violating term explicitly biasing Chern-Simons number. We compare four different dimension six, scalar-gauge CP-violating terms, involving both the Higgs field and another dynamical scalar coupled to SU(2) or U(1) gauge fields. We find that for sensible values of parameters, all implementations can generate a baryon asymmetry consistent with observations, showing that baryogenesis is a generic outcome of a fast tachyonic electroweak transition.
We perform numerical simulations of Cold Electroweak Baryogenesis, including for the first time in the Bosonic sector the full electroweak gauge group SU(2)$times$U(1) and CP-violation. We find that the maximum generated baryon asymmetry is reduced by a factor of three relative to the SU(2)-only model, but that the quench time dependence is very similar. In addition, we compute the magnitude of the helical magnetic fields, and find that it is proportional to the strength of CP-violation and dependent on quench time, but is not proportional to the magnitude of the baryon asymmetry as proposed in the literature. Astrophysical signatures of primordial magnetic helicity can therefore not in general be used as evidence that electroweak baryogenesis has taken place.
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.