the josephson effect has implications going beyond electrical transport @xcite . in 1965 ,
maki and griffin predicted that the interplay between cooper pairs and quasiparticles in tunneling events would provide heat currents with quantum coherence @xcite .
this phenomenon manifests itself as a phase - dependent component of the heat current flowing through a thermally - biased josephson junction ( jj ) , which should , therefore , enable the realization of heat interference . towards this end ,
a _ magnetic flux - controllable _ superconducting heat interferometer was recently theoretically conceived @xcite and subsequently realized experimentally @xcite .
this achievement served to show , on the one hand , how quantum coherence extends also to dissipative observables such as heat current . on the other
, these devices might constitute the building blocks for the implementation of superconducting hybrid _
coherent caloritronic circuits _ consisting of , for instance , thermal modulators , heat transistors , splitters , etc .
the purpose of this paper is to review recent advances done within the emerging field of _ coherent caloritronics_. the paper is organized as follows : we dedicate sec .
[ sec00 ] to distinguish first between coherent caloritronics and thermoelectric effects also present in temperature - biased jj . in sec .
[ sec02 ] we describe the general model used to derive the behavior of the _ heat _ current flowing through the jj . in sec .
[ sec1 ] a superconducting heat interferometer is analyzed theoretically and its operation is demonstrated experimentally . in sec .
[ sec2 ] we discuss the implementation of a double - loop heat modulator an even more versatile device is theoretically analyzed in sec .
it consist of a temperature - biased extended jj that serves to produce coherent diffraction of the thermal currents . in sec .
[ sec5 ] the remarkable heat rectification properties of nis and sis junctions are addressed ( where n stands for a normal metal , i represents a thin insulating layer and s and s refer to two different superconducting electrodes ) .
this review is enriched with realistic calculations showing the achieved modulation of the electronic temperature of a normal metal electrode placed nearby each of the aforementioned structures .
finally , our results are summarized in sec .
[ summary ] along with a few final remarks .
it is convenient to begin by clearly stating the nature of the phenomena we will deal with along this manuscript .
it is our goal to investigate the intriguing quantum behavior of the _ heat _ current @xmath0 flowing through a temperature - biased jj . here
, we will concentrate on jjs consisting of two superconducting electrodes separated by a thin insulating layer , i.e. tunnel sis junctions .
these studies are not to be confused with thermoelectric effects also present in jjs @xcite .
the latter concern the consequences of a thermal bias in the _ electric _ current flowing through the junction and manifests as the generation of a potential difference of thermoelectric nature .
specifically , a thermovoltage @xmath1 might arise across the junction when the electronic temperature difference between the superconducting electrodes exceeds a threshold value @xmath2 @xcite .
these effects can be understood by assuming that the temperature drop across the junction gives rise to a normal electric current ( @xmath3 ) .
as noted by ginzburg @xcite , this current is shorted out in superconducting materials by the generation of a supercurrent ( @xmath4 ) flowing in the opposite direction . as a result ,
the voltage difference generated across the junction remains zero .
if , however , the temperature drop is large enough to generate @xmath5 , @xmath6 being the critical current of the jj , one enters into the dissipative regime leading to the appearance of @xmath1 .
the smallness of this effect has complicated considerably its experimental observation , making it measurable only in certain kind of jjs such as the sns type . in this way ,
kartsovnik _ et al _ observed the generation of @xmath7 v upon temperature differences of the order of @xmath8 k in ta - cu - ta sns sandwiches @xcite .
thermoelectric effects on response to an in - plane magnetic field were also studied by ryazanov _ et al _ in ref .
specifically , @xmath2 was found to depend on the modulus of the cardinal sine function of the magnetic flux piercing the junction leading to a fraunhofer - like dependence .
the influence of thermoelectric phenomena in sis junctions was analyzed in ref .
according to these calculations , a temperature gradient of @xmath9 k across a tunnel junction with normal state resistance @xmath10 @xmath11 will generate a cooper pair thermocurrent @xmath12 a. the latter will imply , in turn , the existence of a tiny phase difference between the superconducting electrodes of only @xmath13 rad that can be neglected for the considerations made here .
we will turn our attention now to the description of the main equations governing heat transport in jjs .
we start by assuming a josephson tunnel junction between two superconducting electrodes @xmath14 and @xmath15 characterized by a phase difference @xmath16 as schematized in fig . [
fig01](a ) . under such circumstances , an _
electric _ supercurrent will flow that depends on @xmath17 .
assume now that the jj is electrically open but a thermal gradient is established by heating up the quasiparticles of @xmath14 up to @xmath18 while keeping quasiparticles in @xmath15 at temperature @xmath19 .
then , a stationary heat current @xmath20 will flow from the hot to the cold reservoir .
and s@xmath21 separated by a thin insulating layer constitute a conventional jj .
if a phase difference @xmath22 arises , an electric supercurrent @xmath23 will flow .
( b ) if s@xmath24 and s@xmath21 are kept at temperatures @xmath18 and @xmath19 ( with @xmath25 ) , respectively , a steady - state total heat current @xmath26 will flow from s@xmath24 to s@xmath21 . ]
the explicit form of @xmath20 was calculated by maki and griffin in 1965 @xcite giving .
@xmath27    equation ( [ heatcurrent ] ) consists of two terms ; @xmath28 is the usual heat flux carried by quasiparticles @xcite and @xmath29 is the phase - dependent part of the heat current , which is peculiar of josephson tunnel junctions . quoting the authors
, the latter term `` gives rise to an oscillatory heat flux and is analogous to a similar term in conventional electric transport through josephson junctions '' .
we stress the fact that , depending on @xmath30 , the phase - dependent component of the heat current can flow backwards , i.e. , opposite to the thermal gradient .
the total heat current goes however from the hot to the cold reservoir , preserving the second principle of thermodynamics @xcite .
it is worthwhile to stress as well the absence of any cooper pair - related energy term .
this is so for , in a static situation , the superconducting condensate shall carry no entropy @xcite .    the quasiparticle term appearing in eq .
( [ heatcurrent ] ) reads @xcite @xmath31 , \label{qqp } \\end{aligned}\ ] ]    where @xmath32 $ ] is the quasiparticle bcs normalized density of states in s@xmath33 at temperature @xmath34 ( @xmath35hot , cold ) and @xmath36 is the energy measured from the condensate chemical potential .
furthermore , @xmath37 is the temperature - dependent superconducting energy gap and @xmath38 , @xmath39 being the critical temperature of the superconductor , @xmath40 , @xmath41 is the heaviside step function , @xmath42 is the boltzmann constant , @xmath43 is the junction normal - state resistance , and @xmath44 is the electron charge .
on the other hand , for the phase - dependent component of the heat current one gets @xcite @xmath45 , \label{qint}\end{aligned}\ ] ] where @xmath46 $ ] is the cooper pair bcs density of states in s@xmath33 at temperature @xmath34 @xcite .
this term originates from those energy - carrying processes involving cooper pair tunneling and recombination or destruction of cooper pairs .
since the phase difference between the annihilated and created pair is relevant in such a process this gives rise to the cos@xmath30 contribution to the transferred heat .
along this manuscript , we will neglect any contribution to thermal transport arising from lattice phonons .
moreover , we will assume the latter to be very well thermalized with the substrate phonons that reside at the bath temperature ( @xmath47 ) .
this is usually the case in the kind of systems of our interest as the kapitza resistance between a metallic thin film ( @xmath48 nm ) and the substrate is typically negligibly small at low temperatures @xcite
. it will be also convenient to provide with the expression for the energy exchanged between electrons and phonons in the superconductor @xcite , @xmath49 , \label{eph}\end{aligned}\ ] ] where @xmath50 $ ] , @xmath51 is the electron - phonon coupling constant , and @xmath52 is the volume of the superconducting electrode
.    the behavior of @xmath0 can be experimentally revealed by probing the electronic temperature of one of the superconducting electrodes or , alternatively , a normal metal electrode tunnel - connected to the junction .
the latter option is more convenient from an experimental point of view as electron thermometry is typically simpler to be performed on normal metals @xcite and will be therefore considered here .
we start by theoretically investigating heat exchange between two normal metal electrodes kept at different temperatures and tunnel - coupled to each other through a thermal ` modulator ' @xcite in the form of a dc - squid .
heat transport in a similar system is subsequently studied experimentally and found to be phase dependent , in agreement with the original prediction @xcite .
the structure we envision is sketched in fig . [
fig11](a ) .
it consists of a dc squid composed of two superconductors s@xmath24 and s@xmath21 in thermal equilibrium kept at different temperatures @xmath18 and @xmath19 , respectively . @xmath53 and @xmath54 denote the normal - state resistance and phase drop of jj@xmath55 , respectively . in this specific case the total heat current given in eq .
( [ heatcurrent ] ) flowing from s@xmath24 to s@xmath21 becomes @xmath56 , where @xmath57 , and @xmath58 .
for definiteness , we assume that @xmath25 so that the squid is only biased with a temperature drop across the junctions , but the voltage across them vanishes .     and
s@xmath21 kept at temperature @xmath18 and @xmath19 , respectively , are tunnel coupled so to implement a thermally - biased dc squid .
( b ) interference heat current @xmath59 vs @xmath60 calculated for a few values of @xmath61 .
( c ) total heat current @xmath62 vs @xmath60 calculated for a few values of @xmath18 at @xmath63 assuming @xmath64 .
, scaledwidth=90.0% ]    by neglecting the geometric inductance of the ring it follows that @xmath65 where @xmath60 is the applied magnetic flux through the loop , @xmath66 is an integer and @xmath67 wb is the flux quantum . for a given @xmath60 ,
the phases @xmath68 and @xmath69 are determined by the equation @xmath70 , which describes conservation of the supercurrent circulating along the loop , where @xmath71 is the ambegaokar - baratoff critical current @xcite of junction jj@xmath55 . by defining @xmath72 ( with @xmath73 ) as the degree of asymmetry of the squid junctions one gets @xcite @xmath74 and @xmath75 . using the latter equations @xmath59
can be rewritten as @xmath76 , which is analogous to the expression for the total josephson _ electric _ critical current in a dc squid with generic junctions asymmetry @xcite . in particular , for a symmetric squid ( @xmath64 ) we get @xmath77 we remark that , in this specific case , the phase - dependent component of the heat current always flows in the direction opposite to the quasiparticle heat flow , i.e. , from the cold to the hot electrode@xcite .
figure [ fig11 ] ( b ) shows @xmath59 vs @xmath60 calculated for a few values of @xmath78 at generic temperatures @xmath18 and @xmath19 such that @xmath79 .
as it can be seen , @xmath59 is a periodic function of @xmath60 maximized at integer values of @xmath80 . by
increasing @xmath78 leads to a suppression of the modulation amplitude combined with a reduction of the average value of the heat current .
eventually , the modulation amplitude is totally suppressed for @xmath81 , as only one junction is driving heat flow through the squid .
therefore high junctions symmetry is desired to maximize heat current modulation in the device .
for that reason , in the following we will restrict our calculations to the case @xmath64 .
figure [ fig11](c ) shows the total heat current @xmath62 vs @xmath60 at @xmath63 for a few values of @xmath18 .
@xmath62 is @xmath80-periodic and minimized for integer values of @xmath80 .      the experimental implementation of our heat interferometer is shown in fig . [ fig13](a ) .
the structure consists of a source and drain copper ( cu ) electrodes tunnel - coupled to a superconducting aluminum ( al ) island defining one half of a dc - squid ( s@xmath24 ) .
the other half of the squid ( s@xmath21 ) extends into a large volume lead to insure proper thermalization of its quasiparticles at the bath temperature @xmath47 .
s@xmath24 is also contacted by an extra al probe ( s@xmath82 ) via a tunnel junction , enabling independent characterization of the squid .
both source and drain are tunnel - coupled to a few external al probes so to realize normal metal - insulator - superconductor ( nis ) junctions , which allow joule heating and thermometry @xcite junctions normal - state resistances are @xmath83 k@xmath11 , @xmath84 k@xmath11 and @xmath85 k@xmath11 , respectively , whereas the resistance of each squid junction is @xmath86 k@xmath11 .
the ring area is @xmath87m@xmath88 .
finally , nis thermometers exhibit @xmath89 k@xmath11 normal - state resistance each . ] .    .
( c ) squid voltage ( @xmath90 ) versus current ( @xmath91 ) characteristics at two representative values of the applied flux and ( d ) @xmath60-dependent experimental pattern of the squid critical current @xmath92 along with the theoretical curve ( solid line ) taken at 240 mk . ]    below the critical temperature of al ( @xmath93 k ) josephson coupling allows dissipationless charge transport through the squid .
the squid voltage - current characteristics at @xmath94 mk for two representative magnetic - flux values are shown in fig . [
fig13](c ) .
in particular , a well - defined josephson current with maximum amplitude of @xmath95 na is observed at @xmath96 . the magnetic - flux pattern of the squid critical current @xmath92 along with the theoretical prediction @xcite is displayed in fig . [ fig13](d ) , and shows a nearly - complete supercurrent modulation , which confirms the good symmetry of the squid .    thermal transport and
, therefore , heat interference in the structure , arises from heating electrons in the source above lattice temperature ( @xmath47 ) so to elevate the quasiparticles temperature in s@xmath24 ( @xmath18 ) .
as the second half of the squid ( s@xmath21 ) is well thermalized at @xmath47 a temperature gradient arises across the squid .
@xmath62 will thus manifest itself leading to a @xmath80-periodic modulation of drain electron temperature ( @xmath97 .     measured for several @xmath98 values .
( b ) modulation amplitude @xmath99 and average temperature @xmath100 vs. @xmath98 .
( c ) flux - to - temperature transfer function @xmath101 versus @xmath60 measured at a few selected @xmath98 .
( d ) flux modulation of @xmath102 recorded at different @xmath47 and @xmath103 mk .
dashed lines in all panels are the theoretical results obtained from the thermal model ( see text ) . , scaledwidth=80.0% ]    investigation of heat transport in our system is performed by stabilizing first the cryostat temperature at a desired @xmath47 and heating the source up to a given temperature @xmath98 . for this purpose ,
one pair of nis junctions in the source is operated as a heater whereas a second pair is used to measure @xmath98 by applying a small dc bias current and recording the corresponding temperature - dependent voltage drop @xcite .
analogously , another pair of nis junctions is used to probe the electronic drain temperature ( @xmath102 ) against a slowly - sweeping external magnetic flux .
figure [ fig14](a ) shows @xmath102 against @xmath60 measured at 235 mk for increasing values of @xmath98 .
notably , @xmath102 is @xmath80-periodic in @xmath60 , as the josephson critical current [ see fig .
[ fig13](d ) ] , and maximized at integer values of @xmath80 .
as we shall argue , such a temperature modulation is of coherent nature , and stems from magnetic flux - control of @xmath62 , which is a hallmark of the josephson effect .
a monotonic enhancement of the average drain temperature over one flux quantum , @xmath100 , is observed when raising @xmath98 , which follows from increased heat flow across the structure . on the other hand ,
the modulation amplitude @xmath99 , defined as the difference between the maximum and minimum values of @xmath102 , turns out to initially increase and then tends to saturate at larger @xmath98 .
in particular , @xmath99 up to @xmath104 mk is observed corresponding to @xmath105 of relative modulation amplitude at 235 mk .
the full @xmath98-dependence of @xmath100 and @xmath99 are displayed in fig .
[ fig14](b ) and confirm the above described behavior . a relevant figure of merit of the heat interferometer is represented by the flux - to - temperature transfer coefficient , @xmath101 , shown in fig .
[ fig14](c ) versus @xmath60 for a few selected @xmath98 .
it turns out that @xmath106 exceeding @xmath107 mk/@xmath80 is obtained at @xmath108 mk .
the dependence on bath temperature is shown in fig .
[ fig14](d ) , which displays @xmath109 at a few increasing @xmath47 for @xmath98 set around 700 mk . besides leading to a monotonic enhancement of @xmath110 , by increasing @xmath47 yields suppression of @xmath99 and smearing of @xmath109 , which can be mainly ascribed to the enhancement of electron - phonon coupling in the drain @xcite as well as to the influence of thermal broadening .
@xmath111 mk is still observable at 450 mk , whereas the modulation disappears for @xmath112 mk .
we emphasize that the latter is substantially smaller than the temperature setting the disappearance of the josephson effect in the squid ( @xmath93 k ) .    to account for our observations
we have elaborated a thermal model sketched in fig . [ fig13](b ) .
we assume s@xmath24 to exchange heat at power @xmath113 and @xmath114 due to quasiparticle heat conduction with source and drain , respectively , at power @xmath62 with s@xmath21 and @xmath115 with s@xmath82 .
both s@xmath21 and s@xmath82 are assumed to be thermalized at @xmath47 . furthermore , drain electrons exchange energy at power @xmath114 with s@xmath24 , and at power @xmath116 with lattice phonons residing at @xmath47 @xcite .
the thermal steady - state of the system may be described by the energy - balance equations @xmath117 and @xmath118 , where first equation accounts for thermal budget in s@xmath24 , while the second one describes heat exchange in the drain . @xmath18 and @xmath102 can be determined under given conditions by numerically solving the aforementioned equations [ see eq .
( [ heatcurrent ] ) ] substituting @xmath119 for @xmath120 , @xmath121 for @xmath122 and @xmath123 for @xmath124 .
furthermore , @xmath113 can be obtained from eq .
( [ qqp ] ) by doing @xmath125 with @xmath126 and substituting @xmath119 for @xmath127 , @xmath114 can be also obtained from eq .
( [ qqp ] ) by setting @xmath128 with @xmath126 and substituting @xmath119 for @xmath129 .
for the normal metal @xmath130 , where @xmath131 wk@xmath132m@xmath133 is the electron - phonon coupling constant for cu @xcite , and @xmath134 m@xmath135 is drain volume . ] .
the model neglects heat exchange with photons due to mismatched impedance @xcite , electron - phonon coupling in s@xmath24 owing to its reduced volume and low experimental @xmath47 @xcite , as well as phonon heat current @xcite .
theoretical curves have been obtained by setting the structure parameters as extracted from the experiment , and varying @xmath120 between @xmath136 and @xmath137 to match measured data .
results from the thermal model are shown in fig .
[ fig14 ] ( dashed lines ) .
although idealized , this model provides reasonable agreement with our observations , and grasps the relevant physical picture at the origin of heat interference in our system .
in this section we envision and theoretically analyze a double - loop heat interferometer @xcite .
the latter allows to finely tune and fully balance the heat flux flowing through the device making it an improved _ heat transistor_. additionally , it is much more robust against fabrication deficiencies such as differences between the normal - state resistances of the jjs .
thanks to the application of two independent magnetic fluxes , even a quite asymmetric device is able of providing a much more robust temperature modulation compared to the simple single - loop heat interferometer or its complete annihilation .
our thermal circuit consists of two superconductors s@xmath24 and s@xmath21 , weak linked forming a double - loop interrupted by three parallel jjs [ see fig .
[ fig21](a ) ] .
@xmath138 and @xmath139 denote the normal - state resistance and phase drop of junction @xmath140 , with @xmath141 .
this structure behaves as a conventional squid pierced by a magnetic flux @xmath142 in which one of the junctions has been replaced by a dc squid .
the characteristics of this second `` junction '' can be tuned thanks to the application of a control magnetic flux @xmath143 .
the system is temperature biased by setting the temperature in s@xmath24 to be @xmath25 , @xmath19 being the temperature in s@xmath21 . furthermore , the voltage drop across the whole structure is set to zero . under these circumstances , a thermal gradient arises across the junctions and a stationary heat current @xmath144 will flow from s@xmath24 to s@xmath21 , which are in steady - state thermal equilibrium @xcite .     is risen up to @xmath145 , which yields a steady - state heat current @xmath144 .
source and drain normal - metal electrodes are tunnel - coupled to s@xmath24 .
superconducting tunnel junctions operated as heaters and thermometers are connected to source and drain .
( b ) idealized thermal diagram accounting for our setup .
arrows indicate the direction of heat flows for @xmath146 .
( c ) density plots showing @xmath147 as a function of the magnetic flux on the main and the control loop , @xmath142 and @xmath143 , respectively .
two representative cases have been considered and are schematized on the top part of each plot . ]
we shall concentrate on the phase - dependent heat current only given by @xmath148 . on the one hand , neglecting the self - induced flux in the loops ,
the fluxoid quantization on both rings imposes @xmath149 and @xmath150 , where @xmath66 and @xmath151 are integers .
the conservation of the circulating supercurrent in both loops , on the other hand , imposes @xmath152 , where @xmath153 . in writing the previous equation
we have established a given current sign convention [ see yellow arrows in fig . [ fig21](a ) ] .
furthermore we define @xmath154 , and @xmath155 .
combining the aforementioned conditions and using simple trigonometric relations one gets @xmath156 , @xmath157 and @xmath158 , where @xmath159 , @xmath160 , @xmath161 and @xmath162 $ ] .
finally , one gets @xmath163 . in fig .
[ fig21](c ) we show two density plots of @xmath147 vs. @xmath142 and @xmath143 for two representative cases . in general , the maximum of @xmath147 is always reduced for the cases in which one resistance is different from the others .
the straightest choice is @xmath164 .
although being the most simple configuration , it enables us to infer most of the characteristics of our thermal interferometer .
when analyzing the behavior of @xmath165 vs. @xmath142 for different @xmath143 values one can distinguish between two regimes ; in the first one , defined by @xmath166 [ fig .
[ fig22](a ) ] , the mean value and the shape of the curves evolves whereas the amplitude of the modulation , denoted @xmath167 , holds unchanged as we increase the amplitude of the control flux @xmath143 . at @xmath168 one
gets @xmath169 |$ ] .
this is to say , apart from a small shift equal to @xmath170 , one recovers the same dependence on @xmath142 obtained for the symmetric single - loop heat interferometer . in the second regime ,
covered by @xmath171 [ fig .
[ fig22](b ) ] , @xmath167 decreases linearly with @xmath143 whereas the mean value of @xmath165 remains constant .
furthermore , at @xmath172 , the modulation disappears completely and @xmath147 becomes independent of @xmath142 .
the aforementioned characteristics are emphasized in figs . [ fig22](c ) and [ fig22](d ) where we plot the transfer function , @xmath173 , for both regimes .
similar curves are obtained in the more general case for which @xmath174 . moreover ,
if the condition @xmath175 is satisfied , it can be demonstrated that @xmath165 can be written as a function of @xmath176 | $ ] where @xmath177 is a shift in @xmath142 .
this is to say it is still possible to suppress completely @xmath165 .
unlikely to a single - loop heat interferometer , even a quite asymmetric double - loop structure offers therefore the possibility of maximizing or suppressing completely @xmath165 through an appropriate choice of @xmath143 .
let us finally analyze how the previously described control over the magnetic flux - to - heat current transfer functions are translated into a realistic situation . for this purpose , a device similar to that of fig .
[ fig13](a ) could be envisioned . in this specific case
, the single squid loop must be substituted by a double - loop consisting of three parallel - connected josephson junctions and with independent magnetic flux controls .
similar devices operating as charge interferometers @xcite have been reported already in the literature proving the feasibility of this structure .
vs. @xmath142 plotted for different values of the control flux @xmath143 for @xmath164 .
the inset shows the case for @xmath178 .
panels ( c ) and ( d ) : transfer function @xmath173 vs. @xmath142 plotted for the same values of @xmath143 as in the top panels .
panels ( e ) and ( f ) : @xmath102 modulation calculated using the thermal model depicted in fig . [
fig21](b ) for the same conditions as in panels ( a ) and ( b ) .
panels ( e ) and ( f ) : corresponding flux - to - temperature transfer coefficient . ]
the device we envision is that schematized in fig . [
fig21](a ) .
as it was argued in section [ subsec12 ] heat transport in this structure can be studied by intentionally heating electrons in the source up to @xmath98 yielding a quasiparticle temperature @xmath179 in s@xmath24 , therefore leading to a finite heat current @xmath144 .
the latter can be modulated by means of the two control magnetic fluxes @xmath142 and @xmath143 and inferred by measuring the drain electron temperature ( @xmath102 ) .
drain temperature can be predicted by solving a couple of thermal balance equations similar to those described in section [ subsec12 ] .
new equations will differ in that no additional probe [ i.e. , @xmath180 electrode in fig . [
fig13](b ) ] is considered this time .
for this reason it will be convenient to include now the electron - phonon interaction in @xmath14 [ @xmath181 , defined in eq . ( [ eph ] ) ] . in this way
, we have now @xmath182 and @xmath118 for s@xmath24 and drain , respectively .
k@xmath11 , @xmath183 @xmath11 , @xmath184 m@xmath133 , @xmath185 wk@xmath132m@xmath133 ( typical of cu ) @xcite , @xmath186 m@xmath133 , @xmath187 wk@xmath132m@xmath133 and @xmath188ev , the last two parameters typical of aluminum ( al ) @xcite . ]
the results of thermal balance equations for drain temperature are shown in fig . [ fig22](e ) and ( f ) for @xmath164 and the same @xmath143 values as in panels ( a ) and ( b ) .
for these calculations we have set @xmath189 mk and @xmath190 mk .
as it can be seen , @xmath102 is @xmath142 periodic reaching its minimum value for @xmath191 , as expected , but the appearance of the oscillation depends strongly on @xmath143 .
we highlight that , for @xmath192 , the @xmath102 modulation is , under equal conditions , more than twice that obtained with the single thermal interferometer discussed in the previous section [ see inset in fig . [ fig22](f ) ] .
this enhancement is related with the overall enhancement of @xmath147 and the fact that @xmath102 does not depend linearly on @xmath144 . in this way , by simultaneously playing with both @xmath142 and @xmath143 , @xmath102 could be modulated from @xmath193 up to @xmath194 mk in overall .
differences between the modulation curves obtained at each @xmath143 values are emphasized in fig . [ fig22](g ) and ( h ) where the flux - to - temperature transfer coefficient @xmath195 is plotted .
in this section we theoretically analyze heat transport in temperature - biased extended jj showing that the phase - dependent component of thermal flux through the weak - link interferes in the presence of an in - plane magnetic field leading to _ heat diffraction _
@xcite , in analogy to what occurs for the josephson _ electric _ critical current @xcite . in particular ,
thermal transport is investigated in three prototypical _ electrically - open _ junctions geometries showing that the quantum phase difference across the junction undergoes @xmath196 _ jumps _ in order to minimize the josephson coupling energy .
we finally propose how to demonstrate thermal diffraction and to prove the existence of @xmath196 jumps in a realistic microstructure .
our system is schematized in fig .
[ fig31](a ) , and consists of an extended jj composed of two superconducting electrodes s@xmath24 and s@xmath21 residing at different temperatures @xmath18 and @xmath19 , respectively .
s@xmath24 and s@xmath21 are characterized by their london penetration depths @xmath197 and @xmath198 , respectively , and separated by an insulating layer with thickness @xmath199 . an external magnetic field ( @xmath200 )
is applied parallel to a symmetry axis of the junction [ see fig . [ fig31](a ) ] penetrating within a distance @xmath201 in the @xmath202 direction ) must be used , where @xmath203 @xcite . ] .
we will consider symmetric jjs in the _ short _ limit and neglect the effect of the edges .
plane in the presence of an in - plane magnetic field @xmath200 along the @xmath204 direction . additionally , source and drain normal - metal electrodes are tunnel - coupled to one of the junction electrodes ( s@xmath24 ) .
superconducting tunnel junctions operated as heaters and thermometers are connected to source and drain . @xmath205vs .
@xmath60 calculated for a rectangular [ ( a ) ] , circular [ ( b ) ] , and annular [ ( c ) ] jj represented as insets .
@xmath206 , @xmath207 , @xmath208 and @xmath78 represent the junctions geometrical parameters .
( e ) thermal model describing the main heat exchange mechanisms existing in the structure shown in ( a ) . ]    we shall focus onto the phase - dependent component only . to this end we need to determine the phase gradient @xmath209 induced by the application of @xmath200 . by choosing the closed integration contour indicated by the dashed line depicted in fig .
[ fig31](a ) it can be shown @xcite that , neglecting screening induced by the josephson current , @xmath210 , where @xmath211 and @xmath212 is the phase difference at @xmath213 .
after integration over @xmath204 the phase - dependent component of the heat current can then be written as @xmath214 , where @xmath215 is the heat current density integrated along @xmath204 .
this equation resembles the expression for the josephson current given by @xmath216 , where @xmath217 is the supercurrent density integrated along @xmath204 @xcite . in the actual configuration of electrically - open junction
, the condition of _ zero _ josephson current for any given value of @xmath200 yields the solution @xmath218 , with @xmath219 .    on the other hand , the josephson coupling energy of the junction ( @xmath220 )
can be expressed as @xmath221 where @xmath222 , @xmath6 is the zero - field critical supercurrent .
minimization of @xmath220 for any applied @xmath200 imposes the second term of @xmath223 to be always negative , so that @xmath212 will undergo a @xmath196 _ jump _ whenever the integral does contribute to @xmath220 with negative sign .
it can be shown that the @xmath196 jumps are present in any junction geometry @xcite .
they have energetic origin only .
this fact makes them very different from those present in low dimensional superconductors , caused by thermal @xcite and quantum @xcite fluctuations .
taking into account the aforementioned @xmath196 _ jumps _ and assuming that the symmetry of the junction is reflected in an analogous symmetry in the heat current , i.e. , @xmath224 , the phase - dependent component of the heat current can be written as @xmath225
.    we can now determine the behavior of @xmath205 for the three prototypical junction geometries sketched in fig .
( [ fig31 ] ) .
in particular , we shall consider two well - known examples such as the _ rectangular _ [ see fig . [ fig31](b ) ] and _ circular _ [ see fig . [ fig31](c ) ] junction , and the more exotic _
one [ see fig . [ fig31](d ) ] .
annular junctions offer the possibility to investigate fluxons dynamics due to the absence of collisions with boundaries ; yet , they provide fluxoid quantization thanks to their geometry , which allows fluxons trapping .
we assume that the total phase - dependent heat current is characterized by a uniform distribution .
@xmath205 can therefore be calculated for the three considered geometries by following , for instance , refs .
@xcite . in particular , for the rectangular junction , the absolute value of the sine cardinal function is obtained , @xmath226 . for the circular geometry one
gets the airy diffraction pattern , @xmath227 , where @xmath228 is the bessel function of the first kind , @xmath229 , and @xmath208 is the junction radius .
finally , for the annular junction @xcite , the phase - dependent component of the heat current takes the form @xmath230 , where @xmath229 , @xmath231 , @xmath232 is the @xmath66th bessel function of integer order , @xmath208 ( @xmath78 ) is the external ( internal ) radius , and @xmath233 is the number of @xmath66 trapped fluxons in the junction barrier .
right panels in fig .
[ fig31 ] illustrate the behavior of @xmath205 for the three geometries .
in particular , the curve displayed in fig .
[ fig31](a ) for the rectangular case shows the well - known fraunhofer diffraction pattern analogous to that produced by light diffraction through a rectangular slit exhibiting minima for integer multiples of @xmath80 . in the case of a circular junction
, the flux values where @xmath234 vanishes do not coincide anymore with multiples of @xmath80 , and @xmath234 falls more rapidly than in the rectangular junction case .
finally , for an annular junction , the heat current diffraction pattern is strongly @xmath66-dependent providing , apparently , enhanced flexibility to tailor the heat current response .
demonstration of diffraction of thermal currents could be achieved in the setup shown in fig . [
fig31](a ) .
this is similar to the device already analyzed in section [ subsec22 ] but the in - plane static magnetic field @xmath200 is applied perpendicular to the josephson weak - link
. the experimental operation mode will be the same and , again , the drain temperature behavior can be predicted by solving a couple of thermal balance equations accounting for the main heat exchange mechanisms existing in the structure , according to the model shown in fig .
[ fig31](e ) . , i.e. , @xmath235 k@xmath11 , @xmath183 @xmath11 , @xmath236 m@xmath133 , @xmath237 wk@xmath132m@xmath133 , @xmath238 m@xmath133 , @xmath239 wk@xmath132m@xmath133 and @xmath188ev . ]
vs @xmath60 calculated at @xmath189 mk for several values of source temperature @xmath98 for a structure based on a a rectangular [ ( a ) ] , circular [ ( c ) ] , and annular with @xmath240 [ ( e ) ] , @xmath241 [ ( g ) ] and @xmath242 [ ( i ) ] .
bottom panels show the corresponding flux - to temperature transfer function @xmath243 vs @xmath60 . for the annular case we set @xmath244 . ]
the results of thermal balance equations for drain temperature and different source temperatures are shown in fig .
[ fig33 ] .
as expected , @xmath102 shows a response to magnetic flux clearly resembling a fraunhofer or an airy - like diffraction patterns for the rectangular and circular cases , respectively . for the annular case the temperature patterns are drastically different depending on the number of fluxons trapped inside the loop . in all cases ,
temperature minima appear at the expected @xmath80 positions being the the unequivocal manifestation of the above - described phase _
jumps_.    the flux - to - temperature transfer coefficient , @xmath245 , corresponding to each of the above - mentioned cases is plotted in the bottom panels of fig .
[ fig33 ] .
as it can be seen there , @xmath243 exceeding @xmath246 mk@xmath247 could be achieved at @xmath248 mk for the annular case . moreover
, the transfer coefficient clearly demonstrates the non - monotonicity of the amplitude of drain temperature lobes as a function of @xmath98 .
in this section we propose and analyze theoretically the performance of two different thermal diodes .
one , consisting of a nis junction , and a second one , consisting of a sis josephson tunnel junction @xcite .
although never considered in the literature so far for such a purpose @xcite , superconducting tunnel junctions appear particularly well suited for the implementation of electron heat rectifiers .
heat transport in such structures is deeply influenced by the strong temperature dependence of the superconducting density of states .
in addition , the sis diode allows for the _ in - situ _ fine tuning of the thermal rectification magnitude and direction .
we shall start , first of all , by defining a heat rectification parameter @xmath249 . to this end , let us consider a tunnel jj made by two different superconductors , s@xmath24 and s@xmath21 characterized by its energy gap @xmath250 and @xmath251 leading to critical temperatures @xmath252 and @xmath253 , respectively .
the electronic temperature in both s@xmath24 and s@xmath21 is kept at fixed @xmath18 and @xmath19 , respectively , and the voltage drop across the junction is set to zero .
additionally , @xmath30 denotes the macroscopic phase difference across the junction with normal - state resistance @xmath119 . in the forward thermal bias configuration ,
a thermal gradient is created by setting @xmath254 , which leads to a total heat current @xmath255 flowing from s@xmath24 to s@xmath21 [ see fig . [
fig41](a ) ] . in the reverse thermal bias configuration ,
the thermal gradient is inverted so that @xmath256 leading to a heat current @xmath257 flowing from s@xmath21 to s@xmath24 [ see fig .
[ fig41](b ) ] . under these hypothesis
we define the rectification coefficient as @xmath258 .
we start by analyzing the case in which one of the two electrodes is a normal metal , i.e. , a nis junction . for this purpose
we can simply set @xmath259 , which leads to the complete suppression of the interference component of the heat current , @xmath260 .
this leads to @xmath261 $ ] and @xmath262 $ ] , where @xmath263 is the tunnel resistance of the nis junction .
we calculate @xmath249 as a function of @xmath18 for different values of @xmath19 , i.e. , for different temperature gradients established across the weak link . as shown in fig .
[ fig42](c ) a maximum positive rectification of @xmath264 is obtained for @xmath265 aproaching @xmath266 .
as @xmath18 increases , heat rectification starts to decrease eventually inverting its sign , which implies that heat flux from n to s becomes preferred .
furthermore , by increasing @xmath19 leads to a reduction of @xmath249 , which reaches its maximum for larger values of @xmath18 .
note that the heat rectification character of the nis diode stems from the presence of two different dos in the materials constituting the junction and the fact that the superconducting dos is strongly temperature - dependent . as a consequence
, @xmath249 depends strongly on temperature as well , inverting its sign in the vicinity of @xmath267 , i.e. , when the superconducting dos becomes energy independent .     vs. @xmath18 for different values of @xmath19 corresponding to a nis thermal diode .
black horizontal line indicates @xmath268 .
( b ) and ( c ) show the results from the thermal model described in fig . [
fig41](d ) corresponding to different fabrication parameters detailed in the legends .
black horizontal line indicates @xmath269 . ]
we conclude by analyzing the theoretical behavior of a realistic nis - based thermal rectifier . for this purpose
we consider the experimental design proposed in fig .
[ fig41](c ) .
two identical normal metal electrodes , source and drain , are weakly connected one each via resistances @xmath270 to both s and n , respectively .
superconducting probes can be tunnel - coupled to these electrodes so to implement sinis thermometers and heaters @xcite . yet
, the forward thermal bias configuration can be realized by intentionally increasing the electronic temperature in source electrode up to @xmath271 and probing the temperature in drain electrode @xmath272 . on the reverse configuration , we set @xmath273 and @xmath274
is measured in a similar way .
the difference @xmath275 for a given @xmath276 can be used to assess experimentally heat rectification .    as we have done in the previous sections , @xmath277 can be computed numerically using the thermal model described in fig . [ fig41](d ) .
thermal balance equations are given by @xmath278 , @xmath279 and @xmath280 for the @xmath281 , @xmath282 and drain electrode , respectively$ ] @xcite .
@xmath283 can be obtained in a similar way for the reverse configuration by simply exchanging the roles of @xmath284 , @xmath285 , @xmath286 and @xmath287 . as representative parameters we set @xmath288 m@xmath135 .
we choose cu for which @xmath289 wk@xmath132m@xmath133 @xcite and al @xmath290 wk@xmath132m@xmath133 @xcite . ] . using realistic parameters and assuming @xmath189 mk
, the computed values of @xmath277 vs. @xmath291 are plotted in fig . [ fig42](b ) for different @xmath263 values .
remarkably , temperature differences exceeding @xmath292 mk can be obtained .
as it can be seen , @xmath263 does not affect much the maximum value of @xmath277 .
the choice of @xmath270 , on the other hand , does influence noticeably the appearance and sign of @xmath277 as shown in fig . [
fig42](c ) .
the expected temperature differences are easily measurable with standard al - based sinis or sns thermometry techniques @xcite .
vs. @xmath18 and @xmath293 calculated for @xmath294 , @xmath295 and @xmath296 , respectively .
panels ( d ) , ( e ) and ( f ) show three selected profiles of @xmath249 vs @xmath293 for the same values of @xmath30 corresponding to the colored straight lines in ( a ) , ( b ) and ( c ) .
notice that the scale is logarithmic above the break in the vertical axis .
in addition , a dashed line indicates @xmath268 .
all curves have been calculated for @xmath297 . ]
we consider finally the case of a sis junction for which we define @xmath298 in this case , the forward and reverse total heat currents flowing through the jj read @xcite @xmath299 and @xmath300 $ ] , where the negative sign before brackets is set so that @xmath257 is positive by definition .    by fixing the temperature of the second electrode to @xmath297 ,
we calculate @xmath249 as a function of @xmath18 and as a function of @xmath293 .
the result is plotted in fig .
[ fig43](a ) , ( b ) and ( c ) for three representative cases , corresponding to @xmath294 , @xmath295 and @xmath296 , respectively .
three selected profiles of @xmath249 vs. @xmath293 for different @xmath18 values , i.e. , for different thermal gradients , are shown as well in fig .
[ fig43](d ) , ( e ) and ( f ) .
the inspection of these graphs reveals , on the first place , how phase biasing across the junction does make a substantial difference .
in particular , the heat rectification coefficient does not only change by almost two orders of magnitude from @xmath294 to @xmath296 but it also switches its sign .
it is worthwhile to emphasize that the sis junction rectifies heat only if @xmath301 . as for the case of the nis diode ,
heat rectification demands the combination of two different dos being ( at least one of them ) strongly temperature - dependent .
let us analyze in more detail the case for which @xmath296 .
in the forward configuration @xmath255 exhibits a local maximum when the temperature - dependent superconducting gaps of both superconducting electrodes coincide , which corresponds to the bright bend curve clearly visible in fig . [
fig43](c ) .
on the other hand , for the reverse configuration , @xmath257 exhibit a local minimum when the @xmath15 electrode reaches its critical temperature , which corresponds to the straight line defined by @xmath302 . as a result , @xmath249 is maximized for @xmath303 and @xmath302 .
in particular , @xmath304 is reached at @xmath305 and @xmath306 , which accounts well for the subgap leackage in realistic sis junctions .
@xmath249 is reduced to @xmath307 when considering @xmath308 and down to @xmath309 for @xmath310 . ] .
the experimental realization of the sis-based thermal rectifier is similar to that of the nis diode . in this case
, we substitute n by a second superconducting electrode s@xmath21 made , for instance , of mn - doped al since the latter allows for fine tuning of the aluminum superconducting gap @xcite .
optimum phase biasing can be achieved by using a rf squid configuration pierced by a control flux @xmath60 [ see fig . [ fig45](a ) ] .
for such a purpose , the thermal diode can be enclosed through clean contacts within a superconducting ring s@xmath82 with energy gap @xmath311 so to suppress heat losses . neglecting the inductance of the loop , the phase - flux relation is given by @xmath312 @xcite enabling the phase drop across the junction to vary within the whole phase space , i.e. , @xmath313 .
the experimental operation procedure will be that described in the previous section by @xmath314 as given in eq .
[ eph ] and setting @xmath315 [ eq . ( [ qqp ] ) ] with @xmath126 and substituting @xmath119 for @xmath129 [ see fig . [
fig45](b ) ] .
as representative parameters we set again @xmath316 m@xmath135 , @xmath317 k@xmath11 and @xmath318 k@xmath11 .
we use cu for which @xmath319 wk@xmath132m@xmath133 @xcite and al and mn - doped al with @xmath290 wk@xmath132m@xmath133 @xcite .
for these materials we set @xmath320 k and @xmath321 . ] . solving the thermal model for @xmath322 mk , the latter would give temperature difference approaching @xmath323 mk as it can be seen in fig .
[ fig45](c ) .
even more interesting , phase - coherence fingerprints are clearly observable as well .
notably , @xmath324 shows the expected @xmath325-periodicity as shown in fig . [ fig45](d ) . these temperature differences are easily measurable using sns thermometry techniques @xcite based on , for instance , vanadium with @xmath326 k.     as a function of @xmath291 and @xmath30 , respectively .
black horizontal line indicates @xmath277 . ]
along this manuscript we have proposed the experimental realization of different josephson - based devices that exploit phase - coherence of heat currents to provide with a robust magnetic flux - to - temperature control .
the successful realization of these architectures relies on the recent experimental demonstration of heat current interference in a simpler device .
the latter consist of a `` thermal '' version of a conventional electric squid that provided with an experimental temperature modulation of @xmath104 mk in amplitude upon the application of an external magnetic field as reviewed in sec .
[ sec1 ] .    in sec .
[ sec2 ] we have theoretically demonstrated how this modulation can be tuned _ in - situ _ thanks to an improved design .
considering realistic parameters , such tunability leads to a much more robust drain temperature oscillations of @xmath327 mk .
this is possible by replacing one of the squid jj by an additional dc squid leading to a double - loop thermal interferometer that allows to maximize / minimize the strength of the phase - dependent component of the heat current .
the existence of two control knobs , i.e. , two externally applied magnetic fluxes , may be exploited to perform non - trivial adiabatic cycles in its control space eventually leading to the realization of a heat pump @xcite .    in sec .
[ sec4 ] we have shown that the heat current through a temperature - biased extended jj under the influence of an in - plane magnetic field displays _ coherent diffraction _ , in full analogy with the josephson _ electric _ critical current @xcite . depending on the junction geometry
, the latter will lead to a plethora of magnetic flux - to - temperature diffraction patterns with modulation amplitudes in the range of @xmath328 mk .
the experimental realization of this device has been reported very recently @xcite setting a complementary and conclusive demonstration of the `` thermal '' josephson effect in weakly - coupled superconductors .
we have finally proposed and analyzed in sec .
[ sec5 ] the concept of a nis - based and a sis josephson thermal rectifier . under appropriate conditions , remarkably large rectification coefficients of @xmath329 and @xmath304 can be obtained for the nis and sis diodes , respectively . combining the diode within a realistic device , temperture differences between the forward and
reverse configurations of the order of @xmath330 and @xmath107 mk could be reached for the nis and sis diodes , respectively .
in addition , thanks to phase - coherence , the rectification character of the sis diode can be maximized or even inverted _ in - situ_. such a device might find a straightforward application of technological interest , e.g. , in the field of electronic refrigeration enabling magnetic - flux dependent heat management and thermal isolation at the nanoscale .
the structures described here can be integrated within , not only superconducting elements , but also hybrid mesosocopic circuits composed of , e.g. , normal metals , two dimensional electron gases and semiconductor nanowires as well . in this way
, many fields of research such as radiation detectors or quantum computing might benefit from our approach @xcite . on the other hand ,
these devices represent the first step towards the realization of coherent caloritronic circuits .
the discipline usually referred to as _
( from _ calor- _ , i.e. , `` heat '' in latin + _ electronics _ ) , deals with the generation and mastering of heat currents through an electronic drive .
yet , our approach provides with the possibility of manipulating heat currents to flow through a series of superconducting components exploiting their coherent character .
the core papers described in this review deal only with static parameter configurations where the superconducting phase and temperature are time - independent .
the natural extension is to control and change the superconducting phase in time by applying a magnetic field . in principle
this opens the possibility to transfer dynamically the heat between the superconductors , cool or heat one them and build a thermal mesoscopic engine .
this new field requires both theoretical and experimental advance and will be the subject of future research .
the fp7 program no .
228464 microkelvin , the italian ministry of defense through the pnrm project terasuper , the marie curie initial training action ( itn ) q - net 264034 and firb - futuro in ricerca 2012 under grant no .
rbfr1236vv hybridnanodev .
are acknowledged for partial financial support