Molecular cavity optomechanics: a theory of plasmonenhanced Raman scattering
Abstract
The exceptional enhancement of Raman scattering crosssection by localized plasmonic resonances in the nearfield of metallic surfaces, nanoparticles or tips has enabled spectroscopic fingerprinting of single molecules and is widely used in material, chemical and biomedical analysis. The conventional explanation attributes the enhancement to the antenna effect focusing the electromagnetic field into subwavelength volumes. Here we introduce a new model that additionally accounts for the dynamical and coherent nature of the plasmonmolecule interaction and thereby reveals an enhancement mechanism not contemplated before: dynamical backaction amplification of molecular vibrations. We first map the problem onto the canonical model of cavity optomechanics, in which the molecular vibration and the plasmon are parametrically coupled. The optomechanical coupling rate, from which we derive the Raman cross section, is computed from the molecules Raman activities and the plasmonic field distribution. When the plasmon decay rate is comparable or smaller than the vibrational frequency and the excitation laser is bluedetuned from the plasmon onto the vibrational sideband, the resulting delayed feedback force can lead to efficient parametric amplification of molecular vibrations. The optomechanical theory provides a quantitative framework for the calculation of enhanced crosssections, recovers known results, and enables the design of novel systems that leverage dynamical backaction to achieve additional, modeselective enhancement. It yields a new understanding of plasmonenhanced Raman scattering and opens a route to molecular quantum optomechanics.
Introduction
In 1973 Fleischmann et al. first reported the dramatic enhancement of the Raman scattering crosssection of molecules on rough metal surfaces (1), an effect confirmed in 1977 by Van Duyne et al. (2). Two decades later, this technique known as surfaceenhanced Raman scattering (SERS) enabled the detection of single molecules (3); (4). Enhancement factors in the range of  have been reported to occur at “hot spots” (5), regions of high electromagnetic fields associated with localized plasmonic resonances. Using the plasmon at the tip of a scanning tunneling microscope (6) has led to a powerful analytical tool for sensitive Raman imaging: TERS (tipenhanced Raman scattering). SERS and TERS both rely on the phenomenon of plasmonenhanced Raman scattering, which is today widely employed in the fields of material and surface science (6), nanotechnology (7), chemistry (3) and even invivo biomedical applications (8).
The generally accepted model for SERS invokes the combined enhancement by the plasmonic hot spot of the incoming electromagnetic field and the Raman scattered field (5); (9). It predicts an enhancement of the Raman crosssection proportional to the fourth power of the field enhancement. Although this “ law” has been verified experimentally, the observation of even larger enhancements and of anomalous Stokes/antiStokes intensity ratio (10) have raised the suspicion that a “vibrational pumping” mechanism was involved (11). Moreover in recent experiments the maximal enhancement was achieved when the laser was bluedetuned from the plasmon resonance by the vibrational frequency (exciting the antiStokes vibrational sideband) (58); (47). Under these conditions large nonlinear effects were also evidenced in (58); (59). These observations are calling for a new theoretical understanding and for further investigations (15).
Here we show for the first time that SERS scenarios can be mapped onto the canonical model of cavity optomechanics (16) (Fig. 1A), in which a dynamic and coherent interaction takes place between two parametrically coupled and nonresonant harmonic oscillators, namely the molecular vibration and the plasmonic cavity. The optomechanical coupling rate can be computed from the Raman activity and the plasmonic field distribution, from which we derive the Raman cross section and recover conventional results.
The novel enhancement mechanism revealed by our approach is dynamical backaction amplification (17) of the vibrational mode due to the nonvanishing plasmon lifetime. In optomechanical microstructures, this effect was first evidenced by the amplification of mechanical breathing modes in silica microtoroids (18) under bluedetuned excitation, leading to a range of new phenomena as reviewed in (42). We find that SERS systems (i) can feature the suitable dissipation and frequency hierarchies (despite short plasmon decay times, i. e. low quality factors), and (ii) exhibit exceptionally large optomechanical vacuum coupling rates (20), so that under suitable conditions dynamical backaction amplification can occur via Raman scattering and lead to rich new physics, such as large nonlinearities and outofequilibrium vibrational occupancies.
Model
The plasmon mode is formally equivalent to an optical cavity and is modeled by a harmonic oscillator of frequency with bosonic creation and annihilation operators and the Hamiltonian (Fig. 1A). Each normal vibrational mode of a molecule is described by an effective mass , a frequency and normal coordinates with the corresponding position operator , where is the zeropoint motion and , are bosonic creation and annihilation operators. The vibrational mode thermal occupancy at temperature is given by ( is the Boltzmann constant) which for high frequency Ramanactive modes (1000 cm corresponding to THz) is well below 1 at room temperature (even under moderate heating of the plasmonic particles). The few lowest energy levels of the molecular vibration are well approximated by a harmonic oscillator () ( Fig. 1C; possible effects of anharmonicity are discussed later).
Since the frequency hierarchy is satisfied , and assuming the molecule has no optically allowed electronic transitions resonant with the plasmon, the coupling between vibrational and plasmonic modes is purely parametric (as in the case of optomechanical systems (21) and in contrast to resonant coupling (22); (23)), and the vibrational displacement leads to a dispersive shift in the plasmon resonance frequency (Fig. 1B) according to
(1) 
Denoting by the polarizability of the molecule, which is dependent on the vibrational mode displacement , the optomechanical coupling rate is given by (cf. Supplementary Material)
(2) 
with the mode volume of the plasmonic cavity (24) and the permittivity of vacuum. The corresponding vacuum optomechanical coupling rate describes the plasmon frequency shift related to the zeropoint motion of the molecular vibration.
The Hamiltonian for the complete system is formally identical to the one obtained in cavity optomechanics, with an interaction term (42)
(3) 
that describes the coherent coupling between a mechanical oscillator (here the molecular vibration) and an electromagnetic cavity mode (here the localized plasmon). The vibration acts on the plasmon via the dispersive plasmon frequency shift of eq. (1) that can change the plasmon occupancy . In turn, the plasmon acts back on the vibration via the timedependent force
(4) 
where is the plasmon occupancy. The delay in the plasmonic cavity response to changes in the resonance frequency leads to a component of the force that is out of phase with respect to the vibrational motion (see Supplementary Material). When the driving laser frequency matches a phonon sideband () the phase shift is , corresponding to a purely viscous force. For blue detuning (Fig. 2A), this leads to an effective mechanical gain of electromagnetic origin and to parametric amplification. The resulting dynamical backaction can be described as a delayed feedback loop (Fig. 2B) (25). In the frequency domain, the plasmonic and vibrational time responses correspond to spectral filtering. The transduction from plasmon to vibration occurs through the force , with an amplification factor proportional to the intracavity plasmonic field (induced by the pump laser) and to the Raman polarizability (contained in ).
Following the conventions in Raman spectroscopy we introduce the massweighted cartesian displacement coordinates(53) . In this notation the vacuum coupling rate is
(5) 
The Raman activity in a simplified one dimensional model satisfies . Even for small molecules, and without considering possible enhancement from electronic resonances, we obtain optomechanical coupling rates in the 10 – 100 GHz range (Supplementary Material, section 4), which are 4 to 5 orders of magnitudes higher than in stateofart microfabricated optomechanical structures (27).
For a realistic value of the plasmon quality factor () (28) in the nearinfrared frequency range ( THz for nm)
Raman crosssection calculations
Assuming that the amplitude of vibration is small, the interaction between the plasmonic and the vibrational mode is linearized and the plasmonenhanced Raman scattered power is then derived by solving classical Langevin equations, which account for both unitary evolution and dissipative terms. When the dissipation hierarchy is satisfied, and neglecting backaction for now (as in conventional SERS theories), the Stokes (resp. antiStokes) power is (see Supplementary Material)
(6a)  
(6b) 
where the incoming photon flux is ( being the laser power) and the average occupation number of the mode in equilibrium with an environment at temperature is . is the fraction of plasmonic energy decaying into radiation coupled to the excitation/detection optics. The Stokes (resp. antiStokes) Raman crosssection is defined by where is the laser intensity at the excitation spot of area . We obtain the crosssections for Stokes and antiStokes scattering without backaction (as in conventional electromagnetic theory)
(7a)  
(7b) 
To see how this formalism (without backaction yet) recovers the conventional theory, we note that the field inside the plasmonic cavity satisfies the energy relation . On the other hand the incoming power flux is given by , so that we can approximate the field enhancement factor by , and because we find that . This expression recovers the accepted electromagnetic enhancement (“ law” (33); (3)) while explicitly showing the contribution of the plasmonic density of states at the laser and Stokes (resp. antiStokes) frequencies (Fig. 3A). Because the optomechanical coupling rate is derived from the Raman tensor, our model can be extended to account for all effects previously studied such as (i) the dependence of cross section on the molecule spatial orientation, (ii) the selection rules corresponding to different light polarization, and (iii) the increase in cross section close to electronic resonances, possibly altered by chemical effects and charge transfers.
A new enhancement mechanism: Dynamical Backaction
The optomechanical theory uncovers a novel enhancement mechanism: dynamical backaction between the molecular vibration and the plasmon resonance, which leads to a modification of the damping rate of the vibrational mode. Under weak coupling condition (), the change in damping rate can be expressed by (34)
(8) 
where is the plasmon occupancy. When pumping the plasmon with bluedetuned light (), the backaction damping rate is negative () leading to amplification (18) of the vibrational mode and to an outofequilibrium vibrational occupancy
(9) 
where is the cooperativity
(10) 
For the Stokes and antiStokes crosssections are enhanced according to
(11a)  
(11b) 
and they become power dependent, leading to a superlinear increase of the Raman signal with pump intensity. Also the antiStokes/Stokes ratio becomes “anomalous” under dynamical backaction (i.e. exhibits values deviating from the equilibrium Boltzmann factor)
(12) 
where the last approximation is valid for low thermal occupancy typical of highfrequency vibrations. A parametric instability occurs when the amplification rate exceeds the damping rate (), corresponding to coherent regenerative oscillations (phonon lasing) (18). By increasing the quality factor of the plasmonic structure, the threshold is lowered and occurs for a laser detuning closer to the phonon sideband (). Close to this threshold the system exhibits a highly nonlinear response, which could provide a mean to achieving superresolution in TERS (58) (see Section 4 & 5 in the Supplementary Material).
Molecular vibrations exhibit anharmonic potentials at higher occupancy (Fig. 1C) and intermode couplings that cause internal vibrational redistribution (IVR) (35), both of which are not included explicitly in our model. IVR would prevent reaching the threshold for regenerative oscillations of a single mode by introducing additional damping channels, and at the same time could lead to the appearance of a broad Raman background originating from other vibrational or rotational modes indirectly excited (58). Anharmonicity of the potential, on the other hand, could lead to frequency shifts and broadening of the Raman peak under high amplification (close to threshold). We note that although the reduced damping rate should manifest as a reduced linewidth of the vibrational mode, this signature could be masked by both IVR and anharmonicity. Observation of the optical spring effect (a static shift in the frequency of the vibration due to the average value of the radiation pressure force) would also be difficult with molecules.
Collective coupling
The innovative fabrication process and careful characterization methods described by Zhu and Crozier in (47) enable the authors to investigate plasmonic cavities with welldefined field confinement properties, reaching down to the onset of quantum tunneling between the two metallic parts (54); (37). Based on the geometric estimate of the field confinement described in (54), we can approximate the mode volume in (47) by (see Supplementary Material, Section 3.2). This volume also sets the approximate number of molecules contributing to the plasmonenhanced Raman signal. Considering monolayers of thiophenol covering the metallic cylinders, we find with the thiophenol surface density, the radius of the cylinder, its height and the size of the gap.
Under these conditions the optomechanical interaction should be coherently driving a collective oscillator mode (38).
Introducing the experimental parameters (at the onset of quantum tunneling) and assuming (39) in our calculations, we estimate that Zhu and Crozier are reaching the parametric amplification instability for an incoming power of W. This compares favorably to the value of W below which the authors operate in order to avoid sample damages and nonlinear effects in the plasmon. Consequently, we predict that their experiment might be showing signatures of dynamical backaction amplification. It is therefore worth comparing our predictions with the conventional theory of SERS.
Setting the plasmon resonance at a wavelength of nm and the incoming power slightly below the threshold, we provide a comparison between the two models (without and with dynamical backaction) and the experimental results of Zhu et al. (Fig. 3). Under the condition of negligible dynamical backaction (red dashed line) maximal Raman signal occurs when the pump wavelength is halfway between the plasmon resonance and the phonon sideband, recovering the results of the electromagnetic theory (40). The sharpening and shift of the Raman excitation spectrum can be qualititavely predicted and explained by our new model that allows dynamical backaction to take place in the plasmonic cavity.
The anomalous antiStokes/Stokes ratio contemplated in the previous section could also be explored with the experimental setup here described. The deviation from the equilibrium Boltzmann factor is demonstrated in Fig. 4 for realistic pump intensity. The signature of dynamical backaction amplification can be experimentally distinguished from thermal effects such as local heating by (i) the sharp sensitivity of the anomaly on the laser detuning from the plasmon; and (ii) its nonlinear power dependence (see Fig. S1B in the Supplementary Material).
Outlook
The theory provides a unified, physically grounded and quantitative framework for the phenomena of SERS. It challenges the common wisdom, based on the “ law”, that broader plasmonic resonances – overlapping with both laser and Stokes wavelengths – achieve the highest possible enhancement. Instead we show that vibrational amplification is more efficient for larger plasmon factors (i.e. narrower resonances) and for blue detuned excitation on the antiStokes vibrational sideband. We note that even higher amplification factors could be obtained for double plasmonic resonances spaced by the vibrational frequency, corresponding to a multimode optomechanical system. These insights are of major relevance for the design of SERS systems leveraging backaction for more sensitive spectroscopy and higher resolution imaging.
More broadly, the theory lays the foundations of molecular cavity optomechanics and opens unforeseen research directions. The rich physics of cavity optomechanics is now accessible in systems of nanometric dimensions featuring coupling rates several orders of magnitude higher than stateoftheart microfabricated devices. One example of observable phenomenon is the counterintuitive reduction of the Stokes crosssection under simultaneous pumping with a second, reddetuned laser, which would lead to damping of the vibrational mode. Other effects include optomechanicallyinduced transparency and frequency conversion. Functional materials could be designed using layers of molecules collectively coupled to twodimensional cavities. By choosing molecules exhibiting vibrational modes that are both infrared and Raman active, one could achieve coherent frequency conversion between THz radiation and visible light. Finally, because of the high frequency of molecular vibrational modes, they are not thermally excited at room temperature, fostering the prospects for molecular quantum optomechanics. With the mechanical oscillator in its quantum ground state, schemes such as the creation of nonclassical states of motion (41) become feasible with simple experimental setups at room temperature.
Acknowledgments The authors would like to thank Dr. E. Verhagen for stimulating discussions and Dr. E. Bremond for his precious help in running the chemical simulations. This work was partially supported via an ERC Advanced Grant, NCCR of Quantum Engineering, QSIT as well as the Swiss National Science Foundation. C. G. acknowledges the support of the Swiss National Science Foundation through an Ambizione Fellowship. N. P. acknowledges the support of a MarieCurie Fellowship.
– Supplementary Material –
Molecular cavity optomechanics: a theory of plasmonenhanced Raman scattering
Appendix A Cavity optomechanical principles
a.1 Symbols and definitions
Notation  Description 

Dimensionless intracavity field (normalized so that photon number)  
Flux of the pump field in units of s ( is the photon flux)  
Frequency of the excitation laser  
Resonance frequency of the plasmonic cavity  
Total energy damping rate of the plasmonic cavity  
Rate of intrinsic losses, including absorption and uncoupled radiation  
Coupling rate between the cavity field and the incoming/outgoing fields  
Total decay rate of the plasmonic cavity  
Coupling ratio ( for critical coupling)  
Energy damping rate of the vibrational mode  
Mechanical frequency of the vibrational mode  
Displacement coordinate of the vibrational mode  
Effective mass of the vibrational mode  
Zeropoint motion of the vibrational mode  
Average occupation number of the vibrational mode  
Average number of plasmons inside the cavity  
Detuning between the laser and the plasmon resonance frequency  
Coupling rate between the plasmon and the vibrational mode  
Optomechanical vacuum coupling rate to the vibrational mode 
a.2 From the formal Hamiltonian to classical equations of motion
We present in the following a simple treatment of the optomechanical system to give an intuitive picture of the important phenomena considered in the manuscript. A complete derivation including all the assumptions made can be found in reviews on cavity optomechanics (42); (43).
We start from the expression describing the interaction between a radiation mode and a vibrational mode as introduced in the main text
(S1) 
(We recall the expression of the zeropoint motion: ).
A first derivation of this optomechanical Hamiltonian was given by Law (44). The interaction term is corresponding to a nonlinear process that involves a product of three operators. The driving laser at frequency is modeled by a coherent field: . It is convenient to switch to a reference frame rotating at the laser frequency, in which the total Hamiltonian writes
(S2) 
Making use of the BakerCampbellHausdorff formula, this change of reference frame allows to eliminate the explicit time dependence of the Hamiltonian. Dissipations can be introduced by writing the quantum Langevin equations. Within the inputoutput formalism (45), defining as the detuning between the laser and the plasmon, we obtain the equations
(S3)  
(S4)  
where we introduced the input noise terms expressed in the rotating frame: the vacuum noise and driving laser’s noise entering the plasmonic cavity, and the thermal noise on the vibrational mode. Assuming a shotnoise limited laser and neglecting the thermal excitation at the plasmon frequency, the correlators associated with these fluctuations are given by:
(S5)  
(S6)  
(S7)  
(S8) 
where is the thermal occupancy of the bath at the vibrational frequency (valid for sufficiently high vibrational quality factor).
The equations of motion for the corresponding creation operators are obtained from the
relations ,
. These
equations are describing the complete evolution of the plasmonic excitations coupled
to the molecular vibrations. It is instructive
to give first the classical version of these equations.
We express the vibrational degreeoffreedom in the position and momentum operators defined by
(S9)  
(S10) 
In this representation the Langevin equations for the vibrational mode become
(S11)  
(S12) 
where and are the quantum noise operators.
After taking expectation values the noise terms average to zero and one obtains
(S13)  
(S14) 
where and .
This force of electromagnetic origin
a.3 Linearization of the cavity field
Setting the time derivatives to zero in the system of equations described above, we find the stationary solutions
(S15) 
We assume that the amplitude of the vibrations is described
by
and is small in comparison to the typical length scale of the plasmonic
cavity
(S16) 
order equation
(S17) 
The homogeneous solution can be neglected if the measurement time is much longer than the other timescales in the system (as is always the case in SERS). We can thus consider only the particular solution, equal to the steadystate solution without optomechanical coupling
(S18) 
order equation
(S19) 
We can neglect the homogeneous solution following the same reasoning as before. We look for a particular solution of the form and obtain the amplitudes of the antiStokes and the Stokes fields, respectively
(S20) 
a.4 Dynamical Backaction Force
We develop the expression of the force to first order in
(S21) 
In the last expression, the first term corresponds to a constant force applied to the vibrational mode . The second term is timedependent and can be expressed as a sum of in and outofphase components
(S22) 
with
(S23)  
(S24) 
We can now insert the expression for the electromagnetic force into the Langevin equation (S14)
(S25)  
(S26) 
The equation is thus describing a mechanical oscillator displaced by a constant force , with a shifted natural frequency and damped at an effective rate
(S27)  
(S28) 
Depending on the laser detuning from the (shifted) plasmonic resonance, three different situations can arise:

: For resonant excitation there is no radiationinduced change of the damping rate

: For a reddetuned laser, the damping rate is increased, corresponding to “cooling”

: For a bluedetuned laser, the damping rate is decreased, corresponding to amplification
Deep in the sidebandresolved regime () the maximal increase (decrease) of the damping rate is obtained when () and the expression for the maximum backaction damping rate due to optomechanical interactions can be approximated by
(S29) 
where . When the damping is increased trough the coupling with the plasmonic field (red detuning), the vibrational mode is losing power into the plasmonic cavity and thus is being optically “cooled”. On the contrary, when the damping rate is decreased (blue detuning), power is transfered from the plasmonic field to the vibrational mode, whose motion is thereby amplified: this is dynamical backaction amplification.
a.5 Collective optomechanical plasmonvibrational coupling
In this section, we consider identical phonon modes – with annihilation operators – coupled to the same plasmonic cavity. Corresponding experimental situations include (i) a layer of molecules filling the gap of a metal dimer (47), where each represents the same Raman mode for each molecule; (ii) the optical phonon of a piece of bulk or 2D material interacting with a plasmonic cavity, where each represent the same vibrational mode of each unit cell (for example the band of graphene in (48)). For simplicity, we consider identical vibrational frequencies , intrinsic damping rates and coupling rates to the cavity . We follow a derivation similar to (49), starting from the Langevin equations in the rotating frame of the molecular vibration and of the Stokes field. Because the plasmon decay rate is much larger than the vibrational damping the cavity field evolution can be adiabatically eliminated and replaced by its steady state solution. This leads to a linear system of differential equations describing the evolution of the phonon modes coupled via the cavity driven by a laser tuned on the antiStokes (blue) vibrational sideband:
(S30) 
The solutions of these equations reveal two different behaviours. On the one hand, we find linearly independent and degenerate eigenmodes () with the same eigenvalue that can be written in the general form:
(S31) 
such that . The shift of the plasmonic resonance caused by each of the collective modes is proportional to its collective optomechanical coupling rate . This shows that these modes are dark, i.e. decoupled from the plasmonic cavity and thus not affected by dynamical backaction.
On the other hand, the eigenmode with eigenvalue (the equivalent of the superradiant mode in cavity QED) is the only collective mode coupled to the cavity and its backaction damping rate is enhanced by a factor , i.e. . This scaling translates to a scaling of (see equation 8 of the main text).
a.6 Description of the Raman crosssection using the optomechanical formalism
Writing the Stokes scattered power as and expressing following eq. (S20) yields
(S32) 
Here, we assume that the vibrational amplitude is related to the thermal energy, i.e. . We describe the amplitude of molecular vibrations with the help of its quantum mechanical description
(S33) 
The Stokes crosssection is defined as with ( is the incident power and the illuminating spot area
(S34) 
where is the average occupation number of the mode in equilibrium with an environment at temperature . The same development can be done for the antiStokes process, leading to the result expressed in the main text.
Appendix B Plasmonvibration vacuum optomechanical coupling rate
The change in resonance frequency of a cavity when a dielectric is inserted in an air gap is described by (50)
(S35) 
where is the energy stored inside the plasmonic cavity (51), the plasmonic cavity field and the induced dipole per unit volume. As long as the plasmonic cavity dimension is sufficiently small in comparison to the incoming wavelength, the quasistatic approximation remains valid and the magnetic energy can be neglected. We can express the contribution of the electric field to the total energy stored in the metal by where and the factor depends on the dielectric function of the metal and the plasmon resonance frequency (52).
We assume that a single molecule is located at the position of maximum electric field. Introducing the molecular dipole moment
(S36) 
The induced dipole is related to the electric field by ; the dependence of the linear polarizability on the molecular displacement is developed to first order
(S37) 
where the gradient of polarizability has units Note that in the expression for the polarizability of the molecule contains the contributions from the internal electonic transitions so that it is frequency dependent.
We can now quantify the sensitivity of the plasmonic frequency to the molecular vibrations to first order in
(S38) 
expressed as a function of the effective mode volume of the cavity . When the plasmon resonance frequency is getting closer to the bulk plasma frequency, the energy stored in the metal grows and the mode volume defined here is getting larger (52).
We have thus shown that the plasmon frequency is coupled to the molecular vibration and that this coupling can be quantified by the optomechanical vacuum coupling rate (also called single photon coupling rate)
(S39) 
Applying the concepts used in optomechanics, we are thus able to find a quantitative way of describing the interaction between the molecular vibrations and the plasmonic resonance.
Appendix C Quantifying the optomechanical coupling
c.1 Calculation of
Following the conventions used in Raman spectroscopy we use instead of the normal coordinates the reduced coordinates with units , refered to as the massweighted cartesian displacement coordinates (53). In this notation the kinetic energy is given by . Expressing the vacuum coupling rate (S39) in the reduced coordinates () leads to
(S40) 
In optomechanical systems the effective mass is defined to verify the equipartition theorem (), with the energy stored in this mode.
In order to compute the optomechanical vacuum coupling rate, we need the value of the Raman polarizability tensor. The Raman activity can be inferred from experimental data and, using density functional theory (DFT) simulations, it is possible to have access to the tensor components. Generally, the Raman intensity, describing the orientationaveraged magnitude of the Raman scattering, is given by two invariant scalars built from the space derivative of the polarizability tensor
(S41) 
where and are the isotropic and anisotropic parts of the tensor.
In a simplified one dimensional model, the Raman activity can be described as with units . In simulations, using gaussian units, this first quantity is however expressed in . The relation between the usual description of the Raman activity and that one is
(S42) 
Consequently, the vacuum optomechanical coupling rate (S39) is given by