44.05.01 · optimization-control / 05-conic-semidefinite

Conic Optimization: Linear, Second-Order Cone, and Semidefinite Programs

shipped3 tiersLean: none

Anchor (Master): Ben-Tal & Nemirovski 2001 Lectures on Modern Convex Optimization (SIAM) Lec. 1-4; Nesterov & Nemirovski 1994 Interior-Point Polynomial Algorithms in Convex Programming (SIAM) Ch. 1-5; Vandenberghe & Boyd 1996 Semidefinite Programming (SIAM Review) §1-§6

Intuition Beginner

Linear programming is the workhorse of optimization: minimize a linear cost while every variable stays nonnegative and a few linear equations hold. The nonnegativity rule says each coordinate lives on the good side of zero. Geometrically, the allowed region is a flat-faced shape — corners and straight edges — and the best point sits at a corner.

Conic optimization keeps the linear cost and the linear equations but swaps the "every coordinate is nonnegative" rule for a richer kind of rule: the variable must lie inside a chosen cone. A cone is a region that is closed under scaling up from the origin — if a point is in it, the whole ray from the origin through that point is in it too. The nonnegative quadrant is one cone, but there are curved ones as well.

Why bother with curved cones? Because two of them unlock problems that flat shapes cannot describe. The ice-cream cone (its formal name is the second-order or Lorentz cone) lets you say "the length of this vector is at most this number" — a smooth, round constraint that captures distances, norms, and uncertainty. The cone of positive-semidefinite matrices lets you constrain the eigenvalues of a matrix all at once, which controls vibration frequencies, stability, and energy.

The beautiful part is that all three live in one family. A linear program is a conic program over the quadrant; a second-order cone program is a conic program over the ice-cream cone; a semidefinite program is a conic program over the matrix cone. Each one contains the previous as a special case, so a single theory — one notion of duality, one kind of solver — handles all of them at once.

Visual Beginner

Figure: three cones side by side in increasing generality. On the left, the nonnegative quadrant of the plane — a flat wedge with straight edges, labelled LP. In the middle, the round ice-cream cone in three dimensions: all points whose horizontal radius is at most their height, labelled SOCP. On the right, a smooth curved cone of matrices labelled SDP, drawn as a rounded solid. Arrows of inclusion point left-to-right: the wedge sits inside the ice-cream cone, which sits inside the matrix cone.

      LP cone                 SOCP cone                  SDP cone
   (nonneg orthant)        (ice-cream / Lorentz)      (PSD matrices)

      t |  . . . .            t  |    /\               curved smooth
        | . region .             |   /  \              solid: matrices
        |. . . . . .             |  / .. \             with all
        +-----------  x          | /......\            eigenvalues >= 0
        |                        +---------- x
   straight edges,          round circular base,    a "bowl" of matrices;
   sharp corners            smooth side wall         its boundary is curved

          LP    is inside   SOCP    is inside    SDP
   (each is a special case of the next: one theory covers all three)

Worked example Beginner

We solve a small second-order cone problem by hand and watch the round constraint do its work. Minimize the cost (so we want as large as possible) subject to the single round rule that the point stays inside the unit disk: .

Step 1. Read the rule as a cone constraint. The rule "length at most " is the ice-cream-cone condition with the height fixed at . It is a smooth, curved constraint, not a flat one.

Step 2. Guess the best direction. We want large. Pushing equally in and uses the budget best, so try the point on the circle in the direction , normalized to length : that is .

Step 3. Check feasibility. Its length is , so the point sits exactly on the boundary of the disk and the rule holds.

Step 4. Compute the cost. Here , so the cost .

What this tells us. Unlike a linear program, the best point is not at a sharp corner — the disk has no corners. It is the point where the boundary circle is tangent to the cost direction. Curved cones move the optimum onto smooth boundary, and that single change is what lets conic programs capture distances and norms that linear programs cannot.

Check your understanding Beginner

Formal definition Intermediate+

Throughout, carries the standard inner product; in matrix space (real symmetric matrices) the inner product is .

Definition (proper cone and generalized inequality). A convex cone is proper if it is closed, has nonempty interior (solid), and contains no line (pointed: ). A proper cone induces a partial order, the generalized inequality

which is reflexive, antisymmetric, transitive, and compatible with addition and nonnegative scaling [Boyd & Vandenberghe §2.4]. Writing means .

Definition (conic program). Given , an affine map specified by and , and a proper cone , the (primal) conic program in standard form is

It is a convex program: the objective is linear and the feasible set is the intersection of an affine subspace with a convex cone, hence convex (the operation calculus of 44.01.01).

Definition (dual cone). The dual cone of is , itself a proper cone whenever is, with for closed convex . A cone is self-dual if .

The three canonical cones. Each is proper and self-dual under its natural inner product:

  • Nonnegative orthant . The conic program is a linear program (LP); is componentwise inequality.
  • Second-order (Lorentz / ice-cream) cone . A conic program with this cone (or a product of such cones) is a second-order cone program (SOCP).
  • Positive-semidefinite cone , i.e. all eigenvalues nonnegative (01.01.13). The conic program over is a semidefinite program (SDP); means is PSD. In standard SDP form is read as and the objective as .

The inclusion LP SOCP SDP. A componentwise inequality is a product of one-dimensional Lorentz constraints , so every LP is an SOCP. A Lorentz constraint is equivalent to a linear matrix inequality via the Schur complement (below): , so every SOCP is an SDP. The expressive ladder is strict: SDP captures eigenvalue constraints no SOCP can.

Definition (conic dual problem). The Lagrange dual of the primal conic program, pairing the equality with a free multiplier and the cone membership with a multiplier in the dual cone, is the conic dual

Here is the dual slack. The pairing is symmetric: dualizing the dual returns the primal, because .

Counterexamples to common slips

  • Not every closed convex cone is proper. A halfspace is a closed convex cone with nonempty interior but contains a line, so it is not pointed and induces no antisymmetric order. The orthant, Lorentz, and PSD cones are proper precisely because they are pointed.
  • The dual cone is not the complement or the polar without a sign flip. uses ; the polar cone uses . Mixing the sign conventions flips the direction of the dual inequality and corrupts the dual program.
  • Generalized inequality is only a partial order. For it can happen that neither nor (their difference is indefinite). Statements that silently assume comparability of all PSD matrices are false.

Key theorem with proof Intermediate+

The organizing result is conic duality: weak duality is unconditional and identical in form across LP, SOCP, and SDP, while strong duality holds under a conic Slater condition, specializing the Lagrangian strong-duality theorem of 44.02.01 to a proper-cone constraint.

Theorem (conic weak and strong duality).

(i) (Weak duality.) For any primal-feasible and dual-feasible with , the primal cost dominates the dual cost: , and the conic duality gap equals . Hence .

(ii) (Strong duality under Slater.) Suppose the primal is strictly feasible: there exists with and (the conic Slater condition), and is finite. Then and the dual optimum is attained.

Proof of (i). Let be primal feasible (, ) and dual feasible (). The gap computation is direct:

Since and , the definition of the dual cone gives . Therefore . Taking the supremum over dual-feasible on the left and the infimum over primal-feasible on the right yields .

Proof of (ii). Write the conic constraint as a generalized inequality , a -convex constraint (the map is -convex because is a convex cone). The problem is then a convex program with a generalized-inequality constraint, and its Lagrangian uses a multiplier paired through the dual cone:

The dual function is , finite only when , in which case it equals . Eliminating reproduces the conic dual .

The conic Slater point with is exactly a strictly feasible point for the generalized-inequality constraint . The strong-duality theorem of 44.02.01, in its generalized-inequality form, then applies: a strictly feasible point in the relative interior of the cone constraint, with finite optimal value and convex data, forces a zero duality gap and attainment of the dual optimum. Concretely, the value set

is convex; is a boundary point with nothing in strictly below it at ; the supporting-hyperplane theorem (finite-dimensional Hahn-Banach) produces a supporting functional whose cone component lands in , and Slater rules out a vertical support exactly as in 44.02.01. Normalizing gives the multiplier and slack with .

Bridge. This theorem builds toward the interior-point algorithms of 44.05.04 and appears again in every optimality certificate for LP, SOCP, and SDP. This is exactly the Lagrangian strong duality of 44.02.01 read through a proper cone: the dual multiplier that was nonnegative for an inequality constraint is now required to lie in the dual cone , so the per-coordinate sign condition generalises to the cone membership , and the complementary-slackness scalar generalises to the conic complementarity . The foundational reason strong duality can fail for SDPs while never failing for the polyhedral LP case is the same vertical-support phenomenon of 44.02.01: the orthant's polyhedral closedness makes LP duality automatic (the Farkas mechanism of 44.02.04), whereas the curved PSD cone needs the interior-point Slater hypothesis to exclude a vertical support. Putting these together, weak duality is dual to primal feasibility and is the single identity across all three cones, while strong duality is the separation theorem applied to the value set over the cone, and the dual slack is the price vector that closes the conic gap.

Exercises Intermediate+

Advanced results Master

Conic duality theory: the gap can be nonzero, finite, or unattained independently

Conic duality refines the Lagrangian picture of 44.02.01 by separating three phenomena the orthant conflates. For a closed convex cone the weak-duality identity holds verbatim, but strong duality splits into independent guarantees: finiteness of , attainment of the dual optimum, and a zero gap. The Ben-Tal-Nemirovski conic-duality theorem [Ben-Tal & Nemirovski Lec. 2] states that if the primal is strictly feasible and bounded below then the dual is solvable with equal optimal value; symmetrically, dual strict feasibility yields primal solvability. When neither problem is strictly feasible the gap may be positive even with both values finite — the curved-cone pathology absent from LP. The polyhedral case is special: when is the nonnegative orthant the cone is finitely generated, the value set is a polyhedron whose closedness is unconditional, and strong duality and the Farkas certificates of 44.02.04 hold with no interior hypothesis at all.

The Schur complement as the universal SDP-modeling primitive

The single most-used SDP-modeling tool is the Schur complement [Vandenberghe & Boyd §6]. For a symmetric block matrix with ,

and iff additionally the complement is strictly positive. This converts a wide class of nonlinear convex constraints into linear matrix inequalities: a quadratic-over-linear bound becomes ; a matrix-fractional or trace-quadratic objective, a Lorentz-cone constraint (the inclusion above), and the convexity of the perspective and matrix-fractional functions all reduce to one Schur-complement LMI. The proof factors through the congruence , and congruence preserves inertia (01.01.13), so the PSD-ness of matches that of the block-diagonal form.

Eigenvalue optimization and the Lovász theta function

The PSD cone's distinctive power is spectral. Beyond minimizing , the sum of the largest eigenvalues is SDP-representable, as is minimizing the spectral norm, the nuclear norm, and the condition number, and bounding the spectrum into an interval. The Lovász theta function [Lovász 1979] — the SDP — sandwiches between the clique and chromatic numbers, , and is computable to arbitrary accuracy by interior-point methods even though the integer quantities it bounds are NP-hard. The same relaxation philosophy yields the Goemans-Williamson max-cut SDP [Goemans & Williamson 1995] with its randomized-rounding guarantee, the first demonstration that an SDP relaxation can beat every prior combinatorial bound and the template for sum-of-squares and Lasserre hierarchies.

Self-dual cones, symmetric cones, and the limits of the framework

The three workhorse cones are not merely self-dual but symmetric: self-dual and homogeneous (their automorphism group acts transitively on the interior). The Koecher-Vinberg theorem identifies symmetric cones with the cones of squares in Euclidean Jordan algebras, and the classification gives exactly five families, of which the orthant, the Lorentz cone, and the real/complex/quaternionic PSD cones are the cases used in practice. Self-concordant barriers exist for these cones with barrier parameter governing interior-point complexity — for each orthant coordinate, for the Lorentz cone, for the PSD cone — which is why all three admit the same polynomial-time path-following solver of 44.05.04. Cones outside this class (the copositive cone, the cone of nonnegative polynomials in several variables) are convex and proper but lack tractable barriers; optimizing over them is NP-hard, marking the boundary between the conic problems that the LP/SOCP/SDP theory solves and the convex problems it can only approximate.

Synthesis. The conic program is exactly the linear program with its sign constraint replaced by membership in a proper cone, and the central insight is that weak duality is the single cone-pairing identity valid across LP, SOCP, and SDP, so the entire duality theory of 44.02.01 transports to the cone with the dual multiplier living in . The foundational reason the three canonical cones admit one unified solver is that each is self-dual and symmetric with a self-concordant barrier, which generalises the per-coordinate logarithm of the LP barrier to the Lorentz and log-determinant barriers and feeds the interior-point machinery of 44.05.04. Putting these together, the Schur complement is the engine that realizes the inclusion LP SOCP SDP, turning each Lorentz constraint into a linear matrix inequality, while eigenvalue optimization and the max-cut and Lovász relaxations show the strict gain in expressive power at the SDP level. This is exactly the same separation engine of 44.02.01 that is dual to primal feasibility, with the polyhedral orthant case is dual to the Farkas certificates of 44.02.04 needing no Slater hypothesis; the bridge is the dual cone , the single object through which generalized inequalities, conic duality, complementary slackness, and the barrier-based interior-point algorithms of 44.05.04 become one theory, and it appears again in robust optimization, combinatorial relaxation, and control via linear matrix inequalities.

Full proof set Master

Proposition 1 (the three canonical cones are proper and self-dual). The nonnegative orthant , the Lorentz cone , and the PSD cone are each closed, pointed, solid convex cones equal to their own dual under the standard (resp. trace) inner product.

Proof. Each is a convex cone closed under nonnegative scaling and addition; closedness follows from being a sublevel/intersection set of continuous functions, solidity from exhibiting an interior point (, , respectively), and pointedness because each contains no line. Self-duality of is Exercise 1. For the self-duality under is proved in 44.01.01 via the spectral decomposition and rank-one test matrices . For the Lorentz cone, let , so . For any with , Cauchy-Schwarz gives , so . Conversely, if satisfies for all , take and (for ): then , optimizing over the choice forces , so . Hence is self-dual.

Proposition 2 (conic weak duality and conic complementary slackness). For primal-feasible and dual-feasible , ; equality (zero gap) holds iff , the conic complementarity condition.

Proof. The identity is the computation in part (i) of the Key theorem, using and . Nonnegativity is the dual-cone definition applied to , . The gap of the pair is exactly , so it vanishes iff . For the PSD cone with , together with forces : writing with , with each term nonnegative forces on the range of , hence .

Proposition 3 (Schur-complement criterion). Let be symmetric with . Then , and .

Proof. Let , which is invertible. The congruence is verified by block multiplication. Congruence by an invertible matrix preserves inertia (Sylvester's law, 01.01.13), so and the block-diagonal matrix have the same numbers of positive, zero, and negative eigenvalues. Since , the block-diagonal matrix is PSD iff and PD iff . Transferring through the inertia-preserving congruence gives the claim for .

Proposition 4 (eigenvalue-sum SDP representation). For symmetric , the sum of the largest eigenvalues satisfies , an SDP in .

Proof. Fix and minimize over . The constraints and mean in the PSD order; the minimizing has eigenvalues in the eigenbasis of , giving . The inner objective becomes . This is convex and piecewise linear in ; its minimum is attained at , where , because the terms with contribute zero and the first terms sum as shown. Hence the SDP value equals . The case recovers of Exercise 5.

Connections Master

  • Lagrangian duality and Slater's condition 44.02.01 is the direct parent: conic duality is its specialization to a single generalized-inequality constraint , with the nonnegative scalar multiplier replaced by a dual-cone multiplier , weak duality becoming the pairing , and the conic Slater condition replacing the strictly-feasible interior point. The strong-duality proof here invokes that unit's separating-hyperplane argument over the value set verbatim.

  • The convex-cone and dual-cone calculus of 44.01.01 supplies every structural fact this unit relies on: that the orthant, Lorentz, and PSD cones are proper and self-dual, that makes the primal-dual conic pairing symmetric, and that an affine preimage of a cone is convex, which is what makes the LP, SOCP, and SDP feasible sets convex and the inclusion chain meaningful.

  • The spectral theorem and PSD characterization of 01.01.13 is the foundation of the SDP layer: the PSD cone is defined by nonnegativity of eigenvalues, eigenvalue optimization is SDP-representable through the order , and the Schur-complement proof rests on Sylvester's inertia law, which is the spectral structure of symmetric matrices applied under congruence.

  • Farkas' lemma and the theorems of the alternative 44.02.04 are the polyhedral, nonnegative-orthant case of conic duality: when the value set is a polyhedron whose closedness is automatic, so strong duality and infeasibility certificates hold without any Slater hypothesis, and the conic theory of this unit generalizes those alternatives to an arbitrary closed convex cone via . The interior-point algorithms of 44.05.04 then solve the conic programs assembled here in polynomial time using self-concordant barriers for each canonical cone.

Historical & philosophical context Master

Linear programming and its duality were established by George Dantzig's simplex method (1947) and the duality theorem formalized by John von Neumann and by David Gale, Harold Kuhn, and Albert Tucker around 1951, the same circle that produced the nonlinear-programming optimality conditions cited in 44.02.01. The decisive generalization to nonpolyhedral cones came with the interior-point revolution: Narendra Karmarkar's 1984 polynomial-time LP algorithm prompted Yurii Nesterov and Arkadi Nemirovski to identify the abstract property — self-concordance of a barrier — that makes polynomial-time interior-point methods possible, published in their 1994 monograph [Nesterov & Nemirovski 1994], which framed LP, SOCP, and SDP as instances of one conic-programming theory solvable by one algorithmic principle.

Semidefinite programming as an explicit object was surveyed and popularized by Lieven Vandenberghe and Stephen Boyd in their 1996 SIAM Review article [Vandenberghe & Boyd 1996], which assembled the eigenvalue, control, and Schur-complement applications; the conic-programming synthesis was codified pedagogically by Aharon Ben-Tal and Nemirovski [Ben-Tal & Nemirovski Lec. 1-4]. On the applications side, László Lovász's 1979 theta function [Lovász 1979] was the first influential combinatorial quantity computed by what is now recognized as an SDP, and Michel Goemans and David Williamson's 1995 max-cut relaxation [Goemans & Williamson 1995] established SDP relaxation with randomized rounding as an approximation-algorithm technique, seeding the sum-of-squares and Lasserre hierarchies that extend the PSD cone to polynomial optimization.

Bibliography Master

@book{boyd2004convex,
  author    = {Boyd, Stephen and Vandenberghe, Lieven},
  title     = {Convex Optimization},
  publisher = {Cambridge University Press},
  year      = {2004}
}

@book{bental2001lectures,
  author    = {Ben-Tal, Aharon and Nemirovski, Arkadi},
  title     = {Lectures on Modern Convex Optimization: Analysis, Algorithms, and Engineering Applications},
  series    = {MPS-SIAM Series on Optimization},
  publisher = {SIAM},
  year      = {2001}
}

@book{nesterov1994interior,
  author    = {Nesterov, Yurii and Nemirovski, Arkadi},
  title     = {Interior-Point Polynomial Algorithms in Convex Programming},
  series    = {SIAM Studies in Applied Mathematics},
  number    = {13},
  publisher = {SIAM},
  year      = {1994}
}

@article{vandenberghe1996semidefinite,
  author  = {Vandenberghe, Lieven and Boyd, Stephen},
  title   = {Semidefinite Programming},
  journal = {SIAM Review},
  volume  = {38},
  number  = {1},
  pages   = {49--95},
  year    = {1996}
}

@article{goemans1995improved,
  author  = {Goemans, Michel X. and Williamson, David P.},
  title   = {Improved Approximation Algorithms for Maximum Cut and Satisfiability Problems Using Semidefinite Programming},
  journal = {Journal of the ACM},
  volume  = {42},
  number  = {6},
  pages   = {1115--1145},
  year    = {1995}
}

@article{lovasz1979shannon,
  author  = {Lov\'asz, L\'aszl\'o},
  title   = {On the Shannon Capacity of a Graph},
  journal = {IEEE Transactions on Information Theory},
  volume  = {25},
  number  = {1},
  pages   = {1--7},
  year    = {1979}
}