Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIsospin-Asymmetry Dependence of the Thermodynamic Nuclear Equation of State in Many-Body Perturbation Theory
82
0
0.0
(
0
)
Authors
Corbinian Wellenhofer
Ask ChatGPT about the research
No Arabic abstract
The computation of the thermodynamic properties of nuclear matter is a central task of theoretical nuclear physics. The nuclear equation of state is an essential quantity in nuclear astrophysics and governs the properties of neutron stars and core-collapse supernovae. The framework of chiral effective field theory provides the basis for the description of nuclear interactions in terms of a systematic low-energy expansion. In this thesis, we apply chiral two- and three-nucleon interactions in perturbative many-body calculations of the thermodynamic equation of state of infinite homogeneous nuclear matter. The conceptual issues that arise concerning the consistent generalization of the standard zero-temperature form of many-body perturbation theory to finite temperatures are investigated in detail. The structure of many-body perturbation theory at higher orders is examined, in particular concerning the role of the so-called anomalous contributions. The first-order nuclear liquid-gas phase transition is analyzed with respect to its dependence on temperature and the neutron-to-proton ratio. Furthermore, the convergence behavior of the expansion of the equation of state in terms of the isospin asymmetry is examined. It is shown that the expansion coefficients beyond the quadratic order diverge in the zero-temperature limit, implying a nonanalytic form of the isospin-asymmetry dependence at low temperatures. This behavior is associated with logarithmic terms in the isospin-asymmetry dependence at zero temperature.
rate research
Read More
The isospin-asymmetry dependence of the nuclear matter equation of state obtained from microscopic chiral two- and three-body interactions in second-order many-body perturbation theory is examined in detail. The quadratic, quartic and sextic coefficients in the Maclaurin expansion of the free energy per particle of infinite homogeneous nuclear matter with respect to the isospin asymmetry are extracted numerically using finite differences, and the resulting polynomial isospin-asymmetry parametrizations are compared to the full isospin-asymmetry dependence of the free energy. It is found that in the low-temperature and high-density regime where the radius of convergence of the expansion is generically zero, the inclusion of higher-order terms beyond the leading quadratic approximation leads overall to a significantly poorer description of the isospin-asymmetry dependence. In contrast, at high temperatures and densities well below nuclear saturation density, the interaction contributions to the higher-order coefficients are negligible and the deviations from the quadratic approximation are predominantly from the noninteracting term in the many-body perturbation series. Furthermore, we extract the leading logarithmic term in the isospin-asymmetry expansion of the equation of state at zero temperature from the analysis of linear combinations of finite differences. It is shown that the logarithmic term leads to a considerably improved description of the isospin-asymmetry dependence at zero temperature.
We investigate the order-by-order convergence behavior of many-body perturbation theory (MBPT) as a simple and efficient tool to approximate the ground-state energy of closed-shell nuclei. To address the convergence properties directly, we explore perturbative corrections up to 30th order and highlight the role of the partitioning for convergence. The use of a simple Hartree-Fock solution to construct the unperturbed basis leads to a convergent MBPT series for soft interactions, in contrast to, e.g., a harmonic oscillator basis. For larger model spaces and heavier nuclei, where a direct high-order MBPT calculation in not feasible, we perform third-order calculation and compare to advanced ab initio coupled-cluster calculations for the same interactions and model spaces. We demonstrate that third-order MBPT provides ground-state energies for nuclei up into tin isotopic chain that are in excellent agreement with the best available coupled-cluster results at a fraction of the computational cost.
We discuss experimental evidence for a nuclear phase transition driven by the different concentration of neutrons to protons. Different ratios of the neutron to proton concentrations lead to different critical points for the phase transition. This is analogous to the phase transitions occurring in 4He-3He liquid mixtures. We present experimental results which reveal the N/A (or Z/A) dependence of the phase transition and discuss possible implications of these observations in terms of the Landau Free Energy description of critical phenomena.
Normalizing flows are a class of machine learning models used to construct a complex distribution through a bijective mapping of a simple base distribution. We demonstrate that normalizing flows are particularly well suited as a Monte Carlo integration framework for quantum many-body calculations that require the repeated evaluation of high-dimensional integrals across smoothly varying integrands and integration regions. As an example, we consider the finite-temperature nuclear equation of state. An important advantage of normalizing flows is the ability to build highly expressive models of the target integrand, which we demonstrate enables precise evaluations of the nuclear free energy and its derivatives. Furthermore, we show that a normalizing flow model trained on one target integrand can be used to efficiently calculate related integrals when the temperature, density, or nuclear force is varied. This work will support future efforts to build microscopic equations of state for numerical simulations of supernovae and neutron star mergers that employ state-of-the-art nuclear forces and many-body methods.
We compute from chiral two- and three-nucleon interactions the energy per particle of symmetric nuclear matter and pure neutron matter at third-order in perturbation theory including self-consistent second-order single-particle energies. Particular attention is paid to the third-order particle-hole ring-diagram, which is often neglected in microscopic calculations of the equation of state. We provide semi-analytic expressions for the direct terms from central and tensor model-type interactions that are useful as theoretical benchmarks. We investigate uncertainties arising from the order-by-order convergence in both many-body perturbation theory and the chiral expansion. Including also variations in the resolution scale at which nuclear forces are resolved, we provide new error bands on the equation of state, the isospin-asymmetry energy, and its slope parameter. We find in particular that the inclusion of third-order diagrams reduces the theoretical uncertainty at low densities, while in general the largest error arises from omitted higher-order terms in the chiral expansion of the nuclear forces.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
