We introduce several technical and analytical extensions to our recent state-averaged orbital-optimized variational quantum eigensolver (SA-OO-VQE) algorithm (see Yalouz et al. Quantum Sci. Technol. .Show moreWe introduce several technical and analytical extensions to our recent state-averaged orbital-optimized variational quantum eigensolver (SA-OO-VQE) algorithm (see Yalouz et al. Quantum Sci. Technol. 2021, 6, 024004). Motivated by the limitations of current quantum computers, the first extension consists of an efficient state-resolution procedure to find the SA-OO-VQE eigenstates, and not just the subspace spanned by them, while remaining in the equi-ensemble framework. This approach avoids expensive intermediate resolutions of the eigenstates by postponing this problem to the very end of the full algorithm. The second extension allows for the estimation of analytical gradients and nonadiabatic couplings, which are crucial in many practical situations ranging from the search of conical intersections to the simulation of quantum dynamics, in, for example, photoisomerization reactions. The accuracy of our new implementations is demonstrated on the formaldimine molecule CH2NH (a minimal Schiff base model relevant for the study of photoisomerization in larger biomolecules), for which we also perform a geometry optimization to locate a conical intersection between the ground and first-excited electronic states of the molecule. Show less
yalouz, S.; Senjean, B.; Miatto, F.; Dunjko, V. 2021
Variational quantum algorithms (VQA) are considered as some of the most promising methods to determine the properties of complex strongly correlated quantum many-body systems, especially from the... Show moreVariational quantum algorithms (VQA) are considered as some of the most promising methods to determine the properties of complex strongly correlated quantum many-body systems, especially from the perspective of devices available in the near term. In this context, the development of efficient quantum circuit ansatze to encode a many-body wavefunction is one of the keys for the success of a VQA. Great efforts have been invested to study the potential of current quantum devices to encode the eigenstates of fermionic systems, but little is known about the encoding of bosonic systems. In this work, we investigate the encoding of the ground state of the (simple but rich) attractive Bose-Hubbard model using a Continuous-Variable (CV) photonic-based quantum circuit. We introduce two different ansatz architectures and demonstrate that the proposed continuous variable quantum circuits can accurately encode (with a fidelity higher than 99%) the strongly correlated many-boson wavefunction with just a few layers, in all many-body regimes and for different number of bosons and initial states. Beyond the study of the suitability of the ansatz to approximate the ground states of many-boson systems, we also perform initial evaluations of the use of the ansatz in a variational quantum eigensolver algorithm to find it through energy minimization. To this end we also introduce a scheme to measure the Hamiltonian energy in an experimental system , and study the effect of sampling noise. Show less
Senjean, B.; Sen, S.; Repisky, M.; Knizia, G.; Visscher, L. 2021
Two (so‐called left and right) variants of N‐centered ensemble density‐functional theory (DFT) are presented. Unlike the original formulation of the theory, these variants allow for the description... Show moreTwo (so‐called left and right) variants of N‐centered ensemble density‐functional theory (DFT) are presented. Unlike the original formulation of the theory, these variants allow for the description of systems with a fractional electron number. While conventional DFT for open systems uses only the true electron density as basic variable, left/right N‐centered ensemble DFT relies instead on (a) a fictitious ensemble density that integrates to a central (integral) number N of electrons, and (b) a grand canonical ensemble weight α which is equal to the deviation of the true electron number from N. Within such a formalism, the infamous derivative discontinuity that appears when crossing an integral number of electrons is described exactly through the dependence in α of the left and right N‐centered ensemble Hartree‐exchange‐correlation density functionals. Incorporating N‐centered ensembles into existing density‐functional embedding theories is expected to pave the way toward the in‐principle‐exact description of an open fragment by means of a pure‐state N‐electron many‐body wavefunction. Work is currently in progress in this direction. Show less
Modeling chemical reactions and complicated molecular systems has been proposed as the “killer application” of a future quantum computer. Accurate calculations of derivatives of molecular... Show moreModeling chemical reactions and complicated molecular systems has been proposed as the “killer application” of a future quantum computer. Accurate calculations of derivatives of molecular eigenenergies are essential toward this end, allowing for geometry optimization, transition state searches, predictions of the response to an applied electric or magnetic field, and molecular dynamics simulations. In this work, we survey methods to calculate energy derivatives, and present two new methods: one based on quantum phase estimation, the other on a low-order response approximation. We calculate asymptotic error bounds and approximate computational scalings for the methods presented. Implementing these methods, we perform geometry optimization on an experimental quantum processor, estimating the equilibrium bond length of the dihydrogen molecule to within 0.0140.014 Å of the full configuration interaction value. Within the same experiment, we estimate the polarizability of the H22 molecule, finding agreement at the equilibrium bond length to within 0.060.06 a.u. (2%2% relative error). Show less