Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Electro-superlubric springs for continuously tunable resonators and oscillators


Resonators and resonator-based oscillators are used in most electronics systems and they are classified as either mechanical or electrical, with fixed or difficult-to-tune resonant frequencies. Here, we propose an electro-superlubric spring, whose restoring force between two contacting sliding solid surfaces in the structural superlubric state is linearly dependent on the sliding displacement from the balanced position. We use theoretical analysis and finite element methods to study the restoring force and stability. The stiffness of this electro-superlubric spring is proportional to the square of the applied electric bias, facilitating continuous tuning from zero to several megahertz or gigahertz for the microscale or nanoscale resonators, respectively. Furthermore, we propose an electro-superlubric oscillator that is easily operated by varying a pair of harmonic voltages. The resonant frequency, resonant amplitude, quality factor, and maximum resonant speed can be continuously tuned via the applied voltage and bias. These results indicate significant potential in the applications of electro-superlubric resonators and oscillators.


Microscale resonators, which can be classified into either electric or acoustic type1,2, have attracted considerable research interest due to their broad applications in electromechanical filters3, gyroscopes4, resonant accelerometers5 and energy harvesting6,7,8. All electric resonators, like L-C resonators, are based on the electronic oscillations in the circuit. Acoustic resonators, which mainly contain microelectromechanical system (MEMS), nanoelectromechanical system (NEMS), and quartz resonators, are based on the interaction of mechanical motion and thermal or acoustic phonons1. Generally, electric resonators have better performance on the integration with IC process than acoustic resonators, while acoustic resonators have much higher quality factor (Q-factor) and lower energy dissipation than electric acoustic resonators, and thus acoustic resonators are widely used in many applications, including Si MEMS clock oscillators in Apple iPhone1.

Large-range tunability of the resonant frequency is among the most desirable functions of resonators. For instance, with emerging wireless communication systems, an extreme resonant frequency tunability enables a single resonator to cover a wide range frequency; such a resonator could replace a series of frequency-fixed resonators, which will reduce costs as well as the size of the device, thereby offering unprecedented configurability9,10. Furthermore, extreme resonant frequency tunability promises an extremely wide range of harvesting energy6, as the resonance frequency of such a device will be able to adapt to the external energy excitation frequency, which means that much more energy can be efficiently harvested than a fixed or a small-range resonant frequency resonator. There are several techniques for tuning the resonant frequencies of resonators, including piezoelectric6,11,12,13, electrothermal3,14,15,16,17 and electrostatic18,19,20,21 techniques. Although these existing resonators can tune a certain range of frequencies around the natural resonant frequency, the tuning range is not extreme and the tuning techniques contain additional structures, which is inconvenient for integrated fabrication.

Structural superlubricity (SSL) refers to a state of ultralow friction and zero wear between two contacting solid surfaces22,23,24. Since its first realisation in microscale graphite22, the SSL state has been observed in contact between microscale graphite and hBN, DLC (diamond like carbon), Au, etc.25,26,27, at high speeds (up to 294 \({{{{\rm{ms}}}}}^{-1}\))28 and for long ranges of sliding (>100 km)25.

In this paper we propose an electro-superlubric spring (ESL-spring), whose restoring force between two contacting sliding solid surfaces in the structural superlubric state is linearly dependent on the sliding displacement from the balanced position. We show that the stiffness of this ESL-spring is proportional to the square of the applied electric bias, which provides a highly effective solution to realise extreme tunable range. Furthermore, we propose an ESL-oscillator that is easily operated by varying a pair of harmonic voltages, indicating significant potential in the applications of ESL-resonators and oscillators.


The principle of ESL-spring

Figure 1a, c illustrate the side and top views of the elementary structure of the ESL-spring, respectively. We consider an elementary setup, where a rectangular electrode (SLIDER) of length, \(L\), and width, \(W\), and a significantly larger substrate are in contact in the SSL state. The substrate is covered or coated by a dielectric thin film of thickness, \(d\), and relative permittivity, \({\varepsilon }_{{{{\rm{r}}}}},\) on the structure of two rectangular electrodes (Electrodes 1 and 2) of the same width \(W\) and a larger and thicker dielectric material. Electrodes 1 and 2 are placed in a line, namely, on the \(x\)-axis (Fig. 1c), and separated by a distance \(s\). In practice, we can transfer a microscale graphite flake, which is conductive, with the self-retracting motion property29 to an atomically smooth hexagonal boron nitride (hBN) film, which is a dielectric for forming a robust SSL contact27. Silicon dioxide (SiO2), which is an insulator and can be grown directly via chemical vapour deposition, provides another choice of dielectric thin film.

Fig. 1: Elementary structure, driving force of ESL-spring, and its numerical validation.
figure 1

a, c Illustration of ESL-spring viewed from side and top, respectively, where key variables and parameters are indicated. b Equivalent circuit diagram of ESL-spring. df are the two-dimensional FEM simulation results for the potential and electric fields as functions of sliding displacement of the SLIDER (the upper electrode), with \(L=W=4{\,{\upmu }}{{{\rm{m}}}}\), \(H=200\,{{{\rm{nm}}}},s=1\,{{\upmu }}{{{\rm{m}}}}\), and \({V}_{{{{\rm{b}}}}}=40{{{\rm{V}}}}\), and taking silicon oxide as Dielectrics 1 and 2 \(({\varepsilon }_{{{{\rm{r}}}}}=3.9)\). d Changes in the potential field (left panel) and electric field (right panel) with five symmetrical displacements for \(d=100\,{{{\rm{nm}}}}.\) e, f Dimensionless equivalent potential energy and driving force as functions of the relative sliding displacement \(\xi =2x/{L}_{{{{\rm{c}}}}}\) with respect to four different width-thickness ratios, \(W/d\); the approximation accuracy is excellent for \(W/d \, > \,40\).

The SLIDER moves on the substrate (STATOR) in the \(x\)-direction in the range of \(|x| < {L}_{{{{\rm{c}}}}}/2\). By applying a constant electric bias \({V}_{{{{\rm{b}}}}}\) to Electrodes 1 and 2, we demonstrate that the slider will meet a linear restoring force, \(F=–kx\), where \(x\) is the displacement from the balanced position. The slider is located on the upper right over the middle of Electrodes 1 and 2, and \({L}_{{{{\rm{c}}}}}=L–s\) is the overlapping length of SLIDER with Electrodes 1 and 2. Hence, the setup shown in Fig. 1a, c corresponds to a linear spring in the displacement range of \(|x| < {L}_{{{{\rm{c}}}}}/2\). To demonstrate the above result, as shown in the equivalent circuit diagram in Fig. 1b, we note that the slider forms capacitances \({C}_{1}={C}_{0}(1–\xi )/2\) and \({C}_{2}={C}_{0}(1+\xi )/2,\) with Electrodes 1 and 2, respectively. Here, \({C}_{0}={C}_{1}+{C}_{2}={\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}{L}_{{{{\rm{c}}}}}W/d\) is the total capacitor, \(\xi =2x/{L}_{{{{\rm{c}}}}}\in (–1,1)\) is the relative displacement, and \({\varepsilon }_{0}\) is the permittivity of vacuum. As detailed in the Supplementary Note 1, the equivalent potential energy, \(U\), of the entire system can be calculated as:

$$\begin{array}{c}U=–\frac{1}{2}\frac{{C}_{1}{C}_{2}}{{C}_{1}+{C}_{2}}{V}_{{{{\rm{b}}}}}^{2}=–(1–{\xi }^{2})\frac{{C}_{0}}{8}{V}_{{{{\rm{b}}}}}^{2}.\end{array}$$

Substituting (1) into the relationship \(F=–\partial U/\partial x\), we obtain the restoring force \(F=–kx\) with the following ESL-spring coefficient:

$$\begin{array}{c}k=\frac{{C}_{0}}{{L}_{{{{\rm{c}}}}}^{2}}{V}_{{{{\rm{b}}}}}^{2}=\frac{{\varepsilon }_{0}{\varepsilon }_{r}W}{{L}_{{{{\rm{c}}}}}d}{V}_{{{{\rm{b}}}}}^{2}.\end{array}$$

Thus, we introduce a new spring type, called an ESL-spring, with a continuously tunable spring coefficient by changing the applied voltage, \({V}_{{{{\rm{b}}}}}\).

To validate the above theoretical relationships, we carried out a finite element simulation on the electric and potential fields, equivalent potential energy, and driving force as functions of the displacement, \(x\), as detailed in the Supplementary Note 2. We assume \(L=W=4{{{\rm{\mu }}}}{{{\rm{m}}}}\), \(d=100{{{\rm{nm}}}},s=1{{{\rm{\mu }}}}{{{\rm{m}}}},\) and \({V}_{{{{\rm{b}}}}}=40{{{\rm{V}}}}\), and set the thicknesses of the slider and substrate at \(200{{{\rm{nm}}}}\) and 4\({{{\rm{\mu }}}}{{{\rm{m}}}}\), respectively. Using SiO2 as the dielectric \(({\varepsilon }_{{{{\rm{r}}}}}=3.9)\), we obtain the corresponding simulation results on the potential field (left panel) and electric field (right panel) for five symmetrical displacements of the ESL-spring, as depicted in Fig. 1d. The electric field (potential difference) is mainly concentrated in the overlap area between the slider and Electrodes 1 and 2, which means that the electrostatic energy is mainly composed of capacitors \({C}_{1}\) and \({C}_{2}\). Figure 1e, f show the simulation results (solid lines) of the dimensionless equivalent potential energy, \(U/{U}_{0}\), and driving force, \(F/{F}_{0}\), with respect to four different width-to-thickness ratios, \(W/d\) (with the same \(W=4{{{\rm{\mu }}}}{{{\rm{m}}}}\)), where \({U}_{0}={C}_{0}{V}_{{{{\rm{b}}}}}^{2}/8\) and \({F}_{0}={C}_{0}{V}_{{{{\rm{b}}}}}^{2}/2{L}_{{{{\rm{c}}}}}\). The linear relationship \(F=–kx\) (black dashed line) is excellent at \(W/d \, > \,40\) because edge effects decrease with increasing \(W/d\). This condition is very easy to achieve through micro-nano-fabrication technology in the SSL state. Further, the linear relationship holds only within the range of \(|x| < {L}_{{{{\rm{c}}}}}/2\).

However, nonzero frictions in SSL contacts are inevitable, although they are ultralow (by definition). In ambient environments, typical friction forces affecting a sliding 4-\({{{\rm{\mu }}}}\)m-sized graphite flake on substrates in the SSL contacts were measured to be in the range of \(0 \sim 1.0{{{\rm{\mu }}}}{{{\rm{N}}}}\)22,23,27,30. Recently, it was revealed that these frictional forces are mainly caused by interactions between the edges of the flake and adsorptions of the substrates and can be mostly removed through an annealing process31. In comparison, the maximum ESL-spring force is equal to \({F}_{{{\max}}}={\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}W{V}_{{{{\rm{cr}}}}}^{2}/2d\) at \(\xi =\pm 1\) whenever \({V}_{{{{\rm{b}}}}}\) takes its maximum allowed voltage, \({V}_{{{{\rm{cr}}}}}={E}_{{{{\rm{cr}}}}}d\), before the electric breakdown of the dielectric film, where \({E}_{{{{\rm{cr}}}}}\) is the breakdown strength of the dielectric film. According to the theory and experiments of the solid dielectric breakdown32, when the thicknesses of dielectric films are thinner than ~\(4{{{\rm{\mu }}}}{{{\rm{m}}}}\), \({E}_{{{{\rm{cr}}}}}\) and the relative permittivity, \({\varepsilon }_{{{{\rm{r}}}}}\), for most dielectric films can be approximately correlated to a unified form of \({E}_{{{{\rm{cr}}}}}\approx {E}_{0}/\sqrt{{\varepsilon }_{{{{\rm{r}}}}}}\) with \({E}_{0}=2.24\times {10}^{9}{{{{\rm{Vm}}}}}^{–1}\). Therefore, the maximum allowed ESL-spring force can be written as \({F}_{{{\max}}}\approx {\varepsilon }_{0}{E}_{0}^{2}Wd/2\), which is interestingly independent of the dielectric film. For example, if we choose \(W=4{{{\rm{\mu }}}}{{{\rm{m}}}}\) and \(d=100{{{\rm{nm}}}}\), we can obtain \({F}_{{{\max}}}\approx 8.9{{{\rm{\mu }}}}{{{\rm{N}}}}\), which is significantly larger than \(1{{{\rm{\mu }}}}{{{\rm{N}}}}\). The corresponding maximum allowed voltage for BaTiO3 is \({V}_{{{\max}}}={E}_{{{{\rm{cr}}}}}d\approx 7.1{{{\rm{V}}}}\), which has a normal value of \({\varepsilon }_{{{{\rm{r}}}}}=1,000\) as the dielectric film. The above analysis verifies that the ESL-spring is not only physically but also materially feasible. Furthermore, there will be a normal electrostatic force while applying a constant electric bias \({V}_{{{{\rm{b}}}}}\) to Electrodes 1 and 2. However, this would not generate enough additional friction to affect the operation of the ESL-resonator due to the extremely low friction coefficient of \( \sim 0.001\) in the SSL contact state (as detailed in Supplementary Note 8).

Using the ESL-spring as a resonator, we obtain the (circular) resonance frequency as follows:

$$\omega =\sqrt{\frac{k}{m}}=\sqrt{\frac{{\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}}{\rho L{L}_{{{{\rm{c}}}}}Hd}}{V}_{{{{\rm{b}}}}},$$

where \(m\), \(\rho \) and \(H\) denote the mass, mass density and thickness of the slider, respectively. By taking the maximum allowed voltage, \({V}_{{{{\rm{cr}}}}}\), before the dielectric breakdown and choosing, for instance, \(H=L/10\) and \(s=L/2\), we can obtain the maximum achievable resonance frequency \({\omega }_{{{\max}}}\) in the form \({\omega }_{{{\max}}}={v}_{0}{(d/{L}^{3})}^{1/2}\), where \({v}_{0}={E}_{0}\sqrt{20{\varepsilon }_{0}/\rho }\approx 628{{{{\rm{ms}}}}}^{-1}\). Therefore, the resonant frequency can be up to \( \sim25{{{\rm{MHz}}}}\) in the microscale (\(L=4{{{\rm{\mu }}}}{{{\rm{m}}}},d=100{{{\rm{nm}}}}\)) and \( \sim2.5{{{\rm{GHz}}}}\) in the nanoscale (\(L=40{{{\rm{nm}}}},d=1{{{\rm{nm}}}}\)). Since the tunable frequency range of ESL-spring resonator is much smaller than the natural frequency of the graphite flake (\(1.7{{{\rm{GHz}}}}\) for micro scale and \(170{{{\rm{GHz}}}}\) for nano scale), we ignored the deformation of graphite in the above analysis (see Supplementary Note 7).

The tunable range comparison between different types of resonators

As demonstrated above, one of the most remarkable features of the ESL-resonator is its widely tunable range of resonant frequencies10. For comparison, Fig. 2a summarises some representative tunable ranges reported in the literature and the tunable range of the ESL-resonator. To better comprehend the advantage of the tunability of the ESL-resonator, we also plot the relative tuning ranges of \(({\omega }_{{{\max}}}–{\omega }_{{{\min }}})/{\omega }_{{{\max}}}\) in Fig. 2b, where \({\omega }_{{{\max}}}\) and \({\omega }_{{{\min }}}\) denote the maximum and minimum tunable resonant frequencies, respectively. These comparison results demonstrate the superiority of the ESL-resonator in terms of tunability.

Fig. 2: Comparison between the ESL-resonator and conventional resonators.
figure 2

a Tunable resonant ranges for five piezoelectric resonators6,11,12,13, six electrostatic resonators18,19,20,21, five thermoelectric resonators3,14,15,16,17 and that of the ESL-resonator (for \(W=L=4{{\mu}}{{{\rm{m}}}},s=2{{\mu}}{{{\rm{m}}}},d=100{{{\rm{nm}}}},{\varepsilon }_{{{{\rm{r}}}}}=3.9,\rho =2.16{{{{\rm{gcm}}}}}^{-3},H=50{{{\rm{nm}}}}.\)) The solid lines of the ranges are given according to their respective theoretical or experimental relations, such as Eq. (3), or the literature3,6,11,12,13,14,15,16,17,18,19,20,21. b A comparison of the relative frequency tuning ranges of \(({\omega }_{{{\max}}}–{\omega }_{{{\min }}})/{\omega }_{{{\max}}}\). Details in the Supplementary Note 3.

The stability of ESL-resonators

In the setup of the ESL-resonator shown in Fig. 1a, c, there is no mechanical constraint to the motions of the slider in either the lateral direction (i.e., the \(y\)-axis) or in the \(x\)O\(y\) rotational plane. The introduction of such a mechanical constraint is undesirable, as it not only increases the structural complexity and cost but also induces additional and often significant friction. Therefore, a fundamental question to ask is whether the kinetics of such a free-constrained slider are stable.

Hereafter, we theoretically demonstrate the kinetic robustness of the ESL-resonator. Figure 3a illustrates the top view of the ESL-resonator; the slider moves to a deviating position characterised by a coordinate in the \((x,y,\theta )\)-phase space with \(y\,\ne\,0\) and/or \(\theta \,\ne\, 0\). To simplify the kinetic stability analysis, we introduce three dimensionless parameters, \(\xi =2x/{L}_{{{{\rm{c}}}}},\eta =y/W,\,{{{\rm{and}}}}\,\gamma =2\theta /\pi ,\) which are all constrained in the range of \((–1,1)\). The two capacitors denoted by \({C}_{1}\) and \({C}_{2}\), formed between the slider and Electrodes 1 and 2, respectively, can be expressed as \({C}_{1}={\varepsilon }_{0}{\varepsilon }_{r}{A}_{1}/d\) and \({C}_{2}={\varepsilon }_{0}{\varepsilon }_{r}{A}_{2}/d\), where \({A}_{1}(\xi ,\eta ,\gamma )\) and \({A}_{2}(\xi ,\eta ,\gamma )\) are the overlapping areas between the slider and Electrodes 1 and 2, respectively. We can express the equivalent potential energy, \(U\), of the system in the form of the first relationship in Eq. (1), yielding:

$$\begin{array}{c}U=g(\xi ,\eta ,\gamma )\frac{{C}_{0}}{8}{V}_{{{{\rm{b}}}}}^{2},\end{array}$$

where the function \(g(\xi ,\eta ,\gamma )\) represents the normalised energy that degenerates \(g(\xi ,0,0)={\xi }^{2}–1\) and is determined by \({A}_{1}\) and \({A}_{2}\). Mathematically, whenever the normalised energy \(g(\xi ,\eta ,\gamma )\) is at a minimum at \((\xi ,0,0)\), the kinetic system of the slider is stable at \((\xi ,0,0).\)

Fig. 3: Stability analysis of the ESL-resonator.
figure 3

a Illustration of the top view of the ESL-resonator. bd Normalised energy surfaces in the \(\{\eta ,\gamma \}\) region for \(\xi =0.0,0.4{{{\rm{and}}}}0.8\) in b, c, and d, respectively. For better visualisation of the results, we have plotted the contour lines of the surfaces in bd. Two typical size ratios of \(W/L=0.5\,{{{\rm{and}}}}\,s/L=0.2\) are used.

We use numerical methods to calculate \(g(\xi ,\eta ,\gamma )\) for different \(\xi ,\eta \), and \(\gamma \), all in the range of \((–1,1)\), as detailed in the Supplementary Note 4. Figure 3b–d show the normalised energy surfaces of \(g=g({\xi }_{0},\eta ,\gamma )\) for \({\xi }_{0}=0.0,0.4,0.8,\) respectively, and two typical size ratios \(W/L=0.5\) and \(s/L=0.2\). We can clearly see that the energy is always at its minimum at \(\eta =\gamma =0\), as the slider is non-rotationally moved along the \(x\)-direction. This means that the slider suffers restrains on the \(y\) and \(\theta \) directions while moving. Hence, the kinetics of such a free-constrained slider are stable. The above-characterised motion robustness is valid for other values of \({\xi }_{0},\) \(W/L\) and \(s/L\), when \(|{\xi }_{0}|\) is not too large, as explained in the Supplementary Note 4.

The application of ESL-resonators to ESL-oscillators

To apply the above-mentioned ESL-resonator as an oscillator, a major challenge is to apply a stimulating load without affecting the SSL state. To solve this problem, we propose another elementary setup, as illustrated in Fig. 4a, where the STATOR is formed by inserting a middle electrode (Electrode m) with length \(a\), width \(W\), and separation distance \(s/2\) between Electrodes 1 and 2. We shall further show that when choosing the pair of voltages \({V}_{1}\) and \({V}_{2}\) between Electrode m and Electrodes 1 and 2, respectively, in the alternative forms \({V}_{1}={V}_{0}\,\sin \,{\omega }_{{{{\rm{e}}}}}t+{V}_{{{{\rm{b}}}}}/2,{V}_{2}={V}_{0}\,\sin \,{\omega }_{{{{\rm{e}}}}}t–{V}_{{{{\rm{b}}}}}/2,\) the setup will become an oscillator stimulated by the harmonic voltage \({V}_{0}\,\sin \,{\omega }_{{{{\rm{e}}}}}t\) with a resonant frequency that is continuously tuned by changing the bias potential \({V}_{{{{\rm{b}}}}}\). The slider will form capacitors \({C}_{1}\), \({C}_{{{{\rm{m}}}}}\) and \({C}_{2}\) with Electrodes 1, m and 2, respectively, and the equivalent circuit diagram is shown in Fig. 4b. As detailed in the Supplementary Note 5, the dynamic equation of the slider can be expressed as \(m\ddot{x}+c\dot{x}+kx={F}_{{{{\rm{ex}}}}}\,\sin \,{\omega }_{{{{\rm{e}}}}}t,\) where \(m\) is the mass of the slider, \(c\) is the mechanical damping (much larger than ohmic loss and dielectric loss; see Supplementary Note 6), and \(k\) is the spring coefficient, which has exactly the same form as in Eq. (2),

$$\begin{array}{c}{F}_{{{{\rm{ex}}}}}=\frac{W{\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}}{{L}_{{{{\rm{c}}}}}d}a{V}_{0}{V}_{{{{\rm{b}}}}},\end{array}$$

denoting the excitation force with \({L}_{{{{\rm{c}}}}}=L–s\).

Fig. 4: Schematic representation of the ESL-oscillator.
figure 4

a Illustration of an elementary ESL-oscillator, where the two alternative voltages \({V}_{1}\) and \({V}_{2}\) are applied to the left and right electrodes, respectively. b Equivalent circuit diagram of the ESL-oscillator. c Relationship of the oscillating amplitude, \({A}_{{{{\rm{os}}}}},\) versus \({V}_{0}\). In this case, the slider displacement must be less than \(({L}_{{{{\rm{c}}}}}–a)/2=1{{{\rm{\mu }}}}{{{\rm{m}}}}\), which implies an excitation voltage of \({V}_{0}\le 3{{{\rm{mV}}}}\). d, e Relationship of the maximum velocity in the oscillation, \({v}_{{{{\rm{os}}}}}\), and Q-factor versus the tunable electric bias, \({V}_{{{{\rm{b}}}}}\), for \({V}_{0}=2{{{\rm{mV}}}}\) in d. The parameters used in the calculations are: \({\varepsilon }_{{{{\rm{r}}}}}=3.9,W=L=4{{{\rm{\mu }}}}{{{\rm{m}}}},s=1{{{\rm{\mu }}}}{{{\rm{m}}}},a=1{{{\rm{\mu }}}}{{{\rm{m}}}},H=200{{{\rm{nm}}}},d=20{{{\rm{nm}}}},\rho =2.16{{{{\rm{gcm}}}}}^{-3},\) and the damping is \(c=1.2\times {10}^{-11}{{{{\rm{kgs}}}}}^{-1}\)33,34. Details are provided in the Supplementary Note 5.

Based on a physically meaningful understanding of Eq. (5), we can reform it to \({F}_{{{{\rm{ex}}}}}={Q}_{{{{\rm{m}}}}}{E}_{{{{\rm{ho}}}}}\), where \({Q}_{{{{\rm{m}}}}}={C}_{{{{\rm{m}}}}}{V}_{0}\) can represent an equivalent charge of the middle capacitor \({C}_{{{{\rm{m}}}}}={\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}Wa/d\) excited by the voltage \({V}_{0}\), and \({E}_{{{{\rm{ho}}}}}={V}_{{{{\rm{b}}}}}/{L}_{{{{\rm{c}}}}}\) can correspond to an equivalent horizontal electric field formed by the bias voltage, \({V}_{{{{\rm{b}}}}}.\). This understanding endows the physical picture, where the force \({F}_{{{{\rm{ex}}}}}\) is electrostatically applied to the charge \({Q}_{{{{\rm{m}}}}}\) in the electric field \({E}_{{{{\rm{ho}}}}}\).

By solving the dynamic equation, we can obtain the resonance amplitude, \({A}_{{{{\rm{os}}}}}\); velocity, \({v}_{{{{\rm{os}}}}}\); and Q-factor, \(Q,\) as:

$$\begin{array}{c}{A}_{{{{\rm{os}}}}}=\frac{{F}_{{{{\rm{ex}}}}}}{{\omega }_{0}c}=\frac{1}{{c}_{1}}\sqrt{\frac{{\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}\rho LH}{d{L}_{{{{\rm{c}}}}}}}a{V}_{0},\\ {v}_{{{{\rm{os}}}}}={\omega }_{0}{A}_{{{{\rm{os}}}}}=\frac{{\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}}{{c}_{1}d{L}_{{{{\rm{c}}}}}}a{V}_{0}{V}_{{{{\rm{b}}}}},\\ Q=\frac{{\omega }_{0}}{c/m}=\frac{1}{{c}_{1}}\sqrt{\frac{{\varepsilon }_{0}{\varepsilon }_{{{{\rm{r}}}}}\rho LH}{d{L}_{{{{\rm{c}}}}}}}{V}_{{{{\rm{b}}}}},\end{array}$$

where \({c}_{1}=c/W\) represents the damping per unit width. To intuitively display the performance of the ESL-oscillator, we plot the above performance parameters for \({\varepsilon }_{{{{\rm{r}}}}}=3.9,W=L=4{{{\rm{\mu }}}}{{{\rm{m}}}},s=1{{{\rm{\mu }}}}{{{\rm{m}}}},a=1{{{\rm{\mu }}}}{{{\rm{m}}}},H=200{{{\rm{nm}}}},d=20{{{\rm{nm}}}},\rho =2.16{{{{\rm{gcm}}}}}^{-3},\;{{{\rm{and}}}}c=1.2\times {10}^{–11}{{{{\rm{kgs}}}}}^{-1}\)33,34, as shown in Fig. 4c–e; the damping coefficient, \(c\), is estimated in the Supplementary Note 5.

The spring coefficient, \(k\), is independent of \({V}_{0}\), and the electric excitation force, \({F}_{{{{\rm{ex}}}}}\), depends on the product \(a{V}_{0}\). However, different choices of \(a\) affect the maximum allowed oscillation amplitude, \(|x| \, < \,({L}_{{{{\rm{c}}}}}–a)/2\). The Q-factor increases with decreasing \(d\) in the form \( \sim {d}^{-1/2}\), which implies a remarkable advantage of using nanometre-thick dielectric films to separate the SLIDER and underlying electrodes.

The advantages of ESL-resonators and ESL-oscillators

We point out two additional remarkable advantages of ESL-resonators and ESL-oscillators. First, MEMS resonators constitute a type of mechanical resonator that has a compact size, high sensitivity, high resolution, and low cost [1,2]. Because most MEMS resonators are based on cantilevers, the allowed maximum relative displacement of the cantilevers is relatively small, typically less than \(1/10\) the length of the cantilever owing to the pull-in phenomenon35,36,37. This largely restricts the scope of MEMS applications in some important devices, such as energy-harvesting devices6 and micro-gyroscopes. In comparison, the ESL-oscillator allows significantly larger relative sliding displacements. In fact, the maximum allowed displacement, \({x}_{{{\max}}}\), is equal to \(({L}_{c}-a)/2\), which is of the same order as the ESL-oscillator size. This unique property promises applications of ESL-oscillators with better performance capabilities than those of MEMS in energy-harvesting devices6 and micro-gyroscopes (as detailed in Supplementary Note 9).

Second, the restoring force of the ESL-spring is strictly linear, which results in good harmonic vibration dynamic characteristics of ESL-oscillators at any voltage. In contrast, the additional restoring forces of the conventual tunable resonators and oscillators are often nonlinear6,18, which leads to non-harmonic vibrations and often affects the performance of the resonators and oscillators.


In summary, using structural superlubricity, a recently introduced technology, we proposed one type of linear spring, called an electro-superlubric spring (ESL-spring), which possesses a unique feature of continuously tunable spring stiffness by alteration of the applied bias voltage. This suggests the future development of ESL-resonators and ESL-oscillators with continuously tunable resonant frequencies in the range from zero to several MHz for microscale ESL-resonators and from zero to several GHz for nanoscale ESL-resonators. Moreover, the relative displacement amplitudes of ESL-resonators can largely exceed those of mechanical ones. Therefore, the commercialisation of ESL-resonators can be expected to be realised in the near future, which will lead to their broad application. However, the experimental results are not yet accomplished; this remains to be addressed in future studies.


Finite element simulations of the ESL-spring

In order to verify the theoretical model of ESL-spring in Supplementary Note 1 and take into account the influence of edge effect, we develop a finite element model to simulate the structure in Fig. 1. At the boundary, a large area of air is set outside, as shown in Supplementary Fig. 2a, while the zero-charge boundary condition is set on the outside boundary of air (See Supplementary Eq. (S13)). Simultaneously, all electrodes are equipotential; and the voltage of Electrode 1 is set to 0 V, and the voltage of Electrode 2 is set to 40 V, which is equivalent to applying a constant voltage \({V}_{{{{\rm{b}}}}}=40{{{\rm{V}}}}\). For the Electrode 0, we give it a zero-total charge condition, as Supplementary Equation (S14) shows. The whole region is solved by the Maxwell equation in electrostatic field, as details in Supplementary Eq. (S15).

Calculating the overlap area of Slider and Electrodes

Here, we need to use a numerical method to calculate \({A}_{1}(\xi ,\eta ,\gamma )\) and \({A}_{2}(\xi ,\eta ,\gamma )\). Firstly, we divide the area of Electrode 0 into \(N\times N\) (here \(N=500\)) elements uniformly, as illustrated in Supplementary Fig. 3. The position of each element is discretized as a point (\(\frac{{x}_{0}^{(i,j)}L}{2},\frac{{y}_{0}^{(i,j)}W}{2}\)), in which \({x}_{0}^{(i,j)}\) and \({y}_{0}^{(i,j)}\) is defined as Supplementary Eq. (S16). Then, when Electrode 0 produces the displacement \((\xi ,\eta ,\gamma )\), we can calculate \((i,j)\)-th element’s position in the global coordinate system \(x\)-O-\(y\) using Supplementary Eq. (S17). Next, we just need to judge each element whether it is in the area where Electrode1&2 are located, to determine whether it belongs to the overlap area \({A}_{1}\) and \({A}_{2}\) by Supplementary Eq. (S18). Then we can use Supplementary Eq. (S19) to obtain the overlap area.

Data availability

The data that support the plots within this paper and other finding of this study are available from the corresponding authors upon reasonable request.


  1. 1.

    Tabrizian, R., Raiszadeh, M. & Ayazi, F. Effect of phonon interactions on limiting the f.Q product of micromechanical resonators. TRANSDUCERS 2009–2009 International Solid-State Sensors, Actuators and Microsystems Conference. p. 2131–2134 (IEEE, 2009).

  2. 2.

    van Beek, J. T. M. & Puers, R. A review of MEMS oscillators for frequency reference and timing applications. J. Micromech. Microeng. (2012).

  3. 3.

    Svilicic, B., Mastropaolo, E. & Cheung, R. Widely tunable MEMS ring resonator with electrothermal actuation and piezoelectric sensing for filtering applications. Sens. Actuat A Phys. 226, 149–153 (2015).

    CAS  Article  Google Scholar 

  4. 4.

    Sonmezoglu, S., Alper, S. E. & Akin, T. An automatically mode-matched MEMS gyroscope with wide and tunable bandwidth. J. Microelectromech. Syst. 23, 284–297 (2014).

    Article  Google Scholar 

  5. 5.

    Roessig, T. A., Howe, R. T., Pisano, A. P. & Smith, J. H. Surface-micromachined Resonant Accelerometer (IEEE, 1997).

  6. 6.

    Peters, C., Maurath, D., Schock, W., Mezger, F. & Manoli, Y. A closed-loop wide-range tunable mechanical resonator for energy harvesting systems. J. Micromech. Microeng. (2009).

  7. 7.

    Harne, R. L. & Wang, K. W. A review of the recent research on vibration energy harvesting via bistable systems. Smart Mater. Struct. (2013).

  8. 8.

    Amri, M. et al. Stiffness control of a nonlinear mechanical folded beam for wideband vibration energy harvesters. TM-Tech. Mess. 85, 553–564 (2018).

    Article  Google Scholar 

  9. 9.

    Entesari, K. & Rebeiz, G. M. A differential 4-bit 6.5-10-GHz RF MEMS tunable filter. IEEE Trans. Microw. Theory Tech. 53, 1103–1110 (2005).

    Article  Google Scholar 

  10. 10.

    Zhang, W.-M., Hu, K.-M., Peng, Z.-K. & Meng, G. Tunable micro- and nanomechanical resonators. Sensors 15, 26478–26566 (2015).

    Article  Google Scholar 

  11. 11.

    Piazza, G., Abdolvand, R., Ho, G. K. & Ayazi, F. Voltage-tunable piezoelectrically-transduced single-crystal silicon micromechanical resonators. Sens. Actuat. A Phys. 111, 71–78 (2004).

    CAS  Article  Google Scholar 

  12. 12.

    Roundy, S. & Zhang, Y. In Smart Structures, Devices, and Systems II, Pt 1 and 2 Vol. 5649 Proceedings of the Society of Photo-Optical Instrumentation Engineers (Spie) (ed. S. F. AlSarawi) 373–384 (2005).

  13. 13.

    Olivares, J. et al. Mems, Moems, and Micromachining Ii Vol. 6186 Proceedings of SPIE (eds H. Urey & A. ElFatatry) (Society of Photo Optical, 2006).

  14. 14.

    Mastropaolo, E., Wood, G. S., Gual, I., Parmiter, P. & Cheung, R. Electrothermally actuated silicon carbide tunable mems resonators. J. Microelectromech. Syst. 21, 811–821 (2012).

    CAS  Article  Google Scholar 

  15. 15.

    Svilicic, B., Mastropaolo, E., Flynn, B. & Cheung, R. Electrothermally actuated and piezoelectrically sensed silicon carbide tunable MEMS resonator. IEEE Electron Dev. Lett. 33, 278–280 (2012).

    CAS  Article  Google Scholar 

  16. 16.

    Svilicic, B., Mastropaolo, E., Zhang, R. & Cheung, R. Tunable MEMS cantilever resonators electrothermally actuated and piezoelectrically sensed. Microelectron. Eng. 145, 38–42 (2015).

    CAS  Article  Google Scholar 

  17. 17.

    Kazmi, S. N. R. et al. Tunable nanoelectromechanical resonator for logic computations. Nanoscale 9, 3449–3457 (2017).

    CAS  Article  Google Scholar 

  18. 18.

    Adams, S. G. et al. Capacitance based tunable resonators. J. Micromech. Microeng. 8, 15–23 (1998).

    Article  Google Scholar 

  19. 19.

    Eriksson, A. et al. Direct transmission detection of tunable mechanical resonance in an individual carbon nanofiber relay. Nano Lett. 8, 1224–1228 (2008).

    CAS  Article  Google Scholar 

  20. 20.

    Morgan, B. & Ghodssi, R. Vertically-shaped tunable MEMS resonators. J. Microelectromech. Syst. 17, 85–92 (2008).

    Article  Google Scholar 

  21. 21.

    Chen, C. et al. Graphene mechanical oscillators with tunable frequency. Nat. Nanotechnol. 8, 923–927 (2013).

    CAS  Article  Google Scholar 

  22. 22.

    Liu, Z. et al. Observation of microscale superlubricity in graphite. Phys. Rev. Lett. (2012).

  23. 23.

    Hod, O., Meyer, E., Zheng, Q. & Urbakh, M. Structural superlubricity and ultralow friction across the length scales. Nature 563, 485–492 (2018).

    CAS  Article  Google Scholar 

  24. 24.

    Zheng, Q. Superlubricity: the world of ‘0’ friction. Sci. Technol. Rev. 34, 12–26 (2016).

    Google Scholar 

  25. 25.

    Peng, D. 100 km wear-free sliding achieved by microscale superlubric graphite/DLC heterojunctions under ambient conditions. Natl. Sci. Rev. (2021).

    Article  Google Scholar 

  26. 26.

    Shigeki Kawai, A. B. et al. Superlubricity of graphene nanoribbons on gold surfaces. Science 351, 957–961 (2016).

    Article  Google Scholar 

  27. 27.

    Song, Y. et al. Robust microscale superlubricity in graphite/hexagonal boron nitride layered heterojunctions. Nat. Mater. 17, 894–899 (2018).

    CAS  Article  Google Scholar 

  28. 28.

    Peng, D. et al. Load-induced dynamical transitions at graphene interfaces. Proc. Natl Acad. Sci. USA 117, 12618–12623 (2020).

    CAS  Article  Google Scholar 

  29. 29.

    Zheng, Q. et al. Self-retracting motion of graphite microflakes. Phys. Rev. Lett. (2008).

  30. 30.

    Yang, J. et al. Observation of high-speed microscale superlubricity in graphite. Phys. Rev. Lett. (2013).

  31. 31.

    Wang, W. et al. Measurement of the cleavage energy of graphite. Nat. Commun. 6, 7 (2015).

    Article  Google Scholar 

  32. 32.

    Neusel, C. & Schneider, G. A. Size-dependence of the dielectric breakdown strength from nano- to millimeter scale. J. Mech. Phys. Solids 63, 201–213 (2014).

    Article  Google Scholar 

  33. 33.

    Wang, K., Qu, C., Wang, J., Quan, B. & Zheng, Q. Characterization of a microscale superlubric graphite interface. Phys. Rev. Lett. (2020).

  34. 34.

    Qu, C. et al. Origin of friction in superlubric graphite contacts. Phys. Rev. Lett. (2020).

  35. 35.

    Dequesnes, M., Rotkin, S. V. & Aluru, N. R. Calculation of pull-in voltages for carbon-nanotube-based nanoelectromechanical switches. Nanotechnology 13, 120–131 (2002).

    Article  Google Scholar 

  36. 36.

    Nayfeh, A. H., Younis, M. I. & Abdel-Rahman, E. M. Dynamic pull-in phenomenon in MEMS resonators. Nonlinear Dyn. 48, 153–163 (2007).

    Article  Google Scholar 

  37. 37.

    Zhang, W.-M., Yan, H., Peng, Z.-K. & Meng, G. Electrostatic pull-in instability in MEMS/NEMS: a review Sens. Actuat. A-Phys. 214, 187–218 (2014).

    CAS  Article  Google Scholar 

Download references


Q.Z. wishes to acknowledge the financial support of the National Natural Science Foundation of China (Nos. 11572173, 11890671 and 51961145304), the National Key Basic Research Program of China (Grant No. 2013CB934200), the Cyrus Tang Foundation (Grant No. 202003), the Beijing Municipal Science & Technology Commission (Program No. Z151100003315008), the Tsinghua University Initiative Scientific Research (Program Nos. 2014z01007 and 2012Z01015) and the State Key Laboratory of Tribology Tsinghua University Initiative Scientific Research (Program No. SKLT2019D02). We would like to thank Editage ( for English language editing.

Author information




Z.W. and X. H. designed and analysed the resonators. Z.W., X.H. and X.X. discussed and proposed the idea. Q.Z. supervised the project. Z.W., X.H. and Q.Z. wrote the paper.

Corresponding authors

Correspondence to Zhanghui Wu, Xuanyu Huang or Quanshui Zheng.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Communications Materials thanks Ghader Rezazadeh and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Aldo Isidori. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wu, Z., Huang, X., Xiang, X. et al. Electro-superlubric springs for continuously tunable resonators and oscillators. Commun Mater 2, 104 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing