A unitary coupled-cluster (UCC) form for the wavefunction 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 how many and what kind of terms are included in the wavefunction 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.

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.

We investigate the role of the geometric phase (GP) in an internal conversion process when the system changes its electronic state by passing through a conical intersection (CI). Local analysis of a two-dimensional linear vibronic coupling (LVC) model Hamiltonian near the CI shows that the role of the GP is twofold. First, it compensates for a repulsion created by the so-called diagonal Born-Oppenheimer correction (DBOC). Second, the GP enhances the non-adiabatic transition probability for a wave-packet part that experiences a central collision with the CI. To assess the significance of both GP contributions we propose two indicators that can be computed from parameters of electronic surfaces and initial conditions. To generalize our analysis to N -dimensional systems we introduce a reduction of a general N -dimensional LVC model to an effective 2D LVC model using a mode transformation that preserves short-time dynamics of the original N -dimensional model. Using examples of the bis(methylene) adamantyl and butatriene cations, and the pyrazine molecule we have demonstrated that their effective 2D models reproduce the short-time dynamics of the corresponding full dimensional models, and the introduced indicators are very reliable in assessing GP effects.