Heat Modulation on Target Thermal Bath via Coherent Auxiliary Bath

We study a scheme of thermal management where a three-qubit system assisted with a coherent auxiliary bath (CAB) is employed to implement heat management on a target thermal bath (TTB). We consider the CAB/TTB being ensemble of coherent/thermal two-level atoms (TLAs), and within the framework of collision model investigate the characteristics of steady heat current (also called target heat current (THC)) between the system and the TTB. It demonstrates that with the help of the quantum coherence of ancillae the magnitude and direction of heat current can be controlled only by adjusting the coupling strength of system-CAB. Meanwhile, we also show that the influences of quantum coherence of ancillae on the heat current strongly depend on the coupling strength of system—CAB, and the THC becomes positively/negatively correlated with the coherence magnitude of ancillae when the coupling strength below/over some critical value. Besides, the system with the CAB could serve as a multifunctional device integrating the thermal functions of heat amplifier, suppressor, switcher and refrigerator, while with thermal auxiliary bath it can only work as a thermal suppressor. Our work provides a new perspective for the design of multifunctional thermal device utilizing the resource of quantum coherence from the CAB.


Introduction
Quantum thermodynamics mainly studies thermodynamic behaviors emerging in systems that are quantum in nature [1][2][3][4]. Compared with the classical systems dominated by the standard laws of thermodynamics, some novel phenomena can emerge in the quantum systems due to the presence of quantum properties, e.g., quantum coherence [5][6][7][8][9], or entanglement [10][11][12]. For example, the efficiency beyond the Carnot cycle [13], the reversion of heat flowing from the hot system to the colder one [14]. These non-intuitive physical behaviors have become an apparent challenge to the standard laws of thermodynamics [15,16]. The rapid progress of quantum technologies has allowed us to characterize the quantum machine [17][18][19][20] and experimentally realized them in various quantum systems [21][22][23][24][25][26][27][28][29][30]. With the aid of these controllable quantum platforms (systems), some studies focus on the redefinitions of some concepts, such as work and heat [31][32][33][34], and the verification and modification of thermodynamics second law [35] in quantum domain. Others concentrate on heat control/management [36][37][38][39][40][41][42][43] in order to design the of system-CAB, seems to be frozen to the THC, i.e., quantum coherence almost has no effects on the THC at the critical value. For the coupling strength below/over the critical value the THC increases/decreases monotonically as the coherence magnitude of ancillae increases. Thus, whether the coherence is positive, negative, or not correlated with the THC strongly depends on the coupling strength of system-CAB. Especially, the reversal-THC, in some parametric regime of the coherence of ancillae and the coupling strength of system-CAB, could appear which is impossible for the auxiliary bath being thermal bath without coherence. Therefore, with the aid of CAB the three-qubit system can serve as a multifunctional thermal device integrating the functions of heat amplifier, suppressor, switcher and refrigerator. Of particular interest, the multiple functions can be switched only by adjusting a single controllable coupling strength of system-CAB.
The remainder of this paper is arranged as follows: In Section 2, we introduce the scheme of heat modulation of three-qubit system as a quantum thermal device and its dynamics. In Section 3, we mainly focus on the effects of ancilla's coherence (including the relative phase and coherence magnitude) of CAB and coupling strength of system-CAB on the THC, and analyze the thermal functions of system with CAB in different parameter regions. We compare the characteristics of THC in CAB with that in thermal auxiliary bath, and demonstrate the role of coherence of CAB as a resource in thermal modulation. In addition, the effects of temperature on the THC are also discussed. We conclude the whole work in Section 4.

Model
We design a scheme of heat modulation on a TTB via a tripartite system as quantum machine assisted with an auxiliary bath as shown in Figure 1. In our model, the system is composed of three qubits (S a,b,c ) with the frequencies ω a,b,c . The CAB (TTB) consists of a series of identical and independent two-level atoms (TLAs) {L 1 , L 2 , . . . , L N } ({R 1 , R 2 , . . . , R N }) with transition frequency ω L (ω R ) and density matrix ρ L (ρ R ). We assume that the ancillae of left auxiliary bath and the thermal TLAs of right TTB synchronously pass through the left box and the right box one by one, respectively, and two processes denoted as M1 (interaction process among three subsystems) and M2 (interaction process between the system and the baths) are implemented alternatively with equal time interval τ. Here, the boxes simulate triggers to control which interaction channels are on/off in the dynamics of system. Specifically, before the atoms reach the boxes, i.e., in M1 process, there are only interactions among the three qubits (S a,b,c ). Further, when the atoms entering the boxes, i.e., in M2 process, the interaction channels in M1 are off, and the subsystems S a (S c ) couples to the ancillae (thermal TLA) in the left (right) box. The energy exchange among the CAB, system and TTB occurs in this process. After the repeated implementations of M1 and M2 processes, a steady heat current between the system and the TTB can be established. With the aid of this model, the influences of coherence of CAB on THC and the thermal functions of system assisted with CAB can be exploited at length. is first implemented in M1 process; then the subsystems S a and S c are coupled to the nth ancilla (prepared in ρ L ) in coherent auxiliary bath (CAB) and the nth thermal atom (prepared in ρ R ) in TTB, respectively, in M2 process. After that, the two steps of M1 and M2 are implemented repeatedly, and ancilla (thermal atom) in CAB (TTB) interacting with S a (S c ) are refreshed by the next one in each round. Thus, the steady heat current between the system and TTB, after many rounds, is established.
Next, we focus on the specific dynamics of system in a single round (i.e., to implement the M1 and M2 process once). Here, we denoteĤ i = ω iσ z i /2, (i = L, R, a, b and c) as the free Hamiltonians of ancillae in CAB and thermal TLAs in TTB and subsystems S a , S b and S c , respectively, where ω i andσ z i = |1 1| − |0 0| are independently the transition frequencies and Pauli matrices of the TLAs with subscript i, and |1 (|0 ) describing a TLA in the excited (ground) state. The free Hamilton of system is denoted asĤ 0 =Ĥ a +Ĥ b +Ĥ c . In M1 process, the interaction Hamilton of system is given as: with g m is the coupling constant among three subsystems.
being the raising (lowering) operator) representing the off-resonant and the resonant terms, respectively. Here, the XXX-type Hamiltonian as a more likely candidate of three-qubit physical Hamiltonian [82] has been exploited widely in spin system [82][83][84][85], and can be implemented experimentally [59]. Further, its dynamics effectively reduces to that of the resonant couplingV R on phenomenological grounds [86,87] when the frequencies are resonant ω b = ω a + ω c and the coupling is weak g m ω i i.e., the rotating wave approximation applies (the off-resonant termV OR in the interaction Hamiltonian (1) can be omitted). In this paper, we also consider that the condition of resonant frequencies is satisfied. According to the scheme in Figure 1, the dynamics of the system in M1 process for t ∈ (t n−1 , t n−1 + τ] is given as: where the mapping Λ 1 describes the unitary evolution of system, and the unitary operator isÛ 1 (τ) = exp(−iH S τ/ ) withĤ S =Ĥ 0 +V 1 being the total Hamilton of system. ρ n−1 S (ρ n S ) represents the state of system S before (after) the (n − 1)th ancilla reaches (leaves) the box, that is, the state of beginning (ending) of system's internal interaction at time t n−1 (t n−1 + τ). One can identify that the free Hamilton ofĤ 0 does not commute with the unitary operatorÛ 1 (τ) for the finite coupling g m , i.e., Ĥ 0 ,Û 1 (τ) = 0 which implies that the energy of system is non-preserving in this process due to the off-resonant term V OR introduced. Here, the off-resonant termV OR characterizes the external driven on the system by an external agent (or work source). Any energy changes of system in M1 process are solely due to energy leaving or entering the work external agent. In the M2 process for t ∈ (t n−1 + τ, t n ] with t n = t n−1 + 2τ (i.e., at the time interval of the nth ancillae (TLA in TTB) kept in the left (right) box), the interaction channel in M1 process is switched-off, and the subsystems S a and S c are coupled to the ancillae in left box and the thermal TLA in right box, respectively. The interactions are expressed as: whereV S,L (V S,R ) and g L (g R ) represent the interaction Hamilton and the coupling constant of S a -ancillae in the left box (S c -TLA in the right box), respectively. The dynamics of system and TTB are unitary and can be described as: where ρ n S ( ρ n R ) denotes the state of the system (nth thermal TLA of TTB) at time t n = 2nτ (i.e., the moment of the nth thermal TLA in TTB just leaving the right box) the mapping Λ 2 corresponds to a Markov process andÛ 2 (τ) =Û S,L (τ)Û S,R (τ) with: Here, it is noted that as usually treated in most works of collision models [35,65,[70][71][72][73] we have considered that the state of system after its interaction with two baths is embodied by the stroboscopic map. That ensures the system is being always independent to the particle units (ancillae/thermal-TLAs in CAB/TTB) that have collided, and the dynamics of system a memoryless Markov process. Physically, the collision model protocol is consistent with the inspiration from Boltzmann's original Stosszahlansatz. For instance, a particle in Brownian motion interacts with only a few water molecules at a time. Moreover, this interaction lasts for an extremely short time, after that the molecule moves on, and never to return [88]. Since the environment is large, the decoherence and dissipation ensure that the next molecule to arrive will be completely uncorrelated from the previous one, so the process repeats anew [35]. The same scenario is also suitable for the multipartite collision models [79,81], and addressing the thermodynamics of engineered reservoirs [37,[89][90][91]. In terms 89 of the exchange interaction ofV S,L(R) given in Equation (3) one can identify that the commutation relation Ĥ 0 S,L(R) ,Û S,L(R) (τ) = 0 withĤ 0 S,L =Ĥ a +Ĥ L andĤ 0 S,R =Ĥ c +Ĥ R holds, which indicates that the total energy of whole system (system plus CAB and TTB) is preserved in this process. Here, it is mentioned that the unitary evolutionÛ S,L(R) (τ) in Equation (5) above corresponds to a swap gate operation, and we can rewrite it as: where λ L(R) = g L(R) τ is the dimensionless coupling strength of system-ancillae in auxiliary bath (system-TLAs in TTB),Î L(R),a(c) is the 4 × 4 identity operator, andŜ sw L(R),a(c) = |11 11| + |00 00| + |01 10| + |10 01| is the two-TLAs swap operator, havingŜ sw L(R),a(c) = ψ a(c) ⊗ ψ L(R) for all ψ a(c) , ψ L(R) ∈ C 2 . Therefore, from Equation (6), the swap strength between two particles is determined by sin(2λ L(R) ) varying with λ L(R) at the period π/2. Further, when λ L(R) = kπ/2 (k =0, 1, 2, . . . , n) the swap strength sin(2λ L(R) ) is zero conrresponding to the weakest swap strength, and for λ L(R) = kπ/2+π/4 it becomes one characterizing the strongest swap strength (complete swap) between the system and the ancillae in CAB (thermal atoms in TTB). It also means that the energy exchange between the system and the baths (CAB and TTB) will be periodical with the same period as that of swap strength. In our model, only the heat exchange between the system and the TTB occurs in the M2 process due to the energy-persevering evolution. In the same spirit as the definition of heat (or heat flow) in [35,61,78,92,93], the amount of heat exchange between the system and the TTB can be quantified by the energy change of thermal TLAs of TTB in each round. In the arbitrary n-th round the amount of heat flowing to TTB reads: where ρ n R ( ρ n R ) is the state of the n-th thermal TLA in TTB at the beginning (end) of the subsystem S c in M2 process. Here, ∆Q n > 0 (∆Q n < 0) represents that the system pumps heat into the TTB (the TTB delivers heat to the system). It is pointed out that the state ρ n R of the n-th thermal TLA is initially a thermal state and, the state, after the mapping M2, ρ n R remains in a diagonal distribution in the eigenbasis ofĤ R in our model (i.e., both of them are the mixed states with no coherence), and the energy spectrum ofĤ R is kept unchanged in the process: ρ n R → ρ n R . Based on the heat exchange in Equation (7) between the system and the TTB, the average heat current can be defined as: The positive (negative) current J n > 0 (J n < 0) in Equation (8) indicates the heat flowing into (out of) the TTB, which also means that the system works as a heat pump (J n > 0) (refrigerator (J n < 0)) to heat (refrigerate) the TTB. In a long-time limit, the steady THC denoted as J SS can be established: It is well known that when a finite system contacts with an infinite heat/nonequilibrium bath it will relax to (or be thermalized into) a steady state as the time increases, and a dynamical equilibrium can be established among the system and the baths, i.e., the steady heat current emerged. This mechanism is also suitable for our collision model. Specifically, the interaction between the three-qubit system and the two baths in each round is used to mimic the thermal contact process in the conventional model, and the state of system, after each round, can be updated once including the populations and quantum correlation/coherence among qubits. As the collision time increases the system's state is modified by the baths less and less gradually, and the system, after many rounds, will reach a steady state associated with the steady heat current. In the following subsections, we are mainly concerned about the behaviors of THC and the thermal functions of the system with a thermal/coherent auxiliary bath. For the sake of brevity, the THC mentioned in the following subsections refers to the steady THC.

Initial States of System and Baths
We consider that the CAB, system and TTB are initially uncorrelated, and the initial state of the whole composite system (system plus CAB and TTB) is: Here ρ S (0) is the initial state of system being a thermal product state of three subsystems given as: where Z Ĥ = tr exp −βĤ is the partition function with the inverse temperature ⊗N represent the initial states of CAB (ρ tot L (0)) and TTB (ρ tot R (0)) where each bath of CAB and TTB is composed of identical units. Further, the initial state of each thermal TLA ρ R (0) in TTB and ρ L (0) for the ancillae in CAB are independently given by: and where |φ L = √ P 00 |0 + e iθ √ P 11 |1 with θ being the relative phase and is a weight parameter determining the proportion of two components, thermal state ρ β and pure coherent state |φ L , in the state ρ L (0). Here, it is noticed that the parameter α only appears in the non-diagonal elements of ρ L (0), and the diagonal elements of ρ L (0) are the same with the ones of thermal state ρ β . According to the coherence measure of l 1 -norm, for an arbitrary state ρ the coherence reads [5] with ρ mn (m = n) the non-diagonal elements of ρ. One can identify that for the state ρ L (0) given in Equation (13) the l 1 -measure of coherence C l 1 (ρ L (0)) = 2α √ P 00 P 11 is proportional to α. That is, when fixing the population P 00 (P 11 = 1 − P 00 ) of ρ L (0) in the energy basis {|0 , |1 } ofĤ L the coherence C l 1 (ρ L (0)) of state ρ L (0) increases as α increases, and for (α = 0) α = 1 the state ρ L (0) reduces to a complete mixture state without coherence C l 1 (ρ L (0)) = 0 (pure state with the maximum of coherence in the range of 0 ≤ α ≤ 1, i.e., max α∈ [0,1] C l 1 (ρ L (0)) = 2 √ P 00 P 11 ). Thus, the weight α can be regarded as an indicator to measure the coherence magnitude of ρ L (0). For simplicity, we, in the next subsections, will take α-indicator instead of C l 1 to measure the coherence magnitude of ancillae of CAB.

Thermal Modulation with Thermal Auxiliary Bath
First, we consider that the auxiliary bath is a thermal bath without coherence, i.e., each ancilla is in a thermal state, and the situation for CAB is provided in the latter subsections. Here, we focus on the modulation of heat current J SS by the dimensionless coupling strength λ L (λ L = g L τ) between the system and the ancilla with state ρ L (0) = ρ β (ρ β given in Equation (13)) being a thermal state. For simplicity, we set the ancillae of auxiliary bath, subsystems S a and S c , and the thermal TLAs of the TTB with the same transition frequencies, ω L = ω a = ω c = ω R = ω and ω b = 2ω, and the Planck constant = 1, throughout the paper. Though some simplified treatments of parameters have been done in our model it is still hard to get the exact analytical solution of THC given in Equation (9) due to the high dimensions of system. Thus, we will investigate the features of THC numerically below when the auxiliary bath is introduced. By numerical calculations, we find that the THC J SS behaves as a cosine-like periodical function of λ L with the fixed period T λ = π/2, and can be fitted with the form: where the parameters A, ϕ 0 and ∆ have a complex relationship with ω, β, λ R and ξ m (ξ m = g m τ representing the coupling strength among three subsystems).
According to Equation (9) we, in Figure 2, plot the variation of THC J SS with the coupling strength λ L , λ L ∈ [0, π] when fixing the other parameters: λ R = 0.4π, ξ m = 0.15π, τ = 0.01, ω = 10 and β = 0.01 (see the red dotted curve). Further, the blue dotted line corresponds to the fitting function J SS with A = 0.567, ϕ 0 = 0.055 and ∆ = 3.445. We can see that numerical results (red curve) of the heat current J SS is basically matched with that of the fitting function J SS (blue curve), i.e., J SS ≈ J SS . Thus, the THC can be modulated in cosine-like form of λ L and satisfied as J SS (λ L ) = J SS (λ L + kT λ ) for k being zero or positive integer. The maximum (minimum) of THC is J Ther kT λ + /4) corresponding to the minimum (maximum) strength of populations (or energy) exchange between the ancilla and the subsystem S a via the swap operation given in Equation (6). Physically, the nonzero steady heat current also implies that nonequilibrium steady state of system is reached. For simplicity, to denote the increment of system's energy as ∆E S injected by the external work source into system in M1 step of each round, the energy increment of system ∆E S always equals to the sum of increasing amount of two baths' energy, ∆Q L + ∆Q R (the subscript L and R represent the left auxiliary bath and the right TTB, respectively), in M2 step, i.e., ∆E S − (∆Q L + ∆Q R ) = 0 implying no net energy accumulation in the nonequilibrium steady state dynamics. The stronger the swap strength |sin(2λ L )| of system-ancilla (i.e., increasing or decreasing λ L in the range of kπ/2 ≤ λ L ≤ kπ/2 + π/4 or kπ/2 + π/4 ≤ λ L ≤ (k + 1)π/2, k = 0, 1, 2, . . .) is, the larger the energy ∆Q L captured by the auxiliary from the system becomes, and that results in the smaller ∆Q R implying smaller THC J SS , vice versa. Therefore, the variety of the THC J SS in Figure 2 with the coupling strength λ L has an opposite trend with that of swap strength, |sin(2λ L )| with fixed period π/2. In addition, we can see that for a finite λ L , λ L ∈ (0, T λ ), the THC J Ther max > J SS (λ L = 0) > J Ther min with J Ther min = ∆ − A > 0 and J Ther max ≈ J SS (λ L = 0) = A cos ϕ 0 + ∆. In order to observe how the auxiliary bath influence the THC in our model, it might be appropriate take the value of THC for no interaction (i.e., λ L = 0) of system-ancilla as a reference THC denoted as J ref , with J ref = J SS (λ L = 0). Using it, the amplification/suppression of THC can be described intuitively, that is J SS > J ref and 0 < J ref < J SS respectively correspond to the heat amplification and suppression. Here, the refrigeration is indicated by the reversal current J SS < 0. Thus, Figure 2 shows that the quantum machine with the thermal auxiliary bath can only serve as a heat suppressor due to the THC being suppressed for a finite λ L , and zero-current or reversal-current cannot emerge, that is, 0 < J SS (λ L = 0)< J ref . Meanwhile, by changing the coupling strength λ L the thermal suppressor can modulate THC in the form of cosine-like function in the range of ∆ + A ≥ J SS ≥ ∆ − A. Here, in order to conveniently describe the capability of quantum machine to modulate the THC by adjusting the controllable coupling strength, we define modulation width of the quantum machine on THC as the difference of the maximum and the minimum of THC, i.e., D := J SS max (λ L ) = J SS min (λ L ) (16) with 0 ≤ λ L ≤ T λ . As for a thermal control device, one usually expects that it could control the THC varying over a wide range as much as possible. That is, the larger the modulation width of THC is, the better the performance of quantum machine is. In terms of characteristics of THC given in Equation (15), the system with the thermal auxiliary bath can only work as a thermal suppressor with the modulation width D = 2A.

Thermal Modulation with CAB
We, in this subsection, mainly focus on the behaviors of THC and the thermal functions of quantum machine assisted with a CAB.

Effects of Relative Phase on THC
Next, we study the effects of the relative phase θ and the coupling strength λ L of ancillae on the THC J SS numerically for a fixed magnitude of coherence below. From numerical calculations we find that when the thermal auxiliary bath is replaced by the CAB the THC J SS is a periodical function of λ L and θ, i.e., J SS (λ L + T λ , θ + T θ )= J ref (λ, θ) with T λ = π/2 and T θ = π as shown in Figure 3, where each ancilla of CAB is prepared in the same state, ρ L (0) = |φ L φ L | given in Equation (13) with the maximum magnitude of coherence α = 1. It is noted that the periods T λ = π/2 and T θ = π are independent of the other parameters in our model. In order to demonstrate the characteristics of J SS and the thermal functions of quantum machine clearly, the variation of J SS in a single period with 0 ≤ θ ≤ T θ and 0 ≤ λ L ≤ T λ is shown in Figure 4 where all the other parameters are the same as that given in   indicates the values of THC at some certain working points (θ, λ u ) on the purple dotted line are the same as reference current J ref (λ L = 0) as defined before). One can see that the HPAR only lies in the region of λ L below λ u , (i.e., the region of λ L < λ u ), and other function regions of HPSR, SR and RR for λ L over λ u , (i.e., the region of λ L > λ u ). Further, the SR and RR appear in the middle region of parametric space: θ ∈ [0, T θ ] and λ L ∈ [λ u , T λ ], and the reminder part for the HPSR. Third, for the relative phase θ taken in RR the quantum machine could perform the different thermal functions in turn: thermal amplifier→stabilization→suppressor→switcher→refrigerator→suppressor as λ L increases from zero to the maximum T λ in a period.
From Figure 4b, one can see some specific behaviors of J SS varying with θ and λ L . Firstly, for a finite relative phase θ the THC J SS always behaves as a sine-like function with respect to λ L . That is, the THC J SS first increases from the initial value J ref to the maximum then decreases to the minimum below J ref , and monotonically returns to its initial value again. Secondly, for the coupling strength λ L being about 0.24π denoted as λ 0 = 0.24π the values of J SS with different θ approach to the same value of J ref , i.e., J SS (θ, λ u ) ≈ J ref which means that for the critical coupling strength λ 0 = 0.24π the THC is almost independent of the relative phase. Thirdly, in the region of J SS > J ref (HPAR) and fixing λ L the closer the relative phase θ is about to 0.45π denoted as θ S = 0.45π, the larger the THC J SS becomes (except for the larger θ such as θ/π = 0.75, 0.9, 1.0 where the three curves of J SS are basically overlap, i.e., the less influence the larger relative phases have), but for the region of J SS < J ref the opposite is true (where J SS decreases as θ closes to θ S , and it can be less than or equal to zero seen the segments below the black solid line with J SS = 0). Based on the numerical simulations above, we show that the relative phase of ancillae is related to the THC and could be regarded as a useful resource to modulate the THC well. Meanwhile, at the suitable relative phase of ancillae the quantum machine can integrate multiple functions, such as thermal amplifier, stabilizer, suppressor, switcher and refrigerator, and these functions can be switched only by adjusting the coupling of system-CAB.

Effects of Coherence Magnitude on THC
Next, we mainly concern the influences of coherence magnitude α on the THC J SS for a fixed relative phase θ S (here, θ S = 0.45π corresponds to the largest modulation width of J SS (seen in Figure 4a or Figure 4b) where D(λ L , θ) has been given in Equation (16)).
In Figure 5a Figure 5a it can be seen that for small coherence magnitude α (about α < 0.2) the HPAR with λ L < λ 0 becomes very narrow (i.e., for α < 0.2 the coupling strength λ 0 in HPIR is also small), and the rest region of λ L is for the HPSR with λ L < λ 0 . It means that for the small α the quantum machine could serve as a thermal amplifier in the regime of λ L < λ 0 , and the thermal suppressor in the regime of λ L > λ 0 . However, as α increases the domain of HPAR/HPSR enlarges/shrinks rapidly, and when α is beyond a certain value (about 0.6 with min J SS (λ 0 , α = 0.6) ≈ 0 seen the Figure 5b, the SR and RR can appear in the regime λ L > λ 0 , and both of SR and RR enlarge as α increases. To demonstrate the dependence of J SS on α and λ L clearly, we, in Figure 5b, plot the variations of J SS with respect to λ L for some selected α, α = {0, 0.15, 0.3, 0.45, 0.6, 0.75, 0.9, 1.0}. Further, some specific features of J SS have been shown in Figure 5b. Firstly, for a fixed α the THC J SS always rises at first, then falls down, and rises up its initial value again.
Meanwhile, the modulation width D(λ L , α) of the THC becomes larger with the increase of α. Secondly, we notice that for λ L = λ 0 with λ 0 = 0.25π, the THCs with different α have the same value as the one for thermal auxiliary bath (α = 0), i.e., J SS λ 0 , α = J SS λ 0 , 0 with α ∈ [0, 1], which means that at the critical coupling λ L = λ 0 the effects of coherence on THC can be frozen, i.e., the THC is independent of the coherence magnitude of ancillae. For 0 < λ L < λ 0 the THC can be amplified by the coherence of CAB, i.e., J SS (λ L , α) > J SS (λ L , 0) (seen for 0 < λ L < λ 0 all the lines of J SS with different non-zero α are above the black solid line with α = 0), and for a fixed λ L the THC J SS always increases monotonously with the increasing of coherence magnitude α, which implies that the THC is positively correlated with the coherence magnitude of ancillae. However, for λ L > λ 0 the situation is opposite where the THC J SS is suppressed, J SS (λ L , α) = J SS (λ L , 0) (seen all the lines with different α are below the one with α = 0 for λ L > λ 0 ), and the value of J SS (λ L , α) with fixed λ L decreases monotonously with increasing α, which demonstrates that the THC is negatively correlated with the coherence magnitude of ancillae in the regime of λ L > λ 0 . Especially, for the ancillae with strong coherence (α > 0.6) the zero-or reversal-current, J SS ≤ 0, can appear (each curve of J SS with different α, (α = 0.75, 0.9, 1.0), has one segment below the black solid line with J SS = 0) which means that for the CAB with strong coherence the quantum machine can also serve as a switcher or a refrigerator.
Based on the analysis above, one can find that the THC is related to not only the coupling of system-ancillae but also the coherence magnitude of ancillae. Meanwhile, the influences of coherence magnitude of ancillae on the THC, such as amplification, suppression and reverse of THC, strongly depend on the coupling strength of systemancillae. This can be understood that when thermal auxiliary bath is replaced with the CAB the system will reach a new steady state associated with the ancilla's coherence magnitude and relative phase. Based on some previous researches on the thermalization problem including the thermalization of the TLA/micro-cavity as system by a coherent TLAs/threelevel atoms (or atomic-pairs) bath (seen in [94][95][96]) it is known that except for the population of bath units (the coherent TLAs/three-level atoms (or atomic-pairs)) the coherence in bath units and the coupling strength between the system and the units have a nonlinear effect on the coherence and the population of system, and the respective contributions of the coherence (coherence magnitude and relative phase) and the coupling strength to the population of system at steady state cannot be separated though the coherence and the coupling strength are independent parameters. Meanwhile, compared with the thermal bath without coherence the excite populations of system at steady state cannot only increase but also decrease which is determined by the parameters of coherence and coupling strength together. This might be why the THC can exhibit rich behaviours when the TTB is introduced. Thus, based on the characteristics of THC one can see that for a prepared CAB the THC can be modulated well by adjusting the controllable coupling of system-ancilla. Especially, for the ancillae with strong coherence (large α) the THC can go through all the regions: HPAR, HPIR, HPSR, SR and RR orderly only by increasing the coupling strength λ L from zero to T λ .

Maximum and Minimum of THC and Modulation Width
We have shown the characteristics of THC and the thermal functions of quantum machine at fixing coherence magnitude α = 1 (or relative phase θ = 0.45π) before. It demonstrates that quantum machine could behave the different functions in different parametric space of λ L , α and θ. Now, we observe which thermal functions the machine can perform when the coherence parameters, α and θ (0 ≤ α ≤ 1 and 0 ≤ θ ≤ π), of ancillae are arbitrary.
Due to the continuous variation of THC with the parameters λ L , α and θ, the thermal functions of quantum machine are determined by the maximum J SS max and the minimum J SS min of THC, i.e., to judge which functional regions of HPAR, HPIR, HPSR, SR and RR are covered in the range of J SS min~J SS max . Therefore, we plot the maximum and the minimum of THCs, J SS max and J SS max , in Figure 6a and the corresponding modulation width of THC, Figure 6b in the full parametric space of λ L , α and θ with 0 ≤ λ L ≤ T λ , 0 ≤ α ≤ 1 and 0 ≤ θ ≤ π. The lower and upper colorful surfaces in Figure 6a correspond to the variations of J SS max and J SS min , respectively, and the sliver gray ellipsoid-like surface for the working points of switcher with J SS = 0. It is noticed that the values of J SS max,min are characterized by the colors on surfaces not the height of the surfaces (the height of surface is for the vertical axis λ L ). From Figure 6a, several features of J SS max and J SS min have been shown. Firstly, it can be seen that the surface of J SS max is always above the one of J SS min which indicates that for the fixed coherence parameters of (α, θ) the maximum of THC J SS max corresponds to the small λ L being about 0 < λ L < 0.15π, and the minimum of THC J SS min for the large λ L with 0.25π < λ L < 0.45π. Secondly, J SS max and J SS min behave as the positive and the negative correlated to the coherence magnitude α, respectively, i.e., the value of J SS max (J SS min ) increases (decreases) as α increases, which implies that the modulation width, D = J SS max − J SS min , of THC is an increasing function of α as shown in Figure 6b. Thirdly, the maximum (minimum) values of THC are always larger (smaller) than the reference THC, i.e., J SS . Therefore, J SS max always remains in HPAR, and J SS min could distribute into three regions: HPSR (0 < J SS min < J ref ), RR (J SS min = 0) corresponding the parts outside and inside the gray ellipsoid-like surface, and SR (J SS min = 0). In addition, from Figure 6b one can see that the larger of the coherence magnitude α and the closer to relative phase θ is to 0.45π, the wider the modulation width of THC becomes, and the better the performance of quantum machine is. With respect to the characteristics of J SS max,min above, one can infer that the quantum machine assisted with CAB can always work as a heat amplifier or a suppressor by tuning the coupling strength λ L , λ L ∈ [0, π/2], due to both regions of HPAR and HPSR being covered in the full parametric space, (α, θ), of coherence. Meanwhile, for some regimes of (α, θ) with J SS min = 0 (or J SS min < 0) the quantum machine can also work as a switcher (a switcher or a refrigerator) due to the SR (or RR) being involved. Especially, in full coherent parametric space of 0 < α ≤ 1 and 0 ≤ θ ≤ π the maximum (minimum) of J SS max (J SS min ) exists at α = 1 and θ = 0.45π, i.e., max J SS max (α, θ) = J SS max (1, 0.45π) (min J SS min (α, θ) = J SS min (1, 0.45π) < 0) which also corresponds to the largest modulation width of THC (seen Figure 6b). Thus, the CAB consisting of the ancillae with α = 1 and θ = 0.45π can be regarded as the optimal CAB, in which the quantum machine not only integrates all thermal functions as a heat amplifier, suppressor, switcher and refrigerator, but also could perform the strongest capabilities in heat amplification and refrigeration.

Effect of Temperature on THC
From numerical simulation we investigate the effects of temperature on THC. We, in Figure 7, plot the variation of target heat current (THC) J SS with the temperature of TTB, T R , and the coupling strength, λ L , for fixed β L = 0.01. Further, Figure 7a,b respectively depict the behaviors of THC for the auxiliary being the CAB and the thermal auxiliary bath. From Figure 7a one can see that two obvious features of THC J SS have been shown. First, for the fixed λ L the THC J SS always decreases as the temperature T R increases, and approaches to the steady value for the high enough T R . Second, for different T R features of the THC J SS varying with λ L are similar. However, for different λ L the THC J SS can exhibit different behaviors as T R increases. For example, in the range of 0.25π ≤ λ L ≤ 0.47π the THC J SS always decreases from a relative large positive value to zero, and then increases in the opposite direction. This means that the quantum machine in this coupling region can perform a heat pump (J SS > 0 region outside the white curve), switching (J SS = 0 the working points on the white curve) or refrigerator (J SS < 0 region inside the white curve) which depends on the temperature of TTB. It is noted that for the coupling strength about λ L ≤ 0.25π and 0.47π ≤ λ L ≤ 0.5π the heat always flow into the TTB J SS > 0 no matter T R is higher or lower than T L . Thus, the machine assisted with the coherent auxiliary bath can also exhibit thermal-diode-like action in the certain parametric regimes. Further, the amplification (refrigeration) of the machine with CAB can be enhanced in the heating (cooling) region with T R < T L (the left side of black solid line) (T R > T L (the right side of black solid line)). That is, the different temperatures could boost the performance of the machine in the certain region of λ L . Compared with the situation of CAB for the thermal auxiliary bath shown in Figure 7b, though the THC also decreases as the temperature T R it cannot be less than or equal to zero in whole parametric space of 0 ≤ λ L ≤ 0.5π and 0 ≤ T R ≤ 250, which implies that the device assisted with the thermal auxiliary bath cannot work as a thermal switch or refrigerator.

Conclusions
In summary, we have proposed a scheme of heat modulation via a three-partite system assisted with a CAB to control the magnitude and the direction of heat current between the system and the TTB. We have analyzed the influences of quantum coherence of ancillae, the coupling strength of system-ancillae and the temperature of two baths on the THC at length. It is shown that for the thermal auxiliary without coherence, the THC behaves as a cosinelike variation with the coupling strength of system-ancillae and a suitable fitting function has been given. Meanwhile, under the thermal auxiliary bath, no matter the high or the low temperature of TTB the heat always flows into the TTB implying the device assisted with the thermal auxiliary bath only working as a heat pump. However, replacing the thermal auxiliary bath with the CAB, due to the influence of coherence (including the coherence magnitude of and the relative phase) of ancillae the THC, in certain parametric regimes of coherence and the coupling of system-ancillae, could exhibit rich behaviors, such as heat amplification, heat suppression, zero-and reversal-current. Therefore, the three-partite system assisted with the CAB could serve as a multifunctional thermal device integrating with heat amplifier, suppressor, switcher and refrigerator. Via the analysis of the maximum and the minimum THC in full coherence parametric space of coherence magnitude and relative phase, the optimal CAB can be suggested, in which the modulation width of THC is largest, and the machine could perform the strongest capabilities in heat amplification and refrigeration. Besides, it has been demonstrated that the different thermal functions can be switched flexibly only by adjusting the coupling strength of system-ancillae, which is convenient for practical application.
Our research might shed some light on the role of resource of quantum coherence outside the system, could boost the deep understanding of thermodynamic properties of quantum coherence, and provide a new perspective for the design of multifunctional thermal management device with the aid of a non-equilibrium auxiliary bath.