{"id":"ff1ebc80-db8e-479c-8951-6ce5de6be036","arxiv_id":"2511.21059","paper_version":3,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new universal evolution equation, with an Integration Hamiltonian and band–wave correspondence, models nonlinear frequency combs under arbitrary electro-optic modulation in strong-coupling regimes where the standard Lugiato–Lefever equation fails.","lead":"The authors propose a new general equation to model light pulses in electro-optic cavities, claiming it works for both weak and strong microwave modulation, and they test it on fiber and chip-based lithium-niobate setups. The framework is aimed at giving engineers programmable control over comb light sources used in precision measurement and communications.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The band–wave correspondence rests on the excitation condition Ĥ_ig=0, which is asserted rather than derived; if this condition is not the correct resonance condition, the central programmable-spectral-control claim is not established.","rationale":"The reader's weakest assumption identifies the same load-bearing point: the BWC's excitation condition is deferred to the missing SI and is therefore unverifiable. I agree and sharpen it: the condition Ĥ_ig=0 is not only underived but could be wrong, because it is an operator-level statement converted into a c-number resonance condition. The entire programmable spectral control claim depends on this condition. If the condition is correct, the BWC is a genuine consequence of the GEE; if not, the central claim collapses to an additional postulate. The paper's experimental Fig. 2 provides some independent support for the Integration Hamiltonian via measured spectral shapes and coupling strengths, and Fig. 3 claims experimental confirmation of BWC, but the quantitative comparison is not testable from the submitted text. The concrete test I propose—solving the discrete-time Heisenberg equation and comparing peak locations to the BWC formula—would settle whether the excitation condition is a real consequence of the model or an asserted ansatz. Since this is an addressable, checkable concern rather than a demonstrated contradiction, I do not move the verdict away from CONDITIONAL; the reader's conditional verdict already captures this risk, so I mark the verdict as UNCHANGED.","tokens_in":9914,"tokens_out":10276,"duration_ms":112756,"concrete_test":"Numerically solve the discrete-time Heisenberg equation (5) with the pump term √κ_e a_in δ_{0,μ} for an asymmetric triangular modulation V(t), using the Hamiltonian of Eq. (4) (or equivalently the round-trip map of the simplified GEE). Locate the spectral peaks of the stationary state and compare their mode indices/frequencies with the BWC prediction E(φ)/ℏ = ω0 − π f_R V(t=−φ/ω_R)/Vπ − nω_R. If the excited modes do not coincide with the modes satisfying Ĥ_ig=0, the excitation condition is invalid and the BWC must be revised. Alternatively, analytically rederive the resonance condition from the stationary-phase points of the round-trip transfer function T(φ)=i f_R(e^{∫(−iΔ+iΩ)dt}−1) without invoking the Ĥ=0 ansatz, and check whether it yields the same band formula.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of programmable spectral control via the band–wave correspondence (BWC) hinges on the statement in the 'Band-wave correspondence' section: 'The excitation condition of mode φ is Ĥ_ig(φ,t)=0', leading to E(φ)/ℏ = ω0 − π f_R V(t=−φ/ω_R)/Vπ − nω_R. This is not derived in the submitted text; the derivation is relegated to the missing SI. The issue is not merely missing derivations: the condition may be an arbitrary ansatz. In Eq. (3), Ĥ_ig is an operator; setting it to zero is a c-number condition on the round-trip phase integral. It is not shown that the modes satisfying Ĥ_ig=0 are exactly the modes resonantly excited by the periodic pump in the discrete-time Heisenberg equation (5). If the excitation condition is chosen ad hoc, the BWC, the triangular-wave sideband engineering in Figs. 4d–4e, and the claimed 'universality' of the framework lose their quantitative foundation. The reader's missing-SI concern is therefore not only a verifiability deficit: the excitation condition could be incorrect. A direct check is needed to decide whether the correspondence is a consequence of the GEE or an added postulate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a general evolution equation (GEE) for nonlinear frequency combs under arbitrary electro-optic modulation, intended to replace the mean-field Lugiato–Lefever equation in strong-coupling regimes. From the GEE the authors derive a discrete-time Integration Hamiltonian, an energy-band picture, and a 'band–wave correspondence' that maps modulation waveforms to synthetic band structures. They report experimental validation in fiber and thin-film lithium niobate platforms, use triangular-wave modulation to demonstrate directional coupling and single-sideband comb generation, and extend the framework to Kerr soliton dynamics with a 'soliton band-drifting' theory and addressing scheme. The central claims are that the GEE is universal, that the Integration Hamiltonian provides a rigorous frequency-domain formalism, and that the band–wave correspondence enables programmable spectral control.","tokens_in":10248,"tokens_out":3464,"duration_ms":40145,"significance":"If the derivations are correct, this would be a substantial advance: it offers a unified framework for strong-coupling electro-optic and hybrid EO–Kerr combs, connects frequency-domain EO comb physics to synthetic dimensions, and proposes a concrete mechanism for waveform-programmable spectral control. The experimental work is a real asset: two complementary platforms, triangular-wave sideband engineering, and single-sideband generation are demonstrated, and the soliton addressing simulations produce falsifiable predictions. However, the manuscript as submitted defers all key derivations to a missing SI and asserts the central excitation condition, so the validity of the main claims cannot currently be assessed from the text alone. The paper's significance is therefore conditional on the supplied derivation and verification.","major_comments":[{"comment":"The central reduction from the GEE to the discrete-time Heisenberg equation and the Integration Hamiltonian is asserted, with details deferred to the SI, which is not included. In particular, the transformation from the operator Eq. (3) to the longitudinal-mode basis Eq. (4), the Fourier-series expansion of T(φ), and the neglect of dispersion and Kerr terms are not shown. Since the spectra in Fig. 2b and the coupling strengths |W_m| in Fig. 2c are computed from these equations, the quantitative experimental agreement cannot be checked. Please include the full derivation and state all assumptions.","section":"Formalism of Integration Hamiltonian, Eqs. (2)–(5)"},{"comment":"The statement 'The excitation condition of mode φ is Ĥ_ig(φ,t)=0' is asserted without derivation. Ĥ_ig is an operator, and setting it to zero is a c-number condition on the round-trip phase integral; it is not shown to be equivalent to the resonance condition of the discrete-time Heisenberg equation (5) or to the response peak of the GEE. The band–wave correspondence E(φ)/ℏ = ω0 − π f_R V(t=−φ/ω_R)/Vπ − nω_R is the load-bearing step for the programmable-spectral-control claim and for the band-overlap interpretation in Figs. 3 and 4. Either derive this condition from the GEE or state it explicitly as an additional postulate, and validate it numerically (e.g., weakly probe each azimuthal mode in the GEE and compare the resonant frequency with the predicted band).","section":"Band-wave correspondence"},{"comment":"The experimental sections do not provide the full parameters needed to reproduce or quantitatively assess the claims: fiber cavity length/finesse, TFLN quality factors (loaded and intrinsic), coupling rates, fiber coupling efficiency, detuning Δ, pump power, dispersion D2, and modulation waveform parameters are not tabulated. Agreement between theory and experiment is described qualitatively as 'agree well' or 'consistent', without error bars, residuals, or a fitting metric. Given the universality claim, please provide a full parameter table and a quantitative comparison for the spectra and coupling strengths.","section":"Experimental validation, Figs. 2–4"},{"comment":"The derivation of the soliton drifting velocity v_d = −d2 k/3 and the soliton-addressing condition is deferred to the 'Supplement'/'SI', which is missing. The phase diagrams and the claim of EO pulse–Kerr soliton co-excitation rely on these results. Please provide the perturbative Lagrangian derivation, the simulation parameters (Q_i, D2, pumping and sweeping rates, noise model), and the definition of the 'probability of successful soliton addressing' used in Fig. 5c.","section":"Soliton band-drifting, Eq. (7) and Fig. 5"}],"minor_comments":[{"comment":"Several symbols are used before definition or are ambiguous: t_R, α, V_local(φ,t) = V(t−φ/ω_R) vs. V(t) later in the BWC expression, and the distinction between Ω0 and Ω(φ,t). Please define all symbols at first use.","section":"Eq. (1)"},{"comment":"The operator ordering in the integrand a†(φ,t)a(φ,t) and the transition from the integral over φ to the Bessel-function series should be shown explicitly; also check the m=0 term sign convention.","section":"Eq. (6)"},{"comment":"The phrase 'multiple ω_R-separated bands precisely replicate the shape of the modulation wave' is descriptive; please clarify what is plotted in Fig. 3 (axes, units, and which experimental waveform was used).","section":"Section 'Band-wave correspondence'"},{"comment":"Figure axes are often unlabeled or lack units (e.g., Fig. 3 band structure, Fig. 5 phase diagrams). Figure 5e should explain the red horizontal lines and the meaning of each panel in the caption.","section":"Figures"},{"comment":"Reference [34] is a Tidy3D simulation notebook; it should be cited as a software/online resource with proper author and access information. References to 'SI' and 'Supplement' must be included in the submission for review.","section":"References"},{"comment":"There are minor typographical inconsistencies in subscripts and Hamiltonian notation (e.g., Ĥ_ig vs. H_ig, Ω vs. Ω(φ,t)). These should be harmonized.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has the potential to be an important contribution, and the experimental demonstrations are valuable. However, the central derivation is deferred to a missing SI, and the excitation condition behind the band–wave correspondence is asserted rather than derived. The stress-test concern is real: unless the excitation condition is shown to follow from the GEE or is explicitly tested as a postulate, the programmable-spectral-control claim is not established. I recommend major revision with a request for the full derivation and a numerical falsification test, rather than rejection, because these issues are addressable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this paper is worth reading seriously and deserves a proper referee. The GEE is a genuine attempt to move beyond the mean-field LLE for strong-coupling electro-optic combs, and the Integration Hamiltonian plus band–wave correspondence give a plausible frequency-domain picture. If it holds up, it would unify a chunk of recent strong-coupling EO work and give a practical tool for programming combs. But the load-bearing parts of the derivation are in the SI, which is not part of this submission, and the excitation condition that anchors the BWC is asserted rather than derived. That is the core concern.\n\nWhat is genuinely new: the GEE with a periodic delta pump is a real departure from the LLE, and the paper argues convincingly that the mean-field approximation artificially reinjects energy. The reduction to the Integration Hamiltonian, with coupling coefficients as Fourier coefficients of the transmission function, is elegant and testable. The BWC is a nice way to tie modulation waveforms to band structures, and the triangular-wave sideband engineering is a clever application. The two experiments, fiber and TFLN, give some direct support for the coupling spectra and the band overlap features.\n\nSoft spots: the main text compresses Eqs. (2)–(5) and the excitation condition to 'see SI', so as submitted the central claims are unverifiable. The stress-test raises a fair question: why is the excitation condition Ĥ_ig=0? It is not obvious that this is the correct resonance condition rather than an ansatz. That said, the experimental BWC demonstration in Fig. 3 suggests the condition works in practice for the tested cases, so I would call it an addressable gap rather than a fatal flaw. The experimental agreement is described qualitatively, with no error bars or full parameters (finesse, Q, fiber length), so the reader cannot check the claimed comparisons. The 'universal' framing is perhaps over-strong given the limited experimental test space, though the theory itself is general.\n\nOn balance, the paper shows clear thinking and honest engagement with the literature. The missing SI and the asserted excitation condition are addressable rather than fatal. If the SI actually contains the derivation and the condition emerges from the GEE, this is a solid contribution. The referee should require the SI, a derivation of the excitation condition, and quantitative experiment details before publication.\n\nRecommendation: send to peer review. I would not cite it yet, but I would want to see the SI.","headline":"A serious and genuinely new framework for strong-coupling EO combs, but the load-bearing derivation is in the missing SI and the excitation condition is asserted rather than proved.","tokens_in":10734,"tokens_out":2375,"would_cite":false,"duration_ms":25531,"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 evolution equation with a genuinely periodic pump supersedes the mean-field model for electro-optic combs and reduces to a discrete-time Hamiltonian whose bands mirror the drive waveform.","keywords":["nonlinear frequency combs","electro-optic modulation","strong-coupling regime","Lugiato-Lefever equation","Integration Hamiltonian","synthetic frequency dimension","band-wave correspondence","Kerr solitons"],"falsifier":"Simulate the printed GEE and Ĥig side by side: numerically integrate the GEE for a triangular drive with Ω0 between one and two free spectral ranges, and compare the steady-state spectrum with the solution of the discrete-time Heisenberg equation from Eqs. (4)–(5); any mismatch refutes the reduction. Independently, measure the band structure in a fiber EO cavity by sweeping detuning Δ at fixed drive and check that transmission resonances lie on E(φ)/ħ = ω0 − πfR V(−φ/ωR)/Vπ − nωR for the appropriate n; the paper asserts this universal relation, so a measured band that is not a scaled mirror im","tokens_in":9821,"feed_emoji":"⚡","tokens_out":19579,"duration_ms":180459,"temperature":0.7,"pith_summary":"This paper tries to establish a universal starting point, the general evolution equation (GEE), for nonlinear frequency combs in cavities driven by arbitrary electro-optic modulation — including the strong-coupling regime where modulation strength exceeds the cavity's free spectral range and the standard mean-field Lugiato–Lefever equation (LLE) breaks down. The load-bearing move is the pump: the GEE injects the periodic round-trip pump, whereas the LLE's mean-field treatment, the paper argues, artificially reinjects energy and imposes a wrong repetition rate. In the electro-optic-only limit the GEE reduces to a discrete-time Integration Hamiltonian, and the condition that singles out excited modes yields the band–wave correspondence: the synthetic band structure is a scaled, mirrored copy of the applied microwave waveform. The paper reports experimental confirmation on both a fiber cavity and an integrated thin-film lithium niobate chip, and uses the formalism to demonstrate triangular-wave sideband engineering, single-sideband comb control, and a soliton band-drifting mechanism that makes Kerr-soliton addressing deterministic. A sympathetic reader would care because the framework turns a microwave waveform into a designed comb spectrum and, if correct, supersedes the LLE as the foundation for strong-coupling electro-optic and hybrid EO–Kerr comb modeling.","feed_headline":"One equation replaces the flawed model for strong-coupling combs","feed_subtitle":"The synthetic band structure mirrors the microwave drive, so comb spectra become programmable.","key_machinery":"The central object is the Integration Hamiltonian Ĥig, a discrete-time operator obtained by integrating the GEE over one cavity round trip and transforming to the longitudinal-mode basis; its couplings Wm are Fourier coefficients of a round-trip transmission function, reducing for sinusoidal drive to Bessel-function couplings. It converts a continuously pumped cavity into a stationary band problem and makes the synthetic frequency dimension explicit without a rotating-wave approximation. Its companion is the band–wave correspondence E(φ)/ħ = ω0 − πfR V(−φ/ωR)/Vπ − nωR, mapping the local modulation voltage onto ωR-separated synthetic bands.","core_discovery":"The paper claims the mean-field Lugiato–Lefever equation fails for electro-optic combs once modulation strength Ω0 exceeds the free spectral range ωR, since its pump reinjects energy. The replacement general evolution equation (GEE) uses a genuinely round-trip-periodic pump and reduces, in the pure-EO limit, to a discrete-time Heisenberg equation with an Integration Hamiltonian Ĥig whose couplings Wm are Fourier coefficients of a round-trip transmission function. The excitation condition Ĥig = 0 yields the band–wave correspondence E(φ)/ħ = ω0 − πfR V(−φ/ωR)/Vπ − nωR: synthetic bands mirror the drive waveform. Fiber and chip experiments match, and the formalism yields single-sideband control","pith_inferences":["Read as an inverse-design rule, the band–wave correspondence suggests that any target band shape can be reached by synthesizing V(t) from the inverted relation V(−φ/ωR) = (Vπ/πfR)(ω0 − nωR − E(φ)/ħ); the paper demonstrates one corner of this space (triangular waves) but does not test arbitrary-waveform synthesis.","Because the GEE's pump is strictly round-trip-periodic, the cavity is a time-periodic (Floquet) system; computing the topological invariants of the synthetic bands and testing their survival under Kerr nonlinearity is a natural next step the paper points toward but does not perform.","The soliton-drift law vd = −d2k/3 is a clean quantitative prediction that could be tested in isolation: track a soliton's arrival-time shift during a detuning sweep in a fiber cavity and compare with the predicted velocity, separating the drift mechanism from the spectral fits reported.","The paper's outlook extends the GEE to χ(2) effects such as second-harmonic generation and parametric down-conversion; if the periodic-pump structure survives those additions, the framework would cover microcombs well beyond the EO–Kerr systems demonstrated here."],"forward_implications":["Strong-coupling EO and EO–Kerr simulations built on the mean-field LLE misrepresent the dynamics through artificial pump reinjection; GEE-based simulation restores the natural round-trip pump and, the paper argues, resolves the long-standing discrepancy between Kerr theory, simulation, and experiment.","Comb spectra in pure-EO cavities can be computed by solving a discrete-time Heisenberg equation with Ĥig — no rotating-wave approximation needed — and the paper's measured spectra and coupling strengths |Wm| on fiber and lithium-niobate chips match those solutions.","Because the synthetic bands mirror the drive waveform, tailoring the waveform tailors the spectrum: symmetric triangular waves concentrate coupling at a single dominant order m0 ≈ 4Ω0/ωR, asymmetric triangular waves give directional coupling, and slope-dominant dynamics enable single-sideband EO comb generation.","Soliton addressing becomes deterministic when the band slope is steep enough: extra Kerr solitons drift with velocity vd = −d2k/3 and collide with the contracting modulation-instability boundary, relaxing the required dispersion and pump power; stronger EO modulation enlarges this addressing regime.","Strong-coupling band overlap lets EO pulses and Kerr solitons coexist in one cavity, a regime the paper reports observing for the first time, and the GEE provides the model for further hybrid EO–Kerr comb exploration on thin-film lithium niobate."],"fun_headline_variants":["New equation overthrows flawed mean-field model for EO combs","Synthetic band structure mirrors drive waveform, making combs tunable","Strong-coupling comb spectra become programmable via GEE","Old comb model fails; new universal equation steps in"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"Everything depends on two derivations that appear only in the missing supplementary information: the reduction of the round-trip-integrated GEE to the discrete-time Heisenberg equation with Ĥig, and the choice of Ĥig = 0 as the excitation condition that yields the band–wave correspondence; a reader of the main text alone cannot check either step.","fun_headline_variants_meta":{"raw":{"variants":["New equation overthrows flawed mean-field model for EO combs","Synthetic band structure mirrors drive waveform, making combs tunable","Strong-coupling comb spectra become programmable via GEE","Old comb model fails; new universal equation steps in"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000231,"raw_usage":{"total_tokens":1315,"prompt_tokens":727,"completion_tokens":588,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":471,"completion_tokens_details":{"reasoning_tokens":520}},"tokens_in":471,"tokens_out":588,"duration_ms":6682,"temperature":1.0,"reasoning_tokens":520,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T20:04:26.571115+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the printed GEE and Ĥig side by side: numerically integrate the GEE for a triangular drive with Ω0 between one and two free spectral ranges, and compare the steady-state spectrum with the solution of the discrete-time Heisenberg equation from Eqs. (4)–(5); any mismatch refutes the reduction. Independently, measure the band structure in a fiber EO cavity by sweeping detuning Δ at fixed drive and check that transmission resonances lie on E(φ)/ħ = ω0 − πfR V(−φ/ωR)/Vπ − nωR for the appropriate n; the paper asserts this universal relation, so a measured band that is not a scaled mirror im","supporting_citations":[],"review_version":1}