## Silicon Double-Disk Optomechanical Resonators from Wafer-Scale Double-Layered Silicon-on-Insulator

AMY NAVARATHNA, 1,2 BENJAMIN J. CAREY, 1,2 JAMES S. BENNETT, 3 SOROUSH KHADEMI, 1,2 AND WARWICK P. BOWEN 1,2,4*

1 School of Mathematics and Physics, The University of Queensland, St Lucia, Queensland 4072, Australia.

2 ARC Centre of Excellence for Engineered Quantum Systems, St Lucia, Queensland 4072, Australia. 3 Centre for Quantum Dynamics, Griffith University, Nathan, Queensland 4222, Australia. 4 ARC Centre of Excellence in Quantum Biotechnology, St. Lucia, Queensland 4072, Australia

* w.bowen@uq.edu.au

Abstract: Whispering Gallery Mode (WGM) optomechanical resonators are a promising technology for the simultaneous control and measurement of optical and mechanical degrees of freedom at the nanoscale. They offer potential for use across a wide range of applications such as sensors and quantum transducers. Double-disk WGM resonators, which host strongly interacting mechanical and optical modes co-localized around their circumference, are particularly attractive due to their high optomechanical coupling. Large-scale integrated fabrication of silicon double-disk WGM resonators has not previously been demonstrated. In this work we present a process for the fabrication of double-layer silicon-on-insulator wafers, which we then use to fabricate functional optomechanical double silicon disk resonators with on-chip optical coupling. The integrated devices present an experimentally observed optical quality factors of the order of 10 5 and a single-photon optomechanical coupling of approximately 15 kHz.

## 1. Introduction

Optomechanical WGM resonators have applications in many fields including precision sensors [1-4], biosensors [1], memory [5], communications [6] and for quantum information devices [7-9]. Such resonators support optical modes confined to their periphery by the refractive index contrast at their outer boundary. They also support mechanical degrees of freedom that strongly overlap spatially with the confined optical field [10, 11], leading to significant optomechanical coupling. WGM optomechanical resonators have been experimentally demonstrated in structures such as microspheres [12,13], microtoroids [14-18], (single) microdisks [19,20], and doubledisks [6, 21, 22].

A double-disk WGM resonator is made by fabricating two coaxial disks separated by a narrow air gap. Evanescent coupling between the two disks results in hybridization of the WGMs, forming modes that have a high field intensity in the gap [21]. This is illustrated in Fig. 1(a). Modulation of the air gap by out-of-plane mechanical motion results in large optomechanical coupling; particularly for the symmetrical out-of-plane oscillation (as shown in Fig. 1(b)). This is compounded by the tendency for out-of-plane deformations to be much larger than in-plane deformations (of equal energy), due to bending typically requiring less energy than stretching or compressing a material [23]. As such, double-disk systems achieve some of the largest optomechanical coupling rates in WGM optomechanical devices demonstrated to date [22].

The optomechanical coupling can also be increased by improving the confinement of the optical mode and ensuring strong co-localization of the optical and mechanical modes [7,24-28]. Both objectives can be achieved by moving to smaller resonators fabricated from materials with higher refractive indices. Double-disk optomechanical resonators have previously been fabricated from materials such as silica (SiO2, 𝑛 ≈ 1 44) [6, 17, 21, 22], silicon nitride (Si3N4, . 𝑛 ≈ 1 98) [29, 30], and lithium niobate (LiNO3, . 𝑛 ≈ 2 21) [31]. . There is much interest in developing silicon (Si) optomechanical devices for many applications [32-34]. Silicon's larger refractive index ( 𝑛 = 3 48) permits tighter confinement of optical modes and hence higher . coupling rates and further miniaturization of devices. This potential is further bolstered by silicon's compatibility with electronics fabrication processes and operation with light in the telecommunications C-band, leading to opportunities for integration with a range of technologies and systems. In particular, Si optomechanical resonators present ready compatibility with silicon photonic, semiconducting, and superconducting circuits. Further, recent studies have demonstrated control of defect-induced single photon emissions in silicon substrates [35, 36] offering the potential for integration with quantum information systems utilizing such photon sources.

Previous reports have demonstrated chip-scale bonding to produce high-purity Si doublelayers with the necessary oxide spacing layer ( ∼ 100 nm) for production of free-standing Si double-disks [32] or oxide-clad ring resonators [37]. However, due to the complicated nature of their production, demonstration of large-scale Si double-disk WGM resonator fabrication and experimental investigations of their optomechanical performance remained elusive.

In this work we present a truly wafer -scale approach for double layer Silicon-on-Insulator (SOI) substrate fabrication, with subsequent chip-scale fabrication and characterization of free-standing double-silicon-disk WGM optomechanical resonators. We achieve this by bonding two entire 4-inch SOI wafers prior to device fabrication. This allows for both chip- and wafer-scale double-layer SOI photonics, which are compatible with the aforementioned technologies. Further, we present Finite Element Method (FEM) modeling and comparison to experimental results of the fabricated devices. We observe optical and mechanical resonances of the disks under vacuum, and characterize the single-photon optomechanical coupling rate ( 𝑔 0) through measurements of the optomechanical spring effect.

## 2. Device Design and Fabrication

As described above, double-disk resonators consist of two coaxial Si microdisks separated vertically by a silica (SiO2) sandwich layer of thickness 𝑧 𝑔 , acting as a spacer. This SiO2 layer is selectively etched around the periphery, resulting in an air gap between the disks. When the gap is sufficiently small, evanescent coupling results in the emergence of a hybridized optical mode that exists in both disks and the gap simultaneously; this is depicted in Fig. 1(a). As expected, our optical FEM simulation (COMSOL Multiphysics) shows that the fundamental (transverse) WGMs are confined to the outer ∼ 2 μ mof the disks, which is commensurate with the region of greatest deformation caused by mechanical modes. The fundamental symmetric flapping mode is shown in Fig. 1(b).

To characterize the impact of the motion on the optical WGM, we first consider the effective refractive index ( 𝑛 eff ), which can be determined from the weighted integral method [38]:

$$n _ { \text{eff} } = \left ( \frac { \int _ { V } E ( \vec { r } ) n ^ { 2 } ( \vec { r } ) \, \mathrm d V } { \int _ { V } E ( \vec { r } ) \, \mathrm d V } \right ) ^ { \frac { 1 } { 2 } } \,.$$

Here, 𝐸 and 𝑛 are the electric field amplitude and material refractive index, respectively, at position fi 𝑟 within the simulation volume 𝑉 . The effective index of the WGMs with varying disk separation is presented in Fig. 1(c). When there is no air-gap ( i.e. 𝑧 𝑔 = 0) the WGM behaves as if contained within a single 400-nm disk, with the field well-confined inside the Si; thus, 𝑛 eff approaches 𝑛 Si . Conversely, when the gap is very large ( e.g. 𝑧 𝑔 ≥ 300 nm) the mode resembles two independent WGMs with negligible shared field (see Fig. 1(c) insets). The effective refractive index is lower than that of the 0-nm gap because a larger fraction of the optical field extends evanescently beyond the Si.

Theangular resonance frequencies of an optical WGM ( 𝜔 𝑚 ) with a radius of 𝑟 are approximately given by 𝜔 𝑚 = 𝑚𝑐 𝑟𝑛 / eff , where 𝑚 is the mode number associated with the number of wavelengths within the cavity and 𝑐 is the speed of light in vacuum. The dependence of 𝜔 𝑚 on 𝑛 eff , which is modulated by mechanical motion, ensures that changes in the disk separation will result in a shift in optical resonance frequency. For small deformations, the dependence can be summarised by the optomechanical coupling (we suppress the subscript for the 𝑚 th mode),

$$G _ { \text{OM} } = \frac { \delta \omega _ { m } } { \delta z _ { g } }.$$

A high 𝐺 OM value is advantageous for applications such as sensing and transduction [39, 40]. Our FEM model is presented in Fig. 1(d). Here, it can be seen that the optomechanical coupling diminishes with increasing separation. The maximum coupling rate occurs as 𝑧 𝑔 approaches zero; however, this regime is experimentally impracticable due to large van der Waal's and Casimir forces that can promote collapse of the disks, effectively destroying the mechanical modes of interest. For these reasons, and to ensure fabrication consistency, a separation of 60 nm was chosen. With a disk separation of 60 nm we predict an optomechanical coupling of 𝐺 OM / 2 𝜋 ≈ 150 GHz/nm, which is among the highest reported in optomechanical cavities [39,41].

Fig. 1. Axisymmetric eigenmode FEM simulation results of the double-disk resonators. (a) The optical field distribution of a hybridized (transverse electric) whispering gallery mode.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00219v1-figure-67edc2de-608c-4423-8727-c79e1852a689.jpg)

- (b) Mechanical displacement profile of the symmetric mechanical resonance (inset: asymmetric mode-profile).
- (c) Effective refractive index of the optical mode ( cf . Eq. 1) with differing disk separation (insets: mode profiles with 0 nm &amp; 300 nm spacings).
- (d) Optomechanical coupling rate of the disks for varying separations (the vertical line corresponds to the 60 nm used in our fabricated devices).

Asecond key figure of merit is the single-photon optomechanical coupling rate, which is given by:

$$g _ { 0 } = x _ { z p f } G _ { \text{OM} },$$

where 𝑥 zpf = √︁ ℏ / 2 𝑚 eff Ω is the zero point motion of the oscillation with resonance frequency of Ω / 2 𝜋 , ℏ is the reduced Planck's constant, and

$$m _ { \text{eff} } = \int _ { V } \rho ( \vec { r } ) \frac { \left | \vec { X } ( \vec { r } ) \right | ^ { 2 } } { \left | \vec { X } _ { \max } \right | ^ { 2 } } \, \mathrm d V$$

is the effective mass of the mechanical oscillation. Here 𝜌 is the material's density, fi 𝑋 is the displacement, and fi 𝑋 max is the maximum displacement amplitude of the excitation. The single-photon coupling rate is equal to the frequency shift per zero-point displacement of the oscillator, and also determines the optical force per photon [11].

Using the data extracted from both optical and mechanical FEM simulations (see S.1 for the mechanical FEM), a 120 μ m radius Si double-disk resonator with an air gap of 60 nm and undercut of 3.35 μ myields a single photon coupling rate of 𝑔 0 / 2 𝜋 = 116 kHz.

## 2.1. Further design considerations

The design of our double-disk devices broadly follows the concepts of those presented in [32]. However, unlike previous reports-which demonstrated only chip-scale double-layer SOI-we present a wafer -scale fabrication process utilizing two whole wafers of commercially available SOI (see Sec. 3 &amp; S.2.1). Furthermore, unlike many previous works that have coupled light directly from tapered optical fibers into the WGM cavity [17,21,22], our work utilizes coupling via integrated on-chip waveguides, enabling superior flexibility and scalability.

Our wafer-scale double-layer SOI approach allows for rapid production of many devices, allowing for a range of devices of varying radii and waveguide-to-disk coupling distance to be fabricated and explored on the same chip. The radii of our devices range from 15 μ m to 120 μ m. A minimum radius of 15 μ mis used to ensure radiative optical losses in the devices are negligible [42]. This range of radii presents free spectral ranges (FSRs) in the range of 1-7 nm, allowing for expedition of experimental probing.

The coupling distance ranges from 100-400 nm. The on-chip waveguides are 500 nm wide, which ensures single mode operation around 1550 nm wavelength. The waveguides' two right-angled bends form an S-shaped configuration, minimizing the potential for stray light ( e.g. , light reflected off the substrate) to couple into the output. The waveguides are tapered down to 100 nm over the course of 50 μ m, allowing for mode-matched coupling to and from tapered optical fiber tips in a similar manner to that presented in [43]. Thin tethers along the length of the waveguides minimize buckling due to stress in the Si films. These tethers are narrower than the waveguide to suppress their scattering cross-section.

## 3. Fabrication

## 3.1. Double-layer silicon-on-insulator wafer fabrication

Wafer-scale fabrication of double-disk devices requires double-layer SOI wafers. However, such wafers are not commercially available and so were produced in-house in a similar-though scaled-up-manner to that which is presented in [44]. Here SOI wafers (4 inch) consisting of a 215-nm upper layer of crystalline Si, a 3μ mburied oxide (SiO2) layer, and a 750μ mSi carrier layer were purchased from WaferPro. The wafers were thermally oxidized in a dry oxidation furnace (HiTech benchtop oxidation furnace) at 900 ◦ C for two hours, forming an SiO2 layer at the surface of the wafers. This oxide layer forms one half of the spacer layer in the final double-disk devices, as depicted in Fig. 2(b &amp; c). Thin-film reflectometry thickness analysis (SCI FilmTek 2000m) revealed an oxide thickness of 30 nm and underlying silicon layer thickness reduced to 200 nm. This increase in overall thickness is consistent with prior work [45], and

Fig. 2. Wafer scale production of the Si double-disks from two standard silicon-oninsulator (SOI) wafers (not to scale).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00219v1-figure-d5e63337-de21-44b9-b6af-f9c19a96bf64.jpg)

- (a) Initial SOI wafer stacks.
- (b) Oxidation of the top layers prior to wafer bonding, followed by backside etching from one side to remove one carrier and buried oxide layer (detail in Fig. S2).
- (c) The resultant double-layer SOI wafer.
- (d) The wafers are then diced into 15 mm × 15 mm chips and patterned with electron beam lithography, etched and undercut to create the devices (detail in Fig. S4).
- (e) A released double-disk resonator.
- Blue, Si; magenta, silica; green, electron-beam resist.

occurs because of silica's lower density than Si. The oxidized SOI wafers were then inspected for cleanliness with darkfield microscopy; pristine surfaces are a requirement for successful wafer bonding. After determining the surfaces were free of contaminants, the wafers were surface activated by exposure to oxygen plasma (200 sccm, 50 W, Oxford Instruments PlasmaPro80). The bonding process was then performed using an EVG520 IS Hot Embosser. The process consists of two pistons being pushed together. Mounted onto the top piston is a piece of quartz and the bottom piston is flat steel. When the process starts the temperature is slowly ramped up to 600 ◦ C and a force of 3.1 kN is applied. The high temperature reduces the buildup of stress between the wafers, and increases the bonding strength [46]. Bond strength was then assessed by a razor blade test, as described by [47].

Following the bonding, one side of the stack is exposed to hydroflouric acid (HF) Vapor Phase Etching (VPE) in order to remove any SiO2 grown on the back of the wafer during the thermal oxidation process (as the underside of the wafers were also exposed in the oxidation furnace). Afterwards, Deep Reactive Ion Etching (DRIE) (Plasma-Therm Versaline DRIE) is used to remove the 750μ msilicon handle layer on one side of the bonded stack. Here a two-step Bosch process is used, consisting of C4F8 gas, which forms a passivating layer, followed by an isotropic Si etch with SF6. This process was repeated and the progress of the etch was monitored periodically by removing the wafer stack and visually inspecting. Photographs of the etch progression throughout the DRIE process can be found in the Supplementary Material (Sec. S.2, Fig. S3), showing non-homogeneous removal of material. As a result, a final liquid HF etch, used to remove the 3μ mburied oxide (BOX) layer that is now partially exposed, leaves approximately 60% of the wafer area usable for fabrication of devices. Finally, the double-layer SOI substrate was diced into chips of 20 mm × 20 mm for the subsequent fabrication of optomechanical silicon double-disks.

Fig. 3. (a) Optical micrograph of a fabricated device (a wide-view photograph is shown in Fig. S5).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00219v1-figure-b4332477-de75-47a4-a1d4-bad64e5aa9bf.jpg)

Scanning Electron Microscope (SEM) images of a (b) double-disk resonator; (c) suspended waveguide; and (d) edge of the released area.

## 3.2. Device fabrication

The double-disks are fabricated using standard Electron-Beam Lithography (EBL) and ReactiveIon Etching (RIE). A schematic of the fabrication process is depicted in Fig. S4.

Electron-resist AR-P 6200.13 (CSAR, AllResist) was spincoated onto the double-layer SOI chips (1500 RPM), followed by 180 ◦ C bake producing a ∼ 580-nm coating (confirmed with thin-film reflectometry). As the total thickness of material that needs to be etched is relatively thick (2 × 200 nm Si and 60 nm SiO2), and the selectivity of the RIE is quite low (approaching 1:1) a thicker resist than typical SOI fabrication was found to be desirable to protect the top device layer. To compensate for the thicker resist and heterogeneous stack, electron scattering simulations and resultant Proximity Effect Correction (PEC) were performed with GenISys Tracer and GenISys Beamer software respectively. The chips were then patterned with a Raith EBPG-5150 electron beam pattern generator (100 kV accelerating voltage, beam current 1-100 nA) for a total dosage of approximately 300 μ C/cm 2 , followed by developing with AR600-546 developer, O-xylene and isopropyl alcohol (IPA) for 90, 5, and 15 seconds respectively.

After patterning, the chips were then etched with RIE in a three-stage process to remove the Si, SiO2 and Si device layers sequentially. Here, the Si etches consist of SF6/CHF3 20/35 sccm (82 W RF power) for 450 s, and SiO2 etch of CHF3/Ar 25/45 sccm (200 W) for 120 s. Afterwards, any residual resist was removed by soaking the chips in Remover-PG at 60 ◦ C until the resist was visibly removed upon inspection with an optical microscope. The final step in the fabrication process is release of the design with HF VPE (at 35 ◦ C).

The devices were inspected with optical microscopy, due to concerns surrounding degradation of their optical Q-factors when subjected to scanning electron microscopy (SEM). Any devices imaged with SEM were not used for subsequent measurements. SEM pictures of a fabricated double-disk device and waveguide are included in Fig. 3(b-d). The pictures are taken at an oblique angle of 45 ◦ . The high contrast between the edges of the waveguides and disks indicate the existence of an air gap.

## 4. Device Characterization

## 4.1. Optical properties of Si double-disk resonators

Optical transmission spectra were utilized to confirm the functionality of the waveguides and access the optical interaction with the double-disk resonators. For this, a fiber-coupled laser (EXFO T100S-HP model CU) was connected to a single sided taper, which was aligned with the waveguide using a 6-axis positioning stage. Another single sided taper was similarly aligned on the other end of the S-bend waveguide and connects to a photodector (PD) (New focus 1811-FC-AC,

Newport Inc.). The PD is connected to an oscilloscope (Tektronix DPO 3034 Digital Phosphor Oscilloscope), a computer (through a National Instruments Data Acquisition (DAQ) board) and a spectrum analyser (Liquid instruments, Moku:Pro), as depicted in Supplementary Fig. S6. The emission wavelength of the laser was swept whilst monitoring the transmitted intensity, as is presented in Fig. 4. Here, a high density of sharp dips are observed from the spectrum. These correspond to a Free Spectral Range (FSR) of ∼ 2 8 nm. The FSR ( . Δ 𝜆 ) can be approximated as

$$\Delta \lambda \approx \frac { \lambda ^ { 2 } } { 2 \pi r n _ { \text{eff} } },$$

which for an 𝑟 = 45μ m, 𝑧 𝑔 = 60-nm double-disk results in an FSR of 2.88 nm, in good agreement with experiment.

The fundamental WGMs of the double-disks are accompanied by higher order resonances, as observed from the additional families of sharp dips with similar FSRs. The broader features of the transmission spectrum are attributed to the polarization-dependent coupling efficiencies between the fiber tapers and on-chip waveguides [48]. These can be compensated for by manipulating the Polarization Controller (PC). However, it is evident from Fig. 4 (inset) that these polarization effects are negligible over the line-width of individual resonance dips, so we will neglect them in the following characterization.

Fig. 4. Optical transmission spectrum in ambient conditions. The sharp dips are optical resonances of the double-disks. The inset shows the transmission of a narrow span (0.2 THz).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00219v1-figure-a4924040-7986-447f-93a8-8b893956e731.jpg)

Characterization of the optomechanical interaction was performed using a slow wavelength sweep with the devices held in vacuum (see section 4.2 for details). An exemplar transmission trace is shown in Fig. 5. The asymmetrical lineshape is unexpected, as most WGMs exhibit symmetrical Lorentzian profiles [8]. This asymmetry is independent of the PC, optical power, and sweep rate, which suggests that it is not caused by thermal broadening [26]. Rather, the asymmetry appears to be due to interference between the high-Q optical mode of the resonator and low-Q optical resonances supported by the waveguide, i.e. , a Fano resonance. The Fano-type lineshape is given by [49]:

$$T ( \omega ) = T _ { 0 } + ( T _ { \infty } - T _ { 0 } ) \frac { ( q - \epsilon ) ^ { 2 } } { 1 + \epsilon ^ { 2 } },$$

where 𝑇 0 and 𝑇 ∞ are the transmissions at resonance and far off-resonance respectively, 𝑞 is the Fano asymmetry parameter, and 𝜖 = 2 Δ / 𝜅 is the detuning Δ = 𝜔 0 -𝜔 𝐿 normalised by the cavity (amplitude) loss rate, 𝜅 / 2. Here, 𝜔 0 is the cavity resonance frequency and 𝜔 𝐿 is the laser frequency. This lineshape was found to be in excellent agreement with the measured data (Fig. 5(b)). The extracted energy decay rate of the cavity is 𝜅 = 2 𝜋 × 1 71 GHz, which is in . good agreement to the estimate based upon the the Full Width at Half Maximum (FWHM) of the transmission dip. This yields an optical Q-factor of (Q = 𝜔 0 / 𝜅 ) ∼ 2 × 10 5 , which is comparable to reports of other double-disk photonic resonators (which range from 2 5 . × 10 4 [32] to ∼ 10 6 [21] for Si and SiN devices respectively). The fitted Fano parameter is 𝑞 = + 0 1744. .

## 4.2. Optomechanical properties of Si double-disk resonators

Initial investigations of the mechanical resonances of the disks were performed by detuning the laser around the Half Width at Half Maximum (HWHM) point (sideband detuning) of one optical mode and monitoring the Power Spectrum Density (PSD) of the transmitted light with a spectrum analyzer. The mechanical displacement is transduced into intensity fluctuations of the detected signal. The resulting power spectrum shows, as expected, that squeeze-film damping greatly degrades the mechanical Q-factor under ambient pressures, resulting in mechanical Q-factors of less than 10 (see S.3.1, &amp; Fig. S7). This is consistent with prior works [50, 51]. Despite this, we also see the existence of a mechanical resonance at 24 MHz, as predicted by our FEM simulation for double-disks with an undercut of ∼ 3.35 μ m(Fig. S1).

To overcome the effects of squeeze-film damping and study the intrinsic optomechanical properties of the disks, further measurements were performed under vacuum ( ∼ 10 -6 mbar). The vacuum apparatus is depicted in Fig. 5. Here, a free-space laser (New Focus 6300 High-Power Velocity Tunable Diode Laser) is coupled into a single-mode fiber. A 99:1 beamsplitter is used to monitor the input power during the experiment prior to insertion into the chamber, and the output fiber is connected to a PD (New Focus 1811-FC-AC, Newport Inc.). An oscilloscope (Tektronix, DPO3034) is used for data acquisition.

Single-sided tapered fibers are inserted into the vacuum chamber using feedthroughs. The fiber tips are glued on aluminum supports and mounted on positioning stages (SmarAct Linear Stages). Two windows (top and side view) allow for monitoring of the coupling of the fiber tips to the waveguide. Fig. 5 shows the vacuum chamber experimental set-up with a mounted chip. The chamber is pumped down with a roughing-pump-backed turbomolecular pump.

Fig. 6(a) depicts the PSD of the transmitted signal when sideband-detuned (with detuning Δ = -0 744 GHz). . Here, what was a very broad and nondescript feature in the spectrum at ambient pressure is revealed to be four distinct mechanical modes. Fitting Lorentzian lineshapes to each peak shows that the mechanical Q-factors of the four modes are on the order of ∼ 10 3 , which is comparable to previous reports on silica double-disks in vacuum conditions [22]; this is, however, the first measurement of silicon double-disk resonators under vacuum. The complete fit model and parameters are summarised in S.3.2 and Table S.1.

Determining the optomechanical coupling rate is less straightforward, but can be achieved by exploiting the optical spring effect, wherein a detuned laser drive causes a shift of the mechanical resonance frequency [8]. The expected shift of angular frequency 𝛿 Ω 𝑖 is given by

$$\frac { \delta \Omega _ { i } } { \Omega _ { i } } = \left [ 1 - \frac { C _ { i } \left ( 0 \right ) } { Q _ { i } } \cdot \frac { 2 \epsilon } { \left ( 1 + \epsilon ^ { 2 } \right ) ^ { 2 } } \right ] ^ { 1 / 2 } - 1,$$

where 𝑄 𝑖 is the zero-detuning mechanical Q-factor , and 1 𝜅 is the optical linewidth (as discussed in section 4.1). We have introduced the (resonant) optomechanical cooperativity, 𝐶 𝑖 ( 0 , which ) is related to the single-photon optomechanical coupling rate by

$$g _ { 0 } = \sqrt { \frac { C _ { i } \left ( 0 \right ) \kappa \Gamma _ { i } } { 4 n \left ( 0 \right ) } },$$

1 We use the optomechanically shifted Q-factor in the calculation of 𝐶 , which introduces an error on the order of 0.5%. This is far less than the uncertainty in the fit, and therefore negligible. See S.3.2 for details.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00219v1-figure-dc4af069-2433-4445-8dec-f34292cd27e3.jpg)

Fig. 5. (a) Schematic of vacuum chamber measurement set-up for transmission measurements of silicon double-disks. A laser is coupled into a fiber. The light is guided through a polarization controller and split by a 99:1 beamsplitter; 1% is directly connected to a power monitor, whilst the remainder is fed into the chamber and coupled to the on-chip waveguides using a single-sided tapered fiber. Transmitted light is collected by a second taper and guided to the detection photodetector (PD). The photocurrent is monitored on an oscilloscope and post-processed on a connected computer. A waveform generator is used to control the laser sweep, with the control voltage monitored at the oscilloscope to allow tracking of the output wavelength.

(b) Transmission spectrum of double-disks under vacuum. The data (red dots) is fitted with a Fano lineshape, Eq. (5) (black line). The resonance frequency is approximately 193.040 THz (1554.082 nm) and the Fano parameter is 𝑞 = + 0 1744.. This resonance . is used for the investigations presented in section 4.2. Note that the detuning is defined so that positive detuning is to the red side of the resonance.

(c) Pictures of measurement set-up for measurement of double-disks. The chip is situated on a stage, and fibers are glued onto aluminum supports mounted on linear positioners. When the vacuum chamber is closed, coupling to the disks is monitored through a top and a side window.

where 𝑛 ( Δ ) is the number of intracavity photons when the cavity is pumped at detuning Δ . Thus, we can recover 𝑔 0 by fitting multiple mechanical spectra recorded at different detunings.

Fig. 6. (a) Sideband-detuned PSD of the double-disk motion (blue) and fit according to Eq. S.1 (green), with vertical bars marking the center frequencies of each resonance. -.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00219v1-figure-3d23be5f-651b-415c-8465-4b6f30a09209.jpg)

(b) The PSDs as a function of optical detuning. Here the dotted white lines correspond Ω 𝑖

(c) Optical spring effect (Eq. 6) frequency shifts extracted from (b). The color of each curve corresponds to those allocated in (a). Experimental points are shown as circles

Data were acquired with a detuning of 0 744 GHz. to the centre of the fitted Lorentzian ( i.e. ). and fits are shown as dashed/dotted lines.

Detuning-dependent power spectra are shown in Fig. 6(b). The corresponding optical spring effect fittings (Eq. 6) are shown in Fig. 6(c), with numerical values quoted in Table 1. Two of these modes (1 and 4) achieve cooperativities greater than unity, which is an important threshold for onset of instability in optomechanical amplification [8].

To convert 𝐶 𝑖 ( 0 ) to 𝑔 0 we require the resonant intracavity photon number ( 𝑛 ( 0 ) ≃ 5 × 10 4 ), which can be calculated from the approximate waveguide photon number and the coupling strength into the cavity (see S.3.3 &amp; S.3.4 for details). We find that the optomechanical coupling rates lie between 𝑔 0 / 2 𝜋 = 7 78 kHz and . 𝑔 0 / 2 𝜋 = 15 47 kHz for the modes presented in this . study (see Table 1). These are consistent with reports for hybridized silica WGM resonators such as those presented in [22, 52]. Note that we were unable to observe optical resonances in our low-radius devices (likely due to fabrication imperfections), which were expected to have the largest optomechanical coupling rates. Our devices can therefore be further optimized to improve these 𝑔 0 values beyond those obtained in silica. The maximum experimentally determined value of 𝑔 0 is ≈ 1 8 of the expected valued based on the FEM (Sec. (2)), which is reasonably congruent. / The disparity likely resulting from fabrication inconsistencies or imperfections.

Table 1. Cooperativities and coupling rates of the investigated mechanical modes extracted from the curve fitting presented in Fig. 6. All of the uncertainties presented here derive from the fitting of experimental data to Eq. (6).

|   Mode number | 𝐶 ( 0 )            | 𝐶 0 × 10 + 5       | 𝑔 0 / 2 𝜋 (kHz)     |
|---------------|--------------------|--------------------|---------------------|
|             1 | 4 . 11 (± 0 . 15 ) | 9 . 05 (± 0 . 32 ) | 15 . 47 (± 0 . 12 ) |
|             2 | 0 . 82 (± 0 . 06 ) | 1 . 80 (± 0 . 13 ) | 7 . 86 (± 0 . 08 )  |
|             3 | 0 . 87 (± 0 . 17 ) | 1 . 92 (± 0 . 37 ) | 7 . 78 (± 0 . 13 )  |
|             4 | 2 . 31 (± 0 . 07 ) | 5 . 08 (± 0 . 15 ) | 12 . 18 (± 0 . 08 ) |

## 5. Conclusion

We have demonstrated the first wafer-scale fabrication of double-layer silicon-on-insulator wafers. These custom wafers were then utilized in the fabrication of freestanding Si double-disk WGM optomechanical resonators with associated on-chip silicon waveguides, enabling integration with an enormous variety of existing Si-based technologies.

Further, we presented the first observation of the optomechanical spring-effect in silicon double-disks. Four mechanical modes were identified that show a strong change in mechanical resonance frequency as a function of detuning. The optomechanical cooperativity for each mode is extracted from fitting to the experimental data. Using this, we determine the single-photon coupling rates in these devices.

Future devices could be improved by increasing the undercut between the disks and reducing the disks' diameters. Increasing the undercut makes the mechanical modes less sensitive to fabrication imperfections, and also tends to improve the mechanical quality factor. Reducing the disk radius is also predicted to lead to higher optomechanical coupling rates.

Further research with these devices could involve the integration of these devices with a wide range of applications including single-photon defect emitters, opto-electronics or microwave resonators which could prove useful for the fields of quantum computing and quantum state transfer.

Funding. This work was financially supported by the Australian Research Council (ARC) Centres of excellence for Engineered Quantum systems (EQUS, Grant No. CE170100009) and Quantum Biotechnology (QUBIC, Grant No. CE230100021). This project was in part funded under the Next Generation Technologies Fund, which is now administered through the Advanced Strategic Capabilities Accelerator.

Acknowledgments. The authors acknowledge the facilities, and the scientific and technical assistance, of the Australian Microscopy &amp; Microanalysis Research Facility at the Centre for Microscopy and Microanalysis, The University of Queensland. This work was performed in part at the Queensland node of the Australian National Fabrication Facility, a company established under the National Collaborative Research Infrastructure Strategy to provide nano- and microfabrication facilities for Australia's researchers. This work was financially supported by the Australian Research Council (ARC) Centres of excellence for Engineered Quantum systems (EQUS, Grant No. CE170100009) and Quantum Biotechnology (QUBIC, Grant No. CE230100021). This project was in part funded under the Next Generation Technologies Fund, which is now administered through the Advanced Strategic Capabilities Accelerator.

Disclosures. The authors declare no conflicts of interest.

Supplemental document. See Supplemental Material for supporting content, including further FEM modeling, greater detail about the double-layer SOI wafer production process and the device fabrication, results on optomechanics at ambient conditions, and investigations into the cavity photon number and photon coupling rates.

## References

- 1. S. Frustaci and F. Vollmer, 'Whispering-gallery mode (wgm) sensors: review of established and wgm-based techniques to study protein conformational dynamics,' Curr. Opin. Chem. Biol. 51 , 66-73 (2019).
- 2. X. Liu, W. Liu, Z. Ren, Y. Ma, B. Dong, G. Zhou, and C. Lee, 'Progress of optomechanical micro/nano sensors: a review,' Int. J. Optomechatronics 15 , 120-159 (2021).
- 3. F. Gotardo, B. J. Carey, H. Greenall, G. I. Harris, E. Romero, D. Bulla, E. M. Bridge, J. S. Bennett, S. Foster, and W. P. Bowen, 'Waveguide-integrated chip-scale optomechanical magnetometer,' Opt. Express 31 , 37663-37672 (2023).
- 4. S. Basiri-Esfahani, A. Armin, S. Forstner, and W. P. Bowen, 'Precision ultrasound sensing on a chip,' Nat. Commun. 10 , 132 (2019).
- 5. M. B. Kristensen, N. Kralj, E. Langman, and A. Schliesser, 'A long-lived and efficient optomechanical memory for light,' (2023).
- 6. Y. Lei, Z.-G. Hu, M. Wang, Y.-M. Gao, Z. Zuo, X. Xu, and B.-B. Li, 'Fully reconfigurable optomechanical add-drop filters,' Appl. Phys. Lett. 121 , 181110 (2022).
- 7. J. Chan, M. Eichenfield, R. Camacho, and O. Painter, 'Optical and mechanical design of a 'zipper' photonic crystal optomechanical cavity,' Opt. Express 17 , 3802 (2009).
- 8. W. P. Bowen and G. J. Milburn, Quantum optomechanics (CRC Press, Boca Raton London New York, 2015).
- 9. J.-D. Tang, Q.-Z. Cai, Z.-D. Cheng, N. Xu, G.-Y. Peng, P.-Q. Chen, D.-G. Wang, Z.-W. Xia, Y. Wang, H.-Z. Song, Q. Zhou, and G.-W. Deng, 'A perspective on quantum entanglement in optomechanical systems,' Phys. Lett. A 429 , 127966 (2022).
- 10. S. Mancini, D. Vitali, and P. Tombesi, 'Optomechanical Cooling of a Macroscopic Oscillator by Homodyne Feedback,' Phys. Rev. Lett. 80 , 688-691 (1998).
- 11. M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, 'Cavity optomechanics,' Rev. Mod. Phys. 86 , 1391-1452 (2014).
- 12. R. Ma, A. Schliesser, P. Del'Haye, A. Dabirian, G. Anetsberger, and T. J. Kippenberg, 'Radiation-pressure-driven vibrational modes in ultrahigh-Q silica microspheres,' Opt. Lett. 32 , 2200 (2007).

13. Z. Shen, Y.-L. Zhang, Y. Chen, F.-W. Sun, X.-B. Zou, G.-C. Guo, C.-L. Zou, and C.-H. Dong, 'Reconfigurable optomechanical circulator and directional amplifier,' Nat. Commun.

9

, 1797 (2018).

- 14. K. H. Lee, T. G. McRae, G. I. Harris, J. Knittel, and W. P. Bowen, 'Cooling and Control of a Cavity Optoelectromechanical System,' Phys. Rev. Lett. 104 , 123604 (2010).
- 15. E. Verhagen, S. Deléglise, S. Weis, A. Schliesser, and T. J. Kippenberg, 'Quantum-coherent coupling of a mechanical oscillator to an optical cavity mode,' Nature 482 , 63-67 (2012).
- 16. F. Ruesink, J. P. Mathew, M.-A. Miri, A. Alù, and E. Verhagen, 'Optical circulation in a multimode optomechanical resonator,' Nat. Commun. 9 , 1798 (2018).
- 17. C. Bekker, C. G. Baker, R. Kalra, H.-H. Cheng, B.-B. Li, V. Prakash, and W. P. Bowen, 'Free spectral range electrical tuning of a high quality on-chip microcavity,' Opt. Express 26 , 33649 (2018).
- 18. C. Bekker, R. Kalra, C. Baker, and W. P. Bowen, 'Injection locking of an electro-optomechanical device,' Optica 4 , 1196-1204 (2017).
- 19. X. Lu, J. Y. Lee, and Q. Lin, 'High-frequency and high-quality silicon carbide optomechanical microresonators,' Sci. Reports 5 , 17005 (2015).
- 20. W. C. Jiang and Q. Lin, 'Chip-scale cavity optomechanics in lithium niobate,' Sci. Reports 6 , 36920 (2016).
- 21. X. Jiang, Q. Lin, J. Rosenberg, K. Vahala, and O. Painter, 'High-Q double-disk microcavities for cavity optomechanics,' Opt. Express 17 , 20911 (2009).
- 22. C. Meng, G. A. Brawley, S. Khademi, E. M. Bridge, J. S. Bennett, and W. P. Bowen, 'Measurement-based preparation of multimode mechanical states,' Sci. Adv. 8 , 7585 (2022).
- 23. J. Ward and O. Benson, 'Wgm microresonators: sensing, lasing and fundamental optics with microspheres,' (2011).
- 24. J. Chan, T. P. M. Alegre, A. H. Safavi-Naeini, J. T. Hill, A. Krause, S. Gröblacher, M. Aspelmeyer, and O. Painter, 'Laser cooling of a nanomechanical oscillator into its quantum ground state,' Nature 478 , 89-92 (2011).
- 25. E. Gavartin, R. Braive, I. Sagnes, O. Arcizet, A. Beveratos, T. J. Kippenberg, and I. Robert-Philip, 'Optomechanical Coupling in a Two-Dimensional Photonic Crystal Defect Cavity,' Phys. Rev. Lett. 106 , 203902 (2011).
- 26. R. Leĳssen, G. R. La Gala, L. Freisem, J. T. Muhonen, and E. Verhagen, 'Nonlinear cavity optomechanics with nanomechanical thermal fluctuations,' Nat. Commun. 8 , ncomms16024 (2017).
- 27. H. Ren, M. H. Matheny, G. S. MacCabe, J. Luo, H. Pfeifer, M. Mirhosseini, and O. Painter, 'Two-dimensional optomechanical crystal cavity with high quantum cooperativity,' Nat. Commun. 11 , 3373 (2020).
- 28. X. Lu, J. Y. Lee, and Q. Lin, 'Silicon carbide zipper photonic crystal optomechanical cavities,' Appl. Phys. Lett. 116 , 221104 (2020).
- 29. S. Lee, S. C. Eom, J. S. Chang, C. Huh, G. Y. Sung, and J. H. Shin, 'A silicon nitride microdisk resonator with a 40-nm-thin horizontal air slot,' Opt. express 18 , 11209-11215 (2010).
- 30. G. S. Wiederhecker, S. Manipatruni, S. Lee, and M. Lipson, 'Broadband tuning of optomechanical cavities,' Opt. Express 19 , 2782-2790 (2011).
- 31. Y. Zheng, Z. Fang, S. Liu, Y. Cheng, and X. Chen, 'High-q exterior whispering-gallery modes in a double-layer crystalline microdisk resonator,' Phys. Rev. Lett. 122 , 253902 (2019).
- 32. H. Moradinejad, A. H. Atabaki, A. H. Hosseinnia, A. A. Eftekhar, and A. Adibi, 'Double-Layer Crystalline Silicon on Insulator Material Platform for Integrated Photonic Applications,' IEEE Photonics J. 6 , 1-8 (2014).
- 33. Y. A. V. Espinel, F. G. S. Santos, G. O. Luiz, T. P. M. Alegre, and G. S. Wiederhecker, 'Brillouin Optomechanics in Coupled Silicon Microcavities,' Sci. Reports 7 , 43423 (2017).
- 34. R. Dehghannasiri, H. Moradinejad, T. Fan, A. H. Hosseinnia, A. A. Eftekhar, and A. Adibi, 'Integrated Optomechanical Resonators in Double-Layer Crystalline Silicon Platforms,' in 2018 IEEE Photonics Conference (IPC), (IEEE, Reston, VA, 2018), pp. 1-2.
- 35. M. Hollenbach, Y. Berencén, U. Kentsch, M. Helm, and G. V. Astakhov, 'Engineering telecom single-photon emitters in silicon for scalable quantum photonics,' Opt. Express 28 , 26111-26121 (2020).
- 36. A. Durand, Y. Baron, W. Redjem, T. Herzig, A. Benali, S. Pezzagna, J. Meĳer, A. Y. Kuznetsov, J.-M. Gérard, I. Robert-Philip, M. Abbarchi, V. Jacques, G. Cassabois, and A. Dréau, 'Broad diversity of near-infrared single-photon emitters in silicon,' Phys. Rev. Lett. 126 , 083602 (2021).
- 37. R. M. Briggs, M. Shearn, A. Scherer, and H. A. Atwater, 'Wafer-bonded single-crystal silicon slot waveguides and ring resonators,' Appl. Phys. Lett. 94 , 021106 (2009).
- 38. S. Ohke, Y. Satomura, T. Umeda, and Y. Cho, 'A new integral expression for the effective refractive index of optical waveguides and its application for a fast numerical solution finder for the effective refractive index,' Opt. communications 118 , 227-234 (1995).
- 39. J. Xia, Q. Qiao, G. Zhou, F. S. Chau, and G. Zhou, 'Opto-mechanical photonic crystal cavities for sensing application,' Appl. Sci. 10 (2020).
- 40. N. Lauk, N. Sinclair, S. Barzanjeh, J. P. Covey, M. Saffman, M. Spiropulu, and C. Simon, 'Perspectives on quantum transduction,' Quantum Sci. Technol. 5 , 020501 (2020).
- 41. M. Winger, T. D. Blasius, T. P. M. Alegre, A. H. Safavi-Naeini, S. Meenehan, J. Cohen, S. Stobbe, and O. Painter, 'A chip-scale integrated cavity-electro-optomechanics platform,' Opt. Express 19 , 24905-24921 (2011).
- 42. J. R. Buck, 'Cavity QED in microsphere and Fabry-Perot cavities,' Ph.D. thesis, California Institute of Technology (2003). Medium: PDF Version Number: Final.
- 43. T. G. Tiecke, K. P. Nayak, J. D. Thompson, T. Peyronel, N. P. de Leon, V. Vuletić, and M. D. Lukin, 'Efficient fiber-optical interface for nanophotonic devices,' Optica 2 , 70-75 (2015).
- 44. H. Moradinejad, 'Hybrid multi-layer CMOS-compatible material platform for high-performance integrated nanophotonics,' Ph.D. thesis, Georgia Institute of Technology (2017).
- 45. J. X. Zhang and K. Hoshino, 'Fundamentals of nano/microfabrication and scale effect,' in Molecular Sensors and Nanodevices, (Elsevier, 2019), pp. 43-111.
- 46. H. Moriceau, F. Rieutord, F. Fournel, Y. Le Tiec, L. Di Cioccio, C. Morales, A. M. Charvet, and C. Deguet, 'Overview of recent direct wafer bonding advances and applications,' Adv. Nat. Sci. Nanosci. Nanotechnol. 1 , 043004 (2011).
- 47. K. Pantzas, F. Fournel, A. Talneau, G. Patriarche, and E. Le Bourhis, 'Measuring the surface bonding energy: A comparison between the classical double-cantilever beam experiment and its nanoscale analog,' AIP Adv. 10 , 045006 (2020).
- 48. J. M. Ward, A. Maimaiti, V. H. Le, and S. N. Chormaic, 'Contributed Review: Optical micro- and nanofiber pulling rig,' Rev. Sci. Instruments 85 (2014). 111501.
- 49. M. F. Limonov, M. V. Rybin, A. N. Poddubny, and Y. S. Kivshar, 'Fano resonances in photonics,' Nat. Photonics 11 , 543-554 (2017).
- 50. G. S. Wiederhecker, L. Chen, A. Gondarenko, and M. Lipson, 'Controlling photonic structures using optical forces,' Nature 462 , 633-636 (2009).
- 51. Q. Lin, J. Rosenberg, X. Jiang, K. J. Vahala, and O. Painter, 'Mechanical oscillation and cooling actuated by the optical gradient force,' Phys. review letters 103 , 103601 (2009).
- 52. J. Rosenberg, Q. Lin, and O. Painter, 'Static and dynamic wavelength routing via the gradient optical force,' Nat. Photonics 3 , 478-483 (2009).
Amy Navarathna Benjamin J. Carey James S. Bennett Soroush Khademi Warwick P. Bowen
Silicon Double-Disk Optomechanical Resonators from Wafer-Scale Double-Layered Silicon-on-Insulator
Whispering gallery mode (WGM) optomechanical resonators are a promising
technology for the simultaneous control and measurement of optical and
mechanical degrees of freedom at the nanoscale. They offer potential for use
across a wide range of applications such as sensors and quantum transducers.
Double-disk WGM resonators, which host strongly interacting mechanical and
optical modes co-localized around their circumference, are particularly
attractive due to their high optomechanical coupling. Large-scale integrated
fabrication of silicon double-disk WGM resonators has not previously been
demonstrated. In this work we present a process for the fabrication of
double-layer silicon-on-insulator wafers, which we then use to fabricate
functional optomechanical double silicon disk resonators with on-chip optical
coupling. The integrated devices present an experimentally observed optical
quality factors of the order of 10^5 and a single-photon optomechanical
coupling of approximately 15 kHz.