The recent detection of gravitational waves from merging neutron star events has opened a new window on the many unknown aspects of their internal dynamics. A key role in this context is played by the transition from baryon to quark matter described in the neutron star equation of state (EoS). In particular, the binary pulsar observation of heavy neutron stars requires appropriately stiff dense matter in order to counter gravitational collapse, at variance with the predictions of many phenomenological quark models. On the other side, the LIGO observations favor a softer EoS therefore providing a lower bound to the equation stiffness. We introduce a quantum chromodynamics (QCD) description of the neutron star's high baryon density regime where the pressure and energy density distributions are directly obtained from the matrix elements of the QCD energy momentum tensor. Recent ab initio calculations allow us to evaluate the energy-momentum tensor in a model independent way including both quark and gluon degrees of freedom. Our approach is a first effort to replace quark models and effective gluon interactions with a first principles, fully QCD-based description. Most importantly, the QCD energy momentum tensor matrix elements are connected to the Mellin moments of the generalized parton distributions which can be measured in deeply virtual exclusive scattering experiments. As a consequence, we establish a connection between observables from high energy experiments and from the analysis of gravitational wave events. Both can be used to mutually constrain the respective sets of data. In particular, the emerging QCD-based picture is consistent with the GW170817 neutron star merger event once we allow a first-order phase transition from a low-density nuclear matter EoS to the newly-constructed high-density quark-gluon one.

In this conference-proceedings contribution, we review recent advances in placing model-independent constraints on the properties of cold and dense QCD matter inside neutron stars. In addition to introducing new bounds for the Equation of State, we explain how these results may be used to make robust statements about the physical phase of strongly interacting matter in the centers of neutron stars of different masses. Our findings indicate that the existence of quark-matter cores inside massive neutron stars appears to be a very common feature of the allowed Equations of State, and should not be considered an exotic or unlikely scenario.

We combine results from astrophysical observations and theoretical ab-initio calculations to investigate the presence of quark matter (QM) inside neutron stars (NSs). We find a clear qualitative change in the material properties of NS matter at energy densities comparable to those where quark-gluon plasma is created in heavy-ion collisions. In the low-density phase, the system has characteristics closely resembling those of hadronic matter, while the high-density phase can be clearly identified with nearly conformal QM. We show that QM is never present inside 1.44-M_sol NSs, but resides inside maximally massive NSs barring very specific and extreme conditions. Finally, for the heaviest observed NSs of approximately two solar masses, the fate of QM is found to strongly depend on the behavior of the speed of sound in the matter: if this quantity does not strongly violate the conformal limit c_s^2 < 1/3, these stars host sizable quark cores of at least five kilometers in radius.

At high baryon chemical potential mu_B, the equation of state of QCD allows a weak-coupling expansion in the QCD coupling alpha_s. The result is currently known up to and including the full next-to-next-to-leading order (NNLO) alpha_s^2. Starting at this order, the computations are complicated by the modification of particle propagation in a dense medium, which necessitates non-perturbative treatment of the scale Sqrt(alpha_s)*mu_B. In this work, we apply a Hard-Thermal-Loop scheme for capturing the contributions of this scale to the weak-coupling expansion, and use it to determine the leading-logarithm contribution to NNNLO: alpha_s^3 (ln alpha_s)^2. This result is the first improvement to the equation of state of massless cold quark matter in 40 years. The new term is negligibly small, and thus significantly increases our confidence in the applicability of the weak-coupling expansion.

The LIGO/Virgo detection of gravitational waves originating from a neutron-star merger, GW170817, has recently provided new stringent limits on the tidal deformabilities of the stars involved in the collision. Combining this measurement with the existence of two-solar-mass stars, we generate a generic family of neutron-star-matter Equations of State (EoSs) that interpolate between state-of-the-art theoretical results at low and high baryon density. Comparing the results to ones obtained without the tidal-deformability constraint, we witness a dramatic reduction in the family of allowed EoSs. Based on our analysis, we conclude that the maximal radius of a 1.4-solar-mass neutron star is 13.6 km, and that smallest allowed tidal deformability of a similar-mass star is $\Lambda(1.4 M_\odot) = 120$.

We numerically investigate global properties of rotating neutron stars using the allowed band of QCD equations of state derived by Kurkela et al. This band is constrained by chiral effective theory at low densities and perturbative QCD at high densities, and is thus, in essence, a controlled constraint from first-principles physics. Previously, this band of equations of state was used to investigate non-rotating neutron stars only; in this work, we extend these results to any rotation frequency below the mass-shedding limit. We investigate mass--radius curves, allowed mass--frequency regions, radius--frequency curves for a typical 1.4-solar-mass star, and the values of the moment of inertia of the double pulsar PSR J0737-3039A, a pulsar whose moment of inertia may be constrained observationally in a few years. We present limits on observational data coming from these constraints, and identify values of observationally-relevant parameters that would further constrain the allowed region for the QCD equation of state. We also discuss how much this region would be constrained by a measurement of the moment of inertial of the double pulsar PSR J0737-3039A.

The equation of state (EoS) of quantum chromodynamics (QCD) at zero temperature can be calculated in two different perturbative regimes: for small values of the baryon chemical potential $\mu$, one may use chiral perturbation theory (ChEFT); and for large values of $\mu$, one may use perturbative QCD (pQCD). There is, however, a gap for $\mu \in (0.97\text{ GeV},\, 2.6\text{ GeV})$, where these theories becomes non-perturbative, and where there is currently no known microscopic description of QCD matter. Unfortunately, this interval obscures the values of $\mu$ found within the cores of neutron stars (NSs). In this thesis, we argue that thermodynamic matching of the ChEFT and pQCD EoSs is a legitimate way to obtain quantitative constraints on the non-pertubative QCD EoS. Moreover, we argue that this method is effective, verifiable, and systematically improvable. First, we carry out a simplified matching procedure in QCD-like theories that can be simulated on the lattice without a sign problem. Our calculated pressure band serves as a prediction for lattice-QCD practitioners and will allow them to verify or refute the simplified procedure. Second, we apply the state-of-the-art matched EoS of Kurkela et al. (2014) to rotating NSs. This allows us to obtain bounds on observable NS properties, as well as point towards future observations that would more tightly constrain the current state-of-the-art EoS band. Finally, as evidence of the ability to improve the procedure, we carry out calculations in pQCD to improve the zero-temperature pressure. We calculate the full $\mathcal{O}(g^{6} \ln^{2} g)$ contribution to the pQCD pressure for $n_{f}$ massless quarks, as well as a significant portion of the $\mathcal{O}(g^{6} \ln g)$ piece and even some of the $\mathcal{O}(g^{6})$ piece.