While fermion-to-qubit isomorphisms are mathematically invariant under exact unitary statevector transformations, their physical implementation on noisy intermediate-scale quantum (NISQ) devices introduces severe, practical measurement bottlenecks. In this work, we systematically investigate the measurement efficiency of adaptive derivative-assembled pseudo-trotter variational quantum eigensolver (ADAPT-VQE) algorithms across three critical molecular benchmarks (LiH, BeH2, and H2O) at stretched bond lengths (2xRe) covering distinct point-group symmetries. First, we establish exact representation invariance across all target multi-reference geometries, demonstrating that operator pool gradients and statevector electronic energies remain strictly isomorphic under exact mathematical execution (E(LiH) = -7.691469 Ha, E(BeH2) = -15.110110 Ha, and E(H2O) = -74.566874 Ha). Second, we quantify the hardware measurement overhead by partitioning operator pool observables into qubit-wise commuting (QWC) tensor-product basis (TPB) cliques. We reveal that while both Jordan-Wigner (JW) and Bravyi-Kitaev (BK) encodings yield identical total Pauli term counts across the singlet excitation manifold (128 terms for H2O and BeH2), the hierarchical binary-tree layout of the BK mapping dramatically compresses the required measurement circuits---achieving up to a 45.31% reduction in total hardware execution steps compared to 12.50% under JW. Our findings demonstrate that tree-structured mappers serve as implicit, highly efficient measurement-packing engines for adaptive quantum simulations without altering the underlying variational trajectory.
While fermion-to-qubit isomorphisms are mathematically invariant under exact unitary statevector transformations, their physical implementation on noisy intermediate-scale quantum (NISQ) devices introduces severe, practical measurement bottlenecks. In this work, we systematically investigate the measurement efficiency of adaptive derivative-assembled pseudo-trotter variational quantum eigensolver (ADAPT-VQE) algorithms across three critical molecular benchmarks (LiH, BeH2, and H2O) at stretched bond lengths (2xRe) covering distinct point-group symmetries. First, we establish exact representation invariance across all target multi-reference geometries, demonstrating that operator pool gradients and statevector electronic energies remain strictly isomorphic under exact mathematical execution (E(LiH) = -7.691469 Ha, E(BeH2) = -15.110110 Ha, and E(H2O) = -74.566874 Ha). Second, we quantify the hardware measurement overhead by partitioning operator pool observables into qubit-wise commuting (QWC) tensor-product basis (TPB) cliques. We reveal that while both Jordan-Wigner (JW) and Bravyi-Kitaev (BK) encodings yield identical total Pauli term counts across the singlet excitation manifold (128 terms for H2O and BeH2), the hierarchical binary-tree layout of the BK mapping dramatically compresses the required measurement circuits---achieving up to a 45.31% reduction in total hardware execution steps compared to 12.50% under JW. Our findings demonstrate that tree-structured mappers serve as implicit, highly efficient measurement-packing engines for adaptive quantum simulations without altering the underlying variational trajectory.
The Variational Quantum Eigensolver (VQE) is a leading algorithm for noisy intermediate-scale quantum (NISQ) devices, but its adaptive variants (e.g., ADAPT-VQE) suffer from severe classical simulation bottlenecks during the ansatz growth phase. Representing and exponentiating pool operators for multi-qubit systems constructs massive sparse matrices that quickly scale to millions of elements, choking classical memory bandwidth and CPU/GPU cycle capacity. In this work, we present a resource-efficient software-layer framework that completely bypasses dense matrix exponentiation by evaluating state updates through a deterministic fifth-order (O(5)) Taylor series expansion. This approach reframes the costly unitary evolution into a chained sequence of five lightweight, successive sparse matrix-vector multiplications scaling strictly as O(Nz), where Nz is the number of non-zero elements. We validate our framework using equilibrium BeH2 (14 qubits), equilibrium H2O (12 qubits), and strongly correlated asymmetrically stretched H2O molecular profiles under both Jordan-Wigner (JW) and Bravyi-Kitaev (BK) transformations. The simulation results demonstrate that the O(5) truncation maintains exceptional numerical fidelity, retaining a state fidelity > 0.999999 and matching absolute ground state energies down to subchemical accuracy, while effortlessly navigating matrix spaces exceeding 268 million structural elements. This framework provides a scalable, high-performance pathway for executing deep variational simulations on hardware platforms with constrained computational budgets.
ADAPT-VQE is a leading variational quantum algorithm as it circumvents the choice-of-ansatz conundrum by iteratively growing compact and arbitrarily accurate problem-tailored ansätze. However, for molecular Hamiltonians and hardware-efficient operator pools, the gradient-measurement step of the algorithm requires the estimation of O(N8) observables, which may represent a bottleneck for relevant system sizes on real devices. We present an efficient strategy for measuring the gradients of the three best-performing hardware-efficient operator pools by partitioning the required observables into O(N3)-sized sets of commuting Paulis that can be simultaneously measured with an at most N−3 CNOT overhead. We argue that our approach is robust to shot-noise effects and show that measuring the pool gradients is, in fact, not O(N4), but only about 4N times as expensive as measuring the energy once. Our proposed measurement strategy significantly ameliorates the measurement overhead of ADAPT-VQE and brings us one step closer to practical implementations on real devices.
Fixed-ansatz Variational Quantum Eigensolvers (VQE), such as the Unitary Coupled Cluster with Singles and Doubles (UCCSD) framework, frequently suffer from severe initialization paralyzation and zero-gradient traps when evaluated using the non-local Bravyi-Kitaev (BK) fermion-to-qubit mapping. In this work, we systematically demonstrate how the Adaptive Derivative-Assembled Pseudo-Trotter (ADAPT-VQE) framework shatters these structural limitations across three distinct electronic and geometric molecular pillars: Lithium Hydride (LiH), Hydrogen Fluoride (HF), and Water (H2O), under heavily stretched or asymmetric multi-reference configurations. While conventional UCCSD-VQE flatlines completely at a zero energy shift (0.000000~Ha) due to global phase cancellations inherent to the BK tree structures, our dynamic ADAPT-VQE loop successfully isolates the dominant symmetry-breaking operators using analytical commutator gradients. To bypass the severe O(N3) computational bottlenecks of dense matrix exponentiation and Singular Value Decomposition on larger registers, we implement a highly optimized, vector-based Taylor series expansion state-evolution engine. Our numerical results show that the accelerated ADAPT-VQE framework achieves instant, exact Full Configuration Interaction (FCI) convergence within the very first macro-cycle across all three molecular systems, maintaining absolute numerical stability up to a 12-qubit register space. This study establishes a robust, hardware-efficient path for simulating strongly correlated, highly polarized triatomic chemical environments on near-term local architectures.
Accurate determination of ground-state energies for molecules remains a challenge in quantum chemistry and a cornerstone for progress in fields such as drug discovery and materials design. The Variational Quantum Eigensolver (VQE) represents a leading hybrid quantum-classical paradigm for addressing this challenge; however, its widespread realization is limited by noise and the restricted scalability of current quantum hardware. Achieving efficient simulations on Noisy Intermediate-Scale Quantum (NISQ) devices and forthcoming Fault-Tolerant Application-Scalable Quantum (FASQ) systems demands a detailed understanding of how computational resources scale with molecular complexity and fermion-to-qubit encodings. In this study, resource requirements for VQE implementations employing the Unitary Coupled Cluster Singles and Doubles (UCCSD) ansatz are systematically analyzed. The molecular Hamiltonian is formulated in second quantization and mapped to qubit operators through the Jordan-Wigner (JW), Bravyi-Kitaev (BK), and Parity (Pa) transformations. Hamiltonian reduction strategies, including Z2 tapering and frozen-core approximations, are examined to assess their effect on quantum resource scaling. The analysis reveals that appropriate transformations, when combined with symmetry-based reductions, can substantially reduce qubit counts by up to ≈50% and quantum gate counts by up to ≈27.5× and Hamiltonian Pauli string counts by up to ≈2.75×, relative to the corresponding unreduced Hamiltonian representations within the same active-space configuration for the representative set of molecular systems under study. These findings provide practical circuit-level insights for executing chemically relevant simulations on NISQ hardware, while establishing physical-resource baselines that may inform future logical-level analyses targeting FASQ systems.
Abstract Adaptive variational quantum algorithms arguably offer the best prospects for quantum advantage in the Noisy Intermediate-Scale Quantum era. Since the inception of the first such algorithm, the Adaptive Derivative-Assembled Problem-Tailored Variational Quantum Eigensolver (ADAPT-VQE), many improvements have appeared in the literature. We combine the key improvements along with a novel operator pool—which we term Coupled Exchange Operator (CEO) pool—to assess the cost of running state-of-the-art ADAPT-VQE on hardware in terms of measurement counts and circuit depth. We show a dramatic reduction of these quantum computational resources compared to the early versions of the algorithm: CNOT count, CNOT depth and measurement costs are reduced by up to 88%, 96% and 99.6%, respectively, for molecules represented by 12 to 14 qubits (LiH, H 6 and BeH 2 ). We also find that our state-of-the-art CEO-ADAPT-VQE outperforms the Unitary Coupled Cluster Singles and Doubles ansatz, the most widely used static VQE ansatz, in all relevant metrics, and offers a five order of magnitude decrease in measurement costs as compared to other static ansätze with competitive CNOT counts.
The Variational Quantum Eigensolver (VQE) and its adaptive variants, such as ADAPT-VQE, are central to the study of strongly correlated quantum systems. However, the classical simulation of the ansatz growth process remains constrained by the exponential scaling of operator space and the associated computational cost of unitary evolution. We introduce the Greedy Operator Commutativity Partitioning (GOCP) framework, an analytical methodology designed to optimize both operator selection and state evolution. By reformulating complex unitary rotations as a chained sequence of fifth-order O(5) Taylor series expansions, GOCP bypasses the need for explicit matrix exponentiation, reducing the computational task to a sequence of sparse matrix-vector operations. We evaluate the performance of this framework across diverse molecular systems, including BeH2 and strongly correlated H2O geometries, utilizing both Jordan-Wigner and Bravyi-Kitaev mappings. Our results demonstrate that the GOCP framework maintains exceptional numerical fidelity-exceeding 1 - 10^-6 in state fidelity-while achieving sub-chemical accuracy in ground-state energy calculations. By enabling the simulation of operator manifolds exceeding 2.68 x 10^8 elements with high efficiency, this approach provides a scalable and rigorous pathway for exploring deep variational circuits in complex quantum many-body systems.