TABLE OF CONTENTS
- 1. /baryonPotentialModule
- 1.1. baryonPotentialModule/EQS_Type
- 1.2. baryonPotentialModule/dsymm
- 1.3. baryonPotentialModule/noPerturbativePotential
- 1.4. baryonPotentialModule/DeltaPot
- 1.5. baryonPotentialModule/HypPot
- 1.6. baryonPotentialModule/symmetriePotFlag
- 1.7. baryonPotentialModule/SurfacePotFlag
- 1.8. baryonPotentialModule/rho_0
- 1.9. baryonPotentialModule/p_0
- 1.10. baryonPotentialModule/U_0
- 1.11. baryonPotentialModule/bindingEnergy
- 1.12. baryonPotentialModule/compressibility
- 1.13. baryonPotentialModule/nLoopReAdjust
- 1.14. baryonPotentialModule/baryonPotential
- 1.15. baryonPotentialModule/BaryonPotential
- 1.16. baryonPotentialModule/variableSkyrme
- 1.17. baryonPotentialModule/momentumDependentPart
- 1.18. baryonPotentialModule/rhoLaplace
- 1.19. baryonPotentialModule/LDApotential(teilchen)
- 1.20. baryonPotentialModule/LDApotentialWelke(teilchen)
- 1.21. baryonPotentialModule/SurfacePart
- 1.22. baryonPotentialModule/rearrangementPotential
- 1.23. baryonPotentialModule/HandPotentialToDensityStatic
- 1.24. baryonPotentialModule/getNoPertPot_baryon
- 1.25. baryonPotentialModule/getsymmetriePotFlag_baryon
- 1.26. baryonPotentialModule/getPotentialEQSType
/baryonPotentialModule [ Modules ]
NAME
module baryonPotentialModule
PURPOSE
Includes all information about the baryonic potentials.
baryonPotentialModule/EQS_Type [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Switch for equation of state for nucleon resonances with spin 1/2. Parameters for nucleon potentials:
- 0 = nucleon potential is set to zero
- 1 = soft, momentum dependent, lambda = 2.130
- 2 = hard, momentum dependent, lambda = 2.126
- 3 = soft, momentum independent
- 4 = hard, momentum independent
- 5 = medium, momentum dependent, lambda = 2.130
- 6 = LDA potential (Birger Steinmueller)
- 7 = Deuterium potential Argonne V18 (only for Deuterium)
- 8 = LDA Potential Welke
- 9 = Buss PhD, Set#1
- 10 = Buss PhD, Set#2
- 11 = Buss PhD, Set#3
- 98 = use pre-stored values
- 99 = variable Skyrme : E_bind, p_0, U_O, rho_0 must be defined!
NOTES
Can be set in namelist baryonPotential.
SOURCE
integer,save :: EQS_Type = 5
baryonPotentialModule/dsymm [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
SOURCE
real,save :: dsymm = 0.03
PURPOSE
Parameter for symmetriePotential in GeV.
NOTES
Can be set in namelist baryonPotential.
baryonPotentialModule/noPerturbativePotential [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Switch for potential of perturbative particles. If .true. then perturbative baryons feel no potential.
SOURCE
logical,save :: noPerturbativePotential=.false.
NOTES
Can be set in namelist baryonPotential.
baryonPotentialModule/DeltaPot [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Switch for potential of spin=3/2 resonances:
- 0 = no potential
- 1 = nucleon (spin=1/2) potential times 2/3 [according to Ericson/Weise book]
- 2 = 100 MeV * rho/rhoNull
- 3 = nucleon (spin=1/2) potential
SOURCE
integer, save :: DeltaPot = 1
NOTES
Can be set to different value in jobcard, namelist 'baryonPotential'.
baryonPotentialModule/HypPot [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Switch for potential of hyperons:
- 0 = no potential
- 1 = nucleon (spin=1/2) potential times (3+S)/3 (i.e. according to the share of the light quarks)
- 2 = nucleon (spin=1/2) potential
SOURCE
integer, save :: HypPot = 1
NOTES
Can be set to different value in jobcard, namelist 'baryonPotential'
baryonPotentialModule/symmetriePotFlag [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Switch for the asymmetry term in the nucleon potential.
SOURCE
logical, save :: symmetriePotFlag = .false.
NOTES
Can be set to different value in jobcard, namelist 'baryonPotential'. Needs to be turned on together with 'coulombFlag' (in namelist 'coulomb').
baryonPotentialModule/SurfacePotFlag [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Switch for the surface term in the nucleon potential.
SOURCE
logical, save :: SurfacePotFlag=.false.
NOTES
- Do not use it together with yukawa!
- Can be set to different value in jobcard, namelist 'baryonPotential'.
baryonPotentialModule/rho_0 [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Nuclear matter density for EQS-Type=99
SOURCE
real, save :: rho_0=0.16
NOTES
- Units : fm^{-3}
baryonPotentialModule/p_0 [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
momentum for which U(p_0,rho=rho_0)=0 for EQS-Type=99
SOURCE
real, save :: p_0 =0.8
NOTES
- Units : GeV
baryonPotentialModule/U_0 [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
U(p=0,rho=rho_0) for EQS-Type=99
SOURCE
real, save :: U_0 =0.075
NOTES
- Units : GeV
baryonPotentialModule/bindingEnergy [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Nuclear matter binding energy for EQS-Type=99
SOURCE
real, save :: bindingEnergy=0.016
NOTES
- Units : GeV
baryonPotentialModule/compressibility [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
Nuclear matter compressibility for EQS-Type=99
SOURCE
real, save :: compressibility=0.290
NOTES
- Units : GeV
baryonPotentialModule/nLoopReAdjust [ Global module-variables ]
[ Top ] [ baryonPotentialModule ] [ Global module-variables ]
PURPOSE
number of iterations, if density is readjusted (cf. type(nucleus)%ReAdjustForConstBinding)
SOURCE
integer, save :: nLoopReAdjust = 10
NOTES
It is necessary to reiterate (at least for momentum dependent potentials), since we calculate the potential for a given pF and then calculate for the radjusting a new pF
baryonPotentialModule/baryonPotential [ Namelists ]
[ Top ] [ baryonPotentialModule ] [ Namelists ]
NAME
NAMELIST /baryonPotential/
PURPOSE
Includes the following switches:
- symmetriePotFlag
- EQS_Type
- DeltaPot
- HypPot
- noPerturbativePotential
- rho_0
- p_0
- U_0
- bindingEnergy
- compressibility
- SurfacePotFlag
- nLoopReAdjust
- dsymm
baryonPotentialModule/BaryonPotential [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
function BaryonPotential (teilchen, med, positionNotSet, EQS_in)
INPUTS
- type(particle) :: teilchen -- boosted to LRF
- type(medium) :: med -- medium information
- logical :: positionNotSet -- .true. : %position of particle is not well defined
- integer, OPTIONAL :: EQS_in -- If present, then we use EQS_in as EQS type, if not present then EQS is chosen according to EQS_type
Some routines like Yukawa might need the position of the particle, and not only the densities. Therefore, the positionNotSet-flag is used to check whether the position is actually set. If e.g. Yukawa is used and the position is not set, then the code stops.
NOTES
Baryon potential is defined as 0th component of a vector potential in the LRF.
baryonPotentialModule/variableSkyrme [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
real function variableSkyrme(rho,p)
PURPOSE
- This function evaluates a variable Skyrme mean field potential
- It's variable in the sense, that the potential parameters are no longer fixed by the above tables "alpha", "beta",... They are evaluated based on the input values of rhoNull, p_0, u_0, bindingEnergy and compressibility!
INPUTS
- real, intent(in) :: rho ! Density in GeV**3
- real, intent(in) :: p ! Momentumin GeV
OUTPUT
- Single particle potential in GeV
NOTES
- This function is initializing the potential parameters when its called for the first time.
- See Oliver's Phd thesis appendix A.4
- Stops the code if there is no solution for the parameters.
baryonPotentialModule/momentumDependentPart [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
real function momentumDependentPart(pin,c,lambda,rho,pF_in)
PURPOSE
This function provides the analytical momdep. potential.
INPUTS
- real :: pIn -- absolute momentum in GeV in LRF
- real :: c -- parameter of potential
- real :: lambda -- parameter of potential
- real :: rho -- baryon density in LRF
- real, OPTIONAL :: pF_in -- value of fermi mom to use (in GeV)
NOTES
see effenberger, dr.-thesis, pages 14-16
baryonPotentialModule/rhoLaplace [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
function rhoLaplace(rvec,a)
PURPOSE
Calculates div(grad(rho))
baryonPotentialModule/LDApotential(teilchen) [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
function LDApotential(teilchen)
PURPOSE
Calculates the Potential for a nucleus initialised with the LDA approach
baryonPotentialModule/LDApotentialWelke(teilchen) [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
function LDApotentialWelke(teilchen)
PURPOSE
Calculates the Potential for a nucleus initialised with the LDA approach and additionally a momentum dependent part
baryonPotentialModule/SurfacePart [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
function SurfacePart(teilchen,spar)
PURPOSE
Determines the surface contribution to the total baryon potential
INPUTS
- type(particle) :: teilchen -- particle, in a position of which grad(\nabla\rho) has to be calculated.
- real :: spar -- Parameter of surface part of baryon potential
OUTPUT
real :: SurfacePart
baryonPotentialModule/rearrangementPotential [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
real function rearrangementPotential(teilchen,rhoBaryon,rhoProton,rhoNeutron)
PURPOSE
Returns the value of the rearrangement potential.
INPUTS
- type(particle) :: teilchen -- particle whose rearrangement potential should be calculated. It should be boosted to LRF. The position of the particle must be set!!!
- real :: rhoBaryon,rhoProton,rhoNeutron -- densities in LRF
OUTPUT
function value
NOTES
Notation according to Teis Dr.-thesis. Pages 76-78.
This is *not* the rearrangement potential, but the expression
-(1/2 U_b + U_r)
which is needed to calculate the beinding energy correctly.
baryonPotentialModule/HandPotentialToDensityStatic [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
subroutine HandPotentialToDensityStatic(nuc)
PURPOSE
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- The potential at input is frozen and stored in a r dependent grid
- EQS_Type is set to 98
- The static density is adjusted.
baryonPotentialModule/getNoPertPot_baryon [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
logical function getNoPertPot_baryon
PURPOSE
Returns the flag noPerturbativePotential
baryonPotentialModule/getsymmetriePotFlag_baryon [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
logical function getsymmetriePotFlag_baryon
PURPOSE
Returns the flag symmetriePotFlag
baryonPotentialModule/getPotentialEQSType [ Functions ]
[ Top ] [ baryonPotentialModule ] [ Functions ]
NAME
integer function getPotentialEQSType
PURPOSE
Returns the EQS type of the potential (0 = no potential)