HOME ESR EXPERIMENTS LABS etc
E121 E127 E132 E125 E143 laser_cooling_2021 E142 Ê128 E146 E0052 E0018 E0028 E0038
  Nuclear hyperfine mixing in 229Th89+ 2024  Not logged in ELOG logo
Message ID: 40     Entry time: Mon Jun 10 11:59:12 2024
Author: Ruijiu Chen 
Category: Analysis 
Subject: The ion identification of spectrum recorded on 2022 
The simulated frequency is given in simulation_result.out. The yield is calculated by LISE file. The LISE file is attached here also.
Attachment 1: simulation_result.out  189 kB  Uploaded Mon Jun 10 13:19:25 2024  | Show | Hide all | Show all
Attachment 2: E142_TEline-ESR_229Th.lpp  1.297 MB  Uploaded Mon Jun 10 13:23:19 2024  | Hide | Hide all | Show all
Version 17.5.18

{============================= Main Part ======================================}
[general]
                    File = E:/work/experiment/ESR/2022_E142_229Th/E143_TEline-ESR_229Th.lpp
                    Date = 08-06-2024
                    Time = 23:26:50
           Configuration = GSI/FRS-TA-ESR_2014.lcn
             Optionsfile = GSI_FRS_2021.lopt
                   Title = GSI FRS TA-ESR (2014)
          BlockStructure = DSDSMDASDDD
          NumberOfBlocks = 11

[settings]
                   A,Z,Q = 238U 92+                              ; Mass ElementName Charge+ Beam
                  Energy = 453.9                MeV/u            
               Intensity = 1e+9                 pps              ; enA,pna,pps,kW
            RF frequency = 20                   MHz              
            Bunch length = 1                    ns               
         Settings on A,Z = 229Th                                 ; Mass ElementName Charge+ Beam

[OpticsBeam]
                      BX = 0.472                (+/-)mm          ; one-half the horisontal beam extent (x)
                      BT = 1.59                 (+/-)mrad        ; one-half the horisontal beam divergence(x')
                      BY = 0.245                (+/-)mm          ; one-half the vertical beam extent (y)
                      BF = 2.04                 (+/-)mrad        ; one-half the vertical beam divergence (y')
                      BL = 0                    (+/-)mm          
                      BD = 0.03                 (+/-)%           ; one-half of the momentum spread (dp/p)
                  ShiftX = 0                    mm               ; beam respect to the spectrometer axis
                  AngleX = 0                    mrad             ; beam respect to the spectrometer axis
                  ShiftY = 0                    mm               ; beam respect to the spectrometer axis
                  AngleY = 0                    mrad             ; beam respect to the spectrometer axis
            Scheme Angle = 0                    degrees          
                  ShapeX = 1
                  ShapeT = 1
                  ShapeY = 1
                  ShapeF = 1
                  ShapeL = 1
                  ShapeD = 1
               OptBeam_X = 1                    (+/-)mm          
               OptBeam_T = 30                   (+/-)mrad        
               OptBeam_Y = 1                    (+/-)mm          
               OptBeam_F = 30                   (+/-)mrad        
               OptBeam_L = 0                    (+/-)mm          
               OptBeam_D = 1.5                  (+/-)%           

[options]
               NP simple = 64                                    ; Number of points in distribution
        NP charge states = 16                                    ; Number of points in distribution
                NP wedge = 128                                   ; Number of points in distribution
           Charge states = Yes                                   ; No  & Yes
           CutEdgeEffect = 1                                     ; 1-Yes. Default, 0-no - for extended configurations
       Prim.beam scatter = 0                                     ; 0-without, 1-with
              Delta peak = 0                                     ; 0-without, 1-with
             BrhoMeanMax = 1                                     ; 0-Max, 1-Mean
            BrhoMeMaLeRi = 3                                     ; 0-Max, 1-Mean, 2-Left, 3-Right /for fission/
              CentralCut = 1                                     ; 0-no, 1-moderate, 2-strong

[target]
         Target contents = 0,4,1,9.012                           ; Nomer,Z,Atoms,Mass
        Target thickness = 1,1850,1.85,0,0,0                     ; State,Thickness,density,angle,SolidGas,..
  Target fusion compound = 0
   Targ use for Q-states = 1
           Target Defect = 0,0.1                                 ; [0] choice - % or micron at 0 degree,  [1]=value; 
       Degrader contents = 0,6,1,12.011                          ; Nomer,Z,Atoms,Mass
      Degrader thickness = 0,0,2.26,0,0,0                        ; State,Thickness,density,angle,SolidGas,..
  Degra use for Q-states = 1
         Degrader Defect = 0,0.1                                 ; [0] choice - % or micron at 0 degree,  [1]=value; 

[mechanism]
                Reaction = 0                                     ; 0 - fragm, 1 - fusion-resid, 2 - fusion-fission
               CalcOther = 1                                     ; calculate other reactions
           V calculation = 5                                     ; 0 - constant, 1 - Borrel, 2 - Rami, 3-convolution, 4-two body reaction
             V_opt/Vbeam = 1                                     ; default 1
         Velocity_exceed = 1                                     ; 0 - without, 1-with - two-body recations velocity corrections
Binding Energy for Vf/Vp = 8                    MeV              ; Binding energy for Borrel's expression
    Shift for Vf/Vp calc = 0     
        Prefragment_Rami = 1
                  Sigma0 = 90                   MeV/c            ; default 90
                  SigmaD = 0                    MeV/c            ; default 200
                  SigmaM = 87                   MeV/c            ; default 87
               Asymmetry = 0                    %                ; default 0
          Method v-sigma = 0                                     ; 0 - Goldhaber, 1-Morrissey,2-Friedman,3-Convolution
               G_Surface = 0.95                 MeV/fm^2         
 Symmetry around half_Ab = 1
 Pfaff pickup correction = 0
     ChargeExchangePfaff = 0                                     ; 1 - exclude, 0-forget
            Sigma corr 0 = 0                                     ; Coulomb energy
            Sigma corr 1 = 0                                     ; Projectile mass
           Friedman mode = 2                                     ; 0-Qgg, 1-Surface, 2-Qgg+Surface
         Prefragment_Fri = 1
        Coulomb_Friedman = 1
                  K_Morr = 8                    MeV/A            ; E/A=8MeV/A default; D.Morrissey coef.
               K_MorHalf = 8                    MeV/A            ; E/A=8MeV/A default at Afrag=Aproj/2; D.Morrissey coef.
                 AA_fast = 0
             AA_Modefast = 3
               MethodApf = 2                                     ; Prefragment search method: 0-N/Z direction, 1 -  W*EPAX, 2 -  W*CSgeom,
           ExcitationApf = 1                                     ; Excitation energy for APF search: 0 - dSurface, 1 - dA*Ex
            BarrierShape = 1                                     ; 0-classical, 1-quantum mech.
                 H_omega = 3                    MeV              ; default 3
           Probabilty_CN = 1                                     ; 0/1 use Prbabilty for CN formation
            UseVanishing = 1
              VanishMode = 0                                     ; 0-Sierk, 1-Cohen
           NuclPotential = 1                                     ; 0-Bass, 1-WS
                   WS_V0 = 105                  MeV              
                   WS_R0 = 1.12                 fm               
                    WS_a = 0.75                 fm               
          FusDiffuseness = 1    
              Width Coef = 1                                     ; default 1; for Leon's charge state distribution
          gZt Correction = 1                                     ; default 1;      Leon's C.S.D.
           PowerCoefLeon = 0.477                                 ; default 0.477; Leon's C.S.D.
           Cross section = File                                  ; Fit  & File
           Charge method = 3                                     ; charge calculations method  0-5
      EPAX Cross Section = 4                                     ; cross section calculations method 0-4
        SR Cross Section = 1                                     ; EPAX for SR  0-2
             Energy Loss = 4                                     ; energy loss calculation method 0-3
     Anglular straggling = 1
        AngStragInOptics = 0                                     ; 0-LISE, 1-ATIMA
         StragglingCoef1 = 0.217
         StragglingCoef2 = 1.12 
       Energy straggling = 1                                     ; 0-LISE, 1-ATIMA
       EnergyStragMethod = 1                                     ; 0-integrate, 1-table
        EnergyStragShape = 0                                     ; 0-Gauss, 1-Landau-Vavilov
          EquilThickness = 1                                     ; 0-Charge, 1-Global
              MassMethod = 0                                     ; 0-DB+calcul, 1 + just calcul
            MassDataBase = 0                                     ; 0-A&W, 1-User ME
            Mass formula = 2                                     ; 0-LDM, 1-Myerer, 2: 1+corrections
          MassExcessFile = AME2016.lme
      UseChargeForEnergy = 2                                     ; 0-No, 1-Yes, 2-Auto
         EnergyValueAuto = 30                                    ; default value 30 MeV/u
         EquilibriumMode = 0                                     ; 0-Equil, 1-NonEquil
               UB_Global = 70                                    ; default 70 MeV/u
             MinZ_Global = 29                                    ; default Z>=29
        ChargeStateOptim = 1                                     ; 0-No, 1-Yes
         ZmQ_AfterReactn = 0                                     ; default 0 (full stripped)
             EPAX_p_Norm = 1
               EPAX_p_Un = 1.65
              EPAX_p_Up0 = 1.788
              EPAX_p_Up1 = 0.004721
              EPAX_p_Up2 = -1.303e-5
                EPAX_p_H = 1

[fission]
         FisAngDistShape = 0                                     ; 0-isotropic; 1-anisotropic
     FisMomCutForAngDist = 2                                     ; 0-dont use;  1-use just MatrixKinematics; 2-use for all; (default 2)
      OddEvenCorrections = 1                                     ; 0-dont use;  1-use
 PostScissionEvaporation = 1                                     ; 0-dont use;  1-use
   DeexcitFunctionPoints = 0                                     ; 0- average deexcitation energy;  1- 3 points; 2 - manually
           FisEXmanually = 20                                    ; Excitation energy manually
           FisCSmanually = 1000                                  ; Cross section manually
            FisTXEmethod = 0                                     ; 0-from Edissipated, 1 from Q-value
                   Fis_f = 0.0035                                ; default 0.0045
              FisEXsigma = 5.5                  MeV              ; default 5.5
            FisCS_Global = 1e-12
               FisCS_TKE = 1e-8
                FisBeta1 = 0.625                                 ; deformation of light fragment
                FisBeta2 = 0.625                                 ; deformation of heavy fragment
                FisTKE_d = 2                    fm               ; d-param in Wilkins formula
              FisBetaFit = 1                                     ; 0-manual, 1-fit
                      N0 = 83                                    ; default 82
                     dU0 = -2.65                                 ; default -2.5
                      C0 = 0.7                                   ; default 1.4
                   cpol0 = 0.65                                  ; default 0.65
                  width0 = 0.63                                  ; default 0.63
                      N1 = 90                                    ; default 90
                     dU1 = -3.8                                  ; default -5.5
                      C1 = 0.15                                  ; default 0.16
                   cpol1 = 0.55                                  ; default 0.55
                  width1 = 0.97                                  ; default 0.97

[charge_suppression]
                 FragInd = 1e-3
               FragTotal = 1e-5
                 BeamInd = 1e-12
               BeamTotal = 1e-15

[convolution]
        Convolution mode = 1                                     ; 0-Qgg, 1-Surface, 2-E* per Abraded nucleon
                 Sigma_0 = 91.5   
                 Sigma_1 = 91.5   
                 Sigma_2 = 160     
              CoefConv_0 = 3.344 
              CoefConv_1 = 3     
              CoefConv_2 = 1     
             ShiftConv_0 = 0.1581
             ShiftConv_1 = 0.1487
             ShiftConv_2 = -1     

[evaporation]
          NP evaporation = 32                                    ; Number of points in distribution
               NPevapFis = 8                                     ; Number of points in distribution
              EvapMethod = 2
        StateDensityMode = 2                                     ; 0, 1+pairing, 2+shell
      EvapUnstableNuclei = 1                                     ; 0 - only stable,1 +unstable
              Tunnelling = 1
          AvoidUnboundCS = 1
       ProtectedChannels = 1
           R_Evaporation = 5.7                  fm               ; correction for the effective Coulomb barrier
         Mode_Apf_manual = 0                                     ; 1-manual, 0-auto
             Energy_in_T = 2                                     ; default 2
     EvaporationVelocity = 0                                     ; 0 - quality, 1 -fast
        DeltaOddEvenEvap = 12    
     DeltaOddEvenFission = 14    
   BreakupTemperature250 = 4.7   
   BreakupTemperature150 = 5.9   
   BreakupTemperature050 = 8     
      BreakupDiffuseness = 0.05  
      DissipationKramers = 0
 DissipationStepFunction = 1
         DissipationBeta = 1                                     ; default 2.0
                 mode_1n = 1
                 mode_2n = 0
                 mode_1p = 1
                 mode_2p = 0
                  mode_a = 1
                  mode_d = 0
                  mode_t = 0
                mode_3he = 0
                mode_fis = 1
             mode_brk_up = 1
              mode_gamma = 0

[fission_barrier]
    FissionBarrierFactor = 1    
      FissionBarrierMode = 1                                     ; #0-4
      OddEvenCorrections = 1
        ShellCorrections = 1
             FB_InOutMax = 2                                     ; #0-2 - in/out/max
             ModeForUser = 1                                     ; #0-2
           NdeltaOddEven = 2.5  
           ZdeltaOddEven = 9    

[excitation_energy]
           AbrasionModel = 0                                     ; 0-Geometrical, 1-Exponential
       GeomAA_Correction = 1                                     ; 0 - don't use,1 - use -default
          Thermalization = 0
          ThermaTimeCoef = 3e+0                                  ; 2.1e-22 MeV *s/e(t)
            AbraExpSlope = 0.363 
                Friction = 0                                     ; 0 - off,1 - on
          Ev_A_SigmaCoef = 9.6  
         G_FrictionCoef1 = 6.5  
         G_FrictionCoef2 = 0.5  
           G_FactorCoef1 = 1.5  
           G_FactorCoef2 = 2.5  
               DepthHole = 40     
          EnergyCoef_CB0 = 0       
          EnergyCoef_CB1 = 27       
          EnergyCoef_CB2 = 0       
           SigmaCoef_CB0 = 0       
           SigmaCoef_CB1 = 18       
           SigmaCoef_CB2 = 0       
              D_MeanTemp = 13       
             DN_MeanTemp = 0       
             DZ_MeanTemp = 0       
               LN_median = 20       
             LN_variance = 20       
              LN2_median = 0       
            LN2_variance = 0       
               AA_factor = 1       
          ApplyLimitTemp = 0

[evapauto]
                  tun_a0 = -0.61392
                  tun_a1 = 0.44559
                  tun_a2 = 0.12008
                 A_Bound = 300
               A_Pairing =   1

[plot]
            Start target = Detector                              ; Detector & RF      
            Start of TOF = T1
            Stop  of TOF = M1
         dE-detector-1st = M1
         dE-detector-2nd = M1
            TKE-detector = M1
              X-detector = M1
              Y-detector = M1
                 Tilting = M1
                 Stopper = M2
                RO_Wedge = W1
          ConditionBlock = A0
          Plot threshold = 1e-10                pps              ; minimal value for plot scale
     Shift of TOF for RF = 0                    ns               ; for dE-TOF plot with RF
  Fraction of RF trigger = 1
            UseCondition = 0
         TKE_calibration = 1,1,0,MeV

[cs_file]
              UserDiffCS = 0                                     ; Number of User Diff CS saved in this file
         AppendOverwrite = 1
          AttachedInside = 1
            ShowCSinPlot = 1
                    Chi2 = 1

[sec_reactions]
        NP sec.reactions = 32                                    ; Number of points in distribution
     Secondary reactions = 0
         fiss_FilterUse0 = 1
         fiss_FilterUse1 = 1
         fiss_FilterUse2 = 1
... 3843 more lines ...
Attachment 3: The_ion_identification_of_spectrum_recorded_on_2022_v2.pdf  1021 kB  Uploaded Mon Jun 10 13:30:51 2024  | Hide | Hide all | Show all
The_ion_identification_of_spectrum_recorded_on_2022_v2.pdf The_ion_identification_of_spectrum_recorded_on_2022_v2.pdf The_ion_identification_of_spectrum_recorded_on_2022_v2.pdf The_ion_identification_of_spectrum_recorded_on_2022_v2.pdf The_ion_identification_of_spectrum_recorded_on_2022_v2.pdf The_ion_identification_of_spectrum_recorded_on_2022_v2.pdf
ELOG V3.1.5-fc6679b