Quantum heat transport of a two-qubit system: Interplay between system-bath coherence and qubit-qubit coherence

We consider a system consisting of two interacting qubits that are individually coupled to separate heat baths at different temperatures. The quantum effects in heat transport are investigated in a numerically rigorous manner with a hierarchial equations of motion (HEOM) approach for non-perturbative and non-Markovian system-bath coupling cases under non-equilibrium steady-state conditions. For a weak interqubit interaction, the total system is regarded as two individually thermostatted systems, whereas for a strong interqubit interaction, the two-qubit system is regarded as a single system coupled to two baths. The roles of quantum coherence (or entanglement) between the two qubits (q-q coherence) and between the qubit and bath (q-b coherence) are studied through the heat current calculated for various strengths of the system-bath coupling and interqubit coupling for high and low temperatures. The same current is also studied using the time convolutionless (TCL) Redfield equation and using an expression derived from the Fermi golden rule (FGR). We find that the HEOM results exhibit turnover behavior of the heat current as a function of the system-bath coupling strength for all values of the interqubit coupling strength, while the results obtained with the TCL and FGR approaches do not exhibit such behavior, because they do not possess the capability of treating the q-b and q-q coherences. The maximum current is obtained in the case that the q-q coherence and q-b coherence are balanced in such a manner that coherence of the entire heat transport process is realized. We also find that the heat current does not follow Fourier's law when the temperature difference is very large, due to the non-perturbative system-bath interactions.