Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2610.38193

Barren Plateaus in Quantum Linear Combinations: Interference and Lie-Algebraic Control

Amir Koutahi, Pouya Kananian, Hans-Arno Jacobsen

quant-phcs.ETcs.LGSubmitted Oct 28, 2026• Announced Oct 29, 2026• 0 citations
arXivHTMLPDF

Abstract

A coherent sum of quantum circuits introduces interference that can alter the training signal. We distinguish a single linear combination of unitaries (LCU) from a circuit that repeatedly alternates branch operations and ancilla mixing. For a single-pass LCU, independent branches with zero matrix means cannot erase the differentiated branch's gradient variance: the total variance is at least the branch variance times the fourth power of its coefficient magnitude. This bound holds at any gate position and does not require the circuit unitaries to be Haar-random. Beyond this lower bound, we derive exact loss and gradient formulas that quantify the additional contributions from interference between branches. One of these terms depends on relative coefficient phases and vanishes when the branch ensemble is invariant under multiplication by $i \mathbb{I}$; a four-qubit matchgate example shows that it can be nonzero. For repeated interleaving with full special-unitary ancilla control, we determine the generated Lie algebra from the associative algebra of the branches, including the available central phase directions. This control result predicts concentration only with additional assumptions on circuit sampling. Matchgate examples exhibit polynomial single-pass LCU variances and exponential concentration under joint Haar-random sampling. We also derive coefficient-selection rules and bounds for postselection, measurement error, and approximate moments.

AI Summary

The paper asks when coherent combinations of parameterized circuits preserve branch gradient signals, distinguishing a single-pass linear combination $A=\sum_\ell\alpha_\ell U_\ell$ from repeated interleaving of controlled branches and ancilla mixers. Its main objective is the unnormalized loss $\mathcal L=\operatorname{Tr}(OA\rho A^\dagger)$. For a parameter appearing only in branch $\ell$, independence and zero matrix means of the other branches imply $\operatorname{Var}(\partial_\theta\mathcal L)\geq|\alpha_\ell|^4\operatorname{Var}(g_\ell^{\rm branch})$. Conditioning on the differentiated branch proves this at any gate position without Haar sampling. Exact loss and measurement-adjacent gradient formulas additionally retain interference moments $\mathbb E|m_{\ell k}|^2$ and $\mathbb E(m_{\ell k}^2)$. Haar averaging evaluates these through commutants and invariant bilinear forms; the latter produces coefficient-phase dependence, illustrated by a four-qubit matchgate example with $\mathbb E(m_{12}^2)=-5/32$. For repeated interleaving, full $\mathfrak{su}(L)$ ancilla control generates $\bigoplus_i\mathfrak{su}(Ld_i)\otimes I_{m_i}$ plus explicitly characterized central directions, where $(d_i,m_i)$ describe the branches’ associative envelope. The proof generates branch products through ancilla commutators. Matchgate examples contrast inverse-polynomial single-pass variances with exponential concentration under joint Haar sampling, including the bound $\operatorname{Var}(\mathcal L_{\rm BE})\leq\lambda^4\operatorname{Tr}(\rho^2)\operatorname{Tr}(O^2)/(L^24^{n-1}-1)$. Reachability alone does not establish the sampling assumptions: no convergence depth or finite-depth crossover is determined. The gradient floor does not transfer automatically to normalized postselected objectives, whose gradients include success-probability derivatives. Correlated branches also require separate analysis. Coefficient-selection rules, measurement-cost estimates, and approximate-moment error bounds supplement the results.

Some signs of AI-generated textStrong assumptions

Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.

Related papers

Unflattening by Flattening -- How Input Distributions Shape Output Variance in Angle-Encoded Circuits

Melvin Strobl, Gabriel Mejia, Eileen Kuehn +1 more

Barren plateaus hinder training of parameterized quantum circuits by making loss gradients exponentially small in the number of qubits. For angle-encoded product states, we show how the input distribution affects output variation through algebraic input purity, i.e. the input's overlap with the circuit's dynamical Lie algebra. With the specified readouts and Haar or exact group 2-design sampling, matchgate circuits on $n$ qubits retain output variance of order $1/n$ for every pure product input. The off-diagonal family instead has zero output on computational-basis inputs at any depth and for every parameter choice. Independent uniform angles yield mean variance of the same order as matchgates. A count of diagonal Pauli strings identifies these zero-output endpoints. For any fixed dataset of nonzero inputs with at least $2n$ coordinates, we construct a classical preprocessing certificate. A shared random rotation is accepted only when the dataset's mean purity passes a computable threshold. This certifies mean output variance of order $1/n$ for the off-diagonal family before circuit execution, with a constant expected number of trials. Binary and ternary weighted encodings also recover the independent-angle mean purity from one uniform scalar input. Numerical experiments examine training behavior and extensions to larger algebras. The guarantee concerns output variance averaged over inputs, while successful learning also depends on the task and information preserved by the encoding.

quant-pharXiv:2610.01446•Oct 2, 2026•0 citations

Trainability of Parametrised Linear Combinations of Unitaries

Nikhil Khatri, Stefan Zohren, Gabriel Matos

A principal concern in the optimisation of parametrised quantum circuits is the presence of barren plateaus, which present fundamental challenges to the scalability of applications, such as variational algorithms and quantum machine learning models. Recent proposals for these methods have increasingly used the linear combination of unitaries (LCU) procedure as a core component. In this work, we prove that an LCU of trainable parametrised circuits is still trainable. We do so by analytically deriving the expression for the variance of the expectation when applying the LCU to a set of parametrised circuits, taking into account the postselection probability. These results extend to incoherent superpositions. We support our conclusions with numerical results on linear combinations of fermionic Gaussian unitaries (matchgate circuits). Our work shows that sums of trainable parametrised circuits are still trainable, and thus provides a method to construct new families of more expressive trainable circuits. We argue that there is a scope for a quantum speed-up when evaluating these trainable circuits on a quantum device.

arXiv:2506.22310•Jun 27, 2025•0 citations

Multi-qubit controlled gate synthesis without T-count overhead in the small-error limit

Soichiro Yamazaki, Seiseki Akibue, Kaoru Sano

We study the Clifford+T synthesis of quantum multiplexers $U=\bigoplus_{i=1}^{2^n}U_i$: an $n$-qubit control register selects which single-qubit gate $U_i\in\mathrm{SU}(2)$ acts on the target qubit. For each block $U_i$, consider the minimum even T-count needed to approximate it within diamond distance $\varepsilon$ without ancillae, and let $m$ be the largest of these costs. For an arbitrary multiplexer $U$, we construct an approximation of $U$ within the same distance using at most $m+O(\sqrt{2^n(m+1)})$ $T$ gates. With one control qubit, the cost is $m+O(1)$ and no ancillae are required. We also prove a lower bound that allows the Clifford+T circuit to use any number of clean ancillae initialized to zero and returned exactly to zero. For fixed $n$ and independent Haar-random blocks, the minimum T-count among these circuits lies between $(3-δ)\log_2(1/\varepsilon)$ and $(3+δ)\log_2(1/\varepsilon)$ for every fixed $δ>0$, except with probability $O(\varepsilon^{c_{n,δ}})$ as $\varepsilon\rightarrow0$. The optimal leading coefficient is therefore the same as for a single Haar-random $\mathrm{SU}(2)$ gate. In the lower-bound proof, unitarity forces pairs of off-diagonal blocks describing transitions in opposite directions to satisfy cancellation relations with residuals of order $\varepsilon^2$. We combine these relations with the arithmetic of Clifford+T circuits to bound the Haar measure of targets compatible with each pair of off-diagonal blocks. As an application, we obtain ancilla-free approximations of Haar-random $\mathrm{SU}(4)$ gates within diamond distance $\varepsilon$, using at most $(9+δ)\log_2(1/\varepsilon)$ $T$ gates for every fixed $δ>0$, except with probability $O(\varepsilon^{c_δ})$ with $c_δ>0$.

quant-pharXiv:2603.14202•Mar 17, 2026•0 citations

Lie-geometric trainability of quantum dynamical systems: avoiding barren plateaus via low-dimensional Lie subalgebras

Haijian Shao, Yujie Wu, Xing Deng +1 more

Abstract We present a Lie-geometric framework for analyzing the trainability of quantum neural networks (QNNs) through the algebraic structure of their generator spaces. We prove that QNNs whose dynamics are restricted to polynomial-dimensional Lie subalgebras of su ( 2 n ) admit a polynomial lower bound on gradient variance, thereby avoiding barren plateaus by construction. In contrast, architectures exploring the full Lie algebra su ( 2 n ) exhibit exponentially vanishing gradient variance, O ( 2 − n ) , due to the exponential growth of the associated unitary orbit dimensions. Our results interpret barren plateaus as a geometric concentration phenomenon arising from high-dimensional Lie-group actions, rather than as an intrinsic limitation of parameterized quantum models. We further illustrate the mechanism using a finite-size collective-spin benchmark and a minimal two-qubit sanity check, both of which are consistent with the theoretical predictions.

Physica Scripta•Jun 15, 2026•1 citation

Unified Theoretical Characterization of Barren Plateaus in Variational Quantum Algorithms

Yalla Jnan Devi Satya Prasad

Barren plateaus-regions of exponentially small gradients-are a major obstacle for training variational quantum algorithms on non-trivial problem sizes. Existing results analyze specific causes in isolation, such as random deep circuits or global cost functions, but lack a unified picture that integrates circuit design, data encoding, cost structure, and noise. In this work, we develop a comprehensive theory of barren plateaus for parametrized quantum circuits. We prove that three mechanisms jointly control gradient scaling: (1) ansatz-induced concentration, where highly expressive or sufficiently deep circuits approximate unitary 2-designs causing gradient variance to scale as O(2 −n); (2) data-encoding-induced concentration, where certain feature maps wash out class-conditional structure; and (3) probleminduced concentration, where cost functions effectively average over exponentially many degrees of freedom. For each mechanism, we derive explicit scaling laws for the gradient variance as a function of qubit count n, depth L, and encoding complexity. We show that appropriately constrained expressibility, locality, and symmetry breaking can prevent gradients from vanishing exponentially and obtain conditions under which gradients remain at worst polynomially small. This yields architectural and initialization design principles for plateau-resistant circuits, including shallow, problem-informed ansätze, layer-wise training, and structured data encodings. Numerical experiments on up to 100 qubits for VQE, QAOA, and quantum classifiers confirm our theoretical predictions and demonstrate 10-100× improvements in effective gradient magnitudes over standard architectures. Our framework subsumes existing barren plateau results and provides a practical toolkit for diagnosing and mitigating trainability issues in variational quantum algorithms.

•Dec 17, 2025•0 citations

Exploiting all ancilla outcomes in linear combinations of unitaries: low-rank recovery and quantum trapdoor functions

Ammar Daskin

The linear combination of unitaries (LCU) is a fundamental quantum computing primitive used to embed non-unitary operators into larger systems via ancilla post-selection. In standard LCU implementations, only the zero-state ancilla outcome is retained, while all other measurement branches are discarded as unusable overhead. We study these discarded parts. Every computational-basis outcome of the ancilla projects the system onto a different linear combination of the target unitaries; collecting these outcome states in a matrix gives a factorization $Φ=CX$, where $C$ is determined by the ancilla preparation and the rows of $X$ are the states $U_t|ψ\rangle$, so that $rank(Φ)\le K$. Although this holds for any LCU circuit; to make the analysis easier, we use a simplified model circuit with Hadamard gates on the index register and a single rotation qubit in place of the coefficient preparation unitary, in which $C$ is explicit. The structure suggests two complementary applications: (i) in an explicit amplitude-access model, classical low-rank matrix completion can reconstruct the full output (including the target) from a subset of its amplitude entries, while computational basis counts provide intensities and generally require additional phase information for full complex recovery; (ii) treating $C$ as a secret key suggests a candidate quantum trapdoor function, whose security remains informal and depends on the difficulty of structured phase retrieval. The scheme thus turns the "junk" ancilla outcomes into a structured resource, possibly opening paths for further applications.

quant-pharXiv:2605.02986•May 5, 2026•0 citations
A coherent sum of quantum circuits introduces interference that can alter the training signal. We distinguish a single linear combination of unitaries (LCU) from a circuit that repeatedly alternates branch operations and ancilla mixing. For a single-pass LCU, independent branches with zero matrix means cannot erase the differentiated branch's gradient variance: the total variance is at least the branch variance times the fourth power of its coefficient magnitude. This bound holds at any gate position and does not require the circuit unitaries to be Haar-random. Beyond this lower bound, we derive exact loss and gradient formulas that quantify the additional contributions from interference between branches. One of these terms depends on relative coefficient phases and vanishes when the branch ensemble is invariant under multiplication by iIi \mathbb{I}iI; a four-qubit matchgate example shows that it can be nonzero. For repeated interleaving with full special-unitary ancilla control, we determine the generated Lie algebra from the associative algebra of the branches, including the available central phase directions. This control result predicts concentration only with additional assumptions on circuit sampling. Matchgate examples exhibit polynomial single-pass LCU variances and exponential concentration under joint Haar-random sampling. We also derive coefficient-selection rules and bounds for postselection, measurement error, and approximate moments.
The paper asks when coherent combinations of parameterized circuits preserve branch gradient signals, distinguishing a single-pass linear combination A=∑ℓαℓUℓA=\sum_\ell\alpha_\ell U_\ellA=∑ℓ​αℓ​Uℓ​ from repeated interleaving of controlled branches and ancilla mixers. Its main objective is the unnormalized loss L=Tr⁡(OAρA†)\mathcal L=\operatorname{Tr}(OA\rho A^\dagger)L=Tr(OAρA†).
For a parameter appearing only in branch ℓ\ellℓ, independence and zero matrix means of the other branches imply Var⁡(∂θL)≥∣αℓ∣4Var⁡(gℓbranch)\operatorname{Var}(\partial_\theta\mathcal L)\geq|\alpha_\ell|^4\operatorname{Var}(g_\ell^{\rm branch})Var(∂θ​L)≥∣αℓ​∣4Var(gℓbranch​). Conditioning on the differentiated branch proves this at any gate position without Haar sampling. Exact loss and measurement-adjacent gradient formulas additionally retain interference moments E∣mℓk∣2\mathbb E|m_{\ell k}|^2E∣mℓk​∣2 and E(mℓk2)\mathbb E(m_{\ell k}^2)E(mℓk2​). Haar averaging evaluates these through commutants and invariant bilinear forms; the latter produces coefficient-phase dependence, illustrated by a four-qubit matchgate example with E(m122)=−5/32\mathbb E(m_{12}^2)=-5/32E(m122​)=−5/32.
For repeated interleaving, full su(L)\mathfrak{su}(L)su(L) ancilla control generates ⨁isu(Ldi)⊗Imi\bigoplus_i\mathfrak{su}(Ld_i)\otimes I_{m_i}⨁i​su(Ldi​)⊗Imi​​ plus explicitly characterized central directions, where (di,mi)(d_i,m_i)(di​,mi​) describe the branches’ associative envelope. The proof generates branch products through ancilla commutators. Matchgate examples contrast inverse-polynomial single-pass variances with exponential concentration under joint Haar sampling, including the bound Var⁡(LBE)≤λ4Tr⁡(ρ2)Tr⁡(O2)/(L24n−1−1)\operatorname{Var}(\mathcal L_{\rm BE})\leq\lambda^4\operatorname{Tr}(\rho^2)\operatorname{Tr}(O^2)/(L^24^{n-1}-1)Var(LBE​)≤λ4Tr(ρ2)Tr(O2)/(L24n−1−1).
Reachability alone does not establish the sampling assumptions: no convergence depth or finite-depth crossover is determined. The gradient floor does not transfer automatically to normalized postselected objectives, whose gradients include success-probability derivatives. Correlated branches also require separate analysis. Coefficient-selection rules, measurement-cost estimates, and approximate-moment error bounds supplement the results.
Barren plateaus hinder training of parameterized quantum circuits by making loss gradients exponentially small in the number of qubits. For angle-encoded product states, we show how the input distribution affects output variation through algebraic input purity, i.e. the input's overlap with the circuit's dynamical Lie algebra. With the specified readouts and Haar or exact group 2-design sampling, matchgate circuits on nnn qubits retain output variance of order 1/n1/n1/n for every pure product input. The off-diagonal family instead has zero output on computational-basis inputs at any depth and for every parameter choice. Independent uniform angles yield mean variance of the same order as matchgates. A count of diagonal Pauli strings identifies these zero-output endpoints. For any fixed dataset of nonzero inputs with at least 2n2n2n coordinates, we construct a classical preprocessing certificate. A shared random rotation is accepted only when the dataset's mean purity passes a computable threshold. This certifies mean output variance of order 1/n1/n1/n for the off-diagonal family before circuit execution, with a constant expected number of trials. Binary and ternary weighted encodings also recover the independent-angle mean purity from one uniform scalar input. Numerical experiments examine training behavior and extensions to larger algebras. The guarantee concerns output variance averaged over inputs, while successful learning also depends on the task and information preserved by the encoding.
A principal concern in the optimisation of parametrised quantum circuits is the presence of barren plateaus, which present fundamental challenges to the scalability of applications, such as variational algorithms and quantum machine learning models. Recent proposals for these methods have increasingly used the linear combination of unitaries (LCU) procedure as a core component. In this work, we prove that an LCU of trainable parametrised circuits is still trainable. We do so by analytically deriving the expression for the variance of the expectation when applying the LCU to a set of parametrised circuits, taking into account the postselection probability. These results extend to incoherent superpositions. We support our conclusions with numerical results on linear combinations of fermionic Gaussian unitaries (matchgate circuits). Our work shows that sums of trainable parametrised circuits are still trainable, and thus provides a method to construct new families of more expressive trainable circuits. We argue that there is a scope for a quantum speed-up when evaluating these trainable circuits on a quantum device.
We study the Clifford+T synthesis of quantum multiplexers U=⨁i=12nUiU=\bigoplus_{i=1}^{2^n}U_iU=⨁i=12n​Ui​: an nnn-qubit control register selects which single-qubit gate Ui∈SU(2)U_i\in\mathrm{SU}(2)Ui​∈SU(2) acts on the target qubit. For each block UiU_iUi​, consider the minimum even T-count needed to approximate it within diamond distance ε\varepsilonε without ancillae, and let mmm be the largest of these costs. For an arbitrary multiplexer UUU, we construct an approximation of UUU within the same distance using at most m+O(2n(m+1))m+O(\sqrt{2^n(m+1)})m+O(2n(m+1)​) TTT gates. With one control qubit, the cost is m+O(1)m+O(1)m+O(1) and no ancillae are required. We also prove a lower bound that allows the Clifford+T circuit to use any number of clean ancillae initialized to zero and returned exactly to zero. For fixed nnn and independent Haar-random blocks, the minimum T-count among these circuits lies between (3−δ)log⁡2(1/ε)(3-δ)\log_2(1/\varepsilon)(3−δ)log2​(1/ε) and (3+δ)log⁡2(1/ε)(3+δ)\log_2(1/\varepsilon)(3+δ)log2​(1/ε) for every fixed δ>0δ>0δ>0, except with probability O(εcn,δ)O(\varepsilon^{c_{n,δ}})O(εcn,δ​) as ε→0\varepsilon\rightarrow0ε→0. The optimal leading coefficient is therefore the same as for a single Haar-random SU(2)\mathrm{SU}(2)SU(2) gate. In the lower-bound proof, unitarity forces pairs of off-diagonal blocks describing transitions in opposite directions to satisfy cancellation relations with residuals of order ε2\varepsilon^2ε2. We combine these relations with the arithmetic of Clifford+T circuits to bound the Haar measure of targets compatible with each pair of off-diagonal blocks. As an application, we obtain ancilla-free approximations of Haar-random SU(4)\mathrm{SU}(4)SU(4) gates within diamond distance ε\varepsilonε, using at most (9+δ)log⁡2(1/ε)(9+δ)\log_2(1/\varepsilon)(9+δ)log2​(1/ε) TTT gates for every fixed δ>0δ>0δ>0, except with probability O(εcδ)O(\varepsilon^{c_δ})O(εcδ​) with cδ>0c_δ>0cδ​>0.
Abstract We present a Lie-geometric framework for analyzing the trainability of quantum neural networks (QNNs) through the algebraic structure of their generator spaces. We prove that QNNs whose dynamics are restricted to polynomial-dimensional Lie subalgebras of su ( 2 n ) admit a polynomial lower bound on gradient variance, thereby avoiding barren plateaus by construction. In contrast, architectures exploring the full Lie algebra su ( 2 n ) exhibit exponentially vanishing gradient variance, O ( 2 − n ) , due to the exponential growth of the associated unitary orbit dimensions. Our results interpret barren plateaus as a geometric concentration phenomenon arising from high-dimensional Lie-group actions, rather than as an intrinsic limitation of parameterized quantum models. We further illustrate the mechanism using a finite-size collective-spin benchmark and a minimal two-qubit sanity check, both of which are consistent with the theoretical predictions.
Barren plateaus-regions of exponentially small gradients-are a major obstacle for training variational quantum algorithms on non-trivial problem sizes. Existing results analyze specific causes in isolation, such as random deep circuits or global cost functions, but lack a unified picture that integrates circuit design, data encoding, cost structure, and noise. In this work, we develop a comprehensive theory of barren plateaus for parametrized quantum circuits. We prove that three mechanisms jointly control gradient scaling: (1) ansatz-induced concentration, where highly expressive or sufficiently deep circuits approximate unitary 2-designs causing gradient variance to scale as O(2 −n); (2) data-encoding-induced concentration, where certain feature maps wash out class-conditional structure; and (3) probleminduced concentration, where cost functions effectively average over exponentially many degrees of freedom. For each mechanism, we derive explicit scaling laws for the gradient variance as a function of qubit count n, depth L, and encoding complexity. We show that appropriately constrained expressibility, locality, and symmetry breaking can prevent gradients from vanishing exponentially and obtain conditions under which gradients remain at worst polynomially small. This yields architectural and initialization design principles for plateau-resistant circuits, including shallow, problem-informed ansätze, layer-wise training, and structured data encodings. Numerical experiments on up to 100 qubits for VQE, QAOA, and quantum classifiers confirm our theoretical predictions and demonstrate 10-100× improvements in effective gradient magnitudes over standard architectures. Our framework subsumes existing barren plateau results and provides a practical toolkit for diagnosing and mitigating trainability issues in variational quantum algorithms.
The linear combination of unitaries (LCU) is a fundamental quantum computing primitive used to embed non-unitary operators into larger systems via ancilla post-selection. In standard LCU implementations, only the zero-state ancilla outcome is retained, while all other measurement branches are discarded as unusable overhead. We study these discarded parts. Every computational-basis outcome of the ancilla projects the system onto a different linear combination of the target unitaries; collecting these outcome states in a matrix gives a factorization Φ=CXΦ=CXΦ=CX, where CCC is determined by the ancilla preparation and the rows of XXX are the states Ut∣ψ⟩U_t|ψ\rangleUt​∣ψ⟩, so that rank(Φ)≤Krank(Φ)\le Krank(Φ)≤K. Although this holds for any LCU circuit; to make the analysis easier, we use a simplified model circuit with Hadamard gates on the index register and a single rotation qubit in place of the coefficient preparation unitary, in which CCC is explicit. The structure suggests two complementary applications: (i) in an explicit amplitude-access model, classical low-rank matrix completion can reconstruct the full output (including the target) from a subset of its amplitude entries, while computational basis counts provide intensities and generally require additional phase information for full complex recovery; (ii) treating CCC as a secret key suggests a candidate quantum trapdoor function, whose security remains informal and depends on the difficulty of structured phase retrieval. The scheme thus turns the "junk" ancilla outcomes into a structured resource, possibly opening paths for further applications.