Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMicroscopically-based energy density functionals for nuclei using the density matrix expansion: Full optimization and validation
351
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We seek to obtain a usable form of the nuclear energy density functional that is rooted in the modern theory of nuclear forces. We thus consider a functional obtained from the density matrix expansion of local nuclear potentials from chiral effective field theory. We propose a parametrization of this functional carefully calibrated and validated on selected ground-state properties that is suitable for large-scale calculations of nuclear properties. The first component of this functional is a non-local functional of the density and corresponds to the direct part (Hartree term) of the expectation value of local chiral potentials on a Slater determinant. A second component is a local functional of the density and is obtained by applying the density matrix expansion to the exchange part (Fock term) of the expectation value of the local chiral potential. We apply the UNEDF2 optimization protocol to determine the coupling constants of this energy functional. We obtain a set of microscopically-constrained functionals for local chiral potentials from leading-order up to next-to-next-to-leading order with and without three-body forces and contributions from $Delta$ excitations. These functionals are validated on the calculation of nuclear and neutron matter, nuclear mass tables, single-particle shell structure in closed-shell nuclei and the fission barrier of $^{240}$Pu. Quantitatively, they perform noticeable better than the more phenomenological Skyrme functionals. The inclusion of higher-order terms in the chiral perturbation expansion seems to produce a systematic improvement in predicting nuclear binding energies. This result is especially promising since all the fits have been performed at the single reference level of the energy density functional approach, where important collective correlations such as center-of-mass correction have not been taken into account yet.
rate research
Read More
While ab initio many-body techniques have been able to successfully describe the properties of light and intermediate mass nuclei based on chiral effective field theory interactions, neutron-rich nuclei still remain out of reach for these methods. Conversely, energy density functional approaches can be used to calculate properties of heavy nuclei but rely mostly on phenomenological interactions. A usable form of the nuclear energy density functional that is rooted in the modern theory of nuclear forces was presented recently. The first component of this new set of functionals corresponds to the direct part (Hartree term) of the expectation value of local chiral potentials on a Slater determinant. The exchange term, which is a functional of the non-local density, is transformed into a local functional by applying the density matrix expansion. In order to reduce the computational cost due to the direct implementation of non-separable, local interactions in the Hartree term, we use an approximation to represent the regularized Yukawa functions in terms of a sum of (separable) Gaussian functions. These proceedings analyze the accuracy of such an approximation in terms of the number of Gaussian functions and look for an optimal value that gives an acceptable level of accuracy while maintaining the computational memory requirements in a many-body calculation as low as possible.
Nuclear density functional theory is the prevalent theoretical framework for accurately describing nuclear properties at the scale of the entire chart of nuclides. Given an energy functional and a many-body scheme (e.g., single- or multireference level), the predictive power of the theory depends strongly on how the parameters of the energy functionals have been calibrated with experimental data. Expanded algorithms and computing power have enabled recent optimization protocols to include data in deformed nuclei in order to optimize the coupling constants of the energy functional. The primary motivation of this work is to test the robustness of such protocols with respect to some of the technical and numerical details of the underlying calculations, especially when the calibration explores a large parameter space. To this end, we quantify the effect of these uncertainties on both the optimization and statistical emulation of composite objective functions. We also emphasize that Bayesian calibration can provide better estimates of the theoretical errors used to define objective functions.
Relativistic energy density functionals have become a standard framework for nuclear structure studies of ground-state properties and collective excitations over the entire nuclide chart. We review recent developments in modeling nuclear weak-interaction processes: charge-exchange excitations and the role of isoscalar proton-neutron pairing, charged-current neutrino-nucleus reactions relevant for supernova evolution and neutrino detectors, and calculation of beta-decay rates for r-process nucleosynthesis.
Quadrupole and octupole deformation energy surfaces, low-energy excitation spectra and transition rates in fourteen isotopic chains: Xe, Ba, Ce, Nd, Sm, Gd, Rn, Ra, Th, U, Pu, Cm, Cf, and Fm, are systematically analyzed using a theoretical framework based on a quadrupole-octupole collective Hamiltonian (QOCH), with parameters determined by constrained reflection-asymmetric and axially-symmetric relativistic mean-field calculations. The microscopic QOCH model based on the PC-PK1 energy density functional and $delta$-interaction pairing is shown to accurately describe the empirical trend of low-energy quadrupole and octupole collective states, and predicted spectroscopic properties are consistent with recent microscopic calculations based on both relativistic and non-relativistic energy density functionals. Low-energy negative-parity bands, average octupole deformations, and transition rates show evidence for octupole collectivity in both mass regions, for which a microscopic mechanism is discussed in terms of evolution of single-nucleon orbitals with deformation.
Machine learning is employed to build an energy density functional for self-bound nuclear systems for the first time. By learning the kinetic energy as a functional of the nucleon density alone, a robust and accurate orbital-free density functional for nuclei is established. Self-consistent calculations that bypass the Kohn-Sham equations provide the ground-state densities, total energies, and root-mean-square radii with a high accuracy in comparison with the Kohn-Sham solutions. No existing orbital-free density functional theory comes close to this performance for nuclei. Therefore, it provides a new promising way for future developments of nuclear energy density functionals for the whole nuclear chart.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
