Efficient classical representation and quantum state preparation of complete active space wavefunctions
Pith reviewed 2026-06-26 20:17 UTC · model grok-4.3
The pith
CAS wavefunctions admit an MPS representation with O(d²) bond dimension in the symmetry-adapted basis from the Quantum Paldus Transform, enabling O(d³) state preparation on quantum computers.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
CAS states expanded in the symmetry-adapted basis produced by the Quantum Paldus Transform can be represented as matrix product states with bond dimension scaling as O(d²). The MPS can be loaded efficiently onto a quantum computer and transformed back to the Fock basis via the inverse QPT. Both the second-quantized and first-quantized preparation protocols then achieve gate complexity O(d³).
What carries the argument
The Quantum Paldus Transform, a map from the Fock basis to a symmetry-adapted basis that permits any CAS wavefunction to be written as an MPS of bond dimension O(d²).
If this is right
- State-preparation gate counts for CAS initial states drop from exponential in d to cubic in d.
- The same polynomial scaling applies when the target state is prepared in first quantization.
- Classical storage of CAS wavefunctions becomes feasible for active spaces previously considered intractable.
- The MPS representation can be used directly in hybrid quantum-classical algorithms that require multi-reference starting points.
Where Pith is reading between the lines
- The approach may extend to other point-group or spin symmetries if analogous transforms exist that produce low-bond-dimension bases.
- Polynomial preparation removes a key bottleneck that has limited the active-space sizes treatable on near-term quantum hardware.
- The same structural property could be exploited to design classical tensor-network algorithms that operate entirely inside the QPT basis.
Load-bearing premise
Any CAS wavefunction, once written in the basis produced by the Quantum Paldus Transform, admits a matrix-product-state representation whose bond dimension is bounded by O(d²).
What would settle it
An explicit CAS wavefunction of active-space size d whose minimal bond dimension in the QPT basis exceeds c d² for some constant c independent of d.
read the original abstract
Quantum computers promise to solve the electronic structure problem for a large class of molecules. However, the performance of relevant quantum algorithms hinges on preparing initial states with substantial overlap with the target eigenvector. For classically challenging molecules with strong electron correlation, starting from multi-reference states, such as complete active space (CAS) wavefunctions is necessary. Unfortunately, the most advanced state preparation protocols applied to such states result in a gate complexity that scales exponentially with the active space size $d$. In fact, even encoding a CAS state classically is traditionally believed to be intractable for chemically relevant systems. Here, we draw insights from the recently introduced Quantum Paldus Transform (QPT) to show that there exists an efficient classical representation of CAS states and to design a new state preparation routine outperforming previous ones. The QPT represents a transformation from the Fock basis to a friendlier symmetry-adapted basis. Our main contribution consists in showing that CAS states expanded in this basis can efficiently be represented as a matrix product state (MPS) with a bond dimension scaling as $O(d^2)$. One can then efficiently load the MPS on a quantum computer and use the inverse QPT to transform the state to the Fock basis. Moreover, our method can easily be extended to the efficient preparation of CAS states in first quantisation with similar complexity. Crucially, we demonstrate that the complexity of both state preparation protocols only grows polynomially as $O(d^3)$ , which constitutes to the best of our knowledge an exponential improvement over the state of the art.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that the Quantum Paldus Transform (QPT) maps the Fock basis to a symmetry-adapted basis in which any complete active space (CAS) wavefunction admits an efficient matrix product state (MPS) representation with bond dimension scaling as O(d²). This representation is asserted to enable both classical encoding of CAS states and quantum state-preparation protocols (including an extension to first quantization) whose gate complexity scales only as O(d³), constituting an exponential improvement over prior art that scales exponentially with active-space size d.
Significance. If the asserted O(d²) bond-dimension bound and consequent polynomial preparation cost were rigorously established, the result would address a central bottleneck in preparing multi-reference initial states for quantum algorithms applied to strongly correlated molecules, potentially enabling practical quantum advantage in electronic-structure calculations.
major comments (2)
- [Abstract] Abstract: the claim that 'CAS states expanded in this basis can efficiently be represented as a matrix product state (MPS) with a bond dimension scaling as O(d²)' is stated without derivation, explicit construction of the MPS tensors, or any argument bounding the Schmidt rank across cuts. This structural property is load-bearing for both the classical representation and the O(d³) preparation complexity.
- [Abstract] Abstract: the statement that 'the complexity of both state preparation protocols only grows polynomially as O(d³)' is presented without a proof sketch, complexity analysis, or numerical verification for any value of d, rendering the exponential-improvement claim impossible to assess.
Simulated Author's Rebuttal
We thank the referee for their comments on our manuscript. The abstract summarizes the central results, with full derivations, explicit constructions, and complexity analyses provided in the body of the paper. We respond point-by-point to the major comments below.
read point-by-point responses
-
Referee: [Abstract] Abstract: the claim that 'CAS states expanded in this basis can efficiently be represented as a matrix product state (MPS) with a bond dimension scaling as O(d²)' is stated without derivation, explicit construction of the MPS tensors, or any argument bounding the Schmidt rank across cuts. This structural property is load-bearing for both the classical representation and the O(d³) preparation complexity.
Authors: The abstract is intended as a concise summary. The full manuscript derives the O(d²) bond-dimension bound by explicitly constructing the MPS tensors from the QPT mapping and proving an upper bound on the Schmidt rank across any cut. The argument relies on the fact that the symmetry-adapted basis restricts the possible entanglement structure to at most quadratic growth in d, which is shown by enumerating admissible occupation patterns consistent with the Paldus labels. revision: no
-
Referee: [Abstract] Abstract: the statement that 'the complexity of both state preparation protocols only grows polynomially as O(d³)' is presented without a proof sketch, complexity analysis, or numerical verification for any value of d, rendering the exponential-improvement claim impossible to assess.
Authors: The manuscript contains a gate-count analysis establishing the O(d³) scaling: preparing an MPS of bond dimension O(d²) costs O(d³) gates, and the inverse QPT is implemented with O(d³) gates via a sequence of controlled rotations whose depth is linear in d. A proof sketch appears in the main text. While the paper focuses on the analytical bound rather than numerical benchmarks, the scaling holds for arbitrary d and constitutes the claimed exponential improvement over exponential-in-d methods. revision: no
Circularity Check
No significant circularity; O(d²) MPS bound presented as new contribution
full rationale
The abstract positions the efficient MPS representation (bond dimension O(d²)) of CAS states in the QPT basis as the main contribution, from which the O(d³) state-preparation scaling follows. The QPT is referenced only as the source of the symmetry-adapted basis ('draw insights from the recently introduced Quantum Paldus Transform'), not as the source of the bond-dimension bound itself. No equations, self-definitions, fitted parameters renamed as predictions, or load-bearing self-citations that reduce the central claim to prior inputs by construction appear in the provided text. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption CAS wavefunctions in the symmetry-adapted basis produced by the QPT admit an MPS representation with bond dimension O(d²).
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.