## Abstract

Thermal radiative phenomena can be strongly influenced by the coupling of phonons and long-range electromagnetic fields at infrared frequencies. Typically employed macroscopic descriptions of thermal fluctuations often ignore atomistic effects that become relevant at nanometric scales, whereas purely microscopic treatments ignore long-range, geometry-dependent electromagnetic effects. We describe a mesoscopic framework for modeling thermal fluctuation phenomena among molecules near macroscopic bodies, conjoining atomistic treatments of electronic and vibrational fluctuations obtained from density functional theory in the former with continuum descriptions of electromagnetic scattering in the latter. The interplay of these effects becomes particularly important at mesoscopic scales, where phonon polaritons can be strongly influenced by the objects' finite sizes, shapes, and nonlocal or many-body response to electromagnetic fluctuations. We show that, even in small but especially in elongated low-dimensional molecules, such effects can modify thermal emission and heat transfer by orders of magnitude and produce qualitatively different behavior compared to predictions based on local, dipolar, or pairwise approximations.

Original language | English (US) |
---|---|

Article number | 045901 |

Journal | Physical review letters |

Volume | 121 |

Issue number | 4 |

DOIs | |

State | Published - Jul 27 2018 |

## All Science Journal Classification (ASJC) codes

- Physics and Astronomy(all)

## Access to Document

## Other files and links

## Fingerprint Dive into the research topics of 'Phonon-Polariton Mediated Thermal Radiation and Heat Transfer among Molecules and Macroscopic Bodies: Nonlocal Electromagnetic Response at Mesoscopic Scales'. Together they form a unique fingerprint.

## Cite this

- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver

*Physical review letters*,

*121*(4), [045901]. https://doi.org/10.1103/PhysRevLett.121.045901

**Phonon-Polariton Mediated Thermal Radiation and Heat Transfer among Molecules and Macroscopic Bodies : Nonlocal Electromagnetic Response at Mesoscopic Scales**. In: Physical review letters. 2018 ; Vol. 121, No. 4.

}

*Physical review letters*, vol. 121, no. 4, 045901. https://doi.org/10.1103/PhysRevLett.121.045901

**Phonon-Polariton Mediated Thermal Radiation and Heat Transfer among Molecules and Macroscopic Bodies : Nonlocal Electromagnetic Response at Mesoscopic Scales.** / Venkataram, Prashanth S.; Hermann, Jan; Tkatchenko, Alexandre; Rodriguez, Alejandro W.

Research output: Contribution to journal › Article › peer-review

TY - JOUR

T1 - Phonon-Polariton Mediated Thermal Radiation and Heat Transfer among Molecules and Macroscopic Bodies

T2 - Nonlocal Electromagnetic Response at Mesoscopic Scales

AU - Venkataram, Prashanth S.

AU - Hermann, Jan

AU - Tkatchenko, Alexandre

AU - Rodriguez, Alejandro W.

N1 - Funding Information: Venkataram Prashanth S. 1 Hermann Jan 2 Tkatchenko Alexandre 2 Rodriguez Alejandro W. 1 Department of Electrical Engineering, 1 Princeton University , Princeton, New Jersey 08544, USA Physics and Materials Science Research Unit, 2 University of Luxembourg , L-1511 Luxembourg 27 July 2018 27 July 2018 121 4 045901 19 March 2018 © 2018 American Physical Society 2018 American Physical Society Thermal radiative phenomena can be strongly influenced by the coupling of phonons and long-range electromagnetic fields at infrared frequencies. Typically employed macroscopic descriptions of thermal fluctuations often ignore atomistic effects that become relevant at nanometric scales, whereas purely microscopic treatments ignore long-range, geometry-dependent electromagnetic effects. We describe a mesoscopic framework for modeling thermal fluctuation phenomena among molecules near macroscopic bodies, conjoining atomistic treatments of electronic and vibrational fluctuations obtained from density functional theory in the former with continuum descriptions of electromagnetic scattering in the latter. The interplay of these effects becomes particularly important at mesoscopic scales, where phonon polaritons can be strongly influenced by the objects’ finite sizes, shapes, and nonlocal or many-body response to electromagnetic fluctuations. We show that, even in small but especially in elongated low-dimensional molecules, such effects can modify thermal emission and heat transfer by orders of magnitude and produce qualitatively different behavior compared to predictions based on local, dipolar, or pairwise approximations. National Science Foundation 10.13039/100000001 DMR-1454836 DMR 1420541 DGE 1148900 Fonds National de la Recherche Luxembourg 10.13039/501100001866 FNR-11360857 Radiative heat transfer is important when characterizing molecular structures, including graphene, fullerenes, carbon nanotubes, and other low-dimensional materials, and has applications to thermophotovoltaic power generation [1–4] , lasers [5,6] , and single-molecule junctions [7,8] . Recent theoretical [9–11] and experimental [12–15] works in this area have explored transport at nanometric scales, where the interplay of electromagnetic (EM) and phonon mediated transport can no longer be captured solely by macroscopic electrodynamics [12–14] . We present a framework for modeling thermal fluctuation phenomena among molecules near continuum bodies at mesoscopic scales and apply it to compute heat transfer between molecules in contact with separate thermal reservoirs. Our approach captures the emergence of phonon polaritons [16,17] , which typically arise at infrared frequencies and hence probe the bulk of the Planck distribution near room temperature, and their impact on charge delocalization (nonlocality) and radiation. By extending microscopic ab initio models of electronic response in molecular bodies based on density functional theory (DFT) to include molecular vibrations mediated by short-range interactions (phonons), and combining these with classical EM scattering techniques applicable to arbitrary macroscopic objects [18,19] , we account for their coupling to collective electronic excitations (plasmon polaritons), including retardation, multiple scattering, and geometry-dependent many-body effects to all orders. Our framework, detailed below, assumes the Born-Oppenheimer approximation for the nuclear motion and a harmonic approximation to the chemical bonds effecting molecular phonons and ignores thermally driven changes in molecular conformations, all of which are only expected to be violated at large temperatures. For illustration, we investigate thermal emission and heat transfer among fullerene and carbyne molecules at different temperatures, either in vacuum or above a metallic surface. We find differences of over an order of magnitude in the emitted and transferred powers relative to common dipolar approximations that neglect finite size and many-body EM effects, with the greatest differences obtained for elongated carbyne wires. We show that the proximity of the molecules to the metallic surface and resulting induction of image charges can qualitatively change their mutual heat transfer. More interestingly, we find that molecular phonon polaritons can delocalize the spatial extent of their polarization response, from a few angstroms up to tens of nanometers, depending on the shape of the molecules and their proximity to the surface. Such nonlocality qualitatively changes the heat transfer (eventually saturating) at nanometric scales. Recent theoretical models accounting for atomic-scale and vibrational effects on radiative emission and heat transfer suffer from a number of limitations [6,20–27] . Commonly, atomistic mechanical Green’s functions [10,12,16,17,28,29] are used to model phonon mediated transport across material boundaries. While the harmonic spring constants that lead to collective mechanical oscillations have recently been derived through DFT methods, they have thus far only been applied to bulk media, whereas heat transport across vacuum gaps or in structured media has been largely ignored or approximated through ad hoc pairwise models known to fail in complex geometries [30–34] . Other atomistic methods include molecular dynamics [9,35–37] , which capture the anharmonicity of covalent bonds but suffer from similar pitfalls of heuristically approximating noncovalent interactions. Macroscopic approaches based on fluctuational electrodynamics arise from discrete dipolar and multipolar methods [11,38–41] and other continuum EM treatments based on discretization of bulk susceptibilities and fields [1–3,18,19,42–49] , which discretize bulk susceptibilities using point dipoles or other basis functions. While such methods capture the impact of long-range EM effects, they typically treat material responses at atomic scales via semiempirical models rather than ab initio calculations and can therefore miss important effects arising from the spatial dispersion or nonlocality of electronic wave functions and its interplay with geometry. Even for a material like graphene with highly delocalized electronic and mechanical response [2,50–52] , nonlocality is typically treated phenomenologically through macroscopic approximations of quantum behavior. Motivated by recent predictions [31–34,53] of many-body van der Waals (vdW) interactions in low-dimensional systems and experiments investigating heat transfer at the nanoscale [12–15,54] , both of which probe regimes where continuum electrodynamics may no longer be valid, we propose a framework that can capture such thermal radiative processes at mesoscopic scales. Background .—We consider a collection of N mol molecules labeled k and described by electric susceptibilities V k in the presence of a collection of macroscopic bodies described collectively by a macroscopic susceptibility V env , relating induced polarization and electric fields; the former are maintained at temperatures T k , and the latter is maintained at a different temperature T env (henceforth assumed to be zero for conceptual convenience). The radiative energy transfer from molecule m to n (which may or may not be the same), can be shown to be [55] Φ n ( m ) = - 1 2 π Tr [ Im ( V m - 1 ⋆ ) P m T ⋆ asym ( P n G env ) T P m ] , (1) where T - 1 = ∑ k V k - 1 - G env is known as a scattering transition operator, which describes EM scattering (to all orders) by the collection of molecules in the presence of all bodies [43,57] . G env = ( G 0 - 1 - V env ) - 1 is the classical EM Green’s function [43,53,57] , and P k projects onto the space spanned by the degrees of freedom of molecule k . Here, ⋆ implies complex conjugation, † implies conjugate transposition (adjoint), Im ( A ) = ( A - A ⋆ ) / ( 2 i ) , and asym ( A ) = ( A - A † ) / ( 2 i ) . Given (1) , the far-field thermal emission from molecule m , W ( m ) = ∑ n = 1 N s n m Φ n ( m ) Θ ( ω , T n ) , (2) and the net heat transfer from molecule m to molecule n , W m → n = Φ n ( m ) [ Θ ( ω , T m ) - Θ ( ω , T n ) ] , (3) where Θ ( ω , T ) is the Planck function, while the sign function s n m = 1 – 2 δ n m accounts for the direction of energy flow for the emission of a body. In either case, the net power is the integral over all frequencies, P = ∫ 0 ∞ W ( ω ) d ω . The above formulas require accurate and fast computations of G env and V k . The former can be computed using any number of state-of-the-art classical EM techniques for arbitrary macroscopic geometries, including spectral, boundary element, or finite-difference methods [18,19,46–49] . The molecular susceptibilities V k generally require quantum descriptions, but recent work in the related field of vdW interactions [30,32,33,58] has shown that accurate models of the valence electronic response of insulating or weakly metallic molecules, including carbon allotropes, can be obtained by expressing V k = ∑ p , q α p i , q j | f p e i ⟩ ⟨ f q e j | in a localized basis | f p ⟩ (along Cartesian direction e i ) associated with each atom p by modeling its valence electrons as a quantum oscillator in its ground state. We extend that framework in the following way: each atom p in a given molecule comprises an effective valence electron that couples to long-range EM fields and has charge q p , mass m e p , and damping coefficient b e p = m e p γ e p , as well as a nucleus of mass m I p which is coupled only to its associated valence electron with spring constant k e p and to neighboring nuclei [5,16,17,25] with spring constants K p q . While m I p are obtained from elemental data and γ e p come from empirical data [59] , m e p , q p , k e p , and K p q , along with the molecular geometries, are all obtained from DFT [32,34] (see Supplemental Material [55] for more details), thereby accounting for short-range electrostatic, hybridization, and quantum exchange effects. The nuclear harmonic couplings, in particular, arise from lowest-order expansions of DFT force constants, producing anisotropy (tensorial couplings) due to the directionality of covalent bonds. This leads to the simple frequency-domain equations of motion [ K e - i ω B e - ω 2 M e - K e - K e K e + K I - ω 2 M I ] [ x e x I ] = [ Q e e e 0 ] , (4) where ( Q e , M e , M I , K e , K I , B e ) are 3 N k × 3 N k matrices collecting, respectively, q p , m e p , m I p , k e p , K p q , and b e p . The equations of motion determine the amplitude response of the nuclear positions x I and electronic dipole moments p e = Q e x e in the presence of an electric field, represented by a 3 N k -dimensional vector e e obtained by evaluating | E ⟩ at every atomic position. Solving for p e = α e e yields the electric susceptibility matrix, α = Q e [ K e - i ω B e - ω 2 M e - K e ( K e + K I - ω 2 M I ) - 1 K e ] - 1 Q e , (5) which enters the expansion of V k above. The discrepancies in the masses of electronic and nuclear oscillators ensures that the poles of α separate into ultraviolet resonances, corresponding to predominantly electronic modes, and infrared resonances, corresponding to phononic modes, the latter of which is primarily responsible for thermal emission. The form of K I , coupling neighboring nuclei, ensures that nonlocality enters the polarization response in a microscopic rather than phenomenological manner. Hence, in contrast to point dipoles or macroscopic objects modeled via local susceptibilities [13,14,21,60,61] , which lead to diverging EM fields as two or more bodies approach one another, the molecular susceptibility above will always be finite. In particular, the electronic response is accurately described by localized Gaussian basis functions [32,34,62–66] , f p ( x ) = ( 2 π σ p ) - 3 exp ( - ( x - x p ) 2 2 σ p 2 ) , (6) centered at the locations x p of each atom p , normalized such that ∫ d 3 x f p = 1 , and featuring a width that, rather than being phenomenological [67,68] , depends on the bare atomic polarizability at every frequency via σ p ( ω ) = [ ( | α p ( ω ) | / 72 π 3 ) ] 1 / 3 [32,69] , where α p ( ω ) = ∑ q , j α p j , q j ( ω ) . [Note that the effective polarizability of each atom is defined as a trace over the delocalized susceptibility matrix arising from the solution of the above equations of motion (5) ; physically, it includes both short-range Coulombic electronic interactions and polarization smearing due to short-range ion-ion interactions (described by the spring matrix K I ).] Such basis functions mitigate short-distance EM divergences through the smearing of the response over nonzero widths σ p , leading to finite matrix elements, ⟨ f p | ( ⋯ ) f q ⟩ = ∫ d 3 x d 3 x ′ f p ( x ) ( ⋯ ) f q ( x ′ ) < ∞ , even for coincident atoms p = q . While these widths are typically smaller than 1 Å [32,34] when phonons are ignored, we find that the introduction of phonons delocalizes the effective atomic polarization response in a way that depends strongly on the shapes of the molecules, leading to dramatically larger widths, on a scale ξ ∼ 1 nm . As we show below, this stronger nonlocality has a noticeable impact on heat transfer among molecules at nanometric separations. Furthermore, since α enters T , which includes long-range EM scattering, thermal fluctuations will be largest near the frequency poles of T , corresponding to the dressed, phonon-polariton resonances of the fully coupled hybrid molecular-macroscopic system. Note that the root-mean-square amplitude of quantum-thermal fluctuations in x e and x I themselves is over 2 orders of magnitude smaller than the covalent bond length scale of an angstrom, so the equilibrium atomic positions x p are essentially fixed. Predictions .—We now exploit this framework to present new proof-of-concept predictions of thermal emission and heat transfer. For simplicity, we consider G env as arising from the field response in either vacuum or above a gold plate, modeled as a perfect electrically conducting surface [70] . We consider either a C 500 fullerene of radius 1 nm or a 500-atom-long carbyne wire of length 65 nm in various configurations above the plane. Where appropriate, we explicitly compare results obtained by evaluating the retarded many-body (RMB) transfer rates (1) , which capture the impact of multiple scattering, molecular geometry, and EM screening to all orders, to those obtained from a dipolar approximation (DA) that contracts the bare molecular responses V into a point particle, with the susceptibilities and screened responses replaced by contracted dipolar polarizabilities ( α ∞ ) i j = ∑ p , q α p i , q j . For convenience, the integrated powers are normalized to the emitted power P B = σ B T 4 A of a blackbody of the same temperature T and area A . Within this approximation, the thermal emission rate of a molecule can be approximated (to lowest order in the scattering) by [Note that it is also possible to compare the RMB results against a DA that includes higher-order scattering effects [53] , e.g. multiple scattering that becomes relevant in resonant situations; we find that thermal emission matches the RMB description as the molecules are far enough away from the conducting surface to be well-approximated as screened point dipoles, but near-field heat transfer in such a DA still does not match the RMB description due to neglect of the confluence of finite size and smeared nonlocal response relevant at shorter separations.] Φ ∞ ( 0 ) = - 1 2 π Tr [ Im α ∞ Im G env ] . (7) For a molecule at some position r 0 = z e z above the plate, Im G env can be obtained from image theory; it is a constant in the near field and in the far field ( z ≫ ω / c ) given by Im G env ( ω , z , z ) = ( ω 3 / 6 π c 3 ) I - ( ω 2 / 8 π c 2 z ) sin [ ( 2 ω z / c ) ] ( I - e z ⊗ e z ) . Figure 1 shows the far-field thermal emission of either molecule at T = 300 K as a function of their distance z above the gold plate at zero kelvin, normalized by the emission rate P B (on the order of 10 - 14 W ) of blackbodies of the same effective areas. The C 500 fullerene is a compact molecule whose vibrational modes and resulting polaritons are limited in extent by its small size. As a result, its emission spectrum consists of a small set of resonances in a narrow frequency range and is dominated by a single peak at ω 0 ≈ 7 × 10 13 rad / s , causing the emission to scale as z - 1 sin ( 2 ω 0 z / c ) in the limit of z ≫ λ T = ( ℏ c / k B T ) ≈ 7.6 μ m , where λ T is the thermal wavelength. At small z , the emission approaches a constant value that is attenuated relative to vacuum by screening of the field response near the gold surface. Similar behavior is observed in the case of a carbyne wire parallel to the surface, except that the elongated shape of the molecule enables longer-range vibrational modes and stronger polarization anisotropy, which further suppresses emission at small z . For both molecules, the DA (7) underestimates the far-field emission by almost an order of magnitude for all z , due to the neglect of long-range EM scattering and screening effects. More dramatically, for the carbyne wire in the RMB treatment, as z increases, a set of closely spaced infrared polaritonic resonances begin to contribute, causing a beating pattern in the emission as a function of z . These close resonances do not arise within the DA, which only captures the impact of short-range EM interactions on the polarization response. 1 10.1103/PhysRevLett.121.045901.f1 FIG. 1. Contribution of phonon polaritons to molecular thermal emission. (a) RMB (black) and DA (magenta) emission rates of a C 500 fullerene held at T = 300 K above a zero-temperature gold plate as a function of vertical separation z , normalized by the power emitted by a blackbody of the same area P B = σ B T 4 A , where A = 4 π R 2 and R = 1 nm . (Insets) A schematic of the fullerene above the plate as well as the RMB emission spectrum in the limit z → ∞ , along with a schematic of the dominant polaritonic mode. (b) Same as (a) but for a carbyne wire parallel to the plate, where A = 2 π R L , R = 0.13 nm , and L = 65 nm . The RMB emission spectrum is shown at both z = 0.1 μ m (blue) and z = 100 μ m (red). Figure 2 considers the heat transfer between two molecules, one at T = 300 K and the other at zero temperature, as a function of their mutual horizontal separation d for several different values of their (identical) distances z above the zero-temperature plate. For two fullerenes far from the plate, we find that the flux rate transitions from 1 / d 2 in the far field to 1 / d 6 in the near field ( d ≪ λ T ), leading to flux rates P ≫ P B , consistent with the dipolar approximation [60,61] . However, as d decreases further to become comparable to the phonon-induced nonlocal length scale ξ ∼ 1 nm , the EM response begins to probe the finite sizes and smeared electronic response of the molecules, causing the heat transfer to saturate (in contrast to the diverging flux rate arising in the dipolar picture) and the resulting power-law exponent to approach zero. Hence, (3) seamlessly captures the transition between radiation and conduction. 2 10.1103/PhysRevLett.121.045901.f2 FIG. 2. Impact of nonlocal response, shape, and many-body effects on heat exchange between molecules. (a) Power law of heat transfer, ∂ ( log P ) / ∂ ( log d ) , between two C 500 fullerenes in vacuum held at T = 300 and 0 K, with respect to their surface-surface separation d . (Inset) The transfer versus d , normalized by the emission rate from a corresponding room-temperature blackbody P B . (b) Same as (a) but for two 500-atom-long parallel carbyne wires oriented parallel to a gold plate. The various curves denote different values of their vertical distance z from the plate. In the case of two parallel carbyne wires above a gold plate, the departures from dipolar and pairwise predictions are even more stark. For most separations d , as z decreases, the proximate surface severely attenuates the heat flux relative to vacuum (by over 4 orders of magnitude at z ∼ ξ ), except at very small values of d ≲ ξ ≪ z , where the overlapping nonlocal molecular response functions leads to saturation of the heat transfer independent of z . In the far field d ≫ λ T , as the molecules are much smaller than λ T , the metal plate does not affect the heat transfer much, so the power law attains the far-field dipolar limit of - 2 , and the heat transfer powers converge for different z . As d decreases further, if z ≫ λ T , the power law reaches the near-field dipolar limit of - 6 , but if z < λ T , the image dipoles induced in the proximate surface lead to an effective triple-dipole configuration with a power law of - 9 , akin to the so-called Axilrod-Teller-Muto correction for dipolar vdW interactions [30,32,71] . These power laws exhibit progressively stronger deviations from either the dipolar limit of - 6 or the pairwise limit of - 5 (obtained by summing the dipolar heat transfer between pairs of atoms across these parallel linear molecules) in the near field as d decreases further, due to the confluence of various competing length scales, including z , d , the wire lengths, ξ , and λ T . If z > ξ , the heat transfer begins to saturate for d ≲ 10 nm due to the nonlocal response, and the power laws are dominated by pairwise interactions of proximate dressed carbon atoms, so the power laws at each z match those of fullerenes in vacuum, as they are both carbon allotropes. The situation becomes more complicated when both d , z ≲ ξ [72] , due to the interplay of elongated molecular shapes, finite sizes, and nonlocal response of the wires and their images in the conducting surface. Note that, for z < ξ , the atomism of the conducting surface becomes more relevant and should also be treated, but this requires extending our current framework to handle periodically extended objects. Concluding remarks .—We have demonstrated a mesoscopic approach to computing nonequilibrium thermal emission and radiative heat transfer among molecules in the presence of planar metallic surfaces, accounting for molecular nonlocal response as well as many-body and multiple scattering effects to all orders. Our approach demonstrates significant deviations in the emitted and transferred power from commonly used local or dipolar approximations, particularly in the near field and even at relatively large separations of 10 nm, where the modification to long-range EM interactions due to the proximate surface as well as the nonlocal molecular susceptibility both matter; moreover, it has the virtue of deriving atomic characteristics from ab initio DFT methods, in contrast to phenomenological treatments of nonlocal response [73–75] , which give heat transfer results that depend strongly on the nonlocal material model and boundary conditions considered. Our approach can be generalized to consider extended materials like graphene or organic crystals, which might enable accurate computations of radiative contributions to bulk thermal conductivity beyond the Fourier regime, as well as actively studied experimental situations, including molecular materials in the vicinity of heated atomic force microscopy tips with applications to molecular junctions [7,8] . This work was supported by the National Science Foundation under Grants No. DMR-1454836, DMR 1420541, and DGE 1148900, and the Luxembourg National Research within the FNR-CORE program (No. FNR-11360857). P. S. V. thanks Chinmay Khandekar, Weiliang Jin, Sean Molesky, and Theerachai Chanyaswad for the helpful discussions. [1] 1 R. Messina and P. Ben-Abdallah , Sci. Rep. 3 , 1383 ( 2013 ). 10.1038/srep01383 [2] 2 V. B. Svetovoy and G. Palasantzas , Phys. Rev. Applied 2 , 034006 ( 2014 ). PRAHB2 2331-7019 10.1103/PhysRevApplied.2.034006 [3] 3 O. Ilic , M. Jablan , J. D. Joannopoulos , I. Celanovic , and M. Soljačić , Opt. Express 20 , A366 ( 2012 ). OPEXFF 1094-4087 10.1364/OE.20.00A366 [4] 4 A. Lenert , D. M. Bierman , Y. Nam , W. R. Chan , I. Celanović , M. Soljačić , and E. N. Wang , Nat. Nanotechnol. 9 , 126 ( 2014 ). NNAABX 1748-3387 10.1038/nnano.2013.286 [5] 5 X. L. Ruan and M. Kaviany , Phys. Rev. B 73 , 155422 ( 2006 ). PRBMDO 1098-0121 10.1103/PhysRevB.73.155422 [6] 6 S. V. Boriskina , J. K. Tong , W.-C. Hsu , B. Liao , Y. Huang , V. Chiloyan , and G. Chen , Nanophotonics 5 , 134 ( 2016 ). 10.1515/nanoph-2016-0010 [7] 7 L. Cui , R. Miao , C. Jiang , E. Meyhofer , and P. Reddy , J. Chem. Phys. 146 , 092201 ( 2017 ). JCPSA6 0021-9606 10.1063/1.4976982 [8] 8 L. Cui , W. Jeong , S. Hur , M. Matt , J. C. Klöckner , F. Pauly , P. Nielaba , J. C. Cuevas , E. Meyhofer , and P. Reddy , Science 355 , 1192 ( 2017 ). SCIEAS 0036-8075 10.1126/science.aam6622 [9] 9 L. Cui , Y. Feng , and X. Zhang , J. Phys. Chem. A 119 , 11226 ( 2015 ). JPCAFH 1089-5639 10.1021/acs.jpca.5b07995 [10] 10 J. B. Pendry , K. Sasihithlu , and R. V. Craster , Phys. Rev. B 94 , 075414 ( 2016 ). PRBMDO 2469-9950 10.1103/PhysRevB.94.075414 [11] 11 S. Edalatpour and M. Francoeur , Phys. Rev. B 94 , 045406 ( 2016 ). PRBMDO 2469-9950 10.1103/PhysRevB.94.045406 [12] 12 V. Chiloyan , J. Garg , K. Esfarjani , and G. Chen , Nat. Commun. 6 , 6755 ( 2015 ). NCAOBW 2041-1723 10.1038/ncomms7755 [13] 13 K. Kloppstech , N. Könne , S.-A. Biehs , A. W. Rodriguez , L. Worbes , D. Hellmann , and A. Kittel , Nat. Commun. 8 , 14475 ( 2017 ). NCAOBW 2041-1723 10.1038/ncomms14475 [14] 14 L. Cui , W. Jeong , V. Fernández-Hurtado , J. Feist , F. J. García-Vidal , J. C. Cuevas , E. Meyhofer , and P. Reddy , Nat. Commun. 8 , 14479 ( 2017 ). NCAOBW 2041-1723 [15] 15 R. St-Gelais , B. Guha , L. Zhu , S. Fan , and M. Lipson , Nano Lett. 14 , 6971 ( 2014 ). NALEFD 1530-6984 10.1021/nl503236k [16] 16 Z. Tian , K. Esfarjani , and G. Chen , Phys. Rev. B 89 , 235307 ( 2014 ). PRBMDO 1098-0121 10.1103/PhysRevB.89.235307 [17] 17 Z. Tian , K. Esfarjani , and G. Chen , Phys. Rev. B 86 , 235304 ( 2012 ). PRBMDO 1098-0121 10.1103/PhysRevB.86.235304 [18] 18 A. W. Rodriguez , M. T. H. Reid , and S. G. Johnson , Phys. Rev. B 88 , 054305 ( 2013 ). PRBMDO 1098-0121 10.1103/PhysRevB.88.054305 [19] 19 A. G. Polimeridis , M. T. H. Reid , W. Jin , S. G. Johnson , J. K. White , and A. W. Rodriguez , Phys. Rev. B 92 , 134202 ( 2015 ). PRBMDO 1098-0121 10.1103/PhysRevB.92.134202 [20] 20 H. Hajian , A. Ghobadi , S. A. Dereshgi , B. Butun , and E. Ozbay , J. Opt. Soc. Am. B 34 , D29 ( 2017 ). JOBPDE 0740-3224 10.1364/JOSAB.34.000D29 [21] 21 S. Shen , A. Narayanaswamy , and G. Chen , Nano Lett. 9 , 2909 ( 2009 ). NALEFD 1530-6984 10.1021/nl901208v [22] 22 F. J. García de Abajo , ACS Nano 7 , 11409 ( 2013 ). ANCAC3 1936-0851 10.1021/nn405367e [23] 23 D. N. Basov , M. M. Fogler , and F. J. García de Abajo , Science 354 , aag 1992 ( 2016 ). SCIEAS 0036-8075 10.1126/science.aag1992 [24] 24 T. Luo and G. Chen , Phys. Chem. Chem. Phys. 15 , 3389 ( 2013 ). PPCPFQ 1463-9076 10.1039/c2cp43771f [25] 25 V. P. Carey , G. Chen , C. Grigoropoulos , M. Kaviany , and A. Majumdar , Nanoscale Micro. Thermophys. Eng. 12 , 1 ( 2008 ). 1556-7265 10.1080/15567260801917520 [26] 26 S. V. Boriskina , H. Ghasemi , and G. Chen , Mater. Today 16 , 375 ( 2013 ). MATOBY 1369-7021 10.1016/j.mattod.2013.09.003 [27] 27 D. G. Cahill , P. V. Braun , G. Chen , D. R. Clarke , S. Fan , K. E. Goodson , P. Keblinski , W. P. King , G. D. Mahan , A. Majumdar , H. J. Maris , S. R. Phillpot , E. Pop , and L. Shi , Appl. Phys. Rev. 1 , 011305 ( 2014 ). APRPG5 1931-9401 10.1063/1.4832615 [28] 28 A. Dhar and D. Roy , J. Stat. Phys. 125 , 801 ( 2006 ). JSTPBS 0022-4715 10.1007/s10955-006-9235-3 [29] 29 N. Mingo and L. Yang , Phys. Rev. B 68 , 245406 ( 2003 ). PRBMDO 0163-1829 10.1103/PhysRevB.68.245406 [30] 30 A. Tkatchenko , A. Ambrosetti , and R. A. DiStasio, Jr. , J. Chem. Phys. 138 , 074106 ( 2013 ). JCPSA6 0021-9606 10.1063/1.4789814 [31] 31 V. V. Gobre and A. Tkatchenko , Nat. Commun. 4 , 2341 ( 2013 ). NCAOBW 2041-1723 10.1038/ncomms3341 [32] 32 R. A. DiStasio, Jr. , V. V. Gobre , and A. Tkatchenko , J. Phys. Condens. Matter 26 , 213202 ( 2014 ). JCOMEL 0953-8984 10.1088/0953-8984/26/21/213202 [33] 33 A. Tkatchenko , Adv. Funct. Mater. 25 , 2054 ( 2015 ). AFMDC6 1616-301X 10.1002/adfm.201403029 [34] 34 A. Ambrosetti , N. Ferri , R. A. DiStasio, Jr. , and A. Tkatchenko , Science 351 , 1171 ( 2016 ). SCIEAS 0036-8075 10.1126/science.aae0509 [35] 35 A. Henry and G. Chen , Phys. Rev. Lett. 101 , 235502 ( 2008 ). PRLTAO 0031-9007 10.1103/PhysRevLett.101.235502 [36] 36 K. Esfarjani , G. Chen , and H. T. Stokes , Phys. Rev. B 84 , 085204 ( 2011 ). PRBMDO 1098-0121 10.1103/PhysRevB.84.085204 [37] 37 E. G. Noya , D. Srivastava , L. A. Chernozatonskii , and M. Menon , Phys. Rev. B 70 , 115416 ( 2004 ). PRBMDO 1098-0121 10.1103/PhysRevB.70.115416 [38] 38 R. Messina , M. Tschikin , S.-A. Biehs , and P. Ben-Abdallah , Phys. Rev. B 88 , 104307 ( 2013 ). PRBMDO 1098-0121 10.1103/PhysRevB.88.104307 [39] 39 A. Pérez-Madrid , J. M. Rubí , and L. C. Lapas , Phys. Rev. B 77 , 155417 ( 2008 ). PRBMDO 1098-0121 10.1103/PhysRevB.77.155417 [40] 40 S. Edalatpour and M. Francoeur , J. Quant. Spectrosc. Radiat. Transfer 133 , 364 ( 2014 ). JQSRAE 0022-4073 10.1016/j.jqsrt.2013.08.021 [41] 41 S. Edalatpour , J. DeSutter , and M. Francoeur , J. Quant. Spectrosc. Radiat. Transfer 178 , 14 ( 2016 ). JQSRAE 0022-4073 10.1016/j.jqsrt.2015.12.027 [42] 42 R. Yu , A. Manjavacas , and F. J. G. de Abajo , Nat. Commun. 8 , 2 ( 2017 ). NCAOBW 2041-1723 10.1038/s41467-016-0013-x [43] 43 M. Krüger , G. Bimonte , T. Emig , and M. Kardar , Phys. Rev. B 86 , 115423 ( 2012 ). PRBMDO 1098-0121 10.1103/PhysRevB.86.115423 [44] 44 S. Molesky and Z. Jacob , Phys. Rev. B 91 , 205435 ( 2015 ). PRBMDO 1098-0121 10.1103/PhysRevB.91.205435 [45] 45 B. Zhao , B. Guizal , Z. M. Zhang , S. Fan , and M. Antezza , Phys. Rev. B 95 , 245437 ( 2017 ). PRBMDO 2469-9950 10.1103/PhysRevB.95.245437 [46] 46 C. R. Otey , L. Zhu , S. Sandhu , and S. Fan , J. Quant. Spectrosc. Radiat. Transfer 132 , 3 ( 2014 ). JQSRAE 0022-4073 10.1016/j.jqsrt.2013.04.017 [47] 47 M. T. H. Reid , A. W. Rodriguez , and S. G. Johnson , Proc. IEEE 101 , 531 ( 2013 ). IEEPAD 0018-9219 10.1109/JPROC.2012.2191749 [48] 48 A. W. Rodriguez , O. Ilic , P. Bermel , I. Celanovic , J. D. Joannopoulos , M. Soljačić , and S. G. Johnson , Phys. Rev. Lett. 107 , 114302 ( 2011 ). PRLTAO 0031-9007 10.1103/PhysRevLett.107.114302 [49] 49 C. Luo , A. Narayanaswamy , G. Chen , and J. D. Joannopoulos , Phys. Rev. Lett. 93 , 213905 ( 2004 ). PRLTAO 0031-9007 10.1103/PhysRevLett.93.213905 [50] 50 M. Jablan , M. Soljačić , and H. Buljan , Phys. Rev. B 83 , 161409 ( 2011 ). PRBMDO 1098-0121 10.1103/PhysRevB.83.161409 [51] 51 E. H. Hwang , R. Sensarma , and S. Das Sarma , Phys. Rev. B 82 , 195406 ( 2010 ). PRBMDO 1098-0121 10.1103/PhysRevB.82.195406 [52] 52 S. Xiao , X. Zhu , B.-H. Li , and N. A. Mortensen , Frontiers Phys. 11 , 117801 ( 2016 ). 10.1007/s11467-016-0551-z [53] 53 P. S. Venkataram , J. Hermann , A. Tkatchenko , and A. W. Rodriguez , Phys. Rev. Lett. 118 , 266802 ( 2017 ). PRLTAO 0031-9007 10.1103/PhysRevLett.118.266802 [54] 54 B. Song , Y. Ganjeh , S. Sadat , D. Thompson , A. Fiorino , V. Fernández-Hurtado , J. Feist , F. J. Garcia-Vidal , J. C. Cuevas , P. Reddy , Nat. Nanotechnol. 10 , 253 ( 2015 ). NNAABX 1748-3387 10.1038/nnano.2015.6 [55] 55 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.121.045901 , which also cites Ref. [56], for a derivation of (1) along with further details on the molecular susceptibility parameters. [56] 56 L. Novotny and B. Hecht , Principles of Nano-Optics ( Cambridge University Press , Cambridge, England, 2006 ), p. 335 . [57] 57 S. J. Rahi , T. Emig , N. Graham , R. L. Jaffe , and M. Kardar , Phys. Rev. D 80 , 085021 ( 2009 ). PRVDAQ 1550-7998 10.1103/PhysRevD.80.085021 [58] 58 A. Ambrosetti , A. M. Reilly , R. A. DiStasio , and A. Tkatchenko , J. Chem. Phys. 140 , 18A508 ( 2014 ). JCPSA6 0021-9606 10.1063/1.4865104 [59] 59 S. Y. Buhmann , S. Scheel , S. A. Ellingsen , K. Hornberger , and A. Jacob , Phys. Rev. A 85 , 042513 ( 2012 ). PLRAAN 1050-2947 10.1103/PhysRevA.85.042513 [60] 60 A. I. Volokitin and B. N. J. Persson , Phys. Rev. B 63 , 205404 ( 2001 ). PRBMDO 0163-1829 10.1103/PhysRevB.63.205404 [61] 61 G. Domingues , S. Volz , K. Joulain , and J.-J. Greffet , Phys. Rev. Lett. 94 , 085901 ( 2005 ). PRLTAO 0031-9007 10.1103/PhysRevLett.94.085901 [62] 62 A. G. Donchev , J. Chem. Phys. 125 , 074713 ( 2006 ). JCPSA6 0021-9606 10.1063/1.2337283 [63] 63 A. D. Phan , L. M. Woods , and T.-L. Phan , J. Appl. Phys. 114 , 044308 ( 2013 ). JAPIAU 0021-8979 10.1063/1.4816446 [64] 64 Y. V. Shtogun and L. M. Woods , J. Phys. Chem. Lett. 1 , 1356 ( 2010 ). JPCLCD 1948-7185 10.1021/jz100309m [65] 65 H.-Y. Kim , J. O. Sofo , D. Velegol , M. W. Cole , and A. A. Lucas , Langmuir 23 , 1735 ( 2007 ). LANGD5 0743-7463 10.1021/la061802w [66] 66 M. W. Cole , D. Velegol , H.-Y. Kim , and A. A. Lucas , Mol. Simul. 35 , 849 ( 2009 ). MOSIEA 0892-7022 10.1080/08927020902929794 [67] 67 J. Mahanty and B. W. Ninham , J. Chem. Soc., Faraday Trans. 2 71 , 119 ( 1975 ). JCFTBS 0300-9238 10.1039/f29757100119 [68] 68 M. J. Renne , Physica (Utrecht) 53 , 193 ( 1971 ). PHYSAG 0031-8914 10.1016/0031-8914(71)90070-X [69] 69 A. Mayer , Phys. Rev. B 75 , 045407 ( 2007 ). PRBMDO 1098-0121 10.1103/PhysRevB.75.045407 [70] For computational convenience, and since gold is essentially perfectly reflecting at infrared frequencies where nonequilibrium energy transfer is most significant, we model the gold plate as a perfect electrically conducting surface. [71] 71 B. M. Axilrod and E. Teller , J. Chem. Phys. 11 , 299 ( 1943 ). JCPSA6 0021-9606 10.1063/1.1723844 [72] Note that, in principle, as d attains angstrom length scales, the assumption of the molecules being separate breaks down due to the possible formation of new covalent bonds; while this would require recomputing K I for the coupled molecular system, our formalism for computing radiative heat transfer still holds as long as the two molecular components can be maintained at different temperatures. [73] 73 P.-O. Chapuis , S. Volz , C. Henkel , K. Joulain , and J.-J. Greffet , Phys. Rev. B 77 , 035431 ( 2008 ). PRBMDO 1098-0121 10.1103/PhysRevB.77.035431 [74] 74 F. Singer , Y. Ezzahri , and K. Joulain , J. Quant. Spectrosc. Radiat. Transfer 154 , 55 ( 2015 ). JQSRAE 0022-4073 10.1016/j.jqsrt.2014.11.016 [75] 75 R. Schmidt and S. Scheel , J. Phys. B 51 , 044003 ( 2018 ). JPAPEH 0953-4075 10.1088/1361-6455/aa9f58

PY - 2018/7/27

Y1 - 2018/7/27

N2 - Thermal radiative phenomena can be strongly influenced by the coupling of phonons and long-range electromagnetic fields at infrared frequencies. Typically employed macroscopic descriptions of thermal fluctuations often ignore atomistic effects that become relevant at nanometric scales, whereas purely microscopic treatments ignore long-range, geometry-dependent electromagnetic effects. We describe a mesoscopic framework for modeling thermal fluctuation phenomena among molecules near macroscopic bodies, conjoining atomistic treatments of electronic and vibrational fluctuations obtained from density functional theory in the former with continuum descriptions of electromagnetic scattering in the latter. The interplay of these effects becomes particularly important at mesoscopic scales, where phonon polaritons can be strongly influenced by the objects' finite sizes, shapes, and nonlocal or many-body response to electromagnetic fluctuations. We show that, even in small but especially in elongated low-dimensional molecules, such effects can modify thermal emission and heat transfer by orders of magnitude and produce qualitatively different behavior compared to predictions based on local, dipolar, or pairwise approximations.

AB - Thermal radiative phenomena can be strongly influenced by the coupling of phonons and long-range electromagnetic fields at infrared frequencies. Typically employed macroscopic descriptions of thermal fluctuations often ignore atomistic effects that become relevant at nanometric scales, whereas purely microscopic treatments ignore long-range, geometry-dependent electromagnetic effects. We describe a mesoscopic framework for modeling thermal fluctuation phenomena among molecules near macroscopic bodies, conjoining atomistic treatments of electronic and vibrational fluctuations obtained from density functional theory in the former with continuum descriptions of electromagnetic scattering in the latter. The interplay of these effects becomes particularly important at mesoscopic scales, where phonon polaritons can be strongly influenced by the objects' finite sizes, shapes, and nonlocal or many-body response to electromagnetic fluctuations. We show that, even in small but especially in elongated low-dimensional molecules, such effects can modify thermal emission and heat transfer by orders of magnitude and produce qualitatively different behavior compared to predictions based on local, dipolar, or pairwise approximations.

UR - http://www.scopus.com/inward/record.url?scp=85050659241&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85050659241&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.121.045901

DO - 10.1103/PhysRevLett.121.045901

M3 - Article

C2 - 30095944

AN - SCOPUS:85050659241

VL - 121

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 4

M1 - 045901

ER -