{"id":"73297b72-cabb-405c-b07c-70c55ffa3b4b","arxiv_id":"2412.07903","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using Madelung hydrodynamics, the authors derive analytic nonlocal expressions for magnetoplasmons, magneto-optical conductivity, plasmon-enhanced second-harmonic generation, and a velocity-dependent renormalization of the SPP dispersion.","lead":"This paper applies Madelung's hydrodynamic formulation to derive nonlocal corrections in 2D electron gas plasmons, magneto-optical conductivity, second-harmonic generation, and the self-modulation of surface-plasmon polaritons. It provides closed-form dispersion formulas that could simplify design calculations in 2D-material nonlinear plasmonics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Self-modulation claim operates in a regime where the two-harmonic truncation breaks down: at |vω|=20vF, where the effect becomes significant, Eq. (81) gives nω/n0 not small, so Eq. (82) is unvalidated.","rationale":"The reader's conditional verdict is appropriate. The linear-response and SHG cross-checks against Boltzmann theory are strong evidence that the hydrodynamic framework is sound. The self-modulation section, however, is the only place where the paper advances a new, unaided prediction. The load-bearing assumption is not the algebra leading to Eq. (82) (we spot-checked the small-amplitude limit for β=0 and found it consistent), but the claim that the result applies at amplitudes where the effect is visible. Since the paper itself concedes the breakdown at high amplitude and provides no bridging calculation, the conclusion that nonlinearity amplifies nonlocal corrections is not established in a controlled regime. This does not invalidate the framework; it makes the self-modulation result conditional on either a nonperturbative treatment or a demonstration that physically attainable SPP amplitudes are small enough for the truncation to hold. We therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":17226,"tokens_out":34515,"duration_ms":283256,"concrete_test":"Extend the harmonic ansatz of Sec. IIIE to include the third harmonic (n3ω, v3ω) and solve the resulting algebraic system for the same parameters as Fig. 4. If the frequency shift δω = ω(|vω|)-ω0 at |vω|=20vF changes by more than 10% when the 3ω terms are included, the two-harmonic truncation is not converged and Eq. (82) cannot support the claim. A complementary check: evaluate nω/n0 from Eq. (81) at the plotted k range and |vω|=20vF; if it exceeds ~0.3, the perturbative expansion is invalid by the paper's own criterion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. IIIE derives Eq. (82) by truncating the hydrodynamic variables at the second harmonic (nω, n2ω, vω, v2ω) and treating |vω| as a free, real amplitude. The central claim, that nonlinearity renormalizes the SPP dispersion and acts like a nonlocal contribution similar to Fermi pressure, is evaluated in Fig. 4 for |vω| = 20 and 40 v_F. However, the same figure (panel B) and Eq. (81) show that nω/n0 is not a small parameter in this range: at high k and |vω|=40v_F the density fluctuation approaches n0, and the authors concede the perturbative approach is no longer valid. At |vω|=20v_F, where they claim a considerable difference in the dispersion, nω/n0 is already sizable, so the neglect of higher harmonics and of nω^2 terms in the harmonic balance is uncontrolled. No independent estimate connects vω to a physical SPP excitation amplitude, and no estimate of the neglected 3ω terms is given. Thus the regime that supports the headline conclusion lies outside the domain of validity of the derivation, and Eq. (82) overstates the renormalized spectrum.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the Madelung (quantum) hydrodynamic formalism to several problems in two-dimensional plasmonics and nonlinear optics. It derives the magnetoplasmon spectrum and magneto-optical conductivity with nonlocal corrections from the Fermi pressure, computes second-harmonic generation both in the bare and in the plasmon-assisted self-consistent case, and presents a new self-modulation analysis of surface-plasmon polaritons. The central new claim is that nonlinearity renormalizes the SPP dispersion and acts as an additional source of nonlocal behavior, expressed analytically in Eq. (82) and illustrated in Fig. 4. The manuscript also reviews the derivation of the Madelung equations from the Wigner function and verifies several results against independent Boltzmann kinetic calculations.","tokens_in":17490,"tokens_out":13892,"duration_ms":106047,"significance":"If the central claim holds, the paper offers a simple hydrodynamic derivation of nonlinear nonlocal corrections in 2D materials, with a compact analytic formula for the renormalized SPP dispersion. The cross-checks against Boltzmann kinetic results in Refs. 41, 44, and 48, and the independent derivation of the nonlinear current via the polarization formalism in Sec. IIIC.1, are clear strengths. The self-modulation result is the main novelty, but its demonstration relies on a perturbative regime whose validity is not established for the parameter values used to display the effect.","major_comments":[{"comment":"The matrix in Eq. (14) is inconsistent with the preceding equations. Substituting the continuity solution (12a) and the self-consistent potential (13) into the linearized Euler equation (12b) produces a matrix that contains an additional contribution proportional to q^2 n0/(2 m epsilon0 |k|) in the entries, reflecting the Coulomb interaction. The matrix as printed contains only the Fermi-pressure term and the cyclotron term; its determinant gives omega^2 = beta^2 k^2 + omega_c^2, without the ak term that appears in Eq. (15). The derivation should be corrected so that the matrix reproduces the claimed dispersion relation.","section":"Sec. IIIA, Eq. (14)"},{"comment":"The perturbative expansion underlying Eq. (82) is not controlled in the regime where the self-modulation effect is claimed to be significant. The paper acknowledges at the end of Sec. IIIE that for |v_omega| = 40 v_F the density fluctuation approaches n0 and the approach is no longer valid. However, Fig. 4B shows that already for |v_omega| = 20 v_F the density fluctuation is a sizable fraction of n0 at moderate wavevectors, so the neglect of higher harmonics and of higher-order nonlinearities is not justified. No estimate of the small parameter or of the magnitude of the neglected 3-omega terms is provided. Consequently, Eq. (82) is not validated in the regime that supports the headline conclusion that nonlinearity renormalizes the SPP dispersion; the claim should either be restricted to the |v_omega| << v_F regime (where the effect is negligible) or be supported by a controlled expansion with a quantitative error bound.","section":"Sec. IIIE, Eqs. (81)-(83) and Fig. 4"},{"comment":"The assumption that the harmonic amplitudes n_omega, n_2omega, v_omega, and v_2omega are real is introduced without justification. For a propagating SPP described by a factor exp(i k x - i omega t), the complex amplitudes carry phase information; setting them real imposes a specific phase relation that is not generally valid. This assumption is essential for reducing Eqs. (73) to (74), and the authors should either state the physical conditions under which it holds or redo the harmonic balance with complex amplitudes. Without this, the derivation of Eq. (82) rests on an unverified ansatz.","section":"Sec. IIIE, Eqs. (74a)-(74d)"}],"minor_comments":[{"comment":"The continuity equation (29a) mixes perturbation orders: it combines the first-order terms partial_t n1 + n0 div v1 with the second-order terms div(n1 v1) + n0 div v2. It would be clearer to write separate equations at each order and then collect harmonics.","section":"Sec. IIIC, Eq. (29a)"},{"comment":"Eq. (72) contains a typographical error: E_{l omega}^2 should likely be E_{l omega} (or a different notation), since the electric field is not squared.","section":"Sec. IIIE, Eq. (72)"},{"comment":"The notation for the density harmonics is inconsistent: n1 is used for the total first-order density (containing both omega and 2omega components), while n2 in Eq. (36) denotes the second-order density at 2omega. Standardizing the notation would improve readability.","section":"Secs. IIIC and IIIE"},{"comment":"The paper would benefit from a brief discussion of the physical meaning of the amplitude v_omega and how it relates to an external drive or to the intensity of the SPP, since Eq. (82) is expressed in terms of this input parameter.","section":"Sec. IIIE, Eq. (82)"}],"recommendation":"major_revision","confidential_remarks":"The magnetoplasmon matrix error in Eq. (14) is readily correctable and does not invalidate the known final dispersion, but it must be fixed in revision. The main obstacle is the self-modulation section: the central claim is only demonstrated in a regime where the perturbative expansion is uncontrolled. The authors should provide a systematic expansion with error estimates or clearly restrict the claim to the regime of validity. The paper is otherwise well written and the cross-checks with kinetic calculations are valuable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a useful hydrodynamic review of 2D plasmonics with two new pieces bolted on. The linear-response and SHG parts are solid and match Boltzmann kinetic results. The nonlocal SHG formula in Eq. (69) looks right, and the polarization cross-check in Sec. IIIC.1 is a good sign. The self-modulation formula in Eq. (82) is the advertised new result, but it is not under control.\n\nWhat does the paper actually do well? It gives a compact derivation of the Madelung equations from the Wigner function, then applies them to magnetoplasmons, magneto-optical conductivity, SHG, plasmon-assisted SHG, and SPP self-modulation. The agreements with Refs. 41, 44, and 48 are convincing and no fitted parameters are involved. That's genuine evidence the framework is working.\n\nThe soft spots. First, Eq. (14) as printed cannot produce Eq. (15) — the matrix has no Coulomb potential term, so the determinant will not give the a k term. This looks like a typo, but it's a confusing one. Second, Eq. (29a) mixes first- and second-order terms in the continuity equation. You can sort it out by collecting frequencies, but as written it's not a clean perturbative expansion. Third, the self-modulation result depends on an arbitrary velocity amplitude |v_ω| with no physical estimate, and the two-harmonic truncation fails exactly where the effect is claimed. The authors admit the breakdown at |v_ω| = 40 v_F, but panel 4B shows n_ω/n_0 is already sizable at 20 v_F, where they claim a considerable dispersion change. So Eq. (82) is an interesting suggestion, not a demonstrated result.\n\nWho is this for? Someone who wants a readable hydrodynamic derivation of known 2D plasmon results and a formula for nonlocal SHG. The self-modulation section is not ready for quantitative use. The paper deserves a serious referee because the core framework is sound and the flaws are fixable, but it needs revision: fix the matrix in Eq. (14), present the perturbation theory cleanly, and either connect v_ω to a physical excitation amplitude or state clearly that Eq. (82) is a heuristic extrapolation.\n\nMy recommendation: send it to a competent referee with specific instructions to check Eq. (14) and the validity of Eq. (82). If those are addressed, it's a citeable reference. Right now I wouldn't base any of my own work on the self-modulation claim.","headline":"A useful hydrodynamic review with a couple of questionable new results; the core framework is sound but the self-modulation claim needs more work.","tokens_in":18029,"tokens_out":11026,"would_cite":false,"duration_ms":91638,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["73.20.Mf","42.65.Ky"],"model":"deepseek-v4-flash","headline":"Nonlinearity renormalizes surface-plasmon dispersion in 2D electron gases and acts as an extra nonlocal pressure.","keywords":["Madelung transformation","quantum hydrodynamics","surface plasmon polaritons","nonlocal plasmonics","second-harmonic generation","nonlinear dispersion renormalization","two-dimensional electron gas","magneto-optical conductivity"],"falsifier":"Measure the surface-plasmon dispersion of a gated two-dimensional electron gas by near-field or attenuated-total-reflection spectroscopy while increasing the driving field intensity at fixed wavevector; if the resonance frequency does not shift upward according to Eq. (82), or shifts downward, the central prediction is wrong. A complementary check is a fully nonlinear numerical solution of the hydrodynamic equations without truncating at second order, which should reproduce the analytic spectrum only where the perturbative assumption holds.","tokens_in":17045,"feed_emoji":"⚡","tokens_out":11449,"duration_ms":96060,"temperature":0.7,"pith_summary":"This paper establishes that nonlinearity renormalizes the dispersion of surface-plasmon polaritons (coupled light–charge oscillations at an interface) in a two-dimensional electron gas, acting as an extra nonlocal contribution alongside the Fermi pressure. Using the Madelung hydrodynamic picture, the authors derive an analytic expression, Eq. (82), that ties the nonlinear velocity amplitude $v_\\omega$ to a positive frequency shift of the plasmon spectrum. The same framework reproduces known second-harmonic-generation and magneto-optical conductivity results, and shows that plasmon excitation resonantly boosts second-harmonic emission. The payoff is a simple analytic route to nonlinear nonlocal optical response in 2D materials, where kinetic-equation approaches are heavier.","feed_headline":"Nonlinearity shifts 2D surface-plasmon frequencies upward","feed_subtitle":"A hydrodynamic formula ties field strength to a blue shift of 2D plasmon resonances.","key_machinery":"The central object is the Madelung transformation, $\\Psi = \\sqrt{n}\\,e^{iS}$, which converts the Schrödinger equation into continuity and Euler-type equations for the electron density $n$ and velocity field $v$. Supplemented by the Fermi pressure and Poisson's equation, this hydrodynamic system is the machinery that carries every result: harmonic expansions of $n$ and $v$ at $\\omega$ and $2\\omega$ produce algebraic equations whose solution yields the magnetoplasmon spectrum, the nonlocal conductivity tensor, the second-harmonic current, and finally the renormalized dispersion relation in Eq. (82).","core_discovery":"The central claim is that solving the Madelung equations to second order in the electron velocity field, together with Poisson's equation, gives a self-modulated surface-plasmon-polariton dispersion of the form $\\omega^2 = q^2 n_0 k/(2m\\epsilon_0) + \\beta^2 k^2 + k v_\\omega^2 F(k,\\omega)$, so that as the velocity amplitude $|v_\\omega|$ grows the plasmon frequency increases. The paper argues that nonlinear effects play a role analogous to the Fermi pressure, supplying another nonlocal contribution to the optical response. For the same framework, the paper derives nonlocal magnetoplasmon spectra and magneto-optical conductivities, and shows that second-harmonic generation is resonantly enhanced at the fundamental and second-harmonic plasmon poles, with the hydrodynamic results matching earlier Boltzmann kinetic calculations in the appropriate limits.","pith_inferences":["The intensity-dependent blue shift predicted by Eq. (82) suggests a concrete experimental signature: in graphene, tuning the pump intensity should shift the near-field plasmon resonance peak, and the shift should grow with wavevector $k$.","Because the perturbative expansion breaks down near $|v_\\omega| \\approx 40 v_F$, strongly driven systems likely require retaining higher harmonics or a fully nonlinear treatment, a regime the paper leaves open.","The Madelung framework could be extended to include relaxation and retardation, predicting how the nonlinear blue shift competes with damping in realistic surface-plasmon waveguides; the paper treats the renormalized spectrum in the lossless case.","If the renormalized dispersion is measured, the slope of the frequency shift versus intensity would provide a direct estimate of the nonlinear velocity amplitude $v_\\omega$ in the two-dimensional electron gas."],"forward_implications":["At high driving amplitudes the surface-plasmon frequency rises with $|v_\\omega|$, so intense plasmon fields are predicted to blue-shift the resonant response of 2D electron gases.","Because the nonlinear term enters multiplied by the wavevector $k$ in Eq. (82), nonlinearity is a nonlocal correction and grows in importance at shorter wavelengths.","Second-harmonic generation from a 2D electron gas is resonantly enhanced when either $\\omega$ or $2\\omega$ matches a plasmon pole of the nonlocal dielectric function, and vanishes at normal incidence by symmetry.","The hydrodynamic results for the nonlinear current and the magneto-optical conductivity coincide with previous Boltzmann kinetic calculations in the collisionless limit, so the simpler analytic expressions can be used for modeling."],"supporting_citations":[{"why":"Supplies the original Madelung transformation used throughout to rewrite the Schrödinger equation as hydrodynamic equations.","marker":"[19]"},{"why":"Supplies the many-body and Wigner-function derivation that incorporates statistical pressure into the Madelung framework.","marker":"[32]"},{"why":"Provides the magnetoplasmon spectrum that Eq. (16) reproduces with nonlocal corrections.","marker":"[34]"},{"why":"Provides the Boltzmann kinetic calculation against which the nonlocal magneto-optical conductivity is compared.","marker":"[41]"},{"why":"Provides the Boltzmann kinetic result for the nonlinear current that the hydrodynamic second-harmonic expression matches in the collisionless limit.","marker":"[44]"},{"why":"Supplies the Boltzmann kinetic calculation of plasmon-assisted second-harmonic generation that Eq. (68) agrees with.","marker":"[48]"},{"why":"Establishes that surface-plasmon polaritons are intrinsically nonlinear through self-modulation, motivating the renormalized dispersion calculation.","marker":"[49]"}],"fun_headline_variants":["Nonlinearity blue-shifts 2D plasmon resonances","Madelung hydrodynamics: Nonlinearity renormalizes 2D plasmons","2D plasmon frequencies rise with nonlinearity, model shows","Nonlocal and nonlinear effects shape 2D plasmon spectra","Hydrodynamic model links field amplitude to 2D plasmon shift"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the nonlinear velocity amplitude stays small enough that the first and second harmonics are only weak corrections; the paper itself notes that once $|v_\\omega|$ approaches about $40 v_F$, the density fluctuation reaches the full background density and the perturbative expansion is no longer valid.","fun_headline_variants_meta":{"raw":{"variants":["Nonlinearity blue-shifts 2D plasmon resonances","Madelung hydrodynamics: Nonlinearity renormalizes 2D plasmons","2D plasmon frequencies rise with nonlinearity, model shows","Nonlocal and nonlinear effects shape 2D plasmon spectra","Hydrodynamic model links field amplitude to 2D plasmon shift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000196,"raw_usage":{"total_tokens":1320,"prompt_tokens":865,"completion_tokens":455,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":481,"completion_tokens_details":{"reasoning_tokens":367}},"tokens_in":481,"tokens_out":455,"duration_ms":4388,"temperature":1.0,"reasoning_tokens":367,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:26:21.480099+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the surface-plasmon dispersion of a gated two-dimensional electron gas by near-field or attenuated-total-reflection spectroscopy while increasing the driving field intensity at fixed wavevector; if the resonance frequency does not shift upward according to Eq. (82), or shifts downward, the central prediction is wrong. A complementary check is a fully nonlinear numerical solution of the hydrodynamic equations without truncating at second order, which should reproduce the analytic spectrum only where the perturbative assumption holds.","supporting_citations":[{"cited_title":"Born, Zeitschrift für Physik37, 863 (1926)","cited_arxiv_id":null,"evidence_quote":"Supplies the original Madelung transformation used throughout to rewrite the Schrödinger equation as hydrodynamic equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the many-body and Wigner-function derivation that incorporates statistical pressure into the Madelung framework."},{"cited_title":"Locarno and D","cited_arxiv_id":null,"evidence_quote":"Provides the magnetoplasmon spectrum that Eq. (16) reproduces with nonlocal corrections."},{"cited_title":"Wegner, D.-N","cited_arxiv_id":null,"evidence_quote":"Provides the Boltzmann kinetic calculation against which the nonlocal magneto-optical conductivity is compared."},{"cited_title":"Ferreira, J","cited_arxiv_id":null,"evidence_quote":"Provides the Boltzmann kinetic result for the nonlinear current that the hydrodynamic second-harmonic expression matches in the collisionless limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Boltzmann kinetic calculation of plasmon-assisted second-harmonic generation that Eq. (68) agrees with."},{"cited_title":"Scalora, M","cited_arxiv_id":null,"evidence_quote":"Establishes that surface-plasmon polaritons are intrinsically nonlinear through self-modulation, motivating the renormalized dispersion calculation."}],"review_version":1}