Changes between Version 3 and Version 4 of jobCards/deuterium

Timestamp:

Dec 2, 2021, 3:48:35 PM (2 years ago)

Author:

gallmei

Comment:

Update infos on this page

jobCards/deuterium

@@ -5 +5 @@
 First the target has to be adjusted to Deuterium:
 {{{
-$target
-target_Z=1
-target_A=2
-fermimotion=.true.
-$end
+! file: code/density/nucleus.f90
+&target
+      Target_Z =  1, Target_A =  2
+/
 }}}
 To distribute the nucleons in position and momentum space we can choose between two different wave function models:
 {{{
-$deuteriumFermi
-waveFunction_switch=2
-! 1=Bonn
-! 2=Argonne
-$end
+! file: code/init/deuterium.f90
+&deuteriumFermi
+      waveFunction_switch=2     ! 1=Bonn, 2=Argonne
+/
 }}}
 Next, we need to define a potential to bind the two nucleons. For this we can't use a mean field, because Deuterium represents a too small system. Instead we use a real two-body potential. Using the parallel ensemble technique, the potential ''V'' for each nucleon in the ''j''th ensemble is given by \( V=V_\text{2-body}(r_{1,j}-r_{2,j}) \)
 where \( r_{i,j} \) is the position of the ''i''th nucleon in the ''j''th ensemble. For the full ensemble method, a Deuterium potential is not yet properly implemented. So we choose for the general input and the propagation routines the following switches:
 {{{
+! file: ./inputOutput/input.f90
 $input
-delta_T     = 0.025              ! small time step sizes since the two-body potential is stiff and therefore the propagation is sensitive to too large time steps
-fullensemble=.false.             ! => use parallel ensemble technique
-freezeRealParticles=.false.     
-length_perturbative=1            ! We don't use perturbative particles, see comments below
+      delta_T     = 0.025     ! small time step size since the two-body potential 
+                              ! is stiff and therefore the propagation is sensitive 
+                              ! to too large time steps
+      fullensemble=.false.    ! => use parallel ensemble technique
+      freezeRealParticles=.false.     
 ...
-$end
+/
 
-$initDensity
-densitySwitch=1
-splineExtraPolation=.true. !Switch for linear spline extrapolation for dynamically calculated density: Extrapolates density between 
-gridPoints(1)=100
-gridPoints(2)=100
-gridPoints(3)=100
-gridSize(1)=8.
-gridSize(2)=8.
-gridSize(3)=8.
-$end
+! file: code/density/density.f90
+&initDensity
+      densitySwitch = 2         ! 2=analytic
+/
 
-$propagation
-delta_P=0.01                ! Delta Momentum for derivatives
-DerivativeType=2            ! 1=first order Range-Kutta, 2=second order Range-Kutta
-predictorCorrector=.true.   ! Whether to use a predictor/corrector algorithm to do the propagation
-$end
+! file: code/propagation/propagation.f90
+&propagation
+      RungeKuttaOrder=2         ! 2=second order Runge-Kutta
+/
 
-$baryonPotential
-EQS_Type=7   ! => Two body potential for deuterium
-DeltaPot=1   ! Switch for potential of spin=3/2 resonances
-             ! 1=nucleon (spin=1/2) potential times  3/5   [according to ericson/Weise book]
-             ! 2= 100 MeV *rho/rhoNull
-symmetriePotFlag=.false.   ! Switch for the assymetry term in the nucleon potential
-$end
-
-$Yukawa
-yukawaFlag=.false.  !decides whether Yukawa is switched off(.false.)  or on (.true.)
-$end
+! file: code/potential/baryonPotential.f90
+&baryonPotential
+      EQS_Type=7                ! => Two body potential for deuterium
+/
 }}}
 
 
-Oliver prefers not to use perturbative particles with Deuterium, since there is no unperturbed nucleus left if there is a nuclear reaction in deuterium. So he chooses 
+Up to now is probably all you have to change.
+
+Oliver prefers not to use perturbative particles with Deuterium, since there is no unperturbed nucleus left if there is a nuclear reaction in deuterium. So he chooses for his special case
 {{{
+! file: ./inputOutput/input.f90
+$input
+...
+      length_perturbative=1   ! We don't use perturbative particles
+...
+/
+
 $low_photo_induced
 ...
-realRun=.true. ! => reaction products are set into real particle vector
-$end
+      realRun=.true. ! => reaction products are set into real particle vector
+/
 }}}
+This may not work if you use another initialization.
 
-