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: 46     Entry time: Wed Jun 12 01:51:51 2024
Author: RuiJiu, David 
Category: Analysis 
Subject: ID 
The preliminary identification has been thoroughly checked using the RionID software. We have confirmed that the identification is accurate.

Attached can be found the simulation file with each of the expected ions, their revolutions frequencies in the 125th harmonic and expected yield from the LISE file attached.
Also attached are some figures from RionID that highlight the region of interest for the 124th and 125th harmonics. The center frequency of the 245 MHz resonator is closer to the 124th harmonic for 229Th+89 which is located at 244.769 MHz. 

If necessary, we can slightly shift the resonator curve to higher frequencies to fine-tune the signal-to-noise ratio for 229Th+89.
Attachment 1: E142_TEline-ESR_229Th.lpp  1000 kB  Uploaded Wed Jun 12 02:59:34 2024  | Hide | Hide all | Show all
Version 17.5.18

{============================= Main Part ======================================}
[general]
                    File = //wsl.localhost/Ubuntu-18.04/home/chenruijiu-ubuntu/experiment/ESR/2022_229Th/E142_TEline-ESR_229Th.lpp
                    Date = 12-06-2024
                    Time = 01:29:23
           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 = 555.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
... 3039 more lines ...
Attachment 2: simulation_result.out  141 kB  Uploaded Wed Jun 12 02:59:44 2024  | Show | Hide all | Show all
Attachment 3: 229Th89_124h.png  86 kB  Uploaded Wed Jun 12 03:02:14 2024  | Show | Hide all | Show all
Attachment 4: 229Th89_125h.png  87 kB  Uploaded Wed Jun 12 03:02:20 2024  | Hide | Hide all | Show all
229Th89_125h.png
ELOG V3.1.5-fc6679b