{"id":"a088c9ba-4237-44a9-a067-91c4966978aa","arxiv_id":"1908.02060","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An analytic scattering-matrix calculation yields the spectra and photon-number correlations of spontaneous emission from a moving refractive index front in a dispersive dielectric.","lead":"This paper calculates the quantum light emitted when a sharp moving change in refractive index sweeps through a transparent material, an optical stand-in for a black hole horizon. It gives analytic formulas for the spectrum of emitted photons and for correlations between them, which experiments could look for.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative signatures (shark-fin peaks, near-unity correlations) are computed for an infinitely sharp single step; a Kerr RIF is a finite, smooth, two-boundary pulse, and the paper's robustness study never varies the profile shape.","rationale":"This concern is load-bearing because the abstract and conclusion assert that the calculated spectra and correlations are genuine features of optical analogues, and the reader's strongest claim says they describe what an experiment should observe. If the profile dependence materially changes the scattering matrix, every figure and table changes; varying step height or pulse velocity cannot fix that. The concern is not about internal inconsistency: the step S-matrix is a well-defined exact solution for the stated model, and the paper is transparent about the step being physically difficult to implement. The missing piece is a demonstration that the sharp-step limit is representative of a realistic smooth pulse. A numerical smooth-profile study using the same Lagrangian is feasible and would settle the point. Secondary issues, including that only one kinematic scenario is explicitly derived in the main text with the rest relegated to the first author's thesis, and that the conclusion generalizes to non-optical analogues without a supporting calculation, reinforce the conditional verdict but are less decisive than the profile idealization. The reader's conditional verdict already captures this risk; the check proposed here would either remove the concern or force the paper to present its results explicitly as an idealized step-limit calculation.","tokens_in":23948,"tokens_out":10682,"duration_ms":123220,"concrete_test":"Replace the Heaviside step by a smooth profile n(x) = n_R + delta_n [1 - tanh(x/w)]/2 for a single front, and by a finite pulse n_R + delta_n exp(-x^2/w^2), with delta_n = 2e-6 and u = 2c/3 as in Section V. Solve the same linear Hopfield equations (A1)-(A6) for the asymptotic in/out mode coefficients, for example by transfer-matrix or WKB/shooting, ensuring the result reduces to the analytic S-matrix as w tends to 0. Compute Eq. (3), Eq. (15), and Eq. (9) for w = 0, 10 nm, 100 nm, 1 micrometre, 10 micrometres, and 100 micrometres. If the shark-fin peaks survive and the correlation coefficients remain above about 0.9 for w in the sub-micron to micron range typical of Kerr fronts, the step idealization is benign; if they shift, split, or drop sharply, the quantitative predictions must be labeled step-limit results rather than direct predictions for Kerr-pulse analogues.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing weakness is the external validity of the step idealization on which every quantitative result rests. The paper computes spectra and correlations for a single infinitely sharp boundary (Figs. 4-7, Table I), with matching conditions at x=0 (Appendix B). A Kerr-generated RIF is a finite, smooth, two-edged pulse; the introduction concedes that the step 'is physically hard to implement' and that outputs 'will be further modified by gradual horizons, two-horizon interactions' and nonlinear effects. The claimed robustness (Section V C, Table I) varies only the step height and velocity, never the profile shape or pulse length. In smooth profiles, mode conversion is controlled by the local gradient and adiabaticity rather than by point matching, and a finite pulse creates two stationary horizons whose scattering can interfere; the step result is not automatically the limit of the smooth result. Thus the specific predictions, including the shark-fin spectral shape and the correlation peaks at (397 nm, 210 nm) and (3.6 micrometres, 227 nm), are not yet shown to describe what an actual optical analogue experiment should observe.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytical scattering-matrix framework for spontaneous emission at a moving refractive-index step in a dispersive dielectric described by the Hopfield model. It constructs global modes for all kinematic scenarios, including evanescent waves, derives the S-matrix from the field matching conditions at the step, and uses it to compute photon flux spectra and photon-number correlations in both the comoving and laboratory frames. Numerical results for fused silica predict 'shark-fin' spectral peaks at white-hole and black-hole horizon frequencies and near-unity correlations at specific wavelength pairs, with correlation structures reported to be robust to changes in step height and RIF velocity.","tokens_in":24205,"tokens_out":19899,"duration_ms":185977,"significance":"If correct, the paper provides a complete and explicit analytical route from a microscopic Hopfield model to the two-photon observables of an optical analogue of gravity, including a step-by-step derivation of the covariance and variance formulas and a quasi-unitarity check of the S-matrix. The expressions relating flux and correlations directly to S-matrix elements are general and could be reused for other stationary mode-conversion problems. The main limitation is that the quantitative predictions are computed for an idealized sharp step, which the authors themselves note is hard to realize experimentally; the paper therefore demonstrates a calculational method and identifies qualitative types of emission, rather than delivering directly testable experimental spectra.","major_comments":[{"comment":"The quantitative predictions, including the shark-fin spectral shape and the peak correlations of 0.97 and 0.92 at (397 nm, 210 nm) and (3.6 micrometres, 227 nm), are computed for a single infinitely sharp step with matching at x=0 (Appendix B). The robustness study in Section V C and Table I varies only the step height and the RIF velocity, never the profile shape or the pulse length. The Introduction explicitly states that the step 'is physically hard to implement' and that outputs 'will be further modified by gradual horizons, two-horizon interactions and other nonlinear effects'. Since a Kerr-generated RIF has a smooth finite-width profile, the step result is not automatically the limit of the smooth result, and the paper does not demonstrate that the computed spectra and correlations describe a real experiment. The abstract and Section V B should be tempered to state that these are predictions for the idealized step, or the paper should add a concrete discussion (even a qualitative adiabaticity estimate) of how smooth profiles are expected to modify the signatures.","section":"Section V, Figs. 4-7, Table I"},{"comment":"The central claim that 'We calculate the scattering matrix at all frequencies (from all kinematic scenarios)' is only explicitly demonstrated for mode scenario c (8 propagating modes on either side). Appendix C works out the sigma matrices and the S-matrix for that case and refers to the first author's PhD thesis ([36]) for 'further calculations'. Scenarios with evanescent (complex) modes are described verbally in Appendix A 2, but no explicit construction, matching solution, or quasi-unitarity verification for such a scenario appears in the manuscript. To substantiate the all-frequency claim and make the paper self-contained, the authors should include at least one explicit example with evanescent modes (e.g., scenario d) in an appendix, or provide the general algorithm in sufficient detail to reproduce those cases.","section":"Section IV and Appendix C"},{"comment":"The treatment of unphysical (evanescent) modes assigns them a norm of unity in the metric g (Appendix D) and includes them in the S-matrix construction. The physical observables (3) and (9) sum only over beta not in {alpha}, which excludes the unphysical modes if they are placed in {alpha}, so physical predictions are expected to be independent of this arbitrary normalization, but this independence is not shown. The authors should demonstrate that the physical block of the S-matrix is unaffected by the arbitrary norm assignment, or justify why the quasi-unitarity relation with such an assignment is a meaningful numerical check.","section":"Appendices A 2 and D"}],"minor_comments":[{"comment":"The prefactor 'Delta^2/Delta1 Delta2' in the last line of Eq. (9) is printed in the text as 'Delta2/Delta1 Delta2', which is ambiguous and would reduce to 1/Delta1 if taken literally; this should be corrected to clarify that the numerator contains Delta^2.","section":"Equation (9)"},{"comment":"Near the discussion of Fig. 7, the text reads 'starting at 237 mn' and should read '237 nm'.","section":"Section V B"},{"comment":"The distinction between the scattering matrix S and the norm-adjusted matrix S introduced in Appendix E is not defined in the main text at Eq. (6); a short note clarifying that rows of negative-norm modes are complex-conjugated would improve readability.","section":"Section III and Appendix E"},{"comment":"The abstract states 'We supplement the field theory in curved spacetime for this model', but the Lagrangian (A1) is written in a moving frame and the analogue metric is not explicitly used; a brief clarification of how the analogue spacetime enters the calculation would help readers.","section":"Abstract and Introduction"},{"comment":"The caption of Fig. 5 uses labels uL, lL, nlL, and nuL for non-optical modes without defining them; these labels should be defined or referenced to Fig. 8.","section":"Figure 5"}],"recommendation":"major_revision","confidential_remarks":"The evanescent-mode scattering calculation is partly delegated to the first author's PhD thesis; if the thesis is not readily available, this could be a novelty-disclosure concern. The external-validity issue is a scope problem rather than a correctness problem: the internal derivations appear coherent, but the quantitative observational claims should be framed more cautiously."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is the natural follow-up to Finazzi and Carusotto's step-RIF scattering work, and it is not a rehash: it extends the scattering problem to all kinematic scenarios including evanescent modes, and it derives spectrally resolved photon-number correlations for an optical analogue for the first time. The correlation formulas in Eqs. (7)-(9) are general results for any stationary scattering process, not just this system, and they are derived carefully in Appendix E. The S-matrix construction in Section IV is coherent, the quasi-unitarity check in Appendix D is included, and the free parameters (Sellmeier coefficients, step height, velocity) are physical inputs rather than fitted quantities.\n\nThe paper is honest about its own limits. The introduction concedes that the step is physically hard to implement and that smooth profiles, two-horizon interactions, and nonlinearities will modify the outputs. That candor makes the main caveat less damaging than it might be: the quantitative predictions, including the shark-fin spectra and the near-unity correlation peaks at (397 nm, 210 nm) and (3.6 µm, 227 nm), are computed for a single infinitely sharp boundary, and the robustness studies in Fig. 4 and Table I vary only step height and velocity, never the profile shape or pulse length. A Kerr RIF is finite and smooth, with two edges; mode conversion in that setting is controlled by the local gradient and by interference between the two horizons. So the numbers in Figs. 6 and 7 are not yet shown to be what an actual experiment will observe. The authors mostly say as much, but the abstract's claim that these are genuine features of \"optical and non-optical analogues\" goes beyond what is computed, since only fused silica with a three-resonance Sellmeier response appears.\n\nThe other soft spot is self-containedness. The claim that the S-matrix is derived for all frequencies and kinematic scenarios is not fully backed inside the paper: only scenario (c) is worked out in Appendix C, and the rest is referred to the first author's thesis. A referee should ask for the other scenarios to be included or for an explicit statement that they follow from the same construction. The citation pattern is otherwise fine — Finazzi-Carusotto, the Hopfield model, and the Barnett-Loudon quantization line are all credited; the thesis dependence is the only structural issue.\n\nNone of this is fatal. The derivation holds together, the correlation formulas are genuinely new and useful, and the kinematic classification in Section II is worth having. This paper is for people working on optical analogue gravity — theorists who want the toolkit and experimentalists who want target wavelengths to search for — and it deserves a serious referee. My recommendation: send it to review, and accept with revisions that make the other scenarios self-contained, soften or support the non-optical generalization, and state plainly that the sharp-step predictions are the idealized limit, not the experimental forecast.","headline":"Honest, useful extension of the step-RIF programme to two-photon correlations and lab-frame observables; the sharp-step idealization and the thesis dependency are the real caveats.","tokens_in":24673,"tokens_out":5670,"would_cite":true,"duration_ms":54109,"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":"All spontaneous emission at a moving refractive-index front in a dispersive dielectric follows from a single analytically computed scattering matrix, which yields concrete spectra and correlation maps for a fused-silica analogue of black…","keywords":["analogue gravity","Hawking radiation","Hopfield model","scattering matrix","photon-number correlations","moving refractive index front","dispersive media","spontaneous emission"],"falsifier":"Measure the photon-number correlation coefficient at the predicted black-hole pair (397 nm, 210 nm) using a 400 nm pump pulse in fused silica: if C is far below 0.97, or if the spectral peak at 210 nm lacks the sharp 'shark fin' cut-off at the horizon boundary, the step-like Hopfield model does not capture the emission.","tokens_in":23749,"feed_emoji":"⚛️","tokens_out":7562,"duration_ms":72519,"temperature":0.7,"pith_summary":"This paper aims to show that a step-like moving refractive index front in a dispersive optical medium—a laboratory analogue of an event horizon—can be solved exactly: from the Hopfield model of light coupled to three polarisation resonances, the authors build the complete scattering matrix for all modes and frequencies, including evanescent waves. From that matrix they derive, in closed form, the photon flux and the spectrally resolved photon-number correlations in both the frame co-moving with the front and the laboratory frame. Applying the formulas to bulk fused silica, they predict distinct emission types depending on whether a black-hole or white-hole horizon exists at a given frequency, and show the correlation structure is robust when the index step's height is varied over five orders of magnitude. If correct, the paper gives experimenters exact target wavelengths and correlation strengths—for example a 0.97 correlation at (397 nm, 210 nm)—for observing analogue Hawking emission in optical fibres.","feed_headline":"Gravity-analogue photon pairs predicted to correlate at 0.97","feed_subtitle":"A moving index step in fused silica yields exact Hawking-pair spectra and near-unity correlations at fixed wavelengths.","key_machinery":"The machinery is the Hopfield model—a Lagrangian in which light (a scalar potential) is coupled to three polarisation oscillator fields with Sellmeier dispersion—combined with matching conditions at a step-like refractive index front. The matching conditions, derived in Appendix B, require continuity of the potential, the polarisation fields, and their derivatives across the step, with a discontinuity in the polarisation-derivative set by the elastic constants. The scattering matrix S is obtained as the transformation between global in and out mode bases, S = (σ_out_L)^{-1} σ_in_L, and is shown to be quasi-unitary with respect to the Klein-Gordon norm, which preserves the canonical commutation relations. All observables—flux (3), covariance (7), correlation coefficient (9), and g(2) (11)—are expressed directly in terms of S, so the entire emission phenomenology rests on this matrix.","core_discovery":"The central discovery is that the quantum emission at a moving refractive index front is completely characterised by an analytically constructed scattering matrix. The authors supplement the Hopfield-model field theory in curved spacetime, treat the index change as a sharp step, and match the electromagnetic potential, the three polarisation fields, and their derivatives at the boundary. Solving the resulting linear system yields the scattering matrix in all kinematic scenarios—horizonless, white-hole, black-hole, and mixed—including the complex, evanescent modes. From this matrix the authors derive the output photon flux in Eq. (3) and the normally ordered photon-number covariance in Eq. (7), which is independent of the photon flux, and hence the full correlation coefficient (9) and g(2) functions. Applied to fused silica, the formalism predicts narrow 'shark-fin' spectral peaks over the horizon intervals and near-unity pair correlations—0.97 at the black-hole pair (397 nm, 210 nm) and 0.92 at the white-hole pair (3.6 µm, 227 nm)—with the overall structure stable against changes in step height and pulse velocity.","pith_inferences":["Because the correlation coefficient in Eq. (9) is independent of the photon flux, the correlation map may remain measurable even where the flux is too low for direct detection; this makes the near-unity contours, rather than peak heights, the most decisive experimental signature.","The same S-matrix formalism could be transferred to other linear bosonic platforms with step-like inhomogeneities—polariton condensates, cold atoms, or superconducting circuits—where analogous horizon pairs would display the same shark-fin spectral cut-offs if dispersion is tuned to match the Hopfield-Sellmeier structure.","A smooth-front calculation is the natural next test: if a realistic Kerr pulse only broadens the shark-fin peaks rather than shifting their centres, the analytic step model would still serve as the quantitative backbone for designing analogue-gravity experiments."],"forward_implications":["A fused-silica experiment with a 400 nm pump (RIF velocity u = 2c/3) should observe the black-hole partner at 210 nm and the signal at 398 nm, with a photon-number correlation coefficient of 0.97 at (397 nm, 210 nm).","The white-hole pair at 3.6 µm and 227 nm, with C = 0.92, provides a spectrally well-separated target that avoids pump contamination.","The 'shark fin' spectral line shape, cut off abruptly at the horizon interval boundary, is a robust signature of horizon emission, persisting over five orders of magnitude of index-step height.","The method extends to finite-length pulse profiles, allowing theoretical predictions for the two-horizon configurations that real Kerr pulses implement."],"supporting_citations":[{"why":"Supplies the original step-like refractive-index-front scattering approach and the Hopfield-model framework that this paper generalises to all kinematic scenarios, including evanescent waves.","marker":"[42]"},{"why":"Provides the microscopic Lagrangian of light coupled to polarisation oscillators used to model the dielectric across the front.","marker":"[48]"},{"why":"Adapts the Hopfield model to analogue gravity in dispersive media, the canonical field-theory setup the authors supplement with curved-spacetime methods.","marker":"[47]"},{"why":"Earlier work by the same authors that defined the four kinematic scenarios and computed the output spectral density, which this paper extends to photon-number correlations.","marker":"[38]"},{"why":"Source of the Sellmeier coefficients for bulk fused silica used in the computed spectra and correlation maps.","marker":"[44]"}],"fun_headline_variants":["Exact Hawking spectrum from moving index step","Quantum emission analytics for optical gravity analogues","Gravity-analogue pairs correlate at 0.97","Analytic scattering matrix for optical horizons","Fused silica Hawking pairs predicted precisely"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions rest on modelling the moving refractive-index perturbation as a sharp, lossless step in a three-resonance Hopfield dielectric, ignoring phonon interactions, four-wave mixing, and other nonlinearities; a real Kerr pulse has a smooth index profile, and the paper acknowledges that the step 'is physically hard to implement'.","fun_headline_variants_meta":{"raw":{"variants":["Exact Hawking spectrum from moving index step","Quantum emission analytics for optical gravity analogues","Gravity-analogue pairs correlate at 0.97","Analytic scattering matrix for optical horizons","Fused silica Hawking pairs predicted precisely"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000215,"raw_usage":{"total_tokens":1417,"prompt_tokens":919,"completion_tokens":498,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":430}},"tokens_in":535,"tokens_out":498,"duration_ms":5546,"temperature":1.0,"reasoning_tokens":430,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:55:26.802879+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the photon-number correlation coefficient at the predicted black-hole pair (397 nm, 210 nm) using a 400 nm pump pulse in fused silica: if C is far below 0.97, or if the spectral peak at 210 nm lacks the sharp 'shark fin' cut-off at the horizon boundary, the step-like Hopfield model does not capture the emission.","supporting_citations":[{"cited_title":"Finazzi and I","cited_arxiv_id":null,"evidence_quote":"Supplies the original step-like refractive-index-front scattering approach and the Hopfield-model framework that this paper generalises to all kinematic scenarios, including evanescent waves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the microscopic Lagrangian of light coupled to polarisation oscillators used to model the dielectric across the front."},{"cited_title":"Sch ¨utzhold, G","cited_arxiv_id":null,"evidence_quote":"Adapts the Hopfield model to analogue gravity in dispersive media, the canonical field-theory setup the authors supplement with curved-spacetime methods."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier work by the same authors that defined the four kinematic scenarios and computed the output spectral density, which this paper extends to photon-number correlations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the Sellmeier coefficients for bulk fused silica used in the computed spectra and correlation maps."}],"review_version":1}