Auxiliary Function Approach for Determining Symmetry Energy at Supra-saturation Densities
Abstract
Nuclear symmetry energy Esym() at density is normally expanded or simply parameterized as a function of =(-0)/30 in the form of Esym()≈ S+L+2-1Ksym2+6-1Jsym3+·s using its magnitude S, slope L , curvature Ksym and skewness Jsym 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 supra-saturation densities where is not small enough, hindering an accurate determination of high-density Esym() even if its characteristic parameters at 0 are all well determined by experiments/observations. By expanding the Esym() in terms of a properly chosen auxiliary function sym(,sym) with a parameter sym fixed accurately by an experimental Esym(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 Esym(). In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher -order contributions and converges to the same Esym() much faster than the conventional -expansion at densities 30. Several quantitative demonstrations using Monte Carlo simulations are given.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.