r""" A brief overview of VQE ======================= .. meta:: :property="og:description": Find the ground state of a Hamiltonian using the variational quantum eigensolver algorithm. :property="og:image": https://pennylane.ai/qml/_images/pes_h2.png .. related:: tutorial_quantum_chemistry Building molecular Hamiltonians tutorial_vqe_parallel VQE with parallel QPUs tutorial_vqe_qng Accelerating VQE with the QNG tutorial_vqe_spin_sectors VQE in different spin sectors tutorial_vqt Variational quantum thermalizer *Author: Alain Delgado. Last updated: 25 June 2022.* The Variational Quantum Eigensolver (VQE) is a flagship algorithm for quantum chemistry using near-term quantum computers [#peruzzo2014]_. It is an application of the `Ritz variational principle `_, where a quantum computer is trained to prepare the ground state of a given molecule. The inputs to the VQE algorithm are a molecular Hamiltonian and a parametrized circuit preparing the quantum state of the molecule. Within VQE, the cost function is defined as the expectation value of the Hamiltonian computed in the trial state. The ground state of the target Hamiltonian is obtained by performing an iterative minimization of the cost function. The optimization is carried out by a classical optimizer which leverages a quantum computer to evaluate the cost function and calculate its gradient at each optimization step. In this tutorial you will learn how to implement the VQE algorithm in a few lines of code. As an illustrative example, we use it to find the ground state of the hydrogen molecule, :math:`\mathrm{H}_2`. First, we build the molecular Hamiltonian using a minimal basis set approximation. Next, we design the quantum circuit preparing the trial state of the molecule, and the cost function to evaluate the expectation value of the Hamiltonian. Finally, we select a classical optimizer, initialize the circuit parameters, and run the VQE algorithm using a PennyLane simulator. Let's get started! Building the electronic Hamiltonian ----------------------------------- The first step is to specify the molecule we want to simulate. This is done by providing a list with the symbols of the constituent atoms and a one-dimensional array with the corresponding nuclear coordinates in `atomic units `_. """ ############################################################################## # The molecular structure can also be imported from an external file using # the :func:`~.pennylane.qchem.read_structure` function. # # Now, we can build the electronic Hamiltonian of the hydrogen molecule # using the :func:`~.pennylane.qchem.molecular_hamiltonian` function. ############################################################################## # The outputs of the function are the Hamiltonian, represented as # a linear combination of Pauli operators, and the number of qubits # required for the quantum simulations. For this example, we use a # `minimal basis set `_ # to represent the `molecular orbitals `_. # In this approximation, we have four spin orbitals, which defines the # number of qubits. Furthermore, we use the Jordan-Wigner # transformation [#seeley2012]_ to perform the fermionic-to-qubit mapping of # the Hamiltonian. # # For a more comprehensive discussion on how to build the Hamiltonian of more # complicated molecules, see the tutorial :doc:`tutorial_quantum_chemistry`. # # Implementing the VQE algorithm # ------------------------------ # From here on, we can use PennyLane as usual, employing its entire stack of # algorithms and optimizers. We begin by defining the device, in this case PennyLane’s # standard qubit simulator: ############################################################################## # Next, we need to define the quantum circuit that prepares the trial state of the # molecule. We want to prepare states of the form, # # .. math:: # \vert \Psi(\theta) \rangle = \cos(\theta/2)~|1100\rangle -\sin(\theta/2)~|0011\rangle, # # where :math:`\theta` is the variational parameter to be optimized in order to find # the best approximation to the true ground state. In the Jordan-Wigner [#seeley2012]_ encoding, # the first term :math:`|1100\rangle` represents the `Hartree-Fock (HF) state # `_ where the two electrons in # the molecule occupy the lowest-energy orbitals. The second term :math:`|0011\rangle` # encodes a double excitation of the HF state where the two particles are excited from # qubits 0, 1 to 2, 3. # # The quantum circuit to prepare the trial state :math:`\vert \Psi(\theta) \rangle` is # schematically illustrated in the figure below. # # | # # .. figure:: /demonstrations/variational_quantum_eigensolver/sketch_circuit.png # :width: 50% # :align: center # # | # # In this figure, the gate :math:`G^{(2)}` corresponds to the # :class:`~.pennylane.DoubleExcitation` operation, implemented in PennyLane # as a `Givens rotation `_, which couples # the four-qubit states :math:`\vert 1100 \rangle` and :math:`\vert 0011 \rangle`. # For more details on how to use the excitation operations to build # quantum circuits for quantum chemistry applications see the # tutorial :doc:`tutorial_givens_rotations`. # # Implementing the circuit above using PennyLane is straightforward. First, we use the # :func:`hf_state` function to generate the vector representing the Hartree-Fock state. ############################################################################## # The ``hf`` array is used by the :class:`~.pennylane.BasisState` operation to initialize # the qubit register. Then, we just act with the :class:`~.pennylane.DoubleExcitation` operation # on the four qubits. ############################################################################## # The next step is to define the cost function to compute the expectation value # of the molecular Hamiltonian in the trial state prepared by the circuit. # We do this using the :func:`~.expval` function. The decorator syntax allows us to # run the cost function as an executable QNode with the gate parameter :math:`\theta`: ############################################################################## # Now we proceed to minimize the cost function to find the ground state of # the :math:`\mathrm{H}_2` molecule. To start, we need to define the classical optimizer. # PennyLane offers many different built-in # `optimizers `_. # Here we use a basic gradient-descent optimizer. ############################################################################## # We initialize the circuit parameter :math:`\theta` to zero, meaning that we start # from the Hartree-Fock state. ############################################################################## # We carry out the optimization over a maximum of 100 steps aiming to reach a # convergence tolerance of :math:`10^{-6}` for the value of the cost function. # store the values of the cost function # store the values of the circuit parameter ############################################################################## # Let's plot the values of the ground state energy of the molecule # and the gate parameter :math:`\theta` as a function of the optimization step. # Full configuration interaction (FCI) energy computed classically # Add energy plot on column 1 # Add angle plot on column 2 ############################################################################## # In this case, the VQE algorithm converges after thirteen iterations. The optimal # value of the circuit parameter :math:`\theta^* = 0.208` defines the state # # .. math:: # \vert \Psi(\theta^*) \rangle = 0.994~\vert 1100 \rangle - 0.104~\vert 0011 \rangle, # # which is precisely the ground state of the :math:`\mathrm{H}_2` molecule in a # minimal basis set approximation. # # Conclusion # ---------- # In this tutorial, we have implemented the VQE algorithm to find the ground state # of the hydrogen molecule. We used a simple circuit to prepare quantum states of # the molecule beyond the Hartree-Fock approximation. The ground-state energy # was obtained by minimizing a cost function defined as the expectation value of the # molecular Hamiltonian in the trial state. # # The VQE algorithm can be used to simulate other chemical phenomena. # In the tutorial :doc:`tutorial_vqe_bond_dissociation`, we use VQE to explore the # potential energy surface of molecules to simulate chemical reactions. # Another interesting application is to probe the lowest-lying states of molecules # in specific sectors of the Hilbert space. For example, see the tutorial # :doc:`tutorial_vqe_spin_sectors`. Furthermore, the algorithm presented here can be # generalized to find the equilibrium geometry of a molecule as it is demonstrated in the # tutorial :doc:`tutorial_mol_geo_opt`. # # .. _vqe_references: # # References # ---------- # # .. [#peruzzo2014] # # Alberto Peruzzo, Jarrod McClean *et al.*, "A variational eigenvalue solver on a photonic # quantum processor". `Nature Communications 5, 4213 (2014). # `__ # # .. [#seeley2012] # # Jacob T. Seeley, Martin J. Richard, Peter J. Love. "The Bravyi-Kitaev transformation for # quantum computation of electronic structure". `Journal of Chemical Physics 137, 224109 (2012). # `__ # # About the author # ---------------- ############################################################################## # .. bio:: Alain Delgado # :photo: ../_static/alain-delgado.png # # Alain Delgado Gran is a senior quantum scientist working at Xanadu. He has a background in quantum many-body systems and electronic structure methods. His work focuses on the application of quantum algorithms for simulating advanced materials.