Random Phase Approximation (RPA) is the basic method for calculation of excited states of nuclei over the Hartree-Fock ground state, suitable also for energy density functionals (EDF or DFT). We developed a convenient formalism for expressing densities and currents in a form of reduced matrix elements, which allows fast calculation of spectra for spherical nuclei. All terms of Skyrme functional were taken into account, so it is possible to calculate electric, magnetic and vortical/toroidal/compression transitions and strength functions of any multipolarity. Time-odd (spin) terms in Skyrme functional become important for magnetic and isovector toroidal (i.e. second-order term in long-wave expansion of E1) transitions. It was also found that transition currents in pygmy region (low-lying part of E1 resonance) exhibit isoscalar toroidal flow, so the previously assumed picture of neutron-skin vibration is not the only mechanism present in pygmy transitions.
RPA calculations with heavy axially-symmetric nuclei now become feasible on ordinary PC. Detailed formulae for axial Skyrme RPA are given. Some numerical results are shown in comparison with the approximate approach of separable RPA, previously developed in our group for fast calculation of strength functions.
Skyrme RPA for spherical and axially symmetric nuclei
Anton Repko, Jan Kvasil, V.O. Nesterenko and P.-G. Reinhard
Institute of Particle and Nuclear Physics, Charles University, V Holešovičkách 2, 180 00 Prague 8, Czech Republic
Joint Institute for Nuclear Research, Dubna, Moscow region 141980, Russia
Institut für Theoretische Physik II, Universität Erlangen, D-91058 Erlangen, Germany
Energy density functional in nuclear physics is a self-consistent microscopic phenomenological approach to calculate nuclear properties and structure over the whole periodic table. The method is analogous to Kohn-Sham density functional theory (DFT) used in electronic systems. Three types of functionals are frequently used nowadays: non-relativistic Skyrme functional  with zero-range two-body and density dependent interaction, finite-range Gogny force  and relativistic (covariant) mean-field . Typical approach employs Hartree-Fock-Bogoliubov or HF+BCS calculation scheme to obtain ground state and single-(quasi)particle wavefunctions and energies. These results are then utilized to fit the parameters of the functionals, thus obtaining various parametrizations suitable for specific aims, such as: calculation of mass-table, charge radii, fission barriers, spin-orbit splitting and giant resonances.
Random Phase Approximation (RPA) is a textbook standard  to calculate one-phonon excitations of the nucleus. Increasing computing power has enabled to employ fully self-consistent residual interaction derived from the same density functional as the underlying ground state. While the spherical nuclei can be treated directly (by matrix diagonalization) [5, 6, 7], axially deformed nuclei still pose certain difficulties due to large matrix dimensions . Our group developed a separable RPA (SRPA) approach [9, 10], which greatly reduces the computational cost by utilizing separable residual interaction, entirely derived from the underlying functional by means of multi-dimensional linear response theory.
Present article gives a convenient formalism for rotationally-invariant treatment of spherical full RPA, and also gives detailed expressions for matrix elements in axial symmetry. Both time-even and time-odd terms of Skyrme functional are employed, so the method is suitable for various electric and magnetic multipolarities. Finally, some results are shown by means of strength functions and transition currents.
2 Skyrme functional and density operators
Skyrme interaction is defined as
with parameters and a spin-exchange operator . Since it is a zero-range interaction, the solution of many-body problem by Hartree-Fock method can be equivalently reformulated as a density functional theory , where the complete density functional is
where the last line contains the spin (time-odd) terms, which are usually omitted. However, they have quite important contribution for magnetic excitations , as will be illustrated below, so we use them in all calculations. Parameters depend on the parameters from (1):
The operators corresponding to the densities and currents in (2) are defined as
and they are understood as single-particle operators in a many-body system. The densities are time-even and the currents are time-odd. Spin-orbital current and current have two indices, so we will decompose them into scalar, vector and (symmetric) tensor part
The coefficients related to the (non-relativistic) time-reversal symmetry of the operators are defined as
where factors define the Bogoliubov transformation of particles to quasiparticles . Pairing interaction in the particle-particle channel is defined as ( is the nucleon ground-state density)
with constant parameters and . Contribution to the density functional is then
and the pairing density with the corresponding operator is
with energy-cutoff weight . We use monopole pairing, and we involve only diagonal pairs.
3 RPA formalism in spherical symmetry
Notation of the Clebsch-Gordan coefficients comes from the book of Varshalovich . In the language of spherical tensors, a hermitian operator satisfies
for example a vector (rank-1) operator:
where the symbol in bold font represents the (complex-conjugated) spherical harmonics in its scalar, vector or tensor form (depending on the rank of operator : ):
Expression (14) defines the reduced matrix elements , e.g., a standard density matrix element:
and radial wavefunctions are taken from
with spinors . Remaining reduced matrix elements will be published later.
Excitations of a given multipolarity are treated as the RPA phonons. One-phonon state, , is created by action of an operator on the RPA ground state , and has an excitation energy .
The operator is a two-quasiparticle () operator defined by real coefficients
and satisfies the RPA equation
where the index means that we take only the two-quasiparticle portion of the commutator (after normal ordering). Although all commutators should be evaluated in the RPA ground state, it is common to evaluate them in the HF+BCS ground state (i.e., we use quasi-boson approximation), since the contribution of and higher correlations in the ground state is assumed to be low.
and will be rescaled automatically. Diagonal matrix elements contribute only to electric transitions with even.
The hamiltonian is taken as a sum of mean-field part (HF+BCS) and the second functional derivative of Skyrme density functional (+ Coulomb and pairing interaction).
RPA equation (20) then turns into a matrix equation
which can be recast into half-dimensional symmetric-matrix eigenvalue problem . Matrices and are
The expression is symbolical, and includes the integration of the delta function, yielding . The exchange Coulomb interaction can be treated by Slater approximation as a density functional
however, the direct Coulomb interaction gives rise to a double integral instead
After calculation of the RPA states, yielding and , we are interested in the matrix elements of electric and magnetic transition operators and in the transition densities and currents.
Besides electric () and magnetic () operators in long-wave approximation, we implement also electric vortical, toroidal and compression operators  (toroidal operator constitutes the next-order term in long-wave expansion of the exact electric transition operator)
where are effective charges of the nucleons, are orbital/spin g-factors (we take ; spin g-factors are reduced by a quenching factor ) and is the nuclear current composed of convective and magnetization part
where (convective) current and spin one-body operators are the same as in Skyrme functional (4).
4 RPA formalism in axial symmetry
Axial coordinates are
Calculations in axially deformed nuclei don’t conserve total angular momentum (previously denoted or ), nevertheless, they conserve its -projection () and parity, so it is convenient to preserve part of the formalism from spherical symmetry, namely the convention of -components for vector and tensor operators, and their hermitian conjugation (12). Operators of differentiation are then
Single-particle wavefunction (and its time-reversal conjugate) is expressed as a spinor
and the radial parts of its derivatives will be denoted by a shorthand notation
Radial functions are real, and their spinor-wise products will be denoted by a dot to keep the expressions simple:
Vector currents will be decomposed in the style of rank-1 tensor operators. Vector product in the expression for spin-orbital current leads to (for vector product in -scheme see [15, (1.2.28)])
Matrix elements of densities and currents are then
|Factor will be omitted in the following expressions.|
In the actual calculation, it is necessary to choose projection of angular momentum and parity (together denoted also as , and ). Transition operators have the form of
where contains a function (or even derivatives) not dependent on . Choice of the two-quasiparticle pairs is restricted by (time-reversed states are omitted from the discussion here for simplicity; in short, it is necessary to add pairs with during the duplicate removal similar to (21)). Commutators are then evaluated in quasiparticle vacuum as