Gauge invariance of the natural lineshape and dissipative dynamics of a two-level atom

The calculation of the natural lineshape of an excited two-level atom (TLA) has long been known to be gauge-dependent, with certain experiments in better agreement with the lineshape calculated with the dipole gauge. We show that by using a Coulomb gauge Hamiltonian truncated in a manner consistent with the gauge principle, the correct output spectrum can be obtained. For TLAs undergoing dynamics arising from additional Hamiltonian couplings, we also show that the master equation is gauge-invariant under the same conditions of validity as the Born-Markov approximation, despite different gauges having different spectral densities. These results highlight the importance of using correctly truncated gauge-invariant Hamiltonians in input-output theory for accurate frequency-dependent spectra, even in weak coupling regimes.