{"id":"d817e068-71cf-49c7-9f52-16f4ed50fe30","arxiv_id":"2607.06161","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"The capacitance gradient in heterodyne KPFM is exactly representable by a convergent Taylor series, with explicit effective coefficients for static and first/second-eigenmode components.","lead":"This paper shows that the capacitance gradient in heterodyne Kelvin probe force microscopy can be described by an exact, convergent Taylor series, replacing the usual first-order 'small-amplitude' approximation. It gives explicit formulas for the static and first/second eigenmode coefficients and numerically confirms them, which matters for correctly interpreting KPFM sideband measurements.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bimodal coefficient formulas (Eqs. 39–40) are derived for incommensurate frequencies; for commensurate eigenfrequencies, extra high-order contributions enter the static/ω1/ω2 components, so the claim of exactness needs qualification.","rationale":"The paper's central achievement is a rigorous convergence proof for the non-truncated Taylor expansion of the Hudlet-model CG, and the monomodal Taylor–Fourier correspondence is exact. The bimodal extension is more delicate: Eqs. (39)–(40) are derived by selecting spectral components under the assumption that ω1 and ω2 are incommensurate. For commensurate ratios (including the paper's own f2/f1=6.3), additional integer combinations contribute to the same Fourier lines at high Taylor order. The paper openly states this ('The compact expressions derived for incommensurate frequencies are used here') and restricts its numerical validation to n≤30, where the extra terms (n≥72) do not yet enter. Therefore the work is not internally inconsistent; it is precisely scoped. The concern is that the abstract and the reader's strongest_claim present K0, K1, K2 as exactly giving the spectral components independent of amplitude and frequency relationship, which is only true for incommensurate frequencies or finite truncation orders. This does not change the overall verdict because the TSE representation itself is exact and the coefficient formulas are correct within their stated domain. A concrete numerical test would quantify the missing terms and confirm the paper's caveat. The reader's identified quasi-static assumption is indeed a limitation, but the commensurability issue is a more specific and immediately testable gap in the exactness of the bimodal coefficients.","tokens_in":22397,"tokens_out":16150,"duration_ms":149068,"concrete_test":"For the paper's commensurate ratio f2/f1=6.3, choose parameters where high-order CG derivatives are significant (e.g., zc=0.6 nm, z1,0=0.5 nm, z2,0=0.1 nm). Evaluate the exact CG C^(1)(z(t)) from the Hudlet model along the bimodal trajectory over one super-period; compute the Fourier coefficient at ω1 by projection. Compare with the infinite incommensurate series K1z1,0 from Eq. (40b) evaluated to convergence. If the difference is nonzero, Eq. (41b) is not exact for commensurate frequencies. To isolate the mechanism, also compute the TSE truncated at n=71 and n=200; the difference should vanish for n<72 and appear for n≥72, confirming the missing terms.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations (39)–(40) define K0, K1, K2 from a Taylor–Fourier expansion that selects components at static, ω1, and ω2 assuming these frequencies are incommensurate. For commensurate frequencies (e.g., the paper's f2=6.3f1=63fs, f1=10fs), integer combinations m1ω1+m2ω2 can coincide with ω1 for nonzero m2 (e.g., m1=1−10k, m2=63k). Such terms appear at Taylor orders n=|m1|+m2 ≥ 72 and contribute to the exact Fourier coefficient X_p1, but are absent from the 'incommensurate' Λ1. Thus K1z1,0 = X_p1 (Eq. 41b) is not exact for a non-truncated TSE in the commensurate case; it is valid only below the first such order. The paper acknowledges this in Sec. III.G.2 and states these contributions do not enter 'over the Taylor-order range considered here' (n ≤ 30 in the numerical validation), but the abstract and the reader's strongest_claim present the coefficients as generally exact. This is the least secure point in the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a non-truncated Taylor-series description of the tip–sample capacitance gradient (CG) in heterodyne KPFM, for both monomodal and bimodal cantilever motion. It proves convergence of the Taylor expansion for a Hudlet-based CG model under the physically admissible condition zmin > 0, derives closed-form expressions for the effective CG coefficients K0, K1, K2 that govern the static, first-eigenmode, and second-eigenmode spectral components, and introduces order-truncation regimes (OTRs) with a term-significance criterion. Numerical simulations validate the convergence of Taylor-based coefficients to Fourier coefficients and map the OTRs in the (zc, z1,0) plane.","tokens_in":22709,"tokens_out":9722,"duration_ms":90090,"significance":"If the claims hold, the paper provides a rigorous replacement for the conventional first-order Taylor approximation, with explicit convergence guarantees and quantitative truncation criteria. The convergence proof for the elementary functions 1/(z+ξ) and ln(z+a) is clean and correct, and the closed-form coefficient expressions (Eqs. 25 and 40) are valuable for analyzing heterodyne KPFM observables. The numerical validation is thorough and reproducible from the stated parameters. The main caveat is the treatment of frequency commensurability, which requires qualification in the abstract and conclusions.","major_comments":[{"comment":"The identification of the effective coefficients with Fourier coefficients, Xp1 = K1 z1,0 and Xp2 = K2 z2,0, is stated without qualification, but it is exact only for incommensurate frequencies. For commensurate frequencies, integer combinations m1ω1 + m2ω2 can coincide with ω1 or ω2. For the paper's own ratio f1=10fs, f2=63fs, the first extra contribution to the ω1 component appears at Taylor order n = |1−63k| + 10k = 72 (k=1), and to the ω2 component at n=74. Thus for a non-truncated TSE, Eqs. (41b)–(41c) are not exact for commensurate trajectories; they hold only below the first coincidence order. The main text concedes this only in the sentence about the 'Taylor-order range considered here' and in SI Sec. III.G.2, while the abstract and conclusion present the coefficients as exact without this restriction. Since the abstract's claim of a 'rigorous spectral description' and 'effective","section":"Sec. II.C.2, Eqs. (39)–(41)"},{"comment":"The text states that the framework 'applies to integer-multiple, commensurate, and incommensurate frequency relationships,' but the compact formulas presented in Sec. II.C.2 are derived only for incommensurate frequencies. The general commensurate case is not worked out in the main text; the SI only notes the absence of extra contributions for n ≤ 30. This is misleading: a reader cannot infer the validity range from the main text alone. I recommend either presenting the general coincidence-order condition in the main text or clearly marking Eqs. (39)–(40) as valid for incommensurate frequencies and for commensurate frequencies only below the first coincidence order.","section":"Sec. II.C.1"}],"minor_comments":[{"comment":"The abstract and conclusion should carry the same caveat about commensurate frequencies as the body text. As written, they assert general exactness that is not warranted for all frequency relationships.","section":"Abstract and Conclusion"},{"comment":"The OTR order ℓ* is defined as the maximum q with ri,q ≥ τ. If the significance ratios are non-monotonic, the 'highest' significant order may not form a contiguous regime. Please clarify whether monotonic decay is assumed or whether the definition should be interpreted as a set of significant orders.","section":"Eq. (30)"},{"comment":"The sentence 'The compact expressions derived for incommensurate frequencies are used here' is placed after Eq. (39). It would be clearer to state this before presenting the formulas and to specify that the numerical validation uses n ≤ 30, where the extra commensurate terms are absent.","section":"Sec. II.C.2, text before Eq. (40)"},{"comment":"The color map uses a single hue gradient for HOTR-2 through HOTR->5; distinguishing the boundaries would be easier with a discrete color bar and explicit labels for each regime.","section":"Fig. 10"},{"comment":"The name 'Benjamin Gr´evin' contains a formatting artifact (the accent appears as a combining character). Please ensure the typesetting is correct.","section":"Author list"},{"comment":"The sentence 'The static, f1, and f2 components retained in the Taylor-Fourier comparison remain clearly identifiable' is correct, but it would help to state explicitly that the numerical projection in Eq. (54) captures all components at f1 and f2, including any high-order commensurate terms, which is why the comparison is valid for n ≤ 30.","section":"Sec. III.B.4"}],"recommendation":"major_revision","confidential_remarks":"The central derivation is sound for incommensurate frequencies, and the numerical validation is careful. The main issue is the overstatement of exactness for commensurate frequencies, which appears in the abstract and conclusion. This is correctable by qualification rather than new derivation. After such revision, I would be willing to accept."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this paper actually delivers a rigorous replacement for the first-order small-amplitude approximation in heterodyne KPFM. It gives closed-form effective coefficients K0, K1, K2 for the capacitance-gradient spectral components, proves convergence of the Taylor expansion for a realistic Hudlet-based model, and replaces the vague 'small-amplitude regime' with an explicit order-truncation criterion based on term significance. The math is clean: the convergence proof uses simple singularity arguments and the admissibility condition zmin>0; the coefficient formulas follow from binomial/parity selection. The numerical validation is honest: it projects the same model, demonstrating internal consistency, not fitting to data. That is the right kind of validation for a theoretical framework paper.\n\nThe main soft spot is in the bimodal section. The formulas (39)-(40) are derived for incommensurate frequencies. For commensurate eigenfrequencies — including their own f2=6.3f1 — integer combinations of the two tones can land exactly on ω1 or ω2. For their ratio, those contributions first appear at Taylor order n≥72, so they don't show up in the n≤30 numerics. The paper acknowledges this in Sec. III.G.2 of the SI, and the main text says 'the compact expressions derived for incommensurate frequencies are used here.' But the abstract and some of the claims (and the reader's strongest take) present the bimodal coefficients as generally exact. That overstates it. The correct statement is: exact for incommensurate frequencies; for commensurate ones, exact below the first coinciding combination, which is a high order for typical ratios but should be spelled out in the main text, not buried in the SI.\n\nThe quasi-static assumption is stated explicitly, so no complaint there. The paper is a solid formal contribution. It is not a new physical effect paper; it is a proper mathematical foundation for existing KPFM observables. I'd send it out for review: the derivation deserves careful checking, and the SI should be read in full. The authors are clearly serious and the work is coherent within its own framework.","headline":"A rigorous Taylor-series framework for KPFM capacitance gradients that replaces the small-amplitude approximation, with a real caveat about the bimodal exactness claim.","tokens_in":23192,"tokens_out":2859,"would_cite":true,"duration_ms":28043,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A non-truncated Taylor expansion of the capacitance gradient exactly captures heterodyne KPFM spectral components at any oscillation amplitude.","keywords":["Kelvin probe force microscopy","capacitance gradient","heterodyne detection","Taylor series convergence","bimodal atomic force microscopy","electrostatic force spectroscopy","order-truncation regimes","Fourier series"],"falsifier":"Measure the ωmod sideband amplitude in a heterodyne experiment with a well-characterized sample while sweeping z1,0 over a range where K1 and C'' differ by more than a few percent; a deviation from the predicted K1 scaling that cannot be explained by the Hudlet geometry would falsify the model.","tokens_in":22319,"feed_emoji":"🔬","tokens_out":3769,"duration_ms":37574,"temperature":0.7,"pith_summary":"The paper shows that the usual small-oscillation-amplitude approximation in heterodyne Kelvin probe force microscopy is unnecessary. For a realistic Hudlet-based tip–sample capacitance model, the full Taylor series of the capacitance gradient converges for any physically admissible oscillatory trajectory, both monomodal and bimodal. This yields exact effective coefficients K0, K1, K2 that determine the static, first-mode, and second-mode components of the electrostatic interaction, and provides explicit term-significance criteria for when truncated expansions are accurate. If correct, the standard first-order treatment becomes a special case within a controlled hierarchy, and second-eigenmode observables are shown to be especially sensitive to higher-order capacitance derivatives.","feed_headline":"Capacitance-gradient Taylor series proven exact at any tip amplitude","feed_subtitle":"New K0–K2 coefficients replace the small-amplitude limit in heterodyne Kelvin probe microscopy, with explicit truncation criteria.","key_machinery":"The Hudlet-based capacitance-gradient model, split into cantilever, cone, and apex contributions, each a linear combination of functions 1/(z+ξ) and ln(z+a). The convergence proof maps the trajectory to a real variable x = δz/zc and shows |x| < 1 ≤ β, where β is the distance to the nearest singularity of these elementary functions, so the Taylor series converges term-by-term. The effective coefficients K0, K1, K2 (Eqs. 39–40) then carry the argument: they are infinite sums of higher-order CG derivatives with amplitude-weighted combinatorial factors, and the order-truncation regime criterion (comparing significance ratios ri,q against a threshold τ) decides when the leading first-order term s","core_discovery":"The central claim is that the capacitance gradient C'(z(t)) along a sinusoidal tip trajectory can be represented exactly by its non-truncated Taylor series about the average tip–surface distance zc, provided the tip never touches the surface (zmin > 0). For the Hudlet-based model, the proof relies on the capacitance gradient being a linear combination of rational and logarithmic functions whose nearest singularities lie beyond the sampled distance range; hence the series converges regardless of oscillation amplitude. The static, ω1, and ω2 spectral components are then given exactly by the coefficients K0, K1, K2, which are infinite sums of higher-order spatial derivatives weighted by amplitu","pith_inferences":["If the quasi-static assumption breaks down, for instance due to dielectric relaxation in the sample at the bias-modulation frequency, the position-only Taylor expansion would no longer describe the actual force; discrepancies between K1-based predictions and measured sidebands could serve as a diagnostic for such rate-dependent capacitance.","The same Taylor–Fourier machinery could be applied to other oscillating-probe observables that are nonlinear functions of tip–sample distance, such as van der Waals force gradients or microwave impedance, yielding analogous effective coefficients and truncation criteria.","Because K2 is more sensitive to higher-order derivatives, deliberately tuning z1,0 and z2,0 into the higher-order regime could enhance contrast to subsurface or short-range dielectric variations, turning the small-amplitude approximation's failure into a useful measurement channel."],"forward_implications":["Open-loop amplitude-modulation heterodyne KPFM can be modeled at arbitrary first-eigenmode amplitude by replacing C''(zc) with K1(zc), removing the vague small-amplitude restriction.","The second-eigenmode coefficient K2 enters higher-order truncation regimes earlier than K1, making second-mode amplitude and phase observables more sensitive to short-range capacitance nonlinearity.","Taylor-based coefficients converge to Fourier coefficients with exponentially decreasing error as truncation order increases, giving practical stopping rules for numerical simulation.","The framework extends to incommensurate bimodal frequencies, where no finite Fourier super-period exists, because the Taylor representation is frequency-agnostic.","The static coefficient K0 provides an exact expression for the DC electrostatic-force channel, which matters for interpreting force offsets in Kelvin probe measurements."],"fun_headline_variants":["Exact capacitance gradient Taylor series for any amplitude","Taylor series for capacitance gradient proven convergent","Capacitance gradient exact without small-amplitude limit","New coefficients replace small-amplitude KPFM approximation","Heterodyne KPFM: rigorous spectral description without truncation"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The capacitance and its gradient depend only on the instantaneous tip–surface distance, not on how fast the tip moves or on the bias-modulation frequency; if this quasi-static description fails, the position-only Taylor expansion no longer describes the actual electrostatic force.","fun_headline_variants_meta":{"raw":{"variants":["Exact capacitance gradient Taylor series for any amplitude","Taylor series for capacitance gradient proven convergent","Capacitance gradient exact without small-amplitude limit","New coefficients replace small-amplitude KPFM approximation","Heterodyne KPFM: rigorous spectral description without truncation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000359,"raw_usage":{"total_tokens":1817,"prompt_tokens":820,"completion_tokens":997,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":919}},"tokens_in":564,"tokens_out":997,"duration_ms":7169,"temperature":1.0,"reasoning_tokens":919,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T08:17:01.208573+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the ωmod sideband amplitude in a heterodyne experiment with a well-characterized sample while sweeping z1,0 over a range where K1 and C'' differ by more than a few percent; a deviation from the predicted K1 scaling that cannot be explained by the Hudlet geometry would falsify the model.","supporting_citations":[],"review_version":2}