REVIEW 2 major objections 3 minor 38 references
Tight-and-cheap conic relaxation for the AC optimal power flow problem
T0 review · 2 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper introduces a conic relaxation for AC optimal power flow that is stronger than the second-order cone relaxation, nearly as tight as the semidefinite relaxation, and much faster to solve than the chordal relaxation on large…
desk verdict Solid new conic relaxation for ACOPF with correct proofs, but the abstract overstates tightness and the speed claims rest on a generic CVX baseline for chordal SDPs. 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 per-branch $3\times3$ semidefinite block $$ \begin{bmatrix} 1 & v_k^* & v_m^* \\ v_k & V_{kk} & V_{km} \\ v_m & V_{km}^* & V_{mm} \end{bmatrix} \succeq 0, $$ imposed for every branch $(k,m)$. This block simultaneously encodes the second-order cone relation $|V_{km}|^2 \le V_{kk}V_{mm}$ and couples the complex voltage variables $v$ to $V$. The RLT ingredient is a secant lower bound on the concave square root: Lemma 1 gives $\sqrt{y} \ge (y + \sqrt{\ell u})/(\sqrt{\ell}+\sqrt{u})$ on $\ell \le y \le u$, applied to $V_{11}=|v_1|^2$ to yield the reference-bus inequality $\mathrm{Re}(v_1) \ge (V_{11}+\underline{v}_1\overline{v}_1)/(\underline{v}_1+\overline{v}_1)$. That inequality is what lets the proof recover $v$ from $V$, making the equivalence between nSDR and SDR go through.
What would settle it
Run a timing comparison on the 6515-bus test network using a semidefinite solver that exploits sparsity, including the time to construct the chordal extension; if the sparsity-exploiting relaxation then solves in time comparable to or less than TCR while giving the same bound, the paper's central speed claim is refuted.
Extended reading notes
Core claim
The central discovery is that the full positive semidefinite constraint $V \succeq 0$ in the standard SDP relaxation can be replaced, without losing tightness, by the combination of $V \succeq vv^H$ (where $v$ is the complex voltage vector) and a single RLT inequality at the reference bus. Proposition 2 constructs, for any feasible $V$, the vector $v = Ve_1/\sqrt{V_{11}}$ and uses Schur-complement reasoning to show this pair satisfies the new constraints, so the new semidefinite relaxation nSDR is exactly equivalent to SDR. Replacing the global constraint $V \succeq vv^H$ by per-branch $3\times3$ semidefinite blocks on $(1, v_k, v_m)$ yields TCR, which dominates the second-order cone relaxation and is dominated by SDR. The strengthened STCR, which uses the $3\times3$ principal submatrices on the reference bus together with each branch, is strictly stronger than TCR and is proven equivalent to SDR when the network graph with the reference bus removed has no cycles. On unmodified test cases up to 6515 buses, TCR's optimality gaps are close to SDR's while its solve times are about 30 times lower than the chordal relaxation on large instances.
Load-bearing premise
The speed advantage rests on comparing TCR against a sparsity-exploiting semidefinite relaxation that is solved with a generic solver and timed without counting the cost of building the sparsity structure; a purpose-built solver could shrink that advantage.
Editorial extensions
If this is right
- On large networks, TCR gives a lower bound within a few hundredths of a percent of the SDP bound while solving roughly 30 times faster than the chordal relaxation, making SDP-quality bounds practical at scales where the full SDP was too slow.
- Because TCR uses only per-branch $3\times3$ semidefinite blocks, it avoids storing or factoring an $n\times n$ positive semidefinite matrix, which is the main source of SDP's computational cost.
- STCR is strictly tighter than TCR and is equivalent to SDR for networks whose graph minus the reference bus has no cycles, so SDP-quality bounds can be obtained cheaply on such topologies.
- As a valid relaxation, TCR provides a certificate (lower bound) on the global optimum of ACOPF, so it can be used to bound the suboptimality of any local solution found by nonlinear solvers.
Reading between the lines
- The same construction could be applied at other buses: adding RLT cuts on more voltage magnitudes, or on branch current magnitudes, would likely tighten TCR further while preserving its local structure; the paper only uses the reference bus.
- The equivalence strategy of Proposition 2 — global PSD replaced by $V \succeq vv^H$ plus one scalar RLT cut at an anchor variable — may transfer to other complex quadratically constrained programs with bounded magnitudes, not just power flow.
- The speed comparison uses a generic solver for the sparsity-exploiting relaxation and does not count the time to build its chordal structure; a purpose-built solver could narrow TCR's advantage, so the 'cheap' part should be revalidated in production-grade software.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a new convex relaxation for the AC optimal power flow problem, obtained by combining semidefinite constraints with reformulation-linearization inequalities on voltage magnitudes. The authors first introduce nSDR, a reformulation of the standard SDP relaxation with additional voltage variables and an RLT constraint at the reference bus, and prove in Proposition 2 that nSDR is equivalent to SDR. They then drop the global positive-semidefinite constraint to local 3-by-3 blocks, obtaining the tight-and-cheap relaxation (TCR), and a further strengthened variant (STCR) using 3-by-3 principal submatrices anchored at the reference bus. Propositions 3 and 4 establish that STCR is stronger than TCR and, for networks whose graph becomes acyclic after removing the reference bus, equivalent to SDR. Computational experiments on MATPOWER instances up to 6515 buses compare TCR and STCR with SOCR, CHR, and SDR in terms of optimality gap and solver time. The paper claims that TCR is stronger than SOCR, nearly as tight as SDR, and up to an order of magnitude faster than the chordal relaxation.
Significance. The theoretical hierarchy in the paper is clean and the proofs of Propositions 2-4 are self-contained and correct; the nSDR/TCR construction is a legitimate and potentially useful way to localize the semidefinite constraint in OPF relaxations. The empirical study is also carefully done in one respect: the authors use unmodified MATPOWER cases and report both cost and loss objectives, with no fitted parameters. However, two load-bearing applied claims are not fully supported by the evidence as presented. First, the abstract's 'nearly as tight as SDR' is contradicted by the paper's own tables on instances where the SDR bound is not already near-exact. Second, the speed comparison against the chordal relaxation is made with a generic CVX/MOSEK interface that does not exploit chordal structure, so the reported speedups may not represent a state-of-the-art chordal SDP solver. If these computational claims are hardened or appropriately qualified, the paper would be a solid contribution to the OPF relaxation literature.
major comments (2)
- [Abstract and Section V, Tables I and II] The claim that TCR is 'nearly as tight as the standard semidefinite relaxation' is not supported on instances where SDR itself has a non-negligible gap. In Table I, case5 has a TCR optimality gap of 12.75% versus 5.22% for both SDR and CHR, and case_ACTIV_SG_500 has a TCR gap of 4.39% versus 2.11% for SDR/CHR. The paper acknowledges these discrepancies in Section V, bullet 3, but the abstract and conclusion state the tightness claim without qualification. The authors should either restrict the closeness claim to instances where the SDR bound is already tight, or add a quantitative statement about the worst-case or instance-dependent gap between TCR and SDR.
- [Section V, first paragraphs and Tables I and II] The central speed claim, that TCR is up to an order of magnitude faster than the chordal relaxation, is measured by solving every relaxation with CVX 2.1 and MOSEK 8.0.0.60 under default settings. CVX passes each positive-semidefinite block to the solver as a generic constraint and does not exploit the clique tree that makes chordal SDP relaxations scalable, while the paper cites specialized chordal implementations in [25]-[27]. The CHR times in Tables I and II are therefore likely not representative of a properly implemented chordal relaxation, and the abstract's speedup claim needs to be either benchmarked against a chordal-exploiting solver or explicitly limited to generic SDP solvers.
minor comments (3)
- [Model 2] The initialization line contains a doubled word: 'LG + I|N| ≻ 0, where LG is the Laplacian matrix of G and and I|N| is the identity matrix of size |N|.' Please fix the typo.
- [Section V, paragraph before Tables I and II] The sentence 'Both objective functions of test cases from [36] are the same' is ambiguous. Please clarify that for the RTE/PEGASE cases the generation-cost coefficients already encode active loss, which explains why the large-scale rows of Tables I and II are identical.
- [Section IV, Model 5] It would help the reader to state explicitly that constraint (8) implies the usual 2-by-2 SOCR constraints by taking the leading 2-by-2 principal submatrix, rather than leaving the domination argument entirely to the phrase 'Clearly TCR dominates SOCR.'
Circularity Check
No circularity found: the relaxations are derived and proved from stated assumptions, and the computational claims are benchmarked on external unmodified MATPOWER cases.
full rationale
The paper's derivation chain is self-contained. The new relaxations are built from the SDP relaxation by explicit RLT inequalities and a positive-semidefinite block constraint; Proposition 2 proves equivalence of nSDR and SDR using a constructive argument (Lemma 1 and Lemma 2), and Propositions 3 and 4 prove the dominance and exactness results from the given constraints and chordal graph theory. No parameter is fitted to the test data and no 'prediction' is defined in terms of the quantity it is said to predict. The empirical speed comparisons use unmodified MATPOWER test cases and report solver times from MOSEK, which is an external benchmark. The caveat that CHR is solved via CVX's generic interface rather than a specialized chordal-exploiting implementation affects the fairness of the speed comparison but is not a circularity in the derivation. There are no load-bearing self-citations: the cited completion theorem of Grone et al. is an external mathematical result. Therefore the circularity score is 0.
Assumptions & free parameters
assumptions (4)
- standard math Positive semidefinite matrices and Schur complement properties
- standard math Chordal graph theory: PSD of partial matrices equivalent to clique PSD on chordal extension
- domain assumption Voltage magnitude bounds are strictly positive (vk > 0)
- domain assumption Generation cost functions are convex
Cite this review
Pith. "Pith review of Tight-and-cheap conic relaxation for the AC optimal power flow problem." pith.science (2026). https://pith.science/paper/U544HM32
@misc{pith2026190802319,
author = {Pith},
title = {Pith review of: Tight-and-cheap conic relaxation for the AC optimal power flow problem},
year = {2026},
howpublished = {\url{https://pith.science/paper/U544HM32}},
note = {Machine review of arXiv:1908.02319}
}
read the original abstract
The classical alternating current optimal power flow problem is highly nonconvex and generally hard to solve. Convex relaxations, in particular semidefinite, second-order cone, convex quadratic, and linear relaxations, have recently attracted significant interest. The semidefinite relaxation is the strongest among them and is exact for many cases. However, the computational efficiency for solving large-scale semidefinite optimization is lower than for second-order cone optimization. We propose a conic relaxation obtained by combining semidefinite optimization with the reformulation-linearization technique, commonly known as RLT. The proposed relaxation is stronger than the second-order cone relaxation and nearly as tight as the standard semidefinite relaxation. Computational experiments using standard test cases with up to 6515 buses show that the time to solve the new conic relaxation is up to one order of magnitude lower than for the chordal relaxation, a semidefinite relaxation technique that exploits the sparsity of power networks.
Reference graph
Works this paper leans on
-
[25]
Daniel K. Molzahn, Jesse T. Holzer, Bernard C. Lesieutr e, and Christo- pher L. DeMarco. Implementation of a large-scale optimal po wer flow solver based on semidefinite programming. IEEE Transactions on Power Systems, 28(4):3987–3998, 2013
work page 2013
- [27]
-
[1]
J. Carpentier. Contribution to the economic dispatch pr oblem. Bulletin de la Soci´ et´ e franc ¸aise des ´ electriciens, 3(8):431–447, 1962
work page 1962
-
[2]
Mary B. Cain, Richard P . ONeill, and Anya Castillo. Histo ry of optimal power flow and formulations. Federal Energy Regulatory Commission , pages 1–36, 2012
work page 2012
-
[3]
Optimal power flow: A bibliographic survey I
Stephen Frank, Ingrida Steponavice, and Steffen Rebenn ack. Optimal power flow: A bibliographic survey I. Energy Systems , 3(3):221–258, 2012
work page 2012
-
[4]
Optimal power flow: A bibliographic survey II
Stephen Frank, Ingrida Steponavice, and Steffen Rebenn ack. Optimal power flow: A bibliographic survey II. Energy Systems , 3(3):259–289, 2012
work page 2012
-
[5]
Power Grid Security Analysis: An Optimization Approach
Abhinav V erma. Power Grid Security Analysis: An Optimization Approach. Columbia University, 2010
work page 2010
-
[6]
AC- feasibility on tree networks is NP-hard
Karsten Lehmann, Alban Grastien, and Pascal V an Hentenr yck. AC- feasibility on tree networks is NP-hard. IEEE Transactions on Power Systems, 31(1):798–801, 2016
work page 2016
Show all 38 references
-
[7]
Convex Optimization
Stephen Boyd and Lieven V andenberghe. Convex Optimization. Cam- bridge University Press, 2004
2004
-
[8]
Semidefinite Programming for Combinatorial Optimization
Christoph Helmberg. Semidefinite Programming for Combinatorial Optimization. Konrad-Zuse-Zentrum f¨ ur Informationstechnik Berlin, 2000
2000
-
[9]
R. A. Jabr. Radial distribution load flow using conic prog ramming. IEEE Transactions on Power Systems , 21(3):1458–1459, Aug 2006
2006
-
[10]
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
2008
-
[11]
R. A. Jabr. Exploiting sparsity in SDP relaxations of th e OPF problem. IEEE Transactions on Power Systems , 27(2):1138–1139, May 2012
2012
-
[12]
S. H. Low. Convex relaxation of optimal power flow–Part I : Formu- lations and equivalence. IEEE Transactions on Control of Network Systems, 1(1):15–27, March 2014
2014
-
[13]
S. H. Low. Convex relaxation of optimal power flow–Part I I: Exactness. IEEE Transactions on Control of Network Systems , 1(2):177–189, June 2014
2014
-
[14]
Convex quadratic relaxations for mixed-integer nonlinear progra ms in power systems
Hassan Hijazi, Carleton Coffrin, and Pascal V an Henten ryck. Convex quadratic relaxations for mixed-integer nonlinear progra ms in power systems. Mathematical Programming Computation , pages 1–47, 2014
2014
-
[15]
Polynomial SDP cuts for optimal power flow
Hassan Hijazi, Carleton Coffrin, and Pascal V an Henten ryck. Polynomial SDP cuts for optimal power flow. In Power Systems Computation Conference (PSCC), 2016 , pages 1–7. IEEE, 2016
2016
-
[16]
Application of the moment-SOS approach to global optimizat ion of the OPF problem
C´ edric Josz, Jean Maeght, Patrick Panciatici, and Jea n Charles Gilbert. Application of the moment-SOS approach to global optimizat ion of the OPF problem. IEEE Transactions on Power Systems , 30(1):463–470, 2015
2015
-
[17]
Network flow and copper plate relaxations for AC transmission system s
Carleton Coffrin, Hassan Hijazi, and Pascal V an Henten ryck. Network flow and copper plate relaxations for AC transmission system s. In Power Systems Computation Conference (PSCC), 2016 , pages 1–8. IEEE, 2016
2016
-
[18]
Sleiman Mhanna, Gregor V erbiˇ c, and Archie C. Chapman. Tight LP approximations for the optimal power flow problem. In Power Systems Computation Conference (PSCC), 2016 , pages 1–7. IEEE, 2016
2016
-
[19]
Dey, and X
Burak Kocuk, Santanu S. Dey, and X. Andy Sun. Strong SOCP relaxations for the optimal power flow problem. Operations Research , 64(6):1177–1196, 2016
2016
-
[20]
R. D. Zimmerman, C. E. Murillo-Sanchez, and R. J. Thomas . MA T- POWER: Steady-state operations, planning, and analysis to ols for power systems research and education. IEEE Transactions on Power Systems , 26(1):12–19, Feb 2011
2011
-
[21]
Javad Lavaei and Steven H. Low. Zero duality gap in optim al power flow problem. IEEE Transactions on Power Systems , 27(1):92–107, 2012
2012
-
[22]
Johnson, Eduardo M
Robert Grone, Charles R. Johnson, Eduardo M. S´ a, and He nry Wolkow- icz. Positive definite completions of partial Hermitian mat rices. Linear algebra and its applications , 58:109–124, 1984
1984
-
[23]
Exploiting sparsity in semidefinite programming vi a matrix completion I: General framework
Mituhiro Fukuda, Masakazu Kojima, Kazuo Murota, and Ka zuhide Nakata. Exploiting sparsity in semidefinite programming vi a matrix completion I: General framework. SIAM Journal on Optimization , 11(3):647–674, 2001. 7 TABLE I: Cost minimization Test case υ [$/h] υ [$/h] Optimali...
2001
-
[24]
Exploiting sparsity in semidefinite program- ming via matrix completion II: Implementation and numerica l results
Kazuhide Nakata, Katsuki Fujisawa, Mituhiro Fukuda, M asakazu Ko- jima, and Kazuo Murota. Exploiting sparsity in semidefinite program- ming via matrix completion II: Implementation and numerica l results. Mathematical Programming, 95(2):303–327, 2003
2003
-
[26]
Andersen, Anders Hansson, and Lieven V andenb erghe
Martin S. Andersen, Anders Hansson, and Lieven V andenb erghe. Reduced-complexity semidefinite relaxations of optimal po wer flow problems. IEEE Transactions on Power Systems , 29(4):1855–1863, 2014
2014
-
[28]
R. A. Jabr. A conic quadratic format for the load flow equa tions of meshed networks. IEEE Transactions on Power Systems , 22(4):2285– 2286, Nov 2007
2007
-
[29]
McCormick
Garth P . McCormick. Computability of global solutions to factorable nonconvex programs–Part I: Convex underestimating proble ms. Mathe- matical Programming, 10(1):147–175, 1976
1976
-
[30]
Anstreicher
Kurt M. Anstreicher. Semidefinite programming versus t he reformulation-linearization technique for nonconvex qua dratically con- strained quadratic programming. Journal of Global Optimization , 43(2- 3):471–484, 2009
2009
-
[31]
Anstreicher and Samuel Burer
Kurt M. Anstreicher and Samuel Burer. Computable repre sentations for convex hulls of low-dimensional quadratic forms. Mathematical programming, 124(1):33–43, 2010
2010
-
[32]
Chen Chen, Alper Atamt¨ urk, and Shmuel S. Oren. A spatia l branch- and-cut method for nonconvex QCQP with bounded complex vari ables. Mathematical Programming, pages 1–29, 2017
2017
-
[33]
C´ edric Josz and Daniel K. Molzahn. Moment/sum-of-squ ares hierarchy for complex polynomial optimization. arXiv preprint arXiv:1508.02068, 2015
2015 arXiv
-
[34]
Hijazi, and Pascal V an Hent enryck
Carleton Coffrin, Hassan L. Hijazi, and Pascal V an Hent enryck. Strengthening the SDP relaxation of AC power flows with conve x envelopes, bound tightening, and valid inequalities. IEEE Transactions on Power Systems , 32(5):3549–3558, 2017
2017
-
[35]
A. B. Birchfield, T. Xu, K. M. Gegner, K. S. Shetye, and T. J . Over- bye. Grid structural characteristics as validation criter ia for synthetic networks. IEEE Transactions on Power Systems , 32(4):3258–3265, July 2017
2017
-
[36]
AC power flow data in MA TPOWER and QCQP format: iTesla, RTE snapshots, and PEGASE
C´ edric Josz, St´ ephane Fliscounakis, Jean Maeght, and Patrick Panciatici. AC power flow data in MA TPOWER and QCQP format: iTesla, RTE snapshots, and PEGASE. arXiv preprint arXiv:1603.01533 , 2016
2016 arXiv
-
[37]
CVX: MA TLAB software for discipline d convex programming, version 2.1
CVX Research, Inc. CVX: MA TLAB software for discipline d convex programming, version 2.1. http://cvxr.com/cvx, August 2012. 8 TABLE II: Loss minimization Test case υ [MW] υ [MW] Optimality gap [%] Computation time [s] SOCR TCR STCR CHR SDR SOCR TCR STCR CHR SDR Small-scale in...
2012
-
[38]
Grant and S
M. Grant and S. Boyd. Graph implementations for nonsmoo th convex programs. In V . Blondel, S. Boyd, and H. Kimura, editors, Recent Advances in Learning and Control , Lecture Notes in Control and Information Sciences, pages 95–110. Springer-V erlag Limi ted, 2008. http://stanf...
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.