GCM solver (Ver. 3.0): A mathematica notebook for diagonalization of the Geometric Collective Model (Bohr Hamiltonian) with generalized Gneuss-Greiner potential

Abstract

The program diagonalizes the Geometric Collective Model (Bohr Hamiltonian) with generalized Gneuss-Greiner potential with terms up to the sixth power in b. In nuclear physics, the Bohr-Mottelson model with later extensions into the rotovibrational Collective model is an important theoretical tool with predictive power and it represents a fundamental step in the education of a nuclear physicist. Nuclear spectroscopists might find it useful for fitting experimental data, reproducing spectra, EM transitions and moments and trying theoretical predictions, while students might find it useful for learning about connections between the nuclear shape and its quantum origin. Matrix elements for the kinetic energy operator and for scalar invariants as β2 and β3 cos (3γ) have been calculated in a truncated five-dimensional harmonic oscillator basis with a different program, checked with three different methods and stored in a matrix library for the lowest values of angular momentum. These matrices are called by the program that uses them to write generalized Hamiltonians as linear combinations of certain simple operators. Energy levels and eigenfunctions are obtained as outputs of the diagonalization of these Hamiltonian operators.