{"id":"410fbc01-58fc-4b89-ae5d-08e263c09fff","arxiv_id":"2412.00958","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A frequency- and time-resolved simulation of entanglement-based QKD that includes multi-pair events and detector imperfections reproduces measured key rates and quantum bit error rates on two SPDC-based systems.","lead":"This paper simulates quantum key distribution systems that use highly entangled photon pairs, including the unwanted multi-pair events that real sources produce. The simulation matches measured key rates and error rates on two experimental links, which could help engineers design and optimize entanglement-based QKD networks.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reconstructed JSA phase, not the spectral amplitude, carries the load; eq. (17)'s decoupling ansatz is untested against phase-sensitive data, so the central agreement claim rests on an unvalidated phase model.","rationale":"Good-faith read: the paper builds a careful covariance-formalism simulation and compares it to multiple experimental scans; I found no obvious algebraic error in the main formulas, and the observed correlations in figs. 6, 7 and 10 are real evidence. The weakest point is not the covariance machinery but the reconstruction of the complex JSA phase from intensity data. This is exactly the reader's weakest_assumption, and I agree with it. The phase enters every interference and dispersion prediction, and the decoupling form (eq. 17) plus quadratic expansion (eq. 19) is a fitting ansatz, not a derived consequence. The paper's own text notes the asymmetry of the type-II spectrum is 'most likely' due to waveguide imperfections; it does not rule out spectrograph response, and no independent phase measurement is reported. This is an addressable validation gap rather than a demonstrated error, so the appropriate outcome is the reader's CONDITIONAL status, unchanged. The concrete check above would settle whether the concern actually moves the predicted observables.","tokens_in":20094,"tokens_out":7041,"duration_ms":74007,"concrete_test":"Perform a phase-sensitive experiment on the same type-II source: measure the phase-basis (Franson) coincidence visibility as a function of narrowband spectral filtering or of added fiber dispersion, and compare it quantitatively with the simulation. To make the check decisive, also run the simulation with an alternative JSA phase that reproduces the same measured power spectrum but violates the decoupling ansatz (for example, Delta k(omega_-, z) = Delta k(omega_-) + omega_- delta k(z), or a minimum-phase reconstruction). If the predicted visibility changes by more than the measurement uncertainty, the power-spectrum fit is underdetermined and eq. (17) is load-bearing; if the predicted visibility is unchanged, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the simulations reproduce the measured key rates and QBERs—requires the complex joint spectral amplitude to be correct, not merely its modulus. In Section IV A 1, the phase of the phase-matching function is recovered by fitting eq. (20) to a measured power spectrum under the decoupling ansatz Delta k(omega_-, z) approx Delta k(omega_-) + delta k(z) (eq. 17) and the quadratic perturbation delta k(xi) = delta k' xi + (1/2) delta k'' xi^2 (eq. 19). A power spectrum fixes only |Phi(omega_-)|. Any frequency-position coupling not of this separable form, or any uncalibrated wavelength-dependent spectrograph response, would change the fitted delta k' and delta k'' and therefore the phase of the propagated JSA. Because the Franson-interference terms in eqs. (A8)-(A9) and the dispersion-induced time-basis QBER in Section IV D (fig. 9) are phase-sensitive, the agreement in figs. 6, 7 and 10 cannot by itself certify the phase: fig. 7 mainly shows the expected cosine fringe shape, and the sweeps in fig. 10 combine many simultaneous effects. Thus the strongest validation claim rests on an untested phase model, with no internal error bar for this assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports time- and frequency-resolved simulations of two experimental entanglement-based QKD systems implementing the time-bin BBM92 protocol, using the covariance-expansion formalism introduced in part I. The simulation includes the joint spectral amplitude reconstructed from measured pump and pair spectra, wavelength-division multiplexing for the type-0 source, chromatic dispersion, interferometer mode mismatch, and SPADs with dark counts, dead times, and afterpulses. The main results are comparisons of simulated and measured sifted key rates and time-basis QBERs for both systems over varying mean pair number, fiber lengths, repetition rates, and detector settings, with the claim that the results agree with measurements.","tokens_in":20405,"tokens_out":5709,"duration_ms":50922,"significance":"If the central claim holds, this work provides a practical tool for predicting the performance of entanglement-based QKD systems with realistic multi-pair and spectral effects, going beyond the single-pair approximation. The combination of continuous-mode spectral degrees of freedom with multi-pair statistics, the analytic reduction of the covariance expressions, and the testing against two independent systems are notable strengths. The device parameters are largely characterized independently rather than fitted to the key-rate or QBER data. However, the validation is qualitative, the reconstructed JSA phase is not directly tested against phase-sensitive data, and the N=5 truncation for the type-0 system lacks a convergence check, so the strength of the central agreement claim is currently limited.","major_comments":[{"comment":"The complex phase of the phase-matching function is reconstructed from the measured power spectrum under the decoupling ansatz Δk(ω−, z) ≈ Δk(ω−) + δk(z) and the quadratic expansion δk(ξ) = δk′ξ + (1/2)δk″ξ², with δk′, δk″, Δk0, and a background fitted to |Φ(ω−)|². A power spectrum fixes only the modulus of Φ, so any frequency-position coupling not of this separable form would change the JSA phase and, through eqs. (A8)–(A9), the predicted Franson interference and dispersion-induced time-basis QBER. The comparisons in Figs. 6, 7, and 10 do not isolate this phase dependence: Fig. 7 shows mainly the expected cosine fringe shape, and the μ-sweeps in Fig. 10 combine many simultaneous effects. The simulation error bands also do not include uncertainties in the fitted JSA parameters. The authors should provide a sensitivity analysis over the reconstructed phase parameters or validate the phase with a phase-sensitive measurement (e.g., Franson visibility or time-domain correlation) to support the central agreement claim.","section":"IV A 1, eqs. (17) and (20)"},{"comment":"The type-0 simulations truncate the covariance expansion at order N = 5 and rely on the channel-separation condition eq. (24), but no convergence study is given. Equation (24) only ensures that up to order N no same-channel pair contribution occurs; it does not demonstrate that N = 5 is sufficient for the sifted key rate and QBER predictions at the μ values considered. The authors should show results for increasing N (for example N = 3, 5, 7) or otherwise quantify the truncation error.","section":"IV A 2, eq. (24)"},{"comment":"The agreement between simulation and measurement is assessed visually, with best-/worst-case parameter bands, and no quantitative agreement statistic is reported. Because the bands are wide and the horizontal μ error bars are also large, the statement that the simulations 'closely match' or are 'in agreement' is stronger than the evidence supports. The authors should provide a quantitative comparison metric (e.g., reduced chi-square) or explicitly qualify the comparison as a semi-quantitative validation.","section":"IV D, Fig. 10"}],"minor_comments":[{"comment":"The sentence following eq. (39) contains apparent typographical errors: 'where Ie = ∅, Ic = Ie and Il = Ie ∪ Ic' should likely read something like 'Īe = ∅, Īc = Ie, Īl = Ie ∪ Ic'; please correct the notation.","section":"IV C, after eq. (39)"},{"comment":"The text contains a duplicated phrase, 'In section section IV A'; please revise.","section":"IV A, first paragraph"},{"comment":"In eq. (A9), the argument of ψ* contains 'τ (µ′) B', which appears to be a typographical error for a path label such as y′; please check and correct.","section":"Appendix A, eq. (A9)"},{"comment":"The definition of the characteristic loss distance 'L0 = 103/ln(10α·0.1 km)' is difficult to parse; please clarify the intended expression, for example L0 = 10³/(ln(10) α [dB/km] × 0.1 km).","section":"IV B, fiber loss definition"}],"recommendation":"major_revision","confidential_remarks":"The main risk is that the paper's strongest claim—agreement with measured key rates and QBERs—rests on an untested JSA phase model and a fixed truncation order. If the authors can supply the requested convergence study and a phase-sensitivity analysis, the paper would be much stronger. The wide parameter bands also make the visual agreement less convincing than the abstract suggests, and a quantitative comparison would help the reader assess the predictive power."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a paper worth engaging with. It applies the covariance-expansion formalism from Part I to two real BBM92 QKD setups and shows that a simulation with mostly independently characterized parameters reproduces measured sifted key rates and time-basis QBER over several parameter scans. The WDM channel-restriction method (Section IV A 2) and the analytic detection probability formulas in Appendix A (eqs. A8, A9) are genuinely new and useful.\n\nThe validation against type-0 and type-II systems with different distances, repetition rates, and detector settings is the strongest part. Parameters were determined from component test reports, crystal conversion efficiency measurements, and detector tomography, not fitted to the key-rate or QBER data. The worst-/best-case error bands are an honest representation of parameter uncertainty. The paper also explains physically why the time-basis QBER oscillates with distance and why interleaving can lower QBER – that mechanistic insight is valuable.\n\nThe math in the appendix is clean; the reduction from 24x24 to 2x2 via Sylvester's theorem is transparent. Citation pattern is fine, with the method properly building on Part I and the experimental systems being the authors' own published setups.\n\nThe soft spots are real but not fatal. Validation is visual; I would have liked a quantitative goodness-of-fit and a convergence study for the N=5 truncation. The 'beyond single-pair' claim is not quantified against a single-pair baseline. The stress-test point about the JSA phase is fair: the complex phase is reconstructed from a power-spectrum fit under the decoupling ansatz Δk(ω−, z) ≈ Δk(ω−) + δk(z). Since the Franson-interference terms and the dispersion-induced QBER are phase-sensitive, the agreement in Figs. 6, 7, and 10 does not by itself certify that phase model. This is a validation gap, not a demonstrated error; the same phase model is used consistently across many independent measurements. I'd want a direct test of the reconstructed phase, e.g. a phase-sensitive measurement or alternative reconstruction, before fully trusting the predictive claim.\n\nWho is this for? Engineers and experimentalists designing entanglement-based QKD networks who need multi-pair and spectral effects in rate predictions. It deserves a serious referee; I'd send it out with a request for a convergence study, a single-pair baseline comparison, and sensitivity analysis of the phase-reconstruction assumption. Minor: no code or data released.","headline":"Solid applied QKD simulation paper; the reconstructed JSA phase is the main unvalidated load-bearing assumption, but the methods and validation are worth referee time.","tokens_in":20926,"tokens_out":2970,"would_cite":true,"duration_ms":26971,"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":"This paper establishes that frequency-resolved simulations including multi-pair events reproduce measured QKD key rates and quantum bit error rates.","keywords":["quantum key distribution","BBM92 protocol","biphoton states","covariance formalism","joint spectral amplitude","multi-pair effects","chromatic dispersion","single-photon detector imperfections"],"falsifier":"Measure the non-local two-photon interference visibility of the type-II source as a function of added chromatic dispersion and compare it with the simulation's prediction; a systematic deviation growing with dispersion would indicate that the reconstructed complex phase-matching function is wrong. Alternatively, measure the two-photon spectral phase directly using stimulated emission tomography and compare it with the model's phase.","tokens_in":19879,"feed_emoji":"🔐","tokens_out":5766,"duration_ms":50904,"temperature":0.7,"pith_summary":"This paper argues that entanglement-based QKD systems can be simulated accurately without neglecting multi-pair events, by expanding the covariance matrix and detection probabilities of the biphoton state in powers of its joint spectral amplitude. It applies that expansion, developed in Part I, to two SPDC-based implementations of the time-bin BBM92 protocol: a type-II polarized system and a type-0 system with wavelength-division demultiplexing. The simulation includes the measured photon spectra, chromatic dispersion in optical fibers, interferometer mode mismatch, and detector dark counts, dead times, and afterpulses. Across varied mean pair number, fiber distance, repetition rate, and detector settings, the simulated sifted key rates and time-basis quantum bit error rates agree with measurements. If this holds, the method provides a predictive design tool for choosing optimal operating parameters of entanglement-based QKD networks.","feed_headline":"Full-spectrum QKD simulation matches measured key rates and QBER","feed_subtitle":"Multi-pair events, fiber dispersion, and detector dead times included in a time/frequency-resolved BBM92 model.","key_machinery":"The machinery is the covariance-matrix expansion: the renormalized covariance $\\Gamma$ of a Gaussian biphoton state is approximated by $\\Gamma_N = \\sum_{n=1}^N (2Z)^n/(2n!)$, with $Z$ built from the joint spectral amplitude $\\psi$, and detection probabilities follow from the probability-generating function $\\det(1+W\\Gamma)^{-1/2}$. For the type-0 system, the key object is the projected reduced JSA $\\tilde{\\psi} = P\\psi P$, where $P$ selects only frequencies within $N\\Delta_+/2$ of the wavelength-division channel bounds, so that the Schmidt decomposition of the reduced state can be computed directly. For the type-II system, the bivariate Poisson approximation yields analytic expressions for the detection probabilities in terms of the reconstructed JSA. These objects carry the argument because they turn continuous frequency entanglement and multi-pair events into finite matrix computations whose outputs are directly compared with measured key rates and QBER.","core_discovery":"The central claim is that the full photon statistics of highly entangled biphoton states, including multi-pair emission, spectral-temporal correlations, and realistic detector behavior, can be captured by a truncated series expansion of the renormalized covariance, and that the resulting simulations quantitatively reproduce measured QKD performance. For the type-II system, the complex phase-matching function is reconstructed from measured power spectra through a decoupling approximation. For the type-0 system, the joint spectral amplitude is restricted by a projection onto the frequency components that can reach the wavelength-division channels, making higher-order expansions numerically tractable. The predicted sifted key rates and time-basis QBER match experimental data for both systems.","pith_inferences":["A natural next step the paper does not take is to invert the simulation: with a fast surrogate, one could search for the globally optimal combination of pump power, repetition rate, fiber length, and detector settings rather than scanning best-case and worst-case parameter bounds.","The wavelength-division channel-restriction argument suggests a general rule of thumb: for highly entangled states, only frequencies within roughly $N\\Delta_+/2$ of the channel edges matter at expansion order $N$, which could be used to bound truncation error a priori.","The method's reliance on a reconstructed joint spectral amplitude implies an independent experimental check: measuring non-local two-photon interference visibility as a function of added dispersion would isolate phase-reconstruction errors from other imperfections."],"forward_implications":["The simulation reproduces measured sifted key rates and time-basis QBER for both SPDC systems, so it can predict performance at parameter values that have not yet been measured.","Because multi-pair events, fiber dispersion, and detector imperfections are included, the simulation can identify the dominant error source at each operating point, such as detector saturation from dead time at high mean pair number.","The predicted non-monotonic time-basis QBER as a function of fiber distance, caused by spectral side lobes and chromatic dispersion, is a directly observable signature of the model.","The method extends to other biphoton-based quantum communication protocols, since the covariance expansion and the detection model are not specific to BBM92."],"supporting_citations":[{"why":"Supplies the covariance expansion and bivariate Poisson approximation on which the entire simulation rests.","marker":"[26]"},{"why":"Describes the four-user time-bin BBM92 experimental system whose key-rate and QBER measurements are used for validation.","marker":"[27]"},{"why":"Provides the virtual beam-splitter model used for interferometer mode mismatch in the receiver.","marker":"[15]"},{"why":"Establishes the generating-function and detector-tomography approach that the detection-probability calculations build on.","marker":"[16]"},{"why":"Supplies the measured dark-count, dead-time, and afterpulse parameters of the single-photon detectors via detector tomography.","marker":"[34]"},{"why":"Provides the measured photon-pair spectra used to reconstruct the joint spectral amplitudes of the SPDC sources.","marker":"[10]"},{"why":"Underlies the decoupling approximation used to reconstruct the complex phase-matching function from measured power spectra.","marker":"[40]"}],"fun_headline_variants":["Multi-pair QKD simulation matches measured key rates and QBER","Full photon statistics in QKD sim reproduces experimental results","Time-resolved biphoton QKD model verified by measured QBER","Covariance expansion captures multi-pair effects in QKD"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predicted interference and error rates rely on the assumption that the measured photon-pair spectrum's asymmetry is fully explained by a phase mismatch that separates into a uniform frequency-dependent part and a smoothly varying crystal perturbation; if the true phase couples frequency and position differently, all predicted visibilities and error rates shift.","fun_headline_variants_meta":{"raw":{"variants":["Multi-pair QKD simulation matches measured key rates and QBER","Full photon statistics in QKD sim reproduces experimental results","Time-resolved biphoton QKD model verified by measured QBER","Covariance expansion captures multi-pair effects in QKD"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1312,"prompt_tokens":802,"completion_tokens":510,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":418,"completion_tokens_details":{"reasoning_tokens":437}},"tokens_in":418,"tokens_out":510,"duration_ms":5138,"temperature":1.0,"reasoning_tokens":437,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:49:46.238954+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the non-local two-photon interference visibility of the type-II source as a function of added chromatic dispersion and compare it with the simulation's prediction; a systematic deviation growing with dispersion would indicate that the reconstructed complex phase-matching function is wrong. Alternatively, measure the two-photon spectral phase directly using stimulated emission tomography and compare it with the model's phase.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the covariance expansion and bivariate Poisson approximation on which the entire simulation rests."},{"cited_title":"Dorfman, S","cited_arxiv_id":null,"evidence_quote":"Describes the four-user time-bin BBM92 experimental system whose key-rate and QBER measurements are used for validation."},{"cited_title":"Olivares, Quantum optics in the phase space, Eur","cited_arxiv_id":null,"evidence_quote":"Provides the virtual beam-splitter model used for interferometer mode mismatch in the receiver."},{"cited_title":"Brendel, N","cited_arxiv_id":null,"evidence_quote":"Supplies the measured dark-count, dead-time, and afterpulse parameters of the single-photon detectors via detector tomography."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the measured photon-pair spectra used to reconstruct the joint spectral amplitudes of the SPDC sources."},{"cited_title":"Helmfrid and G","cited_arxiv_id":null,"evidence_quote":"Underlies the decoupling approximation used to reconstruct the complex phase-matching function from measured power spectra."}],"review_version":1}