\begin{gather*}
  \begin{pmatrix}
    % row 1
    -\frac{1-h}{\tau_\mathrm{w2}} - \frac{1}{\tau_\mathrm{wR}} & \frac{\NW{}}{\tau_\mathrm{w2}} + \frac{\NQ}{\tau_\mathrm{2w}} & 0\\
    % row 2
    \frac{1-h}{\NQ\tau_\mathrm{w2}} & -\frac{\NW{}}{\NQ\tau_\mathrm{w2}} - \frac{1}{\tau_\mathrm{2w}} - \frac{1-f}{\tau_{21}} - \frac{f}{\tau_{12}}
    & \frac{h}{\tau_{21}} + \frac{1-h}{\tau_{12}} \\
    % row 3
    0 & \frac{1-f}{\tau_{21}} + \frac{f}{\tau_{12}} &
    -\frac{h}{\tau_{21}} - \frac{1-h}{\tau_{12}} -\frac{2f}{\tau_\mathrm{1R}} - \frac{2g_\mathrm{max}\photon}{\NQ\sigma\hbar\omega}
  \end{pmatrix}
\end{gather*}