REVIEW 3 major objections 5 minor 18 references
Dual Conic Proxy for Semidefinite Relaxation of AC Optimal Power Flow
T0 review · 3 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read This paper claims that a neural network can be trained to output verified lower bounds for AC optimal power flow from the semidefinite relaxation, without solving an SDP at run time.
desk verdict Genuinely new dual-feasible SDP proxy with clean eigenvalue-shift completion, but the main guarantee borrows unproven recovery rules from the SOC paper and the speedup claim lacks timing data. 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 Algorithm 1, the differentiable dual-completion layer for the dual SDP (called DSDP-OPF in the paper). The layer predicts $\lambda^p,\lambda^q$ and the $\nu$ flow variables, then recovers $\vec\lambda$ from the linear dual equalities, recovers the $\nu^s$ components from the second-order-cone constraints, recovers the generation-bound multipliers from the cost equalities, and finally sets $\hat S = -A_R(\lambda,\mu,\nu) - jA_I(\lambda,\mu,\nu)$, $S = \hat S - \delta I$ with $\delta = \min(0,\lambda_{\min}(\hat S))$, and $\bar\mu^w = 0$, $\underline\mu^w = -\delta$. The shift makes $S$ Hermitian positive semidefinite by construction, and the $\mu^w$ terms exactly compensate the shift in the linear equality, so constraints (6h) and (6k) hold by construction. Differentiability is what matters for training: the objective of the dual is used as the loss, and gradients flow through these recovery steps.
What would settle it
Run Algorithm 1 on a fresh test instance and check the output against every constraint of the dual SDP: the equalities (6b)-(6g), nonnegativity of the $\mu$ variables, membership in the second-order cones (6j), $S$ Hermitian positive semidefinite, and $A_R+jA_I+S=0$; a single violation, or a dual bound that exceeds the best-known AC-OPF objective, would refute the dual-feasibility guarantee.
Extended reading notes
Core claim
The paper's central claim is that a dual conic proxy can be built for the semidefinite relaxation of AC-OPF, not just for the weaker second-order-cone relaxation used before. The architecture predicts a subset of the dual variables, recovers the rest in closed form, and handles the positive-semidefinite constraint by computing an intermediate Hermitian matrix, shifting it by the negative part of its smallest eigenvalue, and absorbing that shift into the voltage-bound dual variables. The authors assert that this completion always produces dual-feasible solutions, so each forward pass yields a valid lower bound on the AC-OPF optimum. On the reported test systems the SDP proxy improves the dual gap over the SOC proxy on the 30-, 118-, and 500-bus cases, and on the 118-bus case it gives a bound better than the exact SOC relaxation bound.
Load-bearing premise
The load-bearing premise is that the closed-form recovery rules imported from the earlier SOC-relaxation proxy, which this paper does not re-prove, remain valid for the SDP dual's structure, so that the completion always returns a dual-feasible point.
Editorial extensions
If this is right
- Every trained proxy output is a verified lower bound on the AC-OPF optimum, with no SDP solver run and no post-solve feasibility check at inference time.
- Operators can use the much tighter SDP relaxation at the same latency as a learned heuristic, which matters for repeated real-time dispatches.
- The self-supervised training loop removes the need to pre-compute optimal primal or dual solutions, so datasets cost only load perturbations, not solver calls.
- On systems where the SDP relaxation is far tighter than SOC, the proxy can reach small dual gaps the SOC proxy cannot approach, and on one benchmark it beats the exact SOC bound.
Reading between the lines
- If the dual-feasibility guarantee holds under a rigorous floating-point check, this pattern could make SDP-based certificates practical in market-clearing contexts where a valid bound must accompany every bid or dispatch decision.
- The size of $\delta = \min(0,\lambda_{\min}(\hat S))$ is a natural diagnostic for where the proxy is struggling: it measures how much of the semidefinite constraint is being satisfied by the shift rather than by the learned variables, and one could track it during training.
- On large grids, the global eigenvalue shift will not scale; a chordal-sparsity-aware version would need to complete each clique's PSD block separately, which would change the shift from a scalar to block-local corrections.
- Paired with a learned primal-feasible proxy, this dual proxy would yield an end-to-end learned optimality gap, letting a control system know how far its fast solution is from the true optimum without calling a solver.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a dual conic proxy (DCP) architecture for the semidefinite programming relaxation of AC optimal power flow (SDP-OPF). A neural network predicts a subset of dual variables, and a differentiable completion layer enforces all dual feasibility constraints of the conic dual DSDP-OPF, so that the network output is a certified lower bound on the optimal AC-OPF cost. The model is trained self-supervised by maximizing the dual objective, avoiding generation of ground-truth solutions. Numerical experiments on PGLib systems up to 500 buses compare the SDP proxy with a prior SOC-based proxy and with Mosek, reporting dual optimality gaps and gap-closed metrics. The paper claims several orders of magnitude speedups over interior-point SDP solvers, although timing results are not reported.
Significance. If the dual-feasibility guarantee is established, the paper makes a valuable contribution: it is the first DCP architecture for an SDP relaxation of AC-OPF, and the self-supervised training scheme avoids expensive data generation. The central mechanism — the eigenvalue-shift completion in step 5 of Algorithm 1 — is simple, elegant, and correctly enforces the PSD constraint (6h)-(6k) given a Hermitian completion. The experimental comparison against the SOC proxy is informative and shows that the stronger SDP relaxation can yield better bounds when the proxy is trained successfully. However, the main guarantee rests on recovery rules that are not specified or proved in this manuscript, and the headline speedup claim is not backed by any reported timing data. The paper is therefore promising but needs substantial revision before the central claims can be accepted.
major comments (3)
- [Algorithm 1, steps 3–4 and Section 3] The paper's central claim — that the DCP architecture 'always outputs dual feasible solutions' — depends on the closed-form recovery rules in steps 3 and 4 of Algorithm 1. These steps are not specified in the manuscript; the text only says they 'use the same argument as [6]', where [6] addresses a second-order cone relaxation, not the SDP dual. The manuscript does not provide the formulas for recovering ⃗ν^s and ⃗ν^s from (6j), nor the formulas for recovering ¯µ^pg, ¯µ^pg, ¯µ^qg, ¯µ^qg from (6b)-(6c) together with the nonnegativity constraints (6i). Because (6b)-(6c) are equalities with nonnegative multipliers, there is a family of feasible choices, and the paper does not say which closed-form selection is used or why it always satisfies the sign constraints and the SOC cone constraints. The reader cannot verify the dual-feasibility guarantee from this manuscript alone. The authors should include the explicit recovery formulas and a proof (or a precise reduction to a lemma in [6] that maps every constraint to the corresponding SDP dual constraint) for both forward and reverse cones and for the generation-bound multipliers.
- [Abstract, Section 4.3, and Conclusion] The abstract and conclusion state that the proposed SDP proxy achieves 'orders of magnitude of speedups' compared to interior-point solvers, but Section 4.3 says only that 'timing results, not reported for lack of space' support this claim. No timing measurements, hardware-specific runtime comparisons, or methodology for measuring inference time appear anywhere in the paper. Since the computational speedup is a primary motivation for the DCP approach and for the paper's title and abstract, this claim is unsupported as written. The authors should either include the timing results (e.g., median inference time per instance, including the eigenvalue decomposition and completion layer, versus Mosek solve time) or remove the speedup claim from the abstract and conclusion.
- [Section 2.2, Model 3, constraint (6h)] The completion step 5 uses the matrix ˆS = -A_R(λ,μ,ν) - j A_I(λ,μ,ν) and sets S = ˆS - δ I, with δ = min(0, λ_min(ˆS)) ≤ 0, and then sets ¯µ^w = 0, ¯µ^w = -δ. This correctly enforces S ≽ 0 and constraint (6h) only if the definitions of A_R and A_I are consistent with the signs in (6h). However, Equation (4) contains a notational ambiguity: the term for the voltage-bound multipliers is printed as '¯µ^w_i − ¯µ^w_i', which appears to use the same symbol for the upper and lower bounds; the intended expression is presumably (¯µ^w_i − µ̲^w_i). This ambiguity should be fixed, because the sign of δ and the assignment of ¯µ^w and ¯µ^w depend on which multiplier is which. If the symbols are indeed meant to be identical, the expression is wrong and the completion step would not satisfy (6h).
minor comments (5)
- [Section 3, after Algorithm 1] The sentence 'For a more detailed description of the DCP architecture, the reader may refer to [6, 8]' delegates a substantial part of the method description to previous work; the core novelty of this paper is precisely the SDP completion layer, so the main text should be self-contained for Algorithm 1.
- [Section 4.2, Eq. (8)] The 'dual optimality gap' is defined relative to z*_AC, described as the best-known AC-OPF solution obtained by Ipopt. Since Ipopt may return a local optimum, the terminology 'gap' should be clarified to indicate that it is computed against a reference solution, not a proven global optimum; this does not affect the validity of the dual bound itself.
- [Table 1] The table reports both mean and standard deviation, but for the 'gap closed' metric the arithmetic mean is used while the dual gaps use the geometric mean. This is stated, but the reader should be reminded that negative values of gap closed make geometric means undefined; the reporting choices are reasonable but worth a brief justification.
- [Section 4.3] The discussion of the ieee30 system attributes the weaker performance to 'the higher output dimension of the SDP proxy' and 'the additional difficulty of the PSD constraint (6k)', but no diagnostic experiments (e.g., training curves, loss values, or ablations) are provided to support this explanation. This is a minor point, but a short analysis would strengthen the interpretation.
- [Notation] The symbols ⃗ν and ⃗ν are used for forward and reverse dual SOC variables, and the same arrow convention is used for other flow variables; the paper should explicitly state the convention once in the notation section, since the distinction is critical for the cone constraints (6j).
Circularity Check
No circularity: the dual completion is an algebraic projection onto the dual feasible set and training maximizes the dual objective without fitting solver outputs.
full rationale
The derivation chain is not circular. Algorithm 1's completion layer is a sequence of algebraic projections onto the DSDP-OPF dual constraints: Step 2 substitutes predicted values into equalities (6d)-(6g); Step 3 chooses the minimum nu^s satisfying the second-order cone (6j), i.e. nu^s = ||(nu^p,nu^q)||_2; Step 4 chooses nonnegative generation multipliers as the positive and negative parts satisfying (6b)-(6c); Step 5 shifts hat S by -delta I and sets underline_mu^w = -delta to enforce (6h) and (6k). No step fits a parameter to the quantity it later predicts. Training is self-supervised with the dual objective (6a) as the loss, so maximizing it over the feasible dual set is legitimate and does not require solver targets; Mosek and Ipopt values are used only in test-time evaluation. The reference to the authors' prior work [6] for Steps 3-4 is a dependency on an earlier derivation, not a self-citation that carries the conclusion: the closed-form rules for cone recovery and multiplier splitting are elementary and are expressed in the present paper's own constraint set. Even if the correctness of importing SOC recovery rules to the SDP dual were questioned, that would be a proof-completeness or correctness risk, not a circularity in the derivation.
Assumptions & free parameters
free parameters (1)
- Neural network parameters (weights and biases)
assumptions (3)
- standard math Weak duality between SDP-OPF (Model 2) and its conic dual DSDP-OPF (Model 3)
- domain assumption SDP-OPF is a valid relaxation of AC-OPF via W = V V* and rank relaxation
- ad hoc to paper The dual recovery rules of steps 3-4 from [6] produce feasible dual variables for the SDP dual
Cite this review
Pith. "Pith review of Dual Conic Proxy for Semidefinite Relaxation of AC Optimal Power Flow." pith.science (2026). https://pith.science/paper/CY2HSUTV
@misc{pith2026250206978,
author = {Pith},
title = {Pith review of: Dual Conic Proxy for Semidefinite Relaxation of AC Optimal Power Flow},
year = {2026},
howpublished = {\url{https://pith.science/paper/CY2HSUTV}},
note = {Machine review of arXiv:2502.06978}
}
read the original abstract
The nonlinear, non-convex AC Optimal Power Flow (AC-OPF) problem is fundamental for power systems operations. The intrinsic complexity of AC-OPF has fueled a growing interest in the development of optimization proxies for the problem, i.e., machine learning models that predict high-quality, close-to-optimal solutions. More recently, dual conic proxy architectures have been proposed, which combine machine learning and convex relaxations of AC-OPF, to provide valid certificates of optimality using learning-based methods. Building on this methodology, this paper proposes, for the first time, a dual conic proxy architecture for the semidefinite (SDP) relaxation of AC-OPF problems. Although the SDP relaxation is stronger than the second-order cone relaxation considered in previous work, its practical use has been hindered by its computational cost. The proposed method combines a neural network with a differentiable dual completion strategy that leverages the structure of the dual SDP problem. This approach guarantees dual feasibility, and therefore valid dual bounds, while providing orders of magnitude of speedups compared to interior-point algorithms. The paper also leverages self-supervised learning, which alleviates the need for time-consuming data generation and allows to train the proposed models efficiently. Numerical experiments are presented on several power grid benchmarks with up to 500 buses. The results demonstrate that the proposed SDP-based proxies can outperform weaker conic relaxations, while providing several orders of magnitude speedups compared to a state-of-the-art interior-point SDP solver.
Reference graph
Works this paper leans on
-
[6]
Dual conic proxies for ac optimal power flow
Guancheng Qiu, Mathieu Tanneau, and Pascal V an Hentenry ck. Dual conic proxies for ac optimal power flow. Electric Power Systems Research, 236:110661, 2024
work page 2024
-
[1]
Mak, and Pascal V an H entenryck
Ferdinando Fioretto, Terrence W .K. Mak, and Pascal V an H entenryck. Predicting AC Optimal Power Flows: Combining Deep Learning and Lagrangian Dual Methods. Proceedings of the AAAI Conference on Artificial Intelligence, 34(01):630–637, Apr. 2020. doi: 10.1609/aaai.v34i01.54 03. REFERENCES 6
-
[2]
DC3: A lea rning method for optimization with hard constraints
Priya L Donti, David Rolnick, and J Zico Kolter. DC3: A lea rning method for optimization with hard constraints. arXiv preprint arXiv:2104.12225 , 2021
arXiv 2021
-
[3]
D eepOPF-V: Solving AC-OPF problems effi- ciently
Wanjun Huang, Xiang Pan, Minghua Chen, and Steven H Low. D eepOPF-V: Solving AC-OPF problems effi- ciently. IEEE Transactions on Power Systems, 37(1):800–803, 2021
work page 2021
-
[4]
End-to-End Feasible Optimization Proxies for Large-Scale Economic Dispatch
Wenbo Chen, Mathieu Tanneau, and Pascal V an Hentenryck. End-to-End Feasible Optimization Proxies for Large-Scale Economic Dispatch. IEEE Transactions on Power Systems, 39(2):4723–4734, 2024
work page 2024
-
[5]
Energy and operating reserve markets, 2022
MISO. Energy and operating reserve markets, 2022. Busin ess Practices Manual
work page 2022
-
[7]
R. A. Jabr. Radial distribution load flow using conic prog ramming. IEEE Transactions on Power Systems, 21(3): 1458–1459, Aug 2006. ISSN 0885-8950. doi: 10.1109/TPWRS.2 006.879234
doi:10.1109/tpwrs.2 2006
-
[8]
Dual Lagrang ian Learning for Conic Optimization
Mathieu Tanneau and Pascal V an Hentenryck. Dual Lagrang ian Learning for Conic Optimization. In Advances in Neural Information Processing Systems , volume 37, 2024
work page 2024
Show all 18 references
-
[9]
Molzahn and Ian A
Daniel K. Molzahn and Ian A. Hiskens. A Survey of Relaxati ons and Approximations of the Power Flow Equations. F oundations and Trends® in Electric Energy Systems, 4(1-2):1–221, 2019. ISSN 2332-6557
2019
-
[10]
Strong SOCP re laxations for the optimal power flow problem
Burak Kocuk, Santanu S Dey, and X Andy Sun. Strong SOCP re laxations for the optimal power flow problem. Operations Research, 64(6):1177–1196, 2016
2016
-
[11]
Semidefinite programming for optimal power flow problems
Xiaoqing Bai, Hua Wei, Katsuki Fujisawa, and Y ong Wang. Semidefinite programming for optimal power flow problems. International Journal of Electrical Power & Energy Systems, 30(6):383–392, 2008. ISSN 0142-0615. doi: 10.1016/j.ijepes.2007.12.003
2008 doi
-
[12]
The power grid lib rary for benchmarking ac optimal power flow algorithms
Sogol Babaeinejadsarookolaee et al. The power grid lib rary for benchmarking ac optimal power flow algorithms. arXiv preprint arXiv:1908.02788 , 2019
1908 arXiv
-
[13]
Powermodels.jl: An open-source framework for exploring power flow formulations
Carleton Coffrin, Russell Bent, Kaarthik Sundar, Y eesian Ng, and Miles Lubin. Powermodels.jl: An open-source framework for exploring power flow formulations. In 2018 Power Systems Computation Conference (PSCC) , pages 1–8, June 2018. doi: 10.23919/PSCC.2018.8442948
2018
-
[14]
On the implementa tion of an interior-point filter line-search algorithm for large-scale nonlinear programming
Andreas Wächter and Lorenz T Biegler. On the implementa tion of an interior-point filter line-search algorithm for large-scale nonlinear programming. Mathematical programming, 106:25–57, 2006
2006
-
[15]
MOSEK Optimizer API for Julia 10.1.12 , 2022
MOSEK ApS. MOSEK Optimizer API for Julia 10.1.12 , 2022. URL https://docs.mosek.com/10.1/juliaapi/index.html
2022
-
[16]
PGLearn - An Open-Source Learning Toolkit for Optimal P ower Flow. 2024. URL https://openreview.net/forum?id=cecIf0CKnH
2024
-
[17]
Pytorch: An imperative style, high-p erformance deep learning library
Adam Paszke et al. Pytorch: An imperative style, high-p erformance deep learning library. Advances in neural information processing systems, 32, 2019
2019
-
[18]
Adam: A method for stocha stic optimization
Diederik P Kingma and Jimmy Ba. Adam: A method for stocha stic optimization. arXiv preprint arXiv:1412.6980, 2014
2014 arXiv
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.