---
title: 'QFrCoRT: Convex-Roof Optimization Framework'
url: https://www.emergentmind.com/topics/qfrcort
type: topic
---

# QFrCoRT: Convex-Roof Optimization Framework

QFrCoRT denotes a general numerical framework for evaluating a broad class of quantum resource measures defined by convex-roof extension. Introduced in the context of quantum resource theories, it reformulates the optimization over all pure-state ensemble decompositions of a mixed state as an optimization over a complex Stiefel manifold, and then removes the manifold constraint by a polar-projection trivialization so that standard gradient-based optimizers can be used in unconstrained Euclidean space [2406.19683]. The framework is presented as a unified workflow for entanglement, coherence, and magic-state quantification, and is also described as extendable to other convex-roof quantities in quantum information theory.

## 1. Resource-theoretic setting

In the paper’s formulation, a quantum resource theory is written as $\mathcal{R}=(\mathcal{F},\mathcal{O})$, where $\mathcal{F}(\mathcal{H})$ is the set of free states on Hilbert space $\mathcal{H}$, and $\mathcal{O}$ is the set of free operations [2406.19683]. A resource measure $R(\rho)$ is required to be nonnegative, to vanish exactly on free states, and to be monotone under free operations. QFrCoRT targets the class of resource measures obtained by convex-roof extension from a pure-state quantity.

The defining optimization problem is
$$
R(\rho)=\min_{\{p_i,|\psi_i\rangle\}:\,\rho=\sum_i p_i|\psi_i\rangle\langle\psi_i|}\sum_i p_i\,R(|\psi_i\rangle).
$$
This is the formal problem solved by the framework. The difficulty is that the minimization ranges over all ensemble decompositions of a mixed state. The paper characterizes this as a difficult NP-hard decomposition problem and presents QFrCoRT as a route to converting that problem into a numerically tractable unconstrained optimization task [2406.19683].

The significance of this starting point is structural rather than resource-specific. Entanglement of formation, geometric coherence, and stabilizer-based magic measures differ in their pure-state functionals, but they share the same convex-roof architecture. QFrCoRT is built precisely around that shared architecture.

## 2. Stiefel-manifold reformulation of the convex roof

The central reformulation begins from the spectral decomposition of a rank-$r$ state,
$$
\rho=\sum_{j=1}^r \lambda_j|\lambda_j\rangle\langle \lambda_j|.
$$
Any $n$-term pure-state decomposition $\rho=\sum_{i=1}^n p_i|\psi_i\rangle\langle\psi_i|$ is represented באמצעות unnormalized auxiliary states $|\tilde\psi_i\rangle$ satisfying
$$
p_i=\langle \tilde\psi_i|\tilde\psi_i\rangle,
$$
and
$$
|\tilde\psi_i\rangle=\sum_{j=1}^r \sqrt{\lambda_j}\,X_{ij}\,|\lambda_j\rangle.
$$
Here $X$ is an $n\times r$ complex matrix with
$$
X^\dagger X=I_r,
$$
so $X\in \mathrm{St}(n,r)$, the complex Stiefel manifold [2406.19683].

With this parametrization, the convex-roof optimization is equivalent to an optimization over the Stiefel manifold:
$$
\min_{X\in \mathrm{St}(n,r)}\sum_{i=1}^n p_i(X)\,R(|\psi_i(X)\rangle).
$$
This is the framework’s core geometric step. Instead of searching directly over arbitrary decompositions, one searches over matrices with orthonormal columns. The decomposition variables are therefore encoded in a constrained but standard manifold object.

This reformulation is the main unifying device in QFrCoRT. It does not depend on a particular resource theory beyond the requirement that the mixed-state quantity be a convex roof of a pure-state resource functional. That is why the same machinery can be reused across different resource measures.

## 3. Polar-projection trivialization and optimization workflow

QFrCoRT then eliminates the Stiefel-manifold constraint by a trivialization. The chosen trivialization is the polar projection
$$
g(A)=A\left(A^\dagger A\right)^{-1/2},\qquad g:\mathbb{C}^{n\times r}\to \mathrm{St}(n,r),
$$
which maps an unconstrained full-rank matrix $A$ to its polar factor on the Stiefel manifold [2406.19683]. Since $\mathbb{C}^{n\times r}$ has $2nr$ real parameters, the resulting optimization is unconstrained in Euclidean space.

The paper states that the polar-projection trivialization is numerically better than matrix exponential or Euler–Hurwitz parameterizations. It also notes a practical caveat: when $\mathrm{rank}(A)<r)$, the polar map is ill-defined, but such cases are measure-zero and were not encountered in practice. Matrix square roots are computed by eigendecomposition, while gradients through the square root are obtained via a Sylvester equation.

The algorithmic workflow is deliberately uniform. One builds the objective function for the chosen resource, parameterizes the decomposition through $A\mapsto g(A)$, computes gradients by automatic differentiation in PyTorch, and performs optimization using L-BFGS implemented in SciPy. The reported default choices are typically $n=2d$ pure states in the decomposition, $N=3$ optimization repeats to reduce local-minimum issues, and a stopping tolerance of $10^{-14}$ [2406.19683].

The framework is explicitly contrasted with three established classes of methods. SDP methods relax the convex-roof problem into semidefinite constraints and provide lower bounds but scale poorly. Manually derived gradient methods require case-by-case analytical gradients and are not unified. Seesaw strategies alternate optimization blocks but may struggle for nearly pure multipartite states. QFrCoRT’s stated advantage is that the same optimization machinery applies across different measures, with gradients handled automatically.

## 4. Numerical stabilization and differentiable objective design

A substantial part of QFrCoRT is devoted to numerical stability. For entropy-based objectives, direct evaluation can produce $0\log 0$ or division-by-zero issues. The paper therefore rewrites objectives in terms of the unnormalized auxiliary states and truncates small values below a threshold $\varepsilon$ [2406.19683].

For entanglement of formation, the objective is rewritten as
$$
\mathcal{L}=\sum_i p_i E_f(|\psi_i\rangle)
= -\sum_i p_i\,\mathrm{Tr}[\rho_i\ln\rho_i]
= \sum_i p_i\ln p_i - \mathrm{Tr}[\tilde\rho_i\ln \tilde\rho_i],
$$
where
$$
\tilde\rho_i=\mathrm{Tr}_B\!\left[|\tilde\psi_i\rangle\langle \tilde\psi_i|\right].
$$
The same stabilization logic is reported for linear-entropy and stabilizer-purity objectives.

For geometric coherence, the nondifferentiable $\max$ operation is replaced by a log-sum-exp smoothing:
$$
\max_j \tilde\rho_{i,jj} = \lim_{T\to 0^+} T\log\!\left(\sum_j e^{\tilde\rho_{i,jj}/T}\right).
$$
This preserves compatibility with automatic differentiation while keeping the computation within the same unconstrained optimization pipeline [2406.19683].

These implementation choices are not ancillary. They are part of the framework’s claim to being unified in practice, not only in formal geometry. A plausible implication is that the success of the method depends as much on stable differentiable objective construction as on the manifold reformulation itself.

## 5. Resource measures treated and benchmark problems

The paper demonstrates QFrCoRT on entanglement, coherence, and magic-state quantification [2406.19683].

| Resource | Pure-state quantity | Benchmark state or family |
|---|---|---|
| Entanglement of formation | $E_f(|\psi\rangle)=S\!\left(\mathrm{Tr}_B|\psi\rangle\langle\psi|\right)$ | $d\otimes d$ Werner states |
| Linear entropy of entanglement | $E_l(|\psi\rangle)=1-\mathrm{Tr}(\rho_A^2)$ | $3\otimes 3$ Horodecki bound-entangled states |
| Geometric measure of coherence | $C_g(|\psi\rangle)=1-\max_i|\langle i|\psi\rangle|^2$ | Noisy coherent states |
| Stabilizer entropy / linearized version | $M_\alpha(|\psi\rangle)=\frac{1}{1-\alpha}\log P_\alpha(|\psi\rangle)$ | Qubit states, $\alpha=2$ |

For entanglement, the mixed-state entanglement of formation is defined by the convex roof of the von Neumann entropy of the reduced state,
$$
E_f(\rho)=\min_{\{p_i,|\psi_i\rangle\}}\sum_i p_i\,E_f(|\psi_i\rangle),\qquad
E_f(|\psi\rangle)=S\!\left(\mathrm{Tr}_B|\psi\rangle\langle\psi|\right),
$$
with $S(\rho)=-\mathrm{Tr}(\rho\ln\rho)$. The linear entropy of entanglement is treated analogously, with
$$
E_l(|\psi\rangle)=1-\mathrm{Tr}(\rho_A^2),\qquad \rho_A=\mathrm{Tr}_B|\psi\rangle\langle\psi|.
$$
The reported test families are the Werner states
$$
\rho_W(\alpha)=\frac{1}{d^2-d\alpha}(I_{AB}-\alpha F_{AB}),
$$
where $F_{AB}$ is the swap operator, and the $3\otimes 3$ Horodecki bound-entangled states $\rho_{\mathrm H(\alpha)}$.

For coherence, the paper uses the geometric measure of coherence, defined on pure states by
$$
C_g(|\psi\rangle)=1-\max_i|\langle i|\psi\rangle|^2,
$$
and benchmarks it on the noisy coherent state
$$
\rho=p|\psi^+\rangle\langle\psi^+|+(1-p)\frac{I}{d},\qquad
|\psi^+\rangle=\frac{1}{\sqrt d}\sum_{i=1}^d |i\rangle,
$$
for which the quoted analytic value is
$$
C_g(\rho)=1-\frac{1}{d^2}\Big[(d-1)\sqrt{1-p}+\sqrt{1+(d-1)p}\Big]^2.
$$

For magic, the framework is applied to stabilizer entropy. For a pure $n$-qubit state,
$$
M_\alpha(|\psi\rangle)=\frac{1}{1-\alpha}\log P_\alpha(|\psi\rangle),
$$
with stabilizer purity
$$
P_\alpha(|\psi\rangle)=\frac{1}{d_n}\sum_{P\in\mathbb{P}_n} |\langle\psi|P|\psi\rangle|^{2\alpha},\qquad d_n=2^n.
$$
The mixed-state extension uses a convex roof for $P_\alpha(\rho)$, and the paper also discusses the linearized version $M_\alpha^{\mathrm{lin}}=1-P_\alpha$. The reported demonstrations use qubit states with $\alpha=2$, where stabilizer states form an octahedron inside the Bloch sphere.

## 6. Empirical behavior, scope, and nomenclature

The reported numerical results emphasize both accuracy and efficiency [2406.19683]. For $d\otimes d$ Werner states, the computed entanglement of formation agrees closely with analytical results and reproduces the separable/entangled boundary. For the Horodecki $3\otimes 3$ bound-entangled states, the gradient-based method yields accurate upper bounds on the linear entropy of entanglement, complementing SDP lower bounds. For coherence, the method matches analytic values for noisy coherent states and is reported to be substantially faster than SDP; for $d=20$ to $50$, a fit gives time scaling $t\propto d^\beta$ with $\beta=2$ for gradient-based optimization and $\beta=3.6$ for SDP, while numerical errors stay below $10^{-8}$ even at relatively large dimensions. For linear entropy of entanglement, the paper states that polar projection is much faster than PPT/SDP relaxation, which becomes unavailable in larger dimensions such as $3\otimes 4$. For magic, the heatmaps of linear stabilizer entropy clearly separate stabilizer and magic regions, consistent with the Bloch-sphere octahedron picture.

The paper also reports polynomial scaling with system size in the studied examples. For entanglement of formation on Werner states, the parameter count satisfies $2nr\le 4d^4$. This does not alter the NP-hardness of the original decomposition problem, but it does indicate that the practical optimization variables scale in a controlled way for the cases examined.

The stated overall scope is broader than the three showcased resources. The methodology is described as readily extendable to other convex-roof quantities beyond resource theories, including constrained Holevo capacity and quantum Fisher information. This suggests that the unifying principle is the convex-roof structure itself, not any specific operational interpretation of entanglement, coherence, or magic.

In unrelated natural-language-processing literature, the acronym QFrCoRT is also used for the **Quebec-French Corpus of Regional Terms**, a zero-shot multiple-choice benchmark for dialect understanding composed of 171 single-word Quebec regional terms, each paired with 10 answer options [2510.05026]. That usage is distinct from the quantum-information framework discussed above and reflects an acronym collision rather than a conceptual connection.

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