{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "### Corrected Bond Stretch" ] }, { "cell_type": "code", "execution_count": 65, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "import matplotlib as mpl\n", "mpl.rcParams[\"font.size\"] = 10\n", "mpl.rcParams[\"font.family\"] = \"sans-serif\"\n", "mpl.rcParams[\"axes.edgecolor\"] = \"#eae8e9\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "### H2PLUS" ] }, { "cell_type": "code", "execution_count": 84, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Binding Energy (hartree)')" ] }, "execution_count": 84, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjYAAAFzCAYAAAAty0N3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAABuLUlEQVR4nO3dd3hb5fXA8e/xjncSJ3H2IBACAZIAgTDDLKvsFspMobRQoIxfCwVaCG0pLW3DKFAoUELDLHuPQhkCQsJISEISyN5O7Dje2zq/P+6VLcmSLctDtnw+z6NHuu99771HXjp+7ztEVTHGGGOMiQcJsQ7AGGOMMaazWGJjjDHGmLhhiY0xxhhj4oYlNsYYY4yJG5bYGGOMMSZuWGJjjDHGmLhhiY0xxhhj4oYlNsYYY4yJG5bYGGOMMSZuWGJjjDHGmLjRKxMbEfm5iKwVkRoR+VJEDm2j/uFuvRoRWSMilwbtnykiGuKR1pHrGmOMMaZ79brERkTOAu4CbgOmAB7gTREZFab+WOANt94U4I/APSJyRlDVMmCo/0NVa6K9rjHGGGO6n/S2RTBFZD7wlape5le2HHhJVW8IUf/PwMmqOtGv7AFgH1Wd7m7PBO5S1dzOuq4xxhhjul9SrANoDxFJAfYF/hS06x3goDCHTXf3+3sbuFhEklW13i3LFJH1QCKwCPitqi7swHURkVQg1bc9cuRIHnt0zoBJkyYVhzvGGGOMMWFlAVsGDRkctlWmVyU2QB5O4rEtqHwbkB/mmPww9ZPc820FVgAzgSVANnAV8ImI7KOqK6O8LsANwC2+jfr6eiZNmtRKdWOMMca0YQSwOdzO3pbY+ARnahKirK36TeWq+hnwWdNOkU+Ar4ArgV904Lq3A7N9Gz+5+OIsYFP/gQNISOic7k2fz/4B+9c2hc7SMTOZ9KPfd8q5jTHGmJ7C6/Wyc0cxQHlr9XpbYlMENNKylWQwLVtTfArC1G8AdoQ6QFW9IvI5sGsHrouq1gK1vu3CbdsBSEhI6LTEplZSSaivaL5mxbZOO7cxxhjT2/SqT0BVrQO+BI4J2nUM8GmYw+aFqH8s8IVf/5oAIiLAZJzbVNFet1s0pOYEbCfUlMQmEGOMMaYH6G0tNuDc2pkrIl/gJC0/BUYBDwCIyO3AcFW9wK3/AHCFiMwGHsLpTHwx8CPfCUXkFpxbUStx+tj8AiexuTzS68aKpvUP2E6uK41RJMYYY0zs9brERlWfEZGBwM04880sBU5Q1fVulaE4CYev/loROQG4EydR2QL8QlWf9zttLvBPnFtNpcBC4DBVXdCO68aEpAcmNin1ltgYY4zpu3pdYgOgqvcD94fZNzNE2YfA1FbOdw1wTUeuGyuJ6QMCtvs1lsUoEmOMMSb2elUfG9NSStbAgO2MxlY7ixtjjDFxzRKbXi41Ky9gO5NK8DbGKBpjjDEmtiyx6eXScwITmwQUrbF+NsYYY/omS2x6uczcvBZlVaVFMYjEGGOMiT1LbHq5nOxsajU5oKyipDBG0RhjjDGxZYlNL5eZlkwJGQFl1dZiY4wxpo+yxKaXExHKJSugrLY85EoRxhhjTNyzxCYOVCUGJjb1FZbYGGOM6ZsssYkD1YnZAduNlcUxisQYY4yJLUts4kBdSuBCmFTvjE0gxhhjTIxZYhMHGlJyA7ZthW9jjDF9lSU2ccCblhuwnVRXEpM4jDHGmFizxCYOSL+gFb7rbOZhY4wxfZMlNnEgISNwIcy0Blvh2xhjTN9kiU0cSM4aELCd4bXExhhjTN9kiU0cSM0MbLHJ0grwemMUjTHGGBM7ltjEgX45gwK2E/FCXXmMojHGGGNixxKbOJCZO6hFmS2rYIwxpi+yxCYO5OTkUq+JAWXlO22Fb2OMMX2PJTZxIKtfCqVBK3xXlVpiY4wxpu+xxCYOJCYIZcErfJfZrShjjDF9jyU2caIyITCxqbMVvo0xxvRBltjECVvh2xhjjLHEJm4Er/CttsK3McaYPsgSmzhRH5TYJNRYYmOMMabvscQmTnhTcwO2E2ttIUxjjDF9jyU28cJW+DbGGGN6Z2IjIj8XkbUiUiMiX4rIoW3UP9ytVyMia0Tk0qD9l4iIR0R2uo93RWRaUJ1ZIqJBj4KueH/RkPTAxCatwRIbY4wxfU+vS2xE5CzgLuA2YArgAd4UkVFh6o8F3nDrTQH+CNwjImf4VZsBPAUcAUwHNgDviMjwoNN9Awz1e+zVKW+qEyQHLYSZ3mhrRRljjOl7el1iA1wLPKKqD6vqclW9GtgIXBam/qXABlW92q3/MPAv4Je+Cqp6rqrer6qLVHUFcAnO1+aooHM1qGqB36PV6X1FJFVEsn2P2XfOzmqtfkekZgUmNplaDqpddTljjDGmR+pViY2IpAD7Au8E7XoHOCjMYdND1H8b2E9EksMckw4kA8GTwewqIlvc22BPi8i4NkK+ASj1Pe69775NbdSPWr/svIDtFBqgvqqrLmeMMcb0SL0qsQHygERgW1D5NiA/zDH5YeonuecL5U/AZuBdv7L5wAXA93BadPKBT0VkYMvDm9wO5PgeV1x++YhW6nZIRm7Lt1JnK3wbY4zpY5JiHUCUgu+xSIiytuqHKkdErgN+BMxQ1ZqmE6i+6VdtiYjMA1YDFwKzQ15UtRao9W0XbtveSogdk5WbR6MKidL8lipKChkwMGTXI2OMMSYu9bYWmyKgkZatM4Np2SrjUxCmfgMQ0KQhIr8EbgSOVdXFrQWiqpXAEmDXiCLvYjnpqS1W+K4usxW+jTHG9C29KrFR1TrgS+CYoF3HAJ+GOWxeiPrHAl+oar2vQER+BfwWOE5Vv2grFhFJBSYCWyOLvmulJCVQRmDf5OpSuxVljDGmb+lViY1rNvATEblIRCaKyJ3AKOABABG5XUT+7Vf/AWC0iMx2618EXAz81VfBvf30B+AiYJ2I5LuPTL86f3XnwxkrIgcAzwHZwGNd/H4jFrzCd72t8G2MMaaP6XV9bFT1GbfD7s04c8ksBU5Q1fVulaE4iY6v/loROQG4E7gc2AL8QlWf9zvtz4EUnGTF363ALPf1CJy5bvKAQuAz4EC/68ZcVWKWc4PN1WCJjTHGmD6m1yU2AKp6P3B/mH0zQ5R9CExt5XxjIrjm2ZFHGBu1yTkBiY2t8G2MMaav6Y23okwYwSt8iyU2xhhj+hhLbOJIY4sVvktiEocxxhgTK5bYxJOgFb6TbYVvY4wxfYwlNnEkIWiF79SGshhFYowxxsSGJTZxJClzQMB2eqMlNsYYY/qWqEdFichuwHCgH86MwN+qqt37iKHkzMD1ojK95TGKxBhjjImNdiU2InIgcClwPM0LSPrWafKKyCLgcWCOJTndr19O4HqcadRBfTUk94tRRMYYY0z3iiixEZHJwF3AYcBynInsvgK2AzXAAGAccCBwG/B7EfkjMNtdBsF0g4ycQS3KGqt2kphjiY0xxpi+IdIWm/nAE8C1qvpVaxVFJAM4G7jOPf8fOhShiVhWbl6LssqSQrJzhsUgGmOMMab7RZrYTFLVlZFUdFe9fkRE5gAjow3MtF9uZj/KNJ1sqWoqqywpInt0DIMyxhhjulFEo6IiTWqCjmlU1XXtjshELS05kVIyA8pqyopiFI0xxhjT/aIe7i0iOSLyPRE5V0T6t32E6Q4VQSt815ZbYmOMMabviCqxEZHf4qyS/Sbwb2CsW/6eiPy688Iz7VWVGJjYNFQWxygSY4wxpvu1O7ERkZ8DtwCPACfiDPf2ec0tMzFSkxS4EKa3yhbCNMYY03dEM0HfFTjDuK8TkcSgfSuBXTselolWfUqOMwDfZSt8G2OM6UuiuRU1Dng7zL5yIDfqaEyHNQSt8J1gK3wbY4zpQ6JJbEqBIWH2jcGZtM/EiAav8F1rE0AbY4zpO6JJbN4DrnMn4vNREUkCLiN8a47pBtIveIVvS2yMMcb0HdH0sbkZ+BxYBryIs07UFcAUYBTww06LzrRbYkbgCt/9GmyFb2OMMX1Hu1tsVHUVcDDOmlE/xxkVdQHOCt+HquqGTo3QtEtKVuBCmJleS2yMMcb0HdG02KCqy4DjRCQVGAjsVNXqTo3MRCUtK3C9qHRqoKEOklJiFJExxhjTfaKeedjveAHqOyEW0wnSc1ouhKk25NsYY0wfEe3Mw0eIyDyc4d3rgb3d8vtE5PROjM+0U2aIFb6rSm1ZBWOMMX1DNDMPHwm8A6QBfw06RxEws1MiM1HJzcqgQtMCyipKCmMUjTHGGNO9ommx+R3whqpOAX4TtO9rYHJHgzLRy0xNshW+jTHG9FnRJDZTgAfd1xq0rxAY3KGITIeICOUSmNjUlu+IUTTGGGNM94omsWkAksPsG4zT78bEUFVidsB2fYUlNsYYY/qGaBKbz4Hzw+w7E5gXfTiREZGfi8haEakRkS9F5NA26h/u1qsRkTUicmmIOmeIyDIRqXWfT+vodWOlOmiFbypslQtjjDF9QzSJzZ+A00TkReBknNtRB4jIvTiJzR2dGF8LInIWcBdwG85tMQ/wpoiMClN/LPCGW28K8EfgHhE5w6/OdOAZYC6wj/v8HxE5INrrtsbb2Ihq8108r9dLY2MjXq83oF5jYyONUdStSA1cyiuxYmub51XVpvKeVjfUe+5pddv6HrWnbrivT0+oaz8nHavb076f9r3vvrr2N6Ljdb1B9cKJZubhd4ELgUOB53HmsbkPOAeYqaoft/ec7XQt8IiqPqyqy1X1amAjzjpVoVwKbFDVq936DwP/An7pV+dq4L+qeruqrlDV23HWxLq6A9dFRFJFJNv3mH3n7CyAT+fNo76+eeqfjRs34vF4WLlyZcDxn3zyCR6Ph9ra2qayLVu24PF4WLFiRUDdzz77DI/HQ1VVFfXpgYlNQtmmgO3PP/8cj8dDeXnzXcPt27fj8XhYunRpQN0vv/wSj8dDaWnzmlM7duzA4/Hw9ddfB9RdtGgRHo+H4uLiprKdO3fi8XhYuHBhQN0lS5bg8XgoLGwesVVWVobH4+GLL74IqLt06VI8Hg/btze3PFVWVuLxeJg/f35A3eXLl+PxeNi6tTmZq66uxuPxMG9eYGPid999h8fjYdOm5q9PXV0dHo+Hjz8O/DFetWoVHo+HDRuaJ9ZubGzE4/Hg8XgCfgHXrl2Lx+Nh7dq1TWWq2lTX/5d4w4YNeDweVq1aFXC9jz/+GI/HQ11dXVPZpk2b8Hg8fPfddwF1582bh8fjobq6eY7MrVu34vF4WL58eUDd+fPn4/F4qKysbCoL973/4osv8Hg8lJU1z15dWFiIx+NhyZIlAXUXLlyIx+Nh587mOZOKi4vxeDwsWrQooO7XX3+Nx+Nhx47mW6SlpaV4PB6+/PLLgLqhvvfl5eV4PB4+//zzgLrLli3D4/FQUFDQVFZVVYXH4+Gzzz4LqLtixQo8Hg9btmxpKqutrcXj8fDJJ58E1F25ciUej4eNGzc2ldXX1zd9P/2tWbMGj8fDunXrmsq8Xm9TXf8PhnXr1uHxeFizZk3AOXx1u/JvhE9BQQEej4dly5YF1LW/EQ77G+HoKX8jPp0X2Q2hdiU2IpIoIrsBrwEjgWOA84DjgJGq+kR7ztdeIpIC7Isz3NzfO8BBYQ6bHqL+28B+IpLcRp2DOnBdgBtwVkMvBUrvve++Ta3U7TSSPSxgO6ve+tgYY4zpG8Q/i2yzsrOCdw3wfVV9s8uiCn/9YcBm4GBV/dSv/EbgQlWdEOKY74A5qvpHv7KDgE+AYaq6VUTqcFqbnvSrcw7wqKqmRnNdd38qkOrb/vX112dde821m/oP6E9iUhIiAjj/0akqIkJCQnOu6cvaExIS2lX33bdf4ZjPLmjaV0MqabdsA/ccoc6rqk3/TSYmJrYaQ3fXDfWee1rdtr5H7akb7uvTE+raz0nX/pyEe889oa597+1vRKx/ThobGthZvBMgZ9CQwWEXQmzXWlGq2iAiBXR8KYaOCs7GJERZW/WDyyM5Z7uuq6q1QFMbceE2p5k0ITGx6ZsFBPyg+PP/RrenbmbeyIB9adRC9U5IHxD2vCISsrwn1A31nuO5LkT/ve/quvZz0rV1wb730dTtqd9P+znp3LoJIcpDiWYRzKdxVvN+PYpjO6oIaATyg8oHA9vCHFMQpn4DsKONOr5zRnPdmMkZ0rI/c33JZpLdxMYYY4yJmirUVUB1CdSUQE1p8+vg5xHT4ICfdmt40SQ2i4CzROR/wAvAVoJaLVT1hY6H1pKq1onIlzh9e17023UM8HKYw+YB3w8qOxb4QlXr/eocA9wZVOfTDlw3ZvIHZFOo2QyS5pa60u0byBu2VwyjMsYY02N4vVBb1nZiUu3u9y+rKQVvQ4TXaegVic2/3efhwIwQ+xWIrL0oOrOBuSLyBU5C8lNgFPAAgIjcDgxXVV8nkweAK0RkNvAQTkfhi4Ef+Z3zbuAjEbkeJ1E5BTgaOCTS6/Yk/dOTWcYABtGc2FQUbqDl8pjGGGPiQl0VVBVBZRFU7XCfi5qfq4oDX9eWgXrbPm9HVZd0/TWCRJPYHEnr/Vm6lKo+IyIDgZuBocBS4ARVXe9WGYqTcPjqrxWRE3BaYy4HtgC/UNXn/ep8KiJnA38Afg+sBs5S1fl+ddq6bo8hIpQk5oF3XVNZXXG3DMgyxhjTUapO4hEySdkRVLbDea6vavu8sVBT0u2XbHdio6ofdEEc7Y3hfuD+MPtmhij7EJjaxjmfA56L9ro9TWXqEGiesgBv2ZbwlY0xxnQ9Vae1pHwLlG0N8VzgtqjsgMa6ts/XIwik5UC/XEjLDXzdLxf6j+32iNqd2IjIGuA0Vf06xL5JwCuqOq4zgjPRq8sITGySgmYfNsYY04ka6qB8q/Mo2xL07Je8NNa2fa7uJoktE5IWzzmh96VmQ5hRU7ESza2oMfjNzRIkDRgddTSm02jmUGcsl6tfja0XZYwxUakpg7LN7iNM8lJZ2PZ5ukNCEqTnQfpAyBjovM7Ic5/d7fQB0K9/cwtLalbTPGfxIJrEBsL3sRmHre7dIyTljgjYzqrrIb90xhjT06g6fVZ2roXiNVDsPvu2q2I4e3tSPzcxGegmK0FJStO2WyctJ66SlGhElNiIyIU460P5/EPEbyyxox/OApIfdlJspgP65QUmNtlaBvXVkNwvRhEZY0wMeb1Oy0rxmsCkpXit86jrzv/JBTIGQfZQyB4OWUOd11lDIWNwYNKSktGNccWHSFts0oFBftu5tLwdVYuzQvYtHQ/LdFTO4JaT9GnZVmSgdX8yxsSpxgYo3RCYsPiSmJ3roKGm62NISnMTlWHNz/6vs4ZCVj4kJrd9LhOViBIbVf0H8A8AEVkLnBGq87DpOQblDaJC08iU5l/kiqINZFliY4zp7WpKYds3ULAUir5tTmBKNoA2tn18tNIHQtaw5taVpqTFr6xf/z5/KyjW2tXHRkT6AZ8BWV0Tjuksg7P7sVH7kynNo6HKtq0nK+RyncYY0wN5vU5ry7alThKzzX2UbOia6yUkQ//RMGCc8+g/1n09FnJGQnJa11zXdKr2LoJZLSLfx229MT1XSlICOxIHsos2JzbVNkmfMaanqimD7cugYInTGrNtKWxbBvWVnXud5HQ3YfFLWnxJTM4ISOjKifM7QX01LHoCvnvb+VpVFUNSKmQOgVEHwj5nw5hD2j7PslfgP+cHlk04EX70ZNfE3Y2iGRX1NTAJ+KiTYzGdrCx5MPjN8dRYYpP0GWNizOuFknVuC4ybwBQsgZJOnMQ9LSdEq4ubxGQO6b23ijZ8Bs/+2JkTx19jrTNT8Y6VsHCuk6Cc9gCkZYc/15JnW5at+i9U73Rup/Vi0SQ2v8ZZM+kbd0Zf00PV9hsSkNgkBP8yGGNMV6oth+3L3VYYXyLzjbMydEclJMOgCTBkTxg4vjlx6T/Wmacl3mxcAI99P3BG4ozBMGwy1FbAps/B667r/O3r8O9T4KK3nNacYDVlsPKdluWNdU5Lzr4XttzXi0ST2NwPZAL/E5GdtFzdW1V1n84IznRMY2Y+lDZvp1Zvi10wxpj4pgo7VsH6T2HDPNg43+nQ2xnS8yB/EgyZBPl7Oc95u0FSSuecv6drqHVaavyTmgN/Dsf8rnl01c718ORZULjc2d7yFbz/Rzjm1pbnW/5q4AixhOTmpGjJs30ysdlBwJy2pqdKyBkOm5u3M2tt9mFjTCdpbICCxU4Ss/5T5zZJVQc/GhKSnIRlyCQ3kdkThuwFWUM6J+beavEzUObXR3LUQXDc7YF1+o+Gs+bC/QeCt8Ep+/xhOPRa59acP//bUP0GwJ6nwRePONvrP3FmV84e2vnvo5tEswjmjC6Iw3SB1AGBk/TlNO4Ab2PP7xxnjOl56qpg85fNicymzzt2Syl9oJPANCUxk5xbS6FunfR1374ZuD3tJ6Hr5e0K42bAqned7boKWPcx7H5ic52K7bDWr4vs7ifApDOaExv1wtLn4aArOi387hbtkgqmF8gaFDhJXyJeZz2TrPwYRWSM6TWqip3bSb5bS1sWNd+uaA9JdD5wmxKYvZyWmKz83tuJt7ttDZo2bsS08HVHTGtObHzH+ic2S18InOtnj1Nh1HSnv06l26q/5Nm+mdiISA6wG85SCgFU1UZM9QADBg+nXhNJluYf4trijaRaYmOMCVa6ybmd5Etkti+L7jxZQ50PytEHwYj9YNBEm/+lo4LXqsps5dZc5uDA7cqg24P+t6FSc2Ds4c7q3BNPgi/+5ZRvXQRFqyBvfNQhx1K7ExsRSQIeAC4Awt3TsHsdPcCQ3Ay2k8twmn8pyrZtYNDo/WMYlTEm5lSh6LvmJGb9PGcpgmgM3BVGT3f6fYyeDrmjrSUmpoLWqPb/XhSvhc1fNG9POL65A/YepzQnNuAkQEfc0HVhdqFoWmyuAb4PXAT8G7gcqAcuAXKAX3RadKZDstOSWMuAgMSmcsfGgEW/jDF9RFWxM6nbd286/S6iWbFaEmHo3k6LjO+RaX9Rulz6QCjzGwlSsQ1yR4auW1HY8lifJc8F7tvjlObXYw516vp+LvpYYnM+cBvwFE5iM19VvwIeFpG3gSOAEAPkTXcTEUqSBkHjyqay+p02+7AxfcaO1U7H02/fcFpm1Nu+45PSYMT+7q2l6c7rVFtRp9sN3Scwsdm0IHxis2lBy2N9lgYlNq9dDa//X/N2rd8K58WrYfNXMHxqVCHHUjSJzTic2Yd9vyH+N08fAO4GemeaF4eq0gaD/4zkZTZJnzFxy+t1bjV8+waseMNZILI90nKbk5hRBzkfin1lrpiebLfjnO+pz+f/ckYyBStaBWs+aN5OzmheXmHrYihcEVi/oo25zZY812cSm0ogRVVVRIqB0cCn7r5qYGDYI023a8jID0hskioLYheMMabz1VXB2g9hxevw3VvOyMdIZQ8PTGQG7e50JDU9yz5nw4d/bm61Wf8xvH0THH0rJLof4yUb4D8XNM9hA7D/xc1z2IRaQqEt37wAx/6h1/1MRJPYrADGuq8/Ba4VEQ/O5P3XAe38F8F0Jc0aBn7z8qXX2CR9xvR6FYVOEvPtm7D6f9BQHdlxyemwy5FOC8C4wyF3VNvHmNhLSoUzHnGWVPANuZ93r5OsDJ0MdZXO0Hz/4fhD94EjbnJeq8I3LzbvS82BX60K3Rr32jXNnYjLt8I6j/Oz0otEk9g8gzPMG+AWnMUwfauX1QOnd0JcppOk9B8esJ3dUOj8kNuoBWN6D1UoWumsAfTtm866QcGjX8LJGAwTjnMWRhx3OCS3mKHD9Aajp8MFL8NzF0GF2/JesQ1Wvt2y7m7HwWkPNg+zX/8plG70239s+FuME09uOToq3hMbVb3f7/VCEdkDOA2nz81/VdVabHqQ9LzADmb9tMZZBTZ4im1jTM/ibXT+C/f1lyleHfmxgyY6Q3l3PxGGTe11txJMGGMOhqsWwcLHnRa7gqVQXQyJKc7cNqOmwz5nwdjDAo8L7jS8+0mtXONQZ3Xv6p3O9vJX4MS/9aoZoUU1wqzfdFjhtu3ZQOnAQXkkdNMfmkVrtjL537sHlHkvnUdC/h7dcn1jTDvUVji3lr590/ngqi6O7DhJcPrI7H6Ck9AMGNe1cRoTA16vlx2FRQA5g4YMLgtXr0NLKojIIELPPBzlTE+msw0ZmEuxZjJAmtd0KS/aQI4lNsb0DN5GWPM+fPVv+PYtaKyN7LiUTBh/FEw4AXY9FtIHdG2cxvQS0cw8nAXcCfyIwKHe/mzm4R5iUGYqK3VAYGKzfSN2I8qYGCvZCIuecG4r+Pd/aE3WUKdFZsKJMPbQXnV7wJjuEk2LzV3AOcAjwGIgwn8vOk5E+gP3ACe7Ra8AV6pqSSvHCE4n558C/YH5wOWq+o27fwBwK3AsMBIoAl4CfquqpX7nWYcztN3fn1X11x19X10pKTGB4sQ88GtEqy2O8I+oMaZzNdQ5M/9+9W9Y9R4RdQAeMslplZlwPAybYh3/jWlDNInNicCvVfXuzg4mAk8CI4Dj3O1/AnNxlngI5zrgWmAm8B3wG+C/IjJBVcuBYe7jl8AynOTlAbfszKBz3Qw85LddQS9QkTIoIP30ltokfcZ0q8LvYOG/YdFTUFXUet2EJBh9cHMy0z/4/yljTGuiSWzSgCWdHUhbRGQiTkJzoKrOd8suAea5SUqL0Vhua83VwG2q+oJbdiGwDafV6UFVXQr4T+G4WkRuAh4XkSRV9ZvtiHJV7XUz3NWm5wckNgkVW2MXjDF9RV0lLHvZaZ3ZMK/t+vl7wdQLYa8znVEpxpioRDM05w3g0M4OJALTgVJfUgOgqp8BpcBBYY4ZC+Tjt3aVqtYCH7ZyDDiLeZYFJTUA14vIDhFZJCI3iUirc42LSKqIZPses++cHZNFVrxZQwO206rbmEbbGBMdVdiy0Jnk7G+7w0uXtZ7UpGbDfhfDTz+ESz+GaZdYUmNMB0XUYuP2Q/H5A/CciJQDrwItlohV1QjHKLZLPgFz6DbZ7u4Ldww4LTT+ttGyvwwAIjIQ+C3wYNCuu4GvgJ3ANOB2nMTpJ63EfANO/x4A7r3vPq695tpWqneNxJzASfqy6mz2YWM6VfVOZ12drx6DgggatEcdBFMvcFZXTknv+viM6UMivRVVRGAvNwH+4j5CiXhUlIjMwu/DP4z93edQPe0kTLm/4P0hjxGRbOB1nL42twacQPVOv83FIrITJ8G7XlVbJHeu24HZvo0rLr88C+j25bXTBowI2M72lkJDrY2oMKYjVGH9J86tpmUvQ0NN6/XT82DyOTDlfBi0W+t1jTFRizSx+R0Rz9/dbvcCT7dRZx2wNzAkxL5BtGyR8fH1h8kH/DuWDA4+xh3G/hZOh+DTVLWe1n3mPo8nRKsVNN32aurdUrgtNi0lWYNDNE6Vb4X+Y7o9FmN6vfJt8PWTTkJTvKaNygLjj3ZaZ3Y7zlbKNqYbRJTYqOqsrgpAVYtwWoRaJSLzgBwRmaaqC9yyA3D6w3wa5rC1OMnNMcBC95gU4HDger9zZwNv4yQhJ6tqG/96ATDFfe7xPXEHDRpMtabQT+qaymqLN5NqiY0xkWlsgNXvuZPovQna2Hr9nJFOy8zkcyB3ZOt1jTGdqkMzD3cnVV0uIm8BD4nIz9zifwKv+Y+IEpEVwA2q+qKqqojcBdwoIiuBlcCNQBXO0HFfS807QDpwHpDtJjoAharaKCLTgQOB93E6K++PM0nhK71hluX8nH5s1QGMk+YBXaXb1zN4l4NjGJUxvUBdFXz5KHx6L5S3MU1CQrKzNtPUC2DcDEiweUqNiYVIOw9fB/xdVasjPbGI7Avkq+rr0QYXwrk4E/T5Rjm9AlwRVGcCBEysewfOsg/30zxB37HuHDYA+wIHuK9XBZ1rLM5tsFrgLJy+QKk4q5k/5J67x8tITaJIBjKO5sSmqqjH52PGxE5tBXz+MMy7FyoLW6+bN8FJZvY5GzLyuic+Y0xYkbbY/Bi4RkQeAR5X1RWhKolIGs4EfhcBR7vHdRp3tNV5bdSRoG0FZrmPUPU/wOlM3No5v8Jpsem1ylLywK/XUEPJ5tgFY0xPVVMGC/4J8+5rfQHK5HTY83QnoRk5zWYDNqYHiTSx2Qu4HGd23htEZDvO0OftQA0wANjFrZeEM7Joqm/ZAhN71Wn5AYkNZT2+a5Ax3ad6J8x/ED67H2pKw9cbNtVJZiadAWnZ4esZY2Im0s7DDcDdInIvcApwAs6EeQfh3OYpAlYAvweeVNW2hgqYbtaQkQ/lzdspVb1uAmVjOl9VsdM6s+CfUFsWvt74Y+Dw65zWGWNMj9auzsOq2gi84D5MLyLZQ/HrYkNGrU3SZ/qwikKY93f4/BGoa2XJtwknwGG/guFTuy82Y0yH9JpRUaZjUgYEDjnNbtgBXi8kRLOqhjG9VHkBfOomNA2tjIWYeLKT0Azdu/tiM8Z0Ckts+ojMvMDEJpkGZ5XhzMExisiYblS6GT65G76cA421YSoJTDodDv0lDNmjO6MzxnQiS2z6iP6DR9CoQqI0TyDdWLqZREtsTDwr2QAf3wkLH4fGutB1JAH2+oGT0NhSB8b0epbY9BFD+mewnf4MpXkIa9n2DfQfPqWVo4zppYrXwsezYdGT4G0IXUcSYZ8fwaHXwsBdujc+Y0yXscSmj8jLSGWx9meoNCc2lYUb6B/DmIzpdEWrwPM3WPxM+GUPEpJhyrlwyDW2XpoxccgSmz4iIUEoScoD7+qmstqd3b7QeI82Y8YM1q1bx7p162Idimmv7SvA81dY+jyoN3SdxBSYeiEcfJWt32RMHGt3YiMirwL3qurbXRCP6UIVqUPAbyCIlrax9o0xPd22ZfDRHfDNS4CGrpOUBvv+GA7+BWQP687ojDExEM1Y34nAGyLynYhc5bdgpOnh6tOHBGwnVfasSfpmzZqFiAQ8ZsyY0eoxM2fObHHMzJkzuyVeE0NlW+Hly+EfB8E3LxIyqUlOh4OuhKsWw/F/sqTGmD6i3S02qjpeRE7AWXxyNvAHEXkCuE9Vl3R2gKbzaNZQ2NG8nVa9LXbBGBONukr45B749B6orwpdJyUTpl0C06+wRSmN6YOi6mOjqm/gtNrsgrOG1EzgEhHxAPcCL7qzFJseJCl3RMB2Tv12ULUF/EzP5/XC10/Ce7+HijAtjak5cMDP4MDLIH1A98ZnjOkxOtR5WFVXA9eKyO+B54AjgEOBLSJyB05fnDA3vk13Sxs8LmC7n1ZD2RbIGR6jiGJHWknmQu1bu3YtY8aM6cKITFgbP4c3r4MtX4Xen5zh3HI68DLol9utoRljep4OJTYiMgK4FPgJMAh4E3gGZ6HMu4AJOLesTA+QPWQslZpKhjTPvKqF3yJ9MLG55ZZbWpTNmTOHkpISrr766hb7cnNzuz4oE6i8AN6dBV8/FXq/JMCU8+GImyBrSOg6xpg+J6rERkSOxElYTgJqgMeAe1R1pVvl3yJyFTALS2x6jNF5GazS4ewjzYuvl238hpzxR8YwqtiYNWtWi7IPPviAdevWhdxnulFDLXz2D/joL+EXqBx/NBzze1v6wBjTQjTDvZcDuwFrgeuAf6lqWYiq84GcjoVnOlN+dhqfywj2oTmxqdz8jX2TTM+x6j1441dQvDr0/oG7wnG3w67HdG9cxpheI5oWm804Cc1rbfSf+QoYG1VUpkuICCUZY6Hqo+ayou9iGJExrrKt8PYN7tDtEFKyYMb1MO1nkJTSvbEZY3qVaIZ7Hx1hvTpgfbsjMl2qfsBu4DdKNqsizH/GPcSqVatC9nnxWbBgQfcFYzpfYwN8/jD87w9QVx66zpTz4KhbbCV6Y0xEbEmFPiZ16ETwW0khs6EEKndAxsCYxdSazZs3c/fdd8c6DNMVNn0Jr10NBYtD7x++H5xwBwzft1vDMsb0btH0sfESdu5yFCgBvgD+rKrvRx+a6QqDR02gdkESqeK34nHRt5BxUOyC6iE++OCDWIfQN9SUwXu/c1pqQv0pScuFY26FKRdAQjSToxtj+rJo/mr8DucWUzHOaKg7gLnu9gbgcWAE8F8RsR5+Pcwu+bms1aEBZVWbl8UoGtPnfPc23H8gfP4QIZOayefClV/CvjMtqTHGRCWavxzFQAEwRlUvUtUbVHUmTkfhbTidiycDHuCmTorTdJJRA9JZQ+CaOeUbl8YomrYdfvjhqGrYx4UXXhjrEE0kKovg+Z/Akz+Ess0t9w+aCD9+E06935ZBMMZ0SDSJzS+Av6pqpX+hqlYAfwV+rqoNwAPA1I6HaDpTcmICRWmBg9W8hd/GKBoT91Rh8X/g3v1hybMt9yf1g6NvhUs9MNpuhxpjOi6azsMjgPow+xqAfPf1ViA5mqBM16rtv6vTtuZKL+vZI6NML1W2BV69Cla+E3r/2MPg+/fAAJsVwhjTeaJpsfkWuEpEApIid/sqdz/AUKCwY+GZrpA8ZPeA7Zy6bVAbZqhtHNi8eTN33XUXxx57LKNGjSIlJYX8/HzOOOMM5s+fH+vw4o8qfP2005cmVFKTlgMn3wsXvGJJjTGm00XTYnMz8DywSkRewvnffwhwKjAcOMOtdwwwr+Mhms42YNQeNH4tJIpf582i7+J2WO3f//53/vznP7PLLrtwzDHHMHjwYFauXMlLL73ESy+9xFNPPcUPf/jDWIcZHyq2w6tXw7evh94/8WQ44S+QlR96vzHGdFC7W2xU9WWcNaK2A1cCt7nP24CTVPUVt95PVPWsTowVEekvInNFpNR9zBWR3DaOERGZJSJbRKRaRD4QkT2D6nwgIhr0eLqj1+6pxg0dyAYNnOysrmBFjKLpetOmTeOjjz5i1apVPPLII9x+++0899xzvP/++yQmJnLZZZdRW1vb9olM6755Ee47IHRSkzkEfjgXzpprSY0xpku1K7ERkRQRORZYpqrTgCxgJJClqgeo6ttdEaSfJ3FGXB3nPibjDDVvzXXAtTiLce6PM6LrvyKSFVTvIZzbZ77Hzzrh2j3SuEHOYpj+SjcsiVE0Xe/000/n0EMPbVF+6KGHcsQRR1BcXMySJfH7/rtcVTE8dxE8OxOqi1vu3+uH8PPPYI+Tuz00Y0zf095bUQ3Aa8DxwAZVrSJggv6uIyITcRKKA1V1vlt2CTBPRCaoaouhPSIiwNXAbar6glt2IU7r0jnAg37Vq1S1oLOu3ZOlpyRRkDIaGr9qKmvY1qveQqdJTnb6tycl2STcUVn9Prx0GZRvbbkvPQ9OutMSGmNMt2pXi42qenEm5M/umnBaNR0o9SUWbjyfAaVAuHGiY3FGab3jd0wt8GGIY84VkSIR+UZE/hrUohPNtRGRVBHJ9j1m3zk7uJUoZqpzxgdsp5WsilEksbNhwwbeffdd8vPz2WuvvWIdTu/SUAtv3wRzTw2d1Ez8vrXSGGNiIppRUY8Al4tIYmcH04Z8nH49wbbTPMQ81DEQMLi5adv/mCeAHwEzgN/jdIB+oYPXBrgBJ/kpBUrvve++Ta3U7VYJgycEbOfUbIL6mhhF45g1a1aLCfjaWuZgzpw5LY6ZM2dOm9eqr6/n/PPPp7a2ljvuuIPExO7+ce7Fti+Hh46Cefe23JeWC6c/7PSnyRzU7aEZY0w07e91wARguYi8gjNfjf/c6Kqqd0Z6MhGZBdzSRrX9fecOdYow5f6C9wcco6oP+e1bKiIrgS9EZKqq+u7XRHPt24HZvo0rLr88i4AlKGMne+Qk8OsvnIAXilfDkD3DHxQnvF4vF110ER999BGXXHIJ559/fqxD6h1UYcFD8N/fQkOIJHiXI+GU+yF7aMt9xhjTTaJJbP7s9/raEPsViDixAe4Fnm6jzjpgb5xh5cEG0bJFxsfXZyYfJwHzGdzKMQBf4UxCuKv7uiCKa/tuezUNtyncFqrRJzZGDx3MFh3AMGnu7Nm4bQWJcZ7YqCqXXHIJjz/+OOeddx4PPPBArEPqHSqLnL40oealSUx1Fq2c9jNb38kYE3PRJDadOqOWqhYBRW3VE5F5QI6ITFPVBW7ZAUAO8GmYw9biJCXHAAvdY1KAw4HrW7ncnjizJvuSoWiu3aONH5zJN97hDEtsTmzKNy4ld+8zWjmqd/N6vfzkJz/h0Ucf5Uc/+hFz5swhwT6I27bW46zzVBGib/3gPeCMh/tES58xpndod2Kjquu7IpAIrrtcRN4CHhIR31DsfwKv+Y9KEpEVwA2q+qKqqojcBdzo3l5aCdyIM5LrSbf+LsC5wBs4CdYewN9wEqFP2nPt3mRgZiobE0cCzcOcawqWxy6gLuaf1Jx11lnMnTvX+tW0xdsIH/0FPvwzqLfl/gMug6NnQXJat4dmjDHhRD3GVUR2x2n5yAMeUdUCERkG7FTV6s4KMMi5wD00j3J6BWd+Gn8TcFpSfO4A+gH3A/2B+cCxqupbQ6AOOApnOYhMYCPwOnCrqja289q9SkX2LlDWvJ1cvDJ2wXQhr9fLxRdfzJw5c/jBD37A448/bklNW8q2wguXwDpPy30Zg+HUf8CuR3d/XMYY04Z2JzbuaKh/AjNp7jz7Js4tnwdxWjpu7rwQm6lqMXBeG3UkaFuBWe4jVP2NOAlah6/d6+TtFpDY5FSth8YGSIyvOV1+97vfMWfOHDIzM9ltt934wx/+0KLOqaeeyuTJk7s/uJ5o1bvwws+gKsQd4nFHwOn/hMzBLfcZY0wPEM0n2E04k9v9CngLWOq3702chKdLEhvTudKH7wlrmreTtB5K1sPAXWIXVBdYt24dABUVFdx2220h64wZM8YSG28jvH8beP7Wcp8kwpG/gYOvtg7CxpgeLZrEZibwe1WdHWIum7V0cudi03VGjBjJDs1ioDSv7K2FK5A4S2zmzJkT0dw2fVpFITx/Eaz9qOW+7BFw5iMw6sDuj8sYY9opmn+9hhN+1e4anPWjTC8wflBmizWjKjd9E6NoTMxsmA8PHhY6qdnteLjUY0mNMabXiCax2Q6MC7NvAj1kAjrTtuG5/VjLiICyqi3LYhSN6Xaq8NkDMOcEKN8SuC8hCb73R/jRU5A+IDbxGWNMFKK5FfUGcJM7/Nk3sYWKSA7wC+DVzgrOdK2EBKE0c2zAMqYJRd/FLiDTfWor4JUr4ZsXWu7LGgY/fAxGTuv+uIwxpoOiSWxuxlndexnwPs6oqD8Ck3Bm6/19p0VnulzDgAkBiU1WxRrnP3mR8AeZ3q1oJTx9LhSFmIJp7GFwxr9snSdjTK/V7ltRqroNZ+2mp4B9gUZgH5wRUQe5w6JNL9Fv2MSA7VRvNZRtjlE0pst99zY8dGTopObQ/4PzX7KkxhjTq0U1YYmb3FzaybGYGMgfMY7y+f3IEr85FQu/hZwR4Q8yvY/X6wzjfv82WqzbmpoDpz8IE46PSWjGGNOZbEKKPm78kCxW67CAspqt1oE4rtSWw7MXwPt/oEVSM2Qv+NkHltQYY+JGVC02InIIziR9o3GWK/CnqnpURwMz3WPMwAxe1eFMZnVTWcWmb7DVf+LEjtVOf5rCEOuATToDTr4XUtK7Py5jjOki7W6xEZEfAx8BP8RZe0mCHtYK1IukJCWwo9+YgDIttJFRcWHVe/DQES2TGkmAY34HZzxiSY0xJu5E02JzHfAf4EJVre3keEwM1PbfFbY1b2eUrYpdMKbjVGH+A/D2jS1X5U7LhTP/BeOtUdUYE5+iSWxGA1daUhM/kvMnBiQ26Q2lUFkEGXmxC8pEp6EO3vglfPVYy32D94SzH4cB4ebXNMaY3i+a20bLgSGdHYiJnbyRu1GryYGFhSGGA5ueraoY5p4WOqnZ4xS4+B1LaowxcS+axOZG4NciMrzNmqZXGD8khzU6NKCsfuuSGEVjorJ9hdOfZv3HLffNuAF+8BikZnZ/XMYY082iuRV1OZADfCcii4AdQftVVU/paGCm++wyKIN3dBQT2dBUVr3qY5Kn21RFvcLK/8JzF0FtWWB5Uj847R+w52mxicsYY2IgmsRmb5zZhrcDw9yHP21xhOnRstKSWZGyFzQ2/7efsmmeLa3QG8z/J7x1fctOwllDnQUsh02JTVzGGBMj7U5sVHVMF8RhYqw4bxps+0fTdlptEexYBXm7xjAqE5a30Rn1NP+BlvuGTYWzn4TsoS33mZhRVbwKXlW8qqg6/zv4tr3q1ElNSiQtOQGxfyqMiUpUE/SZ+DN83B4UFPQnX3Y2F6772BKbnqi2Ap7/CXz3Zst9e54Op94PycHzZppwquoa2FFRR2FFLUXlteyorKOovJaiilqK3Nc7Kuuoqm1wkg+ak5CmRMWrfkkKTYmLf9LSHiLQLzmR9JRE+qUkkp6c5Dyn+MqSSE9OJD3VV5YUWD8lye914PH9khNJSLCkycSviBIbEbkAeF1Vd/iVDQO2qWqjX9lw4GeqenOnR2q61LRxecz/eCKnJH7aVNaw9mOS9vtxDKMyLZRthSd/CAWLW+47/Ncw49d9/vahqlJaXe8kJhV1FFXUssN99i/zlVfVNbZ90m6mClV1jV0WW3pKIrn9kslJTyG3XzL9M5LJ6ZdCbnqys52eQo7vdUaKWzeZ1KTELonHmM4UaYvNo8B03I7CIpIIbMRZ5fsrv3ojgJsAS2x6mSmjcrldJ3IKzYlN4xoPSdbPpucoWAJPntVy9fWEZDjlPtjnrNjE1c2q6hpYW1TJmkLnsaG4yi9pcZKVhvY2kfQxvqRpS2lNu47rl5xI//TmhCg3PZnc9OaEqGm7n/PcPyOZAekpJCXahPSm+0Sa2IT6ZLNPuziSkZpE8aBpsPORprLU6m1QvAYG7hLDyAzgjHx6dibUVQSW9+sPZz0BYw6OSVhdqb7Ry5rCSlYUlLGioJwVW8v4tqC83R/GpvNU1zdSXdr+hKh/ejIDMlIYmJnKwIwUBmamMDAjNejZ2Z/bL9lulZkOsT42psmI8XtTuCCHQVLaXLj+E0tsYu3LOfDataBBtyUGjINznoW88TEJqzPVNXj5tqCcRZtK+HpjCd9sKWPV9nLqG7uv5SUlKYFBmankZaaQl+l82DrPTll2mvOBmyCQIIK4z84DJKAsfJ0EvzIRAuoA1DY0NrWoVNc1UlXXQHW9f1lD4P76UGUN7rFddzurPXZW1bOzqp7VhZVt1k0QnCQoI9VNhpqTHv/XAzJSyMtIJbtfknW0NgEssTFNpo0dyPzPJnJS4mdNZY1rPyZx6gUxjKoPU4X3b4OP/tJy36jpTktNxsDuj6sTbC2tZsHaYhZtdBKZpVvKqGvwtn1gO2WlJpGX5SQmAzNSyctqTlYGBSUumanx+QGpqtTUe6lyk5/q+kbKa+opra5nZ2U9JdX1lFbVUVLtJB8lVXWUVtdTUlXPzqo6ymsaujVer+L2g6qLqH5yojhJjt/3Mi/T/3vutBINynKSoWS7LRb3LLExTfYfM4C/egMTm4Y1HhKtn033a6iDV38BXz/Vct9eP3D61CSldn9cUdpcUs38NTv4bM0O5q8tZv2Oqg6dTwSG5fRj3KAMxuZlMCQ7zWltyQr8MEtLts6uIkI/d4RUNGlwQ6OXspoGSqrq2FlVT2l1nZv0NCdEviTIlxCVVNVR1k0JUX2jsq2slm1lkS1fmJuezMAMX/LjJkBBrwe5rUMZqfYR2Ru157s2Q0RGuK8TcCbiO0JExvjV2a2zAjPdLyc9me0D9oXyR5vKUiu3QMl66D8mdoH1NTWl8J8LYM0HLfcdci0cdXOPTzTLaur5ZGURH35XyCeri9hYXB31uYblpLH70Gx2z89iQn4Wuw3JYmxehiUt3SQpMYEBGSkMyEhp13ENjV52VtVTXFnHDnfofHGFO5y+oo7iSqej9w53f3clQk7iFdltsX7JiU23JPP8+wO12E6xTtI9SHsSmz+FKAvRRm4zD/dm+btMZsfCLAZKeXPhuk8ssekuZVvgiR/AtqWB5ZIAJ/4N9rsoNnG1QVVZUVDOB98W8sG32/ly/c6oRiaNzctg8shc9h6Rwx5Ds9k9P5uc9OS2DzQ9TlJiAoOyUhmUlQpktVm/rsHLzqrm4fnFle7ryjqKK+rYUVnrJkROIlTZDX2Hqusb2bSzmk07I0vM+6cnN3WQzssM7Bztaw3y9RHKTovPW589QaSJzRFdGkWERKQ/cA9wslv0CnClqpa0cowAtwA/BfoD84HLVfUbd/8YYG2Yw3+oqs+69dYBo4P2/1lVfx3Ne+mppo3LY8GXu3N84udNZd51H5Mw5dwYRtVHbF8Oj5/Rcjh3cjqc+ShMOC42cYXR0Ohlwbpi3lpawDvfbKOgrH0jZbLTkthvzAAmj8xtSmZy09vXKmDiR0pSAkOy0xiSnRZR/Zr6xqbJFH1Jj/+cRc1zFzmtQ90xA4Cvk/SqCOomJ0pgC5DbKjbA7SA9ICOVARnJznN6inWSboeIEhtV/bCrA4nQkzhz5fj+wv8TmAt8v5VjrgOuBWYC3wG/Af4rIhNUtRxnPp7gued/6h4XPLXrzcBDfttBY297v/3H9ud+78SAxKZhjQf7uOli6z+Fp852bkP5yxgE5zwDw/eNTVxB6hu9fLp6B28t3crb32yjuDKyDp7g9G2YNmYAB44byAHjBrB7fjaJNqzXRCktOZHhuf0Yntv2LNuNXqWkqi5gcsaiCvcWWVASVFRRS20XdGQPVt+oFJTVRPwPQVKCkJvuS3qaH/0znDLfc1N5egopSX3z1liv6RklIhNxEpoDVXW+W3YJMM9NUr4NcYwAVwO3qeoLbtmFwDbgHOBBd+bkgqDjTgOeUdXgxKVcVQuIkIikAk09PH99/fVZ115zbaSHx8TgrDQ2Zu8L1f9uKksp3wglGyF3ZAwji2PLXnGWSGgM6vw4cDyc+xwMGBubuFxer/LZ2h28+NVm3v6mIOK+ECmJCRwwbgCH7zaIQ3bNY7fBWTY/iYmJxARxh4unMqGN22KqSmVdY1NLUGG5k+z49xXaUdHcP2hnVR3aDa1BDV5tSsoilZWaxIBMN9lJD0yG+vtNpuibXTo3TpIh0e74jnQCEbkImK2quUHlJcA1qvpoiGPGAauBqaq60K/8ZaBEVS8Mccy+wBfAwar6qV/5OpwkJQWnledZ4C+qGvZfVhGZhXMbDIDMzEzWrFrNwEF5JCT03B+eXz+7iOuXnkh/8cvrTnsQ9jk7dkHFqwUPwRu/okXXtBHTnJaa9AExCQtg1fYKXly4iZcWbmFzSWR9DEb078cREwYzY8Igpu8ykPSUXvO/kzFR8XWS3lEZeBusebv5dXf1DeqIjJTEptmkfUtr9E9PJrdfc5lvhmlfcpTTL7lbWl+9Xi87CosAcgYNGVwWrl5v+quTD2wPUb7d3RfuGHBaaPxto2V/GZ+LgeX+SY3rbpzlI3YC04DbgbHAT1qJ+XZgtm/jissvzwI2tVK/R9h/XB6fL57AsYlfNpXpuo8RS2w6jyr87/fg+VvLfRNOgDMegZT0bg+rpKqOlxdt4YWvNvH1ptK2DwD2HJbN8ZPy+d6e+YwfnGn9AEyfEthJum3VdY3NiU7QyDBf52j/R3V99yZClXWNVNZVR/zPDDiDNLPTApfU8CU9YwamM/Pg7m11jnliE9yqEcb+7nOo5iUJU+4veH/IY0SkH84tqt+3OIHqnX6bi0VkJ/CciFzvvzho0DG1QFO7YeG2UHlZzzNt7ADmeCcGJDb11s+m8zTWw6tXw6LHW+7bdyac8DdI7L5fTVVl0cYSHv9sA68t3hJR/4LJI3M5flI+x08ayqiB3Z+AGdNb9UtJZERKOiP6R/Z7U13XSHGVMzKsuMrpCF1cWe8+t0yESqrru+XWmD9VKK12Jn0MnqNq0vDsvpfYAPcCT7dRZx2wNzAkxL5BtGyR8fH1h8kHtvqVDw5zzJlAOvDvEPuC+WaxG4+7OGi8GNG/H6vT94H65g/elNJ1zsrS2cH9rE271FU6az6tfKflvhk3wOHXd9scNRW1Dby8aDNPfLaBZVvDtuo2mTAkizP2Hc5Jew9jWAQdNo0xHdcvJZHhKZF1kobmjtL+yc6Oyjp2us/Fbr+g0mpnUsWSynrKa7tuDqH+MRjpGPPERlWLgKK26onIPCBHRKap6gK37AAgBwi+beSzFie5OQZY6B6TAhwOXB+i/sXAK6paGEHoU9znra3W6oVEhP7jplK2Ip1s8cu+138Ce50Zu8B6u6piZ46azV8ElksCnHSn01rTDdYVVfKvT9bywlebqWjjD1peZiqnTh7GaVOHs8fQbLvNZEwP599ROlL1jV531ujmWaWbX/tml65rWoLDty+S22SxmMKh3YmNiLS2cJAXKAG+UtUt0QYViqouF5G3gIdE5Gdu8T+B1/xHRInICuAGVX1RVVVE7gJuFJGVwErgRqAKZ+g4fseNBw4DTgi+tohMBw4E3gdKcW6N3YmTBG3ozPfZU+w3bhALlk3g6MSmPtfoWg9iiU10SjbC46dD0XeB5Un94AePwoTju/TyqsoX63fy0Edr+O/yba02VScnCsfukc+Z+43g0PF5NpuqMXEuOTGhaYmJ9qipb3Rmcq52kp7S6jp3Lp86St3nySP7d1HU4UXTYjOH5v4p/v+++Zd5RWQucImqdmYb17k4E/T52vFfAa4IqjMBpxXH5w6gH3A/zRP0HevOYePvImCz37n91QJn4fQFSgXW48xnc0e0b6SnO2DsAJ7xTgxIbBrWerA5YKOwbZkz8V55UK6flgvn/AdGHdBll25o9PLWNwU85FnL1xtLWq07on8/zjlgFD/Yd2TEHSGNMX1XWnIi+TmJ5OdENqlid4kmsZkGPIOTADyF01clH/gRcCxwGbAvzmR264BbOyNQAFUtBs5ro44EbSswy320dtyNOK05ofZ9hdNi02fsMiiT5al7g7e5YSt552oo3wZZobo6mZDWz4Onzmo58V72cDjvBRi8e5dctrahkf98sYkHP1zd6nTwCQJH7j6E8w4cxWG7DrJ5ZowxvV40ic1VwIuq+ku/sm+BD0Xkb8BFqnqWu/zBuXRiYmO6j4iQPWZfylf3I0v8PhjXfwKTTo9dYL3JijfguR9DQ9DMooN2d5KanOGdfsma+kaeWrCBBz9c0+qMpllpSZwzbRQXHDQm4k6JxpjY2VhcxaF3vN+0fcDYATzzs+mtHnPwn/4XMGxbxLntlJmaRF5mCuMHZ3LI+EGcNmU4/VLaXlS2vtHLtNveZWdVfVNZUoKw4Kaj271IaleKJrE5CfhBmH1vAs+5r/8HXBlNUKZn2G/cIL5cuRszEr9uLrTEJjJfPgavXQ0aNHR65AHwo6c7feK9qroGnvhsAw9+tKbVmUmH5/bjokPGctb+I8lMjfnYAWNMN1J1FhstbnBGR323rYI3lhRwx9sruO3UvThx79ZHvX70XWFAUgPOjMivL97C+dPHdGHk7RPNX7ZEYBfg3RD7xtPc76YOvzlcTO8zbewAXvdODEhsGlb+jyTVbhuS3Ouowkd/gfdva7lvt+OcxSw7ceK96rpG/j1vHQ9+tKbVdZv2GZHDTw4dx/GT8q0zsDF91PGT8mn0Kjsq61i6ubRpzqqSqnouf/IrCsr24OJDws858/Ki0GOCXl7U+xObd4A/iMgaVf2vr1BEvoczsZ2v8+3uOH1sTC81cWg2tyXtjf80Q0kla2Djgi7t8NpreRvhzevg84db7ptyHpx0d6dNvFff6OU/X2zk7ndXsr08/P8P08cN5MqjxjN93EAbqm1MH/eP85oX0y2vqeee91bykGdtU9kf31jO1FG5TBnVciRTdV0j7y5vnv4tOVGob3TGDH25YSebdlZFPOlgV4vmX7ercIZ0vyUiJSLyrbte0xs4yw1c7Vf3rg7GZ2IoMUFIHb0fa71BnYUXRjJ/YR9TX+P0pwmV1Bz6f3DyvZ2S1Hi9yqtfb+GY2R9y04tLwyY1h+6ax7OXTuepnx7IQbvkWVJjjAmQlZbMTScGttA0epW73l0Zsv47ywqo8lvn6tLDd2l6rQqvfN2pM7x0SLsTG1XdDOyD03/mLWCD+3wFMNndj6reF2phStO7TBs3kP80HhFQpktfhNrg0fJ9WE2pM5x72ctBOwSOvwOOurnDt+5UlQ+/K+T7937MlU8tZF3QtOU+R0wYxAs/P4i5Fx/A/mNit4CmMaZ3uPLI8aT6rej98aoiymvqW9Tzvw2VnpLI5UeMZ9SA5haaV8LcpoqFqG62q2qVqt6vqmer6jHu8z9UNfRfW9NrHTB2IM81HkqDNv+oSH0lLH0hhlH1IOUF8OiJsP7jwPLEFDjzX3DAz0If1w4rCso4/5EFXPivBXyzJfTSB4fumscrVxzMoz+extQQzcjGGBNKbnoKew1vnvqt0ass3Rz4d6akqg7PyuYJ+Y/YfTBpyYkcP6l5/ekVBeV8W9Az/uG1XoSmVVNG5pKSO4z3vVMCd3xlt6Mo/A4ePga2LQksT8mCc5/r8OixoopabnhhCSfc7eHjVaFXHdlnZC5P/uQA5l58AHuPyO3Q9YwxfVPwBHs7KgNvcb+2eGtTfxqgKaE5zi+xAXh50eYuirB92p3YiEiyiPxGRJaJSKWINAY9um41LdPtEhKEM6YO5+nGGYE7Nn/hzKjbV234DP51LJQGraiRMRh+/DqMOzzqU9c2NPLAh6uZ8ZcPeGrBBrwhlj8YPziTB87bl5d+fhAHjc+L+lrGGBP8J0YIvHXuf5spLTmBI3cfDMCUUf0D5sHqKf1sounNeDtwDc6cNS9hQ7rj3hn7juC+/01mm+YyREqadyycC8fdHrO4Ymb5q/D8T1pOvNd/LJz/IgwIP1yyNarKW0sL+OOby9lYHHq24PzsNK49ZjdOnzrchm0bYzrFlpLAvzf+k+1tKanm8/XFTduH7zaI9JTm1OF7e+bzr0+ckVWbdlbz5fpi9h0d2/590SQ2PwR+p6o2o3AfMXpgBlPH5PHcpsO4POmV5h1fPwVHz4KkPrSu0IKH4I1f0eJ/nGFTnXWfMgdFddpV2yuY9co3YW859UtO5NLDd+Gnh42LaIZQY4yJxM7KOr7x61OTmCBMGp7dtP3K11sCFs39ZNUODvzje03bwSt8v7xoS8wTm2j+5esPfNTZgZie7cx9R/Cf4NtR1TthxWsxiafbeb3w31vgjV/SIqnZ9ViY+VpUSU1lbQN/enMFx9/9Udik5vSpw3n/lzO46uhdLakxxnSqu99bSV1j8wzph4zPIyutebnj4En5KmobKCiraXqUVgeOoHp98VYaGoNmXO9m0SQ2HwGTOzkO08OdsNdQtiUN49PGPQJ3fDU3NgF1p4ZaeOlS+OSulvumXgBnPwUpGe06pary+uKtHD37Qx74cHVAxzyf/Ub35+XLD2b2Dyf3uNVzjTG9W3lNPX94bRlzPl3XVJaYIFx99K5N26u2l7N8a+iRmOHsqKwL+09ad4nmVtQvgJdFZD3wmqqGn8fdxI2stGSOnzSUZ76ewUGJfp2G17wPO9dD/9Exi61LVRXDM+c5a2QFm3EjHH5du+eoWVNYwc0vh7/tlJ+dxk0nTuSkvYfaxHrGmJBWbq/gsse/DLlvn5G5ARPo+Vz2+Jd4VdlRUccSvyUVfH574sSAWYdfWhjYWnPvOVM4ae9hLc67cMNOTrv/06btVxZtYcaEwe16P50pmsRmEZAMPAuoiATPXaOqmtPiKNPrnbnvCC5aOI1SnUOO/7d90RNwxI2xC6yr7FgNT/wAilcHlksifP8up7WmHWobGvnHB6u5//3VAU2/PsmJwsWHjOPKI8eTYQtUGmNaUVxZx5tLC0Luawg1lBLC1s9NT+aPp+3FCXsFLoL56uLmxCYlKSFssjJ5ZC5Dc9LYWuoMqHhn2TZq6htJS47NrfNo/no+T8vRYaYPmD5uIANzsnmx8hBmJr3TvGPh43D49ZAQR/0/1n8KT5/j9CPyl5wBP5gDux3brtPNW72Dm15cwpqiypD7Dxmfx6yT92T84MwoAzbGmNaJQHJCAumpiQzKTGWXQZkcsmsep08dHjDSCZxWmPV+M5wfMj6PzDD/cIkI39szv+m2VkVtA+8u3xaydac7iKrlKN2lcNv2bKB04KA8EhJ651Ddv73zLe+9/y5vpAa10Jz7HOx6TGyC6mxfPwOvXAGNQXdZs4bCOc/A0H0iPlVxZR23vb6c57/aFHL/0Jw0fnvSHhw/Kd9uOxljTCu8Xi87CosAcgYNGRy284+1d5t2OWPqCP7+vzEs9o5l74TmVWH56t+9P7FRhQ/+BB/+qeW+/L3gR89AzvAIT6U89+Um/vjGcnZWtVx3JTFBuPiQsVx11K5228kYYzpRRH9RReQw4CtVrXBft0pVbTh4nBqTl8H+Y/rzzMYjAhObb9+A8m2QNST8wT1ZXSW8fDl882LLfbsdB2c8AqmR3SZav6OSG15Ywqerd4Tcv8/IXP542iT2HGZd0YwxprNF+q/iB8CBwAL3dbj7V+Lui6POFibYGVNHcNu6g/hN0uP0E/d2jbcB3rkJzng4tsFFo3gNPH0ebP+m5b4DLoXv/TGi/kMNjV4e/ngtd/73uxajDQCyUpO47rgJnHPAaBIT7LaTMcZ0hUgTmyOAZX6vTR92wt5DmfVqJq80HsRZSR8071jyLEw6AyYcH7PY2m3Vu/DcxVBTElguCXDcn+GAn0Z0mqWbS7n++cVhV98+ca+h3Pz9PRiSbfPRGGNMV4oosVHVD0O9Nn1Tdloyx+2Zz18X/YDjEhcEDv1+7RoYNR365cYsvoioOhPuvfc70KDWldQcOPORiPoMVdc1cue73/GwZ03IxSqH5aTx+1MncdTEXnqLzhhjepneOTTHxNyZ+46kkP78rj5oLpfyrfDOb2ITVKRqK+DZC+HdWS2TmkG7w0/fjyip+XhlEd+76yP++VHLpEYEZh40hneuPdySGmOM6UaRdh7+VzvOqap6cZTxmF5i+i4DGZaTxvOlh/L9xnnMSPy6eefCuTDpdNjlyNgFGE7BUnj+Yihc0XLfxJPh1PshNavVU5RU1fGH15fz3Jehh3DvNiSTP52xN1P9ZvA0xhjTPSLtY3MkgR2Gc4EcoAHYAQx0z1UK7Aw+2MSfxAThgoPG8Kc3V3Bj/cW8k3AdmVLTXOGVq+Dn8yIeSdTlvF747H5479aW89MgcNTNcMg1rS6PoKq8tngrt776DUUVLVcSSUlM4Mojx/Ozw3chJckaQ40xJhYi+uurqmNUdayqjgV+CFQA5wL9VHUo0A84D6gEzu6qYE3P8uODxzBuUAZbyONPDT8K3Fm6wUkieoKyLfD4ac6oreCkJi3HmVzw0GtbTWq2lFTzk8e+4MqnFoZMavYf0583rjqUK4/a1ZIaY4yJoXbPPCwiHwHPq+rdIfZdA5ypqgd3UnxxJR5mHg72yaoizn14PoKXp1Ju48CE5YEVfvwmjD4oNsEBLHsFXv1Fy6URwJl074f/hgHjwh7e6FUe+3Qdf3vnWyrrGlvsz0xN4tfH784500aRYEO4jTGmy0Q683A0n677AkvD7FsCTI7inBERkf4iMldESt3HXBHJbeOY00XkbREpEhEVkRbxiUiqiPzdrVMpIq+IyIiOXrsvOHh8HifvMwwlgevrL6FaUwIrvHwF1Fd3f2A1pc61/3N+iKRG4OCr4Cf/azWpWbKplFPv+4TfvbYsZFJz9MQhvHvt4Zx34GhLaowxpoeIJrEpA44Os+9od39XeRIncTrOfUwG5rZxTAbwCfDrVurcBZyGcxvtECATeE1E/Gdli+bafcJvTpxIVmoS6zWfvzb8IHBn8Wpnnpja8u4JxtsIXzwK90x1OjEHyxoGF74Cx/wOklJa7sdZwO3WV7/hlPs+Zsnm0hb78zJTuf/cqTx0wb7k59i8NMYY05NEcyvqDuBa4E6cD/sCIB+nz83VwGxVva5zwwQRmYgzSeCBqjrfLTsQmAfsrqrftnH8GGAtMEVVF/mV5wCFwPmq+oxbNgzYCJygqm9He20RSQVSfdu/vv76rGuvuXZTPN2K8pnzyVpmvbqMBLw8nzKLKQmrAivkTYCzn4C8XbsuiLUeeOsG2LYk9P49ToWT7oT0ASF3qyrvLNvGLS9/Q0FZTcg6Z+03khtO2J3c9NBJkTHGmK7RlbeibgQeB/4P+ALY5D5fi5Po3Bj+0A6ZDpT6EgsAVf0MZyRWRzpx7AskA+/4nXcLzu0233mjvfYNbp1SoPTe++4LPT44Dpx34Gj2HJaNlwR+Vf9TajQ5sELRt/DPI2D5a51/8Z3r4Jnz4bGTQic1KZlw6j/gB3PCJjWrCyv48ZzP+dncL0MmNeMHZ/Kfn03nz2fubUmNMcb0YO1eVlhVG4CZInI7zvIKA3GGfH+gqiEmB+k0+cD2EOXb3X0dOW+dqgZ3xNjmd95or307MNu3ccXll2fhJIJxJykxgT+cOonT//Epq3QEF9f/kvuS7yFXKpsr1ZXDM+fCob+EI26MaP2lVm39Gr6a66ws3lgbus6EE+F7t8GAsSF3l9fUc897K3n0k3U0hJg6ODUpgV8ctSuXHDrORjsZY0wv0O7Exse9/dLq7Z9IiMgs4JY2qu3vu2yoU4Qp76jg87b72qpaCzR94hZuC5UbxY8po/pz9v6jeGrBBj7x7sVJdbfxYPKd7JmwPrCi56+wZSGcch9kD23fRap3wpLnnGSmYHH4eoP3gONuh3EzQu72epXnvtrEHW99S1FF6KTo0F3z+P0pkxiTl9G+GI0xxsRM1ImNiAwGRuPMYRNAVT9qx6nuBZ5uo846YG8g1Nz0g3BaV6JVAKSISP+gVpvBwKd+dbri2nHn+uMm8PY3BRRX1rFJB3NG3Sz+nPIIpyR8HFhx9Xswe3cnARl7OIw7HEYfDGnZzXVUobIQSjdByXrnNtbyV8O3zgD0G+C0Bu37Y0gM/eP95fpifvfacr7eWBJy/+CsVG46cSIn7zMMaWVuG2OMMT1PuxMbERmKMxrIt8q37y+/0tyCEfE9BlUtAooiuO48IEdEpqnqArfsAJwZkD9t9eDWfQnUA8cA/3HPOxSYBPg6QXfVteNObnoKNxy/O796zmlNqSGVq+ouY2HiLtyc8gQJ2hB4wPZlzmP+P0ASYfi+kJQKZZuhdHPrSYw/SYRpP4UZ10O/0EsZLN5Uwuz/fscH3xaG3J+SmMDFh47l8iPGk5kadc5vjDEmhqIZFfU8MAOn/8hi/G61+HTVCuAi8iYwDPiZW/RPYL2qft+vzgrgBlV90d0eAIxyj3sdZ0j3t0CBqha4df4BnATMBIqBv+L0HdpXVRsjvXZb4nGCvlBUlSufWshri7cGlO8vK3gs8z7S63d03sWyh8Pkc2HKedB/dMgqy7aUcee73/HfZeEb146eOJjfnLiH3XYyxpgeKtJRUdH8W3o48EtVfTTa4DrgXOAemkcwvQJcEVRnAk5Lis/JgH+svttetwKz3NfX4Kx79R+cW2vvATN9SU07rm0AEeHus6eQlZbEUws2NpV/rrszo/x3PDLsVSZVfILURTm3TUIy7H4iTD0fxh0RthPyym3l3Pnud7yxpCDsqcYNyuDmk/ZgxoTB0cVijDGmR4mmxaYQ+JGqvts1IcWvvtJi46Oq/PWdb7nv/dUt9p05eQh/OrCepPUeWPMhbFoQYnFKlyQ6LTP9RzsJzV4/hIyBIavWNjTy1tICnl6wkXlrwrcM5WWm8PMZ4znvwNE22skYY3qBSFtsokls7scZHn11hyLsg/paYuPzyMdr+f1ry1qUTxqezQXTx/D9vYfRj1rYMA+2fAWJqZAzHHJGQs4IyBzS5tDwldvKeWrBRl5YuImSqvqw9fqnJ3Pp4btw/vTRpKdYPxpjjOktujKxmQE8BLwJvIozh00AVf2qXSftI/pqYgPw4sJN/OrZxSHnislKS+L0KcP50QGj2D0/O8TRLZXX1LN4UykLN+zk/W8L+XJ9iEUu/WSnJfHTw8Yx8+Cx1jHYGGN6oa5MbLx+m8EHC6Cq2sGZ1+JTX05sAN5fsZ3LnviSmnpv2DpTR+Vy2G6DSE9JJC05kbSkRNJSEklLSqC4so5FG0tYuKGE77aXE8mPbm56MhdMH8PFh4wlp19y2wcYY4zpkboysbmwrTqq+li7TtpH9PXEBuDL9Tu56umFbNrZtSt+Tx83kLOnjeR7e+aTlmx5tjHG9HZdltiY6Fli42ho9PLeiu08tWADH35XGFHLSyTyMlP5wX4jOGu/kTZs2xhj4kxXDvc2pkOSEhP43p75fG/PfDYWV/HM5xt55ouNFJZHOBmfKzlR2GNYDlNG5nLQLgM5YvfBJCf23YTRGGNMhC02InIz8LCqbnFft0ZV9fedEl2csRab8Oobvby3fDvvr9jOzqo6ahq81NQ3UlvfSHV9IzX1XpIShD2GZTNlVH8mj8xlz2HZdpvJGGP6iE69FeV2GD5QVRcEdR4OxToPh2GJjTHGGBOdTr0VpaoJoV4bY4wxxvQklqQYY4wxJm50OLERkWwR2U9EQq9AaIwxxhjTTSJObETkRBF5XEQeFZGj3bKfA1uA+cAaEXlWRGwWNGOMMcbERER9bETkRJzlE2qBOuA8EbkC+DvOateLgSnA6cBlOKtgG2OMMcZ0q0hbbH4JvAf0V9Uc4H7gTuCfqnq8ql6vqscCc4A2ZyY2xhhjjOkKkSY2ewL3qmqNu/0XIA14Maje88DYTorNGGOMMaZdIk1s8oACv23f66KgejuAnI4GZYwxxhgTjfaMitIwr40xxhhjeoT2rBU1QUQa3Ne+mYV3FxH/Ort3SlTGGGOMMVFoT2IzJ0TZ3KBtwVpzjDHGGBMjkSY2P+7SKIwxxhhjOkGka0U91tWBGGOMMcZ0lK0VZYwxxpi4YYmNMcYYY+KGJTbGGGOMiRuW2BhjjDEmblhiY4wxxpi40asSGxHpLyJzRaTUfcwVkdw2jjldRN4WkSIRURGZHLR/gIj8XUS+FZEqEdkgIveISE5QvXXu8f6PP3X+uzTGGGNMtNozQV9P8CQwAjjO3f4nziSB32/lmAzgE+BZ4KEQ+4e5j18Cy4DRwANu2ZlBdW8OOkdF+8I3xhhjTFfqNYmNiEzESWgOVNX5btklwDwRmaCq34Y6TlXnunXHhNm/FDjDr2i1iNwEPC4iSara4LevXFULMMYYY0yP1JtuRU0HSn1JDYCqfgaUAgd18rVygLKgpAbgehHZISKLROQmEUlp7SQikioi2b7H7DtnZ3VynMYYY4zx05sSm3xge4jy7e6+TiEiA4HfAg8G7bobOBs4ArgXuBq4v43T3YCTeJUCpffed9+mzorTGGOMMS3FPLERkVkhOuUGP/Zzq4daYLPTFt4UkWzgdZy+Nrf671PVO1X1Q1VdrKoPA5cCF7uJUDi347T+5AA5V1x++YjOiNMYY4wxofWEPjb3Ak+3UWcdsDcwJMS+QcC2jgYhIlnAWzgdgk9T1fo2DvnMfR4P7AhVQVVrgVrfduG2UA1OxhhjjOksMU9sVLUIKGqrnojMA3JEZJqqLnDLDsBpDfm0IzG4LTVv4yQhJ6tqTQSHTXGft3bk2sYYY4zpPDFPbCKlqstF5C3gIRH5mVv8T+A1/xFRIrICuEFVX3S3BwCjcIZvA0wQEYACVS1wW2reAdKB84BsN9EBKFTVRhGZDhwIvI/TX2Z/4E7gFVXd0HXv2hhjjDHt0WsSG9e5wD04iQjAK8AVQXUm4LTi+JwMPOq37bvtdSswC9gXOMAtWxV0rrE4t8FqgbOAW4BUYD3OfDZ3RPUujDHGGNMlRLVT+t2aCBRu254NlA4clEdCQsz7bRtjjDG9htfrZUdhEUDOoCGDy8LVs09XY4wxxsQNS2yMMcYYEzcssTHGGGNM3LDExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3LDExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3LDExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3LDExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3LDExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3LDExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3OhViY2I9BeRuSJS6j7mikhuG8ecLiJvi0iRiKiITA5R5wN3n//j6Y5e2xhjjDHdq1clNsCTwGTgOPcxGZjbxjEZwCfAr9uo9xAw1O/xs064tjHGGGO6UVKsA4iUiEzESSgOVNX5btklwDwRmaCq34Y6TlXnunXHtHGJKlUt6MxrG2OMMaZ79ZrEBpgOlPoSCwBV/UxESoGDgI4mF+eKyHnANuBN4FZVLe/ItUUkFUj1bf/mppuyfnHlL/B6vR0M1RhjjOlbIv3s7E2JTT6wPUT5dndfRzwBrAUKgEnA7cA+wDEdvPYNwC2+jYcfeYRfXPkLdu4o7mC4xhhjTJ+VBZSF2xnzxEZEZuH34R/G/u6zhjpFmPKIqepDfptLRWQl8IWITFXVrzpw7duB2b6N5ORkli5dOmDSpEmW2QCz75ydde9992264vLLR1x7zbXlbR9hOoN93WPDvu6xYV/32OjCr3sWsKW1CqLaoZygw0QkD8hro9o64BxgtqrmBh1fAlyjqo+2cZ0xOK0yU1R1URt1BagFzlfVZ0Tkoo5c24QmItlAKZCjqmGzb9O57OseG/Z1jw37usdGLL/uMW+xUdUioKiteiIyD8gRkWmqusAtOwDIAT7t5LD2BJKBre52d17bGGOMMVHqNcO9VXU58BbwkIgcKCIH4gzRfs1/VJKIrBCR0/y2B7hz1+zhFk0Qkckiku/u30VEbhaR/URkjIicADwLLMQZJh7xtY0xxhgTW70msXGdCywB3nEfi4Hzg+pMwGlJ8TkZJ0l53d1+2t2+1N2uA44C3sYZ3XSPe+6jVbWxndc27VML3Oo+m+5jX/fYsK97bNjXPTZi9nWPeR8bY4wxxpjO0ttabIwxxhhjwrLExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHdSkRuEJHPRaRcRLaLyEsiMiHWcfU17vdBReSuWMcS70RkuIg8LiI7RKRKRBaJyL6xjiueiUiSiPxBRNaKSLWIrHGn9bDPvE4kIoeJyKsissX9e3Jq0H4RkVnu/moR+UBE9uzquOybbLrb4cB9wIE4a3ElAe+ISEZMo+pDRGR/4Kc4UxaYLiQi/XHmw6oHjseZT+v/gJIYhtUXXI8zpccVwETgOuBXwJWxDCoOZQBf43ydQ7kOuNbdvz/Oeoz/FZGsrgzKhnubmBKRQTiLiR6uqh/FOp54JyKZwFfAz4HfAItU9eqYBhXHRORPwMGqemisY+lLROQ1YJuqXuxX9jxQpao2/1gXEBEFTlPVl9xtwVnT6S5V/bNblgpsA65X1Qe7KhZrsTGx5ptM0RYG7R73Aa+r6ruxDqSPOBlnQd1n3VuvC0XkklgH1Qd8DBwlIrsBiMg+wCHAGzGNqm8ZC+TjTGgLgKrWAh8CB3XlhWO+VpTpu9yMfjbwsaoujXU88U5Ezgam4jQJm+4xDrgM5+f8j8A04B4RqVXVf8c0svj2Z5x/mlaISCOQCNykqk/FNqw+Jd993hZUvg0Y3ZUXtsTGxNK9wN44/0mZLiQiI4G7gWNVtSbW8fQhCcAXqnqju73Q7Tx5GWCJTdc5CzgPOAf4BpgM3CUiW1T1sVgG1gcF93eREGWdyhIbExMi8necZvrDVHVTrOPpA/YFBgNfOg1lgPNf7GEicgWQGrQ2mukcW4FlQWXLgTNiEEtf8hfgT6r6tLu9RERGAzcAlth0jwL3OR/n98BnMC1bcTqV9bEx3cod/ncvcDpwpKqujXVMfcR7wF44/7n6Hl8ATwCTLanpMp/gLMzrbzdgfQxi6UvSAW9QWSP2mded1uIkN8f4CkQkBWdk7KddeWFrsTHd7T6c5uFTgHIR8d2HLVXV6tiFFd9UtRwI6MckIpXADuvf1KXuBD4VkRuB/+D0sfmp+zBd51XgJhHZgHMragrOsON/xTSqOOOOshzvVzRWRCYDxaq6wZ0n60YRWQmsBG4EqoAnuzQuG+5tupM7JDCUH6vqnO6Mpa8TkQ+w4d5dTkROAm4HdsX5L3a2qj4U26jimztPyu+B03BufWwBngJ+p6p1sYwtnojIDOD9ELseU9WZ7gCRW4CfAf2B+cDlXf3PlCU2xhhjjIkbdr/RGGOMMXHDEhtjjDHGxA1LbIwxxhgTNyyxMcYYY0zcsMTGGGOMMXHDEhtjjDHGxA1LbIwxxhgTNyyxMcYYY0zcsMTGmC4gIgeIyIsiskFEakVkm4jME5G/xTq2ziQiJ4jIrDD71onInO6NyJkNVUTU71EnIoUi8omI3OYuhhh8zEy37ph2XutGETm1s2LvCUQkV0SKRORsv7I5IlLRyddZJyKvdfAcF4vIZhHJ6Ky4TO9niY0xnUxETsRZ5C0buA44FrgKZ0HEs2IYWlc4AWfK9FBOw5nWPlZuBKYDRwAXAx8AFwHLReTcoLqvu3W30j43Aqd2KMqe5xacJQieiXUgEXgMqMT5PTMGsEUwjekK1+GsCfQ9VW3wK39aRHr0H2ARSVfVqs44l6ou7IzzdMBKVf3Mb/sVt8XsXWCOiCxW1SUAqloIFMYiyJ5ERAbgrOtzjfbg9XZEpB9Qo6oNIvIg8FsR+XNn/eya3s1abIzpfAOBoqCkBgBV9QaXichZ7m2qShGpEJG3RWRKUJ057r49ReQ9t26hiNwrIulBdS8XkY9EZLtbb4mIXCciyUH1PhCRpSJymIh8KiJVuKsfuzG9IyJbRaRaRJaLyJ/8m/zd20yXu6/9b/2Mccta3IoSkVEi8rgbW6173v8TkQS/OmPc8/xSRK4VkbXue58nIgdG9i0ITVWLcT64k4Br/K7Z4laUiEwRkdf8Yt0iIq+LyAjfewYygAv93vsH7r5BInK/iCxzY98uIv8TkUODvh7teq/i3OJ8VUR2iEiNiKwWZwVl/zq7isiTQV/jyyP8Es10vzYhW2si/PlLE5Hb3fdS594quk9EcsOc8zgR+cr9OVshIhcF7fd9b44VkX+JSCHOCtGpbpUncFpHz8YYLLExpivMAw4QkXvcD6LkcBVF5EacVYeXAT8EzgeyAI+I7BFUPRl4A3gP5/bHvTgf0sEfQrsAT7rnOgl4BPgV8GCIEIYCj7v1TwDud8t3da91MXAccJcb36t+x/4eeM59Pd3vEfJ2jogMwrlFdyzwW+BknNaTv7rvJdjlwDHA1cC5OEnEGyKSE+r8kVLVz90YDwtXx03g/gsMCYpjA873B5z3Wo3zdfK995+7+wa4z7cCJwI/BtYAH4izInKwNt+riHwP8ACjgGuB44E/uDH66uwBfA5MAv4P5/v/OnCPiIS7ZejvRGChqpaE2Nfmz5+ICPAS8Etgrnu+2cCFwP9EJJVA+wB/A+4ETgEWA4+ISKjvzb+Aepyf6zPd16hqAbDCvZYxoKr2sIc9OvGB02LjAdR91OH0r/k1kOlXbyTOH+d7go7PxPngfcavbI57rl8E1b3RLT84TCwJOP+Bnw80AP399n3gHntkG+9H3HMc5tbf22/fvc6fkZDHrQPm+G3f7h4/Laje/YAX2M3dHuPWWwwk+tXb3y0/u414Z7j1zmylzmdAld/2TPeYMe72vu72KW1cq8L/PbZSL9H9Gr4LvOBXHvF7BVa5j7RWrvMWsBHIDir/O04S1r+NOCuBf4Qoj+jnD/ieu/2roHo/dMsvCfr5qAZG+ZWlATuAB0J8bx5rJe7HgYL2/J7aI34f1mJjTCdT1R2qeijOh9OvgZeB3XA+2JeISJ5b9Xs4H3b/FpEk3wOoAT7E+YAO9kTQ9pPu8xG+AvcWyisisgNoxEme/o3z4bpb0PE7VfV/wRcRkXHu7YwCv3N86O6e2NbXIIwjgWWquiCofA5O8nRkUPnrqtrot73YfW4xqikK0sb+VcBO4M8icmmI1rO2L+Ac95WI1OAklfXAUYT++rX6XkVkN5yWuEdUtSbM9dLc878IVAX9TL2BkzSEvZXn3ipKB7a38rba+vnzfQ/nBNV7FidpOiqofJGqbvBtuO/tO0J/j59vJa7twGD3vZo+zhIbY7qIqn6hqn9W1R8Aw3Ca28fQPILDdwvhc5wPPf/HWUAegRpUdUdQWYH7PBCcPiw4rUXDcUZi+RIsXx+LfkHHt7htJCKZ7jkOAH6Dk2DtD5we5hyRGhjqejgjcJreg5+A96qqtR28vr9RftdtQVVLgcOBRcAfgW/cPja3tnZr0UdErgX+AcwHzsBJKPbHaVEJFX9b73WQ+7yplcsOxEmUr6Tlz9Mbbp3gnyl/vmuFTJyI4OfPfW5QpzN2E1VVt26r32NXLaG/Rq2NWKvBSVbTWqlj+gjLbo3pBqpaLyK34nRYneQWF7nPZwLrIzhNkogMDPpwyXeffWWn4vTPOF1Vm84pIpPDhRai7EicRGyGqvpaaXz/0XfEDpw+PcGGuc9FIfZ1OhGZhvN1e6S1euqMmDrb7TeyN84tkZtxbp/8qY3LnAd8oKqXBV07K0z9tvgShRGt1NmJ07o2F7gvTJ21rRzv+xkaEGZ/JD9/O9x6g/yTG/drmI+TxEertVFaA4BaVe3UuXZM72QtNsZ0MhEJ9eENzbcgfC0Fb+PcotjFbd1p8QhxjuD5V85xnz9wn31//H3/8fs+VC5px1tocQ7Xz0LUrXWvEUkrynvAHiIyNaj8Avea77cjxqiIM5z5AZxWjDsjOUYdX6vqNUAJ4B9/uNYFJejrJyJ743QwbjdV/Q5YDVwUogOur04VztdwCrA4zM9UqBYS3/F1OB2cd2kllLZ+/t5zn88LqncGTsL9Hl1jHE4HfGOsxcaYLvC2iGzCGUG0AucfiMk4o1QqgLsBVHWdiNwM3CYi43BuU+zEuUU1DahUVf+RLHXA/7m3ij4HDsK5VfSmqn7s1vmvW+8pEbkDp2n+MqB/O+L/1I3jAbeVqR7nA22fEHWXuM/Xi8ibOC0Gi90PyWB34iQxr7vvez3OSJaf43RY/a4dMUZiV3fIdALOLZADcEZ5ZQMXqOo34Q4UkZPcuF7C+bAXnFtxuThfY58lwAwR+T7OrZJyVf0WeA1nbpVbcfomTcBp7VlL9H93L8f5mfpMRO7EGaE1Cme+JF/CcRXwMc6oun/gdNDNAsYD31fV4H5MwT7AGW0VSqQ/f2/j9E3Kxuk0vzfO6LCFOK1JnUqcqQKm0UYLnOlDYt172R72iLcHzgiQJ3A6QZbjfCCsx+nAOzFE/VOA/wGlOH0F1uF0tjzKr84cnKRoL5z/yqtwmv3vBzKCzncSTt+Qapw+GXfgDNlWnNtLvnofAEvDvIfpOAlOJU7HzIdwWgIUmOlXL8Xdtx1nZJP/yKJ1BI0YwvkgfgLntlMdTuL3SyDBr84Y9zy/DBGXArPa+PrPoHlEmuIkZkXu+7kNGB3imJlBsU/A6Ri7yv1al+D0l7kw6Lh9cBKJSvf4D/y+Ln9xv/7VwJfu93kOsC7a94rTV+cNN54aN77ZQXXG4HzIb3K/xttxEoybIvjZPdK97v5B5XOI/OcvDedW3Tr3+lvcerlB9dYBr4WI4QPf1zHoe7NfGzFPjfXvvj16xkNUW7ttaYzpCcSZ6O5MVc2MdSwmvonIYuATDeof1FOJyFxgnKoeHOtYTM9gfWyMMcb4uw6YKe4Myz2ZiOyCM4Lw+ljHYnoOS2yMMcY0UdW3cGaqHhvrWCIwCrhCm/v4GGO3oowxxhgTP6zFxhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3LDExhhjjDFxwxIbY4wxxsQNS2yMMcYYEzcssTHGGGNM3Ph/bOxsVbjb+IUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "h2plus_distance = np.load('./h2plus_distance.npy')\n", "h2plus_lda = np.load('./h2plus_lda.npy')\n", "h2plus_overlap = np.load('./h2plus_overlap.npy')\n", "h_energy = -0.23839575845846386\n", "\n", "# Figure\n", "fig, ax = plt.subplots(1,1, dpi=100)\n", "\n", "ax.axhline(y=0, alpha=0.5, c=\"grey\", ls=\":\")\n", "ax.plot(h2plus_distance, h2plus_lda - 2*h_energy , lw=3)\n", "ax.plot(h2plus_distance, h2plus_overlap - 2*h_energy, lw=3)\n", "\n", "# Text inside plot\n", "# \"H$_2^+$\", fontsize=15, fontweight=\"bold\"\n", "\n", "ax.text(2,-0.05, \"H$_2^+$\", color=\"black\", fontsize=20, fontweight=\"semibold\")\n", "ax.text(9,-0.02, \"OA\", color=\"tab:orange\", fontsize=15, fontweight=\"semibold\")\n", "ax.text(9,-0.09, \"LDA\", color=\"tab:blue\", fontsize=15, fontweight=\"semibold\")\n", "ax.set_ylim(-0.15,0.05)\n", "\n", "ax.set_xlabel(\"Separation Distance (bohr)\", fontsize=12)\n", "ax.set_ylabel(\"Binding Energy (hartree)\", fontsize=12)\n", "\n", "# ax.legend()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "### H2" ] }, { "cell_type": "code", "execution_count": 83, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(-0.2, 0.12)" ] }, "execution_count": 83, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAi0AAAFuCAYAAABNzSgtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAABpsUlEQVR4nO3dd3hb5fXA8e/R8F7Zzl6QhCxCgBDC3ntTaMsKUH4to5TRQltaCqWM0pZSVikzJezVsjcNK4EA2SSBhOzhOMt7yJLe3x/32taVJQ9ZsizlfJ5Hj33f+96rIy8dv1OMMSillFJKdXeuZAeglFJKKdUemrQopZRSKiVo0qKUUkqplKBJi1JKKaVSgiYtSimllEoJmrQopZRSKiVo0qKUUkqplKBJi1JKKaVSgifZAaSLrVtKBRgAVCY7FqWUUioF5QOb+vTrG3XVW01a4mcAsCHZQSillFIpbBCwMdpJTVripxKgR6+euFyJ73Wrum138qSu6XjF0f9m970PT/jzKqWUUvEWDAbZuX0HtNFboUlLnLlcri5JWur9AQqoajpuqKnokudVSimlkkXf5VJUnWQ5jv11VVFqKqWUUulBk5YU5QtLWgL11UmKRCmllOoamrSkKJ/LmbQENWlRSimV5jRpSVEN7mzHcdCnSYtSSqn0pklLivKHJS1o0qKUUirNadKSogItkpaa5ASilFJKdRFNWlJU0JPjOBa/Ji1KKaXSmyYtKSrodba0uBo0aVFKKZXeNGlJUcbrbGlxBWqTFIlSSinVNTRpSVXeXMehR7uHlFJKpTlNWlKUZIQlLYG6KDWVUkqp9JCSSYuIXCYiq0WkTkS+FpGDWqnbX0SeFpFvRSQoIndHqXeGiCwVkXr742kJewFx4Mp0Ji3eoHYPKaWUSm8pl7SIyNnA3cCtwF7AJ8BbIjIkyiWZwFa7/sIo99wfeA6YCexpf3xeRPaLa/BxFJ60ZAS1pUUppVR6S7mkBbgGeNQY84gxZpkx5ipgPXBppMrGmDXGmF8YY54AyqPc8yrgPWPM7caY5caY24EP7PJuyR2WtGQaTVqUUkqlt5RKWkQkA9gbeDfs1LvAtE7cev8I93ynk/dMKG9WnuM4S5MWpZRSaS6lkhagN+AGtoSVbwGKO3Hf4o7eU0QyRaSg8XHX3+/K78Tzd5gnW5MWpZRSu5ZUS1oambBjiVCW6Hv+Bqu7qRwov+/++zd08vk7JCM8aZEGTMDflSEopZRSXSrVkpZtQICWLSB9adlS0hElMdzzdqCw8XHF5ZcP6sTzd1hmTsuGnbqaqq4MQSmllOpSnmQH0BHGGJ+IfA0cBfwn5NRRwCuduPUc+x5/Dyk7GpjdSiz1QH3j8dYtpZ14+o7LzMlrUVZbU0F2flGXxqGUUkp1lZRKWmx3ATNF5CusZOP/gCHAgwAicjsw0BhzfuMFIjLJ/jQP6GMf+4wxS+3yfwAfi8j1WMnPKcCRwIEJfzUxys4taFFWV60tLUoppdJXyiUtxpjnRKQXcCPQH1gCHG+MWWtX6Y+VxISaH/L53sCPgbXAMPues0Xkh8CfgFuA74GzjTFfJOp1dFZ2dh5BI7ikediNr64yiREppZRSiZVySQuAMeYB4IEo56ZHKJN23PNF4MVOB9dF3G4X1WSQ29xDhU/HtCillEpjqTYQV4WokyzHcUOdJi1KKaXSlyYtKSw8afFr0qKUUiqNadKSwnxhSUugrjpJkSillFKJp0lLCvO5nElLsF6TFqWUUulLk5YU5nNlO46D9do9pJRSKn1p0pLC/G5nS4tpqElSJEoppVTiadKSwgLuHGeBT5MWpZRS6UuTlhQW8Di7h9CWFqWUUmlMk5YUFvQ6W1pcmrQopZRKY5q0pDATlrS4/Zq0KKWUSl+atKSy8KQlUJukQJRSSqnE06QlhUlGruPYo0mLUkqpNKZJSwqTDGdLizeoSYtSSqn0pUlLCnNlOltavIG6JEWilFJKJZ4mLSnMk5XnOM40mrQopZRKX5q0pDC3Ji1KKaV2IZq0pDBvWNKSrUmLUkqpNKZJSwoLT1qy8GGCwSRFo5RSSiWWJi0pLDPHmbS4xOCrr05SNEoppVRiadKSwrJyClqU1VVXJiESpZRSKvE0aUlh4S0toEmLUkqp9KVJSwrLyc1vUVZfW5WESJRSSqnE06QlhXm9XuqM11Hmq9GkRSmlVHrSpCXF1UqW47ihTruHlFJKpSdPsgNQnVNHJtCcqDTUaUuLUkqptgWDhtqGANU+P9X1Aarr/dT4rOMa+7jaZ5fV++3jADV2/aG9cvjjKeO7NOaYkxYRGQUMBLKBbcC3xpjyeAWm2qdessA0HwfqdMqzUkqlo4ZAsClxqK73U9WYSNT7qapvLquJkIRY9eyEw9d8rjPGDWg5gzXROpS0iMhU4GfAcUDvxmKst82giCwAngRmaALTNXyubAj5uQvoOi1KKdVt+ANBKuv8VNb5qahroKKugco6Z7JR4wtNQAJNn1fZLR2NZT5/91o8tLre3+XP2a6kRUQmAXcDBwPLgBeBeUApUAf0BEYAU4FbgVtE5DbgLmOML+5RqyYNrixH0hKs1+4hpZSKB2MMNb5AU8JRWddARW1j8uGn0k5AKmobQuo4jzvbmtGdVSfhtbW3peUL4CngGmPMvNYqikgu8EPgOvv+f+pUhKpVDe5saGg+Nr6a5AWjlFLdVCBoKK9tYEd1PTuqG9hR7WNHtY+dNT62V1kfG8vKa5tbRAJB0/bN04TbJeRmuMnN9JCT4SYv00NOhofczMYyD7kZbnIyPeRluinKzujyGNubtIw3xqxoT0VjTDXwqIjMAAbHGphqn4A723FsfNo9pJRKb8ZYA0jDk43GJCT02CprYGeND5Nm+Ue210om8uykIjckwQhPOHIz3PaxVdZ0zi7LyXCT6XEhIsl+Wa1qV9LS3oQl7JoAsKaj16mOCXicSYs0aEuLUip1+QNBtlbVs7m8jpLyOjaX17G5rJbNFdZxSXkd26rqqe9m4zvaw+0SCrI85GV5yMv0NicbdlLRmGw0luXZSUVepqepdaO5vge3q3snGInQmdlDhVhjWHoDbxpjdsYtKtVuQU+O49jl16RFKdU9NQSCbKloTkaaPlbUsqnMOi6trKO79sjkZrgpyPaSn+UhP8tLQePH7MZj61xjnQK7TuNxttfd7VsyuruYkhYR+T3wa6zpzgbYF9gpIh8A7xlj7ohfiKo1QW940lKbpEiUUrsyYwzbqnys3lbN5vLakKSktik52VpV3y26aPIzPfTIzaBHbga9cjPokZNBrzzrY89cL0U5GU0JSKGdcORlevC4dT3WZOtw0iIilwF/AB4A3gLeCDn9OnA6oElLVwlLWtza0qKUSiCfP8i6HdWsLK1m1bYqvi+t5vutVXy/tYrKuq6fAut1i51sWI/QRKSxrGdIYlKU4yXT4+7yOFV8xNLScgXWVObrRCT8O78C2L3zYal2y8h1HHoC2tKilOq8ndW+FknJqq3VrN1R0yUzajLcLooLsyguzKK//XFAYTbFhVn0yc+0EpPcDPIzPdrlsguJJWkZAbwT5VwlUBRzNKrDXJnOlhZPoC5JkSilUo0/EGTDztqmhKQxOfl+azU7qhO3xFamx8WAomyKC5oTEutjNv3tz3vmZuxyycj6HTUcdOf/mo73G96T5366f6vXHHDHh2wsa/5nVQS8bhd5mR5652WwW988DtytD6ftNZDsjLZbmBoCQabc+j47a5rX0vC4hLk3HEnP3K6f4hwulqSlHOgX5dwwrAXnVBeRsJaWjKC2tCilWtpcXsv8dWV8s6m8KUFZs60GXyC+s3AyPC4G9bCSj+KCbAYUhSQlBVZ5UY53l0tIuooxVhfeDr813fu7LVW8ubiEO99Zzq2nTuCEif1bvf7j77Y6EhYAf9DwxqJNnLf/sARG3j6xJC0fANeJyCtYq+ECGBHxAJcSvRUmbuxxNb8C+gPfAFcZYz5ppf4hwF3AOGATcKcx5sGQ89OBxyNcmm2M6dZNF57MPMdxRrBbh6uU6gJ1DQG+2VTOvLVlzF+/k/nrythcHt+/DX3yMxnRO5eRffMY2SePkX1yGdknj4FF2bh2wam43cVx44sJBA3bq30s2VjeNDW8rKaBy5+eR0nFWC4+cHjU619ZsClqeaomLTcCXwJLgf9gzR66AtgLGAKcFbfoIhCRs7G2FLgM+Az4KfCWiIw1xqyLUH848CbwMHAucADwgIhsNca8FFK1Ahgdem13T1gAPFnOpCWz+4eslIojYwwbdtYyf30Z89buZP76MpZuKqch0PlxJ163MLRXblNCMsJOTkb0yaMw2xuH6FW8/fPcvZs+r6xr4J4PVvDwJ6ubym57cxmThxSx15AeLa6t9QV4f9mWpmOvW5p+jr5et5MNO2sY1COnxXVdqcNJizFmpYgcgNVycRnWhonnA/8DzomUOMTZNcCjxphH7OOrROQYrFae30So/zNgnTHmKvt4mYjsA/wSCE1ajDGmJEExJ4wny9k9lIUmLUqls1pfgEUbyhxJytbK+k7dsyjH62gtsRKUXAb3zMGr03xTVn6WlxtOGEvQwKOfWolLIGi4+/0V/PuiKS3qv7u0xLFX0s8OGcm9H64ErG6nVxdu4rJDd+ua4KOIaZ0WY8xS4FgRyQR6ATuNMQkfTCEiGcDetJxS/S4wLcpl+9vnQ70DXCwiXmNMY+ddnoisBdzAAuD3xpj5rcSSCWQ2Hv/6+uvzr7n6mva+lLjxZjtbWrJNvfXTpf3FSqU8Ywxrt9cwf/3Opq6eZZsrY569k5/lYdLgIvboX+Do2ukOAyxV4vz88N148vO1TV1Fn67cRmVdA/lZztay0K6hnAw3lx+2G68s2MS6HdZSGq8uSNGkJYQLq6Wloa2KcdIbK6nYEla+BSiOck1xlPoe+36bgeXAdGAxUAD8AvhMRPZsZQuD32CtVwPAffffTzKSlsycfMexR4L4fXV4MrOjXKGU6q6q6/0sXO9sRYl1Fo8IjO6Xz152V8DkIUWM6J2n4012QUU5GUwYWMhXa62F6wNBw5KNFew/sldTnbIaH5+s2Np0fNiYvmR53Rw3vph/fbwKgOUllXxbUsnoYuf7TleKdUXcw4DbsFbCBZgCzBOR+4EPjDEvxym+aML/zZAIZW3Vbyo3xnwOfN50UuQzYB7wc+DKKPe8HauLDIArLr88H9jQVuDxlpHT8oenpqaSAk1alEoJpZV1vL+0lHeXljB75faYZ/P0yPGy15Ae7DW4iMlDezBxUGGL/6TVrqu4MMtxvL3a2aX4+qLNjnFQx4232gGODUlaAF5ZsJHrjh2TwEhbF8uKuIdjda8sAf4KXBdyehtWi0WikpZtQICWrSp9adma0qgkSn0/sD3SBcaYoIh8SSsL5Rlj6oGm7/rWLcmZ6Z0TIWmpr6mEHn2TEI1Sqj2+31rFu99s4d2lJSxYX9bhpe3dLmFMsdWKMnlID/Ya0oNhvXJ0GrGKquV/7s6flVdDuoayvC4OH2O9h+w1pAcDi7Kb1oJ5deGm1EpagD9ibZB4ij3NOTRpWQhcGJfIIjDG+ETka+AorJlLjY4CXoly2RzgpLCyo4GvQsazOIj1mz8Jq7uoW8vKjZK0KKW6jWDQsHBDGe8u3cK735Tw/dbqDl3fOy/DakWxk5SJgwrJyehs777alWwqcw47DR3HtKmsli/X7mg6PmRUH8fP1zHjinnsM2sg74adtXy9dgd7D+2Z4Igji+Wnfi/gB/bn4cnbVqxWjES6C5gpIl9hJST/hzXV+kEAEbkdGGiMOd+u/yBwhYjchTXteX/gYuBHjTcUkT9gdQ+twBrTciVW0nJ5gl9Lp2VkZOIzbjKkecS3Ji1KJZ/PH2TOqu28+00J7y3dQmk7Z/h4XMK4AQWOJGVQj2xtRVEx21nt45uNFU3HbpcwfmBB0/GrCzc5Wvs+W7mdqbd90HRc29D8/gLWgN1USlr8QLSO0r5YS/knjDHmORHphbVeTH+sbqrjjTFr7Sr9sZKYxvqrReR44O9YScgm4MqwNVqKgIewupHKgfnAwcaYuYl8LfEgItSRRQbN/7k11FYlMSKldl0VdQ3M+nYr735Twqxvt1JV374NBPMyPRw2pi9Hje3HoaP7UKBjUVQc/eODFY6xUgfu1tsx3il8Qbmqen+rP7tvLNrMjSeOTcqu17EkLV8C5xG5O+ZMrNaPhDLGPIC1y3Skc9MjlH0ETG7lflcDV8crvq5WK5kUhCYtdR1relZKxW5LRR3vLd3Cu0u3MOf7be1e1K1vfiZHje3H0eOKmTqip+48rOKusq6Bf7y/ghmz1zSVuV3CVUc2D9dcWVrJss0VEa6Obnu1j09XbuPQ0V0/djKWpOUO4B0R+Q/wBFYX0X4ichFW0nJYHONT7VAvWY6OOn+ddg8plUgrS6t4d2kJ736zhQXry9p93cg+uRw9rpijx/Zjz0FFOv1YtWpFaRWXPvl1xHN7Di7iZ4eMbFF+6ZNfEzSG7VU+Focs49/o9yfs4VgN97/zna0s9/14L06cOKDFfeev28lpD8xuOn51wabUSFqMMe+LyAVYS+mfYhffD5QB040xn8YtOtUu9a5sa06VLaAtLUrF3Zpt1Tz75Xre/aaEVdva/zs2eUgRR48r5qix/RjZJ6/tC5Sy7aj28daSyAu1+6MsMBitflGOl9tOm8DxE5wbJr62qDlpyfC4oiYikwYX0b8wq2kPq3eXbqGuIUCWt2tbCDuUtIiIGxgJvI61BP40rB2ftwGfGWP03TIJGsKSlqBPvw1KxYMxhjmrtvPYp2v4YPmWdk1NznC7mLZbL44eW8yRe/Slb0FW2xcpFSci4HW5yMl00ycvk5F98jhw996cPnlgixln89ftZO32mqbjA3frTV5m5LRARDhmXHFTV1NVvZ/3l22J2CqTSGI6sECAPcW5DjjJGPNWwqJKQVu3lBYA5b369Mbl6trBSQvvOII9675qOv5i5FXsd97NXRqDUumkriHAqws38dinq1le0nZ3a749kPbocf04ZFQfXdRNqQ4KBoNs37oNoLBPv75RB9l0qKXFGOMXkRKs5ftVN+F3O1e/NdrSolRMSivrePLzdTz1+Vq2t7F8fr+CTI4eW8zR4/qx3/BeZHj0z6JSiRbLQNxnsXZ1fiPOsagYBdxhW4X7ayJXVEpFtGRjOY99tprXFm5qdfZPfpaHs/cZzEl7DmDCwEIdSKtUF4slaVkAnC0iH2It17+ZsEXmumDvIRUi6HW2tLgaNGlRqi2BoOH9ZVt49NPVzF29o9W6w3vncuEBwzhj8iByo/T5K6USL5bfvifsjwOBQyOcN1g7MasuEvQ4W1o0aVEqusq6Bp7/agMzZq9m/Y7aVusesFsvLj5wOIeO6qutKkp1A7EkLYfT+o7Kqqt5nUmLO9D6H2KldkVrt1czY/YaXvhqQ6urfWZ4XJw2aSAXHjiMMcUFUesppbpeLOu0zEpAHKozMnIdhx4d06IUYE1Z/nzVDh77bDXvL2t9ynKf/EzOnzqUH+83hF55mV0XpFKq3TqctIjIKuA0Y8zCCOfGA68aY0bEIzjVPpLhbGnxaEuL2sXV+wO8tnAzj326mqVtLFE+fmABFx84nBMmDNAZQEp1c7F0Dw0Dov0bkgUMjTkaFRNXpnOVzYxgXZIiUSq5tlbW89QXa3ny87Vsq4o+ZdklcMy4Yi46cDj7DO2hOygrlSJiHQYfrZF1BAne5Vm15M50dg9lGE1a1K5lS0Ud93ywghe+2uDYzTZcfqaHs/cdzAXThjG4Z07Uekqp7qldSYu919AFIUX/FJHwNtdsYE/gozjFptrJk+VMWjK1pUXtIsprGnjw4+95/LPV1DVET1aG9srhwmnDOHOfwVGXKVdKdX/t/e3NAfqEHBfRsouoHngO+EPnw1Id4clydg9loUmLSm91DQEe/2wN/5y1koq66DOBpo3sxUUHDOewMX1x65RlpVJeu5IWY8w/gX8CiMhq4IxIA3FVcniz8x3HWdo9pNKUMYY3F5dw25vL2FgWecB5htvFKZMGcOEBwxk7QKcsK5VOOrrLczbwOZDfVl3VdTKywwbiSoBAgw+3NyNJESkVf0s2lvPH15Yyd03k1WvdLuGsfQbziyN2p7hQd1ZWKh11dMPEWhE5CbvVRXUPmWFJC0BtTSV5hb2SEI1S8VVaUcdf3vmWF+dtiLrOyokT+3PNUaMY0afl74JSKn3EMiJtITAe+DjOsagYZeW0bPiqq67QpEWlNJ8/yCOfruL+D1dS7QtErDNleE9+f8JYJgwq7OLolFLJEEvS8mtgpoh8Y4zRmULdQHZuy377+hqdea5S15drdvCblxezsrQq4vlBPbK54fg9OHZ8sa6xotQuJJak5QEgD/hQRHbScpdnY4zZMx7BqfbJysrCb1x4pHnKZ31t5D/2SnVn5TUN3PH2Mp6Zuz7i+dwMN5cfvhsXHTCcLK/uy6rUriaWpGU7sC3egajYictFLZnk0zybwqdJi0ohxhheXbiJW15fGnElWxH4wd6D+OXRo+lboINsldpVxbJh4qEJiEN1Up1kOZIWvyYtKkWs31HDDf9dwsffbY14fvzAAm47bQITBxV1bWBKqW5Hl4ZME3WS5eik89dp0qK6N2MMT36+ltveXE5tQ8uBtjkZbq49ejQX7D8Uj1s3MlRKdSJpEZFCYBTW8v0OxhidWdTFfGFJS6C+OnnBKNWGzeW1XPfiIj5ZEbmn+cg9+nLzKeMZWNTiz4tSahfW4aRFRDzAg8D5QLSRcDpCrov5XFkQsvVKUJMW1Q0ZY3hlwSZufGVJxOX3+xVkcvPJ4zhmnM4KUkq1FEtLy9XAScBFwBPA5UADcAlQCFwZt+hUuzW4syHkPSDo06RFdS87qn387r+LeXNxScTz5+w3hOuPG0NBlreLI1NKpYpYkpbzgFuBZ7CSli+MMfOAR0TkHeAw4N34hajaw+8Oa0bXpEV1I/9bXsp1Ly1ia2V9i3P9CjK588w9OWRUnwhXKqVUs1iSlhFYq+I2dkaEzj98EPgH8JtOxqU6KNAiaalJTiBKhfD5g9z59nIe+XR1xPOnThrAzSePpzBHW1eUUm2LJWmpBjKMMUZEdgBDgdn2uVpA145PgqAnx3Esfk1aVHJt2FnDFU/PZ8H6shbninK83HrqBE6Y2L/rA1NKpaxYkpblwHD789nANSLyCeADrgO+jVNsqgOCXmdLi6tBkxaVPO8v3cK1LyykvLahxbnDx/TljtMn6CJxSqkOiyVpeQ5rqjPAH7A2TlxrHzcAp8chLtVBxutsaXEFaqPUVCpxGgJWd9DDn7TsDsrwuPj9iWM5d78hOjNIKRWTWFbEfSDk8/kiMhY4DWuMy3vGGG1pSQZvruPQo91Dqott2FnDz5+Zz/x1ZS3ODe2Vw/0/nsz4gbobs1Iqdp1eEdcYsx64Jw6xqE6QjLCkJVCXpEjUrmj299u47Kl5lNW07A46YUJ/7jhjAvk6lVkp1UmdSlpEpA+RV8Rd15n7qo5zZTqTFm9Qu4dU4hljmPn5Wm5+bSmBoHGcy3C7+P1J2h2klIqfWFbEzQf+DvwI53TnULoibhcLT1oygtrSohLL5w/yh1eX8Mzc9S3OaXeQUioRYmlpuRv4MfAosAhouVpUgonIZcCvgP7AN8BVxphPWql/CHAXMA7YBNxpjHkwrM4ZwC3ASOB74AZjzH8S8wrizx2WtGQaTVpU4myrqufSJ7/myzU7W5w7co9+3HX2nrqyrVIq7mJJWk4Afm2M+Ue8g2kPETkbK3G6DPgM+CnwloiMjdQtJSLDgTeBh4FzgQOAB0RkqzHmJbvO/lizon4P/AdrYPHzInKgMeaLxL+qzvNm5TmOszRpUQnyzaZy/u+Jr9lY1rIL8ueH78bVR47C5dLuIKVSTjAIwQYINDR/bPrcb3/0NX/uzoCBk7s0xFiSlixgcbwD6YBrgEeNMY/Yx1eJyDHApUReifdnwDpjzFX28TIR2Qf4JfBS4z2wZj7dbh/fbrfOXIXVDdZuwUAAEWnqww8GgxhjEBFcLldTvUAgAIDL5YpLXU929KTFGEMwaC1g7HY399xFum9X1430mrtb3ba+Rx2pG+3r0x3qtuf7+daSEq59fiG1DQFCZXld/PUHe3LcuH4YE8QY/TlJ9vcz3t/7tPsbEQhgggGEIC4MBP2YoJ+g3w8mgFuAYABMgKC/ARP0I8auawJW3YAfggFcGMQErboBv7Nu0Kpr7LpCMKSunQyYIC4TgKC/6fmw7yF2uQn6rXiDfrvMegQDdl2CSNCu25RgNCBNSYgfE5KMiP25CUlGxITsutsefcfCZXPi870POP+mRBNL0vImcBDwYQzXdoqIZAB7A3eEnXoXmBblsv1puRfSO8DFIuI1xjTYdf4eoc5VrcSSCWQ2Hv/6+uvzr7n6GmbPmcPUqVPJyMgAYP369axevZr+/fszevTopus/++wzgsEgU6dOJSvLGhq0adMmVq5cSd++fRk7dmxT3c8//5yGhgb23XdfcnOtbqCSkhK+++47evfuzfjx48kIT1rE+mEUt4fS0lKWLVtGjx492HPPPZvqfP3119TU1DBp0iSKiooA2L59O0uWLKGgoIDJk5sz6AULFlBZWcmECRPo1cta9Hjnzp0sWrSIvLw89tlnn6a6ixcvpqysjLFjx9K3b18AKioqmD9/PtnZ2ey3335NdZcsWcKOHTsYM2YMxcXFAFRXV/PVV1+RkZHBtGnN39Zly5axdetWdt99dwYOHAhAbW0tc+fOxePxcOCBBzbV/e677ygpKWHEiBEMGTIEAJ/Px5w5cxARDjnkkKa6K1euZNOmTQwbNoxhw4YB1i/Wp59+CsDBBx/c9Mu1evVq1q9fz+DBgxk5ciRg/cJ98onVO3nggQfi8Vi/VuvWrWPNmjUMGDCAUaNGNT3fp59+ijGG/fffn8xM60dow4YNrFq1iuLiYsaMGdNUd86cOfj9fqZMmUJOjrUWz+bNm1mxYgV9+vRh3LhxTXW/+OILfD4f++yzD3l51s9DaWkpy5cvp2fPnkycOLGp7ldffUVtbS177bUXhYXWuJOtW7eydOlSioqKmDRpUlPd+fPnU1VVxaL63tz1v7WEG1CYxUPn78P4gYXMmzePiooKxo8fT+/evQEoLy9nwYIF5OTkMGXKFMf3fufOneyxxx7069cPgMrKSubNm0dWVhZTp05tqrt06VK2bdvGqFGjGDBgAAA1NTV8+eWXeL1eDjjggKa6y5cvp7S0lN12241BgwYBUF9fz+eff47L5eLggw9uqrtixQo2b97M8OHDGTp0KAANDQ3Mnm0t8n3ooYc21V21ahUbNmxgyJAhjBgxArDeJBu/9wcddFDTH+c1a9awbt06Bg0axG677dZ0j8a606ZN6/K/EY2+/PJL6urqmDx5MgUFBQCR/0YYw7wvv6CuqpyJ48ZQmJsFAR9l27aw8tul5GdnMGb3EeCvh4CP1d8upb6mgiED+5OfnQlBPzVVFWzasI4sr5vBA/vb/6H72bp5I77aanr3LCInKwMCDTTU1bB9WyleF/TuWWS9gQf8VJXvwF9fR25OJpkeNwT9BBt81FRVIGLIy85qSi4a6mox/gY8HhcuARrf8H0+IAghZa4Ib9BC5AGZrghlHa0bqe0xUt143CNaO2ek8k61iQas2YKNfyMmTpxIz549AdixYweLFy8mPz+fvffeu+mShQsXRvwbsWjRIvYYPablc4RpV9IiIj1DDv8EvCgilcBrwPbw+saYHe25bwx6Y/2cbAkr3wIUR7mmOEp9j32/za3UiXZPsFp1/tB4cN/993PN1de0FntCZebktyirq6kiO7+o64NRaSVoDM8sr+edNS0TljG9vMz82YH0yc+McKWKO2OQYAOu+nIoWw8NNeCrpmDHIoyvGvfyzeAKgK+avO0lDNuygbwNwKpsu24NY7ZtAn89ucszgAD4ffTy1bB/XQ0uAvBmsCkR2Rd7RtinzSH0APZtPAgZSdiUmi1tLssFdm88WNZc3vSHNWQMdwbWIEXAGnloK2j8JGT4lBto+otX2Vze9FMYMtJSiMPaHiqyYMslDhJNjDFtVxIJAqEVJezYwRiTkNlDIjIA2AhMM8bMCSm/ATjPGNMiTROR74DHQ7p+EJEDsH4N+xtjSkTEB1xgjHkmpM45WN1QEWdIRWlp2dCjZw/cHk+XN/1u3bKBPv9s/q8bYMeli+nZb8iu2fSr3UMdrhvp61PvD/DL5xfy2qLNhDtr70HcdPJYcjKbB9zqz0kbX3dfHdTuxFVfhtTuhJodBKu3Qc0OpHYHUrcT6srBV4PxVVu7tTfUIHZyYnzVVteAUkliXB5wZyAuLxQNgUs/jcvvfcDvZ+eOnQCFffr1rYj2/O1NQP9IK0lKF9oGBGjZAtKXli0ljUqi1PfT3EoUrU60e2KMqSckn9+6pRQAl9vtWJMi9A9WqNBvYjzqZucWtKhTV10FgIhEvEeksq6uG+k1p3NdiP/3Pl51w19HRV0DP5v5NbO/b9GYyq+OGc1lh45ssf7KLvVz4qtBanfgrtkBNduhdgfU7LCSkhZlO6BmJ+768pb3bVFix9HOMpUE4gJxg8sd8tEVduwGV6R6drnLax27PM0fJfS48XP7Ia6Q4/Zc5wV348cMcHutcrd93Pi5y2uXRfncTlIaP4+05lI8fu9dEcojaVfSYoy5qV13SzBjjE9EvgaOwprl0+go4JUol80BTgorOxr4yh7P0ljnKJzjWo6meffqbi87O4+gEVzSnFv66ipbuUKp6LZU1HHBY3NZXuL8GfK4hDvOmMiZew9KUmRdpKEOytZB2VrYucb+uNYqq95qJST+XXmGnoAnE9yZ4MmwPrq9dlnIm2L4w+11vqmGvvG6Q4+9ON6g3WHH4W/UjoShtbLGhMETIYkIL4uQiDQmD7pYYtKkYlffXcBMEfkKK9n4P2AI8CCAiNwODDTGnG/XfxC4QkTuwpr2vD9wMc5ZQf8APhaR67GSn1OAI4EDSRFut4tqMsgN6cz11VQlMSKVqlaWVnHBY3NbTGnOyXDzwDmTOXR03yRFFkcBP1RssJKQnWtDkhL7Y1VJsiNsP3cmZORY+49l5Eb5PMc69uaCN9tKLhqTjqZkozEBsR+tlbnc+satkqK9A3GvA+41xrR7bXgR2RsoNsa8EWtwkRhjnhORXsCNWOO2lgDHG2MaRwn2x0piGuuvFpHjsVpRLsca4nVl4xotdp3ZIvJDrEHGt2AtLnd2qqzR0qhOshxJS0OdJi2qY77ZVM65j3zBzrA9hHrnZfD49ClMGJQiK9waA1VbwhKSNc3H5RutWSTJlFkA2T0gpydk97Q+5vSCrKI2ko888OY0l7lT8X9PpWLT3p/2C4GrReRR4EljzPJIlUQkC2vxuYuwWioujEuUYeydph+Icm56hLKPgFZXwDHGvAi8GI/4kqVOssA095n7NWlRHbBwfRnnPzaX8lpnwjKsVw5PXLQfQ3rlJCmyVvjrYeu3ULoUtnwDW5fb3Tnruq77RlxW8tGYeGTbyUdOpDL78+weVguGUqpD2pu0TMBqpfgl8BsRKQXmAaVAHdATa/n7CfY93wAmG2O+iXvEKiqfZDmGSwfqqpMXjEopX6/dwfTHvqSy3u8o33NQIY9N35deed1gSnNDLZQsgU3zYfMC6+PWbxPbYpLTC4qGWrMkegy1Pi8YGJKA9LBaRqIMeFZKxVd7B+L6gX+IyH1Y4z2OxxobMg1rl+dtwHKsrpWnjTGrEhOuao3PlQUh6yUF6zVpUW37fNV2LprxJTU+55v/1BE9efSCfcnNTEL3QzAI276DdbNhw9d2grI8/glKRp6ViDQmJI6PQyCz5fpHSqnk6dBfI2NMAHjZfqhuxufKdhwH67V7SLXu0xXb+MkTX1LX4Fwd9KDde/PQefuQnZGQJZda8vus1pO1s2Hd57D+c6htuRljh7kzrOQjUkJSNMxqLdEBpUqlDB3BlUb8buc6eKahJkmRqFTwv+Wl/PTJr/H5nQnL4WP68sA5k8nyJjBhMcYah7LyfVj5Aaz/onNjUIqGQN9x0G8c9N69OTnJK9auG6XSiCYtaSTgDhso6dOkRUX27jclXP70PBoCzjUjjxnXj3t/NJkMTwLe6Gt2wPcfNj8qW66y2yZxQe/RMGAS9J9kfew7FrJaLq6olEo/mrSkkYDH2T2EtrSoCGZ9WxoxYTlpzwHcddaeeN1xTFjKN8Cy12DpK1a3T0cX1i4cDIP3g0H7wIC9oHiCNQVYKbVL0qQljQS9zpYWlyYtKswXq7bz05lft0hYTp88kL+cuSduVxzGd5Sts5KUpa/Ahi87dm3fsTBkKgyZZn0sGtz5eJRKNQ21sOAp+O4dKFlstVJ6MiGvn/V7secPYVg71j5d+io8f56zbPQJ8KOnExN3F9CkJY2YsKTF7dekRTVbuL6Mi//9FfVhY1jO3mcwt58+AVdnEpaaHbDoOVj0PGya1/7rCofAbofDyCOsP8I5Pdu+Rql0tu5zeOFCqNzkLA/UQ30FbF8B82daycdpD7beNbr4hZZlK9+zBrln94hv3F1Ek5Z0Ep60BNq9gLFKc8tLKjj/sblUha3DcuqkAbEnLMEgrJ4F82bC8tch4Gv7Gm+OlZyMPAJ2OwJ67aazd5RqtH4u/Psk5+9Sbl9r7FZ9ldVyGbQXf/z2DXjiFLjobasVJlxdBax4t2V5wGe1wOx9QUJeQqJ1OGkRkdeA+4wx7yQgHtUJEtbX79GkRQGrtlZx7iMtV7o9Zlw//vqDPTuesJRvgPlPwYInra6gtmQVWv8Vjj0FRhwK3qw2L1Fql+Ovt1pYQhOWqZfBUX+09ocCaxuKp8+Grcus403z4H+3wVE3t7zfstecM/Jc3uaEZ/ELu07SAuwBvCki3wP3A48bYyriG5aKhWQ4W1q8QU1adnUbdtZw7iNfsK2q3lF+8Kg+3POjvfB0ZNDt+i/hs7th+Ru0OaA2uweMORHGngrDD9Yl65Vqy6LnrE08Gw2ZBsfe7qzTYyicPRMemApBu9X0y0fgoGusfw5ChXYNZfeEcafBV49ax2s/g4rNUNA//q8jwTo8TcAYsxtwErASa8fljSLyoIhMiHdwqmNcmc6WFm+gi/ZeUd1SaUUd5zzyBZvKnT8HU4b15F/n7k2mpx3rsASD8O3b8Nhx8OiRVjdQtITFnQnjz4RzX4ZfroBT7oPdj9SERan2+PYt5/GUn0Su13t3q8Wyka8K1nzqrFNVCqs/bj4eczyMP6P52ARhyUukopjmNhpj3jTGHA+MAh4GzgIWiMgsETlTRLpoGU0VypOV5zjONJq07Koq6xo4/7G5rN3uHIw9cVAhj05vx0q3fp/VBfTP/eGZs63l9KPpNx6OuxOuXQ5nPmqNVWlszlZKtc/mhc7jQVOi1w0/F37tkpedW16MPRWG7G+Nj2kUaZBuCujUQFxjzPfANSJyC9YOyYcBBwGbROROrLEvHVyYQcXKrUmLAhoCQS57ah7LSyod5aP75fPvC6eQn9VKQhFogK9nwCd3tZy9ECojHyacCZPPt9ZP0cG0SnVOzXbncV6/6HXz+jqPq7c5j0MTksxCGH6ItTL0HifCV49Z5ZsXwLaV0Hu3mENOhk4lLSIyCPgZ8BOgD/AW8BzWpop3A6OBKzoXomqvjBxnn2a+qbKWS9c3lF2GMYbfvLyYT1Y4/4gN65XDzIun0CM3SleNMVbXz/s3wfaV0Z8gfwDsfxlMvkBXoVUqacLaAkL/xu9YDRu/aj4efVxzF+3YU5qTFrCSm8N+k7gwEyCm7iEROVxEXgZWAVcCLwBjjDEnGGOeMMacAVwDnBO/UFVbcnsUO44zJEBtZVlyglFJcff7K3jx6w2Osl65GTxx0X70LYgya2f9XHjsWHju3OgJS58xcOo/4RcLYdrPNWFRKt5yejmPq7ZEr1u1Nfq1i190nht7SvPnww4Kq5t6XUSxTHlehjWWZTVwHfBYlNlDXwCFEcpVghT2ajkSvGz7JrILUnMRIdUxz3+1nn98sMJRluV18ej0fRnSK6flBdu/hw9utlaujWboAXDAL2C3o3TjQaUSqf+eULGx+XjD3OgrQm+Y2/LaRkvCkpbXr4I3rm0+rg/pNt7xPWycBwMnxxRyMsTSPbQRK1l5vY3xKvOA4TFFpWJSUFhIrckgW5rn+VfvKIHh45IYleoKH3+3ld++vNhR5hK490eTmTS4yFm5oRZm3QFz7mueNhluwGQ4+pb2LRWulOq8UcfCt282H3/5mHPGT6NtK2HVrOZjb27z7+nmRbB1ubN+ay02YLXMpFDSEsuU5yONMa+1NcDWGOMzxqyNPTTVUSJCmTib7WvL2viBVSlv6aYKLntqHv6g81fy5pPHcdTYsMF8a+fAgwda661ESliKhsKZj8FPPtCERamutOcPoWBg8/HaT+GdGyAQ8ntatg6eP9/5u7vvxc1rtMTS3fPNy9bSBilC23vTTKW7yHHsKy9NTiBR3HTTTYiI43HooYe2es306dNbXDN9+vQuibe721RWy4UzWi7P/9NDRnDe/sOaC+qr4M1fwePHRR63kt0DjrkdrvjS+u9Ou4KU6lqeTDjjUWvl2kZz7oO/j4WnzoLHT4B7JkPpN83n++8Jh91gfW4MfPOf5nOZhfC7rXBTecvHPhc116vcDGs+Sexri6MO/2USkaCIBKI8/CKyTUTeFpHDEhGwal2N1zl+JRA+YEuljRqfn5/8+yu2VDhXuz1xYn+uP2ZMc8H3H1rrrcx9iBazDlxemHYlXDnfmhUUaQ8TpVTXGLo/nP8K5IVMqqjaAivesVpegiFbcYw6Fs5/tXlbjLWzoXx9yPmjoy/suMfJzuMUGpAby5iWPwLnA/nAa8AWoD9wAlAJvAocCbwnIscZY96LU6yqHeozekLo8iw126LWVanLGMOvXlzE0s3OMfBThvds3k+orhze+S3MfzLyTQbuDafcD3336IKIlVLtMuwA+MUC6/f2u7ehZAnU7gB3hrV2y5D9Yc+zre0xQoUPwB1zYivPcZDVulq70zpe9iqc8LeU+KcllqRlB1ACTDDGVDcWikge8B7WQN1J9uc32B9VFwlk9YSQ9zF37Y7kBaMS5oFZ3/PGos2OshF9cnnovL3J8rqtP3TPnQs7V7e82JMNh/8Opl4KLl28Wqlux5sNUy6xHu114t+tR3u4PXD9mphCS7ZYOq6vBP4amrAAGGOqgL8Clxlj/MCDQOoMSU4TJmyuf2b99ig1Var6YNkW/vrut46y/CwPj5y/D0U5GbDgGXjkyMgJy9AD4dLPYNoVmrAopVJOLC0tg4CGKOf8QGNn3GZANyDpYq585/LO2f6y5ASiEmJlaSW/eHYBoXP3rKnNezGihxdeuwq+frzlhRn51vb1e1+og2yVUikrlqTlW+AXIvKW3aICgIh4gF/Y58Ea56KjQLtYRn4fx3FeoCw5gai4K69t4JInvm4xU+j6Y8dwaL86eOws2DS/5YUD94EfPA5FQ7ooUqWUSoxYkpYbgZeAlSLyX6yBuP2AU4GBQONqOEcBczofouqIrCLnuhxFwQrdfygNBIKGK5+Zz+ptjl5ZTpk0gP8bsBr+9ZPmQXWh9r0Ejrkt+iwCpZRKIR1OWowxr4jIiViziH4OCNY8yq+Anxpj3rHr/SSegar2yevp3H8oUxqoqy4nK68oOQG1w8qVK7nqqquinp87d27Uc7uKO99ZzkffORsuxw8s4K+DPkGeupEWU5k92XDyPTDxrK4LUimlEqxDSYuIZACHAkuNMVNEJAfoAew0xtQkID7VQZH2HyrftrlbJy0bN27kH//4R7LD6LZeXbiJf320ylHWO9fLM8Pewvv+/S0v6DkSzp4J/XT7BqVUeunoiDw/8DqwO4AxpsYYs1ETlu6joKCIOuMc/1y1oyRJ0ajOWllaxa9fWuQoy3IHeWv48+R/HSFhGXMi/N//NGFRSqWlDiUtxpggsAHQfem7KZfbRZk4N9euKUv/pGXjxo3cfffdHH300QwZMoSMjAyKi4s544wz+OKLL5IdXkxqfQEuf2oeNb5AU1kmPt4f+Ah9VkZYwfKIG+HsJ5v3IVFKqTQTy9zHR4HLRUQXeeimKt3ONy1fRfeexHXIIYdgjIn6uOCCC9q8x7333svVV1/NqlWrOOqoo7j22ms58MADeeWVV5g2bRrPP/98F7yS+PrDq0v4dkvzNvIFVPNWz78zqHSWs6K44ZQH4KBrdcC1UiqtxTJ7yAeMBpaJyKtY67GEjgI0xph2LsunEqHGUwTN/5wTqOzeSUs8TJkyhY8//piDDjrIUf7JJ59wxBFHcOmll3LKKaeQmdn9l6kGeOnrDTz/1Yam4z6U8WzOXxhRE7ZgnCcLznwcxhzfxREqpVTXiyVp+XPI59dEOG8ATVqSqD6jJ4TuobcL7D90+umnRyw/6KCDOOyww3j33XdZvHgx++yzTxdH1nErtlTyu/8uaToeJKU8nXE7Q4JbnBUzC+HHz8LQaV0coVJKJUcsScvwuEeh4sqf1cvautLmrt21l/L3eq2ByR5PLD/uXavG5+fSp+ZR22A1lfVlJ097b2WIhLWW5fWDc1+C4glJiFIppZIjlnVa1iYiEBU/Jqe34zijftfdNHHdunW8//77FBcXM2FC936DN8bwu/8sYWVpFQBFVPJkxm0McYUlLD2Gw3n/gZ76/4NSatcS8yYkIjJGRH4qIjeISLFdNkBEsuMXXovn7CEiM0Wk3H7MFJGiNq4REblJRDaJSK2IzBKRcWF1ZomICXs8m6jXkWiufGfSkt0QYaXUXUBDQwPnnXce9fX13Hnnnbjd3Xvs+PNfrefl+RsByKWWGRl/ZpRro7NSv/Fw0TuasCildkkdTlpExC0ijwLfAP/EWhl3gH36X8Bv4hdeC08Dk4Bj7cckYGYb11yHNfbmCmBfoAR4T0Tyw+o9jLVfUuPjp/EKuqtlhG2amBcoT1IkyRMMBrnooov4+OOPueSSSzjvvPOSHVKrlpdUcOMr3wDWtOaHvX9jksu5oBw9R1otLPn9ItxBKaXSXywtLTcAPwZ+BYzHWsa/0VtYyUTcicge9r1/YoyZY4yZA1wCnCgio6NcI8BVwK3GmJeNMUuAC4Ac+zWEqjHGlIQ8UvadPqvImbQUBctwbAuc5owxXHLJJTz55JOce+65PPjgg8kOqVV1DQGufGY+9f4gbgLc672Xae6lzkoFA+H8VyCvb+SbKKXULiCWpGU6cIsx5i6ad3RutJrEDdTdHyg3xjStFGaM+RwoB6JNnxgOFAPvhlxTD3wU4ZpzRGSbiHwjIn+N0BLjICKZIlLQ+Ljr73e1Wr8r5YbtP5QlDdTXViQpmq4VDAa5+OKLeeyxx/jRj37EjBkzcLli7gXtEne+/S3fbalCCHKn9yGOdn/trJDTC877LxQNTkp8SinVXcTy13wg0XdvrgMS9eZdDJRGKC+1z0W7BqydqENtCbvmKeBHWPsq3YK1U/XLbcTzG6yEqRwov+/++ze0Ub/LFPQa0KKsfNvmJETS0k033dRi8bhZs2a1es2MGTNaXDNjxowW9YLBID/5yU94/PHHOfvss5k5c2a3H8fyyYqtPPbZasBwo2cmZ7g/cVbILIBzX4Y+o5ISn1JKdSexJC2lwIgo50ZjLfPfbvYg2fBBsOGPxsU1IvVxSJTyUOHnHdcYYx42xrxvjFlijHkWOBM4UkQmt3LP24HCxscVl18+qI0YukxRYQ/qjXNiWOX28LwtvTS2sDz++OP84Ac/4Mknn+z2CcvOah+/fGEhAD9zv8aFnnecFTxZ8KNnYcCkrg9OKaW6oVgWrngTuEFE3sYa1ApgRKQQuBJ4rYP3uw9oa6bOGmAiEGkEYh9atqQ0aoyvGGvl3kZ9W7kGYB7QgLUx5LxIFexupqYl3LZuidQIlBzW/kMF9KN5qnNtmu8/9Mc//pEZM2aQl5fHqFGj+NOf/tSizqmnnsqkSZO6PrgIjDH89j+L2VJRz6Gu+Vznec5ZweWBs2bCsAOSE6BSSnVDsSQtNwLHAUuB/2G1WNyGNSi3Aat7pd2MMduANpdsFZE5QKGITDHGzLXL9sNq6Zgd5bLVWInLUcB8+5oM4BDg+laebhzgxZnopJRKVyH9gs1Ji6+8+yRVibBmzRoAqqqquPXWWyPWGTZsWLdJWl78egNvLSlhhGziHu99uCS0MVDgtH/BqKOTFp9SSnVHHe4eMsZswZo6/AywN9YuN3tizRyaZoxJyEpmxphlwNvAwyIyVUSmYk1Tft0Y0zQgWESWi8hp9jUGuBv4rYicJiLjgRlADdb0aURkpIjcKCL7iMgwETkeeAEryfksEa+lK1R7eziO/VXpvf9QpHEv4Y/p06cnO0wA1m2v4aZXvyGfGh72/o0CqXVWOOJGmHBmcoJTSqluLKZ1ze3E5WdxjqU9zgHuoXk20KtY66+EGo3V+tLoTiAbeADoAXwBHG2MaVzo3gccAfwCyAPWA28ANxtjQrYdTC0t9h+qTv/9h1KBPxDk6ucXUOtr4GHv/Yx0hTXmjTsdDrw6OcEppVQ31/03Ywlht+Kc20YdCTs2wE32I1L99VjdRWnFn9VT9x/qhv4563u+XruTX3pe4Aj3fOfJ4glwyv0gEvlipZTaxcWUtIjIgViLsw3FasUIZYwxR3Q2MNU5uv9Q97NoQxl3f7CCE1yfc4XnFefJnF7ww6chIyc5wSmlVArocNIiIhcCjwI7gO9wdkKAc4VclSSuXN1/qDup9we49vmFjDar+UvGv5wnXR446wkoGpKc4JRSKkXE0tJyHfA8cIE97Vd1Q94C5+zwvEBZcgJRANz34Uq2lW7itcy7yJGwX5tj74BhByYnMKWUSiGxJC1DgZ9rwtK9ZRY696gpMLvGMv7d0ZKN5TwwayUPef/FIAkbED35fNj3J8kJTCmlUkwsK+IuI/Iib6obye/p/BblUI+vpjJKbZUoPn+QX724iHPknZYDbwdNgeP/qgNvlVKqnWJJWn4L/FpEBsY7GBU/Bb36tygr356ya+WlrAc/+p5AyTfc4HnaeSK3D5w9EzyZyQlMKaVSUCzdQ5djrYPynYgsAMLn0hpjzCmdDUx1TmGP3viMmwxpXmqmcnsJfQbrxntdZXlJBQ99+A0veu8jUxqcJ099EPKj7fOplFIqkliSlolYq+CWAgPsR6i2Ni9UXcDtdrFdCukbsv9QTVl6b5rYnfgDQX71wiKulacZ41rvPDn1Mtj9yOQEppRSKazDSYsxZlgC4lAJUOEqpG/o/kMV6b3/UHfy8Cer6b15FhdmhO3c3G88HPGH5ASllFIpLqVWxFUdU+MpsjYpsPkr03v/oe5iZWkVM9+fy6te53osxpOFnPEoeLOSFJlSSqW2dg3EFZHzRaRXWNkAEXGHlQ0UkT/GM0AVu/oM56aJVGvSkmiBoOG6F+ZzuzxAb3FOM5djboO+Y5IUmVJKpb72zh56HBjZeGAnK+uxdncONQi4IT6hqc5qyOrpOHbp/kMJ98ScNUza9ByHuBc5T4w+Afa5KDlBKaVUmmhv91CkhSR0cYluLpjTx3Gs+w8lVkl5Ha+9+x7PeJ5xlJu8YuTke3U9FqWU6qRY1mlRKcKVF7b/kK8sOYHsIm55dRE3mQfIFH9TmUGQ0/8Fub1auVIppVR7aNKSxjLynS0tubr/UMJ8uHwLg5Y/xkTXake5TPs5jDg0OUEppVSa0dlDaSyzyLmUf5EpS04gaa7G5+ehl99lhudFR3mg12jch/8uSVEppVT66UjScqiIDLI/d2EtIneYiAwLqaPLrXYjuT2cK67mUE9DXTXerNwkRZSe/vH+t1xddx9ZruZVbw2C+9T7dZl+pZSKo44kLXdEKPtLhDJdEbebiLj/0LbN9B60WxKiSU/LSyqom/0w+3mWO0/s9zMYvG9yglJKqTTV3qTlsIRGoRKiqEdvGowbb+j+QztKNGmJk2DQcNcLH3CX2zlbqKFgCN4jfp+kqJRSKn21K2kxxnyU6EBU/Hk8brZKPn0oayrT/Yfi59m56/hx6d/Jc9c5yr2n3gsZ2gWnlFLxprOH0lylq9BxXF+uSUs8bK2sZ8nbD3Goe6Gj3L/nuTpbSCmlEkSTljRX7XEu5e+v0KX84+Efr3zKr8wMR1l9dl88x96anICUUmoXoElLmqsL33+oZltyAkkjs1duY/9v76CHVDnKM0+5G7KLkhKTUkrtCjRpSXMNWc6VWHX/oc5pCAR5++XHOME911FeN/oUGHNCkqJSSqldgyYtac5kO5MWb53uP9QZT3+6nP+r/pejrN5bSNZJf0tSREoptevQpCXNuXLD9h9q2JmkSFLftqp66v73FwaJs4vNe/wdkNcnylVKKaXipcPL+IvI+a2cDgJlwDxjzKZYg1Lx4yno6zjO85clJ5A08Ngr7/EL86pjf/PKflPIn/Sj5AWllFK7kFj2HppB86q3IX++HWVBEZkJXGKM8aOSJits/6ECU56kSFLbgnU72W/5n8l0N/84B3CTf/o/QKSVK5VSSsVLLN1DU4A1wEPAocAeWCvmPgysBY4DbgR+BNwQjyBV7HJ7OJOWPGrx19ckKZrUFAwa3nrxEQ5xL3KU1+19CfQbm6SolFJq1xNLS8svgP8YY34ZUvYt8JGI/A24yBhztoj0AM4Bbo5DnCpG+b2KW5SVbyuh18ARSYgmNf1n7necX/5PR7tidUYfco/WHZyVUqorxdLSciLwdpRzbwHH2J9/CAyJJSgVP0U9++I3zm9z5Y6SJEWTeirqGih/53YGinOqeOYJt0NmfpKiUkqpXVMsSYsbGBnl3G40/z/qA+pjCUrFj9fjoVycb641OzVpaa+Zr73PucFXHWU7+07FM/HMJEWklFK7rliSlneBP4nIUaGFInIMcIt9HmAM1tgXlWTlriLHcV15aXICSTErSirYc/GfyAjZJduPmx4/uEcH3yqlVBLEkrT8Amta89siUiYi34pIGfAmsBO4KqTu3Z2MT8VBjafIcRyo0v2H2mKM4c3nH+RA1xJHefXkn0Kf0UmKSimldm0dHohrjNkoInsC04GDgV7AfOAj4N/GmBq73v1xjFN1Qp23h9VZZzOatLTpg0WrOGu7c/BtRUZfCo/RCXFKKZUsscwewk5MHrAfqptryOoJ1c3Huv9Q63z+IJtfv50jxbnlQeYJf4bMvCRFpZRSKqWW8ReRHiIyU0TK7cdMESlq45rTReQdEdkmIkZEJkWokyki99p1qkXkVREZlKjX0dWC4fsP1ev+Q615+cPZnOX7r6OstO80MieelpyAlFJKATEkLSLiFZHfichS+w0+EPZI5Aq4TwOTgGPtxyRgZhvX5AKfAb9upc7dwGnAD4EDgTzgdRFxdyrabkLC9sXJ8un+Q9HsrPZROPtWMqWhqSyAiz5n3qWDb5VSKsli6R66Hbgaa02W/9JF05pFZA+sRGWqMeYLu+wSYI6IjDbGfBvpOmPMTLvusCj3LQQuBs4zxrxvl50LrAeOBN6J80vpcp58Z9KSGyhLTiAp4OVXXuJiZjvKto89n75990hSREoppRrFkrScBfzRGNPVK93uD5Q3JiwAxpjPRaQcmIa1Km8s9ga8NE/VxhizSUSW2PeNmLSISCaQ2Xj86+uvz7/m6mtiDCGxsgqdS/kXBnX/oUhWbqlg3+V3Otofq1359D3xD8kLSimlVJNYxrT0AD6OdyDtUAxEWmCk1D7Xmfv6jDHhfSZb2rjvb4Dyxsd999+/oRMxJFRO2P5D+dQQ8NUlKZru66MX7mWia5WjrOGg6yCnZ5IiUkopFSqWpOVjrLEkcSEiN9kDZFt77GNXN5FuEaW806G1cd/bgcLGxxWXX95tB+4WRNp/aLuuihtq9tK1nLj1IUfZ9qxhFB18aZIiUkopFS6W7qErgVdEZC3wujHG19YFbbgPeLaNOmuAiUC/COf6YLWKxKoEyBCRHmGtLX0hbHBDCGNMPSHjebZu6b6rzBb17EfACG5pzsEqd5TQs/+w5AXVjQSChjWv3Mo0KXOU55z8Z3B7kxOUUkqpFmJJWhZgjQF5ATAiUhN23hhjCtt7M2PMNmBbW/VEZA5QKCJTjDFz7bL9sFo6oiYX7fA10AAcBTxv37c/MB64rhP37TYyMrzskHx6UtFUVq2bJjZ5/ZMvOL3uZcdCcpv6HMCAsccmLyillFItxJK0vERiumNaZYxZJiJvAw+LyE/t4oewWnuaBuGKyHLgN8aY/9jHPbF2mx5gVxkt1tTVEmNMiTGmXEQeBf4mItuBHcBfgcXA+13x2rpChauQnsHmpKVe9x8CoKreT/asm8kKm+Lc78y/JTEqpZRSkcSyjP/0BMTRXucA99A80+dV4IqwOqOxWl8anQw8HnLc2BV1M3CT/fnVgB+rpSUb+ACYbowJkCaq3EUQXN907K/UpAXglVdf4hzjbKjbMuocBvTTKc5KKdXdxLSMf7IYY3YA57ZRR8KOZwAz2rimDvi5/UhLdRk9rE4wW7C6zR65tLdxZzUTl9zh6BaqduUz4NQ/Ji8opZRSUbUraRGRg4F5xpgq+/NWGWOSMSVataIhs5fuPxRm1gv3c444pzjXTPsVuTrFWSmluqX2trTMAqYCc+3Po41paZwmnBbL36eTYE4va7SOzVu3a+8/9M3aLRy68V+OVpatmUPoc9hlyQtKKaVUq9qbtBwGLA35XKUYye3tOM727dpJy+KX7+SH4uwiyzz+Vp3irJRS3Vi7khZjzEeRPlepw5Pf13Fc5N+apEiSb/bi7zi+7GlHK8vGwr0ZOPGk5AWllFKqTbGsiKtSUN7A0Y7jfmYrvpqKKLXTVyBoKH39FgrClhfqdfqduouzUkp1c+0diPtYB+5pjDEXxxiPSpCBIyYSNIIrZFXcku8XMWTCgUmMquu99+lsTqh7w9HKsnbACQwduk/0i5RSSnUL7R3TcjjOwbdFWGuh+IHtQC/7XuVA+MaDqhsoLCxgvfRjMM0r4e5cu3iXSlrqGgJkzLoFrzQvv+PDw6Azb0tiVEoppdqrXd1DxphhxpjhxpjhwFlAFdZCb9nGmP5YC7KdizWp9oeJClZ1TmnWMMdxQ8ny5ASSJG+99QqHB+c4ykrGXIC757DkBKSUUqpDYhnT8jfgr8aYZxpXjDXGBIwxT9vn/h7PAFX81BTs5jjOKFuRpEi63s6qeobNu8NRVin5DD7590mKSCmlVEfFkrTsDSyJcm4xMCnmaFRCufo5B+P2qlkVpWb6ee/lR9iLbx1llftdjeT0SFJESimlOiqWpKUCODLKuSPt86obKhg83nFcHCgh4KtNUjRdZ11pGVO+v8dRts07gAFHpu2uDUoplZZiSVpmAr8Skb+IyF4i0t/++Ffgl/Z51Q0N2G1Px7FbDCWrojWapY8vXvwbw6TEWXjEH8CTkZyAlFJKxSSWDRN/C/QFrgWuCSkX4En7vOqGevXsxWZ605/mlWC3r1nEwDH7JjGqxFqyah1HbHnMMcV5Q+44Bu13dvKCUkopFZMOJy3GGD8wXURux1rSvxfWtOdZxphdazpKCtqSMZT+vuakpX7z0lZqpzZjDN+//CfGS5WjvPCUP+tCckoplYJiaWkBwBjzLYSNbFTdXlXBbrDt66bjjB3pO4Po84WLOabyZUcry6rehzFi1EHJC0oppVTMYk5aRKQvMBRrjRYHY8zHnQlKJVCf0YT0DtGjZnXyYkmgYNBQ8eYfyZKGpjI/LgaeeUcrVymllOrOOpy0iEh/rMG2jbs9N/4fa+zPDeCOS3Qq7vIHj4dlzcf9/RsJNvhwedNrUOqHn8ziyPr3Ha0sa4aeyW7FY5IWk1JKqc6JpaXlPmAv4HpgEVAf14hUQhWPdM4g8kqAknVLKR45KTkBJUC9P0D2R7fgDtlnqZYshp9xSxKjUkop1VmxJC2HAL80xjwe72BU4vXt24+tpog+UtZUtm31orRKWt5/8yVOCH7tKNs87ieMKChOUkRKKaXiIZZ1WgywPt6BqK4hImzOGOooq9u0LErt1FNR62No2HL9Za4eDD/p+iRFpJRSKl5iSVpeAE6MdyCq61Tmj3Qcu7d/l6RI4u9/Lz/EeL53lFXsdw2SVZCkiJRSSsVLLN1DzwMPi4gLeA1rjRYHY8y8zgamEsf0Hg07mo+LqtNjD6KSHRXs9d0/HINvS7yDGHLkpckLSimlVNzEkrR8aH+8Arg87JzOHkoBuQPHQUjjygD/ekzAj7hjngHfLcx98a+cLKXOwiP+AG5vcgJSSikVV7G8S10Y9yhUl+o7ck/4X/NxJg1s27iS3kNSdzrw9xs2ccBG53L9a3PGM3S/HyQvKKWUUnEVyzL+/05EIKrr9O8/iJ0mnx5S2VRWumphSicty1/8EyeEvB6AgpNu0+X6lVIqjcQyEFelOJfbxUbvEEdZzYZvkhRN5y34ZimH73zBUbay58H02OOQJEWklFIqEdrV0iIiNwKPGGM22Z+3xhhjdBWvbq4ibwSUNScqru2puY2UMYbtr91Itviayvy46H/Gn5MYlVJKqURob/fQTcDbwCb789YYQJOWbi7QaxSUNR8XVKXmDKJPP53FYbXO5fq/H3g6oweOTV5QSimlEqJdSYsxxhXpc5W6sgeOI3Q5kwEN68CYlBoD4vMHyfrfH3CFLNdfQxbDf/CnJEallFIqUTQB2UX1HeHcgyiHOnZuTq3Wlv+98Qz7Bhc6ytaP/SkZRf2TFJFSSqlE6nTSIiIFIrKPiAxtu7bqLgYMGk6VyXaUbfl+YZTa3U9FTR0j5t/uKNvu6sWoU3W5fqWUSlftTlpE5AQReVJEHheRI+2yy7DGuXwBrBKRF0REV/JKAR6Pm/Ue5wyiqo1LkxRNx3364j3sHrYFVsX+1yMZuUmKSCmlVKK1d/bQCVhL9tcDPuBcEbkCuBd4F1gE7AWcDlwK3JOQaFVcleUOh4rmWUOydXkSo2m/zVu3sc/3DzgG327IGM7wI36SvKCUUkolXHtbWn4JfAD0MMYUAg8AfwceMsYcZ4y53hhzNDADuCAhkaq48/cc7TjOr/w+Ss3uZeHzt9JXdjrKXMfcCi7dPUIppdJZe5OWccB9xpg6+/gvQBbwn7B6LwHD4xSbSrCsAXs4jot9a60ZRN3YtytXclDpU46yFfn7MWDvE5IUkVJKqa7S3qSlN1ASctz4+bawetuBws4GpbpG7+ETHccFVFOxbWOSommfjf/5PblS33QcNELv0+9IYkRKKaW6SkdmD5kon3cZEekhIjNFpNx+zBSRojauOV1E3hGRbSJiRGRShDqz7HOhj2cT9Tq6i4HDRlNrMhxlJSu77wyir+Z+xiFVbznKlhWfRI/hk5MUkVJKqa7UkQ0TR4uI3/68cfDAGHEuRpboHfeeBgYBx9rHDwEzgZNauSYX+Ax4AXi4lXoPA6FbFNTGHmZqyPB6+M49mFHB5rEsVRuWAN2vq8UfCBJ890bcIQvJ1ZLJiLNuS2JUSimlulJHkpYZEcpmhh0LCWqFEZE9sJKVqcaYL+yyS4A5IjLaGBNx8xxjzEy77rA2nqLGGFPSRp20szNnOFQ1Jy2mtHvOIPr4jac43P+Vo2zVbtMZ12twkiJSSinV1dqbtFyY0CjaZ3+gvDFhATDGfC4i5cA0oLM7/p0jIucCW4C3gJuNMZWdvGe35+u5O1S933ScU9H9ZhCVVVSy2zzn0vw7pIgxZ/wuSREppZRKhvbuPfTvRAfSDsVAaYTyUvtcZzwFrMYaYDweuB3YEzgq2gUikglkNh7/+vrr86+5+ppOhtH1Mor3gHXNx/3q1yQtlmjmPnsrR+NsBNux/w30zC5IUkRKKaWSIel7D4nITREGwYY/9rGrR+p66nSXlDHmYWPM+8aYJcaYZ4EzgSNFpLURnr8Byhsf991//4bOxJAsPYc5ZxD1pJyasi1JiqallSu/5YCNjznKVmWNZbcjdSE5pZTa1SQ9aQHuA/Zo47EEqxWkX4Tr+2B16cTTPKAB2L2VOrdjTe8uBAqvuPzyQXGOoUsMHjmWeuNscNvcTWYQGWMoeenXLaY45576N3B1hx9dpZRSXakjA3ETwhizjZbrvbQgInOAQhGZYoyZa5fth5U0zI5zWOMAL7A5WgVjTD3WtgYAbN0Sqeeq+8vKzGSVawAjTHMfUfm6JbDP0UmMyvLFrDc4sPZDR9k3xSczYcy0JEWklFIqmVLm31VjzDLgbeBhEZkqIlOxpim/HjpzSESWi8hpIcc97bVZxtpFo0VkkogU2+dHisiN9k7Vw0TkeKzp0fOxpkqnve05zkWMg91gBlFdvY+eHzsH2laSw24/vDNJESmllEq2lElabOcAi7E2aWzcqPG8sDqjca7KezJWAvKGffysffwz+9gHHAG8gzUD6R773kcaYwLxfwndT12Rsxesx44FyQkkxOzn/8Yos9pRtnbiVWT36OyYa6WUUqkq6d1DHWGM2QGc20YdCTueQeQ1ZhrPrwcOiUN4KcszdD/Y+EjT8Ujft+xcv4weg/do5arE2VSyickr73Ps4rzeM5RxJ1+dlHiUUkp1D6nW0qISYMKBJ7HdOKcPr/3fjOQEA6x85tcUSZWjLHDMnxFPRpQrlFJK7Qo0aVHk5WSzuIdzSZp+a15Jyo7Pi7/6lAPKXnWWFR7KsH2P6/JYlFJKdS+atCgA8vb9seO4f3AzG5d81KUx1Pt8uN66Nmx/oQwGn31Xl8ahlFKqe9KkRQEwab/DWMMAR9nWz8K3lkqs2U/ezLiAc+bS8pEXUzRgZJfGoZRSqnvSpEUB4PG4WdXfubvzsJJ3MP76KFfE1/JFc5m29kFH2SZXfyb84Pdd8vxKKaW6P01aVJMBB53vOC6ikpWz/5vw562rq0P++zMyxd9UFjRCw0n348nKTfjzK6WUSg2atKgmo/eYwGK3c5pz3dfPJPx55878PaODzt2lFw4+h6F7HZHw51ZKKZU6NGlRTUSEbcNPcZSNKv+U+qqdCXvOpfM+Zf8NjzrK1rsHMeFcXflWKaWUkyYtymHUYefhM+6m40waWDHryYQ8V11tDZmvX45Xmhce9hsX5pR/areQUkqpFjRpUQ4DBw5iQdZ+jjLvkhcT8lxfP/FrRgbXOMoWDL2QIRMPTsjzKaWUSm2atKgW6see6TgeXbeAipLVUWrHZtmXHzJ10xOOstXu4Uw697a4Po9SSqn0oUmLamHiYWdRYXIcZd9/OCNu96+triLnzZ87FpFrMG7cZzyIJyMrbs+jlFIqvWjSolooLMhnYeFhjrJeq/4Tl2X9g4EAix+6mKFmg6N8wYj/Y8jYqZ2+v1JKqfSlSYuKKHPyjxzHQ/xrKfnuy07d0wSDzP3XZUwpf9tRvsKzO5N//MdO3VsppVT606RFRTTpgOPYRG9H2eZP/t2pe37x798wtfRZR1md8ZL5g4dwe3UHZ6WUUq3TpEVFlOH18F1f587Kgze+gQn4o1zRurnP3c7UsGX6/cbFdwffy5DRk2OOUyml1K5DkxYVVZ8DznMc9zY7Wf3pcx2+z/zXHmDKsjscZUEjLNzndiYe8aMoVymllFJOmrSoqMZOnMK3MsJRNvh/P+f7t+5p96DcJR8+zYSvbmhRPneP37D3ST+LS5xKKaV2DZq0qKhEhI0jnGu2eAkw8ovfs/yRizANda1e/92c1xn10c/xSNBR/tnQy5j6w+vjHq9SSqn0pkmLatV+Z/6SOZ4pLcrHbHyZNXcdTu2Ojc4TwSAb5r3Dwvt+zJC3p5MhzjEwn/X9EdMuuDWRISullEpTYuKw9oaCrVtKC4DyXn1643KlVy5YUVPHJ//6BSeUP9vi3Hbpif/MJ8gp6s2aDx+j7+r/0i9YGvE+swtPYOqVT+Jyp9fXRymlVOcEg0G2b90GUNinX9+KaPU0aYmTdE5aAIJBw1vPPcBhy28mR+od5/y48BCMcqVlbs4h7HX1S3i93kSGqZRSKgW1N2lJv3dXlRAul3DCjy5n4THPs8H0cZxrLWEJGOGTwpOZ8PPnNGFRSinVKZ5kB6BSy/7TDmXVwA+Z98Q5TA4silpvBUNYP+RUhh82nYOGj+zCCJVSSqUrTVpUh40YOoTKa9/lvX9dzlHlLzWVbzcFLCg6ivyp5zF534PZ3eNOYpRKKaXSjY5piZN0H9MSSTBo+OSDV6ld8QmeAROYfPiZ9CzITXZYSimlUowOxO1iu2LSopRSSsWDDsRVSimlVFrRpEUppZRSKUGTFqWUUkqlBE1alFJKKZUSNGlRSimlVErQpEUppZRSKUGTFqWUUkqlBE1alFJKKZUSNGlRSimlVErQpEUppZRSKSGlkhYR6SEiM0Wk3H7MFJGiVup7ReTPIrJYRKpFZJOIPCEiA8LqZYrIvSKyza73qogMSvgLUkoppVS7pVTSAjwNTAKOtR+TgJmt1M8BJgO32B9PB0YBr4bVuxs4DfghcCCQB7wuIrpNsVJKKdVNpMyGiSKyB7AUmGqM+cIumwrMAcYYY75t5332BeYCQ40x60SkENgKnGeMec6uMwBYDxxvjHmnPffVDROVUkqp2LR3w0RP14XUafsD5Y0JC4Ax5nMRKQemAe1KWoBCwABl9vHegBd4N+S+m0RkiX3fiEmLiGQCmY3Hv7vhhvwrf34lwWCw3S9IKaWUUrT7vTOVkpZioDRCeal9rk0ikgXcATxtjGnM5IoBnzFmZ1j1LW3c9zfAHxoPHnn0Ua78+ZXs3L6jPaEopZRSqqV8oPu2tIjITYS8+Uexr/0xUl+WRCkPfx4v8CzWOJ7L2hNaG/e9Hbir8cDr9bJkyZKe48eP16wlAe76+135991//4YrLr980DVXX1OZ7HjSkX6NE0+/xomnX+PES+DXOB/Y1FqFpI9pEZHeQO82qq0BfgzcZYwpCru+DLjaGPN4K8/hBZ4HRgCHG2O2h5w7HPgA6Bna2iIiC4H/GmPaSqhUFxCRAqAcKAxpJVNxpF/jxNOvceLp1zjxkvk1TnpLizFmG7CtrXoiMgcoFJEpxpi5dtl+WGNUZrdyXWPCsjtwWGjCYvsaaACOsushIv2B8cB1HX5BSimllEqIlJnmYoxZBrwNPCwiU+2ZQw8Dr4fOHBKR5SJymv25B3gR2Ac4B3CLSLH9yLDvWw48CvxNRI4Qkb2AJ4HFwPtd+BKVUkop1Yqkt7R00DnAPTTP9HkVuCKszmis1heAQcDJ9ucLwuodBsyyP78a8GO1tGRjdRdNN8YE4hS36rx64Gb7o0oM/Ronnn6NE0+/xomXtK9x0se0KKWUUkq1R8p0DymllFJq16ZJi1JKKaVSgiYtSimllEoJmrQopZRSKiVo0qK6LRH5jYh8KSKVIlIqIv8VkdHJjiud2V9zIyJ3JzuWdCIiA0XkSRHZLiI1IrJARPZOdlzpQkQ8IvInEVktIrUiskpEbhQRfY+LkYgcLCKvicgm+2/CqWHnRURuss/XisgsERmX6Lj0G6q6s0OA+4GpWIv/eYB3RSQ3qVGlKXsH9P8DFiU7lnQiIj2Az7AWsTwOGAtcS/Omrarzrgd+hrUExh5YC4P+Cvh5MoNKcbnAQlouK9LoOuAa+/y+QAnwnojkJzIonfKsUoaI9MHaIPMQY8zHyY4nnYhIHjAPa1+u3wELjDFXJTWoNCEidwAHGGMOSnYs6UpEXge2GGMuDil7CagxxpyXvMjSg4gY4DRjzH/tY8HaI+huY8yf7bJMrI2GrzfG/CtRsWhLi0oljYsG6qaU8Xc/8IYxRleBjr+Tga9E5AW7m3O+iFyS7KDSzKfAESIyCkBE9gQOBN5MalTpazhQTPNCrxhj6oGPgGmJfOJUWxFX7aLszP4u4FNjzJJkx5NOROSHwGSad1NX8TUCuBTr5/c2YApwj4jUG2OeSGpk6ePPWP/ULBeRAOAGbjDGPJPcsNJWsf1xS1j5FmBoIp9YkxaVKu4DJmL996TiREQGA/8AjjbG1CU7njTlAr4yxvzWPp5vD1i8FNCkJT7OBs4Ffgx8A0wC7haRTcaYfyczsDQXPr5EIpTFlSYtqtsTkXuxmtgPNsZsSHY8aWZvoC/wtdWYBVj/pR4sIlcAmboHV6dtBpaGlS0DzkhCLOnqL8Adxphn7ePFIjIU+A2gSUv8ldgfi7F+vhv1pWXrS1zpmBbVbdlT6u4DTgcON8asTnZMaegDYALWf6aNj6+Ap4BJmrDExWdYG7mGGgWsTUIs6SoHCIaVBdD3uERZjZW4HNVYICIZWDM+ZyfyibWlRXVn92M1954CVIpIYz9quTGmNnlhpQ9jTCXgGCMkItXAdh07FDd/B2aLyG+xdpKfgjW1/P+SGlV6eQ24QUTWYXUP7YU1HfexpEaVwuwZhbuFFA0XkUnADmPMOnstp9+KyApgBfBboAZ4OqFx6ZRn1V3Z0+wiudAYM6MrY9mViMgsdMpzXInIicDtwO5Y/6XeZYx5OLlRpQ97bZBbgNOwuig2Ac8AfzTG+JIZW6oSkUOB/0U49W9jzHR7csQfgJ8CPYAvgMsT/c+OJi1KKaWUSgna36eUUkqplKBJi1JKKaVSgiYtSimllEoJmrQopZRSKiVo0qKUUkqplKBJi1JKKaVSgiYtSimllEoJmrQo1UEisp+I/EdE1olIvYhsEZE5IvK3ZMcWTyJyvIjcFOXcGhGZ0bURWQteiYgJefhEZKuIfCYit9r7zYRfM92uO6yDz/VbETk1XrF3ByJSJCLb7J29G8tmiEhVnJ9njYi83sl7XCwiG0UkN15xqdSnSYtSHSAiJ2DtrVEAXAccDfwCa3+Zs5MYWiIcj7XiZSSnYa1Amiy/BfYHDgMuBmYBFwHLROScsLpv2HU30zG/BU7tVJTdzx+wVot9LtmBtMO/gWqs3zOlAN17SKmOug5rGfZjjDH+kPJnRaRb/3EVkRxjTE087mWMmR+P+3TCCmPM5yHHr9otXe8DM0RkkTFmMYAxZiuwNRlBdici0hNryfWrTTdeCl1EsoE6Y4xfRP4F/F5E/hyvn12V2rSlRamO6QVsC0tYADDGhO8yi4icbXcdVYtIlYi8IyJ7hdWZYZ8bJyIf2HW3ish9IpITVvdyEflYRErteotF5DoR8YbVmyUiS0TkYBGZLSI12JvH2TG9KyKbRaRWRJaJyB2hzfB218/l9ueh3THD7LIW3UMiMkREnrRjq7fve62IuELqDLPv80sRuUZEVtuvfY6ITG3ftyAyY8wOrDdlD3B1yHO26B4Skb1E5PWQWDeJyBsiMqjxNQO5wAUhr32Wfa6PiDwgIkvt2EtF5EMROSjs69Gh1ypWt+NrIrJdROpE5HuxNqULrbO7iDwd9jW+vJ1foun21yZiK0s7f/6yROR2+7X47O6b+0WkKMo9jxWRefbP2XIRuSjsfOP35mgReUxEtmJtupdpV3kKq1XzhyiFJi1KddQcYD8Rucd+k/FGqyjWrr7PAEuBs4DzgHzgExEZG1bdC7wJfIDVJXEf1htw+BvMSKxdVM8DTgQeBX4F/CtCCP2BJ+36xwMP2OW72891MXAscLcd32sh194CvGh/vn/II2IXi4j0weo2Oxr4PXAyVqvHX+3XEu5yrG3trwLOwUoQ3hSRwkj3by9jzJd2jAdHq2MnZ+8B/cLiWIf1/QHrtdZifZ0aX/tl9rme9sebgROAC4FVwCyxNpkL1+ZrFZFjgE+AIVi7Ex8H/MmOsbHOWOBLYDxwLdb3/w3gHhGJ1o0X6gRgvjGmLMK5Nn/+RESA/wK/BGba97sLuAD4UEQycdoT+BvWLtenAIuAR0Uk0vfmMaAB6+f6TPtzjDElwHL7uZQCY4w+9KGPdj6wWlo+AYz98GGNZ/k1kBdSbzDWH957wq7Pw3pTfS6kbIZ9ryvD6v7WLj8gSiwurP+czwP8QI+Qc7Psaw9v4/WIfY+D7foTQ87dZ/2JiHjdGmBGyPHt9vVTwuo9AASBUfbxMLveIsAdUm9fu/yHbcR7qF3vzFbqfA7UhBxPt68ZZh/vbR+f0sZzVYW+xlbque2v4fvAyyHl7X6twEr7kdXK87wNrAcKwsrvxUqwerQRZzXwzwjl7fr5A46xj38VVu8su/ySsJ+PWmBISFkWsB14MML35t+txP0kUNKR31N9pO9DW1qU6gBjzHZjzEFYbzy/Bl4BRmG9aS8Wkd521WOw3sieEBFP4wOoAz7CevMN91TY8dP2x8MaC+xujVdFZDsQwEqMnsB64xwVdv1OY8yH4U8iIiPsLoaSkHt8ZJ/eo62vQRSHA0uNMXPDymdgJUaHh5W/YYwJhBwvsj+2mP0TA2nj/EpgJ/BnEflZhFavtp/Aum6eiNRhJYwNwBFE/vq1+lpFZBRWC9qjxpi6KM+XZd//P0BN2M/Um1gJQdTuNbv7JgcobeVltfXz1/g9nBFW7wWshOiIsPIFxph1jQf2a/uOyN/jl1qJqxToa79WtYvTpEWpGBhjvjLG/NkY8wNgAFYT+DCaZzo0Nut/ifWGFvo4G+iNk98Ysz2srMT+2AusMSNYrTwDsWYsNSZPjWMassOub9GVIyJ59j32A36HlTztC5we5R7t1SvS82HNVGl6DSEcr9UYU9/J5w81JOR5WzDGlAOHAAuA24Bv7DEtN7fW3ddIRK4B/gl8AZyBlSzsi9USEin+tl5rH/vjhlaethdWEvxzWv48vWnXCf+ZCtX4XBGTItrx82d/9BtrYHMTY4yx67b6PbbVE/lr1NrMrjqsRDSrlTpqF6GZq1KdZIxpEJGbsQZ/jreLt9kfzwTWtuM2HhHpFfbGUWx/bCw7FWs8xOnGmKZ7isikaKFFKDscK8k61BjT2LrS+J94Z2zHGkMTboD9cVuEc3EnIlOwvm6PtlbPWDOLfmiP05iI1U1xI1aXxh1tPM25wCxjzKVhz50fpX5bGpOAQa3U2YnVKjYTuD9KndWtXN/4M9Qzyvn2/Pxtt+v1CU1c7K9hMVaCHqvWZjP1BOqNMXFdS0alJm1pUaoDRCTSGzM0dws0/of/Dla3wUi7VabFI8I9wtcX+bH9cZb9sfEPe+N/6o1vGJd04CW0uIftpxHq1tvP0Z7Wjw+AsSIyOaz8fPs5/9eBGGMi1pTeB7FaH/7enmuMZaEx5mqgDAiNP1qrgCHs6yciE7EG63aYMeY74HvgogiDWRvr1GB9DfcCFkX5mYrUstF4vQ9rsPDIVkJp6+fvA/vjuWH1zsBKpj8gMUZgDWZXSltalOqgd0RkA9ZMm+VYif8krNkcVcA/AIwxa0TkRuBWERmB1XWwE6vbaApQbYwJnfHhA661u2++BKZhdd+8ZYz51K7znl3vGRG5E6u5/FKgRwfin23H8aDdOtSA9Wa1Z4S6i+2P14vIW1j/6S+y3wDD/R0rQXnDft1rsWZ8XIY1+PO7DsTYHrvb04ZdWN0S+2HNhioAzjfGfBPtQhE50Y7rv1hv5ILVPVaE9TVutBg4VEROwuq+qDTGfAu8jrV2yM1YY4FGY7XSrCb2v6mXY/1MfS4if8eayTQEaz2gxmTiF8CnWLPP/ok12DUf2A04yRgTPm4o3CysWUmRtPfn7x2ssUAFWAPQJ2LNopqP1QoUV2JNl59CGy1naheS7JHA+tBHKj2wZko8hTWgsBLrj/1arMGwe0SofwrwIVCO1Te/Bmvg4hEhdWZgJTwTsP6brsFqin8AyA2734lYYzFqscZA3Ik1bdlgdfk01psFLInyGvbHSl6qsQY5Poz1H7wBpofUy7DPlWLNAAqdgbOGsJk1WG+yT2F1BfmwkrpfAq6QOsPs+/wyQlwGuKmNr/+hNM/cMlhJ1zb79dwKDI1wzfSw2EdjDTJdaX+ty7DGp1wQdt2eWElCtX39rJCvy1/sr38t8LX9fZ4BrIn1tWKNjXnTjqfOju+usDrDsN7AN9hf41Ks5OGGdvzsHm4/775h5TNo/89fFlb32Rr7+TfZ9YrC6q0BXo8Qw6zGr2PY92afNmKenOzffX10j4cY01pXolIq0cRapO1MY0xesmNR6U1EFgGfmbDxON2ViMwERhhjDkh2LKp70DEtSim167gOmC72yr/dmYiMxJppd32yY1HdhyYtSim1izDGvI21gvLwZMfSDkOAK0zzmBqltHtIKaWUUqlBW1qUUkoplRI0aVFKKaVUStCkRSmllFIpQZMWpZRSSqUETVqUUkoplRI0aVFKKaVUStCkRSmllFIpQZMWpZRSSqUETVqUUkoplRL+HxxFSAZFay7DAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "h2_distance = np.load('./h2_distance.npy')\n", "h2_lda = np.load('./h2_lda.npy')\n", "h2_overlap = np.load('./h2_overlap.npy')\n", "\n", "h_energy = -0.47790086199331666\n", "\n", "fig, ax = plt.subplots(1,1, dpi=100)\n", "\n", "ax.axhline(y=0, alpha=0.5, c=\"grey\", ls=\":\")\n", "ax.plot(h2_distance, h2_lda - 2*h_energy, lw=3)\n", "ax.plot(h2_distance, h2_overlap - 2*h_energy, lw=3)\n", "\n", "ax.text(1.6,-0.05, \"H$_2$\", color=\"black\", fontsize=20, fontweight=\"semibold\")\n", "ax.text(9,-0.03, \"OA\", color=\"tab:orange\", fontsize=15, fontweight=\"semibold\")\n", "ax.text(9, 0.04, \"LDA\", color=\"tab:blue\", fontsize=15, fontweight=\"semibold\")\n", "ax.set_ylim(-0.15,0.05)\n", "\n", "ax.set_xlabel(\"Separation Distance (bohr)\", fontsize=12)\n", "ax.set_ylabel(\"Binding Energy (hartree)\", fontsize=12)\n", "\n", "ax.set_ylim(-0.2,0.12)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.5" } }, "nbformat": 4, "nbformat_minor": 4 }