In this series of simulations we assumed continuous a µw irradiation, which is part of model (B), with γ n,ee = 20.0 Hz. Besides, the radiation damping time constant and the transverse relaxation rate γ 2 were set to the fixed value γ rd = 2.0 kHz and γ 2 = 650 Hz, so that the ratio γ rd /γ 2 was above unity, the theoretical maser threshold. Note that changing the value of m o also affects the value of the radiation damping constant, so keeping γ rd fixed in these simulations implied changing the parameter G in Equation 3 accordingly. The value of γ n,ee = 20.0 Hz was chosen so as to mimic a buildup process that occurs on a much longer time scale than individual maser pulses, and therefore likely not interfering with the latter.