Install

KMFR was originally written in Fortran targeting the Compaq Visual Fortran compiler for Windows and relied on some of its non-standard extensions. Since the CVF was last released in 1999 and is no longer supported, in 2019 the model was ported to pure Fortran 2003. Testing has been performed using GFortran 8 and 9 but any recent compiler should work.

System requirements

Fortran 2003 compiler. GCC gfortran 8 & 9 tested.
A Unix-like operating system. KMFR has minor dependencies on Unix conventions such as directory path syntax and shell escape commands. Mac OS X and GNU/Linux are known to work.
CERNLIB, a once-widespread mathematical library. Specifically, the MATHLIB and KERNLIB parts of it are used. On Debian-based Linux distributions (including Ubuntu) it is enough to install the cernlib package:
```
# apt-get install cernlib
```
Unfortunately, Red Hat family of distributions (CentOS, Fedora, Scientific Linux) stopped supporting CERNLIB packages. Generic instructions on how to install CERNLIB on any GNU/Linux can be fond on the DESY website.
On Mac OS X CERNLIB can be conveniently installed using the Fink package manager by executing a shell command
```
$ sudo apt-get install cernlib2006-dev 
```

Build instructions

Uncompress the distribution file and cd to the kmfr-VERSION directory.
On Linux execute the command
```
$ make -f makefile.linux
```
On Mac OS X execute the command
```
$ make -f makefile.mac
```
Note: certain edits to the makefile might be necessary to set up paths to the Fortran compiler and libraries.

Usage

Configuration parameters of the model are specified in a text configuration file using the Fortran namelist format. A sample configuration file with comments is provided (sample-kmfr.cfg). Refer to it for user-modifiable settings.

The calculation is started by running either

$ /path/to/kmfr sample-kmfr.cfg

$ /path/to/kmfr < sample-kmfr.cfg

Calculation results can be found in the automatically created output directory.

Calculated cross sections can be found in the sct1.dat, sct2.dat,... files, where the numeric index corresponds to the reaction code in the text file key. For each reaction this file lists, respectively, the reaction code, its integrated cross section, threshold, and the number of emitted neutrons and protons. All energy values are specified in [MeV] and cross sections in [mb] units.

Sample configuration file

!A sample configuration file for KMFR.
!Fortran namelist format is used.
!Comments are denoted by '!'.
&kmfr_config
!These are the only mandatory parameters.
  element = 'Sm',
  mass_number = 154,
  deformation = 0.319, ! quadrupole deformation parameter β₂ (in ground state)
  ! The provided file 'data/deform.dat' has a reference
  ! also see http://cdfe.sinp.msu.ru/services/radchart/radmain.html
  ! By controlling the amount of deformation splitting of the GDR and GQR
  ! this parameter affects the width of the photoabsorption cross section 

!Other parameters are optional and can be left in their default values
  density_model = 3, ! nuclear level density model
  ! 1 == Back-shifted Fermi gas Model (BFM)
  ! 2 == Generalized Superfluid Model (GSM)
  ! 3 == (recommended) Gilbert-Cameron Constant Temperature Model (CTM)
  isospin_split_model = 1, ! isospin splitting of the GDR
  ! 1 == use the KMFR model (works better for heavy nuclei)
  ! 2 == use the Fallieros expressions (likewise for medium nuclei)
  qd_model = 1, ! effect of correlation in the quasideuteron doorway state
  ! (affects resulting cross sections at higher energies)
  ! 0 == no correlation, damping starts from 2p2h
  ! 1 == energies of particles and holes are correlated, start with 2p1h
  energy_max = 44, ! maximum photon energy (must be <= 44 MeV in this version)
  use_delay = .FALSE., ! effect of delayed decay of the doorway state
  ! generally results in a larger yield of (gamma,1n)
  ! should not be used for soft-surface nuclei
  use_dcm_model = .FALSE., ! use the DCM model of photoabsorption
  ! use the Dynamic Collective Model instead of the KMFR photoabsorption model
  dcm_hw2 = 8, ! DCM parameters should be set by hand
  dcm_bt0 = 0.2,
  use_rotator_model = .FALSE., ! use the rigid rotator model
  ! use the rigid rotator model of photoabsorption instead
  rot_hw2 = 8, ! rotator parameters should be set by hand
  rot_bt0 = 0.2,
  rot_gm0 = 0.1,
  dynamic_deformation = 0.1,
  gqr_nu = 0.0, ! enhancement parameter for the GQR strength
  irb0 = 2, ! this is numbering one of 1.5, 1.55, 1.60, 1.65, 1.70
  ! this is the parameterization used for calculation of the 
  ! damping width of the GDR
  ! recommended value: 2
  ! generally, irb0 should be 1..5 with smaller values for light nuclei,
  ! and larger values for heavy nuclei
  exciton_model_GM = 160, ! parameter of the squared matrix element
  ! of intranuclear transitions of the exciton model
  ! the larger its value, the faster the damping is and the more
  ! statistic are the cross sections
  kw = 0.1, ! smearing parameter
  drw = .TRUE., ! respect collective structure of the GR doorway state
  gdr_energy_shift = 0.1, ! arbitrary energy shift
  meson_exchange_alpha = 0.3 ! meson exchange currents correction
  ! this parameter can be used to tune the photoabsorption cross section
/