Figure 1.

Schematic of OMS in which the optical cavity contains a photon mode, whose creation (annihilation) operator is a†(a), coupling to a phonon mode, whose creation (annihilation) operator is b†(b). The two cavity modes decay into the reservoirs respectively."

Figure 2.

The numerical results of the normalized correlation function ${g}_{\beta }^{(2)}(x)={G}_{\beta }^{(2)}(x)/\max \{{G}_{\beta }^{(2)}(x^{\prime} )\}$. The parameters are ωb = 1, ω = 10−7 and (a) g = 10−5, γa = 0.13, γb = 10−5; (b) g = 10−5, γa = 0.13, γb = 0.1; (c) g = 10−5, γa = 1, γb = 10−5; (d) g = 0.1, γa = 0.13, γb = 10−5."

Figure 3.

The relative differences between the results in figure 2 and those based on the scattering theory. The parameters are the same as those in figure 2."

Figure 4.

Representative transitions among the energy levels of the system corresponding to the three scattering processes, where the subfigures (a)–(c) correspond to processes (a)–(c) respectively. We have omitted the higher-order transitions of g due to the weak coupling condition."

Figure 5.

The relative proportion of ${\int }_{0}^{+\infty }{\rm{d}}{{xG}}_{1\sim 3}(x)$ to their sum. The proportion of process (a)–(c) are given by ${\int }_{0}^{+\infty }{\rm{d}}{{xG}}_{1\sim 3}(x)/{\sum }_{n=1}^{3}{\int }_{0}^{+\infty }{\rm{d}}{{xG}}_{n}(x)$, respectively. The parameters are ωb = 1, g = 10−5, γa = 0.13, γb = 10−5."

Figure B1.

Some other numerical results of the normalized correlation function ${g}_{\beta }^{(2)}(x)={G}_{\beta }^{(2)}(x)/\max \{{G}_{\beta }^{(2)}(x^{\prime} )\}$ with non-experimental parameters. The parameters are ωb = 1, ω = 10−7 and (a) g = 0.01, γa = 0.02, γb = 0.002; (b) g = 0.3, γa = 20, γb = 5; (c) g = 0.05, γa = 1, γb = 0.05; (d) g = 0.3, γa = 2, γb = 5."