vqe_experiment.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from qiskit.aqua.operators import WeightedPauliOperator\n",
    "from qiskit.quantum_info import Pauli\n",
    "\n",
    "def convert_pauli_str(gates):\n",
    "    pauli_string = list('IIII')\n",
    "\n",
    "    if gates == '':\n",
    "        return ''.join(pauli_string)\n",
    "    \n",
    "    else:\n",
    "\n",
    "        for pauli in gates.split(' '):\n",
    "            gate = pauli[0]\n",
    "            qubit = int(pauli[1])\n",
    "\n",
    "            pauli_string[qubit] = gate\n",
    "\n",
    "        pauli_string = ''.join(pauli_string)\n",
    "    \n",
    "\n",
    "    return pauli_string\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "-0.3276081896748089 IIII\n",
      "-0.04919764587136759 XXYY\n",
      "0.04919764587136759 XYYX\n",
      "0.04919764587136759 YXXY\n",
      "-0.04919764587136759 YYXX\n",
      "0.13716572937099494 ZIII\n",
      "0.15660062488237958 ZZII\n",
      "0.10622904490856085 ZIZI\n",
      "0.15542669077992843 ZIIZ\n",
      "0.13716572937099494 IZII\n",
      "0.15542669077992843 IZZI\n",
      "0.10622904490856085 IZIZ\n",
      "-0.1303629205710914 IIZI\n",
      "0.1632676867356435 IIZZ\n",
      "-0.13036292057109133 IIIZ\n"
     ]
    }
   ],
   "source": [
    "H2 = ['-0.3276081896748089 []  ', '-0.04919764587136759 [X0 X1 Y2 Y3]  ','0.04919764587136759 [X0 Y1 Y2 X3]  ','0.04919764587136759 [Y0 X1 X2 Y3]  ','-0.04919764587136759 [Y0 Y1 X2 X3]  ','0.13716572937099494 [Z0]  ', '0.15660062488237958 [Z0 Z1]  ','0.10622904490856085 [Z0 Z2]  ', '0.15542669077992843 [Z0 Z3]  ', '0.13716572937099494 [Z1]  ', '0.15542669077992843 [Z1 Z2]  ','0.10622904490856085 [Z1 Z3]  ', '-0.1303629205710914 [Z2]  ', '0.1632676867356435 [Z2 Z3]  ', '-0.13036292057109133 [Z3]  ' ]\n",
    "op_list = []\n",
    "for weighted_pauli in H2:\n",
    "    line = weighted_pauli.replace(' ', '').replace(']','').split('[')    \n",
    "    weight  = float(line[0])\n",
    "    pauli_str = line[1][:-1]   \n",
    "    \n",
    "    final_str = convert_pauli_str(pauli_str)\n",
    "    \n",
    "    print(weight,final_str)\n",
    "    \n",
    "    op_list.append((weight, Pauli(label = final_str)))\n",
    "    \n",
    "    \n",
    "\n",
    "op = WeightedPauliOperator(op_list)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Result generated from Codi\n",
    "\n",
    "H2 \n",
    "\n",
    "-0.3276081896748089 []  \n",
    "-0.04919764587136759 [X0 X1 Y2 Y3]  \n",
    "0.04919764587136759 [X0 Y1 Y2 X3]  \n",
    "0.04919764587136759 [Y0 X1 X2 Y3]  \n",
    "-0.04919764587136759 [Y0 Y1 X2 X3]  \n",
    "0.13716572937099494 [Z0]  \n",
    "0.15660062488237958 [Z0 Z1]  \n",
    "0.10622904490856085 [Z0 Z2]  \n",
    "0.15542669077992843 [Z0 Z3]  \n",
    "0.13716572937099494 [Z1]  \n",
    "0.15542669077992843 [Z1 Z2]  \n",
    "0.10622904490856085 [Z1 Z3]  \n",
    "-0.1303629205710914 [Z2]  \n",
    "0.1632676867356435 [Z2 Z3]  \n",
    "-0.13036292057109133 [Z3]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAz4AAADWCAYAAAAQNFeqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1 /AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy j8jraAAAgAElEQVR4nOzde1yUdfr/8dfMcD7KQUVBFEVIUDx1Eg oWWm1qZW2Zm4Hd3PFrTQtW612yzJ/lq61u2p2sm8HbXWtNU9lpnhOzCOeUFEBxTOoIIrMzO PWxAUhhmYmXvu4Xo Hj6Em3tu34zXfbjmvu/PrTObzWaEEEIIIYQQwo3p1Q4ghBBCCCGEEI4mjY8QQgghhBDC7UnjI4QQQgghhHB70vgIIYQQQggh3J40PkIIIYQQQgi3J42PEEIIIYQQwu1J4yOEEEIIIYRwe9L4CCGEEEIIIdyeND5CCCGEEEIItyeNjxBCCCGEEMLtSeMjhBBCCCGEcHvS AghhBBCCCHcnjQ QgghhBBCCLcnjY8QQgghhBDC7UnjI4QQQgghhHB70vgIIYQQQggh3J40PkIIIYQQQgi3J42PEEIIIYQQwu1J4yOEEEIIIYRwe9L4CCGEEEIIIdyeND5CCCGEEEIItyeNjxBCCCGEEMLtSeMjhBBCCCGEcHvS AghhBBCCCHcnjQ QgghhBBCCLcnjY8QQgghhBDC7UnjI4QQQgghhHB7HmoHEM7x8/vq/Lt9xtXt9WrkrmtmrdLie63VuhbOI3VtPalrbdBqfWhxXRTuR874CCGEEEIIIdyeND5CCCGEEEIItyeNjxBCCCGEEMLtSeMjhBBCCCGEcHsyuIGo5PFJTSm8nI9eb0CvN5AQ3YUxgz6mUYNmakezSKu5tUar77NWcwvn0Gp9aDW3cA6t1odWcwttkDM otzZC8c5fzGPmaO38cM7hXz5ahYFRaf5ZOl4taNZpNXcWqPV91mruYVzaLU tJpbOIdW60OruYV2SOMjyh3IScfPJ4joRrcBEOQfRmxkJwoKT6uczDKt5tYarb7PWs0tnEOr9aHV3MI5tFofWs0ttEMaH1HuQE468VF3oNPpMBpLSd /grSd39K74xNqR7NIq7m1Rqvvs1ZzC fQan1oNbdwDq3Wh1ZzC 2Qe3xEuQM56WQcXc A1xtQXFKIj6cfYwd9Skr7wWpHs0iruW9mNit/63Tq5qiOVt9nreZ2F1LXjqHV3O7CbHbdmgbt1odWcwvtkManApPJxPTp0/noo4/IyckhPj6eDz/8kOeee46UlBTmzJmjdkSHOpi7lQlPzKNbu4Gcu5jHKx/dQ/bp/WrHqpFWc5fZexzWHoCDp5SdaWQIdIuD22PA4ELnZLX6Pms1t5YZTbD1CKw7ACcKlAPE1o0h5TZo01TtdJVptT60mlvLiktg40Hlz/ki8PSA9s2Uuo4KVTtdZVqtD63mFtrhQodV6hs fDiTJk1ixIgRLF  nMGDBzNkyBCysrLo3Lmz2vEc6vjZQ1wqzqd1VCcAwoKa8FjKWJZsno3JZAIg48gGpswbVv6ad78Zyv7sLarkLWNNbqPJyMsf3cPYWT0ZMb0Df5rWTs3I5cxm GE7zFkDB04qB4smM Tmw7zN8NlaZZorkPoQ1io1widpSg0fz1dq2mhSavyj1bB0h9oJb5C6Fta6VAzTV8APO BcEZiBklKlwZ  ArYdVTvhDVLXQlRPGp/r5s2bx9y5c1m8eDHjxo2jV69eTJw4kS5dulBaWkqnTp3UjuhQB3LSCfQLpXFI8/JpyYn9Kbh0ioyj6wGIi pM1omdAOzOWoenhze3Rd pSt4y1uQ26A28N2IVE4bOIyK0BW89s1ituJXszIFVe5Wvyy4Hqvj1nuPw427n56qK1Iew1vJdsO E8nWFsi6v65V7YFeO02NVSepaWOvLjXD20q3TzSjN/Vcb4cxFp8eqktS1ENWTxue6yZMn07dvX1JSUipNj42NxdPTk6SkJACOHj1KSkoKcXFxtGvXjnXr1qkR1 4yc9KJbdqx0rRg/3ASWiSzbtdCALw8ffD1DuBi0Tk //E1hvd7V42olViTGyDv/BH uSiVFx ZTZPQGGfHrFLavpqvEV fCdeMzsljidSHsEZJKWw4aHkenQ7SXOTKFalrYY2TFyDzZOVG/mYmc8217yxS10JUT2c2my2ty/VCbm4uzZo149NPP XZZ5 t9LMhQ4awf/9 tm/fDsD9999P//79SU1NZePGjQwaNIgjR47g5eVllyw6B90tufI9 /w3z178EqcLsmnTvAuDUsbWOP 9L9ft97FH7mOn9vLFj28w5rGPCfQLqXH uma2hqe3P6mfFlo173/e7ErewY0OTmSf97o 1Ac4p0a0qGl8dwa9vtaqef/1jA/Ga1cdnEjqWuq67jrc/yIpw2bUON/5Ewf48pXbHJ5H9ueutT8XrsHadkbO KA0PgARERGVphcXF5OWllZ mdvZs2dZv349w4cPByA5OZmmTZuyevVq5wZWUUKLZI6ezGBgtxfUjmK1Vz  j1P5x/j7FwMZO6snFy fVzsSeoOn9fN6WD v2qQ 6jeDh/UfANkyr9qkrus3a2vV4Kmdmgapa1E/yahuQHh4OACZmZk88MAD5dOnTp1KXl5e cAG2dnZNG7cGG9v7/J5YmJiOHbsmN2yOOoE3M/v22c5Ow vIbX/h3hYeeBe19/HHrnnvZZr0/zOOAlqMsPfFsGlK5bn0 lgz9Y1BPs5PJJd3uv6UB/gnBrRooLL8OZ3li8JAgj2hStFF50yHLDUtfWkrqu25zh8vMbyPDodJHeI4XMnvIeyP7ee1LS4mTQ QMuWLUlKSmLy5MmEhoYSGRnJwoULWbZsGYDbj hmjTMFufzr 78Q3TiB2 PvUzuO5ul10LU1rLAweIEOaBeFU5qeupL6EAAN/CAhUhmi3dLhRtc4134GShmpawHQponSrF8srr6uzWZlm64FUteiPpPGB9Dr9SxYsIARI0YwcuRIwsLCeOqppxg1ahQTJkwoH9ggOjqaU6dOcfXq1fKzPkeOHKF58 aWFu8WGjaI4s2nv1c7hlvp2QZ2ZCs3zt5MpwNfT/hdx1t/5oqkPkSZ/p0g6wxcuVZ5tMIyTRpAj3jn56oNqWsBoNfD7   ftbHXHXz0zEabnOxZ1RVR pa1Gdyj891cXFxrF69mqKiIrKzs5k0aRK7d 8mISEBX19fQLkkrmvXrnz66acAbNy4kePHj9OrVy81owuN8vGE5  Fzi2UM0AVxTaC0fdDw0BVoglRa42ClNqNbVR5ul4Ht7eA5/sotS ElrRpCiN6QePgytO9PaBPIjzZ9dbtuBDC9cgZHwu2bt3K3XffXWna7Nmzefrpp5kxYwZeXl7MmzfPbiO6ifrH3xuGdYUBneD1Rcq0ib DhkHq5hKiLhoHwag cPoiTP5BmfbmIxDoo24uIeoivgmMfxByzisPLQV461Gl RFCaIOsrtUoLCwkMzOT1NTUStNbtmzJ2rXWDdeqNZ8u yt7jm4gsUVXohrGM3/1u4x dA7tW6XwnzXvsXHP/2gc0pyXH5/LtdKrvDKnD5Fhsbz6xFcumRlg3e5FzPrfi3zzWg7FVwtdInNVAn1vfO2qTU9173W7mO5Mnf8HThdk42HwYuKT8/Hy8FH1va4ua2KLZF6a2YMjJ3cze8wOIsNjyTt/hKnz/4AOHeHBUYwf8iUGvYHXPnuIwuICZoxa7/T87qJRhVp2pabHlvoA C1zJfN/eReT2cSI302jWcN4Veq7utyNQqKrrOFV275m8cZ/E gXyl f AZ/nyCp6zrS6SA67Mb3rtT02FofAAdzt5H6QWdWTLmGweDh9PqwNfPHS14h4 gG9Ho9Ywd9RlTD1lLTwiZyqVs1AgICMBqNPP/882pHcYojJzMounKR6alruXj5HFdKihiU8jLtW6WQX3iaHYdXM2PUemKaJLEh43t8vQOYOHS y2Yus27XQho2aAbgEpm1ytJ7ffjEDjw8vJieupb773iGVdu VvW9tpTVoPfgzae/p3u7x8rnD/BpwNvPLGF66loiQmPYsl8Z1OTtZ5eokl84lq31cfVaMUs3f8SU51YybeQa4qI6q1LflnJXVcOlxmss2Tyb6SPX0qfTMJZu/giQunZXttZHmcWbZtI6slP5986sD1szX7x8nszcrXzwlw0M7zeFHzbNdHpmoX3S AgAMo6s5/Y4ZXSXTq3vRX/90yCAzJyttG/Z8/rP rDv2CYVEt7KUmaAX/cto1PrPuh0UuZ1Zem9Dg OxGQyAlBYXECQX1iVy3AWS1l1Oh0hgY0rzR/oF4K/r3LhvofBE72uch0J92Jrfew9tgmdTs ET/oxZd4wikuKnJq3jKXcVdXw8bMHiYloh8HgQafWfdjrIttt4Ri21gfA0ZN7aBgcha 3OjeT2prZ1ysAf59gjCYjRVfU39cIbZIjQgHApcvn eKnvzF2Vk  WfUOlyo8FKzoSgF Pso1K/4 wRReKVArZiWWMgOs/O0L7un0pErp3Iul9zrIP5yr14p59r02LNk0i27tHlExac11UZ2zF07wW bK8h2xcE 21kf pVOcv5TH5D8uJ7F5Mks3feSkpJVZk7tiDRcWV95uFxW7xnZbOIat9QGwaN0M nf9i7OjlrM1s6eHFxGhMTw7NZ5/ffcX t45XIXUQutc6OpUoaZAv1Ceuv8tkhMfZvPeJZy5cOMhYf4 wZwpUL6/fOUiAT4N1IpZiaXM2w/9QkLzLnhq6OnwrszSe/1b5k8E zfks5f3sXbXQhakvc we99wyazVKSm9ynvfPsVLgz7GYJDNojuztT78fYJp26IbBr2BDrG9WZBmp6dH2qim3DfXsL9PMJevXASg6OpF/H1dY7stHMPW sg9cxA/nyCC/cNVSmx75mOn9nH87EE fyWTg8e38fmKiYwb/JlK6YVWyRkfAUDbmG7szlIGbdh5eE35pUsAcc3uYFdWGgDbDv5Mm Z3V7kMZ7OU ejJDDbtWcxfP 7LsVN7 HzFa2rFdAuW3muz2UyQXygAwf7hFF2p4sFETmQpa3VmLHyOh5NH0bxxgqPjCZXZWh/xze4g /Q AA6f2EFEaIzDM1alptw313BUwziOnszAaDKy/eDPtIl2je22cAxb6 PIyd1k5qTz14/7ciRvFzMW/dnlM4MZf98G6PV6ZV9TrO6 RmiTfLQpAIiJaIuHwZOxs3qS0LwLPl7 GE2lAIQENKJdyx6M/nc3GjWI5pHuo1VOq7CUeWC3FxjY7QUARv 7G8/0fVvNqJpn6b2 Pe4 fkz/jLGzemI2mxg3 HOXzQow6cvBZBxdz/GzB3m85ys0CGjE oxFnMo/xqJ1MxjY7UW6tRuo4m8gHMnW khu25 klim8NLMH3p5 /HXoNy6Xe /RTVXWcL 7/sRLM7sT4BvChCfUyS2cw9b66N7uEbpfvyx57KyejH5ktstn7tZuIH7egYyZ2R2jsZTU/h84PbPQPml8RLnhD7xb/vXaXQuZv3oKkeGtad8qhd/3Gs/ve40v/3nx1UKmzHuS GZ3qBG1nKXMZcqGuHSVzFpl6b1 4w8LK82r9nttKevrw/5zy/yL3750y7TXPnuI0KAmDs0p1GFrfTzaYwyP9hhT/r1a9W0pd1U1fG/nYdzbeVilaVLX7svW igzbeSa8q dXR 2Zn7hkZm3TJOaFrbQmc1ms9ohhOP9rM5l6fQZV7fXq5G7rplra/TXyt8zhqrz72vxvdZqXdcnUte2k7p2fWrWtVbrQ4vronA/0vgI4SLUPkAUwhGkroU7kroWQptkcAMhhBBCCCGE25PGRwghhBBCCOH2pPERQgghhBBCuD1pfIQQQgghhBBuTxofIYQQQgghhNuTxkcIIYQQQgjh9qTxEUIIIYQQQrg9aXyEEEIIIYQQbk8aHyGEEEIIIYTbk8ZHCCGEEEII4fak8RFCCCGEEEK4PWl8hBBCCCGEEG5PGh8hhBBCCCGE25PGRwghhBBCCOH2pPERQgghhBBCuD0PtQMI5/j5fXX 3T7j6vZ6NXLXNbNwHqlr60lda4cW61qLmYVzaXG7J3XtfuSMjxBCCCGEEMLtSeMjhBBCCCGEcHtyqZsQKjpXCAdPQvb5G9O 2QRRodCqEUSGqJdNiNrKPQ HT0Nu/o1pC7ZAs1BoHQFhAeplE6I2TCY4fAaOnYXjFer6f9uU7fVtTcDfW718QgjrSOMjKnl8UlMKL ej1xvQ6w0kRHdhzKCPadSgmdrRLNJa7qzTsHIP7D8B5pt tiVL QMQHQa920D7aNDpnB7TbWitPkB7mc1m2H4MVu DnPO3/nzDQeVvHdCmKdzbFmIaOjWi29FajYD2MpcaYX0mrMtUPqi62ep9yt8eeujYAu5rCw0DnRrRrWitPspoNXd9JJe6iXJnLxzn/MU8Zo7exg/vFPLlq1kUFJ3mk6Xj1Y5mkZZyl5TCf7fChythXxVNz82yz8Hc9fDpWrhU7JSIbkdL9VFGa5kvFMMnafB/G6pueioyA3tPwIc/waKtyjohbKe1GgHtZT6RD9N/hO 3Vd30VFRqgvQsmLoU0vYrHwQI22itPspoNXd9JY2PKHcgJx0/nyCiG90GQJB/GLGRnSgoPK1yMsu0kvvqNZj9C6w7YPtrM3Jhxk9wvoadr7iVVuqjIi1lPlcIM36EPcdte50ZWHsA5qyBq9L82ExLNVJGS5kPnVK2uSfya563omtG O43WJAuzY ttFQfFWk1d30ljY8odyAnnfioO9DpdBiNpaTvX0Hazm/p3fEJtaNZpIXcJjN8vg6yztR GecKYdYvcOWa/XLVB1qoj5tpJXNxCcxaBflFtV/GoVPwxTo5SLSVVmqkIq1kziuAj9fU7WzkxoOwfJfdItULWqmPm2k1d30l9/iIcgdy0sk4up4BrzeguKQQH08/xg76lJT2g9WOZpEWcm88CPvzLM8zY6jy9 ivq5/nzCVYvA0G32W/bO5OC/VxM61k/n4bnK3hLKQ1db33BGw6BMmt7ZfN3WmlRirSQmajSRlgpqazkNbU9coMSIiEFuH2y fOtFAfVdFq7vpKzvhUYDKZeP/992ndujU Pj60b9 etLQ04uPjee6559SO53AHc7cy4Yl5fD pgG8m5hAeHEX26f1qx6qRq e fBV 2G6/5W08BDnn7Le82rpQrNzYu3gb/LzHdS/Dc/X6qIoWMh89C78ett/y/rdNOYOktpMX4KfdSl2vO6Csv65ICzVyMy1k3nSo5vvUrGUGFm5R/2ym2ayMsrhsp7Iv2npEuSTP1WihPqqi1dz1lTQ FQwfPpxJkyYxYsQIli9fzuDBgxkyZAhZWVl07txZ7XgOdfzsIS4V59M6qhMAYUFNeCxlLEs2z8ZkMgGQcWQDU YNK3/Nu98MZX/2FlXylrEmt9Fk5OWP7mHsrJ6MmN6BP01r59SMW47Y/x6GdZn2XZ4tjCblpvQ3v1MOVn/ZB0t2wFv/Uz4pdaUdqhbrWgs1DbC FveqWXK1VDkgU0txiXJ505QlsGyXUtf/3QqvL4KfMtQ/eK1I6toxzGZlBDd7ys1XPiRQy5lLMHUZ/HOlUser9sJXG GNRbAzW71cN9NiTYM26lpUJo3PdfPmzWPu3LksXryYcePG0atXLyZOnEiXLl0oLS2lU6dOakd0qAM56QT6hdI4pHn5tOTE/hRcOkXG0fUAxEV1JuvETgB2Z63D08Ob26LvVCVvGWtyG/QG3huxiglD5xER2oK3nlns1Ixbs y/zG3HlGFW1fDdb8pN6aYqDgS3ZCk7VVehxbrWQk1fM8IOBxw0pavU BhNMGd11QM0GE3KJ Wr9jo/V3Wkrh0j97xyxs/e0h2wD7DGpWKl4anqd7pSAnPXwV4bByVxFC3WNGijrkVl0vhcN3nyZPr27UtKSkql6bGxsXh6epKUlATAG2 8QVxcHHq9noULF6oR1SEyc9KJbdqx0rRg/3ASWiSzbpfye3p5 uDrHcDFonN8/uNrDO/3rhpRK7EmN0De SP8c1EqLz4ymyahMU7LV1IKJwrsv9xSI Q5YAddk7OXav5EdGe2Mgy3K9BiXbt6TYMy0lWpyf7LPZ6vTkOfkQtHavhUfsVu17gUD6SuHeWYg7Zbam0P1x6Ai8VVn60sm7R4u2uczdRiTYM26lpUJoMbALm5uWRkZDBmzJhbfpadnU1iYiLe3sojmfv27cvTTz/Ns88 6 yYDvXnh6dXOf0fqesqfd8m m5m/HcEXRIeJiSwsTOiWWRN7mOn9vLFj28wdvBnBPqFOCsaAKcuVn1mxB5O5EOzUMcsuzpbrPjkUgdsPqQ8fFVtWqxrV69pcEwzD8rZldMXoamTf6XNh5UHBFs6ACw1Kg9odYUBGKSuHSPPQXV9okDZD id/BDqTYcs/9yMcjYo57z622st1jRoo65FZdL4oDQ ABEREZWmFxcXk5aWRr9 /cqnJScnOzSLTueYLePK9 xz9J3QIpktK5YxYeg8q av6 9jj9yvfnwfoYFN PsXAwH421OLCPKrvmOw5/9B0/juDHp9baVpZaMBVae6n988etCIkc zc W/6pDOdvc 9zm3dX0SvaH6TYfRZOSbhT/y N0POjxPfa1rW2sa7FvXHe5/kZRhMypNs1dd335XV/IOOvd6ySenZBAWlWhxHrPZzLgJ77Bp4esOz6PFulZjWw32rev7RnxBm 5/qDTNHnVtMoOXtw/Ga84bKUOnN/DC/1l3c2nvvo9weOt3Dk5knxrR2rYa1K/r sJs5alLaXyA8HBlrMnMzEweeOCB8ulTp04lLy/P7Qc2sMXOw2tI7f8hHgZPtaNYbd5ruar926ZSx10bYyx1/nBTV4svKB NW2I2K/NpiNbqWs2aBjAZ3ayui/Ixm0zo9NVf/a3T6aSuHUztujY6oK7NZjM6nQ6T0blP6TWbjJSWFOPh5VvjvCUaqmut1TSoX9eiMml8gJYtW5KUlMTkyZMJDQ0lMjKShQsXsmzZMgCnNj7Wdqy2 vn9ur3 TEEu//r L0Q3TuD2 Pusfl1df5 65q4Ne/4fFF2FiTfdClbdcx seS5ERd/Nm0Ns4zm1D1cLWafhw5WW59EbPPh/44ew/F9DHJ5H6tp69qzrzJMwc1Xlafaq66P7t LrVftstbH2gDJSYU1 mv8eYQHvOTyPFutajZoG 9b1qj3ww47K0 pa1zqdjvBAnN74gDLKZnrWjft5qhLgDVk7V FhcHyeutRIfdpWg OOBYU0PgDo9XoWLFjAiBEjGDlyJGFhYTz11FOMGjWKCRMmlA9sUJ81bBDFm09/r3YMzfH3hrAAOOeAZ9xEOfn HoCYhsqfo2eq3pnqgIZBkBjp7GS1I3VdO1EOuky9YSBOb3oA7ohRhvotulr9fT4dmyvrshZIXddOMwfd5 LsezHLpNx2fYh4c/XNT 8EnNL01JXUtLAXGdXturi4OFavXk1RURHZ2dlMmjSJ3bt3k5CQgK9vzaeKhahOUjP7LzO CfiocKZfp4PhPW7cfH7zRW9hgTCyNxhky LW/LyhtQPuK3bEumINXy9I7a18 l1R2VWdcRHw 7udn0s4V8uGt9aAPbRXqa4jQ DZHrduj8u22z3ioVcbp8cSQlVyeGLB1q1bb7nM7fXXXycqKopNmzYxYsQIoqKiOHzYjo8vF26na tbGwR7LFMtAT7wUl94ujskVDizM7QLjH8QQvzVyyacp2ucfZenQ90R05qGwMSH4bE7bkxLioIRveDPvcFbro9wex4GuDvWvssM8oF2KjU AG2j4PX 0LfChSt3tVK24Y/cXvMtm0K4G9mUV6OwsJDMzExSU1MrTZ80aRKTJk1SKZVjfbrsr w5uoHEFl2JahjP/NXvMvrROSS2SOalmT04cnI3s8fsIDI8luKrhbwypw RYbG8 sRXLpe5fasU r8eXD6 ftkoKq999hCFxQXMGLXeaRnDA5UDug0H7bO8lg2VnZmaDHroEK38KbvG/Y6W6maqiqX6 Nf3z3MkbzdNwloy5rGPMegNqtRHTVmrWv8uFJ3ljc8fxmDwxN8nmNee/BaTyejUdTIpCmLCa37 jbW6xat/KZmPJ3SLg4XpyvfP9FA3z81sqZFDx3cwa/FoAE7nH2Ng9xd5pPtop9d4dZkbhUQzdf4f0KEjPDiK8UO xKA38Mbn/dmVlcYbwxbSKa4PgFMz92oDvx6GS1fss7wHO6h/BjzYD/q2gxW7lO9d6eylrfXxzNR4QgObAPDCIzNp3jiBaQv yK7Da/ji1RrG71Yx9/xfprB53xJCAhrzypD/w9fLX5XcQs74VCsgIACj0cjzzz vdhSnOHIyg6IrF5meupaLl89xpaSIQSkv075VCga9B28 /T3d2z1WPr vdwATh85XMbHlzAAxEe2YNnIN00auKR868u1nl6iS9eGO9jmo8zTAkC7Ofx6EFlmqjwM56ZSWljBt5BqaN07k171KXahVH7aufwG IfwjdT3TR6YRF9mZzXuXOH2d1OuVWvS0w/0B4QHwUIe6L8ed2VojsZEdyrd/MU2SuKvNQ4Bza9xS5gCfBrz9zBKmp64lIjSGLfuVwYRefHQ2j3QfXWk5zszs7w2D77TPshKawp0u IGQq6hNfQT7Nyyv6 aNEwAYO gTQgIjLP1TquY dzGP3UfWMWPUenp3fILlv36iSm6hkMZHAJBxZD23xykjpXRqfS96/Y2jGZ1O5xIPCruZpcwA2af3MWZmdz5Z9qrqI6R4e8JzvSDQp/p5Rn9teYQgDz0MT1FuABc1s1QfeeeyiGmiXPvRqmkH9hxz7nNjbmbr mfQG9BfH3rZaDYSGa7ONWKNguCZ7pY/0a6proN8lcvJ5FIyy2q7jS4uKSL/0kkiw 18DZcVLGUO9AvB3zcYAA DJ3qd8rOwoCZOz3mzds2gfyfL89RU181C4Q/d5FIyS2pTH5cun elmT2YsXAEJdfsdFrORrbmPl2QXd6ktWragb0q72/qO2l8BKBsTL746W MndWTb1a9w6XL59WOVKOaMs8df5DpI9kNwZ0AACAASURBVNdSeDmfTXt/UCnlDY2D4IX7ajciVgM/ZdCA29Q/JtAMS/UR1TCeXVlpAOw49AtFxQ56ZLuVarP 7c/eQuoHt7Pj0C80CY1xQsqqJUQq98AE12IMmGah8MK9ykiAwrLabqPT9y/n9vi Dk5XNWsyn71wgt8yV5YfSLqKXm2US8K8atGQt4uCUX3UGYBGS2pTH/8YtZ7pqWtpFNKcpb8693EOZWzN3SS0JftztmA0lrLj8C8Uqry/qe/kMzYBQKBfKE/d/xbJiQ zee8Szlxw/Qdu1ZS57PK25LYDOHR8O8mJD6sRs5KGgTCmL6zaC6v3QXENz8vzMMBdLZXLgNQY5lfLLNVHbGQHWkS0ZdzsXrSIaEsDlc9o1mb9uy36Tma uJUFadNYkf4Zj/YY44SkVWvdGMY/BEu2w5YsKDVZnt/PSzmw7J2g/v0PWlHbbfSGjO8Y3PMVB6erWk2ZS0qv8t63T/HSoI8xGFzvcOTuVhDbCL77DfYet/w8HFAuZ34gCTq1kDM91qhNfZTt17u2Hciidf9wemawPXeDgIb06fQkr8zpw23N7iQkwPWuoKlPZJcjAGgb043dWWsB5cnIJpNR5UQ1s5S5uKQI4/Xv9xzdQNOwVqpkrIpBD/e1hTcHwhNd4PYWyiVD3h7gZYBQf2VY3wGdlXkG3SlNT23UVNPD7n2D9/ 8miC/MO5q86AaEcvZuv5dK73RMfv7BOHlqf6Q 35eMPgu PtAGNBJqeEQf6WmvT2UM563t1Bq/u8D4d620vTYojbb6FLjNbJP76NV0/aOjlelmjLPWPgcDyePKr8MyBWFB8Kfeioj/t3fTjnrHuSj3Nvm46mctewSe32e30HnGGl6rGVrfVwrLaGk9Cqg7NebqLRfr01dP3DXn5g2cg3RjRNU39/Ud673EYtQRUxEWzwMnoyd1ZOE5l3w8fLHaLrxpOlJXw4m4 h6jp89yOM9XyG5bX8V0yosZT5 5iDTFjyLr1cAEaEx/OG N1VOeysvD XGV7n51TEs1YfJZOLlj3qj1xvoGHsPbaLvctmscOv6FxrUhDlLX0av0xPoG8r4IV qmL6yAB/o2Ub5I ynNtvo7Yd oUOr3i6Zee/RTazPWMSp/GMsWjeDgd1epFu7gfz7 xfYvG8Jm/Yu5qHzf bBu59TLX9F4YHQT55lble21kdii2QmfNoPX68AAnxDeHWIOiPK1qau3/pyEJcun6dlkyRGPDRNldxCIY2PKDf8gXfLv167ayHzV08hMrw17Vul8Pqw/1Sat/hqIVPmPUl8sztuXoxTWco8a/S2W Z/7bOHCHWBm2eFc1iqj2kj19wyv5r1Ycv6BzB9ZFql711lnRSOY2uN3BF/P3fE319pmrNr3FLmxW9fumX UQM ZNSADytNk 22 7K1Pqrar09b8Ed0Tj7NZmvuN4YtuGWaGrkF6MxqD3clnOLn99X5d/uMq9vr1chd18z1TdnIRjOGOv/flrq2ntS19dSsadBmXWsxc31TH tai9tqkLp2JGl8hBB1ovbOVAh7k5oW7kjqWggZ3EAIIYQQQghRD0jjI4QQQgghhHB70vgIIYQQQggh3J40PkIIIYQQQgi3J42PEEIIIYQQwu1J4yOEEEIIIYRwe9L4CCGEEEIIIdyeND5CCCGEEEIItyeNjxBCCCGEEMLtSeMjhBBCCCGEcHvS AghhBBCCCHcnjQ QgghhBBCCLcnjY8QQgghhBDC7UnjI4QQQgghhHB70vgIIYQQQggh3J6H2gGEc/z8vjr/bp9xdXu9GrnrmlmImmixrrW6DRHOodX60OK6KJxH6tp6WqlrOeMjhBBCCCGEcHvS AghhBBCCCHcnjQ QohauXINjp698f3lq plEcJeLl258XX2OSgpVS LEPZgNsO5whvfn8gHo0m9PEKoSe7xEUJY7UIxbD4Evx2F0xcr/2zCQggLgI7NITkWQgNUiSiEzU5egA0HYXcOFFy MX36CtDpoEkDuCMG7moJft7q5RTCWmYzZJ2BDZmwPw8ul9z42dRl4KGH6DDoEgsdmoOnQb2sQjiTND6ikscnNaXwcj56vQG93kBCdBfGDPqYRg2aqR3NIq3m1gqjCX7eAz9lWP6k8FyhMt qPdCzDfRLAi/ZytSKVmtaS7kvl8D3v8GWrOrnMZuVT8j/lw8rdsHvOkJya9DrnJfTnWipPirSUu5zhTB/Mxw8Vf08pSalMco6A0t3wuN3QZumzsvobrRUHxVpNXddyKVuotzZC8c5fzGPmaO38cM7hXz5ahYFRaf5ZOl4taNZpNXcWlF4BT74CZbvsv7yCDOwep/yiXl kUPjuSWt1rSWcp/Ih6lLLTc9N7taCgvT4ZM1ytfCNlqqj4q0lDsjF/7fUstNz80KLsNHq5UPAUxmx2VzV1qqj4q0mruupPER5Q7kpOPnE0R0o9sACPIPIzayEwWFp1VOZplWc2tB0VX498/KvQ61cfIC/HNl5cuHRM20WtNayX3yAvzr59rX5d4TMGc1XDPaN5e700p93EwruTNy4bO1tb8vbc1 WJSunOUU1tNKfdxMq7nrShofUe5ATjrxUXeg0 kwGktJ37 CtJ3f0rvjE2pHs0irubVgwRbIu2B5nhlDlT/VOV8EX2 UTxJtodWa1kLua0aYu67yPQ83q6mmAQ6fhqU77JvN3WmhPqqihdz5RfDlBsvbWWvqev1B5R5OYT0t1EdVtJq7ruTq wpMJhPTp0/no48 Iicnh/j4eD788EOee 45UlJSmDNnjtoRHepATjoZR9cz4PUGFJcU4uPpx9hBn5LSfrDa0SzSam5XtzMbdmTbZ1kHT8Gmg9A1zj7Lqwuz cZBr5 XcvO6q9FqTWsh90 7lTM 9pC2X7kxvEW4fZZXF0aTUtceevD1UjtN1bRQH1XRQu7/bLHf5Zf/3QrxTSDQxz7Lq4tSIxSXgJcneLvoEasW6qMqWs1dVy5aRuoYPnw4ixYt4vXXX6dz585s3LiRIUOGcObMGV566SW14zncwdytTHhiHt3aDeTcxTxe egesk/vVztWjbSa25WZzbAyw77L/HmvMoKQXqXzzEYTrM EdQfg7PWhXcMDoXscdIsDgwud/9ZqTbt67ivXYO0B y3PjDKYxx9T7LdMWxVdhTX7YOMh5WtQGrGebaBDtHq5quLq9VEdV8 dcx72nbDf8opLlFEO 7az3zJtdfYSrNoLW48oZ2l1QNso6J0AMQ3Vy1UVV6 P6mg1d1250K5eXfPmzWPu3LksXryYcePG0atXLyZOnEiXLl0oLS2lU6dOakd0qONnD3GpOJ/WUcrvGRbUhMdSxrJk82xMJuWO9owjG5gyb1j5a979Zij7s7eokreMNbmNJiMvf3QPY2f1ZMT0Dvxpmopbc43IPge5 fZdZn6RfXfOtig1widp8N1vlZ9nce6SMu3Tta7zXAtZFx3ntyP2H5RgT656A3hcugIzfoSVe240PQDHzimX87nSpXhS146z8aD9l7npIJhU2ibmnof3lyuPTii7j84MZByHD1fCtqPq5KqK1LX2SONz3eTJk nbty8pKZU/uouNjcXT05OkpCTy8/N56KGHiIuLo3379tx3330cOnRIpcT2dSAnnUC/UBqHNC flpzYn4JLp8g4uh6AuKjOZJ3YCcDurHV4enhzW/SdquQtY01ug97AeyNWMWHoPCJCW/DWM4vViqsZmScdtFwbRhqyp1/23Wi6Kl4CX/b13uPKKHSuQNZFx3FEXZuBQyrV9YItyifjNyu7OX3lHjiQ59xM1ZG6dhxH1PWFYjhdRW05msmkfBB1tbTythqUujablXtGL7jIgDlS19ojjQ Qm5tLRkYGgwYNuuVn2dnZJCYm4u3tjU6nY/To0WRmZrJz504eeughnnnmGRUS219mTjqxTTtWmhbsH05Ci2TW7VoIgJenD77eAVwsOsfnP77G8H7vqhG1EmtyA SdP8I/F6Xy4iOzaRIa4 yYmpNz3kHLreXocHVhNCmXt9Vk3QH1PuGsSNZFx3FYXTtouZbkFykPXLU0ZohOZ99L  pC6toxiq5WPottT7UdzbMu9uUptW1pZDmjGTYddl4mS6SutUdnNsvAhZs3b6ZLly4sXbqUBx54oHx6cXExrVq1ol /fnz66ae3vG7r1q0MGDCA3Nxcu2XROehO65Xv2ee/efbilzhdkE2b5l0YlDK2xvnvfbluv489ch87tZcvfnyDMY99TKBfSI3z1zWzOxj89500ib270rSaRgOqyuivK39fcOowX4yNrUMy2zVoEsdT71l39Df3pVZcOG3Dg11qyR51XR/WRbDv viXuVcxeNy48782NQ231vXBXxew7J/OvSE4vssQ o76psb5Sq5cYtYfgxyeR/Yx6uxjwqLa8uSU3ZWm2auuN3z7V7b MKWWyWqn25D36NTvJXQWbgY1m03k7ktj0eTeDs8jda2dYydr2xk54wOEhytD8mRmZlaaPnXqVPLy8ujcuXOVr5sxYwYDBgxweD5XktAimaMnMxjY7QW1o1jt1Y/v41T Mf7 xUDGzurJxcsqfDyrMTocswFzVGNviV5nsHpend76edUm62JtOKj VKhra2tVZ0P9uwKpa9s4cpuqyvZab8DyeUwA3fX5tEPq2nXIqG5Ay5YtSUpKYvLkyYSGhhIZGcnChQtZtmwZQJWNz5tvvsmhQ4f45Zdf7JrFUSfgfn7fPsvZeXgNqf0/xMPgadX8df197JF73mu2nZGTk6DKwxn33jQQwc2fBpYp 3Sxup9X1K5NS6e/vyWl8MYiZUQvS3w94UxuJh5O2J/ao67rw7oI9l0f31gEF4tvfG PmgZ4auhjLP3AuXV98gJMWWJ5Hh0QH XnlHVO9jHWs f/x8Vipa4rslddfzBtMl3 N7n24Wph6xH4aqPleXQ6HU8M6MGCSVLX1dF6XTuSnPEB9Ho9CxYsIDExkZEjR/LMM88QHh7OqFGjMBgMJCUlVZr/7bffZsmSJaxYsQI/Pz VUjvXmYJc/jZ3AH4 Qdwef5/acYSDRYVqa7mWeHnAXa1qnq9LLE5peupK1sXai7Lu6jqbNVOhriOCoWVDyyebzChDtWuB1HXtBPkqfxxBje11  jrz1erYb7k1k6JU2dS165HzvhcFxcXx rVqytNGzZsGAkJCfj63tiqvPnmmyxbtoyVK1fSoEEDZ8dUTcMGUbz59PdqxxBO0qqRg5ar0vMX7m njG5V3YMrmzSA zQyUqesi7XXqtGtZzLtoaVKdf34XfDBT8pzV6r6rLVDtPKAVS2Quq69Vo1g zH7LtPXE5oE23eZ1vA0wNBk DQNMFeua50yiYc6QGPH37ZmF1LXrkfO FiwdevWSpe57dmzh7///e cO3eOnj170qFDBzp06KBiQiEco3UEhAXYd5n 3pCk0gMV/bzghfuUTwk9K5zV8TRAt9bwwr3gY90VCELD7mxp/wfVtmoEjVU4QATl3x1zP7RrVvnMT4A3PNAehnUFvYzV4va6OGC8mDtbqXcGPDESRt1z64NKwwNhaBfok6hOLuEe5IxPNQoLC8nMzCQ1NbV8WmJiomauYRSiLvQ66NUGFqbbb5nd4yo3Hc7m5wWD74SHO8Kr/1GmTXpUGp76JNAX7oiBzXYcCrdXG/stqzYaBsGzPZTnrvzt r0ebz5i/wZPuK7WjZXL0nLtdO 5h17ZXqupVWPlw6qzl Dt64 PmfA7VcYREW5GGp9qBAQEYDQa1Y7hNJ8u yt7jm4gsUVXohrGM3/1u4x dA6NQqKZOv8P6NARHhzF CFfYtAbeO2zhygsLmDGqPUumbt9qxR y1zJ/F/exWQ2MeJ302jWMJ5X5vQhMiyWV5/4StXcWpDcWnlCdtaZui rSbDrfEpXsdFxpaanulpObJHMSzN7cOTkbmaP2UFkuPLx7n/WvMfGPf jcUhzXn58LtdKr6pS37bmBli3exGz/vci37yWQ/HVQqfmfrijcrlbxUEOaqtjc2gbVffl2ENwhfs8XKnpsaU TCYTU f/gdMF2XgYvJj45Hy8PHycXte27A vllzmb18MwGi8hp9PEBOemIefTyDTFvyRXYfX8MWrjn/IuU4HQ 6GacvBZIfPZvslKWdXXEHFHK7U9Nh6zHTo HY XvoKRlMpg1LGcVebB51 HFVd5hYRibzx cMYDJ74 wTz2pPf4u3p6zL7GHtzoc2jUMuRkxkUXbnI9NS1XLx8jislRQxKeZn2rVII8GnA288sYXrqWiJCY9iyXxnp7u1naxhOyAks5b56rZilmz9iynMrmTZyDXFRnfH1DmDi0Plqx9YMvU651jrQx/J8o7 2PEqQj6dyyY0WBg5Qi6VaNug9ePPp7 ne7rHy fMLT7Pj8GpmjFpPTJMkNmR8r0p925q7zLpdC2nYoBmA03P7ecOwZMvNQU01DdAoCB67w77Z3I2t9XH4xA48PLyYnrqW  94hlXbvnZ6fdi6P/QwePLqkK YnrqW5MT /LR1LgBjB31CSGCE03JHhsDAqp 8Uc6auk5oqv5ZTFdXm2Omr36exJtP/4/3/7yau9o8CDj3OMpiZt8Q/pG6nukj04iL7MzmvUtcZh/jCNL4CDKOrOf2OGW0kU6t7600Pn6gXwj vsoF7B4GT5ueieJolnLvPbYJnU7PhE/6MWXeMIpLitSKqWlhATCqDzSo5eCFfl6Qeg80ddBoWu7CUi3rdDpCAhtXmj8zZyvtW/a8Pn8f9h3b5KyoldiaG DXfcvo1LoPOp16u5/WEfDHlNpfehkRrNS1v7d9c7kbW sjPDgSk0m50qKwuIAgvzDnhb3O1v2hl6cPYUFNADDoPVV9vkz3eBjQqfavT4yEp7uDhWeHCmyvkbxzWZSUXuGtLx/jb3MHkH/plEtlNugN6K//pxvNRiLDW7vMPsYRpLwFly6f54uf/sbYWT35ZtU7XKriIVVnL5zgt8yV5SuOK7CUO//SKc5fymPyH5eT2DyZpZs UjGptkUEw8sPQOcWtr0uMRLGPwTRzj920Rxr1sGKiq4U4OejDGvk7xNM4ZUCZ8S8ha25AVb 9gX3dHrSCeksa9MUXnng1huoLdEBPeLhpb61/zCgPrG1PoL8w7l6rZhn32vDkk2z6NbuESclvaG2 8Piq4Us3fwRvTs 4cy4t jZBp7vY9vgNJ4G5WzR8BRl H9hma01kn/pFMfPZPLGsIU8ePcIvln1jstl3p 9hdQPbmfHoV9oEhrjMvsYR5ASFwT6hfLU/W RnPgwm/cu4cyFyg tKim9ynvfPsVLgz7GYHCdkrGU298nmLYtumHQG gQ25sFaXZ6Clk95e tXK7WLQ7WZ8KObDCabp1Pr1PueegWp9xw60rXZLuymtbBm/n7BHOmQJnn8pWLBPioM7S rbm3H/qFhOZd8PTwclJCyxoGKQeJe47D oOQmVf1sNBeHkrj3z1Ozl7awtb6 C3zJ4L9G/LZy/tYu2shC9LeZ9i9bzgpraI2 0Oz2cz7/3mWZ/q9Q4Cv o 5aNUYxj8Ivx2FDZmQm1/1fIE cHcr5X7OEH nRtQ0W2vE3zeYuGZ34OPlR4fY3vx37XSXy3xb9J3MfHErC9KmsSL9M6IaxrnEPsYRXOcoVqimbUw3fkqfS3Liw w8vIaI0BgM hulMWPhczycPIrmjRNUTHkrS7njm93Bsl/nAMp14xGhMWpGdRsxDZU/Q 6GEwXKc3FKSpX7dyKCoWkD cSwNmpaB28W1 wOFm cyeO9XmHbwZ9p0/xuJ6a9wdbcR09msGnPYtL3r DYqT18vuI1nun7thMT30qvV4aDbtcMrlyD4 fh9CWlsffxVGq6cbBrDRigFbbWh9lsJshPeWpmsH84RVeqefCWA9Vmf/jFj2 Q2KIrHWN7Oz1vdbw8lGGuu8RC4RXIOQ/5RWA2K/e5RYUqZ4VkuHPb2VojkeGtKSg8jdFkVO14xFLma6Ul5R9G fsEYTQZXWYf4whyiCKIiWiLh8GTsbN6ktC8Cz5e/hhNpQDsPbqJ9RmLOJV/jEXrZjCw24t0azdQ5cQKS7kbBDQkqWUKL83sgbenH38d o3Kad2Lh0G5hE0uY7MPS7UMMOnLwWQcXc/xswd5vOcrJLftT7uWPRj97240ahDNI91HayL3wG4vMLDbCwCM/nc31Zuem/l4Kp Wt7r11iRRC7bWx11tHuTH9M8YO6snZrOJcYM/d6nMVe0Pb4u i2/X/D8SmiezIeM7erZ/nN8lj3R6bksCfJRLO4V91OaY6YG7/sS42T3R6fS8/Phcl8p8 MQO5ix9Gb1OT6BvKOOHfImPl59L7GMcQRofAcDwB94t/3rtroXMXz2FyPDWtG VwuK3L90y/2ufPUTo9Rs61WQp96M9xvBojzHlPy  WsiUeU8S30yGYhKux1Itvz7sP7fM//te4/l9r/Hl36tV37bmLlM2hKusl 7N1vp44w8LK32vRn3Yuj9cPqXklmnTFvwRnVzr67ZsrZFeHX5Prw6/rzTN2cdRljJPH5l2y/yuso xN51ZnshZL/ys0i0ufcbV7fVq5K5rZuH6yoZ0nTFUnX9fi3Wt1W1IfaJmXWu1PrS4LtY3Ute2k7qunjQ Qoh6R 3GRwhHkLoW7kjqWtiT3K4phBBCCCGEcHvS AghhBBCCCHcnjQ QgghhBBCCLcnjY8QQgghhBDC7UnjI4QQQgghhHB70vgIIYQQQggh3J40PkIIIYQQQgi3J42PEEIIIYQQwu1J4yOEEEIIIYRwe9L4CCGEEEIIIdyeND5CCCGEEEIItyeNjxBCCCGEEMLtSeMjhBBCCCGEcHvS AghhBBCCCHcnjQ QgghhBBCCLfnoXYA4Rw/v6/Ov9tnXN1er0buumYWoiZarGutbkOEc2i1PrS4Lgrnkbq2nlbqWs74CCGEEEIIIdyeND5CCCGEEEIItyeXugkh3J7ZDIdOwaHTkHPuxvQ5q6FZGMQ2hthGoNOpl1EIW5UaYc9xOHoWjuffmD53HUSFQkJTaBqiXj4haqO4BHblQPY5yCu4Mf3rTRAdCm2jIMRfvXxC26TxEZU8PqkphZfz0esN6PUGEqK7MGbQxzRq0EztaBZpNbdwLLMZfs2CVXvgzKVbf773hPLnx93QKAjuSYA7W7pGA6TVmtZqbi25ZoRVe2F9JhReufXnO7KVP0t2QMuG0DcJ4iKcn7MqWq0PrebWksIrsGwXbM2CEuOtP0/PUv4s g3aRsKDHSAi2Pk5q6LV tBq7rqQS91EubMXjnP Yh4zR2/jh3cK fLVLAqKTvPJ0vFqR7NIq7mFYxVchlm/wPzNVTc9Nzt9EeZthtmr4cJlx ezRKs1rdXcWpJzHt5fBit2Vd303CzrDMxcBf/5FUpKHZ/PEq3Wh1Zza8muHJiyBDYerLrpqchsht25ynrwy17lezVptT60mruupPER5Q7kpOPnE0R0o9sACPIPIzayEwWFp1VOZplWcwvHOVcIH/wEmSdtf 2BPPhwJZwvtH8uqzNotKa1mlsrDp2Cf66EUxdtf 3GQ/DRariqYvOj1frQam6t2HgQPlsLhVdte12pCRZvhwXp6jY/Wq0PreauK2l8RLkDOenER92BTqfDaCwlff8K0nZ S  OT6gdzSKt5haOcfUazP4F8otqv4xzhcqZH7U IddqTWs1txacvggfr6lbTR4 DV9uUO8gUav1odXcWrAnFxZsqdsyNh6EnzLsk6c2tFofWs1dV3KPjyh3ICedjKPrGfB6A4pLCvHx9GPsoE9JaT9Y7WgWaTW3cIwfdtR8aduMocrfo7 ufp7TF5V7JB653X7ZrKXVmtZqbldnMsG8TTWfrbGmrjNyIf2Ici bs2m1PrSa29UVXYX5v0JNfbg1df3jbkiMVAb1cDat1odWc9eVnPGpwGQy8f7779O6dWt8fHxo3749aWlpxMfH89xzz6kdz EO5m5lwhPz H5SAd9MzCE8OIrs0/vVjlUjreYW9nfygnLDt72sO1C7y4rqSqs1rdXcrm7bMThy1n7L 982ZYAEZ9NqfWg1t6tbmQGXrLhPzRomM3z3m32WZSut1odWc9eVND4VDB8 nEmTJjFixAiWL1/O4MGDGTJkCFlZWXTu3FnteA51/OwhLhXn0zqqEwBhQU14LGUsSzbPxmQyAZBxZANT5g0rf8273wxlf3Ydz1HXkTW5jSYjL390D2Nn9WTE9A78aVo7NSMLB9pgx6YHlE8iNx607zJrIuuiuJk9m3lQPmnfccy y6yJ1LWoqKQUfj1s32UePl15 GtnkLrWHml8rps3bx5z585l8eLFjBs3jl69ejFx4kS6dOlCaWkpnTp1UjuiQx3ISSfQL5TGIc3LpyUn9qfg0ikyjq4HIC6qM1kndgKwO2sdnh7e3BZ9pyp5y1iT26A38N6IVUwYOo I0Ba89cxiteIKBzKb4bej9l/ub0ece0 ErIuionOFynN67M0R64olUteion0noPia/ZcrdW2d lzX0vhcN3nyZPr27UtKSkql6bGxsXh6epKUlATAgAEDSEpKomPHjtx55538/PPPasS1u8ycdGKbdqw0Ldg/nIQWyazbtRAAL08ffL0DuFh0js9/fI3h/d5VI2ol1uQGyDt/hH8uSuXFR2bTJDTG2TGFE5wrhMsl9l9u4VVlaGxnkXVRVFTxgbt2Xe555zb0UteiomxH1bWDllsdqWvtkcENgNzcXDIyMhgzZswtP8vOziYxMRFvb28A5s6dS4MGDQDYvn07PXv25Pz58xgMBqdmtrc/Pzy9yun/SF1X6fs20Xcz478j6JLwMCGBjZ0RzSJrch87tZcvfnyDsYM/I9BPHmPurk5ecNyy8wqc96RwWRdFRXkOquuiq8r9FUG jln zaSuRUWO2l47an2pmtLN9AAAC8dJREFUjtS19kjjg9L4AEREVH60dXFxMWlpafTr1698WlnTA3DhwgV0Oh1mO35spnPQI NXvmefjAktktmyYhkThs6zav66/j72yP3qx/cRGtiEv38xEIC/PbWIIL/qh35x1P BcKzWdw3mgee/rTStbDSg6lT385tHD o/8DEOpf 3Dukqs0dd14d1EWR97Pr4FG7/XeUHCtqrrps1j HimaO1D1eB7GNkH2OLgeN/IrrdvZWm2aOuT5/NR6ez39BuUtfaqWtrj8Wl8QHCw8MByMzM5IEHHiifPnXqVPLy8m4Z2GDUqFEsX76cCxcu8N///hcPj/rzNu48vIbU/h/iYfBUO4rV5r2Wq3YE4QTGa3YaHqgKpdeKHbbs2pJ1sX5wbF07btm1JXVdPziq9hy5vtSF1LXrqD9H7Ba0bNmSpKQkJk eTGhoKJGRkSxcuJBly5YB3NL4/Pvf/wYgLS2NMWPGsHbtWgICAuySxZ5njyr6 f26vf5MQS7/ v4vRDdO4Pb4 6x XV1/n7rmrg1H/R8IxzpzEd75ofK06p77YM1zISraum4pYfZZxYG61XV9WhdB1sdtR H/NlSeZo 69vOCwvN52OtDWtnHWK  1zTA4m3wy77K0 xR1x3im9j1/ZW6tp5W6loaH0Cv17NgwQJGjBjByJEjCQsL46mnnmLUqFFMmDChfGCDm6WkpKDX69mwYQP333 /k1M7V8MGUbz59PdqxxCiWmGB4OMJV w8UpCfF4Q66f4ea8i6WL80C3PMcqNCsVvTYw9S1/WLo raUcutLalr1yONz3VxcXGsXr260rRhw4aRkJCAr69y92dhYSHnzp2jeXNl L/t27dz PBh2rRp4/S8QojK9Dro2Bw2HbLvcjs2d60DRFG/hAcoTUruefsut2PzmucRwlHaNAVvD7haat/lSl2LmkjjY8HWrVu5  67y78vKiri8ccfp7CwEA8PD3x8fPjqq6 Ijo5WMaUQokzX1vZvfLq2tu/yhLCFTgfdWsP8X 23TB9P6NTCfssTwlY nnBHDKy34wOim4dBM/uNayDclDQ 1SgsLCQzM5PU1NTyaY0bN2bz5s0qpnKsT5f9lT1HN5DYoitRDeOZv/pdRj86hxYRibzx cMYDJ74 wTz2pPfYjIZeWVOHyLDYnn1ia9cLrO/TzCzFo8G4HT MQZ2f5FHuo/mtc8eorC4gBmj1quWWThOVCjc2RK2ZNlneV1ioamTRvGsrpYTWyTz0sweHDm5m9ljdhAZHkup8dot04qvFqqyTtqSG6D/68Hlz48oGyVI1kvL7mipHCDa66zP7zoon7Y7g6318a/vn dI3m6ahLVkzGMfY9AbnF4f1WVuFBLN1Pl/QIeO8OAoxg/5EoPewEsze4BOh0HvwYSh8wgJaCQ1bYX728H2Y1Bkh ev6YABnWuczW5sqZELhWd45 vfA5BfeIrb4 4ntf8Mpi34I7sOr GLV 38aZ2Nmas6xjPoPVxmH2Nv8gDTagQEBGA0Gnn  efVjuIUR05mUHTlItNT13Lx8jmulBQxKOVl2rdKIcA3hH krmf6yDTiIjuzee8SfL0DmDh0vstmjo3swLSRa5g2cg0xTZK4q81DALz97BJVMwvHG9AZQvwszzP665pvlA31h/6d7JfLEku1bNB78ObT39O93WPl81c1TY110tbcADER7crXzbKhUWW9tMygh6FdwKOGx8VZU9fxEZDspLOYttbHgZx0SktLmDZyDc0bJ/LrXqUunFkfFveFPg14 5klTE9dS0RoDFv2KwMgTR2xiukj07i38x9YufULp2fWqkBfGHRnzfNZU9e9EyCmoX1y1cTWGgkNiijf5nWOu4 7rx PjB30CSGBETX8a07IXMUxnqvsYxxBGh8BQMaR9dwep4w40qn1vej1N/awBr0BvV4pFaPZSGS4a1z7YylzmeKSIvIvnSz/NFG4Pz8v HNvCPSp/TKCfJVl Dhp5FFLtazT6W554F1V09Rga26A7NP7GDOzO58se1UzowC5giYN4Nnu4FGHvXZ0GDzd3Xn3rNlaH3nnsohpogwm1KppB/Yc2 icoBVYyhzoF4K/bzAAHgZP9DpD dcAV68V07xxopMTa1uH5vBwx7oto3MLeLC9XeJYpTY1UmZ31lrat rptKxlbD3Gc5V9jCNI4yMAuHT5PF/89DfGzurJN6ve4dLlytdU7M/eQuoHt7Pj0C80CY1RKWVlNWUGSN /nNvj 6qQTqipcTC8eB80D7f9tTHhymsbBdk/V3WsqWVXVJvcc8cfZPrItRRezmfT3h9qnF/ckBAJI  BkFqMMtixOaTeA75e9s9VHVvrI6phPLuy0gDYcegXiooLnBGzEmsyn71wgt8yV5YfSJ7Oz aFf3Zh8YZ/EdOknbMja17vBHiii 0fNOl1cG icjZU78Sj2drUCMCBnK20bJKEweD8u0y0eIznKHKPjwAg0C Up 5/i TEh9m8dwlnLlR cNVt0Xcy88WtLEibxor0z3i0xxiVkt5QU2aADRnfMbjnKyqkE2oLD4QX74W1B2D1PrhQwzNIG/hBrzbQPc65O1GwrpZdUW1yl13eltx2AIeObyc58WFHx3QrrRrB Afhx93KQB41Dd/epAH0S4KkZs7JV5Gt9REb2YEWEW0ZN7sXLSLa0kCFT5xrylxSepX3vn2KlwZ9XH4A2ygkmg f38S6Xf9lQdr7jOr/gdNza92dLaF1Y1iyA3Zk///27i kyjuO4/jHxx3/ZB2mVpZlYmfupDhaiFhNUlx/wP6ouyiC9KJJFEEblI3BckRBXuTwTuoiBl5MUroY7GZEbc7R1dZorlitsI0FSdBFhaw61sWvHHbinD3WeZ7n/PZ gfDw8yfne/F9/P0 5znneaTYZOL55UXS5hXmKqbXZtIjktmP1FV94HW5ktJzj5cqXPGBJKmqrE6/3hiWZJ4wPDkZm/rdo8f/fvMwLyesrFCu5/W9TKKaJelx7JH HL iSLGH18ARKI4jNVRIXS3SzjXSmqi0dJ4JRfPmmOP6ZVJHvXSo2Rx7HXqk5L0cVG7rnnj4QLFnc34b 1HFhZGU12ijnJD5/tnhVmn7SnMTjpICqXC26etlC8074fvWSQeb/Ak90sz6um1dl47vPq/wrELVVmxMdYlxktXcO7RLW1bvVWlRpSSzzjz/yOasnLCyA7I pqP8PKntPenzFqm12nyErfhN09dFYdPHTculTzdJe9f6E3ok9z3y3E9Xv1W1i4eYvk7puMdLFa74QJJUtqBKb2SGtL vQZWlq5STlafYpLnB/vVbv jkN51yMhzNyS3QJ9v7fa7WSFSzJF3845zejTT6WCGCItMxi6ZfG8BkkvXykf6tGh0b0d93rmlbw0Gtrmp 6VjQ656fX6qewZ3KzZqtBQVlal9/2POabZIdkmoj5ieI3PbHysrN6jzRKMfJ1Iq33lfFktpA1Xx57IJGRs/o9t2bOvNDr1rrPtLbi6vVPdAmJ8NRKDNbndu 9Lxm24RzzZtQQeW2R readVf47 rKL/Ut2A8kz1eENaYVCD4YMqHTcemjocvDWngfLcWzS3X8ki9vtjz/bS5E//cV/dXOxQtqfG6zGkS1VwT3aCa6IZp8z87tUkF4YVelwkklaiXD7Wdjpv/4phf56Tbuvs /jlujPPSXm77o2fPd3FjXvdHopq/Pnovbv6L66NET9vObY UzI qq31o2ljPYIcyPHw6tps9nhScNeZ1y3jCbXX F84e9 d11x54tb/3o 5XrRlIJh37Ol3/h8Ab6dof6Xguwjv09X XLn1N8AEAAABgPW5uAAAAAMB6BB8AAAAA1iP4AAAAALAewQcAAACA9Qg AAAAAKxH8AEAAABgPYIPAAAAAOsRfAAAAABYj ADAAAAwHoEHwAAAADWI/gAAAAAsB7BBwAAAID1CD4AAAAArEfwAQAAAGA9gg8AAAAA6xF8AAAAAFiP4AMAAADAegQfAAAAANYj AAAAACw3lOPJI2lLlzQrQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1081.19x264.88 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "VQE Result : -1.0731275404157\n"
     ]
    }
   ],
   "source": [
    "from qiskit.aqua.algorithms import VQE, NumPyEigensolver\n",
    "from qiskit.aqua.components.optimizers import COBYLA\n",
    "from qiskit.circuit.library import EfficientSU2\n",
    "\n",
    "from qiskit import IBMQ, BasicAer, Aer\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "#Define Backend\n",
    "backend = BasicAer.get_backend(\"statevector_simulator\")\n",
    "\n",
    "#Define Optimizer\n",
    "optimizer = COBYLA(maxiter=500, tol=0.0001)\n",
    "\n",
    "# Distance \n",
    "Distance = []\n",
    "Energy = []\n",
    "\n",
    "Distance.append(0.74)\n",
    "\n",
    "\n",
    "#Define Number of qubits\n",
    "n = 4\n",
    "\n",
    "#Define Ansatz\n",
    "ansatz = EfficientSU2(num_qubits=n, entanglement='linear')\n",
    "display(ansatz.draw('mpl'))\n",
    "\n",
    "#Operators\n",
    "\n",
    "operators = [op]\n",
    "\n",
    "# Run VQE \n",
    "\n",
    "for ops in operators:\n",
    "    experiment = VQE(ops,ansatz,optimizer)\n",
    "    energy = np.real(experiment.run(backend)['eigenvalue']) # energy shift?\n",
    "    \n",
    "    print('VQE Result : {}'.format(energy))\n",
    "    Energy.append(energy)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['']"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "''"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    ",\n",
    ",\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['-0.04919764587136759 ', 'X0 X1 Y2 Y3 ']\n",
      "['0.04919764587136759 ', 'X0 Y1 Y2 X3 ']\n",
      "['0.04919764587136759 ', 'Y0 X1 X2 Y3 ']\n",
      "['-0.04919764587136759 ', 'Y0 Y1 X2 X3 ']\n",
      "['0.13716572937099494 ', 'Z0 ']\n",
      "['0.15660062488237958 ', 'Z0 Z1 ']\n",
      "['0.10622904490856085 ', 'Z0 Z2 ']\n",
      "['0.15542669077992843 ', 'Z0 Z3 ']\n",
      "['0.13716572937099494 ', 'Z1 ']\n",
      "['0.15542669077992843 ', 'Z1 Z2 ']\n",
      "['0.10622904490856085 ', 'Z1 Z3 ']\n",
      "['-0.1303629205710914 ', 'Z2 ']\n",
      "['0.1632676867356435 ', 'Z2 Z3 ']\n",
      "['-0.13036292057109133 ', 'Z3 ']\n"
     ]
    }
   ],
   "source": [
    "op = WeightedPauliOperator([(1,Pauli(label = ''))])\n",
    "\n",
    "operator = []\n",
    "\n",
    "\n",
    "#H2 = ['-0.04919764587136759 [X0 X1 Y2 Y3]  ', '0.04919764587136759 [X0 Y1 Y2 X3]  ']\n",
    "\n",
    "\n",
    "    \n",
    "    \n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "test = '[X0 X1 Y2 Y3]'\n",
    "test = test.replace('[','').replace(']','')\n",
    "test.split(' ')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "identity = list('IIII')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['X0', 'X1', 'Y2', 'Y3']"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pauli_str[:-1].split(' ')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['X0', 'X1', 'Y2', 'Y3', '']"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pauli_str.split(' ')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "a = 10\n",
    "b  = a"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "b = b *2 "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "a"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "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.8.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}