Pulsar Nulling and Vacuum Radio Emission from Axion Clouds

Nonrelativistic axions can be efficiently produced in the polar caps of pulsars, resulting in the formation of a dense cloud of gravitationally bound axions. Here, we investigate the interplay between such an axion cloud and the electrodynamics in the pulsar magnetosphere, focusing specifically on the dynamics in the polar caps, where the impact of the axion cloud is expected to be most pronounced. For sufficiently light axions m_{a}≲10^{-7}  eV, we show that the axion cloud can occasionally screen the local electric field responsible for particle acceleration and pair production, inducing a periodic nulling of the pulsar's intrinsic radio emission. At larger axion masses, the small-scale fluctuations in the axion field tend to suppress the backreaction of the axion on the electrodynamics; however, we point out that the incoherent oscillations of the axion in short-lived regions of vacuum near the neutron star surface can produce a narrow radio line, which provides a complementary source of radio emission to the plasma-resonant emission processes identified in previous work. While this Letter focuses on the leading order correction to pair production in the magnetosphere, we speculate that there can exist dramatic deviations in the electrodynamics of these systems when the axion backreaction becomes nonlinear.

Astrophysical Equation-of-State Constraints on the Color-Superconducting Gap

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J0348+0432, PSR J1624-2230, PSR J0740+6620, GW170817). At baryon number chemical potential μ=2.6  GeV we find a 95% upper limit on the color-flavor locked pairing gap Δ of 457 MeV using overly conservative assumptions and 216 MeV with more reasonable assumptions. This constraint may be strengthened by future astrophysical measurements as well as by future advances in high-density QCD calculations.

Momentum Shell in Quarkyonic Matter from Explicit Duality: A Dual Model for Cold, Dense QCD

We present a model of cold QCD matter that bridges nuclear and quark matter through the duality relation between quarks and baryons. The baryon number and energy densities are expressed as functionals of either the baryon momentum distribution, f_{B}, or the quark distribution, f_{Q}, which are subject to the constraints on fermions, 0≤f_{B,Q}≤1. The theory is ideal in the sense that the confinement of quarks into baryons is reflected in the duality relation between f_{Q} and f_{B}, while other possible interactions among quarks and baryons are all neglected. The variational problem with the duality constraints is formulated and we explicitly construct analytic solutions, finding two distinct regimes: a nuclear matter regime at low density and a quarkyonic regime at high density. In the quarkyonic regime, baryons underoccupy states at low momenta but form a momentum shell with f_{B}=1 on top of a quark Fermi sea. Such a theory describes a rapid transition from a soft nuclear equation of state to a stiff quarkyonic equation of state. At this transition, there is a rapid increase in the pressure.

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M⊙ neutron star to be [Formula: see text] km.

Strongly interacting matter exhibits deconfined behavior in massive neutron stars

Neutron-star cores contain matter at the highest densities in our Universe. This highly compressed matter may undergo a phase transition where nuclear matter melts into deconfined quark matter, liberating its constituent quarks and gluons. Quark matter exhibits an approximate conformal symmetry, predicting a specific form for its equation of state (EoS), but it is currently unknown whether the transition takes place inside at least some physical neutron stars. Here, we quantify this likelihood by combining information from astrophysical observations and theoretical calculations. Using Bayesian inference, we demonstrate that in the cores of maximally massive stars, the EoS is consistent with quark matter. We do this by establishing approximate conformal symmetry restoration with high credence at the highest densities probed and demonstrating that the number of active degrees of freedom is consistent with deconfined matter. The remaining likelihood is observed to correspond to EoSs exhibiting phase-transition-like behavior, treated as arbitrarily rapid crossovers in our framework.

Degeneracy in the Inference of Phase Transitions in the Neutron Star Equation of State from Gravitational Wave Data

Gravitational wave (GW) detections of binary neutron star inspirals will be crucial for constraining the dense matter equation of state (EOS). We demonstrate a new degeneracy in the mapping from tidal deformability data to the EOS, which occurs for models with strong phase transitions. We find that there exists a new family of EOS with phase transitions that set in at different densities and that predict neutron star radii that differ by up to ∼500  m but that produce nearly identical tidal deformabilities for all neutron star masses. Next-generation GW detectors and advances in nuclear theory may be needed to resolve this degeneracy.

Plausible presence of new state in neutron stars with masses above 0.98MTOV

We investigate the neutron star (NS) equation of state (EOS) by incorporating multi-messenger data of GW170817, PSR J0030 + 0451, PSR J0740 + 6620, and state-of-the-art theoretical progresses, including the information from chiral effective field theory (χEFT) and perturbative quantum chromodynamics (pQCD) calculation. Taking advantage of the various structures sampling by a single-layer feed-forward neural network model embedded in the Bayesian nonparametric inference, the structure of NS matter's sound speed c_s is explored in a model-agnostic way. It is found that a peak structure is common in the c_s² posterior, locating at 2.4-4.8ρ_sat (nuclear saturation density) and c_s² exceeds c²/3 at 90% credibility. The non-monotonic behavior suggests evidence of the state deviating from the hadronic matter inside the very massive NSs. Assuming the new/exotic state is featured as it is softer than typical hadronic models or even with hyperons, we find that a sizable (⩾10^-3M_⊙) exotic core, likely made of quark matter, is plausible for the NS with a gravitational mass above about 0.98M_TOV, where M_TOV represents the maximum gravitational mass of a non-rotating cold NS. The inferred M_TOV=2.18_-0.13^+0.27M_⊙ (90% credibility) is well consistent with the value of 2.17_-0.12^+0.15M_⊙ estimated independently with GW170817/GRB 170817A/AT2017gfo assuming a temporary supramassive NS remnant formed after the merger. PSR J0740 + 6620, the most massive NS detected so far, may host a sizable exotic core with a probability of ≈0.36.

Energy-Density Modeling of Strongly Interacting Matter: Atomic Nuclei and Dense Stars

We seek a simple but physically motivated model of strongly interacting matter applicable in atomic nuclei and the dense matter in the core of neutron stars. For densities below and somewhat above normal nuclear density, energy density functional (EDF) theory based on nucleonic degrees of freedom is the ideal candidate. We have explored that direction within the KIDS (Korea-IBS-Daegu-SKKU) framework, which we review in this contribution. The formalism for the KIDS-EoS and microscopic KIDS-EDF and optimization options for the EDF are described in a practical way to facilitate further applications. At densities higher than one nucleon per single-nucleon volume, i.e., roughly 0.4 fm−3, nucleonic degrees of freedom are no longer appropriate. The pseudo-conformal symmetry emergent in dense, topologically altered nuclear matter provides a simple expression for the energy per baryon in terms of the baryonic density. Besides resembling a simple EDF for dense matter, the expression has the appeal that it predicts a converged speed of sound at high densities. It can thus be implemented as a special case of the constant speed of sound (CSS) model. Here we consider a matching between representative nucleonic KIDS-EoSs and the CSS model, including the pseudo-conformal EoS, and apply the unified model to describe the mass–radius relation of neutron stars and examine the compatibility of CSS cores with heavy neutron stars. Although an abrupt transition to the pseudo-conformal regime at low densities does not favor heavy neutron stars, intermediate scenarios including a cusp in the speed of sound are not ruled out, while some appear more favorable to heavy stars than purely nucleonic matter.

A physics-based neural network reconstruction of the dense matter equation of state from neutron star observables

We introduce a novel technique that utilizes a physics-driven deep learning method to reconstruct the dense matter equation of state from neutron star observables, particularly the masses and radii. The proposed framework involves two neural networks: one to optimize the EoS using Automatic Differentiation in the unsupervised learning scheme; and a pre-trained network to solve the Tolman–Oppenheimer–Volkoff (TOV) equations. The gradient-based optimization process incorporates a Bayesian picture into the proposed framework. The reconstructed EoS is proven to be consistent with the results from conventional methods. Furthermore, the resulting tidal deformation is in agreement with the limits obtained from the gravitational wave event, GW170817.

Cold Quark–Gluon Plasma EOS Applied to a Magnetically Deformed Quark Star with an Anomalous Magnetic Moment

We consider a QCD cold-plasma-motivated Equation of State (EOS) to examine the impact of an Anomalous Magnetic Moment (AMM) coupling and small shape deformations on the static oblate and prolate core shapes of quark stars. Using the Fogaça QCD-motivated EOS, which shifts from the high-temperature, low-chemical-potential quark–gluon plasma environment to the low-temperature, high-chemical-potential quark stellar core environment, we consider the impact of an AMM coupling with a metric-induced shape deformation parameter in the Tolman–Oppenheimer–Volkov (TOV) equations. The AMM coupling includes a phenomenological scaling that accounts for the weak and strong field characteristics in dense matter. The EOS is developed using a hard gluon and soft gluon decomposition of the gluon field tensor and using a mean-field effective mass for the gluons. The AMM is considered using the Dirac spin tensor coupled to the EM field tensor with quark-flavor-based magnetic moments. The shape parameter is introduced in a metric ansatz that represents oblate and prolate static stellar cores for modified TOV equations. These equations are numerically solved for the final mass and radius states, representing the core collapse of a massive star with a phase transition leading to an unbound quark–gluon plasma. We find that the combined shape parameter and AMM effects can alter the coupled EOS–TOV equations, resulting in an increase in the final mass and a decrease in the final equatorial radius without collapsing the core into a black hole and without violating causality constraints; we find maximum mass values in the range 1.6 Mʘ < M < 2.5 Mʘ. These states are consistent with some astrophysical, high-mass magnetar/pulsar and gravity wave systems and may provide evidence for a core that has undergone a quark–gluon phase transition such as PSR 0943 + 10 and the secondary from the GW 190814 event.

Constraining neutron-star matter with microscopic and macroscopic collisions

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not probed only in astrophysical observations, but also in terrestrial heavy-ion collision experiments. Here we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars1-9 and from heavy-ion collisions of gold nuclei at relativistic energies10,11 with microscopic nuclear theory calculations12-17 to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent observations by the Neutron Star Interior Composition Explorer mission5-8,18. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

Status and Perspectives of Continuous Gravitational Wave Searches

The birth of gravitational wave astronomy was triggered by the first detection of a signal produced by the merger of two compact objects (also known as a compact binary coalescence event). The following detections made by the Earth-based network of advanced interferometers had a significant impact in many fields of science: astrophysics, cosmology, nuclear physics and fundamental physics. However, compact binary coalescence signals are not the only type of gravitational waves potentially detectable by LIGO, Virgo, and KAGRA. An interesting family of still undetected signals, and the ones that are considered in this review, are the so-called continuous waves, paradigmatically exemplified by the gravitational radiation emitted by galactic, fast-spinning isolated neutron stars with a certain degree of asymmetry in their mass distribution. In this work, I will review the status and the latest results from the analyses of advanced detector data.

Oscillating Magnetized Color Superconducting Quark Stars

The main objective of this work is to study the structure, composition, and oscillation modes of color superconducting quark stars with intense magnetic fields. We adopted the MIT bag model within the color superconductivity CFL framework, and we included the effects of strong magnetic fields to construct the equation of state of stable quark matter. We calculated observable quantities, such as the mass, radius, frequency, and damping time of the oscillation fundamental f mode of quark stars, taking into account current astrophysical constraints. The results obtained show that color superconducting magnetized quark stars satisfy the constraints imposed by the observations of massive pulsars and gravitational wave events. Furthermore, the quantities associated with the oscillation f mode of these objects fit the universal relationships for compact objects. In the context of the new multi-messenger gravitational wave astronomy era and the future asteroseismology of neutron stars, we hope that our results contribute to the understanding of the behavior of dense matter and compact objects.

Finding Structure in the Speed of Sound of Supranuclear Matter from Binary Love Relations

Analyses that connect observations of neutron stars with nuclear-matter properties can rely on equation-of-state insensitive relations. We show that the slope of the binary Love relations (between the tidal deformabilities of binary neutron stars) encodes the baryon density at which the speed of sound rapidly changes. Twin stars lead to relations that present a signature "hill," "drop," and "swoosh" due to the second (mass-radius) stable branch, requiring a new description of the binary Love relations. Together, these features can reveal new properties and phases of nuclear matter.

Ranking Love Numbers for the Neutron Star Equation of State: The Need for Third-Generation Detectors

Gravitational-wave measurements of the tidal deformability in neutron-star binary coalescences can be used to infer the still unknown equation of state (EOS) of dense matter above the nuclear saturation density. By employing a Bayesian-ranking test, we quantify the ability of current and future gravitational-wave observations to discriminate among families of nuclear-physics based EOS which differ in particle content and ab initio microscopic calculations. While the constraining power of GW170817 is limited, we show that even twenty coalescences detected by LIGO-Virgo at design sensitivity are not enough to discriminate between EOS with similar softness but distinct microphysics. However, just a single detection with a third-generation detector such as the Einstein Telescope or Cosmic Explorer will rule out several families of EOS with very strong statistical significance and can discriminate among models which feature similar softness, hence, constraining the properties of nuclear matter to unprecedented levels.

Neutron Stars with Baryon Number Violation, Probing Dark Sectors

The neutron lifetime anomaly has been used to motivate the introduction of new physics with hidden-sector particles coupled to baryon number, and on which neutron stars provide powerful constraints. Although the neutron lifetime anomaly may eventually prove to be of mundane origin, we use it as motivation for a broader review of the ways that baryon number violation, be it real or apparent, and dark sectors can intertwine and how neutron star observables, both present and future, can constrain them.

Can we distinguish quark stars from neutron stars with measurements of global properties?

The phase state of the dense stellar matter is an exciting topic in the area of nuclear astrophysics. It may be probed by observed properties of neutron stars from, for example, the currently operating satellites (NICER, Neutron star Interior Composition Explorer) and the gravitational-wave laser interferometers (Advanced LIGO, Virgo, and KAGRA). Based on our recent constrained parameter spaces of the equation of states of neutron stars and quark stars from LIGO/Virgo and NICER, we discuss the important role of an even-accurate determination of the stellar radius for distinguishing possible quark stars from neutron stars and our understanding of the QCD phase transition at finite density.

The Equation of State of Neutron-Rich Matter at Fourth Order of Chiral Effective Field Theory and the Radius of a Medium-Mass Neutron Star

We report neutron star predictions based on our most recent equations of state. These are derived from chiral effective field theory, which allows for a systematic development of nuclear forces, order by order. We utilize high-quality two-nucleon interactions and include all three-nucleon forces up to fourth order in the chiral expansion. Our ab initio predictions are restricted to the domain of applicability of chiral effective field theory. However, stellar matter in the interior of neutron stars can be up to several times denser than normal nuclear matter at saturation, and its composition is essentially unknown. Following established practices, we extend our microscopic predictions to higher densities matching piecewise polytropes. The radius of the average-size neutron star, about 1.4 solar masses, is sensitive to the pressure at normal densities, and thus it is suitable to constrain ab initio theories of the equation of state. For this reason, we focus on the radius of medium-mass stars. We compare our results with other theoretical predictions and recent constraints.

Magnetohydrodynamic Simulations of Self-Consistent Rotating Neutron Stars with Mixed Poloidal and Toroidal Magnetic Fields

We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the development of various instabilities known from analytic perturbative studies, all our models finally succumb to them. Differential rotation is developed spontaneously in the cores of our magnetars which, after sufficient time, is converted back to uniform rotation. The rapidly rotating magnetars show a significant amount of ejecta, which can be responsible for transient kilonova signatures. However, no highly collimated, helical magnetic fields or incipient jets, which are necessary for γ-ray bursts, arise at the poles of these magnetars by the time our simulations are terminated.

Hyperon-nucleon femtoscopy, nuclear production and bearing on astrophysics

The presence of hyperons in the inner core of neutron stars is still an open issue in astrophysics and the so-called hyperon puzzle is far away from being solved. A complete understanding of how the interaction between hyperons and nucleons, at the level of two- and three-body, behaves in vacuum and at finite baryon density is a key requirement in order to build a realistic Equation of State for matter inside the core of neutron stars. Recently, measurements of hadron-hadron correlation functions with the femtoscopy technique in small colliding systems delivered a large amount of precise data on several hyperonnucleon pairs. Several interactions relevant for the physics of neutron stars, such as pΛ and pΞ, have been compared to the femtoscopic data and the effects on the corresponding Equation of State have been investigated. The extension to the three-body forces with the measurement of correlations involving triplets and proton-deuteron pairs will provide for the first time access to the genuine three-body interaction and shed light on the hyperon-puzzle.

Probing Compressed Baryonic Matter

The objective of the compressed baryonic matter (CBM) experiment at the future Facility for Antiproton and Ion Research (FAIR) in Darmstadt, Germany, is the investigation of the fundamental properties of strongly interacting matter. Of particular interest for our understanding of compact stellar objects is the determination of the equation-of-state (EOS) at high baryon densities and the exploration of the microscopic degrees-of-freedom under these conditions. The results of these laboratory experiments will complement astronomical observations, which also constrain the high-density EOS. Recent results of QCD-based calculations suggest that a possible first-order chiral phase transition should be observable in heavy-ion collisions at FAIR energies. This article reviews relevant observables from heavy-ion collisions and describes the detector configuration and the physics performance of the CBM experiment.

Dark Matter Effects on the Compact Star Properties

The neutron star properties are generally determined by the equation of state of β-equilibrated dense matter. In this work, we consider the interaction of fermionic dark matter (DM) particles with nucleons via Higgs exchange and investigate the effect on the neutron star properties with the relativistic mean-field model equation of state coupled with DM. We deduce that DM significantly affects the neutron star properties, such as considerably reducing the maximum mass of the star, which depends on the percentage of the DM considered inside the neutron star. The tidal Love numbers both for electric and magnetic cases and surficial Love numbers are also studied for DM admixed NS. We observed that the magnitude of tidal and surficial Love numbers increases with a greater DM percentage. Further, we present post-Newtonian tidal corrections to gravitational waves decreased by increasing the DM percentage. The DM effect on the GW signal is significant during the late inspiral and merger stages of binary evolution for GW frequencies >500 Hz.

Translating Neutron Star Observations to Nuclear Symmetry Energy via Deep Neural Networks

One of the most significant challenges involved in efforts to understand the equation of state of dense neutron-rich matter is the uncertain density dependence of the nuclear symmetry energy. In particular, the nuclear symmetry energy is still rather poorly constrained, especially at high densities. On the other hand, detailed knowledge of the equation of state is critical for our understanding of many important phenomena in the nuclear terrestrial laboratories and the cosmos. Because of its broad impact, pinning down the density dependence of the nuclear symmetry energy has been a long-standing goal of both nuclear physics and astrophysics. Recent observations of neutron stars, in both electromagnetic and gravitational-wave spectra, have already constrained significantly the nuclear symmetry energy at high densities. The next generation of telescopes and gravitational-wave observatories will provide an unprecedented wealth of detailed observations of neutron stars, which will improve further our knowledge of the density dependence of nuclear symmetry energy, and the underlying equation of state of dense neutron-rich matter. Training deep neural networks to learn a computationally efficient representation of the mapping between astrophysical observables of neutron stars, such as masses, radii, and tidal deformabilities, and the nuclear symmetry energy allows its density dependence to be determined reliably and accurately. In this work, we use a deep learning approach to determine the nuclear symmetry energy as a function of density directly from observational neutron star data. We show, for the first time, that artificial neural networks can precisely reconstruct the nuclear symmetry energy from a set of available neutron star observables, such as masses and radii as measured by, e.g., the NICER mission, or masses and tidal deformabilities as measured by the LIGO/VIRGO/KAGRA gravitational-wave detectors. These results demonstrate the potential of artificial neural networks to reconstruct the symmetry energy and the equation of state directly from neutron star observational data, and emphasize the importance of the deep learning approach in the era of multi-messenger astrophysics.

Equations of State for Hadronic Matter and Mass-Radius Relations of Neutron Stars with Strong Magnetic Fields

Neutron star is an important object for us to verify the equation of state of hadronic matter. For a specific choice of equations of state, mass and radius of a neutron star are determined, for which there are constraints from observations. According to some previous studies, since the strong magnetic field acts as a repulsive force, there is a possibility that neutron stars with strong magnetic fields may have relatively heavier masses than other non-magnetized neutron stars. In this paper, the structure of a neutron star with a strong internal magnetic field is investigated by changing its internal functional form to see how much the neutron star can be massive and also how radius of a neutron star can be within a certain range.

Confining density functional approach for color superconducting quark matter and mesonic correlations

We present a novel relativistic density-functional approach to modeling quark matter with a mechanism to mimic confinement. The quasiparticle treatment of quarks provides their suppression due to large quark selfenergy already at the mean-field level. We demonstrate that our approach is equivalent to a chiral quark model with medium-dependent couplings. The dynamical restoration of the chiral symmetry is ensured by construction of the density functional. Beyond the mean field, quark correlations in the pseudoscalar channel are described within the Gaussian approximation. This explicitly introduces pionic states into the model. Their contribution to the thermodynamic potential is analyzed within the Beth–Uhlenbeck framework. The modification of the meson mass spectrum in the vicinity of thee (de)confinement transition is interpreted as the Mott transition. Supplemented with the vector repulsion and diquark pairing the model is applied to construct a hybrid quark-hadron EoS of cold compact-star matter. We study the connection of such a hybrid EoS with the stellar mass-radius relation and tidal deformability. The model results are compared to various observational constraints including the NICER radius measurement of PSR J0740+6620 and the tidal deformability constraint from GW170817. The model is shown to be consistent with the constraints, still allowing for further improvement by adjusting the vector repulsion and diquark pairing couplings.

Hybrid stars from a constrained equation of state

We determine, within a meta-model, the properties of the nuclear matter equation of state (EoS) that allow for a phase transition to deconfinement matter. It is shown that the properties that define the isoscalar channel are the ones that are affected, in particular, a phase transition implies much larger values of the skewness and kurtosis. The effect of multi-quark interaction channels in the description of the quark phase in hybrid stars is also studied. NS properties, such as the mass and radius of the quark core, show an interplay dependence between the 8-quark vector and the 4-quark isovector-vector interactions. We show that low mass NS, M ~ 1.4M⊙, may already contain a quark core, and satisfy all existing NS observational constraints. We discuss the strangeness content of the quark core and its influence on the speed of sound.

Conformality and percolation threshold in neutron stars

Speed of sound is given attention in multi-messenger astronomy as it encodes information of the dense matter equation of state. Recently the trace anomaly was proposed as a more informative quantity. In this work, we statistically determine the speed of sound and trace anomaly and show that they are driven to their conformal values at the centers of maximally massive neutron stars. We show that the local peak in the speed of sound can be associated deconfinement along with percolation conditions in QCD matter.

Exploring the role of hyperons in high density matter in the Quark-Meson-Coupling model

Microscopic composition and properties of matter at super-saturation densities have been a subject of intense investigations for decades. Experimental and observational data and fundamental laws of physics indicate that heavy strange baryons and mesons are essential components of the matter. The Quark-Meson-Coupling-Model (QMC) is well suited for such a study. The model is based on interaction between quarks in individual baryons instead between the baryons as entities without internal structure. This approach significantly increases transparency and reduces the number of variable parameters of the model, thus offering deeper insight into the physics of high density hadronic matter. In this contribution, we review the effect of hyperons on neutron star (NS) properties, the speed of sound and the symmetry energy, both at zero and finite temperature. The QMC results are contrasted with the outcome of the traditional relativistic mean field DD2Y model. The (lack of) the so-called ‘hyperon puzzle’ in both models is discussed.

Exploring the Phase Diagram of V-QCD with Neutron Star Merger Simulations

Determining the phase structure of Quantum Chromodynamics (QCD) and its Equation of State (EOS) at densities and temperatures realized inside neutron stars and their mergers is a long-standing open problem. The holographic V-QCD framework provides a model for the EOS of dense and hot QCD, which describes the deconfinement phase transition between a dense baryonic and a quark matter phase. We use this model in fully general relativistic hydrodynamic (GRHD) simulations to study the formation of quark matter and the emitted gravitational wave signal of binary systems that are similar to the first ever observed neutron star merger event GW170817.

How does dark matter affect compact star properties and high density constraints of strongly interacting matter

We study the impact of asymmetric bosonic dark matter on neutron star properties, including possible changes of tidal deformability, maximum mass, radius, and matter distribution inside the star. The conditions at which dark matter particles tend to condensate in the star’s core or create an extended halo are presented. We show that dark matter condensed in a core leads to a decrease of the total gravitational mass and tidal deformability compared to a pure baryonic star, which we will perceive as an effective softening of the equation of state. On the other hand, the presence of a dark matter halo increases those observable quantities. Thus, observational data on compact stars could be affected by accumulated dark matter and, consequently, constraints we put on strongly interacting matter at high densities. To confirm the presence of dark matter in the compact star’s interior, and to break the degeneracy between the effect of accumulated dark matter and strongly interacting matter properties at high densities, several astrophysical and GW tests are proposed.

Hybrid stars with large strange quark cores

The possible existence of hybrid stars is studied using several multiquark interaction channels. The hadronic phase consists of an equation of state (EoS) with presently accepted nuclear matter properties and the quark model is constrained by the vacuum properties of several light mesons. The dependence of several NS properties on the different quark interactions is analyzed. We show that the present constraints from neutron star observations allow for the existence of hybrid stars with large strangeness content and large quark cores.

Early deconfinement of asymptotically conformal color-superconducting quark matter in neutron stars

We present a relativistic density functional approach to color superconducting quark matter that mimics quark confinement by a fast growth of the quasiparticle selfenergy in the confining region. The approach is shown to be equivalent to a chiral model of quark matter with medium dependent couplings. While the (pseudo)scalar sector of the model is fitted to the vacuum phenomenology of quantum chromodynamics, the strength of interaction in the vector and diquark channels is varied in order to provide the best agreement with the observational constraints on the mass-radius relation and tidal deformability of neutron stars modelled with our approach. In order to recover the conformal behavior of quark matter at asymptotically high densities we introduce a medium dependence of the vector and diquark couplings motivated by the non-perturbative gluon exchange. Our analysis signals that the onset of deconfinement to color superconducting quark matter is likely to occur in neutron stars with masses below 1.0 M⊙.

Finite-temperature equation of state with hyperons

We present the novel finite-temperature FSU2H* equation-of-state model that covers a wide range of temperatures and lepton fractions for the conditions in proto-neutron stars, neutron star mergers and supernovae. The temperature effects on the thermodynamical observables and the composition of the neutron star core are stronger when the hyperonic degrees of freedom are considered. We pay a special attention to the temperature and density dependence of the thermal index in the presence of hyperons and conclude that the true thermal effects cannot be reproduced with the use of a constant Г law

Hyperonic equation of state at finite temperature for neutron stars

We review the composition and the equation of state of the hyperonic core of neutron stars at finite temperature within a relativistic mean-field approach. We make use of the new FSU2H∗ model, which is built upon the FSU2H scheme by improving on the Ξ potential according to the recent analysis on the Ξ atoms, and we extend it to include finite temperature corrections. The calculations are done for a wide range of densities, temperatures and charge fractions, thus exploring the different conditions that can be found in protoneutron stars, binary mergers remnants and supernovae explosions. The inclusion of hyperons has a strong effect on the composition and the equation of state at finite temperature, which consequently would lead to significant changes in the properties and evolution of hot neutron stars.

Phase Conversions in Neutron Stars: Implications for Stellar Stability and Gravitational Wave Astrophysics

We review the properties of hybrid stars with a quark matter core and a hadronic mantle, focusing on the role of key micro-physical properties such as the quark/hadron surface and curvature tensions and the conversion speed at the interface between both phases. We summarize the results of works that have determined the surface and curvature tensions from microscopic calculations. If these quantities are large enough, mixed phases are energetically suppressed and the quark core would be separated from the hadronic mantle by a sharp interface. If the conversion speed at the interface is slow, a new class of dynamically stable hybrid objects is possible. Densities tens of times larger than the nuclear saturation density can be attained at the center of these objects. We discuss possible formation mechanisms for the new class of hybrid stars and smoking guns for their observational identification.

1S0 Pairing Gaps, Chemical Potentials and Entrainment Matrix in Superfluid Neutron-Star Cores for the Brussels–Montreal Functionals

Temperature and velocity-dependent 1S0 pairing gaps, chemical potentials and entrainment matrix in dense homogeneous neutron–proton superfluid mixtures constituting the outer core of neutron stars, are determined fully self-consistently by solving numerically the time-dependent Hartree–Fock–Bogoliubov equations over the whole range of temperatures and flow velocities for which superfluidity can exist. Calculations have been made for npeμ in beta-equilibrium using the Brussels–Montreal functional BSk24. The accuracy of various approximations is assessed and the physical meaning of the different velocities and momentum densities appearing in the theory is clarified. Together with the unified equation of state published earlier, the present results provide consistent microscopic inputs for modeling superfluid neutron-star cores.

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of ^{208}Pb with Minimal Modeling Assumptions

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283±0.071  fm for the neutron-skin thickness of ^{208}Pb, implying a slope parameter L=106±37  MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S_{0}, L, and R_{skin}^{^{208}Pb} directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Combining astrophysical data with PREX-II and chiral effective field theory constraints yields S_{0}=33.0_{-1.8}^{+2.0}  MeV, L=53_{-15}^{+14}  MeV, and R_{skin}^{^{208}Pb}=0.17_{-0.04}^{+0.04}  fm.

Hyperons in Neutron Stars

I review the issues related to the appearance of hyperons in neutron star matter, focusing in particular on the problem of the maximum mass supported by hyperonic equations of state. I discuss the general mechanism that leads to the formation of hyperons in the core of neutron stars and I review the main techniques and many-body methods used to construct an appropriate equation of state to describe the strongly interacting system of hadrons hosted in the core of neutron stars. I outline the consequences on the structure and internal composition of neutron stars and also discuss the possible signatures of the presence of hyperons in astrophysical dynamical systems like supernova explosions and binary neutron star mergers. Finally, I briefly report about the possible important role played by hyperons in the transport properties of neutron star matter and on the consequences of neutron star cooling and gravitational wave instabilities induced by the presence of hyperons.