We present a new equation of state for QCD in which the temperature $T$ and the three chemical potentials for baryon number $\mu_B$, electric charge $\mu_Q$ and strangeness $\mu_S$ can be varied independently. This result is based on a generalization of the $T'$-expansion scheme, thanks to which the diagonal $\mu_B$ extrapolation was pushed up to a baryo-chemical potential $\mu_B/T \sim 3.5$ for the first time. This considerably extended the coverage of the Taylor expansion, limited to \mu_B/T < 2.5-3. As a consequence, we are able to offer a substantially larger coverage of the four-dimensional QCD phase diagram as well, compared to previously available Taylor expansion results. Our results are based on new continuum estimated lattice results on the full set of second and fourth order fluctuations.

We have performed high-statistics lattice simulations using 4HEX improved staggered fermions on $16^3 \times 8$ lattices. We calculated the Taylor expansion coefficients of the pressure with respect to the baryochemical potential to the tenth order at zero, and fourth order at purely imaginary chemical potentials. We used this data to construct rational function approximations of the free energy. We use a rational ansatz that explicitly satisfies the charge conjugation symmetry and the Roberge-Weiss periodicity, which are exact properties of the QCD free energy. We use this ansatz to estimate the position of Lee-Yang zeros in the complex chemical potential plane. The temperature dependence of the imaginary part of the Lee-Yang zeros is then fitted with ansätze motivated by the universal behavior of the free energy near a 3D Ising critical point. In principle, this allows one to estimate the temperature of the critical endpoint. We consider several sources of systematic errors. On this single lattice spacing we find that with $84\%$ probability, the chiral critical endpoint is either below $103$~MeV temperature or it does not exist. We also identify some caveats of the method, which do not disappear even with the extremely high statistics of this present study. We discuss to what extent these can be eliminated by future high statistics lattice analyses.

We present the full results of the Wuppertal-Budapest lattice QCD collaboration on flavor diagonal and non-diagonal quark number susceptibilities with 2+1 staggered quark flavors, in a temperature range between 125 and 400 MeV. The light and strange quark masses are set to their physical values. Lattices with Nt=6, 8, 10, 12, 16 are used. We perform a continuum extrapolation of all observables under study. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized. All results are compared to the Hadron Resonance Gas model predictions: good agreement is found in the temperature region below the transition.

High-energy cosmic rays are observed indirectly by detecting the extensive air showers initiated in Earth's atmosphere. The interpretation of these observations relies on accurate models of air shower physics, which is a challenge and an opportunity to test QCD under extreme conditions. Air showers are hadronic cascades, which eventually decay into muons. The muon number is a key observable to infer the mass composition of cosmic rays. Air shower simulations with state-of-the-art QCD models show a significant muon deficit with respect to measurements; this is called the Muon Puzzle. The origin of this discrepancy has been traced to the composition of secondary particles in hadronic interactions. The muon discrepancy starts at the TeV scale, which suggests that this change in hadron composition is observable at the Large Hadron Collider. An effect that can potentially explain the puzzle has been discovered at the LHC, but needs to be confirmed for forward produced hadrons with LHCb, and with future data on oxygen beams.

We perform a detailed, fully-correlated study of the chiral behavior of the pion mass and decay constant, based on 2+1 flavor lattice QCD simulations. These calculations are implemented using tree-level, O(a)-improved Wilson fermions, at four values of the lattice spacing down to 0.054 fm and all the way down to below the physical value of the pion mass. They allow a sharp comparison with the predictions of SU(2) chiral perturbation theory (\chi PT) and a determination of some of its low energy constants. In particular, we systematically explore the range of applicability of NLO SU(2) \chi PT in two different expansions: the first in quark mass (x-expansion), and the second in pion mass (\xi-expansion). We find that these expansions begin showing signs of failure around M_\pi=300 MeV for the typical percent-level precision of our N_f=2+1 lattice results. We further determine the LO low energy constants (LECs), F=88.0 \pm 1.3\pm 0.3 and B^\msbar(2 GeV)=2.58 \pm 0.07 \pm 0.02 GeV, and the related quark condensate, \Sigma^\msbar(2 GeV)=(271\pm 4\pm 1 MeV)^3, as well as the NLO ones, l_3=2.5 \pm 0.5 \pm 0.4 and l_4=3.8 \pm 0.4 \pm 0.2, with fully controlled uncertainties. We also explore the NNLO expansions and the values of NNLO LECs. In addition, we show that the lattice results favor the presence of chiral logarithms. We further demonstrate how the absence of lattice results with pion masses below 200 MeV can lead to misleading results and conclusions. Our calculations allow a fully controlled, ab initio determination of the pion decay constant with a total 1% error, which is in excellent agreement with experiment.

We calculate second- and fourth-order cumulants of conserved charges in a temperature range stretching from the QCD transition region towards the realm of (resummed) perturbation theory. We perform lattice simulations with staggered quarks; the continuum extrapolation is based on $N_t=10\dots24$ in the crossover-region and $N_t=8\dots16$ at higher temperatures. We find that the Hadron Resonance Gas model predictions describe the lattice data rather well in the confined phase. At high temperatures (above $\sim$250 MeV) we find agreement with the three-loop Hard Thermal Loop results.

We highlight the progress, current status, and open challenges of QCD-driven physics, in theory and in experiment. We discuss how the strong interaction is intimately connected to a broad sweep of physical problems, in settings ranging from astrophysics and cosmology to strongly-coupled, complex systems in particle and condensed-matter physics, as well as to searches for physics beyond the Standard Model. We also discuss how success in describing the strong interaction impacts other fields, and, in turn, how such subjects can impact studies of the strong interaction. In the course of the work we offer a perspective on the many research streams which flow into and out of QCD, as well as a vision for future developments.

Phase transitions in a non-perturbative regime can be studied by ab initio Lattice Field Theory methods. The status and future research directions for LFT investigations of Quantum Chromo-Dynamics under extreme conditions are reviewed, including properties of hadrons and of the hypothesized QCD axion as inferred from QCD topology in different phases. We discuss phase transitions in strong interactions in an extended parameter space, and the possibility of model building for Dark Matter and Electro-Weak Symmetry Breaking. Methodological challenges are addressed as well, including new developments in Artificial Intelligence geared towards the identification of different phases and transitions.

We provide the most accurate results for the QCD transition line so far. We optimize the definition of the crossover temperature $T_c$, allowing for its very precise determination, and extrapolate from imaginary chemical potential up to real $\mu_B \approx 300$ MeV. The definition of $T_c$ adopted in this work is based on the observation that the chiral susceptibility as a function of the condensate is an almost universal curve at zero and imaganiary $\mu_B$. We obtain the parameters $\kappa_2=0.0153(18)$ and $\kappa_4=0.00032(67)$ as a continuum extrapolation based on $N_t=10,12$ and $16$ lattices with physical quark masses. We also extrapolate the peak value of the chiral susceptibility and the width of the chiral transition along the crossover line. In fact, both of these are consistent with a constant function of $\mu_B$. We see no sign of criticality in the explored range.

The QCD cross-over line in the temperature ($T$) -- baryo-chemical potential ($\mu_B$) plane has been computed by several lattice groups by calculating the chiral order parameter and its susceptibility at finite values of $\mu_B$. In this work we focus on the deconfinement aspect of the transition between hadronic and Quark Gluon Plasma (QGP) phases. We define the deconfinement temperature as the peak position of the static quark entropy ($S_Q(T,\mu_B)$) in $T$, which is based on the renormalized Polyakov loop. We extrapolate $S_Q(T,\mu_B)$ based on high statistics finite temperature ensembles on a $16^3\times 8$ lattice to finite density by means of a Taylor expansion to eighth order in $\mu_B$ (NNNLO) along the strangeness neutral line. For the simulations the 4HEX staggered action was used with 2+1 flavors at physical quark masses. In this setup the phase diagram can be drawn up to unprecedentedly high chemical potentials. Our results for the deconfinement temperature are in rough agreement with phenomenological estimates of the freeze-out curve in relativistic heavy ion collisions. In addition, we study the width of the deconfinement crossover. We show that up to $\mu_B \approx 400$ MeV, the deconfinement transition gets broader at higher densities, disfavoring the existence of a deconfinement critical endpoint in this range. Finally, we examine the transition line without the strangeness neutrality condition and observe a hint for the narrowing of the crossover towards large $\mu_B$.

Approaches to finite baryon density lattice QCD usually suffer from uncontrolled systematic uncertainties in addition to the well-known sign problem. We test a method - sign reweighting - that works directly at finite chemical potential and is yet free from any such uncontrolled systematics: with this approach the only problem is the sign problem itself. In practice the approach involves the generation of configurations with the positive fermionic weights given by the absolute value of the real part of the quark determinant, and a reweighting by a sign. There are only two sectors, +1 and -1 and as long as the average $\left\langle \pm \right\rangle \neq 0$ (with respect to the positive weight) this discrete reweighting has no overlap problem - unlike reweighting from $\mu=0$ - and the results are reliable. We also present results based on this algorithm on the phase diagram of lattice QCD with two different actions: as a first test, we apply the method to calculate the position of the critical endpoint with unimproved staggered fermions at $N_\tau=4$; as a second application, we study the phase diagram with 2stout improved staggered fermions at $N_\tau=6$. This second one is already a reasonably fine lattice - relevant for phenomenology. We demonstrate that the method penetrates the region of the phase diagram where the Taylor and imaginary chemical potential methods lose predictive power.

Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science

Axions are one of the most attractive dark matter candidates. The evolution of their number density in the early universe can be determined by calculating the topological susceptibility $\chi(T)$ of QCD as a function of the temperature. Lattice QCD provides an ab initio technique to carry out such a calculation. A full result needs two ingredients: physical quark masses and a controlled continuum extrapolation from non-vanishing to zero lattice spacings. We determine $\chi(T)$ in the quenched framework (infinitely large quark masses) and extrapolate its values to the continuum limit. The results are compared with the prediction of the dilute instanton gas approximation (DIGA). A nice agreement is found for the temperature dependence, whereas the overall normalization of the DIGA result still differs from the non-perturbative continuum extrapolated lattice results by a factor of order ten. We discuss the consequences of our findings for the prediction of the amount of axion dark matter.

The QCD transition is studied on lattices up to $N_t=16$. The chiral condensate is presented as a function of the temperature, and the corresponding transition temperature is extracted. The equation of state is determined on lattices with $N_t=6,8,10$ and at some temperature values with $N_t=12$. The pressure and the trace anomaly are presented as functions of the temperature in the range 100 ...1000 MeV . Using the same configurations we determine the continuum extrapolated phase diagram of QCD on the $\mu-T$ plane for small to moderate chemical potentials. Two transition lines are defined with two quantities, the chiral condensate and the strange quark number susceptibility.

In a previous letter (arXiv:1306.2287) we determined the isospin mass splittings of the baryon octet from a lattice calculation based on quenched QED and $N_f{=}2{+}1$ QCD simulations with 5 lattice spacings down to $0.054~\mathrm{fm}$, lattice sizes up to $6~\mathrm{fm}$ and average up-down quark masses all the way down to their physical value. Using the same data we determine here the corrections to Dashen's theorem and the individual up and down quark masses. For the parameter which quantifies violations to Dashens's theorem, we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4)~\mathrm{MeV}$ and $m_d=4.67(6)(5)(4)~\mathrm{MeV}$ in the $\bar{\mathrm{MS}}$ scheme at $2~\mathrm{GeV}$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than $24$ standard deviations.

We propose a theoretical analysis of the experimental data by Ichikawa et al. (Phys. Rev. Lett. 112, 071101 (2014)) on a spatial distribution of ultracold neutrons in the gravitational field of the Earth above a mirror, projected onto a pixelated detector by scattering by a cylindrical mirror. We calculate a theoretical spatial distribution of a probability density of quantum gravitational states of ultracold neutrons and analyse a spatial distribution of quantum gravitational states in term of the Wigner function. We cannot confirm that the experimental data by Ichikawa et al. (Phys. Rev. Lett. 112, 071101 (2014)) correspond to a spatial distribution of quantum gravitational states of ultracold neutrons.

In this work we study the non-equilibrium Markov state evolution for a spatial population model on the space of locally finite configurations $\Gamma^2 = \Gamma^+ \times \Gamma^-$ over $\mathbb{R}^d$ where particles are marked by spins $\pm$. Particles of type '+' reproduce themselves independently of each other and, moreover, die due to competition either among particles of the same type or particles of different type. Particles of type '-' evolve according to a non-equilibrium Glauber-type dynamics with activity $z$ and potential $\psi$. Let $L^S$ be the Markov operator for '+' -particles and $L^E$ the Markov operator for '-' -particles. The non-equilibrium state evolution $(\mu_t^{\varepsilon})_{t \geq 0}$ is obtained from the Fokker-Planck equation with Markov operator $L^S + \frac{1}{\varepsilon}L^E$, $\varepsilon > 0$, which itself is studied in terms of correlation function evolution on a suitable chosen scale of Banach spaces. We prove that in the limiting regime $\varepsilon \to 0$ the state evolution $\mu_t^{\varepsilon}$ converges weakly to some state evolution $\overline{\mu}_t$ associated to the Fokker-Planck equation with (heuristic) Markov operator obtained from $L^S$ by averaging the interactions of the system with the environment with respect to the unique invariant Gibbs measure of the environment.

We present the new results of the Wuppertal-Budapest lattice QCD collaboration on flavor diagonal and non-diagonal quark number susceptibilities with 2+1 staggered quark flavors, in a temperature regime between 120 and 400 MeV. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized; the light and strange quark masses are set to their physical values. Lattices with $N_t=6, 8, 10, 12$ are used. We perform a continuum extrapolation of those observables for which the scaling regime is reached, and discretization errors are under control.

The Neutrino Mass Ordering (NMO) remains one of the outstanding questions in the field of neutrino physics. One strategy to measure the NMO is to observe matter effects in the oscillation pattern of atmospheric neutrinos above $\sim 1\,\mathrm{GeV}$, as proposed for several next-generation neutrino experiments. Moreover, the existing IceCube DeepCore detector can already explore this type of measurement. We present rthe development and application of two independent analyses to search for the signature of the NMO with three years of DeepCore data. These analyses include a full treatment of systematic uncertainties and a statistically-rigorous method to determine the significance for the NMO from a fit to the data. Both analyses show that the dataset is fully compatible with both mass orderings. For the more sensitive analysis, we observe a preference for Normal Ordering with a $p$-value of $p_\mathrm{IO} = 15.3\%$ and $\mathrm{CL}_\mathrm{s}=53.3\%$ for the Inverted Ordering hypothesis, while the experimental results from both analyses are consistent within their uncertainties. Since the result is independent of the value of $\delta_\mathrm{CP}$ and obtained from energies $E_\nu \gtrsim 5\,\mathrm{GeV}$, it is complementary to recent results from long-baseline experiments. These analyses set the groundwork for the future of this measurement with more capable detectors, such as the IceCube Upgrade and the proposed PINGU detector.

We calculate a resummed equation of state with lattice QCD simulations at imaginary chemical potentials. This work presents a generalization of the scheme introduced in 2102.06660 to the case of non-zero $\mu_S$, focusing on the line of strangeness neutrality. We present results up to $\mu_B/T \leq 3.5$ on the strangeness neutral line $\left\langle S \right\rangle = 0$ in the temperature range $130 \rm{MeV} \leq T \leq 280 \rm{MeV}$. We also extrapolate the finite baryon density equation of state to small non-zero values of the strangeness-to-baryon ratio $R=\left\langle S \right\rangle / \left\langle B \right\rangle$. We perform a continuum extrapolation using lattice simulations of the 4stout-improved staggered action with 8, 10, 12 and 16 timeslices.