期刊
COMPUTER PHYSICS COMMUNICATIONS
卷 253, 期 -, 页码 -出版社
ELSEVIER
DOI: 10.1016/j.cpc.2020.107184
关键词
Relativistic quasiparticle random-phase approximation; Finite amplitude method; Multipole response; Dirac-Hartree-Bogoliubov; Nuclear energy density functional
资金
- Croatian Science Foundation [IP-2014-09-9159]
- QuantiXLie Centre of Excellence
- Croatian Government
- European Union through the European Regional Development Fund -the Competitiveness and Cohesion Operational Programme [KK.01.1.1.01]
The DIRQFAM code 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. Program summary Program Title: DIRQFAM Program Files doi: http://dx.doi.org/10.17632/6gv8mcg4ns.1 Licensing provisions: GPLv3 Programming language: Fortran 90/95, easily downgradable to FORTRAN 77. External routines/libraries: BLAS/LAPACK, version 3.6.0. or higher. 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 meanfield models. The particle-hole channel is described by a zero-range relativistic effective interaction, 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 modified Broyden's method. The QFAM amplitudes are expanded in the simplex-y harmonic oscillator basis. Restrictions: Open-shell even-even nuclei with axially symmetric ground states are considered. An electric multipole operator with J <= 3 is used to calculate the response function. Unusual features: The code should be recompiled before running it with different values of one of the following input parameters: number of oscillator shells in the expansion of nucleon spinors (n0f ), number of Gauss-Hermite (NGH) or Gauss-Laguerre nodes (NGL), the multipolarity values that define the external perturbation operator (J_multipole, K_multipole). (C) 2020 Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据