We consider a two-level systems coupled to an environment consisting of two distinct energy bands. The system is described by the Hamiltonian
![{\displaystyle H=H_{S}+H_{E}+H_{I}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7cbf3223d90e9fe44eead62200c4a89031172ad3)
The system Hamiltonian
is an ordinary two-level system with an energy difference of
, according to
![{\displaystyle H_{S}=\Delta E|e\rangle \langle e|=\Delta E\sigma _{+}\sigma _{-},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56606b60ebbb6346b9e91b177e6c133b62a15c9a)
where we have set the energy of the ground state
to zero. For the energy bands of the environment, we assume a linear spectrum, i.e,
![{\displaystyle H_{E}=\sum \limits _{n_{1}=1}^{N_{1}}{\frac {\delta \varepsilon }{N_{1}}}n_{1}|n_{1}\rangle \langle n_{1}|+\sum \limits _{n_{2}=1}^{N_{2}}\left(\Delta E+{\frac {\delta \varepsilon }{N_{2}}}n_{2}\right)|n_{2}\rangle \langle n_{2}|,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ca4e894c9e96b8f77cab008d24ea9b1dd31cd929)
where
and
are the total number of states in each band. The interaction Hamiltonian is given by
![{\displaystyle H_{I}=\sum \limits _{n_{1}=1}^{N_{1}}\sum \limits _{n_{2}=1}^{N_{2}}\lambda _{n_{1},n_{2}}{\big (}\sigma _{+}\otimes 1_{E}{\big )}(1_{S}\otimes |n_{1}\rangle \langle n_{2}|)+{\text{H.c.}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/781777ff5549a5b28010e7fba6eca96accdc5150)
i.e., whenever the two-level system is de-excited, an excitation is created in the environment, and vice versa. The matrix elements
are random Gaussianly distributed variables with a variance
. In the following, we assume a separation of energy scales according to
.
This model has some very interesting properties. First of all, the size of the environment is clearly finite, meaning that excitations in the environment may not decay fast enough and the Markov approximation may not be applicable. This is even further strengthened when taking into account that the interaction will lead to very strong correlations between system and environment. At the same time, the weak interaction naturally gives rise to a peturbative expansion.
We proceed by defining a correlated projection operator for the relevant part of the dynamics [1], according to
![{\displaystyle {\mathcal {P}}\rho ={\text{Tr}}\left\{\Pi _{1}\rho \right\}{\frac {1}{N_{1}}}\Pi _{1}+{\text{Tr}}\left\{\Pi _{2}\rho \right\}{\frac {1}{N_{2}}}\Pi _{2}=P_{e}{\frac {1}{N_{1}}}\Pi _{1}+P_{g}{\frac {1}{N_{2}}}\Pi _{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b1aff786394e6461316204841aabbe2451c4ea9)
where we have introducted the projectors
![{\displaystyle {\begin{aligned}\Pi _{1}&=\sum \limits _{n_{1}=1}^{N_{1}}1_{S}\otimes |n_{1}\rangle \langle n_{1}|\\\Pi _{2}&=\sum \limits _{n_{2}=1}^{N_{2}}1_{S}\otimes |n_{2}\rangle \langle n_{2}|\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b7f6e9885ad41e95c839ddd9a0bf6ca3a5037e5d)
and
and
denote the probability to find the two-level system in its excited or ground state, respectively. Because of conservation of probability, we have
, meaning that the relevant part of the dynamics consists of a single variable!
In the following, we will derive its equation of motion. The second order TCL master equation is given by
![{\displaystyle {\mathcal {P}}{\frac {d}{dt}}\rho =\int \limits _{0}^{t}ds{\mathcal {P}}{\mathcal {L}}(t){\mathcal {L}}(s){\mathcal {P}}\rho ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8d7b964a1f535c71d6570cf9b4f98fd4743e92dd)
which in our case can be written as
![{\displaystyle {\dot {P}}_{e}{\frac {1}{N_{1}}}\Pi _{1}+(1-{\dot {P}}_{e}){\frac {1}{N_{2}}}\Pi _{2}=\int \limits _{0}^{t}ds{\text{Tr}}\left\{\Pi _{1}{\mathcal {L}}(t){\mathcal {L}}(s){\mathcal {P}}\rho \right\}{\frac {1}{N_{1}}}\Pi _{1}+{\text{Tr}}\left\{\Pi _{2}{\mathcal {L}}(t){\mathcal {L}}(s){\mathcal {P}}\rho \right\}{\frac {1}{N_{2}}}\Pi _{2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c504aee081e8f7679c60ec67398f4d60f0c5ffed)
since the projection and the time derivative commute. Multiplying with
and taking the trace yields
![{\displaystyle {\begin{aligned}{\dot {P}}_{e}&=\int \limits _{0}^{t}ds{\text{Tr}}\left\{\Pi _{1}{\mathcal {L}}(t){\mathcal {L}}(s){\mathcal {P}}\rho \right\}\\&=-\int \limits _{0}^{t}ds{\text{Tr}}\left\{\Pi _{1}\left[H_{I}(t),\left[H_{I}(s),P_{e}{\frac {1}{N_{1}}}\Pi _{1}+(1-P_{e}){\frac {1}{N_{2}}}\Pi _{2}\right]\right]\right\}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aada75c71d80bcc1e76ac6facc44a7281400b8bb)
Expanding the commutator and making use of the relation
results in
![{\displaystyle {\dot {P}}_{e}=-\int \limits _{0}^{t}ds{\text{Tr}}\left\{P_{e}{\frac {1}{N_{1}}}\Pi _{1}H_{I}(s)H_{I}(t)-(1-P_{e}){\frac {1}{N_{2}}}\Pi _{1}H_{I}(t)H_{I}(s)+{\text{H.c.}}\right\}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f46ae35c7fcd116e2ae6fff07f8060c8d213a9c1)
We can now introduce two relaxation rates
describing the transition rates between the two bands. Specifically, we have
![{\displaystyle \gamma _{i}={\frac {1}{N_{i}}}\int \limits _{0}^{t}ds\sum \limits _{n1}\langle n_{1}|H_{I}(s)H_{I}(t)+{\text{H.c}}|n_{1}\rangle .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3195ad56d3ee3a6819cc554e96bc5460da0c4180)
As the interaction Hamiltonian couples does not contain terms that leave the band index unchanged, we may write the matrix elements as
![{\displaystyle \langle n_{1}|H_{I}(s)H_{I}(t)+{\text{H.c}}|n_{1}\rangle =\sum \limits _{n_{2}}\langle n_{1}|H_{I}(s)|n_{2}\rangle \langle n_{2}|H_{I}(t)|n_{1}\rangle +{\text{c.c.}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/731d7463ede8d5e97aac665e1f791bf6117c4bdd)
As we work in the interaction picture, we have
![{\displaystyle \langle n_{1}|H_{I}(s)|n_{2}\rangle \langle n_{2}|H_{I}(t)|n_{1}\rangle =\exp(i\omega _{12}s)\exp(-i\omega _{12}t)|\lambda _{n1,n2}|^{2},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dcf766c70bec584ff020cf8423318503c1ba9a45)
where we have introduced the transitions frequencies
.
Using the statistical properties of the interaction Hamiltonian and performing the integration over
finally yields [2]
![{\displaystyle \gamma _{i}={\frac {2\lambda ^{2}}{N_{i}}}\sum \limits _{n1,n2}{\frac {\sin(\omega _{12}t)}{\omega _{12}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f5c44a0468d1c87b9200f0105939146fcdb9d3ec)
For the next step, we note that the sinc function appearing in the relaxation rates is a representation of the Dirac
distribution, i.e.,
![{\displaystyle \pi \delta _{t}(\omega _{12})=\lim _{t\to \infty }{\frac {\sin(\omega _{12}t)}{\omega _{12}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b39d1aa259cae4f413c7f2b1b839934bb2a09047)
Then, at large enough times and small enough couplings [3], we may approximate the rates by
![{\displaystyle {\begin{aligned}\gamma _{1}&\approx {\frac {2\pi \lambda ^{2}}{N_{1}}}\sum _{k,l}\delta (E_{n_{1}}-E_{n_{2}})=2\pi \lambda ^{2}{\frac {N_{2}}{\delta \varepsilon }}\\\gamma _{2}&\approx 2\pi \lambda ^{2}{\frac {N_{1}}{\delta \varepsilon }}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1759d8b5d79d68f6f968cbe285a48388bf87f4a4)
The equation of motion for the excited state probability then reads
![{\displaystyle {\dot {P}}_{e}=-(\gamma _{2}+\gamma _{1})P_{e}+\gamma _{2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/810fd10bb0620205660d2e319a448919d6dff300)
Assuming the two-level system is initially in its excited state, the stationary state given by
is approached in an exponential decay according to
![{\displaystyle P_{e}(t)=P_{e}^{s}+{\frac {\gamma _{1}}{\gamma _{1}+\gamma _{2}}}\exp[-2(\gamma _{1}+\gamma _{2})t].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d7cf61919b3d70d575044b40be1c2fcbbd1a5fde)
This solution is in excellent agreement with numerical simulations of the full Schrödinger equation [4].
It is important to stress that this efficient description of such a complex system is only possible because our projection operator correctly captures the essential features of the dynamics. If one choses a different projection operator, say, one without correlations between system and environments, according to
![{\displaystyle {\mathcal {P}}\rho ={\text{Tr}}_{E}\{\rho \}\otimes \rho _{E},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/14f0a3c5a5379fed584a54c03184a275170f4b5f)
then the second order TCL master equation no longer correctly describes the dynamics, as higher orders are divergent [4].
- ↑ Breuer, H.-P.; Petruccione, F. (2002). The Theory of Open Quantum Systems. Oxford University Press.
- ↑ Cite error: Invalid
<ref>
tag; no text was provided for refs named Weimer2008
- ↑ Cite error: Invalid
<ref>
tag; no text was provided for refs named Schwabl2010
- ↑ 4.0 4.1 Cite error: Invalid
<ref>
tag; no text was provided for refs named Breuer2006