Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNormalizing flows for microscopic many-body calculations: an application to the nuclear equation of state
115
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
117 -
Marcello Baldo
, Fiorella Burgio (Istituto Nazionale di Fisican Nucleare
, Sezione di Catania
2000
The Bethe-Brueckner-Goldstone many-body theory of the Nuclear Equation of State is reviewed in some details. In the theory, one performs an expansion in terms of the Brueckner two-body scattering matrix and an ordering of the corresponding many-body diagrams according to the number of their hole-lines. Recent results are reported, both for symmetric and for pure neutron matter, based on realistic two-nucleon interactions. It is shown that there is strong evidence of convergence in the expansion. Once three-body forces are introduced, the phenomenological saturation point is reproduced and the theory is applied to the study of neutron star properties. One finds that in the interior of neutron stars the onset of hyperons strongly softens the Nuclear Equation of State. As a consequence, the maximum mass of neutron stars turns out to be at the lower limit of the present phenomenological observation.
A new convenient method to diagonalize the non-relativistic many-body Schroedinger equation with two-body central potentials is derived. It combines kinematic rotations (democracy transformations) and exact calculation of overlap integrals between bases with different sets of mass-scaled Jacobi coordinates, thereby allowing for a great simplification of this formidable problem. We validate our method by obtaining a perfect correspondence with the exactly solvable three-body ($N=3$) Calogero model in 1D.
Starting from realistic nuclear forces, the chiral N$^3$LO and JISP16, we have applied many-body perturbation theory (MBPT) to the structure of closed-shell nuclei, $^4$He and $^{16}$O. The two-body N$^3$LO interaction is softened by a similarity renormalization group transformation while JISP16 is adopted without renormalization. The MBPT calculations are performed within the Hartree-Fock (HF) bases. The angular momentum coupled scheme is used, which can reduce the computational task. Corrections up to the third order in energy and up to the second order in radius are evaluated. Higher-order corrections in the HF basis are small relative to the leading-order perturbative result. Using the anti-symmetrized Goldstone diagram expansions of the wave function, we directly correct the one-body density for the calculation of the radius, rather than calculate corrections to the occupation propabilities of single-particle orbits as found in other treatments. We compare our results with other methods where available and find good agreement. This supports the conclusion that our methods produce reasonably converged results with these interactions. We also compare our results with experimental data.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
