---
title: 'PhaseLift: Convex Approach to Phase Retrieval'
url: https://www.emergentmind.com/topics/phaselift
type: topic
---

# PhaseLift: Convex Approach to Phase Retrieval

Searching arXiv for recent and foundational PhaseLift papers to ground the article.
arxiv_search.query({"search_query":"all:PhaseLift", "start":0, "max_results":10, "sort_by":"relevance", "sort_order":"descending"})
PhaseLift is a convex programming approach to phase retrieval that recasts magnitude-only quadratic measurements as linear observations of a rank-one positive semidefinite matrix. Given measurements of the form \(b_i = |\langle a_i,x\rangle|^2\), the method introduces the lifted variable \(X = x x^*\), so that \(b_i = \mathrm{Tr}(a_i a_i^* X)\). The nonconvex rank-one feasibility problem is then relaxed to a trace-minimization semidefinite program over the PSD cone, and the signal is recovered up to a global phase from a rank-one factorization or, in the noisy case, from the leading eigenvector of the optimizer [1109.4499].

## 1. Measurement model and lifted formulation

In its standard form, PhaseLift addresses recovery of an unknown vector \(x \in \mathbb{C}^n\) or \(\mathbb{R}^n\) from phaseless measurements
\[
b_i = |\langle a_i,x\rangle|^2,\qquad i=1,\dots,m.
\]
The lifting step defines
\[
X = x x^* \succeq 0,
\]
which converts each quadratic measurement into a linear functional,
\[
b_i = \mathrm{Tr}(a_i a_i^* X).
\]
Writing \(\mathcal{A}(X) = (\mathrm{Tr}(a_i a_i^* X))_{i=1}^m\), one obtains the linear observation model \(b = \mathcal{A}(X)\) with the additional structural constraints \(X \succeq 0\) and \(\mathrm{rank}(X)=1\) [1109.4499].

PhaseLift replaces the nonconvex rank constraint by trace minimization on the PSD cone:
\[
\min_{X}\ \mathrm{Tr}(X)
\quad\text{s.t.}\quad
\mathcal{A}(X)=b,\quad X\succeq 0.
\]
In the presence of additive noise, a standard noise-aware variant is
\[
\min_{X\succeq 0}\ \mathrm{Tr}(X)
\quad\text{s.t.}\quad
\|\mathcal{A}(X)-b\|_2\le \epsilon,
\]
and a common alternative is an \(\ell_1\)-fitting formulation that minimizes the residual \(\sum_i |\mathrm{Tr}(a_i a_i^* X)-b_i|\) subject to \(X\succeq 0\) [1109.4499]. When the optimizer is not exactly rank one, the leading eigenvector of \(\widehat X\), scaled by the top eigenvalue and phased appropriately, is used as the estimate of \(x\) [2509.10300].

The method is intrinsically invariant to a global phase factor: if \(x\) is feasible, then \(e^{i\phi}x\) produces the same magnitude data. This ambiguity is the canonical identifiability limit in phase retrieval and is built into both the theoretical guarantees and the reconstruction procedures [1208.6247].

## 2. Foundational recovery theory

The original exact-recovery theorem for PhaseLift established that if the sensing vectors are sampled i.i.d. uniformly on the unit sphere, then there exist absolute constants \(c_0,\gamma>0\) such that
\[
m \ge c_0\, n\log n
\]
implies, with probability at least \(1-3e^{-\gamma m/n}\), that the unique minimizer is \(X^\star = xx^*\) [1109.4499]. The same work also proved stability under additive noise: if \(b_i = |\langle x,z_i\rangle|^2 + \nu_i\) with \(\|\nu\|_2\le \epsilon\), then the solution \(\widehat X\) to the noise-aware SDP obeys
\[
\|\widehat X - xx^*\|_2 \le C_0\,\epsilon,
\]
and the induced vector estimate satisfies
\[
\|\widehat x - e^{i\phi}x\|_2 \le C_0 \min(\|x\|_2,\epsilon/\|x\|_2)
\]
for some phase \(\phi\) [1109.4499].

A subsequent sharp improvement removed the logarithmic oversampling factor for Gaussian and spherical designs. For sensing vectors drawn i.i.d. either from the complex normal law \(\mathcal N(0,I_n/2)+i\mathcal N(0,I_n/2)\) or uniformly from the sphere of radius \(\sqrt n\), there exist absolute constants \(c_0,\gamma>0\) such that
\[
m \ge c_0\, n
\]
implies, with probability at least \(1-Ce^{-\gamma m}\), that the PhaseLift SDP is exact for every \(x_0\in\mathbb C^n\) simultaneously [1208.6247]. In the noisy setting, the same paper analyzed the \(\ell_1\)-fitting program
\[
\min_{X\succeq 0}\ \sum_{i=1}^m \bigl|\mathrm{Tr}(a_i a_i^* X)-b_i\bigr|
\]
and proved
\[
\|\widehat X - x_0x_0^*\|_F \le C_0\,\|w\|_1/m,
\]
together with a corresponding leading-eigenvector error bound. These bounds were described there as information-theoretically optimal up to constants [1208.6247].

The proof architecture that became standard in the PhaseLift literature combines geometric analysis of the tangent space at \(xx^*\), injectivity or \(\ell_1\)-isometry properties of the measurement operator on that cone, and construction of dual certificates in the range of \(\mathcal A^*\) [1109.4499]. In the Gaussian \(m=O(n)\) regime, one formulation of the certificate conditions is
\[
Y_{T^\perp}\preceq -I_{T^\perp},\qquad \|Y_T\|_F \le \tfrac12,
\]
with \(T\) the tangent space at \(x_0x_0^*\); a net argument then yields universality over all signals [1208.6247].

## 3. Structured measurements, coded diffraction, and derandomization

A central development was the extension of PhaseLift from unstructured Gaussian measurements to coded diffraction patterns (CDP), a model tailored to Fourier magnitude measurements with random masks. In this setting,
\[
a_{k,\ell} = D_\ell f_k,\qquad m=nL,\qquad
y_{k,\ell}=|f_k^* D_\ell^* x|^2 + w_{k,\ell},
\]
where \(f_k\) is a Fourier vector and \(D_\ell\) is a diagonal mask [2509.10300]. Under admissible random masks, exact recovery was proved when the number of patterns satisfies
\[
L \ge c\,\log^4 n,
\]
equivalently \(m=O(n\log^4 n)\), with probability at least \(1-n^{-1}\) [1310.3240]. This result showed that physically realistic masked Fourier measurements admit rigorous convex recovery guarantees.

The same coded-diffraction line has now been sharpened in the noisy regime. Under adversarial noise \(w\in\mathbb R^m\) with \(\|w\|_2\le \tau\), and masks whose entries are i.i.d. symmetric, bounded, and satisfy \(\mathbb E|d|^2=\nu\), \(\mathbb E|d|^4=2\nu^2\), Theorem 3.1 in the 2025 analysis proves that if
\[
L\ge C_0\log^2 n,
\]
then
\[
\|\widehat X - xx^*\|_F
=
O\!\Bigl(\log n\;\frac{\|w\|_2}{\sqrt m}\Bigr).
\]
For mean-zero sub-Gaussian noise with \(\|w\|_{\psi_2}\le \sigma\), and a modified PhaseLift with the additional rank bound \(\mathrm{Tr}(X)\le \|x\|_2^2\), the bound becomes
\[
\|\widehat X - xx^*\|_F
=
O\!\Bigl(\sigma\sqrt{\frac{n\log^4 n}{m}}\Bigr)
\]
with probability \(1-e^{-\omega}\) when \(L\ge C_0\,\omega\,\log^2 n\) [2509.10300]. The same work states that these results settle the conjecture of Soltanolkotabi that CDP PhaseLift should achieve the \(\|w\|_2/\sqrt m\) scaling of the Gaussian-design case, up to logarithmic factors [2509.10300].

PhaseLift does not fundamentally require Gaussian measurement vectors. A partial derandomization result replaced Gaussian draws by iid samples from a complex projective \(t\)-design. For any exact \(t\)-design with \(t\ge 3\), exact recovery holds with probability at least \(1-e^{-\omega}\) provided
\[
m \ge C\,\omega\, t\, d^{1+2/t}\,\log^2 d.
\]
If \(t\approx 2\log d\), this becomes \(m=O(d\log^3 d)\), which matches the Gaussian \(O(d)\) scaling up to polylog-factors [1310.2267]. The same work also proves a converse for arbitrary 2-designs: in prime-power dimensions, a generic 2-design can force \(m=\Omega(d^2)\), so \(t\ge 3\) is essential for near-linear scaling [1310.2267].

There are also deterministic constructions. One explicit frame of cardinality \(5n-6\) was shown to be 1-complete, meaning that every pure state \(xx^*\) is the unique minimizer of the PhaseLift SDP from its associated intensity measurements [1508.00522]. This establishes exact deterministic phase retrieval via PhaseLift with an explicit measurement design, rather than probabilistic sensing.

A recurring misconception is that Fourier oversampling alone suffices to make phase retrieval well posed. A contrary result in the structured-illumination literature shows that even heavy Fourier oversampling, up to fivefold in each dimension, can leave an ill-posed problem: both PhaseLift and classical alternating-projection methods fit the oversampled modulus perfectly yet diverge spectrally from the true solution, with relative MSE \(\approx 0.5\) [1109.0573]. The same line of work introduced multiple structured illuminations and showed that three suitably designed diffraction patterns can uniquely determine the phase in 1D and 2D under nonvanishing Fourier and gcd conditions [1109.0573].

## 4. Stability, adversarial noise, and outliers

Noise robustness is a defining theme of the PhaseLift literature, but the precise scaling depends strongly on the sensing model and loss function. In the original random-sphere analysis, the Frobenius or operator-norm error scales linearly with the additive perturbation level in the constrained \(\ell_2\) formulation [1109.4499]. In the Gaussian \(m=O(n)\) regime, the \(\ell_1\)-fitting SDP achieves the sharper bound \(\|\widehat X-x_0x_0^*\|_F \le C\|w\|_1/m\), which improves prior \(\ell_2\)-noise results by an additional \(1/\sqrt m\) factor in the denominator [1208.6247].

For coded diffraction patterns, the recent near-optimal noisy theory makes explicit why logarithmic factors appear. The proof proceeds through robust injectivity on the tangent space \(T\), where
\[
\|\mathcal A(Z)\|_2 \gtrsim \sqrt m\,\|Z\|_F
\quad\text{for all }Z\in T,
\]
and a dual certificate constructed by a golfing scheme. If \(\widehat X-xx^*=Z\), feasibility and optimality imply \(\mathrm{Tr}(Z)\le 0\) and \(\|\mathcal A(Z)\|_2\le 2\tau\); pairing \(Z\) with the dual certificate yields
\[
\|Z\|_F \lesssim (\log n)\,\tau/\sqrt m.
\]
The logarithmic losses arise from union bounds over \(O(\log n)\) golfing steps and concentration inequalities with tail-probability unions over \(nL\) measurements [2509.10300].

Robustness to gross corruptions led to a distinct branch of PhaseLift variants that eliminate the explicit trace term and rely on PSD structure plus \(\ell_1\) residual minimization. In the real Gaussian model
\[
b_i = |\langle a_i,x_0\rangle|^2 + \eta_i + \delta_i,
\]
where \(\eta_i\) is small noise and \(\delta_i\) are arbitrary gross errors on an unknown support, the convex program
\[
\min_{X\succeq 0}\ \sum_{i=1}^m |a_i^T X a_i - b_i|
\]
was shown to tolerate a small, fixed fraction of gross errors with only \(m=O(n)\) measurements. Specifically, there are constants \(\alpha_0,\gamma,C\) such that if \(|\{i:\delta_i\neq 0\}|\le \alpha_0 m\), then with probability at least \(1-e^{-\gamma m}\),
\[
\|\widehat X - x_0x_0^T\|_F \le C\,\|\eta\|_1/m
\]
[1502.04241]. The paper emphasizes that, in this robust formulation, the rank penalization through a trace term is unnecessary and no tuning parameter needs to be chosen [1502.04241].

The adversarial-outlier threshold has since been quantified sharply. For Gaussian measurements perturbed by \(\ell_1\)-bounded noise and an adversarially chosen sparse vector with support size at most \(sm\), the Robust-PhaseLift model
\[
\min_{X\succeq 0}\ \|\mathcal A(X)-b\|_1
\]
succeeds for any \(s<s^*\approx 0.1185\) with \(\mathcal O(n)\) measurements, even when the sparse outliers may depend on the measurement and the observation [2311.13219]. The analysis is based on the robust outlier bound condition on the tangent space and on the distribution of the product of two correlated Gaussian variables. The same work constructs adaptive counterexamples showing failure when \(s>s^*\) with high probability [2311.13219].

## 5. Variants, priors, and related convex formulations

Because the trace surrogate is only an indirect proxy for rank one, a number of variants refine the penalty or incorporate structure beyond low rank. One example is PhaseLiftOff, which replaces the trace-only bias by the nonconvex rank-one proxy
\[
\mathrm{Tr}(X)-\|X\|_F.
\]
On the PSD cone this quantity is nonnegative and vanishes if and only if \(\mathrm{rank}(X)=1\). The corresponding regularized least-squares problem
\[
\min_{X\succeq 0}
\ \frac12\|\mathcal A(X)-y\|_2^2
+\rho\bigl(\mathrm{Tr}(X)-\|X\|_F\bigr)
\]
is equivalent to rank-one least squares under the condition
\[
\|y\|_2>\|e\|_2,
\qquad
\rho>\frac{\|\mathcal A\|\,\|e\|_2}{\sqrt2-1},
\]
and, for i.i.d. complex Gaussian measurements, obeys the stable bound
\[
\|X^{\mathrm{opt}}-\widehat X\|_F
\le
C_\alpha\,\frac{\|e\|_2}{\sqrt m}
\]
when \(m>\theta[\alpha^{-2}\log(1/\alpha)]n\) [1406.6761]. Computation is carried out by a convergent difference-of-convex algorithm, with each convex subproblem solved by ADMM [1406.6761].

When the unknown signal is a superposition of a few Vandermonde atoms, the PhaseLift relaxation can be coupled to atomic norm minimization. In ANM-PhaseLift, one combines the lifted PSD variable with the atomic norm \(\|x\|_{\mathcal A}\) for line spectra. The classical PhaseLift sample requirement is \(M\ge C_0 N\) under i.i.d. Gaussian sensing, while atomic-norm recovery from linear measurements can achieve \(M_{\mathrm{lin}}\ge C_1 L\log N\) under the minimum separation condition
\[
\Delta_{\min}\ge 1/\lfloor (N-1)/4\rfloor.
\]
Simulations in the quadratic setting verify recovery from \(M\) as small as \(O(L\log N)\) samples, substantially below \(N\), although exact RIP-type bounds were left to future work [1808.01036].

PhaseLift has also been generalized from vectors to unitary matrices. In the unitary setting, the unknown is \(U\in\mathbb C^{d\times d}\), the measurements are \(y_i=|\mathrm{Tr}(C_i^\dagger U)|^2+\epsilon_i\), and the lifted variable is \(\Gamma=\mathrm{vec}(U)\mathrm{vec}(U)^\dagger\) with additional Jamiolkowski-state constraints. With unitary 4-design measurements, the noisy PhaseLift program achieves uniform exact recovery with a near-optimal number of samples; with unitary 2-designs, it recovers almost all signals up to a constant additive error under a non-spikiness condition [1510.08887]. This matrix-valued extension is directly motivated by quantum process tomography.

These variants show that “PhaseLift” names both a specific trace-minimization SDP and a broader lifting paradigm. A plausible implication is that many later methods should be understood as structured or geometry-aware descendants of the original convex lifting idea rather than as outright replacements.

## 6. Applications and empirical practice

PhaseLift has been integrated into inverse problems well beyond textbook phase retrieval. In a stochastic time-fractional diffusion inverse problem, the unknown discrete source \(x\in\mathbb C^N\) is recovered from phaseless Fourier data using two masks, \(M_1=I_N\) and \(M_2=\mathrm{diag}(\omega_1,\dots,\omega_N)\) with \(\omega_j\) drawn i.i.d. from Bernoulli\(\{0,1\}\). The implementation uses the TFOCS library and an accelerated proximal scheme of Auslender–Teboulle/Nesterov type, with stopping criterion
\[
\frac{\|\mathcal A(ff^*)-b\|_2}{\|b\|_2}\le 10^{-6}.
\]
For the reported numerical experiments, the typical problem size is \(N=65\), \(M=130\), and reconstructions of \(F(t)=\sin(t)e^{-t/6}\) and \(F(t)=\sin(2t)\cos(3t)\), both corrupted by \(5\%\) uniform noise in the variance estimate, agree closely with the true \(|F(t)|\) for \(\alpha=0.4\) and \(\alpha=0.8\) [2012.11042].

In integrated photonics, PhaseLift was used to reconstruct the transfer matrix of linear-optical networks from classical coherent-light intensity data. For each output port \(j\), the measurement model is
\[
I_j(\alpha)=|\langle M_j,\alpha\rangle|^2+\epsilon_j,
\]
so each row of the transfer matrix reduces to a standard phase-retrieval instance. The paper employs a PSD-constrained \(\ell_1\)-loss formulation without an explicit trace penalty, proves a stable bound
\[
\|Z^\sharp-xx^*\|_2 \le C'\,\|\epsilon\|_1/m,
\]
and reports that for random \(5\times 5\) unitaries the average circuit fidelity between the PhaseLift reconstruction and a two-photon reference reconstruction is \(0.993\) [2010.00517]. The simulated phase transition appears near \(m\approx 4n\) even with \(5\%\) additive Gaussian noise on intensities [2010.00517].

In photoemission orbital tomography, Robust Sparse PhaseLift augments the trace term by an entrywise \(\ell_1\) sparsity penalty and an \(\ell_1\) data-fidelity term,
\[
\lambda_{tr}\,\mathrm{tr}(C)+\lambda_1\|C\|_1+\sum_{m=1}^M |z_m-[\mathcal A(C)]_m|,
\]
subject to \(C\succeq 0\). The method uses known atomic positions and basis functions inside the measurement operator and is reported to identify molecular orbitals, including three-dimensional phases, from a single PMM. The paper further states that adsorption-induced molecular deformations can be perfectly discriminated with an accuracy of \(0.05\) \([\text{angstrom}]\) [2307.12500].

The broader influence of the lifting paradigm extends into quantum time-series estimation. “Quantum Phaselift” estimates the rank-one matrix \(Z=ff^\dagger\) associated with a discrete Loschmidt amplitude rather than estimating \(f\) directly, proves that a \(O(1)\) bandwidth suffices for generic signals, and demonstrates recovery for signals of size exceeding 100 time points using only a few million measurement shots [2602.09119]. Although this is a distinct framework, it preserves the central PhaseLift principle: move the inference task from nonlinear phase variables to a linear measurement model over PSD rank-one matrices.

Across these applications, a consistent pattern emerges. Whenever the forward model naturally produces magnitude-squared linear functionals, and whenever either randomness, structure, or auxiliary masks make the lifted operator sufficiently injective, PhaseLift provides a mathematically explicit route from phaseless data to stable convex recovery.

Source: https://www.emergentmind.com/topics/phaselift