Source author record

Aleksi Kurkela

Aleksi Kurkela 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

34works
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

34 published item(s)

preprint2022arXiv

Cold quark matter at N3LO: soft contributions

High-order perturbative calculations for thermodynamic quantities in QCD are complicated by the physics of dynamical screening that affects the soft, long-wavelength modes of the system. Here, we provide details for the evaluation of this soft contribution to the next-to-next-to-next-to-leading order (N3LO) pressure of high-density, zero-temperature quark matter (QM), complementing our accompanying paper in arXiv:2103.05658. Our calculation requires the determination of the pressure of the hard-thermal-loop (HTL) effective theory to full two-loop order at zero temperature, which we go through in considerable detail. In addition to this, we comprehensively discuss the structure of the weak-coupling expansion of the QM pressure, and lay out a roadmap towards the evaluation of the contributions missing from a full N3LO result for this quantity.

preprint2022arXiv

How perturbative QCD constrains the Equation of State at Neutron-Star densities

We demonstrate in a general and analytic way how high-density information about the equation of state (EoS) of strongly interacting matter obtained using perturbative Quantum Chromodynamics (pQCD) constrains the same EoS at densities reachable in physical neutron stars. Our approach is based on utilizing the full information of the thermodynamic potentials at the high-density limit together with thermodynamic stability and causality. This requires considering the pressure as a function of chemical potential $p(μ)$ instead of the commonly used pressure as a function of energy density $p(ε)$. The results can be used to propagate the pQCD calculations reliable around 40$n_s$ to lower densities in the most conservative way possible. We constrain the EoS starting from only few times the nuclear saturation density $n \gtrsim 2.2 n_s$ and at $n = 5 n_s$ we exclude at least 65% of otherwise allowed area in the $(ε- p)$-plane. This provides information complementary to astrophysical observations that should be taken into account in any complete statistical inference study of the EoS. These purely theoretical results are independent of astrophysical neutron-star input, and hence, they can also be used to test theories of modified gravity and BSM physics in neutron stars.

preprint2022arXiv

Multimessenger constraints for ultra-dense matter

Recent rapid progress in neutron-star (NS) observations offers great potential to constrain the properties of strongly interacting matter under the most extreme conditions. In order to fully exploit the current observational inputs and to study the impact of future observations of NS masses, radii, and tidal deformabilities, we analyze a large ensemble of randomly generated viable NS-matter equations of state (EoSs) and the corresponding rotating stellar structures. We discuss the compatibility and impact of various hypotheses and measurements on the EoS, including those involving the merger product of the gravitational-wave (GW) event GW170817, the binary merger components in GW190814, and radius measurements of the pulsar PSR J0740+6620. We obtain an upper limit for the dimensionless spin of a rigidly rotating NS, $|χ| < 0.81$, an upper limit for the compactness of a NS, $GM/(Rc^2) < 0.33$, and find that the conservative hypothesis that the remnant in GW170817 ultimately collapsed to a black hole strongly constrains the EoS and the maximal mass of NSs, implying $M_\text{TOV} < 2.53 M_\odot$ (or $M_\text{TOV} < 2.19 M_\odot$ if we assume that a hypermassive NS was created). Additionally, we derive a novel lower limit for the tidal deformability as a function of NS mass and provide fitting formulae that can be used to set priors for parameter estimation and to discern whether neutron stars or other compact objects are involved in future low-mass GW events. Finally, we find that the recent NICER results for the radius of the massive NS PSR J0740+6620 place strong constraints for the behavior of the EoS, and that the indicated radius values $R(2M_\odot) \gtrsim 11$ km are compatible with moderate speeds of sound in NS matter and thus with the existence of quark matter cores in massive NSs.

preprint2022arXiv

Non-equilibrium attractor in high-temperature QCD plasmas

We establish the existence of a far-from-equilibrium attractor in weakly-coupled gauge theory undergoing 0+1d Bjorken expansion which goes beyond the energy-momentum tensor to the detailed form of the one-particle distribution function. We then demonstrate that the dynamics can be re-scaled at intermediate times and represented by universal exponents. Finally, we assess different procedures for reconstructing the full one-particle distribution function from the energy-momentum tensor along the attractor and discuss implications for the freeze-out procedure used in the phenomenological analysis of ultra-relativistic nuclear collisions

preprint2022arXiv

Soft interactions in cold quark matter

Accurate knowledge of the thermodynamic properties of zero-temperature, high-density quark matter plays an integral role in attempts to constrain the behavior of the dense QCD matter found inside neutron-star cores, irrespective of the phase realized inside the stars. In this Letter, we consider the weak-coupling expansion of the dense QCD equation of state and compute the next-to-next-to-next-to-leading-order contribution arising from the non-Abelian interactions among long-wavelength, dynamically screened gluonic fields. Accounting for these interactions requires an all-loop resummation, which can be performed using hard-thermal-loop (HTL) kinematic approximations. Concretely, we perform a full two-loop computation using the HTL effective theory, valid for the long-wavelength, or soft, modes. We find that the soft sector is well-behaved within cold quark matter, contrary to the case encountered at high temperatures, and find that the new contribution decreases the renormalization-scale dependence of the equation of state at high density.

preprint2021arXiv

Thermalization of non-abelian gauge theories at next-to-leading order

We provide the first next-to-leading-order (NLO) weak-coupling description of the thermalization process of far-from-equilibrium systems in non-abelian gauge theory. We study isotropic systems starting from either over- or under-occupied initial conditions and follow their time evolution towards thermal equilibrium by numerically solving the QCD effective kinetic theory at NLO accuracy. We find that the NLO corrections remain well under control for a wide range of couplings and that the overall effect of NLO corrections is to reduce the time needed to reach thermal equilibrium in the systems considered.

preprint2020arXiv

Non-equilibrium attractor in high-temperature QCD plasmas

We establish the existence of a far-from-equilibrium attractor in weakly-coupled gauge theory undergoing one-dimensional Bjorken expansion. We demonstrate that the resulting far-from-equilibrium evolution is insensitive to certain features of the initial condition, including both the initial momentum-space anisotropy and initial occupancy. We find that this insensitivity extends beyond the energy-momentum tensor to the detailed form of the one-particle distribution function. Based on our results, we assess different procedures for reconstructing the full one-particle distribution function from the energy-momentum tensor along the attractor and discuss implications for the freeze-out procedure used in the phenomenological analysis of ultra-relativistic nuclear collisions.

preprint2019arXiv

Flow in AA and pA as an interplay of fluid-like and non-fluid like excitations

To study the microscopic structure of quark-gluon plasma, data from hadronic collisions must be confronted with models that go beyond fluid dynamics. Here, we study a simple kinetic theory model that encompasses fluid dynamics but contains also particle-like excitations in a boost invariant setting with no symmetries in the transverse plane and with large initial momentum asymmetries. We determine the relative weight of fluid dynamical and particle like excitations as a function of system size and energy density by comparing kinetic transport to results from the 0th, 1st and 2nd order gradient expansion of viscous fluid dynamics. We then confront this kinetic theory with data on azimuthal flow coefficients over a wide centrality range in PbPb collisions at the LHC, in AuAu collisions at RHIC, and in pPb collisions at the LHC. Evidence is presented that non-hydrodynamic excitations make the dominant contribution to collective flow signals in pPb collisions at the LHC and contribute significantly to flow in peripheral nucleus-nucleus collisions, while fluid-like excitations dominate collectivity in central nucleus-nucleus collisions at collider energies.

preprint2019arXiv

What attracts to attractors?

Whether, how, and to what extent solutions of Bjorken-expanding systems become insensitive to aspects of their initial conditions is of importance for heavy-ion collisions. Here we study 1+1D and phenomenologically relevant boost-invariant 3+1D systems in which initial conditions approach a universal attractor solution. In Israel-Stewart theory (IS) and kinetic theory where the universal attractor extends to arbitrarily early times, we show that all initial conditions approach the attractor at early times by a power-law while their approach is exponential at late times. In these theories, the physical mechanisms of hydrodynamization operational at late times do not drive the approach to the attractor at early times, and the early-time attractor is reached prior to hydrodynamization. In marked contrast, the attractor in strongly coupled systems is realized concurrent with hydrodynamization. This qualitative difference may offer a basis for discriminating weakly and strongly coupled scenarios of heavy-ion collisions.

preprint2016arXiv

Cool quark matter

We generalize the state-of-the-art perturbative Equation of State of cold quark matter to nonzero temperatures, needed in the description of neutron star mergers and core collapse processes. The new result is accurate to order g^5 in the gauge coupling, and is based on a novel framework for dealing with the infrared sensitive soft field modes of the theory. The zero Matsubara mode sector is treated using a dimensionally reduced effective theory, while the soft non-zero modes are resummed using the Hard Thermal Loop approximation. This combination of known effective descriptions offers unprecedented access to small but nonzero temperatures, both in and out of beta equilibrium.

preprint2016arXiv

Initial conditions for hydrodynamics from weakly coupled pre-equilibrium evolution

We use effective kinetic theory, accurate at weak coupling, to simulate the pre-equilibrium evolution of transverse energy and flow perturbations in heavy-ion collisions. We provide a Green function which propagates the initial perturbations to the energy-momentum tensor at a time when hydrodynamics becomes applicable. With this map, the complete pre-thermal evolution from saturated nuclei to hydrodynamics can be modelled in a perturbatively controlled way.

preprint2016arXiv

On high-order perturbative calculations at finite density

We discuss the prospects of performing high-order perturbative calculations in systems characterized by a vanishing temperature but finite density. In particular, we show that the determination of generic Feynman integrals containing fermionic chemical potentials can be reduced to the evaluation of three-dimensional phase space integrals over vacuum on-shell amplitudes - a result reminiscent of a previously proposed "naive real-time formalism" for vacuum diagrams. Applications of these rules are discussed in the context of the thermodynamics of cold and dense QCD, where it is argued that they facilitate an extension of the Equation of State of cold quark matter to higher perturbative orders.

preprint2015arXiv

Isotropization and hydrodynamization in weakly coupled heavy-ion collisions

We numerically solve 2+1D effective kinetic theory of weak coupling QCD under longitudinal expansion relevant for early stages of heavy-ion collisions. We find agreement with viscous hydrodynamics and classical Yang-Mills simulations in the regimes where they are applicable. By choosing initial conditions that are motivated by color-glass-condensate framework we find that for Q=2GeV and $α_s$=0.3 the system is approximately described by viscous hydrodynamics well before $τ\lesssim 1.0$ fm/c.

preprint2015arXiv

Neutron star structure from QCD

In this review article, we argue that our current understanding of the thermodynamic properties of cold QCD matter, originating from first principles calculations at high and low densities, can be used to efficiently constrain the macroscopic properties of neutron stars. In particular, we demonstrate that combining state-of-the-art results from Chiral Effective Theory and perturbative QCD with the current bounds on neutron star masses, the Equation of State of neutron star matter can be obtained to an accuracy better than 30% at all densities.

preprint2015arXiv

Weak and strong coupling equilibration in nonabelian gauge theories

We present a direct comparison studying equilibration through kinetic theory at weak coupling and through holography at strong coupling in the same set-up. The set-up starts with a homogeneous thermal state, which then smoothly transitions through an out-of-equilibrium phase to an expanding system undergoing boost-invariant flow. This first apples-to-apples comparison of equilibration provides a benchmark for similar equilibration processes in heavy-ion collisions, where the equilibration mechanism is still under debate. We find that results at weak and strong coupling can be smoothly connected by simple, empirical power-laws for the viscosity, equilibration time and entropy production of the system.

preprint2014arXiv

Approach to equilibrium in weakly coupled nonabelian plasmas

We follow the time evolution of nonabelian gauge bosons from far-from-equilibrium initial conditions to thermal equilibrium by numerically solving an effective kinetic equation that becomes accurate in the weak coupling limit. We consider initial conditions that are either highly overoccupied or underoccupied. We find that overoccupied systems thermalize through a turbulent cascade reaching equilibrium in multiples of a thermalization time $t\approx 72./ (1-0.12\log λ)/λ^2 T$, whereas underoccupied systems undergo a "bottom-up" thermalization in a time $t\approx (34. +21. \ln(Q/T))/ (1-0.037\log λ)(Q/T)^{1/2}/λ^2 T$, where $Q$ is the characteristic momentum scale of the initial condition. We apply this result to model initial stages of heavy-ion collisions and find rapid thermalization roughly in a time $Qt \lesssim 10$ or $t\lesssim 1$ fm/c.

preprint2014arXiv

Constraining neutron star matter with QCD

In recent years, there have been several successful attempts to constrain the equation of state of neutron star matter using input from low-energy nuclear physics and observational data. We demonstrate that significant further restrictions can be placed by additionally requiring the pressure to approach that of deconfined quark matter at high densities. Remarkably, the new constraints turn out to be highly insensitive to the amount --- or even presence --- of quark matter inside the stars.

preprint2014arXiv

Interacting quark matter equation of state for compact stars

Lattice QCD studies of the thermodynamics of the hot quark-gluon plasma (QGP) demonstrate the importance of accounting for the interactions of quarks and gluons, if one wants to investigate the phase structure of strongly interacting matter. Motivated by this observation and using state-of-the-art results from perturbative QCD, we construct a simple effective equation of state for cold quark matter that consistently incorporates the effects of interactions and furthermore includes a built-in estimate of the inherent systematic uncertainties. This goes beyond the MIT bag model description in a crucial way, yet leads to an equation of state that is equally straightforward to use.

preprint2014arXiv

Renormalization of Null Wilson Lines in EQCD

Radiation and energy loss of a light, high-energy parton in a perturbative Quark-Gluon Plasma is controlled by transverse momentum exchange. The troublesome infrared contributions to transverse momentum exchange can be computed on the lattice using dimensional reduction to EQCD. However a novel extended operator, the Null Wilson Line of EQCD, is involved. We compute the renormalization properties of this object's lattice implementation to next-to-leading order, which should facilitate its efficient calculation on the lattice.

preprint2014arXiv

UV Cascade in Classical Yang-Mills via Kinetic Theory

We show that classical Yang-Mills theory with statistically homogeneous and isotropic initial conditions has a kinetic description and approaches a scaling solution at late times. We find the scaling solution by explicitly solving the Boltzmann equations, including all dominant processes (elastic and number-changing). Above a scale $p_{max} \propto t^{1/7}$ the occupancy falls exponentially in $p$. For asymptotically late times and sufficiently small momenta the occupancy scales as $f(p)\propto 1/p$, but this behavior sets in only at very late time scales. We find quantitative agreement of our results with lattice simulations, for times and momenta within the range of validity of kinetic theory.

preprint2013arXiv

Next-to-leading order thermal photon production in a weakly coupled quark-gluon plasma

We compute the next-to-leading order O(g) correction to the thermal photon production rate in a QCD plasma. The NLO contributions can be expressed in terms of gauge invariant condensates on the light cone, which are amenable to novel sum rules and Euclidean techniques. We expect these technologies to be generalizable to other NLO calculations. For the phenomenologically interesting value of alpha_s=0.3, the NLO correction represents a 20% increase and has a functional form similar to the LO result.

preprint2013arXiv

Thermalization in collisions of large nuclei at high energies

Hydrodynamical analysis of experimental data of ultrarelativistic heavy ion collisions seems to indicate that the hot QCD matter created in the collisions thermalizes very quickly. Theoretically, we have no idea why this should be true. In this proceeding, I will describe how the thermalization takes place in the most theoretically clean limit -- that of large nuclei at asymptotically high energy per nucleon, where the system is described by weak-coupling QCD. In this limit, plasma instabilities dominate the dynamics from immediately after the collision until well after the plasma becomes nearly in equilibrium at time t α^(-5/2)Q^(-1).

preprint2012arXiv

Numerical properties of staggered quarks with a taste-dependent mass term

The numerical properties of staggered Dirac operators with a taste-dependent mass term proposed by Adams [1,2] and by Hoelbling [3] are compared with those of ordinary staggered and Wilson Dirac operators. In the free limit and on (quenched) interacting configurations, we consider their topological properties, their spectrum, and the resulting pion mass. Although we also consider the spectral structure, topological properties, locality, and computational cost of an overlap operator with a staggered kernel, we call attention to the possibility of using the Adams and Hoelbling operators without the overlap construction. In particular, the Hoelbling operator could be used to simulate two degenerate flavors without additive mass renormalization, and thus without fine-tuning in the chiral limit.

preprint2012arXiv

Two-color QCD via dimensional reduction

We study the thermodynamics of two-color QCD at high temperature and/or density using a dimensionally reduced superrenormalizable effective theory, formulated in terms of a coarse grained Wilson line. In the absence of quarks, the theory is required to respect the Z(2) center symmetry, while the effects of quarks of arbitrary masses and chemical potentials are introduced via soft Z(2) breaking operators. Perturbative matching of the effective theory parameters to the full theory is carried out explicitly, and it is argued how the new theory can be used to explore the phase diagram of two-color QCD.

preprint2012arXiv

UV Cascade in Classical Yang-Mills Theory

We study the real-time behavior of classical Yang-Mills theory under initial conditions with nonperturbatively large, infrared field amplitudes. Our lattice study confirms the cascade of energy towards higher momenta and lower occupancy, which occurs via a scaling solution $f[p,t_1] = (t_0/t_1)^{4/7}\, f[p (t_0/t_1)^{1/7},t_0]$. Above a characteristic scale p_{max}, f falls exponentially; below p_{max}, $f[p] \propto p^{-4/3}$. We find no evidence for different infrared exponents or for infrared occupancies in excess of those described by this scaling solution. We also investigate what the fate of large occupancies would be, both in the electric and the magnetic sector.

preprint2011arXiv

Bjorken Flow, Plasma Instabilities, and Thermalization

At asymptotically high energies, thermalization in heavy ion collisions can be described via weak-coupling QCD. We present a complete treatment of how thermalization proceeds, at the parametric weak-coupling level. We show that plasma instabilities dominate the dynamics, from immediately after the collision until well after the plasma becomes nearly in equilibrium. Initially they drive the system close to isotropy, but Bjorken expansion and increasing diluteness makes the system again become more anisotropic. At time τ~ α^(-12/5) Q^(-1) the dynamics become dominated by a nearly-thermal bath; and at time τ~ α^(-5/2) Q^(-1)$ the bath comes to dominate the energy density, completing thermalization. After this time there is a nearly isotropic and thermal Quark-Gluon Plasma.

preprint2011arXiv

Numerical properties of staggered overlap fermions

We report the results of a numerical study of staggered overlap fermions, following the construction of Adams which reduces the number of tastes from 4 to 2 without fine-tuning. We study the sensitivity of the operator to the topology of the gauge field, its locality and its robustness to fluctuations of the gauge field. We make a first estimate of the computing cost of a quark propagator calculation, and compare with Neuberger's overlap.

preprint2011arXiv

Thermalization in Weakly Coupled Nonabelian Plasmas

We investigate how relativistic, nonabelian plasmas approach equilibrium in a general context. Our treatment is entirely parametric and for small Yang-Mills coupling $α$. First we study isotropic systems with an initially nonequilibrium momentum distribution. We consider both the case of initially very high occupancy and initially very low occupancy. Then we consider systems which are anisotropic. We consider both weak anisotropy and large anisotropy, and allow the occupancy to be parametrically large or small. Writing the typical momentum of an initial excitation as Q and the final temperature as T, full equilibration occurs in a time t ~ α^{-2}/T for T > Q, and t ~ α^{-2} Q^{1/2} T^{-3/2} for T < Q, unless the initial system is sufficiently anisotropic and T > α^{2/3} Q, in which case equilibration occurs somewhat faster, t ~ α^{-13/7} Q^{5/7} T^{-12/7} (or α^{-2}/T if that is longer).

preprint2010arXiv

Cold Quark Matter

We perform an O(alpha_s^2) perturbative calculation of the equation of state of cold but dense QCD matter with two massless and one massive quark flavor, finding that perturbation theory converges reasonably well for quark chemical potentials above 1 GeV. Using a running coupling constant and strange quark mass, and allowing for further non-perturbative effects, our results point to a narrow range where absolutely stable strange quark matter may exist. Absent stable strange quark matter, our findings suggest that quark matter in compact star cores becomes confined to hadrons only slightly above the density of atomic nuclei. Finally, we show that equations of state including quark matter lead to hybrid star masses up to M~2M_solar, in agreement with current observations. For strange stars, we find maximal masses of M~2.75M_solar and conclude that confirmed observations of compact stars with M>2M_solar would strongly favor the existence of stable strange quark matter.

preprint2010arXiv

Looking inside neutron stars: Microscopic calculations confront observations

While QCD appears not to be accurately solvable in the regime of interest for neutron star physics, microscopic calculations are feasible at both low and very high densities. In this work, we propose using the most realistic calculations in these two regimes of nuclear physics and perturbative QCD, and construct equations of state by matching the results requiring thermodynamic consistency. We find that the resulting equations of state --- in contrast to several hadronic ones --- are able to reproduce current observational data on neutron stars without any fine tuning, and allow stable hybrid stars with masses up to 2.1M_{sun}. Using recent observations of star radii, we perform a maximum likelihood analysis to further constrain the equation of state, and in addition show that the effects of rotation on radii and masses should not be neglected in future precision studies.

preprint2010arXiv

The phase diagram of Yang-Mills theory with a compact extra dimension

We present a non-perturbative study of the phase diagram of SU(2) Yang-Mills theory in a five-dimensional spacetime with a compact extra dimension. The non-renormalizable theory is regularized on an anisotropic lattice and investigated through numerical simulations in a regime characterized by a hierarchy between the scale of low-energy physics, the inverse compactification radius, and the cutoff scale. We map out the structure of the phase diagram and the pattern of lines corresponding to fixed values of the ratio between the mass of the fifth component of the gauge field and the non-perturbative mass gap of the four-dimensional modes. We discuss different limits of the model, and comment on the implications of our findings.