The thermodynamic properties of high temperature and high density QCD-matter are studied using the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The CMF model provides a proper description of lattice QCD data, heavy-ions physics, and static neutron stars. The behavior of lines of constant pressure with increase of baryon density is discussed. The rapid change of pressure behavior at $\mu_B/T\approx3$ suggests a strong contribution of baryons to thermodynamic properties at this region. The position of this region is very close to the radius of convergence for a Taylor expansion of the QCD pressure. The role of mesons and unstable hadrons in the hydrodynamic expansion of strongly interacting matter is also discussed.

The gravitational waves emitted from a binary neutron star merger, as predicted from general relativistic magneto-hydrodynamics calculations, are sensitive to the appearance of quark matter and the stiffness of the equation of state of QCD matter present in the inner cores of the stars. This is a new messenger observable from outer space, which does provide direct signals for the phase structure of strongly interacting QCD matter at high baryon density and high temperature. These astrophysically created extremes of thermodynamics do match, to within 20\%, the values of densities and temperatures which we find in relativistic hydrodynamics and transport theory of heavy ion collisions at the existing laboratories, if though at quite different rapidity windows, impact parameters and bombarding energies of the heavy nuclear systems. We demonstrate how one unified equation of state can be constructed and used for both neutron star physics and hot QCD matter excited at laboratory facilities. The similarity in underlying QCD physics allows the gravitational wave signals from future advanced LIGO and Virgo events to be combined with the analysis of high multiplicity fluctuations and flow measurements in heavy ion detectors in the lab to pin down the EoS and the phase structure of dense matter.

In this work we study the structure of neutron stars under the effect of a poloidal magnetic field and determine the limiting highest magnetic field intensity which still allows a satisfactory description of magnetic neutron stars in the spherical symmetry regime. We describe different compositions of stars (nucleonic, hyperonic, and hybrid), using three state-of-the-art relativistic mean field models for the microscopic description of matter, which are in agreement with experimental and observational data. The structure of stars is described by the general relativistic solution of both Einstein's field equations assuming a spherical symmetry, and Einstein-Maxwell's field equations assuming an axi-symmetric deformation. We find a limiting magnetic moment of the order of $2\times 10^{31}$Am$^2$, which corresponds to magnetic fields of the order of 10$^{16}$ G at the surface, and $ \sim 10^{17}$ G at the centre of the star, above which the deformation due to the magnetic field is not negligible. We show that the intensity of the magnetic field developed in the star depends on the EoS, and, for a given baryonic mass and fixed magnetic moment, larger fields are attained with softer EoS. We also show that the appearance of exotic degrees of freedom, such as hyperons or a quark core, is disfavored in the presence of a very strong magnetic field. As a consequence, a highly magnetized nucleonic star may suffer an internal conversion due to the decay of the magnetic field, which could be accompanied by a sudden cooling of the star or a gamma ray burst.

In this work, we discuss the dense matter equation of state (EOS) for the extreme range of conditions encountered in neutron stars and their mergers. The calculation of the properties of such an EOS involves modeling different degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking into account different symmetries, and including finite density and temperature effects in a thermodynamically consistent manner. We begin by addressing subnuclear matter consisting of nucleons and a small admixture of light nuclei in the context of the excluded volume approach. We then turn our attention to supranuclear homogeneous matter as described by the Chiral Mean Field (CMF) formalism. Finally, we present results from realistic neutron-star-merger simulations performed using the CMF model that predict signatures for deconfinement to quark matter in gravitational wave signals.

The thermodynamic properties of high temperature and high density QCD-matter are explored within the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The quark sector of the CMF model is tuned to describe the $\mu_B=0$ thermodynamics data of lattice QCD. The resulting lines of constant physical variables as well as the baryon number susceptibilities are studied in some detail in the temperature/chemical potential plane. The CMF model predicts three consecutive transitions, the nuclear first-order liquid-vapor phase transition, chiral symmetry restoration, and the cross-over transition to a quark-dominated phase. All three phenomena are cross-over, for most of the $T-\mu_B$-plane. The deviations from the free ideal hadron gas baseline at $\mu_B=0$ and $T\approx 100-200$ MeV can be attributed to remnants of the liquid-vapor first order phase transition in nuclear matter. The chiral crossing transition determines the baryon fluctuations at much higher $\mu_B\approx1.5$ GeV, and at even higher baryon densities $\mu_B\approx2.4$ GeV, the behavior of fluctuations is controlled by the deconfinement cross-over. The CMF model also describe well the static properties of high $\mu_B$ neutron stars as well as the new neutron star merger observations. The effective EoS presented here describes simultaneously lattice QCD results at $\mu_B=0$, as well as observed physical phenomena (nuclear matter and neutron star matter) at $T\cong0$ and high densities, $\mu_B>1$ GeV.

Merging binaries of neutron stars are not only strong sources of gravitational waves, but also have the potential of revealing states of matter at densities and temperatures not accessible in laboratories. A crucial and long-standing question in this context is whether quarks are deconfined as a result of the dramatic increase in density and temperature following the merger. We present the first fully general-relativistic simulations of merging neutron stars including quarks at finite temperatures that can be switched off consistently in the equation of state. Within our approach, we can determine clearly what signatures a quark-hadron phase transition would leave in the gravitational-wave signal. In particular, we show that if the conditions are met for a phase transition to take place at several times nuclear saturation density, they would lead to a post-merger signal considerably different from the one expected from the inspiral, that can only probe the hadronic part of the equations of state, and to an anticipated collapse of the merged object. We also show that the phase transition leads to a very hot and dense quark core that, when it collapses to a black hole, produces a ringdown signal different from the hadronic one. Finally, in analogy with what is done in heavy-ion collisions, we use the evolution of the temperature and density in the merger remnant to illustrate the properties of the phase transition in a QCD phase diagram.

We analyze the role played by vector-isovector meson interaction in dense matter present in the interior of neutron stars in the light of new measurements made during the double neutron-star merger GW170817. These concern measurements of tidal deformability from gravitational waves and electromagnetic observations. Our study includes three different equations of state that contain different physical assumptions and matter compositions, namely the NL3 family, MBF, and CMF models. Other related quantities/relations analyzed are the neutron matter pressure, symmetry energy slope, stellar masses and radii, and Urca process threshold for stellar cooling.

In this paper we review the most common descriptions for the first order phase transition to deconfined quark matter in the core of neutron stars. We also present a new description of these phase transitions in the core of proto-neutron stars, in which more constraints are enforced so as to include trapped neutrinos. Finally, we calculate the emission of gravitational waves associated with deconfinement phase transitions, discuss the possibility of their detection, and how this would provide information about the equation of state of dense matter.

Sharp phase transitions described by stiff equations of state allow for the existence of a third family of stable compact stars (besides white dwarfs and neutron stars), twin stars. In this work, we investigate for the first time the role of strong magnetic fields on non-magnetic twin stars sequences and the case in which magnetic fields themselves give rise to a third family of stable stars. We use three sets of equations of state to study such effects from a more general point of view: the Quark-Hadron Chiral Parity-Doublet (Q$\chi$P) model for both hadronic and quark phases, and the Many-Body Forces (MBF) model connected to either the MIT Bag model with vector interaction (MIT) or to the Constant-Sound-Speed (CSS) approximation for the quark phase, through a Maxwell construction. Magnetic field effects are introduced in the structure of stars through the solution of the Einstein-Maxwell equations, assuming a poloidal magnetic field configuration and a metric that allows for the description of deformed stars. We show that strong magnetic fields can destabilize twin star sequences, with the threshold intensity being model dependent. On the other hand, magnetic fields can also give rise to twin stars in models that did not predict these sequences, up to some point when they are again destabilized. In this sense, magnetic fields can play an important role on the evolution of neutron stars.

A unified chiral mean field approach is presented for QCD thermodynamics in a wide range of temperatures and densities. The model simultaneously gives a satisfactory description of lattice QCD thermodynamics and fulfills nuclear matter and astrophysical constraints. The resulting equation of state can be incorporated in relativistic fluid-dynamical simulations of heavy-ion collisions and neutron stars mergers. Access to different regions of the QCD phase diagram can be obtained in simulations of heavy-ion data and observations of neutron star mergers.

Particles with strangeness content are predicted to populate dense matter, modifying the equation of state of matter inside neutron stars as well as their structure and evolution. In this work, we show how the modeling of strangeness content in dense matter affects the properties of isolated neutrons stars and the tidal deformation in binary systems. For describing nucleonic and hyperonic stars we use the many-body forces model (MBF) at zero temperature, including the $\phi$ mesons for the description of repulsive hyperon-hyperon interactions. Hybrid stars are modeled using the MIT Bag Model with vector interaction (vMIT) in both Gibbs and Maxwell constructions, for different values of bag constant and vector interaction couplings. A parametrization with a Maxwell construction, which gives rise to third family of compact stars (twin stars), is also investigated. We calculate the tidal contribution that adds to the post-Newtonian point-particle corrections, the associated love number for sequences of stars of different composition (nucleonic, hyperonic, hybrid and twin stars), and determine signatures of the phase transition on the gravitational waves in the accumulated phase correction during the inspirals among different scenarios for binary systems. On the light of the recent results from GW170817 and the implications for the radius of $\sim1.4\,\mathrm{M_{\odot}}$ stars, our results show that hybrid stars can only exist if a phase transition takes place at low densities close to saturation.

Neutron stars are the endpoint of the life of intermediate mass stars and posses in their cores matter in the most extreme conditions in the universe. Besides their extremes of temperature (found in proto-neutron stars) and densities, typical neutron star' magnetic fields can easily reach trillions of times higher the one of the Sun. Among these stars, about $10\%$ are denominated \emph{magnetars} which possess even stronger surface magnetic fields of up to $10^{15}-10^{16}\,\mathrm{G}$. In this conference proceeding, we present a short review of the history and current literature regarding the modeling of magnetic neutron stars. Our goal is to present the results regarding the introduction of magnetic fields in the equation of state of matter using Relativistic Mean Field models (RMF models) and in the solution of Einstein's equations coupled to the Maxwell's equations in order to generate a consistent calculation of magnetic stars structure. We discuss how equation of state modeling affects mass, radius, deformation, composition and magnetic field distribution in stars and also what are the open questions in this field of research.

We study the transport properties of nuclear pasta for a wide range of density, temperature and proton fractions, relevant for different astrophysical scenarios adopting a quantum molecular dynamics model. In particular, we estimate the values of shear viscosity as well as electrical and thermal conductivities by calculating the static structure factor $S(q)$ using simulation data. In the density and temperature range where the pasta phase appears, the static structure factor shows irregular behavior. The presence of a slab phase greatly enhances the peak in $S(q)$. However, the effect of irregularities in $S(q)$ on the transport coefficients is not very dramatic. The values of all three transport coefficients are found to have the same orders of magnitude as found in theoretical calculations for the inner crust matter of neutron stars without the pasta phase and therefore, is in contrast to earlier speculations that a pasta layer might be highly resistive, both thermally and electrically.

The Quark-Hadron Chiral Parity-Doublet model (Q$\chi$P) is applied to calculate compact star properties in the presence of a deconfinement phase transition. Within this model, a consistent description of nuclear matter properties, chiral symmetry restoration, and a transition from hadronic to quark and gluonic degrees of freedom is possible within one unified approach. We find that the equation of state obtained is consistent with recent perturbative quantum chromodynamics (QCD) results and is able to accommodate observational constraints of massive and small neutron stars. Furthermore, we show that important features of the equation of state, such as the symmetry energy and its slope, are well within their observational constraints.

In this work, we expand on a previously reported realistic calculation of the magnetic field profile for the equation of state inside strongly magnetized neutron stars. In addition to showing that magnetic fields increase quadratically with increasing baryon chemical potential of magnetized matter (instead of exponentially, as previously assumed), we show here that the magnetic field increase with baryon number density is more complex and harder to model. We do so by the analysis of several different realistic models for the microscopic description of matter in the star (including hadronic, hybrid and quark models) combined with general relativistic solutions by solving Einstein-Maxwell's field equations in a self-consistent way for stars endowed with a poloidal magnetic field.

In this Letter, we report a realistic calculation of the magnetic field profile for the equation of state inside strongly magnetized neutron stars. Unlike previous estimates, which are widely used in the literature, we find that magnetic fields increase relatively slowly with increasing baryon chemical potential (or baryon density) of magnetized matter. More precisely, the increase is polynomial instead of exponential, as previously assumed. Through the analysis of several different realistic models for the microscopic description of stellar matter (including hadronic, hybrid and quark models) combined with general relativistic solutions endowed with a poloidal magnetic field obtained by solving Einstein-Maxwell's field equations in a self-consistent way, we generate a phenomenological fit for the magnetic field distribution in the stellar polar direction to be used as input in microscopic calculations.

We use an extended version of the hadronic SU(3) non-linear realization of the sigma model that also includes quarks to study hybrid stars. Within this approach, the degrees of freedom change naturally as the temperature/density increases. Different prescriptions of charge neutrality, local and global, are tested and the influence of strong magnetic fields and the anomalous magnetic moment on the particle population is discussed.

As the density of matter increases, atomic nuclei disintegrate into nucleons and, eventually, the nucleons themselves disintegrate into quarks. The phase transitions (PT's) between these phases can vary from steep first order to smooth crossovers, depending on certain conditions. First-order PT's with more than one globally conserved charge, so-called non-congruent PT's, have characteristic differences compared to congruent PT's. In this conference proceeding we discuss the non-congruence of the quark deconfinement PT at high densities and/or temperatures relevant for heavy-ion collisions, neutron stars, proto-neutron stars, supernova explosions, and compact-star mergers.

In this work, we study the properties of magnetized white dwarfs taking into account possible instabilities due to electron capture and pycnonuclear fusion reactions in the cores of such objects. The structure of white dwarfs is obtained by solving the Einstein-Maxwell equations with a poloidal magnetic field in a fully general relativistic approach. The stellar interior is composed of a regular crystal lattice made of carbon ions immersed in a degenerate relativistic electron gas. The onsets of electron capture reactions and pycnonuclear reactions are determined with and without magnetic fields. We find that magnetized white dwarfs violate the standard Chandrasekhar mass limit significantly, even when electron capture and pycnonuclear fusion reactions are present in the stellar interior. We obtain a maximum white dwarf mass of around $2.14\,M_{\odot}$ for a central magnetic field of $\sim 3.85\times 10^{14}$~G, which indicates that magnetized white dwarfs may play a role for the interpretation of superluminous type Ia supernovae. Furthermore, we show that the critical density for pycnonuclear fusion reactions limits the central white dwarf density to $9.35\times 10^9$ g/cm$^3$. As a result, equatorial radii of white dwarfs cannot be smaller than $\sim 1100$~km. Another interesting feature concerns the relationship between the central stellar density and the strength of the magnetic field at the core of a magnetized white dwarf. For high magnetic fields, we find that the central density increases (stellar radius decrease) with magnetic field strength, which makes ultramagnetized white dwarfs more compact. The opposite is the case, however, if the central magnetic field is less than $\sim 10^{13}$~G. In the latter case, the central density decreases (stellar radius increases) with central magnetic field strengths.

In this work, we study the effects of strong magnetic field configurations on the population of neutron stars. The stellar matter is described within a relativistic mean field formalism which considers many-body force contributions in the scalar couplings. We choose the parametrization of the model that reproduces nuclear matter properties at saturation and also describes massive hyperon stars. Hadronic matter is modeled at zero temperature, in beta-equilibrium, charge neutral and populated by the baryonic octet, electrons and muons. Magnetic effects are taken into account in the structure of stars by the solution of the Einstein-Maxwell equations with the assumption of a poloidal magnetic field distribution. Our results show that magnetic neutron stars are populated essencialy by nucleons and leptons, due to the fact that strong magnetic fields decrease the central density of stars and, hence, supress the appearance of exotic particles.

In this work, we study the effects of different magnetic field configurations in neutron stars described by a many-body forces formalism (MBF model). The MBF model is a relativistic mean field formalism that takes into account many-body forces by means of a meson field dependence of the nuclear interaction coupling constants. We choose the best parametrization of the model that reproduces nuclear matter properties at saturation and also describes massive neutron stars. We assume matter to be in beta-equilibrium, charge neutral and at zero temperature. Magnetic fields are taken into account both in the equation of state and in the structure of the stars by the self-consistent solution of the Einstein-Maxwell equations. We assume a poloidal magnetic field distribution and calculate its effects on neutron stars, showing its influence on the gravitational mass and deformation of the stars.

In view of the new recent observation and measurement of the rotating and highly-magnetized white dwarf AR Scorpii \cite{Marsh:2016uhc}, we determine bounds of its moment of inertia, magnetic fields and radius. Moreover, we investigate the possibility of fast rotating and/or magnetized white dwarfs to be source of detectable gravitational wave (GW) emission. Numerical stellar models at different baryon masses are constructed. For each star configuration, we compute self-consistent relativistic solutions for white dwarfs endowed with poloidal magnetic fields by solving the Einstein-Maxwell field equations in a self-consistent way. The magnetic field supplies an anisotropic pressure, leading to the braking of the spherical symmetry of the star. In this case, we compute the quadrupole moment of the mass distribution. Next, we perform an estimate of the GW of such objects. Finally, we show that the new recent observation and measurement pulsar white dwarf AR Scorpii, as well other stellar models, are able to generate gravitational wave radiation that lies in the bandwidth of the discussed next generation of space-based GW detectors DECI-hertz interferometer Gravitational wave Observatory (DECIGO) and Big Bang Observer (BBO).

We perform calculations of white dwarfs endowed with strong magnetic fields. White dwarfs are the progenitors of supernova Type Ia explosions and they are widely used as candles to show that the Universe is expanding and accelerating. However, observations of ultraluminous supernovae have suggested that the progenitor of such an explosion should be a white dwarf with mass above the well-known Chandrasekhar limit $\sim \,$1.4 $\rm{M_{\odot}}$. In corroboration with other works, but by using a fully general relativistic framework, we obtained also strongly magnetized white dwarfs with masses $\rm{M \sim 2.0\,\, M_{\odot}}$.

It is understood that strong magnetic fields affect the structure of neutron stars. Nevertheless, many calculations for magnetized neutron stars are still being performed using symmetric solutions of Einstein's equations. In this conference proceeding, we review why this is not the correct procedure and we also discuss the effects of magnetic fields on the stellar population and temperature profiles.

We study the effects of high magnetic fields on the structure and on the geometry of the crust in neutron stars. We find that the crust geometry is substantially modified by the magnetic field inside the star. We build stationary and axis-symmetric magnetized stellar models by using well-known equations of state to describe the neutron star crust, namely the Skyrme model (Sky) for the inner crust and the Baym, Pethick, and Sutherland (BPS) equation of state for the outer crust. We show that the magnetic field has a dual role, contributing to the crust deformation via the electromagnetic interaction (manifested in this case as the Lorentz force) and by contributing to curvature due to the energy stored in it. We also study a direct consequence of the crust deformation due to the magnetic field: the thermal relaxation time. This quantity, which is of great importance to the thermal evolution of neutron stars is sensitive to the crust properties and, as such, we show that it may be strongly affected by the magnetic field.

We study the effect of isospin-dependent nuclear forces on the pasta phase in the inner crust of neutron stars. To this end we model the crust within the framework of quantum molecular dynamics (QMD). For maximizing the numerical performance, a newly developed code has been implemented on GPU processors. As a first application of the crust studies we investigate the dependence of the particular pasta phases on the isospin dependence of the interaction, including non-linear terms in this sector of the interactions. Our results indicate that in contrast to earlier studies the phase diagram of the pasta phase is not very sensitive to isospin effects. We show that the extraction of the isospin parameters like asymmetry energy and slope from numerical data is affected by higher-order terms in the asymmetry dependence of the energies per particle. Furthermore, a rapid transition from the pasta to a homogeneous phase is observed even for proton-to-neutron ratios typical for a supernova environment.

The presence of quark-hadron phase transitions in neutron stars can be related to several interesting phenomena. In particular, previous calculations have shown that fast rotating neutron stars, when subjected to a quark-hadron phase transition in their interiors, could give rise to the backbending phenomenon characterized by a spin-up era. In this work, we use an equation of state composed of two phases, containing nucleons (and leptons) and quarks. The hadronic phase is described in a relativistic mean field formalism that takes many-body forces into account, and the quark phase is described by the MIT bag model with a vector interaction. Stationary and axi-symmetric stellar models are obtained in a self-consistent way by solving numerically the Einstein-Maxwell equations by means of a pseudo-spectral method. As a result, we obtain the interesting backbending phenomenon for fast spinning neutron stars. More importantly, we show that a magnetic field, which is assumed to be axi-symmetric and poloidal, can also be enhanced due to the phase transition from normal hadronic matter to quark matter on highly magnetized neutron stars. Therefore, in parallel to the spin-up era, classes of neutron stars endowed with strong magnetic fields may go through a "magnetic-up era" in their lives.

In this work, we study the effects of magnetic fields and rotation on the structure and composition of proto-neutron stars (PNS's). A hadronic chiral SU(3) model is applied to cold neutron stars (NS) and proto-neutron stars with trapped neutrinos and at fixed entropy per baryon. We obtain general relativistic solutions for neutron and proto-neutron stars endowed with a poloidal magnetic field by solving Einstein-Maxwell field equations in a self-consistent way. As the neutrino chemical potential decreases in value over time, this alters the chemical equilibrium and the composition inside the star, leading to a change in the structure and in the particle population of these objects. We find that the magnetic field deforms the star and significantly alters the number of trapped neutrinos in the stellar interior, together with strangeness content and temperature in each evolution stage.

In this work we compute models for relativistic white dwarfs in the presence of strong magnetic fields. These models possibly contribute to super-luminous SNIa. With an assumed axi-symmetric and poloidal magnetic field, we study the possibility of existence of super-Chandrasekhar magnetized white dwarfs by solving numerically the Einstein-Maxwell equations, by means of a pseudo-spectral method. We obtain a self-consistent rotating and non-rotating magnetized white dwarf models. According to our results, a maximum mass for a static magnetized white dwarf is 2.13 $\rm{M_{\odot}}$ in the Newtonian case and 2.09 $\rm{M_{\odot}}$ while taking into account general relativistic effects. Furthermore, we present results for rotating magnetized white dwarfs. The maximum magnetic field strength reached at the center of white dwarfs is of the order of $10^{15}\,$G in the static case, whereas for magnetized white dwarfs, rotating with the Keplerian angular velocity, is of the order of $10^{14}\,$G.

In this work we study the effects of strong magnetic fields on hybrid stars by using a full general-relativity approach, solving the coupled Maxwell-Einstein equation in a self-consistent way. The magnetic field is assumed to be axi-symmetric and poloidal. We take into consideration the anisotropy of the energy-momentum tensor due to the magnetic field, magnetic field effects on equation of state, the interaction between matter and the magnetic field (magnetization), and the anomalous magnetic moment of the hadrons. The equation of state used is an extended hadronic and quark SU(3) non-linear realization of the sigma model that describes magnetized hybrid stars containing nucleons, hyperons and quarks. According to our results, the effects of the magnetization and the magnetic field on the EoS do not play an important role on global properties of these stars. On the other hand, the magnetic field causes the central density in these objects to be reduced, inducing major changes in the populated degrees of freedom and, potentially, converting a hybrid star into a hadronic star.

In view of new constraints put forth by recent observations and measurements in the realm of astrophysics and nuclear physics, we update the non-linear realization of the sigma model as to reflect such constraints. By doing this, we obtain new equations of state that may be used to describe neutron stars. Such equations of state are obtained by investigating different ways by which the vector mesons self-interact. Furthermore, we also investigate the role played by the delta mesons in the model. As a result, we are able to develop equations of state that are in better agreement with data, such as nuclear compressibility and slope of the symmetry energy at saturation, star masses, radii, and cooling profiles.

We study the effects of strangeness on the quark sector of a hybrid star equation of state. Since the model we use to describe quarks is the same as the one we use to describe hadrons, we can also study the effects of strangeness on the chiral symmetry restoration and deconfinement phase transitions (first order or crossover). Finally, we analyze combined effects of hyperons and quarks on global properties of hybrid stars, like mass, radius and cooling profiles. It is found that a large amount of strangeness in the core is related to the generation of twin-star solutions, which can have the same mass as the lower or zero strangeness counterpart, but with smaller radii.

The equation of state (EoS) of asymmetric nuclear matter at high densities is a key topic for the description of matter inside neutron stars. The determination of the properties of asymmetric nuclear matter, such as the symmetry energy ($a_{sym}$) and the slope of the symmetry energy ($L_0$) at saturation density, has been exaustively studied in order to better constrain the nuclear matter EoS. However, differently from symmetric matter properties that are reasonably constrained, the symmetry energy and its slope still large uncertainties in their experimental values. Regarding this subject, some studies point towards small values of the slope of the symmetry energy, while others suggest rather higher values. Such a lack of agreement raised a certain debate in the scientific community. In this paper, we aim to analyse the role of these properties on the behavior of asymmetric hyperonic matter. Using the formalism presented in Ref. (R.O. Gomes et al 2014}, which considers many-body forces contributions in the meson-baryon coupling, we calculate the EoS of asymmetric hyperonic matter and apply it to describe hyperonic matter and hyperon stars.

Our main goal in this work is to study magnetized neutron stars by using a fully general$-$relativity approach presented in the LORENE package\footnote{http://www.lorene.obspm.fr}. Here we have adopted a non-uniform magnetic field profile which depends on the baryon density. This profile has been used in many previous works and seems to be a good choice to explore maximum effects of the internal magnetic field in these objects. Equally important, the magnetic field treated here is poloidal and axisymmetric. The preliminary results show that stars endowed with a strong magnetic field will be deformed and the mass somewhat increased.

Neutron stars feature extremely high magnetic fields with deduced field strengths of $10^{15}$ G in the case of magnetars and potentially much higher values inside of the star. In this context we consider the appearance of $\rho^-$ meson condensation taking into account the effect of the magnetic field. The results show that, depending on parameters, such a condensation in magnetized neutron stars might (just) occur.

In this work we introduce an extended version of the formalism proposed originally by Taurines et al. that considers the effects of many-body forces simulated by non-linear self-couplings and meson-meson interaction contributions. In this extended version of the model, we assume that matter is at zero temperature, charge neutral and in beta-equilibrium, considering that the baryon octet interacts by the exchange of scalar-isoscalar ($\sigma$,$\,\sigma^*$), vector-isoscalar ($\omega$,$\,\phi$), vector-isovector ($\varrho$) and scalar-isovector ($\delta$) meson fields. Using nuclear matter properties, we constrain the parameters of the model that describe the intensity of the indirectly density dependent baryon-meson couplings to a small range of possible values. We then investigate asymmetric hyperonic matter properties. We report that the formalism developed in this work is in agreement with experimental data and also allows for the existence of massive hyperon stars (with more than $2M_{\odot}$) with small radii, compatible with astrophysical observations.

The possibility of meson condensation in stars and in heavy-ion collisions has been discussed in the past. Here, we study whether rho meson condensation ($\rho^{-}$) can occur in very dense matter and determine the effect of strong magnetic fields on this condensation. We find that rho meson condensates can appear in the core of the neutron star assuming a rho mass which is reduced due to in-medium effects. We find that the magnetic field has a non-negligible effect in triggering condensation.

Magnetars are compact stars which are observationally determined to have very strong surface magnetic fields of the order of $10^{14}-10^{15}$G. The centre of the star can potentially have a magnetic field several orders of magnitude larger. We study the effect of the field on the mass and shape of such a star. In general, we assume a non-uniform magnetic field inside the star which varies with density. The magnetic energy and pressure as well as the metric are expanded as multipoles in spherical harmonics up to the quadrupole term. Solving the Einstein equations for the gravitational potential, one obtains the correction terms as functions of the magnetic field. Using a nonlinear model for the hadronic EoS the excess mass and change in equatorial radius of the star due to the magnetic field are quite significant if the surface field is $10^{15}$G and the central field is about $10^{18}$ G. For a value of the central magnetic field strength of $1.75\times10^{18}$ G, we find that both the excess mass and the equatorial radius of the star changes by about $3-4\%$ compared to the spherical solution.

The recent measurement of two solar mass pulsars has initiated an intense discussion on its impact on our understanding of the high-density matter in the cores of neutron stars. A task force meeting was held from October 7-10, 2013 at the Frankfurt Institute for Advanced Studies to address the presence of quark matter in these massive stars. During this meeting, the recent oservational astrophysical data and heavy-ion data was reviewed. The possibility of pure quark stars, hybrid stars and the nature of the QCD phase transition were discussed and their observational signals delineated.

Strong interaction physics under extreme conditions of high temperature and/or density is of central interest in modern nuclear physics for experimentalists and theorists alike. In order to investigate such systems, model approaches that include hadrons and quarks in a unified approach, will be discussed. Special attention will be given to high-density matter as it occurs in neutron stars. Given the current observational limits for neutron star masses, the properties of hyperonic and hybrid stars will be determined. In this context especially the question of the extent, to which exotic particles like hyperons and quarks affect star masses, will be discussed.

Shock waves constitute discontinuities in matter which are relevant in studying the plasma behaviour in astrophysical scenarios and in heavy-ion collision. They can produce conical emission in relativistic collisions and are also thought to be the mechanism behind the acceleration of energetic particles in active galactic nuclei and gamma ray bursts. The shocks are mostly hydrodynamic shocks. In a magnetic background they become magnetohydrodynamic (MHD) shocks. For that reason we study the space-like and time-like shock discontinuity in a magnetic plasma. The shocks induce a phase transition in the plasma, here assuming a transition from hadron to quarks. The MHD conservation conditions are derived across the shock. The conservation conditions are solved for downstream velocities and flow angles for given upstream variables. The shock conditions are solved at different baryon densities. For the space-like shocks the anisotropy in the downstream velocity arises due to the magnetic field. The downstream velocity vector always points downward with respect to the shock normal. With the increase in density the anisotropy is somewhat reduced. The magnetic field has effectively no effect on time-like shocks. The slight anisotropy in the downstream flow velocities is caused by the boosting that brings the quantities from the fluid frame to normal incidence (NI) frame.

In this work we consider the thermal evolution of rigidly rotating neutron stars. In order to perform such study we first calculate the structure of rotating objects, which is considerably more complicated than that of spherical objects. The structure of rotating neutron stars is obtained by solving Einstein's equation for a rotationally deformed fluid distributions. The numerical method used is based on the the KEH. The equation of state used for computing the neutron star structure and composition is a simple relativistic mean field model, with parameter set G300. With the structure of rotating neutron stars computed we calculate the thermal evolution of these objects. In order to do so, we re-derive the thermal evolution equations to account for the metric of a rotating object. The cooling of neutron stars with different frequencies is then calculated. We show that the cooling of the star strongly depends on the frequency of the object, with higher frequencies stars showing a substantial temperature difference between the equator and poles.

The study of neutron stars, or more general compact stars, is a topic of central interest in nuclear astrophysics. Furthermore, neutron stars serve as the only physical systems whose properties can be used to infer information on cold and dense matter at several times nuclear saturation density. Therefore, neutron star physics is ideally suited to complement the studies of ultra-relativistic heavy-ion collisions that sample strongly interacting matter at high temperature and relatively small net baryon density. In general, in order to pin down or at least constrain the properties of dense matter, accurate measurements of neutron star properties like masses, radii, rotational frequency, and cooling behavior are needed. Here, in relatively recent times the reliable mass determination of the pulsar PSR J1614-2230 of $M = 1.97 \pm 0.04 M_\odot$ has introduced an important benchmark for modeling stars and strongly interacting matter. It puts constraints on the structure of compact stars and possible exotic phases in the core of the stars as will be discussed in this article. In order to investigate this point we will consider a model for star matter that includes hyperonic and quark degrees of freedom, and present results for compact star properties in the following.

First order phase transitions (PTs) with more than one globally conserved charge, so-called non-congruent PTs, have characteristic differences compared to congruent PTs (e.g., dimensionality of phase diagrams, location and properties of critical points and endpoints). In the present article we investigate the non-congruence of the nuclear liquid-gas PT at sub-saturation densities and the deconfinement PT at high densities and/or temperatures in Coulomb-less models, relevant for heavy-ion collisions and neutron stars. For the first PT, we use the FSUgold relativistic mean-field model and for the second one the relativistic chiral SU(3) model. The chiral SU(3) model is one of the few models for the deconfinement PT, which contains quarks and hadrons in arbitrary proportions (i.e. a "solution") and gives a continuous transition from pure hadronic to pure quark matter above a critical point. The study shows the universality of the applied concept of non-congruence for the two PTs with an upper critical point, and illustrates the different typical scales involved. In addition, we find a principle difference between the liquid-gas and the deconfinement PTs: in contrast to the ordinary Van-der-Waals-like PT, the phase coexistence line of the deconfinement PT has a negative slope in the pressure-temperature plane. As another qualitative difference we find that the non-congruent features of the deconfinement PT become vanishingly small around the critical point.

The unexpected temperature evolution of the neutron star in the Cassiopeia A supernova remnant (Cas A, for short) has renewed tremendous interest in the cooling mechanisms of neutron stars. In particular, the formation of superconducting protons and superfluid neutrons deep inside the cores of neutron stars have become focal points of the discussion. The purpose of this letter is to add a new aspect to this discussion, which focuses on the connection between proton superconductivity and the masses of neutron stars. Assuming (as is currently the case) that the temperature evolution of Cas A is largely controlled by superconducting protons, we study a series of phenomenological proton-pairing models to determine how deep into the stellar core superconducting protons actually penetrate. This allows us to establish a heretofore unknown relationship between the mass of the neutron star in Cas A and the penetration depth of the superconducting proton phase. This relationship can be used to either predict the depth of the superconducting proton phase, or, conversely, determine the mass of Cas A from a reliable calculation of the size of the proton superconducting phase in superdense neutron star matter. We emphasize that the strategy outlined in this paper can be applied to any other neutron star of similar age, whose temperature might be reliably monitored over a several years period. High-mass neutron stars, such as the recently discovered neutron stars J1614-2230 ($1.97 \pm 0.04\, \msun$) and J0348+0432 ($2.01 \pm 0.04 \, \msun$), appear particularly appealing as a significant fraction of the protons in their cores may be superconducting.

We model magnetars as hybrid stars, which have a core of quark matter surrounded by hadronic matter. For this purpose, we use an extended version of the SU(3) non-linear realization of the sigma model in which the degrees of freedom change naturally from hadrons to quarks as the temperature/density increases. The presence of a variable magnetic field allows us to study in detail the influence of Landau quantization and the anomalous magnetic moment on the particle population of the star, more precisely on particles with different spin projections. This allows us to calculate the polarization of the system throughout different phases of the star, hadronic, quark and also a mixed phase.

We apply an extended version of the SU(3) parity model, containing quark degrees of freedom, to study neutron stars. The model successfully reproduces the main thermodynamic features of QCD which allows us to describe the composition of dense matter. Chiral symmetry restoration is realized inside the star and the chiral partners of the baryons appear, their masses becoming degenerate. Furthermore, quark degrees of freedom appear in a transition to a deconfined state. Performing an investigation of the macroscopic properties of neutron stars, we show that observational constraints, like mass and thermal evolution, are satisfied and new predictions can be made.

Driven by the loss of energy, isolated rotating neutron stars (pulsars) are gradually slowing down to lower frequencies, which increases the tremendous compression of the matter inside of them. This increase in compression changes both the global properties of rotating neutron stars as well as their hadronic core compositions. Both effects may register themselves observationally in the thermal evolution of such stars, as demonstrated in this Letter. The rotation-driven particle process which we consider here is the direct Urca (DU) process, which is known to become operative in neutron stars if the number of protons in the stellar core exceeds a critical limit of around 11% to 15%. We find that neutron stars spinning down from moderately high rotation rates of a few hundred Hertz may be creating just the right conditions where the DU process becomes operative, leading to an observable effect (enhanced cooling) in the temperature evolution of such neutron stars. As it turns out, the rotation-driven DU process could explain the unusual temperature evolution observed for the neutron star in Cas A, provided the mass of this neutron star lies in the range of 1.5 to 1.9 \msun and its rotational frequency at birth was between 40 (400 Hz) and 70% (800 Hz) of the Kepler (mass shedding) frequency, respectively.

We study the effects of high magnetic fields on the particle population and equation of state of hybrid stars using an extended hadronic and quark SU(3) non-linear realization of the sigma model. In this model the degrees of freedom change naturally from hadrons to quarks as the density and/or temperature increases. The effects of high magnetic fields and anomalous magnetic moment are visible in the macroscopic properties of the star, such as mass, adiabatic index, moment of inertia, and cooling curves. Moreover, at the same time that the magnetic fields become high enough to modify those properties, they make the star anisotropic.

We study the so called hybrid stars, which are hadronic stars that contain a core of deconfined quarks. For this purpose, we make use of an extended version of the SU(3) chiral model. Within this approach, the degrees of freedom change naturally from hadrons (baryon octet) to quarks (u, d, s) as the temperature and/or density increases. At zero temperature we are still able to reproduce massive stars, even with the inclusion of hyperons.

In this paper we investigate the thermal evolution of hybrid stars, objects composed of a quark matter core, enveloped by ordinary hadronic matter. Our purpose is to investigate how important are the microscopic properties of the quark core to the thermal evolution of the star. In order to do that we use a simple MIT bag model for the quark core, and a relativistic mean field model for the hadronic envelope. By choosing different values for the microscopic parameters (bag constant, strange quark mass, strong coupling constant) we obtain hybrid stars with different quark core properties. We also consider the possibility of color superconductivity in the quark core. With this simple approach, we have found a set of microscopic parameters that lead to a good agreement with observed cooling neutron stars. Our results can be used to obtain clues regarding the properties of the quark core in hybrid stars, and can be used to refine more sophisticated models for the equation of state of quark matter.

There are many factors that contribute to the breaking of the spherical symmetry of a neutron star. Most notably is rotation, magnetic fields, and/or accretion of matter from companion stars. All these phenomena influence the macroscopic structures of neutron stars, but also impact their microscopic compositions. The purpose of this paper is to investigate the cooling of rotationally deformed, two-dimensional (2D) neutron stars in the framework of general relativity theory, with the ultimate goal of better understand the impact of 2D effects on the thermal evolution of such objects. The equations that govern the thermal evolution of rotating neutron stars are presented in this paper. The cooling of neutron stars with different frequencies is computed self-consistently by combining a fully general relativistic 2D rotation code with a general relativistic 2D cooling code. We show that rotation can significantly influence the thermal evolution of rotating neutron stars. Among the major new aspects are the appearances of hot spots on the poles, and an increase of the thermal coupling times between the core and the crust of rotating neutron stars. We show that this increase is independent of the microscopic properties of the stellar core, but depends only on the frequency of the star.

The study of neutron stars is a topic of central interest in the investigation of the properties of strongly compressed hadronic matter. Whereas in heavy-ion collisions the fireball, created in the collision zone, contains very hot matter, with varying density depending on the beam energy, neutron stars largely sample the region of cold and dense matter with the exception of the very short time period of the existence of the proto-neutron star. Therefore, neutron star physics, in addition to its general importance in astrophysics, is a crucial complement to heavy-ion physics in the study of strongly interacting matter. In the following, model approaches will be introduced to calculate properties of neutron stars that incorporate baryons and quarks. These approaches are also able to describe the state of matter over a wide range of temperatures and densities, which is essential if one wants to connect and correlate star observables and results from heavy-ion collisions. The effect of exotic particles and quark cores on neutron star properties will be considered. In addition to the gross properties of the stars like their masses and radii their expected inner composition is quite sensitive to the models used. The effect of the composition can be studied through the analysis of the cooling curve of the star. In addition, we consider the effect of rotation, as in this case the particle composition of the star can be modified quite drastically.

We discuss the properties of neutron stars and their modifications due to the occurrence of hyperons and quarks in the core of the star. More specifically, we consider the general problem of exotic particles inside compact stars in light of the observed two-solar mass pulsar. In addition, we investigate neutron star cooling and a possible explanation of the recently measured cooling curve of the neutron star in the supernova remnant Cas A.

We investigate the properties of dense matter and neutron stars. In particular we discuss model calculations based on the parity doublet picture of hadronic chiral symmetry. In this ansatz the onset of chiral symmetry restoration is reflected by the degeneracy of baryons and their parity partners. In this approach we also incorporate quarks as degrees of freedom to be able to study hybrid stars.

We develop a unified model of hadrons and quarks. Within this approach we investigate the phase structure of the model as function of temperature and chemical potential. Computing the equation of state of cold matter we determine neutron and hybrid star masses and radii. In an extension of the investigation we consider the cooling behavior of the compact stars and derive a general relation between the star's mass and rotation and its cooling behavior. Finally we study the effect of Delta resonances for star matter, especially with respect to possible solutions of stars with small radii.

In this paper we investigate the electrodynamic surface properties of bare strange quark stars. The surfaces of such objects are characterized by the formation of ultra-high electric surface fields which might be as high as $\sim 10^{19}$ V/cm. These fields result from the formation of electric dipole layers at the stellar surfaces. We calculate the increase in gravitational mass associated with the energy stored in the electric dipole field, which turns out to be only significant if the star possesses a sufficiently strong {\em net} electric charge distribution. In part two of the paper, we explore the intriguing possibility of what happens when the electron layer (sphere) rotates with respect to the stellar strange matter body. We find that in this event magnetic fields can be generated which, for moderate effective rotational frequencies between the electron layer and the stellar body, agree with the magnetic fields inferred for several Central Compact Objects (CCOs). These objects could thus be comfortably interpreted as strange stars whose electron atmospheres rotate at frequencies that are moderately different ($\sim 10$ Hz) from the rotational frequencies of the strange star itself.

We study the behavior of hybrid stars using an extended hadronic and quark SU(3) non-linear sigma model. The degrees of freedom change naturally, in this model, from hadrons to quarks as the density/temperature increases. At zero temperature, we reproduce massive neutron stars containing a core of hybrid matter of 2 km for the non-rotating case and 1.18 km and 0.87 km, in the equatorial and polar directions respectively, for stars rotating at the Kepler frequency (physical cases lie in between). The cooling of such stars is also analyzed.

Recent neutron star observations suggest that the masses and radii of neutron stars may be smaller than previously considered, which would disfavor a purely nucleonic equation of state. In our model, we use a the flavor SU(3) sigma model that includes delta resonances and hyperons in the equation of state. We find that if the coupling of the delta resonances to the vector mesons is slightly smaller than that of the nucleons, we can reproduce both the measured mass-radius relationship and the extrapolated equation of state.

We describe an extension of the hadronic SU(3) non-linear sigma model to include quarks. As a result, we obtain an effective model which interpolates between hadronic and quark degrees of freedom. The new parameters and the potential for the Polyakov loop (used as the order parameter for deconfinement) are calibrated in order to fit lattice QCD data and reproduce the QCD phase diagram. Finally, the equation of state provided by the model, combined with gravity through the inclusion of general relativity, is used to make predictions for neutron stars.

We extend the hadronic SU(3) non-linear sigma model to include quark degrees of freedom. The choice of potential for the Polyakov loop as a function of temperature and chemical potential allows us to construct a realistic phase diagram from the analysis of the order parameters of the system. These parameters are the chiral condensate, for the chiral symmetry restoration and the Polyakov loop, for the deconfinement to quark matter. Besides reproducing lattice QCD results, for zero and low chemical potential, we are in agreement with neutron star observations for zero temperature.

We investigate strongly interacting dense matter and neutron stars using a flavor-SU(3) approach based on a non-linear realization of chiral symmetry as well as a hadronic flavor-SU(2) parity-doublet model. We study chiral symmetry restoration and the equation of state of stellar matter and determine neutron star properties using different sets of degrees of freedom. Finally, we include quarks in the model approach. We show the resulting phase diagram as well as hybrid star solutions for this model.

We study different stages of the neutron star cooling by computing neutron star properties at various temperatures and entropies using an effective chiral model including hadronic and quark degrees of freedom. Macroscopic properties of the star such as mass and radius are calculated and compared with observations. It can be seen that the effects of chiral restoration and deconfinement to quark matter in the core of the neutron star at different stages of the evolution can be significant for the evolution of the star and allow insight into the behaviour of matter at extreme densities.

A hadronic chiral SU(3) model is applied to neutron and proto-neutron stars, taking into account trapped neutrinos, finite temperature and entropy. The transition to the chirally restored phase is studied and global properties of the stars like minimum and maximum masses and radii are calculated for different cases. In addition, the effects of rotation on neutron star masses are included and the conservation of baryon number and angular momentum determine the maximum frequencies of rotation during the cooling.

Neutron star matter is investigated in a hadronic chiral model approach using the lowest flavor-SU(3) multiplets for baryons and mesons. The parameters are determined to yield consistent results for saturated nuclear matter as well as for finite nuclei. The influence of baryonic resonances is discussed. The global properties of a neutron star such as its mass and radius are determined. Proto-neutron star properties are studied by taking into account trapped neutrinos, temperature and entropy effects.

The parity doublet model, containing the SU(2) multiplets including the baryons identified as the chiral partners of the nucleons is applied to neutron stars. The maximum mass for the star is calculated for different stages of the cooling taking into account finite temperature/entropy effect, trapped neutrinos and fixed baryon number. Rotation effects are also included.

In these lecture notes, the role of strangeness in relativistic astrophysics of compact stars is addressed. The appearance of strange particles, as hyperons, kaons, and strange quarks, in the core of compact stars is examined and common features as well as differences are presented. Impacts on the global properties of compact stars and signals of the presence of exotic matter are outlined for the various strange phases which can appear in the interior at high densities.

The Parity doublet model containing the SU(2) multiplets including the baryons identified as the chiral partners of the nucleons is applied for neutron star matter. The chiral restoration is analyzed and the maximum mass of the star is calculated.

We investigate various properties of neutron star matter within an effective chiral $SU(3)_L \times SU(3)_R$ model. The predictions of this model are compared with a Walecka-type model. It is demonstrated that the importance of hyperon degrees are strongly depending on the interaction used, even if the equation of state near saturation density is nearly the same in both models. While the Walecka-type model predicts a strange star core with strangeness fraction $f_S \approx 4/3$, the chiral model allows only for $f_S \approx 1/3$ and predicts that $\Sigma^0$, $\Sigma^+$ and $\Xi^0$ will not exist in star, in contrast to the Walecka-type model.

Neutrino-nucleus cross sections for the detection of atmospheric neutrinos are calculated in relativistic impulse and random phase approximations. Pion production is estimated via inclusive Delta-hole nuclear excitations. The number of pion events can confuse the identification of the neutrino flavor. As a further source of pions we calculate coherent pion production (where the nucleus remains in the ground state). We examine how these nuclear structure effects influence atmospheric neutrino experiments and study possible improvements in the present detector simulations.