arXiv:2301.10690 · Journal of Chemical Theory and Computation
Efficient Construction of Involutory Linear Combinations of Anticommuting Pauli Generators for Large-Scale Iterative Qubit Coupled Cluster Calculations
We present an efficient method for construction of a fully anticommutative set of Pauli generators (elements of the Pauli group) from a commutative set of operators that are composed exclusively from Pauli x̂ operators (purely X generators) and sorted by an associated numerical measure, such as absolute energy gradients. Our approach uses the Gauss–Jordan elimination applied to a binary matrix that encodes the set of X generators to bring it to the reduced row-echelon form, followed by the construction of an anticommutative system in a standard basis by means of a modified Jordan–Wigner transformation and returning to the original basis. The algorithm complexity is linear in the size of the X set and quadratic in the number of qubits. The resulting anticommutative sets are used to construct the qubit coupled cluster Ansatz with involutory linear combinations of anticommuting Paulis (QCC-ILCAP) proposed in J. Chem. Theory Comput. 2021, 17 (1), 66–78. We applied the iterative qubit coupled cluster method with the QCC-ILCAP Ansatz to calculations of ground-state potential energy curves for symmetric stretching of the water molecule (36 qubits) and dissociation of N 2 (56 qubits).
We present an efficient method for construction of a fully anticommutative set of Pauli generators (elements of the Pauli group) from a commutative set of operators that are composed exclusively from Pauli x̂ operators (purely X generators) and sorted by an associated numerical measure, such as absolute energy gradients. Our approach uses the Gauss–Jordan elimination applied to a binary matrix that encodes the set of X generators to bring it to the reduced row-echelon form, followed by the construction of an anticommutative system in a standard basis by means of a modified Jordan–Wigner transformation and returning to the original basis. The algorithm complexity is linear in the size of the X set and quadratic in the number of qubits. The resulting anticommutative sets are used to construct the qubit coupled cluster Ansatz with involutory linear combinations of anticommuting Paulis (QCC-ILCAP) proposed in J. Chem. Theory Comput. 2021, 17 (1), 66–78. We applied the iterative qubit coupled cluster method with the QCC-ILCAP Ansatz to calculations of ground-state potential energy curves for symmetric stretching of the water molecule (36 qubits) and dissociation of N 2 (56 qubits).
An iterative version of the qubit coupled cluster (QCC) method [I. G. Ryabinkin et al., J. Chem. Theory Comput. 2019, 14, 6317] is proposed. The new method seeks to find ground electronic energies of molecules on noisy intermediate-scale quantum devices. Each iteration involves a canonical transformation of the Hamiltonian and employs constant-size quantum circuits at the expense of increasing the Hamiltonian size. We numerically studied the convergence of the method on ground-state calculations for LiH, H 2 O, and N 2 molecules and found that the exact ground-state energies can be systematically approached only if the generators of the QCC ansatz are sampled from a specific set of operators. We report an algorithm for constructing this set that scales linearly with the size of the Hamiltonian.
Abstract The iterative qubit coupled cluster (iQCC) method is a systematic variational approach to solve the electronic structure problem on universal quantum computers. It is able to use arbitrarily shallow quantum circuits at expense of iterative canonical transformation of the Hamiltonian and rebuilding a circuit. Here we present a variety of a posteriori corrections to the iQCC energies to reduce the number of iterations to achieve the desired accuracy. Our energy corrections are based on a low-order perturbation theory series that can be efficiently evaluated on a classical computer. Moreover, capturing a part of the total energy perturbatively, allows us to formulate the qubit active-space concept, in which only a subset of all qubits is treated variationally. As a result, further reduction of quantum resource requirements is achieved. We demonstrate the utility and efficiency of our approach numerically on the examples of 10-qubit N 2 molecule dissociation, the 24-qubit H 2 O symmetric stretch, and 56-qubit singlet-triplet gap calculations for the technologically important complex, tris-(2-phenylpyridine)iridium(III) Ir(ppy) 3 .
A unitary coupled cluster (UCC) form for the wave function in the variational quantum eigensolver has been suggested as a systematic way to go beyond the mean-field approximation and include electron correlation in solving quantum chemistry problems on a quantum computer. Although being exact in the limit of including all possible coupled cluster excitations, practically, the accuracy of this approach depends on the number and type of terms are included in the wave function parametrization. Another difficulty of UCC is a growth of the number of simultaneously entangled qubits even at the fixed Fermionic excitation rank. Not all quantum computing architectures can cope with this growth. To address both problems, we introduce a qubit coupled cluster (QCC) method that starts directly in the qubit space and uses energy response estimates for ranking the importance of individual entanglers for the variational energy minimization. Also, we provide an exact factorization of a unitary rotation of more than two qubits to a product of two-qubit unitary rotations. Thus, the QCC method with the factorization technique can be limited to only two-qubit entanglement gates and allows for very efficient use of quantum resources in terms of the number of coupled cluster operators. The method performance is illustrated by calculating ground-state potential energy curves of H 2 and LiH molecules with chemical accuracy, ≤1 kcal/mol, and a symmetric water dissociation curve.
Immense interest in quantum computing has prompted the development of electronic structure methods that are suitable for quantum hardware. However, the slow pace at which quantum hardware progresses forces researchers to implement their ideas on classical computers despite the obvious loss of any “quantum advantage.” As a result, the so-called q uantum-inspired methods emerge. They allow one to look at the electronic structure problem from a different angle, yet to fully exploit their capacity, efficient implementations are highly desirable. Here, we report two schemes for improving the amplitude optimization in the iterative qubit coupled cluster (iQCC) method: a variational quantum eigensolver-type approach, which is based on the qubit coupled cluster (QCC) Ansatz. Our first scheme approximates the QCC unitary as a sum of symmetrical polynomials of generators up to a given order. The resulting energy expression allows for flexible control of computational complexity via the order parameter. It also guarantees smoothness of trial energies and their derivatives, which is important for gradient-based optimization strategies. The second scheme limits the size of the expansion space in which the QCC unitary is generated. It provides better control of memory requirements but, in general, may lead to the nonsmooth variation of energy estimates upon changes in amplitudes. It can be used, however, to extrapolate energies for a given set of amplitudes toward the exact QCC value. Both schemes allow for a larger number of generators to be included in the QCC form compared to the exact formulation. This reduces the number of iterations in the iQCC method and/or leads to higher accuracy. We assess the capabilities of the new schemes to perform QCC amplitudes optimization for a few molecular systems: dinitrogen (N 2, 16 qubits), water (H 2 O, 36 qubits), and tris(2-(2,4-difluorophenyl)pyridine) iridium(III), (Ir(F 2 ppy) 3, 80 qubits).
Recently (Physica Scripta, 100(10):105401, 2025), an algorithm was introduced that deterministically generates a Clifford transformation from the Qubit Coupled Cluster (QCC) algorithm which we call Q-Cliff (QCC+Clifford). There, it was shown that Q-Cliff could be utilized to generate a hardware efficient version of the QCC ansatz. Here, we examine and refine these techniques and show that Q-Cliff can be utilized to generate efficient classical and quantum approximations to the ground states of chemical systems. The algorithm generates an efficient variational method that generally has accuracy between MP2 and CISD with O(N6). Furthermore, we show through DMRG calculations that the entanglement between qubits is reduced significantly and therefore the accuracy for a given bond dimension can be vastly improved (up to an order of magnitude). Finally, we refine the previously reported algorithm to generate low-depth and CNOT efficient circuits that can be optimized with a comparable number of energy evaluations to state-of-the-art VQE algorithms. All these results show that this Hamiltonian derived Clifford transformation should be a tool used for many classical and quantum algorithms.
We propose an efficient method for Hamiltonian simulation of multi-qubit quantum systems with special types of interaction. In our approach, the Hamiltonian of a n-qubit system should be represented as a linear combination of the standard Pauli basis operators, and then decomposed into a sum of partial Hamiltonians, which are, in general, not Pauli operators and satisfy some anticommutation relations. For three types of Hamiltonians, which are invariant with respect to permutations of qubits, the effectiveness of the main algorithm in the three-qubit cluster model is shown by calculating the operator exponentials for these Hamiltonians in an explicit analytical form. We also calculate the density operator, partition function, entropy, and free energy of the cluster weakly coupled to a thermal environment. In our model, the cluster is in the Gibbs state in the temperature interval 0.1−2K, which corresponds to the operating range of modern quantum processors. It follows from our analysis that the thermodynamic properties of such systems strongly depend on the type of internal interaction of qubits in the cluster.