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 N2 (56 qubits).

On the Importance of Well-Defined Thermal Correlation Functions in Simulating Vibronic Spectra.

Two difficulties associated with the computations of thermal vibrational correlation functions are discussed. The first one is the lack of a well-behaved expression that is valid at both high-temperature and T → 0 K limits. Specifically, if the partition function and the propagator are considered separately, then thermal vibrational correlation functions may have an indeterminate form 0/0 in the limit T → 0 K. This difficulty is resolved when the partition function and the propagator are jointly considered in the harmonic approximation, which allows a problematic term that emanates from the zero-point energy to be canceled out, thereby producing a thermal correlation function with a determinate form in the T → 0 K limit. The second difficulty is related to the multivaluedness of the vibrational correlation function. We show numerically that an improper selection of branch leads to discontinuities in the computed correlation function and incorrect vibronic spectra. We propose a phase tracking procedure that ensures continuity of both real and imaginary parts of the correlation function to recover the correct spectra. We support our findings by simulating the UV-vis absorption spectra of pentacene at 4 K and benzene at 298 K. Both are found to be in good agreement with their experimental counterparts.

Estimating Phosphorescent Emission Energies in IrIII Complexes Using Large-Scale Quantum Computing Simulations.

Here we calculate T1 →S0 transition energies in nine phosphorescent iridium complexes using the iterative qubit coupled cluster (iQCC) method to determine if quantum simulations have any advantages over classical methods. These simulations would require a gate-based quantum computer with at least 72 fully-connected logical qubits. Since such devices do not yet exist, we demonstrate the iQCC method using a purpose-built quantum simulator on classical hardware. The results are compared to a selection of common DFT functionals, ab initio methods, and empirical data. iQCC is found to match the accuracy of the best DFT functionals, but with a better correlation coefficient, demonstrating that it is better at predicting the structure-property relationship. Results indicate that the iQCC method has the required accuracy to design organometallic complexes when deployed on emerging quantum hardware and sets an industrially relevant target for demonstrating quantum advantage.

Iterative Qubit Coupled Cluster Approach with Efficient Screening of Generators.

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, H2O, and N2 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.

Constrained Variational Quantum Eigensolver: Quantum Computer Search Engine in the Fock Space.

Variational quantum eigensolver (VQE) is an efficient computational method promising chemical accuracy in electronic structure calculations on a universal-gate quantum computer. However, such a simple task as computing the electronic energy of a hydrogen molecular cation, H2+, is not possible for a general VQE protocol because the calculation will invariably collapse to a lower energy of the corresponding neutral form, H2. The origin of the problem is that VQE effectively performs an unconstrained energy optimization in the Fock space of the original electronic problem. We show how this can be avoided by introducing necessary constraints directing VQE toward the electronic state of interest. The proposed constrained VQE can find an electronic state with a certain number of electrons, a certain spin, or any other property. Moreover, the new algorithm naturally removes unphysical kinks in potential energy surfaces (PESs), which frequently appeared in the regular VQE and required significant additional quantum resources for their removal. We demonstrate the performance of the constrained VQE by simulating PESs of various states of H2 and H2O on Rigetti Computing Inc.'s 19Q-Acorn quantum processor.

Relation between fermionic and qubit mean fields in the electronic structure problem.

For quantum computing applications, the electronic Hamiltonian for the electronic structure problem needs to be unitarily transformed into a qubit form. We found that mean-field procedures on the original electronic Hamiltonian and on its transformed qubit counterpart can give different results. We establish conditions of when fermionic and qubit mean fields provide the same or different energies. In cases when the fermionic mean-field (Hartree-Fock) approach provides an accurate description (electronic correlation effects are small), the choice of molecular orbitals for the electron Hamiltonian representation becomes the determining factor in whether the qubit mean-field energy will be equal to or higher than that of the fermionic counterpart. In strongly correlated cases, the qubit mean-field approach has a higher chance to undergo symmetry breaking and lower its energy below the fermionic counterpart.

Qubit Coupled Cluster Method: A Systematic Approach to Quantum Chemistry on a Quantum Computer.

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 H2 and LiH molecules with chemical accuracy, ≤1 kcal/mol, and a symmetric water dissociation curve.

Design of algal film photobioreactors: material surface energy effects on algal film productivity, colonization and lipid content.

A parallel plate air lift reactor was used to examine the growth kinetics of mixed culture algal biofilms grown on various materials (acrylic, glass, polycarbonate, polystyrene and cellulose acetate). The growth kinetics of the algal biofilms were non-linear overall and their overall productivities ranged from 1.10-2.08g/m(2)day, with those grown on cellulose acetate having the highest productivity. Overall algal biofilm productivity was largely explained by differences in the colonization time which in turn was strongly correlated to the polar surface energy of the material, but weakly correlated to water-material contact angle. When colonization time was taken into account, the productivity for all materials except acrylic was not significantly different at approximately 2g/m(2)/day. Lipid content of the algal biofilms ranged from 6% to 8% (w/w) and was not correlated to water-material contact angle or polar surface energy. The results have potential application for selecting appropriate materials for algal film photobioreactors.