{"id":"67f72857-fc66-4583-a6b8-be5bf6c98203","arxiv_id":"2501.15148","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For CubeSat downlink QKD modeled with an elliptical beam channel, efficient BB84 consistently beats standard BB84 in both finite-size and asymptotic key rates.","lead":"This paper simulates secure key rates for two versions of the BB84 quantum key distribution protocol on a CubeSat-to-ground link, using an elliptical beam model of the turbulent atmosphere. It finds that the efficient BB84 variant generates more key bits than the standard variant across all tested weather and daytime conditions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"As printed, Eq. (15) is not a valid transmittance for a two-dimensionally wandering elliptical beam; it depends only on a scalar centroid coordinate and omits the y0 term, so the quantitative key-rate results in Figs. 2–3 are not supported by the stated model.","rationale":"The reader's strongest claim is the quantitative demonstration that efficient BB84 beats standard BB84, e.g., 8e-5 vs 4e-5 bits per pulse. The condition for that claim is that the transmittance model and the Monte-Carlo averaging are correct. The printed Eq. (15) has a clear internal inconsistency: the integrand should depend on the full centroid (x0,y0), but as written it depends on a single scalar ρ0 and has no y0 term. This makes the relationship between the sampled beam parameters and the resulting η ill-defined, and therefore the numerical rates and PDRs are not reproducible from the manuscript. The unspecified normal distributions in Appendix A are another gap in the same pipeline, but the Eq. (15) issue is more load-bearing because it would invalidate the results even after the distributions are specified. Both are fixable by providing the correct formula and full parameter distributions, so the appropriate editorial action remains conditional acceptance rather than outright rejection: the qualitative efficient-vs-standard comparison is structurally expected from the sifting advantage and unlikely to flip, but the paper's quantitative contribution cannot be assessed until the model is corrected and specified. Thus I do not change the reader's CONDITIONAL verdict; I add a more concrete, higher-priority condition.","tokens_in":17059,"tokens_out":10719,"duration_ms":99798,"concrete_test":"Re-derive Eq. (15) from the elliptical-Gaussian beam intensity and compare with the published formula in Liorni et al. [27] or Vasylyev et al. [28]. Then evaluate the transmittance for (x0=ra/2, y0=0) and (x0=0, y0=ra/2) with W1=W2 and arbitrary α: a correct isotropic circular beam must give identical values by rotational symmetry, while Eq. (15) as printed depends only on x0 (through ρ0) and cannot satisfy this. If the authors' implementation reproduces the symmetry, Eq. (15) needs rewriting; if it does not, all rates and PDRs in Figs. 2 and 3 must be recomputed with the corrected transmittance.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Every key-rate number in Figs. 2 and 3 is obtained by Monte-Carlo sampling of the beam parameter vector v=(x0,y0,W1,W2,α) and evaluating Eq. (15) (see Fig. 2 text and Eq. 16). The printed Eq. (15) contains the centroid only through the scalar ρ0 in (ρ cosθ - ρ0), with no (ρ sinθ - y0) dependence. Expanding the intensity of an elliptical Gaussian centered at (x0,y0) gives -2[A1(ρ cosθ - x0)^2 + A2(ρ sinθ - y0)^2 + A3(ρ cosθ - x0)(ρ sinθ - y0)] in the exponent (with x0=ρ0 cosθ0, y0=ρ0 sinθ0), so Eq. (15) as written cannot be correct for a beam displaced in both transverse directions. If the code follows the printed formula, the simulated transmittance depends on the choice of coordinate origin; if the code uses the correct two-dimensional formula, the manuscript has a serious transcription error in the central model equation. Either way, the specific quantitative claim (e.g., ~8e-5 vs ~4e-5 bits/pulse) is not currently reproducible. The reader's point about the unspecified normal means and variances for x0, y0, Θ1, and Θ2 in Appendix A compounds this: even with a corrected Eq. (15), the sampling distribution is undefined, so the PDT and PDR curves cannot be independently reproduced.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a simulation study of finite and asymptotic secret key rates for the efficient BB84 and standard decoy-state BB84 protocols in a CubeSat-to-ground QKD downlink. A free-space channel is modeled with an elliptical Gaussian beam approximation, including turbulence-induced beam wandering and weather-dependent attenuation. The authors use Monte Carlo sampling of beam parameters to compute the probability distribution of transmittance (PDT) and then average the key-rate expressions from Refs. [17,21] over this distribution. They report that efficient BB84 consistently outperforms standard BB84 across six weather conditions and provide probability distributions of key rates (PDR) at selected zenith angles.","tokens_in":17460,"tokens_out":7761,"duration_ms":61918,"significance":"If the quantitative results are correct, the paper would provide a practical, comparative engineering assessment of two BB84 variants for CubeSat QKD, with explicit weather dependence. A clear strength is that all key-rate formulas and channel parameters are taken from published work; no constants are fit to the paper's own results, so there is no circularity. The qualitative ranking of the two protocols (efficient BB84 > standard BB84) is robust and follows directly from the biased-basis structure of efficient BB84. However, the quantitative key rates are presently not reproducible because the central transmittance formula (Eq. (15)) is not the transmittance of a two-dimensionally wandering elliptical beam, and the Monte Carlo sampling distributions in Appendix A are incompletely specified. These issues affect every number in Figs. 2 and 3.","major_comments":[{"comment":"Equation (15) as printed cannot represent an elliptical Gaussian beam displaced by (x0,y0). The beam centroid enters only through the scalar ρ0 in (ρcosθ - ρ0) and (ρcosθ - ρ0)ρsinθ, and the expression has no dependence on y0 when θ0=0. For a beam centered at (ρ0cosθ0, ρ0sinθ0), the exponent must contain (ρcosθ - x0)^2, (ρsinθ - y0)^2, and the cross term (ρcosθ - x0)(ρsinθ - y0). The printed formula is therefore either a misprint or an incorrect model. Because the text (page 13) states that all points in Figs. 2 and 3 are computed from Eq. (15), the simulated transmittances and all resulting key rates are not reproducible from the manuscript as written.","section":"§2.3, Eq. (15)"},{"comment":"The Monte Carlo sampling distributions are under-specified. The manuscript states that (α - θ0) is uniform on [0,π/2] and that x0, y0, Θ1, Θ2 follow normal distributions, but it does not give the means or variances of those normals. Since Wi is derived from Θi via Θi = ln(Wi^2/W0^2), the lognormal parameters of the semi-axes are also unspecified. Every average key rate in Figs. 2 and 3 is an expectation over these draws (Eqs. (15)-(16)), so the quantitative results cannot be reproduced or cross-checked without these parameter values. The authors should also clarify whether the same distribution parameters are used for all six weather conditions or are scaled with Cn^2 and n0.","section":"Appendix A (after Eq. (18))"}],"minor_comments":[{"comment":"The text states the zenith angle is restricted to [0°,66°], but Fig. 2c,d plot to 80° and the text reports finite-key cutoff at 62° and asymptotic cutoff at 75.2°; please reconcile the range.","section":"Section 2.3 (p.10)"},{"comment":"Equation (3) contains an unbalanced parenthesis and a garbled term '− (12 log2 21 εsec − 2 log2 2 εcorr'; rewrite it in the same style as Eq. (1).","section":"Eq. (3)"},{"comment":"Table 2 lists 'Error correction efficiency' as a parameter but does not give a value, and λec is said to 'depend on block size'; since Eq. (2) determines λec from nX and Q, please state explicitly whether an additional reconciliation-efficiency factor is used and what value it takes.","section":"Table 2"},{"comment":"Fig. 2a,b show zenith angles from 10° upward, while the text reports key rates 'at the zenith position' (0°); please include the 0° point or adjust the text.","section":"Fig. 2"},{"comment":"Since the plotted key rates are Monte Carlo estimates (1,000 samples for Fig. 2, 30,000 for Fig. 3), the absence of error bars or confidence intervals makes it difficult to judge whether the small differences between protocols are significant.","section":"Figs. 2 and 3"},{"comment":"References [10] and [18] are the same article (Ecker et al., npj Quantum Information 7:5, 2021); please merge them.","section":"References [10] and [18]"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely and well-scoped problem, and the qualitative recommendation (efficient BB84 for CubeSat downlinks) is likely sound. My main concern is the reproducibility of the quantitative results, which stems from an apparent error in Eq. (15) and the absence of the sampling distribution parameters in Appendix A. These are fixable in revision. The paper would also benefit from a closer proofread of the equations and a clarification of the zenith-angle range."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is an application paper, not a theory paper. The finite-key formulas come from Sidhu et al. and Lim et al., the elliptical beam model from Liorni et al., and the qualitative conclusion — efficient BB84 beats standard BB84 — follows directly from the biased-basis protocol structure. That conclusion is robust, and the paper does a legitimate job of extending these tools to a CubeSat downlink with six weather conditions and a probability-distribution-of-key-rate analysis.\n\nThe problem is the central quantitative engine. Eq. (15), as printed, does not look like a valid two-dimensional transmittance for an elliptical beam wandering in both x and y. It depends on the centroid only through a scalar ρ0 in (ρ cosθ − ρ0), with no (ρ sinθ − y0) term. A beam displaced in both directions would have both coordinates in the exponent. So either the code implements something different from what is written, or the simulation is using a formula that is not the stated model. Either way, the exact numbers in Figs. 2 and 3 — e.g., the 8e−5 vs 4e−5 bits per pulse figures — are not reproducible from the manuscript. Appendix A compounds this: it says x0, y0, Θ1, Θ2 are normal but never gives their means or variances, so the sampling distribution is undefined. There is also no circular-beam baseline, despite the abstract claiming the elliptical model is more accurate. That is a missing control.\n\nA few smaller issues: no error bars on any plot, and a couple of transcription slips in the equations. These are minor compared to the main problem.\n\nWho is this for? Engineers or mission planners working on CubeSat QKD who want a realistic ballpark. They can take the qualitative message — use efficient BB84 — but they should not trust the quantitative rates as printed. The authors need to correct the transmittance formula, specify the distribution parameters, and re-run the plots before this is publishable.\n\nShould a serious editor send this to review? Yes, because the framework is credible and the flaws are addressable. But a referee should flag the transmittance equation immediately and require major revision. I would not accept it in its current form, and I would not cite the numbers.","headline":"The qualitative takeaway is sound but the quantitative results rest on a suspect transmittance formula and underspecified simulation parameters; worth a referee only if that equation gets fixed.","tokens_in":17954,"tokens_out":2775,"would_cite":false,"duration_ms":27278,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81P94"],"pacs":["03.67.Dd"],"model":"deepseek-v4-flash","headline":"For CubeSat QKD downlinks, the efficient BB84 protocol consistently yields higher secret key rates than standard BB84 across all six modeled weather conditions, in both finite and asymptotic regimes.","keywords":["quantum key distribution","BB84 protocol","decoy states","finite-key analysis","asymptotic key rate","CubeSat","elliptical beam model","atmospheric turbulence"],"falsifier":"Recompute the average key rates using explicit, published normal-distribution parameters for the beam centroid and widths, and compare the predicted zenith dependence (for example, $8 \\times 10^{-5}$ bits per pulse for efficient BB84 at zenith in clear daytime) against measured secret key lengths from an actual CubeSat downlink at 400 km altitude, 785 nm wavelength, and a 50 cm ground aperture; if the measured rates fall below the predicted curves or standard BB84 outperforms efficient BB84, the central claim would be falsified.","tokens_in":16854,"feed_emoji":"🛰️","tokens_out":3131,"duration_ms":29239,"temperature":0.7,"pith_summary":"This paper tries to establish that, for a CubeSat-based quantum key distribution downlink, the efficient BB84 protocol is practically preferable to standard BB84. Using an elliptical beam model of atmospheric transmittance, the authors compute finite-block and asymptotic secret key rates per pulse for both protocols under clear, foggy, and windy day and night conditions. They find that efficient BB84 consistently outperforms standard BB84, roughly doubling the key rate at zenith, and that the advantage persists at larger zenith angles. A sympathetic reader would care because CubeSats are a leading low-cost platform for global QKD, and choosing the right protocol and accurately modeling the channel are key to whether such systems can generate usable keys.","feed_headline":"Biased-basis BB84 doubles CubeSat quantum key rates","feed_subtitle":"Finite and asymptotic simulations across six weather conditions favor the efficient protocol for practical downlinks.","key_machinery":"The load-bearing object is the elliptical beam transmittance model, Eq. (15), which gives the transmittance $\\eta(x_0,y_0,W_1,W_2,\\alpha)$ through a circular receiving aperture as a function of beam centroid, semi-axes, and orientation. This transmittance is averaged over the probability distribution of transmittance (PDT) via Eq. (16), using Monte Carlo samples of the beam parameters. The key rate formulas, Eqs. (1), (3), (13), and (14), then convert transmittance samples into finite and asymptotic secret key rates for efficient and standard BB84, with finite statistics handled by multiplicative Chernoff bounds.","core_discovery":"The central claim is that, in the modeled CubeSat downlink, efficient BB84 with two decoy states yields higher finite and asymptotic average secret key rates than standard BB84 under every weather condition considered, with the gap widening as turbulence and fog increase. The mechanism is the biased basis choice of efficient BB84, which increases the sifting ratio and reduces the statistical uncertainty from two-basis parameter estimation. Quantitatively, at zenith in clear daytime, the finite key rate is about $8 \\times 10^{-5}$ bits per pulse for efficient BB84 versus about $4 \\times 10^{-5}$ for standard BB84, and asymptotic rates are roughly 1.2 to 1.5 times higher than finite rates. The paper also analyzes the probability distribution of key rates across zenith angles, showing broader distributions at low zenith angles and a consistent qualitative advantage for efficient BB84.","pith_inferences":["Editorial inference: the quantitative key rates in the paper depend on the unspecified means and variances of the normal distributions used to sample $x_0$, $y_0$, $\\Theta_1$, and $\\Theta_2$; specifying these would make the results reproducible and testable.","Editorial inference: a natural testable extension is to replace the normal and uniform beam-parameter distributions with log-normal, Gamma-Gamma, or Double Weibull transmittance models and check whether the efficient-BB84 advantage persists, as the paper itself suggests as future work.","Editorial inference: because the efficient protocol's advantage stems from its higher sifting ratio, real hardware with imperfect basis selection or additional background noise may reduce, but likely not eliminate, the gap shown here.","Editorial inference: if real CubeSat missions collect secret key length data, comparing measured rates to the predicted zenith-dependent curves (e.g., $8 \\times 10^{-5}$ bits/pulse at zenith in clear daytime) would provide a direct validation of the elliptical beam model."],"forward_implications":["If the claim holds, efficient BB84 should be the protocol of choice for CubeSat-based QKD downlinks, especially when turbulence and fog limit transmittance.","The modeled advantage implies that a CubeSat QKD mission using efficient BB84 could generate roughly twice the secret key length of standard BB84 in a single overpass at zenith.","The finite-key analysis shows that key generation stops at a lower zenith angle (about $62^\\circ$) than the asymptotic limit (about $75.2^\\circ$), so mission design must account for finite statistics when planning pass geometry.","The consistent weather ordering (Day 1, Night 1, Day 2, Night 2, Day 3, Night 3) suggests that clear daytime conditions are optimal for CubeSat downlink QKD, while foggy and windy conditions degrade rates more severely for standard BB84.","The elliptical beam model, by including beam wandering and elliptical deformation, provides a more realistic transmittance distribution for designing ground-station aperture sizes and link budgets."],"supporting_citations":[{"why":"Supplies the elliptical beam approximation and weather-dependent channel parameters used to model transmittance in the downlink.","marker":"[27]"},{"why":"Provides the finite and asymptotic key rate formulas and the satellite pass time window that the paper adapts to CubeSats.","marker":"[17]"},{"why":"Gives the concise security bounds and phase error rate estimation used in the finite-key analysis.","marker":"[21]"},{"why":"Provides the tight multiplicative Chernoff bounds used for finite statistics corrections.","marker":"[43]"},{"why":"Supplies the error correction leakage bound and its asymptotic form used in the key rate formulas.","marker":"[44]"},{"why":"Defines the efficient BB84 protocol with biased basis selection that is the main protocol under study.","marker":"[34]"},{"why":"Establishes the decoy-state method that underlies the two-decoy-state standard BB84 analysis.","marker":"[35]"},{"why":"Provides the atmospheric quantum channel model with weak and strong turbulence that the elliptical beam approximation builds on.","marker":"[28]"}],"fun_headline_variants":["Efficient BB84 outperforms standard in CubeSat QKD","Biased-basis BB84 doubles CubeSat key rates","CubeSat QKD: efficient BB84 wins in all weather","Foggy skies favor efficient BB84 for CubeSat keys","Elliptical beam model shows efficient BB84 advantage"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical results assume a specific but incompletely specified probability distribution for the random beam parameters: the angle difference $\\alpha - \\theta_0$ is uniform on $[0, \\pi/2]$, while $x_0$, $y_0$, $\\Theta_1$, and $\\Theta_2$ are said to be normal, yet their means and variances are never stated.","fun_headline_variants_meta":{"raw":{"variants":["Efficient BB84 outperforms standard in CubeSat QKD","Biased-basis BB84 doubles CubeSat key rates","CubeSat QKD: efficient BB84 wins in all weather","Foggy skies favor efficient BB84 for CubeSat keys","Elliptical beam model shows efficient BB84 advantage"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000261,"raw_usage":{"total_tokens":1571,"prompt_tokens":901,"completion_tokens":670,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":517,"completion_tokens_details":{"reasoning_tokens":587}},"tokens_in":517,"tokens_out":670,"duration_ms":6203,"temperature":1.0,"reasoning_tokens":587,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:34:05.499788+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the average key rates using explicit, published normal-distribution parameters for the beam centroid and widths, and compare the predicted zenith dependence (for example, $8 \\times 10^{-5}$ bits per pulse for efficient BB84 at zenith in clear daytime) against measured secret key lengths from an actual CubeSat downlink at 400 km altitude, 785 nm wavelength, and a 50 cm ground aperture; if the measured rates fall below the predicted curves or standard BB84 outperforms efficient BB84, the central claim would be falsified.","supporting_citations":[{"cited_title":"Satellite-based links for quantum key distribution: beam effects and weather dependence","cited_arxiv_id":null,"evidence_quote":"Supplies the elliptical beam approximation and weather-dependent channel parameters used to model transmittance in the downlink."},{"cited_title":"Finite key effects in satellite quantum key distribution","cited_arxiv_id":null,"evidence_quote":"Provides the finite and asymptotic key rate formulas and the satellite pass time window that the paper adapts to CubeSats."},{"cited_title":"Con- cise security bounds for practical decoy-state quantum key distribution","cited_arxiv_id":null,"evidence_quote":"Gives the concise security bounds and phase error rate estimation used in the finite-key analysis."},{"cited_title":"Tight security bounds for decoy-state quantum key distribution","cited_arxiv_id":null,"evidence_quote":"Provides the tight multiplicative Chernoff bounds used for finite statistics corrections."},{"cited_title":"Funda- mental finite key limits for one-way information reconciliation in quantum key distribution","cited_arxiv_id":null,"evidence_quote":"Supplies the error correction leakage bound and its asymptotic form used in the key rate formulas."},{"cited_title":"Efficient quantum key distribution scheme and a proof of its unconditional security","cited_arxiv_id":null,"evidence_quote":"Defines the efficient BB84 protocol with biased basis selection that is the main protocol under study."},{"cited_title":"Decoy state quantum key distribution","cited_arxiv_id":null,"evidence_quote":"Establishes the decoy-state method that underlies the two-decoy-state standard BB84 analysis."},{"cited_title":"Atmospheric quantum channels with weak and strong turbulence","cited_arxiv_id":null,"evidence_quote":"Provides the atmospheric quantum channel model with weak and strong turbulence that the elliptical beam approximation builds on."}],"review_version":1}