gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/baryonPotentialModule [ Modules ]

[ Top ] [ 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:


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

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

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)