Collective dynamics of strain-coupled nanomechanical pillar resonators.

Semiconductur nano- and micropillars represent a promising platform for hybrid nanodevices. Their ability to couple to a broad variety of nanomechanical, acoustic, charge, spin, excitonic, polaritonic, or electromagnetic excitations is utilized in fields as diverse as force sensing or optoelectronics. In order to fully exploit the potential of these versatile systems e.g. for metamaterials, synchronization or topologically protected devices an intrinsic coupling mechanism between individual pillars needs to be established. This can be accomplished by taking advantage of the strain field induced by the flexural modes of the pillars. Here, we demonstrate strain-induced, strong coupling between two adjacent nanomechanical pillar resonators. Both mode hybridization and the formation of an avoided level crossing in the response of the nanopillar pair are experimentally observed. The described coupling mechanism is readily scalable, enabling hybrid nanomechanical resonator networks for the investigation of a broad range of collective dynamical phenomena.

Introduction

Nanomechanical pillar resonators[1-5] represent particularly versatile types of nanomechanical resonators. They allow for the convenient integration of multiple functionalities using semiconductor heterostructures[3,6,7], and for the dense integration into large arrays[8,9]. In recent years, the dynamics of both bottom-up and top-down fabricated nanomechanical pillar resonators has been explored[2,4,5,10]. Each pillar exhibits two orthogonally polarized fundamental flexural modes of similar eigenfrequency, which are typically not completely degenerate, even for the case of a circular cross section, as a result of fabrication imperfections[4]. Most experimental investigations focus on a single mechanical mode[2,10] or on the interplay of a single mechanical mode with another degree of freedom, e.g. an embedded quantum dot[3,7], an NV center spin[8], a surface acoustic wave[11], or an optical cavity mode[12]. Some scanning probe applications target the mechanical coupling of the two fundamental flexural modes of a single nanopillar induced by an external force field for vectorial force field sensing[4,5]. To date the intrinsic coupling of the flexural modes, or the coupling between the flexural modes of adjacent nanopillars has not been addressed.

At the same time, the quest of engineering a controlled intrinsic coupling between nanomechanical resonators is receiving an increasing amount of attention. At present, intrinsic coupling of adjacent beam or string resonators has been reported and relies on the strain distribution in a shared clamping point[13-15]. In a related approach, the membranes constituting the building blocks of nanomechanical phonon waveguides exchange energy through physical connections[16,17]. Note that the appearance of collective dynamical effects in a resonator array is enabled by a sufficiently large coupling strength, mitigating the unavoidable, and fabricational disorder-induced detuning of its elements[18].

Here, we translate the concept of strain-induced intrinsic coupling to vertically oriented nanomechanical pillar resonators sharing the same substrate[3-6]. We consider a pair of inverted conical nanopillar resonators like the one displayed in Fig. 1. The pillars are etched into a (100) GaAs substrate using reactive ion etching[2] and feature eigenfrequencies in the range of a few megahertz. In the following, we employ the strong coupling condition as a well-defined and experimentally easy-to-assess smoking gun demonstration of large inter-pillar coupling. Strong coupling, i.e. a coupling rate exceeding the linewidth of the resonances, is not required for most applications targeting the collective dynamics of nanopillar resonators. However, for realistic pillar parameters, fulfilling this condition is indicative of a sufficiently large coupling rate to beat the disorder in a nanopillar array.

Fig. 1

Nanomechanical system. Scanning electron micrograph of a pair of nanopillar resonators with a height of , a foot radius of  and a center-to-center distance of in a tilted view. The taper angle of gives rise to a somewhat larger head radius. Scale bar corresponds to . Inset shows a zoom of the clamping point of the two nanopillars

Results

Mode hybridization

We first investigate a nanopillar pair with bottom radius , height , and center-to-center distance (see inset of Fig. 1). The taper angle of is the same for all nanopillars discussed in this work. The nanopillar pair is driven at frequency using a piezo actuator in a scanning electron microscope, allowing to image the resulting envelope of its mechanical vibration. In total, we find four well separated vibrational modes with no spectral overlap as shown in Supplementary Fig. 1, which are identified with the four eigenmodes of the pillar pair. With respect to the indistinguishable crystal directions, these modes are in the following referred to as ‘horizontal’ (H) and ‘vertical’ (V) polarization of the left (L) and right (R) pillar, and correspondingly labeled LH, LV, RH, and RV (see Supplementary Fig. 1). One of these eigenmodes (LV) is shown in Fig. 2a, b, which display the vibrational envelope of the nanopillar pair driven at , imaged from the top and in a tilted view. A more careful inspection of the mode reveals that in Fig. 2a, b not only the left resonator vibrates with a large amplitude but also the right resonator exhibits a simultaneous vibration, albeit with a much smaller amplitude (see also Supplementary Movies 1 and 2). This is also apparent from Fig. 2c which displays the vibrational amplitude of both pillars extracted from micrographs obtained at different drive frequencies in the vicinity of the resonance . Clearly, the two amplitudes evolve simultaneously with the drive frequency , which implies the existence of a hybridized mode. The underlying coupling between the two pillars is mediated by the joint substrate.

Fig. 2

Hybridized mode. a Scanning electron micrograph of a resonantly driven () pair of nanopillars (, , ) imaged from the top and b, in a tilted view. Red and blue dotted lines indicate the circumference of the undriven left and right resonator, respectively. c Amplitude for different drive frequencies of left (red circles) and right (blue squares) resonator determined from the scanning electron micrographs. Scale bar in a and b corresponds to

Avoided level crossing

To obtain a more thorough understanding of the observed inter-pillar coupling, we measure the response of a second pillar pair with , , and using the optical detection setup described in the Methods section. Thermal tuning, readily implemented by the laser used for optical detection, is employed to sweep the eigenfrequency of the higher-frequency pillar through that of the other pillar which remains largely unaffected (see Supplementary Note 3 for details). Figure 3 shows an avoided crossing of a pair of nanopillar resonators, which is indicative of strong mechanical coupling. A fit of the data using the model of two linearly coupled harmonic oscillators is also included as a black solid line[19]. It yields a coupling strength which exceeds the linewidth of the mechanical resonances . Hence, we demonstrate strong intrinsic coupling between the two nanopillar resonators, for this specific set of geometry parameters. The two modes to the left of the avoided crossing are assigned to the vertical vibration of each resonator (LV, RV) via scanning electron micrographs shown in the insets of Fig. 3, where the tuned pillar with the higher frequency corresponds to the right pillar. Further evidence for the inter-pillar coupling arises from the evolution of the vibration amplitudes in Fig. 3. Since the laser is focused on the frequency-tuned right pillar, mainly the vibration of this one resonator is detected and the vibration of the left pillar is only weakly resolved through the stray field of the laser. Clearly, the transition of the strong signal of the right pillar from the upper to the lower branch of the avoided crossing reflects the hybridization of the two pillars near their resonance at and .

Fig. 3

Strong coupling. Frequency response measurements of the right nanopillar under thermal tuning reveals an avoided crossing as an evidence for strong coupling. The black line shows a fit to the data, yielding a level splitting of , with the linewidth . Insets show respective scanning electron micrographs of the two modes, indicating a vertical oscillation direction when the resonators are resonantly driven near (lower inset) and (upper inset), far from the avoided crossing

The observed strong coupling supports the conclusions drawn from Fig. 2, the data for which was acquired without frequency tuning the right pillar and with variable . An avoided crossing measured for this pair with is included in Supplementary Fig. 1. The untuned situation probed in the scanning electron microscope measurements corresponds to the state of the system towards the leftmost edge of the avoided crossing: far from resonance, the individual eigenmodes of the two pillars dominate the response, however, a slight hybridization is already apparent as a consequence of their strong coupling.

Geometry dependence of the coupling strength

Strain-mediated coupling through the substrate is expected to depend on the geometry of the pillar pair. In particular, the bottom radius of the nanopillars, their height and their center-to-center distance have proven influential. In the following, we investigate the dependence of the nanopillar coupling strength on these parameters in measurements and finite element simulations, using the evaluation procedure detailed above.

Figure 4 shows the measured coupling strength for pillar pairs of different foot radius and height, plotted as a function of their separation. Clearly, the coupling strength increases with decreasing center-to-center distance, while the largest coupling rates are found for nanopillars with a larger foot radius. In order to validate these findings, finite element simulations are performed. An example for the case of identical pillars with realistic clamping conditions is shown in Fig. 5a, b. A significant overlap of the stress distribution is observed, indicating the intrinsic, strain-mediated coupling between the pillars which is also reflected in the appearance of an antisymmetric hybrid mode.

Fig. 4

Measured geometry dependence of coupling strength. Experimentally determined coupling rate of the vertical modes of the two nanopillars over their center-to-center distance for several samples and pillar pair geometries with & (red), & (blue), & (green). Error margins indicate the tolerance of the fitted coupling rate. The inset shows the vibration direction for differently oriented pillar pairs (see Supplementary Note 5 for details)

Fig. 5

Finite element simulation of nanopillar pairs. a Simulation of the stress distribution between two (identical) nanopillar resonators vibrating in the antisymmetric hybrid mode. Red (blue) corresponds to high (low) stress, respectively. b Finite element simulation model, highlighting the realistic clamping area of the pillar pair as well as the narrowing mesh implemented near the pillar feet. The inset shows a zoom of the transition from a nanopillar to the substrate. c–e Finite element simulation of coupling strength of vertical (filled circles) and horizontal (empty circles) modes as a function of c center-to-center distance (, ), d nanopillar bottom radius (, ) and e nanopillar height (, ). The taper angle for all simulations is . We assume an isotropic substrate with Young’s Modulus

Frequency tuning is incorporated by a variation of Young’s Modulus of one of the nanopillars, mimicking thermal tuning. The level splitting is then also obtained by fitting the model of the coupled oscillators to the simulation data (see Supplementary Note 4). We investigate the dependence of on each of the parameters , , and while the other two parameters remain fixed. Only when sweeping the bottom radius , is adapted accordingly to ensure a constant edge-to-edge distance.

Figures 5c–e display the simulated coupling strength as a function of the three parameters under investigation. Figure 5c depicts an increasing coupling when the center-to-center distance of the two nanopillars is decreased. This is a typical signature of strain coupling[20], and in agreement with the measurements in Fig. 4. Furthermore, Fig. 5d predicts an increasing coupling strength with the bottom radius of the two nanopillars, as a consequence of the strain caused at the pillar foot upon its deflection, and the experimental data in Fig. 4 clearly reflects this behavior. Finally, Fig. 5e displays a decrease of the coupling strength with the height of the pillars, again in consequence of the resulting strain profile. The dependence on the height, however, can not clearly be established in Fig. 4 since the differences in coupling for the two investigated pillar heights are smaller than the error bars.

The results shown Fig. 4 have been obtained for pillar pairs with different orientations on the substrate as indicated in the inset of Fig. 4 (see Supplementary Note 5 for details). It is expected that the angular dependence of Young’s modulus of the crystalline GaAs substrate should lead to an angular dependence of the strength of the strain-induced coupling (see Supplementary Note 5). However, differences in coupling strength with pillar orientation are not experimentally resolved, likely because of the fabrication-induced disorder in the pillar geometry governing the uncertainty in our measurements. To validate the influence of the substrate’s Young’s Modulus, further studies should address the vibration polarization of hybridized (symmetric or antisymmetric) modes for arbitrary orientation of the pillar pair.

Discussion

In conclusion, we establish a new nanomechanical platform for the investigation and direct visualization of collective behavior of coupled nanopillars. We reveal strain-induced coupling between two adjacent nanomechanical pillar resonators with coupling rates of up to  kHz, and show mode hybridization as well as the formation of an avoided level crossing under thermal tuning of one of the nanopillars. The coupling strength is found to depend on the center-to-center distance of the nanopillars, as well as their diameter and height. Numerical finite element simulations reproduce the observed scaling, confirming that the coupling between two nanopillar resonators is mediated by strain in the substrate which acts as a shared clamping point.

The established inter-pillar coupling can be readily employed for hybrid nanopillars e.g. based on III–V semiconductor heterostructures[3,6,7,21,22] featuring quantum dots, quantum wells, or integrated cavities. This may open the path towards the direct, mechanically mediated synchronization of single-photon sources, e.g. for applications in quantum information processing[23]. In addition, the scheme is scalable to a large number of pillars in a straightforward way, since vertically oriented nanopillars are ideally suited for dense integration, while a large coupling rate mitigates their unavoidable disorder-induced detuning. Assuming a typical quality factor of and a realistic disorder of ~%, this condition is fulfilled for a coupling rate exceeding the linewidth by about an order of magnitude, a ratio which has been accomplished for several of the investigated pillar pairs. In addition we envisage to implement in-situ thermal tuning of the pillar eigenfrequencies in an array using a spatial light modulator[24]. The resulting nanomechanical network offers inherent, acoustically mediated nearest-neighbor coupling[25]. The large vibration amplitudes of singly-clamped nanopillars in the range of tens or hundreds of nanometers allow for the microscopic imaging of their response and thus for the direct visualization of the many-body dynamics in an all-mechanical array. This promises important insights for the emergent field of collective dynamical phenomena which includes phenomena such as acoustic metamaterials[17,26,27], synchronization[28-31], topologically protected transport[27,32,33], or non-reciprocal signal transduction[34], and may pave the way towards nanomechanical computing[35] or nanomechanical implementations of neural networks [36].

Methods

Fabrication

The circular nanopillar heads are defined by electron-beam lithography with acceleration voltage and aperture on a GaAs substrate with spin-coated polymethyl methacrylate resist. The pattern transfer is performed via electron-beam evaporation of nickel and subsequent lift-off. A final anisotropic inductively coupled plasma reactive ion etching (ICP-RIE) process reveals the inverted conical nanopillars. Anisotropic etching is achieved with a gas flow of 13-sccm silicon tetrachloride and 0.2-sccm nitrogen, a RF power of and an ICP power of .

Measurement details

Measurements are performed in vacuum and at room temperature using piezo-actuation, while the response of the nanopillars is read-out by scanning electron microscope imaging[37] or by optical detection[3], focusing a laser with wavelength onto the head of one nanopillar and detecting the vibration-induced modulation of its reflection.

Supplementary Movie 1

Supplementary Movie 2