TABLE OF CONTENTS
- 1. /deuterium
- 1.1. deuterium/waveFunction_switch
- 1.2. deuterium/iParam
- 1.3. deuterium/pMax
- 1.4. deuterium/scaleMomentum
- 1.5. deuterium/deuteriumFermi
- 1.6. deuterium/initDeuterium
/deuterium [ Modules ]
NAME
module deuterium
PURPOSE
Includes subroutines to model the momentum distribution of nucleons in a Deuterium.
deuterium/waveFunction_switch [ Global module-variables ]
[ Top ] [ deuterium ] [ Global module-variables ]
SOURCE
integer, save :: waveFunction_switch = 1
PURPOSE
Possible values are:
- 0 -- No Wave functions! Pointlike Deuterium
- 1 -- Wave functions according to Bonn potential
- 2 -- Wave functions according to Argonne V18
deuterium/iParam [ Global module-variables ]
[ Top ] [ deuterium ] [ Global module-variables ]
SOURCE
integer, save :: iParam = 1
PURPOSE
Choose parameterization of momentum distribution when using the Bonn potential. Possible values:
- 1 -- Full Bonn (MaH87)
- 2 -- OBEPQ (MaH87)
- 3 -- OBEPQ-A (Mac89)
- 4 -- OBEPQ-B (Mac89)
- 5 -- OBEPQ-! (Mac89)
- 6 -- OBEPR (MaH87) self-made
- 7 -- Paris
References:
- MaH87: R. Machleidt et al. Phys. Rep. 149, 1 (1987)
- Mac89: R. Machleidt, Advances in Nucl. Phys. Vol 19
deuterium/pMax [ Global module-variables ]
[ Top ] [ deuterium ] [ Global module-variables ]
SOURCE
real,save :: pMax = 0.5
PURPOSE
Cut-off parameter for Fermi momentum
deuterium/scaleMomentum [ Global module-variables ]
[ Top ] [ deuterium ] [ Global module-variables ]
SOURCE
real, save :: scaleMomentum = 1.0
PURPOSE
The selected momentum is multiplied by this factor afterwards, i.e. some rescaling is done
deuterium/deuteriumFermi [ Namelists ]
[ Top ] [ deuterium ] [ Namelists ]
NAME
Namelist /deuteriumFermi/
PURPOSE
Includes the parameters:
which are used for the distribution of nucleons in deuterium.
deuterium/initDeuterium [ Subroutines ]
[ Top ] [ deuterium ] [ Subroutines ]
NAME
subroutine initDeuterium(teilchen,nuc,eventNummer,keepNumber)
PURPOSE
Represents nucleus 'nuc' in phase space by testparticles which are stored in vector 'teilchen'. The ordering in the vector teilchen is choosen to be random. The nucleus must be Deuterium.
NOTES
- momentum distribution implemented by Thomas Falter, spatial distribution added later
- The first nucleon is initialized at the origin, while the second nucleon gets its spatial coordinate according the squared wave function.
- The averaged deuterium radius is <r_d>~=2.0 fm. The variable r used here is the distance between the two nucleons. Therefore: <r_d> = 1/2 * sqrt(Int_0^Infty dr r^2 psi^2)
INPUTS
- type(nucleus) :: nuc
- type(particle),dimension(:,:) :: teilchen
- integer :: eventNummer -- eventNummer is given as "%event" to any initialized nucleon.
- logical :: keepNumber -- flag whether each testparticle gets unique number or not
OUTPUT
- type(particle),dimension(:,:) :: teilchen