A systematic, rigorous, and complete investigation of the Bloch equations in time-harmonic driving classical field is performed. Our treatment is unique in that it takes full advantage of the partial fraction decomposition over real number field, which makes it possible to find and classify all analytic solutions. Torrey's analytic solution in the form of exponentially damped harmonic oscillations [Phys. Rev. {\bf 76}, 1059 (1949)] is found to dominate the parameter space, which justifies its use at numerous occasions in magnetic resonance and in quantum optics of atoms, molecules, and quantum dots. The unorthodox solutions of the Bloch equations, which do not have the form of exponentially damped harmonic oscillations, are confined to rather small detunings
δ2≲(γ−γt)2/27 and small field strengths
Ω2≲8(γ−γt)2/27, where
γ and
γt describe decay rates of the excited state (the total population relaxation rate) and of the coherence, respectively. The unorthodox solutions being readily accessible experimentally are characterized by rather featureless time dependence.
The laser gain equation is derived in two ways. First, the optical Bloch equations are derived, based upon a semiclassical treatment of a two-level atomic system and a radiation field. As a visualization aid, there follows a discussion of the formally identical classical Bloch equations, which describe the behavior of a magnetic moment in a combined static and rotating magnetic field. It is shown that laser gain is proportional to the imaginary part of the linear susceptibility, which is obtained from the semiclassical treatment. The second viewpoint uses a slightly modified form of Einstein’s A and B coefficients adapted to describing interactions between atoms and monochromatic fields. Various damping and broadening mechanisms are introduced in a reasonable but non-rigorous fashion, and the crucial distinction between homogeneous and inhomogeneous broadening is introduced.