Background: The quasiparticle random-phase approximation (QRPA) within nuclear density functional theory provides a powerful framework for describing collective excitations. Although the finite-amplitude method (FAM) efficiently solves the QRPA problem, repeated calculations for different external-field parameters remain computationally demanding. Purpose: We construct a reduced basis method (RBM)-based emulator for the FAM that treats the complex energy of the external field as a model parameter to efficiently reproduce FAM amplitudes and QRPA eigenmodes. Methods: High-fidelity FAM calculations are performed at a small set of training points in the complex-energy plane. The resulting FAM amplitudes form a non-orthogonal reduced basis. A variational equation yields an emulator that can predict the response at arbitrary complex energies and QRPA eigensolutions without additional full FAM calculations. Results: The RBM emulator accurately reproduces FAM strength distributions in both giant-resonance and low-energy regions when the relevant energy domain is covered by the training set. It also reproduces imaginary QRPA modes associated with shape instabilities of the HFB state. Applied to the
Kπ=0+ mode of rare-earth Dy isotopes in a realistic model space, the emulator reproduces strength distributions and the lowest
0+ collective states with precision comparable to full FAM calculations, reducing the computational cost by more than an order of magnitude. Conclusions: The RBM provides an efficient and accurate FAM emulator. Its ability to reproduce giant-resonance, low-energy, and imaginary-energy modes at drastically reduced computational cost makes it promising for density-functional optimization, calculations of collective inertia, and large-scale surveys of nuclear collective excitations.
A nonlinear generalization of the quasiparticle random phase approximation is constructed with the aim of describing anharmonic effects in the low-lying spectra of vibrational nuclei. Angular momentum zero pairing correlations of the ground state are described by a generalized BCS calculation based on a variational method. The dynamics of the vibrational states are described by the Heisenberg equations of motion for fermion pair operators; matrix elements of products of pair operators are treated by sum rule methods. Kinematical constraints are expressed by the algebra of the pair operators and by boson-like polynomial expansions for the quasiparticle conserving multipole operators in terms of the quasiparticle pair changing operators, these expansions enforcing the content of the Pauli principle. The physical basis for the choice of approximations is described. The final algorithm is organized so as to involve little more than repeated BCS-like and RPA-like calculations, converging to a self-consistent solution of the full nonlinear scheme. Detailed discussion is carried out for a modified form of the pairing plus quadrupole interaction Hamiltonian, but the method is not confined to this model.NUCLEAR STRUCTURE Approximate many-body methods applied to modified pairing plus quadrupole model. Description of vibrational spectrum of single closed shell nuclei. Equations of motion and algebraic methods.
For a description of large-amplitude collective motion associated with nuclear pairing, requantization of time-dependent mean-field dynamics is performed using the stationary-phase approximation (SPA) to the path integral. We overcome the difficulty of the SPA, which is known to be applicable to integrable systems only, by developing a requantization approach combining the SPA with the adiabatic self-consistent collective coordinate method (ASCC+SPA). We apply the ASCC+SPA to multilevel pairing models, which are nonintegrable systems, to study the nuclear pairing dynamics. The ASCC+SPA gives a reasonable description of low-lying excited
0+ states in nonintegrable pairing systems.