# Phonon heat transport in cavity-mediated optomechanical nanoresonators

## Abstract

The understanding of heat transport in nonequilibrium thermodynamics is an important research frontier, which is crucial for implementing novel thermodynamic devices, such as heat engines and refrigerators. The convection, conduction, and radiation are the well-known basic ways to transfer thermal energy. Here, we demonstrate a different mechanism of phonon heat transport between two spatially separated nanomechanical resonators coupled by the cavity-enhanced long-range interactions. The single trajectory for thermalization and non-equilibrium dynamics is monitored in real-time. In the strong coupling regime, the instant heat flux spontaneously oscillates back and forth in the nonequilibrium steady states. The universal bound on the precision of nonequilibrium steady-state heat flux, i.e. the thermodynamic uncertainty relation, is verified in such a temperature gradient driven far-off equilibrium system. Our results give more insight into the heat transfer with nanomechanical oscillators, and provide a playground for testing fundamental theories in non-equilibrium thermodynamics.

## Introduction

The energy transfer and heat flux play the most important role in thermodynamics, which have been well understood on macroscopic scales in the framework of classical thermodynamics. However, as the size of the system is reduced, thermal or quantum fluctuations becoming increasingly relevant, some striking phenomena and properties have been revealed recently1,2,3,4. To study heat transport in systems of nano/microscales and at single-atom levels has remained elusive, although it has received great interest for the fundamental physics in many fields5,6,7,8 and the applications in various advanced microscopic devices9,10,11.

In parallel, remarkable accomplishments have been achieved in cavity optomechanics, such as ground-state cooling of mechanical resonators12,13,14, ultrasensitive motion detections15,16,17, nonreciprocal control of photons and phonons18,19,20, and so on. Optomechanical system owing to its high controllability has been proposed as a novel mechanism for nonlocal heat transfer and thermal management21,22,23,24,25,26, however, it has not been experimentally demonstrated because of the technical challenges in precisely manipulating optomechanical arrays inside an optical cavity.

Here we study the transport of thermal energy between two spatially separated nanomechanical resonators mediated by the cavity light field. Compared to other mechanisms of heat transport, the optomechanical coupling can be coherently manipulated and it could be, in principle, infinitely long-range, which provides a new mechanism for heat transfer with high tunability. Importantly, in the strong coupling regime, we find that the heat flux flows back and forth coherently in nonequilibrium steady-state at short time scales, and the oscillation period is fully determined by the eigenmodes of such a composite system. This means that the coherence is dominated over the dissipations in such an inherently random process. In addition, we test the constraints on the precision of steady-state heat flux, i.e., the thermodynamic uncertainty relation27,28,29,30,31,32,33,34,35,36,37, in such a thermal gradient driven far-off equilibrium system at various coupling strengths, including the strong coupling regime. Therefore, our observations provide new insights for nonequilibrium thermodynamics and nontrivial extension beyond the weak coupling, in contrast to most previous systems of stochastic thermodynamics that are overdamped3,4.

## Results

### Experimental setup

The multimode optomechanical system offers an exceptional degree of freedom to investigate phonon heat transport. The thermal energy is transported between two nanomechanical resonators mediated by a cavity field, as shown in Fig. 1a. The resonators are two SiN membranes with dimensions 1 mm × 1 mm × 50 nm. The membranes separated by 6 cm are placed inside a Fabry-Perot cavity38,39. The motions of membranes are monitored by two weak probe laser fields separately (not shown in Fig. 1a). We focus on the vibrational (1,1) modes, which are nearly degenerate with eigenfrequencies ω1,2 ~ 2π × 400 kHz. Piezos are used to precisely control the eigenfrequencies of each membrane40. The cavity is driven by a red-detuned laser field, which interacts with both membranes simultaneously due to the dynamical backaction. Therefore, two individual modes of membranes are effectively coupled mediated by the cavity field. The phonon heat transport is realized by contacting two membranes with independent thermal baths. Such a system can be equivalently viewed as two mutually coupled harmonic oscillators, as shown in Fig. 1b, with one in contact with a room temperature reservoir and the other with a high temperature reservoir (realized by exciting the membrane with a Gaussian white noise). The temperature gradient drives the system out of equilibrium with a mean heat flux from the hot reservoir to the cold reservoir under the effective optomechanical coupling. When the coupling between the membranes is zero by turning off the cavity field, the membranes have the equivalent temperatures as their own reservoirs.

### Effective Hamiltonian

The interaction Hamiltonian of such a composite system of two membranes interacting with a common cavity field is $$H_{int} = - \hbar \mathop {\sum}\nolimits_{i = 1,2} {g_i} \hat a^\dagger \hat a(\hat b_i^\dagger + \hat b_i)$$, where $${\hat{\mathrm{a}}}$$ and $$\hat b_i$$ are the annihilation operators for cavity and mechanical modes, respectively. gi is the optomechanical coupling rate of each membrane. After eliminating the cavity mode, the system can be effectively described by the following Langevin equations (see Supplementary Note 1)

$$i\frac{\partial }{{\partial {\mathrm{t}}}}\left( {\begin{array}{*{20}{c}} {\hat b_1} \\ {\hat b_2} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {\omega _1 - \frac{{i\gamma _1}}{2} + {\mathrm{{\Lambda}}}} \\ {\mathrm{{\Lambda}}} \end{array}\begin{array}{*{20}{c}} {\mathrm{{\Lambda}}} \\ {\omega _2 - \frac{{i\gamma _2}}{2} + {\mathrm{{\Lambda}}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {\hat b_1} \\ {\hat b_2} \end{array}} \right) + \left( {\begin{array}{*{20}{c}} {\sqrt {\gamma _1} \hat \eta _1} \\ {\sqrt {\gamma _2} \hat \eta _2} \end{array}} \right)$$
(1)

Here $$\hat \eta _i$$ (i = 1, 2) is the thermal Langevin noise operator with$$\langle {\hat \eta _{1,2}^\dagger \left( t \right)\hat \eta _{1,2}\left( {t^{\prime} } \right)} \rangle = \delta \left( {t - t^{\prime} } \right)k_BT_{H,L}/\hbar \omega _{1,2}$$. The mechanical resonant frequencies are assumed to be equivalent, i.e., ω1 = ω2 = ω0, which leads to the maximum heat transfer efficiency, unless it is specifically described. Ʌ = g2 χm is the effective coupling rate between membranes, where the membranes are assumed to have the same optomechanical coupling rate g, and χm is the effective mechanical susceptibility introduced by the intracavity field41, which is proportional to the intracavity photon number. Generally, the effective coupling between membranes can be both dispersive and dissipative. Here, the laser frequency is detuned far-off the cavity resonance, and consequently the interacting of membranes is dominated by a conservative coupling.

### Thermomechanical noise spectrum

To demonstrate the thermal energy transportation, we first investigate the thermomechanical noise spectral properties. Figure 2a, b illuminate the noise spectra of membranes in the weak and strong coupling regimes, respectively. The critical point that separates the strong and weak coupling regimes is $${\mathrm{{\Lambda}}}^2 - \frac{{( {\gamma _1 - \gamma _2} )^2}}{{16}} = 0$$. Notably, the mechanical spectrum shows two peaks in the strong coupling regime, and the resonant frequencies are given by the eigenvalues of the effective Hamiltonian, i.e., ω+ = ω0 and ω = ω0 + 2Ʌ, which is called the normal mode splitting. This splitting is caused by the strong phonon-phonon interaction, and it has a different mechanism with the previous observations where one cavity mode and one mechanical mode interact42,43. The frequency shifts of normal modes as a function of the intracavity photon number are plotted in Fig. 2c. Due to the effective spring constant is negative and the eigenfrequencies of membranes are also modified by the optomechanical coupling, the frequency of the breathing mode is fixed and the frequency of the center-of-mass mode reduces as the cavity photon number increases, as shown in Fig. 2c. By integrating the noise power spectra, one can obtain the effective temperatures of membranes. The root mean square vibrational amplitude of membrane $$\sqrt {k_BT/m\omega _0^2}$$ at room temperature is utilized to calibrate the effective temperature (see Supplementary Note 4). Please note that the effective temperature is defined by the motion of mechanical oscillator, instead of the motions of atoms which constitute the membrane. When the optomechanical coupling is turned on, the heat energy is transferred from the high temperature membrane to the low temperature membrane. The effective temperatures in steady states are shown in Fig. 2d. When the coupling is weak, the temperature difference between two membranes is large. The temperatures are close to their respective reservoirs’ temperatures in equilibrium. As the cavity photon number increases, the temperatures of membranes tend to be equal, which implicates that the thermal energy has been redistributed in the two membranes under the optomechanical coupling.

### Instant heat flux

To fully characterize the process of heat transfer, we investigate the heat flux in the nonequilibrium steady-state. The instant heat flux from membrane 1 to 2 is $$j_\tau = - \mathop {{\lim }}\nolimits_{t_0 \to \infty } \frac{{2m}}{\tau }(\omega _0{\mathrm{{\Lambda}}} + {\mathrm{{\Lambda}}}^2){\int}_{t_0}^{t_0 + \tau } {u_1} \dot u_2dt$$, where u1,2 is the membrane motion, and τ is the integration time44,45. Since the integration time of the measurements (typical on the order of 100 μs) is much smaller than the damping rates of membranes, then the heat flux directly obtained from the data of lock-in amplifier can be viewed as the instant heat flux. The time evolution of instant heat flux in the strong coupling regime is plotted in Fig. 3a. Remarkably, the instant heat flux oscillates back and forth between the membranes in the nonequilibrium steady-state with strong coupling (see the inset of Fig. 3a for clarity), which is not previously observed in the thermalization and in contrast with the observations in the macroscopic world. By performing Fourier transformation on the heat flux trace, the oscillation frequency as a function of cavity photon number is presented in Fig. 3b, which is coincident with the twice of the effective coupling strength between membranes. Such oscillations disappear as the coupling strength is smaller than the critical point. The inset of Fig. 3b illuminates such a transition near the critical point.

The probability density functions for various integration time are plotted in Fig. 3d. Although the instant heat flux traces exhibit completely different behaviors in the strong coupling regime, the probability density functions have exponential tails and positive mean values, which is consistent with the observations in overdamped systems4. Based on the probability density function, a symmetry function can be defined as the logarithm of the ratio of the probability finding a positive heat flux with respect to the corresponding negative one, i.e., ln[P(j)/P(−j)]40. The symmetry function at different coupling strengths is shown in Fig. 3c, which has a linear dependence on the instant current value j. The green dots and brown squares are the experimental data corresponding to the cases when the cavity photon numbers are 5 × 108 and 109, respectively. A larger optomechanical coupling strength leads to a smaller slope of the symmetry function, and indicates a more symmetric probability density function, which depends on the difference of the effective temperatures of membranes.

### Thermodynamics uncertainty relation

According to the measurements of instant heat flux for the strong coupling (Fig. 3) and weak coupling cases (see Supplementary Note 5), one can find that the dynamical fluctuation is an essential quantity, which indicates that such an optomechanical system driven by a thermal bias provides an ideal platform for studying the theories in nonequilibrium thermodynamics, the efficiency of stochastic heat engine, and the connection between information and thermodynamics3,4. The instant heat flux in the strong coupling regime could violate the second law of thermodynamics in the small time scales. Then we intend to investigate whether the fundamental constraints of stochastic thermodynamics hold in this distinct regime. The thermodynamics uncertainty relation (TUR) is a recently developed thermodynamic inequality, which offers a precise bound on the heat flux fluctuations in terms of the entropy production. Although the TUR has been well accepted and the most interesting applications take a step beyond its validation, there is rare proof in the experiment30,37.

The statement of the TUR in the nonequilibrium steady-state limit under a constant thermal bias reads27,28,29,30,31,32,33,34,35,36,37

$$\frac{{{\mathrm{Var}}(J_\tau )}}{{J_\tau ^2}} \ge \frac{{2k_B}}{{\sum _\tau }}$$
(2)

where the average steady-state entropy production is $$\sum _\tau = \langle J_\tau \rangle ( {\frac{1}{{T_L}} - \frac{1}{{T_H}}} )$$ and the integrated heat current is $$J_\tau = - \mathop {{\lim }}\nolimits_{t_0 \to \infty } 2m(\omega _0{\mathrm{{\Lambda}}} + {\mathrm{{\Lambda}}}^2){\int}_{t_0}^{t_0 + \tau } {u_1} \dot u_2dt$$. The trajectories of integrated heat current as a function of time at different coupling strengths are plotted in Fig. 4a. In the case of TL << TH, Eq. (2) simplifies to Var(Jτ)/(<Jτ > kBTL) ≥ 2. Figure 4b exhibits the TUR quantity Var(Jτ)/(<Jτ > kBTL) as a function of τ at various coupling strengths. As one can see in Fig. 4b, the TUR quantity is above the TUR limit (a value of 2) at an arbitrary time, which approaches to a constant as τ is relatively large. However, the TUR decreases as τ goes to zero, owing to the possibility that the inertia suppress the fluctuations when the time is smaller than the relaxation time of membranes. A characteristic behavior for the strong coupling case (see the inset of Fig. 4b for clarity) is that the TUR quantity shows oscillation and the period is coincident with the oscillation period of instant heat flux.

### Mean heat flux

Lastly, we demonstrate that the mean heat flux can be flexibly manipulated by tuning the laser power and other parameters, in comparison with many other systems where the thermal conductivity is fixed and the heat flux is difficult to tune. The value of mean heat flux is obtained by averaging the instant heat flux trace for a long time (dozens of seconds in contrast to hundreds of μs for the instant heat flux), which can also be extracted from the effective temperatures of membranes in the nonequilibrium steady-state (see Supplementary Note 2). The mean heat flux is zero when the laser field is off. As the intracavity photon number increases, the mean heat flux grows up and tends to be saturated, as shown in Fig. 5a. The maximum mean heat flux is limited by $$\gamma _1\gamma _2k_B(T_H - T_L)/(\gamma _1 + \gamma _2)$$. The mean heat flux is linearly related to the temperature difference of two reservoirs, i.e. TH − TL, which is depicted in Fig. 5b. Previously, we have studied the situations of resonant heat transfer. Figure 5c presents the dependence on the frequency detuning of membranes (∆ = (ω1 − ω2)/2π). The resonance behavior can be clearly seen in Fig. 5c. Interestingly, the linewidth of the resonance peak is broadened compared to the natural linewidth of mechanical resonators, which can be analogous to the power broadening in atomic physics and can be utilized to relax the condition of mechanical frequency match in the strong coupling regime.

## Discussion

The systems of optomechanical arrays offer many unique possibilities to study nonequilibrium dynamics, as they allow for a full real-time control of almost all relevant parameters. We have demonstrated the phonon heat transportation between two spatially separated nanomechanical membranes with cavity-enhanced long-rang and flexible optomechanical coupling. The spontaneous generation of coherent energy oscillation during the process of heat transfer with strong coupling is experimentally observed. The universal bound on the precision of nonequilibrium steady-state heat flux is verified. In contrast to most previous stochastic thermodynamic systems, which are overdamped, our observations might open an unexplored region. Although one eigenmode of each membrane is utilized in this work, the multimode nature of the mechanical resonator implies the possibility of using this system as a multi-channel phonon-photon interface46. Such a system provides a fertile ground for studying stochastic and quantum thermodynamics, and is a candidate for the photon controlled phononic devices and hybrid quantum networks.

## Data availability

The data that supports the results within this paper are available from the corresponding authors upon reasonable request.

## Code availability

The custom code used to obtain the results within this paper are available from the corresponding authors upon reasonable request.

## References

1. 1.

Fong, K. Y. et al. Phonon heat transfer across a vacuum through quantum fluctuations. Nature 576, 243–247 (2019).

2. 2.

Cui, L. et al. Thermal conductance of single-molecule junctions. Nature 572, 628–633 (2019).

3. 3.

Seifert, U. Stochastic thermodynamics, fluctuation theorems and molecular machines. Rep. Prog. Phys. 75, 126001 (2012).

4. 4.

Ciliberto, S. Experiments in stochastic thermodynamics: Short history and perspectives. Phys. Rev. X 7, 021051 (2017).

5. 5.

Gemmer, J., Michel, M. & Mahler, G. Quantum Thermodynamics: Emergence of Thermodynamic Behavior within Composite quantum Systems, 2nd edn, (Springer, 2009).

6. 6.

Pekola, J. P. Towards quantum thermodynamics in electronic circuits. Nat. Phys. 11, 118–123 (2015).

7. 7.

Dhar, A. Heat transport in low-dimensional systems. Adv. Phys. 57, 457 (2008).

8. 8.

Luckyanova, M. N. et al. Coherent phonon heat conduction in superlattices. Science 338, 936–939 (2012).

9. 9.

Peterson, J. P. S. et al. Experimental characterization of a spin quantum heat engine. Phys. Rev. Lett. 123, 240601 (2019).

10. 10.

Klatzow, J. et al. Experimental demonstration of quantum effects in the operation of microscopic heat engines. Phys. Rev. Lett. 122, 110601 (2019).

11. 11.

Li, N. et al. Manipulating heat flow with electronic analogs and beyond. Rev. Mod. Phys. 84, 1045 (2012).

12. 12.

Teufel, J. et al. Sideband cooling of micromechanical motion to the quantum ground state. Nature 475, 359–363 (2011).

13. 13.

Chan, J. et al. Laser cooling of a nanomechanical oscillator into its quantum ground state. Nature 478, 89–92 (2011).

14. 14.

Peterson, R. et al. Laser cooling of a micromechanical membrane to the quantum backaction limit. Phys. Rev. Lett. 116, 063601 (2016).

15. 15.

Teufel, J., Donner, T., Castellanos-Beltran, M., Harlow, J. & Lehnert, K. Nanomechanical motion measured with an imprecision below that at the standard quantum limit. Nat. Nanotechnol. 4, 820–823 (2009).

16. 16.

Schliesser, A., Arcizet, O., Riviere, R., Anetsberger, G. & Kippenberg, T. Resolved-sideband cooling and position measurement of a micromechanical oscillator close to the Heisenberg uncertainty limit. Nat. Phys. 5, 509–514 (2009).

17. 17.

Rocheleau, T. et al. Preparation and detection of a mechanical resonator near the ground state of motion. Nature 463, 72–75 (2010).

18. 18.

Xu, H. et al. Nonreciprocal control and cooling of phonon modes in an optomechanical system. Nature 568, 65–69 (2019).

19. 19.

Bernier, N. R. et al. Nonreciprocal reconfigurable microwave optomechanical circuit. Nat. Commun. 8, 604 (2017).

20. 20.

Ruesink, F., Miri, M.-A., Alu, A. & Verhagen, E. Nonreciprocity and magnetic-free isolation based on optomechanical interactions. Nat. Commun. 7, 13662 (2016).

21. 21.

Farman, F. & Bahrampour, A. R. Heat transfer between micro- and nano-mechanical systems through optical channels. J. Opt. Soc. Am. B 31, 1525–1532 (2014).

22. 22.

Xuereb, A., Genes, C., Pupillo, G., Paternostro, M. & Dantan, A. Reconfigurable long-range phonon dynamics in optomechanical arrays. Phys. Rev. Lett. 112, 133604 (2014).

23. 23.

Xuereb, A., Imparato, A. & Dantan, A. Heat transport in harmonic oscillator systems with thermal baths: application to optomechanical arrays. N. J. Phys. 17, 055013 (2015).

24. 24.

Barzanjeh, S., Aquilina, M. & Xuereb, A. Manipulating the flow of thermal noise in quantum devices. Phys. Rev. Lett. 120, 060601 (2018).

25. 25.

Seif, A., DeGottardi, W., Esfarjani, K. & Hafezi, M. Thermal management and non-reciprocal control of phonon flow via optomechanics. Nat. Commun. 9, 1207 (2018).

26. 26.

Mahmoud Ashrafi, S., Malekfar, R., Bahrampour, A. R. & Feist, J. Optomechanical heat transfer between molecules in a nanoplasmonic cavity. Phys. Rev. A 100, 013826 (2019).

27. 27.

Barato, A. C. & Seifert, U. Thermodynamic uncertainty relation for biomolecular processes. Phys. Rev. Lett. 114, 158101 (2015).

28. 28.

Gingrich, T. R., Horowitz, J. M., Perunov, N. & England, J. L. Dissipation bounds all steady state current fluctuations. Phys. Rev. Lett. 116, 120601 (2016).

29. 29.

Polettini, M., Lazarescu, A. & Esposito, M. Tightening the uncertainty principle for stochastic currents. Phys. Rev. E 94, 052104 (2016).

30. 30.

Pietzonka, P., Ritort, F. & Seifert, U. Finite-time generalization of the thermodynamic uncertainty relation. Phys. Rev. E 96, 012101 (2017).

31. 31.

Horowitz, J. M. & Gingrich, T. R. Thermodynamic uncertainty relations constrain non-equilibrium fluctuations. Nat. Phys. 16, 15–20 (2020).

32. 32.

Saryal, S., Friedman, H. M., Segal, D. & Agarwalla, B. K. Thermodynamic uncertainty relation in thermal transport. Phys. Rev. E 100, 042101 (2019).

33. 33.

Timpanaro, A. M., Guarnieri, G., Goold, J. & Landi, G. T. Thermodynamic uncertainty relations from exchange fluctuation theorems. Phys. Rev. Lett. 123, 090604 (2019).

34. 34.

Hasegawa, Y. & Van Vu, T. Fluctuation theorem uncertainty relation. Phys. Rev. Lett. 123, 110602 (2019).

35. 35.

Li, J., Horowitz, J. M., Gingrich, T. R. & Fakhri, N. Quantifying dissipation using fluctuating currents. Nat. Commun. 10, 1666 (2019).

36. 36.

Dechant, A. & Sasa, S.-I. Fluctuation–response inequality out of equilibrium. Proc. Natl Acad. Sci. USA 117, 6430 (2020).

37. 37.

Pal, S., Saryal, S., Segal, D., Mahesh, T. S. & Agarwalla, B. K. Experimental study of the thermodynamic uncertainty relation. Phys. Rev. Res. 2, 022044(R) (2020).

38. 38.

Wei, X., Sheng, J., Yang, C., Wu, Y. & Wu, H. Controllable two-membrane-in-the-middle cavity optomechanical system. Phys. Rev. A 99, 023851 (2019).

39. 39.

Sheng, J., Wei, X., Yang, C. & Wu, H. Self-organized synchronization of phonon lasers. Phys. Rev. Lett. 124, 053604 (2020).

40. 40.

Wu, S., Sheng, J., Zhang, X., Wu, Y. & Wu, H. Parametric excitation of a SiN membrane via piezoelectricity. AIP Adv. 8, 015209 (2018).

41. 41.

Aspelmeyer, M., Kippenberg, T. J. & Marquardt, F. Cavity optomechanics. Rev. Mod. Phys. 86, 1391 (2014).

42. 42.

Marquardt, F., Chen, J. P., Clerk, A. A. & Girvin, S. M. Quantum theory of cavity-assisted sideband cooling of mechanical motion. Phys. Rev. Lett. 99, 093902 (2007).

43. 43.

Groblacher, S., Hammerer, K., Vanner, M. R. & Aspelmeyer, M. Observation of strong coupling between a micromechanical resonator and an optical cavity field. Nature 460, 724–727 (2009).

44. 44.

Barton, G. Classical van der Waals heat flow between oscillators and between half-spaces. J. Phys. Condens. Matter. 27, 214005 (2015).

45. 45.

Bérut, A., Imparato, A., Petrosyan, A. & Ciliberto, S. Stationary and transient fluctuation theorems for effective heat fluxes between hydrodynamically coupled particles in optical traps. Phys. Rev. Lett. 116, 068301 (2016).

46. 46.

Nielsen, W. H. P., Tsaturyan, Y., Møller, C. B., Polzik, E. S. & Schliesser, A. Multimode optomechanical system in the quantum regime. Proc. Natl Acad. Sci. USA 114, 62 (2017).

## Acknowledgements

We thank Jordan M. Horowitz for helpful discussions. This research was supported by the National Key R&D Program of China (No. 2017YFA0304201), NSFC (No. 11925401, No. 11734008, No. 11974115, and No. 11704126, No. 11621404), the Shanghai Committee of Science and Technology (No. 17JC1400500), the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning.

## Author information

Authors

### Contributions

C.Y., X.W., J.S., and H.W. carried out the experiment, analyzed the data, and developed the theory. J.S. and H.W. wrote the paper with input from all authors.

### Corresponding authors

Correspondence to Jiteng Sheng or Haibin Wu.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Peer review information Nature Communications thanks Seunghwi Kim and the other anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Yang, C., Wei, X., Sheng, J. et al. Phonon heat transport in cavity-mediated optomechanical nanoresonators. Nat Commun 11, 4656 (2020). https://doi.org/10.1038/s41467-020-18426-4

• Accepted:

• Published: