Finite Element-based Nonlinear Dynamic Optimization of Nanomechanical Resonators
Abstract
Nonlinear dynamic simulations of mechanical resonators have been facilitated by the advent of computational techniques that generate nonlinear reduced order models (ROMs) using the finite element (FE) method. However, designing devices with specific nonlinear characteristics remains inefficient since it requires manual adjustment of the design parameters and can result in suboptimal designs. Here, we integrate an FE-based nonlinear ROM technique with a derivative-free optimization algorithm to enable the design of nonlinear mechanical resonators. The resulting methodology is used to optimize the support design of high-stress nanomechanical \ceSi3N4 string resonators, in the presence of conflicting objectives such as simultaneous enhancement of -factor and nonlinear Duffing constant. To that end, we generate Pareto frontiers that highlight the trade-offs between optimization objectives and validate the results both numerically and experimentally. To further demonstrate the capability of multi-objective optimization for practical design challenges, we simultaneously optimize the design of nanoresonators for three key figure-of-merits in resonant sensing: power consumption, sensitivity and response time. The presented methodology can facilitate and accelerate designing (nano)mechanical resonators with optimized performance for a wide variety of applications.
Introduction
Design of mechanical structures that move or vibrate in a predictable and desirable manner is a central challenge in many engineering disciplines. This task becomes more complicated when these structures experience large-amplitude vibrations, since linear analysis methods fail and nonlinear effects need to be accounted for. This is particularly important at the nanoscale, where forces on the order of only a few pN can already yield a wealth of nonlinear dynamic phenomena worth exploiting 1, 2, 3, 4, 5.
Although design optimization of micro and nanomechanical resonators in the linear regime is well-established 6, the use of design optimization for engineering nonlinear resonances has received less attention 7. This is because designers tend to avoid the nonlinear regime, and optimizing structures’ nonlinear dynamics is more complex, which requires extensive computational resources. As a result, available literature on nonlinear dynamic optimization is limited, although some recent advances have been made that combine analytical methods with gradient-based shape optimization, to optimize nonlinearities in micro beams 8, 9. For nonlinear modeling of more complex structures, several approaches have been developed based on nonlinear reduced order modeling (ROM) of finite element (FE) simulations 10, 11, 12. A particularly attractive class known as STEP (STiffness Evaluation Procedure) 13 can determine nonlinear coefficients of an arbitrary mechanical structure and can be implemented in virtually any commercial finite element method (FEM) package. This, for instance, has been recently shown by using COMSOL to model the nonlinear dynamics of high-stress \ceSi3N4 string 14 as well as graphene nanoresonators 15. Since the number of degrees of freedom in the ROM is much smaller than that in the full FE model, the nonlinear dynamics of the structure can be simulated much more rapidly using numerical continuation packages 16.
In this work, we present a route for nonlinear dynamic optimization that is based on an FE-based ROM. The methodology, which is a combination of Particle Swarm Optimization (PSO) with STEP 13 (OPTSTEP), has several beneficial features. First of all, because it uses a derivative-free optimization routine for approaching the optimal design, it can be implemented and combined with FEM packages that are not able to obtain gradients easily. Secondly, the ROM parameters generated in OPTSTEP can facilitate explicitly expressing the optimization goals. Finally, as will be shown, the developed procedure allows using multiple objective functions to approximate a Pareto front, which can help designers in decision-making processes when having to balance performance trade-offs among different objectives. Considering the outstanding performance as ultrasensitive mechanical detectors and the mature fabrication procedure 17, 18, we select high-stress \ceSi3N4 for the experimental validation of our methodology.
The manuscript is structured as follows. We first introduce and describe the general OPTSTEP methodology. Then we demonstrate the method on the specific challenge of the optimization of the support structure for a high-stress \ceSi3N4 nano string, while taking the maximization of its -factor and nonlinear Duffing constant as examples of linear and nonlinear objectives. By comparing the PSO results to the and values that result from a brute-force simulation of a large number of designs that span the design space, we validate that OPTSTEP finds the optimum designs much faster with the same computational resources.
Subsequently, we turn to the problem of dealing with multiple objective functions and focus on simultaneously maximizing both and , demonstrated by a Pareto front. For validation, the results are compared to experimental measurements of fabricated devices. We conclude by demonstrating the potential of OPTSTEP for optimizing the performance of resonant sensors by using more complex objective functions that are relevant for engineering their response time, sensitivity, and power consumption.
OPTSTEP methodology
An overview of the OPTSTEP method is schematically shown in Fig. 1. In the current work, we use it for engineering a parameterized geometry. We use nanomechanical string resonators with compliant supports, which is shown in Fig. 1a, to demonstrate the methodology. We keep the length and width of the central string constant, while varying the width , length and angle of the supports, as well as the thickness of the device. It is noted that the OPTSTEP methodology might be used with a larger number of parameters, or even might be extended towards shape or topology optimization of nonlinear dynamic structures. However, such extension is out of the scope of the current work.
For a certain set of geometrical parameters, a ROM for the parameterized structure is generated using the STEP method 13, which we implemented with shell elements in COMSOL 14. Besides geometric parameters and boundary conditions (see Fig. 1a), the COMSOL simulation contains material parameters (see Methods), and the initial pre-stress distribution is calculated using a static analysis 14. We conduct this static analysis assuming the material is isotropic and pre-stressed (). We then calculate the stress redistribution during the sacrificial layer underetching process, whereby the high-stress \ceSi3N4 layer releases from the silicon substrate. Note that in the present study we only consider , such that the central string is always in tension (in contrast to Ref. 14). After the static analysis, an eigenfrequency analysis is performed to obtain the out-of-plane eigenmodes (see Fig. 1b). These eigenmodes, together with the redistributed stress field obtained from the static analysis, are then used to determine the -factor, resonance frequency , and effective mass , following the procedure outlined in Ref. 19.
As indicated in Fig. 1b the STEP method generates a set of coupled nonlinear differential equations13, 14, 15, where the effective nonlinear elastic force acting on the th mode is given by the function that depends on the quadratic , cubic coupling coefficients, and the generalized coordinates . describes the instantaneous contribution of the corresponding mode shapes to the deflection of the structure. Thus, the finite element model with several thousand or even millions of degrees of freedom (DOFs) is reduced to a condensed ROM, that can usually describe the nonlinear dynamics to a good approximation with less than ten degrees of freedom. We can visualize the resulting frequency response curves for different harmonic drive levels by numerical continuation 16, as shown in Fig. 1c.
The resulting ROM parameters, including effective mass , -factor, linear stiffness and nonlinear stiffness terms , , are passed to the PSO optimizer (see Fig. 1d). The algorithm randomly generates many different initial designs by varying the geometric parameters, as shown in Fig. 1a. For each of these designs, known as a ’particle’ in PSO, a ROM is generated by STEP and the corresponding objective functions are computed accordingly and passed to the optimizer. The optimizer then generates a next generation of particles based on the designs from the current generation, the objective functions, and the constraints, with the aim of improving their design parameters to optimize the objectives (see Supplementary Note 1). The optimization loop will iterate until it reaches the predefined maximum generation. If multiple objective functions are selected to be optimized, there is an additional step that selects the nondominated particles according to Pareto dominance 20. Because each particle is evaluated independently, PSO enables efficient parallel computing to evaluate all particles in one generation on a high-performance computing cluster.
OPTSTEP implementation and validation
Single objective optimization with OPTSTEP
We implement the presented OPTSTEP methodology to optimize the support geometry of the string resonator shown in Fig. 1a. The motion of the fundamental mode of the resonator can be described with the following nonlinear equation of motion:
(1) |
where is the displacement at the string center, is the resonance frequency, is the -factor, is the mass-normalized Duffing constant, and is the mass-normalized harmonic drive force. We present results of the OPTSTEP methodology for two optimization objectives, respectively: maximizing the -factor (shown in Fig. 2a,c,d) or maximizing the mass-normalized Duffing constant (shown in Fig. 2b,e,f) of the fundamental mode. As design parameters, we use the support parameters (, , and in Fig. 1a). The PSO algorithm can freely initialize and vary these variables between preset constraints , , , and .
We initialize the PSO algorithm with 10 randomly generated particles, as indicated by the blue circles at the first generation in Fig. 2a-b. The and values of the best performing particle per generation are highlighted by the red line, which converges towards an optimum. Simulated response curves at different drive levels of the initial design (median performance of the initialized particles) and the optimized design are shown in Fig. 2c, d for and Fig. 2e, f for . It is obvious that the resonance peaks become narrower from Fig. 2c to Fig. 2d, indicative of an increase in -factor. From the backbone curves shown in Fig. 2e, f, we see that the resonance frequency of the optimized device shifts more at the same vibration amplitude, which suggests a larger, optimized value of .
Numerical validation
In order to validate the PSO results, we compare them to a brute-force parametric study where we simulate a large number of designs that span the full design parameter space, and plot the resulting values of and in the contour plots in Fig. 2g, h. Each of these subfigures consists of 16 small contour plots, each of which has a different combination of and , while along the axes the parameters and are varied. The red-colored regions in the plots contain the optimal values of and , which are indicated by a triangle and a star. In Supplementary Table S1, we compare the optimized design parameters from the OPTSTEP method to the best devices from the parametric study. The close agreement between both approaches provides evidence that the OPTSTEP method is able to optimize both linear () and nonlinear () parameters of the ROM. The results in Fig. 2a are obtained in 30 minutes using a high performance computing cluster, while the parametric study in Fig. 2g takes over 325 hours on the same cluster with the same amount of nodes. This illustrates the advantage in computation time that can be realized with OPTSTEP, although it is noted that these times strongly depends on the resolution of the parameter grid and other simulation parameters.
Experimental characterization
To compare the OPTSTEP method to experimental results, we also perform an experimental parametric study on 15 string resonators with varying support design parameters. For this we fabricated a set of devices with and , while keeping and fixed. Fig. 3a shows a Scanning Electron Microscope (SEM) image of an array of nanomechanical resonators with varying support designs made of high-stress \ceSi3N4 (see “Methods” for more details). To characterize the nonlinear dynamics of the devices, as shown in Fig. 3b, we fix the chip to a piezo actuator that drives the resonator by an out-of-plane harmonic base actuation in the out-of-plane direction. We use a Zurich Instruments HF2LI lock-in amplifier, connected to an MSA400 Polytec Laser Doppler Vibrometer, to measure the out-of-plane velocity at the center of the string resonator as a function of driving frequency (see Fig. 3c). We use a velocity decoder with a calibration factor of 200 mm/s/V. We perform all measurements in a vacuum chamber with a pressure below at room temperature.
Fig. 3c shows the frequency response at the center of the string at various drive levels for a device with , , and . We estimate the linear resonator parameters of all devices by fitting the measured frequency response curves at various drive levels with the following harmonic oscillator function14:
(2) |
where is the measured amplitude, is the value of when driving at the natural resonance frequency , and is the drive frequency. To determine the nonlinear stiffness, we measure the resonator’s frequency response at increasing drive levels, construct the backbone curve, and use the relation between the peak amplitude and the peak frequency to estimate the mass-normalized Duffing constant from21, 22:
(3) |
To compensate for small drifts in during the experiments, before fitting with Eq. (3), we shift and align the frequency response curves to match their values14.
In Fig. 3d-f, we compare the dynamical properties between FE-based ROMs (dots) and measurements on 15 string resonators (diamonds) as a function of and . It is evident that the fundamental resonance frequency , -factor, and the mass-normalized Duffing constant of the fabricated devices, are all well predicted by FE-based ROMs. It can also be seen that for short support lengths the device performance is similar, whereas increasing allows tuning , and as we studied in more detail earlier19, 14. In the next section we will compare these experimental results to multi-objective optimization as further validation of OPTSTEP.
Multi-objective optimization with OPTSTEP
For actual device design there are often multiple performance specifications that need to be met. It might sometimes be possible to condense these performance specifications into a single figure of merit, like the product for nanomechanical resonators. However, to make the best design decisions, it is preferred that the optimizer works with two (or more) objective functions like enhancing and , simultaneously. To enable this, we implement OPTSTEP with a multi-objective particle swarm optimization (MOPSO), which is an extension of single-objective PSO. After multi-objective optimization, the nondominated particles in the swarm are used to determine an approximation of the Pareto front, which is the set of designs for which improving one of the objectives will always lead to a deterioration of the other objective(s). By performing MOPSO, we aim at finding the Pareto front in the design space for multiple objectives, that represents the boundary on which all optimized designs reside for the chosen variables. As the red dots show in Fig. 1d illustrate, the Pareto front represents the boundary between feasible and unfeasible combinations of objectives and thus allows the designer to make the best trade-off among different objectives.
To demonstrate that multi-objective optimization can be combined with OPTSTEP, we use it to simultaneously maximize and . Devices with high quality factor and nonlinear stiffness can be of interest in cases where we are looking for designs that can drive a string into the nonlinear regime with a minimum driving force and power consumption.
The resulting Pareto fronts are shown in Fig. 4a. Since we are also interested in the effect of the constraints on the optimum solutions, we include Pareto fronts with: no constraint (purple), a thickness constraint of (grey), and with thickness and support width constraint (multi-coloured). These 3 Pareto fronts show that there is a clear trade-off between and , with higher -factor leading to lower nonlinearity . The experimental devices share the same constraints ( and ) as the multi-colored Pareto and are plotted as the hollow diamonds with error bars in Fig. 4a (see Supplementary Table 2). We observe that all experimental points reside in the region on the left hand side of the Pareto front, confirming the area enclosed by the Pareto front indeed captures the feasible devices, and experimentally strengthening the confidence in the OPTSTEP approach for multi-objective designs. The color of the points links the points in the graph in Fig. 4a to the corresponding design parameters in Fig. 4b. In Fig. 4b the schematic support geometries are shown as insets for both maximum (dark blue) and maximum (dark red). We choose some of the fabricated devices close to the Pareto front to show typical measured frequency response curves and microscopic images in Fig. 4c-f, which correspond to the star, triangle, circle and square data markers in Fig. 4a and b. Together with the microscopic images, it is apparent that with minor alterations in the support region, the response of the string resonators can be largely tuned. To further explore the effect of other design parameters numerically, we release the constraint on , keeping only constrained, and conduct MOPSO (see the grey Pareto front). We can see from the comparison between the grey and multicolored fronts that the performance gain from changing is not very large. In contrast, if we further relax the constraint on , which shares the same design space in Fig. 2g-h, we obtain the purple Pareto front. The thinner pushes the Pareto front to have much higher . The long plateau at fixed is mainly attributed to the increase in that results from the dependence of the intrinsic quality factor on (see Methods).
Besides validating the MOPSO approach by comparing with experimental data, we also use the data from the parametric study in Fig. 2 to extract and generate reference Pareto fronts that are shown as black solid, dotted, and dashed lines in Fig. 4a (see Supplementary Note 2), with constraints that match those from the MOPSO optimization.
Discussion
The OPTSTEP methodology that is presented in this work enables the optimization of the nonlinear dynamic properties of resonant structures using standard FEM software, since it is based on the STEP and uses a derivative-free optimization method. We note that although derivative-free techniques like PSO are able to efficiently find near-optimal values of design parameters, optimality guarantees can typically not be given, and the techniques are therefore also called metaheuristic optimization techniques. Here, in order to validate the OPTSTEP methodology numerically and experimentally, we have focused on and maximization of the fundamental mode of a string resonator by geometric support design. After having established the methodology, it is now of interest to apply it to explore performance parameters that are more relevant to applications. For example, as shown in Fig. 5, our methodology can directly be extended to optimize the power consumption , sensitivity and response time of resonant sensors 23, 24, since these figure-of-merits can be directly expressed in terms of , , and (see Supplementary Note 3). In Fig. 5, 1000 nondominated particles are found by OPTSTEP to form a 3D surface that approaches the Pareto frontier with the objective of minimizing , and simultaneously. The particles have the same design constraints as in the example in Fig.2 and the purple Pareto front in Fig. 4a, which are and and membrane thickness . The competing design trade-offs between these three objective functions are obtained from OPTSTEP, and are visualized in Fig. 5 by showing five typical designs near the Pareto frontier. As demonstrated by the designs at the upper right corner of the Pareto frontier, we can conclude that the devices with shorter response time are more likely to have thicker supports, which lead to a higher resonance frequency combined with a low , thus resulting in a smaller ratio. At the same time, these thicker supports also contribute to a larger onset of nonlinearity 14, so the resonators are able to work at much larger amplitudes in the linear regime, which provides a better sensitivity . However, the larger and will require more energy to sustain the oscillation at resonance that causes higher power consumption . In contrast, the devices with higher sensitivity , which are shown at the lower left corner in Fig. 5, are equipped with more slender supports. With only a slight increase of support angle from 0, the low torsional stiffness of supports is maintained while the stress in the central string can be significantly increased 19, leading to a higher , which can be confirmed by Fig. 2g. Consequently, when aiming at designing a resonant sensor with relatively low power consumption , high sensitivity and short response time with compliant supports, a pair of slender and slightly angled supports, together with a medium thickness of \ceSi3N4 layer is generally favored.
In other cases, like approaching the quantum regime with a nonlinear nanomechanical resonator 25, it is beneficial to maximize and simultaneously. The OPTSTEP methodology can also be used for more complex design problems that involve multiple modes 8, 26, 5, 14, for avoiding or taking advantage of mode coupling, for instance by optimizing nonlinear coupling coefficients ( and in Fig. 1b) and resonance frequency ratios. Since OPTSTEP generates the ROM parameters at each generation, it is particularly suited for dealing with cases where the device specifications can be expressed in terms of these parameters. Interesting challenges include increasing frequency stability by coherent energy transfer 27, 28, signal amplification 29 and stochastic sensing 4, 30. Moreover, intriguing paths for further research involve inclusion of nonlinear damping or extension to full topology optimization 6. Also the use of alternative optimization strategies, like binary particle swarm optimization (BPSO) 31, that could generate radically new geometries, is an interesting direction.
Conclusions
To sum up, we presented a methodology (OPTSTEP) for optimizing the nonlinear dynamics of mechanical structures by combining an FE-based ROM method with a derivative-free optimization technique (PSO). We demonstrated and validated the methodology by optimizing the support design of high-stress \ceSi3N4 nanomechanical resonators. The method was verified numerically by comparing its results to a brute-force parametric study, for both single- and multi-objective optimization. Experimental data on the -factor and Duffing nonlinearity were in correspondence with the OPTSTEP results. The capability of the method was also demonstrated by multi-objective optimization of the support for the nanomechanical resonator, targeting improvements in power consumption, sensitivity and response time in resonant sensing.
We thus concluded that the method can be applied to a wide range of complex design challenges including nonlinear dynamics, and is expected to be compatible to most FE codes and derivative-free optimization routines. It holds the potential to facilitate and revolutionize the way (nano)dynamical systems are designed, thus pushing the ultimate performance limits of sensors, mechanisms and actuators for scientific, industrial, and consumer applications.
Methods
Sample fabrication. We produce our nanomechanical resonators using electron beam lithography and reactive ion etching techniques on high-stress \ceSi3N4 layers, chosen for their reliability and precision in achieving design specifications 32. These layers are deposited via low pressure chemical vapor deposition (LPCVD) onto a silicon substrate. Following this, the devices undergo suspension through a fluorine-based deep reactive ion underetching process. The mechanical properties of the high-stress \ceSi3N4 are characterized in our previous works 14, with an initial isotropic stress , Young’s modulus , Poisson’s ratio , mass density . The intrinsic quality factor is a function of thickness 33, which is .
Data availability
The data that support the findings of this study are available from the corresponding authors upon reasonable request.
Acknowledgements
Funded/Co-funded by the European Union (ERC Consolidator, NCANTO, 101125458). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. Z.L. acknowledges financial support from China Scholarship Council, the assistance on the FE reduced-order modeling from Vincent Bos, and the instruction about using the high
performance computing cluster from Binbin Zhang. This work is also part of the project, Probing the physics of exotic superconductors with microchip Casimir experiments (740.018.020) of the research program NWO Start-up which is partly financed by the Dutch Research Council (NWO). M.X. and R.A.N. acknowledge valuable support from the Kavli Nanolab Delft.
Author contributions
Z.L., F.A., P.G.S. and A.M.A. conceived the experiments and methods; M.X. and R.A.N. fabricated the \ceSi3N4 samples; Z.L. conducted the measurements and analysed the experimental data; Z.L. and F.A. built the theoretical model; Z.L. performed the reduced-order modelling of the finite element model; Z.L. and A.S. set up the optimization on high performance cluster; F.A. and P.G.S. supervised the project; and the manuscript was written by Z.L. and P.G.S. with inputs from all authors.
Competing interests
The authors declare no competing interests.
References
- Erbil et al. [2020] S. O. Erbil, U. Hatipoglu, C. Yanik, M. Ghavami, A. B. Ari, M. Yuksel, and M. S. Hanay, Physical review letters 124, 046101 (2020).
- Yuksel et al. [2019] M. Yuksel, E. Orhan, C. Yanik, A. B. Ari, A. Demir, and M. S. Hanay, Nano letters 19, 3583 (2019).
- Bayram et al. [2022] F. Bayram, D. Gajula, D. Khan, and G. Koley, Microsystems & nanoengineering 8, 1 (2022).
- Venstra et al. [2013] W. J. Venstra, H. J. Westra, and H. S. Van Der Zant, Nature communications 4, 2624 (2013).
- Miller et al. [2021] J. M. Miller, A. Gomez-Franco, D. D. Shin, H.-K. Kwon, and T. W. Kenny, Physical Review Research 3, 033268 (2021).
- Høj et al. [2021] D. Høj, F. Wang, W. Gao, U. B. Hoff, O. Sigmund, and U. L. Andersen, Nature communications 12, 5766 (2021).
- Schiwietz et al. [2024] D. Schiwietz, M. Hörsting, E. M. Weig, M. Wenzel, and P. Degenfeld-Schonburg, arXiv preprint arXiv:2403.17679 (2024).
- Dou et al. [2015] S. Dou, B. S. Strachan, S. W. Shaw, and J. S. Jensen, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 373, 20140408 (2015).
- Li et al. [2017] L. L. Li, P. M. Polunin, S. Dou, O. Shoshani, B. Scott Strachan, J. S. Jensen, S. W. Shaw, and K. L. Turner, Applied Physics Letters 110, 081902 (2017).
- Mignolet et al. [2013] M. P. Mignolet, A. Przekop, S. A. Rizzi, and S. M. Spottswood, Journal of Sound and Vibration 332, 2437 (2013).
- Touzé et al. [2021] C. Touzé, A. Vizzaccaro, and O. Thomas, Nonlinear Dynamics 105, 1141 (2021).
- Cenedese et al. [2022] M. Cenedese, J. Axås, B. Bäuerlein, K. Avila, and G. Haller, Nature communications 13, 872 (2022).
- Muravyov and Rizzi [2003] A. A. Muravyov and S. A. Rizzi, Computers & Structures 81, 1513 (2003).
- Li et al. [2024] Z. Li, M. Xu, R. A. Norte, A. M. Aragón, P. G. Steeneken, and F. Alijani, Communications Physics 7, 53 (2024).
- Keşkekler et al. [2023] A. Keşkekler, V. Bos, A. M. Aragón, P. G. Steeneken, and F. Alijani, Phys. Rev. Appl. 20, 064020 (2023).
- Dhooge et al. [2008] A. Dhooge, W. Govaerts, Y. A. Kuznetsov, H. G. E. Meijer, and B. Sautois, Mathematical and Computer Modelling of Dynamical Systems 14, 147 (2008).
- Xu et al. [2024] M. Xu, D. Shin, P. M. Sberna, R. van der Kolk, A. Cupertino, M. A. Bessa, and R. A. Norte, Advanced Materials 36, 2306513 (2024).
- Cupertino et al. [2024] A. Cupertino, D. Shin, L. Guo, P. G. Steeneken, M. A. Bessa, and R. A. Norte, Nature Communications 15, 4255 (2024).
- Li et al. [2023] Z. Li, M. Xu, R. A. Norte, A. M. Aragón, F. Van Keulen, F. Alijani, and P. G. Steeneken, Applied Physics Letters 122, 013501 (2023).
- Coello et al. [2004] C. A. C. Coello, G. T. Pulido, and M. S. Lechuga, IEEE Transactions on evolutionary computation 8, 256 (2004).
- Nayfeh and Mook [2008] A. H. Nayfeh and D. T. Mook, Nonlinear oscillations (John Wiley & Sons, 2008).
- Schmid et al. [2016] S. Schmid, L. G. Villanueva, and M. L. Roukes, Fundamentals of nanomechanical resonators, Vol. 49 (Springer, 2016).
- Demir and Hanay [2019] A. Demir and M. S. Hanay, IEEE Sensors Journal 20, 1947 (2019).
- Manzaneque et al. [2023] T. Manzaneque, M. K. Ghatkesar, F. Alijani, M. Xu, R. A. Norte, and P. G. Steeneken, Physical Review Applied 19, 054074 (2023).
- Samanta et al. [2023] C. Samanta, S. De Bonis, C. Møller, R. Tormo-Queralt, W. Yang, C. Urgell, B. Stamenic, B. Thibeault, Y. Jin, D. Czaplewski, et al., Nature Physics , 1 (2023).
- Foster et al. [2016] A. Foster, J. Maguire, J. Bradley, T. Lyons, A. Krysa, A. Fox, M. Skolnick, and L. Wilson, Nano Letters 16, 7414 (2016).
- Antonio et al. [2012] D. Antonio, D. H. Zanette, and D. López, Nature communications 3, 1 (2012).
- Chen et al. [2017] C. Chen, D. H. Zanette, D. A. Czaplewski, S. Shaw, and D. López, Nature communications 8, 1 (2017).
- Badzey and Mohanty [2005] R. L. Badzey and P. Mohanty, Nature 437, 995 (2005).
- Belardinelli et al. [2023] P. Belardinelli, W. Yang, A. Bachtold, M. Dykman, and F. Alijani, arXiv preprint arXiv:2312.14034 (2023).
- Lake et al. [2013] J. J. Lake, A. E. Duwel, and R. N. Candler, Journal of Microelectromechanical Systems 23, 364 (2013).
- Shin et al. [2022] D. Shin, A. Cupertino, M. H. de Jong, P. G. Steeneken, M. A. Bessa, and R. A. Norte, Advanced Materials , 2106248 (2022).
- Villanueva and Schmid [2014] L. G. Villanueva and S. Schmid, Physical review letters 113, 227201 (2014).