there is presently a great interest in the study of the radiative heat transfer between closely placed objects @xcite .
the reason for that is the experimental verification of the prediction that the radiative heat transfer between two bodies can be greatly enhanced when they are brought sufficiently close to each other @xcite .
this enhancement can lead to overcome the far - field limit set by the stefan - boltzmann law for black bodies and it takes place when the gap between the two objects is smaller than the thermal wavelength ( 9.6 @xmath0 m at room temperature ) .
such an enhanced thermal radiation exchange originates from the contribution of evanescent waves that dominate the near - field regime and by now it has been observed in numerous experiments @xcite .
these experiments have in turn triggered off the hope that near - field radiative heat transfer ( nfrht ) could lead to new applications in the context of thermal technologies such as thermophotovoltaics @xcite , heat - assisted magnetic recording @xcite , scanning thermal microscopy @xcite , nanolithography @xcite or thermal management @xcite , just to mention a few .    the theoretical description of nfrht is usually done in the framework of fluctuational electrodynamics ( fe ) introduced by rytov in the 1950s @xcite . in this framework , it is assumed that the thermal radiation is generated by random , thermally activated electric currents in the interior of a material .
these currents vanish on average , but their correlations are given by the fluctuation - dissipation theorem @xcite .
thus , the technical problem involved in the description of the radiative heat transfer between two finite objects is to find the solution of the stochastic maxwell equations with random electric current as radiation sources
. this problem can be quite challenging and an analytical solution can only be obtained in a few cases with simple geometries like , for instance , two parallel plates @xcite , two spheres @xcite or a sphere in front of a plate @xcite . in general , in order to solve this problem for complex geometries , which is usually necessary for a comparison with the experiments , one has to resort to numerical methods . in this respect
, a lot of progress has been done in recent years and standard numerical methods in electromagnetism have already been combined with fe to describe nfrht between objects of arbitrary size and shape ( for a review see ref .
those methods include , among others , the scattering matrix approach @xcite , the finite - difference time- and frequency - domain methods @xcite , the boundary element method @xcite , and volumen - integral - equation methods @xcite .
recently francoeur and coworkers adapted the discrete dipole approximation ( dda ) to describe the nfrht between two optically isotropic objects of arbitrary shape @xcite .
this new approach has been termed as thermal discrete dipole approximation ( tdda ) .
it describes the nfrht between two bodies in the framework of fe by discretizing the objects in terms of point dipoles in the spirit of the dda method that is widely used for describing the scattering and absorption of light by small particles @xcite .
the goal of this work is to generalize the tdda method to describe the radiative heat transfer and the thermal emission in magneto - optical ( mo ) systems .
mo objects are of great interest in the context of nfrht .
thus for instance , in a recent work by some of us we showed that the nfrht between two parallel plates made of doped semiconductors can be largely tuned by applying a static external magnetic field @xcite .
doped semiconductors under a magnetic field present a sizable mo activity , that can be controlled by changing the magnitude and the direction of the field .
more recently , it has been predicted that the nfrht between several mo particles , also made of doped semiconductors , can lead to the appearance of striking phenomena such a near - field thermal hall effect @xcite or the existence of a persistent directional heat current @xcite .    motivated by these recent developments in the context of nfrht of mo systems , we present here a generalization of the tdda method to describe the radiative heat transfer between mo objects of arbitrary shape , something that is still missing in the literature .
to be precise , we present a tdda method that is valid for optically anisotropic systems that can be described by an arbitrary electric permittivity tensor ( with @xmath1 ) .
we also extend the tdda approach to describe the thermal emission of a finite object ( both mo and non - mo ) .
moreover , we correct some basic errors of previous formulations of the tdda method and use our generalized approach to provide a rigorous demonstration of the kirchhoff law relating the emissivity and absorptivity of a finite mo object . finally , we also provide the correct expression for the radiative heat transfer between mo dipoles , which is the basis to study many - body effects in systems of non - reciprocal particles .
our work is focused on the analysis of homogeneous mo objects with constant temperatures , but it can be straightforwardly generalized to more complicated situations including the description of nonhomogeneous materials or situations with complicated temperature profiles .
our work provides a practical method to investigate many interesting effects in the context of the thermal emission and nfrht between mo systems .
the rest of the paper is organized as follows . in section [ sec - dda ]
we review the dda method for mo systems that we use as a starting point for our generalized tdda approach .
section [ sec - dipoles ] is devoted to the description of the radiative heat transfer between mo dipoles . in section [ sec - tdda ]
we describe in detail our generalized tdda approach for the description of the radiative heat transfer between two mo bodies . in section
[ sec - te ] we show how the tdda can be used to describe the total thermal emission of an arbitrary finite object . in section [ sec - kirchhoff ]
we address the issue of the kirchhoff law of thermal radiation of mo objects .
in particular , we derive here the definition of the directional , polarization - dependent emissivity that has to be compared with the corresponding absorption cross section to prove kirchhoff s law . in section [ sec - results ]
we present the numerical results obtained for the thermal radiation and radiative heat transfer between different objects of arbitrary size and shape .
this section illustrates both the validity and the capabilities of the formalisms developed in the previous sections .
we conclude the paper in section [ sec - conclusions ] with some additional remarks about our tdda approach and the main conclusions of our work . on the other hand ,
we have included two appendixes where we present some convergence tests and an explanation of the main differences between our tdda approach and previous formulations ( appendix a ) and an alternative derivation of the main result of section [ sec - tdda ] ( appendix b ) .
our formulation of the tdda is based on a dda extension to describe the electromagnetic response of mo systems that has been recently put forward by one of us @xcite . to make this paper more
self - contained , we present here a brief description of this method .
let us consider a nonhomogeneous finite object characterized by a spatially - dependent dielectric permittivity tensor @xmath2 that is embedded in a homogeneous medium , which hereafter is assumed to be vacuum . in the absence of currents inside the object ,
the electric field is given by the solution of the volume - integral equation @xcite @xmath3 \mathbf e(\mathbf r^{\prime } ) d\mathbf r^{\prime } .\ ] ] here , @xmath4 is source electric field or the field in the absence of the object , @xmath5 is the magnitude of the vacuum wave vector , @xmath6 is the volume of the object , and @xmath7 is the vacuum dyadic green tensor given by @xcite @xmath8 , \end{aligned}\ ] ] where @xmath9 , @xmath10 , and @xmath11 denotes the exterior product .    in the dda approach ,
the previous integral equation is solved by discretizing the volume @xmath6 as @xmath12 , where @xmath13 is the volumen of an homogeneous region where the electric field is assumed to be constant .
thus , eq .  ( [ eq - vie ] ) now reads @xmath14 \mathbf e(\mathbf r_n ) v_n .\ ] ] defining the dipole moments as @xmath15 \mathbf e(\mathbf r_n ) , \ ] ] we can rewrite eq .
( [ eq - vie - dda ] ) as @xmath16 where @xmath17 it can be shown that @xcite @xmath18 and @xmath19 here , @xmath20 is the so - called electrostatic depolarization dyadic that depends on the shape of the volume element @xmath13 @xcite . for the case of a parallelepiped of volume @xmath21
, @xmath20 adopts the form @xcite @xmath22_{ij } = \delta_{ij } \frac{2}{\pi }   \arctan \left ( \frac{1}{l^2_{n , i } } \frac{v_n}{\sqrt{l^2_{n , x } + l^2_{n , y } + l^2_{n , z } } } \right ) .
\nonumber\ ] ] thus , we can now rewrite eq .
( [ eq - vie - dda - p ] ) for the internal field , @xmath23 , as follows @xmath24 \right ] \mathbf e_n & = &   \mathbf e_{0,n } + \nonumber \\
k^2_0 \sum_{m \neq n } \hat g_{nm } [ \hat \epsilon(\mathbf r_m ) - \hat 1 ] v_m   \mathbf e_m , & & \end{aligned}\ ] ] where @xmath25 , @xmath26 and @xmath27 .    the left hand side of eq .
( [ eq - en ] ) can be defined as the exciting field @xmath28 , i.e. , the field that excites the @xmath29-volumen element .
now , defining the polarizability tensor of the @xmath29-volumen element , @xmath30 , as @xmath31 where @xmath32^{-1 }   \right)\ ] ] is the quasistatic polarizability tensor , eq .
( [ eq - en ] ) can be rewritten as a set of coupled dipole equations for the exciting fields at each element @xmath33 it is worth stressing that this dda formulation includes automatically the so - called radiative corrections @xcite , which are related to the imaginary part of the green tensor , and it is thus fully consistent with the optical theorem . on the other hand , from the solution of eq .
( [ eq - exc ] ) , which constitutes a set of 3@xmath34 coupled linear equations for the exciting fields , one can get the dipole moments and the total internal fields as follows @xmath35^{-1 }   \mathbf p_n .\end{aligned}\ ] ] let us conclude this section with a few useful remarks .
first , for cubic volumen elements , the depolarization tensor is diagonal : @xmath36 . for volumen elements of spherical or cubic , optically isotropic materials , the polarizability tensor is diagonal : @xmath37 with @xmath38 finally , from the knowledge of the dipole moments and the internal fields , one can easily compute the different cross sections ( scattering , absorption , and extinction ) .
in particular , the absorption cross section , which will play an important role later on , can be obtained as follows . assuming a plane - wave illumination , @xmath39 , the absorption cross section is given by @xcite @xmath40^{\ast } \right\ } .\ ] ]
before presenting our generalized formulation of the tdda for arbitrary bodies , it is instructive to first discuss the radiative heat transfer between magneto - optical ( mo ) dipoles .
this discussion will allow us to highlight the peculiarities of mo systems and to pinpoint some basic errors in both the existent formulation of the tdda @xcite and in the discussion of the radiative heat transfer between mo particles @xcite .
the problem that we address in this section is the radiative heat transfer ( both in the near and in the far field ) between two mo particles that are small compared to their thermal wavelengths such that they can treated as point electrical dipoles .
let us assume that these two dipoles ( or dipolar particles ) are located in positions @xmath41 and @xmath42 , they have arbitrary polarizability tensors @xmath43 and @xmath44 , and they are at temperatures @xmath45 and @xmath46 . to compute the net power exchanged between the two dipoles , we first compute the power dissipated in dipole 2 due to the emission from dipole 1 , @xmath47 , assuming that dipole 2 does not emit @xcite .
this power is given by @xmath48 where @xmath49 is the local electric current density in the volume @xmath50 and @xmath51 is the local electric field at position @xmath52 .
moreover , @xmath53 denotes the statistical average that takes into account the stochastic nature of the dipoles .
let us recall that the thermal emission originates from the fluctuating part of the dipole moments .
now , we can express @xmath54 and @xmath55 in terms of their fourier transforms @xmath56 using the fact that these two functions are real , one can easily show that @xmath57 since the fluctuation - dissipation theorem ( ftd ) will introduce a @xmath58-function of the type @xmath59 , we focus on the calculation of @xmath60 and in most cases we shall drop the argument @xmath61 to alleviate the notation .    in order to determine the dipole moment and field
, we need to solve the dda equations , see eq .
( [ eq - exc ] ) .
the exciting field at the position of dipole 2 is given by @xmath62 to close this equation , we need the corresponding equation for @xmath63 , which reads @xmath64 notice that there is no source term in this case . introducing eq .
( [ eq - exc1 ] ) in eq .
( [ eq - exc2 ] ) , we arrive at @xmath65 where @xmath66^{-1 } .\ ] ] the source field @xmath67 , due to the emission of the fluctuating part of dipole 1 , @xmath68 , is given by @xmath69 .
thus , @xmath70 where we have defined @xmath71 .
now , using eqs .
( [ eq - pn ] ) and ( [ eq - en - pn ] ) , we can express the dipole moment @xmath72 and the internal field @xmath73 in terms of @xmath74 as @xmath75^{-1 } \hat \alpha_2 \mathbf e_{{\rm exc},2 } .\end{aligned}\ ] ] to compute @xmath76 it is convenient to use a matrix notation where column vectors like @xmath72 are understood as ( @xmath77 ) matrices and row vectors like @xmath78 are ( @xmath79 ) matrices . with this notation
, we get rid of the dot product ( or scalar product ) in favor of matrix multiplications as follows @xmath80^{-1 } \hat \alpha_2   \right\ } .
& & \end{aligned}\ ] ] now , using eqs .
( [ eq - alpha ] ) and ( [ eq - alpha0 ] ) it is easy to show that @xmath81^{-1 } \hat \alpha_2 = v_2 \left [ \hat \alpha_2 - \hat \alpha^{\dagger}_2 \left ( \hat l^{\dagger}_2/v + ik^3_0/(6\pi ) \hat 1 \right ) \hat \alpha_2 \right ] .
\nonumber\ ] ] introducing this relation into eq .
( [ eq - mat ] ) and after several simple algebraic manipulations , we arrive at @xmath82 where @xmath83 notice that @xmath84 .
now , making use of eq .
( [ eq - eexc2-p1 ] ) , we obtain @xmath85 the statistical average appearing in the previous equation can be determined with the help of the fluctuation - dissipation theorem ( fdt ) , which in its most general form reads @xcite @xmath86 \hat \chi_1 , \ ] ] where @xmath87^{-1}$ ] is the bose function .
several remarks are pertinent at this stage .
first , notice that the fdt of eq .
( [ eq - fdt - p1 ] ) involves @xmath88 , which contains two terms , see eq .
( [ eq - chi ] ) .
the first one is the standard contribution , while the second one is related to the radiative correction and its contribution avoids the violation of the optical theorem @xcite .
the origin of this term has been nicely explained by messina _
et al . _
@xcite and the result above is a generalization of the arguments of messina _
et al . _
to the mo case .
it is worth stressing that the authors of ref .
[ ] in the original formulation of the tdda ignored this term , which leads to severe errors in many situations .
on the other hand , notice that in the first term in the expression of @xmath88 the combination @xmath89 appears , which for mo systems differs from @xmath90 .
this latter combination was used in the work of ref .
[ ] for the description of radiative heat transfer between anisotropic particles and therefore , that formulation is not valid for mo systems , as it was assumed in the work on the photon thermal hall effect of ref .  [ ] .
now , we can combine eqs .
( [ eq - fdt - p1 ] ) , ( [ eq - im2 ] ) , and ( [ eq - im1 ] ) to write the power dissipated in dipole 2 due to the emission from dipole 1 as follows @xmath91 k^4_0 \mbox{tr } \left\{\hat c_{21 } \hat \chi_1   \hat c^{\dagger}_{21 } \hat \chi_2 \right\ } .\ ] ] following the same line of reasoning , one can show that the power absorbed by the particle 1 due to the emission of particle 2 is given by @xmath92 k^4_0 \mbox{tr } \left\{\hat c_{12 } \hat \chi_2   \hat c^{\dagger}_{12 } \hat \chi_1 \right\ } , \ ] ] where @xmath93 can be obtained from @xmath94 by interchanging the indexes 1 and 2 .
it is straightforward to show that @xmath95 thus , the net power exchange , @xmath96 , between the two particles is given by the landauer - like formula @xmath97 { \cal t}(\omega ) , \ ] ] where @xmath98 and @xmath99 is the frequency - dependent transmission function given by @xmath100 we emphasize that this result differs from that reported in refs .  [ ] and it only coincides with those results for the case of non - mo particles .
it is also important to remark that it can be shown that in the case of isotropic particles , it reduces to the result reported in the literature in which radiative corrections have been taken into account @xcite .
in this section we present our generalized tdda for the description of the radiative heat transfer between mo objects of arbitrary shape and for arbitrary separations . in particular , we focus here on the case of two finite objects assumed to be at fixed temperatures @xmath45 and @xmath46 and do not consider the interaction with any thermal bath .    in the spirit of the dda , we assume that these two bodies are described by a collection of @xmath101 ( body 1 ) and @xmath102 ( body 2 ) electrical point dipoles .
each dipole is characterized by a volumen @xmath103 and a polarizability tensor @xmath104 , where @xmath105 indicates to which body the dipole belongs and @xmath106 if @xmath107 and @xmath108 if @xmath109 . the information about the individual dipole moments and the internal electrical field can be grouped into column supervectors ( denoted with a bar ) as follows @xmath110 @xmath111 the same notation will be used for the exciting and source electric fields .
the calculation of the net radiative power exchanged by arbitrary objects is a straightforward generalization of the calculation for the case of two dipoles presented in the previous section and we shall follow the same strategy .
thus , we first compute the power absorbed by object 2 due to the thermal emission of object 1 , @xmath47 , which is given by @xmath112 following the previous section , the goal is now to compute the exciting field @xmath113 . for this purpose , we use eq .
( [ eq - exc ] ) , which can be rewritten in a more compact form using the notation above @xmath114 here , @xmath115 with @xmath116 where @xmath117_{ij}$ ] is the vacuum dyadic green tensor connecting dipole @xmath118 in body @xmath119 with dipole @xmath120 in body @xmath121 and @xmath122 with @xmath105 .    solving now for @xmath123 , we have @xmath124 where @xmath125^{-1 } = \\   & & \left ( \begin{array}{cc } \bar 1 - k^2_0 \delta \bar g_{11 } \bar \alpha_1   & - k^2_0 \delta \bar g_{12 } \bar \alpha_2 \\
-k^2_0 \delta \bar g_{21 }    \bar \alpha_1 & \bar 1 - k^2_0 \delta \bar g_{22 } \bar \alpha_2 \end{array } \right)^{-1 } .
\nonumber\end{aligned}\ ] ] the source field , @xmath126 , due to the emission of the fluctuating dipoles in body 1 , @xmath127 , is given by @xmath128 thus , @xmath129 from eq .
( [ eq - dbar ] ) it is easy to show that @xmath130^{-1 } \delta \bar g_{12 } \bar \alpha_2 \right]^{-1 } , \nonumber \\ \bar d_{21 } & = & k^2_0 \bar d_{22 } \delta \bar g_{21 } \bar \alpha_1 \left [ \bar 1 - k^2_0 \delta   \bar g_{11 } \bar \alpha_1 \right]^{-1 } .\end{aligned}\ ] ] therefore , @xmath131 where @xmath132^{-1 } .\ ] ] now , we use eqs .
( [ eq - pn ] ) and ( [ eq - en - pn ] ) to express the dipole moment @xmath133 and the internal field @xmath134 in terms of @xmath135 as @xmath136 where @xmath137^{-1 } , \dots , [ \epsilon_0 v_{n_2,2}(\hat \epsilon_{n_2,2 } - \hat 1)]^{-1})$ ] .
following the same steps as in the previous section , it is straightforward to show that @xmath138 where @xmath139 with @xmath105 .
let us recall that @xmath140 is given by eq .
( [ eq - chi ] ) .
again , the statistical average appearing in the previous equation can be computed with the fdt , which now reads @xmath141 \bar \chi_1 .\ ] ] now , we can combine the last two equations to write the power dissipated in body 2 due to the emission from body 1 as follows @xmath91 k^4_0 \mbox{tr } \left\{\bar c_{21 } \bar \chi_1   \bar c^{\dagger}_{21 } \bar \chi_2 \right\ } .\ ] ] analogously , one can show that the power absorbed by body 1 due to the emission of body 2 is given by @xmath92 k^4_0 \mbox{tr } \left\{\bar c_{12 } \bar \chi_2   \bar c^{\dagger}_{12 } \bar \chi_1 \right\ } , \ ] ] where @xmath142 can be obtained from @xmath143 by interchanging the indexes 1 and 2 .
on the other hand , it can be shown that @xmath144 thus , the net power exchange , @xmath96 , between the two bodies is given by the landauer - like formula @xmath145 { \cal t}(\omega ) , \ ] ] where let us recall that @xmath98 and @xmath99 is the transmission function given by @xmath146 eqs .
( [ eq - pnet - bodies ] ) and ( [ eq - t - bodies ] ) are our central result for the description of the radiative heat transfer between anisotropic objects of arbitrary shape , which is valid in particular for mo systems .
obviously , the result summarized in these equations reduces to that of eqs .
( [ eq - pnet - dip ] ) and ( [ eq - t - dip ] ) for the case of two dipoles .
the result of eqs .
( [ eq - pnet - bodies ] ) and ( [ eq - t - bodies ] ) is valid that arbitrary separations between the two bodies , i.e. , it includes both the near - field and the far - field contributions . in the far - field ,
the formula can be simplified by neglecting the multiple scattering between the two bodies .
this can be done by approximating the matrix @xmath143 by @xmath147^{-1 } \delta \bar g_{21 } \left[\bar 1 - k^2_0 \bar \alpha_1 \delta \bar g_{11 } \right]^{-1 } .\ ] ] let us conclude this section by saying that our derivation of the result for the radiative heat transfer between two finite objects is based on an intuitive division of the problem into two subproblems in which one object is assumed to be the emitter and the other one the receiver , and in each subproblem the receiver is assumed to not radiate .
this is indeed the strategy followed by polder and van hove in their seminal paper in which they computed the radiative heat transfer between two parallel plates @xcite .
anyway , one might wonder if this _
ad hoc _ division is fully justified in the problem addressed in this section . to demonstrate that this is indeed the case , we present in appendix [ sec - derivation ] an alternative derivation of the central result of this section starting from a different point of view in which both objects radiate  simultaneously " and the net exchanged power
is directly calculated .
this alternative derivation detailed in appendix [ sec - derivation ] confirms the validity of the central result of this section summarized in eqs .
( [ eq - pnet - bodies ] ) and ( [ eq - t - bodies ] ) .
the thermal emission of a finite object is another important issue that can be described with the help tdda , but it has not been done so far .
the goal of this section , and the next one , is to fill this gap .
we divide our discussion of the formulation of thermal radiation within tdda into two parts .
the first one will be addressed in this section , where we present a convenient way to compute the total radiation power emitted by a finite object .
the analysis of the angular and polarization dependence of the thermal emission will be carried out in the next section , where , in particular , we generalize kirchhoff s law of thermal radiation to the case of mo objects of arbitrary size and shape .
the problem we want to address in this section is the calculation of the total radiation power emitted by a finite body at temperature @xmath148 . in this case
, we model this body as a collection of @xmath34 point dipoles that interact with the electric field of a thermal bath at temperature @xmath148 .
the total power emitted by the body , @xmath149 , must be equal to the total power absorbed from the bath when they are at the same temperature .
thus , we can follow the previous section and write the emitted power as @xmath150 where we are using the same type of notation as in the previous section , i.e. , @xmath151 in this problem , the source field is the field of the bath , i.e. , @xmath152 , where we have used the notation of previous equation .
thus , the exciting field @xmath153 is simply given by @xmath154 where @xmath155^{-1}$ ] .
here , @xmath115 , where @xmath156_{ij}$ ] is the vacuum dyadic green tensor connecting dipole @xmath118 with dipole @xmath120 inside the body and @xmath157 .    as in the previous section ,
it is straightforward to show that @xmath158 where @xmath159 .
now , to compute the statistical average appearing in the previous equation , we make use of the fdt theorem for the field correlations of the bath , which reads @xcite @xmath160   \mbox{im } \ { \bar g \ } .\ ] ] now , we can combine the last two equations to write the emitted power as follows @xcite @xmath161 where @xmath162 is planck distribution function . in the previous equation
one can identify @xmath163 as a quantity that , when divided by the geometrical cross section of the body , plays the role of an effective angular - averaged frequency - dependent emissivity . in the case of a single dipole , @xmath164 ,
which for a spherical , optically isotropic dipole ( @xmath165 ) reduces to @xmath166 $ ] , where @xmath167 is given by eq .
( [ eq - alpha - sphere ] ) .
another fundamental question that we want to address in this work is the validity of kirchhoff s law for finite mo objects .
this thermal radiation law states that the emissivity is equal to the absorptivity .
while this law is well established in the case of extended systems , it is only very recently that a rigorous proof for arbitrary finite objects has been reported @xcite . in the case of a finite object
, the absorptivity equals the absorption cross section .
this proof was provided for reciprocal objects and our goal is to extend the analysis of the validity of this law to the case of mo ( non - reciprocal ) objects .
in its most general form , kirchhoff s law involves the thermal emission in a given direction and for a given polarization .
thus , our first task is to find a proper expression for the polarization - dependent directional emissivity . for this purpose
, it is convenient to write the total emitted power in a way different than above , i.e. , in terms of poynting s vector as @xmath168 which describes the integrated flux across a differential section @xmath169 perpendicular to a radial unit vector @xmath170 , performed over a sphere of radius @xmath171 enclosing the emitting object . since this quantity does not depend on the actual choice of @xmath171 , for convenience
, we will evaluate poynting s vector in the far field .
now , we need to compute the statistical average of poynting s vector resulting from the thermal emission of an arbitrary body , i.e. , @xmath172 , where @xmath52 is the point of observation ( outside the object ) and @xmath51 and @xmath173 are the electric and magnetic field , respectively . as usual
, we can express this average in terms of the fourier transforms of the fields as @xmath174 anticipating that the fdt theorem will introduce a @xmath58-function of the type @xmath59 , we shall focus on the calculation of @xmath175 and we shall not write explicitly the argument @xmath61 .
it is worth noting that @xmath176 , in virtue of the superposition principle , can always be described in terms of two linearly polarized fields that are perpendicular to each other , both lying on a plane perpendicular to @xmath177 , as @xmath178 .    in the spirit of the tdda
, we assume that our body is described by a collection of @xmath34 point dipoles located at positions @xmath179 . using the notation of previous sections , the electric field generated by these dipoles
is given by @xmath180 where we have defined the following row supervector of green tensors @xmath181 as usual , @xmath4 is the source field that originates from the fluctuating dipoles inside that body and is given by @xmath182 the total dipole moments @xmath183 can be obtained from the exciting field inside the body as @xmath184 , while @xmath153 fulfills the dda equation @xmath185 where @xmath186 . solving for @xmath153 , we obtain that @xmath187 from which it is straightforward to show that the total field is given by @xmath188 where @xmath189 in the far - field region , the magnetic field is simply related to the electric field via @xmath190 where @xmath191 .
now , to select the contribution of a given polarization , we write the electric field as @xmath192 , where @xmath193 is a unit vector defining the polarization state and lies in a plane perpendicular to the radial direction @xmath177 . in our matrix notation ,
the field amplitude reads @xmath194 and we have @xmath195 which with the help of eq .
( [ eq - etotal ] ) can be written as @xmath196 now , we can make use of the fdt to write @xmath197 \bar \chi , \ ] ] which leads to the following expression for the radial component of poynting s vector for a polarization given by the unit vector @xmath193 @xmath198 let us stress that in the previous expression one must use the far - field expression of the green tensors , i.e. , @xmath199 where @xmath200 , @xmath201 , and @xmath202 .    taking into account the standard definition of the total emitted power @xmath203 we can write @xmath204 equation  ( [ eq - totempowerangpol ] )
can be written as @xmath205 where @xmath206 is the polarization - dependent directional emissivity given by @xmath207 let us stress that , as defined here , this emissivity has dimensions of area like the absorption cross section . on the other hand ,
it is worth noting that eq .
( [ eq - totemissionda ] ) permits establishing the following relationship between the total emissivity @xmath208 and the polarization - dependent directional emissivity @xmath206 @xmath209 where it is easy to verify that in the case of isotropic @xmath210 it reduces to @xmath211 .    the quantity @xmath206 has to be compared with the absorption cross section for the same polarization and incident angle to verify kirchhoff s law .
the expression for the absorption cross section is given by eq .
( [ eq - cabs1 ] ) , which can be written as @xmath212 where @xmath213 is a linearly polarized plane wave @xmath214 impinging as @xmath215 .    a general demonstration of kirchhoff s law from eqs .
( [ eq - zeta - emissivity ] ) and ( [ eq - cabs2 ] ) is cumbersome .
in what follows we will provide evidence of this fulfillment in the case of single dipoles whereas a numerical verification for finite objects will be provided in section [ sec - results ] .
in the case of a single dipole we have that @xmath216 in eq .
( [ eq - zeta - emissivity ] ) and @xmath217 in eq .
( [ eq - cabs2 ] ) . then , using @xmath218 , the absorption cross section can be written as @xmath219 .
on the other hand , if we assume that the dipole is located at the origin of the spherical coordinate system , the green tensor in far - field approximation is given by @xmath220 and eq .
( [ eq - zeta - emissivity ] ) becomes @xmath221 finally , for an isotropic system this leads to the well - known @xmath222 kirchhoff s result .
in this section we present a series of results related to the different aspects of radiative heat transfer and thermal emission discussed in previous sections .
our goal is to illustrate the capabilities of the tdda method put forward in this work , as well as to demonstrate its validity by comparing with established results in the literature . for these purposes
, we shall consider below two materials .
first , as an example of a non - mo material we choose sio@xmath223 , which is polar dielectric that has been amply studied in the context of nfrht .
we use the dielectric function of sio@xmath223 reported in palik s book @xcite . as an example of a mo material we shall consider @xmath29-doped insb , a polar semiconductor , that when subjected to an external magnetic field becomes mo . for the sake of concreteness
, we shall assume that magnetic field is applied in the @xmath224-direction , @xmath225 . in this case
, the permittivity tensor of insb adopts the following form @xcite @xmath226 where @xmath227 } \right ) , \nonumber \\ \label{eq - epsilons } \epsilon_2(h ) & = & \frac{\epsilon_{\infty } \omega^2_p \omega_c}{\omega [ ( \omega + i \gamma)^2 - \omega^2_c ] } , \\
\epsilon_3 & = & \epsilon_{\infty } \left ( 1 + \frac{\omega^2_l - \omega^2_t}{\omega^2_t - \omega^2 - i \gamma \omega } - \frac{\omega^2_p}{\omega ( \omega + i \gamma ) } \right ) .
\nonumber\end{aligned}\ ] ] here , @xmath228 is the high - frequency dielectric constant , @xmath229 is the longitudinal optical phonon frequency , @xmath230 is the transverse optical phonon frequency , @xmath231 defines the plasma frequency of free carriers of density @xmath29 and effective mass @xmath232 , @xmath233 is the phonon damping constant , and @xmath234 is the free - carrier damping constant . finally , the magnetic field enters in these expressions via the cyclotron frequency @xmath235 .
let us point out that in this case the magneto - optics is induced by the magnetic field , which modifies the diagonal elements of the permittivity tensor and introduces off - diagonal terms .
there are two contributions to the diagonal components , namely , optical phonons and free carriers .
these latter ones are responsible for the magneto - optics at finite magnetic field .
notice that we have neglected the contribution from intra - band transitions because we are interested in thermal properties at room temperature , where this contribution does not play any role . in all calculations below
, we consider the particular case taken from ref .
[ ] , where @xmath236 , @xmath237 rad / s , @xmath238 rad / s , @xmath239 rad / s , @xmath240 rad / s , @xmath241 @xmath242 , @xmath243 , and @xmath244 rad / s . as a reference ,
let us say that with these parameters @xmath245 rad / s for a field of 1 t.    the explicit form for the dielectric tensor given by eq .
( [ perm - tensor - hz ] ) gives rise to a @xmath246 , see eq .
( [ eq - chi ] ) , having the form @xmath247 for simplicity , the results presented below are for homogeneous bodies , i.e. , with a spatially constant permittivity tensor , and for constant temperatures inside the bodies
. moreover , all the results for finite systems were obtained by discretizing the objects in terms of cubes of equal size .
the corresponding polarizability tensors of the cubes were computed with eqs .
( [ eq - alpha ] ) and ( [ eq - alpha0 ] ) with @xmath36 , which in the case of isotropic materials reduce to eq .
( [ eq - alpha - sphere ] ) .
unless stated otherwise , the calculations of the radiative thermal conductance , power emission , etc . , were carried out at @xmath248 k.      the issue that we want to address now is the description of the radiative heat transfer between two finite objects using the formalism detailed in section [ sec - tdda ] . to avoid known problems related to the so - called shape error , an error due to the description of objects that can not be exactly represented by cubical lattices @xcite , we consider the heat transfer between two identical cubes of side @xmath249 . in fig .
[ fig - nfrht - comparison1 ] we show as solid lines the tdda results for the spectral radiative thermal conductance ( conductance per unit of energy ) as a function of the photon energy for two sio@xmath223 cubes of side @xmath250 @xmath0 m and various separations ranging from 100 nm to 5 @xmath0 m . to demonstrate the validity of our tdda approach , we compare our results with those obtained with the code scuff - em @xcite , which are shown in fig .
[ fig - nfrht - comparison1 ] as open symbols .
the scuff - em solver implements the fluctuating - surface - current ( fsc ) formulation of the heat transfer problem put forward in refs .
[ ] in combination with the boundary element method ( bem ) .
the fsc - bem combination used in scuff - em enables the description of the radiative heat transfer between homogeneous bodies of arbitrary shape and provide numerically exact results within the framework of fluctuational electrodynamics in the local approximation , i.e. , assuming that the dielectric function only depends on frequency .
cubes of side @xmath250 @xmath0 m and various gaps , @xmath251 , at @xmath248 k. the solid lines correspond to the results computed with the tdda approach with 6859 dipoles per cube , while the open symbols correspond to the results obtained with scuff - em . ]    as it can be seen in fig .
[ fig - nfrht - comparison1 ] , our tdda method is able to reproduce the numerically exact results obtained with scuff - em . let us say that the tdda results shown here were obtained modeling each cube with 6859 cubic dipoles ( the analysis of the convergence of these results with the number of dipoles used in the simulations is provided in appendix  [ sec - comparison ] ) .
it is worth stressing that it becomes progressively more demanding to converge the tdda results upon reducing the gap between the cubes .
thus for instance , in this example the difference between the tdda and the scuff - em results for the total conductance ( integrated over energy ) is of 2.5% for a gap of @xmath252 @xmath0 m , while it monotonically increases up to 6.2% for @xmath253 nm . the physical reason for this behavior is the fact that the nfrht in this case is dominated by surface phonon polaritons @xcite that have a penetration depth comparable to the gap size @xcite .
this implies that the electric field inside the cubes varies on a length scale smaller than the gap and therefore , one needs an increasing number of dipoles to properly describe the radiative heat transfer as the gap diminishes .
the important thing is that the error in the tdda calculations can be systematically reduced by refining the dipole grid , which can be done by increasing the number of dipoles or using adaptive meshes .
let us stress that in this work we are not interested in presenting a detailed analysis of the convergence of our tdda approach , but rather in establishing its fundamental validity .
we present another example of the comparison between our tdda results and those obtained with scuff - em for the spectral conductance of two sio@xmath223 cubes of varying size and gap @xmath254 @xmath0 m in fig .
[ fig - nfrht - comparison2 ] .
again , we find a very good agreement between both types of results , which becomes progressively worse ( for a fixed number of dipoles ) as the side of the cube increases , as expected .
more importantly , we always find that this agreement can be systematically improved by increasing the number of dipoles in our cubic lattices ( for more details on the convergence of these results , see appendix  [ sec - comparison ] ) .
so in short , we conclude that our tdda method produces the correct results for the radiative heat transfer provided that a sufficiently large number of dipoles is employed in the discretization of the objects .
cubes of different side values and separated by a gap of 0.5 @xmath0 m at @xmath248 k. the solid lines correspond to the results computed with the tdda approach with 6859 dipoles per cube , while the open symbols correspond to the results obtained with scuff - em .
the error in the tdda calculations in the total conductance is of 3.1% for @xmath250 @xmath0 m and it increases monotonically up to 11.3% for @xmath255 @xmath0 m . ]     @xmath0 m separated by a gap of 500 nm , at @xmath248 k , and for various values of the magnetic field , @xmath256 , applied along the @xmath224-direction . the inset shows the discretization geometry , where the number of dipoles per cube is 4913 ( each one has an edge of 59 nm ) . ]    we now proceed to illustrate our approach in the case of mo objects . for this purpose ,
we consider the radiative heat transfer between two identical insb cubic particles of side @xmath257 @xmath0 m separated @xmath258 nm using 4913 cubic dipoles for each particle .
the results for the spectral conductance for different values of the external magnetic field , @xmath256 , are shown in fig .
[ fig - rht_2cubesinsb ] . in this case
, the magnetic field is oriented as shown in the inset of this figure .
as one can see , the spectral conductance in the absence of field ( @xmath259 ) is dominated by the contribution of two peaks that , as explained in ref .
@xcite , are due to the contribution of surface plasmon polaritons ( low - energy peak ) and surface phonon polaritons ( high - energy peak ) .
notice that the magnetic field induces a splitting of the low - energy peak , while it introduces a new peak at high energies that appears at the cyclotron frequency and therefore , its position is blue - shifted linearly with the external field .
this behavior is very similar to what it was found in ref .
@xcite for the case of two semi - infinite insb plates .
these results will be discussed in more detail in a forthcoming publication and we show them here just to illustrate the capabilities of our tdda method .    to conclude this subsection , we want to illustrate the validity of the far - field approximation of eq .
( [ eq - ffapprox ] ) . for this purpose ,
we show in fig .
[ fig - ff - sio2 ] a comparison between the exact results and the far - field approximation for the room temperature spectral conductance for two identical cubes of sio@xmath223 of various sizes separated a distance of 20 @xmath0 m , which is larger than the thermal wavelength .
the results were computed with 2197 dipoles in each cube , which was sufficient to converge the results with an accuracy of better than 1% . as it is clear from fig .
[ fig - ff - sio2 ] , the approximation of eq .
( [ eq - ffapprox ] ) accurately reproduces the exact results .
cubes with a gap of 20 @xmath0 m , @xmath248 k , and for various values of the cube side , @xmath249 .
the solid lines correspond to the exact results and the dashed lines to the results computed with the far - field approximation of eq .
( [ eq - ffapprox ] ) . ]      in this subsection we want to discuss the total thermal emission of a finite object . to test the validity of the formalism described in section [ sec - te ] we first consider the total thermal emission of a sphere of a non - mo material .
a well - known result , first derived by kattawar and eisner @xcite , is that the emissivity of a non - mo sphere of arbitrary radius obtained within fluctuational electrodynamics is equal to the corresponding absorption efficiency , i.e. , to the absorption cross section normalized by the geometrical cross section .
let us recall that we demonstrated this result in section [ sec - te ] for the case of a single spherical dipole .
since the absorption efficiency can be computed exactly with the help of mie theory @xcite , the result of kattawar and eisner provides a very stringent test for our theory of thermal emission .
we have computed both the emissivity and absorption efficiency using eqs .
( [ eq - emissivity ] ) and ( [ eq - cabs2 ] ) , respectively , for spheres of different sizes and materials and in all cases we have found that they are indeed identical . in fig .
[ fig - emissivity - sio2 ] we show the tdda results for the total emissivity of a sio@xmath223 sphere for different radii ( solid lines ) as a function of the energy .
let us clarify that we are showing the quantity defined in eq .
( [ eq - emissivity ] ) and normalized by the geometrical cross section to make it dimensionless .
we also show the corresponding results for the absorption efficiency calculated with the exact mie theory @xcite .
as one can see , there is a very good agreement between our tdda results and the exact results , which demonstrates the validity of our formalism .
again , it is worth remarking that as one increases the size of the sphere , a larger number of dipoles is required to satisfactorily converge the results . to obtain the results of fig .
[ fig - emissivity - sio2 ] we used about 12000 dipoles , which is enough to reproduce the exact results for spheres of radius up to approximately 10 @xmath0 m . in any case , it is important to emphasize there is no fundamental limitation to obtain the exact result within our tdda approach , although the practical problem to converge the results may become in some cases very difficult to overcome .
sphere as a function of the energy for different values of its radius @xmath171 at @xmath248 k. the solid lines correspond to the results computed with tdda with @xmath260 dipoles , while the open symbols correspond to the results obtained with mie theory . ]    for completeness , we show in fig .
[ fig - total - emission - sio2 ] the total power emitted at room temperature by a sphere and a cube of sio@xmath223 of varying size .
the results are normalized by the corresponding results for a black body ( stefan - boltzmann law ) .
as one can see , for small spheres and cubes the emitted power is proportional to the volumen of the finite objects .
this is a well - known result @xcite , which is due to the fact that in this regime the skin depth of this material at the relevant wavelengths is larger than the characteristic size of the object .
this means in practice that the whole object contributes to the thermal emission .
however , as the size increases , the emitted power becomes proportional the area of the object , which reflects the fact that the size becomes larger than the skin depth and then only the surface of the object significantly contributes to the thermal emission .
notice also that in both cases ( spheres and cubes ) , when the characteristic size becomes on the order of a few microns , the emitted power becomes comparable to that of a black body of the same size .
let us say that the thermal emission of a sio@xmath223 sphere was studied by krger _
et al . _
@xcite and our results agree with those reported in those references .
sphere as a function of its radius @xmath171 at @xmath248 k. the power is normalized with the black body result , @xmath261 , where @xmath262 is the total area of the sphere .
( b ) the same as in panel ( a ) , but for a sio@xmath223 cube of side @xmath249 . ]
the goal of this subsection is to illustrate the validity of kirchhoff s law in the case of finite mo objects .
but before presenting the full numerical analysis , it is instructive to consider a mo dipole ( where we still consider mo activity is due to a magnetic field pointing along the @xmath224-direction ) since it is possible to provide an analytical solution for an arbitrary polarization vector @xmath193 . using the coordinate system depicted in fig .
[ fig - coordinatesys ] the spectral , polarization - dependent directional emissivity for a single dipole with @xmath88 given by eq .
( [ chi - tensor - hz ] ) can be written for any polarization angle @xmath263 as [ see eq .
( [ eq - kirchhoff1dip ] ) ] @xmath264.\end{aligned}\ ] ] notice that there is no dependence on the azimuthal angle @xmath265 due to the symmetry introduced by the fact that @xmath266 .
moreover , since we have to choose two orthogonal vectors , i.e. , @xmath263 and @xmath267 , the sum required to obtain the total emitted power is equal to @xmath268 $ ] , that , as expected , does not depend on the choice of @xmath263 .
since the direction of the external magnetic field imposes a preferential direction , it is convenient to choose two values of @xmath263 so that one polarization vector has no component along the direction of the external magnetic field , @xmath269 rad ( let us call it @xmath270 ) , while the other is collinear with the magnetic field , @xmath271 rad ( let us call it @xmath272 ) .
it is straightforward to see that @xmath273 , i.e. , independent of any incidence angle @xmath274 or @xmath265 , while @xmath275 $ ] presenting thus a @xmath276 dependence , varying from @xmath277 to @xmath278 .    to illustrate the fulfillment of kirchhoff s law for a body of finite size
, we consider an homogeneous insb sphere of radius @xmath279 @xmath0 m , using discretization unit cubes with an edge of @xmath280 nm ( @xmath281 2550 discretization elements ) . in fig .
[ fig - kirchhofflaw](a ) we present spectra of the polarized emissivity and of the absorption cross section for a direction defined by @xmath282 rad and polarization angle by @xmath283 rad .
two values of magnetic field , @xmath259 t and @xmath284 t , are shown .
for both cases @xmath285 and @xmath286 ( for the same angles and polarization ) are identical .    in order to analyze the same system as a function of the geometrical angles ,
we verify the @xmath287 dependence obtained for a single dipole in eq .
( [ eq - emissivmodipolesimple ] ) .
this is clearly visible in fig .
[ fig - kirchhofflaw](b ) , where we plot @xmath288 and @xmath289 for @xmath290 rad , @xmath291 rad and @xmath292 ev , as a function of @xmath263 for an applied field of @xmath284 t normalized to its value at @xmath259 t ( let us point out that , since the geometry is fully isotropic , for @xmath259 t it happens that @xmath293 and thus emissivity and absorption cross section are polarization insensitive ) .
the next step in the verification is to compare the angular profile of the emissivity and of the absorption cross section for fixed polarizations . as it can be understood from eq .
( [ eq - emissivmodipolesimple ] ) , the most interesting choice is the @xmath272-polarization ( @xmath271 rad ) . in fig .
[ fig - kirchhofflaw](c ) , we present a colormap of @xmath294 and @xmath289 ( for the same polarization ) for @xmath284 t normalized to @xmath259 t , as a function of @xmath265 and @xmath274 showing that they are independent of @xmath265 and that the dependence on @xmath274 is of the form of @xmath276 , as evidenced by the associated polar plot . the same plot , but for @xmath270-polarization is presented in fig .
[ fig - kirchhofflaw](d ) demonstrating that for this configuration and polarization , the emissivity and absorption cross section is isotropic , within a numerical error smaller than 2% .
let us conclude this subsection by saying that we have analyzed the polarization - dependent , directional emissivity and the corresponding absorption cross section for different sizes and shapes of insb particles under a magnetic field finding that in all cases these two quantities are identical .
in this work we have focused on the description of the thermal radiation and radiative heat transfer of homogeneous objects with constant temperatures .
however , all the basic formulas derived in this work are directly applicable to the case of inhomogeneous objects , with spatially dependent dielectric functions , and they can be straightforwardly generalized to consider the case of arbitrary temperature profiles in the interior of the objects . on the other hand , the tdda method presented here can easily be generalized in a number of additional ways to treat , for instance , the radiative heat transfer between surfaces and finite objects @xcite or between periodic systems @xcite .    with respect to the improvement of the efficiency and accuracy of the tdda method presented here , there are a number of obvious improvements that one can implement without the need to modify the formulation detailed in this work .
thus for instance , one can employ adaptive dipole lattices to better describe the non - uniform field profiles inside the objects @xcite , one can discretize the objects in terms of non - cubic dipoles with shapes better adapted to the geometries of the objects under study @xcite , or one can use a dda formulation based on the integration of the green s tensor ( igt ) in the spirit of ref .
@xcite , rather than on the standard point - dipole interaction , as done in this work .    in summary ,
we have presented in this work a new formulation of the tdda approach to describe the radiative heat between finite objects of arbitrary size and shape .
this formulation corrects the deficiencies of previous tdda methods and allows us to describe the radiative heat transfer between mo objects and , more generally , between optically anisotropic objects with arbitrary permittivity tensors .
we have shown how this tdda approach can also be used to describe the thermal emission of a finite object . moreover ,
we have provided very compact and transparent formulas for different physical quantities that can be directly employed in generalizations of the method presented here . in particular , we have corrected the existent formulas for the radiative heat transfer between non - reciprocal dipoles , which can be used to analyze different many - body effects in ensembles of mo particles . on the other hand ,
we have used our tdda approach to demonstrate kirchhoff s law relating the emission and absorption of non - reciprocal objects .
our work opens the way for the rigorous description of different thermal radiation phenomena involving finite mo objects of arbitrary shape .
we thank v. fernndez - hurtado for providing the numerical results obtained scuff - em shown in this paper , and j.j .
senz and n. de sousa for stimulating discussions .
this work was financially supported by the comunidad de madrid ( contract no .
s2013/mit-2740 ) .
j.c.c .  also acknowledges financial support from the spanish mineco ( contract no .  fis2014 - 53488-p ) and thanks the dfg and sfb 767 for sponsoring his stay at the university of konstanz as mercator fellow .
m . also acknowledges funding from the spanish mineco ( contract no .
mat2014 - 58860-p ) .
in this appendix we present a brief analysis of the convergence of the results shown in section [ sec - results - rht ] as well as a comparison of our approach with previous formulations of the tdda @xcite .    in fig .
[ fig - convergence ] we show some convergence tests for the radiative heat transfer between two sio@xmath223 cubes of two different sizes . in this figure , together with the exact numerical results obtained with scuff - em , we show the tdda results obtained with different number of dipoles ( per cube ) .
as one can see , upon increasing the number of dipoles , the tdda results first converge rapidly to approximately the correct result and then , the convergence improves slowly , but monotonically until a very satisfactory agreement with the scuff - em results is achieved . of course
, as the size of the objects increases , the required number of dipoles increases accordingly and a proper convergence becomes much more demanding .
the same occurs when the gap size is reduced ( not shown here ) . in that case
, the electric field inside the sio@xmath223 particles varies more rapidly in space due to the smaller penetration depth of the surface phonon polaritons that dominate the nfrht in this polar dielectric @xcite .
cubes of side @xmath250 @xmath0 m and a gap size of @xmath295 @xmath0 m .
the different lines correspond to the results computed with the tdda approach with different number of dipoles per cube , as indicated in the legend .
the open symbols correspond to the converged result obtained with scuff - em .
( b ) the same as in panel ( a ) , but for two sio@xmath223 cubes of side @xmath296 @xmath0 m . ]    let us discuss now the relationship with the formalism presented in ref .
we would like to focus on eqs .
( 3.16 ) and ( 3.22 ) in ref .
@xcite and its consequences .
those equations indicate that the thermally generated incident field is nil for all elemental dipoles in the emitter , and are only affected by the multiple scattering with the absorbing object ( or receiver ) .
this is at variance with our eq .
( [ eq - excitingfield ] ) where all the dipoles in the emitter feel the thermally activated field of all the constituent dipoles but itself and so there is also multiple scattering inside the emitter .
this multiple scattering in the emitting element accounts for the fact that the thermal wave activated deep inside the emitter has to first travel across it to reach the receiver . for small objects , neglecting that multiple scattering results in small deviations from the correct result , precisely due to such small size .
however , for moderately large objects such a multiple scattering plays a more important role , especially in the frequency regions in which absorption is particularly large due to a resonant behavior .
this is illustrated in fig .
[ fig - differences ] where we compare our results for the radiative heat transfer between two sio@xmath223 cubes of side 2 @xmath0 m separated by a distance of @xmath295 @xmath0 m with those obtained ignoring the multiple scattering in the emitter . in particular , in this comparison
we have made sure that both the results based on the approach of ref .
@xcite and our results are properly converged , see fig .
[ fig - differences ] and the corresponding caption .
thus , ignoring this multiple scattering in the emitter prevents recovering the correct result and the error becomes progressively more severe as the size of the particles increases .    on the other hand , we would like to recall that another important difference between our tdda approach and the formulations of refs .
@xcite resides in the form of the fluctuation - dissipation theorem employed in those references .
as one can see in eq .
( [ eq - fdt - p1 ] ) , in our case this theorem involves the quantity @xmath88 , which is given by eq .
( [ eq - chi ] ) .
this quantity contains the so - called radiative corrections , whose absence , as in the expressions of refs .
@xcite , leads to wrong results since the optical theorem is not fulfilled ( energy is not conserved ) .
our result in eq .
( [ eq - chi ] ) is a generalization to the mo ( tensorial ) case , which is valid for non - reciprocal systems .
cube of @xmath296 @xmath0 m at @xmath297 k and a gap distance of @xmath295 @xmath0 m , for a number of elemental dipoles ensuring that convergency is achieved .
the solid lines are the results of eqs .
( 3.16 ) and ( 3.22 ) in ref .
@xcite and the open circles are the converged result of our formalism . ]
given the discrepancies found with the original formulation of tdda of ref .
[ ] , it is reassuring to derive the central result of eq .
( [ eq - pnet - bodies ] ) in an alternative way . for this purpose ,
we have made use of the formalism put forward by messina _
et al . _
@xcite for the description of heat transfer in systems of multiple dipoles .
this formalism was developed to study many - body effects in the heat transfer between dipoles , but as we show below , it can be adapted in the tdda spirit to describe the heat transfer between bodies of arbitrary size .    as in section [ sec - tdda ] , we consider two bodies described by a collection of @xmath101 ( body 1 ) and @xmath102 ( body 2 ) electrical point dipoles .
again , each dipole has a volumen @xmath103 and it is characterized by a polarizability tensor @xmath104 , where @xmath105 . in this case , rather than separating the problem into two problems of the emission of one body and the absorption of the other , as we did in section [ sec - tdda ] , we now consider that the fluctuating dipoles in both bodies are radiating at the same time . in this case
, the total dipole moments can be grouped into column supervectors as follows @xmath298 where there are two contributions to the total dipole moments , one coming from the fluctuating dipoles , @xmath299 , and the other is an induced contribution arising from the interaction with the other dipoles .
this second contribution is given by @xcite @xmath300 here , the indexes @xmath118 and @xmath120 refer to dipoles in both bodies . from this equation
it is easy to show that the total dipole moments are given by @xmath301 here , we have used the notation of section [ sec - tdda ] .    now
the goal is compute the net power balance in , let us say , body 2 .
this net power is given by @xmath302 the dipole moment @xmath133 can be obtained from eq .
( [ eq - t - alt ] ) : @xmath303 while @xmath134 can be related to @xmath133 via eq .
( [ eq - en - pn ] ) . then , with the usual algebraic manipulations , one can show that @xmath304 which with the help of the fdt theorem leads to @xmath305 \nonumber .\end{aligned}\ ] ] in thermal equilibrium ( @xmath306 ) , @xmath307 should vanish .
therefore , they following relation must hold @xmath308 thus , @xmath307 can be rewritten as @xmath309 { \cal t}(\omega ) , \ ] ] where @xmath310 now , the remaining task is to show that the previous expression for the transmission reduces to eq .
( [ eq - t - bodies ] ) . for this purpose
, one can use eq .
( [ eq - t - alt ] ) to show that @xmath311^{-1 } \bar \alpha_1 \delta \bar g_{12 } \right]^{-1 } , \nonumber \\
\bar t^{-1}_{21 } & = & k^2_0 \bar t^{-1}_{22 } \bar \alpha_2 \delta \bar g_{21 }   \left [ \bar 1 - k^2_0 \bar \alpha_1 \delta \bar g_{11 } \right]^{-1 } .\end{aligned}\ ] ] from these equations it is straightforward to show that @xmath312 from which it is obvious that eq .
( [ eq - t - new ] ) reduces to eq .
( [ eq - t - bodies ] ) .
this completes our alternative derivation of eq .
( [ eq - pnet - bodies ] ) , which confirms the validity of our result .
w.  a. challener , c.  b. peng , a.  v. itagi , d. karns , w. peng , y.  y. peng , x.  m. yang , x.  b. zhu , n.  j. gokemeijer , y.  t. hsia , g. ju , r.  e. rottmayer , m.  a. seigler , and e.  c. gage , nat . photon . * 3 * , 220 ( 2009 ) .    b.  c. stipe , t.  c. strand , c.  c. poon , h. balamane , t.  d. boone , j.  a. katine , j.  l. li , v. rawat , h. nemoto , a. hirotsune , o. hellwig , r. ruiz , e. dobisz , d.  s. kercher , n. robertson , t.  r. albrecht , and b.  d. terris , nat . photon . * 4 * , 484 ( 2010 ) .