Auxiliary Function Approach for Determining Symmetry Energy at Suprasaturation Densities

Bao-Jun Cai,Bao-An Li
DOI: https://doi.org/10.1103/physrevc.103.054611
IF: 3.199
2021-01-01
Physical Review C
Abstract:Nuclear symmetry energy E sym ( ρ ) at density ρ is normally expanded or simply parameterized as a function of χ = ( ρ − ρ 0 ) / 3 ρ 0 in the form of E sym ( ρ ) ≈ S + L χ + 2 − 1 K sym χ 2 + 6 − 1 J sym χ 3 + ··· using its magnitude S , slope L , curvature K sym and skewness J sym at the saturation density ρ 0 of nuclear matter. Much progress has been made in recent years in constraining especially the S and L parameters using various terrestrial experiments and astrophysical observations. However, such kind of expansions/parameterizations do not converge at suprasaturation densities where χ is not small enough, hindering an accurate determination of high-density E sym ( ρ ) even if its characteristic parameters at ρ 0 are all well determined by experiments/observations. By expanding the E sym ( ρ ) in terms of a properly chosen auxiliary function Π sym ( χ, Θ sym ) with a parameter Θ sym fixed accurately by an experimental E sym ( ρ r ) value at a reference density ρ r , we show that the shortcomings of the χ -expansion can be completely removed or significantly reduced in determining the high-density behavior of E sym ( ρ ) . In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher χ -order contributions and converges to the same E sym ( ρ ) much faster than the conventional χ -expansion at densities . 3 ρ 0 . Moreover, the still poorly constrained skewness J sym plays a small role in determining the E sym ( ρ ) at these densities in the auxiliary function approach. The new approach thus provides a nearly model-independent constraint on the E sym ( ρ ) at densities . 3 ρ 0 . Several quantitative demonstrations using Monte Carlo simulations are given.
What problem does this paper attempt to address?