We prove several results about the rate of convergence to stationarity, that is, the spectral gap, for the M/M/n queue in the Halfin-Whitt regime. We identify the limiting rate of convergence to steady-state, and discover an asymptotic phase transition that occurs w.r.t. this rate. In particular, we demonstrate the existence of a constant B * ≈ 1.85772 s.t. when a certain excess parameter B ∈ (0, B * ], the error in the steady-state approximation converges exponentially fast to zero at rate B 2 4 . For B > B * , the error in the steady-state approximation converges exponentially fast to zero at a different rate, which is the solution to an explicit equation given in terms of special functions. This result may be interpreted as an asymptotic version of a phase transition proven to occur for any fixed n by van Doorn [Stochastic Monotonicity and Queueing Applications of Birth-death Processes (1981) Springer].We also prove explicit bounds on the distance to stationarity for the M/M/n queue in the Halfin-Whitt regime, when B < B * . Our bounds scale independently of n in the Halfin-Whitt regime, and do not follow from the weak-convergence theory.and Whitt [18], who studied the GI/M/n system for large n when the traffic intensity scales like 1− Bn −1/2 for some strictly positive excess parameter B. They proved that, under minor technical assumptions on the inter-arrival distribution, this sequence of GI/M/n queueing models has the following properties:(i) the steady-state probability that an arriving job has to wait for service has a nontrivial limit;(ii) the sequence of queueing processes, normalized by n 1/2 , converges weakly to a nontrivial positive recurrent diffusion, a.k.a. the HW diffusion;(iii) the sequence of steady-state queue length distributions, normalized by n 1/2 , is tight and converges distributionally to the mixture of a point mass at zero and an exponential distribution.Since the steady-state behavior of the M/M/n queue in the HW regime is quite simple [18], while the transient dynamics are more complicated [18], it is common to use the steady-state approximation to the transient distribution [16]. Thus it is important to understand the quality of the steady-state approximation. The only work along these lines seems to be the recent papers [38,39], in which the authors study the Laplace transform of the HW and related diffusions, and prove several results analogous to our own for these diffusions. The key difference is that in this paper we study the prelimit diffusion-scaled M/M/n queue, not the limiting diffusion. We note that the relevant transform functions were also studied in [1], although in a different context. Also, similar questions were studied for the associated sequence of fluid-scaled queues in [23].The question of how quickly the positive recurrent M/M/n queue approaches stationarity has a rich history in the queueing literature. In [28], Morse derives an explicit solution for the transient M/M/1 queue, and discusses implications for the exponential rate of convergence to stationarity...