gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS

/densityStatic [ Modules ]

[ Top ] [ Modules ]

NAME

module densityStatic

PURPOSE

Collect routines for STATIC density calculations.

densityStatic/useCentroids [ Global module-variables ]

[ Top ] [ densityStatic ] [ Global module-variables ]

SOURCE

  logical, save :: useCentroids = .false.

PURPOSE

If this switch is 'true', then the density of the proton and neutron centers will be tabulated that is different from the matter density.

NOTES

presently relevant only for densitySwitch_static=2 (Luis routine)

densityStatic/staticDensityInit [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine staticDensityInit (nuc)

PURPOSE

decide, which density parametrisation is used. Then tabulate this and also set the extreme values for the MC decision.

INPUTS

OUTPUT

densityStatic/TabulateZero [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateZero(nuc)

PURPOSE

Set the density table to zero

INPUTS

OUTPUT

densityStatic/TabulateSphere [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateFermiGas(nuc)

PURPOSE

Tabulate a sphere with constant density.

Tabulates the static density to make it available faster for later use

INPUTS

OUTPUT

densityStatic/TabulateDensityWoodsSaxon [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateDensityWoodsSaxon (nuc)

PURPOSE

Tabulate the Woods-Saxon distribution. Tabulates the static density to make it available faster for later use.

INPUTS

OUTPUT

densityStatic/TabulateDensityHarmOsc [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateDensityHarmOsc(nuc)

PURPOSE

Tabulate the density distribution according harmonic oscillator shell modell.

Tabulates the static density to make it available faster for later use

INPUTS

OUTPUT

NOTES

The parameter are taken from the FRITIOF package.

densityStatic/TabulateDensityLenske [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateDensityLenske (nuc)

PURPOSE

Tabulate the density distribution according to Woods-Saxon distribution but with refined charge radii for proton and neutron according to H. Lenske. Tabulates the static density to make it available faster for later use.

INPUTS

OUTPUT

NOTES

Everything in fm.

densityStatic/TabulateDensityLuis [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateDensityLuis(nuc)

PURPOSE

Tabulate the density distribution of matter (p and n) and the density of centers ( p and n number densities ) following J. Nieves et al., Pionic atoms... NPA554

Tabulates the static density to make it available faster for later use

INPUTS

OUTPUT

NOTES

everything in fm

densityStatic/TabulateDensityExRTF [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateDensityExRTF(nuc)

PURPOSE

Tabulate the density distribution according to Relativistic Thomas-Fermi model code from Horst Lenske.

Tabulates the static density to make it available faster for later use

INPUTS

OUTPUT

NOTES

everything in fm

densityStatic/SearchMaxVals [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine SearchMaxVals(nuc)

PURPOSE

go through the tabulated distributions to search for the extrema

densityStatic/staticDensity [ Functions ]

[ Top ] [ densityStatic ] [ Functions ]

NAME

type(dichte) function staticDensity(r,nucl)

PURPOSE

gives density in the restframe of the nucleus "nucl" at position "r"

INPUTS

  • real, dimension(1:3) :: r -- position where density should be calculated
  • type(tNucleus),pointer :: nucl -- nucleus which is regarded

USAGE

(dichte)=staticDensity(...)

densityStatic/densityLuis [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine densityLuis(r,z,a,rhop,rhon,rp,ap,rho0p,rn,an,rho0n)

PURPOSE

This routine calculates the proton and neutron densities following J. Nieves et al., Pionic atoms... NPA554, 509 (1993)

returns per default the density of matter, set useCentroids=.true. to switch to density of centers

INPUTS

  • real :: r -- radius (fm)
  • integer :: z -- charge of the nucleus
  • integer :: a -- atomic number

RESULT

  • real :: rhop,rhon -- Proton and neutron densities (fm^-3) at r
  • real :: rp,ap,rho0p,rn,an,rho0n -- parameters of the density distributions

densityStatic/TabulateDensityBirger [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateDensityBirger(nuc)

PURPOSE

Tabulate the density distribution based on a local density approximation first described by Brueckner et al.

INPUTS

OUTPUT

densityStatic/TabulateDensityBirgerWelke [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine TabulateDensityBirgerWelke(nuc)

PURPOSE

Tabulate the density distribution based on a local density approximation first described by Brueckner et al. and a momentum-dependent potential

INPUTS

OUTPUT

densityStatic/ReAdjust [ Subroutines ]

[ Top ] [ densityStatic ] [ Subroutines ]

NAME

subroutine ReAdjust(nuc, potP, potN, potC)

PURPOSE

This routine recalculates the density distributions for protons and neutrons by considering the given potentials as static and fulfill the condition 

   sqrt(p_F^2+m_N^2) + U - m_N == E_sep ~ -8MeV

With the Local-Thomas-Fermi, we connect the resulting fermi momentum to a density, 

   rho = p_F^3/(3pi^2)

Since the potentials are given as function of r, we calculate rho(r).

Thus, given proton and nucleon baryon potential (for fixed momentum) and the coulomb potential, the parametrization of the nuclear density is readjusted.

This routine is called by baryonPotentialMain/HandPotentialToDensityStatic

INPUTS

  • type(tNucleus),pointer :: nuc -- the nucleus to consider
  • real, dimension(0:) :: potP, potN -- the proton,neutron potentials with p=pF. The dimension has to be identical to nuc%densTab(0: ,1:2).
  • real, dimension(0:) :: potC -- The Coulomb potential (>0, in GeV)

OUTPUT

  • nuc%densTab(0: ,1:2) is changed