Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #1.PHYSICALL CONSTANTS:
- c = 2.99792458e8 # speed of light [m/s]
- h = 6.62607004e-34 # Planck constant [J*s]
- e_0 = 8.854187817e-12 # vacuum constant [F/m]
- u_B = 927.4009994e-26 # Bohr constant [J/T]
- m_p = 1.66053904e-27 # mass of protone [kg]
- k_B = 1.38064852e-23 # Boltzman constant [J/K]
- g = 9.81 # gravity constant [m/s^2]
- P = 1 # power of dipole trap [W]
- P_lev = 2.4 # power of levitation potential [W]
- o_0 = 50e-6 # omega_0 for dipole trap
- dB = 31.1e-8 # magnetic field
- C = 7e-6 # trap potential constants [um]
- A = 50e-6 # trap potential constants [um]
- omega_z = 23e-6 # trap potential omega [um]
- T = 0.00025 # 1/4kHz in seconds-1 #period [Hz]
- l_lev = 1.064e-6 # wave length to levitating potential [um]
- l_dip = 1.064e-6 # wave length to dipole potential [um]
- w_lev = 2*math.pi*c/l_lev
- w_dip = 2*math.pi*c/l_dip
- x = 0
- #------------------------------------------------------------------------------------------------------------------------------------------------------
- #PARTIClE PARAMETERS (omega is ok)
- type_particle = '133CS'
- if type_particle == '133CS':
- w_x = np.array([2*math.pi*c*894.59296e-9, 2*math.pi*c*852.3472758e-9]) # omega
- g_x = np.array([2*math.pi*4575000, 2*math.pi*5234000]) # decay rate
- g_j = 2.00254032 # total electronic g-factor
- m_particle = 132.905429*m_p # mass of particle
- aHF=h*2.2981579425e-9 # magnetic dipole
- I=7/2
- j0=1/2
- j1=1/2
- j2=3/2
- elif type_particle == '41K':
- w_x = np.array([2*math.pi*770.107919e-9,2*math.pi*766.7004587e-9]) # omega
- g_x = np.array([2*math.pi*5.956e6,2*math.pi*6.035e6]) # decay rate
- g_j = 2.00229421 # total electronic g-factor
- m_particle = 40.961825*m_p # mass of particle
- aHF=h*127.006935e-6 # magnetic dipole
- I=3/2
- j0=1/2
- j1=1/2
- j2=3/2
- #------------------------------------------------------------------------------------------------------------------------------------------------------
- #PLOT PARAMETERS
- #------------------------------------------------------------------------------------------------------------------------------------------------------
- start = -100e-6
- end = 100e-6
- step = 400
- #---------------------------------------------------------------------------------------------------------------------------
- def Re_alpha(w,w_x,g_x): # is okej as mathematical equation
- return 3*math.pi*e_0*c**3*np.sum(1/w_x**3*(g_x/(w_x-w)+g_x/(w_x+w)))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement