REVIEW 4 major objections 3 minor 28 references
Actuation and mapping of SAW-induced high-frequency wavefields on suspended graphene membranes
T0 review · 4 major / 3 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Suspended graphene membranes can be actuated at 375 MHz by surface acoustic waves, and atomic-force acoustic microscopy maps the resulting standing wavefields, showing a wavelength reduction from 10 µm to about 1.8–2.3 µm.
desk verdict Useful new method for driving and imaging high-frequency wavefields on suspended graphene, but the dispersion fit has an internal inconsistency that needs fixing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the pre-tensioned plate dispersion relation $\omega_k = \sqrt{(D/\rho h)\,k^4 + (T/\rho h)\,k^2}$, with $D$ the bending rigidity, $T$ the tension, $\rho h$ the areal mass density, and $k$ the wavenumber. It connects the low-frequency clamped-circular modes (wavenumbers $k_{01}=3.1962/R$, $k_{11}=4.6109/R$, $k_{21}=5.9057/R$) measured by laser interferometry to the high-frequency wavenumber extracted from AFAM fringe spacing, so that effective $(D/\rho h)$ and $(T/\rho h)$ can be fitted. The enabling readout is AFAM: the cantilever's nonlinear tip-sample interaction converts the 375 MHz, amplitude-modulated vibration into a 10 kHz demodulated signal, producing a spatial map of the standing wavefield with sub-micron resolution.
What would settle it
Repeat the AFAM measurement on the same drum at several tip setpoint forces and with different cantilevers, and compare the extracted wavelengths to an independent sub-micron displacement probe (for example a stroboscopic electron-microscopy or short-wavelength interferometric measurement) at the same 375 MHz drive. If the fringe spacing changes with contact conditions or disagrees with the independent wavelength, the half-wavelength interpretation and the fitted $(D/\rho h)$ and $(T/\rho h)$ values do not hold.
Extended reading notes
Core claim
The central claim is that a Rayleigh surface wave arriving at a suspended graphene drum is converted into a flexural wave with a much shorter wavelength, and that this conversion can be quantified. Using 375 MHz SAWs generated by interdigital transducers on lithium niobate, the authors excite five multilayer graphene drums of 9.4–37 nm thickness and map the standing wavefields with AFAM, whose nonlinear tip-sample interaction demodulates the amplitude-modulated drive into time-averaged spatial fringes. Fringe spacings give wavelengths of 1.83–2.28 µm on suspended graphene versus 10 µm on supported graphene. Combining these points with the first three clamped-circular-drum resonance frequencies measured by laser interferometry, the authors fit the dispersion relation $\omega_k = \sqrt{(D/\rho h)\,k^4 + (T/\rho h)\,k^2}$ and find phase velocities from roughly 160 m/s to 700 m/s, with fits at $R^2>0.99$. Finite-element simulations using the fitted stiffness parameters reproduce the measured wavefield geometry and dimensions.
Load-bearing premise
The quantitative interpretation rests on assuming that every bright AFAM fringe on the suspended graphene marks half a local flexural wavelength, and that the flake behaves as a smooth, uniformly stretched, continuously clamped circular plate with the standard mode wavenumbers; if tip-sample contact averaging, topography, or a wrong mode assignment changes the fringe-to-wavelength relation, the extracted wavelengths and the fitted bending rigidity and tension are systematically wrong.
Editorial extensions
If this is right
- SAW actuation should work on any suspended two-dimensional membrane, conducting or not, because the coupling is purely mechanical, and interdigital transducer design can push the drive into the gigahertz range.
- Suspended membranes act as local acoustic lenses: a fixed 10 µm SAW wavelength is compressed to about 2 µm over the drum, so cavity size, shape, and flake thickness can pattern nodal lines and localized wavefields.
- The dispersion fit means the effective bending rigidity and tension of a multilayer flake can be extracted from two measurement bands—megahertz drum resonances and the 375 MHz AFAM point—without needing optical resolution to see the short wavelength.
- Because the high-frequency branch is bending-dominated ($\omega \propto k^2$), thicker and stiffer membranes will compress the SAW wavelength more, making the compression tunable through flake thickness.
- Probe-based mapping is required for these short wavelengths; non-contact optical Doppler vibrometry with micrometer-scale resolution cannot resolve the 1.83–2.28 µm features.
Reading between the lines
- If the half-wavelength fringe interpretation carries over to other materials, the same platform should map frequency-dependent flexural fields on suspended hexagonal boron nitride or MoS2 membranes, and the compression factor should grow as thickness decreases because $D$ scales as $h^3$.
- Because AFAM contrast depends on nonlinear tip-sample contact, a controlled setpoint-force sweep would reveal how much of the fitted stiffness is intrinsic to the membrane rather than induced by the measurement; the paper does not perform that control.
- Detuning or removing the reflector IDT would turn the standing-wave measurement into a propagating-wave measurement, separating intrinsic flexural dispersion from resonator boundary effects.
- The sub-2 µm acoustic standing patterns are natural templates for trapping or moving nanoparticles on a chip, and could be extended to reconfigurable sorting, though the paper stops short of demonstrating transport.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a combined experimental platform in which 375 MHz surface acoustic waves (SAWs) on LiNbO3 actuate suspended multilayer graphene membranes, and atomic force acoustic microscopy (AFAM) maps the resulting standing wavefields. The central observation is that the acoustic wavelength on the suspended graphene is reduced from the 10 µm SAW wavelength to 1.83–2.28 µm, with phase velocities increasing from ~160 m/s at low frequencies to ~700 m/s at 375 MHz. The authors interpret the data with a pre-tensioned plate dispersion model, fit effective bending rigidity and tension parameters per device, and compare the measured wavefields with COMSOL simulations. They also show wavefield modification for cavities of different size and shape.
Significance. If the quantitative interpretation is sound, this work would demonstrate a useful on-chip, material-agnostic route to high-frequency actuation of suspended 2D membranes, with accessible nanoscale readout via AFAM. The direct mapping of the wavelength reduction on suspended graphene is an advance over optical techniques with micron-scale resolution, and the demonstration of geometry-controlled wavefield patterning is a strength. However, the load-bearing dispersion analysis has internal inconsistencies that currently undermine the fitted parameters and the claimed consistency with pre-tensioned plate theory. The experimental phenomenology, especially the wavelength reduction and the circular fringe patterns, remains credible and valuable; the quantitative modeling needs substantial revision.
major comments (4)
- [Pre-tensioned plate model] The wavenumbers k01 = 3.1962/R, k11 = 4.6109/R, and k21 = 5.9057/R are the roots of the clamped circular plate with zero tension. The model in Eq. (1) includes a tension term T k^2, so for a pre-tensioned plate these constants are not the mode wavenumbers. In the reported tension-dominated regime (phase velocities of ~160 m/s at low frequencies imply T/rho_h >> D/(rho_h R^2)), the true wavenumbers approach the membrane Bessel zeros (j01 = 2.405, j11 = 3.832, j21 = 5.136), i.e., they are 25–30% smaller than the zero-tension plate roots used in the fit. This internal inconsistency biases the fitted (D/rho_h)_eff and (T/rho_h)_eff and invalidates the reported R^2 > 0.99 as evidence for the pre-tensioned plate model.
- [Figure 3 and Table 1] The dispersion curve is fitted with two free parameters, (D/rho_h)_eff and (T/rho_h)_eff, to only three or four data points per device (three low-frequency resonances and one AFAM point at 375 MHz). With this number of points, a two-parameter fit can achieve R^2 > 0.99 almost by construction; this does not provide independent confirmation of the plate model. The paper should report confidence intervals on the fitted parameters and, ideally, cross-validate using more than one high-frequency wavenumber per device.
- [Simulation of SAW-induced wavefields on suspended graphene membranes] The COMSOL model for device D4 uses the fitted (D/rho_h)_eff and (T/rho_h)_eff as inputs, so the good agreement between the simulated and experimental wavefields in Figure 4 is partly by construction. This agreement is therefore not an independent validation of the pre-tensioned plate model. The qualitative pattern geometry and the scaling with cavity size/shape remain informative, but the quantitative match should be presented as a consistency check, not as predictive confirmation.
- [AFAM wavelength extraction] The quantitative analysis assumes that the spacing between adjacent bright fringes in the AFAM maps equals half the local flexural wavelength on the suspended graphene. Tip-sample nonlinearities, topography coupling, or a misinterpretation of the demodulated signal could systematically affect this spacing. The paper should provide an explicit calibration or a control experiment (e.g., comparing the fringe spacing on supported graphene to the known SAW wavelength) to justify that the extracted wavenumber k_exp is quantitatively accurate to the level needed for the dispersion fit.
minor comments (3)
- [Figure 3 caption] The caption states 'D1, D3 and D4' but the text refers to devices with 9.4 nm, 14 nm, and 34 nm thickness, which are D1, D2, and D4. The caption should be corrected.
- [Pre-tensioned plate model text] The statement 'for thin flakes (<20 layers), tension dominates over bending rigidity' is not directly tied to the measured thickness range of 9.4–37 nm (about 28–110 layers); the transition also depends on frequency and wavenumber through Eq. (1). Please rephrase to avoid an apparent inconsistency.
- [Abstract] The abstract says '~2 µm' and later the text gives the range 1.83–2.28 µm; please keep the values consistent throughout, including the abstract.
Circularity Check
Dispersion 'prediction' is an in-sample two-parameter fit; D4 COMSOL match re-inserts fitted E and T, making the main wavefield 'prediction' partly by construction, though the direct AFAM wavelength measurement and generic-parameter simulations retain independent content.
-
fitted input called prediction
[Pre-tensioned plate model, Figure 3 (pp. 11-13)]
"A fitting procedure is performed with the laser interferometer and AFAM data based on equation (1). ... Experimental data points fit well with the estimated dispersion curves, with the goodness parameter R2>0.99 for all devices."
The four points used in the fit are the same points displayed as agreeing with the fitted curve. Because (D/rho_h)_eff and (T/rho_h)_eff are chosen to minimize residuals of exactly these points, R2>0.99 quantifies in-sample fit quality, not an out-of-sample prediction. The abstract's statement that the phase-velocity change is 'consistent with ... plate theory' is therefore weaker than claimed: the curve is not an independent theoretical prediction but a two-parameter least-squares description of the same data.
-
fitted input called prediction
[Simulation of SAW-induced wavefields on suspended graphene membranes (pp. 14-15), device D4]
"The Young's modulus (E) and pretension (T) values obtained from the fitted bending rigidity were inserted into the model, and the vibration response was simulated. ... From the figure, we clearly observe that numerical estimation and experimental measurements align well, and thus the model can be used to predict wavefields on suspended graphene layers"
For D4, the simulated 375 MHz wavefield is produced with E and T fitted from the same D4 data, including the AFAM wavenumber at 375 MHz that is one of the fitted points. The simulation therefore reproduces a fitted input, so its agreement with the experimental map is partly in-sample and the word 'predict' overstates independence. The wavefield geometry is an additional matching element, and the later square/smaller-cavity simulations with generic E=1 TPa provide a genuinely independent check, making the circularity partial rather than total.
full rationale
The direct experimental content—SAW actuation at 375 MHz, AFAM fringe mapping, and the measured wavelength reduction from 10 um to 1.83-2.28 um—is self-contained and not circular. The circularity concerns sit in the interpretative layer. First, the dispersion 'agreement' in Figure 3 is an in-sample two-parameter fit of Eq. (1) to the same laser-interferometry and AFAM points, so R2>0.99 only describes residuals and does not independently confirm the plate-theory prediction claimed in the abstract. Second, the D4 COMSOL comparison inserts E and T obtained from that fit into the simulation and compares the simulated 375 MHz wavefield to the same device whose 375 MHz wavenumber was one of the fitted points; the wavelength agreement is partly by construction, although the spatial pattern comparison adds non-fitted content. The additional square/smaller-cavity simulations with generic E=1 TPa are genuinely predictive and provide partial independent support, which prevents the score from being higher. A separate correctness issue—the low-frequency wavenumbers are taken from zero-tension clamped-plate roots kR=3.1962, 4.6109, 5.9057 while Eq. (1) includes tension—is a model-application inconsistency rather than a circularity and does not drive the score. No load-bearing self-citation chain is present.
Assumptions & free parameters
free parameters (2)
- Effective bending rigidity parameter (D/rho h)_eff per device =
Not reported in the main text; referenced to Supporting Information S4
- Effective tension parameter (T/rho h)_eff per device =
Not reported in the main text; referenced to Supporting Information S4
assumptions (4)
- domain assumption The flexural wave dispersion relation for a pre-tensioned plate, Eq. (1), applies to suspended multilayer graphene.
- domain assumption Clamped circular plate mode wavenumbers k01=3.1962/R, k11=4.6109/R, k21=5.9057/R describe the measured resonance modes.
- domain assumption AFAM fringe spacing on the suspended graphene equals half the local flexural wavelength.
- ad hoc to paper The COMSOL model uses linear elastic material properties with fitted rigidity and tension values and ideal clamped boundary conditions.
Cite this review
Pith. "Pith review of Actuation and mapping of SAW-induced high-frequency wavefields on suspended graphene membranes." pith.science (2026). https://pith.science/paper/CLHIQC4A
@misc{pith2026241218310,
author = {Pith},
title = {Pith review of: Actuation and mapping of SAW-induced high-frequency wavefields on suspended graphene membranes},
year = {2026},
howpublished = {\url{https://pith.science/paper/CLHIQC4A}},
note = {Machine review of arXiv:2412.18310}
}
read the original abstract
High frequency acoustic devices based on two-dimensional (2D) materials are unique platforms to design and manipulate the spatiotemporal response of acoustic waves for next-generation sensing and contactless actuation applications. Conventional methods for actuating suspended membranes, however, cannot be applied to all 2D materials, or are limited in frequency. There is, therefore, a need for a universal high-frequency, on-chip actuation technique that can be applied to all types of membranes. Here, we demonstrate that surface acoustic waves (SAWs) can be used to efficiently actuate suspended 2D materials by exciting suspended graphene membranes with high-frequency (375 MHz) Rayleigh surface waves and mapping the resulting vibration field with atomic force acoustic microscopy (AFAM). Acoustic waves travelling from supported to suspended graphene experience a reduction in acoustic wavelength from 10 \mu m to ~2 \mu um due to the decrease in effective bending rigidity, leading to a decrease in wave velocity on suspended graphene. By varying the excitation frequency, we observed a change in phase velocity from ~160 m/s to ~700 m/s. This behavior is consistent with the nonlinear dispersion of acoustic waves, as predicted by plate theory, in suspended graphene membranes. The geometry and bending rigidity of the membrane thus play key roles in modulating the acoustic wave pattern and wavelength. This combined SAW actuation and AFAM visualization scheme can give new insights into the fundamentals of acoustic transport at the nanoscale limit and provides a route towards the manipulation of localized wavefields for on-chip patterning and transport over 2D materials surfaces.
Figures
Reference graph
Works this paper leans on
-
[1]
Lemme, M. C. et al. Nanoelectromechanical Sensors Based on Suspended 2D Materials. Research 2020, (2020)
work page 2020
-
[2]
Hill, E. W., Vijayaragahvan, A. & Novoselov, K. Graphene sensors. IEEE Sens J 11, 3161– 3170 (2011)
work page 2011
-
[3]
Schedin, F. et al. Detection of individual gas molecules adsorbed on graphene. Nature Materials 2007 6:9 6, 652–655 (2007)
work page 2007
- [4]
-
[5]
Steeneken, P. G., Dolleman, R. J., Davidovikj, D., Alijani, F. & J van der Zant, H. S. Dynamics of 2D material membranes. 2d Mater 8, 042001 (2021)
work page 2021
-
[6]
Jung, M. et al. GHz nanomechanical resonator in an ultraclean suspended graphene p –n junction. Nanoscale 11, 4355–4361 (2019)
work page 2019
-
[7]
Zalalutdinov, M. K. et al. Acoustic cavities in 2D heterostructures. Nature Communications 2021 12:1 12, 1–11 (2021). 22
work page 2021
-
[8]
Soubelet, P. et al. The lifetime of interlayer breathing modes of few -layer 2H -MoSe2 membranes. Nanoscale 11, 10446–10453 (2019)
work page 2019
Show all 28 references
-
[9]
& Bhaskaran, H
Kumar, M. & Bhaskaran, H. Ultrasensitive Room-Temperature Piezoresistive Transduction in Graphene-Based Nanoelectromechanical Systems. Nano Lett 15, 2562–2567 (2015)
2015
-
[10]
Fan, X. et al. Graphene ribbons with suspended masses as transducers in ultra -small nanoelectromechanical accelerometers. Nature Electronics 2019 2:9 2, 394–404 (2019)
2019
-
[11]
Verbiest, G. J. et al. Detecting Ultrasound Vibrations with Graphene Resonators. Nano Lett 18, 5132–5137 (2018)
2018
-
[12]
Riaud, A. et al. Anisotropic Swirling Surface Acoustic Waves from Inverse Filtering for On-Chip Generation of Acoustic Vortices. Phys Rev Appl 4, 034004 (2015)
2015
-
[13]
Nie, X. et al. Surface acoustic wave induced phenomena in two-dimensional materials. This journal is Cite this: Nanoscale Horiz 8, 158 (2023)
2023
-
[14]
Rezk, A. R. et al. Acoustically-Driven Trion and Exciton Modulation in Piezoelectric Two- Dimensional MoS2. Nano Lett 16, 849–855 (2016)
2016
-
[15]
& Santos, P
Rudolph, J., Hey, R. & Santos, P. V. Exciton transport by surface acoustic waves. Superlattices Microstruct 41, 293–296 (2007)
2007
-
[16]
Fandan, R. et al. Dynamic Local Strain in Graphene Generated by Surface Acoustic Waves. Nano Lett 20, 402–409 (2020)
2020
-
[17]
Laitinen, A. et al. A graphene resonator as an ultrasound detector for generalized Love waves in a polymer film with two level states. J Phys D Appl Phys 52, 24LT02 (2019)
2019
-
[18]
Castellanos-Gomez, A., Singh, V., Van Der Zant, H. S. J. & Steele, G. A. Mechanics of freely-suspended ultrathin layered materials. Ann Phys 527, 27–44 (2015)
2015
-
[19]
Chen, C. et al. Performance of monolayer graphene nanomechanical resonators with electrical readout. Nature Nanotechnology 2009 4:12 4, 861–867 (2009)
2009
-
[20]
Eichler, A. et al. Nonlinear damping in mechanical resonators made from carbon nanotubes and graphene. Nat Nanotechnol 6, 339–342 (2011)
2011
-
[21]
Davidovikj, D. et al. Visualizing the Motion of Graphene Nanodrums. Nano Lett 16, 2768– 2773 (2016)
2016
-
[22]
& Santos, P
Pitanti, A., Yuan, M., Zanotto, S. & Santos, P. V. High -resolution acoustic field mapping of gigahertz phononic crystals with atomic force microscopy. Phys Rev Appl 20, 054054 (2023)
2023
-
[23]
Garcia-Sanchez, D. et al. Imaging mechanical vibrations in suspended graphene sheets. Nano Lett 8, 1399–1403 (2008)
2008
-
[24]
Lee, D. et al. Direct Visualization of Gigahertz Acoustic Wave Propagation in Suspended Phononic Circuits. Phys Rev Appl 16, 034047 (2021). 23
2021
-
[25]
Mucientes, M. et al. Mapping nanoscale dynamic properties of suspended and supported multi-layer graphene membranes via contact resonance and ultrasonic scanning probe microscopies. Nanotechnology 31, 415702 (2020)
2020
-
[26]
Keşkekler, A., Arjmandi -Tash, H., Steeneken, P. G. & Alijani, F. Symmetry -Breaking- Induced Frequency Combs in Graphene Resonators. Nano Lett 22, 6048–6054 (2022)
2022
-
[27]
M., Alijani, F
Keşkekler, A., Bos, V., Aragón, A. M., Alijani, F. & Steeneken, P. G. Multimode Nonlinear Dynamics of Graphene Resonators. Phys Rev Appl 20, 064020 (2023)
2023
-
[28]
Castellanos-Gomez, A. et al. Deterministic transfer of two-dimensional materials by all-dry viscoelastic stamping. 2d Mater 1, 011002 (2014). Table of content figure
2014
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.