Abstract
Here P is the density operator of the system under consideration, and σ ± and σ 3 are the usual Pauli matrices, acting on atom i whose states are |1 > or |0 >, representing, respectively, the atom being in an excited state or in the ground state. B and C are appropriate decay constants and s has been called the pumping parameter [1]. It varies from s = 0 for pure damping to s = 1 for full laser action. To solve the corresponding quantum master equations, three approaches have been taken: First, one focuses on the case of one atom. Second, one truncates eq. (1) and derives semi-classical models. Third, one employs numerical simulation methods such as the quantum trajectory method. While the latter method is very popular, it should be noted that the numerical complexity of the problem increases exponentially with the number of atoms, and so numerical methods soon become unfeasible.