def _destroy_excitons_2(alpha, m, beta, n, alpha_p, m_p, beta_p, n_p):
dx = kron_delta(alpha, alpha_p) * kron_delta(beta, 0) * kron_delta(beta_p, 1) \
* kron_delta(m, m_p) * kron_delta(n, n_p)
return dx
def _destroy_excitons(alpha, m, beta, n):
dx = kron_delta(alpha, 0) * kron_delta(beta, 1) * kron_delta(m, n)
return dx
def _destroy_photons_2(alpha, m, beta, n, alpha_p, m_p, beta_p, n_p):
dp = np.sqrt(n_p) * kron_delta(alpha, alpha_p) * kron_delta(beta, beta_p) \
* kron_delta(m, m_p) * kron_delta(n, n_p-1)
return dp
def _destroy_photons(alpha, m, beta, n):
dp = np.sqrt(n) * kron_delta(alpha, beta) * kron_delta(m, n - 1)
return dp