REVIEW 2 major objections 5 minor 32 references
Tighter Asymptotic Key Rates for Intensity-Correlated Decoy-State QKD via Nonlinear Programming
T0 review · 2 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read Solving the full nonlinear decoy-state parameter-estimation problems before linearising yields tighter key-rate bounds than canonical reference points, and coincidence of the two stages certifies optimality.
desk verdict A solid, reproducible improvement for correlated-source decoy-state QKD; the relaxation worry is answerable, and the main fix is presentational. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the two-stage optimisation workflow. Stage one solves the full Cauchy-Schwarz-constrained nonlinear programs (P1)-(P3) with a smooth interior-point solver; stage two replaces each CS inequality by its first-order Taylor expansion around that solution and solves the resulting linear program. The CS constraints themselves are the mechanism that couple n-photon yields across different intensity settings: for each pair of intensities a,b and photon number n they bound Y_{n,b} (and H_{n,b}) in terms of Y_{n,a} through square-root functions G±(y,τ), where τ encodes the maximum intensity-information leakage over the correlation range. The work that these constraints do is
What would settle it
Find a tuple of yields that satisfies the nonlinear Cauchy-Schwarz inequalities and the linear observation constraints but violates one of the tangent-line inequalities at the reference point produced by the nonlinear solve; that would show the outer LP is cutting into the feasible set and the claimed lower bound is not guaranteed. A cleaner version: solve the same instances with a certified exact method (e.g., a mixed-integer second-order cone formulation) and check whether its optimum matches the two-stage value; a mismatch would mean the certificate is not reliable.
Extended reading notes
Core claim
The central claim is that the quality of an outer-linearisation bound depends on where you linearise, and the best linearisation point is the solution of the actual nonlinear problem. The paper formulates the three decoy-state parameter-estimation tasks—lower bounds on the Z- and X-basis single-photon yields and an upper bound on the single-photon error rate—as nonlinear programs whose only nonlinearities are the Cauchy-Schwarz square-root constraints. Solving these with an interior-point nonlinear solver and then linearising the CS inequalities around that solution gives valid asymptotic key-rate lower bounds. The simulations show these bounds dominate the canonical channel-model-based refe
Load-bearing premise
The outer linearisation must be a valid outer relaxation—the tangent lines must never cut into the set of yield values allowed by the Cauchy-Schwarz constraints—but the paper never proves this and is inconsistent about whether those constraints are convex, so the final lower bound would not be guaranteed if they are not.
Editorial extensions
If this is right
- The two-stage workflow can be dropped into existing CS-based decoy-state security analyses as a replacement for canonical linearisation points, improving key rates without changing the security proof.
- When the nonlinear value and the outer-LP value coincide, the paper's construction certifies that the computed key-rate lower bound is optimal for that parameter-estimation problem.
- For fine-grained truncated-Gaussian correlation models inspired by real hardware, the gains over canonical points are larger exactly when measured statistics deviate from the simple channel model.
- The approach removes the need for iterative re-linearisation heuristics, making the computation faster and reproducible.
- The improvement holds for both the coarse-grained model-independent formulation and the fine-grained model, across the tested distances and correlation ranges.
Reading between the lines
- The certificate logic only works when the nonlinear and linearised optima coincide; when they differ, the nonlinear value is merely a candidate and the outer bound, while valid, may be loose, so a user must check for agreement rather than assume it.
- The same two-stage idea—solve the true nonconvex problem, then linearise around its solution—should transfer to other QKD parameter-estimation settings where the feasible set admits convex outer approximations, such as source flaws beyond intensity drift or measurement-device imperfections.
- A natural testable extension is to apply the workflow to measured correlation data from real hardware rather than simulated data, to see whether the tighter bounds survive in practice; the paper only simulates the Gaussian model based on one hardware system.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper addresses decoy-state BB84 with phase-randomized weak coherent pulses whose intensities are correlated across rounds. It reviews the Cauchy--Schwarz (CS) constraints of Zapatero et al. [10] that couple n-photon yields across intensity settings, and proposes a two-stage numerical workflow: (i) solve the full nonlinear CS-constrained estimation problems P1--P3 with the interior-point solver IPOPT, and (ii) linearize the CS constraints around the IPOPT solution and solve the resulting LP to obtain a certified lower bound on the asymptotic key rate. The paper reports simulations for a coarse-grained model-independent formulation and for a fine-grained truncated-Gaussian model, showing tighter rates than the canonical channel-model-based linearization points, and argues that when the two optimisation stages coincide the rate is certified optimal. A public code repository is provided for reproducibility.
Significance. The proposed two-stage workflow is a simple, practical improvement over canonical linearization: it preserves the security guarantee of an outer LP relaxation while potentially yielding tighter bounds. Public code availability is a genuine strength, and the optimality-certificate idea is logically sound provided the outer linearization is a true relaxation and the IPOPT point is feasible. The paper's central claim is defensible, but the manuscript leaves the relaxation property unproved and contains an inconsistent and incorrect statement that the CS constraints are non-convex; the relevant feasible set is in fact convex. With fixes to these points, this is a useful, though methodologically incremental, contribution to decoy-state QKD with correlated sources.
major comments (2)
- [Sec. II.I, Eq. (27)] The assertion that the outer linearisation 'always leads to valid key rates' is load-bearing but never proved. Validity requires that tangent lines to G- underestimate G- and tangent lines to G+ overestimate G+, i.e. that G- is convex and G+ concave. This is true (for fixed tau, g-'' = +sqrt(tau(1-tau))/(2[y(1-y)]^{3/2}) > 0 and g+'' = -... < 0, and the branch definitions in Eq. (21) glue smoothly), but the paper should state and prove this explicitly. The current text is inconsistent: the Introduction calls the CS constraints 'convex functions', while Sec. III.A calls them 'non-convex square-root nonlinearities'. Without a correct and explicit proof of the relaxation property, the 'certified lower bound' claim is unsupported.
- [Sec. III, optimality certificate] The claim that coincidence of the two optimisation stages certifies optimality depends on the IPOPT solution being feasible for the original nonlinear problem. If the IPOPT point is only approximately feasible, its objective value is not a valid upper bound, and the sandwich argument (relaxed LP value <= true optimum <= IPOPT value) fails. The paper does not report IPOPT feasibility tolerances or maximum constraint violations. Please report the IPOPT tolerances used and, for the displayed curves, verify and state the maximum constraint violation of the reference points, e.g. max violation < 1e-8.
minor comments (5)
- [Figs. 1 and 2] The captions and text do not make clear whether the curves labeled 'IPOPT' are the raw IPOPT objective values or the final certified outer-LP values after linearizing at the IPOPT point. Since the optimality certificate rests on comparing these two values, the figures should show both (or state that they are numerically indistinguishable and provide the max difference).
- [Sec. II.H title] 'Cauchy-Schwartz' should be 'Cauchy-Schwarz'.
- [Eq. (29)] The derivative g'+- has a singularity at y=0 and y=1. The text should state that reference points are taken in (0,1) or describe how boundary cases are handled, especially because IPOPT could return values at the bounds.
- [Eq. (44)] The notation 'min ... (= Z_{1,mu})' inside the optimization is confusing. It should be written as an objective equal to Z_{1,mu} with the equality displayed after the minimization, not inside the min.
- [References] Reference [3] is given as 'AIAA Journal 61, 1895 (2023)'; this journal title appears incorrect for a quantum-key-distribution paper. Please verify the source.
Circularity Check
No significant circularity: the IPOPT-assisted outer-linearization workflow is a genuine relaxation-and-certificate scheme, not a fit or self-citation chain.
full rationale
The claimed derivation chain is: keep the nonlinear CS constraints (Eqs. 20-23) from [10]; solve the resulting NLPs (P1)-(P3) with IPOPT; use the returned point only as the reference point for the first-order outer linearization (Eq. 27); solve the resulting LP; report the LP value as a certified lower bound. The IPOPT point is not a parameter fitted to make the key-rate prediction come out; it is a candidate feasible point for the same estimation problem that defines the bounds feeding Eq. (6). The final LP bound is not forced to equal the IPOPT objective by construction - the LP can return a strictly smaller value, and optimality is claimed only when the two values coincide, which is the standard relaxation-gap certificate (provided Eq. (27) is a true outer approximation). The outer-relaxation property is not proven in the paper and Sec. III.A even calls the constraints "non-convex square-root nonlinearities", but this is a correctness/clarity risk, not a circular reduction: the required convexity/concavity of G-/G+ can be checked directly from Eq. (22), independent of the paper's conclusions. The only self-citation is [7] in a general background list; it is not load-bearing. No step defines a predicted quantity in terms of an input, no fitted parameter is renamed a prediction, and no uniqueness or ansatz is imported from the authors' own prior work.
Assumptions & free parameters
assumptions (7)
- domain assumption Poisson photon-number statistics conditional on actual intensity (Eq. 1)
- domain assumption Bounded relative deviation |1 - α_k/a_k| ≤ δ_max (Eq. 2)
- domain assumption Finite correlation range ξ (Assumption 3)
- domain assumption Perfect phase randomization, no state-preparation flaws beyond intensity drift, no side-channels (Sec. II.A)
- domain assumption Basis-independent detection efficiency (M^{Z,f}=M^{X,f}, Sec. II.A)
- standard math Outer linearization of CS constraints is a valid relaxation (Sec. II.I, Eq. 27): G- convex, G+ concave
- domain assumption Key-rate formula Eq. (6) is valid
Cite this review
Pith. "Pith review of Tighter Asymptotic Key Rates for Intensity-Correlated Decoy-State QKD via Nonlinear Programming." pith.science (2026). https://pith.science/paper/3GC2CWY7
@misc{pith2026260204966,
author = {Pith},
title = {Pith review of: Tighter Asymptotic Key Rates for Intensity-Correlated Decoy-State QKD via Nonlinear Programming},
year = {2026},
howpublished = {\url{https://pith.science/paper/3GC2CWY7}},
note = {Machine review of arXiv:2602.04966}
}
abstract
Decoy-state QKD with phase-randomized weak coherent pulses is typically analyzed assuming independent, precisely prepared intensities. Real sources, however, can exhibit correlated intensity drift across rounds, potentially leaking intensity information and breaking the standard decoy-state reduction to linear programs. Cauchy--Schwarz (CS) constraints can restore security by coupling $n$-photon yields across intensities, but they introduce nonlinear square-root constraints that are commonly handled via outer linearisation around channel-model-based reference points. We propose a reproducible alternative: first solve the full CS-constrained parameter-estimation problems using the interior-point nonlinear solver IPOPT, then use the resulting candidate solution as the linearisation point for the outer optimisation that certifies a valid lower bound on the asymptotic key rate. Simulations for both coarse-grained model-independent correlations and fine-grained truncated-Gaussian models show consistently tighter key-rate bounds than canonical reference points, and in some cases allow certifying optimality when both optimisation stages coincide.
Figures
Reference graph
Works this paper leans on
-
[10]
V. Zapatero, ´A. Navarrete, K. Tamaki, and M. Curty, Security of quantum key distribution with intensity correlations, Quantum5, 602 (2021), arXiv:2105.11165 [quant-ph]
arXiv 2021
-
[28]
Trefilov, X
D. Trefilov, X. Sixto, V. Zapatero, A. Huang, M. Curty, and V. Makarov, Intensity correlations in decoy-state bb84 quantum key distribution systems, Optica Quantum3, 417 (2025)
2025
-
[1]
G. S. Vernam, Cipher printing telegraph systems: For secret wire and radio telegraphic communications, Journal of the A.I.E.E.45, 109 (1926)
1926
-
[2]
M. N. Wegman and J. L. Carter, New hash functions and their use in authentication and set equality, Journal of Computer and System Sciences22, 265 (1981)
1981
-
[3]
Renner and R
R. Renner and R. Wolf, Quantum advantage in cryptography, AIAA Journal61, 1895 (2023)
2023
-
[4]
Wang, X.-L
X.-B. Wang, X.-L. Hu, and Z.-W. Yu, Practical Long-Distance Side-Channel-Free Quantum Key Distribution, Physical Review Applied12, 054034 (2019)
2019
-
[5]
Nahar, T
S. Nahar, T. Upadhyaya, and N. L¨ utkenhaus, Imperfect phase randomization and generalized decoy-state quantum key distribution, Physical Review Applied20, 064031 (2023)
2023
-
[6]
G. Curr´ as-Lorenzo, M. Pereira, G. Kato, M. Curty, and K. Tamaki, A security framework for quantum key distribution implementations (2023), arXiv:2305.05930 [quant-ph]
arXiv 2023
Show all 32 references
-
[7]
Pereira, G
M. Pereira, G. Curr´ as-Lorenzo, and M. Ara´ ujo, Optimal key rates for quantum key distribution with partial source characterization (2025), arXiv:2510.13085 [quant-ph]
2025
-
[8]
Nagamatsu, A
Y. Nagamatsu, A. Mizutani, R. Ikuta, T. Yamamoto, N. Imoto, and K. Tamaki, Security of quantum key distribution with light sources that are not independently and identically distributed, Physical Review A93, 042325 (2016)
2016
-
[9]
Pereira, G
M. Pereira, G. Kato, A. Mizutani, M. Curty, and K. Tamaki, Quantum key distribution with correlated sources, Science Advances6, eaaz4487 (2020), arXiv:1908.08261 [quant-ph]
2020 arXiv
-
[11]
Curr´ as-Lorenzo, S
G. Curr´ as-Lorenzo, S. Nahar, N. L¨ utkenhaus, K. Tamaki, and M. Curty, Security of quantum key distribution with imperfect phase randomisation, Quantum Science and Technology9, 015025 (2023)
2023
-
[12]
Pereira, G
M. Pereira, G. Curr´ as-Lorenzo, A. Mizutani, D. Rusca, M. Curty, and K. Tamaki, Quantum key distribution with un- bounded pulse correlations, Quantum Science and Technology10, 015001 (2024)
2024
-
[13]
Marwah and F
A. Marwah and F. Dupuis, Proving security of BB84 under source correlations (2024), arXiv:2402.12346 [quant-ph]
2024 arXiv
-
[14]
Li, Y.-G
J.-X. Li, Y.-G. Shan, R. Wang, F.-Y. Lu, Z.-Q. Yin, S. Wang, W. Chen, D.-Y. He, G.-C. Guo, and Z.-F. Han, Secure quantum key distribution against correlated leakage source (2025), arXiv:2507.11251 [quant-ph]
2025 arXiv
-
[15]
Mizutani, G
A. Mizutani, G. Kato, K. Azuma, M. Curty, R. Ikuta, T. Yamamoto, N. Imoto, H.-K. Lo, and K. Tamaki, Quantum key distribution with setting-choice-independently correlated light sources, npj Quantum Information5, 1 (2019)
2019
-
[16]
Sixto, V
X. Sixto, V. Zapatero, and M. Curty, Security of Decoy-State Quantum Key Distribution with Correlated Intensity Fluc- tuations, Physical Review Applied18, 044069 (2022)
2022
-
[17]
Li, F.-Y
J.-X. Li, F.-Y. Lu, Z.-H. Wang, V. Zapatero, M. Curty, S. Wang, Z.-Q. Yin, W. Chen, D.-Y. He, G.-C. Guo, and Z.-F. Han, Quantum key distribution overcoming practical correlated intensity fluctuations (2025), arXiv:2501.13482 [quant-ph]
2025 arXiv
-
[18]
Burenkov, B
V. Burenkov, B. Qi, B. Fortescue, and H.-K. Lo, Security of high speed quantum key distribution with finite detector dead time (2010), arXiv:1005.0272 [quant-ph]
2010 arXiv
-
[19]
Kamin and N
L. Kamin and N. L¨ utkenhaus, Improved Decoy-state and Flag-state Squashing Methods (2024), arXiv:2405.05069 [quant- ph]
2024 arXiv
-
[20]
Tupkary, S
D. Tupkary, S. Nahar, P. Sinha, and N. L¨ utkenhaus, Phase error rate estimation with basis-efficiency mismatch for decoy- state BB84 (2024), arXiv:2408.17349 [quant-ph]
2024
-
[21]
Marcomini, A
A. Marcomini, A. Mizutani, F. Gr¨ unenfelder, M. Curty, and K. Tamaki, Loss-tolerant quantum key distribution with detection efficiency mismatch, Quantum Science and Technology10, 035002 (2025)
2025
-
[22]
Grasselli, G
F. Grasselli, G. Chesi, N. Walk, H. Kampermann, A. Widomski, M. Ogrodnik, M. Karpi´ nski, C. Macchiavello, D. Bruß, and N. Wyderka, Quantum key distribution with basis-dependent detection probability, Physical Review Applied23, 044011 (2025)
2025
-
[23]
Nahar, D
S. Nahar, D. Tupkary, and N. L¨ utkenhaus, Imperfect detectors for adversarial tasks with applications to quantum key distribution (2025), arXiv:2503.06328 [quant-ph]
2025
-
[24]
Z. Wang, D. Tupkary, and S. Nahar, Phase error estimation for passive detection setups with imperfections and memory effects (2025), arXiv:2508.21486 [quant-ph]
2025
-
[25]
H.-K. Lo, X. Ma, and K. Chen, Decoy State Quantum Key Distribution, Physical Review Letters94, 230504 (2005)
2005
-
[26]
X. Ma, B. Qi, Y. Zhao, and H.-K. Lo, Practical decoy state for quantum key distribution, Phys. Rev. A72, 012326 (2005)
2005
-
[27]
Rice and J
P. Rice and J. Harrington, Numerical analysis of decoy state quantum key distribution protocols (2009), arXiv:0901.0013 15 [quant-ph]
2009 arXiv
-
[29]
W¨ achter and L
A. W¨ achter and L. T. Biegler, On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming, Mathematical Programming106, 25 (2006)
2006
-
[30]
Yoshino, M
K.-i. Yoshino, M. Fujiwara, K. Nakata, T. Sumiya, T. Sasaki, M. Takeoka, M. Sasaki, A. Tajima, M. Koashi, and A. Tomita, Quantum key distribution with an efficient countermeasure against correlated intensity fluctuations in optical pulses, npj Quantum Information4, 1 (2018)
2018
-
[31]
Gurobi Optimization, LLC, Gurobi Optimizer Reference Manual (2026)
2026
-
[32]
C. Coey, M. Lubin, and J. P. Vielma, Outer approximation with conic certificates for mixed-integer convex problems, Mathematical Programming Computation12, 249 (2020), arXiv:1808.05290 [math.OC]
2020 arXiv
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.