“…The nuclear symmetry energy E sym (n, α) at baryonic density n that appears quadratically in α = (N − P )/(N + P ) in the energy per nucleon E(n, α) where P (N ) is the number of protons (neutrons) in many nucleon systems given by E(n, α) = E(n, α = 0) + E sym (n)α 2 + O(α 4 ) + · · · , (1) plays the most important role in the equation of state (EoS)) for compact-stars [1,2]. In standard nuclear physics approaches (SNPAs) anchored on the effective density functional (EDF) such as the Skyrme potential, relativistic mean field (RMF) and varieties thereof ( [3] for a recent review) and on standard chiral perturbation (SχPT) expansion up to manageable chiral order (e.g., [4]), equipped with a certain number of parameters fit to available empirical data, the E(n, α) can be more or less reliably determined in the vicinity of the nuclear matter equilibrium density n 0 ∼ 0.16 fm −3 . It has also been extended, with a rather broad range of uncertainty, up to slightly above n 0 from heavy-ion collision experiments.…”