Pith. sign in

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 →

arxiv 1908.02319 v1 pith:U544HM32 submitted 2019-08-06 math.OC

classification math.OC MSC 90C2290C26
keywords ACoptimalpowerflowconicrelaxationsemidefiniteprogrammingsecond-orderconereformulation-linearizationtechniquechordalconvexsystems
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper proposes a new convex relaxation, the tight-and-cheap relaxation (TCR), for the alternating-current optimal power flow problem, a nonconvex and NP-hard problem central to power-system operations. TCR is built by combining semidefinite constraints with reformulation-linearization (RLT) inequalities on voltage magnitudes, and it is shown to be stronger than the standard second-order cone relaxation and nearly as tight as the full semidefinite relaxation. The authors prove that a closely related formulation is exactly equivalent to the standard SDP relaxation, so TCR inherits a rigorous lower-bounding property. On standard test networks with up to 6515 buses, TCR's optimality gap is close to SDR's, and its solve time is about 30 times lower than the chordal relaxation on large instances and 55 times lower on the largest ones. A strengthened variant, STCR, is tighter still and coincides with SDR when the network minus the reference bus is acyclic.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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.
  3. [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

0 steps flagged · score 0.0 of 10

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 0 free parameters · 4 assumptions · 0 invented entities

No fitted parameters; all constraints are derived from first principles. Standard convex analysis and chordal graph theory are used as axioms. The only domain assumptions are the standard ACOPF model assumptions: positive voltage bounds, convex generation costs, and the reference bus angle fixed.

assumptions (4)
  • standard math Positive semidefinite matrices and Schur complement properties
    Used throughout the proofs of Propositions 2 and 3.
  • standard math Chordal graph theory: PSD of partial matrices equivalent to clique PSD on chordal extension
    Invoked in Section III-B and in the proof of Proposition 4.
  • domain assumption Voltage magnitude bounds are strictly positive (vk > 0)
    Stated in Section II after constraint (1h); needed for the construction of v in Proposition 2.
  • domain assumption Generation cost functions are convex
    Assumed in Section II for the ACOPF model.

how reviews work

0 comments
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.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

38 extracted references · 37 canonical work pages

  1. [25]

    Molzahn, Jesse T

    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

  2. [27]

    Madani, M

    R. Madani, M. Ashraphijuo, and J. Lavaei. OPF solver. https://lavaei.ieor.berkeley.edu/Software.html, 2014

  3. [1]

    Carpentier

    J. Carpentier. Contribution to the economic dispatch pr oblem. Bulletin de la Soci´ et´ e franc ¸aise des ´ electriciens, 3(8):431–447, 1962

  4. [2]

    Cain, Richard P

    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

  5. [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

  6. [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

  7. [5]

    Power Grid Security Analysis: An Optimization Approach

    Abhinav V erma. Power Grid Security Analysis: An Optimization Approach. Columbia University, 2010

  8. [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

Show all 38 references
  1. [7]

    Convex Optimization

    Stephen Boyd and Lieven V andenberghe. Convex Optimization. Cam- bridge University Press, 2004

  2. [8]

    Semidefinite Programming for Combinatorial Optimization

    Christoph Helmberg. Semidefinite Programming for Combinatorial Optimization. Konrad-Zuse-Zentrum f¨ ur Informationstechnik Berlin, 2000

  3. [9]

    R. A. Jabr. Radial distribution load flow using conic prog ramming. IEEE Transactions on Power Systems , 21(3):1458–1459, Aug 2006

  4. [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

  5. [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

  6. [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

  7. [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

  8. [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

  9. [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

  10. [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

  11. [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

  12. [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

  13. [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

  14. [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

  15. [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

  16. [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

  17. [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...

  18. [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

  19. [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

  20. [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

  21. [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

  22. [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

  23. [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

  24. [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

  25. [33]

    C´ edric Josz and Daniel K. Molzahn. Moment/sum-of-squ ares hierarchy for complex polynomial optimization. arXiv preprint arXiv:1508.02068, 2015

  26. [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

  27. [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

  28. [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

  29. [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...

  30. [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...

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.