{"id":"8d64aa99-08a7-4079-a01a-91dd569d7a04","arxiv_id":"1908.05230","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"3D PF-SNOM tomography of an h-11BN microdisk reveals quantized phonon-polariton momenta and tip-distance-dependent momentum tuning.","lead":"A new near-field microscopy variant records the optical signal at hundreds of tip heights while scanning, producing a three-dimensional map of infrared near-fields around a tiny boron nitride disk and its edge. The map shows the disk traps infrared waves into discrete momentum states, and that the microscope tip distance can tune them, a capability ordinary near-field microscopes lack.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Momentum-quantization claim rests on parameters fitted to the same three resonances; the Table S1 alignment is not yet distinguished from a post hoc three-point fit.","rationale":"The reader's weakest assumption and my concern are the same: the quantitative quantization claim is calibrated on the same dataset used to test it. The paper's qualitative finding, that q_p is d-independent at certain frequencies for the disk but not for the edge, is supported by the presented images and the second disk. However, the central physical claim goes further and says these locked momenta are the specific Bessel roots of Eq. (1). That step depends on phi fitted to -0.28 pi, on post hoc (s,n) assignment, and on q_p values without error bars. A blind cross-check on the second disk would settle whether the model has predictive power or only descriptive flexibility. Since the reader already requested conditional acceptance pending such checks, I do not move the verdict.","tokens_in":14219,"tokens_out":4976,"duration_ms":53453,"concrete_test":"Use the second disk (Supplementary Fig. S3) as a blind test: fix phi = -0.28 pi and epsilon_s = 2.05 from the first disk, compute all k_sn roots in the measured frequency window, and check whether the second disk's d-independent q_p values fall within propagated uncertainties of any root without reassigning modes or refitting. Also perturb phi between -0.25 pi and -0.30 pi and epsilon_s between 1.2 and 11.7; if any alternative (s,n) root comes within the q_p uncertainty of all three resonances, the mode identification is not unique and the quantization claim should be weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing step is the identification of the three d-independent momenta with Bessel-mode roots. In Supplementary Note 1 the procedure is explicitly: (1) read q_p off the resonant images, (2) solve Eq. S4 for k_sn, and (3) assign (s,n) if they match. The matching uses phi = -0.28 pi, chosen as 'best fit' on these data, and the dispersion curves use epsilon_s = 2.05, chosen 'to match experimental results' (Methods). With only three resonances, residuals of 0.03-0.06 um^-1 (~1-2%), and no stated uncertainty on q_p, the agreement could be a product of the fit rather than a test of it. The adjacent-root spacing of J_s' is on the order of pi/r0 ~ 0.7 um^-1, larger than the residuals, but the freedom in choosing s and in adjusting phi makes a coincidental match possible. The second disk is shown but is not used as a blind prediction with fixed parameters. Thus the experimental observation of d-independent momenta is solid, but the identification with the specific (1,3), (0,4), (1,4) Bessel modes is underdetermined.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a PF-SNOM (peak force scattering-type near-field optical microscopy) method that records a 3D near-field response cube with explicit tip-sample distance dependence, and applies it to an isotopically pure h-11BN microdisk and to an h-11BN edge. The central experimental claims are: (i) the polariton momentum in the circular microdisk takes fixed, frequency-specific values over a small range of infrared frequencies, independent of tip-sample distance, while between these values the momentum can be tuned by the tip-sample distance; (ii) these fixed momenta are identified with the roots k_sn of the Bessel standing-wave condition J_s'(k_sn r0 + phi)=0, corresponding to geometrically quantized polariton modes; and (iii) the edge structure, lacking circular confinement, does not show such quantization but does show distance-dependent momentum tuning. The paper includes tomographic images, extracted dispersion relations, comparison with a second microdisk, simulations of J_s^2 fringe patterns, and vertical sectional analysis that shows a center dip interpreted as destructive interference between adjacent Bessel modes at non-resonant frequencies.","tokens_in":14486,"tokens_out":4194,"duration_ms":42929,"significance":"If the central quantization claim is validated, the work is significant for two reasons. First, it demonstrates a new measurement capability: PF-SNOM provides a 3D near-field cube with explicit tip-sample distance resolution, which standard s-SNOM cannot deliver; this is a genuine methodological advance and is convincingly supported by the data. Second, it reports a potentially important physical observation: circular confinement of hyperbolic phonon polaritons appears to produce discrete, d-independent momentum plateaus, in contrast to the continuous d-dependent momentum in an unstructured edge. The experimental observation of d-independent momenta is reproducible (a second disk shows the same locking behavior), and the vertical sectional analysis provides a novel way to visualize mode superposition. However, the identification of the plateaus with specific Bessel modes (s,n) relies on a fitting procedure in which the phase phi is fitted to the same data and mode numbers are assigned post hoc, with no error bars for the extracted momenta. This circularity weakens the central quantitative claim and must be addressed before the quantization is established.","major_comments":[{"comment":"The assignment of the three d-independent momenta to specific (s,n) Bessel modes is post hoc: the procedure reads q_p from the experimental images, solves Eq. S4 for k_sn, and then assigns (s,n) if the values match, with phi = -0.28π chosen as the best fit on these same data. With only three resonances and no stated uncertainty on q_p (or on the FFT peak position), the agreement (residuals of 0.03-0.06 μm^-1) does not distinguish the quantization hypothesis from a three-point fit. Please provide the uncertainty on each q_p, list all k_sn in the investigated q_p range, and include at least one predictive check with fixed parameters, such as using the phi and r0 from one disk to predict the resonance frequencies of the second disk in Supplementary Fig. S3 or a leave-one-out test on the three observed resonances.","section":"Supplementary Note 1, Table S1"},{"comment":"The anomalous phase phi is introduced as a fitting parameter in the quantization condition J_s'(k_sn r0 + phi)=0, yet in the fringe-pattern simulation (Eq. S5) the phase is omitted 'to achieve optimal agreement between simulated fringes and PF-SNOM fringes.' This dual role is internally inconsistent: if phi is a physical edge correction it should appear in both the eigencondition and the simulated spatial pattern; if it is only a fitting knob, its use in the quantization condition is not independently constrained. Please clarify how phi is determined, whether it is transferred between the two uses, and how the assigned (s,n) change over the plausible range of phi.","section":"Main text Eq. (1) and Supplementary Note 1"},{"comment":"The dispersion curves in Figs. 2e-f are computed using an effective substrate permittivity epsilon_s = 2.05 that is explicitly chosen to match the experimental results, and the justification (a heterogeneous SiO2/Si substrate) does not pin down this value. Consequently the overlay of experimental data and the brown dispersion curve in Figs. 2e-f is not a parameter-free comparison. This does not invalidate the measured plateaus themselves, but the figure should clearly state that epsilon_s is a fitted parameter, and the sensitivity of the comparison to the plausible range of epsilon_s (1.2 to 11.7) should be shown.","section":"Methods, Dispersion relation simulation"},{"comment":"The interpretation of the center dip in the vertical interaction curve at 1425 cm^-1 as destructive interference between the (0,4) and (1,4) Bessel modes (arising from alternating s between adjacent resonances) depends directly on the correctness of the mode assignments in Table S1. Because those assignments are not independently established (see first major comment), this physical interpretation is also not yet fully supported. A quantitative simulation of the expected vertical profile near the center, using the superposition of the assigned modes with distance-dependent coupling weights, would strengthen the argument.","section":"Main text, Fig. 5d and related discussion"}],"minor_comments":[{"comment":"The caption states that at non-resonant conditions 'q_p increases as d increases,' but panel (a) shows the fringe wavelength lambda_f increasing with d, which implies q_p = pi/lambda_f decreases with d; this is also consistent with the main-text statement that larger d gives lower spatial frequency. Please correct this typo.","section":"Supplementary Fig. S4 caption"},{"comment":"The abbreviation 'h-11BN' is used without a definition at first occurrence; please define it (isotopically pure hexagonal boron nitride with boron-11) when it is introduced.","section":"Introduction, abstract"},{"comment":"The text says the AFM is operated with a peak force tapping frequency of 4 kHz and a peak-to-peak sample oscillation amplitude of 300 nm; it would be helpful to specify the cantilever stiffness and the oscillation amplitude at the tip (or at least the typical peak force setpoint) for reproducibility.","section":"Methods, PF-SNOM setup"},{"comment":"In the derivation, the Neumann boundary condition is applied at r = r0 without including the phase phi; the phase is then added to the argument of the Bessel function. This is a nonstandard step that deserves a more explicit justification, even when citing Ref. 29; a reader should be able to see why the same boundary condition is not applied to the phi-shifted argument.","section":"Supplementary Note 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is in scope for a physics-optics venue and the experimental methodology is a clear advance. The main risk is not the PF-SNOM technique but the strength of the quantization claim: the mode assignment is underdetermined with only three resonances, two free parameters (phi and the (s,n) choice), and no error bars. I would be willing to reconsider after the authors provide a predictive test and quantify uncertainties. The second disk is not currently used as a blind validation, which is a missed opportunity. No concerns about citation behavior or novelty disclosure beyond the need to more explicitly separate fitted from predicted elements."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the experiment is genuinely new and the method is worth attention. PF-SNOM gives distance-resolved near-field cubes, and on a 9-micron h-11BN disk they see something s-SNOM cannot: at certain frequencies the extracted in-plane momentum q_p does not move when tip-sample distance d changes, while between those frequencies q_p shifts with d. The edge control shows no locking. A second disk replicates the behavior. That is a solid qualitative result.\n\nThe paper is also honest about its modeling. Supplementary Note 1 lays out the workflow: read q_p from FFT, compute Bessel roots, assign (s,n) if they match. With phi fitted to -0.28 pi and epsilon_s=2.05 chosen to match the same dispersion curves, the agreement is partly a calibration. Three resonances, residuals around 1-2%, and no error bars on q_p mean the specific labels (1,3), (0,4), (1,4) are not strongly constrained. The adjacent-root spacing is bigger than the residuals, which helps, but the freedom in s and phi still leaves room for coincidence. So the momentum-locking observation is robust; the geometric-mode identification is plausible but underdetermined.\n\nThe vertical analysis is interesting too. The d-dependent behavior at 1425 cm-1 between resonances, explained by mode superposition with alternating s parity, is a nice use of the 3D data. It is more qualitative than quantitative.\n\nWho is this for? People building polaritonic resonators and near-field microscopists. It deserves a serious referee because the method and dataset are worth engaging. But it needs revision: error bars or bootstrapped uncertainties for q_p, a genuinely blind test (fix phi and epsilon_s on disk 1 and predict disk 2), and at least a discussion of how many mode assignments are consistent with the data. I would accept it for review with expectation of major revision, and I would bring it to a reading group. I might cite the 3D cube method, though I would wait for the mode assignment to firm up.","headline":"PF-SNOM delivers real 3D near-field data and a plausible momentum-quantization story, but the Bessel-mode assignment is fit to the same three resonances, so treat the specific (s,n) labels as tentative.","tokens_in":15007,"tokens_out":1521,"would_cite":true,"duration_ms":14375,"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 round h-11BN disk locks polariton momentum to discrete values, and the tip-sample distance fine-tunes it between resonances.","keywords":["peak force scattering-type near-field optical microscopy","phonon polaritons","hexagonal boron nitride","momentum quantization","near-field tomography","tip-sample distance","Bessel functions","mid-infrared nanophotonics"],"falsifier":"Measure h-11BN microdisks with deliberately different radii and check whether every resonant frequency and locked momentum obeys $J'_s(k_{sn}r_0-0.28\\pi)=0$ with the same phase shift and no refitting; a single resonance off the predicted ladder would falsify the quantization claim.","tokens_in":14043,"feed_emoji":"🔬","tokens_out":8556,"duration_ms":79506,"temperature":0.7,"pith_summary":"This paper introduces three-dimensional near-field analysis using peak force scattering-type near-field optical microscopy (PF-SNOM), which records the near-field signal at every lateral position and every tip-sample distance at once. Applied to a circular microdisk of isotopically pure hexagonal boron nitride (h-11BN), the method reveals that phonon-polariton momentum takes fixed, discrete values at particular infrared frequencies, matching the standing-wave solutions of a circular resonator. Between those resonances, the same data show that increasing the tip-sample distance continuously lowers the polariton momentum. The paper argues that a reflective edge, lacking circular symmetry, shows no such quantization, confirming that the effect comes from the disk's geometry. The payoff is a complete vertical-plus-lateral map of near fields, which ordinary s-SNOM cannot provide.","feed_headline":"A round h-11BN disk locks polariton momentum to discrete values","feed_subtitle":"PF-SNOM maps the 3D near field and shows tip-sample distance tunes momentum between resonances.","key_machinery":"The central object is the 3D near-field response cube $S_{\\mathrm{NF}}(x,y,d)$, assembled from vertical near-field interaction curves collected by PF-SNOM at every lateral point; it supplies tomographic images at fixed $d$ and sectional images at fixed $x$ or $y$. The quantization argument rests on the standing-wave solution $\\rho = J_s(k_{sn}r+\\phi)(A\\cos(s\\theta)+B\\sin(s\\theta))e^{i\\omega t}$, with resonant spatial frequencies fixed by $J'_s(k_{sn}r_0+\\phi)=0$, where $\\phi=-0.28\\pi$ is the anomalous edge phase shift and $(s,n)$ are mode indices. The tip-sample distance $d$ acts as a controllable momentum filter: small $d$ couples strongly to high spatial frequencies, large $d$ to low spatial frequencies, which explains both the $d$-dependent dispersion and the mode-superposition behavior between resonances.","core_discovery":"The central claim is that the circular boundary of an h-11BN microdisk quantizes the in-plane momentum of hyperbolic phonon polaritons: at resonant frequencies of 1410, 1421, and 1428 $\\mathrm{cm}^{-1}$, the measured polariton momentum $q_p$ stays locked at $2.11$, $2.48$, and $2.84\\ \\mu\\mathrm{m}^{-1}$ regardless of tip-sample distance, and these values coincide with the roots $k_{sn}$ of $J'_s(k_{sn}r_0+\\phi)=0$ with $\\phi=-0.28\\pi$. At non-resonant frequencies the momentum is not locked; instead larger tip-sample separation produces lower spatial frequency, because a more distant tip supports a more loosely confined, lower-momentum field. The vertical dimension of the data cube also distinguishes the two regimes: at resonance the near-field signal grows monotonically as the tip approaches the surface, while between resonances the center of the disk develops a dip caused by destructive interference between two adjacent standing-wave modes. The paper presents the same analysis on an h-11BN flake edge to show that without circular confinement the momentum is continuously tunable and no discrete resonances appear.","pith_inferences":["If the quantization is genuine, a disk of different radius should show resonances at frequencies predicted by the same Bessel condition with unchanged phase shift; a systematic radius series would separate true geometry quantization from coincidence.","The fitted effective substrate permittivity $\\varepsilon_s=2.05$ could be tested by measuring the same disk on substrates with different oxide thickness or on a metal backplane, which should shift the non-resonant dispersion but not the locked resonant momenta.","The vertical dip between resonances implies that conventional tapping-mode s-SNOM, which averages over a range of tip-sample distances, partially washes out the mode-superposition signature; PF-SNOM's distance-resolved data could be used to deconvolve that averaging.","Because the Bessel modes alternate between bright and dark centers with $s=0$ and $s=1$, the vertical interaction curve at the disk center may serve as a fast readout of the mode order without full imaging."],"forward_implications":["PF-SNOM can map the vertical decay of near fields over less than 20 nm, providing a direct check of near-field models rather than a single lateral image.","Circular h-11BN resonators can be used as in situ momentum filters that select polaritons with discrete momenta when integrated into heterostructures.","At the 1420 $\\mathrm{cm}^{-1}$ resonance, the disk center concentrates near-field intensity into a roughly 500 nm hotspot, a location usable for mid-infrared chemical sensing or nonlinear mixing.","The same 3D tomography should reveal geometric resonances in other polaritonic and plasmonic microstructures, since PF-SNOM and s-SNOM share the same scattering mechanism."],"supporting_citations":[{"why":"Establishes the PF-SNOM method for converting time-domain scattering and cantilever-deflection data into vertical near-field interaction curves and 3D response cubes.","marker":"17"},{"why":"Supplies the 2D wave-equation derivation and Bessel standing-wave solutions used to predict quantized spatial frequencies of the circular resonator.","marker":"29"},{"why":"Provides the hyperbolic-phonon-polariton dispersion relation for the h-BN/SiO2 structure that the paper compares against the measured $q_p$ values.","marker":"3"},{"why":"Underpins the choice of isotopically pure h-11BN, whose reduced damping gives long polariton propagation and sharp resonances.","marker":"24"},{"why":"Is the reference the Methods cites for the h-11BN dielectric parameters used in the calculated dispersion curves.","marker":"22"},{"why":"Supplies the isotopically pure h-11BN material used to fabricate the microdisk and flake samples.","marker":"18"}],"fun_headline_variants":["Polariton momentum locks to discrete values in a round h-11BN disk","3D near-field maps reveal quantized polariton momentum in h-11BN","Tip distance tunes polariton momentum between resonances in h-11BN","Circular h-11BN disk quantizes phonon polariton momentum","PF-SNOM images 3D near fields: momentum quantization in h-11BN"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The momentum-quantization claim rests on fitting the edge phase shift to $-0.28\\pi$ and then assigning whole-number mode labels to the observed resonances; if those fitting choices are wrong, the locked momenta could be a coincidence rather than true geometric resonances.","fun_headline_variants_meta":{"raw":{"variants":["Polariton momentum locks to discrete values in a round h-11BN disk","3D near-field maps reveal quantized polariton momentum in h-11BN","Tip distance tunes polariton momentum between resonances in h-11BN","Circular h-11BN disk quantizes phonon polariton momentum","PF-SNOM images 3D near fields: momentum quantization in h-11BN"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000223,"raw_usage":{"total_tokens":1492,"prompt_tokens":1012,"completion_tokens":480,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":628,"completion_tokens_details":{"reasoning_tokens":376}},"tokens_in":628,"tokens_out":480,"duration_ms":4536,"temperature":1.0,"reasoning_tokens":376,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:19:12.879557+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure h-11BN microdisks with deliberately different radii and check whether every resonant frequency and locked momentum obeys $J'_s(k_{sn}r_0-0.28\\pi)=0$ with the same phase shift and no refitting; a single resonance off the predicted ladder would falsify the quantization claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the hyperbolic-phonon-polariton dispersion relation for the h-BN/SiO2 structure that the paper compares against the measured $q_p$ values."}],"review_version":1}