The dynamic heat transfer between two half-spaces separated by a vacuum gap due to coupling of their surface modes is modelled using the theory that describes the dynamic energy transfer between two coupled harmonic oscillators each separately connected to a heat bath and with the heat baths maintained at different temperatures. The theory is applied for the case when the two surfaces are made up of a polar crystal which supports surface polaritons that can be excited at room temperature and the predicted heat transfer is compared with the steady state heat transfer value calculated from standard fluctuational electrodynamics theory. It is observed that for small time intervals the value of heat flux can be significantly higher than that of steady state value.The theory of photon mediated heat transfer between macroscopic objects in close proximity to each other and separated by a vacuum gap has traditionally been treated using the macroscopic Rytov's fluctuational electrodynamics theory which assumes local thermodynamic equilibrium in the bodies in question [1][2][3][4]. This heat transfer comprises of contributions from both long-range radiative modes as well as near-field evanescent and surface modes [4,5]. When thermally excited, contributions from surface modes -which are electromagnetic eigenmodes of the surface and are characterized by their field decaying exponentially on either side of the interface -dominate the heat transfer between the surfaces when the gap is less than the thermal wavelength of radiation. This is primarily due to a peak in the density of electromagnetic states at such frequencies where these surface modes are resonantly activated as evidenced from the dispersion relation for these modes [4]. In particular, for this effect to be prominent at room temperature the surfaces should be made up of a polar crystal such as silicon carbide or silica which supports surface-phonon polariton modes in the infra-red wavelength around 10 µm and can thus be thermally excited at these temperatures.In general, resonant excitation of surface modes plays an important role in several phenomena and applications including: decreased lifetime of molecules close to metal surfaces [6], surface enhanced raman spectroscopy [7], thermal near-field spectroscopy [8,9], concept of perfect lens [10] and thermal rectification [11]. The study of coupling of surface modes across surfaces is significant as it not only plays an important role in heat transfer but also in the van der Waals and Casimir force between them [12,13]. A coupled harmonic oscillator description for the heat transfer between nanoparticles due to coupling of surface modes was arrived at by Biehs and Agarwal [14] and estimates for both dynamic and steady state heat transfer values were arrived at. Barton [15] has considered the heat flow between two harmonic oscillators using Langevin dynamics and has extended this model to planar surfaces but has limited his description to the steady-state heat flow. Yu et. al.,[16] have recently analysed the dynam...