Abstract We present a fluctuation theorem for Floquet quantum master equations. This is a detailed version of the famous Gallavotti–Cohen theorem. In contrast to the latter theorem, which involves the probability distribution of the total heat current, the former involves the joint probability distribution of positive and negative heat currents and can be used to derive the latter. A quantum two-level system driven by a periodic external field is used to verify this result. Keywords:Floquet quantum master equations;Gallavotti–Cohen fluctuation theorem;large deviation principle;positive and negative heat currents
The Gallavotti–Cohen (GC) fluctuation theorem states that the ratio of the probability distribution p(Σ) of having an average total entropy production rate Σ to that of having −Σ approaches $\exp (t\sigma )$ as the time interval t increases; that is,$ \begin{eqnarray}\displaystyle \frac{p(\sigma )}{p(-\sigma )}\asymp {{\rm{e}}}^{t\sigma },\end{eqnarray}$where the Boltzmann constant kB is set to 1 throughout this paper, and the symbol ≍denotes asymptotic change as $t\to \infty $ [1]. The original GC fluctuation theorem was inspired by a relationship between the probabilities of fluctuations in the shear stress of fluids in nonequilibrium steady states [2] and was proved in modern dynamical system theory [3, 4]. Complicated techniques were used; see the latest review [5]. In contrast, its proof in stochastic systems, e.g. in classical Langevin systems or discrete jump systems, is simple [6–9]. Hence, stochastic dynamical systems are suitable to explore new fluctuation theorems.
The motivation of this paper is as follows. Let us imagine that a system contacts a heat bath having an inverse temperature β and is in a nonequilibrium steady state due to some external force. The GC theorem can be reexpressed in terms of the probability distribution of the total heat current j, the total released heat averaged over the time interval t. Then, the exponent on the right-hand side of equation (1) is replaced by $t\beta j$. We know that the total heat is composed of positive and negative parts. Accordingly, the total heat current j can be divided into ${j}_{+}$ and ${j}_{-}$. Apparently, the individual current does not satisfy the GC fluctuation theorem. However, is the same true of the joined currents? Here, a positive answer is presented for the driven quantum systems described by Floquet quantum master equations [10–14].
The rest of this paper is organized as follows. In section 2, we review the Floquet quantum master equation and its stochastic thermodynamics. In section 3, we prove a fluctuation theorem. In section 4, a two-level quantum system is used to concretely verify this theorem. Section 5 concludes the paper.
2. Floquet quantum master equation and its stochastic thermodynamics
Given the Hamiltonian of a quantum system driven by periodic external forces, denoted as H(t), we have$ \begin{eqnarray}H\left(t+{ \mathcal T }\right)=H(t),\end{eqnarray}$where ${ \mathcal T }=2\pi /{\rm{\Omega }}$ is the periodicity and ω is the driving frequency. According to the Floquet theorem [15, 16], this periodic Hamiltonian satisfies an eigenvalue equation:$ \begin{eqnarray}(H(t)-{\rm{i}}{\partial }_{t})| {u}_{n}(t)\rangle ={\epsilon }_{n}| {u}_{n}(t)\rangle ,\end{eqnarray}$where εn and $| {u}_{n}(t)\rangle $ (n=1, ⋯, N) are quasi-energies and Floquet bases, respectively, and we set ℏ=1. Note that the Floquet bases are orthonormal and periodic. In addition, we emphasize that these quasi-energies are restricted in a zone with a size of ω. The heat bath that interacts with the quantum system has an inverse temperature β. Under the weak system-bath coupling conditions and time scale separation assumptions, the evolution of the reduced density matrix of the quantum system ρ(t) can be described by the Floquet quantum master equation [10, 17, 18]:$ \begin{eqnarray}{\partial }_{t}\rho (t)=-{\rm{i}}[H(t),\rho (t)]+D(t)[\rho (t)].\end{eqnarray}$The D-term represents the dissipation induced by the interaction between the system and the heat bath and is$ \begin{eqnarray}\begin{array}{rcl}D(t)[\rho ] & = & \sum _{\omega }r(\omega )\left[A(\omega ,t)\rho {A}^{\dagger }(\omega ,t)\right.\\ & & \left.-\displaystyle \frac{1}{2}\left\{{A}^{\dagger }(\omega ,t)A(\omega ,t),\rho \right\}\right].\end{array}\end{eqnarray}$In the above equation, ω are Bohr frequencies and are equal to ${\epsilon }_{n}-{\epsilon }_{m}+q{\rm{\Omega }}$, where q are certain integers. The numbers may be positive or negative but always appear in pairs. Additionally, in the same equation, A(ω, t) and ${A}^{\dagger }(\omega ,t)$ are called the Lindblad operators and are related by$ \begin{eqnarray}A(-\omega ,t)={A}^{\dagger }(\omega ,t).\end{eqnarray}$The interaction operator of the quantum system and the heat bath is given as A⨂B, where A and B are the system and heat bath components, respectively. These Lindblad operators are obtained by performing a Fourier-like expansion of the interaction-picture operator of A [17]:$ \begin{eqnarray}A(\omega ,t)=\sum _{m,n,q}{\delta }_{\omega ,{\epsilon }_{n}-{\epsilon }_{m}+q{\rm{\Omega }}}\langle \langle {u}_{m}| A| {u}_{n}{ \rangle\rangle }_{q}| {u}_{m}(t)\rangle \langle {u}_{n}(t)| {{\rm{e}}}^{-{\rm{i}}q{\rm{\Omega }}t},\end{eqnarray}$where δ is the Kronecker symbol, and the time-independent coefficient $\langle \langle {u}_{m}| A| {u}_{n}{ \rangle\rangle }_{q}$ is the qth harmonic of the transition amplitude $\langle {u}_{m}(t)| A| {u}_{n}(t)\rangle ;$ that is,$ \begin{eqnarray}\langle \langle {u}_{m}| A| {u}_{n}{ \rangle\rangle }_{q}=\displaystyle \frac{1}{{ \mathcal T }}{\int }_{0}^{{ \mathcal T }}{\rm{d}}t\langle {u}_{m}(t)| A| {u}_{n}(t)\rangle {{\rm{e}}}^{{\rm{i}}q{\rm{\Omega }}t}.\end{eqnarray}$The last ingredient of the Floquet quantum master equation is the assumption that the heat bath is always in a thermal state. Then, r(ω), the Fourier transformation of the correlation function of the operator B, satisfies the important Kubo–Martin–Schwinger condition [19, 20]:$ \begin{eqnarray}r(-\omega )=r(\omega ){{\rm{e}}}^{-\beta \omega }.\end{eqnarray}$
Stochastic thermodynamics can be established for the Floquet quantum master equation [14, 21, 22]. Roughly, equation (4) is unraveled into the dynamics of individual quantum systems [23–26]. The evolution of each system is composed of a continuous process alternating with discrete random jumps. Assume that these jumps occur at time points ti (i=1, ⋯). Each jump indicates that a quantum ωi is released to the heat bath [21, 22, 27–31]. The subscript i represents the time points of these energy exchanges. When the evolution of a quantum system ends at time t, a quantum jump trajectory is generated and is marked as $\overrightarrow{\omega }=\{{\omega }_{1},\,\cdots \}$. If the density matrixes of these individual quantum systems are $\widetilde{\rho }(\overrightarrow{\omega },t)$, their average weighted by the probabilities of all possible quantum jump trajectories is just the reduced density matrix ρ(t) of equation (4). From a thermodynamic point of view, the quanta are the heat released to the heat bath. Hence, given a quantum jump trajectory $\overrightarrow{\omega }$, we define the total heat along it as$ \begin{eqnarray}Q\{\overrightarrow{\omega }\}=\sum _{i=1}{\omega }_{i}={Q}_{+}\{\overrightarrow{\omega }\}+{Q}_{-}\{\overrightarrow{\omega }\}.\end{eqnarray}$In the second equation, we specifically define the positive and negative heat, ${Q}_{+}\{\vec{\omega }\}$ and ${Q}_{-}\{\vec{\omega }\}$, respectively; clearly, they are simply equal to the sums of positive and negative Bohr frequencies.
3. A fluctuation theorem
Because the occurrences of quantum jump trajectories and time points of quantum jumps are random events, all three types of heat are stochastic quantities. Let the joint probability distribution of the positive and negative heat, ${Q}_{+}$ and ${Q}_{-}$, respectively, be $p({Q}_{+},{Q}_{-})$. We can construct its histogram by directly simulating quantum jump trajectories [21, 22]. Because we are interested in the statistics over long time limits, a more practical approach is to compute its moment generation function,$ \begin{eqnarray}\begin{array}{rcl}{\rm{\Phi }}({\chi }_{+},{\chi }_{-}) & = & \displaystyle \int {\rm{d}}{Q}_{+}{\rm{d}}{Q}_{-}p({Q}_{+},{Q}_{-}){{\rm{e}}}^{{\chi }_{+}{Q}_{+}+{\chi }_{-}{Q}_{-}}\\ & = & \mathrm{Tr}[\hat{\rho }(t)].\end{array}\end{eqnarray}$We introduce a characteristic operator $\hat{\rho }(t)$ in the second equation above and find that it satisfies an evolution equation:$ \begin{eqnarray}{\partial }_{t}\hat{\rho }(t)=-{\rm{i}}[H(t),\hat{\rho }]+D(t,{\chi }_{+},{\chi }_{-})[\hat{\rho }(t)],\end{eqnarray}$where the super-operator$ \begin{eqnarray}\begin{array}{l}D(t,{\chi }_{+},{\chi }_{-})[\hat{\rho }]\\ \quad =\,\displaystyle \sum _{\omega \gt 0}r(\omega )\left[{{\rm{e}}}^{{\rm{i}}{\chi }_{+}\omega }A(\omega ,t)\hat{\rho }{A}^{\dagger }(\omega ,t)\right.\\ \qquad \left.-\,\displaystyle \frac{1}{2}\left\{{A}^{\dagger }(\omega ,t)A(\omega ,t),\hat{\rho }\right\}\right]\\ \qquad +\,\displaystyle \sum _{\omega \lt 0}r(\omega )\left[{{\rm{e}}}^{{\rm{i}}{\chi }_{-}\omega }A({\omega }_{,}t)\hat{\rho }{A}^{\dagger }(\omega ,t)\right.\\ \qquad \left.-\,\displaystyle \frac{1}{2}\left\{{A}^{\dagger }(\omega ,t)A(\omega ,t),\hat{\rho }\right\}\right].\end{array}\end{eqnarray}$If χ±=0, equation (13) reduces to equation (4). This result comes from a simple extension of the previous equation (equation (19) in [21]), which concerned the moment-generating function of the total heat, and we can reobtain it by letting ${\chi }_{+}={\chi }_{-}$ in equation (13). Because the derivation is the same, we do not repeat it here.
The abstract equation (12) is not the most convenient to use in analyses. According to equation (11), ${\rm{\Phi }}({\chi }_{+},{\chi }_{-})$ is equal to a sum of diagonal elements of $\hat{\rho }(t)$. Hence, we write the evolution equations for ${P}_{n}(t)=\langle {u}_{n}(t)| \hat{\rho }(t)| {u}_{n}(t)\rangle $ in the Floquet bases:$ \begin{eqnarray}\displaystyle \frac{{\rm{d}}}{{\rm{d}}t}{\boldsymbol{P}}(t)={\boldsymbol{R}}({\chi }_{+},{\chi }_{-}){\boldsymbol{P}}(t),\end{eqnarray}$where the vectors ${\boldsymbol{P}}{(t)=({P}_{1}(t),\cdots ,{P}_{N}(t))}^{{\rm{T}}}$ and ${\rm{T}}$ represents the transpose, the nondiagonal elements of the ${\boldsymbol{R}}$-matrix are$ \begin{eqnarray}\begin{array}{l}{\left[{\boldsymbol{R}}({\chi }_{+},{\chi }_{-})\right]}_{{nm}}\\ \quad =\,\displaystyle \sum _{\omega \gt 0}{{\rm{e}}}^{{\chi }_{+}\omega }r(\omega )| \langle {u}_{n}(t)| A(\omega ,t)| {u}_{m}(t)\rangle {| }^{2}\\ \qquad +\,\displaystyle \sum _{\omega \lt 0}{{\rm{e}}}^{{\chi }_{-}\omega }r(\omega )| \langle {u}_{n}(t)| A(\omega ,t)| {u}_{m}(t)\rangle {| }^{2},\end{array}\end{eqnarray}$($m\ne n$), and the diagonal elements are$ \begin{eqnarray}\begin{array}{l}{\left[{\boldsymbol{R}}({\chi }_{+},{\chi }_{-})\right]}_{{nn}}\\ =\,\displaystyle \sum _{\omega \gt 0}{{\rm{e}}}^{{\chi }_{+}\omega }r(\omega =q{\rm{\Omega }})| \langle {u}_{n}(t)| A(\omega ,t)| {u}_{n}(t)\rangle {| }^{2}\\ \ \ +\,\displaystyle \sum _{\omega \lt 0}{{\rm{e}}}^{{\chi }_{-}\omega }r(\omega =q{\rm{\Omega }})| \langle {u}_{n}(t)| A(\omega ,t)| {u}_{n}(t)\rangle {| }^{2}\\ \ \ -\,\displaystyle \sum _{\omega }r(\omega )\displaystyle \sum _{m}| \langle {u}_{m}(t)| A(\omega ,t)| {u}_{n}(t)\rangle {| }^{2}.\end{array}\end{eqnarray}$Note that equation (7) reminds us that ${\boldsymbol{R}}({\lambda }_{+},{\lambda }_{-})$ is in fact a constant matrix. We can easily prove that this matrix possesses symmetry:$ \begin{eqnarray}{\left[{\boldsymbol{R}}({\chi }_{+},{\chi }_{-})\right]}^{{\rm{T}}}={\boldsymbol{R}}(-\beta -{\chi }_{-},-\beta -{\chi }_{+}).\end{eqnarray}$Using this property, we obtain a fluctuation theorem by simply following a standard procedure, e.g. that presented by Lebowitz and Spohn [7]. First, because the transpose matrix has the same eigenvalues as the original matrix, equation (17) implies that the scaled cumulant generating function [1], $\phi ({\chi }_{+},{\chi }_{-})$, or the maximal eigenvalue of the ${\boldsymbol{R}}$-matrix has the same symmetry:$ \begin{eqnarray}\phi ({\chi }_{+},{\chi }_{-})=\phi (-\beta -{\chi }_{-},-\beta -{\chi }_{+}).\end{eqnarray}$Given the large deviation function $I({j}_{+},{j}_{-})$ of the distribution $p({j}_{+},{j}_{-})$ for the positive heat current ${j}_{+}={Q}_{+}/t$ and the negative heat current ${j}_{-}={Q}_{-}/t$, and because the function is a Legendre transform of the scaled cumulant generating function, equation (18) immediately leads to$ \begin{eqnarray}I({j}_{+},{j}_{-})=I(-{j}_{-},-{j}_{+})-\beta ({j}_{+}+{j}_{-}).\end{eqnarray}$Then, the probability distribution for these two heat currents satisfies the fluctuation theorem$ \begin{eqnarray}\displaystyle \frac{p({j}_{+},{j}_{-})}{p(-{j}_{-},-{j}_{+})}\asymp {{\rm{e}}}^{t\beta ({j}_{+}+{j}_{-})}.\end{eqnarray}$The conventional GC fluctuation theorem (1) can be easily derived from equation (20).
Fluctuation theorem (20) has a time-reversal explanation. Analogous to that of classical stochastic processes [7], the ratio of the probability distribution ${ \mathcal P }\{\overrightarrow{\omega }\}$ of observing a quantum jump trajectory $\overrightarrow{\omega }=\{{\omega }_{1},\,\cdots ,\,{\omega }_{M}\}$ to ${ \mathcal P }\{\overleftarrow{\omega }\}$ of observing its reversed trajectory $\overleftarrow{\omega }=\{-{\omega }_{M},\,\cdots ,\,-{\omega }_{1}\}$ approaches $\exp (\beta Q\{\overrightarrow{\omega }\})$ as the time interval t increases [22]. Note that not only is the time order of the quantum jumps reversed in the revered trajectory, but the signs of these Bohr frequencies are reversed, we obtain ${Q}_{\pm }\{\overleftarrow{\omega }\}=-{Q}_{\mp }\{\overrightarrow{\omega }\}$ and$ \begin{eqnarray}\begin{array}{rcl}p({j}_{+},{j}_{-}) & = & \displaystyle \int { \mathcal D }\overrightarrow{\omega }{ \mathcal P }\{\overrightarrow{\omega }\}\delta ({Q}_{+}\{\overrightarrow{\omega }\}-{Q}_{+})\\ & & \times \delta ({Q}_{-}\{\overrightarrow{\omega }\}-{Q}_{-})\\ & \asymp & \displaystyle \int { \mathcal D }\overrightarrow{\omega }{ \mathcal P }\{\overleftarrow{\omega }\}{{\rm{e}}}^{\beta ({Q}_{+}\{\overrightarrow{\omega }\}+{Q}_{-}\{\overrightarrow{\omega }\})}\\ & & \times \delta ({Q}_{+}\{\overrightarrow{\omega }\}-{Q}_{+})\delta ({Q}_{-}\{\overrightarrow{\omega }\}-{Q}_{-})\\ & \asymp & \displaystyle \int { \mathcal D }\overleftarrow{\omega }{ \mathcal P }\{\overleftarrow{\omega }\}{{\rm{e}}}^{-\beta ({Q}_{-}\{\overleftarrow{\omega }\}+{Q}_{+}\{\overleftarrow{\omega }\})}\\ & & \times \delta ({Q}_{+}\{\overleftarrow{\omega }\}+{Q}_{-})\delta ({Q}_{-}\{\overleftarrow{\omega }\}+{Q}_{+})\\ & = & {{\rm{e}}}^{t\beta ({j}_{+}+{j}_{-})}p(-{j}_{-},-{j}_{+}).\end{array}\end{eqnarray}$Some unimportant constants are ignored here. This rough proof explains why the positions of ${j}_{+}$ and ${j}_{-}$ in the probability distribution on the right-hand side are exchanged and minus signs are added simultaneously; they are just the consequences of the time-reversal. Equations (19) and(20) are the central results of this paper. In the next section, we use a two-level quantum system to concretely show these results.
4. Two-level quantum system
Consider the Hamiltonian of a two-level quantum system [12–14, 17, 32]$ \begin{eqnarray}H(t)=\displaystyle \frac{1}{2}{\omega }_{0}{\sigma }_{z}+\displaystyle \frac{1}{2}{{\rm{\Omega }}}_{R}\left({\sigma }_{+}{{\rm{e}}}^{-{\rm{i}}{\rm{\Omega }}t}+{\sigma }_{-}{{\rm{e}}}^{{\rm{i}}{\rm{\Omega }}t}\right),\end{eqnarray}$where ω0 is the transition frequency of the bare system, ωR is the Rabi frequency, and ω is the frequency of the periodic external field. The Floquet bases and the quasi-energies of this system are$ \begin{eqnarray}| {u}_{\pm }(t)\rangle =\displaystyle \frac{1}{\sqrt{2{\rm{\Omega }}^{\prime} }}\left(\begin{array}{c}\pm \sqrt{{\rm{\Omega }}^{\prime} \pm \delta }\\ {{\rm{e}}}^{{\rm{i}}{\rm{\Omega }}t}\sqrt{{\rm{\Omega }}^{\prime} \mp \delta },\end{array}\right),\end{eqnarray}$and ${\epsilon }_{\pm }=({\rm{\Omega }}\pm {\rm{\Omega }}^{\prime} )/2$, respectively, where ${\rm{\Omega }}^{\prime} =\sqrt{{\delta }^{2}+{{\rm{\Omega }}}_{{\rm{R}}}^{2}}$ and the detuning parameter $\delta ={\omega }_{0}-{\rm{\Omega }}$. Here, we additionally set ${\rm{\Omega }}\gt {\rm{\Omega }}^{\prime} $. Assume that the coupling between the two-level system and the heat bath is Σx-coupling. There are six Lindblad operators, and the Bohr frequencies are $\pm {\rm{\Omega }}$, $\pm ({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )$, and $\pm ({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )$. By performing some simple derivations, the ${\boldsymbol{R}}$-matrix is obtained:$ \begin{eqnarray}\begin{array}{rcl}{\left({\boldsymbol{R}}\right)}_{11} & = & ({{\rm{e}}}^{{\rm{i}}{\chi }_{+}{\rm{\Omega }}}-1){{\rm{\Gamma }}}_{+{\rm{\Omega }}}+({{\rm{e}}}^{-{\rm{i}}{\chi }_{-}{\rm{\Omega }}}-1){{\rm{\Gamma }}}_{-{\rm{\Omega }}}\\ & & -{{\rm{\Gamma }}}_{-({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )}-{{\rm{\Gamma }}}_{+({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )},\\ {\left({\boldsymbol{R}}\right)}_{12} & = & {{\rm{e}}}^{{\rm{i}}{\chi }_{+}({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )}{{\rm{\Gamma }}}_{+({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )}+{{\rm{e}}}^{-{\rm{i}}{\chi }_{-}({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )}{{\rm{\Gamma }}}_{-({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )},\\ {\left({\boldsymbol{R}}\right)}_{21} & = & {{\rm{e}}}^{-{\rm{i}}{\chi }_{-}({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )}{{\rm{\Gamma }}}_{-({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )}+{{\rm{e}}}^{{\rm{i}}{\chi }_{+}({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )}{{\rm{\Gamma }}}_{+({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )},\\ {\left({\boldsymbol{R}}\right)}_{22} & = & ({{\rm{e}}}^{{\rm{i}}{\chi }_{+}{\rm{\Omega }}}-1){{\rm{\Gamma }}}_{+{\rm{\Omega }}}+({{\rm{e}}}^{-{\rm{i}}{\chi }_{1}^{-}{\rm{\Omega }}}-1){{\rm{\Gamma }}}_{-{\rm{\Omega }}}\\ & & -{{\rm{\Gamma }}}_{+({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )}-{{\rm{\Gamma }}}_{-({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )},\end{array}\end{eqnarray}$where we do not explicitly write ${\chi }_{\pm }$ on the left-hand side and the coefficients are$ \begin{eqnarray}\begin{array}{rcl}{{\rm{\Gamma }}}_{\pm {\rm{\Omega }}} & = & {\left(\displaystyle \frac{{{\rm{\Omega }}}_{{\rm{R}}}}{2{\rm{\Omega }}^{\prime} }\right)}^{2}r(\pm {\rm{\Omega }}),\\ {{\rm{\Gamma }}}_{\pm ({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )} & = & {\left(\displaystyle \frac{\delta -{{\rm{\Omega }}}_{{\rm{R}}}}{2{\rm{\Omega }}^{\prime} }\right)}^{2}r(\pm ({\rm{\Omega }}-{\rm{\Omega }}^{\prime} )),\\ {{\rm{\Gamma }}}_{\pm ({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )} & = & {\left(\displaystyle \frac{\delta +{{\rm{\Omega }}}_{{\rm{R}}}}{2{\rm{\Omega }}^{\prime} }\right)}^{2}r(\pm ({\rm{\Omega }}+{\rm{\Omega }}^{\prime} )).\end{array}\end{eqnarray}$In the matrix, symmetry(17) is apparent. Because this is a simple 2×2 matrix, the analytical expression of the scaled cumulant generating function is$ \begin{eqnarray}\phi ({\chi }_{+},{\chi }_{-})=\displaystyle \frac{1}{2}[{\left({\boldsymbol{R}}\right)}_{11}+{\left({\boldsymbol{R}}\right)}_{22}\pm B],\end{eqnarray}$where$ \begin{eqnarray}B=\sqrt{{[{\left({\boldsymbol{R}}\right)}_{11}-{\left({\boldsymbol{R}}\right)}_{22}]}^{2}+4{\left({\boldsymbol{R}}\right)}_{12}{\left({\boldsymbol{R}}\right)}_{21}}.\end{eqnarray}$We clearly see that symmetry(17) is inherited in equation (26).
We compute the large deviation function $I({j}_{+},{j}_{-})$ by numerically performing a Legendre transformation on equation (26) and depict it in figure 1. Equation (19) can be easily verified (data not shown). On the other hand, it is interesting to approximately compute the large deviation function by simulating quantum jump trajectories with finite times; see the spheres in the same figure. We find that even if the simulation time is short, the simulated data roughly exhibit the profile of the function. Although longer simulation times lead to better results, the regimes of the sampled data decrease around the mean currents $({\overline{j}}_{+},{\overline{j}}_{-})$. This observation reminds us that simply simulating quantum jump trajectories is not enough to solve the full large deviation function.
Figure 1.
New window|Download| PPT slide Figure 1.The large deviation functions (LDF) $I({j}_{+},{j}_{-})$ for a two-level Floquet quantum system. The meshed surface is obtained by performing a Legendre transformation on equation (26). The spheres are the data of simulating quantum jump trajectories, where the simulation times for the green and blue data are 50 and 500, respectively. The Fourier transformation of the correlation function is set to be $r(\omega )={ \mathcal A }| \omega {| }^{3}{{ \mathcal N }}_{k}(\omega )$ for ω<0; otherwise, $r(\omega )={ \mathcal A }| \omega {| }^{3}[{{ \mathcal N }}_{k}(\omega )+1]$, where ${ \mathcal N }(\omega )=1/[\exp (\beta | \omega | )-1]$, and the coefficient A is related to the coupling strength between the system and the heat bath [19]. The parameters used are ω0=1, ωR=0.8, ω=1.1, ${ \mathcal A }=1$, and β=1/3. The red dashed line indicates the location of the two mean currents.
5. Conclusion
In this paper, we present a detailed GC fluctuation theorem for the quantum systems described by the Floquet quantum master equations. Although this theorem is proved for systems that interact with one heat bath, its generalization to the case of multiple heat baths is straightforward. For instance, for the case of two heat baths, we have [33]$ \begin{eqnarray}\displaystyle \frac{p({j}_{1}^{+},{j}_{1}^{-},{j}_{2}^{+},{j}_{2}^{-})}{p(-{j}_{1}^{-},-{j}_{1}^{+},-{j}_{2}^{-},-{j}_{2}^{+})}\asymp {{\rm{e}}}^{t{\beta }_{1}({j}_{1}^{+}+{j}_{1}^{-})+t{\beta }_{2}({j}_{2}^{+}+{j}_{2}^{-})},\end{eqnarray}$where ${j}_{k}^{+}$ and ${j}_{k}^{-}$ are the positive and negative heat currents of the quantum system released to the heat bath with the inverse temperature ${\beta }_{k}$ (k = 1, 2). Finally, if a quantum system contacts two heat baths at two different temperatures and can be described by a Lindblad quantum master equation, the system is able to evolve into a nonequilibrium steady state without external driving fields. In such a situation, by carrying out the same argument presented here, we can prove that the fluctuation theorem (28) is still true.
Acknowledgments
This work was supported by the National Science Foundation of China under Grants No. 11 174 025 and No. 11 575 016. We also appreciate the support of the CAS Interdisciplinary Innovation Team, No. 2060 299.