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:2304.13772 · Journal of Chemical Theory and Computation

Block-Invariant Symmetry Shift: Preprocessing Technique for Second-Quantized Hamiltonians to Improve Their Decompositions to Linear Combination of Unitaries

Ignacio Loaiza, Artur F. Izmaylov

Published Nov 8, 2023• 24 citations
arXivHTMLPDF

Abstract

Computational cost of energy estimation for molecular electronic Hamiltonians via quantum phase estimation (QPE) grows with the difference between the largest and smallest eigenvalues of the Hamiltonian. In this work, we propose a preprocessing procedure that reduces the norm of the Hamiltonian without changing its eigenspectrum for the target states of a particular symmetry. The new procedure, block-invariant symmetry shift (BLISS), builds an operator T̂ such that the cost of implementing H ^ − T ^ is reduced compared to that of Ĥ, yet H ^ − T ^ acts on the subspaces of interest the same way as Ĥ does. BLISS performance is demonstrated for a linear combination of unitaries (LCU)-based QPE approaches on a set of small molecules. Using the number of electrons as the symmetry specifying the target set of states, BLISS provided a factor of 2 reduction of 1-norm for several LCU decompositions compared to their unshifted versions.

DOI
10.1021/acs.jctc.3c00912

Related papers

Efficient Block Encodings of Vibrational Hamiltonians

Shashank G. Mehendale, Ignacio Loaiza, Shreyas Malpathak +3 more

Simulation of molecular vibrations is emerging as a promising application of quantum computing in quantum chemistry. In this work, we develop new strategies for block-encoding the second-quantized $n$-mode representation of the vibrational Hamiltonian, whose algebraic structure closely parallels that of the second-quantized electronic Hamiltonian. We derive linear-combination-of-unitaries (LCU) representations of this Hamiltonian based on Pauli decomposition and tensor hypercontraction (THC) and construct efficient circuits for the corresponding block encodings. We also extend the Block-Invariant Symmetry Shift (BLISS) technique to vibrational Hamiltonians, yielding a lower LCU 1-norm and, consequently, a lower overall circuit cost. For the molecules considered here, BLISS-THC reduces the Toffoli cost of quantum phase estimation by more than two orders of magnitude relative to existing block-encoding approaches. In comparison with recent Trotter-based methods, it shows up to a ninefold reduction in the T-gate count, with the relative savings tending to increase with system size. Asymptotically, and up to polylogarithmic factors, both the Pauli and THC LCUs trade additional ancillas for a quadratic improvement in scaling with the number of modes, while THC provides an additional quadratic improvement in scaling with the number of modals. These results highlight the importance of Hamiltonian factorizations and Toffoli-qubit tradeoffs for reducing the computational cost of quantum simulations of molecular vibrations.

quant-pharXiv:2610.33800•Oct 27, 2026•0 citations

Global Minimization of Electronic Hamiltonian 1-Norm via Linear Programming in the Block Invariant Symmetry Shift (BLISS) Method

Smik Patel, Aritra Sankar Brahmachari, Joshua T. Cantin +2 more

The cost of encoding a system Hamiltonian in a digital quantum computer as a linear combination of unitaries (LCU) grows with the 1-norm of the LCU expansion. The Block Invariant Symmetry Shift (BLISS) technique reduces this 1-norm by modifying the Hamiltonian action on only the undesired electron-number subspaces. Previously, BLISS required a computationally expensive nonlinear optimization that was not guaranteed to find the global minimum. Here, we introduce various reformulations of this optimization as a linear programming problem, which guarantees optimality and significantly reduces the computational cost. We apply BLISS to industrially relevant homogeneous catalysts in active spaces of up to 76 orbitals, finding substantial reductions in both the spectral range of the modified Hamiltonian and the 1-norms of Pauli and fermionic LCUs. Our linear programming techniques for obtaining the BLISS operator enable more efficient Hamiltonian simulation and, by reducing the Hamiltonian's spectral range, offer opportunities for improved LCU groupings to further reduce the 1-norm.

arXiv:2409.18277•Journal of Chemical Theory and Computation•Jan 13, 2025•6 citations

Simultaneously Optimizing Symmetry Shifts and Tensor Factorizations for Cost-Efficient Fault-Tolerant Quantum Simulations of Electronic Hamiltonians

Konrad Deka, E. J. Zak

In fault-tolerant quantum computing, the cost of calculating Hamiltonian eigenvalues using the quantum phase estimation algorithm is proportional to the constant scaling the Hamiltonian matrix block-encoded in a unitary circuit. We present a method to reduce this scaling constant for the electronic Hamiltonians represented as a linear combination of unitaries. Our approach combines the double tensor-factorization method of Burg et al. with the block-invariant symmetry shift method of Loaiza and Izmaylov. By extending the electronic Hamiltonian with appropriately parametrized symmetry operators and optimizing the tensor factorization parameters, our method achieves a 25% reduction in the block-encoding scaling constant compared to previous best techniques for industrially and biologically relevant molecules, including the nitrogenase cofactor (FeMoCo) and cytochrome P450. Benchmark calculations comparing our method with the original double-factorization method show a 75% reduction in the block-encoding scaling constant. The resulting savings in the number of non-Clifford T-gates, which are an essential resource for fault-tolerant quantum computation, are expected to accelerate the journey to the feasibility of practical Hamiltonian simulations.

arXiv:2412.01338•Journal of Chemical Theory and Computation•Apr 22, 2025•4 citations

Assessing the query complexity limits of quantum phase estimation using symmetry-aware spectral bounds

Cristian L. Cortes, Dario Rocca, Jérôme F. Gonthier +7 more

The computational cost of quantum algorithms for physics and chemistry is closely linked to the spectrum of the Hamiltonian, a property that manifests in the necessary rescaling of its eigenvalues. The typical approach of using the 1-norm as an upper bound to the spectral norm to rescale the Hamiltonian suits the most general case of bounded Hermitian operators but neglects the influence of symmetries commonly found in chemical systems. In this work, we introduce a hierarchy of symmetry-aware spectral bounds that provide a unified understanding of the performance of quantum phase estimation algorithms using block-encoded electronic structure Hamiltonians. We present a variational and numerically tractable method for computing these bounds, based on orbital optimization, to demonstrate that the computed bounds are smaller than conventional spectral bounds for a variety of molecular benchmark systems. We also highlight the unique analytical and numerical scaling behavior of these bounds in the thermodynamic and complete basis set limits. Our work shows that there is room for improvement in reducing the 1-norm, not yet achieved through methods like double factorization and tensor hypercontraction, but highlights the potential challenges in improving the performance of current quantum algorithms beyond small constant factors through 1-norm reduction techniques alone.

arXiv:2403.04737•Physical Review A•Aug 14, 2024•5 citations

Reducing molecular electronic Hamiltonian simulation cost for linear combination of unitaries approaches

Ignacio Loaiza, Alireza Marefat Khah, Nathan Wiebe +1 more

Abstract We consider different linear combination of unitaries (LCU) decompositions for molecular electronic structure Hamiltonians. Using these LCU decompositions for Hamiltonian simulation on a quantum computer, the main figure of merit is the 1-norm of their coefficients, which is associated with the quantum circuit complexity. It is derived that the lowest possible LCU 1-norm for a given Hamiltonian is half of its spectral range. This lowest norm decomposition is practically unattainable for general Hamiltonians; therefore, multiple practical techniques to generate LCU decompositions are proposed and assessed. A technique using symmetries to reduce the 1-norm further is also introduced. In addition to considering LCU in the Schrödinger picture, we extend it to the interaction picture, which substantially further reduces the 1-norm.

arXiv:2208.08272•Quantum Science and Technology•May 15, 2023•46 citations

Algebraic Tools for Quantum Simulations of Electronic Structure

Smik Patel

Accurately predicting the properties of molecules and materials from first principles requires solving the electronic structure problem. However, the computational cost of solving the electronic structure problem scales prohibitively in the number of atoms, and the usual approximate methods designed to reduce the overhead are often inaccurate for strongly correlated systems. For this reason, quantum algorithms have been developed as a promising alternative. In this work, we develop algebraic tools to improve the efficiency of three quantum algorithms for the electronic structure problem, addressing key bottlenecks for both near-term and fault-tolerant quantum computers. The first is the variational quantum eigensolver (VQE), which suffers from a large measurement overhead needed to estimate the Hamiltonian expectation value. To this end, we introduce a general class of measurable fragments of the electronic Hamiltonian, which enable partitioning schemes that significantly reduce the sampling cost for energy estimation. The second are quantum subspace methods (QSMs). We introduce a new QSM, the quantum seniority-based subspace expansion algorithm, which is a scalable alternative to VQE that alleviates VQE's state preparation, measurement, and optimization challenges while reducing qubit requirements. The third is quantum phase estimation (QPE), which requires an encoding of the electronic Hamiltonian directly on the quantum computer, the cost of which scales indirectly with the Hamiltonian spectral range. To address this, we develop linear programming approaches to implement the Block Invariant Symmetry Shift protocol, a method which reduces the spectral range of the Hamiltonian, improving its block-encoding cost for QPE. Collectively, these methods demonstrate that one can exploit the Hamiltonian's algebraic structure to mitigate computational costs, accelerating the path toward practical quantum advantage for electronic structure.

•Mar 1, 2026•0 citations
Computational cost of energy estimation for molecular electronic Hamiltonians via quantum phase estimation (QPE) grows with the difference between the largest and smallest eigenvalues of the Hamiltonian. In this work, we propose a preprocessing procedure that reduces the norm of the Hamiltonian without changing its eigenspectrum for the target states of a particular symmetry. The new procedure, block-invariant symmetry shift (BLISS), builds an operator T̂ such that the cost of implementing H ^ − T ^ is reduced compared to that of Ĥ, yet H ^ − T ^ acts on the subspaces of interest the same way as Ĥ does. BLISS performance is demonstrated for a linear combination of unitaries (LCU)-based QPE approaches on a set of small molecules. Using the number of electrons as the symmetry specifying the target set of states, BLISS provided a factor of 2 reduction of 1-norm for several LCU decompositions compared to their unshifted versions.
Simulation of molecular vibrations is emerging as a promising application of quantum computing in quantum chemistry. In this work, we develop new strategies for block-encoding the second-quantized nnn-mode representation of the vibrational Hamiltonian, whose algebraic structure closely parallels that of the second-quantized electronic Hamiltonian. We derive linear-combination-of-unitaries (LCU) representations of this Hamiltonian based on Pauli decomposition and tensor hypercontraction (THC) and construct efficient circuits for the corresponding block encodings. We also extend the Block-Invariant Symmetry Shift (BLISS) technique to vibrational Hamiltonians, yielding a lower LCU 1-norm and, consequently, a lower overall circuit cost. For the molecules considered here, BLISS-THC reduces the Toffoli cost of quantum phase estimation by more than two orders of magnitude relative to existing block-encoding approaches. In comparison with recent Trotter-based methods, it shows up to a ninefold reduction in the T-gate count, with the relative savings tending to increase with system size. Asymptotically, and up to polylogarithmic factors, both the Pauli and THC LCUs trade additional ancillas for a quadratic improvement in scaling with the number of modes, while THC provides an additional quadratic improvement in scaling with the number of modals. These results highlight the importance of Hamiltonian factorizations and Toffoli-qubit tradeoffs for reducing the computational cost of quantum simulations of molecular vibrations.
The cost of encoding a system Hamiltonian in a digital quantum computer as a linear combination of unitaries (LCU) grows with the 1-norm of the LCU expansion. The Block Invariant Symmetry Shift (BLISS) technique reduces this 1-norm by modifying the Hamiltonian action on only the undesired electron-number subspaces. Previously, BLISS required a computationally expensive nonlinear optimization that was not guaranteed to find the global minimum. Here, we introduce various reformulations of this optimization as a linear programming problem, which guarantees optimality and significantly reduces the computational cost. We apply BLISS to industrially relevant homogeneous catalysts in active spaces of up to 76 orbitals, finding substantial reductions in both the spectral range of the modified Hamiltonian and the 1-norms of Pauli and fermionic LCUs. Our linear programming techniques for obtaining the BLISS operator enable more efficient Hamiltonian simulation and, by reducing the Hamiltonian's spectral range, offer opportunities for improved LCU groupings to further reduce the 1-norm.
In fault-tolerant quantum computing, the cost of calculating Hamiltonian eigenvalues using the quantum phase estimation algorithm is proportional to the constant scaling the Hamiltonian matrix block-encoded in a unitary circuit. We present a method to reduce this scaling constant for the electronic Hamiltonians represented as a linear combination of unitaries. Our approach combines the double tensor-factorization method of Burg et al. with the block-invariant symmetry shift method of Loaiza and Izmaylov. By extending the electronic Hamiltonian with appropriately parametrized symmetry operators and optimizing the tensor factorization parameters, our method achieves a 25% reduction in the block-encoding scaling constant compared to previous best techniques for industrially and biologically relevant molecules, including the nitrogenase cofactor (FeMoCo) and cytochrome P450. Benchmark calculations comparing our method with the original double-factorization method show a 75% reduction in the block-encoding scaling constant. The resulting savings in the number of non-Clifford T-gates, which are an essential resource for fault-tolerant quantum computation, are expected to accelerate the journey to the feasibility of practical Hamiltonian simulations.
The computational cost of quantum algorithms for physics and chemistry is closely linked to the spectrum of the Hamiltonian, a property that manifests in the necessary rescaling of its eigenvalues. The typical approach of using the 1-norm as an upper bound to the spectral norm to rescale the Hamiltonian suits the most general case of bounded Hermitian operators but neglects the influence of symmetries commonly found in chemical systems. In this work, we introduce a hierarchy of symmetry-aware spectral bounds that provide a unified understanding of the performance of quantum phase estimation algorithms using block-encoded electronic structure Hamiltonians. We present a variational and numerically tractable method for computing these bounds, based on orbital optimization, to demonstrate that the computed bounds are smaller than conventional spectral bounds for a variety of molecular benchmark systems. We also highlight the unique analytical and numerical scaling behavior of these bounds in the thermodynamic and complete basis set limits. Our work shows that there is room for improvement in reducing the 1-norm, not yet achieved through methods like double factorization and tensor hypercontraction, but highlights the potential challenges in improving the performance of current quantum algorithms beyond small constant factors through 1-norm reduction techniques alone.
Abstract We consider different linear combination of unitaries (LCU) decompositions for molecular electronic structure Hamiltonians. Using these LCU decompositions for Hamiltonian simulation on a quantum computer, the main figure of merit is the 1-norm of their coefficients, which is associated with the quantum circuit complexity. It is derived that the lowest possible LCU 1-norm for a given Hamiltonian is half of its spectral range. This lowest norm decomposition is practically unattainable for general Hamiltonians; therefore, multiple practical techniques to generate LCU decompositions are proposed and assessed. A technique using symmetries to reduce the 1-norm further is also introduced. In addition to considering LCU in the Schrödinger picture, we extend it to the interaction picture, which substantially further reduces the 1-norm.
Accurately predicting the properties of molecules and materials from first principles requires solving the electronic structure problem. However, the computational cost of solving the electronic structure problem scales prohibitively in the number of atoms, and the usual approximate methods designed to reduce the overhead are often inaccurate for strongly correlated systems. For this reason, quantum algorithms have been developed as a promising alternative. In this work, we develop algebraic tools to improve the efficiency of three quantum algorithms for the electronic structure problem, addressing key bottlenecks for both near-term and fault-tolerant quantum computers. The first is the variational quantum eigensolver (VQE), which suffers from a large measurement overhead needed to estimate the Hamiltonian expectation value. To this end, we introduce a general class of measurable fragments of the electronic Hamiltonian, which enable partitioning schemes that significantly reduce the sampling cost for energy estimation. The second are quantum subspace methods (QSMs). We introduce a new QSM, the quantum seniority-based subspace expansion algorithm, which is a scalable alternative to VQE that alleviates VQE's state preparation, measurement, and optimization challenges while reducing qubit requirements. The third is quantum phase estimation (QPE), which requires an encoding of the electronic Hamiltonian directly on the quantum computer, the cost of which scales indirectly with the Hamiltonian spectral range. To address this, we develop linear programming approaches to implement the Block Invariant Symmetry Shift protocol, a method which reduces the spectral range of the Hamiltonian, improving its block-encoding cost for QPE. Collectively, these methods demonstrate that one can exploit the Hamiltonian's algebraic structure to mitigate computational costs, accelerating the path toward practical quantum advantage for electronic structure.