reference/approximate/quantum_chemistry.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<img src=\"../../images/QISKit.gif\" alt=\"Note: In order for images to show up in this jupyter notebook you need to select File => Trusted Notebook\" width=\"250 px\" align=\"left\">"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## _*VQE algorithm: Application to quantum chemistry*_ \n",
    "\n",
    "The latest version of this notebook is available on https://github.com/QISKit/qiskit-tutorial.\n",
    "\n",
    "***\n",
    "### Contributors\n",
    "Antonio Mezzacapo, Jay Gambetta"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Introduction\n",
    "\n",
    "One of the most compelling possibilities of quantum computation is the the simulation of other quantum systems. Quantum simulation of quantum systems encompasses a wide range of tasks, including most significantly:\n",
    "    \n",
    "1. Simulation of the time evolution of quantum systems.\n",
    "\n",
    "2. Computation of ground state properties. \n",
    "\n",
    "These applications are especially useful when considering systems of interacting fermions, such as molecules and strongly correlated materials. The computation of ground state properties of fermionic systems is the starting point for mapping out the phase diagram of condensed matter Hamiltonians. It also gives access to the key question of electronic structure problems in quantum chemistry - namely, reaction rates. The focus of this notebook is on molecular systems, which are considered to be the ideal bench test for early-stage quantum computers, due to their relevance in chemical applications despite relatively modest sizes. Formally, the ground state problem asks the following:\n",
    "\n",
    "For some physical Hamiltonian *H*, find the smallest eigenvalue $E_G$, such that $H|\\psi_G\\rangle=E_G|\\psi_G\\rangle$, where $|\\Psi_G\\rangle$ is the eigenvector corresponding to $E_G$.\n",
    "    \n",
    "It is known that in general this problem is intractable, even on a quantum computer. This means that we cannot expect an efficient quantum algorithm that prepares the ground state of general local Hamiltonians. Despite this limitation, for specific Hamiltonians of interest it might be possible, given physical constraints on the interactions, to solve the above problem efficiently. Currently, at least four different methods exist to approach this problem:\n",
    "\n",
    "1. Quantum phase estimation: Assuming that we can approximately prepare the state $|\\psi_G\\rangle$, this routine uses controlled implementations of the Hamiltonian to find its smallest eigenvalue. \n",
    "2. Adiabatic theorem of quantum mechanics: The quantum system is adiabatically dragged from being the ground state of a trivial Hamiltonian to the one of the target problem, via slow modulation of the Hamiltonian terms. \n",
    "3. Dissipative (non-unitary) quantum operation: The ground state of the target system is a fixed point. The non-trivial assumption here is the implementation of the dissipation map on quantum hardware. \n",
    "4. Variational quantum eigensolvers: Here we assume that the ground state can be represented by a parameterization containing a relatively small number of parameters.\n",
    "\n",
    "In this notebook we focus on the last method, as this is most likely the simplest to be realized on near-term devices. \n",
    "\n",
    "The general idea is to define a parameterization $|\\psi(\\boldsymbol\\theta)\\rangle$ of quantum states, and minimize the energy \n",
    "\n",
    "$$E(\\boldsymbol\\theta) = \\langle \\psi(\\boldsymbol\\theta)| H |\\psi(\\boldsymbol\\theta)\\rangle,$$ \n",
    "\n",
    "The key ansatz is that the number of parameters $|\\boldsymbol\\theta^*|$ that minimizes the energy function scales polynomially with the size (e.g., number of qubits) of the target problem.  \n",
    "\n",
    "Then, any local fermionic Hamiltonian can be mapped into a sum over Pauli operators $P_i$,  \n",
    "\n",
    "$$H\\rightarrow H_P = \\sum_i^M w_i P_i,$$ \n",
    "\n",
    "and the energy corresponding to the state $|\\psi(\\boldsymbol\\theta\\rangle$, $E(\\boldsymbol\\theta)$, can be estimated by sampling the individual Pauli terms $P_i$ (or sets of them that can be measured at the same time) on a quantum computer: \n",
    "\n",
    "$$E(\\boldsymbol\\theta) = \\sum_i^M w_i \\langle \\psi(\\boldsymbol\\theta)| P_i |\\psi(\\boldsymbol\\theta)\\rangle.$$ \n",
    "\n",
    "Last, some optimization technique must be devised in order to find the optimal value of parameters $\\boldsymbol\\theta^*$, such that $|\\psi(\\boldsymbol\\theta^*)\\rangle\\equiv|\\psi_G\\rangle$.\n",
    "\n",
    "\n",
    "\n",
    "### Fermionic Hamiltonians\n",
    "\n",
    "The Hamiltonians describing systems of interacting fermions can be expressed in second quantization language, considering fermionic creation (annihilation) operators $a^\\dagger_\\alpha(a_\\alpha)$, relative to the $\\alpha$-th fermionic mode. In the case of molecules, the $\\alpha$ labels stand for the different atomic or molecular orbitals. Within the second-quantization framework, a generic molecular Hamiltonian with $M$ orbitals can be written as \n",
    "$$H =H_1+H_2=\\sum_{\\alpha, \\beta=0}^{M-1} t_{\\alpha \\beta} \\, a^\\dagger_{\\alpha} a_{\\beta} +\\frac{1}{2}  \\sum_{\\alpha, \\beta, \\gamma, \\delta = 0}^{M-1} u_{\\alpha \\beta \\gamma \\delta}\\, a^\\dagger_{\\alpha} a^\\dagger_{\\gamma} a_{\\delta} a_{\\beta},$$\n",
    "with the one-body terms representing the kinetic energy of the electrons and the potential energy that they experience in the presence of the nuclei, \n",
    "$$ t_{\\alpha\\beta}=\\int d\\boldsymbol x_1\\Psi_\\alpha(\\boldsymbol{x}_1) \\left(-\\frac{\\boldsymbol\\nabla_1^2}{2}+\\sum_{i} \\frac{Z_i}{|\\boldsymbol{r}_{1i}|}\\right)\\Psi_\\beta (\\boldsymbol{x}_1),$$\n",
    "and their interactions via Coulomb forces \n",
    "$$ u_{\\alpha\\beta\\gamma\\delta}=\\int\\int d \\boldsymbol{x}_1 d \\boldsymbol{x}_2 \\Psi_\\alpha^*(\\boldsymbol{x}_1)\\Psi_\\beta(\\boldsymbol{x}_1)\\frac{1}{|\\boldsymbol{r}_{12}|}\\Psi_\\gamma^*(\\boldsymbol{x}_2)\\Psi_\\delta(\\boldsymbol{x}_2),$$\n",
    "where we have defined the nuclei charges $Z_i$, the nuclei-electron and electron-electron separations $\\boldsymbol{r}_{1i}$ and $\\boldsymbol{r}_{12}$, the $\\alpha$-th orbital wavefunction $\\Psi_\\alpha(\\boldsymbol{x}_1)$, and we have assumed that the spin is conserved in the spin-orbital indices $\\alpha,\\beta$ and $\\alpha,\\beta,\\gamma,\\delta$.\n",
    "\n",
    "### Molecules considered in this notebook and mapping to qubits\n",
    "\n",
    "We consider in this notebook the optimization of two potential energy surfaces, for the hydrogen and lithium hydride molecules, obtained using the STO-3G basis. The molecular Hamiltonians are computed as a function of their interatomic distance, then mapped to two-(H$_2$) and four-(LiH$_2$) qubit problems, via elimination of core and high-energy orbitals and removal of $Z_2$ symmetries.\n",
    "\n",
    "### Approximate universal quantum computing for quantum chemisty problems\n",
    "\n",
    "In order to find the optimal parameters $\\boldsymbol\\theta^*$, we set up a closed optimization loop with a quantum computer, based on some stochastic optimization routine. Our choice for the variational ansatz is a deformation of the one used for the optimization of classical combinatorial problems, with the inclusion of $Z$ rotation together with the $Y$ ones. The optimization algorithm for fermionic Hamiltonians is similar to the one for combinatorial problems, and can be summarized as follows: \n",
    "\n",
    "\n",
    "1. Map the fermionic Hamiltonian $H$ to a qubit Hamiltonian $H_P$.\n",
    "2. Choose the maximum depth of the quantum circuit (this could be done adaptively).\n",
    "3. Choose a set of controls $\\boldsymbol\\theta$ and make a trial function $|\\psi(\\boldsymbol\\theta)\\rangle$. The difference with the combinatorial problems is the insertion of additional parametrized $Z$ single-qubit rotations.\n",
    "4. Evaluate the energy $E(\\boldsymbol\\theta) = \\langle\\psi(\\boldsymbol\\theta)~|H_P|~\\psi(\\boldsymbol\\theta)\\rangle$ by sampling each Pauli term individually, or sets of Pauli terms that can be measured in the same tensor product basis.\n",
    "5. Use a classical optimizer to choose a new set of controls.\n",
    "6. Continue until the energy has converged, hopefully close to the real solution $\\boldsymbol\\theta^*$, and return the last value of $E(\\boldsymbol\\theta)$.\n",
    "    \n",
    "\n",
    "Note that, as opposed to the classical case, in the case of a quantum chemistry Hamiltonian one has to sample over non-computational states that are superpositions, and therefore take advantage of using a quantum computer in the sampling part of the algorithm. Motivated by the quantum nature of the answer, we also define a variational trial ansatz in this way: \n",
    "\n",
    "$$|\\psi(\\boldsymbol\\theta)\\rangle  = [U_\\mathrm{single}(\\boldsymbol\\theta) U_\\mathrm{entangler}]^m |+\\rangle$$\n",
    "\n",
    "where $U_\\mathrm{entangler}$ is a collection of cPhase gates (fully entangling gates), $U_\\mathrm{single}(\\boldsymbol\\theta) = \\prod_{i=1}^n Y(\\theta_{i})Z(\\theta_{n+i})$ are single-qubit $Y$ and $Z$ rotation, $n$ is the number of qubits and $m$ is the depth of the quantum circuit. \n",
    "\n",
    "\n",
    "\n",
    "References and additional details:\n",
    "\n",
    "[1] A. Kandala, A. Mezzacapo, K. Temme, M. Takita, M. Brink, J. M. Chow, and J. M. Gambetta, *Hardware-efficient Variational Quantum Eigensolver for Small Molecules and Quantum Magnets*, Nature 549, 242 (2017), and references therein.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# useful additional packages \n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "import numpy as np\n",
    "from scipy import linalg as la\n",
    "from functools import partial\n",
    "\n",
    "# importing the QISKit\n",
    "from qiskit import QuantumProgram\n",
    "import Qconfig\n",
    "\n",
    "# import basic plot tools\n",
    "from qiskit.tools.visualization import plot_histogram\n",
    "\n",
    "\n",
    "# import optimization tools\n",
    "from qiskit.tools.apps.optimization import trial_circuit_ryrz, SPSA_optimization, SPSA_calibration\n",
    "from qiskit.tools.apps.optimization import Hamiltonian_from_file, make_Hamiltonian\n",
    "from qiskit.tools.apps.optimization import eval_hamiltonian, group_paulis\n",
    "\n",
    "# Ignore warnings due to chopping of small imaginary part of the energy \n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')\n",
    "                "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Optimization of H$_2$ at bond length\n",
    "\n",
    "In this first part of the notebook we show the optimization of the H$_2$ Hamiltonian in the STO-3G basis at the bond length of 0.735 Angstrom. After mapping it to a four-qubit system with a binary-tree transformation, two spin-parity symmetries are modded out, leading to a two-qubit Hamiltonian. The energy of the mapped Hamiltonian obtained is then minimized using the variational ansatz described in the introduction, and a stochastic perturbation simultaneous approximation (SPSA) gradient descent method. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The exact ground state energy is:\n",
      "-1.8572746705\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 1.23587343356\n",
      "objective function at theta+ for step # 0\n",
      "-1.12184334801\n",
      "objective function at theta- for step # 0\n",
      "-1.13210224717\n",
      "objective function at theta+ for step # 20\n",
      "-1.35606481784\n",
      "objective function at theta- for step # 20\n",
      "-1.36663521809\n",
      "objective function at theta+ for step # 40\n",
      "-1.8133659706\n",
      "objective function at theta- for step # 40\n",
      "-1.81514145739\n",
      "objective function at theta+ for step # 60\n",
      "-1.85052903197\n",
      "objective function at theta- for step # 60\n",
      "-1.84626633243\n",
      "objective function at theta+ for step # 80\n",
      "-1.85681616062\n",
      "objective function at theta- for step # 80\n",
      "-1.85686130312\n",
      "Final objective function is: -1.85727338303\n"
     ]
    }
   ],
   "source": [
    "n=2\n",
    "m=6\n",
    "device='local_qasm_simulator'\n",
    "\n",
    "initial_theta=np.random.randn(2*n*m)\n",
    "entangler_map = Q_program.get_backend_configuration(device)['coupling_map'] # the map of two-qubit gates with control at key and target at values\n",
    "if entangler_map == 'all-to-all':\n",
    "    entangler_map = {i: [j for j in range(n) if j != i] for i in range(n)}\n",
    "shots=1\n",
    "max_trials=100\n",
    "ham_name='H2/H2Equilibrium.txt\n",
    "#ham_name=\"LIH/LiHEquilibrium.txt\" #For optimization of LiH at bond length\n",
    "\n",
    "# Exact Energy\n",
    "pauli_list=Hamiltonian_from_file(ham_name)\n",
    "H=make_Hamiltonian(pauli_list)\n",
    "exact=np.amin(la.eig(H)[0]).real\n",
    "print('The exact ground state energy is:')\n",
    "print(exact)\n",
    "pauli_list_grouped=group_paulis(pauli_list)\n",
    "\n",
    "# Optimization\n",
    "Q_program = QuantumProgram()\n",
    "Q_program.set_api(Qconfig.APItoken,Qconfig.config[\"url\"])\n",
    "\n",
    "\n",
    "\n",
    "def cost_function(Q_program,H,n,m,entangler_map,shots,device,theta):\n",
    "    \n",
    "    return eval_hamiltonian(Q_program,H,trial_circuit_ryrz(n,m,theta,entangler_map,None,False),shots,device).real\n",
    "\n",
    "\n",
    "initial_c=0.01\n",
    "target_update=2*np.pi*0.1\n",
    "save_step = 20\n",
    "\n",
    "if shots ==1:\n",
    "    SPSA_params=SPSA_calibration(partial(cost_function,Q_program,H,n,m,entangler_map,\n",
    "                                         shots,device),initial_theta,initial_c,target_update,25)\n",
    "    output=SPSA_optimization(partial(cost_function,Q_program,H,n,m,entangler_map,shots,device),\n",
    "                         initial_theta,SPSA_params,max_trials,save_step,1);\n",
    "else:\n",
    "    SPSA_params=SPSA_calibration(partial(cost_function,Q_program,pauli_list_grouped,n,m,entangler_map,\n",
    "                                         shots,device),initial_theta,initial_c,target_update,25)\n",
    "    output=SPSA_optimization(partial(cost_function,Q_program,pauli_list_grouped,n,m,entangler_map,shots,device),\n",
    "                         initial_theta,SPSA_params,max_trials,save_step,1);\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x110e91278>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VFX6x/HPk5l0SEIg1KABQXoPRREVQSlSBKTKCljY\nRXd1XV3FLipWRCmKgtgR6VJVRGFVQCBAgNAFI4QaEgiBEFLm/P6YSX4BQnpyJ8nzfr3ui5k7Z+79\nJhPy5J577zlijEEppZTKKw+rAyillCpdtHAopZTKFy0cSiml8kULh1JKqXzRwqGUUipftHAopZTK\nFy0cSiml8kULh1JKqXzRwqGUUipf7FYHKA5VqlQxYWFhVsdQSqlSY/PmzaeMMSF5aWtJ4RCRgcBL\nQCOgnTEm4irtPgF6ASeNMU3zuv2wsDAiIrLdpFJKqWyIyF95bWtVV1UU0B/4JZd2nwHdiz2NUkqp\nPLPkiMMYsxtARHJr94uIhJVAJKWUUnlUZk6Oi8hoEYkQkYjY2Fir4yilVJlVbEccIrIKqJ7NS88a\nYxYX9f6MMdOB6QDh4eE6VrxSxSw1NZWYmBiSk5OtjqLywcfHh9DQUDw9PQu8jWIrHMaYrsW1baWU\n9WJiYqhYsSJhYWG5djsr92CMIS4ujpiYGOrUqVPg7ZSZriqlVMlKTk6mcuXKWjRKERGhcuXKhT5K\ntKRwiEg/EYkBbgCWi8gPrvU1RWRFlnazgfVAAxGJEZH7rcirlMqeFo3Spyg+M6uuqloELMpm/VGg\nZ5bnQ0sy15Sf9tO5YVWa1gosyd0qpVSpol1VLmeSUpi98RADP1zPdzuOWR1HKaXclhYOlyA/L759\n+EYa1qjImFlbmLRqP8boxVlKuTObzUbLli0zlzfeeCPztbvvvpuDBw8C8Nprr2Wuj46OpmnTPA9E\nAcC3337Lrl27iia0y5o1a+jVq1eB3jtkyBD2799fpHnyQwtHhrQUqi4exty2++nfuhbvrtrHP2dv\n5UJKutXJlFJX4evrS2RkZOYyduxYAHbu3El6ejp169YFLi0cBVEchaMwxowZw1tvvWXZ/svkIIcF\nkn4RAM/lj/JO+zE07DaK11f+wV9x55lxbzg1An0tDqiU+xq3dCe7jp4t0m02rhnAi72bFOi9s2bN\nom/fvgCMHTuWCxcu0LJlS5o0acL48eNJT0/nwQcfZN26ddSqVYvFixfj6+vLgQMHePjhh4mNjcXP\nz48ZM2YQHx/PkiVL+N///serr77KggUL+Pnnn5k+fTopKSnUq1ePL7/8Ej8/v2yzjBw5Eh8fHyIi\nIjh79iwTJ0684kjjpZdeokKFCjzxxBMANG3alGXLlhESEsKgQYOIiYkhPT2d559/nsGDB9OpUydG\njhxJWloadnvJ/xrXI44M3hVh2DxoPwbZMI3RMU/z2ZDriT6VRJ+pa9l66LTVCZVSl8koCBnLnDlz\nAFi7di1t2rQB4I033sg8Mpk1axYA+/fv5+GHH2bnzp0EBQWxYMECAEaPHs2UKVPYvHkzEyZM4KGH\nHuLGG2+kT58+vP3220RGRnLdddfRv39/Nm3axLZt22jUqBEzZ87MMWd0dDQbN25k+fLl/OMf/8jz\n5bDff/89NWvWZNu2bURFRdG9u3PoPg8PD+rVq8e2bdsK9H0rLD3iyMpmhx5vQNWGsPxxbjkzlKXD\nZjJicTyDp//OmwOa0a9VqNUplXI7BT0yKKyMgnC5Y8eOERJy9RHC69SpQ8uWLQFo06YN0dHRnDt3\njnXr1jFw4MDMdhcvXsz2/VFRUTz33HOcOXOGc+fO0a1btxxzDho0CA8PD+rXr0/dunXZs2dPXr48\nmjVrxuOPP85TTz1Fr1696NSpU+ZrVatW5ejRo5kFsiTpEUd22oyEexdDUhx1FvVh2Z1ptKodxGNz\ntvHm93twOPSkuVLuzNfXN8e/6r29vTMf22w20tLScDgcBAUFXXLOZPfu3dm+f+TIkUydOpUdO3bw\n4osv5noEcfm9E5c/t9vtOByOzOcZ27v++uvZsmULzZo147nnnuPll1++pI2vrzVd6Fo4ribsJhi9\nGirWIGD+YGa12M7Qdtcwbc0BRn8ZwbmLaVYnVEpdRaNGjfjjjz8yn3t6epKamprjewICAqhTpw7z\n5s0DnMNzZHQFVaxYkcTExMy2iYmJ1KhRg9TU1Mzur5zMmzcPh8PBgQMHOHjwIA0aNLjk9bCwMLZs\n2QLAli1b+PPPPwE4evQofn5+DB8+nP/+97+ZbQD27duX76vDiooWjpxUCoP7V0L927F//ySveX3K\ny72uZ/XeWAZ8sI7D8UlWJ1SqXLv8HEfGVVV33nkna9asyWw3evRomjdvzj333JPj9mbNmsXMmTNp\n0aIFTZo0YfFi53isQ4YM4e2336ZVq1YcOHCAV155hfbt29OxY0caNmyYa85rrrmGdu3a0aNHDz78\n8EN8fHwueX3AgAHEx8fTpEkTpk6dyvXXXw/Ajh07aNeuHS1btmTcuHE899xzAJw4cQJfX1+qV89u\nHNniJ2XxXoXw8HBTpDMAOtLhp5dh7XsQ1onfw99l9PyD2DyED4e3oX3dykW3L6VKid27d9OoUSOr\nY2TrwoULdO7cmbVr12Kz2SzNMnLkSHr16sXdd99dZNt89913CQgI4P77CzYKU3afnYhsNsaE5+X9\nesSRFx42uH0c9PsIDm+gw08DWT60KpX8vbjn4w18s/GQ1QmVUln4+voybtw4jhw5YnWUYhEUFMSI\nESMs279eVZUfLYZAcF345h5qL+zD0t4fMWZjFcYu3MHeE4k827MRdpvWYqXcQW5XOhW18ePHZ54f\nyTBw4EA+++yzIt/XqFGjinyb+aFdVQVx5jB8MxSOR5He9WXGx3fhk3XRdKpfhalDWxPoV/AJUpQq\nLdy5q0rlTLuqrBBUG+77ARr1xrbqeV5wvM/bdzXg94Nx9PtgLQdjz1mdUCmlio0WjoLy8oeBn8Mt\nT0HkLAbufIg599TjzIVU7np/Lb/u13nPlVJlkxaOwvDwgM7PwN2fwLFttP6hPysGBVEzyJeRn27i\ns7V/6gi7SqkyRwtHUWg6AEZ9B440qs/vy6Iup+ncoCovLd3FM4t2kJLmyH0bSilVSmjhKCq1WsOD\nqyGkAb4L7mV62BoeuqUuszceZvjMDcSfT7E6oVJlTmmej+NyPXv25MyZM0W2vZSUFG6++WbS0op+\nlAur5hwfKCI7RcQhItmexReR2iKyWkR2udo+WtI58y2gBoxaAc0G4rH6FZ48P4EpAxsSefgMfd//\njb3HE3PfhlIqz8rSfBwrVqwgKCioyLbn5eVFly5dMkcMLkpW3ccRBfQHPsqhTRrwuDFmi4hUBDaL\nyI/GGPeZTSU7nr7QfwZUbQQ/vUzv+IPUGf4hoxYcpv8Ha5k0pBVdG1ezOqVSReu7sXB8R9Fus3oz\n52jVBeBu83H4+vqydetWTp48ySeffMIXX3zB+vXrad++feZ9HmFhYURERHDu3Dl69OjBTTfddEW+\nW2+9lQkTJhAeHs6pU6cIDw8nOjqanTt3MmrUKFJSUnA4HCxYsID69etz11138fTTT+c61Ep+WXLE\nYYzZbYzZm0ubY8aYLa7HicBuoFZJ5Cs0Eej0OAyeBbF7abq8L9/d7U/dkAo8+GUE09Yc0JPmShWB\n0jIfx+nTp1m/fj3vvvsuffr04bHHHmPnzp3s2LEj22Hhr5bvaj788EMeffRRIiMjiYiIIDTUOf1D\n06ZN2bRpU/6+qXlQKu4cF5EwoBWwIYc2o4HR4BxQzC006gX3/wCzh1Jl3l0s6DWV/+y6jje/38O+\nE4m83r8ZPp7WjqOjVJEo4JFBYZWW+Th69+6NiNCsWTOqVatGs2bNAGjSpAnR0dGZWXLKl5MbbriB\n8ePHExMTQ//+/alfvz7gPAfk5eVFYmIiFStWzHEb+VFsRxwiskpEorJZ+uZzOxWABcC/jTFXnZvS\nGDPdGBNujAnP6QemxFVv5jxpXqMlXt8+wJTqK3i8az0WbT3CkOm/czIxbzOBKaXyzt3m48jYn4eH\nxyX79vDwyPbkdXb54NJ5O7Luc9iwYSxZsgRfX1969uzJzz//nPnaxYsXrxiNt7CKrXAYY7oaY5pm\nsyzO6zZExBNn0ZhljFlYXFmLXYUQGLEEWg5Hfnmbf8W9wowhDdl7PJG+U9cSdSTB6oRKlSnuNh9H\nUQkLC2Pz5s0AzJ8/P3P9wYMHqVu3Lo888gh9+/Zl+/btAMTFxVGlShU8PYt2GCS3vRxXnFNkzQR2\nG2MmWp2n0Oze0Hcq3DEe9izn9t9HsPhv1yDA3R+uY/n2Y1YnVKrUKS3zcRSVJ554gmnTptGqVStO\nnTqVuX7u3Lk0bdqUli1bEhUVxb333gvA6tWrufPOO4s8hyWDHIpIP2AKEAKcASKNMd1EpCbwsTGm\np4jcBPwK7AAy7qB7xhizIrftF/sgh4W1byUsuB/sPpzu8yn3/yRsOXSGR7vU59Eu9fHwkNy3oZTF\n3HmQQ3eaj8NK/fv354033sicGCpDqRzk0BizyBgTaozxNsZUM8Z0c60/aozp6Xr8mzFGjDHNjTEt\nXUuuRaNUuP4OuP9H8PKn0tx+zOkQzYDWoUz6aT//nL2FpBSdllapwijr83HkRUpKCnfdddcVRaMo\nuG1XVZlXtSE8+DPUbo/nkoeYELSA53pcz3dRxxn44XqOnrlgdUKlSrVu3bqV6BWW48ePv6TbrGXL\nlowfP77E9n85Ly+vzC6roqbzcVgtPRW+exIiPoHru/NLs9d4aP4f+HjamH5vG1pfU8nqhEply527\nqlTOSmVXlcrC5gm93oWeE2D/j9z8yz0sHV4Lf28bQz76nYVbYqxOqJRSl9DC4S7aPQjDF0DiUeos\n7M2y3tDm2kr8Z+42Xv9uN+mOsndkqJQqnbRwuJPrOjtvFvSrTMW5d/Nlq93c0/4aPvrfQUZ/EUFi\ncs7XoSulVEnQwuFuKl8HD6yCOrdgX/5vxvt8xau9G7BmXywDpq3jUFyS1QmVcguXD6keHR1NREQE\njzzySIG3GRYWdsn9EVnXN2vWLHNfhdlHWVAqxqoqd3yDYNhc+PF5+P0DhtfdS7173uHv8w/Q9/3f\nmDa8DR3qVrY6pVKWym6cqrCwMMLD83R+N99Wr15NlSpVimx7aWlp2O2l81ewHnG4K5sdur8OvSdD\n9G90+HkQy4ZVJ9jfi+Efb+DrDYesTqiU21mzZg29evUC4KWXXuK+++7j1ltvpW7dukyePDmz3V13\n3UWbNm1o0qQJ06dPL/D+br31Vp566inatWvH9ddfz6+//gpAeno6//3vf2nbti3Nmzfno48+yszX\nqVMn+vTpQ+PGjQF45ZVXaNCgATfddBNDhw5lwoQJHDhwgNatW2fuZ//+/Zc8t1rpLHflSZsRULke\nzP0btRf0YknfmTy0vgrPLNrBvhOJPHdnI+w2rf/KWm9ufJM98XuKdJsNgxvyVLunrvp6xnAj4BxN\ndtGiRVe02bNnD6tXryYxMZEGDRowZswYPD09+eSTTwgODubChQu0bduWAQMGULlyzkfxnTt3zrwL\nfcSIETz22GOA88hh48aNrFixgnHjxrFq1SpmzpxJYGAgmzZt4uLFi3Ts2JE77rgDgC1bthAVFUWd\nOnXYtGkTCxYsYNu2baSmptK6dWvatGnDddddR2BgIJGRkbRs2ZJPP/2UUaNGFej7WBy0cJQGYR2d\nNwvOHor/3MF82u11Xgu5iY/XRvPHyXO8P6w1gX5FO4iZUu7uakOqZ3XnnXfi7e2Nt7c3VatW5cSJ\nE4SGhjJ58uTMQnP48GH279+fa+G4WldV//79gUuHP1+5ciXbt2/PHIgwISGB/fv34+XlRbt27ahT\npw7gnDekb9+++Pj44OPjQ+/evTO3+8ADD/Dpp58yceJE5syZw8aNG/P2jSkBWjhKi0phcP9KWPAg\nHt8/yXNtRtGg3yM8s2QPd32wlo9HhHNdSAWrU6pyKqcjAytlNzz5mjVrWLVqFevXr8fPz49bb701\n12HR87KPrMOfG2OYMmXKFfN0rFmzBn9//zxtd8CAAYwbN47bbruNNm3a5FrYSpL2cZQm3hVhyCy4\n6THY/CkDdz/CnL814OyFVO56fy2/7Iu1OqFSbi8hIYFKlSrh5+fHnj17+P3334t8H926dWPatGmZ\nQ7nv27eP8+fPX9GuY8eOLF26lOTkZM6dO8eyZcsyX/Px8aFbt26MGTPGrbqpQAtH6eNhg64vQb+P\n4PAGWv/Qn+VDq1AryJeRn27kk9/+1GlplcpB9+7dSUtLo1GjRowdO5YOHTrk6X2dO3fOvBw3tzGg\nHnjgARo3bkzr1q1p2rQpf//737OdsKlt27b06dOH5s2b06NHD5o1a0ZgYGDm6/fccw8eHh6Z50fc\nhY5VVZod3gTfDIPUCyT3nc4jm6uyctcJhrStzct9m+Jl178LVPHRsaqKxrlz56hQoQJJSUncfPPN\nTJ8+PfMKqgkTJpCQkMArr7xSpPss7FhVeo6jNKvdFkavhtlD8Zk3jA+7vszEqncwdc0BDsaeZ9rw\n1lSu4J37dpRSlhk9ejS7du0iOTmZESNGZBaNfv36ceDAgUumgXUXesRRFqSch0X/gN1LoMUwll7z\nJE8s2kNIRW8+HhFOw+oBVidUZZAecZReOjquAi9/GPg53DIWtn1N78i/s/DeeqSkORjwwTp+3HXC\n6oRKqTLEksIhIgNFZKeIOEQk2wonIj4islFEtrnajivpnKWKhwd0fhru/hSObafJ8n6sGFyJ66pW\nYPSXEXyw5g89aa6UKhJWHXFEAf2BX3JocxG4zRjTAmgJdBeRvF3+UJ417Q/3fQeOdKrM6cX8W+Lo\n1bwmb32/l8fmRJKcmm51QqVUKWfVnOO7jTF7c2ljjDHnXE89XYv+yZwXNVs5T5pXbYTXgnuZXHMV\nT9xen28jjzJ4+u+cPFvwm52UUsqtz3GIiE1EIoGTwI/GmA1WZyo1KlaHkcuh2SBk9av88/QbTB/a\nhP0nEukzdS07YhKsTqhUoVw+rPobb7xRZNuOjIxkxYoV2b62Zs0aAgMDL9n3qlWrimzfpUGxXY4r\nIquA6tm89KwxZnFetmGMSQdaikgQsEhEmhpjoq6yv9HAaKBEJ6h3a56+0H86VG0IP73MHfEH+fZv\n0xm1IIaBH61jwsAW9Gpe0+qUShVIXsaqKqjIyEgiIiLo2bNntq936tTpkru8i0JpGma92I44jDFd\njTFNs1nyVDQu29YZYDXQPYc2040x4caY8JCQkMJEL1tEoNPjMORriN3H9Uv6sHyAH01rBvLPr7cy\n8cd9OHRaWlVGJCQk0KBBA/budfaEDx06lBkzZgAwZswYwsPDadKkCS+++GLmezZt2sSNN95IixYt\naNeuHQkJCbzwwgvMmTOHli1bMmfOnDztOzo6mkaNGvHggw/SpEkT7rjjDi5cuADAgQMH6N69O23a\ntKFTp07s2eMcSXjkyJH84x//oH379jz55JPExsZy++2306RJEx544AGuvfZaTp06xQsvvMB7772X\nua9nn32WSZMmFcn3rCDctryJSAiQaow5IyK+wO3AmxbHKr0a3ukcJHH2UILm9GV27yk8U6UBk3/a\nz/4TibwzqAV+Xm7746Dc3PHXXuPi7qIdVt27UUOqP/PMVV/POqw6wNNPP83gwYOZOnUqI0eO5NFH\nH+X06dM8+OCDAIwfP57g4GDS09Pp0qUL27dvp2HDhgwePJg5c+bQtm1bzp49i5+fHy+//DIRERFM\nnTo1233/+uuvl+x7wYIF2Gw29u/fz+zZs5kxYwaDBg1iwYIFDB8+nNGjR/Phhx9Sv359NmzYwEMP\nPZR5Y19MTAzr1q3DZrPxz3/+k9tuu42nn36a77//npkzZwJw33330b9/f/7973/jcDj45ptvLB0t\n15LfFCLSD5gChADLRSTSGNNNRGoCHxtjegI1gM9FxIbzyGiuMaZojw3Lm+pNncOzz/0bnose5K1O\nT9Cg51Be+24vf01LYsaIcGoF+VqdUqk8uVpX1e233868efN4+OGH2bZtW+b6uXPnMn36dNLS0jh2\n7Bi7du1CRKhRowZt27YFICAgbzfLZtdVFR0dTZ06dTILSsYw6+fOnWPdunUMHDgws+3FixczHw8c\nODBzno/ffvstc7j37t27U6lSJcA5s2HlypXZunUrJ06coFWrVpaOlmtJ4TDGLAKumHXFGHMU6Ol6\nvB1oVcLRyr4KIXDvYlj2H+TXCTzQaC/1h7/KP+fupe/U3/job21oc22w1SlVKZPTkUFJczgc7N69\nGz8/P06fPk1oaCh//vknEyZMYNOmTVSqVImRI0cWaij1q7l8GPcLFy7gcDgICgq66vmYvA6z/sAD\nD/DZZ59x/Phx7rvvviLJW1BufVWVKiZ2b+g7Fbq9BnuWc8svw1l677X4e9sZOn0D8zfHWJ1QqQJ7\n9913adSoEV9//TWjRo0iNTWVs2fP4u/vT2BgICdOnOC7774DoEGDBhw7doxNmzYBkJiYSFpaGhUr\nViQxMbFI8gQEBFCnTh3mzZsHOOfqyHoklFXHjh2ZO3cu4JwM6vTp05mv9evXj++//55NmzZdMc9H\nSdPCUV6JwA0Pw7C5cOYvwhb2YmlfT8LDKvHEvG28vmI36XrSXLmxjHMcGcvYsWPZu3cvH3/8Me+8\n8w6dOnXi5ptv5tVXX6VFixa0atWKhg0bMmzYMDp27AiAl5cXc+bM4V//+hctWrTg9ttvJzk5mc6d\nO7Nr166rnhzPOMeRsWTM9Hc1s2bNYubMmbRo0YImTZqweHH21wi9+OKLrFy5kqZNmzJv3jyqV69O\nxYoVM7N27tyZQYMGZXZtWUUHOVQQuxe+Hgxnj5B253uMO9SCL3//i9saVmXSkJZU9NFpadWVdJDD\nonfx4kVsNht2u53169czZsyYzC4uh8NB69atmTdvHvXr1y/UfnSQQ1V4IQ2cJ82v6YB9yUO84jeH\nV/s24n/7Yun/wTr+irty5jKlVNE7dOgQbdu2pUWLFjzyyCOZlxLv2rWLevXq0aVLl0IXjaKgRxzq\n/6WnwndPQcRMqN+NDa3eZPS8/YjAB/e05sbrqlidULkRPeIovfSIQxUdmyf0mgg9J8Afq2i/egjL\nh4dSpYI3987cyKwNf1mdULmZsviHZ1lXFJ+ZFg51pXYPwt8WQuJxQuffyeI7HXSqX4VnF0XxwuIo\nUtMdVidUbsDHx4e4uDgtHqWIMYa4uDh8fHwKtR3tqlJXF3fAedL89J84ekzg9ZPtmfHrn3SsV5n3\nh7UmyM/L6oTKQqmpqcTExBTL/RCq+Pj4+BAaGoqn56UXveSnq0oLh8rZhTMw/z448BO0+zvzq4zh\nmcV7qBnkw8cj2lKvagWrEyqlioCe41BFxzfIea9Hh4dh40fcvecx5tzbkHMX0+j3wVrWH4izOqFS\nqoRp4VC5s9mh+2vQZwpE/0arlXezbFh1qlb05tFvtnI2OdXqhEqpEqSFQ+Vd63thxBK4cJrqc+9k\nesdEYs9dZOLKfVYnU0qVIC0cKn+uvREeXA0BtbjuhxFMuH43n6+PZnvMGauTKaVKiBYOlX+VrnXO\n7XFNB/ofm0gT/0SeWbRDx7ZSqpzQwqEKxrsi3PUB4nDwcZVviDqSwBfro61OpZQqAVo4VMFVCoPO\nz1D9+Gr+E7qXd1bu43iCXtOvVFmnhUMVToeHoHozHrrwEd7p53h52U6rEymlipkWDlU4Njv0noQ9\nKZbPaq9gxY7jrN5z0upUSqliZEnhEJGBIrJTRBwikuOdiiJiE5GtIqLzjburWm2g3d9pemwBvYMP\n8fziKC6kpFudSilVTKw64ogC+gO/5KHto8Du4o2jCu22Z5GAWrzp9QknTicy5ef9VidSShUTSwqH\nMWa3MWZvbu1EJBS4E/i4+FOpQvGuCHdOwO/MPibV/oXpvxxk34mimbNZKeVe3P0cx3vAk4CO410a\nNOgBjfvSI/5LGnvH8uyiHTj03g6lypxiKxwiskpEorJZ+ubx/b2Ak8aYzXlsP1pEIkQkIjY2tlDZ\nVSH0eAux+zAj+Cs2Rcczf3OM1YmUUkXMXlwbNsZ0LeQmOgJ9RKQn4AMEiMhXxpjhV9nfdGA6OIdV\nL+S+VUFVrA5dX6Ta8v/wRLXNvPadF10bVyPYX+fuUKqscNuuKmPM08aYUGNMGDAE+PlqRUO5mTaj\noHZ7xiR/gldyPK+t0GsblCpLrLoct5+IxAA3AMtF5AfX+poissKKTKoIeXhA70nYUs/zSc1FzN8c\nw+8Hdd4OpcoKq66qWuQ6mvA2xlQzxnRzrT9qjOmZTfs1xpheJZ9UFVjVRnDTv2l66nvuCtjLc99G\nkZKm1zgoVRa4bVeVKgM6PQHB1/G65yfEnIxjxq8HrU6klCoCWjhU8fH0gd7v4Xv+MO9VX8nkn/Zz\nKC7J6lRKqULSwqGKV52boeVwuiXMpbGHczgSY/SiN6VKMy0cqvjd8QriG8T0oC/4dd8JVuw4bnUi\npVQh5KlwiMhCEblTRLTQqPzzC4bubxByNoong39l3NKdnE1OtTqVUqqA8loIPgCGAftF5A0RaVCM\nmVRZ1GwgXHcbD6Z+he3cUSau3Gd1IqVUAeWpcBhjVhlj7gFaA9HAKhFZJyKjRMSzOAOqMkIE7pyI\nzTj4pOo8Pl8fzfaYM1anUkoVQJ67nkSkMjASeADYCkzCWUh+LJZkquwJrgO3jqVRwi8M9IvkmUU7\nSNdBEJUqdfJ6jmMR8CvgB/Q2xvQxxswxxvwLqFCcAVUZc8PDUK0pL3t+RvSR43yxPtriQEqp/Mrr\nEcdkY0xjY8zrxphjWV8wxuQ4g59Sl7B5Qu/JeCfH8m6VJbyzch/HE5KtTqWUyoe8Fo5KItL/sqWL\niFQt1nSqbAptg7QbTddzS2mSvoeXl+20OpFSKh/yWjjuxzkL3z2uZQbwFLBWRP5WTNlUWdbleSSg\nJh8EfM7KHTGs3nPS6kRKqTzKa+HwBBoZYwYYYwYAjQEDtMdZQJTKH++K0HMClZMOMDbgR55fHMWF\nlHSrUyneLqraAAAaE0lEQVSl8iCvhSPUGHMiy/OTQG1jTDygd3KpgmnYExr15r70udjO/MmUn/db\nnUgplQd5LRxrRGSZiIwQkRHAYtc6f0AvxlcF1+MtPOzezAiexfRfDrDvRKLViZRSuchr4XgY+BRo\n6Vq+AB42xpw3xnQurnCqHAioCV1e4PrzmxnivY5nF+3Aofd2KOXWci0cImLDOW3rAmPMY65lvtEh\nTlVRCb8fQtvxvP0r/oj+i/mbY6xOpJTKQa6FwxiTDjhEJLAE8qjyyDXVrFf6Od4Nmsdr3+0m/nyK\n1amUUleR166qc8AOEZkpIpMzloLuVEQGishOEXGIyFVvIBSRaBHZISKRIhJR0P2pUqBaY6Tjo9ya\n/BPNLkby2ordVidSSl2FPY/tFrqWohIF9Ac+ykPbzsaYU0W4b+Wubv4v7FzE5POf02Fzfe5uE0qH\nupWtTqWUukyeCocx5nMR8QWuMcbsLexOjTG7AUSksJtSZYmnL/R6l0pf9OXZCkt57ttgVjzSCS+7\nTgOjlDvJ6yCHvYFI4HvX85YisqQ4g7kYYKWIbBaR0SWwP2W1urdCi2EMT1+MLXYXM349aHUipdRl\n8vqn3EtAO1z3bBhjIoG6Ob1BRFaJSFQ2S9985LvJGNMa6AE8LCI357C/0SISISIRsbGx+diFcjt3\nvIqHTwDTAr5gyk97ORSXZHUipVQWeS0cqcaYhMvWOXJ6gzGmqzGmaTbL4ryGM8Yccf17EliEs3hd\nre10Y0y4MSY8JCQkr7tQ7si/MnR/nboXd3GPxyqeXxyFXv2tlPvIa+HYKSLDAJuI1BeRKcC6YsyF\niPiLSMWMx8AdOE+qq/Kg+WCoeytPec5hz769rNhx3OpESimXvBaOfwFNgIvAbOAs8O+C7lRE+olI\nDHADsFxEfnCtrykiK1zNqgG/icg2YCOw3BjzfUH3qUoZEej1Lp6SzsSArxm3dCdnk3VYNKXcgZTF\nLoDw8HATEaG3fZQJv06En8YxOvUxarYfyEt9mlidSKkySUQ253VivrxeVXW9iEwXkZUi8nPGUriY\nSuXBjf+Cqk142+8rFqzfxfYYHVNTKavltatqHrAVeA74b5ZFqeJl84TekwhIPcVzvgt5ZtEO0nUQ\nRKUsldfCkWaMmWaM2WiM2ZyxFGsypTLUbou0e5BBju+wH93CF+ujrU6kVLmW18KxVEQeEpEaIhKc\nsRRrMqWyuu15qFiDyRU+Y9LKXRxPSLY6kVLlVl4LxwicXVPrgM2uRc8+q5LjE4D0fJtrUg8y3LGU\nl5fttDqRUuVWngqHMaZONkuOd44rVeQa9YKGvXjUcyFRUdtYveek1YmUKpdyLBwi8mSWxwMve+21\n4gql1FX1eAu73ZOJfp/z/Lc7uJCSbnUipcqd3I44hmR5/PRlr3Uv4ixK5S6wFtL1JcLTIwk/u4op\nP++3OpFS5U5uhUOu8ji750qVjPD7oFY4r/rOYu4v29h3ItHqREqVK7kVDnOVx9k9V6pkeNig9yT8\nzXme9/qaZxftwKH3dihVYnIrHC1E5KyIJALNXY8znjcrgXxKZa96U+TGf9GXNdgP/cb8zTFWJ1Kq\n3MixcBhjbMaYAGNMRWOM3fU447lnSYVUKlu3PIWpVIeJvp/yzopI4s+nWJ1IqXJB5+RUpZenL9Lr\nXWqkH+XetPm8vmK31YmUKhe0cKjS7brO0HwI/7AvZduW9Ww4GGd1IqXKPC0cqvTrNh4Pn0Am+H7K\nc4u2k5KW4+SUSqlC0sKhSj//Kki38TR37KFd/BJm/HrQ6kRKlWlaOFTZ0GIo1LmZ57zn8M1PGzgU\nl2R1IqXKLC0cqmwQgV7v4SNpPG/7nOcXR1EWZ7dUyh1YUjhEZKCI7BQRh4hcdapCEQkSkfkiskdE\ndovIDSWZU5Uyla9DbnmSO2QDnn98z4odx61OpFSZZNURRxTQH/gll3aTgO+NMQ2BFoBeb6lyduMj\nmJBGvO7zOW8vieBscqrViZQqcywpHMaY3caYvTm1EZFA4GZgpus9KcYYnXBa5czuhfSZTBVHHPcm\nf8XElfusTqRUmePO5zjqALHApyKyVUQ+FhF/q0OpUqB2O6Tt/Yy0/8DW339ie4z+vaFUUSq2wiEi\nq0QkKpulbx43YQdaA9OMMa2A88DYHPY3WkQiRCQiNja2CL4CVap1eQEqVOVtr5k8v3Ar6ToIolJF\nptgKhzGmqzGmaTbL4jxuIgaIMcZscD2fj7OQXG1/040x4caY8JCQkMLGV6WdTyAePSdwPdG0PzGH\nL9ZHW51IqTLDbbuqjDHHgcMi0sC1qguwy8JIqrRp1BvToAePey3km5W/cTwh2epESpUJVl2O209E\nYoAbgOUi8oNrfU0RWZGl6b+AWSKyHWgJ6HS1Ku9EkJ4T8LTbec7M4OWlUVYnUqpMsOqqqkXGmFBj\njLcxppoxpptr/VFjTM8s7SJd3U/NjTF3GWNOW5FXlWKBoXh0fZFOHtuw71rI6j0nrU6kVKnntl1V\nShWZtg/gqNmacd5f8da367mQkm51IqVKNS0cquzzsOHRZzJBnGPkuU+Y8vN+qxMpVapp4VDlQ/Vm\nyI3/ZLB9Ddt+Xca+E4lWJ1Kq1NLCocqPW8aSHngN4z0/5qWFm3HovR1KFYgWDlV+ePlh6/0eYRyj\n/ZHPmL8lxupESpVKWjhU+VKvC6bZIB62L2XO8pXEn0+xOpFSpY4WDlXuSLfXEO8KPJ3+EW8s32l1\nHKVKHS0cqvypEIKt+3jCPfZi2/YlGw7GWZ1IqVJFC4cqn1reQ/o1N/GM52zeWfgLKWkOqxMpVWpo\n4VDlkwi2PpPw80jj3oRpzPj1oNWJlCo1tHCo8qtKPWy3PEkv2wa2//wNh+KSrE6kVKmghUOVbx0f\nJTW4AS96fMqrizZhjN7boVRutHCo8s3uheddk6kpp2gXPY0VO45bnUgpt6eFQ6lrOuBoM4pR9h+Y\nu3gxZ5NTrU6klFvTwqEU4NH1JRy+lXky9QPe+0HnC1MqJ1o4lALwDcKz1wSaePyFbdNHbI85Y3Ui\npdyWFg6lMjTuS2q9bvzHPp/35q8iXQdBVCpbWjiUyiCCZ6938LTZuDduMl+s+9PqREq5JavmHB8o\nIjtFxCEi4Vdp00BEIrMsZ0Xk3yWdVZUzQbWxdX2eW23b2PXjZxxPSLY6kVJux6ojjiigP/DL1RoY\nY/YaY1oaY1oCbYAkYFEJ5VPlmLT/OxertuApPuPtxb9bHUcpt2NJ4TDG7DbG7M3HW7oAB4wxfxVX\nJqUyedjw7jeFYDlH+L73WL3npNWJlHIrpeUcxxBgttUhVDlSowWOGx5iqH018xfO4UJKutWJlHIb\nxVY4RGSViERls/TN53a8gD7AvFzajRaRCBGJiI2NLUx0pQCwd36a5Aqh/Cf5fT5YpfN2KJWh2AqH\nMaarMaZpNsvifG6qB7DFGHMil/1NN8aEG2PCQ0JCCh5cqQxe/vj0ncR1HsfwXP8e+04kWp1IKbdQ\nGrqqhqLdVMoq9btysWF//mFbzLR5K3DovR1KWXY5bj8RiQFuAJaLyA+u9TVFZEWWdv7A7cBCK3Iq\nBeDd602Mpx9DTk5k/uZDVsdRynJWXVW1yBgTaozxNsZUM8Z0c60/aozpmaXdeWNMZWNMghU5lQKg\nQlU8u4+nvcce9qx4n/jzKVYnUspSpaGrSinLebT+G0k1O/Co4yumLv7N6jhKWUoLh1J5IYJf/6n4\ne6TQavebbDgYZ3UipSyjhUOpvKpSH3PT4/S2/c7ieZ+RkuawOpFSltDCoVQ+eN7yH84HXMdDSR/w\n2Zooq+MoZQktHErlh90b/wHvEyqnsP/yBofikqxOpFSJ08KhVH5dewNJzf7GCFnBx/MWYYze26HK\nFy0cShWAX89XuehdmYFH3+K7bTFWx1GqRGnhUKogfIPw7vUWzTyi2btkAmeTU61OpFSJ0cKhVAHZ\nmvXnbO3b+Hv6bGYu/Z/VcZQqMVo4lCooEQIGTMJuE5pvf5Xth09bnUipEqGFQ6nCCLoGx63P0sW2\nlRVzppGugyCqckALh1KF5NPxIc4ENeH+xA/55pftVsdRqthp4VCqsGx2Age9T2VJxHP1yxxPSLY6\nkVLFSguHUkVAarYisdWDDJJVzJqn08eosk0Lh1JFJLDHi5z1rkGfQ2+zZudhq+MoVWy0cChVVLz8\n8e33HvU9jnDg21e5kJJudSKlioUWDqWKkGfD7pwK68XwlPnMWv6j1XGUKhZaOJQqYlUGTCTd5kPz\nrS+y77hOXqnKHqvmHB8oIjtFxCEi4Tm0e8zVLkpEZouIT0nmVKpAKlbD0fVl2nns4aevJ+LQeztU\nGWPVEUcU0B/45WoNRKQW8AgQboxpCtiAISUTT6nCqdBhFCeDWzM0YQZL10daHUepImVJ4TDG7DbG\n7M1DUzvgKyJ2wA84WrzJlCoiHh5UGTwNf7mI14/PEn8+xepEShUZtz3HYYw5AkwADgHHgARjzEpr\nUymVdx7VGnI2/F/0YC0L53xmdRylikyxFQ4RWeU6N3H50jeP768E9AXqADUBfxEZnkP70SISISIR\nsbGxRfNFKFVIlbs/TZzPtXSLfotNe/XeDlU2FFvhMMZ0NcY0zWZZnMdNdAX+NMbEGmNSgYXAjTns\nb7oxJtwYEx4SElIUX4JShWf3psLd71PbI5Y/FzxPSprD6kRKFZrd6gA5OAR0EBE/4ALQBYiwNpJS\n+eddrxNH6w6i/4H5vPPJ51RpeBOhlStyTbAftYP9qODtzv8NlbqSWDFfsoj0A6YAIcAZINIY001E\nagIfG2N6utqNAwYDacBW4AFjzMXcth8eHm4iIrTGKDdy4TTn3mlNhbR40owHJwniuAnmmAnmtL0q\nF/2qQ0BNvINrU7HqtYTUuIbQyoHUCPLB0+a2pyJVGSIim40xV7094pK2VhSO4qaFQ7ml09GYP37m\nYvxhkk4dJj3hCPZzR/FLPoG348IlTR1GiCWQY6YyZzxDSPapjqNiDeyVQvGrcg3BNcKoVqsOwQEV\nEBGLviBVluSncOgxslIlpVIY0vY+fIBL7mQ1BpIT4OxR0s/EkHDyL87HHiL1dAyBZ49Q7cJxApK2\n438+CY5fuslTJpB4ewjnvauS5l8Dj8Ba+FSpTWC1MEJq1cW7Uih46n2zqmhp4VDKaiLgGwS+Qdiq\nNSa4AQRn1y75LElxh4k7Gs3Zk39xMe4QJuEIXknHCUw+QvD5bQTFnoc/Ln1bggSQ6OXsDjMBtfAO\nDqVCyLUEVrsWj8BQCKgJXn4l8ZWqMkILh1KlhU8AfrWa4FerSbYvG2OIPX2aEzEHSTgRTVLsYdLO\nxGA/dwzf5BNUijtMjfhIgv86d8V7k2wVueBbnfQKNbEH1cKvyjX4VK7tLCoBruLiXaG4v0JVSmjh\nUKqMEBFCgoMJCQ4GruyqTklzcPTMBXadjCfueDTnYv8iNf4IkngE76TjVE44RfWzf1H92FZ85OyV\n77dXJNW/Oh6BtfAKro0tsBYE1nIVF9e/3gHOIyhVpmnhUKqc8LJ7EFbFn7Aq/tC4NtDpktcTLqRy\nOD6JzaeTiDl1hoQTh0iOi8EkxOCVdJyQtDhqXIyn+ulD1PgrkhBJwINLL65Jt/tjAmpiC6qFZByp\nZC0sATXBt5IWl1JOC4dSCoBAX08CawXStFYgUANolPmaw2E4kZjMobgk/jh9gZ/jkzgWl0DiqRjS\nTh/B58Jxqks8NdLiqX4xjtC4GGp67KCyOY0Hl970aDz9kOwKSkCt/1/8grW4uDEtHEqpXHl4CDUC\nfakR6Ev7S15pC0Byajoxp5M4FJ/E4fgLbIlP4nB8EjFxiSSfPkpgaizVJZ6aEkf1tHiuTTtD7bOx\nVDN7CUo7hQeXzZZo83YWE/8qYPcBTz/w9P3/xZ7xOMtrdt+rtLlsnd1bi1IhaeFQShWaj6eNelUr\nUq9qxSteM8ZwOinVVVScxeXn+CQOuwrN8aQkghxnnEVF4qnlcZr6nglcm5xA8MVEvEjA25zAy6Tg\n6UjG01zEM/0idkfyFUczeWEQjN0HY/cFTz+Mpw/iKizi6Yd4+iCefq6C5HOVIuSXpaBd3sbv/wua\nzbMovr1uRwuHUqpYiQjB/l4E+3vRsnbQFa+npTs4lpB8SWFZd/oC38Sd53RSKmnpDlLSDWkOB6lp\nDlLTDSnpDsDgSTq+XMSbFHwlBR9S8OUiPqTiKxfxwbXO9diXFLwlBd+0jHbO93mTgi/n8ZXTmet8\nScEnyzYLIh0baR7ezsXmQ7qHj/Nfuw/G5kO63Rdj88F4/v+/GUVJvPwQL188PP0QLx88vPyxefvh\n4eWL3dv52O7th6e3H3Zvf8RWcr/OtXBkseyR/ngeiLE6hlLlViXX0iIf7zFZHphL1pts14Pznsus\n7ZKBZAPxAHgBXpdt1yCutmKy7tVcskHJsl4uWZ/RNhUhBeHsJe0l63sKIamy0G/u7kJuJXdaOJRS\npZpkeSDZvZLrb+OSP99xebHLeO64pMoZMBnlL2O9ueyx89+MYpbqWTI3cmrhyKLX5IVWR1BKKben\nw24qpZTKFy0cSiml8kULh1JKqXzRwqGUUipftHAopZTKFy0cSiml8kULh1JKqXzRwqGUUipfxJjL\nb8Yv/UQkFvirgG+vApwqwjhFRXPlj+bKH82VP2Ux17XGmJC8NCyThaMwRCTCGHPl9GkW01z5o7ny\nR3PlT3nPpV1VSiml8kULh1JKqXzRwnGl6VYHuArNlT+aK380V/6U61x6jkMppVS+6BGHUkqpfNHC\n4SIi3UVkr4j8ISJjLc7yiYicFJGoLOuCReRHEdnv+rdSCWeqLSKrRWSXiOwUkUfdJJePiGwUkW2u\nXONc6+uIyAbX5zlHRLxKMleWfDYR2Soiy9wsV7SI7BCRSBGJcK2z9LN0ZQgSkfkiskdEdovIDVbn\nEpEGru9TxnJWRP5tdS5XtsdcP/dRIjLb9f+h2H/GtHDg/M8NvA/0ABoDQ0WksYWRPgO6X7ZuLPCT\nMaY+8JPreUlKAx43xjQGOgAPu75HVue6CNxmjGkBtAS6i0gH4E3gXWNMPeA0cH8J58rwKJB1Lk93\nyQXQ2RjTMsvlm1Z/lgCTgO+NMQ1xziC72+pcxpi9ru9TS6ANkAQssjqXiNQCHgHCjTFNARswhJL4\nGTPGlPsFuAH4Icvzp4GnLc4UBkRleb4XqOF6XAPYa3G+xcDt7pQL8AO2AO1x3gRlz+7zLcE8oTh/\nodwGLMM5R6nluVz7jgaqXLbO0s8SCAT+xHXu1V1yXZblDmCtO+QCagGHgWCcs7kuA7qVxM+YHnE4\nZXwAGWJc69xJNWPMMdfj40A1q4KISBjQCtiAG+RydQdFAieBH4EDwBljTJqriVWf53vAk4DD9byy\nm+QC54TVK0Vks4iMdq2z+rOsA8QCn7q69z4WEX83yJXVEGC267GluYwxR4AJwCHgGJAAbKYEfsa0\ncJRCxvmnhCWXw4lIBWAB8G9jzFl3yGWMSTfOboRQoB3QsKQzXE5EegEnjTGbrc5yFTcZY1rj7J59\nWERuzvqiRZ+lHWgNTDPGtALOc1n3j8U/+15AH2De5a9Zkct1TqUvzoJbE/Dnyi7uYqGFw+kIUDvL\n81DXOndyQkRqALj+PVnSAUTEE2fRmGWMWeguuTIYY84Aq3EengeJiN31khWfZ0egj4hEA9/g7K6a\n5Aa5gMy/VjHGnMTZX98O6z/LGCDGGLPB9Xw+zkJida4MPYAtxpgTrudW5+oK/GmMiTXGpAILcf7c\nFfvPmBYOp01AfdfVCF44D0eXWJzpckuAEa7HI3CeYygxIiLATGC3MWaiG+UKEZEg12NfnOddduMs\nIHdblcsY87QxJtQYE4bz5+lnY8w9VucCEBF/EamY8Rhnv30UFn+WxpjjwGERaeBa1QXYZXWuLIby\n/91UYH2uQ0AHEfFz/f/M+H4V/8+YVSeZ3G0BegL7cPaPP2txltk4+yxTcf4Vdj/O/vGfgP3AKiC4\nhDPdhPNQfDsQ6Vp6ukGu5sBWV64o4AXX+rrARuAPnF0L3hZ+nrcCy9wllyvDNteyM+Pn3erP0pWh\nJRDh+jy/BSq5SS5/IA4IzLLOHXKNA/a4fva/BLxL4mdM7xxXSimVL9pVpZRSKl+0cCillMoXLRxK\nKaXyRQuHUkqpfNHCoZRSKl+0cCjlIiKVs4yAelxEjmR57nVZ2x8y7oXIYXsxGfeY5GHf94lI9aJq\np1RxsufeRKnywRgTh/M+AkTkJeCcMWZC1jauG63EGNOtiHd/H84BGo8XUTulio0ecSiVCxGpJ855\nSGbhvGGuRtajCRFZ6hoscKeIPJDLtuwi8qVrLowoEXlERAbjLFhzMo5uRGSciGxytflQnLJr11ZE\n/ufa/3ciYuUAgKqc0MKhVN40xDnHQWPjGucpixHGmDZAW+A/uUzo0wbncObNjHMOhS+MMXNw3ok/\n2DjnfUgBJhlj2gLNcA433v3ydjiHaZ8EDHDt/yvglaL7kpXKnnZVKZU3B4wxEVd57TER6eN6HApc\nh3PYjOz8ATQQkcnAcmDlVdp1EZH/Aj5AFZzDZX93WZtGQBNglbMHDRvOIWqUKlZaOJTKm/PZrRSR\nrsDNQAdjzAUR+Q3nL/tsGWPiRKQ5ruHMgQHA6KxtRMQPmAq0NsYcEZFXr7JNAbYbYzoV5AtSqqC0\nq0qpwgkE4l1FownO7qqrEpEQnCfX5wEv4Bw2HCARyLhKyxfn5E+nXFduDciyiaztdgG1RKSda9te\nrgxKFSs94lCqcJYDo0VkF86pRDfk0r42MNN1dZYBnnKt/xT4WEQu4Jwb43OcheHYZdu8vN3dwGQR\nCcDZVfUOzhP4ShUbHR1XKaVUvmhXlVJKqXzRwqGUUipftHAopZTKFy0cSiml8kULh1JKqXzRwqGU\nUipftHAopZTKFy0cSiml8uX/AMQe9/M0Z0sdAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10d2a5128>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.arange(0, max_trials,save_step),output[2],label='E(theta_plus)')\n",
    "plt.plot(np.arange(0, max_trials,save_step),output[3],label='E(theta_minus)')\n",
    "plt.plot(np.arange(0, max_trials,save_step),np.ones(max_trials//save_step)*output[0],label='Final Energy')\n",
    "plt.plot(np.arange(0, max_trials,save_step),np.ones(max_trials//save_step)*exact,label='Exact Energy')\n",
    "plt.legend()\n",
    "plt.xlabel('Trial state')\n",
    "plt.ylabel('Energy')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Optimizing the potential energy surface \n",
    "\n",
    "The optimization considered previously is now performed for two molecules, H$_2$ and LiH, for different interatomic distances, and the correspoding nuclei Coulomb repulsion is added in order to obtain a potential energy surface. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "OPTIMIZING HAMILTONIAN # 0 AT INTERATOMIC DISTANCE 0.2 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 1.13570872082\n",
      "objective function at theta+ for step # 0\n",
      "-0.515095690817\n",
      "objective function at theta- for step # 0\n",
      "-0.525349069681\n",
      "Final objective function is: -2.48741365446\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 1 AT INTERATOMIC DISTANCE 0.3 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 1.22663644505\n",
      "objective function at theta+ for step # 0\n",
      "-0.638388485733\n",
      "objective function at theta- for step # 0\n",
      "-0.644002581478\n",
      "Final objective function is: -2.36572817518\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 2 AT INTERATOMIC DISTANCE 0.4 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 1.61608737683\n",
      "objective function at theta+ for step # 0\n",
      "-0.756439332424\n",
      "objective function at theta- for step # 0\n",
      "-0.74663071353\n",
      "Final objective function is: -2.23709196923\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 3 AT INTERATOMIC DISTANCE 0.5 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 1.819928223\n",
      "objective function at theta+ for step # 0\n",
      "-0.830646452243\n",
      "objective function at theta- for step # 0\n",
      "-0.842640993766\n",
      "Final objective function is: -2.11351314334\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 4 AT INTERATOMIC DISTANCE 0.6 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 1.80872757334\n",
      "objective function at theta+ for step # 0\n",
      "-0.897497811088\n",
      "objective function at theta- for step # 0\n",
      "-0.893957030994\n",
      "Final objective function is: -1.99824502825\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 5 AT INTERATOMIC DISTANCE 0.7 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 1.98451681439\n",
      "objective function at theta+ for step # 0\n",
      "-0.932078419519\n",
      "objective function at theta- for step # 0\n",
      "-0.930205473868\n",
      "Final objective function is: -1.89215600477\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 6 AT INTERATOMIC DISTANCE 0.8 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.57512177917\n",
      "objective function at theta+ for step # 0\n",
      "-0.969907525698\n",
      "objective function at theta- for step # 0\n",
      "-0.972625640183\n",
      "Final objective function is: -1.79557049647\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 7 AT INTERATOMIC DISTANCE 0.9 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.38985906954\n",
      "objective function at theta+ for step # 0\n",
      "-0.946758401073\n",
      "objective function at theta- for step # 0\n",
      "-0.95900230368\n",
      "Final objective function is: -1.70849986845\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 8 AT INTERATOMIC DISTANCE 1.0 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.46098643244\n",
      "objective function at theta+ for step # 0\n",
      "-0.94639751431\n",
      "objective function at theta- for step # 0\n",
      "-0.950753259598\n",
      "Final objective function is: -1.63032695925\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 9 AT INTERATOMIC DISTANCE 1.1 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.26184086219\n",
      "objective function at theta+ for step # 0\n",
      "-0.939014720652\n",
      "objective function at theta- for step # 0\n",
      "-0.938080792348\n",
      "Final objective function is: -1.56020779157\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 10 AT INTERATOMIC DISTANCE 1.2 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.379128495\n",
      "objective function at theta+ for step # 0\n",
      "-0.920838239693\n",
      "objective function at theta- for step # 0\n",
      "-0.928489938075\n",
      "Final objective function is: -1.49772193121\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 11 AT INTERATOMIC DISTANCE 1.3 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.72098550207\n",
      "objective function at theta+ for step # 0\n",
      "-0.91997737694\n",
      "objective function at theta- for step # 0\n",
      "-0.917621631894\n",
      "Final objective function is: -1.44095794626\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 12 AT INTERATOMIC DISTANCE 1.4 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.63721973599\n",
      "objective function at theta+ for step # 0\n",
      "-0.887616754087\n",
      "objective function at theta- for step # 0\n",
      "-0.891145479395\n",
      "Final objective function is: -1.39344522831\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 13 AT INTERATOMIC DISTANCE 1.5 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.68225305556\n",
      "objective function at theta+ for step # 0\n",
      "-0.870369638497\n",
      "objective function at theta- for step # 0\n",
      "-0.868861894965\n",
      "Final objective function is: -1.34079916361\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 14 AT INTERATOMIC DISTANCE 1.6 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 4.03252195352\n",
      "objective function at theta+ for step # 0\n",
      "-0.85162104715\n",
      "objective function at theta- for step # 0\n",
      "-0.847122925931\n",
      "Final objective function is: -1.31313863003\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 15 AT INTERATOMIC DISTANCE 1.7 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.12739716655\n",
      "objective function at theta+ for step # 0\n",
      "-0.828434799826\n",
      "objective function at theta- for step # 0\n",
      "-0.829380227189\n",
      "Final objective function is: -1.28242444982\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 16 AT INTERATOMIC DISTANCE 1.8 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 4.55099555839\n",
      "objective function at theta+ for step # 0\n",
      "-0.812462659254\n",
      "objective function at theta- for step # 0\n",
      "-0.805433925014\n",
      "Final objective function is: -1.25546732201\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 17 AT INTERATOMIC DISTANCE 1.9 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.53659315179\n",
      "objective function at theta+ for step # 0\n",
      "-0.792485083358\n",
      "objective function at theta- for step # 0\n",
      "-0.795344022206\n",
      "Final objective function is: -1.21406175007\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 18 AT INTERATOMIC DISTANCE 2.0 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.41274611943\n",
      "objective function at theta+ for step # 0\n",
      "-0.7708230171\n",
      "objective function at theta- for step # 0\n",
      "-0.770055964926\n",
      "Final objective function is: -1.20929805853\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 19 AT INTERATOMIC DISTANCE 2.1 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 4.13277576587\n",
      "objective function at theta+ for step # 0\n",
      "-0.747356918779\n",
      "objective function at theta- for step # 0\n",
      "-0.757540461584\n",
      "Final objective function is: -1.17927507419\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 20 AT INTERATOMIC DISTANCE 2.2 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.62223550753\n",
      "objective function at theta+ for step # 0\n",
      "-0.736926995216\n",
      "objective function at theta- for step # 0\n",
      "-0.733305040755\n",
      "Final objective function is: -1.17827133761\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 21 AT INTERATOMIC DISTANCE 2.3 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.23246307832\n",
      "objective function at theta+ for step # 0\n",
      "-0.717016529988\n",
      "objective function at theta- for step # 0\n",
      "-0.720431849214\n",
      "Final objective function is: -1.16666728759\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 22 AT INTERATOMIC DISTANCE 2.4 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.91416605273\n",
      "objective function at theta+ for step # 0\n",
      "-0.703400953031\n",
      "objective function at theta- for step # 0\n",
      "-0.703441424565\n",
      "Final objective function is: -1.15181747651\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 23 AT INTERATOMIC DISTANCE 2.5 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.86736564913\n",
      "objective function at theta+ for step # 0\n",
      "-0.688432377068\n",
      "objective function at theta- for step # 0\n",
      "-0.689141036095\n",
      "Final objective function is: -1.14419628722\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 24 AT INTERATOMIC DISTANCE 2.6 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.26286828093\n",
      "objective function at theta+ for step # 0\n",
      "-0.675247137487\n",
      "objective function at theta- for step # 0\n",
      "-0.675094364082\n",
      "Final objective function is: -1.13596019713\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 25 AT INTERATOMIC DISTANCE 2.7 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.02301177432\n",
      "objective function at theta+ for step # 0\n",
      "-0.667471609944\n",
      "objective function at theta- for step # 0\n",
      "-0.660031295288\n",
      "Final objective function is: -1.12912804046\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 26 AT INTERATOMIC DISTANCE 2.8 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 4.01251601605\n",
      "objective function at theta+ for step # 0\n",
      "-0.648070325063\n",
      "objective function at theta- for step # 0\n",
      "-0.652094905423\n",
      "Final objective function is: -1.12310301446\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 27 AT INTERATOMIC DISTANCE 2.9 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.85817709969\n",
      "objective function at theta+ for step # 0\n",
      "-0.638517404808\n",
      "objective function at theta- for step # 0\n",
      "-0.639180392944\n",
      "Final objective function is: -1.11597163699\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 28 AT INTERATOMIC DISTANCE 3.0 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.36154883198\n",
      "objective function at theta+ for step # 0\n",
      "-0.629276321916\n",
      "objective function at theta- for step # 0\n",
      "-0.626974483015\n",
      "Final objective function is: -1.10967583843\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 29 AT INTERATOMIC DISTANCE 3.1 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.25176964204\n",
      "objective function at theta+ for step # 0\n",
      "-0.615439382399\n",
      "objective function at theta- for step # 0\n",
      "-0.62052612576\n",
      "Final objective function is: -1.10402988097\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 30 AT INTERATOMIC DISTANCE 3.2 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.83832469266\n",
      "objective function at theta+ for step # 0\n",
      "-0.611893972462\n",
      "objective function at theta- for step # 0\n",
      "-0.605140723078\n",
      "Final objective function is: -1.0987190963\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 31 AT INTERATOMIC DISTANCE 3.3 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.4291337116\n",
      "objective function at theta+ for step # 0\n",
      "-0.604163623749\n",
      "objective function at theta- for step # 0\n",
      "-0.595813826549\n",
      "Final objective function is: -1.09347842994\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 32 AT INTERATOMIC DISTANCE 3.4 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.51632119157\n",
      "objective function at theta+ for step # 0\n",
      "-0.588956870413\n",
      "objective function at theta- for step # 0\n",
      "-0.593335625752\n",
      "Final objective function is: -1.08881826203\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 33 AT INTERATOMIC DISTANCE 3.5 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.76839646695\n",
      "objective function at theta+ for step # 0\n",
      "-0.582251831886\n",
      "objective function at theta- for step # 0\n",
      "-0.584233208559\n",
      "Final objective function is: -1.08431038112\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 34 AT INTERATOMIC DISTANCE 3.6 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.56374749227\n",
      "objective function at theta+ for step # 0\n",
      "-0.569588461376\n",
      "objective function at theta- for step # 0\n",
      "-0.582082604603\n",
      "Final objective function is: -1.08016497774\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 35 AT INTERATOMIC DISTANCE 3.7 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.35452535399\n",
      "objective function at theta+ for step # 0\n",
      "-0.570542376934\n",
      "objective function at theta- for step # 0\n",
      "-0.56678356374\n",
      "Final objective function is: -1.07601671526\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 36 AT INTERATOMIC DISTANCE 3.8 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 3.15744908508\n",
      "objective function at theta+ for step # 0\n",
      "-0.563182641504\n",
      "objective function at theta- for step # 0\n",
      "-0.560666090809\n",
      "Final objective function is: -1.07239314989\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 37 AT INTERATOMIC DISTANCE 3.9 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.71137294313\n",
      "objective function at theta+ for step # 0\n",
      "-0.557238673045\n",
      "objective function at theta- for step # 0\n",
      "-0.553647165638\n",
      "Final objective function is: -1.06884523397\n",
      "\n",
      "OPTIMIZING HAMILTONIAN # 38 AT INTERATOMIC DISTANCE 4.0 ANGSTROM\n",
      "\n",
      "calibration step # 0 of 25\n",
      "calibration step # 5 of 25\n",
      "calibration step # 10 of 25\n",
      "calibration step # 15 of 25\n",
      "calibration step # 20 of 25\n",
      "calibrated SPSA_parameters[0] is 2.16430434131\n",
      "objective function at theta+ for step # 0\n",
      "-0.552862629432\n",
      "objective function at theta- for step # 0\n",
      "-0.54599750462\n",
      "Final objective function is: -1.06546477498\n"
     ]
    }
   ],
   "source": [
    "# MOLECULE PARAMETERS\n",
    "\n",
    "molecule='H2' # name of the molecule: options are H2 or LiH\n",
    "\n",
    "if molecule=='H2':\n",
    "    n=2 # qubits\n",
    "    Z1=1\n",
    "    Z2=1\n",
    "    min_distance=.2\n",
    "    max_distance=4\n",
    "    number_of_points=39\n",
    "elif molecule=='LiH':\n",
    "    n=4 # qubits\n",
    "    Z1=1\n",
    "    Z2=3\n",
    "    min_distance=.5\n",
    "    max_distance=5\n",
    "    number_of_points=46\n",
    "    \n",
    "# OPTIMIZATION PARAMETERS \n",
    "\n",
    "run_optimization=True # Choose True or False. False just plots the exact potential energy surface\n",
    "\n",
    "def cost_function(Q_program,H,n,m,entangler_map,shots,device,theta):\n",
    "    \n",
    "    return eval_hamiltonian(Q_program,H,trial_circuit_ryrz(n,m,theta,entangler_map,None,False),shots,device).real\n",
    "\n",
    "\n",
    "m=3 # depth (number of layers of sq gates - # entanglers is this number -1)\n",
    "device='local_qasm_simulator'\n",
    "initial_theta=np.random.randn(2*n*m)  # initial angles \n",
    "entangler_map = Q_program.get_backend_configuration(device)['coupling_map'] # the map of two-qubit gates with control at key and target at values\n",
    "if entangler_map == 'all-to-all':\n",
    "    entangler_map = {i: [j for j in range(n) if j != i] for i in range(n)}\n",
    "shots=1\n",
    "max_trials=200\n",
    "\n",
    "initial_c = 0.01\n",
    "target_update = 2*np.pi*0.1\n",
    "save_step = 250 #Setting larger than max_trials to suppress the output at every function call\n",
    "\n",
    "####################     COMPUTING THE POTENTIAL ENERGY SURFACE      ###################\n",
    "\n",
    "mol_distance=np.zeros(number_of_points)\n",
    "coulomb_repulsion=np.zeros(number_of_points)\n",
    "electr_energy=np.zeros(number_of_points)\n",
    "electr_energy_optimized=np.zeros(number_of_points)\n",
    "\n",
    "Q_program = QuantumProgram()\n",
    "Q_program.set_api(Qconfig.APItoken,Qconfig.config[\"url\"])\n",
    "\n",
    "for i in range(number_of_points):\n",
    "\n",
    "    # %%%%%%%%%%%%% Coulomb Repulsion For Diatomic Molecules %%%%%%%%%%%%%%%%%%%%%%\n",
    "        mol_distance[i]=np.around((min_distance+(max_distance-min_distance)*i/(number_of_points-1)),2)\n",
    "        \n",
    "        distance=mol_distance[i]/0.529177\n",
    "        coulomb_repulsion[i]=Z1*Z2/distance\n",
    "        \n",
    "    # exact diagonalization\n",
    "        ham_text=molecule+'/PESMap'+str(i)+'atdistance'+str(mol_distance[i])+'.txt'\n",
    "        pauli_list=Hamiltonian_from_file(ham_text)\n",
    "        H=make_Hamiltonian(pauli_list)\n",
    "        pauli_list_grouped = group_paulis(pauli_list)\n",
    "        eigen=la.eig(H)\n",
    "        electr_energy[i]=np.amin(eigen[0])\n",
    "           \n",
    "    # optimization\n",
    "    \n",
    "        if run_optimization:\n",
    "            print('\\nOPTIMIZING HAMILTONIAN # '+str(i)+' AT INTERATOMIC DISTANCE ' + str(mol_distance[i]) + ' ANGSTROM\\n')\n",
    "            \n",
    "            \n",
    "            if shots !=1:\n",
    "                H=group_paulis(pauli_list)\n",
    "         \n",
    "        \n",
    "            SPSA_params = SPSA_calibration(partial(cost_function,Q_program,H,n,m,entangler_map,\n",
    "                                           shots,device),initial_theta,initial_c,target_update,25)\n",
    "            \n",
    "            electr_energy_optimized[i] = SPSA_optimization(partial(cost_function,Q_program,H,n,m,entangler_map,shots,device),\n",
    "                                                           initial_theta,SPSA_params,max_trials,save_step,1)[0];\n",
    "\n",
    "    \n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x1112bfc50>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8FfW9//HX55zsZGPfMSg7KKgBpVbrAq219WJdq9dW\nbC2ttb1281avv4e23t5bettbb1vrz6L2h1p7i0td2qIWXGoXF1CRVQRUZN8EkhASkpzP74+Z4CGc\nrORkEvJ+Ph6HM8t3Zj5nQs4nM9+Zz5i7IyIi0lKxqAMQEZGuRYlDRERaRYlDRERaRYlDRERaRYlD\nRERaRYlDRERaRYlDRERaRYlDRERaRYlDRERaJSPKjZvZucDPgDhwj7vPbjD/W8A1QC2wA/iCu69v\nbr19+vTxkpKS9g9YROQo9dprr+10974taRtZ4jCzOPBLYDqwEVhkZk+6+8qkZm8Ape5eaWbXAv8F\nXNbcuktKSli8eHE6whYROSqZWbN/lNeL8lTVFGCtu7/j7geA3wEzkhu4+/PuXhmOvgwM6eAYRUSk\ngSgTx2BgQ9L4xnBaY74IPJXWiEREpFmR9nG0lJldCZQCH2uizSxgFsCwYcM6KDIRke4nysSxCRia\nND4knHYIM5sG3Ax8zN2rG1uZu88B5gCUlpaqVrzIUaCmpoaNGzdSVVUVdShHjZycHIYMGUJmZmab\n1xFl4lgEjDSz4QQJ47PAFckNzOxE4FfAue6+veNDFJEobdy4kYKCAkpKSjCzqMPp8tydXbt2sXHj\nRoYPH97m9UTWx+HutcDXgGeAVcBD7r7CzG4zs38Km/0YyAceNrMlZvZkROGKSASqqqro3bu3kkY7\nMTN69+59xEdwkfZxuPt8YH6DabckDU/r8KBEpFNR0mhf7bE/ded4yBMJ3pr3f3j35T9EHYqISKem\nxBGyWIwhK+9h++tKHCISiMfjTJo06eBr9uzZzS/UQkuWLGH+/PnNN+yEusTluB2l3PKJV++JOgwR\n6SRyc3NZsmRJWta9ZMkSFi9ezHnnnZeW9aeTjjiS7IsXkFmzN+owRKQT27t3L6NHj2b16tUAXH75\n5dx9990AXHvttZSWljJ+/HhuvfXWg8ssWrSIj3zkI0ycOJEpU6awd+9ebrnlFubNm8ekSZOYN29e\nJJ+lrXTEkaQqXkh2TVnUYYhICt//wwpWbm7f389xgwq59fzxjc7fv38/kyZNOjh+0003cdlll3HH\nHXcwc+ZMrr/+enbv3s2XvvQlAP7jP/6DXr16UVdXxznnnMPSpUsZM2YMl112GfPmzWPy5MmUlZWR\nl5fHbbfdxuLFi7njjjva9TN1BCWOJNWZhRTsfzfqMESkk2jsVNX06dN5+OGHue6663jzzTcPTn/o\noYeYM2cOtbW1bNmyhZUrV2JmDBw4kMmTJwNQWFjYYfGnixJHktqsYvL3lUcdhoik0NSRQUdLJBKs\nWrWKvLw8du/ezZAhQ3j33Xf5yU9+wqJFi+jZsyczZ848au94Vx9HkkROEQVeAa6KJSLSuNtvv52x\nY8fy29/+lquvvpqamhrKysro0aMHRUVFbNu2jaeeCmqyjh49mi1btrBo0SIAysvLqa2tpaCggPLy\nrvmHqo44knhOT7KslqrKcnJ6dP3DSRE5Mg37OM4991yuvvpq7rnnHl599VUKCgo444wz+MEPfsD3\nv/99TjzxRMaMGcPQoUM57bTTAMjKymLevHl8/etfZ//+/eTm5rJw4ULOOussZs+ezaRJkw72nXQV\nShxJYj16AlC+Z4cSh4hQV1eXcvqqVasODv/0pz89ODx37tyU7SdPnszLL7982PT6o5CuRqeqkmT2\n6AXAvj07I45ERKTzUuJIkpXfG4DKvUocIiKNUeJIklvUB4AD5bsijkREpPNS4kiSVxQccdTs+yDi\nSEREOi8ljiT5PfsCkNi3O+JIREQ6LyWOJAX5RRzwOOxX4hARaYwSRxKLxSi3fEwVckUktHHjRmbM\nmMHIkSM57rjjuP766zlw4ECj7ffs2cOdd955cHzz5s1cfPHFrdrmLbfcwsKFC9scc738/PwjXkcq\nShwNVFgBGUocIkLwjO4LL7yQCy64gDVr1vD2229TUVHBzTff3OgyDRPHoEGDeOSRR1q13dtuu41p\n0zrvA1B1A2ADlfECslQhV0SA5557jpycHK6++mogeLDT7bffzvDhwxk+fDjPPPMMe/fuZdOmTVx5\n5ZXceuut3Hjjjaxbt45JkyYxffp0rrvuOj796U+zfPly5s6dy+OPP86+fftYs2YN3/nOdzhw4AAP\nPPAA2dnZzJ8/n169ejFz5kw+/elPU1JSwjXXXAMENyMuX74cd2fdunVcd9117Nixg7y8PO6++27G\njBnDu+++yxVXXEFFRQUzZsxI235R4migOrOQ/AO6j0Ok03nqRti6rH3XOeB4+GTjT/VbsWIFJ598\n8iHTCgsLGTZsGLW1tbz66qssX76cvLw8Jk+ezKc+9Slmz57N8uXLD1bVfe+99w5Zfvny5bzxxhtU\nVVUxYsQIfvSjH/HGG2/wzW9+k/vvv59vfOMbB9uWlpYeXM8NN9zAueeeC8CsWbO46667GDlyJK+8\n8gpf/epXee6557j++uu59tpr+fznP88vf/nL9thDKSlxNFCTWUSPKpVWF5HmTZ8+nd69g8v4L7zw\nQv72t79xwQUXNLnMWWedRUFBAQUFBRQVFXH++ecDcPzxx7N06dKUy8ybN4/XX3+dP//5z1RUVPCP\nf/yDSy655OD86upqAP7+97/z6KOPAvC5z32O7373u0f8GVNR4migNruI/LKuWbFS5KjWxJFBuowb\nN+6w/omysjLef/99MjIyMLND5jUcTyU7O/vgcCwWOzgei8Wora09rP3y5cv53ve+x4svvkg8HieR\nSFBcXNzoI21bEsORirRz3MzONbPVZrbWzG5MMT/bzOaF818xs5J0x+Q5PSlgP3W1NenelIh0cuec\ncw6VlZXcf//9QNDP8O1vf5uZM2eSl5fHggUL+OCDD9i/fz+PP/44p512WruWS9+zZw+XX345999/\nP337BveZFRYWMnz4cB5++GEg6MCvf5jUaaedxu9+9zsAHnzwwXaJIZXIEoeZxYFfAp8ExgGXm9m4\nBs2+COx29xHA7cCP0h5XblAht0KFDkW6PTPjscce4+GHH2bkyJGMGjWKnJwc/vM//xOAKVOmcNFF\nF3HCCSdw0UUXUVpaSu/evTnttNOYMGECN9xwwxFt/4knnmD9+vV86UtfYtKkSQdLvD/44IPce++9\nTJw4kfHjx/PEE08A8LOf/Yxf/vKXHH/88WzatOnIPnwTzCN6aJGZTQW+5+6fCMdvAnD3Hya1eSZs\n85KZZQBbgb7eTNClpaW+ePHiNsX16hN3MeWN77Lpyr8yeMQJbVqHiLSPVatWMXbs2KjDSGnu3Lld\n9pnhqfarmb3m7qUtWT7KU1WDgQ1J4xvDaSnbuHstsBfonWplZjbLzBab2eIdO3a0Oais/OCIo3Kv\nCh2KiKRy1NwA6O5z3L3U3UvrzwW2RU5hUCG3ShVyRaQJM2fO7JJHG+0hysSxCRiaND4knJayTXiq\nqghI6zf6h6XV1cch0hlEdTr9aNUe+zPKxLEIGGlmw80sC/gs8GSDNk8CV4XDFwPPNde/caR6hImj\nThVyRSKXk5PDrl27lDzaibuza9cucnJyjmg9kd3H4e61ZvY14BkgDvza3VeY2W3AYnd/ErgXeMDM\n1gIfECSXtCqsL61eqcQhErUhQ4awceNGjqTfUg6Vk5PDkCFDjmgdkd4A6O7zgfkNpt2SNFwFXNJw\nuXTKysqi3HOxKiUOkahlZmYyfPjwqMOQBo6azvH2VGH5xKv3Rh2GiEinpMSRwr54IZkHlDhERFJR\n4khhf7yA7FqVVhcRSUWJI4UDmYXk1qnQoYhIKkocKdRmF5OfUOIQEUlFiSOFuuxiCrwcdO24iMhh\nlDhSyS0my+qoqtRRh4hIQ0ocKcTzegFQtls3HYmINKTEkUJGfpA49u1R4hARaUiJI4XsMHHsL1OF\nXBGRhpQ4UqgvrV6t0uoiIodR4kihR3GQOGorlDhERBpS4kghP0wcdaqQKyJyGCWOFPLzizjgcdiv\nxCEi0pASRwoWi1FmBcSq9kQdiohIp6PE0Yh9sXwyVCFXROQwShyNqIwXkFWjCrkiIg0pcTSiOqOI\nHJVWFxE5jBJHI2qyCslThVwRkcMocTSiLiuskCsiIodQ4mhEIqeYAvZTV1sTdSgiIp1KJInDzHqZ\n2QIzWxO+90zRZpKZvWRmK8xsqZld1pExxvKCkMp37+zIzYqIdHpRHXHcCDzr7iOBZ8PxhiqBz7v7\neOBc4H/MrLijAoz3CAodVqhCrojIIaJKHDOA+8Lh+4ALGjZw97fdfU04vBnYDvTtqAAzwwq5laqQ\nKyJyiKgSR3933xIObwX6N9XYzKYAWcC6dAdWL6egNwBVZTpVJSKSLCNdKzazhcCAFLNuTh5xdzez\nRh/ubWYDgQeAq9w90US7WcAsgGHDhrUp5mS5RUGhwxpVyBUROUTaEoe7T2tsnpltM7OB7r4lTAzb\nG2lXCPwJuNndX25me3OAOQClpaWNJqKWyi8OzorV7lOhQxGRZFGdqnoSuCocvgp4omEDM8sCHgPu\nd/dHOjA2AAp7BkccrtLqIiKHiCpxzAamm9kaYFo4jpmVmtk9YZtLgTOAmWa2JHxN6qgAMzOzKPdc\nrEqJQ0QkWdpOVTXF3XcB56SYvhi4Jhz+DfCbDg7tEOWxAmLVqpArIpJMd443oTJWQJZKq4uIHEKJ\nowlVGQVkq0KuiMghlDiaUJ1ZRG6dCh2KiCRT4mhCbVYR+SqtLiJyCCWOJiSyg9Lqnmj0vkMRkW5H\niaMpuT3JsjqqKiuijkREpNNQ4mhCrEdQWr1sT8ob20VEuiUljiZkhqXV9+1RoUMRkXpKHE3Iyq+v\nkKtChyIi9ZQ4mpBTFCSO6nIlDhGRekocTThYIVel1UVEDlLiaEJBmDjqVCFXROQgJY4m5PUo4IDH\nYb8Sh4hIPSWOJlgsRpkVEKvaE3UoIiKdhhJHM/bF8slQhVwRkYOUOJqxP15AVo0q5IqI1FPiaEZV\nRhE5Kq0uInKQEkczarKKyFOFXBGRg5Q4mlGXVUSBK3GIiNRT4miG5xZTwH5qaw5EHYqISKegxNEM\nyw0KHZbv0d3jIiKgxNGseFhavWLPjogjERHpHCJLHGbWy8wWmNma8L1nE20LzWyjmd3RkTHChxVy\nK/cqcYiIQLRHHDcCz7r7SODZcLwx/w682CFRNZBdGJZWL/8gis2LiHQ6LUocZvZ7M/uUmbVnopkB\n3BcO3wdc0Mi2Twb6A39ux223WF5YWl0VckVEAi1NBHcCVwBrzGy2mY1uh233d/ct4fBWguRwiDBR\n/TfwnXbYXpvkF4Wl1ffpiENEBCCjJY3cfSGw0MyKgMvD4Q3A3cBv3L0m1XJmthAYkGLWzQ3W72bm\nKdp9FZjv7hvNrMkYzWwWMAtg2LBhzXyilivs2QeARKUKHYqIQAsTB4CZ9QauBD4HvAE8CHwUuAo4\nM9Uy7j6tifVtM7OB7r7FzAYC21M0mwqcbmZfBfKBLDOrcPfD+kPcfQ4wB6C0tDRVEmqTjMwsyj0X\nq1JpdRERaGHiMLPHgNHAA8D5SaeY5pnZ4jZu+0mCpDM7fH+iYQN3/+ekGGYCpamSRrqVxwqIV+uI\nQ0QEWn7E8XN3fz7VDHcvbeO2ZwMPmdkXgfXApQBmVgp8xd2vaeN6211lrIDMAyp0KCICLU8cPc3s\nwgbT9gLL3D3VKaZmufsu4JwU0xcDhyUNd58LzG3Lto5UVUaBKuSKiIRamji+SNDfUH/UcSbwGjDc\nzG5z9wfSEFunUZ1ZREHluqjDEBHpFFqaODKBse6+DcDM+gP3A6cQ3Jh3VCeO2qxielRURB2GiEin\n0NL7OIbUJ43QdmCou38ApLwU92iSyCmi0MvxRCLqUEREItfSI44XzOyPwMPh+EXhtB7A0X+5UW5P\nsqyOysoy8vKLo45GRCRSLU0c1wEXEty3AcFpqkfd3YGz0hFYZxLPC0qrl+3eqcQhIt1es4nDzOLA\nQnc/C3g0/SF1Ppn5QeHefXt2wtAREUcjIhKtZvs43L0OSITlRrqlzLC0elXZzogjERGJXktPVVUA\ny8xsAbCvfqK7/0taoupkcguDQocHVCFXRKTFieP34atb6lEcHHHUVKhCrohIS6vj3mdmucAwd1+d\n5pg6ncKewRFHQqXVRURa/CCn84ElwNPh+CQzezKdgXUmuXkFHPAMvOrov/JYRKQ5Lb0B8HvAFMJ7\nNtx9CXBsmmLqdCwWo8zyiSlxiIi0OHHUuPveBtO61W3U+2IFZB5ouAtERLqflnaOrzCzK4C4mY0E\n/gX4R/rC6nwq4wVk1ShxiIi09Ijj68B4oBr4X6AM+Ea6guqMqjMKyaktjzoMEZHItfSqqkqC54Tf\n3Fzbo1VNVhE9qlRaXUSkpY+OHQV8ByhJXsbdz05PWJ1PXXYx+QmVVhcRaWkfx8PAXcA9QF36wum8\nPKeYAttP7YFqMrKyow5HRCQyLU0cte7+f9MaSScXywsKHZbv3UXPvoMijkZEJDot7Rz/g5l91cwG\nmlmv+ldaI+tk4j2Cj1uxZ0fEkYiIRKulRxxXhe83JE1zutFNgJn5QeLYt1cVckWke2vpVVXD23Oj\n4dHKPILO9veAS919d4p2wwj6VYYSJKrz3P299oylpXIKgkKH1WWqkCsi3VuTp6rM7F+Thi9pMO8/\nj2C7NwLPuvtI4NlwPJX7gR+7+1iCkifbj2CbRyS3qA+gCrkiIs31cXw2afimBvPOPYLtzgDuC4fv\nAy5o2MDMxgEZ7r4AwN0rwvtJIlFQHFTIrVWFXBHp5ppLHNbIcKrx1ujv7lvC4a1A/xRtRgF7zOz3\nZvaGmf04fIxtJArCZ3L4/sPOqImIdCvN9XF4I8Opxg9hZguBASlmHXL3ubu7maVaVwZwOnAi8D5B\nn8hM4N5GtjcLmAUwbNiwpkJrk4zMLHZTQGzftnZft4hIV9Jc4phoZmUERxe54TDheE5TC7r7tMbm\nmdk2Mxvo7lvMbCCp+y42Akvc/Z1wmceBU2kkcbj7HGAOQGlpaZNJra22Zg6loExlR0Ske2vyVJW7\nx9290N0L3D0jHK4fzzyC7T7Jh5f4XgU8kaLNIqDYzPqG42cDK49gm0esrGAEg2rewxPdqqK8iMgh\nWnoDYHubDUw3szXAtHAcMys1s3sA3L2OoD7Ws2a2jOAo5+6I4gXA+46hmAp2bd8YZRgiIpFq6Q2A\n7crddwHnpJi+GLgmaXwBcEIHhtak/CETYDVsXfMGfQa0fz+KiEhXENURR5c0YOSJAFRsXB5xJCIi\n0VHiaIXe/Yawh3xsx1tRhyIiEhkljlawWIzNmSUUlq+NOhQRkcgocbRSeeEIBtes15VVItJtKXG0\nVt8xFLKPnVvfjzoSEZFIKHG0Uo+hEwDYuvaNiCMREYmGEkcrDRwRXFm1T1dWiUg3pcTRSr37D2E3\nhcR0ZZWIdFNKHG2wOesYCstVs0pEuicljjaoKBzJYNWsEpFuSomjLfqOocD2s33zu1FHIiLS4ZQ4\n2iB/6PEAbFu7JOJIREQ6nhJHGwwKa1ZVbtKVVSLS/ShxtEHPvgPZRZGurBKRbkmJo422ZpVQVKEr\nq0Sk+1HiaKOKwhEMrnlfV1aJSLejxNFW/caSb/vZtlFHHSLSvShxtFHBsPDKqnW6skpEuhcljjYa\nHF5ZtV9XVolIN6PE0UZFvfuzk2LiO1dHHYqISIdS4jgCW7N1ZZWIdD9KHEegonAkQ2rWk6irizoU\nEZEOE1niMLNeZrbAzNaE7z0bafdfZrbCzFaZ2c/NzDo61sZYv7HkWTVbN+gZ5CLSfUR5xHEj8Ky7\njwSeDccPYWYfAU4DTgAmAJOBj3VkkE0pCq+s2vGOngYoIt1HlIljBnBfOHwfcEGKNg7kAFlANpAJ\nbOuQ6Fpg4KiTANi/cUXEkYiIdJwoE0d/d98SDm8F+jds4O4vAc8DW8LXM+6+quNCbFpRzz5spxfx\nXbqySkS6j4x0rtzMFgIDUsy6OXnE3d3MPMXyI4CxwJBw0gIzO93d/5qi7SxgFsCwYcOONPQW25Zd\nQrGurBKRbiSticPdpzU2z8y2mdlAd99iZgOB7SmafQZ42d0rwmWeAqYChyUOd58DzAEoLS09LAml\ny76ikRy37XESdXXE4vGO2qyISGSiPFX1JHBVOHwV8ESKNu8DHzOzDDPLJOgY7zSnqgBi/YMrq7as\nfzvqUEREOkSUiWM2MN3M1gDTwnHMrNTM7gnbPAKsA5YBbwJvuvsfogi2MYUHr6xSzSoR6R7Seqqq\nKe6+CzgnxfTFwDXhcB3w5Q4OrVUGqWaViHQzunP8CBUW92YbvcnUlVUi0k0ocbSDbTnDKd73TtRh\niIh0CCWOdlBZNJIhte9TV1sbdSgiImmnxNEO4v3HkmM1bHnvrahDERFJOyWOdqArq0SkO1HiaAeD\nRwVXVlVvXhlxJCIi6afE0Q7yC3uylb5kfKArq0Tk6KfE0U625ZTQc59qVonI0U+Jo53sLx7FkNqN\nurJKRI56ShztJNZ/LNlWw+Z39WwOETm6KXG0k+JjTgBgxztLI45ERCS9lDjayeCREwGo3qyaVSJy\ndFPiaCc9CopZGz+OPpuejToUEZG0iqw67tFo53Gf4dS3f8L6Va9xzNiTow5HOovaA/DuX2DVk1BT\nBQMnwqBJMOAEyClsn224w/7dUFsF8SyIZ4bvWRDTA8akfSlxtKOR51xNzerb2fzir5U4jjaJBLz9\nFFTugn7joO8YyM5vvH19sljxGLz1J6jaA9mFkJUPyx4KGxn0Pg4GTgoSycBJUDAQ6g5Aogbq6l8H\ngvdEDezfQ6JsC4nyLXjZVrx8C7GKbcT2bSNWV50yFCdGIpZJIpZB4uBJhvAhme4N2oITJ2FxEsQO\ne6+zOHVkBO+WQR1x6ohTGw47hnmCGHXBuycw6oiRwDx81W8bxzzYogGGk/LRnU7SdA/226Gf4hBG\n4uD68A/XfXhrS729VrbypHjq29cv8+Fwqm3bIcs3fG98e42rjBcx7ua/N7l8e1DiaEe9+w9hSY8p\nHLdlPnW1tcQztHu7vERd8OX/4o9hR4NaZD1LoN946DcW+o8LEsqeDbDycXjrj1C1F7KLYMx51I2d\nwe4Bp7HnQIz9uzfD5jfJ2v4mebuW03PNX8lf/kiLQ4oBFZ7Ldu/Jdi9mG0PZ5iew3XuynywyqCOL\nWjLrX1ZLFnVkUkuMxMH1HPZFZUYcJ8MSB19xEmRaIpxeR0aQNsgMU0YwrYYMryOLOsycBLEgWVmM\nYA0x3DLC9xgHv4rNSP5ax4LhenbI96eF/3qDKYcywBus9+Bw0voPXU/yV/GHickPtm38y/zwpOKH\nf66D6cMOGatvX7+11OtrbnuHqs1spyPYZuibrZ0lTvgs/V6+nqV/fZwTzro46nCkrepqYcXvg4Sx\n8+3gCOOie2HQibB9FWxfSWLbShJbVxB/+2nM6w4uuj+Wz5K8qbxQ/FFeqB3PtmXO3ldrcP9r0gby\ngKnhC/raXiZnb6B/ZiWxjCziGVnEMrPIyMwmIzOLjKxssjKzSeQUUJvbn3hOAdkZMbIzY2RnxBmc\nEePYjBiZGTGy48F7ZjxGZtzIisfICIcz4zEyYkZGLEY8buGwEY8ZZk3/pStST4mjnY0/81L2vvxv\n1Lz+IChxdD11tbDs4SBhfLCOmt5jWHf6L3itx+ms31DFpqV72bK3J1v3nsy28gnUJZxsDnCsbWGU\nbaDC8lmZcxIF8Tx65WUxokc2U3pk0bNHFr17ZFGcl0lhTiYFORkU5gbvBTmZ9MiK64tbugwljnaW\nnZPHkj6fYOKOP1C2ZxeFxb2jDkmasX/PDjatfpXyd1/jmHd/R6/qTayLDef2um/xp00n4ZtiwEqy\nMmIMKc5lYHEOU4/rw4CibAYU5TKwMIcBRTkMLMqhZ14WsZgSgBzdlDjSoOdHZpLz5O9ZuvA+plz8\nrajDkZDX1bJ9/Sp2rF1E9YalZH+wkv6Va+nruxgRtlnux/LzHjezuf9ZlPTN5z9696Ckdx4lfXow\noDBHSUEEJY60GDnpdNb/cSiFqx8GlDiiUr17IxuX/ZWytS+Tt/0NhlStpj9V9AdqPM77scGs6zGJ\nVb3HkztsIgNGTWbc4GF8T8lBpEmRJA4zuwT4HjAWmOLuixtpdy7wMyAO3OPuszssyCNgsRibSz7D\n1Hd+zoa1yxg64vioQzp6uQdXL5VvoXz7era8vZjEhsX02buMPomdHAcc8DjvxIfzRs9PwqATKR5+\nEsPGnMhx+fkcF3X8Il1QVEccy4ELgV811sDM4sAvgenARmCRmT3p7l3iaUnHnfMF6tb9go0v/Jqh\nI26POpyuL5GA1fNh02Io2wJlm/CyLSTKNhOvrQSgIHy97/1YmTOe/f1OpGjEKRx7wkcY07M40vBF\njiaRJA53XwU0dxXJFGCtu78Ttv0dMAPoEomj3+DhLM09meEb/0Ci7ifE4rp7t03cYc2f4dl/h23L\n8FgGldn92Oo9WVPVl421I9lGL7J7DWHoMSMYOe5Exo04lmGZ2t8i6dKZ+zgGAxuSxjcCp0QUS5sc\nmHAZAxbfwPKX5zPhtPOjDqfree9v8OxtsOEV9uYM5ld53+auD04kURmjf2E2Z4zvy8dG9+XiEX0o\nzsuKOlqRbiNticPMFgIDUsy62d2fSMP2ZgGzAIYNG9beq2+TCWdfQfmiW9j/6gOgxNFym15n39Pf\no8eGv7DTevHTmi/yUNXHmHRMX757Sn8+Nrovo/sX6L4HkYikLXG4+7QjXMUmYGjS+JBwWmPbmwPM\nASgtLW1JCZq0y8nLZ2nvaUzY9Wf2le+hR4HOszdl/YpXOPDcDxm563mqPZ//qftnVg6+lOknHMPf\nJgxkQFFO1CGKCJ37VNUiYKSZDSdIGJ8Frog2pNYrPOXz5D31BxY9+xsmX/C1qMPpdPZuWce65+bS\n650nKammJvreAAAP1UlEQVR7j3LP5eHCK6mdci1fmjiCfoVKFiKdTVSX434G+AXQF/iTmS1x90+Y\n2SCCy27Pc/daM/sa8AzB5bi/dvcu91zW0ZOnsfHpgeStfAiUOACoLdvOuhcfJGP5IxxXtZyTgBXx\nsfxj9I2MOmcml/QbGHWIItKEqK6qegx4LMX0zcB5SePzgfkdGFq7s1iMDcNmMHX9XWx+bzWDSkZH\nHVI06mrY+vJD7Hv1AY7Z+wqjSbCWoTw76CsMPeNKxo/RvS4iXUVnPlV11DjmrC/A3LtY//yvGXT1\nj6IOp0MdKNvBO0/fQb+3HmBAYhebvA9/Lr6MolMuZ/KU0xmRoYdQinQ1ShwdYFDJaFZkTWTo+0/g\niR9isaP/y3L72tfZ+uf/YdT2pxjDAV6JTWLH8bdyyvRLOa8wN+rwROQIKHF0kMpxlzJ+yc28tWgh\nY075eNThpEWiro4Vf3mI+KtzGFf1OgWexcsFHyf39K9SOvkjxFUDSuSooMTRQcadcyX73vh3Egu/\nz4FJZ5CVffRcLVS2Zycr/3QnQ9Y+yPG+la305i/DvsrIc6/jzEFDog5PRNqZEkcH6VFQzOLS2yh9\n7V9ZdOfnKb3+d13+lNX6t15j64JfcPzO+Zxq1byVOY5FE29g4vTP8bHs7KjDE5E0UeLoQKXnf5mX\ndr7D1PV38dJ9NzG1C3aU19XWsvyFh4gvmsOE6jcY4Jks7TmNnmd9jTETPxp1eCLSAZQ4OtipV/2Q\nRT97l6nr72LxH46l9PwvRx1Si5Tt2s6qp+5k6LrfMtG3sY3evFRyHaPOu47J/QZHHZ6IdCAljg5m\nsRgnXDuXFT/9OCcs/jdW9Sth7CmfiDqsRq1982/sfv5Ojt/9Z06xGlZkTmDLSTdxwrR/pn+mCguK\ndEdKHBHIzsljyJcfZdsdZzLgqS+yofdTnephT9VVlSxfcD893pzLmNpVVHo2S3qfR5+zr2P8hC5V\noFhE0kCJIyJFvftT8blH4P5PwIOXsue65ynuk6qYcMfZumEt7z59B6M2/Z6T2csGG8RLI29g3Ce/\nwqm9+kQam4h0HkocERp87Hje+uSvOXb+5az91UXkfnsB2Tl5HRpDZcVeVjz7IDkrH2J81RL6AUvz\nTmXTqbOY8NEZDNUDqESkASWOiI055eMs3vFDShffwOI7P8/J33go7ZfpJurqWPXyU1QueoDxu59n\nslWz2frxyrAvcszZX2LS8DFp3b6IdG1KHJ1A6adn8dKOdUxdfxdv/vhcvPQLTDjjQjLasfPZEwne\nf/sNtvztAUo2/pHx7KDCc1neaxr5Uz7HmCkfZ5COLkSkBcy9UzzzqF2Vlpb64sWLow6jVTyR4OX7\nbmL0+t/SizJ2UszaAZ9i4JnXcMyYk9q0zq3rV7Phtaex9/7KMWWL6ctu6txYnnsyNRMuY/xZl5Pb\no6CdP4mIdEVm9pq7l7aorRJH53KguooVf3kEljzI8fteJsMSrM4Yw57RlzJ2+kwKi3sf0r7mQBVV\nlRVU76+gqqKMbatfoe6dvzB496sM9m0A7KSY9wpOou6YMxg+9TP0G1wSwScTkc5MiaMLJ45kO7du\nYO3CexnwziOUJDaw37P4INaTbK8m26vJ4QCZVnfYcmXeg7U9JlE95KP0nzSd4WNO7vLlTUQkvZQ4\njpLEUc8TCdYseZHdLz1AvHovdRm5ePgiMxfLysMy87CsPIqGjuO4E04jIzMz6rBFpAtpTeJQ53gX\nYLEYo046E046M+pQRETQ+QsREWkVJQ4REWkVJQ4REWkVJQ4REWmVSBKHmV1iZivMLGFmKXvxzWyo\nmT1vZivDttd3dJwiInK4qI44lgMXAi820aYW+La7jwNOBa4zs3EdEZyIiDQukstx3X0VgJk11WYL\nsCUcLjezVcBgYGVHxCgiIql1iT4OMysBTgReiTYSERFJ2xGHmS0EUj2Z6GZ3f6IV68kHHgW+4e5l\nTbSbBcwKRyvMbHVr4u1AfYCdUQfRBMV3ZBTfkVF8R+ZI4jumpQ0jLTliZi8A33H3lPVBzCwT+CPw\njLv/tCNjSxczW9zS2/qjoPiOjOI7MorvyHRUfJ32VJUFHSD3AquOlqQhInI0iOpy3M+Y2UZgKvAn\nM3smnD7IzOaHzU4DPgecbWZLwtd5UcQrIiIfiuqqqseAx1JM3wycFw7/DWj8squua07UATRD8R0Z\nxXdkFN+R6ZD4jsqy6iIikj6dto9DREQ6JyWONDGzc81stZmtNbMbU8yfaWY7kvpvrunA2H5tZtvN\nbHkj883Mfh7GvtTM2vbQ8/TFd6aZ7U3ad7d0cHzNlsOJch+2ML7I9qGZ5ZjZq2b2Zhjf91O0yTaz\neeH+eyW8l6szxRfZ729SDHEze8PM/phiXnr3n7vr1c4vIA6sA44FsoA3gXEN2swE7ogovjOAk4Dl\njcw/D3iKoI/pVOCVThbfmcAfI/z5DgROCocLgLdT/Hwj24ctjC+yfRjuk/xwOJPgxt5TG7T5KnBX\nOPxZYF4niy+y39+kGL4F/DbVzzHd+09HHOkxBVjr7u+4+wHgd8CMiGM6yN1fBD5ooskM4H4PvAwU\nm9nAjomuRfFFyt23uPvr4XA5UF8OJ1lk+7CF8UUm3CcV4Whm+GrY2ToDuC8cfgQ4x5qqUdTx8UXK\nzIYAnwLuaaRJWvefEkd6DAY2JI1vJPUv7kXhaYxHzGxox4TWIi2NP0pTw1MJT5nZ+KiCaKIcTqfY\nh82U64lsH4anWZYA24EF7t7o/nP3WmAv0LsTxQfR/v7+D/CvQKKR+Wndf0oc0fkDUOLuJwAL+PCv\nA2ne68Ax7j4R+AXweBRBtLQcTlSaiS/Sfejude4+CRgCTDGzCR25/ea0IL7Ifn/N7NPAdnd/raO2\n2ZASR3psApL/AhkSTjvI3Xe5e3U4eg9wcgfF1hLNxh8ldy+rP5Xg7vOBTDPr05ExWFAO51HgQXf/\nfYomke7D5uLrDPsw3PYe4Hng3AazDu4/M8sAioBdHRtd4/FF/Pt7GvBPZvYewWnws83sNw3apHX/\nKXGkxyJgpJkNN7Msgs6pJ5MbNDjf/U8E56E7iyeBz4dXBp0K7PWgzH2nYGYD6s/XmtkUgv/HHfal\nEm67uXI4ke3DlsQX5T40s75mVhwO5wLTgbcaNHsSuCocvhh4zsOe3s4QX5S/v+5+k7sPcfcSgu+W\n59z9ygbN0rr/Irlz/Gjn7rVm9jXgGYIrrH7t7ivM7DZgsbs/CfyLmf0TwQOrPiC4SqNDmNn/ElxV\n08eC0i+3EnQA4u53AfMJrgpaC1QCV3dUbC2M72LgWjOrBfYDn+2oL5VQfTmcZeF5cIB/A4YlxRjl\nPmxJfFHuw4HAfWYWJ0hYD7n7Hxv8ftwLPGBmawl+Pz7bQbG1NL7Ifn8b05H7T3eOi4hIq+hUlYiI\ntIoSh4iItIoSh4iItIoSh4iItIoSh4iItIoSh6SVmV1gZm5mY5KmlZjZFe24jdvMbFobl51rZheH\nw/eY2bgm2s40s0FtjfNImdk3zOzzSeMZYYXW2WncXl6a1t3XzJ5Ox7ol/ZQ4JN0uB/4WvtcrAdot\ncbj7Le6+sB3Wc427r2yiyUwgksQR3v37BYJqqPWmE1S+vSRNBQC/AaRMHOE9Dm3m7juALWZ22pGs\nR6KhxCFpE9ZK+ijwRQ69AWk2cLoFzzH4pgXPP/h/ZrbMgucLnBUuP9PMHjezBWb2npl9zcy+FbZ5\n2cx6he2Sjxomm9k/wuJ9r5pZQYOYzMzusOBZKQuBfknzXjCzUgsK3M01s+VhTN8M118KPBjGnWtm\nt5jZorDdnKQ7sV8wsx+F23/bzE4Pp8fN7Cdh+6Vm9vVw+slm9hcze83MnrHUVXTPBl4PC9bVuxz4\nGfA+MDXpc7xnZt83s9fD+MeE0/uG+3JFeHS13sz6mFkPM/tTuM+Wm9llZvYvBEnyeTN7Ply+wsz+\n28zeJCiQeE74s1hmwTNUspO2/8NwPy02s5PCz7XOzL6SFP/jwD83899IOqP2rNGul17JL4IvhXvD\n4X8AJ4fDZ5L0DAHg2wR31wOMIfgizCH4C38twTMl+hJU+PxK2O52guJ9AHMJ7oTOAt4BJofTC4GM\nBjFdSFCULk7wxbgHuDic9wJBcjiZoCJq/TLFyfOTpvdKGn4AOD+p3X+Hw+cBC8PhawlKXGfUL09w\nR/w/gL7htMvq90WDuL8PfD1pPAfYDOQCs4BfJM17r74twXMZ7gmH7wBuCofPJSgV3ge4CLg7afmi\npPX0SZruwKVJ298AjArH70/6ebwHXJv0c1qa9DPclrS+wcCyqP+f6tX6l444JJ0uJyjCRvh+eSPt\nPgr8BsDd3wLWA6PCec+7e7kHpzb2ElQlBVhGcMor2Whgi7svCtdV5of+hQ7BQ6L+14Pqp5uB51LE\n8w5wrJn9wszOBRqrfHuWBU9XW0ZwRJBcmry+sOBrSXFOA35VH5O7fxDGPAFYYEF5kP9DUBCxoYHA\njqTxTxPsm/0ExQwvaHD6KNX2P0r483D3p4Hd4fRlwPTwKOl0d9/byOetC7dFGPe77v52OH4fwb6t\nV1+bbRnBQ6zqf4bVFtaBIihZHlmfkbSdalVJWoSnkc4GjjczJ/gL383shlauqjppOJE0niBN/3/d\nfbeZTQQ+AXwFuJSgf+EgM8sB7iQ4AtlgZt8j+Cu8Ydx1zcRpwAp3n9pEGwjqSSWv/3LgoxZUSIXg\nWQtnExxNtWb7uPvbFjza9jzgB2b2rLvflqJplbvXNRNnveSfU8OfYX08OQSfS7oYHXFIulwMPODu\nx7h7ibsPBd4FTgfKCU5d1Psr4bluMxtFUIxvdRu2uRoYaGaTw3UVhJ3KyV4ELgv7GwYCZzVciQXl\nxWPu/ijBEUD988KT467/Et8Z9uVc3IL4FgBfro8pTK6rgb5mNjWclmmpH6q0ChgRtikk2I/Dwn1b\nAlxH40d09f5OkAQxs48DPcPhQUClu/8G+HEjn7eh1UCJmY0Ixz8H/KWZ7Tc0Ckj5XHnp3JQ4JF0u\nBx5rMO3RcPpSoC7sjP0mwV/usfCUzzxgpn/4rIMW8+AxvZcBvwg7cBdw6F/phDGtAVYSnJd/KcWq\nBgMvhKeOfgPcFE6fC9wVTq8G7ib44nuGoJR+c+4h6L9ZGsZ3RRjzxcCPwmlLgI+kWPYpPjwV9BmC\nMtnJ++gJ4Pz6DupGfB/4uJktBy4BthIkh+OBV8PPdSvwg7D9HODp+s7xZO5eRVDx9+Hw55YA7mpu\nBzRwFvCnVi4jnYCq44p0EWb2GPCv7r6mjctnA3UelP2fCvxfD55yFwkzexGY4e67m20snYr6OES6\njhsJOsnblDgITgE+ZGYx4ADwpfYKrLXMrC/wUyWNrklHHCIi0irq4xARkVZR4hARkVZR4hARkVZR\n4hARkVZR4hARkVZR4hARkVb5/4U1ULle7Ob0AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x111715550>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(mol_distance,electr_energy+coulomb_repulsion,label='Exact')\n",
    "if run_optimization:\n",
    "    plt.plot(mol_distance,electr_energy_optimized+coulomb_repulsion,label='Optimized') \n",
    "plt.xlabel('Atomic distance (Angstrom)')\n",
    "plt.ylabel('Energy')\n",
    "plt.legend()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<h2>Version information</h2>\n",
       "<p>Please note that this tutorial is targeted to the <b>stable</b> version of the QISKit SDK. The following versions of the packages are recommended:</p>\n",
       "<table>\n",
       "<tr><th>Package</th><th colspan=\"2\">Version</th></tr>\n",
       "<tr><td>QISKit</td><td> 0.4.8</td></tr>\n",
       "<tr><td>IBMQuantumExperience</td><td>&gt;= 1.8.26</td></tr>\n",
       "<tr><td>numpy</td><td>&gt;= 1.13, &lt; 1.14</td></tr>\n",
       "<tr><td>scipy</td><td>&gt;= 0.19, &lt; 0.20</td></tr>\n",
       "<tr><td>matplotlib</td><td>&gt;= 2.0, &lt; 2.1</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%run \"../version.ipynb\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "celltoolbar": "Slideshow",
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}