اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدContinuity equation and local gauge invariance for the N3LO nuclear Energy Density Functionals
37
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
Background: The next-to-next-to-next-to-leading order (N3LO) nuclear energy density functional extends the standard Skyrme functional with new terms depending on higher-order derivatives of densities, introduced to gain better precision in the nuclear many-body calculations. A thorough study of the transformation properties of the functional with respect to different symmetries is required, as a step preliminary to the adjustment of the coupling constants. Purpose: Determine to which extent the presence of higher-order derivatives in the functional can be compatible with the continuity equation. In particular, to study the relations between the validity of the continuity equation and invariance of the functional under gauge transformations. Methods: Derive conditions for the validity of the continuity equation in the framework of time-dependent density functional theory. The conditions apply separately to the four spin-isospin channels of the one-body density matrix. Results: We obtained four sets of constraints on the coupling constants of the N3LO energy density functional that guarantee the validity of the continuity equation in all spin-isospin channels. In particular, for the scalar-isoscalar channel, the constraints are the same as those resulting from imposing the standard U(1) local-gauge-invariance conditions. Conclusions: Validity of the continuity equation in the four spin-isospin channels is equivalent to the local-gauge invariance of the energy density functional. For vector and isovector channels, such validity requires the invariance of the functional under local rotations in the spin and isospin spaces.
قيم البحث
اقرأ أيضاً
We introduce a finite-range pseudopotential built as an expansion in derivatives up to next-to-next-to-next-to-leading order (N$^3$LO) and we calculate the corresponding nonlocal energy density functional (EDF). The coupling constants of the nonlocal
EDF, for both finite nuclei and infinite nuclear matter, are expressed through the parameters of the pseudopotential. All central, spin-orbit, and tensor terms of the pseudopotential are derived both in the spherical-tensor and Cartesian representation. At next-to-leading order (NLO), we also derive relations between the nonlocal EDF expressed in the spherical-tensor and Cartesian formalism. Finally, a simplified version of the finite-range pseudopotential is considered, which generates the EDF identical to that generated by a local potential.
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 f
or 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.
The Coulomb exchange and correlation energy density functionals for electron systems are applied to nuclear systems. It is found that the exchange functionals in the generalized gradient approximation provide agreements with the exact-Fock energy wit
h one adjustable parameter within a few dozen $ mathrm{keV} $ accuracy, whereas the correlation functionals are not directly applicable to nuclear systems due to the existence of the nuclear force.
We address the question of how to improve the agreement between theoretical nuclear single-particle energies (SPEs) and experiment. Empirically, in doubly magic nuclei, the SPEs can be deduced from spectroscopic properties of odd nuclei that have one
more, or one less neutron or proton. Theoretically, bare SPEs, before being confronted with experiment, must be corrected for the effects of the particle-vibration-coupling (PVC). In the present work, we determine the PVC corrections in a fully self-consistent way. Then, we adjust the SPEs, with PVC corrections included, to empirical data. In this way, the agreement with experiment, on average, improves; nevertheless, large discrepancies still remain. We conclude that the main source of disagreement is still in the underlying mean fields, and not in including or neglecting the PVC corrections.
It is known that some well-established parametrizations of the EDF do not always provide converged results for nuclei and a qualitative link between this finding and the appearance of finite-size instabilities of SNM near saturation density when comp
uted within the RPA has been pointed out. We seek for a quantitative and systematic connection between the impossibility to converge self-consistent calculations of nuclei and the occurrence of finite-size instabilities in SNM for the example of scalar-isovector (S=0, T=1) instabilities of the standard Skyrme EDF. We aim to establish a stability criterion based on computationally-friendly RPA calculations of SNM that is independent on the functional form of the EDF and that can be utilized during the adjustment of its coupling constants. Tuning the coupling constant $C^{rho Deltarho}_{1}$ of the gradient term that triggers scalar-isovector instabilities of the standard Skyrme EDF, we find that the occurrence of instabilities in finite nuclei depends strongly on the numerical scheme used to solve the self-consistent mean-field equations. The link to instabilities of SNM is made by extracting the lowest density $rho_{text{crit}}$ at which a pole appears at zero energy in the RPA response function when employing the critical value of the coupling constant $C^{rho Deltarho}_{1}$ extracted in nuclei. Our analysis suggests a two-fold stability criterion to avoid scalar-isovector instabilities: (i) The density $rho_{text{min}}$ corresponding to the lowest pole in the RPA response function should be larger than about 1.2 times the saturation density; (ii) one needs to verify that $rho_{p}(q_{text{pq}})$ exhibits a distinct global minimum and is not a decreasing function for large transferred momenta.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
