Sympathetic Cooling of Levitated Optomechanics through Nonreciprocal Coupling

TL;DR

The paper addresses the cooling limits in levitated optomechanics imposed by cavity dissipation and environmental noise. It introduces a non-Hermitian cooling scheme based on nonreciprocal coupling between two nanoparticles, where only one particle is directly coupled to a cavity while the other is cooled sympathetically via unidirectional energy transfer. analytically derives an effective bath for the cavity-coupled particle and a reduced two-mode non-Hermitian model, providing steady-state solutions that show deeper cooling of the target particle as nonreciprocity increases, validated by full master-equation simulations. The work demonstrates a new energy-flow engineering approach with potential for scalable quantum control and ultra-low-noise sensing in levitated systems.

Abstract

Optomechanical cooling of levitated nanoparticles has become an essential topic in modern quantum physics, providing a platform for exploring macroscopic quantum phenomena and high-precision sensing. However, conventional cavity-assisted cooling is fundamentally constrained by cavity dissipation and environmental noise, limiting the attainable minimum temperature. In this work, we propose a non-Hermitian optomechanical cooling scheme through nonreciprocal coupling between two levitated nanoparticles, where one particle is directly cooled by an optical cavity and the other is cooled indirectly through a non-Hermitian interaction. Both analytical solutions and numerical simulations reveal that increasing nonreciprocity enhances directional energy transfer, enabling the target particle to reach a lower phonon occupation than is achievable in conventional cavity cooling. This study demonstrates a new cooling mechanism driven by non-Hermitian interactions, offering theoretical guidance for realizing controllable energy flow and deep cooling in levitated optomechanical systems, and paving the way for future developments in quantum control and sensing technologies.

Figures (4)

• Figure 1: (a) Schematic diagram of the nonreciprocal optomechanical cooling system. particle $A$ is optically trapped inside the cavity and coupled to the cavity mode through radiation pressure, while particle $B$ is levitated outside the cavity. Each particle supports a quantized harmonic mode with frequencies $\omega_a$ and $\omega_b$, respectively. (b) Illustration of the nonreciprocal coupling between particles $A$ and $B$ and their mechanical dissipation processes. The coupling strengths and mechanical decay rates are explicitly labeled.
• Figure 2: Steady--state phonon numbers of particles $A$ and $B$ versus the nonreciprocal coupling parameter $K_g$. All parameters are expressed in dimensionless units with $\omega_a=\omega_b=\Delta=1$. The remaining parameters are fixed as $\gamma = 1$, $\overline{n}_a = \overline{n}_b = 20$, $\kappa_a = \kappa_b = 10^{-4}$, and $g_2 = 0.1$. The reciprocal coupling coefficient is chosen as $g_1 = g_2 (K_g + 1)/(K_g - 1)$ so that $K_g>1$ corresponds to an increasingly nonreciprocal interaction. Particle $A$ and $B$ correspond to the red and blue lines, respectively. (a) Results for cavity-$A$ coupling $g=0.05$. Here $n_A$ remains almost constant, while $n_B$ decreases significantly with $K_g$ due to efficient directional energy transfer into the strongly cooled particle $A$. (b) Results for weaker coupling $g=0.01$. In this regime $n_B$ continues to decrease exponentially with $K_g$, whereas $n_A$ increases slowly and reaches values one order of magnitude larger than in panel (a), reflecting phonon accumulation at particle $A$. In both panels, $n_A$ and $n_B$ approach each other as $K_g\to1^{+}$, consistent with the reciprocal coupling limit.
• Figure 3: Time evolution of the phonon number $n$ for particles $A$ and $B$. particle $B$ is the target to be cooled, while particle $A$ acts as the auxiliary system mediating the cooling process. The initial phonon numbers of particles $A$ and $B$ are $0$ and $2$, respectively. The parameters are $g = 0.1$, $g_2 = 0.1$, $\gamma = 1.0$, $\kappa = \kappa_a = \kappa_b = 0.01$, and $\overline{n}_a = \overline{n}_b = 1$. Particle $A$, $B$, and the cavity correspond to the red solid, blue solid, and green dashed lines, respectively. (a) $g_1 = 0.11$ ($K_g = 1.1$). The steady--state phonon numbers are $n_B = 0.32$ and $n_A = 0.38$. (b) $g_1 = 0.40$ ($K_g = 1.4$). The steady--state phonon numbers are $n_B = 0.22$ and $n_A = 1.69$. As $K_g$ increases, the nonreciprocity strengthens, leading to lower $n_B$ and higher $n_A$, which highlights the directional energy transfer and enhanced non-Hermitian cooling efficiency.
• Figure 4: Comparison between the analytical and numerical steady--state phonon numbers of particle $B$ as a function of the nonreciprocal coupling parameter $K_g$. The red dots represent the results obtained from the sparse-Liouvillian steady--state solver, and the blue solid line shows the analytical solution based on Eq. \ref{['eq:adjoint']}. $g_2 = 0.1$, $g_1 = g_2(K_g + 1)/(K_g - 1)$, $g = 0.05$, $\kappa_a = \kappa_b = 10^{-4}$, $\gamma = 2$, and $\overline{n}_a = \overline{n}_b = 20$.