Constraints on the symmetry energy using the mass-radius relation of neutron stars

The nuclear symmetry energy is intimately connected with nuclear astrophysics. This contribution focuses on the estimation of the symmetry energy from experiment and how it is related to the structure of neutron stars. The most important connection is between the radii of neutron stars and the pressure of neutron star matter in the vicinity of the nuclear saturation density n_s. This pressure is essentially controlled by the nuclear symmetry energy parameters S_v and L , the first two coefficients of a Taylor expansion of the symmetry energy around n_s. We discuss constraints on these parameters that can be found from nuclear experiments. We demonstrate that these constraints are largely model-independent by deriving them qualitatively from a simple nuclear model. We also summarize how recent theoretical studies of pure neutron matter can reinforce these constraints. To date, several different astrophysical measurements of neutron star radii have been attempted. Attention is focused on photospheric radius expansion bursts and on thermal emissions from quiescent low-mass X-ray binaries. While none of these observations can, at the present time, determine individual neutron star radii to better than 20% accuracy, the body of observations can be used with Bayesian techniques to effectively constrain them to higher precision. These techniques invert the structure equations and obtain estimates of the pressure-density relation of neutron star matter, not only near n_s, but up to the highest densities found in neutron star interiors. The estimates we derive for neutron star radii are in concordance with predictions from nuclear experiment and theory.     

中子星物质中的K介子凝聚与三体核力效应

利用Brueckner-Hartree-Fock方法，计算了β稳定中子星物质的状态方程以及三体核力的影响，特别是研究了三体核力对中子星物质中K介子凝聚的影响. 结果表明三体核力对β稳定中子星物质中出现K介子凝聚的临界密度以及中子星物质中各种粒子所占的比例均有重要影响. 三体核力的主要作用是降低了中子星物质中出现K介子凝聚的临界密度并使K凝聚相中的核物质更加接近于对称核物质.     

Crust-core properties of neutron stars in the Nambu–Jona-Lasinio model

We adopt the Nambu–Jona-Lasinio(NJL) model to study the crust-core transition properties in neutron stars(NSs). For a given momentum cutoff and symmetry energy of saturation density in the NJL model, decreasing the slope of the symmetry energy gives rise to an increase in the crust-core transition density and transition pressure.Given the slope of the symmetry energy at saturation density, the transition density and corresponding transition pressure increase with increasing symmetry energy. The increasing trend between the fraction of the crustal moment of inertia and the slope of symmetry energy at saturation density indicates that a relatively large momentum cutoff of the NJL model is preferred. For a momentum cutoff of 500 Me V, the fraction of the crustal moment of inertia clearly increases with the slope of symmetry energy at saturation density. Thus, at the required fraction(7%) of the crustal moment of inertia, the NJL model with momentum cutoff of 500 Me V and a large slope of the symmetry energy of saturation density can give the upper limit of the mass of the Vela pulsar to be above 1.40 M_⊙.     

原子核对称能对中子星壳层结构的影响

中子星内壳层中存在原子核、中子、电子等非均匀分布的物质。在Wigner-Seitz近似下,共存相方法和自洽Thomas-Fermi近似方法是描述这种非均匀物质的有效方法。中子在非均匀物质所占的比例远远大于其他组分,因此原子核的对称能对非均匀物质的性质会产生十分重要的影响,而原子核对称能的密度依赖关系在核物质饱和密度附近有较大的不确定性。采用相对论平均场理论描述核子间相互作用,研究原子核对称能对中子星内壳层的密度范围、pasta相结构、壳核相变密度等性质的影响,探寻其中可能存在的关联。计算结果表明,原子核对称能及其密度依赖性在决定中子星内壳层非均匀物质的性质中起着重要作用,这与之前相关研究中得到的结论基本相符。Within Wigner-Seitz approximation, both the coexisting phases method and the self-consistent Thomas-Fermi approximation can be used to describe the nonuniform matter consisting of nuclei, neutrons, and electrons, which may coexist in the inner crust of neutron star. Since the neutron fraction is very large, nuclear symmetry energy may have an important impact on the properties of nonuniform matter. However, the density dependence of nuclear symmetry energy around saturation density is still rather uncertain. This paper focuses on the influence of nuclear symmetry energy on the density range of inner crust, pasta phase structure, and crust-core transition density of neutron star, where the relativistic mean field theory is adopted to describe the nucleon-nucleon interaction. It is turned out that the nuclear symmetry energy and its density dependence play an import role in determining the properties of nonuniform matter in the inner crust of neutron star, which is consistent with the former related studies.     

Constraints on the symmetry energy from observational probes of the neutron star crust

A number of observed phenomena associated with individual neutron star systems or neutron star populations find explanations in models in which the neutron star crust plays an important role. We review recent work examining the sensitivity to the slope of the symmetry energy L of such models, and constraints extracted on L from confronting them with observations. We focus on six sets of observations and proposed explanations: i) The cooling rate of the neutron star in Cassiopeia A, confronting cooling models which include enhanced cooling in the nuclear pasta regions of the inner crust; ii) the upper limit of the observed periods of young X-ray pulsars, confronting models of magnetic field decay in the crust caused by the high resistivity of the nuclear pasta layer; iii) glitches from the Vela pulsar, confronting the paradigm that they arise due to a sudden recoupling of the crustal neutron superfluid to the crustal lattice after a period during which they were decoupled due to vortex pinning; iv) the frequencies of quasi-periodic oscillations in the X-ray tail of light curves from giant flares from soft gamma-ray repeaters, confronting models of torsional crust oscillations; v) the upper limit on the frequency to which millisecond pulsars can be spun-up due to accretion from a binary companion, confronting models of the r-mode instability arising above a threshold frequency determined in part by the viscous dissipation timescale at the crust-core boundary; and vi) the observations of precursor electromagnetic flares a few seconds before short gamma-ray bursts, confronting a model of crust shattering caused by resonant excitation of a crustal oscillation mode by the tidal gravitational field of a companion neutron star just before merger.     

Properties of nuclear pastas

In this review we study the nuclear pastas as they are expected to be formed in neutron star crusts. We start with a study of the pastas formed in nuclear matter (composed of protons and neutrons), we follow with the role of the electron gas on the formation of pastas, and we then investigate the pastas in neutron star matter (nuclear matter embedded in an electron gas).Nuclear matter (NM) at intermediate temperatures (1 MeV ≲ T ≲ 15 MeV), at saturation and sub-saturation densities, and with proton content ranging from 30% to 50% was found to have liquid, gaseous and liquid–gas mixed phases. The isospin-dependent phase diagram was obtained along with the critical points, and the symmetry energy was calculated and compared to experimental data and other theories. At low temperatures (T ≲ 1 MeV) NM produces crystal-like structures around saturation densities, and pasta-like structures at sub-saturation densities. Properties of the pasta structures were studied with cluster-recognition algorithms, caloric curve, the radial distribution function, the Lindemann coefficient, Kolmogorov statistics, Minkowski functionals; the symmetry energy of the pasta showed a connection with its morphology.Neutron star matter (NSM) is nuclear matter embedded in an electron gas. The electron gas is included in the calculation by the inclusion of an screened Coulomb potential. To connect the NM pastas with those in neutron star matter (NSM), the role the strength and screening length of the Coulomb interaction have on the formation of the pastas in NM was investigated. Pasta was found to exist even without the presence of the electron gas, but the effect of the Coulomb interaction is to form more defined pasta structures, among other effects. Likewise, it was determined that there is a minimal screening length for the developed structures to be independent of the cell size.Neutron star matter was found to have similar phases as NM, phase transitions, symmetry energy, structure function and thermal conductivity. Like in NM, pasta forms at around T ≈ 1.5 MeV, and liquid-to-solid phase changes were detected at T ≈ 0.5 MeV. The structure function and the symmetry energy were also found to depend on the pasta structures.     

Neutron star matter superfluidity and superconductivity

The possibility of superfluidity or superconductivity in neutron or proton subsystems in the nuclear-matter region in neutron stars is investigated. The energy gap and corresponding critical temperature and critical magnetic field is calculated or estimated as function of density or Fermi momentum. In the calculations are used reaction matrix elements calculated earlier by means of Brueckner theory by the author. The final results indicate that neutron superfluidity, corresponding specifically toS-state pairing, may exist in a low-density shell in the nuclear-matter region of a neutron star. There is probably anisotropic neutron superfluidity, corresponding to the³ P ₂ or the singletD state, for higher densities. Superfluidity or superconductivity, corresponding toS-state pairing for the proton subsystem, is quite likely in most of the nuclear-matter region. The expected temperatures and magnetic fields in neutron stars seem to be well below the estimated critical temperatures or critical magnetic fields corresponding to the calculated values of the energy gap. However, similar methods have earlier predicted a much too high critical temperature for liquid³He.     

基于中能重离子碰撞研究高密对称能

用于描述核物质中子质子单粒子能量之差的对称能，最近20年得到了核物理学界的广泛关注。在饱和核密度附近，对称能的数值及斜率基本得到约束，然而其高密行为至今仍具有很大的不确定性。当前，探测对称能的研究计划正在世界范围内能够提供放射性束流的实验室展开。伴随着对称能相关实验的规划、进行，发展更加先进的同位旋依赖的输运理论模型变得非常必要。我们将核子-核子短程关联及介质中同位旋依赖的重子-重子非弹性散射截面等研究的新进展融入到了同位旋依赖的输运模型里面；探索了新的敏感于高密对称能的可观测量，比如挤出核子的中质比、光子、轻碎片，以及包含奇异隐奇异夸克的介子产生等；提出了高密对称能探测盲点的问题并给出解决办法；对于常见的对称能敏感观测量的模型预言的不确定性进行了彻底细致的研究；提出采用定性观测量，比如高能出射粒子的中质比，将高密对称能进行定性约束；率先提出并研究了对称能敏感观测量的探测密度区间问题，指出对称能敏感观测量的探测密度往往小于核反应最大压缩密度；发现核子-核子短程关联明显削弱观测量的对称能效应；考虑到饱和点处对称能斜率的约束范围，基于输运模型，提出通过探测对称能的曲率来约束高密对称能。除了利用重离子碰撞约束高密对称能之外，人们也可以通过与中子星相关的大量天体观测来间接约束高密对称能。     

Nuclear symmetry energy effects in finite nuclei and neutron star

Within a relativistic mean-field model with nonlinear isoscalar–isovector coupling, we explore the possibility of constraining the density dependence of nuclear symmetry energy from a systematic study of the neutron skin thickness of finite nuclei and neutron star properties. We find the present skin data supports a rather stiff symmetry energy at subsaturation densities that corresponds to a soft symmetry energy at supranormal densities. Correlation between the skin of ²⁰⁸Pb and the neutron star masses and radii with kaon condensation has been studied. We find that ²⁰⁸Pb skin estimate suggest star radii that reveals considerable model dependence. Thus precise measurements of neutron star radii in conjunction with skin thickness of heavy nuclei could provide significant constraint on the density dependence of symmetry energy.     

Unified neutron star EOSs and neutron star structures in RMF models

In the framework of the Thomas-Fermi approximation, we systematically study the EOSs and microscopic structures of neutron star matter in a vast density range with n_b ≈ 10⁻¹⁰-2 fm⁻³, where various covariant density functionals are adopted, i.e., those with nonlinear self couplings (NL3, PK1, TM1, GM1, MTVTC) and density-dependent couplings (DD-LZ1, DDME-X, PKDD, DD-ME2, DD2, TW99). It is found that the EOSs generally coincide with each other at n_b ≲ 10⁻⁴ fm⁻³ and 0.1 fm⁻³ ≲ n_b ≲ 0.3 fm⁻³, while in other density regions they are sensitive to the effective interactions between nucleons. By adopting functionals with a larger slope of symmetry energy L, the curvature parameter K_sym and neutron drip density generally increases, while the droplet size, proton number of nucleus, core-crust transition density, and onset density of non-spherical nuclei, decrease. All functionals predict neutron stars with maximum masses exceeding the two-solar-mass limit, while those of DD2, DD-LZ1, DD-ME2, and DDME-X predict optimum neutron star radii according to the observational constraints. Nevertheless, the corresponding skewness coefficients J are much larger than expected, while only the functionals MTVTC and TW99 meet the start-of-art constraints on J. More accurate measurements on the radius of PSR J0740 + 6620 and the maximum mass of neutron stars are thus essential to identify the functional that satisfies all constraints from nuclear physics and astrophysical observations. Approximate linear correlations between neutron stars' radii at M = 1.4M_⊙ and 2M_⊙, the slope L and curvature parameter K_sym of symmetry energy are observed as well, which are mainly attributed to the curvature-slope correlations in the functionals adopted here. The results presented here are applicable for investigations of the structures and evolutions of compact stars in a unified manner.     

Links between heavy ion and astrophysics

Heavy-ion experiments provide important data to test astrophysical models. The high-density equation of state can be probed in HI collisions and applied to the hot protoneutron star formed in core collapse supernovae. The parity radius experiment (PREX) aims to accurately measure the neutron radius of ²⁰⁸Pb with parity-violating electron scattering. This determines the pressure of neutron-rich matter and the density dependence of the symmetry energy. Competition between nuclear attraction and Coulomb repulsion can form exotic shapes called nuclear pasta in neutron star crusts and supernovae. This competition can be probed with multifragmentation HI reactions. We use large-scale semiclassical simulations to study nonuniform neutron-rich matter in supernovae. We find that the Coulomb interactions in astrophysical systems suppress density fluctuations. As a result, there is no first-order liquid-vapor phase transition. Finally, the virial expansion for low-density matter shows that the nuclear vapor phase is complex with significant concentrations of alpha particles and other light nuclei in addition to free nucleons.     

Effects of the density dependence of the symmetry energy on neutron stars

In this paper,we include the density dependence behavior of the symmetry energy in the improved quark mass density dependent (IQMDD) model.Under the mean field approximation,this model is applied to investigate neutron star matter and neutron stars successfully.Effects of the density dependence of the symmetry energy on neutron stars are described.     

Constraints on Non-Radial Oscillation Modes by Symmetry Energy Slope

The effect of the nuclear symmetry energy slope on the non-radial oscillation f-modes in neutron stars is calculated and discussed. Based on a conservative range of the symmetry energy slope constrained by the experiment and theoretical analysis, a constraint on the frequency and damping time of the gravitational radiation from the f-mode in neutron star is obtained. It is also shown that a higher symmetry energy slope corresponds with a smaller frequency and a longer damping time. Meanwhile, a new set of parameters is given to present the universal properties of the scaled frequency and damping time.     

中子星内部强子-夸克相变的有限尺度效应研究

在致密星体内部极高密度条件下,强子物质可能发生退禁闭相变成为夸克物质,即强子-夸克相变。这种相变过程对于中子星的性质有着重要影响。考虑库仑能和表面能的影响,即有限尺度效应,相变过程中的混杂相包含了被称为pasta相的几何结构。强子-夸克共存相的平衡条件是通过求总能量的最小值得到的。采用相对论平均场（RMF）模型来描述强子物质相,采用Nambu-Jona-Lasinio（NJL）模型来描述夸克物质相。有限尺度效应一定程度上增加了中子星的最大质量,增加幅度取决于强子-夸克表面张力的大小。有限尺度效应能够降低混杂相的范围,其结果介于Gibbs结构和Maxwell结构的结果之间。研究结果表明,中子星中可能包含一个混杂相的核心部分,其大小受到表面张力等参数的影响。It is generally considered that hadron matter may undergo a deconfinement phase transition becoming quark matter at very high density in massive neutron stars. This hadron-quark phase transition has important impact on neutron stars, which has received much attention. We consider finite-size effect in this phase transition process, which contains the impact of Coulomb energy and surface energy. By including this effect, the mixed phase forms the pasta structures. The equilibrium conditions for coexisting hadronic and quark phases are derived by minimizing the total energy including the surface and Coulomb contributions. We employ the relativistic mean-field(RMF) model to describe the hadronic phase, while the Nambu-Jona-Lasinio(NJL) model is used for the quark phase. We conclude that the finite-size effect will raise the stiffness of EOS, and then increase the maximum mass of neutron stars, which depend on the value of surface tension. Our results show that finite-size effects can significantly reduce the region of the mixed phase, and the results lie between those from the Gibbs and Maxwell constructions. We show that a massive star may contain a mixed phase core and its size depends on the surface tension of the hadron-quark interface.     

Nuclear symmetry energy in terms of single-nucleon potential and its effect on the proton fraction of <Emphasis Type="Italic">β</Emphasis>-stable <Emphasis Type="Italic">npeμ</Emphasis> matter

Momentum and density dependence of single-nucleon potential u_τ (k, ρ, β) is analyzed using a density dependent finite range effective interaction of the Yukawa form. Depending on the choice of the strength parameters of exchange interaction, two different trends of the momentum dependence of nuclear symmetry potential are noticed which lead to two opposite types of neutron and proton effective mass splitting. The 2nd-order and 4th-order symmetry energy of isospin asymmetric nuclear matter are expressed analytically in terms of the single-nucleon potential. Two distinct behavior of the density dependence of 2nd-order and 4th-order symmetry energy are observed depending on neutron and proton effective mass splitting. It is also found that the 4th-order symmetry energy has a significant contribution towards the proton fraction of β-stable npeμ matter at high densities.     

Effects of the Symmetry Energy and its Slope on Neutron Star Properties

In this work, we study the influence of the symmetry energy and its slope on three major properties of neutron stars: the maximum mass, the radii of the canonical 1.4 M _⊙ and the minimum mass that enables the direct Urca effect. We utilize four parametrizations of the relativistic quantum hadrodynamics and vary the symmetry energy within accepted values. We see that although the maximum mass is almost independent of it, the radius of the canonical 1.4M _⊙ and the mass that enables the direct Urca effect is strongly correlated with the symmetry energy and its slope. Also, since we expect that the radius grows with the slope, a theoretical limit arises when we increase this quantity above certain values.     

Neutron stars

Neutron stars are laboratories for dense matter and gravitational physics. Observations of neutron stars from sources such as radio pulsars, low-mass X-ray binaries, X-ray bursts and thermally-emitting neutron stars are setting bounds to neutron star masses, radii, rotation rates, temperatures and ages. Mass measurements constrain the equation of state at the highest densities and set firm bounds to the highest possible density of cold matter. Radii constrain the equation of state in the vicinity of the nuclear saturation density and yield information about the density dependence of the nuclear symmetry energy. Laboratory measurements and theoretical studies of pure neutron matter are in remarkable agreement with observational bounds.     

The lead nucleus as a miniature surrogate for a neutron star

The nucleus of ²⁰⁸Pb, a system 18 orders of magnitude smaller and 55 orders of magnitude lighter than a neutron star, may be used as a miniature surrogate to establish important correlations between its neutron skin and several neutron-star properties. Indeed, models with a thicker neutron skin in ²⁰⁸Pb generate larger neutron stars that have a lower liquid-to-solid transition density. Further, we illustrate how the correlation between the neutron skin in ²⁰⁸Pb and the radius of a 1.4 solar-mass neutron star may be used to place important constraints on the equation of state of neutron-rich matter and how it may help elucidate the existence of a phase transition at the core of the star.