Implementation of the quasiparticle finite amplitude method within the relativistic self-consistent mean-field framework (II): The program DIRQFAM v2.0.0

We describe the new version 2.0.0 of the code DIRQFAM that calculates the multipole response of even-even axially symmetric deformed nuclei using the framework of relativistic self-consistent mean-field models. The response is calculated by implementing the finite amplitude method for relativistic quasiparticle random phase approximation. In the new version, we have implemented the following features: (i) meson-exchange interactions, (ii) low-rank method for calculating induced densities and currents, (iii) generalized minimal residual method (GMRES), (iv) evaluation of the nucleon localization function, (v) extraction of the QRPA transition matrix elements and eigenfrequencies, (vi) evaluation of the H11 part of the induced Hamiltonian, (vii) multipole excitation operators up to J <= 5 are now included.Program summary Program title: DIRQFAM v2.0.0 CPC Library link to program files: https://doi .org /10 .17632 /6gv8nxg4ns .2 Licensing provisions: GPLv3 Programming language: Fortran 90. External routines/libraries: BLAS/LAPACK, version 3.1. or higher. Journal reference of previous version: A. Bjelcic, T. Niksic, Comp. Phys. Comm. 253, 107184 (2020). Does the new version supersede the previous version?: Yes Nature of problem: Multipole response of deformed even-even open-shell nuclei can be calculated using the quasiparticle finite amplitude method (QFAM), based on the relativistic self-consistent mean-field models. The particle-hole channel is described by a zero-range relativistic effective interaction or interaction based on meson exchange, while the particle-particle channel of the effective inter-nucleon interaction is described by a separable finite-range pairing force. The method can be applied to perform systematic studies of collective modes even in heavy deformed nuclei. Solution method: The current implementation computes the multipole response by solving the QFAM equations in a self-consistent iteration scheme. At each iteration the QFAM solutions are updated using the GMRES method. The QFAM amplitudes are expanded in the simplex-y harmonic oscillator basis. Summary of revisions:1. meson-exchange interaction are included in the code; 2. calculation of the nucleon localization function; 3. extraction of the QRPA transition matrix elements and eigenfrequencies; 4. evaluation of the H11 part of the induced Hamiltonian; 5. multipoles excitation operators up to J <= 5 are included; 6. low-rank method for calculating induced densities and currents; 7. generalized minimal residual method; 8. dynamic memory allocation; 9. reduced memory requirements;Additional comments including restrictions and unusual features: Open-shell even-even nuclei with parity conserved axially symmetric ground states are considered. An electric multipole operator with J <= 5 is used to calculate the response function.(c) 2023 Elsevier B.V. All rights reserved.

Automated summary

We present version 2.0.0 of the DIRQFAM code, which calculates the multipole response of deformed nuclei using the relativistic self-consistent mean-field models framework. The response is computed using the finite amplitude method for relativistic quasiparticle random phase approximation. In this new version, various additional features have been implemented, such as meson-exchange interactions and evaluation of the nucleon localization function.