A state-averaged orbital-optimized hybrid quantum-classical algorithm for a democratic description of ground and excited states Saad Yalouz,1, 2, ∗ Bruno Senjean,2, 1, y Jakob G¨unther,1 Francesco Buda,3 Thomas E. O'Brien,4, 2 and Lucas Visscher1 1Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands 2Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 3Leiden Institute of Chemistry, Leiden University, Einsteinweg 55, P.O. Box 9502, 2300 RA Leiden, The Netherlands. 4Google Research, 80636 Munich, Germany In the Noisy Intermediate-Scale Quantum (NISQ) era, solving the electronic structure problem from chemistry is considered as the \killer application" for near-term quantum devices. In spite of the success of variational hybrid quantum/classical algorithms in providing accurate energy pro- files for small molecules, careful considerations are still required for the description of complicated features of potential energy surfaces. Because the current quantum resources are very limited, it is common to focus on a restricted part of the Hilbert space (determined by the set of active orbitals). While physically motivated, this approximation can severely impact the description of these complicated features. A perfect example is that of conical intersections (i.e. a singular point of degeneracy between electronic states), which are of primary importance to understand many prominent reactions. Designing active spaces so that the improved accuracy from a quantum com- puter is not rendered useless is key to finding useful applications of these promising devices within the field of chemistry. To answer this issue, we introduce a NISQ-friendly method called \State- Averaged Orbital-Optimized Variational Quantum Eigensolver" (SA-OO-VQE) which combines two algorithms: (1) a state-averaged orbital-optimizer, and (2) a state-averaged VQE. To demonstrate the success of the method, we classically simulate it on a minimal Schiff base model (namely the formaldimine molecule CH2NH) relevant also for the photoisomerization in rhodopsin | a crucial step in the process of vision mediated by the presence of a conical intersection. We show that merg- ing both algorithms fulfil the necessary condition to describe the molecule's conical intersection, i.e. the ability to treat degenerate (or quasi-degenerate) states on the same footing. I. INTRODUCTION a chemical reaction), several chemical reactions also de- pend on excited states (see for instance Ref. 21 and refer- Quantum computing promises to solve problems that ences therein). To date, most variational quantum algo- are beyond the capacities of classical devices. Technolog- rithms have been specifically designed to prepare ground ical advances in quantum computers occur rapidly [1{8], states [1, 22{37] while excited state calculations have re- and it is now of primary importance to develop quantum ceived less attention until only recently [23, 38{49]. It algorithms dedicated to specific tasks of high industrial should be noted that some excited states can also be ob- and societal impact. The electronic structure problem in tained as the lowest energy solutions of a given set of chemistry is considered as one of the killer applications symmetries [50{55]. of quantum computers [9{14], and it remains an open The landscape of ground- and excited-state potential question whether these applications will be achieved in energy surfaces (PES) can be quite complicated, making the noisy intermediate-scale quantum (NISQ) era [15]. even qualitative numerical description challenging. As While quantum phase estimation algorithms [16{18] can atomic nuclei move, ground- and excited-state energies in principle solve the electronic structure problem in its may tend towards each other, leading to either avoided entirety, they have received limited focus in the NISQ era crossings or true inversions in the ordering of the elec- due to their apparent longer circuit depth requirements. tronic states. These inversions are referred to as conical arXiv:2009.11417v2 [quant-ph] 14 Jan 2021 Instead, variational hybrid quantum/classical algorithms intersections and are ubiquitous in photochemistry and require much shallower circuits at the cost of additional photobiology [56, 57]. Their presence mediates impor- measurements, with as a spearhead the variational quan- tant reactions, for instance the photoisomerisation of the tum eigensolver (VQE) [7, 19, 20]. Those NISQ-friendly chromophore of the visual photoreceptor rhodopsin, en- algorithms are hoped to solve quantum chemical prob- suing absorption of light. Photoisomerisation is a crucial lems of moderate-large size in the near future. step in the process of vision [58{62]. Conical intersections While computing the ground state of a molecule gives also find their importance in photosynthesis [63], photo- information about the equilibrium geometries and tran- stabilization of DNA [64{67] and excitation-energy trans- sition states (essential to determine the energy barrier of fer [68] such as in supramolecular light-harvesting anten- nae [69]. Despite their indisputable importance, their description remains a challenge in the chemistry commu- nity [56, 57, 70]. Diagonalization of the electronic Hamil- ∗ [email protected] tonian matrix being intractable on classical devices, it is y [email protected] 2 of common use to reduce the Hilbert space size by consid- II. THEORY ering only a subspace of (so-called active) orbitals of the full problem. However, such a truncated Hilbert space is A. Electronic structure problem not guaranteed to be accurate enough to represent a con- ical intersection anymore. A possible solution is provided In quantum chemistry, the Born-Oppenheimer approx- by state-averaged orbital-optimization, which consists in imation assumes that molecular structure and dynamics modifying the truncated Hilbert space (without modify- can be treated in two parts: considering first the atomic ing its size) to hopefully recover the conical intersection. nuclei, and then the motion of the electrons around them. On classical computers, this orbital-optimization is part Atomic nuclei, being many times heavier than their elec- of the so-called state-averaged complete active space self- tronic counterparts, are often treated as classical point consistent field (SA-CASSCF) method [71, 72], known for particles with definite position RA and momentum PA its democratic treatment of multiple (possibly degener- (with A the index of the nucleus). As classical parti- ate) eigenstates. cles have no zero-point energy, the ground state of the For quantum computers in the NISQ era, the concept nuclei is motionless, i.e. PA = 0. By contrast, the of active space also represents an important issue. In- light electrons behave as fundamentally quantum parti- deed, the limitations of quantum resources (e.g. max- cles, and form a collective many-body wavefunction Ψ(r), imum number of qubits) constrain algorithms to focus where r = (r ; r ;:::; r ) is the 3N -dimensional global on small-sized active spaces, thus preventing any calcu- 1 2 Ne e position vector of the Ne electrons. This wavefunction lation on large systems. Designing Hilbert spaces so that evolves by the Schr¨odingerequation with the fixed nuclei the improved accuracy from a quantum computer is not providing a background potential field. The quantum op- rendered useless by the active space approximation itself erator describing this dynamics is the so-called electronic is key to finding useful applications of these promising structure Hamiltonian (in atomic units) devices within the field of chemistry. In this context, getting access to complicated chemical features such as Ne Ne Ne Na ^ 1 X 1 X 1 X X ZA avoided crossing or conical intersections can appear as a H = − rri + − ; (1) 2 2 rij riA real challenge for near-term quantum computers. i=1 i6=j i=1 A=1 Motivated by this issue, we introduce in this paper a NISQ-friendly variational hybrid quantum/classical al- where rij = jri − rjj, riA = jri − RAj and ZA are re- gorithm inspired by the SA-CASSCF method, so-called spectively the distance between electrons i and j, the state-averaged orbital-optimized VQE (SA-OO-VQE). distance between electron i and nuclei A, and the atomic The ability of SA-OO-VQE to capture a conical inter- number of atom A. In this context, solving the electronic section is demonstrated on a minimal Schiff base model structure problem means solving the time-independent of the rhodopsin chromophore, i.e. the formaldimine Schr¨odingerequation molecule. ^ The paper is organized as follows. After a brief descrip- H jΨ`i = E` jΨ`i (2) tion of the electronic structure problem in Sec. IIA, the concept of conical intersection is introduced in Sec.IIB. which yields a set of electronic eigenstates jΨ`i ≡ Ψ`(r) A particular attention is paid to the efficiency of classical with corresponding energies E`. These eigenstates and algorithms to represent this characteristic in chemistry energies depend on the atomic positions themselves | (Sec. IIC), as shown on the formaldimine molecule in E` = E`(R) and jΨ`i = jΨ`(R)i, where R is the 3Na- Sec. IID. Then, the SA-OO-VQE hybrid quantum/clas- dimensional position vector of the Na nuclei. (Note the sical algorithm is detailed in Sec. III, where each building change of notation: Ψ`(r) is a single number from evalu- blocks of the algorithm is described separately. Following ating the function Ψ` at the fixed electronic co-ordinates a brief introduction on quantum variational algorithms in r, while jΨ`(R)i is an entire function defined by a fixed Sec.IIIA, we detail how SA-OO-VQE mimics its classical set of nuclear co-ordinates R). In practice, it is conve- analog SA-CASSCF by combining a state-averaged VQE nient to work in a finite basis of No orthonormal molec- (Sec. IIIB1) with a state-averaged orbital-optimization ular orbitals (MO) fφp(r)g. The latter are usually so- procedure (Sec. IIIB2). The importance of choosing a lutions to the mean-field single-particle problem, deter- flexible wavefunction ansatz such as the generalized uni- mined by the Hartree{Fock (HF) method [72].