examples/pyscf_demo.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Running pyscf to populate the MolecularData class\n",
    "\n",
    "The module run_pyscf.py provides a user-friendly way of running pyscf calculations in FermiLib. The basic idea is that once one generates a MolecularData instance, one can then call pyscf with a specification of certain options (for instance, how much memory to use and what calculations to do) in order to compute things about the molecule, update the MolecularData object, and save results of the calculation. The most common calculations users will want in FermiLib would probably be self-consistent field (aka Hartree-Fock calculations). Our code uses these \"SCF\" calculations compute orbitals, integrals, Hartree-Fock energy, and more. Other common calculations are CISD and FCI calculations which also compute the 1-RDM and 2-RDM associated with their answers, CCSD calculations which also compute the CCSD amplitudes (useful for UCC) and MP2 calculations which can also be used to UCCSD initial guesses. \n",
    "Note that the \"delete_input\" and \"delete_output\" options indicate whether to save automatically generated pyscf input files or not. To use this plugin, you will need to personally download pyscf.\n",
    "\n",
    "Warnings: electronic structure calculations are finicky. They sometimes fail for surprising reasons. See the pyscf documentation for more information or consult and electronic structure theory expert."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 0.275 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.465763692022 Hartree.\n",
      "MP2 energy of -0.471483763171 Hartree.\n",
      "FCI energy of -0.473370003213 Hartree.\n",
      "Nuclear repulsion energy between protons is 1.92428076698 Hartree.\n",
      "Spatial orbital 0 has energy of -0.817340000936 Hartree.\n",
      "Spatial orbital 1 has energy of 1.42093584882 Hartree.\n",
      "\n",
      "At bond length of 0.35 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.780454902072 Hartree.\n",
      "MP2 energy of -0.786933019271 Hartree.\n",
      "FCI energy of -0.789269392404 Hartree.\n",
      "Nuclear repulsion energy between protons is 1.51193488834 Hartree.\n",
      "Spatial orbital 0 has energy of -0.773812219107 Hartree.\n",
      "Spatial orbital 1 has energy of 1.26613647236 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 0.425 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.9500662942 Hartree.\n",
      "MP2 energy of -0.957467753563 Hartree.\n",
      "FCI energy of -0.960393757432 Hartree.\n",
      "Nuclear repulsion energy between protons is 1.24512284922 Hartree.\n",
      "Spatial orbital 0 has energy of -0.731224742484 Hartree.\n",
      "Spatial orbital 1 has energy of 1.1202060156 Hartree.\n",
      "\n",
      "At bond length of 0.5 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.04299627454 Hartree.\n",
      "MP2 energy of -1.05148606962 Hartree.\n",
      "FCI energy of -1.05515979447 Hartree.\n",
      "Nuclear repulsion energy between protons is 1.05835442184 Hartree.\n",
      "Spatial orbital 0 has energy of -0.690822328662 Hartree.\n",
      "Spatial orbital 1 has energy of 0.988673673002 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 0.575 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.09157049678 Hartree.\n",
      "MP2 energy of -1.10131707528 Hartree.\n",
      "FCI energy of -1.10591805325 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.920308192904 Hartree.\n",
      "Spatial orbital 0 has energy of -0.652936838321 Hartree.\n",
      "Spatial orbital 1 has energy of 0.873105751443 Hartree.\n",
      "\n",
      "At bond length of 0.65 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.11299654567 Hartree.\n",
      "MP2 energy of -1.12417429634 Hartree.\n",
      "FCI energy of -1.12990478432 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.814118786031 Hartree.\n",
      "Spatial orbital 0 has energy of -0.617605463235 Hartree.\n",
      "Spatial orbital 1 has energy of 0.773022995439 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 0.725 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.11734326912 Hartree.\n",
      "MP2 energy of -1.13013540113 Hartree.\n",
      "FCI energy of -1.13722137707 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.729899601269 Hartree.\n",
      "Spatial orbital 0 has energy of -0.584812015544 Hartree.\n",
      "Spatial orbital 1 has energy of 0.686877089471 Hartree.\n",
      "\n",
      "At bond length of 0.8 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.11085039748 Hartree.\n",
      "MP2 energy of -1.12545353399 Hartree.\n",
      "FCI energy of -1.13414766668 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.66147151365 Hartree.\n",
      "Spatial orbital 0 has energy of -0.554495881131 Hartree.\n",
      "Spatial orbital 1 has energy of 0.612618086735 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 0.875 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.0974543197 Hartree.\n",
      "MP2 energy of -1.11408436316 Hartree.\n",
      "FCI energy of -1.12467175357 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.604773955337 Hartree.\n",
      "Spatial orbital 0 has energy of -0.526513679939 Hartree.\n",
      "Spatial orbital 1 has energy of 0.54813298342 Hartree.\n",
      "\n",
      "At bond length of 0.95 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.07963692821 Hartree.\n",
      "MP2 energy of -1.09853448444 Hartree.\n",
      "FCI energy of -1.11133941774 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.557028643074 Hartree.\n",
      "Spatial orbital 0 has energy of -0.500643053408 Hartree.\n",
      "Spatial orbital 1 has energy of 0.491560510155 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 1.025 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.05898071013 Hartree.\n",
      "MP2 energy of -1.08041447981 Hartree.\n",
      "FCI energy of -1.09580487601 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.516270449678 Hartree.\n",
      "Spatial orbital 0 has energy of -0.476626271643 Hartree.\n",
      "Spatial orbital 1 has energy of 0.441437629683 Hartree.\n",
      "\n",
      "At bond length of 1.1 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.03653887503 Hartree.\n",
      "MP2 energy of -1.0608064532 Hartree.\n",
      "FCI energy of -1.07919294497 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.481070191745 Hartree.\n",
      "Spatial orbital 0 has energy of -0.454218694119 Hartree.\n",
      "Spatial orbital 1 has energy of 0.396695911145 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 1.175 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -1.01306099773 Hartree.\n",
      "MP2 energy of -1.04048753024 Hartree.\n",
      "FCI energy of -1.06231541668 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.450363583762 Hartree.\n",
      "Spatial orbital 0 has energy of -0.433218017521 Hartree.\n",
      "Spatial orbital 1 has energy of 0.356574452957 Hartree.\n",
      "\n",
      "At bond length of 1.25 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.989113814091 Hartree.\n",
      "MP2 energy of -1.02004941564 Hartree.\n",
      "FCI energy of -1.04578314455 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.423341768736 Hartree.\n",
      "Spatial orbital 0 has energy of -0.413471036649 Hartree.\n",
      "Spatial orbital 1 has energy of 0.320515554372 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 1.325 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.965138979586 Hartree.\n",
      "MP2 energy of -0.999955738808 Hartree.\n",
      "FCI energy of -1.03005661113 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.399379027109 Hartree.\n",
      "Spatial orbital 0 has energy of -0.394866384636 Hartree.\n",
      "Spatial orbital 1 has energy of 0.288081397752 Hartree.\n",
      "\n",
      "At bond length of 1.4 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.941480654708 Hartree.\n",
      "MP2 energy of -0.980569709566 Hartree.\n",
      "FCI energy of -1.01546824929 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.377983722086 Hartree.\n",
      "Spatial orbital 0 has energy of -0.377322825369 Hartree.\n",
      "Spatial orbital 1 has energy of 0.258901974587 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 1.475 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.918402055716 Hartree.\n",
      "MP2 energy of -0.962170870967 Hartree.\n",
      "FCI energy of -1.0022359731 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.358764210793 Hartree.\n",
      "Spatial orbital 0 has energy of -0.360778808505 Hartree.\n",
      "Spatial orbital 1 has energy of 0.232648634323 Hartree.\n",
      "\n",
      "At bond length of 1.55 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.896098958325 Hartree.\n",
      "MP2 energy of -0.944968709166 Hartree.\n",
      "FCI energy of -0.990476340937 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.341404652206 Hartree.\n",
      "Spatial orbital 0 has energy of -0.345185065649 Hartree.\n",
      "Spatial orbital 1 has energy of 0.209023237141 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 1.625 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.874712020072 Hartree.\n",
      "MP2 energy of -0.929114904673 Hartree.\n",
      "FCI energy of -0.980219639774 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.325647514412 Hartree.\n",
      "Spatial orbital 0 has energy of -0.330499875199 Hartree.\n",
      "Spatial orbital 1 has energy of 0.187754608434 Hartree.\n",
      "\n",
      "At bond length of 1.7 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.854337626951 Hartree.\n",
      "MP2 energy of -0.914713958611 Hartree.\n",
      "FCI energy of -0.971426688458 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.311280712306 Hartree.\n",
      "Spatial orbital 0 has energy of -0.316686045438 Hartree.\n",
      "Spatial orbital 1 has energy of 0.168597322543 Hartree.\n",
      "\n",
      "At bond length of 1.775 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.835036819455 Hartree.\n",
      "MP2 energy of -0.901831833626 Hartree.\n",
      "FCI energy of -0.964006402031 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.298128006152 Hartree.\n",
      "Spatial orbital 0 has energy of -0.303708862 Hartree.\n",
      "Spatial orbital 1 has energy of 0.151330576934 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 1.85 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.816842292251 Hartree.\n",
      "MP2 energy of -0.890502698872 Hartree.\n",
      "FCI energy of -0.957832967869 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.286041735632 Hartree.\n",
      "Spatial orbital 0 has energy of -0.291534608794 Hartree.\n",
      "Spatial orbital 1 has energy of 0.135756534497 Hartree.\n",
      "\n",
      "At bond length of 1.925 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.799763828147 Hartree.\n",
      "MP2 energy of -0.880734213989 Hartree.\n",
      "FCI energy of -0.952761485739 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.274897252426 Hartree.\n",
      "Spatial orbital 0 has energy of -0.280129531526 Hartree.\n",
      "Spatial orbital 1 has energy of 0.121698234354 Hartree.\n",
      "\n",
      "At bond length of 2.0 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.783792654277 Hartree.\n",
      "MP2 energy of -0.872511876467 Hartree.\n",
      "FCI energy of -0.948641112176 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.26458860546 Hartree.\n",
      "Spatial orbital 0 has energy of -0.269459223656 Hartree.\n",
      "Spatial orbital 1 has energy of 0.108997369549 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 2.075 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.76890514493 Hartree.\n",
      "MP2 energy of -0.865802852139 Hartree.\n",
      "FCI energy of -0.945325124834 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.255025161889 Hartree.\n",
      "Spatial orbital 0 has energy of -0.259488420123 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0975121876015 Hartree.\n",
      "\n",
      "At bond length of 2.15 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.755066140849 Hartree.\n",
      "MP2 energy of -0.860559520937 Hartree.\n",
      "FCI energy of -0.942677785696 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.246128935312 Hartree.\n",
      "Spatial orbital 0 has energy of -0.25018114306 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0871156542694 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 2.225 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.742231993159 Hartree.\n",
      "MP2 energy of -0.856722791678 Hartree.\n",
      "FCI energy of -0.940578294125 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.237832454346 Hartree.\n",
      "Spatial orbital 0 has energy of -0.241501099769 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0776939126166 Hartree.\n",
      "\n",
      "At bond length of 2.3 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.730353321355 Hartree.\n",
      "MP2 energy of -0.854225123201 Hartree.\n",
      "FCI energy of -0.938922385987 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.230077048226 Hartree.\n",
      "Spatial orbital 0 has energy of -0.23341221005 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0691449985036 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 2.375 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.71937741672 Hartree.\n",
      "MP2 energy of -0.852993149483 Hartree.\n",
      "FCI energy of -0.937622232412 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.222811457229 Hartree.\n",
      "Spatial orbital 0 has energy of -0.225879143757 Hartree.\n",
      "Spatial orbital 1 has energy of 0.061377742722 Hartree.\n",
      "\n",
      "At bond length of 2.45 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.709250219239 Hartree.\n",
      "MP2 energy of -0.852949830064 Hartree.\n",
      "FCI energy of -0.936605256416 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.215990698335 Hartree.\n",
      "Spatial orbital 0 has energy of -0.218867775452 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0543107898459 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 2.525 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.699917829065 Hartree.\n",
      "MP2 energy of -0.854016106823 Hartree.\n",
      "FCI energy of -0.935812374184 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.209575133038 Hartree.\n",
      "Spatial orbital 0 has energy of -0.212345500745 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0478716815187 Hartree.\n",
      "\n",
      "At bond length of 2.6 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.691327561197 Hartree.\n",
      "MP2 energy of -0.856112114321 Hartree.\n",
      "FCI energy of -0.935196030847 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.203529696508 Hartree.\n",
      "Spatial orbital 0 has energy of -0.2062813977 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0419959755719 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 2.675 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.683428596297 Hartree.\n",
      "MP2 energy of -0.859158041179 Hartree.\n",
      "FCI energy of -0.934718271951 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.197823256419 Hartree.\n",
      "Spatial orbital 0 has energy of -0.200646248214 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0366263937864 Hartree.\n",
      "\n",
      "At bond length of 2.75 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.676172310463 Hartree.\n",
      "MP2 energy of -0.863074762874 Hartree.\n",
      "FCI energy of -0.934348987747 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.192428076698 Hartree.\n",
      "Spatial orbital 0 has energy of -0.195412454328 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0317120058829 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 2.825 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.669512378614 Hartree.\n",
      "MP2 energy of -0.867784361057 Hartree.\n",
      "FCI energy of -0.934064392349 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.187319366697 Hartree.\n",
      "Spatial orbital 0 has energy of -0.190553892381 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0272074643943 Hartree.\n",
      "\n",
      "At bond length of 2.9 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.663404741599 Hartree.\n",
      "MP2 energy of -0.873210617868 Hartree.\n",
      "FCI energy of -0.933845750779 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.182474900317 Hartree.\n",
      "Spatial orbital 0 has energy of -0.186045745933 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0230723055488 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 2.975 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.657807511304 Hartree.\n",
      "MP2 energy of -0.87927953613 Hartree.\n",
      "FCI energy of -0.933678338578 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.177874692746 Hartree.\n",
      "Spatial orbital 0 has energy of -0.18186434986 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0192703273508 Hartree.\n",
      "\n",
      "At bond length of 3.05 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.652680866726 Hartree.\n",
      "MP2 energy of -0.885919898164 Hartree.\n",
      "FCI energy of -0.933550604948 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.173500724892 Hartree.\n",
      "Spatial orbital 0 has energy of -0.177987066597 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0157690500134 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n",
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "At bond length of 3.125 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.647986972272 Hartree.\n",
      "MP2 energy of -0.893063845078 Hartree.\n",
      "FCI energy of -0.933453506176 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.169336707494 Hartree.\n",
      "Spatial orbital 0 has energy of -0.174392204073 Hartree.\n",
      "Spatial orbital 1 has energy of 0.0125392576825 Hartree.\n",
      "\n",
      "At bond length of 3.2 angstrom, molecular hydrogen has:\n",
      "Hartree-Fock energy of -0.643689931023 Hartree.\n",
      "MP2 energy of -0.900647439214 Hartree.\n",
      "FCI energy of -0.933379977349 Hartree.\n",
      "Nuclear repulsion energy between protons is 0.165367878413 Hartree.\n",
      "Spatial orbital 0 has energy of -0.171058975462 Hartree.\n",
      "Spatial orbital 1 has energy of 0.00955461527119 Hartree.\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warn: Ipython shell catchs sys.args\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPX1+PH3yb5vJCwBwhJBVATBSBUVFLW1WkWtIu62\nfkXbKnbR1i62tv3Zalu1VaoWd+tScau4CyhuVDaBiLJG2cMSyL4nc35/3EkYhplkQjKZmeS8nidP\nZu69M/dMBubMZ7nnI6qKMcYYE6ioUAdgjDEmsljiMMYY0yGWOIwxxnSIJQ5jjDEdYonDGGNMh1ji\nMMYY0yGWOIwxxnSIJQ5jjDEdYonDGGNMh8SEOoBgyM7O1qFDh4Y6DGOMiRjLly8vUdWcQI7tkYlj\n6NChLFu2LNRhGGNMxBCRzYEea11VxhhjOsQShzHGmA6xxGGMMaZDLHEYY4zpEEscxhhjOsQSR4vC\nOXDvaLg9w/ldOCfUERljTFjqkdNxO6xwDrw2ExprnfvlW537AGOmhS4uY4wJQ9biACrf+O3+pNGi\nsdbZbowx5gCWOICU+p0d2m6MMb2ZJQ5A0gd1aLsxxvRmljgATvstxCYeuC020dlujDHmAJY4AMZM\nY/2EO6jVOBTYQTbrJ9xhA+PGGOODJQ5gUVEJ0z/NY0HsKZRHZbHpisVM/zSPRUUloQ7NGGPCjiUO\noHBbObMuHYcm9SHVVcbEYVnMunQchdvKQx2aMcaEHbuOA7h+cj4A/03MJrrCBXVlTMzPZmJ+dogj\nM8aY8BOSFoeIZInIPBHZ4P6d6ee4TSLyuYisFJGgL7AhyU6i0Oo9wT6VMcZErFB1Vd0KLFDVEcAC\n931/TlXVY1S1INhBRaf2BaCufHewT2WMMRErVIljKvCk+/aTwHkhiuMA8Rn9AKjaWxziSIwxJnyF\nKnH0U9WWT+edQD8/xykwX0SWi8iMtp5QRGaIyDIRWbZnz6F1NSW5E0dt+a5DerwxxvQGQRscF5H5\nQH8fu37teUdVVUTUz9OcpKrbRaQvME9E1qrqh74OVNXZwGyAgoICf8/XptQ+TuJotK4qY4zxK2iJ\nQ1VP97dPRHaJyABVLRaRAYDPT2pV3e7+vVtEXgEmAD4TR1fISk2mTJNxVdnguDHG+BOqrqq5wFXu\n21cBr3ofICLJIpLachv4JrA6mEH1SY5nr6YhNXbhnzHG+BOqxHEncIaIbABOd99HRHJF5E33Mf2A\nj0VkFbAEeENV3w5mUIlx0ZRJOjF1e4N5GmOMiWghuQBQVfcCp/nYvgM4y337K2BsN4dGZUwGuQ02\nOG6MMf5YyREvtbFZJDeWhjoMY4wJW5Y4vDTGZ5HsqgBXc6hDMcaYsGSJw0tzYh+icUGttTqMMcYX\nSxxeJMVdr6rKruUwxhhfLHF4iUl1Xz1eZgPkxhjjiyUOL/FpTqHDqn07QxyJMcaEJ0scXpKynCop\nViHXGGN8s8ThJS2zHy4VGiusq8oYY3yxxOGlT1oipaTgqrKyI8YY44slDi9ZyXFWr8oYY9pgicNL\nQmw05ZJOrNWrMsYYnyxx+FAZk0F8w75Qh2GMMWHJEocP9XGZJDfZlePGGOOLJQ4fGuL7kOqqhOam\nUIdijDFhxxKHD65Ep+wINTbOYYwx3ixx+NBar6ralpA1xhhvljh8iEnNAaC61C4CNMYYb5Y4fEjI\ncAodVu8rDnEkxhgTfixx+JCUOQCA+nJrcRhjjDdLHD6kZebQpFE0VlihQ2OM8RaSxCEiWSIyT0Q2\nuH9n+jkuQ0ReFJG1IrJGRE7ojviyUxMpJRVXtZUdMcYYb6FqcdwKLFDVEcAC931f/gG8raqjgLHA\nmu4ILjM5lr2aRpTVqzLGmIOEKnFMBZ50334SOM/7ABFJByYBjwKoaoOqlnVHcPEx0ZRFpRNbZ2VH\njDHGW6gSRz9VbZmytBPo5+OYYcAe4HERWSEij4hIcncFWBOdQYLVqzLGmIMELXGIyHwRWe3jZ6rn\ncaqqgPp4ihhgPPCgqo4DqvHfpYWIzBCRZSKybM+ezl+4VxuXZfWqjDHGh5hgPbGqnu5vn4jsEpEB\nqlosIgMAX9OXtgHbVHWx+/6LtJE4VHU2MBugoKDAVyLqkMaELJJrq6GpAWLiOvt0xhjTY4Sqq2ou\ncJX79lXAq94HqOpOYKuIHO7edBrwZfeE51mvygbIjTHGU6gSx53AGSKyATjdfR8RyRWRNz2OuxF4\nRkQKgWOAP3VXgFEpTtkRV5XVqzLGGE9B66pqi6ruxWlBeG/fAZzlcX8lUNCNobWKSesLQE3pTlJy\nx4YiBGOMCUt25bgfLfWqqqzQoTHGHMAShx/J7sRRX2aJwxhjPFni8CM9M4dGjaap0hKHMcZ4ssTh\nR5/UePaRilq9KmOMOYAlDj8yk+LYq+lWr8oYY7xY4vAjLibK6lUZYyJD4Ry4dzTcnuH8LpwT1NNZ\n4mhDTUwGCY1WdsQYE2JtJYbCOfDaTCjfCqjz+7WZQU0eIbmOI1LUxWWRXNstBXmNMb1d4RxY8Aco\n3wbpg+C038KYafsTQ2Otc1z5Vnj1Btj0EcSnwdJHoKnuwOdqrHWea8y0oIRqiaMNTQlZJNXUOG9C\nbGKowzHG9FS+ksPcmVC6GT59YP/2Fs318NlTEB3v3PalfFvQwrWuqjY0JzplR7CZVcaYzmqru2ne\n7w5ODk218P7/g1p/46wCvy6G9MG+d6cP6pKwfbHE0YboVKfQoavKEocxph0dHYf47w/ggRPhryOg\ncof/503p73t7+iCIina6tLx7RGITne1BEnDiEJGkoEURpmJSnXpV1aXF7RxpjOnV/A1QL3kY1r0F\nb/zs4BaFqwlK1sGIMyAhw/fzpg+Gb/6x7cQwZhqcc5+75SHO73PuC9r4BgQwxiEiE4FHgBQgT0TG\nAtep6g+DFlWYSHSXHakp3UVqiGMxxoSYv8FrcLZ7J4bGWnjz5raf09UE5z1w8BgH7E8OnufwdW5w\nbgcxUXgLZHD8XuBbOGtooKqrRGRSUKMKE8mZThOxrtzKjhjTq/kcvL4Rti8HiXa3NPz4/jvw4jVQ\n4WOwumUcor3k0M2JoT0BzapS1a0i4rmpOTjhhJeMzCzqNZamSluTw5ger6MtiqY6WPyQM7MpOg6a\nGw5+zvTBkHc8nP47/y2KFmGWHNoSSOLY6u6uUhGJBW4C1gQ3rPCQlRLPXlLRKl8r2xpjegx/LYpN\nH0FTfRstCoFfboUvX207MQTS3RRBAkkc1wP/AAYC24F3gR8FM6hwkZUUxxpNI6V2b6hDMcZ0BX+t\nCn8tis+eguQcJwl47wfnOWLiw3IcIpjaTRyqWgJc1g2xhJ2Y6CgqotLJrLd6VcZEPH9XYBfOabtF\ncfMG+PyFHtXV1FntTscVkZEiskBEVrvvjxGR3wQ/tPBQHZNJYoMlDmMinq9WRXM9bJznjFP4kj4I\nREIy5TWcBdJV9TBwC/AvAFUtFJFngf8XzMDCRb3VqzImcnh3RZ36K0jLhXVvt92qmDrLWhQdEEji\nSFLVJV6zqpqCFE/YaUzoQ0JNPTRUQ1xyqMMxxvjjqyvqvz9wbkfHQ0zCwcUAwUkwPWzwOtgCSRwl\nIpIPKICIXAh06lJqEckCngeGApuAaapa6nXM4e5jWgwHfquqf+/MuTssqQ/sw6lXZYnDmNDyN7hd\nsQPevMX3AHZSNvy4ENa+0f7MJ0sUAQmk5MiPcLqpRonIduDHODOtOuNWYIGqjgAWuO8fQFXXqeox\nqnoMcCxQA7zSyfN2mKQ4hQ6bq+xaDmNCyme9px/CP46Be46AOj9dyjV7nS99Nk7RZdpscYhIFFCg\nqqeLSDIQpaqVXXDeqcAp7ttPAguBX7Rx/GlAkapu7oJzd0hculN2pGpfsd8ilMaYbuBrcNvV6LQ+\npvwGljwCVTsPfpxnlVhrVXSJNlscquoCfu6+Xd1FSQOgn6q2dHftBPq1c/x04Lm2DhCRGSKyTESW\n7dnTda2DeHfiqC2zsiPGBJ2vCrNVe5xigf4Gt11NMOmW9osBmi4TyBjHfBG5GWe8obplo6q2OUdV\nROYDvuoB/9rzjqqqiGgbzxMHnAv8sq3zqepsYDZAQUGB3+frqJQs5yXUW+IwJrh8DW6/ch2oAgpR\nsU4Lw1ug9Z5MlwkkcVzs/u15tbjiDFb7paqn+9snIrtEZICqFovIAKCtmh7fBj5T1ZB8cmdmZFKr\ncTTZGIcxwbXg9wd3RakL4lOdQoG7vrAps2EikMRxhKoeMIdNRBI6ed65wFXAne7fr7Zx7CW0000V\nTH1S4thLmtNcNsZ0jq9ZUYMnwMpn/S91Wl8F/Y5yfsBaFGEgkMSxCBgfwLaOuBOYIyLXAJuBaQAi\nkgs8oqpnue8nA2cA13XiXJ2SmRTHak0jzepVGdM5fruiXIA4NZ+afKyfbYPbYcdv4hCR/jiFDRNF\nZBzQcgVgGtCp1QBVdS/OTCnv7TuAszzuVwN9OnOuzoqOEiqiMsiut8RhTKf4mhWlLohPgx8sgi3/\na78ryoSFtloc3wKuBgYBd7M/cVQAvwpuWOGlOjaTxMY2Fmoxxux3UNmPX0Nsgv9ZUfWVkDHY+QHr\niooAfhOHqj4pIv8GLlHVZ7oxprDTEJdJSnWZM7vjwNIrxhhPPst+uK8XlmhQH2vAWVdUxAnkOo6f\ndFMsYaspsQ9xNDrfjIwx/vnqjgKn7Md5D9h1Fj1EICVH5ovIzSIyWESyWn6CHlkY0aRs50ZNSWgD\nMSZcNTfB6pf8d0fV7IWx063kRw8RtOs4epIod72qpsrdxGT1mpdtjG+eYxhpuZA3EbZ+6iSNqBjn\nSm5vnhfpWaKIeIGsADisOwIJZ3FpfQGo2reTjCEhDsaYUPIew6jYDqtfgD4jYfpzTnfu6zfZzKge\nLpAWByIyGjgSaL3wT1WfClZQ4SYhwyk7Ulu2i4wQx2JMSM27zfcYRlMtjHLPpBexmVE9XLuJQ0R+\nh1PJ9kjgTZwSIB8DvSZxtNSraii3elWmF/B1dXe/o+Cju6HSR/VZOPCqb+uO6vECaXFcCIwFVqjq\n90SkH/B0cMMKL1npaVRpAk2VVnbE9HA+r+6+3plGG5fi1I3yNbvQc0qt6fECmVVV656W2yQiaTgF\nCXvVyhRZyXHs1TSotsRhejifV3c3O1d3//hzOPsem1JrAmpxLBORDOBhYDlQBfwvqFGFmYykOLaQ\nRmZtm5XkjYl8fgsNVkJSlpUuN0Bgs6p+6L75kIi8DaSpamFwwwovLfWq+lm9KtMT+BrDyB0H79+B\nM9PeB7u623hoq8ih3+q3IjJeVT8LTkjhqTY2g8TGr0MdhjGd094YxuHfgaIFziypFtYVZby01eK4\n2+P2sTjdVC0UmBKUiMJUXXwWqVVWr8pEOH9jGHGpcNNKSM723SKxFobx0FaRw1NbbovICs/7vVFz\nYh9iqpqhrgwSM0MdjjGHxt8YRkOVkzTAuqJMuwKZVQV+Oz57j0Eu9/z1u4bBvaOdb2XGRAqXy/k3\nK37+y9t0WtMBAV053usVzuHYfa+776jTL/zaTOeufTMz4ca7q2nsJbBxPuz4DNLzoGoXNHustGdj\nGKaD2hocv5/9LY1BInKf535VnRnMwMJJ5Ru/JVUbD9zYWOtst8Rhwomvwe8P/wLxGXD+v+DoabD6\nRRvDMJ3SVotjmcft5X6P6gVS6n2XWfC33ZiQ8bceRnyKU9YcbAzDdFqbKwB2ZyDhTNIH+VxnQKxf\n2IQTVf/rYVRs795YTI8W6OB4l3IvBjVPRDa4f/ucpiQiPxGRL0RktYg8JyIJvo4LutN+iyvGyiyY\nMFayEZ6+wP9++5JjulBIEgdwK7BAVUcAC9z3DyAiA4GZQIGqjgaigendGmWLMdNYd9wdVGs8Cuwg\nm/UT7rDmvgmNwjnOzL7bM+CeI+G5S+CB42HbMhgz3WpJmaALVeKYCrR0hT0JnOfnuBggUURigCRg\nRzfEdpBFRSVctiSPR2IuQYBtF73F9E/zWFRkS8mabtYy+F2+FVCnC2rdmzCwAG5cDhf8y5ZnNUEX\nyHocOcC1wFDP41X1+504bz9VLXbf3gn08z5AVbeLyN+ALUAt8K6qvtuJcx6ywm3lzLp0HEveXQ+7\nYELKXmZdOo7CbeVMzM8ORUimt/I3+F2xDVKclSpt8NsEWyDXcbwKfATMB5oDfWIRmQ/097Hr1553\nVFVF5KALDN3jHlOBYUAZ8IKIXK6qPtcCEZEZwAyAvLy8QMMMyPWT8wFY0X8U7ALXnvVMLDjBkobp\nXi6X/8Fvf1eEGxMEgSSOJFX9RUefWFVP97dPRHaJyABVLRaRAThrfHg7HfhaVfe4H/MyMBE/i0ip\n6mxgNkBBQUFQrnTPGTicupWxNGz/krSCYJzBGD9KNu6/6NQXG/w23SiQxPG6iJylqm924XnnAlcB\nd7p/v+rjmC3A8SKShNNVdRoHXlvS7Yb3TedrHUDOrrWhDMP0Bp5XfyekQX2VU712/JXw+QsHdlfZ\n4LfpZoEMjt+EkzxqRaRCRCpFpKKT570TOENENuC0LO4EEJFcEXkTQFUXAy8CnwGfu2Od3cnzdsrw\nnBSKNJfY0qJQhmF6Ou8B8LpyZ/tpt8G599vgtwm5QBZySu3qk6rqXpwWhPf2HcBZHvd/B/yuq89/\nqLKS49gePYjU2iXQWAexobmsxPRwC37vu/T5J/+ACdfa4LcJubZqVY1S1bX+FnTqbQs5tahJyyeq\nwgX7voJ+R4Y6HNPT7F7rf6DbBsBNmGirxfFTnFlKd/vY1+sWcmohOSOgAihZb4nDdB2XCz59wBnX\nkChQ18HH2AC4CRNt1aqa4f7dqxdw8pY68AgogrritSQcFepoTMTyHPxO7e+swLd3PRx+NuRPgXm/\nsQFwE7ZsPY4OyuufwzbNJql4DTbCYQ6Jd+nzymKgGMZfDef83VmaOCHNSp+bsGWJo4Pyc5L5yjWA\n0SUbQh2KiVT+rv4uWrB/PXsbADdhLFS1qiLW4KwkvmYgKZVfOWWsjekoG/w2Ea7dxCEiL4vI2SL+\nFivuXWKjoyhNGkqcqxYqQlJz0USq5kZY8Ef2L6zpxQa/TYQIJBk8AFwKbBCRO0Xk8CDHFPaasw5z\nbpSsD20gJnKUbobHz4KP/gZDTgRb38VEsEAuAJwPzBeRdOAS9+2twMPA06rei3H3fPEDjoBiaN6z\njuh8m3RmfPCcNZWUBfU1EBMLFz4Go7974H4b/DYRJqDBcRHpA1wOXAGsAJ4BTsKpM3VKsIILV/0G\n5FGhiej2NaSHOhgTfrxnTdXsda7NOOU2J2mADX6biBbIGMcrOGXVk4BzVPVcVX1eVW8EUoIdYDjK\n75fKV5pL0651oQ7FhCNfs6bUBZ8+GJp4jOligbQ47lPV933tUNVeWVw8PzuF+ZrLYeWWOIwPtmaG\n6eECSRyZInKB17Zy4HNV9bWORo+XnhRLcexgUuo/gvpKiO/yOpAmEjU1wLzb/O+3WVOmhwgkcVwD\nnAC0tDpOAZYDw0TkD6r67yDFFtZq0/OhFCjZAAN91oE0vUn5Nnjhati2FPJPhy2fWMkQ02MFkjhi\ngSNUdReAiPQDngK+AXwI9MrEEZ1zuCWO3sxzVlRyNjTUOFd9X/QEHHW+zZoyPVogiWNQS9Jw2w0M\nVtV9ItLrpuK2yBg4kqZ1UTQWryFxbKijMd3Ke9ZU9R5A4Iw/OEkDbNaU6dECuQBwoYi8LiJXiUjL\nMq8LRSQZKAtueOFrWL9MNms/aovXhDoU09181ppSWBLSBSqN6TaBtDh+BFyAc90GON1UL6mqAr32\n6rf8nBTWaS4T9m4MdSimu1mtKdPLtZk4RCQamO9ek+Ol7gkpMgzMTORtyeXU6kJoboJoKzTcK6x7\ny/8+mzVleok2u6pUtRlwucuNGA/RUUJl8nBitBHKNoc6HBNsLhd88Fd47hJIHwwxXqux2Kwp04sE\n8jW5CvhcROYB1S0bVXXmoZ5URLKA54GhwCZgmqqW+jjuJuBaQICHVfXvh3rOYHD1OQy24cys6pMf\n6nBMV/KcFZWWCyl9YccKGHMxnPMPWPOazZoyvVYgieNl909XuhVYoKp3isit7vu/8DxAREbjJI0J\nQAPwtoi8rqphM6iQ2H8UbIOm3euIOfzMUIdjuor3rKmK7c7P6Ivg/H85025t1pTpxQKpjvukiCQC\nearaVTU2prK/OOKTwEK8EgdwBLBYVWsAROQDnEH6v3RRDJ02cGAue5amEbfjSyt22JP4W6Fv66f7\nV+gzphcLpMjhOcBK4G33/WNEZG4nz9tPVYvdt3cC/Xwcsxo4WUT6iEgScBYwuJPn7VL5OSl8pbm4\ndtu6HD2KzZoypk2BdFXdjtNdtBBAVVeKyPD2HiQi84H+Pnb92vOOqqqIHLQkmqquEZG7gHdxxlZW\nAs1tnG8GMAMgLy+vvfC6xLDsZF5z5TKmfHm3nM90A1czxCVDQ9XB+2zWlDFAYImjUVXL5cAmuqu9\nB6nq6f72icguERmgqsUiMgDnanRfz/Eo8Kj7MX/CGYr2d77ZwGyAgoKCblkMPDUhll3xeSQ2vQfV\neyG5T3ec1gRLXQW8dI2TNKJiwNW0f5/NmjKmVSCJ4wsRuRSIFpERwExgUSfPOxdnEag73b9f9XWQ\niPRV1d0ikoczvnF8J8/b5RoyDoO9OMvIJp8Q6nBMR3jOnErtDwhU7YKz73EqHtusKWN8CiRx3IjT\nvVQPPAe8A/yxk+e9E5gjItcAm4FpACKSCzyiqme5j3vJvfpgI/AjVQ27EicxfUfCXtCS9cgQSxwR\nw3vmVKV7yO3En8Bx1zi3LVEY41Mgs6pqcBLHr9s7NlCquhc4zcf2HTiD4C33T+6qcwZLn4H51H0Z\ni2vHGpKODXU0JmD+Zk6tfhHOuL3bwzEmkrSbOERkJHAzzsV6rcer6pTghRU5hvdN52sdQP+da0kK\ndTAmcDZzyphDFkhX1QvAQ8AjtDGrqbcanpPMSs1lUGnYXJdo2tPU4Ax2N9YcvM9mThnTrkDKqjep\n6oOqukRVl7f8BD2yCJGbnshmGUhyzXZorAt1OKY9deXw7EVO0ojy+t5kM6ci2kMfFLGoqOSAbYuK\nSnjog6J293fmscHe39nnDoZAEsdrIvJDERkgIlktP0GLKMJERQlVqcOJwgX7vgp1OKYt5dvgsTNh\n08cw9QE470GnYCHi/D7nPhsQD6Jgf/iOGZTODc+uaD1mUVEJNzy7gjGD0tvd35nHBnt/Z587GMRZ\nVqONA0S+9rFZVbXdiwBDpaCgQJctW9Zt57vzsee5dcsMuOhJOOq8bjuvaYfndNuUvvsHw6c9Bfm9\ndimZTnnogyLGDEpnYn5267ZFRSUUbivn+sn5be5v+YCbdek4JuZnt37AtXv/knEcOzSTjzaUcPML\nq/jD1NGMGZjO0k37+OPrX3LLtw5n1IA0GptcrNpWxj/fL2LSiGw+2FDClScMYVh2Mk3NSpPLxYZd\nVbz02TaOGZzBii1lnH30AHIzEmhyKVv21TD/y10c3j+VtTsrOXlENn3TEnC5lGaXsrOijsVf72No\nnyQ2ldQwLi+DrOQ4XKq4FPZW1fP59nIGpCdSXF7LqP6ppCc6+1WhtKaBjbur6Jsaz+7KeobnJJMc\nH4MqKFBZ28jmfTVkJceyr7qRQRmJzn6gqq6RHWV1ZCTFUlbTwICMBBJjnX0ANfVN7Kyoo19aAvVN\nrta/YUeIyHJVLQjo2PYSRyTq7sRx/9uruPHTSTRO/hWxp3qX3DIh4T3dFgBxuqJO/mnIwgq1znzw\nXz85v80P+2OHZPLemt3c+vLn/OybIxmencJnW0r514dFTCsYTL+0BNbtrOCtz3eS3zeFDbuqGDM4\nnaS4GOoamqltbGZvdT07y+uIj4mirtFFdJTQ5Ar+Z1RMlBAdJTS7lCaXEh8TRXJ8DFEiRImzjEKU\nCJV1jVTUNZGRGEt2ajxRAlEiiPu4PZX17K6sp19aPLkZic4+nGMQ2F5aw/ayOgZlJjK0T3Jr6TPn\nOeDrkmo2761haJ8kDuubglMY3CmRVrS7iq9KqsnPSWZkv9TW7eI+Zt2uSjburmLmlMP46TcP7/Df\noEsSh4j8XFX/4r59kaq+4LHvT6r6qw5H1k26O3HMXbWD8S+fTOrIk0m/7IluO69pw72joXzrwdvT\nB8NPVnd/PN2kMx/8nvfvuuBohuWk8OH6Pdwzbz3TJwwmMymOfdUNrN9Zyadf76VPchy7K+tJTYil\nrrGZ+qZ2C0og4nxINzYr6YkxDEhPJCE2moTYKBJjo0mMi+arPdWs3VnJ2EHpnJCfTVxMFPEePwvX\n7WHB2t18e3R/zh83kNjoKPePsHZXJXe/s47vjB3AG4U7uf2cIzluWBYxUVHERAsrNpfy85cKuWRC\nHs8t2cKsS8Zz4ojs1r/TDc+u4PJv5PH04i0HfWsP5f7OPncguipxfKaq471v+7ofbro7cWz94Aly\n3ruZeGlE0gfbVcbh4PYMwNe/bYHbw+460i7TXmJoaHLx5ufF3Pbqar4xrA+fbCxh0shsokTYVVHH\n7krnG7+/b/nJcdFkJsfR0ORid2U9h/VN4bihmaQlxJKWGEtaQgxpibHMX7OL11YVM/24wVw3OZ/k\n+GiS42JYuaWMG/8T3A/XDneFXToO4JAfG+z7nY0tUF2VOFao6jjv277uh5tuTRyFc9C5M5Emjy6R\n2EQbaA2l5ka4cwg0Vh+8L8JbHO21KFSVd77YyS0vFjJhaBafFJVQMCSLhiYX20pr2FlRh3dOSEuI\noW9aAv3S4umbmkDftHi+3FHBRxtKOO+YXK6dNJys5Dgyk+JIiI0+5A/3YH/YdqYbDuhUF14w93c2\ntkBZi6M7E0cv7RIJWw3V8MLVsOFd34UKIyChBzLAfN/0ceRmJPDG58U88H4RE4ZlUlnXRNGeaspr\nGw94voEZiQzKTGRQZhKDMhOpa2zm2SVbOH/cQF5btYN/Xja+W77Vt8QfrA9f0zldlTiaccqZC5AI\ntFwtJUCCqsZ2QaxB0Z2JQ2/PQHx0iSiC9OAukbBUXQLPXATFK51ChXHJEVmo0PvDd+G63cx8bgUX\nHzeYukZXCsxCAAAalElEQVQXn35VwobdB7amclLjyc9JJj8nhSgR/rtyOxeMG8hcr8QQym/19uEe\n3mxWVTcmjrq/HkFC9Y6DtyfnknDLmm6Jodfyrm7raoL6KrjwMRh1VvuPD5G2PlyvmzScoj1V/Gfp\nVp7632ZS42PYW93QelxqfAxH5KZR39jMqm3lTCsYxG++cyRpCbGtz2Mf/OZQWOLo5jGO5ldnEt28\nf4yjOTqB6Kn3R8S324jlc7otMPlWOPWXoYkpQN5TWJ9bvIW73l7HEQNS+bqkmtIap6spIdaZklow\nJJPvnzSMo3LTGJyZxKdf7/XblWQf/OZQdSRxBFKryrRlzDSigdLXfkNGwy5EIPrEmZY0gs1fdduV\nz4Q8cbT14T3j5OGkJcRyxhH9uPLRJa0XjwGU1TRyxpH9KBiSRXQU3PHmWmac7CSHjKRYhvRJPqgF\ncXx+nwPu+0oOE/OzOzw105i2WOLoAouSp3BD8z+JcZWzKOpadu6rxErlBVkYV7f1vkL67dXF/OyF\nVYwfnMkjH31FSZXT9ZSTEseeqgbOPnoAv596FNkp8cDB3U2eyaFwW/kBLYyJ+dmt2y05mO5iiaOT\nWv+TXzaOOUu3Urh2JImr32bRuJvtP3IwJaQ5BQu9hUF124n52fzy26O45ollpCbEsLuyHoAviis4\neUQ2k0fmkBgbza//u5qZUw7j6cVbWL+rsjVxtJUcrEVhwoEljk7y/E++s7yO+Z+P4eexc3iyqMj+\nMweDywXzbnOShkSDelT676bqtj67ojaWsGDtblLiY3jz82I27K4CoLaxmeOHZfGrs49gdG46UVFi\n3U0m4gVSHde04frJ+a3/oU8ekcMHrmMAuCrH1ufocs2N8N/r4X+zYMKMkFW3ba1GurGEtTsr+Mnz\nK7j80cU8+vHX3PfeBrKS47h64lAyEmOZOeUw1u+uoqq+iagop6ZQWy0KYyKBtTi6UE5qPDLgaErL\nMsncOA/GXRbqkCKb53TbtFxISIfdX8KU2+DknzmFj8Ze3O1hHTkgje+MGcAVjy6h2T0r8YgBqVw6\nIY9vje7Pxt1V3PDsCh64fLy1KEyPZImji00a2ZcFn4zhuxvfQ5qbINr+xIfEe7ptxXbnZ/yVMOnm\noJ7aV1fUJxtKmFu4g+r6Jt79YhcNza7W8tjXnDSM275zZOuxL3+23QawTY9mXVVdbPLIHN5vHovU\nl8P27iu02OP4m25b9H7QT+25MM72slpunrOKKx5bzPNLt/LRhhIu/UYefz7/aJpcyswph/HKiu0H\nLDDk2X3Zwl9Lw5hIFJKvwyJyEXA7cAQwQVV9fsKKyJnAP4Bo4BFVvbPbgjxE44dk8pPYY3ARRdSG\neZB3fKhDikwhnG57wvA+3HDqYVz12BIam52uqKNy05gxaTjfOqo/n20pbXNw25ieLlT9KKuBC4B/\n+TtARKKBfwJnANuApSIyV1W/7J4QD01sdBSj84fw+dcjGbNxHnLabaEOKTIlZ0P1noO3d9F0W1/d\nUR9t2MMLy7ayaW8NhdvKiY+JApSrThjC76eObj3OrqUwvV1IuqpUdY2qrmvnsAnARlX9SlUbgP8A\nU4MfXedNPjyHdxvGIMWroHJXqMOJPF++CjX7aFn9rFUXTrf17I4qr2nkly8XctVjS5i7qpiquia+\nf+JQkuKimTnlMF4rLLauKGM8hPPI7UDAs175NuAbIYqlQyaNyOF611huYQ4ULYBjLg11SJFj6SPw\nxs0weAKMnQ4f3ROU6rYT87P5/blH8b3Hl7YuF3pUbho/++ZI4qOjufE/K1qrylpXlDEHClriEJH5\nQH8fu36tqq8G4XwzgBkAeXl5Xf30HTI4K4m6PkdSVpNJxoZ5ljj88Zxumz4I+o+Fda/DyG87FW7j\nkqDg+11+2h1ltTywcCNzlm6jyeXCpXDJhMH8+YIxgNONZV1RxvgXtMShqqd38im2A4M97g9yb/N3\nvtnAbHCq43by3J026fB+vLf0aM4vsmm5PnlPty3f6vwMOREufrrTfy9fYxivrdrBY598zertzoV2\nk0bksHxzKVeeMISnF2/hnLEldp2FMQEI5+m4S4ERIjJMROKA6cDcEMcUsEkjc3ivaSxSVwbbl4c6\nnPDjb7pt2ZYuSbKeYxg7ymr5vyeXcuNzKyjcVsa0gsHcPe0YVmwt44HLx/PTbx7OrEvHtR5vjGlb\nSBKHiJwvItuAE4A3ROQd9/ZcEXkTQFWbgBuAd4A1wBxV/SIU8R6K44f1YXHUWFxEwcZ5oQ4n/AR5\nuu3E/GzuOG8033t8KSfd9R7z1+zmtFF9+fDnU7jj/KPZUVZrZT+MOUQh6T9R1VeAV3xs3wGc5XH/\nTeDNbgytyyTGRTNqWB5f7hjJ6A3zYMpvQh1SeEntD5XFB2/vgum2NQ1NPPzh18z+sIiGJhcKXD1x\nKLefe1TrMdYdZcyhC+euqog3eWQOb9cd7ayBXbU71OGEjx0rnCVevXVguu1DHxQd1K308YYSfvD0\nck7560Lunb+eI3PTSHMXGpy7aod1QxnTRSxxBNHkkTksdI117mxcENpgwsWG+fD42ZCYCWf84ZCr\n23qOYagqD7y/kaseW8xbq3cyKDOR2885kqI91TxoYxjGdDlbczyIVJUT/zyft5uvJe2IKc4U097G\nc8ptYgbUlkH/0XDZi053VScsKirh+n8vJzUhlu1ltfRLi+d35xzFt0f3518ffmVrbxvTAbbmeJgQ\nESYd3o/3C4/m3KL3EFczREWHOqzu4z3ltrYUJMpZS6OTSaO0uoE3CoupqGuioq6JSSOyeeSq44iL\ncRrRNoZhTPBYV1WQTR6ZQ2lzPFJbCn/oA/eOdj5QewNfU27VBR/85ZCfstmlPLN4M6fevZDnlmwh\nPiaKGZOGs3pHBcs27+tkwMaYQFiLI8gm1y9Eoj9w31PnIrfXZjp3u2G1upDqxJRbXxfwPbFoEw8u\n3MiuinpG9U9FVXnw8mOZmJ/NKYfnWFkQY7qJtTiCLOmjO0iUhgM3NtY638Z7sn1f+++WC2DKrefg\n957Keq58bDG3z/2ChiYX910yjqnH5LYmDbDrMIzpTtbiCLYQrisRMpv/B/+5FKLjQKKhuX7/vgCn\n3E7Mz+a+6eP4vyeX0dTsoqFZOXdsLn++4GiS433/s7UxDGO6h7U4gqwy3vcgsL/tEW/Vf+CpcyEp\nC67/GKbOOqQpt6u3l/OXd9ZS09BMQ7NyxfF53HfJOL9JwxjTfex/YZAVF9xCzCe3ksj+7qoGYigu\nuIXUEMbVZTyn28anQn0FDD0Zpj3lJI8++R0ay6mub+Keeet5/JOvSYmPJTk+hu+fOJRnFm/h20cP\nsBaFMWHAEkeQjTzjGtYDyR//iQHspVmiID6dkVOuDHVonec93ba+wumaOuZSJ2m0wdfg9/3vbWD2\nh19RWdfEaaP68tmWUh66zBnHOMHWxDAmbFhXVTcYecY1/O3Ilxhe/wz/yL6duPq9sPyJUIfVeT6n\n2zbD+39q96He1WsvfHARd7+7nozEWF76wQkcNyyrdSElsMFvY8KJXTneDRYVlXDDsytIiotme2kN\nn+XdR2ZVEdy00uneiVS3ZwC+/v0I3F7W7sM/3lDCjH8vo7HJRaNLmX7cYP543mhio+37jDHdrSNX\njtv/0CBrSRqzLh3H/ZeMQxFmFJ8LNSXwyX2hDu/QuFzw4d/wnTQIaLrt2p0V/O3dddQ0NNPoUq46\nYQh3fneMJQ1jIoD9Lw2ywm3lrf3y4/IyOf2Ifqwmn3XZZ8D/ZkGFj9Li4ay2DJ6/DN77Iww6DmIS\nD9zfznTbusZm/vrOWr5z38cU7akiOT6aG6ccxmuFxVaA0JgIYYkjyK6fnH/AYO7PvjmSuiYXCwdd\nD82NsPDPIYwuAIVznDIpt2fA30bCrALY8C6ceRdcMw/Ovc/vdFvv0uf/K9rLKX9dyD/fL+KE/D5E\nRwkPX1nAz6x6rTERxcY4QuCm/6zgnS928tlxC0ha8Sj88FPIOTzUYR3Me9ZUi8m/gFN/1e7DW7rp\n7rpgDPPX7OL5ZVuJEvjFmaNQsOq1xoSRjoxxWOIIgU0l1Zx2zwdcV5DOz9dNh6EnwSXPhTqsg907\n2qmt5S19MPxkdbsPV1XumbeeWe9tBCA+NoqHLjuWU0b17epIjTGdZIPjYW5odjLTCgbx8PJyyo/9\nEax7E/46wukOCqfquZ0ol7J1Xw1XP76U+9/bSE5qPArMOHm4JQ1jegBLHCFy45QRiAhvb3IBAtW7\nOaB6biiTR0M1vPlzDmXWVGOzi4c+KOKMez9g2aZ9XHnCEJqalZlTDuPpxVtsDMOYHiAkiUNELhKR\nL0TEJSJ+m0Yi8piI7BaR9vtFIkxuRiJXHD+EE7c/ykEf0N1ZPddz8Pve0fDeHfDgibDkX5B/Wpuz\nprwHvz/bUsqUvy3kzrfWMmlEDnddOIbXC4uZddk4W77VmB4kVC2O1cAFwIftHPcEcGbQowmRH5yS\nTy57fe/sjuq5LYPf5Vtpbe18+Beoq4Cr34ArXm5z1lTL1d/z1+zitv+u5oIHFrG9rJafnjGS2VcW\nsK209oASIXb1tzE9Q0hqVanqGnCWVm3nuA9FZGg3hBQS2SnxlCf0I71+58E7A7iIrtN8lQwBiE1w\nBuzBSRJ+ihQeP6wPFxcM5tqnlqEK8TFRPHj5eKaM6gfY8q3G9FQ2xhFii4b8kFqNO2j7Lle6M9YQ\nTP5aNRU72n1o4bYyvvvQIh78oIi+qfEAXDdpeGvSMMb0XEFLHCIyX0RW+/iZGqTzzRCRZSKybM+e\nPcE4RVCkf+MybtMZbHNlowh1ybm8yYn0rfoSHv0mlG7q3Am8xzAK58Cmj+Hxs2lv8Nt7DAPg7dXF\nnHv/x0z95yds3VfL9ZPzaWyywW9jepOgdVWp6unBem4/55sNzAbnOo7uPHdnTMzPpv7SmZz8xEQy\nk2Nx1cADV4xHdCW8+H2YfSocezV8/oLTQkgf5AxOB7LGhfcFfOVb4ZXrQF2Q0g/GTIcvX4Umj+4q\nj8HvljGMWZeO47ihWfzhtS95+tPNiMC1Jw/nG8OyuOXFQmZd5oxjHG+lz43pFWw9jjBw6qh+TD1m\nIP9duZ3UhBgGpCdC9ulw7fvwxNnw8T37D26Zrgv7k4fnYkqeicVn2XMXJGTATaucJHHYab4fi3sw\n+5JxzHhqObHRQmlNI2MGpXPPtGM4rG8KD31Q5Hfw2xKHMT1XSK4cF5HzgfuBHKAMWKmq3xKRXOAR\nVT3LfdxzwClANrAL+J2qPtre84f7lePeWkpznH5EX15Yto3UhBie/r9vMGZQBtxzJFRsP/hBiVlw\n+Yuw6wt46+cHJoioWBgwFrb7+xu0X/ZcVfloQwn3zl/Pii3Osd8ZM4D7LxnX7qQGY0zksZIjEZQ4\nPMuuT8zP5qXlW7nlxUJio6N4+MoCJj07Ar9jEW2RaIiJh8aag/e5S4b4WoVv0cYSXl21g427q1i+\nuZQ+ybHUNrr43sShPLd0q3VDGdNDWcmRCOJZdh3gu8cOZtYl40hLiOH7TyylJrG/7wem9oeLn/H/\nxOqCc/7hdEd58jGGsaioBFVl9odFXPnYEp5fupUdZbV8/8ShuBQeuaqAW84cZRfwGWMAa3GErYq6\nRq59chn9Ns/lb/GPEqf1rfuaoxN5f+SvOf3iG9svROhv/MPtg/W7+dEzK0hLiGFHeR2ZSbH89IyR\nTDtuMI9/sskq2BrTS1hXVQ9IHOAsenTlo4vpv+U1bk96icym3dQnD+APNRfynctvcj7QfZU+j02E\nc+7jodJjfX7wr9paxri8TF5duZ03CoupqGsCYPLIHGZfeSzxMdHd/VKNMSHWkcRhs6rCWEJsNM/N\nOIHr/h3L+DUnkZkUS215Mz8+YyTj8zIBeKj0WKZMuIORq+9tbVWsH/0T3nMnDc/xk+eXbuV3c1eT\nHBfN3upGkuKiGZ+XyaptZVx1whCeXbKV5ZtLbQzDGNMma3FEAFXl/55cxoK1u4kScLnLe0wYlkVe\nVhKvFxbz4GXjmXhYNvPX7OJnc1Zy3aR80hJjWfr1Pt76Yiep8THsrW4gSpyWxXnjBpKWEMPPXihs\nTSzeA/XGmN7Duqp6WOJo+UC//Bt5/PvTzcyYNJw9lQ18uGEPG3dXASBAbLTQ0Hzg+xkTJSTGRlNZ\n38TkkTncPW0s2SlOiRCfs6psDMOYXsm6qnoQ71aA59XZvz3nSHaU1fLRhj08sWgTa4orOemwbKZP\nGExuRiK56Yls3F3FzP+sYOaJTkmQ9bsqWxOHFSE0xhwKm44b5ryn63qXJs/NSGRwVhK7KuqZOeUw\nviyuICs5jvF5mXxV4iSNWZfaehjGmK5jXVURzrtF4nm/cFu5dUUZYwJiYxy9KHHYOIUxpitY4uhF\nicMYY7qClRwxxhgTNJY4jDHGdIglDmOMMR1iicMYY0yHWOIwxhjTIT1yVpWI7AE2hzoOH7KBnnT1\nnb2e8GavJ7yF2+sZoqo5gRzYIxNHuBKRZYFOd4sE9nrCm72e8BbJr8e6qowxxnSIJQ5jjDEdYomj\ne80OdQBdzF5PeLPXE94i9vXYGIcxxpgOsRaHMcaYDrHEEQQicqaIrBORjSJyq4/9p4hIuYisdP/8\nNhRxBkJEHhOR3SKy2s9+EZH73K+1UETGd3eMHRHA64mY9wZARAaLyPsi8qWIfCEiN/k4JmLeowBf\nT8S8RyKSICJLRGSV+/X83scxEfP+tFJV++nCHyAaKAKGA3HAKuBIr2NOAV4PdawBvp5JwHhgtZ/9\nZwFv4axeezywONQxd/L1RMx74453ADDefTsVWO/j31vEvEcBvp6IeY/cf/MU9+1YYDFwfKS+Py0/\n1uLoehOAjar6lao2AP8BpoY4pkOmqh8C+9o4ZCrwlDo+BTJEZED3RNdxAbyeiKKqxar6mft2JbAG\nGOh1WMS8RwG+nojh/ptXue/Gun+8B5Yj5v1pYYmj6w0Etnrc34bvf/gT3c3St0TkqO4JLSgCfb2R\nJCLfGxEZCozD+VbrKSLfozZeD0TQeyQi0SKyEtgNzFPViH9/YkIdQC/1GZCnqlUichbwX2BEiGMy\njoh8b0QkBXgJ+LGqVoQ6ns5q5/VE1Hukqs3AMSKSAbwiIqNV1ecYW6SwFkfX2w4M9rg/yL2tlapW\ntDRfVfVNIFZEsolM7b7eSBKJ742IxOJ8yD6jqi/7OCSi3qP2Xk8kvkcAqloGvA+c6bUrot4fsMQR\nDEuBESIyTETigOnAXM8DRKS/iIj79gSc92Fvt0faNeYCV7pnhhwPlKtqcaiDOlSR9t64Y30UWKOq\n9/g5LGLeo0BeTyS9RyKS425pICKJwBnAWq/DIub9aWFdVV1MVZtE5AbgHZwZVo+p6hcicr17/0PA\nhcAPRKQJqAWmq3t6RbgRkedwZrFki8g24Hc4A3wtr+VNnFkhG4Ea4HuhiTQwAbyeiHlv3E4ErgA+\nd/ejA/wKyIOIfI8CeT2R9B4NAJ4UkWicBDdHVV/3+jyIpPcHsCvHjTHGdJB1VRljjOkQSxzGGGM6\nxBKHMcaYDrHEYYwxpkMscRhjjOkQSxwmZESk2V3ddJWIfCYiE7voeU8RkdcD3d4F5ztPRI70uL9Q\nRNpcS1pEckXkxa6OpTNE5GoRyQ11HCb8WeIwoVSrqseo6ljgl8CfQx3QIToPOLLdozyo6g5VvTBI\n8RyqqwGficN9HYIxgCUOEz7SgFJoXZ/gryKyWkQ+F5GL3dtPcX+bf1FE1orIMx5XEJ/p3vYZcEF7\nJxORZHHW5lgiIitEZKp7+9Ui8rKIvC0iG0TkLx6PuUZE1rsf87CIzHK3ks4F/upuPeW7D7/Ifdx6\nETnZx/mHintNkLbO6fWY34rIUvffZbbHa18oInd5n09EkkRkjjhrW7wiIotFpECcontPePx9fyIi\nFwIFwDPu15EoIpvcz/uZ+/UcIyKfilNc8BURyfQ4/70iskxE1ojIce7Xs0FE/l/7b72JOKGu624/\nvfcHaAZW4pRgKAeOdW//LjAP58r7fsAWnCtwT3EfNwjnS8//gJOABJzqoiNw1jSYg4/1GvBYxwH4\nE3C5+3YGzroPyTjfur8C0t3PuxmnjlAusAnIwrnS/CNglvvxTwAXepxnIXC3+/ZZwHwfsQzFvSaI\nv3P6eEyWx+1/A+e0dT7gZuBf7tujgSac5HAsTpXWlufK8HieAo/tm4Cfe9wvBCa7b/8B+LvH4+5y\n374J2OF+v+JxKr32CfW/Nfvp2h9rcZhQaumqGoVT+O0p97fok4DnVLVZVXcBHwDHuR+zRFW3qaoL\nJ+kMBUYBX6vqBnU+vZ4O4NzfBG51l7VYiPOBnefet0BVy1W1DvgSGIKzzsoHqrpPVRuBF9p5/pbi\nfMvdMbbH1zm9nepuNXwOTAE8y4n7Ot9JOOvBoE411kL39q+A4SJyv4icCbRVTfd5ABFJx0kwH7i3\nP4mzKFaLlnpsnwNfqLOuRr37XJ4F/EwPYInDhAVV/R+QDeS0c2i9x+1mDr3emgDfdSeuY1Q1T1XX\ndOE5Wp4j0Me3eU4RSQAewGnZHA08jJPsOnw+VS0FxuIkzOuBR9o4vDqA2D3P7+LA1+JqLx4TeSxx\nmLAgIqNwuqb24nQDXezui8/B+Wa7pI2HrwWGeowvXBLAKd8BbvQYJxjXzvFLgckikikiMTjdaS0q\ncZY5DaaWJFEizloVgQysfwJMA3DP+jrafTsbiFLVl4Df4CylC228DlUtB0o9xmuuwGkJml7IvgmY\nUEqU/RVQBbhKVZtF5BXgBJz12hWnn32nO7kcRFXrRGQG8IaI1OAknvY+yP8I/B0oFJEo4GvgO/4O\nVtXtIvInnAS2j/3jMuB0Bz0sIjMJ7AO9w1S1TEQeBlYDO3ESWXsewKnM+qU73i9wYh4IPO5+3eDM\naANnrOYhEanF+ft7u8q9PwmnCyrsq7ia4LDquMYESERS1Fl1LgZ4Badk/iuhjssf9xTaWHdizQfm\nA4erakOIQzMRzlocxgTudhE5Hafb6F2cJUvDWRLwvjgr6gnwQ0sapitYi8MYY0yH2OC4McaYDrHE\nYYwxpkMscRhjjOkQSxzGGGM6xBKHMcaYDrHEYYwxpkP+P3NMWEbESQ4VAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7efb3436cb10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from fermilib.utils import MolecularData\n",
    "from fermilibpluginpyscf import run_pyscf\n",
    "\n",
    "# Set molecule parameters.\n",
    "basis = 'sto-3g'\n",
    "multiplicity = 1\n",
    "n_points = 40\n",
    "bond_length_interval = 3.0 / n_points\n",
    "\n",
    "# Set calculation parameters.\n",
    "run_scf = 1\n",
    "run_mp2 = 1\n",
    "run_cisd = 0\n",
    "run_ccsd = 0\n",
    "run_fci = 1\n",
    "delete_input = True\n",
    "delete_output = True\n",
    "\n",
    "# Generate molecule at different bond lengths.\n",
    "hf_energies = []\n",
    "fci_energies = []\n",
    "bond_lengths = []\n",
    "for point in range(1, n_points + 1):\n",
    "    bond_length = bond_length_interval * float(point) + 0.2\n",
    "    bond_lengths += [bond_length]\n",
    "    geometry = [('H', (0., 0., 0.)), ('H', (0., 0., bond_length))]\n",
    "    molecule = MolecularData(\n",
    "        geometry, basis, multiplicity,\n",
    "        description=str(round(bond_length, 2)))\n",
    "    \n",
    "    # Run pyscf.\n",
    "    molecule = run_pyscf(molecule,\n",
    "                         run_scf=run_scf,\n",
    "                         run_mp2=run_mp2,\n",
    "                         run_cisd=run_cisd,\n",
    "                         run_ccsd=run_ccsd,\n",
    "                         run_fci=run_fci)\n",
    "\n",
    "    # Print out some results of calculation.\n",
    "    print('\\nAt bond length of {} angstrom, molecular hydrogen has:'.format(\n",
    "        bond_length))\n",
    "    print('Hartree-Fock energy of {} Hartree.'.format(molecule.hf_energy))\n",
    "    print('MP2 energy of {} Hartree.'.format(molecule.mp2_energy))\n",
    "    print('FCI energy of {} Hartree.'.format(molecule.fci_energy))\n",
    "    print('Nuclear repulsion energy between protons is {} Hartree.'.format(\n",
    "        molecule.nuclear_repulsion))\n",
    "    for orbital in range(molecule.n_orbitals):\n",
    "        print('Spatial orbital {} has energy of {} Hartree.'.format(\n",
    "            orbital, molecule.orbital_energies[orbital]))\n",
    "    hf_energies += [molecule.hf_energy]\n",
    "    fci_energies += [molecule.fci_energy]\n",
    "\n",
    "# Plot.\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "\n",
    "plt.figure(0)\n",
    "plt.plot(bond_lengths, fci_energies, 'x-')\n",
    "plt.plot(bond_lengths, hf_energies, 'o-')\n",
    "plt.ylabel('Energy in Hartree')\n",
    "plt.xlabel('Bond length in angstrom')\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}