• API
• FAQ
• Tools
• Archive
daily pastebin goal
29%
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.

Top