No Arabic abstract
We propose a new method to solve the eigen-value problem with a two-center single-particle potential. This method combines the usual matrix diagonalization with the method of separable representation of a two-center potential, that is, an expansion of the two-center potential with a finite basis set. To this end, we expand the potential on a harmonic oscillator basis, while single-particle wave functions on a combined basis with a harmonic oscillator and eigen-functions of a one-dimensional two-center potential. In order to demonstrate its efficiency, we apply this method to a system with two $^{16}$O nuclei, in which the potential is given as a sum of two Woods-Saxon potentials.
Shell corrections of the finite deformed Woods-Saxon potential are calculated using the Greens function method and the generalized Strutinsky smoothing procedure. They are compared with the results of the standard prescription which are affected by the spurious contribution from the unphysical particle gas. In the new method, the shell correction approaches the exact limit provided that the dimension of the single-particle (harmonic oscillator) basis is sufficiently large. For spherical potentials, the present method is faster than the exact one in which the contribution from the particle continuum states is explicitly calculated. For deformed potentials, the Greens function method offers a practical and reliable way of calculating shell corrections for weakly bound nuclei.
We apply a new calculation scheme of a finite element method (FEM) for solving an elliptic boundary-value problem describing a quadrupole vibration collective nuclear model with tetrahedral symmetry. We use of shape functions constructed with interpolation Lagrange polynomials on a triangle finite element grid and compare the FEM results with obtained early by a finite difference method.
We present an approach to derive effective shell-model interactions from microscopic nuclear forces. The similarity-transformed coupled-cluster Hamiltonian decouples the single-reference state of a closed-shell nucleus and provides us with a core for the shell model. We use a second similarity transformation to decouple a shell-model space from the excluded space. We show that the three-body terms induced by both similarity transformations are crucial for an accurate computation of ground and excited states. As a proof of principle we use a nucleon-nucleon interaction from chiral effective field theory, employ a $^4$He core, and compute low-lying states of $^{6-8}$He and $^{6-8}$Li in $p$-shell model spaces. Our results agree with benchmarks from full configuration interaction.
We present the code HF-SHELL for solving the self-consistent mean-field equations for configuration-interaction shell model Hamiltonians in the proton-neutron formalism. The code can calculate both ground-state and finite-temperature properties in the Hartree-Fock (HF), HF+Bardeen-Cooper-Schrieffer (HF+BCS), and the Hartree-Fock-Bogoliubov (HFB) mean-field approximations. Particle-number projection after variation is incorporated to reduce the grand-canonical ensemble to the canonical ensemble, making the code particularly suitable for the calculation of nuclear state densities. The code does not impose axial symmetry and allows for triaxial quadrupole deformations. The self-consistency cycle is particularly robust through the use of the heavy-ball optimization technique and the implementation of different options to constrain the quadrupole degrees of freedom.
Fourth-order differential equations play an important role in many applications in science and engineering. In this paper, we present a three-field mixed finite-element formulation for fourth-order problems, with a focus on the effective treatment of the different boundary conditions that arise naturally in a variational formulation. Our formulation is based on introducing the gradient of the solution as an explicit variable, constrained using a Lagrange multiplier. The essential boundary conditions are enforced weakly, using Nitsches method where required. As a result, the problem is rewritten as a saddle-point system, requiring analysis of the resulting finite-element discretization and the construction of optimal linear solvers. Here, we discuss the analysis of the well-posedness and accuracy of the finite-element formulation. Moreover, we develop monolithic multigrid solvers for the resulting linear systems. Two and three-dimensional numerical results are presented to demonstrate the accuracy of the discretization and efficiency of the multigrid solvers proposed.