This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Showcasing extractable work of quantum battery via entropic uncertainty relations

Meng-Long Song School of Physics & Optoelectronic Engineering, Anhui University, Hefei 230601, People’s Republic of China    Xue-Ke Song School of Physics & Optoelectronic Engineering, Anhui University, Hefei 230601, People’s Republic of China    Liu Ye School of Physics & Optoelectronic Engineering, Anhui University, Hefei 230601, People’s Republic of China    Dong Wang dwang@ahu.edu.cn School of Physics & Optoelectronic Engineering, Anhui University, Hefei 230601, People’s Republic of China
Abstract

In this study, we investigate the effectiveness of entropic uncertainty relations (EURs) in discerning the energy variation in quantum batteries (QBs) modelled by battery-charger-field in the presence of bosonic and fermionic reservoirs. Our results suggest that the extractable works (exergy and ergotropy) have versatile characteristics in different scenarios, resulting in a complex relationship between tightness and extractable work. It is worth noting that the tightness of the lower bound of entropic uncertainty can be a good indicator for energy conversion efficiency in charging QBs. Furthermore, we disclose how the EUR including uncertainty and lower bound contributes to energy conversion efficiency in the QB system. It is believed that these findings will be beneficial for better understanding the role of quantum uncertainty in evaluating quantum battery performance.

I Introduction

The storage and utilisation of energy have always been an important topic worldwide and can be achieved by devices, that is, so-called batteries. In contrast to the conversion of chemical energy into electric energy in the classical battery scheme, quantum batteries (QBs) review are regarded as a new and promising energy storage process, pioneeringly presented by Alicki and Fannes propose . The core principle of quantum batteries is using quantum resources (e.g., quantum entanglement and coherence) for energy storage resource1 ; resource2 ; resource3 . With the development of quantum information theory, quantum batteries have received considerable attention, and progress has been made in theory, including the construction of charging models model1 ; model2 ; model3 ; model4 ; model5 ; model6 ; model7 , open QB systems open1 ; open2 ; open3 ; open4 ; open5 ; open6 ; open7 ; open8 ; Dou-F-Q1 , many-body QBs (e.g., Dicke quantum battery and spin-based batteries) model4 ; r1 ; cells1 ; cells2 ; cells3 ; cells4 ; Dou-F-Q2 , and so on. Several promising experiments have been conducted for QBs in organic microcavities experiment1 , superconducting experiment2 , quantum dots experiment4 and NMR experiment3 platforms.

The great contribution of quantum resources is to allow QBs to have superextensive scaling of the charging power in many-body QB systems, such as entanglement of battery subspace resulting in quantum advantages over classical schemes review ; model1 ; model2 . However, many investigations review ; open8 ; Dou-F-Q2 ; x18 ; x18++ ; capacity ; measurement have revealed that quantum sources cannot independently and adequately predict the energy variation in QBs. This is due to several reasons, such as non-equilibrium thermodynamic processes in which multiple quantum resources alternately or simultaneously act on the battery. Therefore, it is argued that the analysis of QB’s performance needs more indicators in the specific charging protocol. Conversely, Gu¨\rm\ddot{u}hne and Lewenstein have studied the intrinsic relationship between entanglement and entropic uncertainty relations (EURs) ee . With this in mind, it is natural to ask whether EURs can be used to predict the change in quantum battery’s energy. Upon studying the relationship between QB’s energy and EURs, we found the answer to be positive.

The uncertainty relation, originally proposed by Heisenberg, is the cornerstone of quantum information theory. It proposes that one cannot accurately predict the measurement outcomes with respect to two arbitrary incompatible observables 1 . Afterwards, Kennard 2 and Robertson 3 generalized into the inequality for arbitrary two operators: ΔM^1ΔM^2|[M^1,M^2]|/2\Delta{\hat{M}_{1}}\cdot\Delta{\hat{M}_{2}}\geq|\langle{[{\hat{M}_{1}},{\hat{M}_{2}}]}\rangle|/2 in terms of standard deviation, with Δ()\Delta(\bullet) denoting the standard deviation and \langle\circ\rangle denoting the expected value. Eversince, there have been various quantifications of uncertainty relations, including Shannon entropy 4 , Rényi entropy 6 , and von Neumann entropy JR . Furthermore, Berta et al. 7 presented quantum-memory-assisted entropic uncertainty relations (QMA-EURs). Liu etet al.al. novelly derived an expression for QMA-EURs as multiple observables liu

x=1mS(M^x|B)log21f+(m1)S(A|B),\displaystyle\sum\limits_{x=1}^{m}{S\left({{{\hat{M}}_{x}}\left|B\right.}\right)}\geq{\log_{2}}\frac{1}{f}+\left({m-1}\right)S\left({A\left|B\right.}\right), (1)

herein, the von Neumann conditional entropy S(A|B)=S(ρ^AB)S(ρ^B)S({A|B})=S({{{\hat{\rho}}_{AB}}})-S\left({{{\hat{\rho}}_{B}}}\right) m ; x ; wu-lin with S(ρ^)=tr(ρ^log2ρ^)S\left({\hat{\rho}}\right)=-tr\left({\hat{\rho}{{\log}_{2}}\hat{\rho}}\right), and f=maxiN{i2in1maxi1[c(ui11,ui22)]x=2n1[c(uixx,uix+1x+1)]}f=\mathop{\max}\limits_{{i_{N}}}\left\{{\sum\limits_{{i_{2}}\sim{i_{n-1}}}{\mathop{\max}\limits_{{i_{1}}}}\left[{c\left({u_{{i_{1}}}^{1},u_{{i_{2}}}^{2}}\right)}\right]\prod\limits_{x=2}^{n-1}{\left[{c\left({u_{{i_{x}}}^{x},u_{{i_{x+1}}}^{x+1}}\right)}\right]}}\right\}, where uixxu_{i_{x}}^{x} is the iith eigenvector of the operator M^x{{{\hat{M}}_{x}}} and c(uixx,uix+1x+1){c\left({u_{{i_{x}}}^{x},u_{{i_{x+1}}}^{x+1}}\right)} denotes the maximal overlap of the operators uixx{u_{{i_{x}}}^{x}} and uix+1x+1u_{{i_{x+1}}}^{x+1}.

In this study, we focused on the effectiveness of EURs in discerning the energy variation for QBs modelled by a charger-mediated protocol with an external field and reservoir. We explored the relationship between the extractable work and the EUR of the QB. Compared with previous results reference1 ; x18 ; x18++ ; capacity , the merits of the current work can be summarised as follows: (i) Improving the extractable work (i.e., exergy and ergotropy) and conversion efficiency of the battery was revealed. (ii) We found that EUR tightness could serve as an excellent indictor for predicting extractable work and conversion efficiency. (iii) We revealed how the corresponding EURs contribute to the conversion efficiency of QB.

The remainder of this paper is organised as follows. Section II introduces a typical quantum-battery model with a battery-charger-field in the presence of reservoirs. In Section III, in addition to the performance of the QB, the relationship between extractable work and tightness is investigated. In Section IV, the mechanism for improving extractable work via EUR is provided. Finally, our article ends with a brief discussion and conclusions in Section V.

II The model of quantum battery

The diagrammatic charging model of a QB consists of a battery and a two-level charger. Simultaneously, the charger is powered by a driving laser field and coupled to a bosonic or fermionic reservoir. As illustrated in Fig. 1, the energy can be transferred to the battery owing to the interaction between the charger and the battery. Consequently, the energy storage of the battery was achieved. In the rotating frame, the Hamiltonian of the charging system s{{\cal H}_{s}} can be written as reference1 ; reference2 ; reference3

Refer to caption
Figure 1: Schematic diagram of the charging process of the quantum battery. The charger (AA) is driven by a classical laser field with frequency ωd{\omega_{d}} and amplitude FF, and is coupled with the battery (BB) and interacts with its surrounding reservoir (RR). As a result, the energy is not only transferred from the charger to the battery (BB) but is also dissipated into the reservoir (RR) mediated by the charger. The two-level charger and the battery are characterized by transition frequencies ωA{\omega_{A}} and ωB{\omega_{B}}, respectively.
s=\displaystyle{{\cal H}_{s}}= ΔA2σ^zA+ΔB2σ^zB+g(σ^+Aσ^B+σ^Aσ^+B)\displaystyle\frac{{{\Delta_{A}}}}{2}\hat{\sigma}_{z}^{A}+\frac{{{\Delta_{B}}}}{2}\hat{\sigma}_{z}^{B}+{g}\left({\hat{\sigma}_{+}^{A}\hat{\sigma}_{-}^{B}+\hat{\sigma}_{-}^{A}\hat{\sigma}_{+}^{B}}\right)
+F2(σ^+A+σ^A),\displaystyle+\frac{F}{2}\left({\hat{\sigma}_{+}^{A}+\hat{\sigma}_{-}^{A}}\right), (2)

here, s{{\cal H}_{s}} denotes the Hamiltonian of the charger-battery system (the Boltzmann constant kB==1{{k_{B}}=\hbar=1} hereafter, for simplicity). Moreover, σ^jA,B(j=x,y,z)\hat{\sigma}_{j}^{A,B}\left({j=x,y,z}\right) are the Pauli operators and σ^kA,B(k=+,)\hat{\sigma}_{k}^{A,B}\left({k=+,-}\right) are the upper and down operators, Δl=ωlωd(l=A,B){\Delta_{l}}={\omega_{l}}-{\omega_{d}}\left({l=A,B}\right) is the detuning of a two-level system with a transition frequency of ωl{\omega_{l}} to the driving field with frequency ωd{\omega_{d}}, assume the charger has the same frequency as the battery, i.e., ωA=ωB=ω0{\omega_{A}}={\omega_{B}}={\omega_{0}}. g{g} is the coupling strength between the charger and battery, and FF is the amplitude of the driving field.

The Markovian master equation describing the charger-battery system can be obtained experiment1 ; reference1 ; reference2 ; measurement ; reference3

ρ^˙AB\displaystyle{\dot{\hat{\rho}}_{AB}} =i[s,ρ^AB]\displaystyle=-i\left[{{{\cal H}_{s}},{\hat{\rho}_{AB}}}\right]
+J{N(T)Dσ^A[ρ^AB]+n(T)Dσ^+A[ρ^AB]},\displaystyle+J\left\{{N\left(T\right){D_{\hat{\sigma}_{-}^{A}}}\left[{{\hat{\rho}_{AB}}}\right]+n\left(T\right){D_{\hat{\sigma}_{+}^{A}}}\left[{{\hat{\rho}_{AB}}}\right]}\right\}, (3)

where DΛ[ρ^]=2Λρ^ΛΛΛρ^ρ^ΛΛ{D_{\Lambda}}\left[\hat{\rho}\right]=2\Lambda\hat{\rho}{\Lambda^{\dagger}}-{\Lambda^{\dagger}}\Lambda\hat{\rho}-\hat{\rho}{\Lambda^{\dagger}}\Lambda denotes the dissipator, and JJ is the dissipative rate. Regarding the different types of reservoirs, the average particle number n(T)n(T) is different at frequency ωk{\omega_{k}} and temperature TT. Specifically, n(T)=nb=1/[exp(ωk/T)1]n\left(T\right)={n_{b}}=1/\left[{\exp\left({\omega_{k}/T}\right)-1}\right] and N(T)=1+n(T)N\left(T\right)=1+n\left(T\right) in the bosonic reservoir. For the fermionic one, n(T)=nf=1/[exp(ωkμ)/T+1]n\left(T\right)={n_{f}}=1/\left[{\exp\left({\omega_{k}-\mu}\right)/T+1}\right] and N(T)=1n(T)N\left(T\right)=1-n\left(T\right), where μ\mu denotes the chemical potential of the reservoir.

Suppose that the initial state of the charger-battery is prepared in the ground state |g\left|g\right\rangle (henceforth, the excited state is written as |e\left|e\right\rangle), that is,

ρ^AB(t0)=|ggABgg|.\displaystyle{\hat{\rho}_{AB}}\left({t\leq 0}\right)={\left|{gg}\right\rangle_{AB}}\left\langle{gg}\right|. (4)

Once the charging process begins (t0)\left({t\geq 0}\right), the driving field continues to power the charger, and the charger is connected to the reservoir because it interacts with the charger, whereas energy is transferred through the charger to the battery by interaction. The internal energy of the charger (AA) or battery (BB) at any time tt can be expressed as

Ei(t)=tr[0iρ^i(t)],\displaystyle{E_{i}}\left(t\right)={\rm tr}\left[{{{\cal H}_{0}^{i}}{\hat{\rho}_{i}}\left(t\right)}\right], (5)

where 0i=ω0σ^+iσ^i{{\cal H}_{0}^{i}}={\omega_{0}}{\hat{\sigma}_{+}^{i}}{\hat{\sigma}_{-}^{i}} and ρ^i(t){{\hat{\rho}_{i}}\left(t\right)} (i=Ai=A or BB) is the reduced density matrix of ρ^AB(t){{\hat{\rho}_{AB}}\left(t\right)}. In addition to the storage capacity of the QB, the extractable work is also an important factor in determining the QB’s performance. First, the extraction of useful work from a battery by thermalization (a nonunitary extraction process) can be quantified by the variations in the Helmholtz free energy of the battery. In other words, the battery is connected to a thermal bath with an inverse temperature β\beta. Consequently, the work can be extracted as exergy x18++ ; exergy1 ; exergy2

Wf(t)=(ρ^B)(ρ^β),\displaystyle{W_{f}}\left(t\right)=\mathcal{F}\left({{{\hat{\rho}}_{B}}}\right)-\mathcal{F}\left({{{\hat{\rho}}_{\beta}}}\right), (6)

where (ρ^B)=tr(0ρ^B)S(ρ^B)/β\mathcal{F}\left({\hat{\rho}_{B}}\right)={\rm tr}\left({{\mathcal{H}_{0}}\hat{\rho}_{B}}\right)-S\left({\hat{\rho}_{B}}\right)/\beta is Helmholtz free energy, and ρ^β=eβ0B/tr(eβ0B){{\hat{\rho}}_{\beta}}={e^{-\beta{\mathcal{H}_{0}^{B}}}}/{\rm tr}\left({{e^{-\beta{\mathcal{H}_{0}^{B}}}}}\right) is the thermal equilibrium state of the battery. To characterize the maximal amount of energy that can be extracted from a QB in the charging process under the cyclic unitary operations, the ergotropy reference4 ; reference5 is introduced as We(t)=tr[0Bρ^B(t)]tr[0Bρ^p]{W_{e}\left(t\right)}={\rm tr}\left[{{{\cal H}_{0}^{B}}{\hat{\rho}_{B}}\left(t\right)}\right]-{\rm tr}\left[{{{\cal H}_{0}^{B}}{\hat{\rho}_{p}}}\right] where ρ^p{{\hat{\rho}_{p}}} is called the passive state of ρ^B(t){{\hat{\rho}_{B}}}(t) with the zero extractable energy by cyclic unitary operations. In spectral decomposition, the Hamiltonian HH and density matrix ρ^\hat{\rho} are written as H=jej|ejej|H=\sum\nolimits_{j}{{e_{j}}}\left|{{e_{j}}}\right\rangle\left\langle{{e_{j}}}\right| and ρ^=jrj|rjrj|\hat{\rho}=\sum\nolimits_{j}{{r_{j}}}\left|{{r_{j}}}\right\rangle\left\langle{{r_{j}}}\right| respectively, where rj{{r_{j}}} and ej{{e_{j}}} are the eigenvalues of ρ^\hat{\rho} and HH with corresponding eigenstates |rj\left|{{r_{j}}}\right\rangle and |ej\left|{{e_{j}}}\right\rangle. After ordering eigenvalues (r0r1r2{r_{0}}\geq{r_{1}}\geq{r_{2}}\geq\cdots and e0e1e2{e_{0}}\leq{e_{1}}\leq{e_{2}}\cdots), the passive state can be derived as ρ^p=jrj|ejej|{\hat{\rho}_{p}}=\sum\nolimits_{j}{{r_{j}}}\left|{{e_{j}}}\right\rangle\left\langle{{e_{j}}}\right|. The ergotropy can be reexpressed as

We(t)=EB(t)jrjej.\displaystyle{W_{e}\left(t\right)}={E_{B}}\left(t\right)-\sum\limits_{j}{{r_{j}}}{e_{j}}. (7)

Furthermore, combining the density matrix of the charger-battery system ρ^AB(t){\hat{\rho}_{AB}}\left(t\right), EB(t){E_{B}}\left(t\right) and W(t)W\left(t\right) can be described as follows:

EB(t)\displaystyle{E_{B}}\left(t\right) =ω0[ρ11(t)+ρ33(t)],\displaystyle={\omega_{0}}\left[{{\rho_{11}}\left(t\right)+{\rho_{33}}\left(t\right)}\right],

and

We(t)\displaystyle{W_{e}\left(t\right)} =ω02{4|ρ12(t)+ρ34(t)|2+κ2+κ},\displaystyle=\frac{{\omega_{0}}}{2}\left\{{\sqrt{4{{\left|{{\rho_{12}}\left(t\right)+{\rho_{34}}\left(t\right)}\right|}^{2}}+{\kappa^{2}}}+\kappa}\right\}, (9)

where ρij(t)(i,j=1,2,3,4){\rho_{ij}}\left(t\right)\left({i,j=1,2,3,4}\right) is the ii-th row and jj-th column matrix entry of ρ^AB(t){\hat{\rho}_{AB}}\left(t\right) with κ=2[ρ11(t)+ρ33(t)]1\kappa=2\left[{{\rho_{11}}\left(t\right)+{\rho_{33}}\left(t\right)}\right]-1. Note that when charging for a sufficiently long period (i.e., tt\to\infty), the charger-battery system becomes steady, and the system’s energy does not vary over time tt. Hence, we have Ei(t)Ei(){E_{i}}\left(t\right)\to{E_{i}}\left(\infty\right), Wf(t)Wf(){W_{f}\left(t\right)}\to{W_{f}\left(\infty\right)} and We(t)We(){W_{e}\left(t\right)}\to{W_{e}\left(\infty\right)} (ρ^˙AB=0{\dot{\hat{\rho}}_{AB}}=0) in the limit. In fact, when S(ρ^β)S(ρ^p)>β[tr(0Bρ^β)tr(0Bρ^p)]S\left({{{\hat{\rho}}_{\beta}}}\right)-S\left({{{\hat{\rho}}_{p}}}\right)>\beta\left[{{\rm tr}\left({{\mathcal{H}_{0}^{B}}{{\hat{\rho}}_{\beta}}}\right)-{\rm tr}\left({{\mathcal{H}_{0}^{B}}{{\hat{\rho}}_{p}}}\right)}\right], the battery always stores more exergy than ergotropy. The left of the above inequality is the entropy production of the battery during the thermalization process, and the right is the change in battery energy caused by thermalization exergy2 ; exergy3 . In addition, the work that can be extracted Wf(t)[0,ΔEB(t)]{W_{f}\left(t\right)}\in\left[{0,\Delta{E_{B}}\left(t\right)}\right] is hold, where ΔEB(t)=EB(t)EB(0)\Delta{E_{B}}\left(t\right)={E_{B}}\left(t\right)-{E_{B}}\left(0\right) represents the stored energy of the battery.

As a matter of fact, the conversion rate of energy is termed as another nontrivial factor to evaluate the performance of QB, in addition to the amount of energy. Therefore, the conversion efficiency of the battery is given by

η1(t)=Wf(t)ΔEB(t),\displaystyle{\eta_{1}}\left(t\right)=\frac{{{W_{f}\left(t\right)}}}{{\Delta{E_{B}}\left(t\right)}},
η2(t)=We(t)ΔEB(t).\displaystyle{\eta_{2}}\left(t\right)=\frac{{{W_{e}\left(t\right)}}}{{\Delta{E_{B}}\left(t\right)}}. (10)

It quantifies the proportion of useful work in QB’s energy storage and deals with the energy cost review ; open1 ; cost3 of charging QBs (higher efficiency means lower energy costs). Therefore, a high-performance QB essentially needs to its higher efficiency to pursue more extractable work.

III Performance of quantum batteries

In general, quantum batteries with excellent performances have two significant advantages: stability and high power. In this consideration, we investigated the energy changes for the considered charging scenarios while the state of the charger-battery system was steady. Importantly, the role of the uncertainty’s tightness of QMA-EUR is discussed in detail while determining the acquirement and extraction of QB’s energy.

To probe the dynamics of entropic uncertainty, we use Pauli operators as the incompatible observables to be measured, which leads to post-measurement states (i=x,y,z)\left({i=x,y,z}\right)

ρ^(σ^i|B)=12[ρ^AB(t)+(σ^i𝕀)ρ^AB(t)(σ^i𝕀)],\displaystyle\hat{\rho}\left({{\hat{\sigma}_{i}}\left|B\right.}\right)=\frac{1}{2}\left[{{\hat{\rho}_{AB}}\left(t\right)+\left({{\hat{\sigma}_{i}}\otimes\mathbb{I}}\right){\hat{\rho}_{AB}}\left(t\right)\left({{\hat{\sigma}_{i}}\otimes\mathbb{I}}\right)}\right], (11)

where 𝕀\mathbb{I} denotes an identity matrix. It can be observed that different measurement choices for AA lead to different post-measurement states. Subsequently, we explored the entropic uncertainty relationship of the system under different measurement choices. When two measurement operators (σ^x,σ^z)\left({{\hat{\sigma}_{x}},{\hat{\sigma}_{z}}}\right) are selected, the entropic uncertainty (i.e., left-hand side of Eq. (1)) is expressed as

Ulxz(t)=j=x,zS[ρ^(σ^j|B)]2S[ρ^B(t)],\displaystyle U_{l}^{xz}\left(t\right)=\sum\limits_{j=x,z}{S\left[{\hat{\rho}\left({{{\hat{\sigma}}_{j}}\left|B\right.}\right)}\right]}-2S\left[{{{\hat{\rho}}_{B}}\left(t\right)}\right], (12)

and the corresponding lower bound (right-hand side of Eq.   (1)) can be expressed as

Urxz(t)=S[ρ^AB(t)]S[ρ^B(t)]+1.\displaystyle U_{r}^{xz}\left(t\right)=S\left[{{\hat{\rho}_{AB}}\left(t\right)}\right]-S\left[{{\hat{\rho}_{B}}\left(t\right)}\right]+1. (13)

Consequently, the uncertainty’s tightness is

Δxz(t)=Ulxz(t)Urxz(t).\displaystyle{\Delta^{xz}}\left(t\right)=U_{l}^{xz}\left(t\right)-U_{r}^{xz}\left(t\right). (14)

Incidentally, a smaller value of Δxz(t){\Delta^{xz}}\left(t\right), i.e., stronger tightness-implies a tighter lower bound, and vice versa. For the steady state, Ulxz(t)Ulxz()U_{l}^{xz}\left(t\right)\to U_{l}^{xz}\left(\infty\right) and Urxz(t)Urxz()U_{r}^{xz}\left(t\right)\to U_{r}^{xz}\left(\infty\right); thus, we have Δxz(t)Δxz(){\Delta^{xz}}\left(t\right)\to{\Delta^{xz}}\left(\infty\right) in terms of Eq. (14).

III.1 Resonance

To probe the energy-transfer issues of QBs, we first consider the case of eliminating detuning between the charger (battery) and driving field regarding the system’s energy variation, that is, ΔA=ΔB=0{\Delta_{A}}={\Delta_{B}}=0. A zero-temperature reservoir where T0T\to 0 leads to nb,nf0{n_{b}},{n_{f}}\to 0. Consequently, the state of the charging system ρAB(){\rho_{AB}}\left(\infty\right) can be described as

ρAB()|=T=ΔA=ΔB=01k4+2(2+k2)l2(k440ik3l2ik2l0k440ik3l2ik3l20k4+4k2l242kl2ik2lik3l22kl2k4+4(4+k2)l24),\displaystyle{\rho_{AB}}\left(\infty\right)\left|{{}_{T={\Delta_{A}}={\Delta_{B}}=0}}\right.=\frac{1}{{{k^{4}}+2\left({2+{k^{2}}}\right){l^{2}}}}\left({\begin{array}[]{*{20}{c}}{\frac{{{k^{4}}}}{4}}&0&{-\frac{{i{k^{3}}l}}{2}}&{i{k^{2}}l}\\ 0&{\frac{{{k^{4}}}}{4}}&0&{-\frac{{i{k^{3}}l}}{2}}\\ {\frac{{i{k^{3}}l}}{2}}&0&{\frac{{{k^{4}}+4{k^{2}}{l^{2}}}}{4}}&{-2k{l^{2}}}\\ {-i{k^{2}}l}&{\frac{{i{k^{3}}l}}{2}}&{-2k{l^{2}}}&{\frac{{{k^{4}}+4\left({4+{k^{2}}}\right){l^{2}}}}{4}}\end{array}}\right), (19)

where k=F/gk=F/g, l=J/gl=J/g and ii is an imaginary unit. As a result, the analytical solutions for energy and conversion efficiency in the resonance case can be obtained by combining Eqs. (6)-(10) and (19).

Refer to caption
Refer to caption
Figure 2: The extractable work, efficiency and tightness as a function of J/gJ/g or F/gF/g for the steady state, where exergy Wf(){{W_{f}\left(\infty\right)}} (the dashed red line), ergotropy We(){{W_{e}\left(\infty\right)}} (the dot-dashed blue line), efficiency η1(){\eta_{1}}\left(\infty\right) (the red rectangle), efficiency η2(){\eta_{2}}\left(\infty\right) (the blue triangle) and bound’s tightness Δxz()\Delta^{xz}(\infty) (black solid circle). Graphs (a): F=2gF=2g. Graphs (b): J=2gJ=2g. And the thin dashed gray line represents J/g=F/gJ/g=F/g. So J/g<F/gJ/g<F/g refers to the left area of Graphs (a) and the right area of Graphs (b); J/g>F/gJ/g>F/g refers to the right area of Graphs (a) and the left area of Graphs (b). All is plotted with β=100/ω0\beta=100/{\omega_{0}} and T=ΔA=ΔB=0T={\Delta_{A}}={\Delta_{B}}=0.

Fig. 2(a) has plotted the extractable work and efficiency of battery with different J/gJ/g. Following the figure, one can obtain: the greater J/gJ/g is in favor of ergotropy and η2(){\eta_{2}}\left(\infty\right), but against exergy. That is, more ergotropy can be obtained while reducing the cost of charging. In addition, η1(){\eta_{1}}\left(\infty\right) has been stable at a high level without being disturbed by J/gJ/g, i.e., Wf()ΔEB(){{W_{f}\left(\infty\right)}}\approx{{\Delta{E_{B}}\left(\infty\right)}}, which suggests exergy represents the variation of stored energy, and almost all of the stored energy can be extracted into exergy. Specifically, it can be seen that strong tightness (i.e., the decrease of Δxz(){\Delta^{xz}}\left(\infty\right)) can enhance exergy (for J/g<F/gJ/g<F/g) or enhance ergotropy (for J/g>F/gJ/g>F/g) to a certain extent.

Fig. 2(b) examines the extractable work and efficiency for different F/gF/g. It is found that, the extractable work may increase (e.g., exergy) or decreases (e.g., ergotropy) with the increase of F/gF/g, where both the efficiency is decreasing, which means a stronger drive will increase the cost of charging. In addition, the strong tightness can increase the extractable work within a certain range.

More importantly, both of Figs. 2(a) and 2(b) indicate that strong tightness always corresponds to higher energy conversion efficiency when J/g>F/gJ/g>F/g. To prove generality, this result is fully validated in Fig. 7 of the Appendix A, which proves the result is hold for large-scale range of J/g[0,1000]J/g\in[0,1000]. And the strong tightness stems from the fact that the uncertainty continues to decline (i.e., Ulxz()U_{l}^{xz}\left(\infty\right) is decreasing) when J/g>F/gJ/g>F/g (see Fig. 7).

In short, when the resonant drive and zero-temperature reservoir co-exist, it is concluded that: (i) the larger dissipation can enhance the extractable work at a lower energy cost (e.g., ergotropy), while a stronger drive cannot; (ii) The variation in tightness always indicates the energy conversion efficiency with respect to J/g>F/gJ/g>F/g.

Refer to caption
Refer to caption
Figure 3: The extractable work, efficiency and tightness as a function of Δ/g\Delta/g. In addition, we also consider the relative intensity of dissipation and drive, where F=2J=6gF=2J=6g for Graphs (a) and J=2F=4gJ=2F=4g for Graphs (b). And the thin dashed gray line represents Δ=g\Delta=g. For all plots, β=100/ω0\beta=100/{\omega_{0}}, T=0T=0 and ΔA=ΔB=Δ{\Delta_{A}}={\Delta_{B}}=\Delta are set. And all the markers are same as those in Fig. 2.

III.2 Detuning

In this subsection, we further investigate the effect of detuning drive on extractable work and the relationship between tightness and energy conversion efficiency. For simplicity, we assume detuning ΔA=ΔB=Δ\Delta_{A}=\Delta_{B}=\Delta. Therefore, we can attain the state of the charging system ρAB()|T=0{\rho_{AB}}\left(\infty\right)\left|{{}_{T=0}}\right., which solution is a 4x4 matrix given in the Appendix B.

As shown in Fig. 3(a), increase of the detuning reduces the extractable work and efficiency η2(){\eta_{2}}\left(\infty\right) for J/g<F/gJ/g<F/g. Moreover, the weak tightness (i.e., the increase of Δxz(){\Delta^{xz}}\left(\infty\right)) promotes the accumulation of extractable work. Since the efficiency η1(){\eta_{1}}\left(\infty\right) is nearly invariable regarding Δ/g\Delta/g, we will focus on investigating η2(){\eta_{2}}\left(\infty\right) in the following.

Fig. 3(b) illustrates the growth of detuning still inhibits extractable work for J/g>F/gJ/g>F/g. While efficiency may be improved, it comes at the cost of losing extractable work. Combining with Fig. 3(a), detuning is always detrimental to charging of QB. Compared with Fig. 2, it finds that the resonant driving field is beneficial to boost the system¡¯s energy,which is a key factor in the course of battery energy storage and extraction. In addition, when Δ/g<1\Delta/g<1, the stronger tightness indicates more extractable work, while the stronger tightness shows less extractable work in the case of Δ/g>1\Delta/g>1.

Furthermore, regardless of the relative strength of dissipation and drive, weak tightness frequently represents a higher energy conversion efficiency. Which suggests that tightness has great potential to indicates the efficiency in detuned charging protocol.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 4: For the bosonic reservoir, the extractable work, efficiency and tightness as a function of nbn_{b}. Graphs (a) and (c): F=2J=4gF=2J=4g. Graphs (b) and (d): J=2F=4gJ=2F=4g. And ΔA=ΔB=0{\Delta_{A}}={\Delta_{B}}=0 for (a) and (b); ΔA=ΔB=0.1g{\Delta_{A}}={\Delta_{B}}=0.1g for (c) and (d). For all plots, β=100/ω0\beta=100/{\omega_{0}} are set. And all the markers are same as those in Fig. 2.

III.3 Bosonic and fermionic reservoirs

Next, we turn to consider the effects of different reservoirs on extractable work and conversion efficiency. Specifically, the temperature of the reservoirs is in charge of the average particle number. To be explicit, a higher temperature generally leads to an increase in the average particle number. As a matter of fact, it is largely relevant to battery performance. Herein, we examine the connection between the tightness and energy of the QB system of interest in nonzero-temperature bosonic and fermionic reservoirs. Regarding bosonic and fermionic reservoirs, the charging system can be depicted as ρAB()|ΔA=ΔB=0{\rho_{AB}}\left(\infty\right)\left|{{}_{{\Delta_{A}}={\Delta_{B}}=0}}\right. for resonance-driven and non-zero temperature reservoirs and the solutions have been provided in the Appendix B.

Figs. 4(a) and 4(b) show the effect of increasing average particle number nb{n_{b}} on QB in the bosonic reservoir for resonance drive. It can be seen that regardless of the relative magnitude of dissipation and drive, an initial increase in nb{n_{b}} may cause a transient mutation in extractable work or efficiency (e.g., nb[0,0.4]{n_{b}}\in\left[{0,0.4}\right] for 4(a), nb[0,0.2]{n_{b}}\in\left[{0,0.2}\right] for 4(b)); However, for a greater range of nb{n_{b}}, the increasing nb{n_{b}} may result in decreased ergotropy, efficiency and improved exergy. This manifests that the different average particle number nb{n_{b}} has versatile influence on extractable work and efficiency of QBs.

Generally speaking, the high-performance QB requires high efficiency while retaining as much extractable work as possible, and detuning is destructive to energy, as mentioned before. Therefore, it is necessary to consider a small detuning case while we observe the non-resonance scenario. Figs. 4(c)-(d) show the influence of nb{n_{b}} on work and efficiency in the presence of a detuning drive. It can be found that the detuning cannot change the variation trend of extractable work and efficiency compared to the case of resonance. Combining with Figs. 4(a)-(d), we have that the non-zero temperature bosonic reservoir (nb0n_{b}\neq 0) may benefit exergy while inhibit ergotropy and efficiency compared with zero temperature (nb=0n_{b}=0).

In the meantime, Δxz(){\Delta^{xz}}\left(\infty\right) decreases with the increase of nb{n_{b}}, which suggest weak tightness often improved efficiency in the non-zero temperature bosonic reservoir.

Figs. 5(a)-(d) show the effect of the average particle number nf{n_{f}} of the fermionic reservoir on extractable work and efficiency in the presence of resonant or detuning drive. The figures shows that non-zero temperature reservoir (i.e., nf0{n_{f}}\neq 0) have the potential to improve extractable work and efficiency compared to zero temperature reservoir (i.e., nf=0{n_{f}}=0), since the work and efficiency have larger peaks at non-zero temperature. In addition, the variation of nf{n_{f}} leads to a phase transition of battery’s energy, and we clearly reveal the existence of phase transition point (i.e., nf=0.5n_{f}=0.5) in the Appendix C (see Fig. 8).

In brief, in the non-zero temperature fermionic reservoir, we obtain: (i) weak tightness is accompanied by improvements in ergotropy and efficiency; Raising exergy requires either strong tightness (before the phase transition) or weak tightness (after the phase transition). (ii) Tightness can predict the energy phase transition of the QB, e.g., the optimal tightness (Δxz()=0{\Delta^{xz}}\left(\infty\right)=0) corresponds to the phase transition point in Fig. 5.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 5: For the fermionic reservoir, all settings are same as those in Fig. 4.

IV Contribution of EUR to the conversion efficiency of QB

It was found that the tightness of the entropic uncertainty relation was closely associated with the battery’s conversion efficiency. However, the tightness is dependent on the uncertainty and its lower bound. Thus, the mechanism by which the EUR affects the efficiency of QB deserves further investigation. Motivated by this, we investigated the dynamics of the lower bound and the uncertainty.

Figs. 6(a)-(d) display the dynamics of uncertainty and lower bound. It can be seen that in the case of different parameters regimes, smaller uncertainty and lower bound often represent more efficient energy conversion. In particular, higher efficiency is accompanied by the fact that the distance of the uncertainty from the lower bound (i.e., Δxz(){\Delta^{xz}}\left(\infty\right)) is always getting smaller (for Fig. 6(a) when J/g>F/gJ/g>F/g) or larger (for Figs. 6(b)-(d)). Therefore, the tightness can be an excellent indicator of efficiency.

Owing to the contribution of tightness to efficiency, it concludes that the entropic uncertainty relation is strongly related to the conversion efficiency of extractable work in QBs. The current result provides a new reliable indicator that predicts energy and efficiency of QBs.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 6: Uncertainty Ulxz()U_{l}^{xz}\left(\infty\right) (red diamond), lower bound Urxz()U_{r}^{xz}\left(\infty\right) (orange pentagon), efficiency η2(){\eta_{2}}\left(\infty\right) (blue triangle) as a function of different parameters. Graph (a): J=2gJ=2g and T=ΔA=ΔB=0T={\Delta_{A}}={\Delta_{B}}=0; Graph (b): F=2J=6gF=2J=6g, T=0T=0 and ΔA=ΔB=Δ{\Delta_{A}}={\Delta_{B}}=\Delta; Graphs (c) and (d): J=2F=4gJ=2F=4g and ΔA=ΔB=0{\Delta_{A}}={\Delta_{B}}=0. All is plotted for β=100/ω0\beta=100/{\omega_{0}}.

V Discussion and conclusion

In summary, we investigated the extractable work and the conversion efficiency of a QB model with a battery-charger-field in the presence of bosonic or fermionic reservoir using entropic uncertainty relations. We discussed in detail the effect of resonance and detuning between the charger-battery and driving field on the charging process of the QBs and examined how the different reservoirs coupling to the charger affect the extractable work and the efficiency. Several interesting results can be outlined: (1) resonance drive and higher-temperature reservoirs (both bosonic and fermionic) often induce battery to store more exergy. By contrast, lower-temperature bosonic reservoir favors ergotropy, and more ergotropy leads to more efficiency in most cases. Particularly, stronger dissipation or weaker drive will also improve efficiency. (2) On the one hand, in non-zero temperature reservoirs, strong tightness may increase exergy while suppressing ergotropy. In addition, energy variation induced by dissipation, drive strength and detuning complicates the relationship between tightness and extractable work. On the other hand, the tightness maintains a unified relationship with conversion efficiency compared with extractable work. The stronger tightness always indicates higher efficiency for dissipation greater than drive, when resonant drive and zero-temperature reservoir coexist. Meanwhile the weak tightness is often accompanied by higher efficiency when they do not coexist. (3) In various parametric regimes, smaller uncertainty and its bound often indicate higher efficiency. Therefore, the current investigations reveal that EURs can be applied to predict the QB’s extractable work and its conversion efficiency, which is of fundamental importance in the pursuit of high-performance practical QBs.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 7: The efficiency η1(){\eta_{1}}\left(\infty\right) (solid red line), η2(){\eta_{2}}\left(\infty\right) (dashed blue line), uncertainty Ulxz()U_{l}^{xz}\left(\infty\right) (dotted magenta line) and tightness Δxz(){\Delta^{xz}}\left(\infty\right) (dot-dashed black line) as a function of F/gF/g where Graph (a): J=4gJ=4g; Graph (b): J=10gJ=10g; Graph (c): J=100gJ=100g and Graph (d): J=1000gJ=1000g. In addition, the parameters are the same as in Fig. 2.
Acknowledgements.
This study was supported by the National Natural Science Foundation of China (Grant Nos. 12075001, 61601002, 12004006 and 12175001). Anhui Provincial Key Research and Development Plan (Grant No. 2022b13020004), and Anhui Provincial Natural Science Foundation (Grant No. 1508085QF139).

Appendix A

Fig. 2 has shown that the strong tightness is beneficial for high efficiency of energy conversion in QBs when J/g>F/gJ/g>F/g and J/g[0,4]J/g\in\left[{0,4}\right], and we investigate their relationship over a wider range of J/g>F/gJ/g>F/g with J/g[0,1000]J/g\in\left[{0,1000}\right] in Fig. 7, which reveals the close relationship between tightness and conversion efficiency. That is, in the charging protocol where resonant drive coexist with zero-temperature reservoirs, strong tightness always means higher efficiency if the dissipation of the charger is stronger than the drive field. And the strong tightness is related to the attenuation of uncertainty.

Appendix B

For the detuning drive and zero-temperature reservoir, the stable charging system is in state of

ρAB()|=T=0(a11a14a41a44),\displaystyle{\rho_{AB}}\left(\infty\right)\left|{{}_{T=0}}\right.=\left({\begin{array}[]{*{20}{c}}{a_{11}}&\cdots&{a_{14}}\\ \vdots&\ddots&\vdots\\ {a_{41}}&\cdots&{a_{44}}\end{array}}\right), (23)

where the matrix elements aij(i,j=1,2,3,4){a_{ij}}\left({i,j=1,2,3,4}\right) are given as

a11=\displaystyle a_{11}= (k4(k4m2+k2(2+4l2+m2)(1+2m2)+2((1+2l2)2+(2+5l2)m2+m4)))/4α,\displaystyle({k^{4}}({k^{4}}{m^{2}}+{k^{2}}(2+4{l^{2}}+{m^{2}})(1+2{m^{2}})+2({(1+2{l^{2}})^{2}}+(-2+5{l^{2}}){m^{2}}+{m^{4}})))/4\alpha,\hfill (24)
a21=\displaystyle a_{21}= {(k3m(2il(1+2l2)2+2k2(2+l2)m+4(1+2l2)2m+il(4+3k210l2)m2\displaystyle\{({k^{3}}m(-2il{(1+2{l^{2}})^{2}}+2{k^{2}}(2+{l^{2}})m+4{(1+2{l^{2}})^{2}}m+il(4+3{k^{2}}-10{l^{2}}){m^{2}}\hfill
+2(4+k2+10l2)m32ilm4+4m5))}/2α,\displaystyle+2(-4+{k^{2}}+10{l^{2}}){m^{3}}-2il{m^{4}}+4{m^{5}}))\}/2\alpha,\hfill (25)
a31=\displaystyle{a_{31}}= {(k3(2i(1+k2+2l2)(l+2l3)2(1+2l2)(2+k2+4l2)m+il(4+k4+10l2+k2(5+8l2))m2\displaystyle\{({k^{3}}(2i(1+{k^{2}}+2{l^{2}})(l+2{l^{3}})-2(1+2{l^{2}})(2+{k^{2}}+4{l^{2}})m+il(-4+{k^{4}}+10{l^{2}}+{k^{2}}(5+8{l^{2}})){m^{2}}\hfill
(8+k4+20l2+k2(2+8l2))m3+2i(1+k2)lm42(2+k2)m5))}/2α,\displaystyle-(-8+{k^{4}}+20{l^{2}}+{k^{2}}(2+8{l^{2}})){m^{3}}+2i(1+{k^{2}})l{m^{4}}-2(2+{k^{2}}){m^{5}}))\}/2\alpha,\hfill (26)
a41=\displaystyle{a_{41}}= {(k2(4i(1+k2+2l2)(l+2l3)+2(1+2l2)(4+3k2+10l2+4l4)m\displaystyle\{({k^{2}}(-4i(1+{k^{2}}+2{l^{2}})(l+2{l^{3}})+2(1+2{l^{2}})(4+3{k^{2}}+10{l^{2}}+4{l^{4}})m\hfill
+2il(10+3k2k4+2(72k2)l2+24l4)m2+3(k4+4k2l24(2+l4))m3\displaystyle+2il(10+3{k^{2}}-{k^{4}}+2(7-2{k^{2}}){l^{2}}+24{l^{4}}){m^{2}}+3({k^{4}}+4{k^{2}}{l^{2}}-4(2+{l^{4}})){m^{3}}\hfill
+2il(14+5k2+30l2)m46(4+k2+6l2)m5+12ilm68m7))}/2α,\displaystyle+2il(-14+5{k^{2}}+30{l^{2}}){m^{4}}-6(-4+{k^{2}}+6{l^{2}}){m^{5}}+12il{m^{6}}-8{m^{7}}))\}/2\alpha,\hfill (27)
a12=\displaystyle{a_{12}}= (k3m(2il(1+2l2)2+2k2(2+l2)m+4(1+2l2)2m+il(43k2+10l2)m2\displaystyle({k^{3}}m(2il{(1+2{l^{2}})^{2}}+2{k^{2}}(2+{l^{2}})m+4{(1+2{l^{2}})^{2}}m+il(-4-3{k^{2}}+10{l^{2}}){m^{2}}\hfill
+2(4+k2+10l2)m3+2ilm4+4m5))/2α,\displaystyle+2(-4+{k^{2}}+10{l^{2}}){m^{3}}+2il{m^{4}}+4{m^{5}}))/2\alpha,\hfill (28)
a22=\displaystyle{a_{22}}= {(k2(k6m2+8m2(l2+4m2)((1+2l2)2+(2+5l2)m2+m4)+2k2((1+2l2)2+(2+9l2+8l4)m2\displaystyle\{({k^{2}}({k^{6}}{m^{2}}+8{m^{2}}({l^{2}}+4{m^{2}})({(1+2{l^{2}})^{2}}+(-2+5{l^{2}}){m^{2}}+{m^{4}})+2{k^{2}}({(1+2{l^{2}})^{2}}+(-2+9{l^{2}}+8{l^{4}}){m^{2}}\hfill
+3(3+2l2)m4+16m6)+k4(2+5m26m4+l2(4+8m2))))}/4α,\displaystyle+3(3+2{l^{2}}){m^{4}}+16{m^{6}})+{k^{4}}(2+5{m^{2}}-6{m^{4}}+{l^{2}}(4+8{m^{2}}))))\}/4\alpha,\hfill (29)
a32=\displaystyle{a_{32}}= (k2m(k4m2+2k2(1+2l27m2)(1+m2)4(l2+4m2)((1+2l2)2+(2+5l2)m2+m4)))/2α,\displaystyle({k^{2}}m({k^{4}}{m^{2}}+2{k^{2}}(1+2{l^{2}}-7{m^{2}})(1+{m^{2}})-4({l^{2}}+4{m^{2}})({(1+2{l^{2}})^{2}}+(-2+5{l^{2}}){m^{2}}+{m^{4}})))/2\alpha,\hfill (30)
a42=\displaystyle{a_{42}}= {(k(ik6(l+im)m2+8(l+2im)m(il+2m)(1+2l23ilmm2)(1+2l2+3ilmm2)(lm+i(1+m2))\displaystyle\{(k(i{k^{6}}(l+im){m^{2}}+8(l+2im)m(il+2m)(1+2{l^{2}}-3ilm-{m^{2}})(1+2{l^{2}}+3ilm-{m^{2}})(lm+i(-1+{m^{2}}))\hfill
+k4(4il3(1+2m2)4l2(m+2m3)+il(2+5m26m4)+2m(15m2+3m4))\displaystyle+{k^{4}}(4i{l^{3}}(1+2{m^{2}})-4{l^{2}}(m+2{m^{3}})+il(2+5{m^{2}}-6{m^{4}})+2m(-1-5{m^{2}}+3{m^{4}}))\hfill
+2ik2(8il4m3+l5(4+8m2)+2il2m(2+4m2+3m4)+l3(4+9m2+6m4)\displaystyle+2i{k^{2}}(8i{l^{4}}{m^{3}}+{l^{5}}(4+8{m^{2}})+2i{l^{2}}m(2+4{m^{2}}+3{m^{4}})+{l^{3}}(4+9{m^{2}}+6{m^{4}})\hfill
+l(12m2+9m4+16m6)+2i(m6m33m5+8m7))))}/2α,\displaystyle+l(1-2{m^{2}}+9{m^{4}}+16{m^{6}})+2i(m-6{m^{3}}-3{m^{5}}+8{m^{7}}))))\}/2\alpha,\hfill (31)
a13=\displaystyle{a_{13}}= {(k3(2i(1+k2+2l2)(l+2l3)2(1+2l2)(2+k2+4l2)mil(4+k4+10l2+k2(5+8l2))m2\displaystyle\{({k^{3}}(-2i(1+{k^{2}}+2{l^{2}})(l+2{l^{3}})-2(1+2{l^{2}})(2+{k^{2}}+4{l^{2}})m-il(-4+{k^{4}}+10{l^{2}}+{k^{2}}(5+8{l^{2}})){m^{2}}\hfill
(8+k4+20l2+k2(2+8l2))m32i(1+k2)lm42(2+k2)m5))}/2α,\displaystyle-(-8+{k^{4}}+20{l^{2}}+{k^{2}}(2+8{l^{2}})){m^{3}}-2i(1+{k^{2}})l{m^{4}}-2(2+{k^{2}}){m^{5}}))\}/2\alpha, (32)
a23=\displaystyle{a_{23}}= (k2m(k4m2+2k2(1+2l27m2)(1+m2)4(l2+4m2)((1+2l2)2+(2+5l2)m2+m4)))/2α,\displaystyle({k^{2}}m({k^{4}}{m^{2}}+2{k^{2}}(1+2{l^{2}}-7{m^{2}})(1+{m^{2}})-4({l^{2}}+4{m^{2}})({(1+2{l^{2}})^{2}}+(-2+5{l^{2}}){m^{2}}+{m^{4}})))/2\alpha,\hfill (33)
a33=\displaystyle{a_{33}}= {(k2(2(1+2l2)(1+k2+2l2)(k2+4l2)+(32+20k2+5k4+k6+2(56+27k2+6k4)l2\displaystyle\{({k^{2}}(2(1+2{l^{2}})(1+{k^{2}}+2{l^{2}})({k^{2}}+4{l^{2}})+(32+20{k^{2}}+5{k^{4}}+{k^{6}}+2(56+27{k^{2}}+6{k^{4}}){l^{2}}\hfill
+8(21+4k2)l4)m2+2(32+3k4+84l2+k2(9+20l2))m4+8(4+k2)m6))}/4α,\displaystyle+8(21+4{k^{2}}){l^{4}}){m^{2}}+2(-32+3{k^{4}}+84{l^{2}}+{k^{2}}(9+20{l^{2}})){m^{4}}+8(4+{k^{2}}){m^{6}}))\}/4\alpha,\hfill (34)
a43=\displaystyle{a_{43}}= {(k(8l2(1+2l2)(1+k2+2l2)2i(l+2l3)(k2+2(2+k2)l2+8l4)m2(16+52l2+68l4\displaystyle\{(k(-8{l^{2}}(1+2{l^{2}})(1+{k^{2}}+2{l^{2}})-2i(l+2{l^{3}})(-{k^{2}}+2(2+{k^{2}}){l^{2}}+8{l^{4}})m-2(16+52{l^{2}}+68{l^{4}}\hfill
16l6+k4(2+l2)+10k2(1+l2))m2il(32+20k25k4+2(56+k2)l2+168l4)m3\displaystyle-16{l^{6}}+{k^{4}}(-2+{l^{2}})+10{k^{2}}(1+{l^{2}})){m^{2}}-il(32+20{k^{2}}-5{k^{4}}+2(56+{k^{2}}){l^{2}}+168{l^{4}}){m^{3}}\hfill
2(48+4k2+5k4+4(7+k2)l284l4)m42il(32+15k2+84l2)m5\displaystyle-2(-48+4{k^{2}}+5{k^{4}}+4(7+{k^{2}}){l^{2}}-84{l^{4}}){m^{4}}-2il(-32+15{k^{2}}+84{l^{2}}){m^{5}}\hfill
+4(24+7k2+42l2)m632ilm7+32m8))}/2α,\displaystyle+4(-24+7{k^{2}}+42{l^{2}}){m^{6}}-32il{m^{7}}+32{m^{8}}))\}/2\alpha,\hfill (35)
a14=\displaystyle{a_{14}}= {(k2(4i(1+k2+2l2)(l+2l3)+2(1+2l2)(4+3k2+10l2+4l4)m+2il(k4+k2(3+4l2)\displaystyle\{({k^{2}}(4i(1+{k^{2}}+2{l^{2}})(l+2{l^{3}})+2(1+2{l^{2}})(4+3{k^{2}}+10{l^{2}}+4{l^{4}})m+2il({k^{4}}+{k^{2}}(-3+4{l^{2}})\hfill
2(5+7l2+12l4))m2+3(k4+4k2l24(2+l4))m32il(14+5k2+30l2)m4\displaystyle-2(5+7{l^{2}}+12{l^{4}})){m^{2}}+3({k^{4}}+4{k^{2}}{l^{2}}-4(2+{l^{4}})){m^{3}}-2il(-14+5{k^{2}}+30{l^{2}}){m^{4}}\hfill
6(4+k2+6l2)m512ilm68m7))}/2α,\displaystyle-6(-4+{k^{2}}+6{l^{2}}){m^{5}}-12il{m^{6}}-8{m^{7}}))\}/2\alpha,\hfill (36)
a24=\displaystyle{a_{24}}= {(k(k6(ilm)m2+8(l2im)(l+2im)m(1ilmm2)(1+2l23ilmm2)(1+2l2+3ilmm2)\displaystyle\{(k({k^{6}}(-il-m){m^{2}}+8(l-2im)(l+2im)m(1-ilm-{m^{2}})(1+2{l^{2}}-3ilm-{m^{2}})(1+2{l^{2}}+3ilm-{m^{2}})\hfill
+k4(4il3(1+2m2)4l2(m+2m3)+2m(15m2+3m4)+il(25m2+6m4))\displaystyle+{k^{4}}(-4i{l^{3}}(1+2{m^{2}})-4{l^{2}}(m+2{m^{3}})+2m(-1-5{m^{2}}+3{m^{4}})+il(-2-5{m^{2}}+6{m^{4}}))\hfill
2ik2(8il4m3+l5(4+8m2)2il2m(2+4m2+3m4)+l3(4+9m2+6m4)\displaystyle-2i{k^{2}}(-8i{l^{4}}{m^{3}}+{l^{5}}(4+8{m^{2}})-2i{l^{2}}m(2+4{m^{2}}+3{m^{4}})+{l^{3}}(4+9{m^{2}}+6{m^{4}})\hfill
+l(12m2+9m4+16m6)2i(m6m33m5+8m7))))}/2α,\displaystyle+l(1-2{m^{2}}+9{m^{4}}+16{m^{6}})-2i(m-6{m^{3}}-3{m^{5}}+8{m^{7}}))))\}/2\alpha,\hfill (37)
a34=\displaystyle{a_{34}}= {(k(8l2(1+2l2)(1+k2+2l2)+2i(l+2l3)(k2+2(2+k2)l2+8l4)m2(16+52l2+68l416l6\displaystyle\{(k(-8{l^{2}}(1+2{l^{2}})(1+{k^{2}}+2{l^{2}})+2i(l+2{l^{3}})(-{k^{2}}+2(2+{k^{2}}){l^{2}}+8{l^{4}})m-2(16+52{l^{2}}+68{l^{4}}-16{l^{6}}\hfill
+k4(2+l2)+10k2(1+l2))m2+il(32+20k25k4+2(56+k2)l2+168l4)m3\displaystyle+{k^{4}}(-2+{l^{2}})+10{k^{2}}(1+{l^{2}})){m^{2}}+il(32+20{k^{2}}-5{k^{4}}+2(56+{k^{2}}){l^{2}}+168{l^{4}}){m^{3}}\hfill
2(48+4k2+5k4+4(7+k2)l284l4)m4+2il(32+15k2+84l2)m5\displaystyle-2(-48+4{k^{2}}+5{k^{4}}+4(7+{k^{2}}){l^{2}}-84{l^{4}}){m^{4}}+2il(-32+15{k^{2}}+84{l^{2}}){m^{5}}\hfill
+4(24+7k2+42l2)m6+32ilm7+32m8))}/2α,\displaystyle+4(-24+7{k^{2}}+42{l^{2}}){m^{6}}+32il{m^{7}}+32{m^{8}}))\}/2\alpha,\hfill (38)
a44=\displaystyle{a_{44}}= {((k8m2+32(l2+4m2)(1+(2+l2)m2+m4)((1+2l2)2+(2+5l2)m2+m4)+k6(2+5m22m4\displaystyle\{(({k^{8}}{m^{2}}+32({l^{2}}+4{m^{2}})(1+(-2+{l^{2}}){m^{2}}+{m^{4}})({(1+2{l^{2}})^{2}}+(-2+5{l^{2}}){m^{2}}+{m^{4}})+{k^{6}}(2+5{m^{2}}-2{m^{4}}\hfill
+4l2(1+3m2))+2k4(16m2+49m4+4m6+12l4(1+2m2)+l2(8+31m2+10m4))\displaystyle+4{l^{2}}(1+3{m^{2}}))+2{k^{4}}(1-6{m^{2}}+49{m^{4}}+4{m^{6}}+12{l^{4}}(1+2{m^{2}})+{l^{2}}(8+31{m^{2}}+10{m^{4}}))\hfill
+8k2(4l6(1+3m2)+4(3+5m2)(mm3)2+l4(12+13m2+35m4)\displaystyle+8{k^{2}}(4{l^{6}}(1+3{m^{2}})+4(3+5{m^{2}}){(m-{m^{3}})^{2}}+{l^{4}}(12+13{m^{2}}+35{m^{4}})\hfill
+l2(5+13m2+35m4+43m6))))}/4α,\displaystyle+{l^{2}}(5+13{m^{2}}+35{m^{4}}+43{m^{6}}))))\}/4\alpha, (39)

with α=k8m2+k6(1+2l2)(2+5m2)+8(l2+4m2)(1+(2+l2)m2+m4)((1+2l2)2+(2+5l2)m2+m4)+2k4(1+17m4+6m6+4l4(2+3m2)+3l2(2+6m2+3m4))+4k2(4m2(1+m2)2(2+3m2)+l6(4+8m2)+l4(8+19m2+28m4)+l2(3+14m2+35m4+32m6)))\alpha={k^{8}}{m^{2}}+{k^{6}}(1+2{l^{2}})(2+5{m^{2}})+8({l^{2}}+4{m^{2}})(1+(-2+{l^{2}}){m^{2}}+{m^{4}})({(1+2{l^{2}})^{2}}+(-2+5{l^{2}}){m^{2}}+{m^{4}})+2{k^{4}}(1+17{m^{4}}+6{m^{6}}+4{l^{4}}(2+3{m^{2}})+3{l^{2}}(2+6{m^{2}}+3{m^{4}}))+4{k^{2}}(4{m^{2}}{(-1+{m^{2}})^{2}}(2+3{m^{2}})+{l^{6}}(4+8{m^{2}})+{l^{4}}(8+19{m^{2}}+28{m^{4}})+{l^{2}}(3+14{m^{2}}+35{m^{4}}+32{m^{6}}))), k=F/gk=F/g, l=J/gl=J/g, m=Δ/gm=\Delta/g, and ii is an imaginary unit.

In non-zero temperature reservoirs, the charging system can be expressed as

ρAB()|=ΔA=ΔB=0(b11b14b41b44),\displaystyle{\rho_{AB}}\left(\infty\right)\left|{{}_{{\Delta_{A}}={\Delta_{B}}=0}}\right.=\left({\begin{array}[]{*{20}{c}}{b_{11}}&\cdots&{b_{14}}\\ \vdots&\ddots&\vdots\\ {b_{41}}&\cdots&{b_{44}}\end{array}}\right), (43)

where, the matrix elements bij(i,j=1,2,3,4){b_{ij}}\left({i,j=1,2,3,4}\right) of ρAB()|ΔA=ΔB=0{\rho_{AB}}\left(\infty\right)\left|{{}_{{\Delta_{A}}={\Delta_{B}}=0}}\right. in bosonic and fermionic reservoirs corresponds to

b11=\displaystyle{b_{11}}= k4(1+2l2)+4l2nb(k2(5+2l2)+2nb(2+3k2+2(2+k2)l2+16l2nb(1+nb)))/4β\displaystyle{k^{4}}(1+2{l^{2}})+4{l^{2}}{n_{b}}({k^{2}}(5+2{l^{2}})+2{n_{b}}(2+3{k^{2}}+2(2+{k^{2}}){l^{2}}+16{l^{2}}{n_{b}}(1+{n_{b}})))/4\beta\hfill
=\displaystyle= k4(1+2l2)+4l2nf(k2(52k2+2l2)+2nf(2+k4+4l22k2(1+2l2)+4k2l2nf))/4γ,\displaystyle{k^{4}}(1+2{l^{2}})+4{l^{2}}{n_{f}}({k^{2}}(5-2{k^{2}}+2{l^{2}})+2{n_{f}}(2+{k^{4}}+4{l^{2}}-2{k^{2}}(1+2{l^{2}})+4{k^{2}}{l^{2}}{n_{f}}))/4\gamma,\hfill (44)
b21=\displaystyle{b_{21}}= 2kl2nb(1+2l2+4l2nb)/β\displaystyle 2k{l^{2}}{n_{b}}(-1+2{l^{2}}+4{l^{2}}{n_{b}})/-\beta\hfill
=\displaystyle= 2kl2nf(1+2nf)(12l2+4l2nf)/γ,\displaystyle 2k{l^{2}}{n_{f}}(-1+2{n_{f}})(1-2{l^{2}}+4{l^{2}}{n_{f}})/-\gamma, (45)
b31=\displaystyle{b_{31}}= ikl(k2(1+2l2)+16l2nb(1+nb))/2β\displaystyle ikl({k^{2}}(1+2{l^{2}})+16{l^{2}}{n_{b}}(1+{n_{b}}))/2\beta\hfill
=\displaystyle= ikl(1+2nf)(k2(1+2l2)+8(2+k2)l2(1+nf)nf)/2γ,\displaystyle ikl(-1+2{n_{f}})({k^{2}}(1+2{l^{2}})+8(-2+{k^{2}}){l^{2}}(-1+{n_{f}}){n_{f}})/-2\gamma,\hfill (46)
b41=\displaystyle{b_{41}}= ik2l(1+2l2)/β\displaystyle i{k^{2}}l(1+2{l^{2}})/-\beta\hfill
=\displaystyle= ik2l(1+2l2)/γ,\displaystyle i{k^{2}}l(1+2{l^{2}})/\gamma,\hfill (47)
b12=\displaystyle{b_{12}}= 2kl2nb(1+2l2+4l2nb)/β\displaystyle 2k{l^{2}}{n_{b}}(-1+2{l^{2}}+4{l^{2}}{n_{b}})/-\beta\hfill
=\displaystyle= 2kl2nf(1+2nf)(12l2+4l2nf)/γ,\displaystyle 2k{l^{2}}{n_{f}}(-1+2{n_{f}})(1-2{l^{2}}+4{l^{2}}{n_{f}})/-\gamma,\hfill (48)
b22=\displaystyle{b_{22}}= k4(1+2l2)+4l2nb(4+5k2+2(4+k2)l2+2nb(2+3k2+2(10+k2)l2+16l2nb(2+nb)))/4β\displaystyle{k^{4}}(1+2{l^{2}})+4{l^{2}}{n_{b}}(4+5{k^{2}}+2(4+{k^{2}}){l^{2}}+2{n_{b}}(2+3{k^{2}}+2(10+{k^{2}}){l^{2}}+16{l^{2}}{n_{b}}(2+{n_{b}})))/4\beta\hfill
=\displaystyle= k4(1+2l2)+4l2nf(4+5k22k4+2(4+k2)l2+2nf(22k2+k44(1+k2)l2+4k2l2nf))/4γ,\displaystyle{k^{4}}(1+2{l^{2}})+4{l^{2}}{n_{f}}(4+5{k^{2}}-2{k^{4}}+2(4+{k^{2}}){l^{2}}+2{n_{f}}(-2-2{k^{2}}+{k^{4}}-4(1+{k^{2}}){l^{2}}+4{k^{2}}{l^{2}}{n_{f}}))/4\gamma,\hfill (49)
b32=\displaystyle{b_{32}}= 0,\displaystyle 0,\hfill (50)
b42=\displaystyle{b_{42}}= ikl(k2(1+2l2)+16l2nb(1+nb))/2β\displaystyle ikl({k^{2}}(1+2{l^{2}})+16{l^{2}}{n_{b}}(1+{n_{b}}))/2\beta\hfill
=\displaystyle= ikl(1+2nf)(k2(1+2l2)+8(2+k2)l2(1+nf)nf)/2γ,\displaystyle ikl(-1+2{n_{f}})({k^{2}}(1+2{l^{2}})+8(-2+{k^{2}}){l^{2}}(-1+{n_{f}}){n_{f}})/-2\gamma,\hfill (51)
b13=\displaystyle{b_{13}}= ikl(k2(1+2l2)+16l2nb(1+nb))/2β\displaystyle ikl({k^{2}}(1+2{l^{2}})+16{l^{2}}{n_{b}}(1+{n_{b}}))/-2\beta\hfill
=\displaystyle= ikl(1+2nf)(k2(1+2l2)+8(2+k2)l2(1+nf)nf)/2γ,\displaystyle ikl(-1+2{n_{f}})({k^{2}}(1+2{l^{2}})+8(-2+{k^{2}}){l^{2}}(-1+{n_{f}}){n_{f}})/2\gamma,\hfill (52)
b23=\displaystyle{b_{23}}= 0,\displaystyle 0,\hfill (53)
b33=\displaystyle{b_{33}}= k2(1+2l2)(k2+4l2)+4l2nb(4+8l2+k2(7+6l2)+2nb(2+3k2+2(10+k2)l2+16l2nb(2+nb)))/4β\displaystyle{k^{2}}(1+2{l^{2}})({k^{2}}+4{l^{2}})+4{l^{2}}{n_{b}}(4+8{l^{2}}+{k^{2}}(7+6{l^{2}})+2{n_{b}}(2+3{k^{2}}+2(10+{k^{2}}){l^{2}}+16{l^{2}}{n_{b}}(2+{n_{b}})))/4\beta\hfill
=\displaystyle= k2(1+2l2)(k2+4l2)+4l2nf(42k4+8l2+k2(310l2)+2nf(2+k44l2+k2(2+8l2)4k2l2nf))/4γ,\displaystyle{k^{2}}(1+2{l^{2}})({k^{2}}+4{l^{2}})+4{l^{2}}{n_{f}}(4-2{k^{4}}+8{l^{2}}+{k^{2}}(3-10{l^{2}})+2{n_{f}}(-2+{k^{4}}-4{l^{2}}+{k^{2}}(-2+8{l^{2}})-4{k^{2}}{l^{2}}{n_{f}}))/4\gamma,\hfill (54)
b43=\displaystyle{b_{43}}= 2kl2(1+nb)(1+2l2+4l2nb)/β\displaystyle 2k{l^{2}}(1+{n_{b}})(1+2{l^{2}}+4{l^{2}}{n_{b}})/-\beta\hfill
=\displaystyle= 2kl2(1+nf)(1+2nf)(12l2+4l2nf)/γ,\displaystyle 2k{l^{2}}(-1+{n_{f}})(-1+2{n_{f}})(-1-2{l^{2}}+4{l^{2}}{n_{f}})/\gamma,\hfill (55)
b14=\displaystyle{b_{14}}= ik2l(1+2l2)/β\displaystyle i{k^{2}}l(1+2{l^{2}})/\beta\hfill
=\displaystyle= ik2l(1+2nf)(1+2l2+8l2(1+nf)nf)/γ,\displaystyle i{k^{2}}l(-1+2{n_{f}})(1+2{l^{2}}+8{l^{2}}(-1+{n_{f}}){n_{f}})/-\gamma,\hfill (56)
b24=\displaystyle{b_{24}}= ikl(k2(1+2l2)+16l2nb(1+nb))/2β\displaystyle ikl({k^{2}}(1+2{l^{2}})+16{l^{2}}{n_{b}}(1+{n_{b}}))/-2\beta\hfill
=\displaystyle= ikl(1+2nf)(k2(1+2l2)+8(2+k2)l2(1+nf)nf)/2γ,\displaystyle ikl(-1+2{n_{f}})({k^{2}}(1+2{l^{2}})+8(-2+{k^{2}}){l^{2}}(-1+{n_{f}}){n_{f}})/2\gamma,\hfill (57)
b34=\displaystyle{b_{34}}= 2kl2(1+nb)(1+2l2+4l2nb)/β\displaystyle 2k{l^{2}}(1+{n_{b}})(1+2{l^{2}}+4{l^{2}}{n_{b}})/-\beta\hfill
=\displaystyle= 2kl2(1+nf)(1+2nf)(12l2+4l2nf)/γ,\displaystyle 2k{l^{2}}(-1+{n_{f}})(-1+2{n_{f}})(-1-2{l^{2}}+4{l^{2}}{n_{f}})/\gamma,\hfill (58)
b44=\displaystyle{b_{44}}= (1+2l2)(k4+4(4+k2)l2)+4l2nb(8+7k2+6(8+k2)l2+2nb(2+3k2+2(26+k2)l2+16l2nb(3+nb)))/4β\displaystyle(1+2{l^{2}})({k^{4}}+4(4+{k^{2}}){l^{2}})+4{l^{2}}{n_{b}}(8+7{k^{2}}+6(8+{k^{2}}){l^{2}}+2{n_{b}}(2+3{k^{2}}+2(26+{k^{2}}){l^{2}}+16{l^{2}}{n_{b}}(3+{n_{b}})))/4\beta\hfill
=\displaystyle= (1+2l2)(k4+4(4+k2)l2)+4l2nf(8+3k22k42(8+5k2)l2+2nf(2+k4+4l2+k2(2+8l2)4k2l2nf))/4γ,\displaystyle(1+2{l^{2}})({k^{4}}+4(4+{k^{2}}){l^{2}})+4{l^{2}}{n_{f}}(-8+3{k^{2}}-2{k^{4}}-2(8+5{k^{2}}){l^{2}}+2{n_{f}}(2+{k^{4}}+4{l^{2}}+{k^{2}}(-2+8{l^{2}})-4{k^{2}}{l^{2}}{n_{f}}))/4\gamma, (59)

with β=(1+2l2)(k4+2(2+k2)l2)+8l2nb(1+nb)(2+3k2+2(4+k2)l2+16l2nb(1+nb))\beta=(1+2{l^{2}})({k^{4}}+2(2+{k^{2}}){l^{2}})+8{l^{2}}{n_{b}}(1+{n_{b}})(2+3{k^{2}}+2(4+{k^{2}}){l^{2}}+16{l^{2}}{n_{b}}(1+{n_{b}})), γ=(1+2l2)(k4+2(2+k2)l2)+8k2l2(2+k2+2l2)(1+nf)nf\gamma=(1+2{l^{2}})({k^{4}}+2(2+{k^{2}}){l^{2}})+8{k^{2}}{l^{2}}(-2+{k^{2}}+2{l^{2}})(-1+{n_{f}}){n_{f}}, and nb{n_{b}} (nf{n_{f}}) being the average particle number in the bosonic (fermionic) reservoir.

Refer to caption
Refer to caption
Figure 8: The first-order partial derivatives of extractable work and efficiency as a function of nfn_{f}, which includes Wf/nf\partial{W_{f}}/\partial{n_{f}} (dashed red line), We/nf\partial{W_{e}}/\partial{n_{f}} (dot-dashed blue line), η1/nf\partial{\eta_{1}}/\partial{n_{f}} (red rectangle) and η2/nf\partial{\eta_{2}}/\partial{n_{f}} (blue triangle). Graphs (a): F=2J=4gF=2J=4g; Graphs (b): J=2F=4gJ=2F=4g. For all plots, ΔA=ΔB=0{\Delta_{A}}={\Delta_{B}}=0 and β=100/ω0\beta=100/{\omega_{0}} are set.

Appendix C

Recent studies phase1 ; Dou-F-Q2 ; cells2 have shown that energy phase transitions play an important role in different charging schemes, and the energy phase transitions that may occur in batteries lead to drastic changes in energy. Quantum phase transition can even lead to adiabatic excitation that increases the energy of the system phase2 . Therefore, it is necessary to investigate the possible phase transition points in the charging protocol. Shown as Fig. 5, in a fermionic reservoir, the changing average particle number results in a continuous but non-differentiable singularity in the battery’s energy. Fig. 8 has depicted the existence of this singularity as a phase transition point (for We(){W_{e}}\left(\infty\right) and η2(){\eta_{2}}\left(\infty\right)) by calculating the first-order partial derivatives of extractable work and efficiency.

Figs. 8(a)-(b) directly show the presence of phase transition point (i.e., nf=0.5n_{f}=0.5). The partial derivative reaches its peak after the phase transition, i.e., the extractable work and efficiency surge after the phase transition. In combination with Fig. 5, it can be seen that the extractable work and efficiency can be rapidly improved after the phase transition and exceed the peak before the phase transition, which supports the positive role of phase transition in the charging protocol.

References