Extremely low loss phonon-trapping cryogenic acoustic cavities for future physical experiments.

Low loss Bulk Acoustic Wave devices are considered from the point of view of the solid state approach as phonon-confining cavities. We demonstrate effective design of such acoustic cavities with phonon-trapping techniques exhibiting extremely high quality factors for trapped longitudinally-polarized phonons of various wavelengths. Quality factors of observed modes exceed 1 billion, with a maximum Q-factor of 8 billion and Q × f product of 1.6 · 10(18) at liquid helium temperatures. Such high sensitivities allow analysis of intrinsic material losses in resonant phonon systems. Various mechanisms of phonon losses are discussed and estimated.

A Fabry-Pérot interferometer is an optical resonant system that typically consists of two parallel highly reflecting mirrors. Due to their high quality factors, a Fabry-Pérot cavity is often used to store photons for almost milliseconds while they bounce between the mirrors. Low loss is necessary in order to increase interaction between light and matter or other physical substances under study. Such systems are widely used in telecommunication, lasers, gravitational wave detection, spectroscopy, etc.

Similar principles can be applied to the field of acoustics, bulk phonons may be resonantly confined in a thin film or a crystal lattice to build a phonon cavity. In this case, resonant phonons are stored between borders that separate different media playing a role of mirrors.

Losses in acoustic resonators can be classified into two groups: first, losses attributed to intrinsic material processes, second, losses due to resonator design. The latter are mostly due to energy leakage through clamping, substrates, electrodes, etc1. The goal of the cavity design is to make “engineering losses” negligible in order to achieve a regime where losses are due to the material itself. If this regime is achieved the cavity quality factor can be increased by orders of magnitude by cooling to cryogenic temperatures.

The thin film approach (realised in Thin Film Bulk Acoustic Resonators and High Overtone Bulk Acoustic Resonators) provides technology that can guarantee smaller distances between the borders, and thus higher resonant frequencies. The main disadvantages of these systems are the need for a substrate and the inability to efficiently control the film geometry due to limitations of the deposition process. Both result in additional losses due to phonon leakage into the environment.

The crystalline lattice approach utilises structures built from pieces of single crystal material (e.g. quartz Bulk Acoustic Wave (BAW) resonator). Due to a completely different fabrication process, this technology can provide a more sophisticated design that guarantees sufficient limitation of phonon leakage. In addition, extremely low losses of such devices extend the possible frequency range to comparatively high frequencies. In this work, we show that this approach allows Q-factors of above 1 billion with a maximum of 8 billion corresponding to a Q × f product of 1.6 · 1018 at liquid helium temperatures. This is 310 times higher than the previous value achieved at this temperature2. Such high sensitivity should allow future applications of acoustical cavities to be broadened and be used as an alternative (to optical cavities) experimental base for broad range of physical experiments, including material characterisation3, quantum cavity electromechanics45, quantum metrology67 etc.

Results

Phonon excitation and trapping

Vibrations of ideal plate resonators have been studied for many decades8. Unlike their idealised counterparts, actual cavities pose two almost contrary requirements. First, they have to be coupled to the environment for effective excitation and measurements, second, they must be uncoupled from the environment in order to avoid phonon leakage. The solution to this corundrum is to excite vibrations with the electro-magnetic field while mechanically separating the crystal from the environment.

Coupling to electromagnetic measuring and excitation fields is introduced by making a vibrating body from piezoelectric material. Imposing a potential difference between two surfaces of the plate creates a gradient of the electric potential along the plate thickness. The spatial Fourier harmonics of this gradient are coupled to mechanical vibrations through the piezoelectric effect. Frequencies of these vibrations are set by plate dimensions, material properties and the order of the harmonic. This type of vibration is referred as a thickness vibration that has to be distinguish from parasitic flexure vibration modes. Detection of the plate vibration is performed by measuring the displacement current between the surfaces of the plate.

Due to the surface boundary conditions, such a cavity can sustain resonant phonons of discrete wavelengths. When describing an acoustic mode in a piezoelectric plate, the number of half-waves in the plate thickness is referred as the overtone (OT) number, and the resonance itself as an overtone. The cavity geometry allows only anti-symmetric types of vibration to be excited piezoelectrically, thus limiting experimentally observed OTs only to odd numbers.

In order to avoid any energy leakage, the vibrating part of the plate has to be mechanically separated from the environment. This is achieved by two means: first, the excitation electrodes are separated from the crystal plate by a tiny vacuum gap; second, acoustic energy is trapped in the centre of the disk plate while the disk is clamped from the sides. The trapping is done by introducing curvature into one of the plate surfaces by manufacturing a plano-convex shape910. A schematic cross-section of the cavity is shown in Fig. 1. This variable plate thickness increases the local phonon frequency along the plate disk radius closer to the edges. The corresponding energy difference along the radial direction due to this curvature can be considered as an effective potential well. So, phonons are trapped if their total transverse energy (corresponding to in-plane wave numbers) do not exceed the energy difference between the edge and the centre. This situation significantly decreases the probability to find a resonant phonon near the clamping points. Theoretical treatment of the energy trapping problem from the point of view of linear elastodynamics is given in9. The calculations predict the mode shape in the form of polynomial-normal distribution. This result has been confirmed by means of X-ray spectroscopy1112.

Figure 1

Cross section in the r − z plane of the quartz Bulk Acoustic Wave (BAW) resonator of cylindrical geometry with typical distribution of trapped phonons in the radial direction.

Details of cavity design

Tested acoustic cavities are made from Top High Quality (THQ) alpha synthetic quartz produced by Gemma Quartz & Crystal according to a hydrothermal growth process. Impurities are removed from the heart of the quartz bars before the latter are cut up in oriented blanks. Disks come from slices cut out from a synthetic quartz crystal bar according to specific angles with respect to the crystalline axis (so-called singly rotated cuts or doubly rotated cuts). The tested cavities are doubly-rotated Stress Compensated (SC) cut plates13, which is one of the most commonly used cuts. Doubly rotated cut cavities support three possible phonon polarizations with different velocities: the quasi-longitudinal (V = 6757 m/s), fast quasi-shear (V = 3966 m/s) and slow quasi-shear (V = 3611 m/s). All three types of polarisation are not pure due to crystal anisotropy. According to the contoured-resonator axis [x, y, z]14, the y-axis being normal to the plate, their polarizations are respectively: [0.226, 0.968, 0.111], [−0.198, −0.0657, 0.979], [0.957, −0.243, 0.177].

After the slices are cut they are then rounded, lapped and polished resulting in plano-convex contoured plates. Polishing is a key step of the process: it consists of a mixing of mechanical polishing and chemical etching. This last operation is mandatory to remove the “dead layer”, i.e. a damaged layer - a blend of amorphous and polycrystalline randomly-oriented matter resulting from the first process. Chemical etching is also one efficient way to control precisely the thickness of the resonator, a key parameter for defining the frequency of waves propagating along the thickness axis. At least, mechanical strains induced by machining should also be removed by an additional reheat.

The cavity design for which both electrodes are not deposited on the vibrating plate itself is known as (BVA)15. Electrodes are deposited on two quartz supports placed on each side of the resonator, and in such a way that the gap between electrodes and resonator surfaces is typically between 5 and 10 μm210. A photograph of an assembled cavity is sown in Fig. 2.

Figure 2

Photograph of the electrode and cavity mounting arrangement of the quartz BAW cavities.

In this work we compare performance of the following types of acoustical cavities:

1.08 mm thick, 13 mm diameter electrode-separated disk cavities initially designed to sustain shear vibration of 5 MHz at room temperature (manufactured by BVA Industrie)2;

0.54 mm thick, 13 mm diameter electrode-separated disk cavities initially designed to sustain shear vibration of 10 MHz at room temperature (manufactured by BVA Industrie)10;

1 mm thick, 30 mm diameter electrode-separated disk cavities with higher grade surface polishing initially designed to sustain shear vibration of 5 MHz at room temperature (manufactured by Oscilloquartz SA);

0.54 mm thick, 13 mm diameter electrode-deposited disk cavities initially designed to sustain shear vibration of 10 MHz at room temperature (manufactured at FEMTO-ST).

It has to be noted that all cavities are initially designed for room temperature operation with slow-speed thickness-shear vibrations and not for longitudinally-polarized ones. This is due to the fact that the frequency-temperature characteristic of the former exhibits a turnover point, which is not the case of the latter. Moreover at liquid helium temperatures, longitudinal phonons exhibit significantly lower losses than any shear phonons2. Direct cavity optimisation for liquid helium operation is still impossible due to unknown temperature coefficients of the elastic constants at these temperatures.

It has to be noted that two parallel excitation electrodes form a parasitic capacitor (typically close to 3 pF) creating a parallel photon branch. At higher frequencies, the number of these non-resonant photons exceeds the number of acoustic phonons decreasing the mode contrast. Hence, the parasitic capacitance uncouples the acoustic motion from the electromagnetic environment. Nevertheless, the excitation of phonons of shorter wavelengths is possible with effective electrical compensation of this shunt capacitance.

Discussion

The four cavities described above were characterised for resonant longitudinal phonons of half wavelengths varying from to of the plate thickness. Q-factors were extracted by fitting measured electric admittance of the resonator around each resonance frequency. Fig. 3 shows an example of experimental results for the 3rd type of acoustic cavities measured for the 5th OT at 3.75 K.

Figure 3

Measured impedance (absolute value and phase) of the 5th OT of the longitudinal mode (15.7 MHz) at 3.75 K for the cavity number 3.

The results (see Fig. 4) show that the cavity with the largest diameter and better polishing exhibits significantly lower losses. This fact could be attributed to better trapping (due to larger potential well) and lower scattering on surface imperfections. In addition, all four samples exhibited an increase of Q as a function of OT number, at least for a certain range of frequencies. This is explained by improving of the phonon trapping with increasing OT number as follows from the theory9. Temperature dependence of the selected overtones for this cavity is shown in Fig. 5.

Figure 4

Measured quality factors for different OTs of the (quasi-) longitudinal vibration of the four cavities under investigation at 3.75 K. Error bars for standard deviations lower than ±2.5% are not represented.

Both higher values in brackets are close to the measurement set-up limitation.

Figure 5

Temperature dependance of cavity losses for different OTs of the (quasi-)longitudinal vibration in the cavity number 3.

Total intrinsic losses of an acoustic cavity apart from lack of trapping are a sum of losses due to different mechanisms: where the Q-values represent different mechanisms. These mechanisms are explained below combining the results of frequency (Fig. 4) and temperature (Fig. 5) dependences.

First, the most important mechanism of energy loss in acoustic cavities at low temperatures and relatively high frequencies is loss mechanism due to phonon-phonon scattering over crystal lattice anharmonicity (Q in equation (1)). This mechanism was first analytically considered by Landau for the case of shear-polarized phonons1617. It was also pointed out that longitudinal phonons do not exhibit this type of losses in the case of two phonon collision, due to significantly higher sound velocity. Later the theory was extended to the case of three colliding phonons18 predicting the following law: . This law is experimentally observed for the large acoustic cavity between between 6 K and 13 K (Fig. 5). According to the theory18, the acoustic wave attenuation α is proportional to the resonance angular frequency ω as long as , and where τ denotes the average lifetime of thermal phonons, ℏ the Planck constant, and k the Boltzmann constant. As a consequence, the Q-factor of an acoustic wave propagating with velocity V would no longer depend on frequency ω. Since these conditions are met in the present experiment, the Q-factor tends to be independent of f (Q = const law) as shown in Fig. 4.

Second, for temperatures below 6 K the Q(T) dependence changes the slope tending to T−0.3 scaling law (Fig. 5). Such dependence have been observed in a various of electromechanical systems1920212223, and acoustic cavities at lower temperatures10, in particular. The T−1/3 dependence is theoretically explained by coupling to so-called “two level systems” (TLS)24 (Q in equation (1)). In synthetic quartz these TLSs are attributed to impurity ions, typically Al3+, Na+, Li+, Si4+ and H+, etc25262728 studied by electron-spin resonance, infrared spectroscopy, dielectric relaxation spectroscopy and thermoluminescence. Although the utilised material is premium pure, concentration of these ions could reach 1–3 ppb.

Third, the increase of cavity losses at higher frequencies occurs for the most of resonators (see Fig. 4), and can be explained by phonon scattering on the cavity surfaces (Q in equation (1)). Scattering originates from the surface roughness as well as from the presence of a dead layer. Indeed, the reflection coefficient on a rough surface with a roughness variance σ2, assuming a gaussian probability density, and can be expressed as29 , where k is the wave number. At the same time, propagating through the thickness h, the acoustic wave is attenuated by e−. Thus, the reflection is equivalent to a modified path causing an attenuation equal to α′ = α + ln(R/h). Consequently, for an OT number n such that kh = nπ, the associated Q-factor becomes: Since the resonant frequency ω is almost proportional to the overtone number n, this mechanism results in Q × ω = const law. The surface scattering mechanism strongly depends on the polishing quality denoted by the roughness variance σ. This explains the better performance of cavity number 3 (Fig. 5), which has a better surface polishing grading. The largest observed values of quality factor were at n = 65 for cavity number 3, which has a σ of one nm (consistent with specification of Oscilloquartz fabrication process). For this type of cavity, Q is estimated to be , which gives an upper limit of 3 × 1010 for the highest measured OT n = 65 at 204 MHz. For other cavities the polishing process is inferior with σ varying from 4 to 7 nm. For example, for cavity type 2 with a specified σ = 4 nm, Q is found to be and gives a value of 4.2 × 108 at n = 19 (118.6 MHz), which is consistent with the experimentally measured value of 3.3 × 108. Thus, we can conclude that the highest measured Q-values are close to estimated values of Q, demonstrating that the cavity roughness could be a major factor limiting the Q-factor for high frequencies.

Fourth, although thermoelastic losses do not exist for pure shear vibration, since there is no volume change, they are present for the longitudinal modes. So, the Q for longitudinal vibration can be estimated based on the description of the wave propagation including the coupled equations of motion of the medium and thermal conduction30: with λ = c, and where u is the displacement vector, c the tensor of isothermal elastic constants, α the thermal expansion tensor, Γ the thermal conductivity tensor, ρ the volumic mass, and C the specific heat at constant strain, δT the temperature change around the reference temperature T0. Although, most of these coefficients for the lightly anisotropic quartz are unknown at cryogenic conditions, simplified calculations can be made assuming an isotropic medium and propagation along the Z-axis. The optical Z-axis is the most experimentally studied case. This analysis implies calculation of a turning-point frequency f * that shows the regime of acoustic wave propagation - adiabatic or isothermal: where is the material diffusivity, and for a propagation along the Z-axis. The former is calculated based on data on thermal conductivity and heat capacity at cryogenic temperatures3132 giving 1 ≤ f * ≤ 5.5 MHz for 3 ≤ T ≤ 10 K. Over this range, all the resonant frequencies are always greater than f *, and acoustic wave propagates in the quasi-isothermal regime. For this regime, the maximum attenuation is: with , where λ33 = 2c12α11 + c33α33, because α11 = α22 and c13 = c23 in case of crystalline quartz. Whereas the thermal expansion coefficients have been measured at low temperatures31, the elastic constants on the other hand can just be roughly approximated by their values at room temperature. Nevertheless, the corresponding thermoelastic losses estimated from equation (2) remain always at least two orders of magnitude lower than those measured below 10 K. Since α∞ is almost independent of frequency for a given material, these type of losses result in Q/ω = const law.

We demonstrate that acoustic cavities with effective trapping of longitudinally polarised phonons can exhibit quality factors above 1 Billion at liquid helium temperatures and frequencies 10–200 MHz. This shows that phonon cavities can be considered as an alternative to photon cavities in many low-frequency, low-power experiments. The advantages of acoustic devices is the ease of operation and integration into higher-order electronic systems, compact size, natural way of combination of mechanical motion and electromagnetic environment, as well as artificial lattice impurities. In particular, these systems can be used as parts of quantum hybrid systems in miliKelvin temperatures integrating quantum circuit electrodynamics and cavity elastodynamics. Moreover, unlike nano-scale mechanical resonators, acoustic cavities are a macroscopic objects of considerable effective mass varying in the range of 0.1–10 milligrams for modes between 100 and 10 MHz. This, can be a promising experiment on unifying quantum mechanics and general relativity in a single experiment7.

Apart from fundamental physics, cryogenic acoustic cavities can be used for building ultrastable frequency sources for variety of applications.

Methods

The device under test is embedded inside an oxygen-free copper block fixed on the end of the cold finger of a Helium pulse-tube cryocooler2. The cold finger is temperature-controlled, and residual fluctuations of temperature are within ±2 mK around the setting temperature.

The device characterisation is performed using the impedance analyzis method: measurements of the displacement current through the resonator thickness while applying alternating potential difference between electrodes. Their ratio is the complex impedance that is measured for a range of frequencies around a resonance. This characteristic is used to calculate a phonon resonance bandwidth. For accurate measurements of intrinsic properties of a device under test (DUT), the method requires a calibration procedure in order to compensate any influence of a long connecting cables. For this purpose, three identical extra cables follow the same path to be ended respectively by an open-circuit, a short circuit and a 50 Ohms load. Measurements of these etalons helps to cancel cable influence for certain values of the analyzer output power, frequency span and sweep time.

The excitation power is kept as low as possible in order to avoid nonlinear and thermal effects. Power dissipation inside the resonator is estimated to vary from 0.1 nW to 3 nW, depending on a mode coupling and quality factor. Nonlinear effects of the investigated devices are due to lattice anharmonicity. They strongly depend on quality factor of a mode and can be noticeable at powers as low as 10 nW.

Special attention should be paid to the accuracy when characterising modes with extremely high quality factors due to extreme narrowness of bandwidths (below 100 mHz). The frequency span of the analyzer is adjusted to get the maximum benefit of the number N of measurement points. The minimum available span of 0.5 Hz for 401 measured points provides resolution of about 1.25 mHz. The minimum sweep time required to avoid all transient effects is (N − 1) × τ where τ is phonon relaxation time of a mode of frequency f. It is given as where Q is a mode loaded quality factor, f is the resonant frequency.

Based on the specifications of the measurement equipment and the post-processing procedure, the measurement uncertainty of the intrinsic (unloaded) quality-factor Q is evaluated as: which demonstrates importance of the fractional frequency stability of both the measurement apparatus frequency source and the DUT itself. The latter stability is improved by locking to an ultra-stable frequency source, such as a Hydrogen maser. The long term (due to requirement for large sweep time) frequency stability of the resonator mode is determined predominantly by its temperature sensitivity and does not typically exceed 10−10 providing measurement uncertainty below a few percent for all characterised modes.

Author Contributions

S.G. and M.G. made cavity characterisation. P.A. participated in organising cryogenic experiments. J.P.A. prepared cavity 4 and consulted on its design. R.B. consulted on cavity design and manufacturing processes. S.G., M.G., R.B. and M.T. discussed loss mechanisms. MG proposed the terminology. S.G., M.G. and M.T. wrote the main manuscript text. M.T. supervised the work. All authors reviewed the manuscript.