The Kerr Parametric Oscillator (KPO) is a nonlinear resonator system that is often described as a synthetic two-level system. In the presence of noise, the system switches between two states via a fluctuating trajectory in phase space, instead of following a straight path. The presence of such fluctuating trajectories makes it hard to establish a precise count or even a useful definition, of the “lifetime” of the state. Addressing this issue, we compare several rate counting methods that allow to estimate a lifetime for the levels. In particular, we establish that a peak in the Allan variance of fluctuations can also be used to determine the levels' lifetime. Our work provides a basis for characterizing KPO networks for simulated annealing where an accurate determination of the state lifetime is of fundamental importance.

Synthetic two level systems (TLSs) generated in driven nonlinear resonators have recently caught a significant attention in the physics community.^{1,2} A particularly prominent example is the Kerr Parametric Oscillator (KPO, also known as parametron)^{3–15} whose potential energy is pumped at frequency *f _{p}* close to twice its resonance frequency

*f*

_{0}, i.e., at $fp\u22482f0$. If the modulation strength

*λ*exceeds a threshold $\lambda th$, the device responds with oscillations locked to $fd\u2261fp/2$ within a certain detuning range. This well-known “period doubling” of the response relative to the pump gives rise to two stable “phase states” with the same amplitude but separated by a phase difference of

*π*. The phase states can be used to encode the two polarization states (up/down) of a classical spin. This analogy leads to the idea of using networks of coupled KPOs to build noisy intermediate-scale quantum (NISQ) machines.

^{16,17}These machines can simulate the dynamics of mathematical problems that overwhelm traditional computers, such as the ground state of an Ising Hamiltonian,

^{18–27}or of other complex systems that can be mapped onto the same framework.

^{28–32}

An important quantity for many applications of TLSs is their lifetime *τ.*^{33} It is the typical time spent on a level before the interaction with an environment induces a (seemingly) spontaneous “jump” from one state to the other. The rates of environmental noise-induced switching have previously been investigated for different systems, such as trapped electrons,^{34} cold atoms,^{35} micromechanical systems,^{36–38} and analog electronic circuits.^{39}

The situation is more subtle for a KPO. Here, the synthetic levels are formed by coherent bosonic states forming attractors in phase space. These attractors are not separated by an energy gap but by a phase gap.^{11} When switches occur on a slow timescale (relative to the resonator relaxation time) and follow narrow channels in phase space, the fluctuations are termed as “weak.” Such a setting allows for situations with negligible back action where the fluctuations during a single switch can be observed. Currently, however, there exist very few studies of the fascinating physics unfolding during individual switches.^{37,40–42}

In this work, we study a classical micromechanical KPO and investigate its switching rates in the presence of weak fluctuations. We invoke and compare several methods previously used to characterize the rates of charge and parity state switching in cooper pair boxes and superconducting qubits.^{43,44} Furthermore, we propose a method to calculate the switching rate that is based on the Allan variance of the resonator displacement.^{45} In the final part of the paper, we compare all methods and find good agreement between several (but not all) of them.

Our KPO consists of a micro-electromechanical resonator (MEMS) in a room-temperature setup schematically shown in Fig. 1(a). The resonator is a doubly clamped beam, with the length of 200 *μ*m, width 3 *μ*m, and 60 *μ*m in thickness with a lumped mass of 25.4 ng made from highly doped single crystal silicon and fabricated in a wafer-scale encapsulation process.^{46} Electrodes on both sides separated from the conducting beam with a gap $\u2248$1 *μ*m enable capacitive driving and sensitive detection of oscillations in the presence of a bias voltage, $Vbias=10\u2009V$.^{47} We use a Zurich Instruments HF2LI lock-in amplifier to apply the driving voltage $Vin$ and to readout the resonator displacement $x\u221dVout=u\u2009cos(\omega t)\u2212v\u2009sin\u2009(\omega t)$ with quadrature amplitudes *u* and *v*. For convenience, we drop the proportionality factor between *x* and $Vout$ and identify in the following *x* $\u2261Vout$.^{10}

Our mechanical resonator can be described by the nonlinear equation of motion (in units of the measured electrical signal)

Here, dots indicate time derivatives, $\omega 0/2\pi =f0=439.56\u2009kHz$ is the resonance frequency, $\alpha =1.47\xd71018\u2009V\u22122\u2009s\u22122$ is the coefficient of the Duffing nonlinearity, $\gamma =\omega 0/Q=770\u2009Hz$ is the resonator relaxation rate, and *Q* = 3580 is the quality factor of the resonator. The potential energy term ($\u221dx$) is pumped with the parametric modulation depth $\lambda =2Vin/(VthQ)$ at the angular frequency $2\omega d=4\pi fd$ where $Vth=320\u2009\mu V$ is the voltage threshold for parametric oscillations for the case $fd=f0$ (demodulation frequency). The potential modulation arises because the electrostatic force due to $Vin$ pulls the beam closer toward one electrode. The force is nonlinear, i.e., it grows stronger for small beam-electrode distances, which corresponds to a change in the overall spring constant that the beam experiences. As a consequence, the drive generates small frequency variations $\delta f0\u221dVin$. The term *ξ* in Eq. (1) represents a fluctuating thermal bath (see the supplementary material for details.).

Figure 1(b) shows the *v*-quadrature response of the resonator during two sweeps of *f _{d}* from positive to negative detuning $\Delta \u2261fd\u2212f0$. Close to $\Delta =50\u2009Hz$, the response jumps from

*v*= 0 to $v=\xb150\u2009\mu V$, marking a bifurcation point of the underlying nonlinear system. At the bifurcation, the resonator experiences a spontaneous $Z2$ symmetry breaking, also known as a period-doubling bifurcation or a discrete time-translation breaking.

^{48,49}At this point, the resonator jumps to a positive or negative response with equal probability. The two responses belong to stable attractors (1 and 2) with opposite phases, i.e., $v1=\u2212v2$ (and $u1=\u2212u2$).

^{27,50–52}

To study switching between the phase states of our KPO, we apply white electrical noise *ξ* characterized by a standard deviation *σ _{V}* (over a bandwidth of 30 MHz) that causes the state of the resonator to fluctuate around its initial solution. If the fluctuations are large enough, they will occasionally carry the resonator across the threshold in the middle between the phase states. The resonator is then captured by the opposite attractor, corresponding to a switch of the synthetic TLS. Several such processes can be observed in Fig. 1(c). From this observation, it appears natural to attribute a lifetime to the inverse switching rate, $\tau =\Gamma \u22121$. However, calculating the switching rate is not straightforward due to the fluctuating trajectory.

For a deeper understanding of the system's transient behavior during switching events, we perform measurements with a high temporal resolution. In Figs. 2(a) and 2(b), we display a narrow time segment before, during, and after a single switch. We find many data points in the unstable zone between the two phase states. A ten-point moving average filter helps to visualize the trajectory of the system during the transition. The total switching time is roughly 10 ms, much longer than the lock-in integration time of 15 *μ*s and the moving-average filter time of 700 *μ*s. The measurement error of each data point is 3.7 *μ*V, in agreement with the measured point-to-point fluctuations, but significantly smaller than the $\u223c10\u2009\mu V$ fluctuations visible on the 5 ms scale.

Our observation depicted in Fig. 2 demonstrates that activated switches between the phase states are not deterministic but include prominent random elements. For instance, in the phase-space representation of the switch in Fig. 2(b), we can clearly see that the system performs a winding path close to the origin. In our device, the fluctuations generally have a slight preference for counterclockwise rotations around the phase states and clockwise ones around the origin. This can be explained by the combination of the drive and the nonlinearity, which leads to an effective detuning of the fluctuations from the lock-in amplifier clock.^{53} In the corresponding Fokker–Planck steady-state calculation presented in Fig. 2(c), we, therefore, find a channel with a significant probability density between the phase states.

These visualizations of the fluctuating trajectories expose a fundamental problem in estimating the lifetime *τ*: since transitions follow no straight lines, they can cross any point in phase space multiple times during a single switching event. An example of this can be observed in Fig. 2(b), where the averaged (dark) trajectory crosses the dotted threshold line from bottom to top, describes a clockwise winding that traverses back across the threshold, and finally crosses the line a third time before completing the switch. A simple counting algorithm will in this case register three crossing events during a single switch. In general, any counting method based on a simple threshold (such as a line) will, therefore, overestimate the switching number $Nswitch$ during the full time *T*, and therefore, also $\Gamma =Nswitch/T$. This problem has been known since a long time.

The problem of overestimating the switching count can be reduced by defining multiple thresholds that have to be crossed in a particular order to constitute an event. In Fig. 3(a), we demonstrate this with the example of two circles in phase space. The count is increased by one each time a circular threshold is left and the opposing circle is entered. This method is less sensitive to small fluctuations, but it requires a subjective measure that impacts the estimated Γ, in our case, the radii of the circular thresholds. Calibrating the measured switching rate Γ as a function of the radius helps to reduce this degree of arbitrariness (see the supplementary material), but it cannot be removed entirely.

To avoid overcounting and subjective dependencies, it is desirable to extract Γ from a method that does not require thresholds at all. Interestingly, the parity lifetime of superconducting qubits can be determined via their charge-parity power spectral density (PSD).^{43,54,55} Assuming that the switching is dominated by telegraph noise, the PSD of *v* of our KPO can be fitted to a Lorentzian function,

where the lifetime *τ* corresponds to the characteristic timescale between level switching events, and *F* is a constant related to the measurement fidelity.^{55} In this case, the lifetime or the switching rate is related to the width of the spectral peak in the frequency domain.^{56} To make the estimate quantitative in Fig. 3(b), we fit the measured displacement power spectral density with Eq. (2), yielding a third estimation for $\Gamma =1/\tau $ lifetime. The method can also be applied after a Fourier transform by fitting the sliding average autocorrelation with the function $AC(\Delta t)=Ae\u22122\Delta t\Gamma $ under the assumption of stationarity and ergodicity (not shown).

Crucially, the autocorrelation is intimately related to the Allan variance (see the supplementary material for the derivation). Originally invented to characterize the fidelity of clocks, the Allan variance measures the frequency fluctuations of a resonator as a function of integration time *τ _{A}*. As we are interested in the time

*τ*over which the typical fluctuations of

*u*(or

*v*) of our KPO are maximal, we apply the Allan variance formalism

^{57}to the measured values,

In this notation,

are sums over the measured *v* values (or *u* values) and $\u27e8\u20db\u27e9i$ denotes the mean over *i*, running from *i* = 1 to $i=N\u22122\tau A/\delta t$, where *N* is the total number of data points and *δt* is the sampling time. Assuming that the signal is dominated by telegraph-like switching with lifetime *τ* and amplitude *B*, we obtain^{45}

Hence, the maximum of $\sigma Allan2(\tau A)$ should occur around the value $\tau A\u2248\tau =\Gamma \u22121$. In Fig. 3(c), we, indeed, find a peak at the expected value, yielding $\Gamma \u22484\u2009Hz$. In contrast to the PSD method, the Allan variance method does not necessarily require a fitting process, as the peak can be read off directly and is easy to interpret even in the presence of noise.

We compare the results of the different methods in Fig. 4. We find excellent agreement between four out of five of the methods for values of Γ varying over more than two orders of magnitude. The only method that we wish to discard from this comparison is the simple line threshold approach, which consistently overestimates the count rate as expected from the discussion above. The method using two circles for thresholding overestimates Γ slightly for $Vin<0.4\u2009V$, where the separation between the attractors is small and the “clouds” start to overlap significantly, cf. the example in Fig. 3(a). Additional comparison as a function of the noise strength *σ _{V}* can be found in the supplementary material.

We emphasize that there is no fundamental reason why the estimators we obtain should be identical at all. The surprisingly good agreement between most of the estimators confirms that the notion of a lifetime *τ* is useful to characterize the switching between phase states in a KPO, where a parametric pump generates a synthetic potential landscape.^{49} This approach may be useful in other systems where multi-stable potentials in dimensions higher than one are present, such as three-dimensional protein folding or other chemical reactions. For advanced applications in the future, the resonator networks could be realized through bilinear, resonant coupling of several KPOs^{25,53} (see the supplementary material for details). For MEMS, such as those studied here, bilinear coupling can be achieved in multiple ways, such as pairwise capacitive, inductive, optical, or mechanical coupling, or indirect all-to-all coupling through a separate radio frequency cavity.

See the supplementary material for theory derivations for the probability density, the Allan variance, and the autocorrelation of telegraph noise; experimental demonstrations of the dependence of the extracted switching rate on the circle threshold radius and on the noise strength; and a short summary of various coupling methods for parametric oscillators.

Fabrication was performed in nano@Stanford labs, which are supported by the National Science Foundation (NSF) as part of the National Nanotechnology Coordinated Infrastructure under Award No. ECCS-1542152, with the support from the Defense Advanced Research Projects Agency's Precise Robust Inertial Guidance for Munitions (PRIGM) Program, managed by Ron Polcawich and Robert Lutwak. This work was further supported by the Swiss National Science Foundation through Grant Nos. CRSII5_177198/1 and PP00P2 190078 and a Deutsche Forschungsgemeinschaft (DFG) - Project No. 449653034.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

### Author Contributions

**Gabriel Margiani:** Data curation (equal); Investigation (equal); Resources (equal); Software (equal); Visualization (equal); Writing – original draft (equal). **Nicholas E. Bousse:** Resources (equal); Writing – review & editing (equal). **Thomas W. Kenny:** Funding acquisition (equal); Supervision (equal); Writing – review & editing (equal). **Oded Zilberberg:** Formal analysis (equal); Supervision (equal); Writing – original draft (equal). **Deividas Sabonis:** Supervision (equal); Writing – original draft (equal). **Alexander Eichler:** Conceptualization (equal); Project administration (equal); Supervision (equal); Writing – original draft (equal). **Sebastián Guerrero:** Resources (equal); Software (equal); Writing – review & editing (equal). **Toni Louis Heugel:** Data curation (equal); Formal analysis (equal); Writing – review & editing (equal). **Christian Marty:** Resources (equal); Writing – review & editing (equal). **Raphael Pachlatko:** Resources (equal); Writing – review & editing (equal). **Thomas Gisler:** Investigation (equal); Methodology (equal); Writing – review & editing (equal). **Gabrielle D. Vukasin:** Resources (equal); Writing – review & editing (equal). **Hyun-Keun Kwon:** Resources (equal); Writing – review & editing (equal). **James M Lehto Miller:** Conceptualization (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.