bandObject.discretize_FS()
bandObject.dos_k_func()
bandObject.doping()
# # bandObject.figMultipleFS2D()
# # bandObject.figDiscretizeFS2D()
## Conductivity Object
condObject = Conductivity(bandObject,
Bamp=Bmin,
Bphi=0,
Btheta=0,
gamma_0=1,
gamma_k=0,
power=0)
condObject.Ntime = 2000
## Empty arrays
rhoxx_array = np.empty_like(B_array, dtype=np.float64)
rhoxy_array = np.empty_like(B_array, dtype=np.float64)
RH_array = np.empty_like(B_array, dtype=np.float64)
for i, B in enumerate(B_array):
condObject.Bamp = B
condObject.solveMovementFunc()
condObject.chambersFunc(0, 0)
condObject.chambersFunc(0, 1)
sigma_xx = condObject.sigma[0, 0]
sigma_xy = condObject.sigma[0, 1]
hPocket.dos_k_func()
hPocket.doping()
ePocket.discretize_FS(mesh_xy_rough=2001)
ePocket.dos_k_func()
ePocket.doping()
## Conductivity Object ---------------
hPocketCondObject = Conductivity(hPocket,
Bamp=45,
gamma_0=15,
gamma_k=0,
power=12,
gamma_dos_max=0,
factor_arcs=1)
hPocketCondObject.Ntime = 1000 # better for high magnetic field values
## Transport coeffcients -------------
## Empty arrays
rhoxx_array = np.empty_like(B_array, dtype=np.float64)
rhoxy_array = np.empty_like(B_array, dtype=np.float64)
rhozz_array = np.empty_like(B_array, dtype=np.float64)
RH_array = np.empty_like(B_array, dtype=np.float64)
for i, B in enumerate(B_array):
hPocketCondObject.Bamp = B
hPocketCondObject.solveMovementFunc()
hPocketCondObject.chambersFunc(0, 0)
hPocketCondObject.chambersFunc(0, 1)