# Near-field radiative energy transfer between sphere-sphere and sphere-flat surfaces

### From Thermal-FluidsPedia

m |
|||

Line 25: | Line 25: | ||

- | where | + | where ''d'' = |x2-x1|, <font style="vertical-align:25%;"> <math>\rm \overline {\overline I}</math></font> is the identity tensor, and <font style="vertical-align:25%;"> <math>\rm {{\hat u} _r}{{\hat u} _r}</math></font> the dyadic notation of unit vectors. Because of thermal fluctuations, particle 1 has a random electric dipole that yields the correlation function of the dipole: |

+ | <div style="text-align: right;"> | ||

+ | <div style="float: left; width: 1115px;">[[File:Sphere4.jpg|center|EQ4]]</div> | ||

+ | (4) | ||

+ | </div> | ||

+ | |||

+ | |||

+ | |||

+ | where < > represents ensemble averaging, and * denotes the complex conjugate. Equation (4) is in fact the modified form of the ensemble average derivation provided by [[Fluctuation-dissipation_theorem|fluctuation-dissipation theorem]]. By combining the above equations, we finally obtain the heat exchange between two dipoles at temperatures ''T<sub>1</sub>'' and ''T<sub>2</sub>'', written as | ||

+ | <div style="text-align: right;"> | ||

+ | <div style="float: left; width: 1115px;">[[File:Sphere5.jpg|center|EQ5]]</div> | ||

+ | (5) | ||

+ | </div> | ||

+ | |||

+ | |||

+ | |||

+ | |||

+ | The conductance between two dipoles due to the near-field radiative heat transfer can thus be expresses as [3] | ||

+ | <div style="text-align: right;"> | ||

+ | <div style="float: left; width: 1115px;">[[File:Sphere6.jpg|center|EQ6]]</div> | ||

+ | (6) | ||

+ | </div> | ||

+ | |||

+ | |||

+ | |||

+ | |||

+ | It should be noted that radiative heat transfer between two spheres has the 1/''d''6 spatial dependence, which is typical of the dipole-dipole interactions. The thermally fluctuating dipole at one nanoparticle induces electromagnetic field on the other nanoparticle to cause the second dipole to fluctuate. Equation (6) suggests that near-field heat transfer between two nanoparticles have a similar resonant behavior at the surface polariton resonance when the dielectric constant approaches –2. In such case, the polarizability ''α'' has a resonance. Provided that the surface polariton resonance occurs when the dielectric constant approaches –1 in case the material is interfaced with the vacuum [4], this resonant behavior is not directly because of the surface polariton resonance: instead, is named as the localized surface polariton resonance – grouped oscillations of the charge density confined to nanostructures [5]. The localized surface polariton resonance appears in the visible range for metals and in the infrared for polar materials. | ||

+ | |||

+ | While the dipole approximation elucidates the ''d''<sup>-6</sup> dependence of the near-field radiative heat transfer between two spheres, this dependence is valid only when ''R << λ<sub>T</sub>'' and ''d >> R<sub>1</sub> + R<sub>2</sub>'', where ''λ<sub>T</sub>'' is the characteristic wavelength defined from Wien’s displacement law. For cases when the dipole approximation is not valid, calculation of near-field thermal radiation between two spheres becomes computationally challenging, mainly due to the difficulty in determining the dyadic Green’s function. Domingues ''et al.'' [3] attempted to address this challenge by using the molecular dynamics (MD) scheme. After computing all the atomic positions and velocities as function of time using the second Newton’s law, <math>\sum _{j} f_{ij} = m_i \ddot x _i</math>, where <big>''m<sub>i</sub>''</big> and <math>\ddot x_i</math> are the atomic mass and acceleration and <big>''f<sub>ij</sub> </big> is the interatomic force exerted by atom ''j'' on atom ''i'', the power exchange between two nanoparticles (NP1 and NP2) is computed as the net work done by a particle on the ions of the other particle: | ||

+ | <div style="text-align: right;"> | ||

+ | <div style="float: left; width: 1115px;">[[File:Sphere7.jpg|center|EQ7]]</div> | ||

+ | (7) | ||

+ | </div> | ||

+ | |||

+ | |||

+ | |||

+ | |||

+ | The interatomic force ''f<sub>ij</sub> is derived from the van Beest, Karmer, and van Santen (BKS) interaction potential [6], in which a Coulomb and a Buckingham potentials are included. On the other hand, the dyadic Green’s function for the two-sphere configurations was directly determined by using vector spherical wave functions and the recurrence relations [7]. They investigated the scattering between two spheres by expanding the electromagnetic field in terms of the vector spherical waves at each sphere and re-expanding the vector spherical waves of one sphere with the vector spherical waves of the second sphere to satisfy the boundary conditions. Reoccurrence relations for vector spherical waves are also required to reduce the computational demands in determining translation coefficients of each spherical wave function term. Due to complexities in the formulations, we do not include equations of the dyadic Green’s function for two spheres: see [7] for more details. | ||

+ | |||

+ | Near-field radiative heat transfer between two spheres computed in both Refs. [3] and [7] agree well with the dipole approximation for the interparticle distance larger than the radii of spheres. However, when the distance is smaller than the particle radius, two studies show different trends: Ref. [3] shows a higher gap dependence than ''d''<sup>-6</sup> in the thermal conductance, but Ref. [7] has the ''d''<sup>-1</sup> gap dependence as the gap distance further decreases. In Ref. [3], a deviation between the MD and the dipole approximation appears when ''d < 8R '', reaching four orders of magnitude higher than the dipole approximation. This enhanced heat transfer appears to be due to the contribution of multipolar Coulomb interactions. However, the thermal conductance in Ref. [7] asymptotically approaches a ''d''<sup>-1</sup> slope when the interparticle gap distance is much smaller than the particle radius. In fact, this slope change is consistent with the result of the proximity force approximation or the Derjaguin approximation [8], in which the flux between two curved surfaces can be locally described as a flux between two parallel plates separated by a distance ''d(r)''. For the near-field radiative heat transfer between two flat surfaces, the corresponding thermal conductance per unit area, or the radiative heat transfer coefficient ''h<sub>r</sub>'', varies with the ''d''<sup>-2</sup> dependence in the short distance regime [9-11]. When the spheres are separated by a minimum gap distance ''d<sub>0</sub>'' much smaller than the sphere radii, the thermal conductance between two spheres can be approximated as |

## Revision as of 10:03, 3 March 2012

Near-field thermal radiation between two spheres was first reported by Volokitin and Persson [1] by assuming them as dipoles. When there are two spherical nanoparticles whose dielectric constants are *ε _{1}* and

*ε*, the spectral power dissipated in particle 2 by the electromagnetic field induced by particle 1 can be written using the dipolar approximation as

_{2}

where x_{2} is the position of the particle 2 and α_{2} = 4π*R*^{3}(ε_{2} − 1) / (ε_{2} + 2) is the polarizability of a sphere of radius *R* having the relative permittivity of ε_{2}. The electric field incident on the particle 2, E_{inc}(x_{2},*ω*), is created by the thermal fluctuating dipole of particle 1 at Temperature *T _{1}* :

where is the electric dyadic Green’s function between two dipoles in vacuum and expressed as [2]

where *d* = |x2-x1|, is the identity tensor, and the dyadic notation of unit vectors. Because of thermal fluctuations, particle 1 has a random electric dipole that yields the correlation function of the dipole:

where < > represents ensemble averaging, and * denotes the complex conjugate. Equation (4) is in fact the modified form of the ensemble average derivation provided by fluctuation-dissipation theorem. By combining the above equations, we finally obtain the heat exchange between two dipoles at temperatures *T _{1}* and

*T*, written as

_{2}

The conductance between two dipoles due to the near-field radiative heat transfer can thus be expresses as [3]

It should be noted that radiative heat transfer between two spheres has the 1/*d*6 spatial dependence, which is typical of the dipole-dipole interactions. The thermally fluctuating dipole at one nanoparticle induces electromagnetic field on the other nanoparticle to cause the second dipole to fluctuate. Equation (6) suggests that near-field heat transfer between two nanoparticles have a similar resonant behavior at the surface polariton resonance when the dielectric constant approaches –2. In such case, the polarizability *α* has a resonance. Provided that the surface polariton resonance occurs when the dielectric constant approaches –1 in case the material is interfaced with the vacuum [4], this resonant behavior is not directly because of the surface polariton resonance: instead, is named as the localized surface polariton resonance – grouped oscillations of the charge density confined to nanostructures [5]. The localized surface polariton resonance appears in the visible range for metals and in the infrared for polar materials.

While the dipole approximation elucidates the *d*^{-6} dependence of the near-field radiative heat transfer between two spheres, this dependence is valid only when *R << λ _{T}* and

*d >> R*, where

_{1}+ R_{2}*λ*is the characteristic wavelength defined from Wien’s displacement law. For cases when the dipole approximation is not valid, calculation of near-field thermal radiation between two spheres becomes computationally challenging, mainly due to the difficulty in determining the dyadic Green’s function. Domingues

_{T}*et al.*[3] attempted to address this challenge by using the molecular dynamics (MD) scheme. After computing all the atomic positions and velocities as function of time using the second Newton’s law, , where

*m*and are the atomic mass and acceleration and

_{i}*f*

_{ij}*is the interatomic force exerted by atom*j

*on atom*i

*, the power exchange between two nanoparticles (NP1 and NP2) is computed as the net work done by a particle on the ions of the other particle:*

The interatomic force *f _{ij} is derived from the van Beest, Karmer, and van Santen (BKS) interaction potential [6], in which a Coulomb and a Buckingham potentials are included. On the other hand, the dyadic Green’s function for the two-sphere configurations was directly determined by using vector spherical wave functions and the recurrence relations [7]. They investigated the scattering between two spheres by expanding the electromagnetic field in terms of the vector spherical waves at each sphere and re-expanding the vector spherical waves of one sphere with the vector spherical waves of the second sphere to satisfy the boundary conditions. Reoccurrence relations for vector spherical waves are also required to reduce the computational demands in determining translation coefficients of each spherical wave function term. Due to complexities in the formulations, we do not include equations of the dyadic Green’s function for two spheres: see [7] for more details. *

Near-field radiative heat transfer between two spheres computed in both Refs. [3] and [7] agree well with the dipole approximation for the interparticle distance larger than the radii of spheres. However, when the distance is smaller than the particle radius, two studies show different trends: Ref. [3] shows a higher gap dependence than *d*^{-6} in the thermal conductance, but Ref. [7] has the *d*^{-1} gap dependence as the gap distance further decreases. In Ref. [3], a deviation between the MD and the dipole approximation appears when *d < 8R *, reaching four orders of magnitude higher than the dipole approximation. This enhanced heat transfer appears to be due to the contribution of multipolar Coulomb interactions. However, the thermal conductance in Ref. [7] asymptotically approaches a *d*^{-1} slope when the interparticle gap distance is much smaller than the particle radius. In fact, this slope change is consistent with the result of the proximity force approximation or the Derjaguin approximation [8], in which the flux between two curved surfaces can be locally described as a flux between two parallel plates separated by a distance *d(r)*. For the near-field radiative heat transfer between two flat surfaces, the corresponding thermal conductance per unit area, or the radiative heat transfer coefficient *h _{r}*, varies with the

*d*

^{-2}dependence in the short distance regime [9-11]. When the spheres are separated by a minimum gap distance

*d*much smaller than the sphere radii, the thermal conductance between two spheres can be approximated as

_{0}