TABLE OF CONTENTS
- 1. /skyrme
- 1.1. skyrme/evaluate_skyrmeParameters
- 1.2. skyrme/U
- 1.3. skyrme/E_div_A
- 1.4. skyrme/dRho_E_div_A
- 1.5. skyrme/dRhoSquared_E_div_A
- 1.6. skyrme/pf
- 1.7. skyrme/ekin
- 1.8. skyrme/d_drho_ekin
- 1.9. skyrme/d_drhoSquared_ekin
- 1.10. skyrme/getParameters
- 1.11. skyrme/check
/skyrme [ Modules ]
NAME
module skyrme
PURPOSE
- This module includes various routines for the Skyrme mean field potential
NOTES
- For documentation see appendix A.4 of Oliver's thesis
skyrme/evaluate_skyrmeParameters [ Subroutines ]
[ Top ] [ skyrme ] [ Subroutines ]
NAME
subroutine evaluate_skyrmeParameters (rhoNull,pNull, uNull,bindingEnergy, compressibility,A,B,C,tau,lambda,success)
PURPOSE
Evaluates the Skyrme parameters A,B,C,tau,lambda for given input * rhoNull -- saturation density in GeV^3 * pNull -- U(rho=rhoNull,p=p_0)=0 -> p_0 is root of potential in GeV * uNull -- U(rho=rhoNull,p=0 )=uNull -> depth of potential at p=0 in GeV * bindingEnergy -- Nuclear matter binding energy in GeV * compressibility -- Nuclear matter compressibility in GeV
OUTPUT
- real, intent(out) :: A,B,C,lambda ! Potential parameters, units: GeV
- real, intent(out) :: tau ! Potential parameter , no units
- logical, intent(out) :: success ! true if parameters could be determined
NOTES
See appendix A.4 of Oliver's thesis for details
skyrme/U [ Functions ]
[ Top ] [ skyrme ] [ Functions ]
NAME
real function U(rho, p, rhoNull, A,B,C, tau, lambda)
PURPOSE
Returns the single-particle potential in units of GeV.
INPUTS
- rhoNull -- saturation density in GeV^3
- rho -- baryon density in GeV^3
- p -- momentum in GeV
- A,B,C,lambda -- potential parameters in GeV
- tau -- potential parameters (no units)
NOTES
See appendix A.4 of Oliver's thesis for details
skyrme/E_div_A [ Functions ]
[ Top ] [ skyrme ] [ Functions ]
NAME
real function E_div_A(rho, rhoNull, A,B,C, tau, lambda)
PURPOSE
Returns the nuclear binding energy in GeV
INPUTS
- rhoNull -- saturation density in GeV^3
- rho -- baryon density in GeV^3
- A,B,C,lambda -- potential parameters in GeV
- tau -- potential parameters (no units)
NOTES
See appendix A.4 of Oliver's thesis for details
skyrme/dRho_E_div_A [ Functions ]
[ Top ] [ skyrme ] [ Functions ]
NAME
real function dRho_E_div_A(rho, rhoNull, A,B,C, tau, lambda)
PURPOSE
Returns the rho-derivative of the nuclear binding energy
INPUTS
- rhoNull -- saturation density in GeV^3
- rho -- baryon density in GeV^3
- A,B,C,lambda -- potential parameters in GeV
- tau -- potential parameters (no units)
NOTES
See appendix A.4 of Oliver's thesis for details
skyrme/dRhoSquared_E_div_A [ Functions ]
[ Top ] [ skyrme ] [ Functions ]
NAME
real function dRhoSquared_E_div_A(rho, rhoNull, A,B,C, tau, lambda)
PURPOSE
Returns the second rho-derivative of the nuclear binding energy
INPUTS
- rhoNull -- saturation density in GeV^3
- rho -- baryon density in GeV^3
- A,B,C,lambda -- potential parameters in GeV
- tau -- potential parameters (no units)
NOTES
See appendix A.4 of Oliver's thesis for details
skyrme/pf [ Functions ]
[ Top ] [ skyrme ] [ Functions ]
NAME
real function pf(rho)
PURPOSE
Returns the Fermi momentum in GeV
INPUTS
- rho -- baryon density in GeV^3
skyrme/ekin [ Subroutines ]
[ Top ] [ skyrme ] [ Subroutines ]
NAME
real function ekin(rho)
PURPOSE
Evaluates the kinetic energy density
INPUTS
real, intent(in) :: rho ! density in GeV^3
NOTES
- See appendix A.4 of Oliver's thesis for details
skyrme/d_drho_ekin [ Subroutines ]
[ Top ] [ skyrme ] [ Subroutines ]
NAME
real function d_drho_ekin(rho)
PURPOSE
Evaluates the rho-derivative of the kinetic energy density
INPUTS
real, intent(in) :: rho ! density in GeV^3
NOTES
- See appendix A.4 of Oliver's thesis for details
skyrme/d_drhoSquared_ekin [ Subroutines ]
[ Top ] [ skyrme ] [ Subroutines ]
NAME
real function d_drhoSquared_ekin(rho)
PURPOSE
Evaluates the 2nd rho-derivative of the kinetic energy density
INPUTS
real, intent(in) :: rho ! density in GeV^3
NOTES
- See appendix A.4 of Oliver's thesis for details
skyrme/getParameters [ Subroutines ]
[ Top ] [ skyrme ] [ Subroutines ]
NAME
subroutine getParameters(rhoNull,pNull, uNull, compressibility,A,B,C,tau,lambda)
PURPOSE
Evaluates A,B,C and tau if lambda and rhoNull, pNull, uNull, compressibility are given
INPUTS
real, intent(in) :: lambda real, intent(in) :: rhoNull,pNull, uNull, compressibility
OUTPUT
real, intent(out) :: A,B,C,tau
NOTES
- All inputs and outputs in multiples of GeV
- See appendix A.4 of Oliver's thesis for details
skyrme/check [ Subroutines ]
[ Top ] [ skyrme ] [ Subroutines ]
NAME
subroutine check(rhoNull, pNull, bindingEnergy, compressibility, A,B,C, tau, lambda)
PURPOSE
Checks whether A,B,C, tau, lambda are valid parameters for given input rhoNull, pNull, bindingEnergy, compressibility
INPUTS
- real, intent(in) :: rhoNull, pNull, bindingEnergy, compressibility, A,B,C, tau, lambda
- all units in multiples of GeV!
NOTES
See appendix A.4 of Oliver's thesis for details