Solomon equations for qubit and two-level systems

We model and measure the combined relaxation of a qubit, a.k.a. central spin, coupled to a discrete two-level system (TLS) environment. We present a derivation of the Solomon equations starting from a general Lindblad equation for the qubit and an arbitrary number of TLSs. If the TLSs are much longer lived than the qubit, the relaxation becomes non-exponential. In the limit of large numbers of TLSs the populations are likely to follow a power law, which we illustrate by measuring the relaxation of a superconducting fluxonium qubit. Moreover, we show that the Solomon equations predict non-Poissonian quantum jump statistics, which we confirm experimentally.