daily pastebin goal
71%
SHARE
TWEET

Untitled

a guest May 16th, 2018 116 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1.  
  2. #1.PHYSICALL CONSTANTS:
  3.  
  4.  
  5. c = 2.99792458e8                  # speed of light    [m/s]
  6. h = 6.62607004e-34                # Planck constant   [J*s]
  7. e_0 = 8.854187817e-12             # vacuum constant   [F/m]
  8. u_B = 927.4009994e-26             # Bohr constant     [J/T]
  9. m_p = 1.66053904e-27              # mass of protone   [kg]
  10. k_B = 1.38064852e-23              # Boltzman constant [J/K]
  11. g = 9.81                          # gravity constant  [m/s^2]
  12.  
  13.  
  14. P = 1                             # power of dipole trap                 [W]
  15. P_lev = 2.4                       # power of levitation potential        [W]
  16. o_0 = 50e-6                       # omega_0 for dipole trap
  17. dB = 31.1e-8                      # magnetic field
  18. C = 7e-6                          # trap potential constants             [um]
  19. A = 50e-6                         # trap potential constants             [um]
  20. omega_z = 23e-6                   # trap potential omega                 [um]
  21. T = 0.00025                       # 1/4kHz in seconds-1 #period          [Hz]
  22. l_lev = 1.064e-6                  # wave length to levitating potential  [um]
  23. l_dip = 1.064e-6                  # wave length to dipole potential      [um]
  24. w_lev = 2*math.pi*c/l_lev
  25. w_dip = 2*math.pi*c/l_dip
  26. x = 0
  27. #------------------------------------------------------------------------------------------------------------------------------------------------------
  28.  
  29. #PARTIClE PARAMETERS (omega is ok)
  30. type_particle = '133CS'
  31.  
  32. if type_particle == '133CS':
  33.     w_x = np.array([2*math.pi*c*894.59296e-9, 2*math.pi*c*852.3472758e-9])     # omega
  34.     g_x = np.array([2*math.pi*4575000, 2*math.pi*5234000])                     # decay rate
  35.     g_j = 2.00254032                                                           # total electronic g-factor
  36.     m_particle = 132.905429*m_p                                                # mass of particle
  37.     aHF=h*2.2981579425e-9                                                      # magnetic dipole
  38.     I=7/2
  39.     j0=1/2
  40.     j1=1/2
  41.     j2=3/2
  42.    
  43.  
  44. elif type_particle == '41K':
  45.     w_x = np.array([2*math.pi*770.107919e-9,2*math.pi*766.7004587e-9])         # omega
  46.     g_x = np.array([2*math.pi*5.956e6,2*math.pi*6.035e6])                      # decay rate
  47.     g_j = 2.00229421                                                           # total electronic g-factor
  48.     m_particle = 40.961825*m_p                                                 # mass of particle
  49.     aHF=h*127.006935e-6                                                        # magnetic dipole
  50.     I=3/2
  51.     j0=1/2
  52.     j1=1/2
  53.     j2=3/2
  54.  
  55. #------------------------------------------------------------------------------------------------------------------------------------------------------
  56.  
  57. #PLOT PARAMETERS
  58.  
  59. #------------------------------------------------------------------------------------------------------------------------------------------------------
  60. start = -100e-6
  61. end = 100e-6
  62. step = 400
  63. #---------------------------------------------------------------------------------------------------------------------------
  64.  
  65. def Re_alpha(w,w_x,g_x): # is okej as mathematical equation
  66.     return 3*math.pi*e_0*c**3*np.sum(1/w_x**3*(g_x/(w_x-w)+g_x/(w_x+w)))
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top