{"id":"80754661-01f7-470b-9540-f16692107473","arxiv_id":"2504.17657","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new mean-field model, a Fabry-Perot Lugiato-Lefever equation with a nonlocal Brillouin response, reproduces coupled-wave simulations of Kerr-Brillouin combs while speeding up computation by up to four orders of magnitude.","lead":"Physicists derive a single averaged equation that captures both Kerr and Brillouin light interactions in a Fabry-Perot fiber cavity. The equation matches full simulations while running up to about four orders of magnitude faster, and it yields a simple formula for when the comb forms.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model derivation and numerical comparisons are internally consistent; the load-bearing gap is that no simulation code or data are released, so the headline speedup and accuracy figures cannot be independently reproduced.","rationale":"The reader's weakest-assumption analysis focused on the slow-variation condition, but a quantitative check shows that condition is well satisfied for the simulated parameters; it is therefore not the most load-bearing risk. The paper's derivation of Eq. (48) from the CWE is internally consistent: the modal expansion, rotating-wave averaging, and Brillouin response all check out algebraically, and the numerical agreement in Figs. 3-5 is strong evidence for the model's accuracy. The remaining concern is reproducibility of the speedup claims, which is exactly the basis for the reader's CONDITIONAL verdict. Since the reader already flagged the absence of released code and data and set the verdict accordingly, no change to the verdict is needed. The stress-test pass therefore reports no scientific objection beyond the reproducibility gap.","tokens_in":16149,"tokens_out":27426,"duration_ms":263745,"concrete_test":"Release the exact simulation code and scripted parameters used for Figs. 3-5, including both the CWE solver and the mean-field split-step solver, then independently rerun all three cases on a neutral platform to confirm the reported runtime ratios (2000, 8000, 9600) and the spectral agreement down to -80 dB (Fig. 5d) and -150 dB (Fig. 3d).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has two parts: (i) Eq. (48) accurately reproduces the coupled-wave-equation (CWE) dynamics, and (ii) the mean-field solver runs 2000-9600 times faster. Part (i) is supported by careful comparisons for the three operating points in Figs. 3-5, including a 20 THz comb, and the slow-variation assumption invoked in the derivation is comfortably satisfied for the stated parameters: the fastest modal dynamics are of order 10^7 s^{-1}, while even the first cavity mode has ω_n ≈ 7.4×10^9 rad/s, so |a_dot/a| ≪ ω_n. Part (ii) is the more load-bearing element. The speedup figures rest entirely on the authors' implementations, for which neither code nor input data are released. Without access to the CWE solver (finite-difference predictor-corrector with the characteristic condition) and the mean-field split-step solver, a reader cannot verify that the comparison is fair, that the time steps and spatial grids are equivalently converged, or that the runtime measurements are reproducible. This is not a claim of misconduct; it is an evidentiary gap for the quantitative performance claim, which is a headline result of the paper.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a mean-field model for Fabry-Perot resonators with both Kerr and Brillouin nonlinearities. Starting from the coupled wave equations (CWE) for the forward and backward optical fields and the acoustic wave, the authors perform a modal expansion, project onto cavity modes, and under the good-cavity/slow-amplitude assumptions obtain a single integro-differential equation, Eq. (48), with a nonlocal Brillouin response. They derive steady-state curves, a Brillouin gain spectrum, and a linear stability analysis yielding a perturbation growth rate. The model is then validated numerically against the original CWE for three operating points, including a 20 THz comb, with reported speedups between 2000 and 9600 times.","tokens_in":16334,"tokens_out":15739,"duration_ms":155685,"significance":"If the result holds, this is a useful and potentially influential contribution: it provides a compact Lugiato-Lefever-type description of Kerr-Brillouin comb generation that is amenable to standard split-step Fourier solvers and to analytical stability analysis. The derivation from the prior coupled-wave framework of Dong and Winful is a genuine strength, as is the numerical comparison over a wide dynamic range (agreement to about 150 dB below the pump line) and the reduction of the growth-rate formula to the pure Brillouin gain in the appropriate limit. The claimed computational speedup is very attractive, but it currently rests on unpublished implementation details rather than on a released code or reproducible benchmark. The central physical derivation is plausible, but the printed equations contain an internal inconsistency in the dispersion term that must be resolved before the manuscript can be considered self-consistent.","major_comments":[{"comment":"The GVD term is not self-consistent across the printed equations. From Eq. (43), the linear dispersive term is -i beta2/(2 beta1) omega_n^2 a_n. Multiplying by Tr = 2 beta1 L and summing against e^{-i beta1 omega_n z}, and using d^2 psi/dz^2 = -beta1^2 sum omega_n^2 a_n e^{-i beta1 omega_n z}, gives a contribution +i L beta2/beta1^2 d^2 psi/dz^2 to Eq. (48). The printed Eq. (48) instead has -i L beta2/beta1^2 d^2 psi/dz^2; the sign would be consistent only if Eq. (43) contained +i beta2/(2 beta1) omega_n^2 a_n. There is also a factor-of-two discrepancy: Eq. (55), with d2 = beta2/beta1^3, yields a modal contribution D = (d2/2) Tr k_n^2 = L beta2 omega_n^2, whereas Eq. (59) contains 2 L beta2 omega_n^2. All of Eqs. (43), (48), (55), and (59) cannot be simultaneously correct; at least one is misprinted. Since the numerical comparisons in Figs. 3-5 are very close, it appears that the implemented equation differs from at least one printed version, and the manuscript should be corrected so that the derivation is internally consistent and reproducible from the displayed equations alone.","section":"III.B, Eq. (48) vs. Eq. (43); IV, Eqs. (55) and (59)"},{"comment":"The headline quantitative claim, namely speedups of about 2000, 8000, and 9600 relative to the CWE solver, rests entirely on the authors' implementations. No simulation code, input parameter files, or benchmark scripts are provided, and the only implementation details are the brief description of the split-step/FFT method for Eq. (48) and the predictor-corrector characteristic scheme for the CWE. A reader cannot verify that the two solvers are equally converged, what time steps and spatial grids were used for each case, how the runtime comparison was made, or what hardware and compiler were used. This is an evidentiary gap for a central performance claim. I ask the authors to provide a reproducibility statement with either released code or, at minimum, a detailed benchmark report including convergence tests in dt and dz for both solvers and the criterion used to declare the spectra matched.","section":"Section V, Figs. 3-5 and the speedup claims"}],"minor_comments":[{"comment":"There is a typo in 'unfiorm- or mean-field limit' in the introduction; it should read 'uniform-field or mean-field limit'.","section":"Introduction"},{"comment":"The text and caption report '8092 spatial modes' for Fig. 4; this is presumably a typo for 8192 (2^13).","section":"Section V, Fig. 4"},{"comment":"The sentence 'Equation (60) represents a double Lorentzian, with maximum gain (absorption) at negative (postitive) frequency shift ∓Omega_B' is confusingly worded; it should state explicitly that the gain maximum is at omega = -Omega_B and the absorption maximum at omega = +Omega_B.","section":"Eq. (60)"},{"comment":"The notation in the Brillouin term of Eq. (48) is ambiguous because the symbol * is used both for complex conjugation and for convolution; the two should be clearly distinguished, and the argument of the convolution should be written explicitly.","section":"Eq. (48) and Eq. (49)"},{"comment":"The phrase 'as seldom done in nonlinear fiber optics devices' in the conclusion is awkward; consider rewording to 'which is rarely done' or similar.","section":"Section VI"}],"recommendation":"major_revision","confidential_remarks":"I see no reason to question the authors' integrity; the numerical agreement is strong enough that the equation inconsistencies are likely typographical errors rather than fundamental flaws. However, the printed derivation must be made internally consistent, and the speedup claim needs to be backed by at least a detailed benchmark description or released code. If those issues are addressed, I would expect the paper to be publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — read it. The main thing to know: Eq. (48) is a genuine generalization of the Fabry-Perot Lugiato-Lefever equation. It folds SBS into a nonlocal term and handles Brillouin cascades to all orders from a single set of modal amplitudes, with no auxiliary acoustic variables. The derivation from Dong-Winful is a standard mean-field projection, carefully laid out, and it passes the right checks: zero Brillouin gain reduces to the known FPLLE, and the undepleted-pump limit recovers the double-Lorentzian Brillouin gain. That is what a model paper needs.\n\nThe numerical section does the heavy lifting. They compare mean-field and CWE at three operating points, the largest running 2^15 modes and spanning a 20 THz comb. Spectra overlap even at components 150 dB below the main line; temporal profiles nearly superimpose. That is convincing evidence that Eq. (48) actually replaces the coupled equations in the regimes it claims.\n\nThe slow-variation worry flagged in the stress-test does not bite for these parameters: fastest modal dynamics are around 10^7 per second against mode frequencies around 10^10 rad/s, so the separation is huge. The lack of a formal validity bound is a minor gap.\n\nThe real soft spot is reproducibility. No code or input data are released. The speed-up claims—2000x, 8000x, 9600x—rest entirely on the authors' implementations. I have no reason to doubt the order of magnitude; CWE with an explicit CFL constraint is much more expensive than split-step. But “no reason to doubt” is not the same as checking convergence and fairness. This is an evidentiary gap in a headline claim, easily fixed by releasing both solvers and runtimes. It is not a scientific flaw.\n\nCitations are fine. Refs. [49,50,51] get proper credit, and [27] is an extension source, not a substitute for evidence.\n\nThis paper is for fiber-photonics and microcomb modelers, particularly anyone studying the Kerr-Brillouin interplay in Fabry-Perot cavities. I would cite it, and I would bring it to a reading group. Send it to peer review; ask the authors for code and data as a condition of acceptance.","headline":"A clean extension of the Fabry-Perot LLE with a nonlocal Brillouin response that matches coupled-wave simulations across extreme comb regimes; the main gap is missing code, not the physics.","tokens_in":16927,"tokens_out":2815,"would_cite":true,"duration_ms":27404,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.65.Es","42.65.Ky","42.60.Da"],"model":"deepseek-v4-flash","headline":"Starting from three coupled wave equations, this paper derives one mean-field envelope equation that reproduces Kerr-Brillouin frequency combs down to 150 dB below the carrier and runs up to four orders of magnitude faster.","keywords":["Kerr-Brillouin frequency combs","Fabry-Perot resonators","stimulated Brillouin scattering","Kerr nonlinearity","mean-field model","nonlocal Brillouin response","coupled wave equations","split-step Fourier method"],"falsifier":"A concrete check: simulate the same high-finesse cavity of Table I but reduce the mirror reflectivities so that the finesse drops to about 20, keep the same pump power and detuning, and compare the spectra produced by Eq. (48) and by the full coupled wave equations. If the comb lines beyond the first Brillouin sideband no longer match, the slow-amplitude / good-cavity premise is the controlling assumption; if they still match, the equation is more robust than the derivation suggests.","tokens_in":15923,"feed_emoji":"💡","tokens_out":11357,"duration_ms":106394,"temperature":0.7,"pith_summary":"The paper derives a single mean-field partial differential equation, Eq. (48), that describes a Fabry-Perot resonator whose medium responds with both Kerr and Brillouin nonlinearities. It starts from the three coupled wave equations for the forward and backward optical fields and the acoustic wave, and eliminates the acoustic variable through a nonlocal response function. The authors show that the resulting equation reproduces broadband Kerr-Brillouin frequency combs generated by cascaded stimulated Brillouin scattering and four-wave mixing, with spectral agreement down to 150 dB below the main line and computation times up to four orders of magnitude shorter than the coupled-wave benchmark. The same equation yields an explicit growth-rate formula for harmonic perturbations, which explains how the comb's line spacing is selected by the combined Kerr-Brillouin gain rather than simply by the cavity free-spectral range. If correct, this gives researchers a fast and interpretable tool for designing combs in which Brillouin scattering is used as a feature rather than avoided.","feed_headline":"A single equation models Kerr-Brillouin combs up to 10,000x faster","feed_subtitle":"It matches coupled-wave spectra down to 150 dB below the carrier and explains the comb's line spacing.","key_machinery":"The load-bearing object is Eq. (48), a single mean-field envelope equation for the backward field in a Fabry-Perot resonator, written as a driven, damped nonlinear Schr\\\"odinger-type equation with a nonlocal Brillouin convolution. The derivation rests on a modal expansion: forward and backward fields share one set of cavity-mode amplitudes $a_n(t)$, and the acoustic wave is eliminated through the frequency-domain response $H_B(\\omega)=\\Omega_B\\Gamma_B/(\\Omega_B^2-\\omega^2-i\\omega\\Gamma_B)$, whose inverse Fourier transform, periodically replicated, forms the spatial kernel $h_B(z)$ of the convolution in Eq. (48). This kernel does the work of the acoustic equation: it converts the Brillouin cascade into a single nonlocal term that can be evaluated by FFT within a split-step Fourier integrator, and it makes the stability calculation a $2\\times2$ matrix problem.","core_discovery":"The central claim is that the full coupled-wave description of stimulated Brillouin scattering plus the Kerr effect in a high-finesse Fabry-Perot cavity—two counter-propagating optical envelopes together with a damped acoustic oscillator—reduces, under the mean-field limit, to a single evolution equation for one envelope, Eq. (48). The Brillouin contribution appears as a nonlocal convolution with a kernel $h_B$ built from the response $H_B(\\omega)=\\Omega_B\\Gamma_B/(\\Omega_B^2-\\omega^2-i\\omega\\Gamma_B)$, so every Stokes and anti-Stokes order of the SBS cascade is generated without tracking the acoustic wave as a separate variable. The authors demonstrate numerically that this one equation reproduces the coupled-wave spectra in fiber Fabry-Perot resonators—down to 150 dB below the primary line, including combs spanning more than 20 THz—and that it runs 2,000 to 9,600 times faster than the coupled-wave solver in their examples. They also derive from Eq. (48) a compact linear-stability growth rate for harmonic perturbations, Eq. (59), that identifies which cavity modes are excited by the combined Kerr-Brillouin gain.","pith_inferences":["The same elimination-by-convolution should extend to any delayed nonlinearity in a Fabry-Perot cavity—Raman scattering, thermal nonlinearities—by substituting the appropriate response kernel; the paper itself treats only the Brillouin case.","The growth-rate formula suggests a practical tuning knob: by choosing detuning and pump power, the Kerr-broadened Brillouin lobe can be moved across cavity modes, making the comb repetition rate selectable in integer steps of the free-spectral range.","The slow-amplitude assumption is the boundary of validity; a quantitative criterion (for instance, finesse times comb bandwidth relative to the Brillouin frequency) would be the natural next test, and can be checked by low-finesse comparisons against the coupled-wave equations."],"forward_implications":["The reported examples reduce runtime from 140 minutes to 3.7 seconds (about 2000x), from 1327 minutes to 10.4 seconds (about 8000x), and from 6716 minutes to 42 seconds (about 9600x).","The linear-stability formula gives the gain of each cavity mode directly, so the dominant sideband spacing (nine times the free-spectral range in the examples) can be predicted from parameters, not only from full simulation.","Because the derivation includes group-velocity dispersion of either sign, the model covers both normal-dispersion switching-wave combs and the anomalous-dispersion soliton combs that are standard in Kerr resonators.","The equation resolves SBS cascades to all orders in a single convolution, so it can describe combs where the pure Brillouin gain peak would fall between cavity resonances yet comb generation still occurs through the joint Kerr-Brillouin gain."],"supporting_citations":[{"why":"Supplies the Fabry-Perot Kerr-comb mean-field equation that Eq. (48) generalizes with the Brillouin response.","marker":"[50]"},{"why":"Provides the three coupled wave equations for forward and backward optical fields and the acoustic wave that are the derivation's starting point.","marker":"[51]"},{"why":"Gives the mean-field technique: transformed variables that make boundary conditions periodic, followed by the good-cavity expansion.","marker":"[43]"},{"why":"The authors' earlier Kerr-only Fabry-Perot mean-field and modulation-instability treatment, whose modal expansion is extended here to include SBS to all orders.","marker":"[49]"},{"why":"Reports the experimental Brillouin-induced Kerr comb in a normal-dispersion fiber Fabry-Perot cavity whose parameters the simulations adopt.","marker":"[27]"},{"why":"Supplies the stable numerical solver for the coupled wave equations used as the accuracy benchmark.","marker":"[53]"},{"why":"Identifies the switching-wave mechanism in normal-dispersion fiber Fabry-Perot resonators reproduced in the 20 THz comb simulation.","marker":"[6]"}],"fun_headline_variants":["One equation models Kerr-Brillouin combs 10,000x faster","Single equation for Kerr-Brillouin combs in Fabry-Perot","Kerr-Brillouin comb model: one equation, 10k speedup","Fast single-equation model for Kerr-Brillouin combs","One equation unifies Kerr-Brillouin comb simulation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument rests on the premise that the cavity-mode amplitudes change slowly over one round trip (high mirror reflectivity, weak nonlinearity, weak dispersion), so that second time derivatives and rapidly oscillating terms can be dropped from the modal equations.","fun_headline_variants_meta":{"raw":{"variants":["One equation models Kerr-Brillouin combs 10,000x faster","Single equation for Kerr-Brillouin combs in Fabry-Perot","Kerr-Brillouin comb model: one equation, 10k speedup","Fast single-equation model for Kerr-Brillouin combs","One equation unifies Kerr-Brillouin comb simulation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000552,"raw_usage":{"total_tokens":2646,"prompt_tokens":970,"completion_tokens":1676,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":586,"completion_tokens_details":{"reasoning_tokens":1583}},"tokens_in":586,"tokens_out":1676,"duration_ms":10151,"temperature":1.0,"reasoning_tokens":1583,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:34:07.154716+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check: simulate the same high-finesse cavity of Table I but reduce the mirror reflectivities so that the finesse drops to about 20, keep the same pump power and detuning, and compare the spectra produced by Eq. (48) and by the full coupled wave equations. If the comb lines beyond the first Brillouin sideband no longer match, the slow-amplitude / good-cavity premise is the controlling assumption; if they still match, the equation is more robust than the derivation suggests.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Fabry-Perot Kerr-comb mean-field equation that Eq. (48) generalizes with the Brillouin response."},{"cited_title":"Kartashov, O","cited_arxiv_id":null,"evidence_quote":"Gives the mean-field technique: transformed variables that make boundary conditions periodic, followed by the good-cavity expansion."},{"cited_title":"Firth, Stability of nonlinear Fabry-Perot resonators, Optics Communications 39, 343 (1981)","cited_arxiv_id":null,"evidence_quote":"The authors' earlier Kerr-only Fabry-Perot mean-field and modulation-instability treatment, whose modal expansion is extended here to include SBS to all orders."},{"cited_title":"Dynamic Interplay Between Kerr Combs and Brillouin Lasing in Fiber Cavities","cited_arxiv_id":"2212.08534","evidence_quote":"Reports the experimental Brillouin-induced Kerr comb in a normal-dispersion fiber Fabry-Perot cavity whose parameters the simulations adopt."},{"cited_title":"Dong and H","cited_arxiv_id":null,"evidence_quote":"Supplies the stable numerical solver for the coupled wave equations used as the accuracy benchmark."}],"review_version":1}