Abstract
Engineering phonon transport in physical systems is a subject of interest in the study of materials, and has a crucial role in controlling energy and heat transfer. Of particular interest are non-reciprocal phononic systems, which in direct analogy to electric diodes, provide a directional flow of energy. Here, we propose an engineered nanostructured material, in which tunable non-reciprocal phonon transport is achieved through optomechanical coupling. Our scheme relies on breaking time-reversal symmetry by a spatially varying laser drive, which manipulates low-energy acoustic phonons. Furthermore, we take advantage of developments in the manipulation of high-energy phonons through controlled scattering mechanisms, such as using alloys and introducing disorder. These combined approaches allow us to design an acoustic isolator and a thermal diode. Our proposed device will have potential impact in phonon-based information processing, and heat management in low temperatures.
Similar content being viewed by others
Introduction
Controlling the flow of heat is important for several fields including thermoelectrics, thermal management, and information processing. For example, suppressing thermal conductivity can improve the performance of thermoelectrics, and can also isolate circuit elements from external heat. The thermal analog of an electric diode is of fundamental importance to efforts in managing heat. Thermal diodes have numerous applications, including blocking unwanted back scattering in phonon-based information processing as well as managing heat and maximizing efficiency in nanostructures. The operation of a thermal diode requires a nonlinear material or broken time-reversal symmetry1; most implementations have exploited the former2,3,4,5,6.
Theoretical and experimental advances in our understanding of the contribution of coherent phonons to heat transport has provided new insights that allow for enhanced control of heat flow in nanostructured materials7. Specifically, due to the very long mean free paths and coherence of these phonons8,9, periodic structures can modify their dispersion and transport properties10. Thus, engineered bandgaps, which have been used to manipulate sound11, can also be used to alter the thermal properties of a material12. Moreover, adding impurities have been proven useful in modifying thermal transport properties of a material by manipulating high-energy phonons13,14.
At the same time, there have been remarkable advances in cavity optomechanics15, where interactions between photons and acoustic phonons confined in an optomechanical cavity can be controlled at the single phonon level16, with potential applications in quantum information processing17,18. More recently, non-reciprocal optical transport was proposed in ring resonators, where the directional laser pump selects one circulation direction19. This scheme and approaches based on stimulated Brillouin scattering in photonics20,21 were experimentally demonstrated in multiple optomechanical systems22,23,24,25,26. Meanwhile, the resulting chirality for phonons in such ring resonators has been investigated27,28,29. Moreover, there have been intriguing proposals to synthesize gauge fields in optomechanical systems, from photonic crystals30,31, to quantum wells32, and superconducting circuits33, and to investigate their associated topological properties34.
In this article, we combine the physics of heat transport in nanostructures and optomechanics to develop a new platform to manipulate both low-energy and high-energy phonons. We propose a method to engineer a tunable non-reciprocal bandgap for acoustic phonons, where a laser field with a phase gradient optically drives an array of optomechanical cavities and induces the non-reciprocal transport by breaking the time-reversal symmetry. We propose an experimental implementation of the scheme in optomechanical crystals. We discuss two applications of such a system, one as an acoustic isolator, and the other as a thermal diode. For the latter, the introduction of alloy and nanoparticle disorder suppresses the transport of high-energy phonons, leaving the low-energy acoustic band as the dominant channel for heat conduction10,35, thus enhancing the overall optomechanically induced non-reciprocity. Our proposed device works in the linear regime and introduces an alternative approach to previous works.
Results
Tight-binding model
To illustrate the basic concepts, we study non-reciprocal transport in a system connected to heat baths. The system considered here is described by a tight-binding model of an array of coupled optomechanical cavities36,37,38. An optomechanical cavity supports localized electromagnetic and mechanical modes. Owing to the radiation pressure, changes to the shape of the cavity change its electromagnetic resonance frequency, effectively coupling mechanical vibrations to electromagnetic excitations. The Hamiltonian describing a collection of isolated cavities (setting ħ = 1) is
where \(\hat a_n(\hat b_n)\) is a bosonic operator that destroys a photonic (phononic) excitation with energy ωcav(mech) at site n, and g is the vacuum coupling rate. In addition, there is a loss rate γcav(γmech) associated with the optical (mechanical) mode of the cavity (see Methods).
In a linear array of cavities, nearest-neighbor couplings dominates due to the tunneling of excitations between adjacent sites. The Hamiltonian describing these processes is
where h.c. denotes Hermitian conjugate, and t and J are tunneling strengths of phononic and photonic excitations, respectively. The system Hamiltonian is then given by
The vacuum coupling rate, g, is typically small, and can be enhanced by means of an external laser drive. To break reciprocity in a spatially dependent manner, in contrast to the usual setup in which the optical mode on each site is excited by a laser with a uniform phase15, we consider a phase gradient in the laser field. This phase, which breaks time-reversal symmetry introduces a position dependent phase in the effective coupling between phonons and photons31. The breaking of time-reversal symmetry is crucial for the effects we consider in this work. The laser frequency, ωd = ωcav + Δ, is detuned from the resonance frequency of the cavity by Δ, and the phase offset between adjacent sites is θ, as shown in Fig. 1b.
To bring electromagnetic and mechanical excitations on resonance, the driving laser is red-detuned from the cavity resonance frequency by Δ ≈ −ωmech. In a rotating frame of photons with angular frequency ωd, after making the rotating-wave approximation (RWA) in the resolved sideband regime (\(\omega _{{\mathrm{mech}}} \gg \gamma _{{\mathrm{mech}}}\)), linearizing and displacing the cavity field, the effective Hamiltonian is39
where G = αg, the enhanced optomechanical coupling strength is large compared to g by an order of α, the square root of the number of photons in the cavity.
The propagating modes of the system are polaritons, which are superpositions of electromagnetic and mechanical quanta. To diagonalize the Hamiltonian in Eq. (4), we write the Fourier transform for \(\hat a_n\) and \(\hat b_n\) as
where d0 is the lattice constant. Then Eq. (4) in the Fourier basis is
which shows the phase gradient of the laser field acts as a momentum shift, and leads to the coupling of phonons and photons with different momenta.
The band structure is shown in Fig. 1c. In the absence of coupling the phononic and photonic dispersions intersect. For G ≠ 0, two bandgaps develop. The asymmetry (under k → −k) of the band structure is controlled by θ. The eigenmodes are polaritons, \(\alpha _{k,j}\hat a_{k - \theta } + \beta _{k,j}\hat b_k\), where j \(\in\) {1, 2} is the band index. The quantities |αk,j|2 and |βk,j|2 indicate the relative weights of photons and phonons composing the polaritons, respectively. The effect of cavity loss is to broaden the bands, and for the gaps to be effective, we need to be in the high cooperativity regime, i.e., \(G^2{\mathrm{/}}\gamma _{{\mathrm{cav}}}\gamma _{{\mathrm{mech}}} \gg 1\).
The asymmetry of the gaps controls the non-reciprocal transport properties of the system, as shown Fig. 1d. To study thermal transport properties of this model, we consider connecting the system to two thermal contacts. The contacts are impedance matched to the non-driven (G = 0) system. Thus, the dispersion of phonons in the contacts is ωcontact(k) = ωmech − 2t cos(d0k). The transmission probabilities can be calculated by mode matching at the boundaries of the system; see Fig. 1d and the Supplementary Note 1. This continuum picture remains valid for a finite number of lattice sites and the transmission probabilities are close to zero for phonons with energies in the gap. These phonons are converted to photons and reflected. The probabilities exhibit Fabry–Perot oscillations whose period is proportional to the inverse of the number of sites in the system. The direction dependent phonon transmission probability for θ close to π can be approximated by
where H(ω) is the Heaviside step function.
The phonon thermal current, assuming that the photon contacts are at zero temperature, can be calculated in the Landauer–Büttiker formalism40,41, and is given by
where ΘL(R) is the temperature of the left(right) contact, and nB(Θ, ω) = 1/(exp(ħω/kBΘ) − 1). We introduce an alternative set of variables to denote the mean temperature, Θ0, and the temperature bias, ΔΘ, such that
with
This relationship implies that if Θ0 is fixed and only the sign of ΔΘ is changed, the temperatures of the two contacts are swapped. In a reciprocal system, there is no distinction between left and right, and taking ΔΘ → −ΔΘ only changes the sign of the current and leaves the magnitude unchanged. However, due to the broken time-reversal symmetry in our system, transport is non-reciprocal \({\cal T}_{{\mathrm{L}} \to {\mathrm{R}}}(\omega ) \ne {\cal T}_{{\mathrm{L}} \leftarrow {\mathrm{R}}}(\omega )\), and the current magnitudes are different. Figure 1e shows the current I as a function of temperature bias ΔΘ. The base temperature Θ0 is chosen close to the energy scale ωmech of the system, so that the non-reciprocal effect is enhanced. In this plot, the non-zero intercept (I0 ∝ G2) is a measure of the non-reciprocity. This is different from an electrical diode mechanism, where the slope changes as the bias is reversed.
To quantify the non-reciprocity, we introduce the contrast C, defined as
which is non-zero in a non-reciprocal system. In the shaded region in Fig. 1e and f, the current does not change its direction when the bias is reversed. In this case, the contrast is maximized and C = 1.
The relation between I and ΔΘ for \(\hbar \omega _{{\mathrm{mech}}}{\mathrm{/}}k_{\mathrm{B}}{\Theta} _0 \ll 1\), and \({\mathrm{\Delta }}{\Theta} \ll {\Theta} _0\) is well described by
In the same regime, the contrast is given by
where \({\mathrm{\Delta }}{\Theta} _{\mathrm{c}} = \hbar G^2{\mathrm{/}}[k_{\mathrm{B}}(2t - G)]\). These approximations are compared with the exact values in Fig. 1e and f.
This non-reciprocal model can be implemented in an optomechanical crystal42,43. An optomechanical crystal is an engineered dielectric, which supports localized phononic and photonic excitations with energies in the bandgaps. Given a uniform dielectric, bandgaps can be introduced by drilling a periodic array of identical holes. Deforming these holes to form a superlattice introduces defect cavities which co-localize phonons and photons, thereby enhancing their mutual couplings. In Fig. 2, we show the correspondence between a cavity and its implementation in the actual optomechanical crystal. Each unit cell is a few microns in size, and a total system size of hundreds of microns leads to the non-reciprocity shown in Fig. 1d. The bare optomechanical coupling g varies between several kilohertz and tens of megahertz in various materials such as Si or GaAs26,44,45. The tunneling strengths depend on the structure design, and values of a few megahertz for the mechanical tunneling strength t, and hundreds of megahertz for its optical counterpart, have been realized in the experiments26. In Methods and Supplementary Note 2, we present more details, and specifically show how to engineer the non-reciprocal band in an optomechanical crystal. The tight-binding model is applicable not only to optomechanical crystal arrays, but also to other optomechanical systems such as coupled ring resonators that have been realized in experiments46.
Applications
Now that we have established that non-reciprocal transport for a continuum band of phonons can be achieved in an array of optomechanical cavities, we further discuss two applications in an optomechanical crystal (see Methods): (1) An acoustic isolator, and (2) a thermal diode in low temperatures.
An acoustic isolator is a device that only permits propagation of coherent monochromatic phonons in one direction. Such a device can be realized using a nonlinear medium attached to a phononic crystal47,48, or through spatio-temporal modulation of system properties in a transmission line49. Our optomechanical crystal operates as an isolator for frequencies which lie in the bandgap. An appropriate figure of merit analogous to the contrast in Eq. (12) for monochromatic waves is
The device discussed in the previous section, with the same red-detuned laser drive acts as an isolator for phonons with frequencies close to ωmech. Specifically, in the high cooperativity regime, the effect of loss is negligible, and we can use the transmission probabilities shown in Fig. 1d. We see that for propagating elastic waves with frequencies inside the bandgap, Ciso(ω) approaches unity, and otherwise, is very close to zero, thus realizing an isolator with a bandwidth of 2G. The non-reciprocity in our scheme is tunable, and the frequency range of the gap is also controllable and depends on the phase gradient of the laser field. Furthermore, as the Hamiltonian in Eq. (4) is linear in both the optical and mechanical fields, in principle the device works at the quantum limit. It is therefore useful for quantum information routing18, and may find new applications to hybrid devices such as superconducting qubits50 coupled to optomechanical crystals17, as they both work in the same energy regime.
As a second application, we show that our system can serve as a thermal diode. A perfect thermal diode would allow heat transport in only one direction. Our system relies on the modification of the material properties in a narrow frequency range at low energies, whereas a major part of heat current is carried out by high-frequency phonons. To suppress the contribution of these high-energy phonons, we introduce various scattering mechanisms to shorten their mean free paths10.
Specifically, to evaluate the figure of merit C in a realistic material, we analyze frequency-dependent phonon scattering processes characterized by a length scale λph(ω). The total transmission at a given frequency, summed over all bands, can be approximated as
where the factor \(\frac{{\lambda _{{\mathrm{ph}}}}}{{\lambda _{{\mathrm{ph}}} + L_{\mathrm{s}}}}\) is the probability of transmittance51,52, and ML⇄R(ω) is the number of conducting bands at a given energy for right (L → R), or left (L ← R) moving phonons, and ϕ is the sample’s porosity. The function \(f(\phi ) = \frac{{1 - \phi }}{{1 + \phi }}\) comes from Maxwell–Garnett effective medium approach, and takes the effect of holes on the number of modes into account53,54. The parameter Ls is the length of the sample, and λph is the mean free path of back scattering, which is related to the mean free path of scattering, lph, by λph = 2lph in 1D, and λph = 4/3lph in 3D. Note that because of the dependence of λph on frequency, the performance of the device gets better for larger samples. Specifically, for larger samples the overall transmission decreases, however, as the mean free path is smaller for higher frequencies, the contrast improves. In our calculations, we considered N = 100 sites.
Following refs. 10,55, we propose using alloy and nano-paricle impurities to modify λph. In this case, the optomechanical crystal is made of an alloy into which nanoparticles are embedded. These impurities lead to mass-difference scatterings, and their respective rates \(\tau _{{\mathrm{alloy}}}^{ - 1}\) and \(\tau _{{\mathrm{np}}}^{ - 1}\) scales with ω4, thus lowering the contribution of high-energy phonons and increasing the contribution of low-energy phonons to thermal transport. To obtain a realistic estimate of λph, it is necessary to consider two additional scattering mechanisms: intrinsic anharmonicity and boundary scattering, characterized by rates \(\tau _{{\mathrm{an}}}^{ - 1}\) and \(\tau _{\mathrm{b}}^{ - 1}\), respectively. The corresponding mean free paths l i ’s, are obtained from the scattering rates by l i = vgτ i , where vg is the group velocity. The total mean free path is given by Matthiessen’s rule, i.e.,
The effect of these scattering mechanisms on the cumulative thermal current, shown in Fig. 3a, is evaluated by a hybrid method using the bulk silicon dispersion (ω ∝ k) for high-energy phonons with short mean free paths and the superlattice dispersion for lower energy phonons with mean free paths longer than several lattice constants of the superlattice54 (see Methods). It can be seen that phonons with frequencies below 25 GHz contribute more to heat transport in a Si90Ge10 optomechanical crystal with 10 nm nanoparticles and a filling factor of 5% than to a nanobeam of the same dimensions composed of nonporous silicon. Moreover, the ratio of the current carried by the optomechanically coupled band to the total current in the engineered optomechanical crystal is close to 9% at 4 K and about 22% at 0.4 K with a bias of ΔΘ/Θ0 = 10−3 (see Supplementary Note 3).
Finally, we calculate the contrast for a driven optomechanical crystal, as displayed in Fig. 3b. As the base temperature decreases and approaches the energy of the optomechanical band (kBΘ0 ≈ ħωmech), the contrast and the non-reciprocal effect increase. Although a significant constrast can be achieved at sub-Kelvin temperatures, the generalization of this scheme to room temperature requires a significant improvement in the material properties such as the optomechanical coupling strength (see Supplementary Note 3 for room temperature). Although there is intrinsic photon loss in the silicon beam, such loss does not directly lead to the generation of phonons in the beam56, and therefore, the implicit assumptions of our approach remain valid (see Methods and Supplementary Note 3).
To measure the contrast, we envision a setup similar to ref. 57, where a pair superconductor/insulator/normal metal/insulator/superconductor (SINIS) tunnel junction58 are mounted at the two ends of the device and serve as both a sensitive thermometer and a heater. The device is heated from one side and the change in temperature is measured at the other side. By interchanging the role of the heater and thermometer and comparing the measurements, the existence of non-reciprocal thermal current in the system can be verified.
Discussion
In this work, we have shown that phase-modulated driven optomechanical systems can be utilized to engineer a non-reciprocal phonon band in a material. We discussed two possible applications of our scheme, as an acoustic isolator and a thermal diode. Although we considered a specific silicon-based optomechanical crystal, these methods can be readily generalized to other materials and designs. These devices may find application to on-chip heat management and quantum information processing, both to increase coherence time and to exploit phonons as information carriers59.
In our approach, we proposed a coherent dynamical control of low-energy phonons by using optomechanical structures, and combined it with incoherent control of high-energy phonons by designing bulk material properties through the introduction of disorder. This strategy of combing low-energy and high-energy phononic physics could be generalized to designing of other thermal technologies such as thermoelectrics, thermal insulation, and the development of new metamaterials.
Methods
Implementation with optomechanical crystals
In an optomechanical cavity electromagnetic and mechanical modes are co-localized. The coupling between these modes arises due to radiation pressure, which changes the cavity’s electromagnetic resonance frequency. More rigorously, the energy \(\hbar \omega (\hat x)\) of a cavity depends on its shape, where \(\hat x = x_{{\mathrm{ZPF}}}(\hat b + \hat b^{\mathrm{\dagger }})\) is the quantized mechanical displacement, and xZPF denotes zero point fluctuations. The Hamiltonian describing a single defect cavity is
The optomechanical coupling arises as the shape of the dielectric boundary changes. Expanding \(\omega (\hat x)\) to first order in \(\hat x\) results in
where g can be calculated from moving boundaries perturbation theory for Maxwell’s equations60,61. Equation (19) reproduces the form of \(\hat H_{{\mathrm{sites}}}\) (1).
Imperfect localization within each cavity leads to nearest-neighbor hopping of the phonons and photons as captured by Htunneling (4). Using finite-element simulations (FEM), we find that to a very good approximation these bands follow the dispersion ω ∝ cos(k) as predicted by the tight-binding model (see Supplementary Note 2). Typical values of ωcav, and ωmech for a silicon optomechanical crystal are 100 THz, and 10 GHz, respectively26. The laser frequency ωd should be ωcav − ωmech ≈ 100 THz to bring it in resonance with the mechanical mode.
To linearize the Hamiltonian \(\hat H_{{\mathrm{sites}}} + \hat H_{{\mathrm{tunneling}}}\) with the laser drive \(\mathop {\sum}\nolimits_n {\kern 1pt} \epsilon _{\mathrm{d}}e^{i\theta n}{\kern 1pt} {\mathrm{cos}}(\omega _{\mathrm{d}}t)\,(a_n^{\mathrm{\dagger }} + a_n)\), we solve for the steady-steady of the cavity in the absence of the optomechanical coupling (g = 0)15, and find that the steady-state is given by
where |α|2 is the number of photons in the cavity. Consequently, displacing the cavity field by \(\hat a_n \to \hat a_n + \alpha e^{i\theta n}\), and using RWA results in \(\hat H_{{\mathrm{eff}}}\) (4).
The phases (einθ) can be tuned off the chip by using stretchable fiber phase shifters26. An on-chip implementation is possible by using a 1 × N multi-mode interferometer to divide the power, and meandered waveguides or zero-loss resonators to tune the phase (see Supplementary Note 4).
Scattering rates
The scattering rate associated with the alloy disorder of Si x Ge1−x is described by an effective mass-difference Rayleigh scattering62,63
where A = 3.01 × 10−41 s5 10,55 is a constant that depends on the alloy properties. Scattering due to nanoparticles can be described by interpolating between the long- and short-wavelength scattering regimes55,64,
where f and V = 4πr3/3 are the filling fraction and the volume of nanoparticles, vg is the magnitude of the group velocity of phonons, and
Here, r is the radius of nanoparticles, ΔD is the difference between particle and alloy densities, and D is the alloy’s density. In our calculations, we consider a filling fraction of f = 0.05, σs = 6.28 × 10−16 m2, and σl = 2.2 × 10−48 m6 × (ω/vg)4, corresponding to r = 10 nm germanium nanoparticles in Si90Ge10.
The anharmonic scattering rate, which takes both the normal and umklapp processes into account, is given by
where B(T) = 3.28 × 10−19 s K−1 and C = 140 K for Si90Ge1010,55. Scattering from the boundaries of a thin film can be modeled by \(l_{\mathrm{b}} = \frac{{1 + p}}{{1 - p}}t\)65,66, where t, the thickness of the sample, is the mean free path in the diffusive limit. The parameter p is the probability that the scattering is specular. It takes the effect of surface roughness into account and depends on the phonon’s wavelength. It is given by65,67
where Λ is the wavelength of the phonons, and η is the surface roughness of the sample, which is taken to be 1 nm in our calculations, and corresponds to an estimated surface roughness achieved in silicon thin film fabrications51.
Thermal current
In Fig. 3, we compared the operation of our optomechanical crystal with a nanobeam of nonporous silicon. The thermal current in the latter can be calculated as follows. At the temperatures considered, only phonons with frequencies smaller than 3 THz are relevant for thermal transport. In this frequency range, only the acoustic branch contributes to thermal conductivity and the dispersion is approximately linear. Specifically, we employed the Debye dispersion, i.e., ω = cs|k|, where cs is the sound velocity, and k is the wavevector. The density of modes is given by \(M(\omega ) = S{\kern 1pt} 3\omega ^2{\mathrm{/}}4\pi c_{\mathrm{s}}^2\), where S is the cross sectional area. In addition, because the sample in this case is nonporous, f(ϕ) = 1, and the only scattering mechanisms are scatterings due to surface roughness and crystal anharmonicities.
The thermal current for the proposed Si90Ge10 optomechanical crystal with nanoparticles is calculated using a hybrid method, depending on the frequency of phonons54. Specifically, The mean free path of phonons depends on their frequency. As the frequency of phonons is lowered, their mean free path becomes comparable with the superlattice spacing, and therefore, the bulk dispersion is no longer a good description. In our system this threshold frequency corresponds to 25 GHz. To get this number, using the silicon bulk dispersion and Eqs. (17) and (21–26), we find the frequency for which the mean free path is comparable to several lattice spacings.
For phonons with frequencies above the threshold, we use the bulk Debye dispersion and associated group velocities and density of states. For phonons with frequencies lower than this threshold, we use the superlattice dispersion, calculated by FEM simulations, which gives group velocities and density of states different from those of the bulk. The scattering rates are then calculated using the superlattice dispersion in this regime. Taken together, these account for the total reciprocal phonon contribution to the thermal current. Finally, we add the contribution of the single non-reciprocal optomechanically coupled band to the calculated current. For this single band, we have used θ = 1.3π and parameters ωmech/2π = 4.3 GHz, J/2π = 0.5 GHz, t/2π = 0.2 GHz, and G/2π = 0.1 GHz, which leads to asymmetric gaps in 3.96 GHz < ω/2π < 4.13 GHz, and 4.47 GHz < ω/2π < 4.64 GHz for right-going, and left-going phonons, respectively.
Data availability
The data that support the findings of this study are available from the authors upon reasonable request.
References
Maznev, A., Every, A. & Wright, O. Reciprocity in reflection and transmission, what is a ‘phonon diode’? Wave Motion 50, 776–784 (2013).
Terraneo, M., Peyrard, M. & Casati, G. Controlling the energy flow in nonlinear lattices: a model for a thermal rectifier. Phys. Rev. Lett. 88, 094302 (2002).
Li, B., Wang, L. & Casati, G. Thermal diode: rectification of heat flux. Phys. Rev. Lett. 93, 184301 (2004).
Segal, D. & Nitzan, A. Spin-boson thermal rectifier. Phys. Rev. Lett. 94, 034301 (2005).
Chang, C., Okawa, D., Majumdar, A. & Zettl, A. Solid-state thermal rectifier. Science 314, 1121–1124 (2006).
Yang, N., Li, N., Wang, L. & Li, B. Thermal rectification and negative differential thermal resistance in lattices with mass gradient. Phys. Rev. B 76, 020301 (2007).
Maldovan, M. Phonon wave interference and thermal bandgap materials. Nat. Mater. 14, 667–674 (2015).
Luckyanova, M. N. et al. Coherent phonon heat conduction in superlattices. Science 338, 936–939 (2012).
Esfarjani, K., Chen, G. & Stokes, H. T. Heat transport in silicon from first-principles calculations. Phys. Rev. B 84, 085204 (2011).
Maldovan, M. Narrow low-frequency spectrum and heat management by thermocrystals. Phys. Rev. Lett. 110, 025902 (2013).
Li, N. et al. Colloquium: phononics: manipulating heat flow with electronic analogs and beyond. Rev. Mod. Phys. 84, 1045–1066 (2012).
Hopkins, P. E. et al. Reduction in the thermal conductivity of single crystalline silicon by phononic crystal patterning. NanoLetters 11, 107–112 (2010).
Kim, W. et al. Thermal conductivity reduction and thermoelectric figure of merit increase by embedding nanoparticles in crystalline semiconductors. Phys. Rev. Lett. 96, 045901 (2006).
Snyder, G. J. & Toberer, E. S. Complex thermoelectric materials. Nat. Mater. 7, 105–114 (2008).
Aspelmeyer, M., Kippenberg, T. J. & Marquardt, F. Cavity optomechanics. Rev. Mod. Phys. 86, 1391–1452 (2014).
Chan, J. et al. Laser cooling of a nanomechanical oscillator into its quantum ground state. Nature 478, 89–92 (2011).
Stannigel, K. et al. Optomechanical quantum information processing with photons and phonons. Phys. Rev. Lett. 109, 013603 (2012).
Habraken, S., Stannigel, K., Lukin, M. D., Zoller, P. & Rabl, P. Continuous mode cooling and phonon routers for phononic quantum networks. New J. Phys. 14, 115004 (2012).
Hafezi, M. & Rabl, P. Optomechanically induced non-reciprocity in microring resonators. Opt. Express 20, 7672–7684 (2012).
Huang, X. & Fan, S. Complete all-optical silica fiber isolator via stimulated Brillouin scattering. J. Light Technol. 29, 2267–2275 (2011).
Kang, M. S., Butsch, A. & Russell, P. S. J. Reconfigurable light-driven opto-acoustic isolators in photonic crystal fibre. Nat. Photonics 5, 549–553 (2011).
Kim, J., Kuzyk, M. C., Han, K., Wang, H. & Bahl, G. Non-reciprocal Brillouin scattering induced transparency. Nat. Phys. 11, 275–280 (2015).
Dong, C.-H. et al. Brillouin-scattering-induced transparency and non-reciprocal light storage. Nat. Commun. 6, 6193 (2015).
Shen, Z. et al. Experimental realization of optomechanically induced non-reciprocity. Nat. Photonics 10, 657–661 (2016).
Ruesink, F., Miri, M.-A., Alù, A. & Verhagen, E. Nonreciprocity and magnetic-free isolation based on optomechanical interactions. Nat. Commun. 7, 13662 (2016).
Fang, K. et al. Generalized non-reciprocity in an optomechanical circuit via synthetic magnetism and reservoir engineering. Nat. Phys. 13, 465–471 (2017).
Fleury, R., Sounas, D. L., Sieck, C. F., Haberman, M. R. & Alù, A. Sound isolation and giant linear non-reciprocity in a compact acoustic circulator. Science 343, 516–519 (2014).
Kim, S., Xu, X., Taylor, J. M. & Bahl, G. Dynamically induced robust phonon transport and chiral cooling in an optomechanical system. Nat. Commun. 8, 205 (2017).
Xu, X. & Taylor, J. M. Optomechanically induced chiral transport of phonons in one dimension. Preprint at http://arxiv.org/abs/1701.02699 (2017).
Ludwig, M. & Marquardt, F. Quantum many-body dynamics in optomechanical arrays. Phys. Rev. Lett. 111, 073603 (2013).
Schmidt, M., Kessler, S., Peano, V., Painter, O. & Marquardt, F. Optomechanical creation of magnetic fields for photons on a lattice. Optica 2, 635–641 (2015).
Poshakinskiy, A. V. & Poddubny, A. N. Phonoritonic crystals with a synthetic magnetic field for an acoustic diode. Phys. Rev. Lett. 118, 156801 (2017).
Barzanjeh, S., Aquilina, M. & Xuereb, A. Manipulating the flow of thermal noise in quantum devices. Phys. Rev. Lett. 120, 060601 (2018).
Peano, V., Brendel, C., Schmidt, M. & Marquardt, F. Topological phases of sound and light. Phys. Rev. X 5, 031011 (2015).
Garg, J., Bonini, N., Kozinsky, B. & Marzari, N. Role of disorder and anharmonicity in the thermal conductivity of silicon-germanium alloys: a first-principles study. Phys. Rev. Lett. 106, 045901 (2011).
Chang, D., Safavi-Naeini, A. H., Hafezi, M. & Painter, O. Slowing and stopping light using an optomechanical crystal array. New J. Phys. 13, 023003 (2011).
Chen, W. & Clerk, A. A. Photon propagation in a one-dimensional optomechanical lattice. Phys. Rev. A 89, 033854 (2014).
Marquardt, F. & Rakich, P. T. Quantum theory of continuum optomechanics. New J. Phys. https://doi.org/10.1088/1367-2630/aaac4f (2018).
Aspelmeyer, M., Kippenberg, T. & Marquardt, F. Cavity Optomechanics: Nano- and Micromechanical Resonators Interacting with Light, 1st edn (Springer, Berlin, 2014).
Yang, N., Xu, X., Zhang, G. & Li, B. Thermal transport in nanostructures. AIP Adv. 2, 041410 (2012).
Datta, S., Bagwell, P. F. & Anantram, M. Scattering theory of transport for mesoscopic superconductors. Phys. Chem. Mater. Low Dimens. Struct. 3, 1 (1996).
Chan, J. Laser Cooling of An Optomechanical Crystal Resonator to its Quantum Ground State of Motion. PhD Dissertation, California Institute of Technology (2012).
Safavi-Naeini, A. H. & Painter, O (eds) in book, Cavity Optomechanics 195–231 (Springer, Berlin, 2014).
Balram, K. C., Davanço, M., Lim, J. Y., Song, J. D. & Srinivasan, K. Moving boundary and photoelastic coupling in GaAs optomechanical resonators. Optica 1, 414–420 (2014).
Leijssen, R., La Gala, G. R., Freisem, L., Muhonen, J. T. & Verhagen, E. Nonlinear cavity optomechanics with nanomechanical thermal fluctuations. Nat. Commun. 8, 16024 (2017).
Gil-Santos, E. et al. Light-mediated cascaded locking of multiple nano-optomechanical oscillators. Phys. Rev. Lett. 118, 063605 (2017).
Liang, B., Yuan, B. & Cheng, J.-C. Acoustic diode: rectification of acoustic energy flux in one-dimensional systems. Phys. Rev. Lett. 103, 104301 (2009).
Liang, B., Guo, X., Tu, J., Zhang, D. & Cheng, J. An acoustic rectifier. Nat. Mater. 9, 989–992 (2010).
Zanjani, M. B., Davoyan, A. R., Engheta, N. & Lukes, J. R. NEMS with broken T symmetry: graphene based unidirectional acoustic transmission lines. Sci. Rep. 5, 9926 (2015).
Clarke, J. & Wilhelm, F. K. Superconducting quantum bits. Nature 453, 1031–1042 (2008).
Jeong, C., Datta, S. & Lundstrom, M. Thermal conductivity of bulk and thin-film silicon: a landauer approach. J. Appl. Phys. 111, 093708 (2012).
Datta, S. Electronic Transport in Mesoscopic Systems 1st edn (Cambridge University Press, Cambridge, 1997).
Nan, C.-W., Birringer, R., Clarke, D. R. & Gleiter, H. Effective thermal conductivity of particulate composites with interfacial thermal resistance. J. Appl. Phys. 81, 6692–6699 (1997).
Alaie, S. et al. Thermal transport in phononic crystals and the observation of coherent phonon scattering at room temperature. Nat. Commun. 6, 7228 (2015).
Mingo, N., Hauser, D., Kobayashi, N., Plissonnier, M. & Shakouri, A. “Nanoparticle-in-alloy” approach to efficient thermoelectrics: silicides in SiGe. NanoLetters 9, 711–715 (2009).
Meenehan, S. M. et al. Silicon optomechanical crystal resonator at millikelvin temperatures. Phys. Rev. A 90, 011803 (2014).
Zen, N., Puurtinen, T. A., Isotalo, T. J., Chaudhuri, S. & Maasilta, I. J. Engineering thermal conductance using a two-dimensional phononic crystal. Nat. Commun. 5, 3435 (2014).
Giazotto, F., Heikkilä, T. T., Luukanen, A., Savin, A. M. & Pekola, J. P. Opportunities for mesoscopics in thermometry and refrigeration: physics and applications. Rev. Mod. Phys. 78, 217–274 (2006).
Chu, Y. et al. Quantum acoustics with superconducting qubits. Science 358, 199–202 (2017).
Johnson, S. G. et al. Perturbation theory for Maxwell’s equations with shifting material boundaries. Phys. Rev. E 65, 066611 (2002).
Eichenfield, M., Chan, J., Camacho, R. M., Vahala, K. J. & Painter, O. Optomechanical crystals. Nature 462, 78–82 (2009).
Klemens, P. The scattering of low-frequency lattice waves by static imperfections. Proc. Phys. Soc. Lond. Sect. A 68, 1113–1128 (1955).
Abeles, B. Lattice thermal conductivity of disordered semiconductor alloys at high temperatures. Phys. Rev. 131, 1906–1911 (1963).
Kim, W. & Majumdar, A. Phonon scattering cross section of polydispersed spherical nanoparticles. J. Appl. Phys. 99, 084306 (2006).
Ziman, J. M. Electrons and Phonons: The Theory of Transport Phenomena in Solids 1st edn (Oxford University Press, Oxford, 2001).
Sondheimer, E. H. The mean free path of electrons in metals. Adv. Phys. 1, 1–42 (1952).
Zhang, Z. Nano/Microscale Heat Transfer 1st edn (McGraw-Hill Education, New York, 2007).
Acknowledgements
We thank Krishna Balram for providing the initial FEM simulations, Hirokazu Miyake for helping with the simulations, and Reza Ghodssi for providing access to computational resources. We also thank Sunil Mittal, Raphaël Van Laer, Amir Safavi-Naeini, and Oskar Painter for helpful discussions. The work was partially supported by Sloan Fellowship, YIP-ONR, the NSF PFC at the JQI.
Author information
Authors and Affiliations
Contributions
A.S. and W.D. formulated the transmission probabilities in the tight-binding model. A.S. performed the calculation and numerical simulations. K.E. and M.H. directed the study. All authors discussed and analyzed the results and contributed and commented on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Seif, A., DeGottardi, W., Esfarjani, K. et al. Thermal management and non-reciprocal control of phonon flow via optomechanics. Nat Commun 9, 1207 (2018). https://doi.org/10.1038/s41467-018-03624-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-018-03624-y
This article is cited by
-
Deducing Phonon Scattering from Normal Mode Excitations
Scientific Reports (2019)
-
Perfect optical nonreciprocity in a double-cavity optomechanical system
Frontiers of Physics (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.