A novel numerical method called the Thermal Discrete Dipole Approximation (T-DDA) isproposed for modeling near-field radiative heat transfer in three-dimensional arbitrary geometries. The T-DDA is conceptually similar to the Discrete Dipole Approximation, except that the incident field originates from thermal oscillations of dipoles. The T-DDA is described in details in the paper, and the method is tested against exact results of radiative conductance between two spheres separated by a sub-wavelength vacuum gap. For all cases considered, the results calculated from the T-DDA are in good agreement with those from the analytical solution. When considering frequency-independent dielectric functions, it is observed that the number of sub-volumes required for convergence increases as the sphere permittivity increases.Additionally, simulations performed for two silica spheres of 0.5 µm-diameter show that the resonant modes are predicted accurately via the T-DDA. For separation gaps of 0.5 µm and 0.2 µm, the relative differences between the T-DDA and the exact results are 0.35% and 6.4%, respectively, when 552 sub-volumes are used to discretize a sphere. This work suggests that the T-DDA is a robust numerical approach that can be employed for solving a wide variety of near-field thermal radiation problems in three-dimensional geometries.

The thermal discrete dipole approximation (T-DDA) is a numerical approach for modeling nearfield radiative heat transfer in complex three-dimensional geometries. In this work, the convergence of the T-DDA is investigated by comparison against the exact results for two spheres separated by a vacuum gap. The error associated with the T-DDA is reported for various sphere sizes, refractive indices and vacuum gap thicknesses. The results reveal that for a fixed number of subvolumes, the accuracy of the T-DDA degrades as the refractive index and the sphere diameter to gap ratio increase. A converging trend is observed as the number of subvolumes increases. The large computational requirements associated with increasing the number of subvolumes, and the shape error induced by large sphere diameter to gap ratios, are mitigated by using a nonuniform discretization scheme. Nonuniform discretization is shown to significantly accelerate the convergence of the T-DDA, and is thus recommended for near-field thermal radiation simulations. Errors less than 5% are obtained in 74% of the cases studied by using up to 82712 subvolumes. Additionally, the convergence analysis demonstrates that the T- † Corresponding authors. Tel.: +1 801 581 5721, Fax: +1 801 585 9825 E-mail addresses: mfrancoeur@mech.utah.edu (M. Francoeur), sheila.edalatpour@utah.edu (S. Edalatpour) 2 DDA is very accurate when dealing with surface polariton resonant modes dominating radiative heat transfer in the near field.

A fluctuational electrodynamics-based formalism for calculating near-field radiative heat transfer between objects of arbitrary size and shape and an infinite surface is presented. The surface interactions are treated analytically via Sommerfeld's theory of electric dipole radiation above an infinite plane. The volume integral equation for the electric field is discretized using the thermal discrete dipole approximation (T-DDA). The framework is verified against exact results in the sphere-surface configuration, and is applied to analyze near-field radiative heat transfer between a complex-shaped probe and an infinite plane both made of silica. It is found that when the probe tip size is approximately equal to or smaller than the gap d separating the probe and the surface, coupled localized surface phonon (LSPh)-surface phonon-polariton (SPhP) mediated heat transfer occurs. In this regime, the net spectral heat rate exhibits four resonant modes due to LSPhs along the minor axis of the probe while the net total heat rate in the near field follows a d-0.3 power law. Conversely, when the probe tip size is much larger than the separation gap d, heat transfer is mediated by SPhPs resulting in two resonant modes in the net spectral heat rate corresponding to those of a single emitting silica surface while the net total heat rate approaches a d-2 power law. It is also demonstrated that a complex-shaped probe can be approximated by a

The physics underlying spectral redshift of thermally generated surface phonon-polaritons (SPhPs) observed in near-field thermal spectroscopy is investigated. Numerically exact fluctuational electrodynamics simulations of the thermal near field emitted by a silicon carbide surface scattered in the far zone by an intrinsic silicon probe show that SPhP resonance redshift is a physical phenomenon. A maximum SPhP redshift of 19 cm -1 is predicted for a 200-nmdiameter hemispherical probing tip and a vacuum gap of 10 nm. Resonance redshift is mediated by electromagnetic gap modes excited in the vacuum gap separating the probe and the surface when the probing tip is much larger than the gap size. The impact of gap modes on the scattered field can be mitigated with a probing tip size approximately equal to or smaller than the vacuum gap. However, sharp probing tips induce important spectral broadening of the scattered field. It is also demonstrated that a dipole approximation with multiple reflections cannot be used for explaining the physics and predicting the amount of redshift in near-field thermal spectroscopy.This work shows that the scattered field in the far zone is a combination of the thermal near field † Corresponding author (S. Edalatpour).

No abstract