def setup_dqd_system(bias, T_c, Gamma_L, Gamma_R, reduce_dim):
'''Instantiates LindbladSystem for double quantum dot model'''
H = np.array([[0, 0, 0], [0, bias / 2., T_c], [0, T_c, -bias / 2.]])
ops = np.array([
np.array([[0, 0, 0], [1., 0, 0], [0, 0, 0]]),
np.array([[0, 0, 1.], [0, 0, 0], [0, 0, 0]])
])
rates = np.array([Gamma_L, Gamma_R])
return LindbladSystem(H, ops, rates, reduce_dim=reduce_dim)
def setup_srl_system(Gamma_L, Gamma_R, reduce_dim):
'''Instantiates LindbladSystem for single resonant level model'''
H = np.array([[0, 0], [0, 0]])
ops = np.array([np.array([[0, 0], [1, 0]]), np.array([[0, 1], [0, 0]])])
rates = np.array([Gamma_L, Gamma_R])
return LindbladSystem(H, ops, rates, reduce_dim=reduce_dim)