Temperature-dependent Equation Of State Research Articles

Hyperonic uncertainties in neutron stars, mergers and supernovae

Abstract 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 and FSU2H*U. 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.

Impact of extreme spins and mass ratios on the post-merger observables of high-mass binary neutron stars

ABSTRACTThe gravitational-wave events GW170817 and GW190425 have led to a number of important insights on the equation of state of dense matter and the properties of neutron stars, such as their radii and the maximum mass. Some of these conclusions have been drawn on the basis of numerical-relativity simulations of binary neutron-star mergers with vanishing initial spins. While this may be a reasonable assumption in equal-mass systems, it may be violated in the presence of large mass asymmetries accompanied by the presence of high spins. To quantify the impact of high spins on multimessenger gravitational-wave events, we have carried out a series of high-mass binary neutron-star mergers with a highly spinning primary star and large mass asymmetries that have been modelled self-consistently using two temperature-dependent equations of state. We show that, when compared with equal-mass, irrotational binaries, these systems can lead to significant differences in the remnant lifetime, in the dynamical ejecta, in the remnant disc masses, in the secular ejecta, and on the bulk kilonova properties. These differences could be exploited to remove the degeneracy between low- and high-spin priors in the detection of gravitational waves from binary neutron-star mergers.

Spritz: a new fully general-relativistic magnetohydrodynamic code

The new era of multimessenger astrophysics requires the capability of studying different aspects of the evolution of compact objects. In particular, the merger of neutron star binaries is a strong source of gravitational waves and electromagnetic radiation, from radio to γ-rays, as demonstrated by the detection of GW170817 and its electromagnetic counterparts. In order to understand the physical mechanisms involved in such systems, it is necessary to employ fully general relativistic magnetohydrodynamic (GRMHD) simulations able to include the effects of a composition and temperature dependent equation of state describing neutron star matter as well as neutrino emission and reabsorption. Here, we present our new code named Spritz that solves the GRMHD equations in 3D Cartesian coordinates and on a dynamical spacetime. The code can support tabulated equations of state, taking into account finite temperature effects and allowing for the inclusion of neutrino radiation. In this first paper, we present the general features of the code and a series of tests performed in special and general relativity to assess the robustness of the basic GRMHD algorithms implemented. Among these tests, we also present the first comparison between a non-staggered and a staggered formulation of the vector potential evolution, which is used to guarantee the divergence-less character of the magnetic field. With respect to other publicly available GRMHD codes, Spritz combines the robust approach of a staggered formulation of the vector potential together with the use of an equation of state driver (EOS_Omni) that can allow the code to use finite temperature equations of state. A next version of the code will fully test the EOS_Omni driver by coupling it with a neutrino leakage scheme.

Numerical simulations of neutron star-black hole binaries in the near-equal-mass regime

Simulations of neutron star-black hole (NSBH) binaries generally consider black holes with masses in the range $(5-10)M_\odot$, where we expect to find most stellar mass black holes. The existence of lower mass black holes, however, cannot be theoretically ruled out. Low-mass black holes in binary systems with a neutron star companion could mimic neutron star-neutron (NSNS) binaries, as they power similar gravitational wave (GW) and electromagnetic (EM) signals. To understand the differences and similarities between NSNS mergers and low-mass NSBH mergers, numerical simulations are required. Here, we perform a set of simulations of low-mass NSBH mergers, including systems compatible with GW170817. Our simulations use a composition and temperature dependent equation of state (DD2) and approximate neutrino transport, but no magnetic fields. We find that low-mass NSBH mergers produce remnant disks significantly less massive than previously expected, and consistent with the post-merger outflow mass inferred from GW170817 for moderately asymmetric mass ratio. The dynamical ejecta produced by systems compatible with GW170817 is negligible except if the mass ratio and black hole spin are at the edge of the allowed parameter space. That dynamical ejecta is cold, neutron-rich, and surprisingly slow for ejecta produced during the tidal disruption of a neutron star : $v\sim (0.1-0.15)c$. We also find that the final mass of the remnant black hole is consistent with existing analytical predictions, while the final spin of that black hole is noticeably larger than expected -- up to $\chi_{\rm BH}=0.84$ for our equal mass case.

CO2 ice structure and density under Martian atmospheric conditions

Clouds composed of CO2 ice form throughout the Martian atmosphere. In the mesosphere, CO2 ice clouds are thought to form via heterogeneous ice nucleation on nanoparticles of meteoric origin at temperatures often below 100K. Lower altitude CO2 ice clouds in the wintertime polar regions form up to around 145K and lead to the build-up of the polar ice caps. However, the crystal structure and related fundamental properties of CO2 ice under Martian conditions are poorly characterised. Here we present X-ray diffraction (XRD) measurements of CO2 ice, grown via deposition from the vapour phase under temperature and pressure conditions analogous to the Martian mesosphere. A crystalline cubic structure was determined, consistent with the low-pressure polymorph (CO2-I, space group Pa-3 (No. 205)). CO2 deposited at temperatures of 80-130K and pressures of 0.01–1 mbar was consistent with dry ice and previous literature measurements, thus removing the possibility of a more complicated phase diagram for CO2 in this region. At 80K, a lattice parameter of 5.578 ± 0.002Å, cell volume of 173.554 ± 0.19Å3 and density of 1.684 ± 0.002gcm−3 was determined. Using these measurements, we determined the thermal expansion of CO2 across 80–130K that allowed for a fit of CO2 ice density measurements across a larger temperature range (80–195K) when combined with literature data (CO2 density= 1.72391-2.53 × 10−4T-2.87 × 10−6T2). Temperature-dependent CO2 density values are used to estimate sedimentation velocities and heterogeneous ice nucleation rates, showing an increase in nucleation rate of up to a factor of 1000 when compared to commonly used literature values. This temperature-dependent equation of state is therefore suggested for use in future studies of Martian mesospheric CO2 clouds. Finally, we discuss the possible shapes of crystals of CO2 ice in the Martian atmosphere and show that a range of shapes including cubes and octahedra as well as a combination of the two in the form of cubo-octahedra are likely.

Spacetime singularity resolution in Snyder noncommutative space

Inspired by quantum gravity proposal, we construct a deformed phase space which supports the UV and IR cutoffs. We show that the Liouville theorem is satisfied in the deformed phase space which allows us to formulate the thermodynamics of the early Universe in the semiclassical regime. Applying the proposed method to the Snyder noncommutative space, we find a temperature dependent equation of state which opens a new window for natural realization of inflation as a phase transition from quantum gravity regime to the standard radiation dominated era. Also we obtain finite energy and entropy densities for the Universe, when at least the Weak Energy Condition is satisfied. We show that there is a minimum size for the Universe which is proportional to the Planck length and consequently the Big Bang singularity is removed.

A new temperature-dependent equation of state of solids

In the present paper, a temperature-dependent equation of state (EOS) of solids is discussed which is found to be applicable in high-pressure and high-temperature range. Present equation of state has been applied in 18 solids. The calculated data are found in very good agreement with the data available from other sources.

Cosmological evolution of finite temperature Bose-Einstein condensate dark matter

Once the temperature of a bosonic gas is smaller than the critical, density dependent, transition temperature, a Bose - Einstein Condensation process can take place during the cosmological evolution of the Universe. Bose - Einstein Condensates are very strong candidates for dark matter, since they can solve some major issues in observational astrophysics, like, for example, the galactic core/cusp problem. The presence of the dark matter condensates also drastically affects the cosmic history of the Universe. In the present paper we analyze the effects of the finite dark matter temperature on the cosmological evolution of the Bose-Einstein Condensate dark matter systems. We formulate the basic equations describing the finite temperature condensate, representing a generalized Gross-Pitaevskii equation that takes into account the presence of the thermal cloud in thermodynamic equilibrium with the condensate. The temperature dependent equations of state of the thermal cloud and of the condensate are explicitly obtained in an analytical form. By assuming a flat Friedmann-Robertson-Walker (FRW) geometry, the cosmological evolution of the finite temperature dark matter filled Universe is considered in detail in the framework of a two interacting fluid dark matter model, describing the transition from the initial thermal cloud to the low temperature condensate state. The dynamics of the cosmological parameters during the finite temperature dominated phase of the dark matter evolution are investigated in detail, and it is shown that the presence of the thermal excitations leads to an overall increase in the expansion rate of the Universe.

Numerical simulation of the hydrodynamical combustion to strange quark matter

We present results from a numerical solution to the burning of neutron matter inside a cold neutron star into stable (u,d,s) quark matter. Our method solves hydrodynamical flow equations in 1D with neutrino emission from weak equilibrating reactions, and strange quark diffusion across the burning front. We also include entropy change due to heat released in forming the stable quark phase. Our numerical results suggest burning front laminar speeds of 0.002-0.04 times the speed of light, much faster than previous estimates derived using only a reactive-diffusive description. Analytic solutions to hydrodynamical jump conditions with a temperature dependent equation of state agree very well with our numerical findings for fluid velocities. The most important effect of neutrino cooling is that the conversion front stalls at lower density (below approximately 2 times saturation density). In a 2-dimensional setting, such rapid speeds and neutrino cooling may allow for a flame wrinkle instability to develop, possibly leading to detonation.

Thermal effects on nuclear symmetry energy with a momentum-dependent effective interaction

The knowledge of the nuclear symmetry energy of hot neutron-rich matter is important for understanding the dynamical evolution of massive stars and the supernova explosion mechanisms. In particular, the electron capture rate on nuclei and/or free protons in presupernova explosions is especially sensitive to the symmetry energy at finite temperature. In view of the above, in the present work we calculate the symmetry energy as a function of the temperature for various values of the baryon density by applying a momentum-dependent effective interaction. In addition to a previous work, the thermal effects are studied separately both in the kinetic part and the interaction part of the symmetry energy. We focus also on the calculations of the mean-field potential, employed extensively in heavy-ion reaction research, both for nuclear and pure neutron matter. The proton fraction and the electron chemical potential, which are crucial quantities for representing the thermal evolution of supernova and neutron stars, are calculated for various values of the temperature. Finally, we construct a temperature dependent equation of state of {beta}-stable nuclear matter, the basic ingredient for the evaluation of the neutron star properties.

From microscales to macroscales in 3D: selfconsistent equation of state for supernova and neutron star models

First results from a fully self-consistent, temperature-dependent equation of state that spans the density range of neutron stars and supernova cores above neutron drip density are presented. The equation of state (EoS) is calculated using a mean-field Hartree-Fock method in three dimensions (3D). The nuclear interaction is represented by the phenomenological Skyrme model in this work, but the EoS can be obtained in our framework for any suitable form of the nucleon-nucleon effective interaction. The scheme we employ naturally allows effects such as (i) neutron drip, which results in an external neutron gas, (ii) the variety of exotic nuclear shapes expected for extremely neutron heavy nuclei, and (iii) the subsequent dissolution of these nuclei into nuclear matter. In this way, the equation of state is calculated across phase transitions without recourse to interpolation techniques between density regimes described by different physical models. EoS tables are calculated in the wide range of densities, temperature and proton/neutron ratios on the ORNL NCCS XT3, using up to 2000 processors simultaneously.

High-resolution calculations of merging neutron stars - II. Neutrino emission

The remnant resulting from the merger of two neutron stars produces neutrinos in copious amounts. In this paper we present the neutrino emission results obtained via Newtonian, high-resolution simulations of the coalescence event. These simulations use three-dimensional smoothed particle hydrodynamics together with a nuclear, temperature dependent equation of state and a multi-flavour neutrino leakage scheme. We present the details of our scheme, discuss the neutrino emission results from a neutron star coalescence and compare them to the core-collapse supernova case where neutrino emission has been studied for several decades. The average neutrino energies are similar to those in the supernova case, but contrary to the latter, the luminosities are dominated by electron-type antineutrinos which are produced in the hot, neutron-rich, thick disk of the merger remnant. The cooler parts of this disk contain substantial fractions of heavy nuclei, which, however, do not influence the overall neutrino emission results significantly. Our total neutrino luminosities from the merger event are considerably lower than those found in previous investigations. This has serious consequences for the ability to produce a gamma-ray burst via neutrino annihilation. The neutrinos are emitted preferentially along the initial binary rotation axis, an event seen ``pole-on'' would appear much brighter in neutrinos than a similar event seen ``edge-on''.

Comparison of Three Temperature-Dependent Equations of State of Condensed Matter

Considering temperature effects, we compare the equation of Vinet et al., the Birch equation and the K-D equation with experimental data. The results show that the Birch equation is the best.