Code for the calculation of Pairing Vibration with phenomenological RPA
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Type Name Latest commit message Commit time
Failed to load latest commit information.


Code for the calculation of Pairing Vibration with phenomenological Random Phase Approximation (RPA). Eventually reads the Atomic Mass Evaluation Database for an automatic usage.

Code is not (and never will be) complete, releasing as process, not as product. For suggestions, contact me.

Andrea Idini, Univ. of Surrey

How to Run:

Use Makefile to compile. Change nucleus and proton/neutron vibration in

If you use the option 2, the program will scan the Atomic Mass Evaluation 2012 for masses, automatically selecting pairing correlation energy. Consequently, it fixes nuclear particle-hole gap. It calculates levels using the standard parametrization of Bohr and Mottelson for Wood Saxon + Coulomb, enjoy the simplicity.

If you use the option 1, you can specify levels and correlation energies following the format in the example.


Factor.dat and RPA.dat, for X and Y RPA amplitudes and energy of the phonons (in two formats). Dispersion.dat, dispersion relation of RPA.


If you don't have much idea what a Pairing Vibration and what a Random Phase Approximation (RPA) is, have a torough look at:

and refs. therein. If you use this code in a scientific work, please consider citing.


The Software is provided "as is" according to the GNUv3 License. I do not offer support. However, feel free to drop me an email with the scientific/didactic/business/life problem you want to discuss, I usually answer.

Read file LICENSE for more information on GNUv3.