a practical figure of merit to compare fuel and battery materials is the specific energy , or energy to mass ratio  @xcite . as a material constant
, it measures the energy that will be harvested by using a unit mass of the material . about a decade ago
, a highly non - traditional fuel , called `` phaseonium '' , which is a three level atom with two lower states in a quantum coherent superposition , was proposed to be used in a photonic carnot engine ( pce )   @xcite .
phaseonium engine could work with a single heat bath and a phaseonium reservoir  @xcite .
this proposal stimulated much interest to quantum heat engines  @xcite .
it was later argued that existing resonator systems can not implement such an engine , due to high cavity losses and atomic dephasing @xcite .
here , we address a fundamental question of how the specific energy of phaseonium fuel is scaled with the number of quantum coherent levels .
a favorable scaling law against decoherence and dephasing could bring the phaseonium engine closer to available practical systems .
we describe multilevel generalization of phaseonium fuel in fig .
the block diagonal density matrix @xmath2 of an @xmath0 level atom is shown in fig .
[ fig1](a ) .
the excited level , denoted by ``@xmath3 '' , and the lower levels , denoted by `` @xmath4 '' are well separated from each other by an energy @xmath5 measured from the central lower level @xmath6 as shown in fig .
[ fig1](b ) .
the lower levels can be degenerate or non - degenerate .
the diagonal elements @xmath7 , and @xmath8 , with @xmath9 , determine the level populations , while the off diagonal elements @xmath10 , with @xmath11 , indicate the coherence between the levels .
coherence can be characterized by the magnitude and phase of the complex number @xmath10 .
though both the amplitude and the phase of coherent superposition states can be controlled in experiments @xcite , the main control variable for the photonic carnot engine is the phase of the coherence as the amplitude is required to be small enough to keep the system only slightly out of thermal equilibrium .
the complete graphs in fig .
[ fig1](a ) have @xmath1 nodes and @xmath12 links , representing the atomic energy levels and the coherences between them , respectively .
the simplest graph has @xmath13 nodes , which is the case of the original phaseonium proposal @xcite .
the interplay between quantum coherence and energy discussed in photon carnot engine @xcite revealed that the energy content of the phaseonium with @xmath13 can be optimized at a certain phase of the coherence .
we could envision as if we are considering more complex , larger , phaseonium molecules with the graphs having @xmath14 , corresponding to @xmath0 level atom phaseonium ( nlap ) .
( color online ) @xmath0 level atom phaseonium ( nlap ) fuel .
a ) density matrix @xmath2 and complete graph representations of nlap .
@xmath2 is @xmath0 dimensional square matrix .
its coherent block can be represented by a complete graph with @xmath1 nodes and @xmath15 links .
graphs were shown up to @xmath16 number of nodes .
b ) nlap for non - degenerate and degenerate atoms .
the excited state is denoted by @xmath3 and the lower levels are denoted by @xmath17 with @xmath18 .
the upper level is well separated from the lower levels by an energy @xmath19 measured from the central lower level @xmath6.,width=326 ]    we can imagine different phaseonium molecules can be possible for a given atom of unit mass and explore how the specific energy of the atom depends on the size of the phaseonium molecule characterized by @xmath1 .
next to the phase of coherence , @xmath1 becomes another control parameter which could favourably contribute to the enhancement of the specific energy of the single atom quantum fuel . for @xmath20 1 , the number of coherences would scale quadratically , @xmath21 .
if the quadratic coherence scaling could be translated into the energy content of the atomic fuel , we could overcome the cavity losses for implementation and boost the performance of quantum carnot engine for applications .
from fundamental point of view , such scaling analysis could reveal profound difference of quantum fuel from a classical resource as such a scaling can not exist without quantum coherence .
complete graphs of phaseonium molecules serve more than a simple counting of coherences .
they emphasize the generality of our question we address in the present contribution .
can we beat decoherence with the scaling advantage of quantum coherent resources ?
( color online ) photonic carnot engine with @xmath0 level atom phaseonium ( nlap ) fuel .
photon gas in a high quality cavity of frequency @xmath5 is the working substance and the mirrors of the cavity play the role of the piston .
nlap leaves the hohlraum at temperature @xmath22 and is subsequently prepared in a state with quantum coherence among its lower levels characterized by @xmath12 phase parameters @xmath23 with @xmath24 .
created nlaps are repeatedly injected into the cavity at a rate @xmath25 in the quantum isothermal expansion process , where heat @xmath26 is transferred to the cavity . the cycle continues with quantum adiabatic expansion and quantum isothermal compression and is completed with a quantum adiabatic compression .
an amount of heat @xmath27 is rejected into the entropy sink in the isothermal compression.,width=307 ]    we have recently numerically examined a superradiant quantum otto engine  @xcite which exhibits a similar scaling law of work with the number of atoms ; though the efficiency is independent of the temperature and the number of atoms .
our present results are fully analytical , and for a carnot cycle which can harvest work from a single heat bath and the scaling laws with the number of coherences are both for the work and the efficiency .
furthermore , specific energy of a superradiant cluster increase linearly with the number of atoms , while here for the single atom phaseonium fuel it is quadratic . quantum advantage in the charging power of quantum batteries with the number of qubits has been examined very recently  @xcite .
linear scaling of work with the number of qubits is reported ; while due to a quantum speed up of the operation time , the charging power scales quadratically  @xcite . the preparation of phaseonium fuel and using it in pce can be compared to charging and discharging a single qubit quantum battery with multiple quantum coherences .
phaseonium fuel or battery allows for quadratic scaling in harvested work , efficiency , and the specific energy with the number of quantum coherences .
we examined the charging or preparation cost of the phaseonium battery and compared it with the harvested work by pce .
our result verified that the second law is obeyed in our system .
the operation of photonic carnot engine is described in fig .
the working fluid of the engine is the photon gas in a high quality single mode cavity of frequency @xmath5 .
the radiation pressure by the cavity photons applies on the cavity mirrors playing the role of the piston of the engine .
the quantum fuel of the engine is an nlap .
the quantum carnot cycle consists of two quantum isothermal and two quantum adiabatic processes .    in the isothermal expansion ,
nlaps are generated and injected into the cavity at a rate @xmath25 .
the interaction time @xmath28 between an nlap and the cavity field is short , @xmath29 , so that only one nlap can be present in the cavity  @xcite .
coherences in nlap are characterized by @xmath12 phase parameters @xmath23 , with @xmath30 .
coherent superposition states in @xmath0 level atom system can be generated by stimulated raman adiabatic passage  @xcite , morris - shore transformation  @xcite , or quantum householder reflection techniques  @xcite .
thermalization of the single atom can be considered relatively fast and hence the injection rate would be limited by the time of coherence preparation .
the choice of specific technique of coherence induction depends on the details of a particular implementation .
if the amplitudes of the coherences are much smaller than the level populations , then nlap can be assumed in an approximate thermal equilibrium with a thermal reservoir ( hohlraum ) at a temperature @xmath22 . during the interaction , the mean number of photons , @xmath31 , and the cavity temperature increases ; while the expansion cools down the cavity when there is no atom inside .
repeated injection of nlaps into the cavity maintains the cavity field at a temperature @xmath32 by transferring a total amount of heat into the cavity as @xmath26 .
@xmath32 is an effective temperature defined in terms of the steady state photon number @xmath33 as @xmath34 , with @xmath35 is the boltzman constant
. it can be higher than @xmath22 in the presence of coherence  @xcite .
the cavity volume , and hence the frequency , change negligibly , @xmath36 .
the cycle continues with an adiabatic expansion where the entropy remains constant and the temperature drops as the @xmath5 changes appreciably . following step is the isothermal compression in which heat @xmath27 is transferred from cavity to a cold reservoir at a temperature @xmath37 .
the cycle is completed by adiabatic compression where the temperature is raised back to @xmath32 .
the net work extracted from the cycle is @xmath38 where @xmath39 and @xmath40 .
the mean photon number @xmath41 and the temperature @xmath42 at the beginning of the @xmath43 stage determine the entropy @xmath44 by @xmath45 .
using @xmath46 , @xmath47 , @xmath48 , and @xmath49 , we write @xmath50 .
the efficiency of the engine is defined as @xmath51 .
it reduces to @xmath52 .
this coincides with the standard definition of thermodynamic efficiency in carnot cycle and used in the original proposal of the phaseonium fuel  @xcite as well as in the arguments against its feasibility in the presence of decoherence channels  @xcite . in order to present results comparable to the previous works
, we calculate the efficiency as defined by these studies . to avoid any misleading impressions however , we emphasize that in practical considerations round trip efficiency can be more relevant figure of merit .
the round trip efficiency of the engine should include the cost of the preparation of the quantum coherent atom ; which would ensure the validity of the second law  @xcite . on the other hand
, it was noted that the cost of quantum fuel can be expensive  @xcite , but it is still appealing as it can be used to harvest work from a single heat reservoir .
our objective here is not to discuss if such pces can be efficient enough for certain applications , but to examine if such devices , proposed in  @xcite , can produce positive work in the presence of decoherence by exploiting a scaling advantage of multiple coherence resources , in contrast to the negative conclusions of earlier studies  @xcite .    during the adiabatic process @xmath53
does not change so that @xmath54 and @xmath55 .
these relations reveal that work and efficiency of the photonic carnot engine can be calculated by determining the @xmath33 at the end of the isothermal expansion stage .
this manuscript organised as follws : in sec .
ii , we describe the n - level phaseonium model , and pce system dynamics .
we review and discuss the analytical and numerical verification of the analytical results in sec .
we also estimate the preparation cost of nlap in this section .
we conclude the results in sec .
in order to find the @xmath33 , we solve @xmath56 where @xmath57 . here
@xmath58 is the fock number state for the cavity photons and @xmath2 is the reduced density matrix of the cavity field .
the equation of motion for @xmath2 can be obtained by tracing the equation of motion of the complete system over atomic degrees of freedom @xmath59_{nn})\label{eq.ofmot},\ ] ] where @xmath60 is the hamiltonian of the arbitrary @xmath61 atom in the interaction picture relative to the cavity photons , with @xmath62 and @xmath63 . here
@xmath64 , @xmath65 are the energies of atomic states @xmath66 and @xmath67 , with @xmath18 , @xmath68 is the coupling rate between the atom and the field , and @xmath69 is the photon annihilation operator .
the model hamiltonian describes a situation where @xmath0 level atom is coupled to a single mode cavity in a fan shaped transition scheme .
a more realistic model requires consideration of multimode cavity coupled to an atom with multiple upper and lower hyperfine levels  @xcite .
such models can be reduced to effective single mode cavity and multilevel atom interactions  @xcite or can be directly described by generalized master equations of micromasers  @xcite . atoms with fan shaped degenerate level schemes are also studied from the perspective of generating large superposition states  @xcite .
the central question for us here is the dependence of work and efficiency on the number of the superposed quantum states and we will only consider single upper level and a set of degenerate or non - degenerate lower levels for simplicity
.    analytically calculating the right hand side of the eq .
( [ eq.ofmot ] ) , we find ( see appendix for details )    @xmath70\nonumber\\ & + ( \sum_{i=1}^n k_{b_i}\rho_{b_ib_i}+\sum_{i < j } k_{ij}^{\phi_{ij}}|\rho_{b_ib_j}|)\nonumber\\ & \times [ n\rho_{nn}-(n+1)\rho_{n+1,n+1}]\},\end{aligned}\ ] ]    where the coefficients @xmath71 and @xmath72 depend on the atomic relaxation rate @xmath73 , atomic dephasing rate @xmath74 , detuning parameter @xmath75 , with @xmath76 , and @xmath77 , as well as the coherence parameters @xmath23 and @xmath78 , by the relations given in the appendix .
thus , we obtain the rate of change of average photon number @xmath79 where @xmath80 and @xmath81 .
the equation of motion for the evolution of population elements in the density matrix coincides with the thermalization dynamics of a resonator coupled to a heat bath .
accordingly , the coarse grained dynamics effectively describes sequence of nlap injected into the resonator as a mesoscopic ensemble of @xmath0 level atoms acting as a heat bath .
the off diagonal elements of the density matrix or the coherences can be kept vanishingly small to describe the steady state approximately as a thermal equilibrium state .
the corresponding effective temperature can be determined by the modified detailed balance condition to reach such a quasi equilibrium state in eq .
( [ eq.mot2 ] ) which gives @xmath82 the detailed balance between the thermal reservoir at @xmath22 and the photon gas in the resonator is broken but there is a modified detailed balance between the coherent atomic ensemble and the resonator photons .
accordingly the resonator can reach a thermal equilibrium at a different temperature @xmath32 than @xmath22 .    the steady state of the eq .
( [ rateeq . ] ) yields the average photon number as @xmath83 where @xmath84 is the average photon number in the absence of coherence . using @xmath85
, we determine the effective cavity temperature as @xmath86 by using high temperature approximations @xmath87 and @xmath88 in eq .
( [ nbarf ] ) .
therefore , the efficiency of the photonic carnot engine in the case of nlap becomes @xmath89 where @xmath90 is the carnot efficiency .
note that for @xmath91 , @xmath92 but @xmath93 could have a positive value for particular values of control parameters @xmath94 . in order to get further analytical results we will make some simplifying assumptions .
we focus on degenerate nlap case to proceed analytically , for which @xmath95 , @xmath96 , @xmath97 , @xmath98 , @xmath99 and @xmath100 .
in addition , we consider phase locked equal amplitude coherences with @xmath101 and @xmath102 . hence the coefficients in eq .
( [ eq.mot2 ] ) become @xmath103 , @xmath104 and @xmath105 , and hence eq .
( [ rateeq . ] ) reduces to @xmath106-\kappa\bar{n}_{\phi } , \label{degrateeq.}\end{aligned}\ ] ] for @xmath107 , @xmath108 , where @xmath109 , @xmath110 , @xmath111 with @xmath112 .
here we introduced @xmath113 and @xmath114 , with @xmath115 , as the decoherence rate due to the dissipation in the cavity and a phenomenological decoherence factor due to atomic dephasing , respectively  @xcite .
while the dephasing factor is phenomenologically introduced in  @xcite we provide its rigorous microscopic derivation in the appendix .
steady state solution of eq .
( [ degrateeq . ] ) yields an effective temperature given by @xmath116 in the high temperature limit where @xmath117 with @xmath118 . for small coherence and decoherence terms in @xmath119
, an approximate expression can be written for the effective temperature @xmath120 this result shows that if the reduction of the magnitude of coherence due to dephasing is slower than the quadratic increase with @xmath1 , then the multilevel coherence could be used to beat the decoherence induced by the cavity dissipation .
the magnitude of coherence @xmath121 is limited by the positivity requirement of the density matrix as well as the thermal equilibrium requirement of the cavity field .
the former condition require @xmath122 so that @xmath123 for @xmath107 as @xmath124 for @xmath107 .
accordingly one can fix the coherence amplitude @xmath121 as a constant as long as it remains smaller than @xmath125 for the range of @xmath1 values . for a realistic number of levels
this is not a very restrictive condition .
more severe limitation on @xmath121 is due to the quasi - equilibrium condition of the photon gas .
we will take @xmath126 and consider @xmath127 in our numerical examinations . in the classical asymptotical limit of @xmath128
then @xmath129 as @xmath125 and hence the quadratic scaling reduces to a linear one for which the specific energy becomes a constant as expected for classical systems .
mesoscopic systems in quantum regime are therefore necessary to exploit the quadratic scaling in the specific energy .    we note that the coarse grained dynamics is designed on purpose to determine the steady state by rapid convergence using a numerically efficient dynamical equations . analytical solution of the mean photon number dynamics for the degenerate case and with @xmath130 , @xmath131 , @xmath132 gives @xmath133 where @xmath134 is the thermalization time .
we will first discuss the predicted analytical state states by the coarse grained master equation in modern resonator settings then examine the exact numerical desription of the dynamics of the system to verify the analytical results in the subsequent subsections .
( color online ) extracted work ( @xmath135 ) and efficiency ( @xmath136 ) of photonic carnot engine , with @xmath0 level atom phaseonium ( nlap ) fuel , depending on the number of degenerate coherent ground state levels @xmath1 , for different decoherence factor models ( a)-(b ) @xmath137 , ( c)-(d ) @xmath138 .
( e)-(f ) @xmath139 , where @xmath140 .
coherence parameter is @xmath141 and the initial thermal coherent atomic temperature is @xmath142 in units of @xmath143 .
@xmath144 values are 0.15 , 0.1 , 0.05 , 0.001 for ( a)-(b ) , 0.14 , 0.12 , 0.1 , 0.08 for ( c)-(d ) and 0.012 , 0.01 , 0.008 , 0.006 for ( e)-(f ) in decreasing order for the lower to upper curves respectively .
the plots are given for the circuit qed parameters in  @xcite .
the quantities @xmath145 , @xmath146 , @xmath147 , and @xmath148 , which are the coupling coefficient to the cavity field , atomic injection rate , cavity loss term , and atomic decay respectively , are dimensionless and scaled with the resonance frequency @xmath149 ghz .
@xmath136 is dimensionless , @xmath135 is dimensionless and scaled with @xmath5.,width=326 ]      in the high temperature limit ( @xmath150 , the entropy change in the isothermal expansion stage is @xmath151 and the heat input becomes @xmath152 .
the work output at @xmath153 is found to be @xmath154 where @xmath155 , respectively . in superconducting circuit , microwave and optical resonators
, it is estimated that @xmath156  @xcite .
@xmath157 should be replaced by @xmath12 for smaller number of levels .
accordingly , by using five or more level quantum phaseonium fuel , the working fluid can beat quantum decoherence to harvest positive work .    in fig .
[ fig : fig3 ] , we plot the work output and efficiency of the photonic carnot engine with degenerate nlap fuel , depending on the number of quantum coherent levels .
we consider @xmath1 independent as well as @xmath1 dependent scaling models  @xcite for the decoherence factor and take @xmath137 in fig .
[ fig : fig3](a)-(b ) , @xmath138 in fig .
[ fig : fig3](c)-(d ) , and @xmath139 in fig .
[ fig : fig3](e)-(f ) , where @xmath140 as shown in the appendix .
the plots are given for the circuit qed parameters in  @xcite .
we consider larger atomic dephasing rates than the typical values to demonstrate its limiting effect on @xmath135 and @xmath136 .
the plots indicate that even when there is large dephasing , which can increase with @xmath1 linearly or quadratically , @xmath135 and @xmath136 can retain their quadratic power law with @xmath1 up to a critical @xmath1 .
( color online ) extracted work ( @xmath135 ) and efficiency ( @xmath136 ) of microwave resonator photonic carnot engine , with @xmath0 level atom phaseonium ( nlap ) fuel , depending on the number of degenerate coherent ground state levels @xmath1 for different decoherence factors ( a)-(b ) @xmath158 , ( c)-(d ) @xmath159 , ( e)-(f ) @xmath139 , where @xmath140 .
coherence parameter is @xmath141 and the initial thermal coherent atomic temperature is @xmath142 in units of @xmath143 .
@xmath144 values are 0.1 , 0.15 , 0.05 , 0.001 for ( a)-(b ) , 0.12 , 0.1 , 0.08 , 0.06 for ( c ) , 0.14 , 0.12 , 0.1 , 0.08 for ( d ) , 0.12 , 0.1 , 0.08 , 0.06 for ( e ) and 0.01 , 0.008 , 0.006 , 0.004 in decreasing order for the lower to upper curves respectively . the parameters @xmath160 , @xmath161 , @xmath162 and @xmath163 are the atom - field coupling coefficient , atomic injection rate , cavity loss term and atomic decay respectively .
they are dimensionless and scaled with the typical resonance frequency is @xmath164 ghz @xcite .
@xmath136 is dimensionless , @xmath135 is dimensionless and scaled with @xmath5.,width=326 ]    similar results are found for the cases of optical and microwave cavities .
we see both in fig .
[ fig : fig4](a)-(b ) and fig .
[ fig : fig5](a)-(b ) that when dephasing is independent of @xmath1 , the work output and efficiency increases quadratically with the number of coherent levels . if the dephasing rate is increasing linearly with @xmath1 as in fig .
[ fig : fig4](c)-(d ) and fig .
[ fig : fig5](c)-(d ) , or if it is increasing quadratically with @xmath1 as in fig .
[ fig : fig4](e)-(f ) and fig .
[ fig : fig5](e)-(f ) , the work output and efficiency of the photonic engine is enhanced quadratically with the number of coherent levels only up to critical @xmath1 .
beyond which the work output and efficiency decays exponentially due to the dominating effect of decoherence .
( color online ) extracted work ( @xmath135 ) and efficiency ( @xmath136 ) of optical resonator photonic carnot engine , with @xmath0 level atom phaseonium ( nlap ) fuel , depending on the number of degenerate coherent ground state levels @xmath1 for different decoherence factors ( a)-(b ) @xmath158 , ( c)-(d ) @xmath159 , ( e)-(f ) @xmath165 , where @xmath140 .
coherence parameter is @xmath141 and the initial thermal coherent atomic temperature is @xmath142 in units of @xmath143 .
@xmath144 values are 0.15 , 0.1 , 0.05 , 0.001 for ( a)-(b ) , 0.14 , 0.12 , 0.1 , 0.08 for ( c)-(d ) and 0.01 , 0.008 , 0.006 , 0.004 for ( e)-(f ) in decreasing order for the lower to upper curves respectively . the parameters @xmath166 , @xmath167 , @xmath168 and @xmath169 are the coupling frequency to the cavity field , atomic injection rate , cavity loss term and atomic decay respectively .
they are dimensionless and scaled with the typical resonance frequency is @xmath170 thz @xcite .
@xmath136 is dimensionless , @xmath135 is dimensionless and scaled with @xmath5.,width=345 ]      in order to perform a faithful numerical simulation of a typical set up described in the theory , we investigate the injection process in detail .
we assume a regular atomic injection of rydberg atoms into a fabry - perot cavity  @xcite with an atomic interaction time @xmath28 with the cavity field and an empty cavity time @xmath171 such that @xmath172 where @xmath25 is the injection rate . during the time interval @xmath28 ,
the hamiltonian is @xmath173 while for the time interval @xmath171 , @xmath174 ( @xmath175 and @xmath176 for degenerate ground state levels ) .
we choose injection time @xmath177 where @xmath113 is the cavity decay rate and @xmath178 is the number of atoms kicking the cavity field in the photon lifetime .
the time elapsed when cavity is empty @xmath171 can be related to the interaction time such that @xmath179 .
thus , we can write @xmath180 . here
, @xmath181 is a factor introduced to measure @xmath171 in terms of @xmath28 so that @xmath182 .
we solve the master equation by numerical methods and compare the results with the developed theory .
we use qutip package @xcite in python software to solve the master equation .
we perform single atomic injection in two steps .
first step is the atom - cavity field interaction ( during @xmath28 ) and the second one is the free cavity field evolution ( during @xmath171 ) .
the master equation for the first step is written under markov and born - markov approximations as @xcite @xmath183+\gamma\sum_m^{n+1 } \mathcal{l}[l_m^{\gamma}]+\frac{\gamma_{\phi}}{2}\sum_n^n \mathcal{l}[l_n^{\phi}]\end{aligned}\ ] ] where last two terms stand for pure spontaneous emission and pure dephasing @xcite , respectively . here , @xmath184=(2x\rho
x^{\dagger}-xx^{\dagger}\rho-\rho x^{\dagger}x)/2 $ ] is a liouvillian superoperator in lindblad form and @xmath185 , @xmath186 .
we include an auxiliary state @xmath187 to the atomic state space to model the decay of the excited and the lower levels .
presence of @xmath187 is not altering the initial phaseonium state .
the auxiliary state is unpopulated and at the lower level energy .
its use allows for faithful simulation of the excited state and the degenerate ground state @xmath188 ) decay equations in eq .
( [ eq15 ] ) .
this decay model is already used in the original master equation developed for the two - level phaseonium engine  @xcite .
different decay models , for example decay of excited level to the lower levels are employed for other systems such as many atom superradiant otto engine  @xcite and similar effect of beating decoherence with scaling up coherence is found .
present contribution discusses the original photo - carnot engine  @xcite as well as the objection to its feasibility due to dephasing , phenomenologically described by factor @xmath114  @xcite .
our introduction of @xmath73 is an additional decoherence channel not included in ref .
we have found that @xmath189 can be analytically expressed in terms of @xmath114 ( see appendix for details ) .
the microscopical master equation approach describes both the original photo - carnot engine coarse - grained master equations with  @xcite and without  @xcite dephasing and generalizes them to the multilevel case ; as our mesoscopic master equation , eq .  [ eq.mot2 ]
, does analytically .    in the second step , cavity decay @xmath190 is present during the time interval @xmath171 in accordance with the key assumptions of micromaser theory  @xcite and the corresponding master equation is @xmath191+\kappa\mathcal{l}[\hat{a}].\end{aligned}\ ] ]    ( color online ) time evolution of average photon number during thermalisation process of pce field while regular injection of 2lap depending on different @xmath178 parameters .
@xmath178 values are 4500 , 1500 , 150 , 100 and 50 in decreasing order for the upper to lower curves respectively .
coherence parameter is @xmath192 , initial field temperature @xmath193 and temperature of thermal coherent atoms is @xmath194 in units of @xmath143 .
the resonant field frequency is @xmath164 ghz , cavity quality factor is @xmath195 , atom cavity field interaction time is @xmath196 @xmath197s , atom decay rate is @xmath198 hz , atom dephasing rate is @xmath199 hz and atom cavity field coupling is @xmath200 khz .
time is dimensionless and scaled with @xmath5.,width=307 ]    in fig .
[ fig : fig6 ] , we present the thermalisation process of the cavity field depending on different @xmath178 values by depicting the photon number versus scaled time
. physical parameters  @xcite are given in the figure caption consistent with the rydberg atoms in a superconducting febry - perot cavity  @xcite .
for low values of @xmath178 which corresponds to large @xmath181 , we have zigzag like curves and for high values of @xmath178 , we have smoother lines .
average photon number @xmath33 converges to the theoretical value for @xmath201 .
thermalisation time is much longer than the convergence rate of coarse grained master equation and the microscopic exact method is much more costly numerically .
( color online ) comparison of time evolution of average photon number in presence of thermal and coherent atom injection with @xmath202 and @xmath203 .
the horizontal dashed and dashed - dotted lines stands for the analytical @xmath33 and @xmath53 values in absence of loss mechanisms .
solid and dotted lines stands for time evolution of @xmath33 and @xmath53 respectively in presence of dissipation channels .
@xmath204 for both subplots corresponding to @xmath205 @xmath197s .
insets magnifies the lines between @xmath206 and @xmath207 .
all the remaining parameters are the same with that of fig .
[ fig : fig6 ] .
time is dimensionless and scaled with @xmath5.,width=307 ]    in fig .
[ fig : fig7](a)-(b ) , we express the effect of the number of degenerate atomic ground state levels @xmath1 of the coherent atoms against the dephasing and the cavity loss mechanisms .
horizontal dotted and dotted - dashed lines are the analitical values of average photon numbers ( @xmath33 , @xmath53 ) of the no loss case for each @xmath1 . in fig .
[ fig : fig7](a ) when decoherence channels are open , average photon number saturates below the analytical values for @xmath202 in accordance with the argument that 2lap phaseonium can not beat decoherence  @xcite .
average photon number exceeds @xmath53 for @xmath203 in fig .
[ fig : fig7](b ) by keeping all the other parameters same .
thus , we show that the decoherence can be beaten by using higher @xmath1 . to make the effect more visible in the figures , we take larger coherence magnitude , @xmath192 .
( color online ) comparison of numerical results ( dashed line ) with developed theory ( solid line ) of effective field temperature t@xmath208 , harvested work @xmath135 and efficiency @xmath136 respectively depending on the number of degenerate ground state levels @xmath1 .
@xmath209 and corresponding @xmath210 @xmath197s .
all the remaining parameters are the same that of fig .
[ fig : fig6 ] .
t@xmath208 , @xmath135 and @xmath136 is dimensionless and scaled with @xmath5.,width=336 ]    we also compare the consistency of effective field temperature t@xmath208 , harvested work @xmath211 and efficiency @xmath212 versus @xmath1 in fig .  [
fig : fig8 ] ( a ) , ( b ) , ( c ) between developed analytical and numerical results when dissipation channels open .
we observe a good consistency between numerical and theoretical results in steady state .
typical methods to generate quantum superposition states , such as pulse area , adiabatic passage , or stirap techniques  @xcite , utilize optical pulses interacting with the atomic system to transfer an initial atomic state to a target one .
the initial and target quantum states are known and hence one can easily determine the required unitary transformation between them .
physical implementation of the required propagator is however a much more challenging problem than the calculation of the transformation matrix .
an efficient strategy to sythesize the transformation matrix is to decompose it into a product of matrices , representing interacting steps that can be implemented by using optical pulses coupled to the atom .
an arbitrary @xmath1 dimensional unitary matrix @xmath213 can be decomposed into so called @xmath1 generalized quantum householder reflection ( qhr ) matrices or @xmath214 standard qhrs and a phase gate   @xcite .
a generalized qhr is defined by @xmath215 where @xmath216 is the identity operator and the @xmath217 is the normalized column vector with dimension @xmath1 , same with the number of the pulses , and @xmath218 is an arbitrary phase factor . the decomposition of @xmath213 in terms of generalized qhrs can be written as @xmath219 for @xmath108 , eq .
( [ genreflct ] ) reduces to @xmath220 which is the standard qhr .
the interaction represented by each householder matrix can be described by a propagator which can be determined by the morris - shore transformation  @xcite .
our n+1 level atom coupled to n optical pulses in a fan shaped transition scheme , or so called n - pod model , is a generic model that is used to discuss generation of arbitrary multilevel superposition states . under the morris - shore transformation
, the lower levels of the atom are grouped into a single bright level coupled to an effective single pulse and n-1 dark levels uncoupled from the optical pulses .
the propagator is then easily determined in this morris - shore basis .
back transformation to the original basis gives the full propagator , or the generalized qhr matrix .
both the number of qhr steps and the number of pulses used in each step are in the order of @xmath1 , therefore the total number of pulses to be used to generate the target state would be in the order of @xmath157 .
this shows that the energetic cost of preparation of the target state scales with @xmath157 , same with the work and efficiency scaling in the corresponding photonic carnot engine .
the preceding discussion is applicable to the case of mixed states as well , for which the normalized vectors of generalized qhrs are defined as @xcite @xmath221 here , @xmath222 is the @xmath43 column of @xmath213 , @xmath223^t$ ] and @xmath218 is an arbitrary phase where @xmath224 . it s shown that for an n - pod system any standard qhr @xmath225 can be realized by single pulses with an rms pulse area @xmath226 @xcite .
the corresponding unitary transformation can only link the mixed states with the identical dynamical invariants . in our case
we consider initial thermal states out of a hohlarum transformed to a final state with small coherences .
the initial and final states would then possess different spectral decompositions so that they can not be unitarily connected .
a resolution to this is suggested to exploit decoherence channels such as spontaneous emission or pure dephasing , in combination with the unitary transformation  @xcite .
we will not follow this route but use an alternative , which allows for a fully unitary procedure to generate desired coherences .
as the coherences contribute additively , to exploit their scaling advantage we do not need an exact state but an approximate one would be sufficient .
accordingly , we can simply consider an approximate approach and do not specify an exact value for the coherences .
we only need to keep them small enough to ensure slightly out of thermal equilibrium final state .
we illustrate our strategy for @xmath202 case , and suggest that in principle larger nlap can be generated by straightforward extension of this technique .
unitarity of our procedure also makes the details of generation process immaterial for the cost estimation
. the cost would be the same for other unitary equivalent processes to generate same states .
initial state of the atom out of the hohlarum at @xmath22 is the thermal density matrix @xmath227 where @xmath228 inverse temperature @xmath229 and @xmath230 is the partition function , with @xmath231 being the atomic hamiltonian .
taking @xmath202 and @xmath194 , we find @xmath232 .
target density matrix @xmath233 is taken to be @xmath234 where the off - diagonal elements between degenerate ground state levels are taken real and much smaller than the diagonal elements that are equal to those of the @xmath235 .
the initial and final density matrices have distinct dynamical invariants in their spectral decompositions and hence they can not be linked by a coherent evolution .
let us assume however an approximate link such that @xmath236 .
the unitary transformation @xmath237 can be determined from the matrix that diagonalize @xmath233 and found to be @xmath238 writing @xmath239 , the fidelity between @xmath233 and @xmath240 is determined by @xmath241 and found be @xmath242 .
thus , @xmath243 can be approximately be used as the unitary operator linking the initial and the target density matrices .
@xmath244 can be synthesized by using two standard qhrs and a phase gate as @xmath245 where@xmath246 and @xmath247 is a one - dimensional phase gate . by using eq .
( [ normalisedvect ] ) one finds the normalized column vectors to be @xmath248^t$ ] , @xmath249^t$ ] .
the example for @xmath202 here illustrates the basic principles to generate arbitrary nlap .
one can use @xmath1 pulses for each @xmath214 standard qhrs and a phase gate to synthesize a unitary transformation matrix , which is the diagonalization matrix of a quasi equilibrium thermal state at @xmath22 with small coherences .
when the pulses applied to the actual atom out of hohlarum at @xmath22 , the final state will be approximately the same with the target state used to determine the properties of the pulses .
we can now estimate the energy cost @xmath250 of generating @xmath233 using @xmath251 where @xmath157 is the total number of pulses used in the qhr technique and @xmath252 is the energy of a single pulse .
@xmath252 can be determined from the pulse area @xmath253 as in  @xcite . for a square pulse of duration @xmath254 and amplitude @xmath255
we have @xmath256 , where @xmath257 is the magnitude of the dipole moment @xmath258 where @xmath259 is the vacuum permittivity and @xmath260 is the speed of light . taking the pulse area
@xmath226 we find @xmath261 the intensity of the pulse is given by @xmath262 .
the pulse energy in a beam of radius @xmath263 can be estimated by @xmath264 . using @xmath265 and @xmath266 , where @xmath121 and @xmath267 are the wavelength of the optical field and the radial beam divergence , respectively
, we find @xmath268 taking @xmath269 and @xmath270  @xcite , we find @xmath271 .
total energy cost to reach steady state per cycle is @xmath272 where @xmath273 . here , @xmath274 is the number of atoms needed for thermalisation , @xmath25 is the injection rate , and @xmath275 is the time elapsed to reach the steady state . in our results
we have found that the harvested work per cycle is much less than the resonance energy , @xmath276 , thus the generation cost of nlap fuel is several orders of magnitude larger than the harvested work @xmath277 , which confirms that the second law of thermodynamics obeyed in photonic carnot engine with nlap .
generation cost is not included in the thermodynamic efficiency but it can be significant figure of merit in the round trip efficiency . to make such photonic carnot engines more appealing for certain applications , it is necessary to increase their round trip efficiency as well .
for that aim one could consider the cases of larger compression ratios ( @xmath278 ) and high operation temperatures ( @xmath279 ) .
our focus here is on the discussion if such engines can operate under decoherence . despite the negative conclusions for two level phaseonium case  @xcite
, we have found that larger phaseonium fuel allows for operational photonic carnot engines .
the question of how to increase their round trip efficiency requires further analysis which is beyond the scope of present contribution .
summarizing , we examined scaling of work and efficiency of a quantum heat engine with the number of quantum resources . specifically , we considered a photonic carnot engine with a multilevel phaseonium quantum fuel .
we derived a generalized master equation for the cavity photons , which forms the working fluid of the engine , and determined the steady state photon number to calculate the work output and thermodynamic efficiency .
we find that they scale quadratically with the number of quantum coherent levels @xmath1 .
we examined the case of degenerate levels to get analytical results and to examine scaling laws against decoherence due to cavity dissipation and atomic dephasing .
we verified our analytical results with detailed numerical methods and have shown consistency of coarse grained analytical results with the microscopical numerical approach .
generation of multilevel phaseonium fuel using morris - shore transformation determined quantum householder reflection technique as well as its cost are examined . using typical parameters in modern resonator systems , such as circuit qed
, our calculations reveal that decoherence due to cavity dissipation could be overcome by the multilevel quantum coherence even in the presence of large dephasing rate . if the dephasing rate increases with @xmath1 , then work and efficiency
can still overcome the decoherence and retain their @xmath157 scaling up to a critical number of coherent levels .
authors warmly thank n. allen , a. imamoglu and i. adagideli for illuminating discussions . authors acknowledge support from ko university and lockheed martin university research agreement .
we generalize micromaser mesoscopic master equation treatment @xcite applied for three - level phaseonium engine @xcite to multi - level case . in @xmath1lap model ,
the hamiltonian of the whole system is @xmath280 where label @xmath35 implies injected @xmath61 arbitrary atom .
we adopt the notation of ref .
@xcite for the hamiltonian and relevant quantities . here , @xmath281 and      in the interaction picture where @xmath283 , @xmath284 are atomic energy levels , @xmath68 is the atom - field coupling coefficient and @xmath5 is the single mode cavity frequency . here , we assume all levels coupled to the excited one with the same coefficient @xmath68 . the equation of motion of overall system is        the equation of motion of the radiation field which is the working substance of the heat engine can be found by tracing out atomic part as @xmath288_{nn})\nonumber\\ + \mbox{tr\scriptsize at}\mathcal{l}_a[\rho]_{nn}+\mbox{tr\scriptsize at}\mathcal{l}_f[\rho]_{nn}\label{eq - mot}\end{aligned}\ ] ] where @xmath57 . here , @xmath289_{nn}={\langle
n|}(h\rho-\rho h){|n\rangle}= \sum_{m}{\langle n|}h{|m\rangle}{\langle m|}\rho{|n\rangle}-\sum_{m}{\langle n|}\rho{|m\rangle}{\langle m|}h{|n\rangle}=\sum_{m}(h_{nm}\rho_{mn}-\rho_{nm}h_{mn})$ ] . in micromaser theory , due to the short atom - cavity interaction time ,
last term , cavity decay is usually ignored when the atom is inside the cavity .
the second term is treated perturbatively and will be considered for zeroth order in g. here , we will keep it but assume it can be treated independently .
one may write the partial trace operation over atomic degrees of freedom for a random single atom as ,    @xmath290_{nn}&=\sum_{\alpha}{\langle \alpha , n|}[h,\rho]{|\alpha , n\rangle}\nonumber\\ & = { \langle a n|}[h,\rho]{|a n\rangle}+\sum_{i=1}^{n}{\langle b_i n|}[h,\rho]{|b_i n\rangle}\label{ptrace}\end{aligned}\ ] ]    where @xmath291 are the atomic basis as expressed at the right hand side .
each term of eq .
( [ ptrace ] ) can be calculated by using the selective rules of the hamiltonian between certain levels @xmath292 and @xmath274 ; for instance inserting for the first term of eq .
( [ ptrace ] ) we have ,    @xmath293{|a n\rangle}&={\langle a n|}(h\rho-\rho h){|a n\rangle}\nonumber\\ & = \sum_{\alpha ' m}\{{\langle a n|}h{|\alpha ' m\rangle}{\langle \alpha ' m|}\rho{|a n\rangle}\nonumber\\ & -{\langle a n|}\rho{|\alpha ' m\rangle}{\langle \alpha ' m|}h{|a n\rangle}\}\nonumber\\ & = \sum_{\alpha ' m}\{h_{an,\alpha'm}\rho_{\alpha'm , an}-\rho_{an,\alpha'm}h_{\alpha'm , an}\}. \label{eq4}\end{aligned}\ ] ]    @xmath294 is the hamiltonian eq .
( [ eq - int ] ) to be inserted into eq .
( [ eq4 ] ) which can be written as , @xmath295 where @xmath296 and @xmath297 . expressing the terms of eq .
( [ eq4 ] ) conveniently , we write @xmath298 and @xmath299 . besides , @xmath300 . substituting these terms into the first part of eq .
( [ eq4 ] ) we have @xmath301 . the second part of eq .
( [ eq4 ] ) is simply the complex conjugate .
hence the first term of eq .
( [ ptrace ] ) is @xmath293{|a n\rangle}&=\hbar g e^{-i\omega t}\sqrt{n+1}\sum_{i=1}^{n}\rho_{b_i n+1,an}-c.c.\nonumber\\   \label{eq7}\end{aligned}\ ] ]    the second term of eq .
( [ ptrace ] ) would be calculated by similar considerations .
we can write @xmath302   { |b_i n\rangle}=\sum_{\alpha ' m}\{h_{b_i n,\alpha'm}\rho_{\alpha'm , b_i n}\}$ ] . in this case
@xmath303 , @xmath304 and @xmath305 .
then , @xmath302{|b_i n\rangle}=\hbar g e^{i\omega t}\sqrt{n}\rho_{an-1,b_in}-c.c.$ ] and finally the second term of eq .
( [ ptrace ] ) becomes , @xmath306{|b_in\rangle}=\hbar g e^{i\omega t}\sqrt{n}\sum_{i=1}^{n}\rho_{an-1,b_in}-c.c.\nonumber\\\label{eq8bi}\end{aligned}\ ] ] inserting these results into eq .
( [ eq - mot ] ) and after some rearrangements we have the field equation of motion , @xmath307_{nn}.\label{9eqmot}\end{aligned}\ ] ] we use @xmath308 hereafter for simplicity . here ,
@xmath309_{nn}&={\langle n|}\frac{\kappa}{2}(2\hat{a}\rho\hat{a}^{\dagger}-\rho\hat{a}^{\dagger}\hat{a}-\hat{a}^{\dagger}\hat{a}\rho){|n\rangle}\nonumber\\ & = \kappa\left\lbrace(n+1)\rho_{n+1,n+1}-n\rho_{nn}\right\rbrace.\label{louvil - rho}\end{aligned}\ ] ]    in order to proceed the calculation , any single term in the summation of eq .  ( [ 9eqmot ] )
should be calculated and inserted therein .
the terms can be obtained by the integration of corresponding equation of motions by using the selecetive rules of the hamiltonian as expressed above .
we evaluate atomic equation of motion @xmath310 order in @xmath68 first , we have      where @xmath312 and @xmath313 .
final two terms of eq .
( [ atomic-eq.mot ] ) corresponds to @xmath314 $ ] with @xmath315 and @xmath316 . here ,
@xmath317 and @xmath318 is taken equal to @xmath73 for simplicity .
note that the atomic part of the master equation is for the case of pure dephasing @xmath319 relaxation and we follow the usual assumption of micromaser theory  @xcite that cavity decay and atomic dynamics can be separately treated . when atom is inside the cavity , decay is not included .
the equation of motion of the @xmath43 term of the first summation of eq .
( [ 9eqmot ] ) for a single atom is @xmath320 here , we have neglected the the matrix element @xmath321{|an\rangle}$ ] in accord with the assumptions indicated above .
the equation of motion for @xmath322 which is the @xmath43 term of second summation of ( [ 9eqmot ] ) could be obtained by simply replacing @xmath323 in eq .
( [ eq10 ] ) , that is @xmath324 therefore we obtain @xmath325 and @xmath326 terms by integrating eq.s  ( [ eq10 ] and [ eq.11 ] ) formally in the following form ,        the terms of eq.s  ( [ eq12 ] and [ eq13 ] ) can be factorized to atomic and field density matrices for zeroth order solution in g. for instance , @xmath329 @xmath330 here , @xmath331 , @xmath332 which are the initial atomic density matrix elements , obey the respective atomic equations of motion @xmath333 in which the solutions are @xmath334 where @xmath335 .
eq.s ( [ eq15 ] ) and ( [ eq16 ] ) imply that excited and ground state levels decay to a lower level .
the off - diagonal elements of the atomic density matrix are equal to @xmath336 .
the @xmath337 is assigned with the coherence preparation .
thus , we can express eq.s  ( [ eq12 ] and [ eq13 ] ) by using zeroth order atomic eq.s  ( [ eq14]-[eq16 ] ) to find first order solutions in g ,          before proceeding , we replace the summation over number of injected atoms by integration over injection time as @xmath341 where @xmath25 is the injection rate .
we also define @xmath342 where @xmath343 .
then , @xmath344+c.c\big\}+\mathcal{l}_f[\rho]_{nn}.\label{eq20}\end{aligned}\ ] ] note that @xmath345 while @xmath346 .
evaluating the integrals in ( [ eq20 ] ) over @xmath347 and @xmath348 after changing integration order as @xmath349 we have    @xmath350-\big[\frac{1}{\gamma}\big(\sum_{i=1}^{n}\frac{1}{-i\delta_i+\gamma}\big)\rho_{b_i b_i}^0\nonumber\\ & + \sum_{i < j}\big(\frac{1}{(-i\delta_i+\gamma)}\frac{1}{(i\omega_{b_i b_j}+\bar{\gamma})}e^{i\phi_{ij}}+\frac{1}{(-i\delta_j+\gamma)}\frac{1}{(-i\omega_{b_i b_j}+\bar{\gamma})}e^{-i\phi_{ij}}\big)|\rho_{b_i , b_j}^{0}|\big]\nonumber\\ & \times\big[(n+1)\rho_{n+1 n+1}-n\rho_{nn}\big]+c.c.\bigg\}+\mathcal{l}_f[\rho]_{nn}\label{eq21}\end{aligned}\ ] ]      @xmath351-\big[\frac{1}{\gamma}\big(\sum_{i=1}^{n}\frac{2\gamma}{\delta_i^2+\gamma^2}\big)\rho_{b_i b_i}^0\nonumber\\ & + \sum_{i < j}\big[\big(\frac{2\cos\phi_{ij}(\delta_i\omega_{b_i b_j}+\gamma\bar{\gamma})+2\sin\phi_{ij}(\omega_{b_i b_j}-\delta_i\bar{\gamma})}{(\delta_i^2+\gamma^2)(\omega_{b_i b_j}^2+\bar{\gamma}^2)}+\frac{2\cos\phi_{ij}(\gamma\bar{\gamma}-\delta_j\omega_{b_i b_j})+2\sin\phi_{ij}(\omega_{b_i b_j}+\delta_j\bar{\gamma})}{(\delta_j^2+\gamma^2)(\omega_{b_i b_j}^2+\bar{\gamma}^2)}\big)|\rho_{b_i b_j}^0|\big]\nonumber\\ & \times\big[(n+1)\rho_{n+1 n+1}-n\rho_{nn}\big]\bigg\}\label{eq22}\end{aligned}\ ] ]            @xmath355 , @xmath356 and @xmath357 .
note that @xmath358 has @xmath1 number of terms and @xmath359 has @xmath15 number of terms in the summation .
since we seek the solutions in the steady state , we obtain the steady state photon number @xmath33 by solving @xmath360 where @xmath361 and we write @xmath362 by using previously obtained @xmath363 as @xmath364                where @xmath371 here , @xmath372 is the average photon number in the absence of atomic coherence in terms of average photon number @xmath373 which is the average photon number in the absence of atomic coherence and in the absence of cavity decay @xmath113 .
finally we look at the degenerate ground state case ( @xmath374 , @xmath375 ) . in this case
, @xmath53 can be simplified to @xmath376 where @xmath377 , @xmath378 for any @xmath379 .
the simplified forms of other parameters are            44ifxundefined [ 1 ] ifx#1 ifnum [ 1 ]
# 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) in  @noop _ _ ,  vol .
( ,  )  pp .   @noop * * ,   ( ) in  @noop _ _ ,  vol .
( ,  )  pp .   @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop _ _ ,  ed .
( ,  ,  ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _ ,  ed .
( ,  ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( )