---
title: Low-Rank Phase Retrieval Overview
url: https://www.emergentmind.com/topics/low-rank-phase-retrieval
type: topic
---

# Low-Rank Phase Retrieval Overview

Searching arXiv for recent and foundational papers on low-rank phase retrieval and lifted low-rank matrix recovery.
Low-rank phase retrieval refers to the recovery of a low-rank matrix, or of a collection of signals arranged as a low-rank matrix, from magnitude-only or intensity-only measurements. Across the literature, the term covers two closely related but distinct viewpoints. In one viewpoint, phase retrieval is lifted to recovery of a rank-one positive semidefinite matrix \(X=xx^*\), and “low-rank phase retrieval” denotes the extension from rank one to general low-rank Hermitian or positive semidefinite matrices under rank-one quadratic measurements [1410.6913], [1610.08070], [1602.02737]. In the other viewpoint, one seeks a matrix \(X=[x_1,\dots,x_q]\) whose columns are multiple signals or frames sharing a low-dimensional subspace, while each column is observed through its own phaseless measurements \(y_{ik}=|\langle a_{ik},x_k\rangle|\) or \(y_{ik}=|\langle a_{ik},x_k\rangle|^2\) [1608.04141], [2006.06198], [1811.01574]. A further extension replaces matrix structure by multilinear tensor structure, using Tucker models or structured tensor lifts when the signal ensemble is naturally multiway or when phaselessness is coupled to additional bilinear structure [2202.08260], [2509.10834]. This suggests that low-rank phase retrieval is best understood as a family of phaseless inverse problems in which rank, subspace, positivity, or tensorial structure supplies the missing information that magnitude-only data alone cannot provide.

## 1. Matrix lifting and the low-rank generalization of phase retrieval

A central formulation begins from Hermitian low-rank matrix recovery under rank-one measurements. The measurements are written as
\[
y_j=\operatorname{tr}(X a_j a_j^*) = a_j^*Xa_j,\qquad j=1,\dots,m,
\]
or, in operator form,
\[
b=\mathcal A(X)+\epsilon,
\qquad
(\mathcal A(Z))_j=\operatorname{tr}(Za_ja_j^*)=a_j^*Za_j.
\]
When \(X\) is rank one and positive semidefinite, \(X=xx^*\), this reduces exactly to phase retrieval because
\[
y_j=\operatorname{tr}(xx^* a_j a_j^*)=|\langle a_j,x\rangle|^2
\]
[1410.6913], [1610.08070], [1612.03108].

This lifted viewpoint is the most direct bridge between ordinary phase retrieval and low-rank recovery. Standard phase retrieval becomes the special case \(r=1\), while higher-rank models allow recovery of general Hermitian or PSD matrices under the same rank-one quadratic sensing architecture [1410.6913]. In this sense, the literature does not merely borrow intuition from phase retrieval; it treats phase retrieval as the rank-one endpoint of a broader low-rank matrix inverse problem.

Several works emphasize that the target may be Hermitian low rank or specifically PSD low rank. The distinction matters algorithmically. For arbitrary Hermitian matrices, the natural convex estimator is nuclear norm minimization,
\[
\min_Z \|Z\|_1
\quad\text{subject to}\quad
\|\mathcal A(Z)-b\|_{\ell_q}\le \eta,
\]
whereas for PSD matrices one may replace the nuclear norm by trace minimization or, in some settings, by a PSD-constrained least-squares or empirical-risk estimator [1410.6913], [1610.08070], [1602.02737]. Since \(xx^*\succeq 0\) automatically, phase retrieval falls into the PSD regime.

A second major formulation treats low-rank phase retrieval directly at the matrix level without lifting a single vector. In that model, the unknown is a matrix
\[
X=[x_1,x_2,\dots,x_q]\in\mathbb C^{n\times q},
\qquad
\operatorname{rank}(X)=r\ll \min(n,q),
\]
and each column \(x_k\) is observed through its own phaseless measurements
\[
y_{i,k}=|\langle a_{i,k},x_k\rangle|,
\qquad
i=1,\dots,m,\;k=1,\dots,q.
\]
Equivalently,
\[
y_k = |A_k^* x_k|,
\qquad k=1,\dots,q,
\]
with a separate sensing matrix \(A_k\) for each column [2202.08260], [1811.01574], [2006.06198], [1608.04141]. This is a multiple-signal or dynamic setting rather than a single lifted rank-one matrix setting, but the rank assumption again couples all columns through a shared low-dimensional subspace.

The literature is careful to distinguish these viewpoints. One line studies low-rank matrix recovery under quadratic rank-one measurements, with phase retrieval as the rank-one PSD special case [1410.6913], [1610.08070], [1612.03108], [1602.02737], [2607.08671]. Another line studies multiple signals or frames that jointly form a low-rank matrix and are measured column by column [1608.04141], [2006.06198], [1811.01574]. Both are properly described as low-rank phase retrieval, but they differ in their unknown object, sensing architecture, and proof methods.

## 2. Measurement models and structural assumptions

The dominant measurement model in the lifted literature uses rank-one sensing matrices \(A_j=a_ja_j^*\). For Gaussian measurements, the vectors \(a_j\) are independent standard Gaussian vectors in \(\mathbb C^n\) or \(\mathbb R^n\), leading to sample complexities of order \(m\gtrsim rn\) for low-rank recovery and \(m\gtrsim n\) in the rank-one phase retrieval case [1410.6913]. Structured alternatives include weighted complex projective \(4\)-designs, random tight frames, Clifford-group orbits, stabilizer states, and heavy-tailed ensembles with only finite \(4+\delta\) moments [1410.6913], [1612.03108], [1610.08070], [2607.08671].

The measurement structure is often decisive. Clifford-orbit measurements form a highly structured family with \(d=2^n\), where the vectors are sampled uniformly and independently from an orbit of the complex Clifford group,
\[
\mathrm{Cl}_n\cdot z=\{Uz:U\in \mathrm{Cl}_n\}/\{e^{i\phi}\}_\phi,
\]
or equivalently from the orbit of projectors
\[
\mathrm{Cl}_n\cdot zz^* = \{Uzz^*U^\dagger:U\in \mathrm{Cl}_n\}.
\]
A particularly important orbit is the orbit of stabilizer states, which includes and generalizes the discrete Fourier basis by replacing linear phases with quadratic ones [1610.08070]. The paper on Clifford orbits argues that stabilizer states are “an ideal model for structured measurements for phase retrieval” because they are highly structured yet still permit recovery guarantees close to Gaussian ones [1610.08070].

In the multi-signal LRPR setting, the measurement model is column-local rather than global over the whole matrix:
\[
y_{ik}=|\langle a_{ik},x_k\rangle|.
\]
The sensing vectors are typically taken i.i.d. Gaussian and independent across both \(i\) and \(k\) [1608.04141], [2006.06198]. This independence across columns is not a minor technical convenience. The literature states explicitly that the gain from low rank critically relies on different measurement vectors for different columns; if the same sensing vectors are reused for all columns, the benefit largely disappears [1608.04141], [2006.06198].

Low-rank phase retrieval also appears in robust and noise-corrupted settings. A prominent PSD model considers
\[
z=\mathcal A(X_0)+\beta+w,
\]
where \(X_0\succeq 0\) is rank \(r\), \(\beta\) is a sparse outlier vector, and \(w\) is bounded dense noise [1602.02737]. Another recent line replaces Gaussian or sub-Gaussian assumptions by heavy-tailed sampling vectors with independent mean-zero variance-one entries and only finite \(4+\delta\) moments, together with the nondegeneracy conditions
\[
\beta:=\min_i \mathbb E|a_i|^4>1,
\qquad
\gamma:=\max_i |\mathbb E a_i^2|<1
\]
[2607.08671]. This substantially weakens the distributional assumptions under which uniform lifted recovery can be proved.

The multi-signal literature imposes additional structure on the unknown matrix. Writing
\[
X^*=U^*\Sigma^*{B^*}^\top,
\]
a standard assumption is right singular vector incoherence, expressed as
\[
\max_k \|b_k^*\|_2 \le \mu \sqrt{\frac{r}{q}}
\]
or, in simplified form under constant condition number,
\[
\max_k \|x_k\|_2^2 \le C \frac{\|X^*\|_F^2}{q}
\]
[2006.06198], [2006.13298]. This prevents a small number of columns from dominating the energy and is analogous to right incoherence in matrix completion.

## 3. Convex recovery methods and uniform guarantees

The foundational convex approach is nuclear norm minimization under the lifted operator \(\mathcal A\). For Gaussian rank-one measurements, if
\[
m\ge C_1 nr,
\]
then with probability at least
\[
1-e^{-C_2 m},
\]
the constrained nuclear-norm program recovers every rank-\(r\) target stably, with error
\[
\|X-X^\sharp\|_2 \le C_3 \frac{\eta}{\sqrt m}
\]
in the PSD case and a corresponding uniform recovery statement for all rank-\(r\) matrices under exact \(4\)-design sampling with
\[
m\ge C_4 nr\log n
\]
[1410.6913]. The same paper emphasizes that these are uniform guarantees: one random draw of the sensing vectors succeeds simultaneously for all rank-\(r\) Hermitian matrices [1410.6913].

For PSD targets, positivity simplifies the program. One may replace nuclear norm minimization by trace minimization over the PSD cone,
\[
\min \operatorname{tr}(Z)
\quad\text{subject to}\quad
Z\succeq 0,\;
\|\mathcal A(Z)-b\|_{\ell_2}\le \eta,
\]
or, in some settings, by the even simpler constrained least-squares program
\[
Z^\sharp=\arg\min_{Z\succeq 0}\|\mathcal A(Z)-y\|_{\ell_q}.
\]
This simplification is especially important for phase retrieval because \(xx^*\succeq 0\) automatically [1410.6913], [1610.08070].

Clifford-orbit recovery gives a uniform, stable theorem for low-rank Hermitian matrices with measurements \(A_k=a_ka_k^*\) drawn from a Clifford orbit. If
\[
m\ge C_1\,\kappa(z,r)\, r\, d\,\log d,
\]
then with probability at least
\[
1-\exp\!\left(-\frac{\gamma m}{\kappa(z,r)}\right),
\]
every Hermitian rank-\(r\) matrix can be recovered stably via nuclear-norm minimization, with
\[
\|Z^\sharp-X\|_2 \le \frac{C_2}{\sqrt r}\,\sigma_r(X) + C_3 \sqrt{\kappa(z,r)}\, d\, m^{-1/q}\,\eta.
\]
For PSD targets, the same sampling rate suffices for the simpler PSD-constrained least-squares estimator, and in the phase retrieval specialization \(r=1\), since \(\kappa(z,1)\le 4\), one obtains
\[
m\ge 4C_1 d\log d
\]
for stable phase retrieval from any Clifford orbit, including stabilizer measurements [1610.08070].

Random tight-frame measurements provide another structured alternative. If the measurement vectors are the rows of a random tight frame multiplied by \(\sqrt m\), and
\[
m \ge C_1 nr,
\]
then with probability at least
\[
1 - 3e^{-C_2 m},
\]
every minimizer of the nuclear norm program obeys
\[
\|X - X^\sharp\|_2 \le \frac{D_1}{\sqrt r}\,\|X_{r,c}\|_1 + D_2 \frac{\eta}{m^{1/q}}
\]
[1612.03108]. In the rank-one case this gives a PhaseLift-type guarantee under random tight-frame sensing.

Heavy-tailed quadratic sampling substantially enlarges the admissible sensing models. Under independent entries with finite \(4+\delta\) moments and the nondegeneracy conditions stated above, both nuclear norm minimization and PSD-constrained empirical risk minimization achieve
\[
m=\mathcal O(rn)
\]
up to moment-dependent constants, with uniform stable and robust recovery guarantees [2607.08671]. This yields lifted phase retrieval guarantees with \(m=\mathcal O(n)\) under similarly weak moment assumptions [2607.08671].

A distinct robust convex formulation addresses sparse arbitrary outliers for PSD targets:
\[
\hat X = \arg\min_{X\succeq 0}\ \|z-\mathcal A(X)\|_1.
\]
If the outlier fraction \(s\) satisfies
\[
s\le \frac{s_0}{r},
\qquad
m\ge c_1 n r^2,
\]
then with probability exceeding
\[
1-\exp(-\gamma m/r^2),
\]
the solution obeys
\[
\|\hat X-X_0\|_{\mathrm F}\le c_2 \frac{r\epsilon}{m}
\]
[1602.02737]. In the rank-one case this is a robust PhaseLift guarantee with sparse outlier corruption.

These convex results establish a broad picture. Low-rank phase retrieval in lifted form admits uniform recovery theorems under Gaussian, \(4\)-design, random tight-frame, Clifford-orbit, and heavy-tailed sensing, with positivity often allowing simplification from nuclear norm minimization to PSD-constrained least squares or empirical risk minimization [1410.6913], [1610.08070], [1612.03108], [1602.02737], [2607.08671].

## 4. Nonconvex formulations, alternating minimization, and landscape analysis

Nonconvex approaches avoid the dimensional cost of lifting or operate directly in a low-rank factorization. In the multi-signal LRPR setting, the principal formulation is
\[
X = U B,
\qquad
x_k = U b_k,
\]
with \(U\in\mathbb R^{n\times r}\) or \(\mathbb C^{n\times r}\) orthonormal and \(b_k\in\mathbb R^r\) or \(\mathbb C^r\) [1608.04141], [2006.06198], [2006.13298]. The most developed algorithmic line is AltMinLowRaP, which alternates between low-dimensional phase retrieval for each \(b_k\), phase estimation, and least-squares updates of the shared subspace \(U\) [2006.06198], [2006.13298].

The improved guarantee for AltMinLowRaP states that, under right incoherence and i.i.d. Gaussian sensing,
\[
mq \gtrsim nr^2 (r+\log(1/\epsilon))
\]
and
\[
m \gtrsim \max(r,\log q,\log n)
\]
suffice for geometric convergence to \(\epsilon\)-accuracy [2006.06198]. More precisely, after
\[
T = C\log(1/\epsilon)
\]
iterations, the subspace error and matrix error satisfy
\[
SubsDist_F(U^*,U^T)\le \epsilon,
\qquad
MatDist(X^*,\hat X^T)\le \epsilon
\]
[2006.06198]. The paper emphasizes that this improves the sample complexity of earlier AltMinLowRaP analysis by a factor of \(r^2\) for the alternating-minimization stage and a factor of \(r\) for initialization [2006.06198].

The earlier LRPR paper proposed LRPR1 and LRPR2 in addition to LRPR-init. LRPR1 performs projected truncated gradient descent: initialize from a spectral estimate, take one truncated-WF-like step per column, then project the resulting matrix onto rank \(r\) [2006.13298]. LRPR2 alternates among phase estimation, coefficient updates, and a least-squares update for the shared subspace under the factorization \(X=UB\) [1608.04141]. Theoretical guarantees in that work focus on initialization, while extensive experiments indicate that LRPR2 has the best empirical sample complexity and LRPR1 offers a favorable speed–accuracy tradeoff [1608.04141].

Phase retrieval of low-rank matrices by anchored regression takes a different nonconvex-to-convex route. Given an anchor matrix \(X_0\) roughly aligned with the target \(X_\sharp\), the paper studies the convex program
\[
\min_X - \operatorname{Re}\,\langle X_0, X \rangle + \lambda \|X\|_*
\]
subject to phaseless convex sublevel constraints. In the rank-1 phaseless blind deconvolution model, if the anchor satisfies
\[
\inf_{\theta\in[0,2\pi)} \|u_0v_0^* - e^{\mathfrak{i}\theta}u_\sharp v_\sharp^*\|_{\mathrm F} \le \delta,
\qquad
\delta\le 0.2,
\]
and
\[
\frac{M}{\log^2 M} \ge C_\delta (d_1+d_2),
\]
then the estimate is stable, with a recovery error proportional to the average noise level divided by \(\|X_\sharp\|_{\mathrm F}\) [1910.11477]. For general-rank real Gaussian sensing, the theorem requires
\[
M \ge C_\delta\, r(d_1+d_2)\log(d_1+d_2)
\]
plus an anchor quality condition and an additional conditioning restriction [1910.11477].

Another nonconvex direction analyzes the quartic Burer–Monteiro objective
\[
f_p(X)=\|y-\mathcal A(XX^*)\|^2
\]
for semidefinite low-rank sensing. In phase retrieval this reduces, for \(p=1\), to
\[
f(x)=\sum_{i=1}^n \bigl(y_i-|\langle a_i,x\rangle|^2\bigr)^2.
\]
A recent landscape analysis shows that it can be helpful to overparametrize the factor rank \(p\), and for phase retrieval with sub-Gaussian measurements, overparametrizing by a factor at most logarithmic in the dimension allows recovery with optimal statistical sample complexity [2505.02636]. More generally, for semidefinite matrix sensing with PSD structure and rank-1 Gaussian measurements, the same framework gives optimal \(rd\)-type sample complexity for higher-rank PSD targets [2505.02636]. This suggests that the PSD lifted structure can make overparametrization beneficial rather than harmful.

A deterministic theory for exact non-convex phase retrieval studies Wirtinger Flow through the lifted rank-one PSD matrix \(xx^H\). It identifies a sufficient condition equivalent to a restricted isometry-type property over rank-1 PSD matrices and shows that, under a uniform concentration bound on the lifted normal operator, spectral initialization and geometric convergence follow [2001.02855]. Although specialized to rank one, this work is conceptually important because it interprets nonconvex phase retrieval as recovery on the manifold of rank-1 PSD matrices and derives a weaker condition than those used in generic low-rank matrix recovery [2001.02855].

The nonconvex literature is therefore heterogeneous. Some methods operate directly on low-rank factors in multi-signal models [2006.06198], [1608.04141]. Others analyze quartic least-squares landscapes in the lifted PSD domain [2505.02636], [2001.02855]. Still others use convex anchored regression over the original matrix variable to avoid higher-order lifting [1910.11477]. Across these lines, initialization quality, phase estimation, and the special geometry of rank-one or PSD structure are recurring themes.

## 5. Bayesian, gauge-dual, and optimization-oriented perspectives

Several works approach low-rank phase retrieval from inference or optimization design rather than from sample-optimal recovery theory alone.

A Bayesian formulation for multi-signal LRPR places a Gaussian–Wishart hierarchical prior on the unknown matrix
\[
\boldsymbol{X}=[\boldsymbol{x}_1,\ldots,\boldsymbol{x}_M]\in\mathbb{C}^{N\times M},
\qquad
\operatorname{rank}(\boldsymbol{X})=r,
\]
with measurement model
\[
y_{p,m} = \big|\boldsymbol{a}_{p,m}^H \boldsymbol{x}_m + w_{p,m}\big|.
\]
The shared precision matrix \(\boldsymbol{\Sigma}\) induces low-rank structure across columns, and a variational EM algorithm alternates among Gaussian posteriors for the columns, a Wishart posterior for the precision, a Gamma posterior for the noise precision, and deterministic updates of the hidden phases [1811.01574]. The paper provides simulation evidence that this variational Bayesian learning approach is less sensitive to initialization and can work even with random initialization, unlike a competing alternating-minimization baseline [1811.01574].

A gauge-duality approach studies the lifted real symmetric phase retrieval problem
\[
a_i^T X a_i = b_i,
\qquad
\operatorname{rank}(X)\le r,
\]
and solves the gauge dual approximately using cheap spectral approximations. The resulting approximate dual information is used to initialize a nonconvex rank-1 least-squares refinement [2006.01014]. The contribution is not a new exact recovery theorem for low-rank phase retrieval, but rather a practical hybrid strategy that extracts the principal eigenspace of the lifted solution from the dual certificate matrix
\[
\mathcal A^*(y)=\sum_i y_i a_i a_i^T
\]
without solving the primal SDP to high accuracy [2006.01014]. In the rank-one phase retrieval case, this acts as an improved initializer for the quartic nonconvex objective.

A related optimization-oriented result studies projected SGD for convex low-rank matrix relaxations over the spectrahedron
\[
\mathcal S_n := \{X \in \mathbb{S}^n \mid X \succeq 0,\ \operatorname{Tr}(X)=1\}.
\]
The paper does not develop a phase-retrieval-specific theorem, but it explicitly lists phase retrieval among the motivating applications. Its main result shows that if a convex lifted relaxation has a low-rank optimum satisfying an eigen-gap condition, then projected SGD with a warm start produces iterates that remain low rank with constant probability, so each projection can be computed using only a low-rank eigendecomposition [2001.11668]. For PhaseLift-type relaxations, this suggests that large-scale convex lifted methods may be computationally viable when the target rank is small.

Another line studies PSD low-rank recovery through nonconvex regularization of the exact rank constraint rather than through nuclear norm minimization. For fixed-rank PSD matrix recovery,
\[
\min_{X\succeq 0,\ \operatorname{rank}(X)\le K}\|\mathcal A(X)-b\|^2,
\]
a nonconvex quadratic-envelope regularization preserves the global minimizers of the original constrained problem when \(\gamma>\|\mathcal A\|^2\), and can be optimized by forward–backward splitting or FISTA [1907.09537]. In the Fourier phase retrieval specialization \(K=1\), the paper shows that oversampling improves the rank of the lifted operator only up to approximately factor \(2\) per dimension and that masked measurements are needed to stabilize the lifted inverse problem beyond that [1907.09537].

These optimization-oriented papers broaden the meaning of low-rank phase retrieval. The topic is not only about identifiability and sample complexity; it also concerns how lifted convex programs, dual certificates, low-rank projections, spectral initializers, and factorized nonconvex objectives can be made computationally practical in the phaseless regime [2006.01014], [2001.11668], [1907.09537].

## 6. Tensor-structured extensions, blind deconvolution couplings, and open directions

Matrix low-rankness is not always the most natural structure. For video or image sequences, reshaping each frame into its spatial dimensions and time axis yields a tensor
\[
\underline{X} \in \mathbb{C}^{n_1\times n_2\times q},
\]
and Tucker structure may reduce the number of unknowns from roughly \(nr+qr\) to
\[
n_1r_1 + n_2r_2 + qr_3 + r_1r_2r_3.
\]
The Tucker-Structured Phase Retrieval method models
\[
\underline{X} = \underline{G}\times_1 D \times_2 E \times_3 F
\]
and alternates among spectral frame initialization, HOSVD, low-dimensional phase retrieval for the temporal factors \(f_k\), phase updates, and CGLS updates for \(D,E,\underline{G}\) [2202.08260]. The paper does not provide formal recovery guarantees, but on real video datasets it reports lower reconstruction error than matrix-based baselines in both underdetermined and over-determined regimes when the Tucker ranks are chosen appropriately [2202.08260]. A plausible implication is that multilinear low-rank structure can regularize low-rank phase retrieval more effectively than matrix rank when the signals are intrinsically multiway.

A different tensor extension arises when phase retrieval is coupled to blind deconvolution. In a partially coherent imaging model,
\[
y_i = |F(x \odot g_i)|^2 \circledast (P_i s),
\]
with \(s=Bh\), the lifted object becomes a third-order structured rank-one tensor
\[
T^\star = x^* \circ x \circ h.
\]
A 2025 paper studies a surrogate tensor sensing problem
\[
\min_{\substack{x\in\mathbb{R}^N,\ \|x\|_2=1\\ h\in\mathbb{R}^K}}
\frac12\|A(x\circ x\circ h)-y\|_2^2
\]
and gives a landscape analysis showing benign population geometry and local linear convergence of Riemannian gradient descent under a tensor restricted isometry property [2509.10834]. This is not a general LRPR theorem, but it shows how phaseless quadratic structure combined with an additional bilinear nuisance parameter naturally leads from low-rank matrices to structured low-rank tensors [2509.10834].

The topic continues to develop in more direct amplitude-based formulations as well. A 2025 paper studies real-valued low-rank phase retrieval from amplitude measurements
\[
\mathbf y = |\mathcal A(\mathbf X_0)|+\eta,
\qquad
\operatorname{rank}(\mathbf X_0)\le r,
\]
and proves that the rank-constrained nonlinear least-squares estimator
\[
\hat{\mathbf X}^r \in \arg\min_{\operatorname{rank}(\mathbf X)\le r}
\big\||\mathcal A(\mathbf X)|-\mathbf y\big\|_2^2
\]
achieves
\[
\min\{\|\hat{\mathbf X}^r-\mathbf X_0\|_F,\|\hat{\mathbf X}^r+\mathbf X_0\|_F\}
\lesssim \frac{\|\eta\|_2}{\sqrt p}
\]
with high probability under Gaussian measurements and
\[
p\gtrsim (m+n)r.
\]
The same paper gives comparable guarantees for constrained and penalized nuclear-norm formulations and introduces a strong restricted isometry property for matrices tailored to amplitude measurements [2509.24699]. This suggests that a more direct, non-lifted matrix LRPR theory can recover near-degree-of-freedom sample complexity for Gaussian linear maps, at least in the real-valued setting [2509.24699].

Several open problems recur across the literature. The multi-signal LRPR survey lists reduction of sample complexity from \((n/q)r^3\) toward \((n/q)r^2\), removal of sample splitting, extension to practical Fourier and masked-Fourier measurements, and treatment of dynamic low-rank evolution as major directions [2006.13298]. Tensor-structured LRPR lacks recovery and convergence theorems and still relies on trial-and-error rank selection [2202.08260]. Anchored regression for general-rank LRPR still requires restrictive anchor and conditioning assumptions [1910.11477]. Overparametrized semidefinite landscape results are currently strongest for PSD sensing and rank-1 sub-Gaussian phase retrieval, not arbitrary complex LRPR models [2505.02636]. Heavy-tailed convex theory covers lifted quadratic sampling, but not the full range of structured non-Gaussian acquisition operators encountered in practice [2607.08671].

Taken together, the literature presents low-rank phase retrieval as a technically diverse area at the intersection of phase retrieval, low-rank matrix recovery, matrix completion, tensor methods, and nonconvex optimization. Its unifying principle is simple: phaseless data are insufficient without structure, and low rank—whether expressed as a PSD lift, a shared subspace, a nuclear norm, a Tucker tensor, or a structured factorization—provides that structure in a form amenable to both analysis and computation [1410.6913], [1608.04141], [2006.06198], [2202.08260], [2509.24699].

Source: https://www.emergentmind.com/topics/low-rank-phase-retrieval