Rigorous solution to a quantum statistical mechanical laser model

A quantum mechanical laser model with relaxation and pumping mechanisms is solved rigorously. A basic equation for the density matrix is derived by the damping theory and is transformed into a corresponding c-number equation for a (quasi-) probability density. This is done with the aid of the quantum phase space method. The probability density is expanded in terms of orthogonal polynomials. The expansion coefficients are solved to give a continued fraction. A complete solution is obtained, namely, time evolution of the probability density is determined as well as that for certain physical quantities. The solution is valid even for strong coupling between photons and atoms: it is free from restriction on system parameters. Detailed studies on dynamics are performed for typical values of the system parameters. This is a prototype of interacting quantum nonequilibrium systems. Relevance to systems other than a laser is briefly mentioned.