# Near-field radiative transfer between two semi-infinite media

### From Thermal-FluidsPedia

Consider the structure shown in Fig. 1, where both the emitter and receiver are *n*-doped silicon. The emitter and receive are assumed to be at 400 and 300 K, respectively. The dielectric function model of doped Si may be modeled as a combination of Drude term and other contributions [1] and the details are described in [2,3]. The spectral energy transfer per unit area can be obtained as [4]

(1)

so that . Note that the integration of *s(ω,β)* over ω gives a weighted function to modify the Planck blackbody distribution function. Expression of *s(ω,β)* is different for propagating (*β < ω/c*) and evanescent (*β > ω/c*) waves,

(2a)

and,

(2b)

In Eqs. (2a,b), the first term on the right-hand side refers to the contribution of *s* polarization or TE wave, while the second term refers to the contribution of *p* polarization or TM wave. Note that and are the Fresnel reflection coefficients for *s* and *p* polarization, respectively, at the interface between vacuum and medium *j* (1 or 2). On the other hand, is the far-field reflectivity at the interface between vacuum and medium *j*. When different doping levels are considered, the location of the peak in *s(ω,β)* shifts towards higher frequencies with increased doping level.

Notice that *s(ω,β)* is independent of temperature and contains all the information about the material properties as well as the geometry of the emitting media. The predicted radiative heat transfer between two doped Si plates is plotted in Fig. 1(a) as a function of the vacuum gap width [3]. Both plates are maintained at the same doping level, which is varied from 10^{18} to 10^{21} cm^{-3}. The dotted line with circles is the radiative heat flux between two blackbodies. The net heat flux at *d* = 1 nm between 10^{19} or 10^{20} cm^{-3} doped Si plates can exceed that between two blackbodies by five orders of magnitude, because of photon tunneling and surface waves. Increase in the doping level of Si does not always enhance the energy transfer. In fact, the radiative heat transfer is the smallest for 10^{21} cm^{-3} doped Si plates as compared with other doping levels considered here. At *d* > 200 nm, doping concentrations between 10^{18} and 10^{19} cm^{-3} yield the largest radiative heat transfer, which is comparable to that between SiC and SiC. A detailed parametric study has been performed to determine the ideal Drude or Lorentz dielectric functions that yield the largest near-field enhancement [5]. Figure 1(b) illustrates the effect of doping concentration on nanoscale energy flux when the vacuum gap width is fixed at *d* = 1 nm. The doping level of medium 1 is represented as *N _{1}* while that for medium 2 is represented by

*N*. Generally speaking, surface waves are better coupled when the two media have similar dielectric functions. As a result there exist peaks when

_{2}*N*≈

_{1}*N*, at doping levels up to 10

_{2}^{20}cm

^{-3}.

The enhancement of near-field radiation in metallic and polar materials can be well explained by surface polaritons. The coupling of SP(h)Ps allows a significant increase in the function given in Eq. (2b) for evanescent waves. Furthermore, for magnetic materials, the enhancement can occur for both *s*- and *p*-polarizations, resulting in multiple spectral peaks in near-field radiative transfer [5-7]. It should be noted that for intrinsic Si or dielectric materials without strong phonon absorption bands, the tunneling is limited and saturate at extremely small distances. Also, for good metals, the SPP excitation frequency is too high that does not result significant enhancement unless the distance is less than 1 nm [5,8]. To illustrate surface wave effect, Figure 2 displays the contour plots of *s(ω,β) / 2π* for SiC plates separated at 100 nm (Figure 6a) and for 10^{20} cm^{-3} *n*-doped Si plates separated at 10 nm for TM waves since the contribution of TE waves is negligibly small. Here *β* is normalized with respect to ω/c. The brightest color represents the peak value at *ω*_{m} = 1.79×10^{14} rad/s, *β*m = 50ω/c