Search

The inner and outer cores of neutron stars are believed to contain type-I and -II proton superconductors, respectively. The type-I superconductor exists in an intermediate state, comprising macroscopic flux-free and flux-containing regions, while the type-II superconductor is flux-free, except for microscopic, quantized flux tubes. Here, we show that the inner and outer cores are coupled magnetically, when the macroscopic flux tubes subdivide dendritically into quantized flux tubes, a phenomenon called flux branching. An important implication is that up to |${\sim} 10^{12} (r_1/10^6 \, {\rm cm}) \, {\rm erg}$| of energy are required to separate a quantized flux tube from its progenitor macroscopic flux tube, where |$r_1$| is the length of the macroscopic flux tube. Approximating the normal-superconducting boundary as sharp, we calculate the magnetic coupling energy between a quantized and macroscopic flux tube due to flux branching as a function of, |$f_1$|⁠ , the radius of the type-I inner core divided by the radius of the type-II outer core. Strong coupling delays magnetic field decay in the type-II superconductor. For an idealized inner core containing only a type-I proton superconductor and poloidal flux, and in the absence of ambipolar diffusion and diamagnetic screening, the low magnetic moments (⁠|${\lesssim} 10^{27} \, {\rm G \, cm^3}$|⁠) of recycled pulsars imply |$f_1 \lesssim 10^{-1.5}$|⁠. [ABSTRACT FROM AUTHOR]

Millisecond magnetars are widely suggested as the central engines powering hydrogen-poor superluminous supernovae (SLSNe). These magnetars primarily lose huge rotational energy through gravitational wave radiation (GWR) and magnetic dipole radiation (MDR), with MDR serving as an energy source for SLSNe. We study the evolution of the magnetar spin, magnetic inclination angle, and the resulting thermal radiative luminosity of the SLSNe, where the impacts of the nucleon–nucleon short-range correlation, the mass and initial spin of the magnetar, and the density-dependent symmetry energy of the dense nuclear matter on the evolution are discussed. The relativistic mean-field theory is employed to calculate the nuclear matter properties, and we particularly concentrate on the time- and space-dependent bulk viscosity, which is crucial for the magnetic inclination angle evolution. It is found that the nucleon–nucleon short-range correlation weakens the damping of bulk viscosity of dense matter and therefore inhibits the growth of magnetic inclination angle, and it reduces the MDR (GWR) peak luminosity of a canonical magnetar by several times while it raises the peak thermal radiation luminosity of SLSNe by several times. For magnetars with non-rotating mass obviously lower than the |$1.4 \, \rm M_\odot$| with slow initial rotation, the magnetic inclination angle is more likely to evolve towards 0 deg quickly, and these magnetars are not suitable as the central engine for SLSNe. Within the 'family' of FSUGarnet interaction, a stiffer symmetry energy gives a lower threshold of direct Urca process and hence gives a much larger bulk viscosity coefficient, and thus it promotes the growth of the magnetic inclination angle and the GWR for canonical stars but reduces the peak brightness of SLSNe significantly. [ABSTRACT FROM AUTHOR]

The accumulation of accreted matter onto the neutron star surface triggers exothermic reactions in the crust. The heat released as a result influences the luminosity exhibited by the X-ray transient. The most common approach to the kinetics of exothermic reactions in the crust of accreting neutron stars is to consider an infinite reaction rate. Here, we investigate accretion-related heat release in the accreted outer crust of a neutron star by including a time-dependent accretion cycle and experimentally based reaction rates in the kinetics of electron captures above the reaction threshold. A simple model was used to compute the zero temperature equation of state of a crust in which two nuclei can coexist. We solved the abundance of parent nuclei as a function of the depth in the star and the time variable using astrophysically motivated features of the accreting system. We calculated the heat release and neutrino loss associated to reactions in the outer crust. We report the existence of layers in the outer crust, which contain both parent and grand-daughter nuclei of electron captures. The reactions can occur deeper in the shell than the reaction threshold, thus releasing more heat per accreted baryon for a given accretion rate. The electron capture layers continue to exist even when the accretion has stopped. The heat sources are time- and pressure-dependent in accreting crusts of neutron stars. The total heat released is a function of astrophysical (active and quiescent time) and microscopic (reaction rate) parameters Therefore, we conclude these parameters should be considered individually and carefully for a range of different sources. [ABSTRACT FROM AUTHOR]

According to the hypothesis that strange quark matter may be the true ground state of matter at extremely high densities, strange quark stars should be stable and could exist in the Universe. It is possible that pulsars may actually be strange stars, but not neutron stars. Here we present a short review on recent progresses in the field of strange quark stars. First, three popular phenomenological models widely used to describe strange quark matter are introduced, with special attention being paid on the corresponding equation of state in each model. Combining the equation of state with the Tolman-Oppenheimer-Volkov equations, the inner structure and mass-radius relation can be obtained for the whole sequence of strange stars. Tidal deformability and oscillations (both radial and non-radial oscillations), which are sensitive to the composition and the equations of state, are then described. Hybrid stars as a special kind of quark stars are discussed. Several other interesting aspects of strange stars are also included. For example, strong gravitational wave emissions may be generated by strange stars through various mechanisms, which may help identify strange stars via observations. Especially, close-in strange quark planets with respect to their hosts may provide a unique test for the existence of strange quark objects. Fierce electromagnetic bursts could also be generated by strange stars. The energy may come from the phase transition of neutron stars to strange stars, or from the merger of binary strange stars. The collapse of the strange star crust can also release a huge amount of energy. It is shown that strange quark stars may be involved in short gamma-ray bursts and fast radio bursts. [ABSTRACT FROM AUTHOR]

Abstract This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided.

With current and planned gravitational-wave (GW) observing runs, coincident multimessenger timing of resonant shattering flares (RSFs) and GWs may soon allow for neutron star (NS) asteroseismology to be used to constrain the nuclear symmetry energy, an important property of fundamental nuclear physics that influences the composition and equation of state of NSs. In this work, we examine the effects of combining multiple RSF detections on these symmetry energy constraints, and consider how realistic uncertainties in the masses of the progenitor NSs may weaken them. We show that the detection of subsequent multimessenger events has the potential to substantially improve constraints beyond those obtained from the first, and that this improvement is insensitive to the mass of the NSs that produce the RSFs and its uncertainty. This sets these asteroseismic constraints apart from bulk NS properties such as radius, for which the NS mass is highly important, meaning that any multimessenger RSF and GW events can equally improve our knowledge of fundamental physics. [ABSTRACT FROM AUTHOR]

We conducted a study on the thermal properties of stellar matter with the nuclear energy density functional BCPM. This functional is based on microscopic Brueckner–Hartree–Fock calculations and has demonstrated success in describing cold neutron stars. To enhance its applicability in astrophysics, we extended the BCPM equation of state to finite temperature for β-stable neutrino-free matter, taking into consideration the hot inner crust. Such an equation of state holds significant importance for hot compact objects, particularly those resulting from a binary neutron star merger event. Our exploration has shown that with increasing temperature, there is a fast decrease in the crust-core transition density, suggesting that for hot stars it is not realistic to assume a fixed value of this density. The microscopic calculations also reveal that the presence of nuclear clusters persists up to T = 7.21 MeV, identified as the limiting temperature of the crust. Above this threshold, the manifestation of clusters is not anticipated. Below this temperature, clusters within the inner crust are surrounded by uniform matter with varying densities, allowing for the distinction between the upper and lower transition density branches. Moreover, we computed mass–radius relations of neutron stars, assuming an isothermal profile for β-stable neutron star matter at various temperature values. Our findings highlight the significant influence of the hot inner crust on the mass–radius relationship, leading to the formation of larger and more inflated neutron stars. Consequently, under our prescription, the final outcome is a unified equation of state at finite temperature. [ABSTRACT FROM AUTHOR]

Neutron stars (NS) that are born in binary systems with a main-sequence star companion can experience mass transfer, resulting in the accumulation of material at the surface of the NS. This, in turn, leads to the continuous growth of the NS mass and the associated steepening of the gravitational potential. Supposing the central density surpasses the onset for the phase transition from nuclear, generally hadronic matter to deconfined quark-gluon plasma, which is a quantity currently constrained solely from an upper limit by asymptotic freedom in quantum chromodynamics (QCD), the system may experience a dynamic response due to the appearance of additional degrees of freedom in the equation of state (EOS). This dynamical response might give rise to a rapid softening of the EOS during the transition in the hadron-quark matter co-existence region. While this phenomenon has long been studied in the context of hydrostatic configurations, the dynamical implications of this problem are still incompletely understood. It is the purpose of the present paper to simulate the dynamics of NSs with previously accreted envelopes caused by the presence of a first-order QCD phase transition. Therefore, we employed the neutrino radiation hydrodynamics treatment based on the fully general relativistic approach in spherical symmetry, implementing a three-flavor Boltzmann neutrino transport and a microscopic model EOS that contains a first-order hadron-quark phase transition. The associated neutrino signal shows a sudden rise in the neutrino fluxes and average energies, becoming observable for the present generation of neutrino detectors for a galactic event, and a gravitational wave mode analysis revealed the behaviors of the dominant f mode and the first and the second gravity g modes that are excited during the NS evolution across the QCD phase transition. [ABSTRACT FROM AUTHOR]

Context. Relating different global neutron-star (NS) properties, such as tidal deformability and radius, or mass and radius, requires an equation of state (EoS). Determining the NS EoS is therefore not only the science goal of a variety of observational projects, but it also enters in the analysis process; for example, to predict a NS radius from a measured tidal deformability via gravitational waves (GW) during the inspiral of a binary NS merger. To this aim, it is important to estimate the theoretical uncertainties on the EoS, one of which is the possible bias coming from an inconsistent treatment of the low-density region; that is, the use of a so called non-unified NS crust. Aims. We propose a numerical tool allowing the user to consistently match a nuclear-physics informed crust to an arbitrary high-density EoS describing the core of the star. Methods. We introduce an inversion procedure of the EoS close to saturation density that allows users to extract nuclear-matter parameters and extend the EoS to lower densities in a consistent way. For the treatment of inhomogeneous matter in the crust, a standard approach based on the compressible liquid-drop (CLD) model approach was used in our work. A Bayesian analysis using a parametric agnostic EoS representation in the high-density region is also presented in order to quantify the uncertainties induced by an inconsistent treatment of the crust. Results. We show that the use of a fixed, realistic-but-inconsistent model for the crust causes small but avoidable errors in the estimation of global NS properties and leads to an underestimation of the uncertainties in the inference of NS properties. Conclusions. Our results highlight the importance of employing a consistent EoS in inference schemes. The numerical tool that we developed to reconstruct such a thermodynamically consistent EoS, CUTER, has been tested and validated for use by the astrophysical community. [ABSTRACT FROM AUTHOR]

Crystallization of dense matter in neutron star crusts and white dwarf cores may be similar to epitaxial crystal growth in terrestrial laboratories. However, in stellar crystals, the spacing between horizontal planes has to gradually increase with the outward movement of the crystallization front, tracing decrease of the electron density. This process produces Coulomb crystals with stretched rather than cubic elementary cells. We extend the analysis of the elastic and breaking properties of such crystals to the face‐centered (fc) lattice. Shear deformations orthogonal to the stretch direction have been studied for 22 crystallographic shear planes. A common property for all these planes is the reduction and eventual nulling of the breaking shear strain with deviation from the unstretched configuration. The effective shear moduli for deformations orthogonal to the stretch direction have been calculated. It is possible that the epitaxial crystallization in compact stars results in a formation of large‐scale crystallites or, at least, in growth of the whole crystallization front perpendicular to particular crystallographic planes. For fc structure growth orthogonal to the {111}$$ \left\{111\right\} $$ planes, we expect that, at any density, ∼5%$$ \sim 5\% $$ (∼0.5%$$ \sim 0.5\% $$) of crystallite height is occupied by layers one (two) orders of magnitude weaker than the bulk of the crystallite. This may be important for realistic modeling of crustquakes on neutron stars. [ABSTRACT FROM AUTHOR]

This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided. [ABSTRACT FROM AUTHOR]

We study cold strange quark stars employing an enhanced version of the quark-mass density-dependent model, which incorporates excluded volume effects to address non-perturbative QCD repulsive interactions. We provide a comparative analysis of our mass formula parametrization with previous models from the literature. We identify the regions within the parameter space where three-flavor quark matter is more stable than the most tightly bound atomic nucleus (stability window). Specifically, we show that excluded volume effects do not change the Gibbs free energy per baryon at zero pressure, rendering the stability window unaffected. The curves of pressure versus energy density exhibit various shapes—convex upward, concave downward, or nearly linear—depending on the mass parametrization. This behavior results in different patterns of increase, decrease, or constancy in the speed of sound as a function of baryon number density. We analyze the mass–radius relationship of strange quark stars, revealing a significant increase in maximum gravitational mass and a shift in the curves toward larger radii as the excluded volume effect intensifies. Excluded volume effects render our models compatible with all modern astrophysical constraints, including the properties of the recently observed low-mass compact object HESSJ1731. [ABSTRACT FROM AUTHOR]

The massive stars end their lives by supernova explosions leaving central compact objects that may evolve into neutron stars. Initially, after birth, the star remains hot and gradually cools down. We explore the matter and star properties during this initial stage of the compact stars considering the possibility of the appearance of deconfined quark matter in the core of the star. At the initial stage after the supernova explosion, the occurrence of non-radial oscillation in the newly born compact object is highly possible. Non-radial oscillations are an important source of gravitational waves (GWs). There is a high chance for GWs from these oscillations, especially the nodeless fundamental (f) mode to be detected by next-generation GW detectors. We study the evolution in frequencies of non-radial oscillation after birth considering phase transition and predicting the possible signature for different possibilities of theoretical compact star models. [ABSTRACT FROM AUTHOR]

Measurements of neutron star mass and radius or tidal deformability deliver unique insight into the equation of state (EOS) of cold dense matter. EOS inference is very often done using generalized parametric or non-parametric models, which deliver no information on composition. In this paper, we consider a microscopic nuclear EOS model based on a field theoretical approach. We show that current measurements from NICER and gravitational wave observations constrain primarily the symmetric nuclear matter EOS. We then explore what could be delivered by measurements of mass and radius at the level anticipated for future large-area X-ray timing telescopes. These should be able to place very strong limits on the symmetric nuclear matter EOS, in addition to constraining the nuclear symmetry energy that determines the proton fraction inside the neutron star. [ABSTRACT FROM AUTHOR]

According to the hypothesis that strange quark matter may be the true ground state of matter at extremely high densities, strange quark stars should be stable and could exist in the Universe. It is possible that pulsars may actually be strange stars, but not neutron stars. Here we present a short review on recent progresses in the field of strange quark stars. First, three popular phenomenological models widely used to describe strange quark matter are introduced, with special attention being paid on the corresponding equation of state in each model. Combining the equation of state with the Tolman-Oppenheimer-Volkov equations, the inner structure and mass-radius relation can be obtained for the whole sequence of strange stars. Tidal deformability and oscillations (both radial and non-radial oscillations), which are sensitive to the composition and the equations of state, are then described. Hybrid stars as a special kind of quark stars are discussed. Several other interesting aspects of strange stars are also included. For example, strong gravitational wave emissions may be generated by strange stars through various mechanisms, which may help identify strange stars via observations. Especially, close-in strange quark planets with respect to their hosts may provide a unique test for the existence of strange quark objects. Fierce electromagnetic bursts could also be generated by strange stars. The energy may come from the phase transition of neutron stars to strange stars, or from the merger of binary strange stars. The collapse of the strange star crust can also release a huge amount of energy. It is shown that strange quark stars may be involved in short gamma-ray bursts and fast radio bursts.

We present microscopic molecular dynamics simulations including the efficient Ewald sum procedure to study warm and dense stellar plasmas consisting of finite-size ion charges immerse in a relativistic neutralizing electron gas. For densities typical of Supernova matter and crust in a proto-neutron star, we select a representative single ion composition and obtain the virialized equation of state (vEoS). We scrutinize the finite-size and screening corrections to the Coulomb potential appearing in the virial coefficients B 2, B 3, and B 4 as a function of temperature. In addition, we study the thermal heat capacity at constant volume, CV , and the generalized Mayer's relation i.e. the difference CP − CV with CP being the heat capacity at constant pressure, obtaining clear features signalling the onset of the liquid–gas phase transition. Our findings show that microscopic simulations reproduce the discontinuity in CV , whose value lies between that of idealized gas and crystallized configurations. We study the pressure isotherms marking the boundary of the metastable region before the gaseous transition takes place. The resulting vEoS displays a behaviour where effective virial coefficients include extra density dependence showing a generalized density–temperature form. As an application we parametrize pressure as a function of density and temperature under the form of an artificial neural network showing the potential of machine learning for future regression analysis in more refined multicomponent approaches. This is of interest to size the importance of these corrections in the liquid–gas phase transition in warm and dense plasma phases contributing to the cooling behaviour of early Supernova phases and proto-neutron stars. [ABSTRACT FROM AUTHOR]

In this work, we delve into the temperature-dependent Equation of State (EoS) of baryonic matter within the framework of the FSU2H* hadronic model, which comprehensively incorporates hyperons and is suitable for relativistic simulations of neutron star mergers and supernovae. To assess the impact of the uncertainties in the hyperonic sector on astrophysical observables, we introduce two additional models, namely FSU2H*L (FSU2H*-Lower) and FSU2H*U(FSU2H*-Upper). These models cover the entire spectrum of variability of hyperonic potentials, as derived from experimental data. Our investigations reveal that these uncertainties extend their influence not only to the relative abundances of various particle species but also to the EoS itself and, consequently, have an impact on the global properties of both cold and hot neutron stars. Notably, their effects become more pronounced at large temperatures, owing to the increased presence of hyperons. These findings have direct implications for the outcomes of relativistic simulations of neutron star mergers and supernovae, emphasizing the need of accounting for hyperonic uncertainties to ensure the accuracy and reliability of such simulations in astrophysical contexts. [ABSTRACT FROM AUTHOR]

We report elastic and inelastic cross-sections for fast superthermal 12C(3 P) and 13C(3 P) atoms scattering on 12CO2. The cross-sections were computed using quantum-mechanical rotationally close-coupling formalism with the electronic interaction described by a newly constructed potential energy surface correlating to the lowest energy asymptote of the complex. State-to-state rotational cross-sections, differential cross-sections, and derived transport properties of interest for energy relaxation are also reported. The computed elastic differential cross-sections are found to be strongly anisotropic, show significant energy dependence, and differ by up to 2 percent between the two isotopes of carbon. [ABSTRACT FROM AUTHOR]

Near-equilibrium bottom-up crystallization of fully ionized neutron star crusts or white dwarf cores is considered. We argue that this process is similar to liquid-phase epitaxial (i.e. preserving order of previous layers) crystal growth or crystal pulling from melt in Earth laboratories whereby lateral positions of newly crystallizing ions are anchored by already solidified layers. Their vertical positions are set by charge neutrality. Consequently, interplane spacing of a growing crystal either gradually increases, tracing n e decrease, as the crystallization front moves away from the stellar centre, or decreases, tracing decrease of 〈 Z 〉, when the crystallization front crosses a boundary between layers of different compositions. This results in the formation of stretched Coulomb crystals, in contrast to the standard assumption of cubic crystal formation, which is based on energetics arguments but does not take into account growth kinetics. Overstretched crystals break, which limits the vertical sizes of growing crystallites. We study breaking shear strain and effective shear modulus of stretched matter and discuss possibility of macrocrystallite formation. The latter has interesting astrophysical implications, for instance, appearance of weak crustal layers, whose strength may increase by a few orders of magnitude upon breaking and refreezing at a late-time event. We also analyse interaction of adjacent Coulomb crystals, having different ion compositions, and estimate the strength of such interfaces. [ABSTRACT FROM AUTHOR]

We report on long-term evolution of gamma-ray flux and spin-down rate of two bright gamma-ray pulsars, PSR J2021+4026 and Vela (PSR J0835−4510). PSR J2021+4026 shows repeated state changes in gamma-ray flux and spin-down rate. We report two new state changes, a first one from a low gamma-ray flux to a high flux that occurred around MJD 58910, and a second one from high to low flux that occurred around MJD 59510. We find that the flux changes associated with these two state changes are smaller than those determined in the previous events, and the waiting time of the new state change from the high gamma-ray flux to low gamma-ray flux is significantly shorter than previous events. Since the waiting time-scale of the quasi-periodic state changes of PSR J2021+4026 is similar to the waiting time-scale of the glitch events of the Vela pulsar, we search for the state change of the gamma-ray emission of the Vela pulsar to investigate the possibility that the glitching process is the trigger of the state change of PSR J2021+4026. For the Vela pulsar, the flux of the radio pulses briefly decreased around the 2016 glitch, suggesting that the glitch may have affected the structure of the magnetosphere. Nevertheless, we could not find any significant change of the gamma-ray emission properties using 15 yr of Fermi -LAT data. Overall, it seems inconclusive that a glitch-like process similar to that occurred to the Vela pulsar triggers the structure change of the global magnetosphere and causes state changes of PSR J2021+4026. Further and deep investigations to clarify the mechanism of the mode change for PSR J2021+4026 are required. [ABSTRACT FROM AUTHOR]

Quantum Chromodynamics predicts phase transition from hadronic matter to quark matter at high density, which is highly probable in astrophysical systems like binary neutron star mergers. To explore the critical density where such phase transition can occur, we performed numerical relativity simulations of binary neutron star mergers with various masses (equal and unequal binaries). We aim to understand the effect of the onset of phase transition on the merger dynamics and gravitational wave spectra. We generated a set of equations of states by agnostically changing the onset of phase transition, having the hadronic matter part and quark matter part fixed. This particular arrangement of the equation of states explores the scenario of mergers where mixed phases of matter are achieved before or during the merger. Under these circumstances, if the matter properties with hadronic and quark degrees differ significantly, it is reflected in the stability of the final merger product for the intermediate mass binary. We performed a case study on mixed species merger, where one of the binary companions is hybrid star. If quark matter appears at low densities, we observe significant change in post-merger gravitational wave analysis in terms of higher peak frequencies and post-merger frequencies in power spectral density. We report indications expressed as spikes in phase difference plots at merger time for mixed mergers. We found that the expression of phase transition in post-merger gravitational wave signals is more significant for unequal mass binary than for equal mass binary having the same total baryonic mass. [ABSTRACT FROM AUTHOR]

The puzzling mechanism of coherent radio emission remains unknown, but fortunately, repeating fast radio bursts (FRBs) provide a precious opportunity, with extremely bright subpulses created in a clear and vacuum‐like pulsar magnetosphere. FRBs are millisecond‐duration signals that are highly dispersed at distant galaxies but with uncertain physical origin(s). Coherent curvature radiation by bunches has already been proposed for repeating FRBs. The charged particles are created during central star's quakes, which can form bunches streaming out along curved magnetic field lines, so as to trigger FRBs. The nature of narrow‐band radiation with time‐frequency drifting can be a natural consequence that bunches could be observed at different times with different curvatures. Additionally, high linear‐polarization can be seen if the line of sight is confined to the beam angle, whereas the emission could be highly circular‐polarized if off‐beam. It is also discussed that pulsar surface may be full of small hills (i.e., zits) which would help producing bulk of energetic bunches for repeating FRBs as well as for rotation‐powered pulsars. [ABSTRACT FROM AUTHOR]

At the high temperatures inside recently formed neutron stars (⁠|$T\gtrsim 5\times 10^{8}\, \text{K}$|⁠), the particles in their cores are in the 'strong-coupling' regime, in which collisional forces make them behave as a single, stably stratified, and thus non-barotropic fluid. In this regime, axially symmetric hydromagnetic quasi-equilibrium states are possible, which are only constrained to have a vanishing azimuthal Lorentz force. In these states, the particle species deviate from chemical (β) equilibrium, which tends to be restored by β decays (Urca reactions), inducing fluid motions that change the magnetic field configuration. If the stars remained hot for a sufficiently long time, this evolution would eventually lead to a chemical equilibrium state, in which the fluid is barotropic and the magnetic field, if axially symmetric, satisfies the non-linear Grad–Shafranov equation. Here, we present a numerical scheme that decouples the magnetic and thermal evolution, enabling to efficiently perform, for the first time, long-term magnetothermal simulations in this regime for different magnetic field strengths and geometries. Our results demonstrate that, even for magnetar-strength fields |$\gtrsim 10^{16} \, \mathrm{G}$|⁠ , the feedback from the magnetic evolution on the thermal evolution is negligible. Thus, as the core passively cools, the Urca reactions quickly become inefficient at restoring chemical equilibrium, so the magnetic field evolves very little, and the Grad–Shafranov state is not attained. Therefore, any substantial evolution of the core magnetic field must occur later, in the 'weak-coupling' regime (⁠|$T\lesssim 5\times 10^8 \, \mathrm{K}$|⁠), when Urca reactions are frozen and ambipolar diffusion becomes relevant. [ABSTRACT FROM AUTHOR]

The dynamics in mergers of binary neutron star (BNS) systems depend sensitively on the equation of state (EOS) of dense matter. This has profound implications on the emission of gravitational waves (GWs) and the ejection of matter in the merger and post-merger phases and is thus of high interest for multimessenger astronomy. Today, a variety of nuclear EOSs are available with various underlying microphysical models. This calls for a study to focus on EOS effects from different physical nuclear matter properties and their influence on BNS mergers. We perform simulations of equal-mass BNS mergers with a set of nine different EOSs based on Skyrme density functionals. In the models, we systematically vary the effective nucleon mass, incompressibility, and symmetry energy at saturation density. This allows us to investigate the influence of specific nuclear matter properties on the dynamics of BNS mergers. We analyse the impact of these properties on the merger dynamics, the fate of the remnant, disc formation, ejection of matter, and GW emission. Our results indicate that some aspects of the merger, such as the frequencies of the post-merger GW spectrum and the shock-heated ejecta mass, are sensitive to the EOS around saturation density while others, such as the contraction of the remnant and the tidal ejecta mass, are sensitive to the behaviour towards higher densities, e.g. characterized by the slope of the pressure versus density. The detailed density dependence of the EOS thus needs to be taken into account to describe its influence on BNS mergers. [ABSTRACT FROM AUTHOR]

We investigate the structure of admixed neutron stars with a regular hadronic component and a fraction of fermionic self-interacting dark matter. Using two limiting equations of state for the dense baryonic interior, constructed from piecewise generalized polytropes, and an asymmetric self-interacting fermionic dark component, we analyse different scenarios of admixed neutron stars depending on the mass of dark fermions m χ, interaction mediators m ϕ, and self-interacting strengths g. We find that the contribution of dark matter to the masses and radii of neutron stars leads to tension with mass estimates of the pulsar J0453+1559, the least massive neutron star, and with the constraints coming from the GW170817 event. We discuss the possibilities of constraining dark matter model parameters g and y ≡ m χ/ m ϕ, using current existing knowledge on neutron star estimations of mass, radius, and tidal deformability, along with the accepted cosmological dark matter freeze-out values and self-interaction cross-section to mass ratio, σSI/ m χ, fitted to explain Bullet, Abell, and dwarf galaxy cluster dynamics. By assuming the most restrictive upper limit, σSI/ m χ < 0.1 cm2 g−1, along with dark matter freeze-out range values, the allowed g – y region is 0.01 ≲ g ≲ 0.1, with 0.5 ≲ y ≲ 200. For the first time, the combination of updated complementary restrictions is used to set constraints on self-interacting dark matter. [ABSTRACT FROM AUTHOR]

Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii, and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in significant biases in the inferred mass and radius, particularly for the high count rates that are nominally required to obtain meaningful constraints. In this follow-on study, we show that the bias can be mitigated by slicing the bursts into shorter segments where variability can be neglected, and jointly fitting the segments. Using this approach, the systematic uncertainties on the mass and radius are brought within the range of the statistical uncertainty. With about 106 source counts, this yields uncertainties of approximately 10 per cent for both the mass and radius. However, this modelling strategy requires substantial computational resources. We also confirm that the posterior distributions of the mass and radius obtained from multiple bursts of the same source can be merged to produce outcomes comparable to that of a single burst with an equivalent total number of counts. [ABSTRACT FROM AUTHOR]

We carefully develop the framework required to model the dynamical tidal response of a spinning neutron star in an inspiralling binary system, in the context of Newtonian gravity, making sure to include all relevant details and connections to the existing literature. The tidal perturbation is decomposed in terms of the normal oscillation modes, used to derive an expression for the effective Love number which is valid for any rotation rate. In contrast to previous work on the problem, our analysis highlights subtle issues relating to the orthogonality condition required for the mode-sum representation of the dynamical tide and shows how the prograde and retrograde modes combine to provide the overall tidal response. Utilizing a slow-rotation expansion, we show that the dynamical tide (the effective Love number) is corrected at first order in rotation, whereas in the case of the static tide (the static Love number) the rotational corrections do not enter until second order. [ABSTRACT FROM AUTHOR]

Magnetars form a special class of neutron stars possessing superstrong magnetic fields and demonstrating power flares triggered by these fields. Observations of such flares reveal the presence of quasi-periodic oscillations (QPOs) at certain frequencies; they are thought to be excited in the flares. QPOs carry potentially important information on magnetar structure, magnetic field, and mechanisms of magnetar activity. We calculate frequencies of torsional (magneto-elastic) oscillations of the magnetar crust treating the magnetic field effects in the first order of perturbation theory. The theory predicts the splitting of non-magnetic oscillation frequencies into Zeeman components. Zeeman splitting of the torsional oscillation spectrum of magnetars was suggested, clearly described and estimated by Shaisultanov and Eichler (2009), but their work has not been given considerable attention. To extend it, we suggest the technique of calculating oscillation frequencies, including Zeeman splitting at not too strong magnetic fields for arbitrary magnetic field configuration. Zeeman splitting enriches the oscillation spectrum and simplifies the theoretical interpretation of observations. We calculate several low-frequency oscillations of magnetars with a pure dipole magnetic field in the crust. The results qualitatively agree with the low-frequency QPOs detected in the hyperflare of SGR 1806–20 and in the giant flare of SGR 1900+14. [ABSTRACT FROM AUTHOR]

We use Bayesian analysis in order to constrain the equation of state for nuclear matter from astrophysical data related to the recent measurements from the NICER mission, LIGO/Virgo collaboration, and probability distributions of mass and radius from other 12 sources, including thermonuclear busters, and quiescent low-mass X-ray binaries. For this purpose, we base our study on a relativistic hadronic mean field model including an ω − ρ interaction. Our results indicate optimal ranges for some bulk parameters at the saturation density, namely, effective mass, incompressibility, and symmetry energy slope (L 0). For instance, we find |$L_0 = 50.79^{+15.16}_{-9.24}$|  MeV (Case 1) and |$L_0 = 75.06^{+8.43}_{-4.43}$|  MeV (Case 2) in a 68 per cent confidence interval for the two cases analysed (different input ranges for L 0 related to the PREX-II data). The respective parametrizations are in agreement with important nuclear matter constraints, as well as observational neutron star data, such as the dimensionless tidal deformability of the GW170817 event. From the mass–radius curves obtained from these best parametrizations, we also find the ranges of 11.97 km ≤ R 1.4 ≤ 12.73 km (Case 1) and 12.34 km ≤ R 1.4 ≤ 13.06 km (Case 2) for the radius of the |$1.4\, \mathrm{M}_\odot$| neutron star. [ABSTRACT FROM AUTHOR]

Motivated by the recent suggestions that very massive pulsar (PSR J0952−0607) and very light compact object (HESS J1731−347) exist, in this article, we revisit the possibility of such objects being strange stars instead of the standard hadronic neutron stars. We study the possible presence of local anisotropy and how it affects the macroscopic properties of strange stars and compare our results with the recent constraints presented in the literature. We found that the presence of anisotropy increases the maximum mass, the radius of the canonical star, and its tidal deformability for positive values of λBL and the opposite for negative values. We also show that although we cannot rule out the possibility of very compact objects being standard hadronic neutron stars, strange stars easily fulfill most of the observational constraints. [ABSTRACT FROM AUTHOR]

Neutron stars (and quark stars) have recently been found to follow certain universality, the I-Love-Q relations. These relations can provide a greater understanding of the structural and macro properties of compact astrophysical objects with knowledge of any one of the observable. The reason behind this is the lack of sensitivity to the relations with the equation of state (EoS) of matter. In this work, we have investigated the consistency of universal relations for a generic family of equations of state. These equations of state are consistent with the prescribed limits of nuclear theory (low densities) and perturbative QCD (asymptotically high densities) and respecting recent astronomical observations. At intermediate densities, the equations of states are parametrized by the adiabatic speed of sound. Although the spread in the EoS is significant, the universal nature of the I-Love-Q trio holds relatively well up to a specific tolerance limit of 5 per cent fractional percentage error. We have compared the result first with a similar set, but not respecting the recent astrophysical constraints. A deviation from universality is seen and the tolerance exceeds the 5 per cent limit. Next, we have compared our result with a family of EoS having monotonic sound speed inside a neutron star. A slight deviation from universality is seen for this case; however, the tolerance remains within the 5  per cent limit. This shows the existence of a correlation between universal relations and observational constraints. [ABSTRACT FROM AUTHOR]

Nuclear many-body theory is based on the tenet that nuclear systems can be accurately described as collections of point-like particles. This picture, while providing a remarkably accurate explanation of a wealth of measured properties of atomic nuclei, is bound to break down in the high-density regime, in which degrees of freedom other than protons and neutrons are expected to come into play. Valuable information on the validity of the description of dense nuclear matter in terms of nucleons, needed to firmly establish its limit of applicability, can be obtained from electron–nucleus scattering data at large momentum transfer and low energy transfer. The emergence of y-scaling in this kinematic region, unambiguously showing that the beam particles couple to high-momentum nucleons belonging to strongly correlated pairs, indicates that at densities as large as five times nuclear density—typical of the neutron star interior—nuclear matter largely behaves as a collection of nucleons.

We study cold strange quark stars employing an enhanced version of the quark-mass density-dependent model, which incorporates excluded volume effects to address non-perturbative QCD repulsive interactions. We provide a comparative analysis of our mass formula parametrization with previous models from the literature. We identify the regions within the parameter space where three-flavor quark matter is more stable than the most tightly bound atomic nucleus (stability window). Specifically, we show that excluded volume effects do not change the Gibbs free energy per baryon at zero pressure, rendering the stability window unaffected. The curves of pressure versus energy density exhibit various shapes—convex upward, concave downward, or nearly linear—depending on the mass parametrization. This behavior results in different patterns of increase, decrease, or constancy in the speed of sound as a function of baryon number density. We analyze the mass–radius relationship of strange quark stars, revealing a significant increase in maximum gravitational mass and a shift in the curves toward larger radii as the excluded volume effect intensifies. Excluded volume effects render our models compatible with all modern astrophysical constraints, including the properties of the recently observed low-mass compact object HESSJ1731.

Context. Proto-neutron stars are born hot, with temperatures exceeding a few times 1010 K. In these conditions, the crust of the proto-neutron star is expected to be made of a Coulomb liquid and composed of an ensemble of different nuclear species. Aims. In this work, we perform a study of the beta-equilibrated proto-neutron-star crust in the liquid phase in a self-consistent multi-component approach. This also allows us to perform a consistent calculation of the impurity parameter, which is often taken as a free parameter in cooling simulations. Methods. To this aim, we developed a self-consistent multi-component approach at finite temperature using a compressible liquid-drop description of the ions, with surface parameters adjusted to reproduce experimental masses. The treatment of the ion centre-of-mass motion was included through a translational free-energy term accounting for in-medium effects. The results of the self-consistent calculations of the multi-component plasma are systematically compared with those performed in a perturbative treatment as well as in the one-component plasma approximation. Results. We show that the inclusion of non-linear mixing terms arising from the ion centre-of-mass motion leads to a breakdown of the ensemble equivalence between the one-component and multi-component approach. Our findings also illustrate that the abundance of light nuclei becomes important and eventually dominates the whole distribution at higher density and temperature in the crust. This is reflected in the impurity parameter, which, in turn, may have a potential impact on neutron-star cooling. For practical application to astrophysical simulations, we also provide a fitting formula for the impurity parameter in the proto-neutron-star inner crust. Conclusions. Our results obtained within a self-consistent multi-component approach show important differences in the prediction of the proto-neutron-star composition with respect to those obtained with a one-component approximation or a perturbative multi-component approximation, particularly in the deeper region of the crust. This highlights the importance of a full, self-consistent multi-component plasma calculation for reliable predictions of the proto-neutron-star crust composition. [ABSTRACT FROM AUTHOR]

We model the nuclear evolution of an accreted matter as it sinks toward the stellar centre, in order to find its composition and equation of state. To this aim, we developed a simplified reaction network that allows for redistribution of free neutrons in the inner crust to satisfy the recently suggested neutron hydrostatic and diffusion equilibrium condition. We analyse the main reaction pathways for the three representative thermonuclear ash compositions: Superburst, Kepler, and Extreme rp. In contrast to the previous results, which neglect redistribution of free (unbound) neutrons in the inner crust, the most significant reactions in our calculations are neutron captures and electron emissions. The pycnonuclear fusion plays some role only for Kepler ashes. For the direct application of our results in astrophysical codes we present profiles of the average charge, 〈 Z 〉, impurity parameter, Q imp, and equation of state for a set of models, parametrized by the pressure at the outer–inner crust interface. Typically, for Superburst ashes Q imp ≈ 1 − 4, while for Kepler ashes Q imp decreases from ≈23 at the outer–inner crust interface to ≈5 at the end of our simulation (the corresponding density equals ρdc ≈ 2 × 1012 g cm−3). At the same time, for Extreme rp ashes Q imp remains large ≈30 − 35 in the considered inner crust region. Our results are important for modelling the thermal relaxation of transiently accreting neutron stars after the end of the outburst. [ABSTRACT FROM AUTHOR]

Strangeon star model has passed various observational tests, such as the massive pulsars and the tidal deformability during binary mergers. Pulsar glitch, as a useful probe for studying the interior structure of pulsars, has also been studied in strangeon star model in our previous papers, including the recovery coefficient, the waiting time of glitches, and glitch activity. In this paper, the recovery process of a glitch is described in the strangeon star model, based on the starquake picture established before (in Paper I). After the starquake, the inner motion of the stellar matter would reduce the tangential pressure in the cracked places at the equatorial plane. The recovery (increase) of the tangential pressure would be achieved by a viscous flow towards the cracked places at equatorial plane, which leads to the exponential recovery of the spin frequency. A uniform viscous flow can reproduce the single exponential decay observed in some glitches, and the viscous time-scale τ and the depth h of the cracking place below the surface can be fitted by the recovery data. It is found that h increases with glitch size Δν/ν, which is expected in the glitch scenario of strangeon stars. The magnitude of the recovery predicted in this recovery model is also consistent with that derived from observations. The single exponential decay reproduced by a uniform viscous flow can be generalized to two or more exponentials by the multicomponent of viscous flows. [ABSTRACT FROM AUTHOR]

Thermal emission of neutron stars in soft X-ray transients (SXTs) in a quiescent state is believed to be powered by the heat deposited in the stellar crust due to nuclear reactions during accretion (deep crustal heating paradigm). Confronting observations of SXTs with simulations helps to verify theoretical models of the dense matter in the neutron stars. Usually, such simulations were carried out assuming that the free neutrons and nuclei in the inner crust move together. A recently proposed thermodynamically consistent approach allows for independent motion of the free neutrons. We simulate the thermal evolution of the SXTs within the thermodynamically consistent approach and compare the results with the traditional approach and with observations. For the latter, we consider a collection of quasi-equilibrium thermal luminosities of the SXTs in quiescence and the observed neutron star crust cooling in SXT MXB 1659−29. We test different models of the equation of state and baryon superfluidity and take into account additional heat sources in the shallow layers of neutron-star crust (the shallow heating). We find that the observed quasi-stationary thermal luminosities of the SXTs can be equally well fitted using the traditional and thermodynamically consistent models, provided that the shallow heat diffusion into the core is taken into account. The observed crust cooling in MXB 1659−29 can also be fitted in the frames of both models, but the choice of the model affects the derived parameters responsible for the thermal conductivity in the crust and for the shallow heating. [ABSTRACT FROM AUTHOR]

We study the effects of the time-variable properties of thermonuclear X-ray bursts on modelling their millisecond-period burst oscillations. We apply the pulse profile modelling technique that is being used in the analysis of rotation-powered millisecond pulsars by the Neutron Star Interior Composition Explorer to infer masses, radii, and geometric parameters of neutron stars. By simulating and analysing a large set of models, we show that overlooking burst time-scale variability in temperatures and sizes of the hot emitting regions can result in substantial bias in the inferred mass and radius. To adequately infer neutron star properties, it is essential to develop a model for the time-variable properties or invest a substantial amount of computational time in segmenting the data into non-varying pieces. We discuss prospects for constraints from proposed future X-ray telescopes. [ABSTRACT FROM AUTHOR]

The coupled, time-dependent Gross–Pitaevskii and Ginzburg–Landau equations are solved simultaneously in three dimensions to investigate the equilibrium state and far-from-equilibrium, spin-down dynamics of an interpenetrating neutron superfluid and proton type-II superconductor, as an idealized description of the outer core of a neutron star. The simulations generalize previous calculations without the time-dependent Ginzburg–Landau equation, where proton feedback is absent. If the angle θ between the rotation and magnetic axes does not equal zero, the equilibrium state consists of geometrically complicated neutron vortex and proton flux-tube tangles, as the topological defects pin to one another locally but align with different axes globally. During spin-down, new types of motion are observed. For θ = 0, entire vortices pair rectilinearly with flux tubes and move together while pinned. For θ ≠ 0, vortex segments pair with segments from one or more flux tubes, and the paired segments move together while pinned. The degree to which proton feedback impedes the deceleration of the crust is evaluated as a function of θ and the pinning strength, η. Key geometric properties of vortex-flux-tube tangles, such as filament length, mean curvature, and polarity are analysed. It is found that proton feedback smooths the deceleration of the crust, reduces the rotational glitch sizes, and stabilizes the vortex tangle dynamics. The dimensionless control parameters in the simulations are mutually ordered to match what is expected in a real neutron star, but their central values and dynamic ranges differ from reality by many orders of magnitude due to computational limitations. [ABSTRACT FROM AUTHOR]

Based on the linear-mixing approach, we calculate the latent heat for crystallizing fully ionized 12C/16O and 16O/20Ne mixtures in white dwarf (WD) cores for two different parametrizations of the corrections to the linear-mixing energies and with account of ion quantum effects. We report noticeable composition-dependent deviations of the excess entropy in both directions from the standard value of 0.77 per ion. Within the same framework, we evaluate the excess entropy and released or absorbed heat accompanying the exsolution process in solidified WD layers. The inclusion of this effect is shown to be important for reliable interpretation of WD cooling data. We also analyse the latent heat of crystallization for the eutectic 12C/22Ne mixture, where we find a qualitative dependence of both the phase diagram and the latent heat behaviour on ion quantum effects. This may be important for the model with 22Ne distillation in cooling C/O/22Ne WD proposed as a solution for the ultramassive WD multi-Gyr cooling anomaly. Astrophysical implications of our findings for crystallizing WD are discussed. [ABSTRACT FROM AUTHOR]

Nuclear many-body theory is based on the tenet that nuclear systems can be accurately described as collections of point-like particles. This picture, while providing a remarkably accurate explanation of a wealth of measured properties of atomic nuclei, is bound to break down in the high-density regime, in which degrees of freedom other than protons and neutrons are expected to come into play. Valuable information on the validity of the description of dense nuclear matter in terms of nucleons, needed to firmly establish its limit of applicability, can be obtained from electron–nucleus scattering data at large momentum transfer and low energy transfer. The emergence of y-scaling in this kinematic region, unambiguously showing that the beam particles couple to high-momentum nucleons belonging to strongly correlated pairs, indicates that at densities as large as five times nuclear density—typical of the neutron star interior—nuclear matter largely behaves as a collection of nucleons. [ABSTRACT FROM AUTHOR]

Neutron stars are the densest objects in the Universe. They have a microscopically homogeneous core and heterogeneous crust. In particular, there may be a specific layer inside neutron stars, the mantle, which consists of substantially non-spherical nuclei immersed in a background of relativistic degenerate electrons and quasi-free neutrons. In this paper, we reconsider the transverse shear modulus for cylindrical phases of the mantle within the framework of the compressible liquid drop model. We demonstrate that transverse shearing affects the shape of nuclear clusters: their cross-section becomes elliptical. This effect reduces the respective elastic constant. Using a simple model, we perform all derivations analytically and obtain the expression for the transverse shear modulus, which can be useful for astrophysical applications. [ABSTRACT FROM AUTHOR]

The presence of nuclear pasta is expected to modify the transport properties in the mantle of neutron stars. The non-spherical geometry of the pasta nuclear clusters leads to anisotropies in the collision frequencies, impacting the thermal and electrical conductivity. We derive analytical expressions for the anisotropic collision frequencies using the Boltzmann equation in the relaxation time approximation. The average parallel, perpendicular, and Hall electrical conductivities are computed in the high-temperature regime above crustal melting, considering incoherent elastic electron-pasta scattering and randomly oriented pasta structures. Numerical values are obtained at different densities and temperatures by using the Indiana University–Florida State University (IUFSU) parametrization of the non-linear Walecka model to determine the crustal structure. We find that the anisotropy of the collision frequencies grows with the length of the pasta structures and, independently of the magnetic field, the presence of rod and slab phases decreases the conductivity by more than one order of magnitude. Our numerical results indicate that, even if the pasta structures might survive above the crustal melting point, no strong anisotropies are to be expected in the conduction properties in this temperature regime, even in the presence of a very high magnetic field. [ABSTRACT FROM AUTHOR]

Understanding the equation of state of dense neutron-rich matter remains a major challenge in modern physics and astrophysics. Neutron star observations from electromagnetic and gravitational wave spectra provide critical insights into the behavior of dense neutron-rich matter. The next generation of telescopes and gravitational wave detectors will offer even more detailed neutron-star observations. Employing deep learning techniques to map neutron star mass and radius observations to the equation of state allows for its accurate and reliable determination. This work demonstrates the feasibility of using deep learning to extract the equation of state directly from observations of neutron stars, and to also obtain related nuclear matter properties such as the slope, curvature, and skewness of nuclear symmetry energy at saturation density. Most importantly, it shows that this deep learning approach is able to reconstruct realistic equations of state and deduce realistic nuclear matter properties. This highlights the potential of artificial neural networks in providing a reliable and efficient means to extract crucial information about the equation of state and related properties of dense neutron-rich matter in the era of multi-messenger astrophysics. [ABSTRACT FROM AUTHOR]

Context. The crust of a neutron star is known to melt at a temperature that increases with increasing matter density, up to about 1010 K. At such high temperatures and beyond, the crustal ions are put into collective motion and the associated entropy contribution can affect both the thermodynamic properties and the composition of matter. Aims. We studied the importance of this effect in different thermodynamic conditions relevant to the inner crust of the proto-neutron star, both at beta equilibrium and in the fixed-proton-fraction regime. Methods. To this aim, we solved the hydrodynamic equations for an ion moving in an incompressible, irrotational, and non-viscous fluid, with different boundary conditions, thus leading to different prescriptions for the ion effective mass. We then employed a compressible liquid-drop approach in the one-component plasma approximation, including the renormalisation of the ion mass to account for the influence of the surrounding medium. Results. We show that the cluster size is determined by the competition between the ion centre-of-mass motion and the interface properties, namely the Coulomb, surface, and curvature energies. In particular, including the translational free energy in the minimisation procedure can significantly reduce the optimal number of nucleons in the clusters and lead to an early dissolution of clusters in dense beta-equilibrated matter. On the other hand, we find that the impact of translational motion is reduced in scenarios where the proton fraction is assumed constant and is almost negligible on the inner-crust equation of state. Conclusions. Our results show that the translational degrees of freedom affect the equilibrium composition of beta-equilibrated matter and the density and pressure of the crust-core transition in a non-negligible way, highlighting the importance of its inclusion when modelling the finite-temperature inner crust of the (proto-)neutron star. [ABSTRACT FROM AUTHOR]

Context. In the upcoming years, present and next-generation gravitational wave observatories will detect a larger number of binary neutron star (BNS) mergers with increasing accuracy. In this context, improving BNS merger numerical simulations is crucial to correctly interpret the data and constrain the equation of state (EOS) of neutron stars (NSs). Aims. State-of-the-art simulations of BNS mergers do not include muons. However, muons are known to be relevant in the microphysics of cold NSs and are expected to have a significant role in mergers, where the typical thermodynamic conditions favour their production. Our work is aimed at investigating the impact of muons on the merger remnant. Methods. We post-process the outcome of four numerical relativity simulations of BNS mergers performed with three different baryonic EOSs and two mass ratios considering the first 15 milliseconds after merger. We compute the abundance of muons in the remnant and analyse how muons affect the trapped neutrino component and the fluid pressure. Results. We find that depending on the baryonic EOS, the net fraction of muons is between 30% and 70% the net fraction of electrons. Muons change the flavour hierarchy of trapped (anti-)neutrinos such that deep inside the remnant, muon anti-neutrinos are the most abundant, followed by electron anti-neutrinos. Finally, muons and trapped neutrinos modify the neutron-to-proton ratio, affecting the remnant pressure by up to 7% when compared with calculations neglecting them. Conclusions. This work demonstrates that muons have a non-negligible effect on the outcome of BNS merger simulations, and they should be included to improve the accuracy of a simulation. [ABSTRACT FROM AUTHOR]

We explored the impact of the latest equation of state (EOS) for dense hydrogen–helium mixtures, which takes into account the interactions between hydrogen and helium species during the evolution of very low-mass stars and brown dwarfs (BDs). These interactions modify the thermodynamic properties of the H/He mixture, notably the entropy, a quantity of prime importance for these fully convective bodies, but also the onset and the development of degeneracy throughout the body. This translates into a faster cooling rate, that is, cooler isentropes for a given mass and age, and thus larger BD masses and smaller radii for a given effective temperature and luminosity than the models based on previous EOSs. This means that objects of a given mass and age in the range M ≲ 0.1 M⊙, τ ≳ 108 yr will have cooler effective temperatures and fainter luminosities. Confronting these new models with several observationally determined BD dynamical masses, we show that this improves the agreement between evolutionary models and observations and resolves at least part of the observed discrepancy between the properties of dynamical mass determinations and evolutionary models. A noticeable consequence of this improvement of the dense H/He EOS is that it yields a larger H-burning minimum mass, now found to be 0.075 M⊙ (78.5 MJup) with the ATMO atmosphere models for solar metallicity. These updated BD models are made publicly available. [ABSTRACT FROM AUTHOR]