Not yet reviewed by Pith; the record is open.

Many applications of nitrogen-vacancy (NV) centers in diamond crucially rely on a spectrally narrow and stable optical zero-phonon line transition. Though many impressive proof-of-principle experiments have been demonstrated, much work remains in engineering NV centers with spectral properties that are sufficiently robust for practical implementation. To elucidate the mechanisms underlying their interactions with the environment, we apply multi-dimensional coherent spectroscopy to an NV center ensemble in bulk diamond at cryogenic temperatures. Our spectra reveal thermal dephasing due to quasi-localized vibrational modes as well as ultrafast spectral diffusion on the picosecond timescale. The intrinsic, ensemble-averaged homogeneous linewidth is found to be in the tens of GHz range by extrapolating to zero temperature. We also observe a temperature-dependent Stark splitting of the excited state manifold, relevant to NV sensing protocols.