ID |
Date |
Author |
Subject |
|
12
|
Wed Nov 8 15:20:17 2023 |
Selin Ezgi Birincioglu | TALYS Parameters from Different Models |
This is the parameters of Ca42+n and Ca50+n acquired from differen models. They were chosen as they have Gaussian distributions which are consistent with TALYS calculations of MACS. Ca42+n has experimental values while Ca50+n does not. |
| Attachment 1: models.xlsx
|
|
11
|
Mon Nov 6 09:20:53 2023 |
Selin Ezgi Birincioglu | TALYS including MACS |
|
| Attachment 1: TALYSparameters.xlsx
|
|
10
|
Sat Nov 4 14:05:19 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of Ne22+n |
1. This graph is Ne22+n with a log scale. Theorotical cross section is around 0.05 while the experimental is 0.061, which matches relatively well.
2. This the the distribution of l = 0 with log scale.
3. This the the distribution of l = 1 with log scale.
4. This is the TALYS output for Ne22+n.
5. mygoe.f input code which consists of TALYS parameters (Gg, D, S). |
| Attachment 1: macs1.pdf
|
|
| Attachment 2: macs01.pdf
|
|
| Attachment 3: macs11.pdf
|
|
| Attachment 4: talys.out
|
TALYS-1.96 (Version: December 30, 2021)
Copyright (C) 2021 A.J. Koning, S. Hilaire and S. Goriely
Dimensions - Cross sections: mb, Energies: MeV, Angles: degrees
########## USER INPUT ##########
USER INPUT FILE
#
# general
#
projectile n
element ne
mass 22
energy 1.
astroe 0.02
racap y
#
# output
#
#channels y
outgamma y
outdensity y
outlevels y
outbasic y
USER INPUT FILE + DEFAULTS
Keyword Value Variable Explanation
#
# Four main keywords
#
projectile n ptype0 type of incident particle
element Ne Starget symbol of target nucleus
mass 22 mass mass number of target nucleus
energy energyfile file with incident energies
#
# Basic physical and numerical parameters
#
ejectiles g n p d t h a outtype outgoing particles
maxz 7 maxZ maximal number of protons away from the initial compound nucleus
maxn 10 maxN maximal number of neutrons away from the initial compound nucleus
bins 40 nbins number of continuum excitation energy bins
equidistant y flagequi flag to use equidistant excitation bins instead of logarithmic bins
equispec n flagequispec flag to use equidistant bins for emission spectra
popmev n flagpopmev flag to use initial population per MeV instead of histograms
segment 1 segment number of segments to divide emission energy grid
maxlevelstar 30 nlevmax maximum number of included discrete levels for target
maxlevelsres 10 nlevmaxres maximum number of included discrete levels for residual nucleus
maxlevelsbin g 10 nlevbin maximum number of included discrete levels for gamma channel
maxlevelsbin n 30 nlevbin maximum number of included discrete levels for neutron channel
maxlevelsbin p 10 nlevbin maximum number of included discrete levels for proton channel
maxlevelsbin d 5 nlevbin maximum number of included discrete levels for deuteron channel
maxlevelsbin t 5 nlevbin maximum number of included discrete levels for triton channel
maxlevelsbin h 5 nlevbin maximum number of included discrete levels for helium-3 channel
maxlevelsbin a 10 nlevbin maximum number of included discrete levels for alpha channel
ltarget 0 ltarget excited level of target
isomer 1.00E+00 isomer definition of isomer in seconds
transpower 15 transpower power for transmission coefficient limit
transeps 1.00E-18 transeps limit for transmission coefficient
xseps 1.00E-17 xseps limit for cross sections
popeps 1.00E-13 popeps limit for population cross section per nucleus
Rfiseps 1.00E-06 Rfiseps ratio for limit for fission cross section per nucleus
elow 1.00E-06 elow minimal incident energy for nuclear model calculations
angles 90 nangle number of angles
anglescont 18 nanglecont number of angles for continuum
anglesrec 1 nanglerec number of recoil angles
maxenrec 20 maxenrec number of recoil energies
channels n flagchannels flag for exclusive channels calculation
maxchannel 4 maxchannel maximal number of outgoing particles in individual channel description
micro n flagmicro flag for completely microscopic Talys calculation
best n flagbest flag to use best set of adjusted parameters
bestbranch y flagbestbr flag to use flag to use only best set of branching ratios
bestend n flagbestend flag to put best set of parameters at end of input file
relativistic y flagrel flag for relativistic kinematics
recoil n flagrecoil flag for calculation of recoils
labddx n flaglabddx flag for calculation of DDX in LAB system
recoilaverage n flagrecoilav flag for average velocity in recoil calculation
channelenergy n flagEchannel flag for channel energy for emission spectrum
reaction y flagreaction flag for calculation of nuclear reactions
ngfit n flagngfit flag for using fitted (n,g) nuclear model parameters
nnfit n flagnnfit flag for using fitted (n,n), (n,2n) and (n,p) nuclear model parameters
nafit n flagnafit flag for using fitted (n,a) nuclear model parameters
astro y flagastro flag for calculation of astrophysics reaction rate
astrogs y flagastrogs flag for calculation of astrophysics reaction rate with target in ground state only
astroex n flagastroex flag for calculation of astrophysics reaction rate to long-lived excited states
nonthermlev -1 nonthermlev excited level non-thermalized in the calculation of astrophysics rate
massmodel 2 massmodel model for theoretical nuclear mass
expmass y flagexpmass flag for using experimental nuclear mass if available
disctable 1 disctable table with discrete levels
production n flagprod flag for isotope production
outfy n flagoutfy flag for output detailed fission yield calculation
gefran 50000 gefran number of random events for GEF calculation
Estop 1000.000 Estop incident energy above which TALYS stops
rpevap n flagrpevap flag for evaporation of residual products at high incident energies
maxZrp 12 maxZrp maximal number of protons away from the initial compound nucleus before residual evaporation
maxNrp 32 maxNrp maximal number of neutrons away from the initial compound nucleus before residual evaporation
#
# Optical model
#
localomp y flaglocalomp flag for local (y) or global (n) optical model
dispersion n flagdisp flag for dispersive optical model
jlmomp n flagjlm flag for using semi-microscopic JLM OMP
riplomp n flagriplomp flag for RIPL OMP
riplrisk n flagriplrisk flag for going outside RIPL mass validity range
optmodall n flagompall flag for new optical model calculation for all residual nuclei
incadjust y flagincadj flag for OMP adjustment on incident channel also
omponly n flagomponly flag to execute ONLY an optical model calculation
autorot n flagautorot flag for automatic rotational coupled channels calculations for A > 150
spherical n flagspher flag to force spherical optical model
soukho y flagsoukho flag for Soukhovitskii OMP for actinides
coulomb y flagcoulomb flag for Coulomb excitation calculation with ECIS
statepot n flagstate flag for optical model potential for each excited state
maxband 0 maxband highest vibrational band added to rotational model
maxrot 2 maxrot number of included excited rotational levels
sysreaction sysreaction particles with reaction cross section from systematics
rotational n rotational particles with possible rotational optical model
core -1 core even-even core for weakcoupling (-1 or 1)
ecissave n flagecissave flag for saving ECIS input and output files
eciscalc y flageciscalc flag for new ECIS calculation for outgoing particles and energy grid
inccalc y flaginccalc flag for new ECIS calculation for incident channel
endfecis y flagendfecis flag for new ECIS calculation for ENDF-6 files
radialmodel 2 radialmodel model for radial matter densities (JLM OMP only)
jlmmode 0 jlmmode option for JLM imaginary potential normalization
alphaomp 6 alphaomp alpha OMP (1=normal, 2= McFadden-Satchler, 3-5= folding potential, 6,8= Avrigeanu, 7=Nolte)
deuteronomp 1 deuteronomp deuteron OMP (1=normal, 2=Daehnick, 3=Bojowald, 4=Han-Shi-Shen, 5=An-Cai)
#
# Compound nucleus
#
widthfluc y flagwidth flag for width fluctuation calculation
widthmode 1 wmode designator for width fluctuation model
compound y flagcomp flag for compound nucleus model
fullhf n flagfullhf flag for full spin dependence of transmission coefficients
eciscompound n flageciscomp flag for compound nucleus calculation by ECIS
cpang n flagcpang flag for compound angular distribution calculation for incident charged particles
urr n flagurr flag for URR calculation
urrnjoy n flagurrnjoy flag for normalization of URR parameters with NJOY method
lurr 2 lurr maximal orbital angular momentum for URR
#
# Gamma emission
#
gammax 2 gammax number of l-values for gamma multipolarity
strength 9 strength model for E1 gamma-ray strength function
strengthM1 3 strengthM1 model for M1 gamma-ray strength function
electronconv y flagelectron flag for application of electron conversion coefficient
racap y flagracap flag for radiative capture model
ldmodelracap 1 ldmodelracap level density model for direct radiative capture
upbend y flagupbend flag for low-energy upbend of photon strength function
psfglobal n flagpsfglobal flag for global photon strength functions only
#
# Pre-equilibrium
#
preequilibrium n flagpreeq flag for pre-equilibrium calculation
preeqmode 2 preeqmode designator for pre-equilibrium model
multipreeq n flagmulpre flag for multiple pre-equilibrium calculation
mpreeqmode 2 mpreeqmode designator for multiple pre-equilibrium model
breakupmodel 1 breakupmodel model for break-up reaction: 1. Kalbach 2. Avrigeanu
phmodel 1 phmodel particle-hole state density model
pairmodel 2 pairmodel designator for pre-equilibrium pairing model
preeqspin 1 pespinmodel model for pre-equilibrium spin distribution
giantresonance y flaggiant flag for collective contribution from giant resonances
preeqsurface y flagsurface flag for surface effects in exciton model
preeqcomplex y flagpecomp flag for Kalbach complex particle emission model
twocomponent y flag2comp flag for two-component pre-equilibrium model
ecisdwba y flagecisdwba flag for new ECIS calculation for DWBA for MSD
onestep n flagonestep flag for continuum one-step direct only
#
# Level densities
#
ldmodel 1 ldmodelall level density model
ldmodelCN 1 ldmodelCN level density model for compound nucleus
shellmodel 1 shellmodel model for shell correction energies
kvibmodel 2 kvibmodel model for vibrational enhancement
spincutmodel 1 spincutmodel model for spin cutoff factor for ground state
asys n flagasys flag for all level density parameters a from systematics
parity n flagparity flag for non-equal parity distribution
colenhance n flagcolall flag for collective enhancement of level density for all nuclides
ctmglobal n flagctmglob flag for global CTM model (no discrete level info)
gshell n flaggshell flag for energy dependence of single particle level density parameter g
colldamp n flagcolldamp flag for damping of collective effects in effective level density
#
# Fission
#
fission n flagfission flag for fission
fismodel 3 fismodel fission model
fismodelalt 3 fismodelalt alternative fission model for default barriers
hbstate y flaghbstate flag for head band states in fission
class2 y flagclass2 flag for class2 states in fission
fispartdamp n flagfispartdamp flag for fission partial damping
massdis n flagmassdis flag for calculation of fission fragment mass yields
ffevaporation y flagffevap flag for calculation of particle evaporation from fission fragment mass yields
fisfeed n flagfisfeed flag for output of fission per excitation bin
fymodel 2 fymodel fission yield model, 1: Brosa 2: GEF
ffmodel 1 ffmodel fission fragment model, 1: GEF 2: HF3D (Okumura) 3: SPY
pfnsmodel 1 pfnsmodel PFNS model, 1: Iwamoto 2: from FF decay
ffspin n flagffspin flag to use spin distribution in initial FF population
#
# Output
#
outmain y flagmain flag for main output
outbasic y flagbasic flag for output of basic information and results
outpopulation y flagpop flag for output of population
outcheck y flagcheck flag for output of numerical checks
outlevels y flaglevels flag for output of discrete level information
outdensity y flagdensity flag for output of level densities
outomp y flagoutomp flag for output of optical model parameters
outkd n flagoutkd flag for output of KD03 OMP parameters
outdirect y flagdirect flag for output of direct reaction results
outinverse y flaginverse flag for output of transmission coefficients and inverse reaction cross sections
outdecay n flagdecay flag for output of decay of each population bin
outtransenergy y flagtransen flag for output of transmission coefficients per energy
outecis n flagoutecis flag for output of ECIS results
outgamma y flaggamma flag for output of gamma-ray information
outpreequilibrium n flagpeout flag for output of pre-equilibrium results
outfission n flagfisout flag for output of fission information
outdiscrete y flagdisc flag for output of discrete state cross sections
outspectra n flagspec flag for output of double-differential cross sections
outbinspectra n flagbinspec flag for output of emission spectrum per excitation bin
resonance n flagres flag for output of low energy resonance cross sections
group n flaggroup flag for output of low energy groupwise cross sections
addiscrete y flagadd flag for addition of discrete states to spectra
addelastic y flagaddel flag for addition of elastic peak to spectra
outangle n flagang flag for output of angular distributions
outlegendre n flaglegendre flag for output of Legendre coefficients
ddxmode 0 ddxmode mode for double-differential cross sections
outdwba n flagoutdwba flag for output of DWBA cross sections for MSD
outgamdis n flaggamdis flag for output of discrete gamma-ray intensities
outexcitation n flagexc flag for output of excitation functions
components n flagcompo flag for output of cross section components
endf n flagendf flag for information for ENDF-6 file
endfdetail y flagendfdet flag for detailed ENDF-6 information per channel
sacs n flagsacs flag for statistical analysis of cross sections
partable n flagpartable flag for output of model parameters on separate file
block n flagblock flag to block spectra, angle and gamma files
########## BASIC REACTION PARAMETERS ##########
Projectile : neutron Mass in a.m.u. : 1.008665
Target : 22Ne Mass in a.m.u. : 21.991385
Included channels:
gamma
neutron
proton
deuteron
triton
helium-3
alpha
73 incident energies (LAB):
2.092E-06
4.183E-06
6.275E-06
8.367E-06
1.046E-05
2.092E-05
3.138E-05
4.183E-05
5.229E-05
6.275E-05
7.321E-05
8.367E-05
9.413E-05
1.046E-04
2.092E-04
3.138E-04
4.183E-04
5.229E-04
6.275E-04
7.321E-04
8.367E-04
9.413E-04
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.010
0.015
0.019
0.023
0.027
0.031
0.036
0.040
0.044
0.048
0.052
0.056
0.061
0.065
... 38605 more lines ...
|
| Attachment 5: Ne23run
|
#!/bin/bash
rm -r Ne23/
mkdir Ne23
for((c=1;c<=1000;c++))
do
./a.out << EOF
22 0 # amass spin
0 1000000 #min max
0 657808.56250 0.8570e-4 #ell D S
0.12333 0.024666 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Ne23/run1_0l_$c
cp out_all2 Ne23/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
22 0 # amass spin
0 1000000 #min max
1 290617.09375 2.3193e-4 #ell D S
0.18086 0.036172 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Ne23/run1_1l_$c
cp out_all2 Ne23/run2_1l_$c
cp out_all3 Ne23/run3_1l_$c
cp out_all4 Ne23/run4_1l_$c
rm out*
done
|
|
9
|
Sat Nov 4 14:01:18 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of Ca43+n |
1. This graph is the Gaussian fitted Ca43+n with rebinned to 10 total. Theorotical cross section is around 110 while the experimental is 66.1.
2. This the the distribution of l = 0 with bins = 10.
3. This the the distribution of l = 1 with bins = 10.
4. This is the TALYS output for Ca43+n.
5. mygoe.f input code which consists of TALYS parameters (Gg, D, S). |
| Attachment 1: macs1wrg.pdf
|
|
| Attachment 2: macs01wbin.pdf
|
|
| Attachment 3: macs11wbin.pdf
|
|
| Attachment 4: talys.out
|
TALYS-1.96 (Version: December 30, 2021)
Copyright (C) 2021 A.J. Koning, S. Hilaire and S. Goriely
Dimensions - Cross sections: mb, Energies: MeV, Angles: degrees
########## USER INPUT ##########
USER INPUT FILE
#
# general
#
projectile n
element ca
mass 43
energy 1.
astroe 0.02
racap y
#
# output
#
#channels y
outgamma y
outdensity y
outlevels y
outbasic y
USER INPUT FILE + DEFAULTS
Keyword Value Variable Explanation
#
# Four main keywords
#
projectile n ptype0 type of incident particle
element Ca Starget symbol of target nucleus
mass 43 mass mass number of target nucleus
energy energyfile file with incident energies
#
# Basic physical and numerical parameters
#
ejectiles g n p d t h a outtype outgoing particles
maxz 12 maxZ maximal number of protons away from the initial compound nucleus
maxn 21 maxN maximal number of neutrons away from the initial compound nucleus
bins 40 nbins number of continuum excitation energy bins
equidistant y flagequi flag to use equidistant excitation bins instead of logarithmic bins
equispec n flagequispec flag to use equidistant bins for emission spectra
popmev n flagpopmev flag to use initial population per MeV instead of histograms
segment 1 segment number of segments to divide emission energy grid
maxlevelstar 30 nlevmax maximum number of included discrete levels for target
maxlevelsres 10 nlevmaxres maximum number of included discrete levels for residual nucleus
maxlevelsbin g 10 nlevbin maximum number of included discrete levels for gamma channel
maxlevelsbin n 30 nlevbin maximum number of included discrete levels for neutron channel
maxlevelsbin p 10 nlevbin maximum number of included discrete levels for proton channel
maxlevelsbin d 5 nlevbin maximum number of included discrete levels for deuteron channel
maxlevelsbin t 5 nlevbin maximum number of included discrete levels for triton channel
maxlevelsbin h 5 nlevbin maximum number of included discrete levels for helium-3 channel
maxlevelsbin a 10 nlevbin maximum number of included discrete levels for alpha channel
ltarget 0 ltarget excited level of target
isomer 1.00E+00 isomer definition of isomer in seconds
transpower 15 transpower power for transmission coefficient limit
transeps 1.00E-18 transeps limit for transmission coefficient
xseps 1.00E-17 xseps limit for cross sections
popeps 1.00E-13 popeps limit for population cross section per nucleus
Rfiseps 1.00E-06 Rfiseps ratio for limit for fission cross section per nucleus
elow 1.00E-06 elow minimal incident energy for nuclear model calculations
angles 90 nangle number of angles
anglescont 18 nanglecont number of angles for continuum
anglesrec 1 nanglerec number of recoil angles
maxenrec 20 maxenrec number of recoil energies
channels n flagchannels flag for exclusive channels calculation
maxchannel 4 maxchannel maximal number of outgoing particles in individual channel description
micro n flagmicro flag for completely microscopic Talys calculation
best n flagbest flag to use best set of adjusted parameters
bestbranch y flagbestbr flag to use flag to use only best set of branching ratios
bestend n flagbestend flag to put best set of parameters at end of input file
relativistic y flagrel flag for relativistic kinematics
recoil n flagrecoil flag for calculation of recoils
labddx n flaglabddx flag for calculation of DDX in LAB system
recoilaverage n flagrecoilav flag for average velocity in recoil calculation
channelenergy n flagEchannel flag for channel energy for emission spectrum
reaction y flagreaction flag for calculation of nuclear reactions
ngfit n flagngfit flag for using fitted (n,g) nuclear model parameters
nnfit n flagnnfit flag for using fitted (n,n), (n,2n) and (n,p) nuclear model parameters
nafit n flagnafit flag for using fitted (n,a) nuclear model parameters
astro y flagastro flag for calculation of astrophysics reaction rate
astrogs y flagastrogs flag for calculation of astrophysics reaction rate with target in ground state only
astroex n flagastroex flag for calculation of astrophysics reaction rate to long-lived excited states
nonthermlev -1 nonthermlev excited level non-thermalized in the calculation of astrophysics rate
massmodel 2 massmodel model for theoretical nuclear mass
expmass y flagexpmass flag for using experimental nuclear mass if available
disctable 1 disctable table with discrete levels
production n flagprod flag for isotope production
outfy n flagoutfy flag for output detailed fission yield calculation
gefran 50000 gefran number of random events for GEF calculation
Estop 1000.000 Estop incident energy above which TALYS stops
rpevap n flagrpevap flag for evaporation of residual products at high incident energies
maxZrp 12 maxZrp maximal number of protons away from the initial compound nucleus before residual evaporation
maxNrp 32 maxNrp maximal number of neutrons away from the initial compound nucleus before residual evaporation
#
# Optical model
#
localomp y flaglocalomp flag for local (y) or global (n) optical model
dispersion n flagdisp flag for dispersive optical model
jlmomp n flagjlm flag for using semi-microscopic JLM OMP
riplomp n flagriplomp flag for RIPL OMP
riplrisk n flagriplrisk flag for going outside RIPL mass validity range
optmodall n flagompall flag for new optical model calculation for all residual nuclei
incadjust y flagincadj flag for OMP adjustment on incident channel also
omponly n flagomponly flag to execute ONLY an optical model calculation
autorot n flagautorot flag for automatic rotational coupled channels calculations for A > 150
spherical n flagspher flag to force spherical optical model
soukho y flagsoukho flag for Soukhovitskii OMP for actinides
coulomb y flagcoulomb flag for Coulomb excitation calculation with ECIS
statepot n flagstate flag for optical model potential for each excited state
maxband 0 maxband highest vibrational band added to rotational model
maxrot 2 maxrot number of included excited rotational levels
sysreaction sysreaction particles with reaction cross section from systematics
rotational n rotational particles with possible rotational optical model
core -1 core even-even core for weakcoupling (-1 or 1)
ecissave n flagecissave flag for saving ECIS input and output files
eciscalc y flageciscalc flag for new ECIS calculation for outgoing particles and energy grid
inccalc y flaginccalc flag for new ECIS calculation for incident channel
endfecis y flagendfecis flag for new ECIS calculation for ENDF-6 files
radialmodel 2 radialmodel model for radial matter densities (JLM OMP only)
jlmmode 0 jlmmode option for JLM imaginary potential normalization
alphaomp 6 alphaomp alpha OMP (1=normal, 2= McFadden-Satchler, 3-5= folding potential, 6,8= Avrigeanu, 7=Nolte)
deuteronomp 1 deuteronomp deuteron OMP (1=normal, 2=Daehnick, 3=Bojowald, 4=Han-Shi-Shen, 5=An-Cai)
#
# Compound nucleus
#
widthfluc y flagwidth flag for width fluctuation calculation
widthmode 1 wmode designator for width fluctuation model
compound y flagcomp flag for compound nucleus model
fullhf n flagfullhf flag for full spin dependence of transmission coefficients
eciscompound n flageciscomp flag for compound nucleus calculation by ECIS
cpang n flagcpang flag for compound angular distribution calculation for incident charged particles
urr n flagurr flag for URR calculation
urrnjoy n flagurrnjoy flag for normalization of URR parameters with NJOY method
lurr 2 lurr maximal orbital angular momentum for URR
#
# Gamma emission
#
gammax 2 gammax number of l-values for gamma multipolarity
strength 9 strength model for E1 gamma-ray strength function
strengthM1 3 strengthM1 model for M1 gamma-ray strength function
electronconv y flagelectron flag for application of electron conversion coefficient
racap y flagracap flag for radiative capture model
ldmodelracap 1 ldmodelracap level density model for direct radiative capture
upbend y flagupbend flag for low-energy upbend of photon strength function
psfglobal n flagpsfglobal flag for global photon strength functions only
#
# Pre-equilibrium
#
preequilibrium n flagpreeq flag for pre-equilibrium calculation
preeqmode 2 preeqmode designator for pre-equilibrium model
multipreeq n flagmulpre flag for multiple pre-equilibrium calculation
mpreeqmode 2 mpreeqmode designator for multiple pre-equilibrium model
breakupmodel 1 breakupmodel model for break-up reaction: 1. Kalbach 2. Avrigeanu
phmodel 1 phmodel particle-hole state density model
pairmodel 2 pairmodel designator for pre-equilibrium pairing model
preeqspin 1 pespinmodel model for pre-equilibrium spin distribution
giantresonance y flaggiant flag for collective contribution from giant resonances
preeqsurface y flagsurface flag for surface effects in exciton model
preeqcomplex y flagpecomp flag for Kalbach complex particle emission model
twocomponent y flag2comp flag for two-component pre-equilibrium model
ecisdwba y flagecisdwba flag for new ECIS calculation for DWBA for MSD
onestep n flagonestep flag for continuum one-step direct only
#
# Level densities
#
ldmodel 1 ldmodelall level density model
ldmodelCN 1 ldmodelCN level density model for compound nucleus
shellmodel 1 shellmodel model for shell correction energies
kvibmodel 2 kvibmodel model for vibrational enhancement
spincutmodel 1 spincutmodel model for spin cutoff factor for ground state
asys n flagasys flag for all level density parameters a from systematics
parity n flagparity flag for non-equal parity distribution
colenhance n flagcolall flag for collective enhancement of level density for all nuclides
ctmglobal n flagctmglob flag for global CTM model (no discrete level info)
gshell n flaggshell flag for energy dependence of single particle level density parameter g
colldamp n flagcolldamp flag for damping of collective effects in effective level density
#
# Fission
#
fission n flagfission flag for fission
fismodel 3 fismodel fission model
fismodelalt 3 fismodelalt alternative fission model for default barriers
hbstate y flaghbstate flag for head band states in fission
class2 y flagclass2 flag for class2 states in fission
fispartdamp n flagfispartdamp flag for fission partial damping
massdis n flagmassdis flag for calculation of fission fragment mass yields
ffevaporation y flagffevap flag for calculation of particle evaporation from fission fragment mass yields
fisfeed n flagfisfeed flag for output of fission per excitation bin
fymodel 2 fymodel fission yield model, 1: Brosa 2: GEF
ffmodel 1 ffmodel fission fragment model, 1: GEF 2: HF3D (Okumura) 3: SPY
pfnsmodel 1 pfnsmodel PFNS model, 1: Iwamoto 2: from FF decay
ffspin n flagffspin flag to use spin distribution in initial FF population
#
# Output
#
outmain y flagmain flag for main output
outbasic y flagbasic flag for output of basic information and results
outpopulation y flagpop flag for output of population
outcheck y flagcheck flag for output of numerical checks
outlevels y flaglevels flag for output of discrete level information
outdensity y flagdensity flag for output of level densities
outomp y flagoutomp flag for output of optical model parameters
outkd n flagoutkd flag for output of KD03 OMP parameters
outdirect y flagdirect flag for output of direct reaction results
outinverse y flaginverse flag for output of transmission coefficients and inverse reaction cross sections
outdecay n flagdecay flag for output of decay of each population bin
outtransenergy y flagtransen flag for output of transmission coefficients per energy
outecis n flagoutecis flag for output of ECIS results
outgamma y flaggamma flag for output of gamma-ray information
outpreequilibrium n flagpeout flag for output of pre-equilibrium results
outfission n flagfisout flag for output of fission information
outdiscrete y flagdisc flag for output of discrete state cross sections
outspectra n flagspec flag for output of double-differential cross sections
outbinspectra n flagbinspec flag for output of emission spectrum per excitation bin
resonance n flagres flag for output of low energy resonance cross sections
group n flaggroup flag for output of low energy groupwise cross sections
addiscrete y flagadd flag for addition of discrete states to spectra
addelastic y flagaddel flag for addition of elastic peak to spectra
outangle n flagang flag for output of angular distributions
outlegendre n flaglegendre flag for output of Legendre coefficients
ddxmode 0 ddxmode mode for double-differential cross sections
outdwba n flagoutdwba flag for output of DWBA cross sections for MSD
outgamdis n flaggamdis flag for output of discrete gamma-ray intensities
outexcitation n flagexc flag for output of excitation functions
components n flagcompo flag for output of cross section components
endf n flagendf flag for information for ENDF-6 file
endfdetail y flagendfdet flag for detailed ENDF-6 information per channel
sacs n flagsacs flag for statistical analysis of cross sections
partable n flagpartable flag for output of model parameters on separate file
block n flagblock flag to block spectra, angle and gamma files
########## BASIC REACTION PARAMETERS ##########
Projectile : neutron Mass in a.m.u. : 1.008665
Target : 43Ca Mass in a.m.u. : 42.958766
Included channels:
gamma
neutron
proton
deuteron
triton
helium-3
alpha
82 incident energies (LAB):
2.047E-06
4.094E-06
6.141E-06
8.188E-06
1.023E-05
2.047E-05
3.070E-05
4.094E-05
5.117E-05
6.141E-05
7.164E-05
8.188E-05
9.211E-05
1.023E-04
2.047E-04
3.070E-04
4.094E-04
5.117E-04
6.141E-04
7.164E-04
8.188E-04
9.211E-04
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.010
0.014
0.018
0.023
0.027
0.031
0.035
0.039
0.043
0.047
0.051
0.055
0.059
0.063
... 57247 more lines ...
|
| Attachment 5: Ca44run
|
#!/bin/bash
rm -r Ca44/
mkdir Ca44
for((c=1;c<=1000;c++))
do
./a.out << EOF
43 3.5 # amass spin
0 1000000 #min max
0 1692.91443 3.0122e-4 #ell D S
1.02046 0.204092 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Ca44/run1_0l_$c
cp out_all2 Ca44/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
43 3.5 # amass spin
0 1000000 #min max
1 908.66852 1.4938e-4 #ell D S
0.77314 0.154628 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Ca44/run1_1l_$c
cp out_all2 Ca44/run2_1l_$c
cp out_all3 Ca44/run3_1l_$c
cp out_all4 Ca44/run4_1l_$c
rm out*
done
|
|
8
|
Sat Nov 4 13:56:02 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of Ca42+n |
1. This graph is the Gaussian fitted Ca42+n with rebinned to 10 total. Theorotical cross section is around 15 while the experimental is 16.5, which matches well.
2. This the the distribution of l = 0 with bins = 10.
3. This the the distribution of l = 1 with bins = 10.
4. This is the TALYS output for Ca42+n.
5. mygoe.f input code which consists of TALYS parameters (Gg, D, S). |
| Attachment 1: macs1wrg.pdf
|
|
| Attachment 2: macs01wbin.pdf
|
|
| Attachment 3: macs11wbin.pdf
|
|
| Attachment 4: talys.out
|
TALYS-1.96 (Version: December 30, 2021)
Copyright (C) 2021 A.J. Koning, S. Hilaire and S. Goriely
Dimensions - Cross sections: mb, Energies: MeV, Angles: degrees
########## USER INPUT ##########
USER INPUT FILE
#
# general
#
projectile n
element ca
mass 42
energy 1.
astroe 0.02
racap y
resonance y
#
# output
#
#channels y
outgamma y
outdensity y
outlevels y
outbasic y
USER INPUT FILE + DEFAULTS
Keyword Value Variable Explanation
#
# Four main keywords
#
projectile n ptype0 type of incident particle
element Ca Starget symbol of target nucleus
mass 42 mass mass number of target nucleus
energy energyfile file with incident energies
#
# Basic physical and numerical parameters
#
ejectiles g n p d t h a outtype outgoing particles
maxz 12 maxZ maximal number of protons away from the initial compound nucleus
maxn 20 maxN maximal number of neutrons away from the initial compound nucleus
bins 40 nbins number of continuum excitation energy bins
equidistant y flagequi flag to use equidistant excitation bins instead of logarithmic bins
equispec n flagequispec flag to use equidistant bins for emission spectra
popmev n flagpopmev flag to use initial population per MeV instead of histograms
segment 1 segment number of segments to divide emission energy grid
maxlevelstar 30 nlevmax maximum number of included discrete levels for target
maxlevelsres 10 nlevmaxres maximum number of included discrete levels for residual nucleus
maxlevelsbin g 10 nlevbin maximum number of included discrete levels for gamma channel
maxlevelsbin n 30 nlevbin maximum number of included discrete levels for neutron channel
maxlevelsbin p 10 nlevbin maximum number of included discrete levels for proton channel
maxlevelsbin d 5 nlevbin maximum number of included discrete levels for deuteron channel
maxlevelsbin t 5 nlevbin maximum number of included discrete levels for triton channel
maxlevelsbin h 5 nlevbin maximum number of included discrete levels for helium-3 channel
maxlevelsbin a 10 nlevbin maximum number of included discrete levels for alpha channel
ltarget 0 ltarget excited level of target
isomer 1.00E+00 isomer definition of isomer in seconds
transpower 15 transpower power for transmission coefficient limit
transeps 1.00E-18 transeps limit for transmission coefficient
xseps 1.00E-17 xseps limit for cross sections
popeps 1.00E-13 popeps limit for population cross section per nucleus
Rfiseps 1.00E-06 Rfiseps ratio for limit for fission cross section per nucleus
elow 1.00E-06 elow minimal incident energy for nuclear model calculations
angles 90 nangle number of angles
anglescont 18 nanglecont number of angles for continuum
anglesrec 1 nanglerec number of recoil angles
maxenrec 20 maxenrec number of recoil energies
channels n flagchannels flag for exclusive channels calculation
maxchannel 4 maxchannel maximal number of outgoing particles in individual channel description
micro n flagmicro flag for completely microscopic Talys calculation
best n flagbest flag to use best set of adjusted parameters
bestbranch y flagbestbr flag to use flag to use only best set of branching ratios
bestend n flagbestend flag to put best set of parameters at end of input file
relativistic y flagrel flag for relativistic kinematics
recoil n flagrecoil flag for calculation of recoils
labddx n flaglabddx flag for calculation of DDX in LAB system
recoilaverage n flagrecoilav flag for average velocity in recoil calculation
channelenergy n flagEchannel flag for channel energy for emission spectrum
reaction y flagreaction flag for calculation of nuclear reactions
ngfit n flagngfit flag for using fitted (n,g) nuclear model parameters
nnfit n flagnnfit flag for using fitted (n,n), (n,2n) and (n,p) nuclear model parameters
nafit n flagnafit flag for using fitted (n,a) nuclear model parameters
astro y flagastro flag for calculation of astrophysics reaction rate
astrogs y flagastrogs flag for calculation of astrophysics reaction rate with target in ground state only
astroex n flagastroex flag for calculation of astrophysics reaction rate to long-lived excited states
nonthermlev -1 nonthermlev excited level non-thermalized in the calculation of astrophysics rate
massmodel 2 massmodel model for theoretical nuclear mass
expmass y flagexpmass flag for using experimental nuclear mass if available
disctable 1 disctable table with discrete levels
production n flagprod flag for isotope production
outfy n flagoutfy flag for output detailed fission yield calculation
gefran 50000 gefran number of random events for GEF calculation
Estop 1000.000 Estop incident energy above which TALYS stops
rpevap n flagrpevap flag for evaporation of residual products at high incident energies
maxZrp 12 maxZrp maximal number of protons away from the initial compound nucleus before residual evaporation
maxNrp 32 maxNrp maximal number of neutrons away from the initial compound nucleus before residual evaporation
#
# Optical model
#
localomp y flaglocalomp flag for local (y) or global (n) optical model
dispersion n flagdisp flag for dispersive optical model
jlmomp n flagjlm flag for using semi-microscopic JLM OMP
riplomp n flagriplomp flag for RIPL OMP
riplrisk n flagriplrisk flag for going outside RIPL mass validity range
optmodall n flagompall flag for new optical model calculation for all residual nuclei
incadjust y flagincadj flag for OMP adjustment on incident channel also
omponly n flagomponly flag to execute ONLY an optical model calculation
autorot n flagautorot flag for automatic rotational coupled channels calculations for A > 150
spherical n flagspher flag to force spherical optical model
soukho y flagsoukho flag for Soukhovitskii OMP for actinides
coulomb y flagcoulomb flag for Coulomb excitation calculation with ECIS
statepot n flagstate flag for optical model potential for each excited state
maxband 0 maxband highest vibrational band added to rotational model
maxrot 2 maxrot number of included excited rotational levels
sysreaction sysreaction particles with reaction cross section from systematics
rotational n rotational particles with possible rotational optical model
core -1 core even-even core for weakcoupling (-1 or 1)
ecissave n flagecissave flag for saving ECIS input and output files
eciscalc y flageciscalc flag for new ECIS calculation for outgoing particles and energy grid
inccalc y flaginccalc flag for new ECIS calculation for incident channel
endfecis y flagendfecis flag for new ECIS calculation for ENDF-6 files
radialmodel 2 radialmodel model for radial matter densities (JLM OMP only)
jlmmode 0 jlmmode option for JLM imaginary potential normalization
alphaomp 6 alphaomp alpha OMP (1=normal, 2= McFadden-Satchler, 3-5= folding potential, 6,8= Avrigeanu, 7=Nolte)
deuteronomp 1 deuteronomp deuteron OMP (1=normal, 2=Daehnick, 3=Bojowald, 4=Han-Shi-Shen, 5=An-Cai)
#
# Compound nucleus
#
widthfluc y flagwidth flag for width fluctuation calculation
widthmode 1 wmode designator for width fluctuation model
compound y flagcomp flag for compound nucleus model
fullhf n flagfullhf flag for full spin dependence of transmission coefficients
eciscompound n flageciscomp flag for compound nucleus calculation by ECIS
cpang n flagcpang flag for compound angular distribution calculation for incident charged particles
urr n flagurr flag for URR calculation
urrnjoy n flagurrnjoy flag for normalization of URR parameters with NJOY method
lurr 2 lurr maximal orbital angular momentum for URR
#
# Gamma emission
#
gammax 2 gammax number of l-values for gamma multipolarity
strength 9 strength model for E1 gamma-ray strength function
strengthM1 3 strengthM1 model for M1 gamma-ray strength function
electronconv y flagelectron flag for application of electron conversion coefficient
racap y flagracap flag for radiative capture model
ldmodelracap 1 ldmodelracap level density model for direct radiative capture
upbend y flagupbend flag for low-energy upbend of photon strength function
psfglobal n flagpsfglobal flag for global photon strength functions only
#
# Pre-equilibrium
#
preequilibrium n flagpreeq flag for pre-equilibrium calculation
preeqmode 2 preeqmode designator for pre-equilibrium model
multipreeq n flagmulpre flag for multiple pre-equilibrium calculation
mpreeqmode 2 mpreeqmode designator for multiple pre-equilibrium model
breakupmodel 1 breakupmodel model for break-up reaction: 1. Kalbach 2. Avrigeanu
phmodel 1 phmodel particle-hole state density model
pairmodel 2 pairmodel designator for pre-equilibrium pairing model
preeqspin 1 pespinmodel model for pre-equilibrium spin distribution
giantresonance y flaggiant flag for collective contribution from giant resonances
preeqsurface y flagsurface flag for surface effects in exciton model
preeqcomplex y flagpecomp flag for Kalbach complex particle emission model
twocomponent y flag2comp flag for two-component pre-equilibrium model
ecisdwba y flagecisdwba flag for new ECIS calculation for DWBA for MSD
onestep n flagonestep flag for continuum one-step direct only
#
# Level densities
#
ldmodel 1 ldmodelall level density model
ldmodelCN 1 ldmodelCN level density model for compound nucleus
shellmodel 1 shellmodel model for shell correction energies
kvibmodel 2 kvibmodel model for vibrational enhancement
spincutmodel 1 spincutmodel model for spin cutoff factor for ground state
asys n flagasys flag for all level density parameters a from systematics
parity n flagparity flag for non-equal parity distribution
colenhance n flagcolall flag for collective enhancement of level density for all nuclides
ctmglobal n flagctmglob flag for global CTM model (no discrete level info)
gshell n flaggshell flag for energy dependence of single particle level density parameter g
colldamp n flagcolldamp flag for damping of collective effects in effective level density
#
# Fission
#
fission n flagfission flag for fission
fismodel 3 fismodel fission model
fismodelalt 3 fismodelalt alternative fission model for default barriers
hbstate y flaghbstate flag for head band states in fission
class2 y flagclass2 flag for class2 states in fission
fispartdamp n flagfispartdamp flag for fission partial damping
massdis n flagmassdis flag for calculation of fission fragment mass yields
ffevaporation y flagffevap flag for calculation of particle evaporation from fission fragment mass yields
fisfeed n flagfisfeed flag for output of fission per excitation bin
fymodel 2 fymodel fission yield model, 1: Brosa 2: GEF
ffmodel 1 ffmodel fission fragment model, 1: GEF 2: HF3D (Okumura) 3: SPY
pfnsmodel 1 pfnsmodel PFNS model, 1: Iwamoto 2: from FF decay
ffspin n flagffspin flag to use spin distribution in initial FF population
#
# Output
#
outmain y flagmain flag for main output
outbasic y flagbasic flag for output of basic information and results
outpopulation y flagpop flag for output of population
outcheck y flagcheck flag for output of numerical checks
outlevels y flaglevels flag for output of discrete level information
outdensity y flagdensity flag for output of level densities
outomp y flagoutomp flag for output of optical model parameters
outkd n flagoutkd flag for output of KD03 OMP parameters
outdirect y flagdirect flag for output of direct reaction results
outinverse y flaginverse flag for output of transmission coefficients and inverse reaction cross sections
outdecay n flagdecay flag for output of decay of each population bin
outtransenergy y flagtransen flag for output of transmission coefficients per energy
outecis n flagoutecis flag for output of ECIS results
outgamma y flaggamma flag for output of gamma-ray information
outpreequilibrium n flagpeout flag for output of pre-equilibrium results
outfission n flagfisout flag for output of fission information
outdiscrete y flagdisc flag for output of discrete state cross sections
outspectra n flagspec flag for output of double-differential cross sections
outbinspectra n flagbinspec flag for output of emission spectrum per excitation bin
resonance y flagres flag for output of low energy resonance cross sections
group n flaggroup flag for output of low energy groupwise cross sections
addiscrete y flagadd flag for addition of discrete states to spectra
addelastic y flagaddel flag for addition of elastic peak to spectra
outangle n flagang flag for output of angular distributions
outlegendre n flaglegendre flag for output of Legendre coefficients
ddxmode 0 ddxmode mode for double-differential cross sections
outdwba n flagoutdwba flag for output of DWBA cross sections for MSD
outgamdis n flaggamdis flag for output of discrete gamma-ray intensities
outexcitation n flagexc flag for output of excitation functions
components n flagcompo flag for output of cross section components
endf n flagendf flag for information for ENDF-6 file
endfdetail y flagendfdet flag for detailed ENDF-6 information per channel
sacs n flagsacs flag for statistical analysis of cross sections
partable n flagpartable flag for output of model parameters on separate file
block n flagblock flag to block spectra, angle and gamma files
########## BASIC REACTION PARAMETERS ##########
Projectile : neutron Mass in a.m.u. : 1.008665
Target : 42Ca Mass in a.m.u. : 41.958618
Included channels:
gamma
neutron
proton
deuteron
triton
helium-3
alpha
73 incident energies (LAB):
2.048E-06
4.096E-06
6.144E-06
8.192E-06
1.024E-05
2.048E-05
3.072E-05
4.096E-05
5.120E-05
6.144E-05
7.168E-05
8.192E-05
9.216E-05
1.024E-04
2.048E-04
3.072E-04
4.096E-04
5.120E-04
6.144E-04
7.168E-04
8.192E-04
9.216E-04
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.010
0.014
0.018
0.023
0.027
0.031
0.035
0.039
0.043
0.047
0.051
0.055
0.059
... 39897 more lines ...
|
| Attachment 5: Ca43run
|
#!/bin/bash
rm -r Ca43/
mkdir Ca43
for((c=1;c<=1000;c++))
do
./a.out << EOF
42 0 # amass spin
0 1000000 #min max
0 12888.49219 2.3936e-4 #ell D S
0.23701 0.047402 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Ca43/run1_0l_$c
cp out_all2 Ca43/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
42 0 # amass spin
0 1000000 #min max
1 4857.46973 1.9670e-4 #ell D S
0.23342 0.046684 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Ca43/run1_1l_$c
cp out_all2 Ca43/run2_1l_$c
cp out_all3 Ca43/run3_1l_$c
cp out_all4 Ca43/run4_1l_$c
rm out*
done
|
|
7
|
Sat Nov 4 13:45:53 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of P31+n |
1. This graph is fitted P31+n with rebinned to 10. Theorotical cross section is around 2 and the experimental value is 1.74, which relatively match.
2. This the the distribution of l = 0 with bins = 10.
3. This the the distribution of l = 1 with bins = 10.
4. This is the TALYS output for P31+n.
5. mygoe.f input code which consists of TALYS parameters (Gg, D, S). |
| Attachment 1: macs1wbg.pdf
|
|
| Attachment 2: macs01.pdf
|
|
| Attachment 3: macs11wbin.pdf
|
|
| Attachment 4: talys.out
|
TALYS-1.96 (Version: December 30, 2021)
Copyright (C) 2021 A.J. Koning, S. Hilaire and S. Goriely
Dimensions - Cross sections: mb, Energies: MeV, Angles: degrees
########## USER INPUT ##########
USER INPUT FILE
#
# general
#
projectile n
element p
mass 31
energy 1.
astroe 0.02
racap y
#
# output
#
#channels y
outgamma y
outdensity y
outlevels y
outbasic y
USER INPUT FILE + DEFAULTS
Keyword Value Variable Explanation
#
# Four main keywords
#
projectile n ptype0 type of incident particle
element P Starget symbol of target nucleus
mass 31 mass mass number of target nucleus
energy energyfile file with incident energies
#
# Basic physical and numerical parameters
#
ejectiles g n p d t h a outtype outgoing particles
maxz 12 maxZ maximal number of protons away from the initial compound nucleus
maxn 14 maxN maximal number of neutrons away from the initial compound nucleus
bins 40 nbins number of continuum excitation energy bins
equidistant y flagequi flag to use equidistant excitation bins instead of logarithmic bins
equispec n flagequispec flag to use equidistant bins for emission spectra
popmev n flagpopmev flag to use initial population per MeV instead of histograms
segment 1 segment number of segments to divide emission energy grid
maxlevelstar 30 nlevmax maximum number of included discrete levels for target
maxlevelsres 10 nlevmaxres maximum number of included discrete levels for residual nucleus
maxlevelsbin g 10 nlevbin maximum number of included discrete levels for gamma channel
maxlevelsbin n 30 nlevbin maximum number of included discrete levels for neutron channel
maxlevelsbin p 10 nlevbin maximum number of included discrete levels for proton channel
maxlevelsbin d 5 nlevbin maximum number of included discrete levels for deuteron channel
maxlevelsbin t 5 nlevbin maximum number of included discrete levels for triton channel
maxlevelsbin h 5 nlevbin maximum number of included discrete levels for helium-3 channel
maxlevelsbin a 10 nlevbin maximum number of included discrete levels for alpha channel
ltarget 0 ltarget excited level of target
isomer 1.00E+00 isomer definition of isomer in seconds
transpower 15 transpower power for transmission coefficient limit
transeps 1.00E-18 transeps limit for transmission coefficient
xseps 1.00E-17 xseps limit for cross sections
popeps 1.00E-13 popeps limit for population cross section per nucleus
Rfiseps 1.00E-06 Rfiseps ratio for limit for fission cross section per nucleus
elow 1.00E-06 elow minimal incident energy for nuclear model calculations
angles 90 nangle number of angles
anglescont 18 nanglecont number of angles for continuum
anglesrec 1 nanglerec number of recoil angles
maxenrec 20 maxenrec number of recoil energies
channels n flagchannels flag for exclusive channels calculation
maxchannel 4 maxchannel maximal number of outgoing particles in individual channel description
micro n flagmicro flag for completely microscopic Talys calculation
best n flagbest flag to use best set of adjusted parameters
bestbranch y flagbestbr flag to use flag to use only best set of branching ratios
bestend n flagbestend flag to put best set of parameters at end of input file
relativistic y flagrel flag for relativistic kinematics
recoil n flagrecoil flag for calculation of recoils
labddx n flaglabddx flag for calculation of DDX in LAB system
recoilaverage n flagrecoilav flag for average velocity in recoil calculation
channelenergy n flagEchannel flag for channel energy for emission spectrum
reaction y flagreaction flag for calculation of nuclear reactions
ngfit n flagngfit flag for using fitted (n,g) nuclear model parameters
nnfit n flagnnfit flag for using fitted (n,n), (n,2n) and (n,p) nuclear model parameters
nafit n flagnafit flag for using fitted (n,a) nuclear model parameters
astro y flagastro flag for calculation of astrophysics reaction rate
astrogs y flagastrogs flag for calculation of astrophysics reaction rate with target in ground state only
astroex n flagastroex flag for calculation of astrophysics reaction rate to long-lived excited states
nonthermlev -1 nonthermlev excited level non-thermalized in the calculation of astrophysics rate
massmodel 2 massmodel model for theoretical nuclear mass
expmass y flagexpmass flag for using experimental nuclear mass if available
disctable 1 disctable table with discrete levels
production n flagprod flag for isotope production
outfy n flagoutfy flag for output detailed fission yield calculation
gefran 50000 gefran number of random events for GEF calculation
Estop 1000.000 Estop incident energy above which TALYS stops
rpevap n flagrpevap flag for evaporation of residual products at high incident energies
maxZrp 12 maxZrp maximal number of protons away from the initial compound nucleus before residual evaporation
maxNrp 32 maxNrp maximal number of neutrons away from the initial compound nucleus before residual evaporation
#
# Optical model
#
localomp y flaglocalomp flag for local (y) or global (n) optical model
dispersion n flagdisp flag for dispersive optical model
jlmomp n flagjlm flag for using semi-microscopic JLM OMP
riplomp n flagriplomp flag for RIPL OMP
riplrisk n flagriplrisk flag for going outside RIPL mass validity range
optmodall n flagompall flag for new optical model calculation for all residual nuclei
incadjust y flagincadj flag for OMP adjustment on incident channel also
omponly n flagomponly flag to execute ONLY an optical model calculation
autorot n flagautorot flag for automatic rotational coupled channels calculations for A > 150
spherical n flagspher flag to force spherical optical model
soukho y flagsoukho flag for Soukhovitskii OMP for actinides
coulomb y flagcoulomb flag for Coulomb excitation calculation with ECIS
statepot n flagstate flag for optical model potential for each excited state
maxband 0 maxband highest vibrational band added to rotational model
maxrot 2 maxrot number of included excited rotational levels
sysreaction sysreaction particles with reaction cross section from systematics
rotational n rotational particles with possible rotational optical model
core -1 core even-even core for weakcoupling (-1 or 1)
ecissave n flagecissave flag for saving ECIS input and output files
eciscalc y flageciscalc flag for new ECIS calculation for outgoing particles and energy grid
inccalc y flaginccalc flag for new ECIS calculation for incident channel
endfecis y flagendfecis flag for new ECIS calculation for ENDF-6 files
radialmodel 2 radialmodel model for radial matter densities (JLM OMP only)
jlmmode 0 jlmmode option for JLM imaginary potential normalization
alphaomp 6 alphaomp alpha OMP (1=normal, 2= McFadden-Satchler, 3-5= folding potential, 6,8= Avrigeanu, 7=Nolte)
deuteronomp 1 deuteronomp deuteron OMP (1=normal, 2=Daehnick, 3=Bojowald, 4=Han-Shi-Shen, 5=An-Cai)
#
# Compound nucleus
#
widthfluc y flagwidth flag for width fluctuation calculation
widthmode 1 wmode designator for width fluctuation model
compound y flagcomp flag for compound nucleus model
fullhf n flagfullhf flag for full spin dependence of transmission coefficients
eciscompound n flageciscomp flag for compound nucleus calculation by ECIS
cpang n flagcpang flag for compound angular distribution calculation for incident charged particles
urr n flagurr flag for URR calculation
urrnjoy n flagurrnjoy flag for normalization of URR parameters with NJOY method
lurr 2 lurr maximal orbital angular momentum for URR
#
# Gamma emission
#
gammax 2 gammax number of l-values for gamma multipolarity
strength 9 strength model for E1 gamma-ray strength function
strengthM1 3 strengthM1 model for M1 gamma-ray strength function
electronconv y flagelectron flag for application of electron conversion coefficient
racap y flagracap flag for radiative capture model
ldmodelracap 1 ldmodelracap level density model for direct radiative capture
upbend y flagupbend flag for low-energy upbend of photon strength function
psfglobal n flagpsfglobal flag for global photon strength functions only
#
# Pre-equilibrium
#
preequilibrium n flagpreeq flag for pre-equilibrium calculation
preeqmode 2 preeqmode designator for pre-equilibrium model
multipreeq n flagmulpre flag for multiple pre-equilibrium calculation
mpreeqmode 2 mpreeqmode designator for multiple pre-equilibrium model
breakupmodel 1 breakupmodel model for break-up reaction: 1. Kalbach 2. Avrigeanu
phmodel 1 phmodel particle-hole state density model
pairmodel 2 pairmodel designator for pre-equilibrium pairing model
preeqspin 1 pespinmodel model for pre-equilibrium spin distribution
giantresonance y flaggiant flag for collective contribution from giant resonances
preeqsurface y flagsurface flag for surface effects in exciton model
preeqcomplex y flagpecomp flag for Kalbach complex particle emission model
twocomponent y flag2comp flag for two-component pre-equilibrium model
ecisdwba y flagecisdwba flag for new ECIS calculation for DWBA for MSD
onestep n flagonestep flag for continuum one-step direct only
#
# Level densities
#
ldmodel 1 ldmodelall level density model
ldmodelCN 1 ldmodelCN level density model for compound nucleus
shellmodel 1 shellmodel model for shell correction energies
kvibmodel 2 kvibmodel model for vibrational enhancement
spincutmodel 1 spincutmodel model for spin cutoff factor for ground state
asys n flagasys flag for all level density parameters a from systematics
parity n flagparity flag for non-equal parity distribution
colenhance n flagcolall flag for collective enhancement of level density for all nuclides
ctmglobal n flagctmglob flag for global CTM model (no discrete level info)
gshell n flaggshell flag for energy dependence of single particle level density parameter g
colldamp n flagcolldamp flag for damping of collective effects in effective level density
#
# Fission
#
fission n flagfission flag for fission
fismodel 3 fismodel fission model
fismodelalt 3 fismodelalt alternative fission model for default barriers
hbstate y flaghbstate flag for head band states in fission
class2 y flagclass2 flag for class2 states in fission
fispartdamp n flagfispartdamp flag for fission partial damping
massdis n flagmassdis flag for calculation of fission fragment mass yields
ffevaporation y flagffevap flag for calculation of particle evaporation from fission fragment mass yields
fisfeed n flagfisfeed flag for output of fission per excitation bin
fymodel 2 fymodel fission yield model, 1: Brosa 2: GEF
ffmodel 1 ffmodel fission fragment model, 1: GEF 2: HF3D (Okumura) 3: SPY
pfnsmodel 1 pfnsmodel PFNS model, 1: Iwamoto 2: from FF decay
ffspin n flagffspin flag to use spin distribution in initial FF population
#
# Output
#
outmain y flagmain flag for main output
outbasic y flagbasic flag for output of basic information and results
outpopulation y flagpop flag for output of population
outcheck y flagcheck flag for output of numerical checks
outlevels y flaglevels flag for output of discrete level information
outdensity y flagdensity flag for output of level densities
outomp y flagoutomp flag for output of optical model parameters
outkd n flagoutkd flag for output of KD03 OMP parameters
outdirect y flagdirect flag for output of direct reaction results
outinverse y flaginverse flag for output of transmission coefficients and inverse reaction cross sections
outdecay n flagdecay flag for output of decay of each population bin
outtransenergy y flagtransen flag for output of transmission coefficients per energy
outecis n flagoutecis flag for output of ECIS results
outgamma y flaggamma flag for output of gamma-ray information
outpreequilibrium n flagpeout flag for output of pre-equilibrium results
outfission n flagfisout flag for output of fission information
outdiscrete y flagdisc flag for output of discrete state cross sections
outspectra n flagspec flag for output of double-differential cross sections
outbinspectra n flagbinspec flag for output of emission spectrum per excitation bin
resonance n flagres flag for output of low energy resonance cross sections
group n flaggroup flag for output of low energy groupwise cross sections
addiscrete y flagadd flag for addition of discrete states to spectra
addelastic y flagaddel flag for addition of elastic peak to spectra
outangle n flagang flag for output of angular distributions
outlegendre n flaglegendre flag for output of Legendre coefficients
ddxmode 0 ddxmode mode for double-differential cross sections
outdwba n flagoutdwba flag for output of DWBA cross sections for MSD
outgamdis n flaggamdis flag for output of discrete gamma-ray intensities
outexcitation n flagexc flag for output of excitation functions
components n flagcompo flag for output of cross section components
endf n flagendf flag for information for ENDF-6 file
endfdetail y flagendfdet flag for detailed ENDF-6 information per channel
sacs n flagsacs flag for statistical analysis of cross sections
partable n flagpartable flag for output of model parameters on separate file
block n flagblock flag to block spectra, angle and gamma files
########## BASIC REACTION PARAMETERS ##########
Projectile : neutron Mass in a.m.u. : 1.008665
Target : 31P Mass in a.m.u. : 30.973762
Included channels:
gamma
neutron
proton
deuteron
triton
helium-3
alpha
102 incident energies (LAB):
2.065E-06
4.130E-06
6.195E-06
8.261E-06
1.033E-05
2.065E-05
3.098E-05
4.130E-05
5.163E-05
6.195E-05
7.228E-05
8.261E-05
9.293E-05
1.033E-04
2.065E-04
3.098E-04
4.130E-04
5.163E-04
6.195E-04
7.228E-04
8.261E-04
9.293E-04
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.010
0.014
0.019
0.023
0.027
0.031
0.035
0.039
0.043
0.047
0.052
0.056
0.060
0.064
... 56676 more lines ...
|
| Attachment 5: P32run
|
#!/bin/bash
rm -r P32/
mkdir P32
for((c=1;c<=1000;c++))
do
./a.out << EOF
31 0.5 # amass spin
0 1000000 #min max
0 62044.46094 0.8401e-4 #ell D S
0.69793 0.139586 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 P32/run1_0l_$c
cp out_all2 P32/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
31 0.5 # amass spin
0 1000000 #min max
1 33466.96484 6.1765e-4 #ell D S
1.25400 0.2508 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 P32/run1_1l_$c
cp out_all2 P32/run2_1l_$c
cp out_all3 P32/run3_1l_$c
cp out_all4 P32/run4_1l_$c
rm out*
done
|
|
6
|
Sat Nov 4 13:40:39 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of Si28+n |
1. This graph is the exponential decay Si28+n total. Theorotical cross section is between 0-1.5 while the experimental should be towards 1.15
2. This the the distribution of l = 0.
3. This the the distribution of l = 1.
4. This is the TALYS output for Si28+n.
5. mygoe.f input code which consists of TALYS parameters (Gg, D, S). |
| Attachment 1: macs1fit.pdf
|
|
| Attachment 2: macs01.pdf
|
|
| Attachment 3: macs11.pdf
|
|
| Attachment 4: talys.out
|
TALYS-1.96 (Version: December 30, 2021)
Copyright (C) 2021 A.J. Koning, S. Hilaire and S. Goriely
Dimensions - Cross sections: mb, Energies: MeV, Angles: degrees
########## USER INPUT ##########
USER INPUT FILE
#
# general
#
projectile n
element si
mass 28
energy 1.
astroe 0.02
racap y
#
# output
#
#channels y
outgamma y
outdensity y
outlevels y
outbasic y
USER INPUT FILE + DEFAULTS
Keyword Value Variable Explanation
#
# Four main keywords
#
projectile n ptype0 type of incident particle
element Si Starget symbol of target nucleus
mass 28 mass mass number of target nucleus
energy energyfile file with incident energies
#
# Basic physical and numerical parameters
#
ejectiles g n p d t h a outtype outgoing particles
maxz 11 maxZ maximal number of protons away from the initial compound nucleus
maxn 12 maxN maximal number of neutrons away from the initial compound nucleus
bins 40 nbins number of continuum excitation energy bins
equidistant y flagequi flag to use equidistant excitation bins instead of logarithmic bins
equispec n flagequispec flag to use equidistant bins for emission spectra
popmev n flagpopmev flag to use initial population per MeV instead of histograms
segment 1 segment number of segments to divide emission energy grid
maxlevelstar 30 nlevmax maximum number of included discrete levels for target
maxlevelsres 10 nlevmaxres maximum number of included discrete levels for residual nucleus
maxlevelsbin g 10 nlevbin maximum number of included discrete levels for gamma channel
maxlevelsbin n 30 nlevbin maximum number of included discrete levels for neutron channel
maxlevelsbin p 10 nlevbin maximum number of included discrete levels for proton channel
maxlevelsbin d 5 nlevbin maximum number of included discrete levels for deuteron channel
maxlevelsbin t 5 nlevbin maximum number of included discrete levels for triton channel
maxlevelsbin h 5 nlevbin maximum number of included discrete levels for helium-3 channel
maxlevelsbin a 10 nlevbin maximum number of included discrete levels for alpha channel
ltarget 0 ltarget excited level of target
isomer 1.00E+00 isomer definition of isomer in seconds
transpower 15 transpower power for transmission coefficient limit
transeps 1.00E-18 transeps limit for transmission coefficient
xseps 1.00E-17 xseps limit for cross sections
popeps 1.00E-13 popeps limit for population cross section per nucleus
Rfiseps 1.00E-06 Rfiseps ratio for limit for fission cross section per nucleus
elow 1.00E-06 elow minimal incident energy for nuclear model calculations
angles 90 nangle number of angles
anglescont 18 nanglecont number of angles for continuum
anglesrec 1 nanglerec number of recoil angles
maxenrec 20 maxenrec number of recoil energies
channels n flagchannels flag for exclusive channels calculation
maxchannel 4 maxchannel maximal number of outgoing particles in individual channel description
micro n flagmicro flag for completely microscopic Talys calculation
best n flagbest flag to use best set of adjusted parameters
bestbranch y flagbestbr flag to use flag to use only best set of branching ratios
bestend n flagbestend flag to put best set of parameters at end of input file
relativistic y flagrel flag for relativistic kinematics
recoil n flagrecoil flag for calculation of recoils
labddx n flaglabddx flag for calculation of DDX in LAB system
recoilaverage n flagrecoilav flag for average velocity in recoil calculation
channelenergy n flagEchannel flag for channel energy for emission spectrum
reaction y flagreaction flag for calculation of nuclear reactions
ngfit n flagngfit flag for using fitted (n,g) nuclear model parameters
nnfit n flagnnfit flag for using fitted (n,n), (n,2n) and (n,p) nuclear model parameters
nafit n flagnafit flag for using fitted (n,a) nuclear model parameters
astro y flagastro flag for calculation of astrophysics reaction rate
astrogs y flagastrogs flag for calculation of astrophysics reaction rate with target in ground state only
astroex n flagastroex flag for calculation of astrophysics reaction rate to long-lived excited states
nonthermlev -1 nonthermlev excited level non-thermalized in the calculation of astrophysics rate
massmodel 2 massmodel model for theoretical nuclear mass
expmass y flagexpmass flag for using experimental nuclear mass if available
disctable 1 disctable table with discrete levels
production n flagprod flag for isotope production
outfy n flagoutfy flag for output detailed fission yield calculation
gefran 50000 gefran number of random events for GEF calculation
Estop 1000.000 Estop incident energy above which TALYS stops
rpevap n flagrpevap flag for evaporation of residual products at high incident energies
maxZrp 12 maxZrp maximal number of protons away from the initial compound nucleus before residual evaporation
maxNrp 32 maxNrp maximal number of neutrons away from the initial compound nucleus before residual evaporation
#
# Optical model
#
localomp y flaglocalomp flag for local (y) or global (n) optical model
dispersion n flagdisp flag for dispersive optical model
jlmomp n flagjlm flag for using semi-microscopic JLM OMP
riplomp n flagriplomp flag for RIPL OMP
riplrisk n flagriplrisk flag for going outside RIPL mass validity range
optmodall n flagompall flag for new optical model calculation for all residual nuclei
incadjust y flagincadj flag for OMP adjustment on incident channel also
omponly n flagomponly flag to execute ONLY an optical model calculation
autorot n flagautorot flag for automatic rotational coupled channels calculations for A > 150
spherical n flagspher flag to force spherical optical model
soukho y flagsoukho flag for Soukhovitskii OMP for actinides
coulomb y flagcoulomb flag for Coulomb excitation calculation with ECIS
statepot n flagstate flag for optical model potential for each excited state
maxband 0 maxband highest vibrational band added to rotational model
maxrot 2 maxrot number of included excited rotational levels
sysreaction sysreaction particles with reaction cross section from systematics
rotational n rotational particles with possible rotational optical model
core -1 core even-even core for weakcoupling (-1 or 1)
ecissave n flagecissave flag for saving ECIS input and output files
eciscalc y flageciscalc flag for new ECIS calculation for outgoing particles and energy grid
inccalc y flaginccalc flag for new ECIS calculation for incident channel
endfecis y flagendfecis flag for new ECIS calculation for ENDF-6 files
radialmodel 2 radialmodel model for radial matter densities (JLM OMP only)
jlmmode 0 jlmmode option for JLM imaginary potential normalization
alphaomp 6 alphaomp alpha OMP (1=normal, 2= McFadden-Satchler, 3-5= folding potential, 6,8= Avrigeanu, 7=Nolte)
deuteronomp 1 deuteronomp deuteron OMP (1=normal, 2=Daehnick, 3=Bojowald, 4=Han-Shi-Shen, 5=An-Cai)
#
# Compound nucleus
#
widthfluc y flagwidth flag for width fluctuation calculation
widthmode 1 wmode designator for width fluctuation model
compound y flagcomp flag for compound nucleus model
fullhf n flagfullhf flag for full spin dependence of transmission coefficients
eciscompound n flageciscomp flag for compound nucleus calculation by ECIS
cpang n flagcpang flag for compound angular distribution calculation for incident charged particles
urr n flagurr flag for URR calculation
urrnjoy n flagurrnjoy flag for normalization of URR parameters with NJOY method
lurr 2 lurr maximal orbital angular momentum for URR
#
# Gamma emission
#
gammax 2 gammax number of l-values for gamma multipolarity
strength 9 strength model for E1 gamma-ray strength function
strengthM1 3 strengthM1 model for M1 gamma-ray strength function
electronconv y flagelectron flag for application of electron conversion coefficient
racap y flagracap flag for radiative capture model
ldmodelracap 1 ldmodelracap level density model for direct radiative capture
upbend y flagupbend flag for low-energy upbend of photon strength function
psfglobal n flagpsfglobal flag for global photon strength functions only
#
# Pre-equilibrium
#
preequilibrium n flagpreeq flag for pre-equilibrium calculation
preeqmode 2 preeqmode designator for pre-equilibrium model
multipreeq n flagmulpre flag for multiple pre-equilibrium calculation
mpreeqmode 2 mpreeqmode designator for multiple pre-equilibrium model
breakupmodel 1 breakupmodel model for break-up reaction: 1. Kalbach 2. Avrigeanu
phmodel 1 phmodel particle-hole state density model
pairmodel 2 pairmodel designator for pre-equilibrium pairing model
preeqspin 1 pespinmodel model for pre-equilibrium spin distribution
giantresonance y flaggiant flag for collective contribution from giant resonances
preeqsurface y flagsurface flag for surface effects in exciton model
preeqcomplex y flagpecomp flag for Kalbach complex particle emission model
twocomponent y flag2comp flag for two-component pre-equilibrium model
ecisdwba y flagecisdwba flag for new ECIS calculation for DWBA for MSD
onestep n flagonestep flag for continuum one-step direct only
#
# Level densities
#
ldmodel 1 ldmodelall level density model
ldmodelCN 1 ldmodelCN level density model for compound nucleus
shellmodel 1 shellmodel model for shell correction energies
kvibmodel 2 kvibmodel model for vibrational enhancement
spincutmodel 1 spincutmodel model for spin cutoff factor for ground state
asys n flagasys flag for all level density parameters a from systematics
parity n flagparity flag for non-equal parity distribution
colenhance n flagcolall flag for collective enhancement of level density for all nuclides
ctmglobal n flagctmglob flag for global CTM model (no discrete level info)
gshell n flaggshell flag for energy dependence of single particle level density parameter g
colldamp n flagcolldamp flag for damping of collective effects in effective level density
#
# Fission
#
fission n flagfission flag for fission
fismodel 3 fismodel fission model
fismodelalt 3 fismodelalt alternative fission model for default barriers
hbstate y flaghbstate flag for head band states in fission
class2 y flagclass2 flag for class2 states in fission
fispartdamp n flagfispartdamp flag for fission partial damping
massdis n flagmassdis flag for calculation of fission fragment mass yields
ffevaporation y flagffevap flag for calculation of particle evaporation from fission fragment mass yields
fisfeed n flagfisfeed flag for output of fission per excitation bin
fymodel 2 fymodel fission yield model, 1: Brosa 2: GEF
ffmodel 1 ffmodel fission fragment model, 1: GEF 2: HF3D (Okumura) 3: SPY
pfnsmodel 1 pfnsmodel PFNS model, 1: Iwamoto 2: from FF decay
ffspin n flagffspin flag to use spin distribution in initial FF population
#
# Output
#
outmain y flagmain flag for main output
outbasic y flagbasic flag for output of basic information and results
outpopulation y flagpop flag for output of population
outcheck y flagcheck flag for output of numerical checks
outlevels y flaglevels flag for output of discrete level information
outdensity y flagdensity flag for output of level densities
outomp y flagoutomp flag for output of optical model parameters
outkd n flagoutkd flag for output of KD03 OMP parameters
outdirect y flagdirect flag for output of direct reaction results
outinverse y flaginverse flag for output of transmission coefficients and inverse reaction cross sections
outdecay n flagdecay flag for output of decay of each population bin
outtransenergy y flagtransen flag for output of transmission coefficients per energy
outecis n flagoutecis flag for output of ECIS results
outgamma y flaggamma flag for output of gamma-ray information
outpreequilibrium n flagpeout flag for output of pre-equilibrium results
outfission n flagfisout flag for output of fission information
outdiscrete y flagdisc flag for output of discrete state cross sections
outspectra n flagspec flag for output of double-differential cross sections
outbinspectra n flagbinspec flag for output of emission spectrum per excitation bin
resonance n flagres flag for output of low energy resonance cross sections
group n flaggroup flag for output of low energy groupwise cross sections
addiscrete y flagadd flag for addition of discrete states to spectra
addelastic y flagaddel flag for addition of elastic peak to spectra
outangle n flagang flag for output of angular distributions
outlegendre n flaglegendre flag for output of Legendre coefficients
ddxmode 0 ddxmode mode for double-differential cross sections
outdwba n flagoutdwba flag for output of DWBA cross sections for MSD
outgamdis n flaggamdis flag for output of discrete gamma-ray intensities
outexcitation n flagexc flag for output of excitation functions
components n flagcompo flag for output of cross section components
endf n flagendf flag for information for ENDF-6 file
endfdetail y flagendfdet flag for detailed ENDF-6 information per channel
sacs n flagsacs flag for statistical analysis of cross sections
partable n flagpartable flag for output of model parameters on separate file
block n flagblock flag to block spectra, angle and gamma files
########## BASIC REACTION PARAMETERS ##########
Projectile : neutron Mass in a.m.u. : 1.008665
Target : 28Si Mass in a.m.u. : 27.976927
Included channels:
gamma
neutron
proton
deuteron
triton
helium-3
alpha
73 incident energies (LAB):
2.072E-06
4.144E-06
6.216E-06
8.288E-06
1.036E-05
2.072E-05
3.108E-05
4.144E-05
5.180E-05
6.216E-05
7.252E-05
8.288E-05
9.324E-05
1.036E-04
2.072E-04
3.108E-04
4.144E-04
5.180E-04
6.216E-04
7.252E-04
8.288E-04
9.324E-04
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.010
0.015
0.019
0.023
0.027
0.031
0.035
0.039
0.044
0.048
0.052
0.056
0.060
0.064
... 38861 more lines ...
|
| Attachment 5: Si29run
|
#!/bin/bash
rm -r Si29/
mkdir Si29
for((c=1;c<=1000;c++))
do
./a.out << EOF
28 0 # amass spin
0 1000000 #min max
0 211419.28125 1.0356e-4 #ell D S
0.56467 0.112934 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Si29/run1_0l_$c
cp out_all2 Si29/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
28 0 # amass spin
0 1000000 #min max
1 84723.37500 5.9809e-4 #ell D S
0.98687 0.197374 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Si29/run1_1l_$c
cp out_all2 Si29/run2_1l_$c
cp out_all3 Si29/run3_1l_$c
cp out_all4 Si29/run4_1l_$c
rm out*
done
|
|
5
|
Sat Nov 4 13:18:53 2023 |
Selin Ezgi Birincioglu | Gaussian fit of Cl35+n with a binning of 10 |
1. This graph is the Gaussian fitted Cl35+n with rebinned total. Experimental cross section is around 12 while the graph shows the maximum probability as around 7.
2. This the the distribution of l = 0 with bins = 10.
3. This the the distribution of l = 1 with bins = 10.
4. This is the TALYS output for Cl35+n.
5. mygoe.f input code which consists of TALYS parameters (Gg, D, S). |
| Attachment 1: macs1wrg.pdf
|
|
| Attachment 2: macs01wbin.pdf
|
|
| Attachment 3: macs11wbin.pdf
|
|
| Attachment 4: talys.out
|
TALYS-1.96 (Version: December 30, 2021)
Copyright (C) 2021 A.J. Koning, S. Hilaire and S. Goriely
Dimensions - Cross sections: mb, Energies: MeV, Angles: degrees
########## USER INPUT ##########
USER INPUT FILE
#
# general
#
projectile n
element cl
mass 35
energy 1.
astroe 0.02
racap y
#
# output
#
#channels y
outgamma y
outdensity y
outlevels y
outbasic y
USER INPUT FILE + DEFAULTS
Keyword Value Variable Explanation
#
# Four main keywords
#
projectile n ptype0 type of incident particle
element Cl Starget symbol of target nucleus
mass 35 mass mass number of target nucleus
energy energyfile file with incident energies
#
# Basic physical and numerical parameters
#
ejectiles g n p d t h a outtype outgoing particles
maxz 12 maxZ maximal number of protons away from the initial compound nucleus
maxn 16 maxN maximal number of neutrons away from the initial compound nucleus
bins 40 nbins number of continuum excitation energy bins
equidistant y flagequi flag to use equidistant excitation bins instead of logarithmic bins
equispec n flagequispec flag to use equidistant bins for emission spectra
popmev n flagpopmev flag to use initial population per MeV instead of histograms
segment 1 segment number of segments to divide emission energy grid
maxlevelstar 30 nlevmax maximum number of included discrete levels for target
maxlevelsres 10 nlevmaxres maximum number of included discrete levels for residual nucleus
maxlevelsbin g 10 nlevbin maximum number of included discrete levels for gamma channel
maxlevelsbin n 30 nlevbin maximum number of included discrete levels for neutron channel
maxlevelsbin p 10 nlevbin maximum number of included discrete levels for proton channel
maxlevelsbin d 5 nlevbin maximum number of included discrete levels for deuteron channel
maxlevelsbin t 5 nlevbin maximum number of included discrete levels for triton channel
maxlevelsbin h 5 nlevbin maximum number of included discrete levels for helium-3 channel
maxlevelsbin a 10 nlevbin maximum number of included discrete levels for alpha channel
ltarget 0 ltarget excited level of target
isomer 1.00E+00 isomer definition of isomer in seconds
transpower 15 transpower power for transmission coefficient limit
transeps 1.00E-18 transeps limit for transmission coefficient
xseps 1.00E-17 xseps limit for cross sections
popeps 1.00E-13 popeps limit for population cross section per nucleus
Rfiseps 1.00E-06 Rfiseps ratio for limit for fission cross section per nucleus
elow 1.00E-06 elow minimal incident energy for nuclear model calculations
angles 90 nangle number of angles
anglescont 18 nanglecont number of angles for continuum
anglesrec 1 nanglerec number of recoil angles
maxenrec 20 maxenrec number of recoil energies
channels n flagchannels flag for exclusive channels calculation
maxchannel 4 maxchannel maximal number of outgoing particles in individual channel description
micro n flagmicro flag for completely microscopic Talys calculation
best n flagbest flag to use best set of adjusted parameters
bestbranch y flagbestbr flag to use flag to use only best set of branching ratios
bestend n flagbestend flag to put best set of parameters at end of input file
relativistic y flagrel flag for relativistic kinematics
recoil n flagrecoil flag for calculation of recoils
labddx n flaglabddx flag for calculation of DDX in LAB system
recoilaverage n flagrecoilav flag for average velocity in recoil calculation
channelenergy n flagEchannel flag for channel energy for emission spectrum
reaction y flagreaction flag for calculation of nuclear reactions
ngfit n flagngfit flag for using fitted (n,g) nuclear model parameters
nnfit n flagnnfit flag for using fitted (n,n), (n,2n) and (n,p) nuclear model parameters
nafit n flagnafit flag for using fitted (n,a) nuclear model parameters
astro y flagastro flag for calculation of astrophysics reaction rate
astrogs y flagastrogs flag for calculation of astrophysics reaction rate with target in ground state only
astroex n flagastroex flag for calculation of astrophysics reaction rate to long-lived excited states
nonthermlev -1 nonthermlev excited level non-thermalized in the calculation of astrophysics rate
massmodel 2 massmodel model for theoretical nuclear mass
expmass y flagexpmass flag for using experimental nuclear mass if available
disctable 1 disctable table with discrete levels
production n flagprod flag for isotope production
outfy n flagoutfy flag for output detailed fission yield calculation
gefran 50000 gefran number of random events for GEF calculation
Estop 1000.000 Estop incident energy above which TALYS stops
rpevap n flagrpevap flag for evaporation of residual products at high incident energies
maxZrp 12 maxZrp maximal number of protons away from the initial compound nucleus before residual evaporation
maxNrp 32 maxNrp maximal number of neutrons away from the initial compound nucleus before residual evaporation
#
# Optical model
#
localomp y flaglocalomp flag for local (y) or global (n) optical model
dispersion n flagdisp flag for dispersive optical model
jlmomp n flagjlm flag for using semi-microscopic JLM OMP
riplomp n flagriplomp flag for RIPL OMP
riplrisk n flagriplrisk flag for going outside RIPL mass validity range
optmodall n flagompall flag for new optical model calculation for all residual nuclei
incadjust y flagincadj flag for OMP adjustment on incident channel also
omponly n flagomponly flag to execute ONLY an optical model calculation
autorot n flagautorot flag for automatic rotational coupled channels calculations for A > 150
spherical n flagspher flag to force spherical optical model
soukho y flagsoukho flag for Soukhovitskii OMP for actinides
coulomb y flagcoulomb flag for Coulomb excitation calculation with ECIS
statepot n flagstate flag for optical model potential for each excited state
maxband 0 maxband highest vibrational band added to rotational model
maxrot 2 maxrot number of included excited rotational levels
sysreaction sysreaction particles with reaction cross section from systematics
rotational n rotational particles with possible rotational optical model
core -1 core even-even core for weakcoupling (-1 or 1)
ecissave n flagecissave flag for saving ECIS input and output files
eciscalc y flageciscalc flag for new ECIS calculation for outgoing particles and energy grid
inccalc y flaginccalc flag for new ECIS calculation for incident channel
endfecis y flagendfecis flag for new ECIS calculation for ENDF-6 files
radialmodel 2 radialmodel model for radial matter densities (JLM OMP only)
jlmmode 0 jlmmode option for JLM imaginary potential normalization
alphaomp 6 alphaomp alpha OMP (1=normal, 2= McFadden-Satchler, 3-5= folding potential, 6,8= Avrigeanu, 7=Nolte)
deuteronomp 1 deuteronomp deuteron OMP (1=normal, 2=Daehnick, 3=Bojowald, 4=Han-Shi-Shen, 5=An-Cai)
#
# Compound nucleus
#
widthfluc y flagwidth flag for width fluctuation calculation
widthmode 1 wmode designator for width fluctuation model
compound y flagcomp flag for compound nucleus model
fullhf n flagfullhf flag for full spin dependence of transmission coefficients
eciscompound n flageciscomp flag for compound nucleus calculation by ECIS
cpang n flagcpang flag for compound angular distribution calculation for incident charged particles
urr n flagurr flag for URR calculation
urrnjoy n flagurrnjoy flag for normalization of URR parameters with NJOY method
lurr 2 lurr maximal orbital angular momentum for URR
#
# Gamma emission
#
gammax 2 gammax number of l-values for gamma multipolarity
strength 9 strength model for E1 gamma-ray strength function
strengthM1 3 strengthM1 model for M1 gamma-ray strength function
electronconv y flagelectron flag for application of electron conversion coefficient
racap y flagracap flag for radiative capture model
ldmodelracap 1 ldmodelracap level density model for direct radiative capture
upbend y flagupbend flag for low-energy upbend of photon strength function
psfglobal n flagpsfglobal flag for global photon strength functions only
#
# Pre-equilibrium
#
preequilibrium n flagpreeq flag for pre-equilibrium calculation
preeqmode 2 preeqmode designator for pre-equilibrium model
multipreeq n flagmulpre flag for multiple pre-equilibrium calculation
mpreeqmode 2 mpreeqmode designator for multiple pre-equilibrium model
breakupmodel 1 breakupmodel model for break-up reaction: 1. Kalbach 2. Avrigeanu
phmodel 1 phmodel particle-hole state density model
pairmodel 2 pairmodel designator for pre-equilibrium pairing model
preeqspin 1 pespinmodel model for pre-equilibrium spin distribution
giantresonance y flaggiant flag for collective contribution from giant resonances
preeqsurface y flagsurface flag for surface effects in exciton model
preeqcomplex y flagpecomp flag for Kalbach complex particle emission model
twocomponent y flag2comp flag for two-component pre-equilibrium model
ecisdwba y flagecisdwba flag for new ECIS calculation for DWBA for MSD
onestep n flagonestep flag for continuum one-step direct only
#
# Level densities
#
ldmodel 1 ldmodelall level density model
ldmodelCN 1 ldmodelCN level density model for compound nucleus
shellmodel 1 shellmodel model for shell correction energies
kvibmodel 2 kvibmodel model for vibrational enhancement
spincutmodel 1 spincutmodel model for spin cutoff factor for ground state
asys n flagasys flag for all level density parameters a from systematics
parity n flagparity flag for non-equal parity distribution
colenhance n flagcolall flag for collective enhancement of level density for all nuclides
ctmglobal n flagctmglob flag for global CTM model (no discrete level info)
gshell n flaggshell flag for energy dependence of single particle level density parameter g
colldamp n flagcolldamp flag for damping of collective effects in effective level density
#
# Fission
#
fission n flagfission flag for fission
fismodel 3 fismodel fission model
fismodelalt 3 fismodelalt alternative fission model for default barriers
hbstate y flaghbstate flag for head band states in fission
class2 y flagclass2 flag for class2 states in fission
fispartdamp n flagfispartdamp flag for fission partial damping
massdis n flagmassdis flag for calculation of fission fragment mass yields
ffevaporation y flagffevap flag for calculation of particle evaporation from fission fragment mass yields
fisfeed n flagfisfeed flag for output of fission per excitation bin
fymodel 2 fymodel fission yield model, 1: Brosa 2: GEF
ffmodel 1 ffmodel fission fragment model, 1: GEF 2: HF3D (Okumura) 3: SPY
pfnsmodel 1 pfnsmodel PFNS model, 1: Iwamoto 2: from FF decay
ffspin n flagffspin flag to use spin distribution in initial FF population
#
# Output
#
outmain y flagmain flag for main output
outbasic y flagbasic flag for output of basic information and results
outpopulation y flagpop flag for output of population
outcheck y flagcheck flag for output of numerical checks
outlevels y flaglevels flag for output of discrete level information
outdensity y flagdensity flag for output of level densities
outomp y flagoutomp flag for output of optical model parameters
outkd n flagoutkd flag for output of KD03 OMP parameters
outdirect y flagdirect flag for output of direct reaction results
outinverse y flaginverse flag for output of transmission coefficients and inverse reaction cross sections
outdecay n flagdecay flag for output of decay of each population bin
outtransenergy y flagtransen flag for output of transmission coefficients per energy
outecis n flagoutecis flag for output of ECIS results
outgamma y flaggamma flag for output of gamma-ray information
outpreequilibrium n flagpeout flag for output of pre-equilibrium results
outfission n flagfisout flag for output of fission information
outdiscrete y flagdisc flag for output of discrete state cross sections
outspectra n flagspec flag for output of double-differential cross sections
outbinspectra n flagbinspec flag for output of emission spectrum per excitation bin
resonance n flagres flag for output of low energy resonance cross sections
group n flaggroup flag for output of low energy groupwise cross sections
addiscrete y flagadd flag for addition of discrete states to spectra
addelastic y flagaddel flag for addition of elastic peak to spectra
outangle n flagang flag for output of angular distributions
outlegendre n flaglegendre flag for output of Legendre coefficients
ddxmode 0 ddxmode mode for double-differential cross sections
outdwba n flagoutdwba flag for output of DWBA cross sections for MSD
outgamdis n flaggamdis flag for output of discrete gamma-ray intensities
outexcitation n flagexc flag for output of excitation functions
components n flagcompo flag for output of cross section components
endf n flagendf flag for information for ENDF-6 file
endfdetail y flagendfdet flag for detailed ENDF-6 information per channel
sacs n flagsacs flag for statistical analysis of cross sections
partable n flagpartable flag for output of model parameters on separate file
block n flagblock flag to block spectra, angle and gamma files
########## BASIC REACTION PARAMETERS ##########
Projectile : neutron Mass in a.m.u. : 1.008665
Target : 35Cl Mass in a.m.u. : 34.968853
Included channels:
gamma
neutron
proton
deuteron
triton
helium-3
alpha
73 incident energies (LAB):
2.058E-06
4.115E-06
6.173E-06
8.231E-06
1.029E-05
2.058E-05
3.087E-05
4.115E-05
5.144E-05
6.173E-05
7.202E-05
8.231E-05
9.260E-05
1.029E-04
2.058E-04
3.087E-04
4.115E-04
5.144E-04
6.173E-04
7.202E-04
8.231E-04
9.260E-04
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.010
0.014
0.019
0.023
0.027
0.031
0.035
0.039
0.043
0.047
0.051
0.056
0.060
0.064
... 52292 more lines ...
|
| Attachment 5: Cl36run
|
#!/bin/bash
rm -r Cl36/
mkdir Cl36
for((c=1;c<=1000;c++))
do
./a.out << EOF
35 1.5 # amass spin
0 1000000 #min max
0 23005.84570 1.1310e-4 #ell D S
0.74567 0.149134 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Cl36/run1_0l_$c
cp out_all2 Cl36/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
35 1.5 # amass spin
0 1000000 #min max
1 13704.58008 3.6277e-4 #ell D S
1.06018 0.212036 #<G> SigmaG
0 0 #Gx sigma
0 #bias
EOF
cp out_all1 Cl36/run1_1l_$c
cp out_all2 Cl36/run2_1l_$c
cp out_all3 Cl36/run3_1l_$c
cp out_all4 Cl36/run4_1l_$c
rm out*
done
|
|
4
|
Wed Nov 1 16:14:37 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of Cl35+n |
Cross Section Distributions for Cl35+n; where macs01 is for l=0, macs11 is for l=1and macs1 is for the total. The most probable value is less than what is expected (around 13 in the EXCEL sheet). |
| Attachment 1: macs1.pdf
|
|
| Attachment 2: macs01.pdf
|
|
| Attachment 3: macs11.pdf
|
|
|
3
|
Thu Oct 26 18:09:12 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of Na23+n |
Cross Section Distributions for Na23+n;where macs01 is for l=0, macs11 is for l=1and macs1 is for the total. This is an example of logarithmic decline due to high level spacing (stable nuclei). |
| Attachment 1: macs1.pdf
|
|
| Attachment 2: macs01.pdf
|
|
| Attachment 3: macs11.pdf
|
|
|
2
|
Thu Oct 26 17:59:17 2023 |
Selin Ezgi Birincioglu | Maxwellian Cross Section Distribution of Cl37+n |
Cross Section Distributions for Cl37+n; where macs01 is for l=0, macs11 is for l=1 and macs1 is for the total. This is an example of Gaussian-like behaviour. |
| Attachment 1: macs1.pdf
|
|
| Attachment 2: macs01.pdf
|
|
| Attachment 3: macs11.pdf
|
|
|
1
|
Wed Oct 18 15:54:26 2023 |
Selin Ezgi Birincioglu | Parameters from TALYS |
These are the parameters taken from TALYS to further utilize them as inputs for MyGOE.f - generating nuclear levels. In the sheet the masses are target + n while the input mass is the target nuclei. |
| Attachment 1: TALYSparameters.xlsx
|