{"id":"f2ef2efa-a16e-49be-ab6b-97ece930542c","arxiv_id":"2412.00258","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A deterministic finite-element ring-down calculation combined with the fluctuation-dissipation theorem reproduces the Brownian noise spectra of the first eleven flexural modes of a doubly-clamped nanobeam in air and water.","lead":"This paper computes the random thermal vibrations of a tiny clamped beam in air or water using one deterministic finite element simulation rather than many random ones. It reports accurate noise spectra for the first eleven bending modes, including the effects of beam tension and a nearby solid wall.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The largest deviations occur exactly at modes 10-11 in air, but no convergence study of the stochastic spectra is reported, so the 'first eleven modes' claim rests on unverified numerical resolution.","rationale":"The paper is a credible application of the standard fluctuation-dissipation route to deterministic finite-element ring-downs. The FDT formalism is sound, there are no fitted parameters, and the natural-frequency validation to within 0.3% is genuine supporting evidence. The water results agree convincingly with the analytical theory, and the floor comparison with Tung et al. is a meaningful additional check. My concern is specifically with the quantitative air claim, which is the quantitative backbone of the abstract: the largest errors in Tables IV and V occur at the highest modes, exactly where temporal integration and fluid-mesh resolution are most stressed. The paper does not report a convergence study of the stochastic spectra themselves, only of natural frequencies and static displacements. A systematic frequency error growing to -1% and height errors of -6% to -9% at mode 11 is consistent with generalized-alpha numerical dissipation or an under-resolved Stokes layer, and the absence of FFT/windowing details leaves spectral processing as another uncontrolled factor. This does not undermine the method's validity, but it does mean the 'first eleven modes' accuracy claim is not yet fully established. I partially agree with the reader: the continuum/no-slip assumption at Wi_11 + Kn = 0.15 is worth noting, but 0.15 is an order of magnitude below the cited <~1 continuum threshold, so it is less immediately load-bearing than the numerical convergence question. If the proposed refinement test shows negligible shifts, then the continuum/molecular regime becomes the leading remaining risk; if it shows large shifts, the current agreement is partly numerical. The verdict should remain conditional pending this check.","tokens_in":18562,"tokens_out":13873,"duration_ms":144629,"concrete_test":"Rerun the air case x0 = 1/4, lf = 10 µm with dt = tau_1/640 (doubled temporal resolution) and a fluid mesh refined so that the mode-11 Stokes layer (delta_s,11 ~ 0.23 µm) is covered by at least 4-5 elements near the beam; compute the FFT with zero-padding to reduce binning artifacts. Compare f_p,10, f_p,11, G_W,p,10, and G_W,p,11 against Tables IV and V. If these values shift by more than roughly half of the current reported errors, the stated 'within 1% / 11%' agreement is not converged and the claim should be revised; if they shift by less, the remaining error is physical or theoretical rather than numerical.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper validates spatial and temporal resolution using natural frequencies and static displacement (Sec. III), but not the noise spectra that constitute the central result. In the air case (Sec. IV.B), Tables IV and V show a systematic degradation with mode number: f_p,11 error reaches -1.03% and G_W,p,11 error reaches -6.1% (lf=10 µm) or -9.2% (lf=2 µm), while modes 10-11 are described as having 'some error evident in the location of the peaks.' Twenty time steps per period with the generalized-alpha method, which has numerical dissipation, is a plausible source of this bias; the Stokes layer for mode 11 in air is about 0.23 µm, and the fluid mesh within it is not described. Because the finite-element spectra are compared only with the same continuum theory used to define the reference, the observed errors could be numerical rather than physical. The reader's continuum/no-slip concern is valid in spirit, but Wi_11 + Kn = 0.15 is an order of magnitude below the quoted <~1 threshold, so the unresolved discretization of the high air modes is the more immediate load-bearing risk.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a computational approach to obtain the multimodal stochastic dynamics of a doubly-clamped nanobeam immersed in a viscous fluid, based on the fluctuation-dissipation theorem. A single deterministic finite-element ring-down simulation is used to compute the autocorrelation and noise spectrum of the beam's displacement for the first eleven flexural modes, including the effects of intrinsic tension and a nearby rigid floor. Results in water and air are compared with analytical and semianalytical predictions (Eqs. 3-6 and Ref. 38), with reported relative errors of about 1% in peak frequencies and up to about 11% in peak amplitudes for the highest modes in air. The authors conclude that the approach can be extended to complex geometries where analytical theories are unavailable.","tokens_in":18791,"tokens_out":6288,"duration_ms":60335,"significance":"The central methodological claim is valuable: if the deterministic ring-down with the fluctuation-dissipation theorem correctly reproduces the noise spectra of many modes simultaneously, it would provide a general numerical route to Brownian dynamics of complex three-dimensional elastic structures in fluids. The study is genuinely parameter-free on the simulation side: no mode expansion, no fitted parameters, and the FEM is an independent first-principles discretization of the coupled elastodynamics and Navier-Stokes equations. The favorable comparison with the experimentally validated theory of Ref. 20 for the lower, well-resolved modes in water is a true strength. The main limitation is the lack of a direct convergence study of the noise spectra for the highest air modes, where the reported deviations are largest.","major_comments":[{"comment":"The claim that the method accurately resolves the first eleven modes is not yet supported for modes 10 and 11 in air, where the peak-frequency error reaches -1.03% and the peak-height error reaches -6.1% (lf=10 µm) or -9.2% (lf=2 µm). The only reported convergence checks are for vacuum natural frequencies and static displacement (Sec. III), which do not control the fluid-loaded spectral peak heights or the numerical dissipation of the generalized-α time integrator. At ~20 time steps per period for mode 11, numerical damping could plausibly lower the high-frequency peaks, and the paper does not describe the fluid mesh inside the air Stokes layer (δ_s,11 ≈ 0.23 µm). A convergence study that halves the time step and refines the fluid mesh, reporting the resulting spectra for modes 8-11 in air, is needed before \"excellent agreement\" can be claimed for all eleven modes.","section":"IV.B / Tables IV and V"},{"comment":"The fluid mesh is essentially not described. The paper states only that the beam mesh uses a maximum length scale on the order of the beam thickness h (Sec. III), but gives no information about the fluid mesh size, element type, boundary-layer resolution, or refinement strategy. Because the Stokes layer thickness for mode 11 in air is about 0.23 µm, the number of elements across this layer directly controls the accuracy of the added mass and damping for the highest modes. Without this information, the numerical resolution of the fluid dynamics cannot be assessed, and the computations cannot be reproduced by other groups.","section":"III"}],"minor_comments":[{"comment":"There are several typographical errors that should be corrected: \"centralto\" in the abstract, \"the the\" in Sec. III, \"noise spectrum spectrum\" in Sec. IV.A, and \"the displacement the displacement is measured\" in the caption of Fig. 5.","section":"Throughout"},{"comment":"Refs. 52 and 54 appear to be the same paper (same journal, volume, article number, and authors); this duplication should be resolved.","section":"References"},{"comment":"The statement that modes 4 and 8 have \"large relative errors\" is imprecise: in Table IV their frequency errors are not conspicuously large, and the large values appear only in the peak-amplitude errors of Table V. The text should clarify that the large relative errors are in the spectral peak magnitudes and are due to the smallness of the signals near nodes.","section":"IV.B / Tables IV and V"},{"comment":"The comparison theory, Eq. (3), is the authors' own framework (Refs. 20, 52-54). This is not a circularity problem because the FEM simulation contains no fitted parameters and the theory was previously checked against experiments, but the paper should explicitly note that for the floor case the green semianalytical curve (Ref. 38) is itself a fit to numerical simulations, so the agreement between the FEM and that curve is partly a comparison between two numerical approaches.","section":"III"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and the approach is promising. The main obstacle is the missing spectral convergence study for the high-mode air results; if the authors can provide it and show that the observed deviations are not numerical artifacts, I would support acceptance. The self-citation pattern is not a concern given that the comparison theory was previously validated experimentally."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it claims: a single deterministic FEM ring-down, run through the fluctuation-dissipation theorem, reproduces the Brownian noise spectra of the first eleven flexural modes of a tensioned, doubly-clamped nanobeam in water and air, with and without a nearby floor. Agreement with the analytical theory is generally excellent, and the method captures overlapping modes and wall effects without modal expansion.\n\nWhat's new is the multimodal generalization. The FDT-FEM idea was already established by this group for the fundamental mode; extending it to eleven modes, including tension and a nearby boundary, is a real step. The comparison against Eq. (3) from their 2023 PR Applied paper is meaningful because that theory was experimentally validated, and the FEM calculation is parameter-free first-principles. Natural-frequency errors under 0.3% and static displacement errors around 0.2% show the mechanics are well resolved. Self-citation is not a real problem here, since the cited theory was externally validated and the numerics are independent.\n\nThe soft spot is exactly where the stress-test note lands. Spatial and temporal resolution is validated on natural frequencies and static displacement, but not on the noise spectra themselves. In air, modes 10 and 11 systematically drift: peak frequency error reaches -1.03% and peak height error -6% to -9%. That pattern is consistent with the ~20 time steps per period plus generalized-alpha numerical dissipation, and the Stokes layer for mode 11 in air (~0.23 µm) gets no mesh description. Without a convergence study on the spectra, we cannot tell whether the high-mode errors are numerical or physical. The reader's continuum/no-slip worry is less pressing: Wi_11+Kn = 0.15 is well below the Kara et al. threshold, so molecular effects are probably not the main culprit. But that shifts the burden to numerical resolution.\n\nAlso, the paper ships no code or data, and Tables IV and V lack error bars. 'Available upon request' is a minor drag, not a fatal one.\n\nWho is this for? Researchers computing or measuring multimode Brownian dynamics of nanomechanical resonators in viscous fluids. They get a validated workflow and a clear statement of where it currently degrades. It deserves a serious referee; the central approach holds up, and the high-mode air discrepancy is a fixable but necessary point for revision.\n\nRecommendation: send it to peer review, with a request for a convergence study of the noise spectra at high modes in air (finer time step and/or mesh, showing the errors shrink), and for the data or a clear reason it cannot be shared.","headline":"A credible multimodal extension of the FDT-FEM approach with strong agreement against an independently validated theory; the main risk is unverified numerical convergence of the high air modes.","tokens_in":19292,"tokens_out":2508,"would_cite":false,"duration_ms":24409,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single deterministic ring-down reproduces the Brownian noise spectrum of eleven nanobeam modes.","keywords":["nanobeam","stochastic dynamics","fluctuation-dissipation theorem","finite element method","multimode noise spectrum","hydrodynamic function","squeeze-film damping","Brownian motion"],"falsifier":"Measure the mode-10 and mode-11 noise peak heights of a comparable tensioned beam in atmospheric air; if the measured peak heights differ from the simulation by more than the reported 11 percent, the continuum damping description is the weak link.","tokens_in":18367,"feed_emoji":"📈","tokens_out":5921,"duration_ms":53784,"temperature":0.7,"pith_summary":"The paper's central claim is that the full Brownian-driven vibration spectrum of a tensioned nanobeam in a viscous fluid can be computed from one deterministic finite-element ring-down, without directly simulating thermal noise. Using the fluctuation-dissipation theorem, the ring-down after removal of a static point force yields the displacement autocorrelation and power spectral density for all modes at once. The authors demonstrate this for the first eleven flexural modes of a doubly-clamped beam in water and air, with and without a nearby rigid floor, matching theory in peak frequency to about one percent and peak height to about eleven percent. If correct, this makes multimode stochastic dynamics computable for experimental geometries where no analytical theory exists.","feed_headline":"One ring-down reproduces 11 nanobeam noise modes","feed_subtitle":"Brownian spectra of a tensioned beam in air or water follow from a single deterministic simulation, to about 1% in frequency.","key_machinery":"The load-bearing object is the fluctuation-dissipation relation $\\langle w(x_0,0)w(x_0,t)\\rangle = (k_B T/F_0)\\, W(x_0,t)$, which converts the deterministic return-to-equilibrium displacement $W$ after removal of a point force $F_0$ into the equilibrium autocorrelation of the stochastic displacement $w$; a cosine transform then gives the noise spectrum. The finite-element computation supplies $W(x_0,t)$ by solving the coupled Navier-Stokes and Euler-Bernoulli problem, including three-dimensional fluid flow, intrinsic tension, and no-slip walls. The hydrodynamic function $\\Gamma(\\omega)$ for a thin blade in unbounded fluid, or the wall-corrected version for a nearby floor, is what the analytical side uses to predict the same spectra, so the comparison tests both the finite-element fidelity and the adequacy of the hydrodynamic function.","core_discovery":"The discovery being argued is that deterministic fluid-solid finite-element computation plus the fluctuation-dissipation theorem, expressed in Eqs. (1)-(2), is sufficient to reproduce the stochastic multimodal response of a nanobeam. No mode expansion is needed: the single computed ring-down retains all modes, including overlapping modes and any fluid-mediated couplings. For the specified silicon-nitride beam, the simulated noise spectra agree with the analytical prediction based on a mode-independent hydrodynamic function for a thin blade, and with the wall-corrected semianalytical prediction when a floor is nearby. The quantitative standard is agreement in peak frequencies within about 1% and peak heights within 0.5% to 11% for the first eleven modes.","pith_inferences":["A natural inverse use is to match computed multimode spectra to measured ones to extract high-frequency fluid properties or effective damping near walls in regimes where simple analytical formulas are untested.","The systematic high-mode errors in air, up to about 11%, hint that a mode-dependent damping correction could refine the continuum description; the same computational setup could test such corrections by varying the ratio of Stokes length to mode wavelength.","Because the simulation is deterministic and parameter-free once geometry, tension, and fluid properties are fixed, it provides a convenient benchmark for calibrating reduced-order models of Brownian nanomechanical sensors.","The method could also compute cross-spectra between different points on the beam or between two nearby beams, quantities relevant to correlation-based force spectroscopy, without any new conceptual steps."],"forward_implications":["The same single-ring-down recipe should give autocorrelations and noise spectra at any axial position on the beam, including points where even and odd modes all contribute and where modal peaks overlap.","The good match with a mode-independent hydrodynamic function means axial-flow corrections are not needed for these eleven modes, and the analytical multimode formula can be trusted for similar slender beams.","The floor's influence is confined mainly to low frequencies: the Stokes length grows as frequency drops, so a nearby wall changes primarily the fundamental-mode peak.","Because the approach is based on linear response, it extends in principle to other modes of motion and other dissipation mechanisms, and to complex three-dimensional structures where no analytical theory exists.","The numerics require no mode expansion, so overlapping and fluid-coupled modes are included automatically, which matters for low-quality-factor cases such as a beam in water."],"supporting_citations":[{"why":"Supplies the fluctuation-dissipation relations, Eqs. (1)-(2), that turn a deterministic ring-down into the stochastic autocorrelation and noise spectrum.","marker":"[5]"},{"why":"Gives the theoretical multimode noise spectrum, Eq. (3), and the experimentally motivated beam and fluid conditions used for comparison.","marker":"[20]"},{"why":"Establishes the deterministic finite-element approach with the fluctuation-dissipation theorem for fundamental-mode stochastic dynamics of elastic cantilevers.","marker":"[32]"},{"why":"Provides the frequency-dependent hydrodynamic function for an oscillating beam used in the analytical predictions.","marker":"[33]"},{"why":"States the fluctuation-dissipation theorem that underlies the entire numerical method.","marker":"[35]"},{"why":"Supplies the semianalytical hydrodynamic function that accounts for a nearby rigid wall, used for the floor-effect comparisons.","marker":"[38]"},{"why":"Provides the natural frequencies and mode shapes of a beam under tensile load, used to validate the solid mechanics.","marker":"[40]"},{"why":"Gives the continuum-validity criterion for unsteady fluid flow that the paper uses to justify applying Navier-Stokes dynamics to the air modes.","marker":"[42]"},{"why":"Provides the detailed theory for the natural frequencies and effective spring constants of a tensioned beam used in the validation.","marker":"[52]"}],"fun_headline_variants":["Single simulation yields all 11 nanobeam noise modes","One deterministic run predicts nanobeam Brownian modes","All 11 nanobeam modes from one force-decay run","One computation gives all 11 flexural noise modes","No eigenmodes needed: one simulation captures thermal noise"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fluid must still behave as a continuous no-slip medium for the fastest air modes, where molecular and timing effects are close to the empirical limit of about 0.15 on the combined scale.","fun_headline_variants_meta":{"raw":{"variants":["Single simulation yields all 11 nanobeam noise modes","One deterministic run predicts nanobeam Brownian modes","All 11 nanobeam modes from one force-decay run","One computation gives all 11 flexural noise modes","No eigenmodes needed: one simulation captures thermal noise"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001095,"raw_usage":{"total_tokens":4550,"prompt_tokens":901,"completion_tokens":3649,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":517,"completion_tokens_details":{"reasoning_tokens":3568}},"tokens_in":517,"tokens_out":3649,"duration_ms":21906,"temperature":1.0,"reasoning_tokens":3568,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:32:47.823283+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the mode-10 and mode-11 noise peak heights of a comparable tensioned beam in atmospheric air; if the measured peak heights differ from the simulation by more than the reported 11 percent, the continuum damping description is the weak link.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the fluctuation-dissipation relations, Eqs. (1)-(2), that turn a deterministic ring-down into the stochastic autocorrelation and noise spectrum."},{"cited_title":"Gress , author J","cited_arxiv_id":null,"evidence_quote":"Gives the theoretical multimode noise spectrum, Eq. (3), and the experimentally motivated beam and fluid conditions used for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the deterministic finite-element approach with the fluctuation-dissipation theorem for fundamental-mode stochastic dynamics of elastic cantilevers."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the semianalytical hydrodynamic function that accounts for a nearby rigid wall, used for the floor-effect comparisons."},{"cited_title":"Bokaian ,\\ title title Natural frequencies of beams under tensile axial loads , \\ https://doi.org/https://doi.org/10.1016/0022-460X(90)90663-K journal journal J","cited_arxiv_id":null,"evidence_quote":"Provides the natural frequencies and mode shapes of a beam under tensile load, used to validate the solid mechanics."},{"cited_title":"Kara , author V","cited_arxiv_id":null,"evidence_quote":"Gives the continuum-validity criterion for unsteady fluid flow that the paper uses to justify applying Navier-Stokes dynamics to the air modes."},{"cited_title":"Barbish , author C","cited_arxiv_id":null,"evidence_quote":"Provides the detailed theory for the natural frequencies and effective spring constants of a tensioned beam used in the validation."}],"review_version":1}