gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/RMF [ Modules ]

[ Top ] [ Modules ]

NAME

module RMF

PURPOSE

Includes all information about relativistic mean-field potential for baryons and mesons.

NOTES

  • When hyperon coupling is scaled by the well known factor of 2/3, the kaons are scaled by the factor 1/3 in order to compensate the missing self energy between incoming and outgoing channel. This is because the threshold condition, e.g. \pi N->YK sqrt(s*)>m*_y+m*_k, in the medium assumes no changes in the self energy between initial and final states (see Alexei's paper on three-body collisions). This prescription should better not be used at energies near the kaon-production threshold, since in this method the kaon potential is not consistent within Chiral Perturbation Theory or One-Boson-Exchange models.
  • The same non-trivial feature appears if the baryon self energies depend on isospin. Presently no isospin-dependent part is included in the baryon fields.
  • Going beyond this simple approximation means to explicitly include different threshold conditions for all channels considered in the collision term.


RMF/RMF_flag [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  logical, save :: RMF_flag = .false.

PURPOSE

If .true. then use relativistic mean fields.


RMF/N_set [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  integer, save :: N_set = 1

PURPOSE

Select which parameter set to use:

  • 1 --- NL1 from G.A. Lalazissis et al., PRC 55, 540 (1997), (K=211.29 MeV, m*/m=0.57)
  • 2 --- NL3 from G.A. Lalazissis et al., PRC 55, 540 (1997), (K=271.76 MeV, m*/m=0.60)
  • 3 --- NL2 set from A. Lang et al., NPA 541, 507 (1992), (K=210 MeV, m*/m=0.83)
  • 4 --- NLZ2 set from M. Bender et al., PRC 60, 34304 (1999), (K=172 MeV, m*/m=0.583)
  • 5 --- NL3* set from G.A. Lalazissis, private communication, (K=258.28 MeV, m*/m=0.594)
  • 6 --- Same as N_set=3, but including the rho meson.
  • 7 --- NL1 set from S.J. Lee et al., PRL 57, 2916 (1986), (K=212 MeV, m*/m=0.57)
  • 8 --- NL2 set from S.J. Lee et al., PRL 57, 2916 (1986), (K=399 MeV, m*/m=0.67)
  • 9 --- Set I from B. Liu et al., PRC 65, 045201 (2002), (K=240 MeV, m*/m=0.75)

*31 --- Parity doublet model Set P3 from D. Zschiesche et al., PRC 75, 055202 (2007)


RMF/grad_flag [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  logical, save, public :: grad_flag = .false.

PURPOSE

If .true. then include space derivatives of the fields.


RMF/lorentz_flag [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  logical, save, public :: lorentz_flag = .true.

PURPOSE

If .false. then the space components of the omega field are put to zero.


RMF/Tens_flag [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  logical, save, public :: Tens_flag = .false.

PURPOSE

If .true. then compute the energy-momentum tensor and four-momentum density field (not used in propagation).


RMF/flagCorThr [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  logical, save, public :: flagCorThr=.false.

PURPOSE

If .true. then the srtfree of colliding particles is corrected to ensure in-medium thresholds of BB -> BB and MB -> B


RMF/fact_pbar [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  real, save :: fact_pbar    = 1.

PURPOSE

Modification factor for the antiproton coupling constants.


RMF/fact_hyp [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  real, save :: fact_hyp     = 1.

PURPOSE

Modification factor for the hyperon coupling constants.


RMF/fact_antihyp [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  real, save :: fact_antihyp = 1.

PURPOSE

Modification factor for the antihyperon coupling constants.


RMF/fact_Xi [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  real, save :: fact_Xi      = 1.

PURPOSE

Modification factor for the Xi and XiStar coupling constants.


RMF/fact_antiXi [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  real, save :: fact_antiXi  = 1.

PURPOSE

Modification factor for the antiXi and antiXiStar coupling constants.


RMF/fact_kaon [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  real, save :: fact_kaon    = 0.

PURPOSE

Modification factor for the Kaon and antikaon coupling constants.


RMF/kaonpot_flag [ Global module-variables ]

[ Top ] [ RMF ] [ Global module-variables ]

SOURCE

  logical, save, public :: kaonpot_flag = .false.

PURPOSE

This switch turns on the Kaon potential in RMF mode.


RMF/ModificationFactor [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function ModificationFactor(Id,antiFlag)

PURPOSE

Returns the modification factor of the RMF coupling constants for a given particle.

INPUTS

  • integer, intent(in) :: Id ! Id of particle
  • logical, intent(in) :: antiFlag ! if .true. the particle is an antiparticle


RMF/init [ Subroutines ]

[ Top ] [ RMF ] [ Subroutines ]

NAME

subroutine init

PURPOSE

Reads input switches. Initializes the mean field parameters.


RMF/RMF_input [ Namelists ]

[ Top ] [ RMF ] [ Namelists ]

NAME

NAMELIST /RMF_input/

PURPOSE

Includes the following input switches:


RMF/walecka [ Subroutines ]

[ Top ] [ RMF ] [ Subroutines ]

NAME

subroutine walecka(rhobar,shift,em0,rhoscalar,endens,S,V,potential)

PURPOSE

Determine the mass shift of the nucleon in equilibrated isospin symmetric nuclear matter at zero temperature within Walecka model with nonlinear sigma-coupling.

INPUTS

  • real, intent(in) :: rhobar ! -- baryon density (fm^-3),
  • real, optional, intent(in) :: em0 ! -- starting value of mass for iterations (GeV),

OUTPUT

  • real, intent(out) :: shift ! = m - m^* -- mass shift (GeV),
  • real, optional, intent(out) :: rhoscalar ! -- scalar density (fm^-3),
  • real, optional, intent(out) :: endens ! -- energy density (GeV/fm^3),
  • real, optional, intent(out) :: pressure ! -- pressure (GeV/fm^3),
  • real, optional, intent(out) :: S ! -- scalar potential (GeV),
  • real, optional, intent(out) :: V ! -- vector potential (GeV),
  • real, optional, intent(out) :: potential ! -- Schroedinger equivalent potential (GeV)


RMF/fshift [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function fshift(rho)

PURPOSE

Fit of the nucleon mass shift m - m* for the various RMF parameter sets.

INPUTS

  • real, intent(in) :: rho -- baryon density (fm**-3)

OUTPUT

  • real :: fshift -- m - m* (GeV)

NOTES

This is a very rough fit which is only good to provide the starting value for iterations in walecka. The density rho must be in the interval from 0 up to 12*rhoNull.


RMF/PD [ Subroutines ]

[ Top ] [ RMF ] [ Subroutines ]

NAME

subroutine PD(mubStar,sigma,shift,flagPlot,sigma_inp,mub,rhoPlus,rhoMinus,rhoscalar,endens,pressure,S,V,potential)

PURPOSE

Determine the scalar field and mass shifts of the nucleon and its negative parity partner in equilibrated isospin symmetric nuclear matter at zero temperature within parity doublet model.

INPUTS

  • real, intent(in) :: mubStar ! = sqrt(pf_pm**2+m_pm**2) -- kinetic part of baryon chemical potential (GeV)
  • logical, optional, intent(in) :: flagPlot ! when = true -- equation for sigma field f(sigma)=0 is not solved, only function f vs sigma is plotted
  • real, optional, intent(in) :: sigma_inp ! -- starting value of sigma field (GeV),

OUTPUT

  • real, intent(out) :: sigma ! -- scalar field (GeV)
  • real, intent(out) :: shift(1:2) ! = m - m^* -- mass shift (GeV), 1 - nucleon, 2 - negative parity partner
  • real, optional, intent(out) :: mub ! -- baryon chemical potential (GeV)
  • real, optional, intent(out) :: rhoPlus ! -- nucleon density (fm^-3),
  • real, optional, intent(out) :: rhoMinus ! -- partner density (fm^-3),
  • real, optional, intent(out) :: rhoscalar ! -- scalar density (fm^-3),
  • real, optional, intent(out) :: endens ! -- energy density (GeV/fm^3),
  • real, optional, intent(out) :: pressure ! -- pressure (GeV/fm^3),
  • real, optional, intent(out) :: S(1:2) ! -- scalar potential (GeV),
  • real, optional, intent(out) :: V(1:2) ! -- vector potential (GeV),
  • real, optional, intent(out) :: potential(1:2) ! -- Schroedinger equivalent potential (GeV)


RMF/mPD [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function mPD

PURPOSE

Computes effective mass (in GeV) in the parity doublet model

INPUTS

  • real, intent(in) :: sigma -- scalar field (GeV)
  • integer, intent(in) :: parity -- +1 for nucleon, -1 for S11_1535


RMF/dmPD [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function dmPD

PURPOSE

Computes derivative of effective mass over sigma field d m/d sigma in the parity doublet model

INPUTS

  • real, intent(in) :: sigma -- scalar field (GeV)
  • integer, intent(in) :: parity -- +1 for nucleon, -1 for S11_1535


RMF/d2mPD [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function d2mPD

PURPOSE

Computes second derivative of effective mass over sigma field d^2 m/d sigma^2 (in GeV^-1) in the parity doublet model

INPUTS

  • real, intent(in) :: sigma -- scalar field (GeV)


RMF/f [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function f

PURPOSE

Computes analytically the expression

  3*a*\int_0^1 dx x^2/\sqrt(x^2+a^2)

INPUTS

  • real, intent(in) :: a -- dimensionless parameter equal to m^*/p_F


RMF/fprime [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function fprime(a)

PURPOSE

Computes analytically the derivative of function f(a) with respect to a.

INPUTS

  • real, intent(in) :: a -- dimensionless parameter equal to m^*/p_F


RMF/g [ Functions ]

[ Top ] [ RMF ] [ Functions ]

NAME

real function g(a)

PURPOSE

Computes analytically the expression

  \int_0^1 dx x^2*\sqrt(x^2+a^2)

INPUTS

  • real, intent(in) :: a -- dimensionless parameter equal to m^*/p_F