Source-linked AI summary
Low rank representations for quantum simulation of electronic structure
Mario Motta, Erika Ye, Jarrod R. McClean, Zhendong Li, Austin J. Minnich, Ryan Babbush, Garnet Kin-Lic Chan
TL;DR
The paper develops low-rank representations and approximate basis transformations for electronic-structure operators. Retaining O(N^2) uCC eigenvalues and using partial basis rotations requiring on the order of ρℓN rotations provides a reduced representation, including applications to iron-sulfur clusters.
Problem
The antisymmetry properties of uCC amplitudes lead to less favorable retained-eigenvalue scaling than for the Hamiltonian.
Method
The approach uses eigenvalue truncation, an extended-threshold procedure, and partial basis rotations based on QR decompositions of reduced matrices.
Results
O(N^2) eigenvalues are retained for the uCC operator as system size increases, while approximate basis transformation requires on the order of ρℓN rotations.
Takeaways & Limitations
The reduced-basis procedure supports linear-depth implementation of the relevant rotations on devices with linear connectivity and was applied to iron-sulfur clusters.
Abstract
from arXiv · showhide
The quantum simulation of quantum chemistry is a promising application of quantum computers. However, for N molecular orbitals, the $\mathcal{O}(N^4)$ gate complexity of performing Hamiltonian and unitary Coupled Cluster Trotter steps makes simulation based on such primitives challenging. We substantially reduce the gate complexity of such primitives through a two-step low-rank factorization of the Hamiltonian and cluster operator, accompanied by truncation of small terms. Using truncations that incur errors below chemical accuracy, we are able to perform Trotter steps of the arbitrary basis electronic structure Hamiltonian with $\mathcal{O}(N^3)$ gate complexity in small simulations, which reduces to $\mathcal{O}(N^2 \log N)$ gate complexity in the asymptotic regime, while our unitary Coupled Cluster Trotter step has $\mathcal{O}(N^3)$ gate complexity as a function of increasing basis size for a given molecule. In the case of the Hamiltonian Trotter step, these circuits have $\mathcal{O}(N^2)$ depth on a linearly connected array, an improvement over the $\mathcal{O}(N^3)$ scaling assuming no truncation. As a practical example, we show that a chemically accurate Hamiltonian Trotter step for a 50 qubit molecular simulation can be carried out in the molecular orbital basis with as few as 4,000 layers of parallel nearest-neighbor two-qubit gates, consisting of fewer than 100,000 non-Clifford rotations. We also apply our algorithm to iron-sulfur clusters relevant for elucidating the mode of action of metalloenzymes.
Low-rank factorization of uCC amplitudes
The uCC cluster operator is factorized into a low-rank form using amplitude symmetries, normal-operator squares, and truncated singular values. The retained rank scales less favorably with molecular size than the Hamiltonian rank.
- Low-rank factorization of uCC amplitudes: The antisymmetric uCC amplitudes are represented through a sparse tensor and singular value decomposition after exploiting their mixed symmetries.The amplitudes satisfy tabij = −tbaij = −tabji, and the sparse tensor is decomposed by SVD.
- Low-rank factorization of uCC amplitudes: O(N^2) retained eigenvalues are required as molecular size increases, making uCC scaling less favorable than Hamiltonian scaling.The less favorable scaling is attributed to the antisymmetry properties of the uCC amplitudes.
- Low-rank factorization of uCC amplitudes: The resulting cluster operator is expressed as a linear combination of squares of normal operators.The construction rewrites iτ as i times a linear combination of squared normal operators suitable for implementation.
- Low-rank factorization of uCC amplitudes: Small eigenvalues and singular values are discarded to approximate the uCC operator with a reduced retained rank.Eigenvalues with |σℓ| < ε are discarded, while the largest singular values ρℓ,µ are retained in the subsequent decomposition.
Partial Basis Rotation
Partial basis rotations retain only the leading eigenvectors of each factorized matrix, reducing the number of Givens rotations. The resulting rotations remain compatible with linear-connectivity implementations.
- Partial Basis Rotation: Givens rotations can be ordered to achieve linear depth on a linearly connected device.The algorithm combines partial rotations with consecutive single-particle basis changes between Cholesky-vector bases.
- Partial Basis Rotation: O(ρℓN) Givens rotations implement the approximate basis transformation, where ρℓ is the number of retained eigenvalues.Because ρℓ ≤ N, the partial rotation uses fewer rotations than a general full basis transformation.
- Partial Basis Rotation: A partial basis rotation retains the first ρℓ eigenvectors associated with the largest-magnitude eigenvalues of L(ℓ).This truncation reduces the dimensions of the matrices involved in the basis transformation.
- Partial Basis Rotation: The QR decomposition of the retained N × ρℓ matrix requires at most Nρℓ − ρℓ(ρℓ+1)/2 rotations.This count equals the number of entries below the diagonal that must be eliminated.
Details of calculations
The calculations cover small molecules, alkane chains, and iron-sulfur clusters using electronic-structure methods and selected basis sets. The iron-sulfur studies use chemically motivated active spaces and DMRG calculations for experimentally identified ground states.
- Details of calculations: The study evaluates nine small molecules across STO-6G, 6-31G*, cc-pVDZ, and cc-pVTZ basis sets.These systems are treated with RHF and RCCSD calculations using PySCF.
- Details of calculations: The alkane-chain set contains ethane through octane, represented as CnH2n+2 with n = 2 ... 8.The chains are studied at experimental equilibrium geometries with RHF and RCCSD methods.
- Details of calculations: The Fe-S set includes [2Fe-2S], [4Fe-4S], and the nitrogenase PN cluster.The active spaces include Fe 3d, S 3p, and ligand σ-bonds, forming minimal chemically meaningful active spaces.
- Details of calculations: DMRG calculations target S = 0 states with bond dimensions 8000, 4000, and 2000 for the three Fe-S clusters.The S = 0 states are identified as the experimentally observed ground states.
- Details of calculations: Broken-symmetry UHF calculations use MS = 0 for [2Fe-2S] and [4Fe-4S], while high-spin UHF is used for [8Fe-7S] because of convergence issues.The choice of high-spin UHF is specific to the [8Fe-7S] calculation.