The coupling between acoustic vibrations in a lithium niobate bulk acoustic wave resonator and microwave photons of a re-entrant microwave cavity was investigated at a temperature close to 4 K. Coupling was achieved by placing the acoustic resonator in the location of the re-entrant cavity electric field maxima in a symmetric “split-post” configuration with a large overlap between the microwave field and the acoustic mode, allowing acoustic modulations of the microwave frequency. We show that the acoustic modes in this setup retain large inherent quality factors of greater than 106. A maximum optomechanical coupling rate was determined to be g0 = 0.014 mHz, four orders of magnitude larger than previous results obtained using a quartz BAW at 4 K in a similar experimental setup but using a single post-re-entrant cavity resonator.
Investigating the fundamental interaction between photonic and phononic systems is crucial for advancing various applications within the domain of quantum technology and quantum electrodynamics.1,2 This understanding has led to the emergence of numerous associated predicted phenomena. Resolved sideband cooling,3,4 parametric amplification,5–7 optomechanically induced transparency/absorption,6,8 and long-range entanglement9,10 are just a few of these phenomena that have had important ramifications for testing fundamental quantum limits,11–15 dark matter detection,16–18 quantum communication, information processing, and storage, along with generation and manipulation of quantum states.19–21 Importantly, for these systems, highly coherent quantum state readouts and increased coupling rates are integral for parametric detection techniques. This is particularly relevant in the search for gravitons22,23 and in the development of parametric transducers for achieving low noise readouts at microwave frequencies.3,24
Ground state preparation of truly macroscopic (gram-scale) mechanical resonators is a highly sought after experimental demonstration.12,25,26 For one, there still remain questions regarding the quantum-to-classical crossover point that will be revealed by such systems, but in addition, the precision readout of such devices may offer insights into quantum gravity and high frequency gravitational wave detection.13,27–31
Studies on these photon–phonon interactions employ various types of resonant photonic device architectures, such as superconducting circuits,32 co-planar resonators,33 3-D microwave cavity resonators,34 and whispering gallery resonators,21,35,36 combined with various types of mechanical architectures displaying high coupling rates, such as membrane-in-the-middle systems,6 trapped ion particles,37 and acoustic-mechanical modes such as bulk acoustic wave (BAW) resonators.6,25,26,35 BAW resonators are low-loss mechanical resonators that exhibit long mechanical coherence times due to their specifically engineered convex surface, which helps trap the majority of phonons in the center of the resonator.27
In this work, a macroscopic centimeter size-scale lithium niobate (LiNbO3) crystal with its axis cut in the z-direction was manufactured with a diameter of 30 mm, a center thickness of 2 mm, and a convex radius of curvature of 100 mm. The crystal exhibited coherent phonon modes trapped in the center of the crystal with frequencies of order MHz. Figure 1 shows a photo of the LiNbO3 crystal along with the two halves of a “split-post” microwave cavity resonator,38 which sandwiches the BAW in the post gap, allowing the maximum electric field to be applied to the region of maximum phonon trapping as shown in Fig. 2. This means the acoustic phonons will be parametrically upconverted to microwave frequencies. The motivation was drawn from similar experiments conducted with quartz BAW resonators.26,39 However, the cavity design differs from before as we have used a symmetric design in which the gap is located in the center of a coaxial cavity. This particular design results in a more highly confined electric field distribution focused into the center of the BAW, and we term this design the “split-post re-entrant cavity.” Moreover, the reason for the emerging popularity of LiNbO3 over quartz is that it exhibits a stronger piezoelectric coefficient in addition to a low loss tangent. For example, LiNbO3-BAW has demonstrated a high quality factor × frequency product of Hz at room temperature.40
FIG. 1.
The convex lithium niobate crystal BAW resonator, along with the two halves of the re-entrant split-post microwave cavity, fabricated from oxygen-free copper.
FIG. 1.
The convex lithium niobate crystal BAW resonator, along with the two halves of the re-entrant split-post microwave cavity, fabricated from oxygen-free copper.
Close modal
FIG. 2.
(a) Not to scale microwave-BAW cavity setup with a frequency discriminator output. The synthesizer was tuned to the microwave cavity and measured at the mixer output with the FFT to detect the upconverted BAW frequency. (b) Microwave transmission measurements (S21) with a vector network analyzer between ports 1 and 2, comparing 300 K with 4 K, showing an increase in microwave Q-factor. (c) To scale electric flux density and magnetic field intensity density plots inside the split-post microwave cavity. The electric flux density is concentrated in the center of the BAW in the vicinity of the MHz phonon modes, which are indicated in (a).
FIG. 2.
(a) Not to scale microwave-BAW cavity setup with a frequency discriminator output. The synthesizer was tuned to the microwave cavity and measured at the mixer output with the FFT to detect the upconverted BAW frequency. (b) Microwave transmission measurements (S21) with a vector network analyzer between ports 1 and 2, comparing 300 K with 4 K, showing an increase in microwave Q-factor. (c) To scale electric flux density and magnetic field intensity density plots inside the split-post microwave cavity. The electric flux density is concentrated in the center of the BAW in the vicinity of the MHz phonon modes, which are indicated in (a).
Close modal
Much smaller LiNbO3 nanostructures have been used previously in various optical applications, such as quantum sensors41,42 and optical transducers.43–48 Our study investigated LiNbO3 microwave-BAW coupling rates, with the goal of creating strongly coupled coherent microwave readouts of large macroscopic phonon modes for fundamental physics applications that require larger masses. Here upconversion opens up possibilities for harnessing effects such as parametric amplification, cooling, and squeezing, while also providing a method for reading out displacement in large non-piezoelectric crystals.
To model the photon–phonon interaction, the dynamics are characterized by the following Hamiltonian:
1 where
a(
a†) and
b(
b†) are the bosonic annihilation (creation) operators of the microwave and the acoustic resonance, at frequencies
ωc and
ωm, respectively, and
g0 is the optomechnical coupling rate.
The Hamiltonian in
(1) describes the interaction between a microwave electric field between the two posts and the acoustic phonon mode of the BAW device, which results in radiation pressure acting on the acoustic modes. This radiation pressure displaces the oscillator (by displacement
x) and leads to parametric modulation of the microwave cavity fields as shown in
(2), where the second term on the RHS represents the linear coupling, which can be implemented as a displacement readout,
1 Furthermore, we can define the following quantities:
where
G, also known as the frequency pull factor, is measured from the experimental data and helps evaluate
g0, which is the photon–phonon coupling rate, also known as the single photon coupling rate. Additionally,
is the zero-point fluctuation of the mechanical resonance, where ℏ is the reduced Planck’s constant and
Meff,m is the effective mode mass.
Figure 2(a) illustrates the experimental setup for measurements conducted at 4 K. The experiment utilizes the split post-microwave cavity, which functions as a frequency discriminator. This phase bridge detects frequency shifts of the cavity in response to mechanical modes by ensuring that the LO and RF signals are in quadrature, making the system insensitive to amplitude fluctuations. The power splitting ratio is set by a −6 dB directional coupler. The signal generator is an analog synthesizer (Keysight E8663D), covering the frequency range of 100 kHz–9 GHz. Therefore, the experiment measures the frequency modulation of the BAW resonator up-converted to microwave frequencies. The output of the mixer is observed on the FFT spectrum analyzer.
The setup in Fig. 2 also shows the two posts of the microwave cavity, one of which was anodized to provide two electrically isolated posts, allowing them to function as electrodes to piezoelectrically excite the BAW mechanical modes with an arbitrary waveform generator (AWG). Since the applied voltage to the electrodes is known, we can estimate the piezoelectric-induced crystal displacement using known LiNbO3 electrostatic material parameters, allowing us to determine the change in displacement dx. The measured change in microwave cavity frequency, dωc, is then used to determine G.
It is known that for piezoelectric materials, the mechanical displacement
x is proportional to the charge on the surface.
1 In our setup, this charge is enhanced due to the split post geometry and is given by
q =
kmx, where
km is the effective piezoelectric coupling constant defined using the Butterworth–Van Dyke (BWD) model, which can be understood from the BAW-resonator lumped circuit element theory.
49 For a mechanical mode,
m, relations between the piezoelectric constant and the corresponding LCR circuit elements of the BWD model are represented as
28 Here,
Qm is the mode quality factor,
Rm is the motional resistance, and
Lm and
Cm are the effective inductance and capacitance, respectively. Similarly, the current from the equivalent BAW circuit model may be presented in terms of mechanical resonance
ωm as
Therefore, the values of motional resistance could be measured at 4 K using an impedance analyzer.
The operational cavity mode was the fundamental split-post-mode [Fig. 2(c)], with a frequency of 6.075 GHz at 4 K. The loaded quality factor, QL, improved significantly from 300 to 4 K to QL = 330–2500, respectively, due to reduced resistive losses from the copper cavity walls as shown in Fig. 2(b). The coupling coefficients were determined from reflection measurements from ports 1 and 2 and calculated to be β1 = 0.8 and β2 = 0.136. Therefore, the unloaded microwave Q-factor can thus be determined to be Q0 = QL(1 + β1 + β2) = 4250. The increased Q-factor helped achieve better discriminator sensitivity for the phase bridge setup. The incident input power to the microwave cavity was 0.01 mW.
The change in the displacement of the mechanical mode is given by the simple relationship between the calculated displacement Δ
x and the change in the output voltage Δ
V from the mixer, as follows:
where
dV/
dωc is the discriminator sensitivity (phase-bridge); the ratio at which frequency fluctuations are converted to synchronous voltage fluctuations in the readout, while
G = −
dωc/
dx is the frequency pull factor introduced previously and is evaluated from the voltage response output from the mixer corresponding to the displacement of mechanical resonance.
From
(1), we can also calculate
in terms of
xzpf,
Driving an applied time varying MHz signal across the split-posts, piezoelectrically exciting the BAW while simultaneously driving the microwave cavity at its resonant frequency, yields a mechanically induced frequency modulation of the microwave pump frequency, producing side bands on the microwave carrier. This is triggered by the induced displacement of the LiNbO
3 BAW when tuned to the acoustic modes. Due to the trigonal symmetry of lithium niobate, two-mode families exist. “A” refers to longitudinal polarization, where bulk displacement occurs in the z-direction, while “B” denotes a shear wave polarized mode with dominant displacement in the x-y plane. The first subscript refers to the fundamental polarization inside the BAW resonator, which has been studied in Ref.
40. More details of the different suffixes in the modes can be found in Ref.
27. In this work, we study the A
3,0,0 and A
5,0,0 longitudinal modes, and B
5,0,0, B
7,0,0, and B
9,0,0 shear modes.
Figures 3(c) and 3(d) show the mixer output voltages for the longitudinal A5,0,0 mode and the shear B7,0,0 modes, respectively. Normalized reflection coefficients S11 are measured by driving the BAW through the anodized post and measuring the reflected power on an impedance analyzer, which is shown as the green trace. Figure 3 illustrates that the mixer output voltage is indeed a result of the mechanical mode being driven.
FIG. 3.
FEM simulations of the displacement field of the (a) 8.3 MHz longitudinal and (b) 5.8 MHz shear modes, with dashed lines indicating the effective mode volumes. (c) and (d) show the corresponding mixer output voltage in blue and S11 reflection measurements in green using an impedance analyzer across the split-post. Loaded Q-factors are determined by curve fits, as shown in orange.
FIG. 3.
FEM simulations of the displacement field of the (a) 8.3 MHz longitudinal and (b) 5.8 MHz shear modes, with dashed lines indicating the effective mode volumes. (c) and (d) show the corresponding mixer output voltage in blue and S11 reflection measurements in green using an impedance analyzer across the split-post. Loaded Q-factors are determined by curve fits, as shown in orange.
Close modal
We can determine the linewidth of the mechanical modes from the mixer output voltage and obtain Q-factors. These values closely agree with the previously measured Q at 4 K using a parallel cathode plate.40 Figure 4 shows measured mechanical Q-factors at 4 K for various A and B modes.
FIG. 4.
Q values for the A3,0,0 and A5,0,0 longitudinal and B5,0,0, B7,0,0, and B9,0,0 shear modes, obtained from mixer trace output. The Q-factors demonstrate insignificant power dependence until the onset of duffing nonlinearity. The Q × f product of the best modes is of order 1014.
FIG. 4.
Q values for the A3,0,0 and A5,0,0 longitudinal and B5,0,0, B7,0,0, and B9,0,0 shear modes, obtained from mixer trace output. The Q-factors demonstrate insignificant power dependence until the onset of duffing nonlinearity. The Q × f product of the best modes is of order 1014.
Close modal
The effective mode mass for various A and B modes of the BAW vary slightly from each other, depending on the value of their lumped circuit parameters. These values were calculated using COMSOL finite element modeling, similar to the analysis for quartz, as detailed in the works of,12,26 and related analytical effective mass calculations,50,51 with details given in the Appendix. Various determined parameters for each of the mechanical modes are presented in Table I, where co-operativity C0 is defined as , Γm is the linewidth of the mechanical resonance and κc is the linewidth of the microwave cavity.1 As observed, different modes show different coupling rates due to different displacement polarizations and different overlaps between the re-entrant post-mode and the mechanical modes.
TABLE I.Comparison of the parameters for some BAW modes. Frequencies, fm, in MHz;40 Meff,m in gram, evaluated from finite element modeling; g0, the coupling rate in Hz; and the single photon cooperativity, C0.
Xn,m,p
. | fm
. | Meff,m
. | g0
. | C0
. |
---|
B5,0,0 | 4.20 | 3.53 × 10−4 | 2.67 × 10−7 | 2.40 × 10−20 |
A3,0,0 | 4.70 | 4.15 × 10−4 | 1.44 × 10−5 | 2.40 × 10−17 |
B7,0,0 | 5.80 | 3.14 × 10−4 | 2.12 × 10−8 | 2.96 × 10−23 |
B9,0,0 | 7.40 | 2.77 × 10−4 | 4.79 × 10−8 | 7.62 × 10−22 |
A5,0,0 | 8.30 | 3.52 × 10−4 | 2.38 × 10−6 | 1.96 × 10−18 |
Xn,m,p
. | fm
. | Meff,m
. | g0
. | C0
. |
---|
B5,0,0 | 4.20 | 3.53 × 10−4 | 2.67 × 10−7 | 2.40 × 10−20 |
A3,0,0 | 4.70 | 4.15 × 10−4 | 1.44 × 10−5 | 2.40 × 10−17 |
B7,0,0 | 5.80 | 3.14 × 10−4 | 2.12 × 10−8 | 2.96 × 10−23 |
B9,0,0 | 7.40 | 2.77 × 10−4 | 4.79 × 10−8 | 7.62 × 10−22 |
A5,0,0 | 8.30 | 3.52 × 10−4 | 2.38 × 10−6 | 1.96 × 10−18 |
The displacement of the crystal structure is calculated from the electrical current arising from the applied piezoelectric voltage using Eqs. (4) and (5). The frequency shift from the phase bridge measurement divided by the calculated piezoelectric induced displacement provides the frequency pull factor, as detailed by Eqs. (6) and (7). By measuring Δfc for different input voltages to the piezoelectric BAW (converted to displacements), coupling rates at 4 K for different modes were calculated from the gradients of dfc/dx. Figure 5 shows coupling rates for A and B modes of LiNbO3-BAW, compared to other mechanical resonators. The inset figure shows dfc/dx for the A5,0,0 mode and is equal to 4.4 × 1013 Hz/m. The coupling rate g0 is derived from this value by multiplying it with xzpf.
FIG. 5.
Coupling rates for various A (blue) and B (orange) modes of the LiNbO3-BAW at 4 K, compared to other acoustic resonator coupling rates.6,26,39 The inset figure shows the change in microwave frequency Δfc (Hz) corresponding to displacement Δx for A5,0,0 mode, the slope of which multiplied by xzpf gives g0.
FIG. 5.
Coupling rates for various A (blue) and B (orange) modes of the LiNbO3-BAW at 4 K, compared to other acoustic resonator coupling rates.6,26,39 The inset figure shows the change in microwave frequency Δfc (Hz) corresponding to displacement Δx for A5,0,0 mode, the slope of which multiplied by xzpf gives g0.
Close modal
The value of these coupling rates is benchmarked against similar experiments using a quartz BAW at 4 K with a single post-re-entrant resonator, as described in Refs. 26 and 39. The coupling rates g0 for LiNbO3-BAW in the split-post re-entrant cavity are observed to be up to four orders of magnitude greater compared to quartz for the longitudinal mode at A3,0,0, demonstrating higher coupling rates and making it a suitable candidate for further studies. Additionally, we compared our results with a Si3N4 membrane inside a re-entrant cavity, which achieved coupling rates of 25 mHz.6 These findings motivate further investigation into the coupling of LiNbO3-BAW with photonic systems.
The degree to which the electromagnetic and mechanical mode volumes overlap plays a key role in determining the optomechanical coupling. The improvement in coupling rates resulting from the use of a split-post vs a single-post cavity can be estimated by the ratios of a so-called overlap factor. This factor can be determined from the percentage of electromagnetic mode that exists within the mechanical mode volume determined by the dashed lines in Fig. 3. The A5,0,0 mechanical mode, for example, has all of its displacement field within a diameter of ∼6.0 mm. Inside this volume, the electromagnetic field distribution for both the split- and single-post cavities is primarily determined by the Ez field, which makes up of the electrical energy, which is only 5% in the case of a single post. The amount of electrical field within the mechanical mode volume is 20 times larger for a split-post design compared to the single-post design, explaining one of the contributions to the observed coupling rate increase of these results compared to the results in quartz.26,39 The remaining contribution stems from the material properties.
In conclusion, we investigated the acoustic modes of a LiNbO3-BAW resonator by directly exciting them with a radio frequency external voltage at mechanical resonance through piezoelectricity at 4 K. Coupling rates of order 0.014 mHz were obtained for the A3,0,0 longitudinal mode, four orders of magnitude larger than previous results obtained using quartz. Therefore, with further development, it could be utilized as an alternate material in quantum sensing, communication and testing of fundamental physics.
ACKNOWLEDGMENTS
This research was supported by the ARC Center of Excellence for Engineered Quantum Systems (EQUS, Grant No. CE170100009) along with support from the Defense Science and Technology Group (DSTG) as part of the EQUS Quantum Clock Flagship program. Additional support was provided by the ARC Center of Excellence for Dark Matter Particle Physics (CDM, Grant No. CE200100008).
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
S. Parashar: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Supervision (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). W. M. Campbell: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Supervision (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). J. Bourhill: Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Supervision (equal); Validation (equal); Writing – review & editing (equal). E. Ivanov: Formal analysis (equal); Investigation (equal); Methodology (equal); Supervision (equal); Validation (equal); Writing – review & editing (equal). M. Goryachev: Conceptualization (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Supervision (equal); Writing – review & editing (equal); M.E. Tobar: Conceptualization (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing – review & editing (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.
APPENDIX: SOME ACOUSTIC MODE PROPERTIES
1. Effective mass of the mechanical modes
Evaluation of mode mass was performed using COMSOL; evaluating the potential energy integral for the 3D/2D model for a given mode is given in terms of the following equations:
The total potential energy of the mode is then given by an integral over the entire volume as represented by
Here material density is uniform and effective mass for a 3D model
This process begins by evaluating the point displacement at the calculation point, as in the case of the odd mode, where the displacement at the point of maximum displacement is used to normalize the displacement function
rm(x). The integral of the normalized displacement over the entire geometry of the resonator is then multiplied by
to obtain the effective mass (our model in COMSOL was a 2D axisymmetrical model). By investigating the potential energy U for any mode with a small volume element dV centered on
x, the potential energy centered around the mass element
dm =
ρ(
x)
dV is defined as below, where m
eff,m is defined as effective modal mass for a given mode m.
2. Power dependence of the mechanical mode
Figure 6 shows the density plot for A(5,0,0) mode where the input voltage to the cavity is varied in the range of 0.1–3.0 V, and the output of the mixer is plotted for different input power. As the power is increased, the amplitude of the sideband increases; however, there is no significant change in linewidth, which indicates no variation in the Q-factor.
FIG. 6.
Mixer output voltage for different input power.
FIG. 6.
Mixer output voltage for different input power.
Close modal
REFERENCES
1.M.
Aspelmeyer
,
T. J.
Kippenberg
, and
F.
Marquardt
, “
Cavity optomechanics
,”
Rev. Mod. Phys.
86
,
1391
(
2014
).
2.A.
Blais
,
A. L.
Grimsmo
,
S. M.
Girvin
, and
A.
Wallraff
, “
Circuit quantum electrodynamics
,”
Rev. Mod. Phys.
93
,
025005
(
2021
).
3.B.
Cuthbertson
,
M.
Tobar
,
E.
Ivanov
, and
D.
Blair
, “
Parametric back-action effects in a high-Q cyrogenic sapphire transducer
,”
Rev. Sci. Instrum.
67
,
2435
–
2442
(
1996
).
4.J. D.
Teufel
,
T.
Donner
,
D.
Li
,
J. W.
Harlow
,
M. S.
Allman
,
K.
Cicak
,
A. J.
Sirois
,
J. D.
Whittaker
,
K. W.
Lehnert
, and
R. W.
Simmonds
, “
Sideband cooling of micromechanical motion to the quantum ground state
,”
Nature
475
,
359
–
363
(
2011
).
5.D. J.
Wilson
,
V.
Sudhir
,
N.
Piro
,
R.
Schilling
,
A.
Ghadimi
, and
T. J.
Kippenberg
, “
Measurement-based control of a mechanical oscillator at its thermal decoherence rate
,”
Nature
524
,
325
–
329
(
2015
).
6.S.
Kumar
,
M.
Kenworthy
,
H.
Ginn
, and
X.
Rojas
, “
Optomechanically induced transparency/absorption in a 3D microwave cavity architecture at ambient temperature
,”
AIP Advances
14
(
3
),
035107
(
2024
).
7.H.
Xiong
and
Y.
Wu
, “
Fundamentals and applications of optomechanically induced transparency
,”
Appl. Phys. Rev.
5
,
031305
(
2018
).
8.G. S.
Agarwal
and
S.
Huang
, “
Electromagnetically induced transparency in mechanical effects of light
,”
Phys. Rev. A
81
,
041803
(
2010
).
9.X.
Yao
,
M. H.
de Jong
,
J.
Li
, and
S.
Gröblacher
, “
Long-range optomechanical interactions in sin membrane arrays
,”
arXiv:2408.03066 (
2024
).
10.T.
Palomaki
,
J.
Teufel
,
R.
Simmonds
, and
K. W.
Lehnert
, “
Entangling mechanical motion with microwave fields
,”
Science
342
,
710
–
713
(
2013
).
11.B.
Schrinski
,
Y.
Yang
,
U.
Von Lüpke
,
M.
Bild
,
Y.
Chu
,
K.
Hornberger
,
S.
Nimmrichter
, and
M.
Fadel
, “
Macroscopic quantum test with bulk acoustic wave resonators
,”
Phys. Rev. Lett.
130
,
133604
(
2023
).
12.P.
Bushev
,
J.
Bourhill
,
M.
Goryachev
,
N.
Kukharchyk
,
E.
Ivanov
,
S.
Galliou
,
M.
Tobar
, and
S.
Danilishin
, “
Testing the generalized uncertainty principle with macroscopic mechanical oscillators and pendulums
,”
Phys. Rev. D
100
,
066020
(
2019
).
13.W. M.
Campbell
,
M. E.
Tobar
,
M.
Goryachev
, and
S.
Galliou
, “
Improved constraints on minimum length models with a macroscopic low loss phonon cavity
,”
Phys. Rev. D
108
,
102006
(
2023
).
14.A.
Lo
,
P.
Haslinger
,
E.
Mizrachi
,
L.
Anderegg
,
H.
Müller
,
M.
Hohensee
,
M.
Goryachev
, and
M. E.
Tobar
, “
Acoustic tests of Lorentz symmetry using quartz oscillators
,”
Phys. Rev. X
6
,
011018
(
2016
).
15.M.
Goryachev
,
Z.
Kuang
,
E. N.
Ivanov
,
P.
Haslinger
,
H.
Müller
, and
M. E.
Tobar
, “
Next generation of phonon tests of Lorentz invariance using quartz BAW resonators
,”
IEEE Trans. Ultrason. Ferroelectrics Freq. Control
65
,
991
–
1000
(
2018
).
16.A.
Arvanitaki
,
S.
Dimopoulos
, and
K.
Van Tilburg
, “
Sound of dark matter: Searching for light scalars with resonant-mass detectors
,”
Phys. Rev. Lett.
116
,
031102
(
2016
).
17.D.
Carney
,
A.
Hook
,
Z.
Liu
,
J. M.
Taylor
, and
Y.
Zhao
, “
Ultralight dark matter detection with mechanical quantum sensors
,”
New J. Phys.
23
,
023041
(
2021
).
18.W. M.
Campbell
,
B. T.
McAllister
,
M.
Goryachev
,
E. N.
Ivanov
, and
M. E.
Tobar
, “
Searching for scalar dark matter via coupling to fundamental constants with photonic, atomic, and mechanical oscillators
,”
Phys. Rev. Lett.
126
,
071301
(
2021
).
19.V. J.
Gokhale
,
B. P.
Downey
,
D. S.
Katzer
,
N.
Nepal
,
A. C.
Lang
,
R. M.
Stroud
, and
D. J.
Meyer
, “
Epitaxial bulk acoustic wave resonators as highly coherent multi-phonon sources for quantum acoustodynamics
,”
Nat. Commun.
11
,
2314
(
2020
).
20.Y.
Chu
,
P.
Kharel
,
W. H.
Renninger
,
L. D.
Burkhart
,
L.
Frunzio
,
P. T.
Rakich
, and
R. J.
Schoelkopf
, “
Quantum acoustics with superconducting qubits
,”
Science
358
,
199
–
202
(
2017
).
21.P.
Salzenstein
,
A.
Kuna
,
L.
Sojdr
, and
J.
Chauvin
, “
Significant step in ultra-high stability quartz crystal oscillators
,”
Electron. Lett.
46
,
1433
–
1434
(
2010
).
22.G.
Tobar
,
S. K.
Manikandan
,
T.
Beitel
et al, “
Detecting single gravitons with quantum sensing
,”
Nat. Commun.
15
,
7229
(
2024
).
23.G.
Tobar
,
I.
Pikovski
, and
M. E.
Tobar
, “
Detecting kHz gravitons from a neutron star merger with a multi-mode resonant bar
,”
arXiv:2406.16898 (
2024
).
24.C.
Locke
,
M.
Tobar
,
E.
Ivanov
, and
D.
Blair
, “
Parametric interaction of the electric and acoustic fields in a sapphire monocrystal transducer with a microwave readout
,”
J. Appl. Phys.
84
,
6523
–
6527
(
1998
).
25.J.
Bourhill
,
E.
Ivanov
, and
M.
Tobar
, “
Precision measurement of a low-loss cylindrical dumbbell-shaped sapphire mechanical oscillator using radiation pressure
,”
Phys. Rev. A
92
,
023817
(
2015
).
26.J.
Bourhill
,
N. C.
Carvalho
,
M.
Goryachev
,
S.
Galliou
, and
M. E.
Tobar
, “
Generation of coherent phonons via a cavity enhanced photonic lambda scheme
,”
Appl. Phys. Lett.
117
,
164001
(
2020
).
27.M.
Goryachev
and
M. E.
Tobar
, “
Gravitational wave detection with high frequency phonon trapping acoustic cavities
,”
Phys. Rev. D
90
,
102005
(
2014
).
28.W. M.
Campbell
,
M.
Goryachev
, and
M. E.
Tobar
, “
The multi-mode acoustic gravitational wave experiment: Mage
,”
Sci. Rep.
13
,
10638
(
2023
).
29.M.
Goryachev
,
W. M.
Campbell
,
I. S.
Heng
,
S.
Galliou
,
E. N.
Ivanov
, and
M. E.
Tobar
, “
Rare events detected with a bulk acoustic wave high frequency gravitational wave antenna
,”
Phys. Rev. Lett.
127
,
071102
(
2021
).
30.S.
Singh
,
L. A. D.
Lorenzo
,
I.
Pikovski
, and
K.
Schwab
, “
Detecting continuous gravitational waves with superfluid 4He
,”
New J. Phys.
19
,
073023
(
2017
).
31.N.
Aggarwal
,
O. D.
Aguiar
,
A.
Bauswein
,
G.
Cella
,
S.
Clesse
,
A. M.
Cruise
,
V.
Domcke
,
D. G.
Figueroa
,
A.
Geraci
,
M.
Goryachev
et al, “
Challenges and opportunities of gravitational-wave searches at MHz to GHz frequencies
,”
Living Rev. Relativ.
24
,
4
(
2021
).
32.A.
Blais
,
S. M.
Girvin
, and
W. D.
Oliver
, “
Quantum information processing and quantum optics with circuit quantum electrodynamics
,”
Nat. Phys.
16
,
247
–
256
(
2020
).
33.R.
Manenti
,
A. F.
Kockum
,
A.
Patterson
,
T.
Behrle
,
J.
Rahamim
,
G.
Tancredi
,
F.
Nori
, and
P. J.
Leek
, “
Circuit quantum acoustodynamics with surface acoustic waves
,”
Nat. Commun.
8
,
975
(
2017
).
34.J.
Kitzman
,
J.
Lane
,
C.
Undershute
,
M.
Drimmer
,
A.
Schleusner
,
N.
Beysengulov
,
C.
Mikolas
, and
J.
Pollanen
, “
Free-space coupling and characterization of transverse bulk phonon modes in lithium niobate in a quantum acoustic device
,”
Appl. Phys. Lett.
123
,
224001
(
2023
).
35.A.
Schliesser
and
T. J.
Kippenberg
, “
Cavity optomechanics with whispering-gallery mode optical micro-resonators
,” in
Advances in Atomic, Molecular, and Optical Physics
(
Elsevier
,
2010
), Vol.
58
, pp.
207
–
323
.
36.C.
Locke
and
M.
Tobar
, “
Measurement of the strain-induced coefficient of permittivity of sapphire using whispering gallery modes excited in a high-Qacoustic sapphire oscillator
,”
Meas. Sci. Technol.
15
,
2145
(
2004
).
37.S.
Kotler
,
R. W.
Simmonds
,
D.
Leibfried
, and
D. J.
Wineland
, “
Hybrid quantum systems with trapped charged particles
,”
Phys. Rev. A
95
,
022327
(
2017
).
38.J.-M.
Le Floch
,
Y.
Fan
,
M.
Aubourg
,
D.
Cros
,
N.
Carvalho
,
Q.
Shan
,
J.
Bourhill
,
E.
Ivanov
,
G.
Humbert
,
V.
Madrangeas
, and
M. E.
Tobar
, “
Rigorous analysis of highly tunable cylindrical transverse magnetic mode re-entrant cavities
,”
Rev. Sci. Instrum.
84
,
125114
(
2013
).
39.N.
Carvalho
,
J.
Bourhill
,
M.
Goryachev
,
S.
Galliou
, and
M.
Tobar
, “
Piezo-optomechanical coupling of a 3D microwave resonator to a bulk acoustic wave crystalline resonator
,”
Appl. Phys. Lett.
115
,
211102
(
2019
).
40.W. M.
Campbell
,
S.
Parashar
,
M. E.
Tobar
, and
M.
Goryachev
, “
Low temperature properties of low-loss macroscopic lithium niobate bulk acoustic wave resonators
,”
arXiv:2407.17693 (
2024
).
41.I.
Marinkovic
,
M.
Drimmer
,
B.
Hensen
, and
S.
Groblacher
, “
Hybrid integration of silicon photonic devices on lithium niobate for optomechanical wavelength conversion
,”
Nano Lett.
21
,
529
–
535
(
2021
).
42.A. Y.
Cleland
,
E. A.
Wollack
, and
A. H.
Safavi-Naeini
, “
Studying phonon coherence with a quantum sensor
,”
Nat. Commun.
15
,
4979
(
2024
).
43.P.
Arrangoiz-Arriola
,
E. A.
Wollack
,
M.
Pechal
,
J. D.
Witmer
,
J. T.
Hill
, and
A. H.
Safavi-Naeini
, “
Coupling a superconducting quantum circuit to a phononic crystal defect cavity
,”
Phys. Rev. X
8
,
031007
(
2018
).
44.K. J.
Satzinger
,
Y.
Zhong
,
H.-S.
Chang
,
G. A.
Peairs
,
A.
Bienfait
,
M.-H.
Chou
,
A.
Cleland
,
C. R.
Conner
,
É.
Dumur
,
J.
Grebel
et al, “
Quantum control of surface acoustic-wave phonons
,”
Nature
563
,
661
–
665
(
2018
).
45.L.
Shao
,
M.
Yu
,
S.
Maity
,
N.
Sinclair
,
L.
Zheng
,
C.
Chia
,
A.
Shams-Ansari
,
C.
Wang
,
M.
Zhang
,
K.
Lai
, and
M.
Lončar
, “
Microwave-to-optical conversion using lithium niobate thin-film acoustic resonators
,”
Optica
6
,
1498
–
1505
(
2019
).
46.M.
Mirhosseini
,
A.
Sipahigil
,
M.
Kalaee
, and
O.
Painter
, “
Superconducting qubit to optical photon transduction
,”
Nature
588
,
599
–
603
(
2020
).
47.W.
Jiang
,
C. J.
Sarabalis
,
Y. D.
Dahmani
,
R. N.
Patel
,
F. M.
Mayor
,
T. P.
McKenna
,
R.
Van Laer
, and
A. H.
Safavi-Naeini
, “
Efficient bidirectional piezo-optomechanical transduction between microwave and optical frequency
,”
Nat. Commun.
11
,
1166
(
2020
).
48.E. A.
Wollack
,
A. Y.
Cleland
,
P.
Arrangoiz-Arriola
,
T. P.
McKenna
,
R. G.
Gruenke
,
R. N.
Patel
,
W.
Jiang
,
C. J.
Sarabalis
, and
A. H.
Safavi-Naeini
, “
Loss channels affecting lithium niobate phononic crystal resonators at cryogenic temperature
,”
Appl. Phys. Lett.
118
,
123501
(
2021
).
49.S.
Bible
,
Crystal Oscillator Basics and Crystal Selection for rfPICTM and PICmicro® Devices
(
Microchip Technology, Inc
,
2002
).
50.D.
Stevens
and
H.
Tiersten
, “
An analysis of doubly rotated quartz resonators utilizing essentially thickness modes with transverse variation
,”
J. Acoust. Soc. Am.
79
,
1811
–
1826
(
1986
).
51.J.
Shi
,
C.
Fan
,
M.
Zhao
, and
J.
Yang
, “
Variational formulation of the Stevens–Tiersten equation and application in the analysis of rectangular trapped-energy quartz resonators
,”
J. Acoust. Soc. Am.
135
,
175
–
181
(
2014
).
© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).