Saturation effects in coherent anti-Stokes Raman scattering (CARS) spectroscopy are discussed. The discussion is limited to Raman-resonant CARS (${\ensuremath{\omega}}_{l}\ensuremath{-}{\ensuremath{\omega}}_{s}\ensuremath{\simeq}{\ensuremath{\omega}}_{21}$, where ${\ensuremath{\omega}}_{l,s}$ are the frequencies of the pump fields with powers
Pl,s, and ${\ensuremath{\omega}}_{21}$ is the frequency of the Raman transition $|1〉\ensuremath{\rightarrow}|2〉$) with the possible addition of a one-photon resonance (${\ensuremath{\omega}}_{l}\ensuremath{\simeq}{\ensuremath{\omega}}_{31}$ is the frequency of the electronic transition $|1〉\ensuremath{\rightarrow}|3〉$). For these cases we show that the CARS polarization
pCARS is proportional to the off-diagonal density-matrix element ${\ensuremath{\rho}}_{21}$. In order to determine ${\ensuremath{\rho}}_{21}$, we use Laplace transforms to solve the Bloch equations for the effective two-level system
∣1〉 and
∣2〉 when ${\ensuremath{\omega}}_{l}$ is far from resonance, or for the three-level system
∣1〉,∣2〉, and
∣3〉, when ${\ensuremath{\omega}}_{l}\ensuremath{\simeq}{\ensuremath{\omega}}_{31}$. The steady-state expression for
pCARS in the former case gives ${\mathcal{p}}_{\mathrm{CARS}}\ensuremath{\propto}{P}_{l}{P}_{s}^{\frac{1}{2}}$ at low powers and ${\mathcal{p}}_{\mathrm{CARS}}\ensuremath{\propto}{P}_{l}^{0}{P}_{s}^{\frac{\ensuremath{-}1}{2}}$ at high powers. In the three-level system, we show that when the pressure is low and at least one field is weak, the slow time dependence of ${\ensuremath{\rho}}_{21}$ must be considered. When one field is strong and the other weak, the CARS spectrum is Stark split. When
Pl is high, for example, ${\mathcal{p}}_{\mathrm{CARS}}\ensuremath{\propto}{P}_{l}^{0}{P}_{s}^{\frac{1}{2}}$ for ${\ensuremath{\omega}}_{s}\ensuremath{\simeq}{\ensuremath{\omega}}_{32}$ and ${\mathcal{p}}_{\mathrm{CARS}}\ensuremath{\propto}{({P}_{l}{P}_{s})}^{\frac{1}{2}}$ when ${\ensuremath{\omega}}_{s}\ensuremath{\simeq}{\ensuremath{\omega}}_{32}\ifmmode\pm\else\textpm\fi{}{V}_{13}$ where
V13 is the one-photon Rabi frequency for the $|1〉\ensuremath{\rightarrow}|3〉$ transition. The Wilcox-Lamb approximation is used to reduce the three-level Bloch equations to rate equations containing one- and two-photon terms. When the fields are so weak that both one- and two-photon terms are small compared to the decay terms, the usual expression for
pCARS is reproduced. If only the direct two-photon processes are important, the effective-two-system results are reproduced. When both fields are intense and nearly resonant, the steady state is rapidly achieved. The results for the case where one field is much stronger than the other are essentially the same as those for one strong and one weak field. When the fields are of comparable strength, ${\mathcal{p}}_{\mathrm{CARS}}\ensuremath{\propto}{P}_{l}^{\frac{1}{2}}{P}_{s}^{0}$ for ${\ensuremath{\omega}}_{l}\ensuremath{\simeq}{\ensuremath{\omega}}_{31}$ and ${\ensuremath{\omega}}_{s}\ensuremath{\simeq}{\ensuremath{\omega}}_{32}$, and the CARS spectrum is split into five components.