Skip to content

Instantly share code, notes, and snippets.

@rnelsonchem
Last active December 23, 2015 17:59
Show Gist options
  • Save rnelsonchem/6672179 to your computer and use it in GitHub Desktop.
Save rnelsonchem/6672179 to your computer and use it in GitHub Desktop.
Aufbau for Cr ion
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Ignore This Stuff"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline\n",
"plt.rc('font', size=16)\n",
"\n",
"npoints = 100 # The number of points in our simulation."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 1-Electron Orbital Energies \n",
"\n",
"Here's where we will set 1-electron energies for our *d* and *s* orbitals. Remember that the *d* orbitals are always lowest in energy, and they drop in energy faster than the *s* orbitals.\n",
"\n",
"(The *linspace* command makes a linear spaced array of numbers that starts from the value on the left and goes to the value on the right. *npoints* is setting the number of points in the array.) "
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"d = np.linspace(-10, -130, npoints)\n",
"s = np.linspace(-5, -50, npoints)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAosAAAF3CAYAAADNUgHoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8VfX9x/HX9yaBJIWEYCCEJcQgeyQoWnCEoVYQxToK\nDsCqKHVRqNqCTIltHUX6K4r6U4baKu5drSjkJ2JrM0CGjCIgchMII4CQ/f39cZI0NwMC5ObcJO/n\n45FHzLnnnnwYwpvvOd/Px1hrERERERGpisftAkREREQkcCksioiIiEi1FBZFREREpFoKiyIiIiJS\nLYVFEREREamWwqKIiIiIVKtRhkVjTAdjzOvGmIPGmBxjzBvGmA5u1yUiIiISaExj67NojAkH1gDH\ngIdKDs8FwoE+1tqjbtUmIiIiEmiC3S7ABbcDnYGzrbXbAIwxa4EtwB3APBdrExEREQkojXFlcTnQ\nxFp7YYXjKwCstUkulCUiIiISkBrjM4s9gXVVHN8A9KjjWkREREQCWmMMi1HAgSqO7y95TURERERK\nNMawKCIiIiI11Bg3uByg6hXEljiriz6MMY3roU4RERGp16y1pjav1xjD4nqgVxXHe+A8t1hJY9sE\n1JDMmjWLWbNmuV2GnAL92tVv+vWr3/TrV38ZU6s5EWict6HfBc43xnQuPWCM6QQMLHlNREREREo0\nxrD4HLAdeMcYc6Ux5krgHWAn8IybhYmIiIgEmkYXFksmtAwBNgMvAi8B/wGGaHpLw5OUlOR2CXKK\n9GtXv+nXr37Tr5+U1+iacp8sY4zVz5GIiIjUB8aYWt/g0uhWFkVERESk5hQWRURERKRaCosiIiIi\nUi2FRRERERGpVmNsyi0iIuIX/miILFKeG5tuFRZr4KMtH5EYm0hMsxi3SxERkQCnDhriL279Y0St\nc07AGGOHLBlCujedsJAwEmMTSWyT6HyOTaR9RHv9S1JERICytiVulyENVE1+f/mjdY7C4gmU9lm0\n1rIjZwdp3jSfj8LiwrLgWPoRFxWHx+hxUBGRxkZhUfxJYTFAnagpt/ewtyw4pnpTSc9M52DuQRLa\nJJAYm1j2uWt0V4I9uusvItKQKSyKPyksBqhTmeCy7+g+0jPT/xsgvensPryb3jG9SWyTSEKsEyB7\ntupJ0+CmfqpcRETqmsKi+JPCYoCqrXF/h/IOkZGZ4bMK+d2B7+gW3c3nFnafmD6Eh4TXQuUiIlLX\nFBbFnxQWA5Q/Z0MfLTjK2qy1pHv/uwr5bfa3xEXFlYXHhDYJ9GvTj8jQSL/UICIitUdhUfxJYTFA\n+TMsViW/KJ/1e9aT5k0ru5W9Nmstsc1jfXZiJ8QmEB0eXWd1iYjIiSksVm3x4sX88pe/ZNGiRYwb\nN87VWpKSkkhJSaG4uLhG59dF7ePHj2fp0qVs376djh07VnueW2FROy4CTJOgJiTEJpAQm1B2rKi4\niE37NpXdwn7ki0dI86YR2TSy0k7s2GaxauUjIiIB6WT/ftq0aRPz5s1j+fLl7N69m+DgYM466yxG\njhzJpEmTiIqK8nsdxpiyj1Lbt28nLi6OcePGsWjRolOq4XRqqmsKi/VAkCeIHq160KNVD27qcxMA\nxbaY7w58VxYg//zPP5PmTSPYE1wpQJ4ZeWZA/yYUERGp6Nlnn+Wuu+4C4JJLLuHaa68lPz+fFStW\n8PDDD7Nw4ULeffddzjvvPL/WcfXVV/PTn/6UNm3alB0r/Tu1Nv9uDeQVaYXFespjPJzV8izOankW\n1/W8DnB+o+06tKssQC7KWMQ9H93DsYJjzg7scs3Eu5zRRb0gRUQkIL3zzjvceeedxMTE8P7779O/\nf3+f1xctWsSECRMYPnw4GRkZdOjQwW+1REREEBER4XOsNNgFcsCrTUoLDYgxhg6RHbiq21XMHjyb\n98a8xw+Tf2DjXRu5f+D9RIVF8cbGN7j85cuJ/EMkFy66kHs/upfFGYtZm7WWgqICt38IIiJSj/34\n449MmTKFtm3bEh4eTv/+/XnzzTdP6hqFhYXcd999GGN45ZVXKgVFgFtuuYWpU6dy4MABpk6d6vNa\nUlISHo+HvLw8pk6dSufOnWnSpAnz58/3OS83N5cpU6bQvn17wsLC6NevHy+99FKl77V48WI8Hg9L\nly4t+zouLg6AJUuW4PF4yj5SUlIA2L17NzNmzGDAgAG0atWK0NBQunTpwv3338+RI0dO6ucjEGhl\nsQaOHYOwMLerOHUxzWL4WfzP+Fn8z8qOHTh2oGwDzSf/+YQ/rvojOw7uoFfrXmWNxBNjE+kd05vQ\n4FAXqxcRkfqguLiYK664gpUrV9K/f3/Gjx/P7t27ueGGGxg2bFiNr/PZZ5+xc+dOBg4cyMUXX1zt\neZMnT+bRRx/ltdde45lnniE83Lft3KhRo1i/fj2XX345ERERPquP1lquvfZa1q9fz+jRo8nLy2PZ\nsmWMHTuW7OxsJk2aVO33TUhI4L777mP+/Pn069ePUaNGlb3WqVMnAFJSUnjyyScZNmwYF1xwAcYY\nvvrqK5544glSUlJYtWoVwcH1J4LVn0pd1KULPPQQ/PKX0KSJ29XUjqiwKIZ0HsKQzkPKjh3JP8Ka\nzDWkelP56oeveOrfT7Fl3xbOPuNsn9vYfdv0pVmTZi5WLyIigWbx4sWsXLmSUaNG+awmjhs3jqFD\nh9b4OqtXrwY44XsiIyNJTExk9erVpKamcuGFF/q8np2dzbp16yrdQi61Y8cO1q9fXxYyp06dSkJC\nAr/73e/4xS9+QWxsbJXv69u3L5MmTSoLizNmzKh0zrBhw8jKyiKswkpTcnIy06dP59VXX+XGG288\n7o8vkOg2dA289Zbz0a0bLF0KRUVuV+QfzZo0Y1DHQdx73r0sumoRa+5cw/4H9/PcyOc4v935rN+7\nnl9//GtiHo+h+4Lu3PjmjTz+5eN89t1nHDh2wO3yRUTERaW3cOfOnetzfPDgwVxyySU1vk5mZiYA\n7du3P+G5pauFWVlZlV6bPXt2tUERYNq0aT6rkbGxsdx3333k5eXx6quvHvf7nuhZxejo6EpBEWDi\nxIkALF++/LjvDzRaWayBc8+Fjz+GlBSYNg3+8AeYMwd+/nPwNPC4HRocyrntzuXcdueWHSsoKmBj\n9kbSvemkelN5Z9M7ZGRm0Cq8lU8z8cTYRGKaxbhYvYhI/VDXDSv8sS9jzZo1tGjRgh49elR6bdCg\nQfzjH/+o/W9aDWMM55xzznHPqbgSCU6dAGvXrj3tGkpvj2dkZHDw4EGfvo5er/e0r1+XFBZPwkUX\nOYHx44+d29KPPAJz58Lll9f9/+huCgkKoU9MH/rE9GFcP6dBaVFxEVv3by3bif3Yl4+RnpnOT0J+\nUmkndvuI9mrlIyJSTkPYVJuTk8NZZ51V5WsxMTVfOChtUfP999+f8Nxdu3YBVHnLuHXr1tW+zxhT\n5euldebk5NSo1uo89thjPPjgg8TExDBixAjatWtHaGgo1lpmz55NXl7eaV2/riksniRj4Gc/g8su\nc25N338/JCc7H0lJblfnniBPEF2ju9I1uitjeo8BnGX67Qe3l02jWZi6kDRvGsW2uNI0mrioOLXy\nERGpxyIjI9m7d2+Vr1V1m7g6AwcOBJxbtXPmzKn2vJycHFJTU2nSpEmVO6aPx1rLnj17aNeuXZV1\nRkae+ojdwsJC5s6dS7t27VizZg0tW7b0uf7s2bNP+dpuUVg8RcY4t6Gvugr+9je47Tbo3NlZafRz\nf9B6wxhD56jOdI7qzDU9rgGc/0G9R7xOgPSm89d1f2XKJ1PIycvx2YWd0CaBbtHdCPIEufyjEBGR\nmujXrx+ff/4569evp2fPnj6vffHFFzW+zpAhQ+jYsSOrV68mJSWFiy66qMrz5s2bR15eHjfffDOh\noSfftSMlJYUxY8ZUWWffvn2P+96gIOfvpqIqNjFkZ2dz+PBhLrvsMp+gCLBq1aqTrjMQaCnnNAUF\nwU03wcaNcN11cO21cOWVUAuPOzRIxhjaNm/LFWdfwfSLp/PWL95i56938p97/8O0C6cR85MY3t/8\nPqNeHUXkHyL56fM/5a4P7uL5tOdJ96aTX5Tv9g9BRESqcNNNzoSxhx56yGcDyGeffcann35a4+sE\nBQXx5JNPAjB69GjS0tIqnbNkyRKSk5Np2bIlycnJp1RvcnIyP/74Y9nXXq+X+fPnExoayvXXX3/c\n95aOGdy5c2el11q3bk1oaCipqank5ub6XL9iT8j6QiuLtSQkBCZMgLFjYeFCuPRS57b07NnQtavb\n1QW+6PBoLjnrEi4567875g7lHSIjM4M0bxord6xk3lfz2HZgG92iu/mMM+wT04fwkPDjXF1ERPxt\n/PjxvPjii7zzzjuce+65DBs2DK/Xy7JlyxgxYgQffPBBja81atQoFixYwL333suAAQO49NJL6d27\nNwUFBaSkpJCWlkarVq149913q9w1XZPJKnFxcfTq1YtrrrmmbAf0vn37+NOf/lRt25xSzZo1Y8CA\nAaSkpDB27Fji4+PxeDyMHTuWjh07MnHiRObNm0dCQgLDhw9n//79fPDBByQlJbF58+Ya/zwEDGut\nPo7z4fwUnbzDh6195BFro6OtveUWa7dvP6XLSAU/5v9oV3+/2i741wJ76zu32sRnEm3Y3DDbc0FP\ne/ObN9t5q+fZldtX2pzcHLdLFZFG6FT/zmgojhw5YidPnmxjY2NtWFiYTUxMtG+++aZdvHix9Xg8\ndsmSJSd1vY0bN9oJEybY+Ph4GxYWZiMiImxCQoKdMWOGPXDgQJXvSUpKsh6Pp9prlr6em5trp0yZ\nYtu3b2+bNm1q+/TpY1988cVK51dX++bNm+2IESNsVFSU9Xg81uPx2JUrV1prrc3Pz7dz5syx8fHx\nNjQ01MbHx9vZs2fb/Px8a4yxgwcP9rnW+PHjrcfjsTt27Djuz0dNfn+VnFOrWcjYhrAFy4+MMfZ0\nfo4OHoTHH4enn4YbboCpU+EE/2CRk5RflM+6PetI9zoTadIy0/gm6xtim8dW2kgTHR7tdrki0oAZ\nYxrNvGCpezX5/VVyTq22HFFYPIHTDYul9uxx+jMuXuxshnnwQTjjjNOvT6pWWFzIpuxNZSMN07xp\nZGRmENE0wucWdmJsIrHNYtXKR0RqhcKi+JPCYoCqrbBYatcuZ8f066/D3XfD5MlwnAbzUouKbTHf\nHfiOVG+qswqZmUbq7lSCPcGVekF2atFJAVJETprCoviTwmKAqu2wWGrbNmfzy0cfwW9+4wTHcO3R\nqHPWWnYd2uUTINO8aRwrOFYpQHY5o4t6QYrIcSksij8pLAYof4XFUhs2wMyZsGqV8zzj7bdD06Z+\n+3ZSQ1lHskjPTCd1d2rZrey9R/fSr00/n36Q3aO7ExIU4na5IhIgFBbFnxQWA5S/w2KptDSYPh3W\nr4cZM5wWPMFqbBRQDhw7QEZmhrMKWRIgd+bspGernj7zsHvH9CY0+OQbxIpI/aewKP6ksBig6ios\nllq1CqZNA6/XuU19/fXg0Z3PgHUk/whrMteQ5k0rC5Fb9m3h7DPO9rmN3bdNX5o1aeZ2uSLiZwqL\n4k8KiwGqrsMiOAPlP/0UHnoIcnPh4Ydh5EhnxKAEvtzCXL7J+qZs9THVm8qGvRvoGNmxrJVPQmwC\nCW0SiAqLcrtcEalFCoviTwqLp8EYczZwDzAU6AAcBr4GpltrKw3eM8bcDkwBOgHbgXnW2mequXad\nh8VS1sJ77zmhMTzc2UU9dKhCY31UUFTAxuyNZb0gU72prMlaQ6vwVj5tfBLaJBDTLMbtckXkFCks\nij8pLJ4GY8zdwJ3AYuDfQAvgAaAfcIG1Nq3cubcDC4FHgE+BYcBU4C5r7cIqru1aWCxVXAzLljnP\nMrZrB8nJMHCgqyVJLSi2xWzZt4U0b5pPP8iwkLBKzcQ7RHRQKx+RekBhUfxJYfE0GGPOsNbuq3As\nAmfV8D1r7biSY8HAbuADa+0t5c59HrgSiLXWFla4juthsVRhISxd6jzL2KuXs9KYkOB2VVKbrLXs\nyNlRFhxLVyGLiot8ViD7x/YnLipOAVIkwCgsij8pLPqBMeafwCFr7SUlX18IrAQusdYuL3deEvAZ\nMMRau6LCNQImLJbKy4PnnoNHHoFBg2DOHOje3e2qxF+stXiPeH0CZHpmOjm5OfRr088nRHY9oytB\nniC3SxZptBQWxZ8UFmuZMaYl8D3wvLX23pJjdwJP4awgZpU7tzWQiXMr+ukK1wm4sFjq6FH4y1+c\n2dOXX+70a4yLc7sqqSvZR7N95mGne9PZfXg3vWN6+zQT79m6J02CmrhdrkijoLAo/qSwWMuMMS8D\nVwF9rLXbSo5NBeYCodba/HLnBgP5OBtikitcJ2DDYqmcHJg3zwmO113nbIhp187tqsQNObk5rMla\n47MKue3ANrpFd/NZgewT04fwEI0MEqltCoviT26FxYDs4GeMGWaMKa7Bx2fVvP93wBjg7tKg2JBF\nRsKsWfDtt9C8OfTpA1OmwN69blcmdS0yNJKLzryISedPYunVS1n3q3VkP5DNUyOeon9sf9K96dz1\n4V1EPxpNz6d6cvNbNzNv9TxWbl9JTm6O2+WLiPhNUlISnpNoXLx48WI8Hg9LlizxW03jx4/H4/Gw\nc+dOv32P2hCoM0JWAd1qcN7RigdKbjUnA9OstYsrvHyg5HMUkFXueMuSz/ur+iazZs0q+++kpCSS\nkpJqUFrdi46GRx+FX//a2THdrRv86ldOcGzRwu3qxC3hIeGc3/58zm9/ftmx/KJ81u9ZXzYTe9mG\nZazNWkvb5m0r7cSODo92sXoRaYw2bdrEvHnzWL58Obt37yY4OJizzjqLkSNHMmnSJKKiTq1H7cls\nCjTGlH2U2r59O3FxcYwbN45FixadUg2nU1NVVqxYwYoVK2qlluo0qNvQxpibcdrnPGGtfaCK1y8C\nVlD9BpfB1tqVFd4T8Lehq7N9u7P55b33nAB5773QTENEpBqFxYVs3rfZZx52emY6LUJb+IwzTIxN\nJLZZrHZii1RBt6FP37PPPstdd90FwCWXXELfvn3Jz89nxYoVpKen06pVK959913OO++8k7puUlIS\n//d//0dRUVGNzj906BCZmZm0adOGiIgIAHbs2EHnzp0ZP348L7zwwsn9wKowfvx4li5dyvbt2+nY\nseMJz3frNnSgriyeNGPM1cALwHNVBcUSXwLZwI3A8nLHbwL24axoNhidOsELL8CmTc7ml/h4+O1v\n4c47IVSji6WCYE8wPVr1oEerHtzc92bA6QX53YHvylr4/M+//ofU3akEe4IrBchOLTopQIrIaXnn\nnXe48847iYmJ4f3336d///4+ry9atIgJEyYwfPhwMjIy6NChg99qiYiIKAuJpUqDWqP7B4G1tt5/\nABcBuTgNuX8KnF/uI6HCuXcARcDDQBIwp+TridVc2zYUGRnWjhxpbfv21j77rLX5+W5XJPVRcXGx\n3Xlwp31749t2+mfT7RV/vcK2faKtjfpDlB2yZIj9zce/sX9d+1f77d5vbVFxkdvlitSphvR3xqlY\ntmyZveCCC2x0dLQNDQ217dq1s8OHD7f/+Mc/TvjegoICe+aZZ1qPx2NXrFhR7XkzZsywxhh70003\n+Ry/+OKLrTHG5ubm2t/97ne2U6dONiQkxD755JNlr3s8Hnvs2DE7efJk265dOxsaGmr79u1rX3zx\nxUrfZ9GiRdYYY5csWeLzdVUfK1eutNZa+8MPP9jp06fbc88910ZHR9umTZva+Ph4+5vf/MYePny4\n0vcYN26cNcbYHTt2nPDnx9qa/f4qOadWc1ZDWVkcDDQBEqi8OrgdKGsoY619xhhjccb93Q/soJrp\nLQ1N377w7rvw1VfOjuk//tFp8D16NASpNZ/UkDGGDpEd6BDZgau6XVV2PPNIJunedNIz03l94+tM\n+2wae4/udXpBlszDToxNpHt0d0KCQlz8EYiIPyxYsIB77rmH+Ph4xowZQ/Pmzfnhhx/48ssv+eST\nTxg2bNhx3//ZZ5+xc+dOBg4cyMUXX1zteZMnT+bRRx/ltdde45lnniE83Lezw6hRo1i/fj2XX345\nERERPquP1lquvfZa1q9fz+jRo8nLy2PZsmWMHTuW7OxsJk2aVO33TUhI4L777mP+/Pn069ePUaNG\nlb3WqVMnAFJSUnjyyScZNmwYF1xwAcYYvvrqK5544glSUlJYtWoVwcH1MHrVdvpsaB804H8lfvaZ\ntT/9qbU9e1r75pvWFhe7XZE0NPuP7rfLty23j6963N7wxg226/90teHJ4fbcZ8+1d7x3h33m38/Y\nf+36lz1WcMztUkVqRUP+O+NEEhISbPv27e2xY5X/f963b98J3z9r1ixrjLEzZsw44bkDBw60xhib\nkpJSdqx0ZfGcc86xOTk5ld5T+nqvXr3sjz/+WHZ89+7dNiYmxoaGhtrdu3eXHa+4smittdu3b7fG\nGHvLLbdUWdfevXvt0aNHKx2fO3euNcbYl156yee4VhYl4A0eDKtWwYcfOiuNjzzijBC89FLQo2dS\nG6LCohjSeQhDOg8pO3Yk/whrMp1ekF/t+ooFXy9gy74tdDmjS9lO7ITYBPq16UezJtqRJY2DmV23\nf+jambX/zJ0xhpCQkCrb07Rs2bKKd/jKzMwEoH379ic8t0OHDqxevZqsrKxKr82ePbvSs4blTZs2\nzWc1MjY2lvvuu49p06bx6quvHnd10cli1YuOrrp7xMSJE5k+fTrLly/nxhtvPO41ApHCYiNnDIwY\n4UyAeeMNmDQJWrVyWu9ceKHb1UlD1KxJMwZ1HMSgjoPKjuUW5rJuz7qyRuIvrn2RdXvWcWaLMyu1\n8mkRqj5Q0vD4I7zVtV/84hf89re/pVevXowePZqkpCTOP//8SreJ/ckYwznnnHPccy6s4i+3QYOc\nP4/Wrl172jWU3h7PyMjg4MGDFBcXl73m9XpP+/puUFgUADweZ/rL1VfDyy/D2LHQtauz0niC/+9E\nTltocCjntD2Hc9r+9zdbQVEBG7M3lgXItze9TUZmBq3CW/lMo0lok0BMsxgXqxcRgAceeICWLVuy\ncOFC5s6dy9y5c2natCnXXnstTzzxBK1btz7u+9u0aQPA999/f8LvtWvXLsBZFazoeN/HGFPl6zEx\nzp8hOTmnN5zgscce48EHHyQmJoYRI0bQrl07QkNDsdYye/Zs8vLyTuv6blFYFB/BwTBuHIwZA88/\nD6NGwYABTr/GXr3crk4ak5CgEPrE9KFPTB/G9xsPQFFxEVv3by0LkI9/+Thp3jTCQsJ8ViATYxNp\nH9FerXxE6thtt93GbbfdRnZ2NitXrmTp0qW8/PLLfP/99ydsHD1w4EAAli9fzpw5c6o9Lycnh9TU\nVJo0aVKptc6JWGvZs2cP7SrMxC29nR0ZGXlS1yuvsLCQuXPn0q5dO9asWeNz6z0rK4vZs2ef8rXd\nprAoVWrSBCZOhPHj4amnYOhQGDbM2T0dH+92ddJYBXmC6Brdla7RXRnTewzg/OG//eD2skbiz6Q+\nQ6o3lWJbXNYHsn9sfxJjE4mLilOAFKkD0dHRXHPNNVxzzTUkJCSQkpLC4cOHad68ebXvGTJkCB07\ndmT16tWkpKRw0UUXVXnevHnzyMvL4+abbyb0FJoGp6SkMGbMGJ9jX3zxBQB9+/Y97nuDSlqHVNXY\nOzs7m8OHD3PZZZdVekZz1ar63cY5IGdDS+AIC3PGBW7dCt27w/nnw+23Qw3uEojUCWMMnaM68/Pu\nP2fukLl8eOOHZE7JJOOODO49717CQ8J5+ZuXGbJ0CC3+2IKkxUlM/ngyL619iQ17N1BUXLNpDiJy\nfCtXrqx07NixY+Tk5BASElIWtKoTFBTEk08+CcDo0aNJS0urdM6SJUtITk6mZcuWJCcnn1KdycnJ\n/Pjjj2Vfe71e5s+fT2hoKNdff/1x31s6ZrCqWc6tW7cmNDSU1NRUcnNzfa4/derUU6o1UGhlUWqk\neXNnx/SvfgWPPeb0bLz5Zpg6FWL0uJgEGGMM7SLa0S6iHVecfUXZ8eyj2c4YQ286721+j9krZ+M9\n7KVPTB+faTQ9W/ekSVATF38EIvXPVVddRVRUFOeddx4dO3YkNzeXjz76iO3bt3PPPffUaKPLqFGj\nWLBgAffeey8DBgzg0ksvpXfv3hQUFJCSkkJaWlrZuL+qdk2faLcyQFxcHL169eKaa64hLy+PV199\nlX379vGnP/2pymcgy2vWrBkDBgwgJSWFsWPHEh8fj8fjYezYsXTs2JGJEycyb948EhISGD58OPv3\n7+eDDz4gKSmJzZs3n7C2gFXbvXga2geNuGfW8Xi91t5zj7VRUdb+9rfW1qCFlkhAOnjsoF25faX9\n05d/sje/ebPtuaCnDZsbZhOfSbS3vnOrXfCvBXb196vtj/k/nvhi0ug15r8znn76aXvllVfaTp06\n2dDQUNu6dWt74YUXVuotWBMbN260EyZMsPHx8TYsLMxGRETYhIQEO2PGDHvgwIEq35OUlGQ9Hk+1\n1yx9PTc3106ZMsW2b9/eNm3a1Pbp06fKCS6LFy+2Ho/Hp8+itdZu3rzZjhgxwkZFRVmPx2M9Hk/Z\nBJf8/Hw7Z84cGx8fb0NDQ218fLydPXu2zc/Pt8YYO3jwYJ9rjR8/3no8noDvs2hsDVJ4Y2aMsfo5\nqt7Onc7ml7ffhvvuc1rvHOeRFJF64WjBUb7J+oZUbyrp3nTSMtPYuHcjcVFxPjOxE2ITiGhafT83\naXyMMTVa3RI5FTX5/VVyTq0+nK2weAIKizWzZQvMnAnLl8MDDzi3q8PC3K5KpPbkF+Wzfs960rxp\nTojMTGdt1lraNm9bqRdkdHjVjXml4VNYFH9SWAxQCosn55tvYMYM+PprmDYNbr3V2Vkt0hAVFhey\nKXsT6ZnppO52AmR6ZjotQlv4TKNJjE0ktlmsdmI3AgqL4k8KiwFKYfHUfP21syGmdMXxppvgBBvh\nRBqEYlvMtgPbnNvX5VYhg0yQTyPxxNhEOrXopADZwCgsij8pLAYohcXTk5LirDBmZzvPNl5zjTMt\nRqQxsdbkNkvZAAAgAElEQVSy69AuZyd25n9D5LGCY87KY7lb2F1adiHIo39Z1VcKi+JPCosBSmHx\n9FkLH3/shMbiYmeE4PDhzlxqkcYs60hWWXgsDZDZR7PpG9PXZ6Rh9+juhASFuF2u1IDCoviTwmKA\nUlisPdbCm286zzS2aAHJyZCU5HZVIoHlwLEDZQGy9PPOnJ30bNXTJ0D2at2L0OCTn14h/qWwKP6k\nsBigFBZrX1ER/O1vMGsWdO7srDSed57bVYkEriP5R1iTucZZgcx0mopv3reZLmd08dmJ3bdNX5o1\naeZ2uY2awqL4k8JigFJY9J+CAli0CB5+GBITnc99+rhdlUj9kFuYy7o968puYad501i3Zx0dIzv6\nrEAmtEkgKizK7XIbDYVF8SeFxQClsOh/ubmwcCH84Q8weDDMng1nn+12VSL1T0FRARv2bvB5DnJN\n1hpahbfyCZCJsYm0/klrt8ttkBQWxZ8UFgOUwmLdOXIE/vxnmDcPrrzSebbxzDPdrkqkfisqLmLL\n/i1lrXxKb2OHh4T77MROjE2kfUR7tfI5TQqL4k8KiwFKYbHuHTwIjz8OTz8NN9wAU6fCCWa7i8hJ\nsNayI2dHWSPx0lXIIltU1gOy9CMuKg6PUb+rmlLYFn9TWAxACovu2bPHuTW9eDHcdhs8+CCccYbb\nVYk0XN7DXp9G4mneNA7mHiShTYJPiOwa3ZVgT7Db5YpIFRQWXaCw6L5du5wd06+/DnffDZMnQ0SE\n21WJNA7ZR7NJ96b7BMjdh3fTJ6aPzy3snq170iRIsz1F3Kaw6AKFxcCxbZuz+eWjj+D+++GuuyA8\n3O2qRBqfnNwc1mSt+e8qpDedbQe20b1Vd59pNH1i+hAeov9JReqSwqILFBYDz4YNzuaX1aud5xlv\nuw2aNnW7KpHG7WjBUb7J+sbnNvbGvRuJi4rzaePTr00/IkMj3S5XpMFSWHSBwmLgSkuDhx76b3gc\nOxaC9RiVSMDIL8pn/Z71PvOwv8n6htjmsT7NxBNiE4gOj3a7XJEGQWHRBQqLgW/VKmfutNfr3Ka+\n/nrwaPOmSEAqKi5i075NPvOwMzIzaBHaotJO7LbN27pdrki9o7DoAoXF+sFa+PRTZ6UxN9eZBjNy\nJKiLhUjgK7bFbDuwrSxApmemk7o7lWBPcKVm4mdGnqn2NCLHobDoAoXF+sVaeO89JzSGhzu7qIcO\nVWgUqW+stew6tMsnQKZ50zhacLRSM/EuZ3RRL0iREgqLLlBYrJ+Ki2HZMudZxnbtIDkZBg50uyoR\nOV1ZR7J8GomnedPYe3QvfWP6loXH/rH96d6qu3pBSqOksOgChcX6rbAQli51nmXs1ctZaUxIcLsq\nEalN+4/tJyMzw2cVcmfOTnq26ulzC7tX616EBoe6Xa6IXyksukBhsWHIy4PnnoNHHoFBg2DOHOje\n3e2qRMRfjuQfISMzw5mJnemEyC37ttDljC4+O7H7tulLsybN3C5XpNYoLLpAYbFhOXoU/vIXZ/b0\n5ZfDzJkQF+d2VSJSF3ILc/km6xuf29jr966nY2RHn53YCW0SiAqLcrtckVOisOgChcWGKScH5s1z\nguN11zkbYtq1c7sqEalrBUUFfJv9bdkkmrTMNDIyM2gV3qrSRpqYZjFulytyQgqLLlBYbNiys+HR\nR+H552H8ePjtb6FVK7erEhE3FRUXsXX/Vp9pNGneNMJDwn1uYSfGJtI+or1a+UhAUVh0gcJi4+D1\nOjum//Y3+NWvYMoUaNHC7apEJFBYa9mRs8MJkLudAJnqTaXYFpfdui4NkHFRcWrlI65RWKwhY8xo\n4K/AD9baDlW8fjswBegEbAfmWWufqeZaCouNyPbtzuaX996DyZPh3nvhJz9xuyoRCUTWWrxHvM7t\n63KrkAdzD5aFx9LPXaO7qpWP1AmFxRowxrQAvgWKgUJrbccKr98OLAQeAT4FhgFTgbustQuruJ7C\nYiO0aZOz+WXlSufW9B13QKg6bohIDew7us9nHna6N53dh3fTO6Y3iW2cWdiJsYn0bNWTpsFN3S5X\nGhiFxRowxjwLdAAygWHlVxaNMcHAbuADa+0t5Y4/D1wJxFprCytcT2GxEVuzBqZPh4wM5/P48RAS\n4nZVIlLfHMo75NMLMtWbyncHvqNbdDefXpB9YvoQHhLudrlSjyksnoAxZhDwCdAbmAEMrRAWLwRW\nApdYa5eXO54EfAYMsdauqHBNhUXhq6+cHdPbtzsNvkePhqAgt6sSkfrsaMFR1matdRqJe51nIL/N\n/pa4qDifNj792vQjMjTS7XKlnlBYPA5jTAiQASyz1s42xiymcli8E3gKZwUxq9zx1jgrkXdZa5+u\ncF2FRSnz+ecwbRocOuQ823j11Zo7LSK1J78on/V71vvMw16btZbY5rE+O7ETYhOIDo92u1wJQAqL\nx2GMeQgYC/Sy1uZXExanAnOBUGttfrnjwUA+MN1am1zhugqL4sNa+PBDZ6UxKMgZIXjZZQqNIuIf\nRcVFbNq3yWcednpmOpFNI31uYSfGJtK2eVu3yxWXNZqwaIwZhnM7+URWWGuHGGPigbXAKGvtJyXX\nWIxzW7lsg4vCotSm4mJ44w2YMcPpzTh3Llx0kdtViUhjUGyL+c/+//hMo0nzphHsCa4UIM+MPFO9\nIBuRxhQWw3A2qZzIUWvtLmPMhzi7n28ESn+CngIuAnoCedbaXGPMRGABJ3kbeubMmWVfJyUlkZSU\ndMo/Nml4Cgvh5Zdh1izo2tUJjeec43ZVItLYWGv5/tD3Za180jKdZyGPFhytNI2myxld1AuygVix\nYgUrVqwo+3r27NmNIyyeLGPMd8CZxznlSWvtZGPMRcAKqt/gMthau7LCtbWyKDWSn+9MgklOhgED\nnGcae/VyuyoRaeyyjmT5PAOZ5k1j79G99I3p67MC2T26OyFBavdQ3zWalcWTZYw5DyjfrMoAvwX6\nA9fiNOf+T8kmmB+A9621vyz3/v8FrkKtc6QWHDsGTz3ljBEcNszZPR0f73ZVIiL/deDYAdIz08vm\nYad509hxcAe9WvfymUbTO6Y3ocFqMlufKCyehKo2uJQcvwPnFvUjwHJgCDANuLviLeiS8xUW5ZQc\nPgxPPgnz5zu7pmfMgA41ebhCRMQFR/KPsCZzjc887C37ttDljC4+O7H7tulLsybN3C5XqqGweBKM\nMYtwwmLHKl6bgDPu70xgB864v0rTW0rOVViU07J/Pzz+ODzzDNx0E0ydCjExblclInJiuYW5rNuz\nrmwedpo3jfV719MxsqPPOMOENglEhUW5Xa6gsOgKhUWpLVlZ8Pvfw4svwoQJcP/90LKl21WJiJyc\ngqICNmZvLGsknp6ZTkZmBq3CW/kEyMTYRGKa6V/GdU1h0QUKi1Lbvv/e2fzy1ltw330waRI0b+52\nVSIip66ouIit+7f6jDNMz0wnPCTc5xZ2Ymwi7SPaq5WPHyksukBhUfxlyxan3c6nn8IDD8CvfgVh\nYW5XJSJSO6y1bD+4vez2dao3lTRvGsW2uFKAjIuKU4CsJQqLLlBYFH/75htn88vXXzujBG+9FZo0\ncbsqEZHaZ63Fe8RbNg87LTON1N2p5OTl+Ny+TmiTQLfobgR5gtwuud5RWHSBwqLUla+/dkYIbtkC\nM2c6m2GC9OekiDQC2UezSfem+6xCeg976R3T22cFsmfrnjQJ0r+mj0dh0QUKi1LXUlKcFcbsbOfZ\nxmuuAY8GLYhII3Mo7xAZmRk+4wy3HdhGt+huPs3E+8T0ITwk3O1yA4bCogsUFsUN1sLHHzsrjUVF\nzgjB4cNBj/SISGN2tOAoa7PWloXH9Mx0Nu7dSFxUnE+A7NemHxFNI9wu1xUKiy5QWBQ3Wevsmp4+\nHSIjnVGCgwe7XZWISODIL8pn3Z51PjOxv8n6htjmsT4baRJiE4gOj3a7XL9TWHSBwqIEgqIi+Nvf\nnN3TnTo5K43nn+92VSIigamwuJBN2Zt85mFnZGYQ0TTCZwUyMTaR2GaxDWontsKiCxQWJZAUFMCi\nRfDww9CvnxMa+/Z1uyoRkcBXbIv57sB3Tg/Icjuxgz3BJMQm+Gyk6dSiU70NkAqLLlBYlECUmwsL\nF8If/gBJSTB7NnTt6nZVIiL1i7WWXYd2+QTING8axwqOVQqQXc7ogscE/m5DhUUXKCxKIDtyBP78\nZ5g3D0aOdFrunHmm21WJiNRvWUeySM9M95mJvffoXvq16efTD7J7dHdCgkLcLteHwqILFBalPjh4\nEB5/HJ5+Gm64AaZOhdhYt6sSEWk4Dhw7QEZmRtkowzRvGjtzdtKzVU+fmdi9Y3oTGhzqWp0Kiy5Q\nWJT6ZM8e59b04sVw223w4INwxhluVyUi0jAdyT/Cmsw1PvOwt+zbwtlnnO1zG7tvm740a9KsTmpS\nWHSBwqLUR7t2OZtfXn8d7r4bJk+GiMbZckxEpE7lFubyTdY3PtNoNuzdQMfIjmWtfBJiE0hok0BU\nWFStf3+FRRcoLEp9tm2bs/nlo4/gN79xgmO4Bh2IiNSpgqICNmZvLOsFmepNZU3WGlqFt/Jp45PQ\nJoGYZjGn9b0UFl2gsCgNwYYNzuaXVauc5xlvvx2aNnW7KhGRxquouIit+7eWTaIp7QcZFhJWqZl4\nh4gONW7lo7DoAoVFaUjS0pwRghs2wIwZMHYsBAe7XZWIiIDTymdHzg6fedip3lSKiosqNROPi4qr\nspWPwqILFBalIVq1CqZNA6/XuU19/fXgCfz2YSIijY61Fu8Rr8887DRvGofyDrF78m7CQsJ8zldY\ndIHCojRU1sLy5U5ozM11psKMHAn1dGiBiEijcuDYgSo3yCgsukBhURo6a+G995zb0+Hhzi7qoUMV\nGkVE6iOFRRcoLEpjUVwMy5Y5zzK2awfJyTBwoNtViYjIyfBHWNRTSiICOM8sjh7tbH65+WYYMwZG\njID0dLcrExERNyksioiP4GD45S9h82a4/HInMF53HWzc6HZlIiLiBoVFEalS06ZOE++tW+Hcc+Hi\ni2HcOKfRt4iINB4KiyJyXOHh8MADsGULdO4MAwbAxInwww9uVyYiInVBYVFEaiQyEmbNgm+/hebN\noXdvmDIF9u51uzIREfEnhUUROSnR0fDoo7B+PeTlQbduMH06HDzodmUiIuIPCosickpiY+Evf4HU\nVOeWdJcu8Pvfw48/ul2ZiIjUJoVFETktnTrBCy/AF1/AmjUQHw/z5ztTYUREpP5TWBSRWtG1K7zy\nCvz9784YwbPPhueeg4ICtysTEZHTobAoIrWqb194911nGsyrr0L37vDyy1BU5HZlIiJyKjTu7wQ0\n7k/k9Hz+OUybBocOwcMPw6hRmjstIuIvmg3tAoVFkdNnLXz4ITz0kDMhZu5cuPRShUYRkdqmsOgC\nhUWR2lNcDG+8ATNmQKtWkJwMF17odlUiIg2HwqILFBZFal9hofMc46xZzsaYuXPhnHPcrkpEpP7z\nR1hsUBtcjDHtjDEvGGO8xphcY8w2Y8wjVZx3uzHm25JzvjXG3OFGvSKNVXCwM2d60ya46irnOcaf\n/xzWrXO7MhERqajBhEVjTCfgX0A8cA9wCTALKKhw3u3AQuA14LKSz08ZY+6su2pFBKBJE2fO9JYt\nMGgQDB0KN94IW7e6XZmIiJRqMLehjTF/B1oAg6y1VTbpMMYEA7uBD6y1t5Q7/jxwJRBrrS2s8B7d\nhhapI4cPw5NPOk29r77aebaxQwe3qxIRqT90G7oaxpizgEuB/6kuKJb4KRANvFTh+IvAGcAF/qlQ\nRGqieXNnzvTmzc4M6n79YNIkyMpyuzIRkcarQYRFYFDJ51xjzD9KnkXcb4xZYoxpWe68niWfKz4Z\ntaHkc3e/VikiNdKypTNnev16p+1Ojx7wu9/B/v1uVyYi0vjUKCwaY1YbY8YaY5r6u6BT1Lbk8wvA\nt8DPgAeBEcDHxpR1cysNjgcqvH9/hddFJAC0aePckk5Ph+xsZ4Tgww87t6tFRKRuBNfwvDxgMTDP\nGLMUeMZa+62/ijLGDAM+qcGpK6y1Q/hv6P3cWntP6WvGmBzgFZyNLH8/1XpmzZpV9t9JSUkkJSWd\n6qVE5BR07OjMmX7wQafdTnw8PPAA/OpXEBbmdnUiIu5ZsWIFK1as8Ov3qPEGF2NMN2ACMA6IAlJw\ndhW/Ya0tON57T7ooY8KAmjzWftRau6uk9c3TwD3W2gXlrhMF7AN+Z639ozFmIrAAZyNLVrnzWgOZ\nwF3W2qcr1KINLiIBZt06Z/PLv/7ljBK89VZnZ7WISGPn6gYXa+231trJQDucwBgM/BXYZYz5ozEm\nrraKstYes9ZursHHrpK31LQ72/qSz70qHO9R8nkDIhLwevWCN9+Et96Ct9+Gbt1gyRIoOt72NhER\nOSUnvcHFWptrrX0RuBf4P6AVcD+wxRjzujGmTS3XWBNf4awM/qzC8dKvvy75vBrIBm6scN5NOCuQ\nq/xVoIjUvnPPhY8/hsWL4X//1wmRr73mjBUUEZHacVJ9Fo0x4cAY4E6gP7AJ5/bv6zibSWYD35Y8\nR1injDFjcZ6rfAZ4C6c591wg3Vo7tNx5dwBPAY8Ay4EhwDTg7oq3oEvO121okXrAWic4PvSQs8I4\ndy4MHw6mVm/GiIgENtdmQxtj+gB34KzIhQPvAE9baz+rcN5I4HVrrSu7po0xN+Hsgu6Cs1L4Os7z\nikcrnDcBmAKcCewA5llrF1ZzTYVFkXrEWuf29PTpEBkJyckweLDbVYmI1A03w2IxzuST54BnrbXe\nas7rASyw1jaYP5oVFkXqp6Ii+NvfYOZM6NzZCY3nned2VSIi/uVmWLwGePsE01EaJIVFkfqtoAAW\nLXL6MyYkOJ/79nW7KhER/3BtN7S19o3GGBRFpP4LCYEJE2DLFhgyBC67DEaPhk2b3K5MRKR+qOnK\n4kyguhOLgRwgzVrb4HYTa2VRpGE5cgT+/GeYNw9GjnRuU595pttViYjUDrefWayJ1cBwa23OaVUV\nQBQWRRqmgwfh8cfh6afhhhtg6lSIjXW7KhGR0+NmU+4ewFb+u4M4DOhESX9FYBBwPdAd+H1tFigi\n4g8tWjjtdTZudG5V9+zpjBDct8/tykREAktNVxaXA59Ya/9YxWsPApdZa4cYYx7AGblXk1F99YJW\nFkUah127nPD42mtw773w619DRITbVYmInBw3VxbPB/5dzWtpJa8DpAIxp1uUiEhda98eFi505k3/\n5z/QpQs89hgcPXri94qINGQ1DYuHgGHVvDYUZ4MLQGjJuSIi9dJZZ8HSpfD55/DPfzqhccECyMtz\nuzIREXfUNCw+DzxgjPmLMeZiY0z3ks8LcJ5bfL7kvPOAb/xRqIhIXerRA15/Hd57Dz74ALp2hRde\ngMJCtysTEalbNX1mMQhn7vMknHF/pX4EngRmWmuLjTEDgCPW2g3+KNYNemZRRABWrYJp08Drhdmz\n4frrwVPTf26LiNQR11rnlCsgCugNxAJe4Btr7YHaLCjQKCyKSClrYflyJzTm5jrTYEaOBFOrfyyL\niJw6V8KiMaYpkAmMs9a+W5vfvD5QWBSRiqx1bk8/9BCEhztzp4cOdbsqERGXdkNba/OAQiC3Nr+x\niEh9ZQxceSVkZMB998HEic4owS+/dLsyEZHaV9Mnbt4GrvVnISIi9Y3HA2PGwIYNcOONzn+PGAHp\n6W5XJiJSe2q6weVq4H+AfwJv4Tyv6PNGa+1n/ijQbboNLSI1lZcHzz0HjzwCgwbBnDnQvbvbVYlI\nYxLIs6GttTaodkoKLAqLInKyjh6Fv/zFmT19+eUwcybExbldlYg0Bm6GxaQTnWOtXVEL9QQchUUR\nOVU5OTBvnhMcr7vO2RDTrp3bVYlIQ+Z665zGSGFRRE5XdjY8+ij87//CLbfAb38LrVq5XZWINERu\nzoYuLSDaGHOFMWacMeaMkmNhJU27RUSkCtHRTlhct855rrFbN5g+HQ4edLsyEZETq1FYNI7HgR+A\nd4EXgDNLXn4bmOaf8kREGo62bZ1b0qmp8MMPztzp3/8efvzR7cpERKpX05XF3wF34Yz8Ow8ov7z5\nHjCilusSEWmwOnVy5kx/8QWsWQPx8TB/vjMVRkQk0NQ0LN4GPGytfQSo2EHsP0B8rVYlItIIdO0K\nr7wCf/+7M0bw7LOd1jsFBW5XJiLyXzUNi+2A1dW8lg/8pHbKERFpfPr2hXffhWXL4NVXnd6ML78M\nRUVuVyYiUvOwuBvoXc1rfYDvaqccEZHG6/zz4dNPndXFBQucEPnWW84sahERt9Q0LC4DZhhjLqDc\n5BZjTFdgCvCKH2oTEWmUBg+GVavgj390psAMGAAff6zQKCLuqGlT7nDgY2AQsANnJ/R3QAfgS+Ay\na22eH+t0jfosioibiovhjTdgxgynN2NyMlx4odtViUigcrUptzEmGBgD/AxoDWQDfwdettYW1mZR\ngURhUUQCQWGh8xzjrFnOxpi5c+Gcc9yuSkQCjSa4uEBhUUQCSX4+PP+8s8J47rnw8MPQq5fbVYlI\noHB9gouIiLirSROYOBG2bIELLoChQ+HGG2HrVrcrE5GGqqYTXJoaY2YZYzYZY44ZY4orfKjBg4hI\nHQoLgylTnJDYrZuzk/r22+H7792uTEQamppucJmPM8HlI2AdUHEzi7XWzq798tyn29AiUh/s3w+P\nPQbPPgs33QRTp0JMjNtViUhdc+2ZRWPMD8DT1tq5tfnN6wOFRRGpTzIznXnTL70EEybA/fdDy5Zu\nVyUidcXNZxab4bTIERGRANamjTNnOj0dsrOdEYIPPwyHD7tdmYjUVzUNi+8DF/mzEBERqT0dOzqT\nYL76CjZtgvh4eOIJOHbM7cpEpL6p6W3o84AXgZeBD4D9Fc+x1m6r9eoCgG5Di0hDsG6d09j7X/+C\nadPg1ludndUi0rC4+cxi8QlOsdbaoNop6dQYY1oBM4HhQBsgEyfYzrbWZlc493acMYWdgO3APGvt\nM9VcV2FRRBqMf/8bHnoINm+GmTOdzTBBrv7pLSK1yc2wOP5E51hrF9dCPafEGGOA1UAcMB3YCPQE\n5gBbrbU/LXfu7cBC4BHgU2AYMBW4y1q7sIprKyyKSIOTkuKsMGZnO/Onr7kGPOq8K1LvBeQEF2NM\nEBBpra10a7quGGO64gTEO6y1z5U7fgfwNNDVWrulZGThbuADa+0t5c57HrgSiK04ulBhUUQaKmvh\nk0+c0FhU5IwQHD4cTK3+NSMidalOd0MbY/YbYxLLfW2MMe8aY+IqnHousLc2izoFpTdRciocL/26\n9Mf5UyAaeKnCeS8CZwAX+KU6EZEAZAxcdhl8/bXzPOODD8KgQfD5525XJiKB5Hg3HVoAweW+DgKu\nKDlekav/DrXWbgA+AaYbY/obY5oZYwYAM4APrbWbSk7tWfJ5XYVLbCj53N3/1YqIBBZj4OqrYc0a\nuOsuZxLMsGHwz3+6XZmIBIKG9ITK1cAO4GvgEPAVsBW4ttw5pa1pD1R47/4Kr4uINDpBQc6c6Y0b\n4Re/gGuvhSuvdEKkiDRewSc+pe4ZY4bhrBSeyApr7RBjjAd4HegH3IHz/GIPYDbwujFm5Ok8eDhr\n1qyy/05KSiIpKelULyUiEvBCQpzVxZtvhoULnVvVSUkwezZ07ep2dSJS3ooVK1ixYoVfv0e1G1xK\n2uWcb639V8nXwUA+cI61Nq3ceecDX1pra22V0hgTBnSowalHrbW7jDFXAW8BQ621ZU/blAudo6y1\n7xpjJgILcDayZJU7rzVOq527rLVPV6hFG1xEpFE7cgT+/GeYNw9GjnSeb+zUye2qRKQqboz7a2+M\niSvZ1BJX8VjJ8Xa1WRCAtfaYtXZzDT52lbylR8nnf1e41Ncln7uVfF5f8rlXhfNK378BERHx0awZ\nTJ0KW7ZA27bQvz/cfTd4vW5XJiJ14URh8XWc5/62At+WHHu73LGtwGt+q67mSkPjuRWOn1fy+YeS\nz6uBbODGCufdBOwDVvmlOhGRBqBFC6e9zsaNzvSXnj3h/vudXo0i0nAd7zb0+JO4jrXWLqmVik6B\nMaYZzqphMPAwsAlnNXEmkAv0sNYeLTn3DuApnKbcy4EhwDTg7oq3oEvO121oEZEq7NrlhMfXXoN7\n7oHJkyEiwu2qRBq3gGzKHSiMMW2BWcBQIBbnGcR/ALOstd4K507AGfd3Js4O6nlVTW8pOVdhUUTk\nOLZtg1mz4O9/d1Ya77oLwsPdrkqkcVJYdIHCoohIzWzY4Gx++fJL5xnH22+Hpk3drkqkcXFjg4uI\niEiN9OgBr78O778PH37otNlZtAgKC0/8XhEJXFpZPAGtLIqInJpVq5y505mZTo/G664Dj5YoRPxK\nt6FdoLAoInLqrIVPP3VCY16esyHmiiucEYMiUvsUFl2gsCgicvqshXffhenTnc0vyckwdKjbVYk0\nPAqLLlBYFBGpPcXF8OqrMHMmtG/vrDQOHOh2VSINhza4iIhIvebxwJgxzs7pG290/nvECEhPd7sy\nEamOwqKIiNS54GC49VbYvBkuv9wJjNdd50yHEZHAorAoIiKuadrUmTO9ZQuccw5cfDGMGwfffed2\nZSJSSmFRRERc95OfwIMPOqGxUycnOE6cCD/84HZlIqKwKCIiASMy0unJuGkTNGsGvXvDlCmwd6/b\nlYk0XgqLIiIScKKj4bHHYN06pz9jt25O252DB92uTKTxUVgUEZGA1bYt/OUvkJrq3JLu0gV+/3v4\n8Ue3KxNpPBQWRUQk4HXqBC+8AF98AWvWQHw8zJ8PubluVybS8CksiohIvdG1K7zyCvz977B8OZx9\nNjz3HBQUuF2ZSMOlsCgiIvVO377O+MBly5yJMD16wMsvQ1GR25WJNDwa93cCGvcnIhL4Pv8cpk2D\nQ4fg4Ydh1CgwtTrwTKR+0GxoFygsiojUD9bChx/CQw85E2LmzoVLL1VolMZFYdEFCosiIvVLcTG8\n8WixLu0AABgqSURBVAbMmAGtWkFyMlx4odtVidQNhUUXKCyKiNRPhYXOc4yzZjkbY+bOdSbDiDRk\n/giL2uAiIiINUnCwM2d60ya46irnOcaf/xzWr3e7MpH6RWFRREQatCZNnDnTW7bAoEEwZAjcdBNs\n3ep2ZSL1g8KiiIg0CmFhzpzprVud29Lnnw+33w7ff+92ZSKBTWFRREQalebNnTnTmzc7M6j79YP7\n7oOsLLcrEwlMCosiItIotWzpzJkufYaxRw+YOhUOHHC3LpFAo7AoIiKNWps2zpzp9HTYuxe6dHEa\nex8+7HZlIoFBYVFERATo2NGZM716NXz7LcTHwxNPwLFjblcm4i6FRRERkXK6dHH6My5fDl984Xz9\n9NOQn+92ZSLuUFgUERGpQq9e8NZbzsc770C3brBkCRQVuV2ZSN3SBJcT0AQXEREBSElx5k7v3Qtz\n5sA114BHSy4SYDTuzwUKiyIiUspa+OQTmDbNWWGcOxeGDwdTq381i5w6hUUXKCyKiEhF1sLbbzv9\nGiMjndA4eLDbVYkoLLpCYVFERKpTVASvvAIzZ0KnTpCcDOed53ZV0pj5IyzqaQsREZFTFBQEN94I\nGzfC9dfDtdfClVfCmjVuVyZSexQWRURETlNICEyYAFu2wJAhcNllMHo0bNrkdmUipy/gw6IxZrIx\n5j1jjNcYU2yMmXmcc283xnxrjMkt+XxHNeeNMsakG2OOGWO2G2OmGWMC/udCREQCW2goTJoEW7dC\nnz5wwQXwy1/C9u1uVyZy6upDQLoNiAbeKvm6ygcIjTG3AwuB14DLSj4/ZYy5s8J5lwGvA/8EfgbM\nBx4CHvFH8SIi0vg0a+bMmd6yBdq2hf794e67wet1uzKRk1dvNrgYY4KAAmCWtXZOhdeCgd3AB9ba\nW8odfx64Eoi11haWHEsHDlprB5c7bzpOYOxorc2qcG1tcBERkdOyZw/84Q+weDH8f3v3HiVVdSd6\n/PvjoWh0UIxZjq/gCIpAgniZ5TtijC+yVKIm0dG5aqIxOskkd5zRKCjIQ0ezNIvEG8kYdWniDXrx\n/UpiwNYrGs11FCMqwvU5xLcYjCAI7PvHPh3Koqq7we4+Vc33s9ZZVXXOPlW/7t/q7l/vc/bep54K\n55wDW21VdlTqiTb0AS5tfeF7k3sff1m1/xfAVsB+ABGxAzCiTru+wOGdEqkkSRU+8xm4/HJ46ilY\nsgR22QUmTszPpUbXTMViW4YVj09X7X+meNytrXYppZeApRXtJEnqdNtvD9Onw2OPwQsvwKBBcOml\nsHRp2ZFJ9fWUYnFA8bi4av+7VcfrtWvdN6DGfkmSOtXOO8P110NLSy4cBw2Cn/wEli8vOzJpbd1a\nLEbEl4oRze1ts7szrtbwSvhMSdIGbOhQmDkT7roL7r0Xdt0VrrkGVq4sOzJpjT7d/HlzgCEdaLeu\nHfKtPYVbApUDVFp7Ct+t0a7aFhXtPmbixIl/fT569GhGjx69juFJklTfHnvAPffAnDl53elLLoEL\nL8wTfffqKdcA1SVaWlpoaWnp0s9optHQfYAV1B4N/QWgBTg4pTSrYv9oYDZwYErpgYjYEXgJOC2l\ndHVFu4HAC8ApKaXrqt7b0dCSpG6TEvzud7loXL4cJk+GI46A8PqXOmBDHw3dloeBt4ETqvafCLxD\n7tEkpfQKMLdOuxXAvV0bpiRJbYuAgw+GRx+FSZNg/HjYe2+YNav9c6Wu0N2XoddZRIwCBrKmsB0W\nEccWz+9OKS1LKa0s5kr8aUQsAmYBXwROAb7TOsdi4TzgroiYDswARgLjgGkppTe7/iuSJKl9EXDU\nUblX8cYb4Ywz8mjqKVNgn33Kjk4bkoa/DB0R1wInFS8TawaiJGCnorewte23gLOAzwIvAz9KKU2v\n8Z5fASaQ7598Hfg5MLXW9WYvQ0uSGsHKlXDddbm3cfjwXDSOHFl2VGo0XXEZuuGLxbJZLEqSGsny\n5XDVVXDRRbDvvrl43M1ZglXwnkVJkjZwG2+c15leuBD+/u/hgAPgpJPgxRfLjkw9lcWiJElNaNNN\n4eyzYcEC2GknGDUq39e4aFHZkamnsViUJKmJ9e+f15mePx822ww+/3k46yx4662yI1NPYbEoSVIP\n8OlPww9/CH/8Y76vccgQOP98eO+9siNTs7NYlCSpB9l2W7jiCnj88XxJevBguPhi+OCDsiNTs7JY\nlCSpBxo4MK8z/dBDMHcuDBoE06bBhx+WHZmajcWiJEk92K67wowZ8Otf51VgdtklT73z0UdlR6Zm\nYbEoSdIGYMQIuOMOuOmmvCLMbrvBDTfAqlVlR6ZG56Tc7XBSbklSTzR7dl53eskSmDwZxo7NSwyq\nubmCSwksFiVJPVVKcM89uWjs0ycvIXjIIRaNzcxisQQWi5Kknm71arj5ZrjgAth6a5g6Ffbfv+yo\ntD4sFktgsShJ2lCsXJnvY5w4MQ+MmTIlrwyj5uHa0JIkqcv06ZPXmZ4/H446Kt/HePTR8PTTZUem\nMlksSpKkj9loo7zO9IIFsO++cNBBcMIJsHBh2ZGpDBaLkiSppk02yetML1yYlw/cay847TR49dWy\nI1N3sliUJElt2nzzvM7088/nNah33x2+9z14442yI1N3sFiUJEkdMmBAXmd63rz8euhQOPdcePfd\ncuNS17JYlCRJ62SbbfI6008+Ce+8k5cQnDwZ3n+/7MjUFSwWJUnSetlhB/iP/4BHHoHnnoNBg+Cy\ny2DZsrIjU2eyWJQkSZ/I4MF5fsZZs2DOnPz6yithxYqyI1NnsFiUJEmdYvhwuOUWuO02uP32PIL6\nuutg1aqyI9Mn4Qou7XAFF0mS1s+DD+Z1p996CyZNgmOOgV52U3Upl/srgcWiJEnrLyX47W9h3Ljc\nwzhlCowZA9Gp5YxaWSyWwGJRkqRPLqV8efr886F//1w0Hnhg2VH1PBaLJbBYlCSp86xaBTNmwIQJ\nMHAgTJ0Ke+5ZdlQ9R1cUi945IEmSuk3v3nmd6Wefha9/HY49Fo48Ep56quzIVI/FoiRJ6nZ9++Z1\nphcsgIMOgkMPheOOg/nzy45M1SwWJUlSafr1y+tML1gAI0bAfvvBN74BL79cdmRqZbEoSZJKt9lm\neZ3pBQtgu+1gjz3gO9+B114rOzJZLEqSpIaxxRZ5nennnoONN4Zhw+Dss/Ma1CqHxaIkSWo4W2+d\n15l+6ilYsgR22QUmTszP1b0sFiVJUsPafnuYPh0eewxeeAEGDYJLL4WlS8uObMNhsShJkhrezjvD\n9ddDS0suHAcNgiuugOXLy46s57NYlCRJTWPoUJg5E+68E+65B3bdFa65BlauLDuynssVXNrhCi6S\nJDWuhx6C8ePzqOkLL4SvfQ16bcBdYRvkCi4R8S8RcWdEvBYRqyNiQo02fxsRl0TEExHxXkS8GRG/\ni4j967zn2KLtsoh4KSLGRUTDfy8kSdLH7bcf3H9/viR9+eUwciTccUdei1qdoxkKpFOBTwO3Fq9r\npf+/AV8r2hwLnAx8CLRExJcrG0bEocBM4FHgMGAaMB64qAtilyRJXSwCDj4YHn0UJk3KPY177QX3\n3WfR2Bma5jJ0RPQGPgImppQmVR3rD/wlpbSqqv084I2U0gEV+58A3kspHVix73xywbhjSumNqvf2\nMrQkSU1k9Wq48UaYMCFP8D11KuyzT9lRdY8N8jJ0hbpfeErpz5WFYrFvFTAX2PavbxCxAzAC+GXV\nW/wC6Asc3mnRSpKkUvTqBccfD888AyeemJ9/+cvwxBNlR9acmqlYXCcRsRGwN/Bsxe5hxePTlW1T\nSi8BS4HduiU4SZLU5fr0gW9+E55/Hg47DMaMga9+FZ59tv1ztUaPLRaBicB2wCUV+wYUj4trtF9c\ncVySJPUQG28M3/0uLFwIo0bBAQfAySfDiy+WHVlz6NZiMSK+VIxobm+b/Qk/5x+Ac4BJKaU5HT3t\nk3ymJElqbJ/6FJxzDixYAJ/9bC4czzgDFi0qO7LG1qebP28OMKQD7dZ7EZ+IOAK4Fvh5SunCqsOt\nPYpb1jh1C+DdWu85ceLEvz4fPXo0o0ePXt/wJElSyfr3z3Myfve7cMkl8LnPwSmnwA9+kNekbiYt\nLS20tLR06Wc002joPsAKaoyGrmhzEHA3cGtK6fgax3cEXgJOSyldXbF/IPACcEpK6bqqcxwNLUlS\nD/anP8FFF8GvfgVnnglnnQVbbFF2VOtnQx8N3aaI2Bu4HbgPOLFWm5TSK+QR0idUHTqRXIje25Ux\nSpKkxrPttnlS78cfz5ekBw+Giy+GDz4oO7LG0PA9ixExChhILmxnAP+72ADuTikti4ghwMPAn8kT\ncn9sWfGU0u8r3u9w4C7gquL9RpIn5P5xSumcGp9vz6IkSRuQ+fPzHI0PPJAvTZ9+OvTrV3ZUHdMV\nPYvNUCxeC5xUvEysGYiSgJ1SSq9ExMnANVXHW6WUUu+q9/wKMIF8/+TrwM+BqbWqQotFSZI2THPn\nwvnnw5NP5seTT4a+fcuOqm0bZLFYNotFSZI2bL//fV5C8OWXYeJEOO446N273dNKYbFYAotFSZIE\ncP/9MG4cLFkCkyfD2LF5XepGYrFYAotFSZLUKiW4557c09inD0yZAocc0jhFo8ViCSwWJUlStdWr\n4eab4YIL8tyMU6fC/vuXHZXFYiksFiVJUj0rV8INN+R7GXfdNfc0jhpVXjzOsyhJktRA+vSBk07K\n0+0cdVTejj4ann667Mg6j8WiJEnSJ7TRRnmd6YULYd994aCD4IQT8utmZ7EoSZLUSTbZJC8XuHAh\nDBkCe+0F3/oWvPpq2ZGtP4tFSZKkTrb55nki7+efh622gt13h+9/H954o+zI1p3FoiRJUhcZMCCv\nMz1vXp52Z+hQOO88WLy47Mg6zmJRkiSpi22zDUybBk88AW+9BYMH54m933+/7MjaZ7EoSZLUTXbc\nEa66Ch55BJ57DgYNgssug2XLyo6sPotFSZKkbjZ4cJ6fcdYsmDMnv77ySlixouzI1maxKEmSVJLh\nw+GWW+C22/I2ZAhcdx2sWlV2ZGu4gks7XMFFkiR1lwcfhHHj4O23YdIkOOYY6LUOXXsu91cCi0VJ\nktSdUoLf/AbGj89rUE+eDGPGQHSgBLRYLIHFoiRJKkNKcOutcMEF0L9/Xnf6wAPbPsdisQQWi5Ik\nqUyrVsGMGTBhAgwcCFOnwp571m5rsVgCi0VJktQIPvoIrr02X5YeOTKPpt5884+3sVgsgcWiJElq\nJB9+mEdQH3/82vcxWiyWwGJRkiQ1i64oFp1nUZIkSXVZLEqSJKkui0VJkiTVZbEoSZKkuiwWJUmS\nVJfFoiRJkuqyWJQkSVJdFouSJEmqy2JRkiRJdVksSpIkqS6LRUmSJNVlsShJkqS6LBYlSZJUl8Wi\nJEmS6rJYlCRJUl0NXyxGxL9ExJ0R8VpErI6ICR04Z5+i7eqIWOtrjIixEfFERCyLiJciYlytdpIk\nSRu6ZiiQTgU+DdxavE5tNY6IvsDPgNdrtY2IQ4GZwKPAYcA0YDxwUeeFLEmS1DP0KTuA9qSUhgJE\nRG/g2x045d/IReI1wHk1jv878H9SSq3v9UBEbAaMj4gfpZTe6ISwJUmSeoRm6FlsFe02iNgZGAec\nCayscXwHYATwy6pDvwD6Aod/8jDVSFpaWsoOQevJ3DU389fczJ8qNVOx2BHTgZtSSg/VOT6seHy6\ncmdK6SVgKbBb14WmMvgLr3mZu+Zm/pqb+VOlhr8M3VERcSKwB3B8G80GFI+LaxxbXHFckiRJdHPP\nYkR8qWKUclvb7HV83wHAZcC5KaW31ze89TxPkiSpx4qU2hxc3LkfFrEJsEMHmi5NKf1X1bl9gBXA\nxJTSpKpjPwX2AQ5kzQjoHwBnk0dSL08pfRARhwN3A3unlB6teo+/AP8zpXRO1f7u+wZJkiR9Qiml\nTu0A69bL0CmlZcDzXfDWuwGfB96pcext4DbgaGBesW84eeocACJiILAp8Ez1yZ39DZckSWomPeWe\nxe8D/av2nQKcBBwEvAGQUnolIuYCJwBXV7Q9kdxreW/XhypJktQ8Gr5YjIhRwEDW3F85LCKOLZ7f\nnVJallKaW+O8LxZPH0gpra44dB5wV0RMB2YAI8nT7UxLKb3ZFV+DJElSs2qGqXP+CbiJXNgl4KvF\n6xuBrds4L1FjBZeU0r3AscBewK+B7wFTyfc4Ank+xoiYGRHvRcSfI+LmYo5GNZCIODYibouIVyJi\naUQ8FxEXFZOsV7bbMiJ+HhFvRcRfIuK+iBheVtyqLSJ+XQxwm1y13/w1qIgYExEPRsT7xe/KP0TE\ngRXHzV0Dioj9i1y8GRFLIuLxiDilqo25K1lEbB8RP4mIR4q/casjYsca7TqUq4joFxE/LJZPXhoR\nD0fE/h2JpeGLxZTSKSmlXsXWu+r5K22cd2HRZnWNY7emlHZPKfVLKQ1MKU1JxUifiNgUmA3sAvx3\n4B+BwcD9xTE1jrOAj8iF/mHAlcAZwH0REQDF453AIcB3gGPIE7DfHxHblRG01hYRx5PvO4aKf/LM\nX+OKiNPJ94P/ARjLmn/kNy2Om7sGFBEjgfvIf/+/CXyFnMOrI+LbRRtz1xgGkX+u3gEerNVgHXN1\nNXkJ5fHAl4HXgN9ExIh2I0kpuVVs5J7GlcDfVewbSC5K/kfZ8bl9LFdb1dj3j8Bq4MDi9VHF6wMq\n2vxN8cM3reyvwS0BbFn80vp6katJFcfMXwNuxe/EZcA/t9HG3DXgBlwMfAhsWrX/YeBhc9c4G8WM\nNcXzU4uc7FjVpkO5Iq9etxo4qWJfb+A54Pb2Ymn4nsUSHAk8klJ6oXVHyiu8zCEnRQ0ipVRr9Pv/\nLR63LR6PBBallB6oOG8J+T8x89kYLgH+mFK6scYx89eYvkH+p3p6G23MXWPqTe78WFa1fwlr5hs2\ndw0gFRVdOzqaqyPJeb+xot0q8i1+h0ZE37Y+xGJxbcOoWg6w8AwwtJtj0bo7oHh8tnhsK587emtB\nuSJiP3Jv8D/VaWL+GtN+wHzgHyLi/0XERxGxICLOrGhj7hrT1cAq4McR8bcRsUVEnAZ8EfhR0cbc\nNY+O5moY8EJK6cMa7TYiX/Kuy2JxbVtSeznAd4tjalDF/RmTgPtSSv9Z7B5A/XyCOS1NRGwE/Az4\nYUppQZ1m5q8xbUu+l/tS4CLgYPJ9cFdExD8XbcxdA0opzQcOJd8Lt4icjyuA01NKNxXNzF3z6Giu\n2mvX5nLHDT91jtQRxQjo28nzZVaO6nMFnsZ1NrAxeTaCesxfY+oFbE6+/+m2Yl9LscDBucCPS4pL\n7ShGyd5FvmXnJ+TL0WOBn0XE8pTS/yozPq2zbvkdabG4tsXU/q9pAGsqcDWQYhnJO8k33R+QUvpT\nxeHF1P6PaUDFcXWzYvqHceTRmJsUOWzVLyL6A3/B/DWqd4Cdyb2Jle4DDouIbTB3jWoy8B5wREpp\nZbHv/ojYCpgWEb/C3DWTjuZqMbDWtDsV7dqsb7wMvbZ55OUAqw2lxnKAKldxU+5MYA9gTEppXlWT\neeR7NaoNBV5OKS3t4hBV29+RexV/Sf4l1boB/Cv5F9twzF+jmseawRBttTF3jWco8FRFodjqD8BW\nwGcwd82ko7maB+wUEf1qtFsBLGzrQywW13YHsFdE7NS6o7i0sk9xTA0iInoBNwCjgbEppcdqNLsD\n2C4ivlBx3t8AR2A+y/QEOW+VW+tkzr8oXi/E/DWqW4rHw6r2Hwa8mlJ6HXPXqP4LGFFj9Oue5EvS\n75Bv6TF3zaGjP2d3kOdf/FpFuz7kKct+k1L6qK0PiY6NzN5wFCOH5pJ/aMYXuycDnwI+739UjSMi\nrgROJ9/zdnfV4VdTSouKCUsfAnYA/o18+eVccq/ViJTSom4MWe2IiNXAlJTSBcVr89egImIWee62\nccCL5AET3wROTildb+4aU0QcBdwK/Bb4KXnOxSOBM4HLU0r/au4aR6xZ3vgg8t+7M4G3gTdTSg+u\nS66KWwwOLdq9RF7EYgywT0rpyTYDKXvSyUbcim/6TODP5LmnbqFqIky38jfyH6hV5IlGq7cLKtpt\nSZ4u4h3gA/J9VZ8rO363mjn92KTc5q9xN/IAlyuA14HlwJPAceau8Tfy6PXZwJvF37j/BL4N9DJ3\njbVV/V2r/Hs3e11zBfQDLiMvgrAMeAT4QkfisGdRkiRJdXnPoiRJkuqyWJQkSVJdFouSJEmqy2JR\nkiRJdVksSpIkqS6LRUmSJNVlsShJkqS6LBYlSZJUl8WiJEmS6vr/Nbvn7RNXE9gAAAAASUVORK5C\nYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x73dc330>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(10,6))\n",
"plt.plot(d, label='d Orbital')\n",
"plt.plot(s, label='s Orbital')\n",
"plt.ylabel('Energy')\n",
"plt.legend()\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Orbital Energy Difference\n",
"\n",
"Here we are calcuating the difference in energy between the *s* and *d* orbitals. Because of the way we've set up our *d* and *s* orbital energies, this difference is getting progressively larger. We'll use this as the basis of our comparison. "
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"ediff = s - d"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Interaction Energies\n",
"\n",
"Here are the electron-electron interaction energies, which the paper calls (4s,4s), (4s,3d), and (3d,3d). The 3*d* orbitals are smaller than the 4*s*, so remember that the ordering of the magnitude of these values is (3d,3d) > (4s,3d) > (4s,4s)."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"dd = 10.\n",
"ds = 7.\n",
"ss = 1."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## [Ar]4$s^{2}$ 3$d^{4}$ Configuration\n",
"\n",
"Here we'll first calculate the energies of the 3*d* and 4*s* orbitals based on equations 1 and 2 from the paper. In order to calculate the energy of this configuration, we must calculate an average energy for each configuration based on the total number of electrons in each orbital. In this example, there are 2 *s* electrons and 4 *d* electrons."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"d4_e = d + 3.*dd + 2.*ds\n",
"s2_e = s + 1.*ss + 4.*ds\n",
"\n",
"s2d4_ave = (2*s2_e + 4.*d4_e)/6."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## [Ar]4$s^1$ 3$d^5$ Configuration\n",
"\n",
"This is the same as above, but we just need to change some of the electron numbers."
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"d5_e = d + 4.*dd + 1.*ds\n",
"s1_e = s + 0*ss + 5.*ds\n",
"\n",
"s1d5_ave = (1.*s1_e + 5.*d5_e)/6."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Plot Configuration energy vs Orbital Energy Difference\n",
"\n",
"Below we're making a plot of the *average* configuration energy versus the energy difference of the *s* and *d* orbitals. Notice that once the orbital energies reach a certain energy gap, the [Ar]4$s^1$ 3$d^5$ configuration becomes lower in energy. (You can ignore all of the plotting code.)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x74a4b50>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAnwAAAGuCAYAAAAK89vvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmcjWX/wPHP9yDL2OMRYUhEJIlHZRtbPNosT3tk6VdJ\nspS02iJJixbtmKIeIRUlkgZjSSqyPDwmmRTJTnZmvr8/rnvGcZwxZ/Yz4/t+vc7rdK77uu/ruu85\nzXxdq6gqxhhjjDEm7/LldAWMMcYYY0zWsoDPGGOMMSaPs4DPGGOMMSaPs4DPGGOMMSaPs4DPGGOM\nMSaPs4DPGGOMMSaPs4DPGGPOISKSKCIxOV2PtBKRWiLyhYjs8O5hZU7XyZjcJH9OV8AYE5yIVAF+\nTSXbQlVtkfW1CX8icj7wIHAdUB0oAuwBfgKmAR+q6omcq2FYyVULsIpIPuBTIBJ4H9gGbM/RShmT\ny1jAZ0z42wh8lMKx+GysR9gSkWuBKUBJYC3wIbAX+AcQBUwA7gJa51AVw0lN4HBOVyKNLgJqAG+p\n6gM5XRljciML+IwJfxtVdXhOVyJcicgVwOfASeBWVZ0WJE974L7srls4UtWNOV2HdCjvvf+Vo7Uw\nJhezMXzG5BEiUsUb2zRRRC4WkU9FZK+IHBSReSJSN4Xzqnnn/CEix0TkdxEZJyJlznL92iLyuYjs\n8dKK+13rExHZ773misjlIhLt5avs5evufX4khTrd4R1/LIRbfxUoCPQOFuwBqOps4OaAMoqKyEgR\niRORo97YsOnBnpOILPDqU1BERnvP6pCILBaRhl6eyiIyRUR2e898hohccJZnWFdE5vg9q89EpEaQ\nslt6+Td61z0gIktE5JYgeUP5GZ0xhk9EKojI6yLyi4gc9s5ZIyKvikj+gLxVReR9EfnT+77Ee/lO\n+774lyUi5bxzdnrXXyYizYP9rIJcIx5Y4H0c4l0zUUTu9stzufe92+n9LDeKyAgRiUjr8zEmrxLb\nS9eY8CSnxvB9oao3piH/QqAOrmvzB+Bi4CZcF2ctVd3hd87VwBzgPGAmsBmoBVzv/XdDVd0bcP0l\nQF1gJfA9UA54ACgNrMB1o34G/BeoD1wDrAKaAVVUdYuIFMaNw/pTVS8Nci/zgaZAZVVNcayWiFQH\n/gdsUdUqqT0jv/MKA7Fe/b4DYoDKwC24lsJ2qhrrl3+BV/+Z3vP5CjgfuA045B37CvjFeyb1gFbA\nIlWN8rtOFdwzjAUu98r+AdfN2hE35vAqVd3kd85XQFVgObDVK/dG3HPvr6qvBLl+0J+Rqh4SkURg\ngaq29M6JANbhWtG+8J5nBK4LNQo4X1UPe3lreXUvhRtT9z+gIa6rfLNX951+9UkEfvaut887txxw\nK5AAXKmq64L8iJKJSF/ved6NC/wWeIc+U9XVXuD4FSDAVO8ZtQT+CfwINFXVoyE8n96qevBsdTEm\nV1NVe9nLXmH4AqoAibg/qkNTeDUKkj8RGBhwreFe+iC/tPOA34BdwCUB+f/t5X8thes/FaS+H3rH\negekP+2lJ+ACuKT01730qwPyV/XSPwvhGd3t5Y1O47Md6p33bkB6Sy89Du8fxF76Ai/9W6CgX/oA\nL30/8FzAtT73jtVP4RkODcjfw0ufGZAeGaT+RXBB9D6gcKg/Iy9PIvCt3+cbvbQ+QfKWSOE53BWQ\nb5iXPiFIWad9jwLu9c0Qf15RXv7BAen5gE24IL1pwLGJgeeE8nzsZa+8+srxCtjLXvYK/gr44xTs\nlQA8FCT/L0GuFekdm+aX1slLeySF8n8AdgS5/h9AvoC8BYGjwBbAF3CsMLCbMwO+ugQPup7x0m8I\n4RkN8vKOTOOz3QwcAf4R5Ngs75pN/dKSAp3GAXkv5FTAVyjg2J3esbuDPMNd+AVqfsc3eMFL6RDu\nISnYbB7Kz8gvT0oB3/+lUl7Sd+inIMcKATtwrZ0FAso6ABQJyJ8POAGsCPHnFUXwgK+5l/5pkHPK\ne9/JTWl5PvayV1592aQNY8JfSF26flYFSdvqvZf0S2vkvdcRkaFBzikMlBGR0qq6xy99taomBOSt\niWsx/EFVE/0PqOoREVkFtAhIXy0i3wO3ikg/dd2NPqAb8CfwZYp3mAHeWK1I7z52BMmyELe0S11c\nF2RylXHdk/6Supvj1Os29JM0waBCkDJWquqRIOlLcV2pdYBFfvV9FNctfxHu5+LvAs4U7GeUkgW4\n+xgnIq1xXfyxqvpLQL7LvfdFgRdQ1aPez7I9cAluOEGSjep1CfvlTxCRvzj9+5geSXVaGKROf4rI\nL0AtEYlQ1UN+h9PyfIzJEyzgMybvORCYoKonRQRcy0qS0t5717NcS3Hjr/wDvmAzJYt57zuDHAPX\n+hPMO8B7uLFzE4G2uFaz5wIDxxT86b1fGELeJEmD81Oa8bk9IF8yDRjj5QUuEOSZ41rqAAoEOZbS\n80iqUwkAETkPF8xcjmtxnYD7WSQAV+CCwIJnuU6qVPWAN5ZzOHAD3uQWEdkIDFPV/3hZQ31uxQLS\ngz0bcM8nXwrHQhVKnS718vkHfDbb15xzbJauMeeupD/E7VTVl8Irn6r+HnBesJlef3vvZVMo6x8p\npE/x6tHD+9zTu/74EO9hqffeXLzIKwRJ910uhePlAvJlhZSeR1LZ+733m3DB3juq+k9V7aOqQ9Qt\n0/PdWa6fptl4qvqbqt4NlAEaAE/hJmZ8KCJJLbPh8NwChVIn5cw62WxFc86xgM+Yc9dy7/3qTLjW\neuA40EDcrgjJvBmx9QjyR9br6vsIaCwiTXHjyRaq3yzVs/G6HRfjumjvOlter7UMVT2AW7D6EhEJ\nFnglLRcSrGs8s1whIkUC6ie4Gc0JwBovuZr3PjPINRpndqVUNVFVf1LVZ4F7veTrvfek59E08DwR\nKYSbFXsEN8kouyRtr3bGEi8iUh7XPf5rQHeuMeckC/iMOXd9BvwODBSRRoEHRaRwsPRgVPU4MAOo\nyJkLHA/AtRal5B3v/T+4YSahtu4leQg3OP8NEekULIOIXA9M90t6H9cV+kxAvihcgPOLqi5JYz3S\nojRuwom/HrgAZbZ6S+HgZlEDNAmoZydOBWIZIm6P2mCBb9LYwKMAqroF1718hYjcHpB3EK51d4qq\nniT7LMYts3KjiAQGwCNx3ekfZGN9jAlbNobPmPB3SQqTKgCOqOro9FxUVY+LyM24NcyWisjXuJa6\n/LjZjM1x65W1D/GSj+HWnnvNG/y/Adeydw1u8kNT3AzJwHqsEpEfcF2J+zg9MAvlPlaJyE3Ax8B0\nEVmDCwT24YKQ5ri9def5nTYaFzD9n4gkTZCoiBtLeJhTXcz+Qu0yDkUs0FdErsKtFXcJbh2+3UB/\nv3yzcDOfB4lIbVzrWW3cWMdPvXMyqi3wvIjE4tYR3IMLPK/Hrd34nl/eXrhnO8kLOuOAK4E2uMAr\nMIjNUqqqItIDN9HkGxGZilvfMQo3KekH4PnsrJMx4coCPmPC38XA4BSO7cMFL+miqt+LSD3cLNB/\n4WbSHsItWxFNGlpH1C2ofA0wBhf4tQKW4QK9Z71sf6dw+iRcwPeRqh5Lx33M8xZhfhAXoN6Bm2yy\nG9ftN8orIyn/UW9s2uO4IG8AbpzXLNxEhTWnl4CSueO+NgG9cc+qt3ftmbj1E3/1q+dBEWkJvIBr\n5Wvh3c+/cLN/O2RCXebgusSbAZ1xa/z9gWtpHa2qSa2MqOoGcTuLDAWuxXXB/4lbU3G4qu5KQ7mZ\n8jxVdZE36WQwbnZ1MVzL6LPAs+n5PhmTF9lOG8aYLOWN6duEW6cu2BIiiMgE3HIs9VU1K8fO5Si/\nnR6iVTVYK6IxxmQJG8NnjMk03kD5QANx25Z9fpZzbge+z8vBnjHG5CTr0jXGZKb1IrIcN34P3MzN\nRrj10Ib6ZxSR63Djv27GLdo8PPuqaYwx5xYL+Iwxmekt3DiqRrixYH/iBv0PV9U/A/L+G7cX7u9A\nf1WdnZ0VNcaYc4mN4TPGGGOMyeNsDJ8xxhhjTB5nAZ8xxhhjTB5nAZ8xxoQhEblbRFaLyCERSRSR\nh7z0RBGJyen6ZSURiReRzTldD2PyEgv4jMlB3h/vUF9h8wdQRCqJyBgR+VlE9ovIURHZIiLTRaSz\nty9sriIibUXkY+8+jojI317A9bKI1MrmujQBJuK2BnsFN8P5O78suXrwtYhEe9/pyilkyeyFro05\n59ksXWNy1jBO/8MmwBDcDhpjA/Luy65KnY2IdMHtf1sQWI4LTA4CFwItgU7ABOCenKpjWohIEdyu\nIv/G3cfXuC3D8gGXAvcCfUSklaouzKZq/ct7v1tVvw84VhO3/Vtud7aArmW21cKYc4QFfMbkIFUd\nFpgmIkOAfaoaduvSiUh74H1gJ3BLYADktezdidsCLLeYiAv2vgK6qupu/4MiUgZ4BiiejXVKWsD6\nr8ADqroxG+uRlYQU9idW1bBpzTYmr7AuXWPCnIh08Lq/7gtIH+Glxwak1/HSXwpIv1xEPhGRnV4X\n7EbvGhEh1iM/bs9UJUiwB24ze1WdDPTyO2+oV5/mItJDRFaKyGER+TSEMn0i8n8islxEDnrdrEtE\npGModQ7h+q1wCz9vADoGBnvePe1S1V7A3IBzm4vIXBHZ693PGhEZ6D0n/3xR3v0PEZEGIjLPu499\nIjJDRCID8+K2mQPYnNSl75cn6Bg+EWklIku9uvwlIhNE5Pxg4+FEZIH/NQOOndHdKiLdvLS7ReQm\nr5yDIrLSO15CRB4TkVgR2S4ix0TkNxF5S0TKBVw/HugaeH/+9xSszl56VRF5X0T+9MqIF5FXvaA8\nMG+iiMSISDnvnJ3es1kmIs2D5K8gIq+LyC9evj3ez/TVwJ+pMbmRfYmNCX8LcUFWFPC2X3pSK1oD\nESmsqke8z1He+4KkjN4fuK9wLSpTga24brMngLYi0lRVj6ZSjxZAFWBxal2bqno8SPIgoCkw06vL\n32e7htda+B9cQLYO1xIHcD3wiYj0U9VXU6lzapL2s30hhTon8z8uIrcBH+Lu4WNgr1ev0bh7vDHI\nJRoCjwLfAm8C9YEOwGUiUkdVjwGbcd38HYDLcd36wbryT+sO9VpeZwJHvHrtAtrjuqfzA8Hu7Wxd\nqikduxVojdsmbyGn/obUwg1FmI/7fh0B6uG6w68VkfqqmnQfL+MC2sD7i0/lHmsBsUAp4FPgf7hn\n+iBwnYhcpao7A65R0jtnH/ABUM67h7kicqWqrvOuHQEsxbWsfuHdQwRQA/g/4DHgZArPxJjcQVXt\nZS97hdELSAR+DUhbCfzp9zkCOAZ84+Vv7XdsOu6PU0nvcz5gk5fWNOC6E73zB4dQryFe3qFpvJ+h\n3nn7gFppOO8+77zX8BaJ99KL4CYwHAXKZ/BZbwYSgCppOKcEsB8X7F3il54fmOfVuatfepSXlgjc\nHHCt9730WwPSo730yil8P771+5wP2OJ9H+r5pftwrZLBvk8LgIQU7u+MsnEBWiJwAmgW5JziSd+3\ngPQ7vfOeDPX+vOPxKdQ5EbgrIH2Ylz4hyHNKBF4LSO/hpb/pl3ajl9YnhZ+3BKunveyVm17WpWtM\n7rAQKCenZos2xs3gfA4X+LSA5FaxZsBqPdWi0gSoCsxS1dO6f3EtfMdxW5yl5gLvfWs67+EdVV2f\nhvy9cUFif1VNbu1R1cPACNz+u53SWZck6bmnDkAx4D1V/Z9fvU7iWoIg+PNcqKrTAtImeO8N0lB+\noKZARWCGqq7yq08i8HQGrhvMp6q6KDBRVQ/4fd/8fYQLjFtlpFCv27sZsErdkAF/o3AtmreKSIGA\nYwdxLcv+3scF+cGe+Rmt3Kq63//7Z0xuZV26xuQOMcBDuNai9d77YVwguIxT3bu1gTKA/x/Fy733\nYGPu/hSRX4BaIhKhqoeyovKeH0LNKG7mbB1cy9VTcuYqL2W990u8/PVwgZi/zar6fvqqelZne54/\nisghoG6Q834MkpYUaJbMhPosC3LsB1xwk1lS/Bl64yH7Af8ESuNaHpOUD3pS6JLuMViweVREvsd1\nYV8CrPU7vNH7B4J//gQR+YvTn/kCYDswTkRaA3OAWFX9JYP1NiZsWMBnTO4QixvT1AI3/isKWKaq\nJ0RkIS4oiiDI+D1OzS49Y8anZztu+ZHiwNkCvu3ee8U01j1JSuUHU8p7rwwMTiGP4rq2wY0XC8y3\nANeaczbbvTIq4rp3Q5Ha8/wLqBQk/UCQtKRxYfmCHAtVMe89cPwaqpooIrsycO1AQe9ZRG7Fjbfc\njwuW4nHj+AQXBBbMYLmhfIfh1LNIEuyZg3vuyc9cVQ+IyNXAcOAG3LhRRGQjMExV/5OeShsTTqxL\n15hcQFX3AGuA5iJSFDdYfYF3eAHuH29NODVezL8lJOmP3mmzJf2UwwVPKf1xTLLEe0/vGmlp6RZL\nqst3qupL4ZVPVXsCqGp0kOOh1HMJLihJyz2F8jxTe5aZKWnyS9nAAyLiA84Pck6i3/FAJc5SVko/\nw8G4FucrVfV2VX1cVYerW3Yoo8EehPbM/fOlmar+pqp341rIGwBP4f7h8aGI5KZlhowJygI+Y3KP\nGNwf9ftxrRNJy1h8h2tNacmZ4/fATfgACLYURXncTMRfQ+jOjcG1gjUWkWZnyygi56VyrbNS1b9x\nS6XU9gLcrJI0hm6AiJw1MPG7p7M9z/q4VsdVgceyUFJZVwc51oDgPTl7cYHuaa21XgB4OWnf5aIa\n8F9V/TXgevWBQkHyJ3Uzh9qymfTMmwYeEJFCuG7kI7iZuxmiqomq+pOqPoubZQxuBrYxuZoFfMbk\nHgu894G4rtfvIXm5kGW4mZRlOL07F2Ax8Ctwo4g0Djg2Ejf544PUClfVBKAPLhiYKiLB/vj6RKQr\n8FYoN5SK13BddG8FC8ZEpLaInNGqlRaq+i1uCY5awHQRKR2knLIi8g7Qzkv6HNeSdI+IVPfLlx83\niQZCeJ6ZKBb4A+gsIsljB0UkH66LMpik3Tu6BaQPwC29k1a/ATX8fx4iUhy3LVwwe7z3lLZWO42q\n/o4bM3mFiNwecHgQ7h9CU7yJM2kmIrVE5B9BDiVN6kltySJjwl6uHMMnInOAa4GRqvq0X3opYAxw\nE1AY90ewv6quDXohY3KXpG7assA3AX/cFnCqW3KB/0mqqiLSAze26hsRmQpsw3X/NsINxH8+lAqo\n6mwRuRu3tdpCEfkOWIGbDVkBNxuzIvBuGu8tWFlvisg1uKU9mnoL827HTQC4DDdu7yqCjF1Lo+64\nf/z+G4gXka+BX3CtT7VwzzU/MMmr1wERuR83MWaFiHyMm018HW4s5CxVnZTBOoXMG6fXC/gMWCIi\n/8EFVP/CjVXbhteF62cibk3Aod6El19xrYG1cYHVGa2XqXgdF9z9JCIzcN247XCTUrZx5o4a84GH\ngXe8/IeA+CAzcP31wv3jZZKIdMJtf3cl0Marf+Bs3LRoCzwvbhHzX3DPrwauZW8v8F4Grm1MeMjp\ndWHS+gJu59QvsOF+6YL7ZbAFt7BmW9wfvp3AhTldb3vZK9QXQdZN8zu2Etcd9kRAemPvvOT194Kc\neznwCW4Ji2PARtzyJkXSUcdKuCBxFW6g/jHgd2AG0Ckg7xCvzmes3xZiWXfgFiveg2tpicct3Hxf\neup+lnLaAlNwrVVHcWPjVuMWB74kSP7muCB6L647cQ2u9TVfQL4oUljrENeaFmwNuYneM0t1HT6/\n9Na4xYMP4yY3jMfNlv0bWBkkf13cuoEHvXuYAVwUrGzcMjMJ+K0vGOR6D+AWyD7sPcOXgKK4YQBn\nfJ+BR3BdsMcC7+ks51TFTcT50zsvHngVKBPqcwp2fdz+xC/jZlLv8u5hI/AGEJlZ3zF72SsnX6Ka\ne5YX8lrw/oub9fUfYISqDvaO3YRbfb2FersAeF0Km4HJqto3Z2ptjDE5Q0Sq4VrCpqrqbTldH2NM\nzsltY/hGA2tU9eMgx24Etqrflk+qegCYheviNcaYPElEIgInt3jjHl/0Pn6e/bUyxoSTXDOGT0Sa\nAF0IvqApuLEnwcbq/RfoKiJFNGABTmOMySNqATEi8hWum7MUbuxhVdzQlik5VjNjTFjIFQGftxzC\n28AYVY1LIVtp3MDdQEmzwUrhxmUYY0xeswU3PrMZbvKI4H4fDgVGa24au2OMyRK5IuDDzSYriFtC\nIiX2C80Yc05S1R2cucSKMcYkC/uAT0QqA08CPYHCIlLY73AhESnBqVlmZ6yh5Ze2N+C6FiAaY4wx\nJtdQ1TM2Fg9Vbpi0cRGudW8yrns26QVuWv9e3Cbr63Dj+AJdCvwWbPxeTk+RDsfXkCFDcrwO4fay\nZ2LPxZ6LPRd7JvZccvqVUWHfwodbdywqIE1w2zxNwq019QswE+guIs1UdREkL8tyAy5YNMYYY4w5\nJ4V9wKeq+zl9I3gARARcy11ScDcTt7PGZBEZiFv5/nHc2L6QdhEwxhhjjMmLckOXbkjUtXdej1s5\n/g3cqvEncAsxb83JuuUmUVFROV2FsGPPJDh7LsHZcwnOnsuZ7JkEZ88la+SqnTYyk4jouXrvxhhj\njMldRATN45M2jDHGGGNMBljAZ4wxxhiTx1nAZ4wxxhiTx1nAZ4wxxhiTx4X9six53ejFo2lfvT2X\nlbssp6tijDF5jreElzFhLTsmkdos3Ry0ZMsSmkxsgk983HflfQyLGkbZiLI5WidjjMlLvJmNOV0N\nY1IU6nfUZunmYjXL1OTBhg8iCG/+8CbVX6vOy8te5njC8ZyumjHGGGPyEGvhCwPrdqxjwNcD+HrT\n1wDUOL8GL177ItdVv866I4wxJgOshc+Eu+xq4bOAL0yoKrPjZjPg6wFs3L0RgGurXcvLbV/m0rKX\n5nDtjDEmd7KAz4Q7C/iyWLgFfEmOJxxn3PfjGLZwGPuP7Sef5KNXg14MjRrK+UXOz+nqGWNMrmIB\nnwl3FvBlsXAN+JLsPLSTIQuG8PaPb5OoiZQqVIqhUUPp1aAXBfIVyOnqGWNMrmABnwl3FvBlsXAP\n+JKs+WsN/ef2Z/7m+QDUKlOLl9q+RLuL2+VwzYwxJvxZwGfCnQV8WSy3BHzgxvfN/N9MHpn3CL/s\n+QWA9tXb89K1L3FJmUtyuHbGGBO+LOAz4c6WZTHJRISbat7E2l5rGdNmDMULFmd23GzqvFmH/nP6\ns/fI3pyuojHGGGOAxMREmjdvntPVOIO18OVCfx38i6djnua9n95DUc4vfD7DWwzn3ivvJb/PNk8x\nxpgk1sKXt82aNYu4uDi2b9/O1q1beeWVVyhTpkyO1ik6OpoePXqQmJgYUn5r4TMpKle0HO/c8A4/\n3fcTzSObs/vIbnrP7k29t+oxb9O8nK6eMcYYk+Xi4+OJi4tjwIABPP/885QrV45u3brlaJ3Wr19P\nuXLlcrQOKbGALxerd0E9Yu6OYfrN06lSsgrrdq7j2snXcuN/biRud1xOV88YY0wOio6OxufzJb+6\nd++epeXFx8efVl7VqlWztLzVq1fzxBNPcOzYMQBatGhBTExMlpZ5NidPnuTrr7/mX//6V47V4Wws\n4MvlRITOl3Zmfe/1jGo1iqLnFWXWxlnUfqM2j3z9CPuO7svpKhpjjMlBHTp0YOjQoXTs2DHFPJ07\nd8bn81GrVq10l1OqVCmGDh3K0KFDKVGiREg7RfXu3Ts5QNyzZ0+ayvvXv/5FbGwsBQsWBGDLli3U\nqFEj1fOOHTtG//79adKkCeXLl6dQoUJUqFCB5s2b8+GHH4bUFXv48GFat25NREQEkZGRAHzwwQc5\n3sJ4NjaGL4dt2QKVKkFm7aD2599/8uS3TxK9KhpFKVukLM+0eIZ76t9DPl++zCnEGGNyiXN5DF/S\nWLLo6Gi6du2aYr5du3ZRoUIFTp48CcDixYu55pprMlR2lSpV8Pl8/PrrrynmiYmJoXXr1hQpUoTD\nhw+zc+dOSpcuna7yjh8/ToMGDXjjjTdo0qTJWfPu2rWLyMhIrrrqKqpXr06ZMmXYsWMHs2bN4q+/\n/uKWW25hypQpIZUbFRVFxYoVGThwILt376Zly5YA+Hw+G8NnTtm3Dy6/HJo3hwULMuea5YuVZ8JN\nE1jxfytoUrkJOw/v5P4v76f+O/WJ2ZxzTd3GGGPC06RJkzh58iR9+/YFYMKECVle5qFDh+jZsycd\nO3akQYMGGQ7KH3vsMUaNGpVqsAdw/vnns3//fubPn89bb73FiBEjeOedd/jll1+49NJLmTp1Kj//\n/HOq1zlx4gQ//PADzZo1IyYmhhUrVjB69GiGDRsGwPPPP88vv/ySofvKVKp6Tr7creesRYtUS5VS\nBfeKinJpmSUxMVE/XvuxRr4cqQxFGYp2nNJRN+3ZlHmFGGNMGAuH3/U5ZeLEiSoi+v77758132WX\nXabFihXTQ4cOac2aNZP/O1BMTIyKiA4dOlRjY2O1VatWWrx4cS1VqtQZeSMjI7Vq1aopltm7d28t\nXbq0bt++XZs3b64+n093796d9ptU1VdeeUVjY2NVVTUuLi5d10jSv39/FRGdM2dOqnmXLVumIqLr\n168/LX3z5s3q9SKGJNTvqJcv3XGPtfDloKZNIT4ehg+HEiVcK1+zZtCmDSxdmvHriwi31L6F9b3X\nM6LFCCIKRPDphk+pNa4Wg+YN4sCxAxkvxBhjTK61YsUK1q5dS6dOnShSpAhdu3bl4MGDTJ06NcVz\nlixZQqtWrShYsCC9evU669jAYBYuXMibb77Jiy++mOEZrR9++CFVq1alevXqbN++nRkzZqT7WkeP\nHuXbb78lIiKCK6+8MtX8sbGxlC1blpo1ayanrVixgiFDhiAi9OvXj7i48JlAaWP4wsS+fTB2LLz8\nMhzw4rC2bWHYMGjUKHPK2Pb3Np6Y/wTv//w+AOUiyjGy5Ui61etm4/uMMXmSjeE7+xi+Xr168fbb\nbzNv3jxatWrF77//TpUqVWjcuDGLFi06Le+CBQuSx6h98MEH3HXXXSmWndIYvsOHD1O3bl2qVavG\n3LlzATcOLjY2NsUxfCtXruSNN96gRIkSFC5cmOrVq9O1a1eWLl1K8+bNSUhISM57/fXXM3PmzJCe\nz5EjRxgrqD3PAAAgAElEQVQ9ejSqyo4dO5g9ezbHjh1j/PjxXHfddWfknzZtGjNmzKBSpUqULVuW\nxYsXU6BAAaZPnx5SeSnJrjF8Od61mlMvwrSZf/du1SefVC1a9FRX73XXqf7wQ+aV8f0f3+s1469J\n7ua94q0rdGH8wswrwBhjwkS4/q7PDql16R45ckRLliypFStWPC29ZcuWKiJndI8mdek2bNgw1bJT\n6tLt06ePFitWTOPj45PTztal+/vvv+vVV1+d3MU8e/ZsjYiI0ISEhFTrkJqdO3eqiKjP51MR0Xz5\n8mmfPn10x44dZ+R97rnntHr16sn1mDZtmp533nn6yiuvZLgeoX5HsS7dvKV0aRgxAjZvhsceg4gI\n+PJLaNAAOnSAVasyXkbDCxuyuPtiPur0EZWKV2Ll9pU0j27OzdNuZvPezRkvwBhjcimR7H/llE8+\n+YT9+/dzxx13nJbepUsXIOXJGw0aNEhXebGxsYwbN44RI0YkL2USSh2LFStGkSJFANdyOHbsWHy+\njIcvZcqUITExkZMnTxIfH8+oUaN49913adSoEfv370/O99133/HEE0/w7rvvJtejTp06nDhxgqZN\nm2a4HtnFAr4wVaYMjBrlAr+BA6FwYfj8c7jiCujcGdasydj1RYTbL7udDQ9uYFjUMArnL8z0/06n\n1rhaPDH/Cf4+9nfm3IgxxpiwlBTQBXbNdu7cmcKFC/PBBx8EXVokPePuTp48SY8ePbjqqqvo06dP\n0DwapFuzdOnSzJs3j2uuuYbBgwdTsGBB7rnnnjSXfzYiQuXKlRk4cCAjR44kPj6e119/Pfn4008/\nTbVq1U7bH3fp0qUUL16cevXqZWpdspIFfGGubFl4/nkX+A0YAIUKwYwZULcu3HILrFuXsesXKVCE\nwc0Hs7HPRu687E6OJRxj1OJR1Hi9BtGroknU0NYRMsaYvODUYJrse+WE+Pj45F0pLr/88tN2yChR\nogRHjhxh27ZtzJkz54xzQ1lQOdDBgwfZtGkTy5YtI1++fKeVt2jRIlSVsmXL4vP52LJlS/J5d955\nJw8//DBxcXGMGDGCunXrsmLFivTfeCpat24NwNq1awHYs2cPMTExtG3b9rR8ixYtonHjxul6Fjkl\nf05XwISmXDl48UV45BF47jl4+22YNg2mT4fbboPBg8FvolCaVSxekcmdJvPgPx+k75y+fL/1e7p/\n3p1xK8Yxtu1YGldunHk3Y4wxJkdNnDgRcBMmLr744jOO79mzhxkzZjBhwgTat2+f4fIKFSpEz549\ngwZIX3zxBdu3b6dLly4ULFiQokWLJh/z+XyMGTOGMWPGsHz5crp06UJ0dDQNGzbMcJ2C2bZtGwDF\nihUDYNOmTSQmJvLPf/7ztHyLFi3ivvvuA1zwXKVKlSypT6bKyADA3Pwilw/k/f131QceUD3vPPdv\nRJ9P9a67VDduzPi1ExITdNLPk7TCixWSJ3bcNv02/W3fbxm/uDHGZKPc/rs+I1KatJGQkKCRkZGa\nP39+3bp1a9BzExIStHz58lqwYEHdtWuXqp6atDFs2LBUy05tHT5/KU3auO6667R9+/anpQ0cOFDH\njBkT0nVTsn79+qDrDO7Zs0fr16+vIqKzZs1SVdUNGzaoiGhMTExyvo0bN6qI6NKlS3XXrl06cuTI\nDNUn1O8oNmnj3FSxIowbB3FxcN994PPB5Mmula9bN9i0Kf3X9omPu+rexcYHN/J0s6cplL8QU9ZO\n4ZLXL2FwzGAOHT+UafdhjDEme82fP58tW7bQpk0bKlSoEDSPz+eja9euHD9+nEmTJmV5nTRI3/a6\ndeuIiopK/rx582a+++67DI/hmzJlCuXLl+f666+nd+/eDBo0iNtvv53IyEhWrlxJjx49uP766wGo\nUaMGderUYfNmN6Fx3759DBw4EBGhWrVqTJs2jU6dOmWoPtnFAr5crnJleOstF/jdc4+b8fX++3DJ\nJe5zfHz6rx1xXgTDWwxnQ+8N3Fr7Vo6ePMozi57hktcvYfLqyTa+zxhjcqEJEyYgInTr1u2s+bp3\n7w649fyykogE7eqNjo7m2LFjPP744/Tp04cXX3yRyZMnU7JkyQyVd8MNN3DrrbcSHx/PRx99xMsv\nv8zChQtp3rw5M2bM4L333jutblOnTmXGjBk8/PDDjBgxgvHjx9O/f3/69evH9u3bT1t4OZzZwst5\nzKZNblmXSZMgIQHy54cePeDJJ11wmBFLtiyh75y+/PjnjwA0urARY9uN5aqKV2VCzY0xJvPZwss9\nmDhxInfffXe2lp3SwsvmTNm18LK18OUx1arBxImwfj106QKJifDOO3DxxdC7N/zxR/qv3bhyY77/\nv++ZeNNELih6Acu3Lufq8Vdz14y7+ONABi5sjDEmy3Tv3h2fz5fcYpdV4uPjk2fe+s+0NeHBWvjy\nuA0b3F69U6a46f8FC8K997pFnVMYuhGSv4/9zajFo3hp2UscSzhG4fyFGdR4EAMbD6RIgSKZdwPG\nGJMB53IL388//8xnn32W3F1ar149brzxxiwrb//+/YwdOza5vJIlS/LQQw9lWXl5RXa18FnAd45Y\nt84Ffkn7YRcqBPffD4MGwQUXpP+6m/du5tFvHmX6f91egpWKV2J069HcVue2XLU+kTEmbzqXAz6T\nO1jAl8XOtYAvyerVMGyYW7wZ3A4evXvDo4+6RZ7Ta2H8QvrN7ceq7W7vt6srXs0r7V6h4YVZs1aS\nMcaEwgI+E+4s4Mti52rAl2TVKhg61G3XBm7P3j593MLO55+fvmsmJCYQvSqaJ759gh2HdgDQ9fKu\njGo1igrFMtB/bIwx6WQBnwl3FvBlsXM94Evy448u8PviC/e5aFHo29dt41a6dPqueeDYAUYuGsnY\n5WM5nnCciAIRPN7kcQZcPYDCBQpnWt2NMSY1FvCZcGcBXxazgO9033/vAr+vvnKfixeH/v2hXz9I\n75JHm/ZsYuC8gXy64VMAIktE8nyb57n50pttfJ8xJltYwGfCnQV8WcwCvuCWLYMhQ2DePPe5ZEnX\n2te3rwsC0yNmcwz95vZj9V+rAWhauSlj242lfvn6mVRrY4wJzgI+E+4s4MtiFvCd3eLFLvD79lv3\nuVQpN76vTx/w9pROk4TEBMavHM9T3z7FzsM7EYRu9brxbKtnuaBoBqYJG2PMWVjAZ8JdWAV8InIv\n8KGq5plNVC3gC82CBTB4MMTGus/nnw8DB7qZvUWLpv16+4/u55lFz/Dq8lc5kXiCoucV5cmmT9Lv\nqn4Uyl8oU+tujDEW8JlwF24BXyLwN/Ah8Laq/pzeAsOFBXyhU4X5812L39KlLq1sWbeGX69eUCQd\n6yzH7Y7j4a8fZtbGWQBULVmVF659gY41O9r4PmNMprGAz4S7cNtarRrwBtAJWCkiy0Skm4hkeZOM\niPxbRD4TkS0iclhENojIsyJSNCBfKRF5T0R2ishBEZknInWyun7nAhFo3dp1886ZA40awc6drov3\nootg7Fg4ciRt16x+fnVm3j6Tr+/6mtpla7N532Y6T+1Myw9aJq/lZ4wxxpjMkaYxfCJSAOgA3A9E\nAfuASbhWv/VZUkGRZcAfwKfe+xXAUGADcI2qqrgmoVigMjDQq9fjQG2gnqpuDXJda+FLJ1U3m3fI\nEPjhB5dWvjw88QTcc4/bxSMtTiae5J0f32FwzGB2H9mNINxT/x5GtBzBPyL+kfk3YIw5Z1gLnwl3\nYdWlm0LB1YH3gKZeUiwwRlW/SG9lUijnfFXdHZDWBXgfaKWqMSJyEy4gbKGqC708xYHNwGRV7Rvk\nuhbwZZCqW79vyBBYudKlVazoAr+ePeG889J2vb1H9jJ84XBeX/E6JxNPUrxgcZ5q+hQPNXqIgvkL\nZv4NGGPyPAv4TLgLty5d/wKLi0hv4BNcsLcSeArID8wUkWfSW5lgAoM9j9euRNL2DTcCW5OCPe+8\nA8As4KbMrI85RQRuuMEt3vzpp1C3LvzxBzzwAFSvDu++CydOhH69UoVL8XK7l1nTaw3tq7fnwLED\nPPrNo9R+ozafb/jcfmkbY4wx6RRywCciDUXkPWAb8AKwCteleqWqPquqjYEhwANZU9XTNPfek7qR\nawNrg+T7L1BZRNIxrcCESgQ6dHCtfNOnQ+3asGUL3Hsv1KgBEyakLfCrWaYmX97xJV/d+RU1y9Rk\n095NdPi4A20mtWHNX2uy7kaMMcaYPCqkgE9EfgKW48btDQMqqmpXVf0uIOs3QKlMreGZdbkQGA7M\nU9WfvOTSwN4g2fd471laJ+P4fNC5M6xeDVOmQM2aEB/vundr1YL334eTJ0O/XruL27H6/tW82u5V\nShUqxfzN86n3dj16fdGLnYd2Ztl9GGOMMWnVo0cP8uXLR4ECBWjYsCHr12fJ1IZ0C7WFbytwHVBd\nVcek0M0K8CNwUabULAhvZu7nwHGgu98h6+sLIz4f3HorrF0LH37ounc3bYJu3eDSS2HyZEhICO1a\nBfIVoE+jPsT1iaPPP/sgCG/9+BbVX6vOy8te5njC8Sy9F2OMMeFt27ZtdOrUKaerQWRkJFu3buW3\n335jxYoV1KpVK6erdJpcs9OGiBQGZgOXAc1VdZ3fse+AfaraLuCcR4HngKKqejjgmA4ZMiT5c1RU\nFFFRUVl3A+ewkydd4PfMMy7wA9f6N2QI3Hwz5MsX+rX+u/O/DJg7gLmb5gJQ4/wavHTtS7Sv3t7W\n7zPGnMEmbeRtb7zxBv/73/+YOXMmmzdvztG6DBs2DP+4IlQpfUcXLFjAggULTrt+diy8XPkshxOB\n/ar6d3orEUL5BYDPgCZAG1X9PuD4eOBaVa0UkB6NCw6rBrmmzdLNZidOuNa94cNdVy+4Fr+hQ11X\nsC/E9mZVZXbcbAZ8PYCNuzcC0LZaW15q+xKXlr00S+pujMmdzuWALzo6mh49eiR/vvvuu5k4cWKW\nlRcfH89FF53q5IuMjMyWIOy3334jKioqxwO+xx57jPLly1O2bFmWLVvGo48+SqVKlVI9L9xm6cbj\nljjZ7P13vF/aFmC/iMR5W7BlKhHx4Xb4iAI6BAZ7npnAhSLSzO+84sAN3jETBgoUgO7d4X//g3fe\ngcqV4b//hVtugXr13EzfUH4viwjX1biONb3W8NK1L1GiYAnmbppL3Tfr8tBXD7HnyJ7UL2KMMeeI\nDh06MHToUDp27Jhins6dO+Pz+TLUDVmqVCmGDh3K0KFDKVGiREi9Lr1798bn8+Hz+dizJ32/u9Ma\n0B87doz+/fvTpEkTypcvT6FChahQoQLNmzfnww8/JDExMdVrHD58mNatWxMREUFkZCQAl112Gb17\n9+aOO+6gY8eO3H333em6n6ySlr10n8RNjJgB/AWUAzoDJYFxQDOgPdBTVTPtnxAi8iZwHzAS+DLg\n8O+qutVbeHkxUInTF16uA1xuCy+Hp+PH3QzekSPdci4AV1wBw4bB9de72b+h2HloJ0MWDOHtH98m\nURMpVagUw6KGcX+D+ymQr0DW3YAxJuxZC18PoqOj6dq1a4r5du3aRYUKFTjpzapbvHgx11xzTYbK\nrlKlCj6fj19//TXFPDExMbRu3ZoiRYpw+PBhdu7cSenSpdNcVnx8PC1atAi5hW/Xrl1ERkZy1VVX\nUb16dcqUKcOOHTuYNWsWf/31F7fccgtTpkwJ6VpRUVFUrFiRyZMnc/LkSfLnz59cp4suuoi9e/dS\nokSJs14ju1r4UNVUX7hlWD4Jki64APBl7/MkYGUo1wz1hWtFTMB1HQe+BvvlKwWMB3YDh4B5wGVn\nua6a8HDkiOprr6mWL6/q2vhUGzRQ/fJL1cTE0K+zevtqbfV+K2UoylC01uu19Ku4r7Ku4saYsHcu\n/66fOHGiioi+//77Z8330ksvqYhov379VES0Z8+eGS47MjJSq1atmuLxgwcPatWqVbVz584aFRWl\nIqK7d+9OV1mbN2/WKlWqhJw/MTFRT5w4EbROtWvXVhHRVatWpXqd48ePa0REhL799tu6bNkyLVGi\nhJ48eVJVVVevXq358uXTgwcPpnqdUL+jXr50x1Ohdul2we2qERgsqpd+h5c0DagZ4jVDoqpVVTWf\nqvqCvIb75durqj1V9XxVjVDVNqpqi7blAoUKwYMPugkdY8dCuXJuy7brroOrr4a5c0Pr6r2s3GXM\n6zKPz279jGqlqrF+13r+9eG/uO6j69iwa0PW34gxxuRCEydOpGjRoowcOZJLLrmEqVOncvjw4TPy\nLViwAJ/Px7Bhw1i8eDGtW7emRIkS6WqVGzRoEPv372fcuHGoarZOuhOR5JY4fxEREVx77bUAbN++\nPdXr/Pjjjxw+fJhmzZpx4YUX8sgjj5DPm4W4ZMkS/v3vfxMREZG5lc+AUAO+okCZFI6VAYp5//03\nrjXOmDQrXBj69oVff4UXXoCyZWH5cmjXDpo0gW++ST3wExFuqnkT6x5Yx5g2YyhesDiz42Zz2ZuX\n0W9OP/YeCbZcozHGnJtWrFjB2rVr6dSpE0WKFKFr164cPHiQqVOnpnjOkiVLaNWqFQULFqRXr15n\nHRsYzMKFC3nzzTd58cUXKVeuXEZvIdMcPXqUb7/9loiICK688spU88fGxlK2bFlq1qxJpUqVaNSo\nEWPGjOGFF15g48aNvPfeGe1kOerMEDe4hcBIEVmvqknbmiEiDXFj62K8pOq4SRzGpFuRIvDww3D/\n/TBuHDz/PCxdCm3aQNOmbpZvaivoFMxfkEeueYQudbvwdMzTvPfTe7yy/BUmr57M8BbDuffKe8nv\nC/Xrb4wxedOECRMA6NKlCwB33XUXTz31FBMmTKBbt25Bz/nmm2/44IMPuOuuu9Jc3uHDh+nZsyet\nW7dO8fqBVq5cyRtvvEGJEiUoXLgw1atXTx6TOHHiRL755hu2b9/Oww8/TPv27WnVqlVI1z1y5Aij\nR49GVdmxYwezZ8/m2LFjTJkyhTJlzmzjmjZtGjNmzKBSpUqULVuWxYsX07Rp0+Tjbdq0oU2bNiGV\nnRNCnbRxEW5MXFXgN2AHbtJGZeBX3JIov4rIAOCoqr6RdVXOHDZpI/f4+2947TXX6rfXa6Br0cJN\n7vD7f+2sVm1fRb85/Vj4m9tuuc4/6vBy25dpfVHrLKq1MSYcpHXShgzL/vU8dUjW/C1KbdLG0aNH\nKV++PEWLFuX3339PTm/VqhUxMTFs3LiRiy++ODl9wYIFtGzZkgYNGvD998EWzDglpUkbDz30ENHR\n0axZsyZ5dmtUVBSxsbFBJ2388ccf3HLLLXzzzTcUKVKEr776iptvvpkDBw7gC3UtrxTs2rWLf/zj\nH8nfEZ/PxwMPPMDTTz9N2bJlT8s7evRoxo8fz6pVqyhSpAjTp0/nzjvvZMyYMTz00EMZqkdYLcui\nqr8CtYD7ca15e4BvcbNna3nHUdWXckOwZ3KXYsXgiSfc2n3PPAMlS0JMDDRr5lr9li1L/Rr1LqhH\nzN0xTL95OlVKVmHtjrW0mdSGG/9zI3G747L8HowxJtx88skn7N+/nzvuuOO09KTWvqTWv0ANGjRI\nV3mxsbGMGzeOESNGJAd7odSxWLFiFClSBHCB5NixYzMc7AGUKVOGxMRETp48SXx8PKNGjeLdd9+l\nUaNG7N+/Pznfd999xxNPPMG7776bXI86depw4sSJ01r4wl5qszqA84CxQMOMzA4Jtxfn8Myt3G7v\nXtUhQ1SLFz81q7ddO9Xly0M7/8iJIzoqdpQWfbaoMhQtMLyAPjz3Yd13ZF+W1tsYk/3O5d/1qc3S\nbdmypYqIrl69+rT0AwcOaJEiRfTCCy/UhISE5PSYmBgVER06dGiqZQfO0j1x4oRefPHFes0112hi\nwPILzZs3V5/Pp7t27TrjOh988IGKiF599dX69NNP66ZNm1ItOyNefPFFFREdMWJEclrr1q21evXq\np+UbP368lihR4ox7SY9Qv6Nk9SxdVT0O3AsUzrKo05g0KFnS7c6xeTM89RQULQpz5kCjRm79vh9/\nPPv5hfIX4rEmj7HxwY30qNeDk4kneXHZi1R/rTrv/PgOCYk278gYk7fFx8cTE+OG319++eXJix/7\nfD5KlCjBkSNH2LZtG3PmzDnj3PTMqD148CCbNm1i2bJl5MuX77TyFi1ahKpStmxZfD4fW7acmgpw\n55138vDDDxMXF8eIESOoW7cuK1asSP+Np6J1azfMZ+3atQDs2bOHmJgY2rZte1q+RYsW0bhx41y1\npWeoo9ZX4fawXZSFdTEmTUqXdl28ffvCiy+6cX5ffuleN97ogsIrrkj5/PLFyjP+pvE80PAB+s3t\nx+Iti7nvi/sYt2IcY9uOpUXVFtl2L8YYk52StliLioo6bZxekj179jBjxgwmTJhA+/btM1xeoUKF\n6NmzZ9AA6YsvvmD79u106dKFggULUrRo0eRjPp+PMWPGMGbMGJYvX06XLl2Ijo6mYcOGGa5TMNu2\nbQOgWDG3+MimTZtITEzkn//852n5Fi1axH333Qe44LlKlSpZUp9MFUozIHA1biu1G/AmeuT2F+dw\nM39etWOH6iOPqBYufKqrt1Mn1YDeiqASExP147Ufa+TLkckLN3ec0lE37cna7gNjTNY6l3/Xp9Sl\nm5CQoJGRkZo/f37dunVr0HMTEhK0fPnyWrBgweSu1qQu3WHDhqVadmoLL/tL6tINXHj5uuuu0/bt\n25+WNnDgQB0zZkxI103J+vXr9dChQ2ek79mzR+vXr68iorNmzVJV1Q0bNqiIaExMTHK+jRs3qojo\n0qVLddeuXTpy5MgM1SfU7yjZtPDyVKA08DlwRER+915bkt4zOQ41Js3KloUxY9w6fv37uwWdZ8yA\nunXdfr3r1qV8rohwS+1bWN97PSNajCCiQASfbviUWuNq8dg3j3Hg2IHsuxFjjMlC8+fPZ8uWLbRp\n04YKFSoEzePz+ejatSvHjx9n0qRJWV4nF8+cbt26dUT5rcG1efNmvvvuO+65554MlTVlyhTKly/P\n9ddfT+/evRk0aBC33347kZGRrFy5kh49enD99dcDUKNGDerUqZO8bdu+ffsYOHAgIkK1atWYNm0a\nnTp1ylB9skuoXbrzUzlu65uYsHHBBfDSSzBwIDz3HLz9NkybBtOnw223weDBUDOF/WAKFyjMk82e\npPsV3Xl8/uN88PMHjF4ymuhV0YxsOZJu9bqRz5cve2/IGGMy0YQJExCRVNfB6969O88//zzR0dH0\n69cvy+ojIkG7eqOjo4mNjeXxxx/n4MGDqCqTJ0+mZMmSGSrvhhtuYNu2bSxdupQlS5Zw6NAhypQp\nQ/PmzenZsycdOnQ4rW5Tp05l4MCBrF27FhFh/PjxjBo1in79+lGjRg1qpvQHJcyEtA5fXmTr8J07\ntm6FZ5+F996D48fB54M77nCBX/XqZz93xdYV9Jvbj6W/LwXgiguuYGy7sTSLbJYNNTfGZFRa1+HL\nS5LW4Zs4cSJ33313tpad0jp85kxhtQ6fMbnZhRe6HTvi4uC++yBfPpg8GWrVgu7dXRdwShpe2JDF\n3RfzUaePqFi8Iiu3r6R5dHNunnYz8fvis+0ejDEmvbp3747P56N79+5ZWk58fHzyzFv/mbYmPITc\nwici9YGngWZASdy6fD+JyChgoaqeOXc7jFkL37krPh5GjoSJEyEhwQWA3bq5JV7ONtHq8InDvLD0\nBZ5b/BxHTh6hYL6CPHz1wzze9HGKnlc05RONMTnmXG7h+/nnn/nss8+Su0vr1avHjTfemGXl7d+/\nn7FjxyaXV7JkyQzvQnEuyK4WvlC3VmsCfIPbRm0+0Bto4AV8I4HaqtrhbNcINxbwmV9/hREj4IMP\nXOCXPz/06AFPPgmVK6d83h8H/uCxbx7jwzUfAnBB0QsY1WoUXS/vik+s0dyYcHIuB3wmdwi3gG8x\nsBvoiOsGPs6pgK8zMFZVK6W3EjnBAj6TJC7Oref34YeQmAgFCsD//R88/jhUrJjyed/98R195/Tl\n+61uT8kGFRowtu1YGldunE01N8akxgI+E+7CbQxffeAtVU0McmwXUDZIujG5QvXqrpVv3Tq4/XY4\neRLeeAMuvhgeegj+/DP4eVdVvIplPZcxqeMkKhSrwA/bfqDJxCbc/sntbNlv41eMMcaEj1ADvqOk\nvLXaBcD+FI4Zk2vUrAkffQRr1rh1+44dc7t3XHQRDBgAf/115jk+8XFX3bvY+OBGnm72NIXyF2LK\n2ilc8volDI4ZzKHjh7L/RowxxpgAoXbpzsRN1GjpJR0HGgArgbnAblW9PasqmRWsS9ekZs0aGDYM\nPvnEfS5cGB580K3vVzaFNu0t+7cw6JtBTFk7BYALi13Ic62f447L7rDxfcbkAOvSNeEu3MbwXQ4s\nxW2vNg0YDLwK1AOuxM3Y3ZDeSuQEC/hMqFatcvvyfv65+xwRAX36wCOPwPnnBz9nyZYl9J3Tlx//\n/BGARhc2Ymy7sVxV8arsqbQxBrCAz4S/sAr4vILqA2Nwy7LkAxKBWGCAqq5MbwVyigV8Jq1+/NEF\nfl984T4XLQp9+7ru3tKlz8yfqIlM+nkSj89/nD8PuoGAd152J8+1fo6Kxc8yG8QYk2ks4DPhLuwC\nPr8CC+P21d2nqrl2gJIFfCa9vv8ehgyBOd7Kk8WLQ79+bv/eYDv+/H3sb55b/BwvLnuRYwnHKJy/\nMIMaD2Jg44EUKVAkeytvzDnGAj4T7sI24MsrLOAzGbVsmQv85s1zn0uUcK19ffu6/w60ee9mHv3m\nUab/dzoAFYtX5PnWz3NbnduC7iNpjMk4C/hMuAu7gE9EqgG3AJWAQoHHVbVHeiuREyzgM5klNtYF\nfjEx7nOpUm58X58+UKzYmfkXxi+k39x+rNq+CoCrK17NK+1eoeGFDbOx1sacG+wfUyY3CJuAT0Q6\n4CZrCLADOOZ/GFBVrZreSuQEC/hMZluwAAYPdgEguAkdAwdC795uvJ+/hMQEJq6ayJPfPsmOQzsA\n6NcXvtYAACAASURBVHp5V0a1GkWFYhWyt+LGGGPCXnYFfGuAP4E7VXVnegsLJxbwmaygCt9+6wK/\npUtdWtmyMGgQ9OoFRQKG7B04doCRi0YydvlYjiccJ6JABI83eZwBVw+gcIGUlr40xhhzrsmugO8Q\n0ElV56a3oHBjAZ/JSqpubN/gwbB8uUsrVw4eewzuu8+t6edv055NDJw3kE83fApAZIlInm/zPDdf\nerN1SRljjMm2gO8n4AVV/Si9BYUbC/hMdlCFr75yY/x++MGllS8PTzwB99wDhQJGw367+Vv6z+3P\n6r9WA9C0clPGthtL/fL1s7nmxhhjwkl2BXytgbHATaq6Kb2FhRML+Ex2UnXr9w0ZAiu9VSsvvBCe\nfBJ69ICCBU/lTUhMYPzK8Tz17VPsPLwTQeherzsjW43kgqIX5MwNGGOMyVHZFfDFAtWA84E4YI//\nYdykjWbprUROsIDP5ARVt2PHkCGw2jXiUbkyPPUUdOsGBQqcyrv/6H6eWfQMry5/lROJJyh6XlGe\nbPok/a7qR6H8Z0yUN8YYk4dlV8C3AFBccBeMqmqL9FYiJ1jAZ3JSYiLMmOF27li3zqVVqQJPPw1d\nupwe+MXtjuPhrx9m1sZZAFQtWZUXrn2BjjU72vg+Y4w5R9jCy+lkAZ8JB4mJMG2aC/w2eLtRX3SR\nm+xx552QP/+pvPM2zaP/3P6s2+kixKgqUbzc9mXqXVAv+ytujDEmW1nAl04W8JlwkpAAH38Mw4bB\nxo0urXp11/V7222QL59LO5l4krd/eJshC4aw+8huBOGe+vcwouUI/hHxj5y7AWOMMVkqowGf7ywX\nfklEKgWkdRKRkgFpNURkZnorYIxxAd0dd7ju3fffh2rVIC4O7roL6tSBKVNca2B+X356/7M3cX3i\n6NuoL/l8+Xj3p3ep/lp1Xlj6AscTjuf0rRhjjAlDKbbwiUgicJWqfu99zg8cBxqo6k9++a4Clqpq\nisFjOLIWPhPOTpyAyZNh+HCIj3dptWu7Fr/OncHn/d+2YdcGHv76YWbHzQbg4tIX80KbF7jxkhtt\nfJ8xxuQhWdbCZ4zJOQUKQPfu8L//wdtvu5m869bBLbdAvXrw6aduxm/NMjX58o4v+erOr6hZpia/\n7PmFDh93oM2kNqz5a01O34YxxpgwYQGfMWHsvPPg3nvduL433nBr961ZA506Qf36MHOmC/zaXdyO\n1fev5pV2r1CqUCnmb55Pvbfr8cCXD7Dz0P+zd+dxOpf7H8df31kJY5mESKSRiKMoIRpZ09GiTZvU\nKRIythjbGLKHKSSlcjqH9vzaHHtDSiklyhE5lmQJ2WKMWT6/P645nWmyjNvc9z0z3s/Hw2N8r/t7\nf6/P/FNv1/W9rqtQnIYoIiJnQYFPpACIjHRn8f74I0ye7E7rWL0abrkFrrkG5s6FsJBwnmjwBBt7\nbKT71d3x8Jj21TRiJscwacUkvd8nInIOO5PAd7IX3vQinEiAFCkC3bvDpk0waZI7n/err+Cmm6Bh\nQ1iwAMoUjWZy28l8+9i3tKrWioOpB+m9oDe1p9Xmww0fondXRUTOPadbtPEtcPi/TUBjYDXwW7Zb\no4A6WrQhEnhHj8Lzz8OYMbAna+a2USO32OOGGwCMuRvn0ntBbzbsc/u9tKrWikmtJ1GzbM2g1S0i\nImfGb/vwZZ2ukVs6aUMkiI4cgalTYdw42LfPtTVt6vb1i42F4xnHee7L5xiWPIyDqQcJ9ULpWr8r\nic0SKVO0TFBrFxGR09PGyz5S4JPC6PBhmDIFxo+H/ftdW7NmLvg1aQJ7j+5l6MdDmb5qOpmWSeki\npUmMTeSx+o8RHhp+6oeLiEjQKPD5SIFPCrNDh+DZZ2HCBDhwwLW1aOGCX6NGsHb3WnrN78XizYsB\nuPz8y5nYeiJtLm0TxKpFRORkFPiyyToZZBLQAvfO4SIgzsx+OsG9CnxS6B044BZ3JCW5EAjQurUL\nftdcY7z/w/v0WdCHTfs3AdA2pi0TWk2gxvk1gli1iIjkpMCXxfO883CLTFKAwVnNTwHn4RaVHM1x\nvwKfnDN+/fV/we+3rCVXbdu64Fe7biqTV05mxLIRHEo95I5vu7obCdcnULpo6eAWLiIigALf7zzP\n6wlMAKqb2X+y2qoAG4EnzWxSjvsV+OScs3evm+adPNkt9AC4+WYX/CpcupshHw9hxtczMIzootEM\nbzaczvU6ExYSFtzCRUTOcQp8WTzPWwxEmFmTHO3JAGYWm6NdgU/OWXv2uIUdU6ZASopra98ehg2D\njLKriZsXx9KtSwGoVbYWk1pPomW1lsErWETkHKfAl8XzvF3AHDPrmqP9OeAOM7sgR7sCn5zzdu+G\nsWNh2jQ4dsy13XUXDB1q/BAyh74L+rL5wGYA2lVvx4RWE4iJjglixSIi56azDXxnvFmy53kXeJ5X\nOecfXwvIQ6WB/Sdo/zXrMxHJoVw5mDgR/vMfeOIJd4Tbm29C7doeb49oz5wW6xjdfDTFI4rzwYYP\nqPVcLfou6MvBYweDXbqIiJyBXAU+z/NKep430/O8FGAXsCXHn83+KU9EAqFCBXjmGXdkW7duEB4O\nr70GV9UpwvfTBzD/po08XPdh0jPTmbBiAjGTY5j+1XQyMjOCXbqIiORCrqZ0Pc/7B3A78BLwHZCa\n8x4zm5nXxZ2J00zp3m5m5XK0a0pX5CS2bYNRo+DllyEtDUJD4f774bbHV/H093Es37YcgDrl6pDU\nOolmVQvUQTsiIgVOQN7h8zxvD5BoZlN87cjfTrNo409Hv3meZwkJCb9fx8bGEhsbG4BKRQqOLVvg\nqadg5kzIyHDBr+ODxlX3v8XT3z7J1oNbAbitxm083eppLil9SVDrFREpLJKTk0lOTv79OjExMWCB\n7x4zW+RrR/6WtS3L07htWTZntVUBNgD9tS2LiO82bXLB79VXITMTwsKg48MplPnrRKatHc2RtCNE\nhEYQ1yCOQU0HERUZFeySRUQKlUCN8E0BMsysp68d+dtJNl4eARRDGy+L5ImNG2HECJg1ywW/8HC4\np/MOjjYayNsb/w5AuWLlGHnDSDrV7URoSGiQKxYRKRwCFfjaAc8AS4GPcCtf/8DMlvhaRF7JdrRa\nS/54tNq2E9yrwCfio/XrYfhweP11MIOICLj18S/ZVL0nq35ZAcCV5a8kqU0STS9uGuRqRUQKvkAF\nvszT3GJmVqD+Ka/AJ3L2vv/eBb8333TXkUWMZj1eZ80FT7LjyHYA7qh5B+NbjqdKqSrBK1REpIAL\nVOCLPd09ZpbsaxHBoMAnknfWrHHHs737rrsuUuIoV3UfzzfFxpKSnkJkaCR9GvYhvkk8xSOKB7dY\nEZECSCdt+EiBTyTvrV7tjmd77z13XbTcdqo+OoB1YbMAKF+8PKObj6bjXzoS4p3xvu8iIuesgAY+\nz/OigWuBMrj3+FaY2Z/e5ysIFPhE/GfVKhf8PvzQXReN+ZzS9/RkR8hKAOpVqMczbZ6hceXGwStS\nRKQACVjg8zxvJNAHiMjWnApMMLPBJ/5W/qXAJ+J/K1dCQgLMmwd4mRS5ZjbhbQZw2PsZgLtr3c24\nluOoXDI/nM4oIpJ/BeodvjhgIu6kjVm449XKA/cBfwN6mdkzvhYRDAp8IoGzYoULfgsXAuFHiGw+\nloxrx5POMYqEFaFfo370b9yfYhHFgl2qiEi+FKjAtx6YZ2ZxJ/hsEnCjmdXwtYhgUOATCbzly13w\nW7IEKLmN8Lb9SbvsdQAqlqjImBZjuLf2vXq/T0Qkh0AFvmPAX0900obneS2BD80s0tcigkGBTyR4\nli6FoUNh2TLgok8J/WtPMsqtAqBBxQYktUni2krXBrdIEZF85GwDX27/Gf0rUPskn9UE9vlagIic\ne66/HpKTYdEiaHRRYzKeXwn/9wohR8rzxc9f0PClhtz/7v1sP7Q92KWKiBQKuQ187wIjPM/r6Hle\nGIDneWGe592LO77sHX8VKCKFk+dB8+Zumnf+vBAaRHYi85kNsGwgpEcya+0sqk+uTmJyIkfTjp7+\ngSIiclK5ndKNwh2p1hjIxI34lcEFxuXATWZ22I915jlN6YrkL2bwr3+5qd5VmzZDyyeh1tsAVCpx\nEeNajqXDFR3wPJ9nNERECqxAbssSArQFmvK/ffiSgX8VxOSkwCeSP5m5/fsSEuCbX5dCmziosBqA\nays24tkbk7i64tVBrlJEJLB00oaPFPhE8jczd2LH0IQM1obOhOYDofgvANxfuyNjW47mwhIXBrdI\nEZEAUeDzkQKfSMGQmenO6B3y1CHWlx0J1yZB2HEivWIMbBpPv8a9KRpeNNhlioj4ld8Cn+d5mcC1\nZrYy6+8GnKwjM7NQX4sIBgU+kYIlMxPeegsGPb2JTdX6weVzADg/7GKebTeODrXv1Pt9IlJo+TPw\nDQNeNLOfs/5+SmZ22nvyEwU+kYIpIwPeeAP6P7+E7bV6Qfk1AFQvch3/uC+JayrVC3KFIiJ5T1O6\nPlLgEynY0tNh1uwM+r42g721B0OxvWAe15fsxOy/jeLCqPLBLlFEJM8EZONlz/Ne9jyv6kk+u9jz\nvJd9LUBExBdhYfBgx1B2ftCFKZdtJOq73pAZxtJDr3DRuBjumzaGo8ePBbtMEZF8Ibf78P3+Pt8J\nPqsPrDSzAnX4pUb4RAqXtDQYO2MDo77qS0rlDwCIOFKVuJpPM7rjbYSE6P0+ESm4AnW02qmUA1Ly\n4DkiIj4LD4fBXatzYNr7PFFqAWG/1uJ4sc2M23o7JZ+4gWff+Bb9G09EzlWnWrRxG3AbbmXufcA8\nYG+O284DmgA/mFlTP9aZ5zTCJ1K4HUlJ529TX+TNvUOwovvAPM7f9gjP3PIU99x8AVrQKyIFiT9X\n6cYBcVmXlYFdQGqO21KBdUC8mf3gaxHBoMAncm7YsX8/904fztKjUyA0HY5FUXnLEKZ2fIKb2kQo\n+IlIgRCQVbqe520BbjWz1b52lN8o8ImcW775aT33/r0P6zPmuoZ9l1Jj2wSe7daOFi08BT8Ryde0\nLYuPFPhEzk3/9908Or/biz223jVsasFVeyYyoX9tYmODWpqIyEkFNPB5nlcaqA5E5vzMzJb5WkQw\nKPCJnLvSMtKYtHwaw5KHkcJ+yAyBVV1odDyRMUPK0qRJsCsUEfmjQE3pFgFeAe7ixMer6Wg1ESlw\n9h3dx8CFw5ixehqZZMCxkpCcwA0luvFUYgQNGwa7QhERJ1DbsgwBYoEHs667AX8DPgE2Ae18LUBE\nJFiiz4tm+i2TWfv4Gppf3BqKHIQ2vVlyWW0adfqQNjcaK/+0+6iISMGT2xG+9cAzwIvAcaC+mX2d\n9dnbwA4ze8KfheY1jfCJSHZmxtyNc+n5r95sOrDBNf7YCuZP4qZrapKYCPV0TK+IBEmgRvgqA98B\nGUAaUCzbZy8Dd/tagIhIfuB5HjdVv4l13dcysdVEoiJKwqULoGsdPsrsQf0m+7jlFlhdaPYqEJFz\nSW4D3z6gVNaQ2HagbrbPooGieV2YiEgwRIRG0KthLzb1/JGu9bsSEmrQYAo8EcP7u5/lynpptG8P\na9YEu1IRkdzL7ZTu28C3ZjbC87yxQBdgHJAO9AFWmVlbv1aaxzSlKyK5sXb3WnrN78XizYsB8PZe\njs2bCD+24c47ISEBatUKcpEiUugFapXu1UBlM3vH87wo3Irdm4FQ4HPgHjPb6msRwaDAJyK5ZWZ8\nsOED+izow4+//ghAyI9tyfzXRLxfL+Puu13wq1EjyIWKSKEVtI2Xs7ZqiTSzg752HkwKfCJyplLT\nU5mycgrDlw3nUOohQiwMVnYn8+OhhBwvzb33wtChEBMT7EpFpLDxe+DzPC8SN4rX38wW+NpRfqPA\nJyK++uXILwxZMoQXv34RwyiSGc3x+cPJ/LIzIYTRsSMMGQKXXBLsSkWksAjUlO5+4HYzW+JrR/mN\nAp+InK1vd31L3Pw4krckA1A6rRYH35hE5o8tCQ2FTp1g8GCoUiWYVYpIYRCowPcWsMnMBvjaUX6j\nwCciecHMmLN+Dn0X9GXzgc0AXHS0HdtfnoDtjSEsDB5+GAYNgsqVg1ysiBRYgQp8TYBZwFvAHGAn\n8Icvmtl/fC0iGBT4RCQvHUs/RtLnSYz8ZCS/Hf+N8JBwLt3Xg/UvDMFSShEeDo8+CvHxUKlSsKsV\nkYImUIEv8zS36CxdERFg5+GdDFoyiJmrZ2IYpSPP59KtT/HlC49AZigREdClCwwYABdeGOxqRaSg\nCFTg63S6e8xspq9FBIMCn4j406odq4ibH8fybcsBqF6yDhd+m0TyzGYAFCkCjz0G/ftD+fLBrFRE\nCoKgbctS0CnwiYi/mRlvrXuLJxc+ydaDbqvS5hfeRtiSp5n/ulvCW7QodOsGTz4JZcsGs1oRyc8U\n+HykwCcigZKSlsLEFRMZvXw0R9KOEBEawb1V4/jlnUHMnRMFQLFi0L079O0L558f5IJFJN8J1JTu\nK+RYpJH9Y9w7fA/7WkQwKPCJSKDtOLyD+MXxvPrtqwCUK1aOR6uN5JuXO/HRh+416OLFoWdP6N0b\nypQJZrUikp8EKvBt4Y+BzwPKAMWBg8ABM6vqaxHBoMAnIsHy5c9fEjc/js9++gyAK8tfSZeLn+G9\nZ5vwr3+5e6KiIC4OevWCUqWCWKyI5AtBndL1PK8p8Dxwv5l97fODgkCBT0SCycx4/bvXeXLRk2w/\ntB2AO2veSYfocTw/pgoLF7r7SpaEPn3cqF9UVBALFpGgCvo7fJ7nPQo8aGbXndWDTvzs6kAPoDlw\nEXAY+BIYYmZrTlJLH6AKsAWYZGbTT/JsBT4RCbqjaUd5+rOnGbN8DCnpKUSGRtKnYR+ahsQzdkRx\nPv7Y3Ve6tHu/r0cPKFEiuDWLSODlh8DXAnjfzM47qwed+NndgceAmcBXQCngSaAucF32UcWssPc8\nMApYBLQABgLdzOz5EzxbgU9E8o3th7YzYNEAZq2dBUD54uUZ3Xw0lfd3JHFYCMuWufuio6FfP7ey\nt3jxIBYsIgEV7CndcOBloIGZVff5QSd/frSZ7cvRFoUbvfvAzB7MagsDdgAfmdlD2e59CbgZqGBm\n6Tmeo8AnIvnO59s/p+e8nqz8eSUA9SrUI6n1M6T+2JiEBPj0U3df2bJuD7+uXeG8PP/ntojkN4Fa\ntPExf16lGwlUB6KBx8zsBV+LOFOe530BHDKzllnXTYClQEszW5ztvlhgCXCDmSXneIYCn4jkS5mW\nyey1sxmwaAA/H/4ZgA5XdGBM87H8sLIyQ4fCF1+4e8uVc6d2dOni9vQTkcLpbANfSG77yfoTku3P\nYeAdXJgKZNgrA1wB/Dtbc62sn9/luH1d1s/L/V2XiEheCfFCuL/O/fzQ/QeGNB1CkbAivP7d69SY\nehnLw4eyaOkR5s6F+vVh9263krdaNZgyBY4dC3b1IpIfFbiNlz3PmwXcAtQxs/9ktQ0EngKKmNnx\nbPeGAcdxizxG5niORvhEpEDYemAr/Rf1543v3wCgYomKjGkxhnuuuJe5H4WQkADffOPurVQJBg6E\nhx+GyMggFi0ieSpQI3x5wvO8Fp7nZebiz5KTfD8euAfo/t+wJyJS2F1c6mJev+N1lj+0nHoV6vHz\n4Z95YM4DNH65EWWv/JxVq+Ddd6FOHdi+HR5/HKpXhxdfhLS0YFcvIvlBbt/he5CTn7SRidt8+Rsz\n236a5xTFba9yOkdzPsvzvMeA54BBZjY6x2ddgam4xRm7s7VfAOzCrdSdluM7lpCQ8Pt1bGwssbGx\nuShNRCR4Mi2TV799lfjF8ez6bRcA99W+jzEtxnBh8Uq8+y4MGwbff+/ur1IFhgyBBx6A8PCglS0i\nZyg5OZnk5OTfrxMTEwOyaCMzF88y4A2gU/Zp1bzged4DuK1ZJpjZkyf4vCmQzMkXbTQzs6U5vqMp\nXREpsA6nHmbM8jFMWDGB1IxUioYVpX/j/vRr3I8ioefx1lsu+K1f7+6vVs0Fv/vug7CwoJYuIj4I\n1CrdRsAs4H3cQo3dQDngTuCvQDegJjAceMbM4n0t6AR93wa8CbxkZo+d5J7/bsvyYfYzfT3Pm4F7\n30/bsohIobR5/2b6L+rPW+veAuCiqIsY22IsHa7oQGamx+uvw/DhsGGDuz8mBoYOhXvugdDQIBYu\nImckUIHvXeCHEwU5z/NGA5eb2a2e540A7jOzS3wtKMezmwILcKtve/DHaeVUM/sm271dcFO+o4DF\nwA3AINz7fn+Yzs26X4FPRAqNZVuXETcvjm92uf8sNrqoEUmtk7i64tWkp8Ps2S74bdrk7q9RAxIS\n4K67ICSgb3OLiC8CFfgOA7dmny7N9llL4F0zK+F5XivcKFuErwXleHYCkIALejl/yS05g6XneZ1x\nR6tdDGzFHa32p1M2su5V4BORQiUjM4OZq2cycMlAfjnyCwAd/9KR0c1Hc2GJC0lPh1dfhREjYMsW\n951atdzUb/v2Cn4i+VmgAt8+YJyZjT3BZ/2BJ80s2vO8NsBrZlba14ICRYFPRAqrQ6mHGLlsJElf\nJHE84zjFwosRf108vRv2pmh4UdLSYOZMeOop2LbNfadOHRf8br0VPJ//lyIi/hKowDcNeAgYDLwF\n/AJcANwFjABeNrPHPc8bANxiZg19LShQFPhEpLDb9Osm+i3sx5z1cwC4uOTFjGs5jjtr3onneaSm\nwssvw8iR8LM70IMrr3TBr107BT+R/CRQge884AXcHnjZOzNgNtDZzFI8z/sr7sizZb4WFCgKfCJy\nrliyeQm95vdize41AFxX+TqSWidR78J6gDud48UXYfRo2LnTfad+fRf82rZV8BPJDwIS+LJ1dhnQ\nAKgA7ARWmtl6XzsPJgU+ETmXZGRmMOPrGQz+eDB7j+7Fw6NT3U6Maj6K8sXLA5CSAtOnw5gx7sg2\ngGuucYs9WrVS8BMJpoAGvsJEgU9EzkUHjh1gxNIRPLvyWdIz0ykeUZxBTQYRd20cRcKKAHD0KEyb\nBmPHwp497nuNGrngd8MNCn4iweC3wOd5XmVgl5kdz/r7KZnZNl+LCAYFPhE5l23Yt4G+C/rywYYP\nAKhaqirjW46n/eXt8bIS3ZEjMHUqjBsH+/a57zVtComJoIOJRALLn4EvE7jWzFbm4qQNM7MCtYWn\nAp+ICCzctJBe83vx/R53FltslVgmtZ5E3fJ1f7/n8GGYMgXGj4f9+11bs2Yu+DVpEoyqRc49/gx8\nnXB76u3N+vspmdlMX4sIBgU+EREnPTOdF1a9wNCPh7IvZR8eHo9c9QhP3fAUFxS74Pf7Dh2CZ56B\niRPhwAHX1qKFm+ptmO/3ZhAp2PwZ+HoCr5vZbs/zLgZ25vUZucGkwCci8kf7U/YzfOlwpnw5hfTM\ndKIioxjcZDBPNHiCyLDI3+87cACSkmDSJBcCAVq3diN+DRoEqXiRQi6QU7rXmtlKXzvKbxT4RERO\nbP3e9fRZ0Ie5G+cCUK10NSa0msDNl938+/t9AL/+6kJfUhL89ptru+kmF/zq1QtG5SKFlz8D3z7g\nYTN7T4FPROTcM+/HefSa34v1e93uW82rNmdS60nULlf7D/ft3QsTJsDkyW6hB8DNN7t9/K68MsBF\nixRS/gx87wNNgdVZP78GDp3oVtyijRt8LSIYFPhERE4vLSONaV9NY1jyMPYf20+IF0KXel1IjE2k\nbLGyf7j3l1/cwo6pU92efgC33eaCX506ga9dpDDxZ+ArDwwFLgeuB9YAv53kOWZmBWqtlgKfiEju\n7Tu6j2HJw5j21TQyLIOSkSVJuD6Bbtd0IyI04g/37trltnKZNs2d4gFw552QkAC1agWheJFCIFBH\nq2UCDc3sC187ym8U+EREzty6PevoNb8XCzYtAKB6dHUmtppI25i2f3i/D9wxbWPGuNM7UlPdhs13\n3+2CX40awahepOAKVOCrAuzQKl0RETEz5m6cS+8FvdmwbwMArau1ZmLridQsW/NP92/f7s7pffFF\nSEuDkBC4914YOhRiYgJdvUjBFPCj1TzPuwAokrNdJ22IiJxbjmccZ+rKqSQuTeRg6kFCvVC61u/K\nsNhhRJ8X/af7t22DkSPh5ZchPR1CQ+H++2HIEKhWLQi/gEgBEqgRvpLAM8DdQOQJbtFJGyIi56g9\nR/aQkJzA9FXTybRMShcpTWJsIo/Vf4zw0PA/3b9lCzz1FMycCRkZLvh16gSDB0OVKgEuXqSACFTg\n+wdwO/AS8B2QmvMenbQhInJuW7t7Lb3m92Lx5sUAXH7+5UxsPZE2l7Y54f2bNrng9+qrkJkJYWHw\n8MMwaBBUPu0J7iLnlkAFvj1AoplN8bWj/EaBT0Qk75kZH2z4gD4L+vDjrz8C0DamLRNbTeSy8y87\n4Xc2boQRI2DWLBf8wsPh0Udh4ECoWDGQ1YvkX4EMfPeY2SJfO8pvFPhERPwnNT2VySsnM2LZCA6l\nHiIsJIzuV3dn6PVDKV209Am/88MP7lze114DM4iMhM6dIT4eKlQI8C8gks8EKvBNATLMrKevHeU3\nCnwiIv73y5FfGLJkCC9+/SKGEV00muHNhtO5XmfCQsJO+J1169zxbG++6a6LFIGuXaF/fyhXLoDF\ni+QjgQp87XCLNpYCHwG/5rzHzJb4WkQwKPCJiATOt7u+JW5+HMlbkgG44oIrmNR6Ei0uaXHS76xd\n64LfO++466JFoVs3ePJJKFv2pF8TKZQCufHyqWiVroiInJKZMWf9HPou6MvmA5sBaFe9HRNaTSAm\n+uQb8q1e7Y5ne+89d12sGHTvDn37wvnnB6BwkXwgUIEv9nT3mFmyr0UEgwKfiEhwHEs/RtLnSYz8\nZCS/Hf+N8JBwelzTgyHXD6FUkVIn/d6qVe6Ujo8+ctfFi0PPntC7N5QpE6DiRYIk4BsvFxYKfCIi\nwbXrt10MWjyIV1a/gmGcf975PNXsKR656hFCQ04+abRypQt+8+a566goiIuDXr2g1MnzokiBjw14\nlAAAIABJREFUFtDA53leNHAtUAb3Ht8KM/vT+3wFgQKfiEj+8PXOr4mbF8cn2z4BoE65OiS1TqJZ\n1Wan/N6KFS74LVzorkuWhD593KhfVJS/qxYJrIAFPs/zRgJ9gIhszanABDMb7GsBwaLAJyKSf5gZ\nb697m34L+7H14FYAbqtxG+NbjqdamVOfu7Z8uQt+S7KWDpYu7d7v69EDSpTwd+UigRGod/jigIm4\nkzZmAbuA8sB9wN+AXmb2jK9FBIMCn4hI/pOSlsLEFRMZvXw0R9KOEBEaQVyDOAY1HURU5KmH7ZKT\nXfBbtsxdR0e7Fb3durmFHiIFWaAC33pgnpnFneCzScCNZlbD1yKCQYFPRCT/2nF4B/GL43n121cB\nKFesHCNvGEmnup1O+X6fGXz8MQwdCp9+6trKlnV7+HXtCuedF4jqRfJeoALfMeCvJzppw/O8lsCH\nZhbpaxHBoMAnIpL/rfx5JXHz4lixfQUAV5a/kqQ2STS9uOkpv2fm3u0bOhS++MK1lSsHAwZAly5u\nTz+RguRsA19ILu/7Fah9ks9qAvt8LUBERORkrql4DZ8+/Cmz28+mUlQlvtn1DdfPvJ4737qTzfs3\nn/R7ngetWrmFHXPnQv36sHu3W8lbrRpMngzHjgXwFxEJsjM5Wq0T8Dgw28zSPc8LA+4Cngf+bmY9\n/FloXtMIn4hIwXI07SjjPx3P2E/HkpKeQmRoJL0b9ib+unhKRJ56dYYZfPihG/Fbvdq1VawIgwbB\nww+7c3tF8rNATelG4Y5Uawxk4kb8yuBGCJcDN5nZYV+LCAYFPhGRgumngz8RvzieWWtnAVC+eHlG\nNx9Nx790JMQ79cSVGfzf/7mTO9ascW0XXQSDB0OnThARcapviwRPILdlCQHaAk353z58ycC/CmJy\nUuATESnYPt/+OT3n9WTlzysBqH9hfZJaJ9G4cuPTfjczE9591wW/7793bVWqwJAh8MADEB7uv7pF\nfKGTNnykwCciUvBlWiaz186m/6L+7Di8A4AOV3RgbIuxVC5Z+fTfz4Q334TERFi/3rVVq+amfu+9\nF8LC/Fm9SO75LfB5nlcOmAbMMLO5J7mnLW4fvi5mttfXIoJBgU9EpPA4cvwIYz8dy/jPxnMs/RhF\nworQr1E/+jfuT7GI02/Cl5EBb7zhgt+GDa4tJsbt69ehA4SefCcYkYDwZ+AbBdwO1DSzjJPcEwas\nBd42syG+FhEMCnwiIoXP1gNb6b+oP298/wYAFUtUZEyLMdxb+97Tvt8HkJ4Os2fD8OGwaZNrq1HD\nBb+77oKQ3O5tIZLH/Bn41uBW3044TQG9gAfNrK6vRQSDAp+ISOG1fNty4ubFsWrnKgAaVGxAUpsk\nrq10ba6+n5YG//gHjBgBW7a4tlq1XPC7/XYFPwk8fwa+I7jNlj8+TQE34DZeLlD7lyvwiYgUbpmW\nyavfvkr84nh2/bYLgHtr38vYFmOpFFUpV884fhz+/nd46inYts211a7tpn5vvdXt9ycSCP7cePlM\n0pCSk4iI5CshXgid6nZiQ/cNDLxuIJGhkcxeO5vqk6uTmJzI0bSjp31GRAQ8+qh7r2/aNKhUCdau\nhfbtoV49eP99t9WLSH53qsC3Faifi2fUA7bkSTUiIiJ5rERkCUY2H8m/u/2bO2reQUp6CsOWDqPG\nlBq8tvY1cjPbExkJjz0GGze6UzoqVIBvvoFbboFrrnGneSj4SX52qind0cD9wJUnW4Hred75wGrg\nH2YW77cq/UBTuiIi56ZlW5cRNy+Ob3Z9A0CjixqR1DqJqytenetnpKTACy/A6NHuyDaABg3cVG+r\nVprqlbznz3f4/hvmjgIDgLlmdizrsyK4TZjHAkWBv5hZgTpPV4FPROTclZGZwczVMxm4ZCC/HPkF\ngAf/8iCjmo/iwhIX5vo5R4+6qd6xY2HPHtfWqJELfs2bK/hJ3vHrxsue59UE3gEuA9KBPYAHlAVC\ngR+A9mb2b18LCBYFPhEROZR6iJHLRpL0RRLHM45TLLwY8dfF07thb4qGF831c44cgalTYdw42Jc1\n/NGkidveJTbWP7XLucXvJ21k7bV3G9ASuCir+SdgIfDuyfbo8wfP8zoAs4GfzeyiE3z+KNAHqIJ7\nr3CSmU0/ybMU+EREBIBNv26i38J+zFk/B4CLS17MuJbjuLPmnXhnMEx3+DBMmQLjx8P+/a6tWTM3\n4tekiT8ql3PFOXO0mud5pYD1QCaQbmaVc3z+KPA8MApYBLQABgLdzOz5EzxPgU9ERP7g480fEzc/\njjW71wDQpHITktokcVWFq87oOYcOwbPPwoQJcOCAa2vRwgW/Ro3yumo5F5xLge8F3AjjLqBF9hG+\nrFHIHcBHZvZQtvaXgJuBCmaWnuN5CnwiIvInGZkZzPh6BkM+HsKeo3vw8Hio7kOMbD6S8sXLn9Gz\nDhyApCSYNMmFQIDWrV3wa9DAD8VLoeXPffjyDc/zGgP3Ad1w7xDm1BA4H/hnjvZ/ANHAdX4tUERE\nCo3QkFC61O/Cxh4b6dOwD2EhYby8+mViJscwZvkYjqUfy/WzSpWCYcPcaR2DBkHx4jB/Plx7Lfz1\nr7Bqld9+DZE/yPeBz/O8cOAFYJyZ/eckt9XK+vldjvZ1WT8v90dtIiJSeJUsUpKnWz3N949/T7vq\n7fjt+G/EL46n5tSavPvvd3O1f99/lS7tTuvYsgXi46FYMfjoI6hf3+3lt3q1/34PESgAgQ/oD4QD\no09xT5msn/tztP+a43MREZEzEhMdw/v3vM+C+xdQq2wtNh/YzO1v3s4Nr97A6l1nltSio2HUKNi8\nGfr1g/POc6d1XHmlO6N37Vo//RJyzgto4PM8r4XneZm5+LMk6/5LcQsvupvZ8WyP0st3IiISUC2r\ntWT1Y6uZ2nYq0UWjSd6SzFXTr6LzB51/38svt8qWdVu4/Oc/0KsXFCkC774LderA3XfDunWnf4bI\nmTijRRue55UFrsWNmH1oZvs8zysKHM/N9ixZ9/5pO5UTOGpm2z3Pm4tblXsf/3t37zmgKW4aN9XM\njnme1xWYilucsTtbfxfgFnl0M7NpOWqxhISE369jY2OJ1WZJIiKSC/tT9pO4NJGpX04lPTOdqMgo\nhjQdwhMNniAiNOKMn7dzJ4wZA9OnQ2qq27C5QwcYOhRq1PDDLyD5XnJyMsnJyb9fJyYm+n+Vruc2\nIRoP9MBNrxpwtZl97XnefOBTMxvuaxGn6HczcPEpbkkys96e5zUFkoGWZrY42/djgSVAMzNbmuPZ\nWqUrIiJnZf3e9fSe35t//fgvAC4tcylPt3yamy+7+Yz27/uv7dvdcW0vvghpaRASAvfdB0OGQExM\nXlcvBUlAtmXxPG8gMAQYgdtw+Qugflbg6w48YGZ5vsDc87wGQGT2Jtwxb/WAO3AbMG/KWtjxM27U\n8eFs358B3IK2ZRERET+a9+M8es3vxfq96wFoXrU5k1pPona52j49b9s2967fSy9BejqEhsIDD7jg\nd8kleVm5FBSBCnz/AWaY2aisPe+O87/AdyPwTzOL9rWIM+F53kygec6TNjzP64Kb7h0FLAZuAAbh\n3v+bdoLnKPCJiEieSctI4/mvnichOYH9x/YT4oXQpV4XEmMTKVusrE/P3LLFre6dORMyMiAsDB58\nEAYPhipV8rJ6ye8CtQ9fRWDFST47DhTztQAfGCdYtJF1hFpX4C5gHnA3J3h3T0RExB/CQ8Pp0aAH\nG3tspPvV3fHwmPbVNGImxzBpxSSOZxw//UNyqFIFZsyAH36ATp3AzI36xcTAY4+5kUCR3MjtCN9m\n3Lm0z55ghK8X0NnMCtRedxrhExERf/r+l+/pvaA3CzYtAKB6dHUmtJrATTE3+fR+H8DGjTBiBMya\nBZmZEB4Ojz4KAwdCxYp5Wb3kN4Ga0h0L/A24FTfSlwbUB47gpk9fNLNEX4sIBgU+ERHxNzPjo40f\n0Xt+bzb+uhGAVtVaMan1JGqWrenzc3/4AYYPh9dec6N+kZHQpQsMGAAVKuRV9ZKfBCrwnQfMBxoD\nW3ErZzfjtlj5DGhtZqm+FhEMCnwiIhIoxzOOM2XlFIYvHc7B1IOEeqF0rd+VYbHDiD7P91fg161z\n5/K++aa7LlIEunaF/v2hXLk8Kl7yhYAEvqyOwoB7gDbABcBe3Ltys3KugC0IFPhERCTQ9hzZQ0Jy\nAtNXTSfTMildpDTDYofRtX5XwkPDfX7u2rUu+L3zjrsuWhS6dYMnn3SbPEvBF7DAV9go8ImISLCs\n3b2WXvN7sXiz2zr28vMvZ2LribS5tM1ZPXf1ahg2DN57z10XKwY9ekDfvu5YNym4FPh8pMAnIiLB\nZGZ8sOEDes/vzab9mwBoG9OWia0mctn5l53Vs1etcsHvww/ddfHi0LMn9O4NZXS6fIEUqHf4NvPH\nrVC8bNeZwEHga+AZM/vO12ICSYFPRETyg9T0VCavnMyIZSM4lHqIsJAwul/dnaHXD6V00dJn9eyV\nKyEhAebNc9dRURAX587vLVUqD4qXgAlU4JsJNAPKAZ8Cv2T9vTHurNpVQCOgJNDCzD71taBAUeAT\nEZH85JcjvzB4yWBmfD0Dw4guGs3wZsPpXK8zYSFhZ/XsFSvcubyLFrnrkiWhTx836hcVlQfFi98F\nKvD9DegJtDKzXdnaKwALgGeB14BFwGEza+lrQYGiwCciIvnR6l2r6TW/F8lbkgG44oIrmNR6Ei0u\naXHWz/7kEzfi9/HH7rp0afd+X48eUKLEWT9e/ChQgW8jMNDM3jrBZ3cBo8zsUs/z7gGmm1m+//eC\nAp+IiORXZsac9XPou6Avmw9sBqBd9XZMaDWBmOiYs35+crIb8fvkE3cdHQ39+rmVvcWLn/XjxQ8C\ndbRaJeBk++wdy/ocYAcQ4WsxIiIi4v7n3v7y9qzrto7RzUdTPKI4H2z4gFrP1aLvgr4cOHbgrJ4f\nGwtLl7op3saNYd8+t2nzJZfA00/D0aN583tI/pHbEb5vgEO4DZaPZWsvituQuYSZXZk1wjfGzC72\nV8F5RSN8IiJSUOz6bReDFg/ildWvYBhlzyvLiGYjeOSqRwgNCT2rZ5vBwoVuxO+LL1xbuXIuAHbp\n4vb0k+AL1JRuC+Aj4AAwl/8t2miLW6hxk5kt8jxvMhBpZp19LShQFPhERKSgWbVjFXHz41i+bTkA\ndcrVIal1Es2qNjvrZ5u51bxDh8JXX7m2ChUgPt6d11ukyFl3IWchkCdt1AQGA9cC5YGdwOfAU2b2\nb18LCBYFPhERKYjMjLfXvU2/hf3YenArALfVuI3xLcdTrUy1PHi+278vIQG++ca1VawIgwbBww+7\nc3sl8LTxso8U+EREpCBLSUth4oqJjF4+miNpR4gIjSCuQRyDmg4iKvLs106auRM7EhJgzRrXVrmy\nC36dOkGE3tgPKAU+HynwiYhIYbDj8A7iF8fz6revAnBBsQsYecNIHqr70Fm/3weQmQnvvutO7vj+\ne9dWpQoMGQIPPADhvh8BLGcgkFO65YB7gOpA9pl8DzAze9jXIoJBgU9ERAqTL3/+krj5cXz202cA\n1C1fl2faPEPTi5vmyfMzM+Gtt1zwW7/etV1yiXvn7777IOzs9oaW0wjUoo3LgBVAGFAc2ANE47Z1\nOQAcNLOqvhYRDAp8IiJS2JgZb3z/Bk8ufJKfDv0EwB0172Bci3FULZ03/5vOyIDXX4fhw2HDBtcW\nE+OC3z33QOjZDyrKCQQq8L2PG9W7FfgNuBpYAzwAJALtzGy1r0UEgwKfiIgUVkfTjvL0Z08zZvkY\nUtJTiAyNpHfD3sRfF0+JyLw5UiM9HWbNghEjYNMm11ajhnvn7847FfzyWqAC307gMeADIB24xsy+\nyvqsP9DGzM5+TXgAKfCJiEhht/3QdgYsGsCstbMAqFC8AqOaj6LjXzoS4uX27IVTS0uDf/7Tjfht\n2eLaatZ0U7+33w4hedPNOS9QJ20UB/abWSZwEDg/22dfAdf4WoCIiIj4R6WoSvyz/T9Z8bcVNKjY\ngJ2/7eSh9x6iwYwGfLrt0zzpIzwcHnoIfvgBXngBLroI1q2Du+6CunVhzhy34leCK7eBbwtQMevv\nG4C7sn12E+49PhEREcmHrq10LZ/97TP+cds/uLDEhXy14yuue+U67nnnHrYd3JYnfUREuA2aN26E\n555ze/etXQvt20O9evD++wp+wZTbKd1JuOPTHvE87y7gdWAjkAHUAEaa2RC/VprHNKUrIiLnoiPH\njzD207GM/2w8x9KPUSSsCH0b9mXAdQMoFlEsz/o5dgxmzIBRo2DnTtdWr56b+r3xRvB8npw8NwXq\nHb5I3JFph7Ku2wEdgPOAfwEvFrT0pMAnIiLnsq0HttJ/UX/e+P4NAC4scSFjmo/hvjr35dn7fQAp\nKW6qd/Ro2L3btTVoAImJ0KqVgl9u+T3weZ4XClwB7DSzX3ztKL9R4BMREYHl25YTNy+OVTtXAdCg\nYgOS2iRxbaVr87Sfo0dh2jQYOxb27HFtDRu6Eb/mzRX8TidQgS8VaGtmC3ztKL9R4BMREXEyLZNX\nv32V+MXx7PptFwD31b6PMS3GUCmqUp72deQITJ0K48bBvn2urUkTF/xiY/O0q0IlUFO6m4HeZjbH\n147yGwU+ERGRPzqcepgxy8cwYcUEUjNSKRpWlP6N+9OvcT/OCz8vb/s6DJMnw9NPw/79rq1ZMxf8\nrrsuT7sqFAIV+AYANwKtzCzV187yEwU+ERGRE9u8fzNPLnqSt9e9DcBFURcxtsVYOlzRAS+P514P\nHYJnnoGJE+FA1p4fLVq44NewYZ52VaAFKvA9BXTKupwH7AT+8EUzG+prEcGgwCciInJqS7csJW5+\nHKt3ucO0Gl3UiKTWSVxd8eo87+vAAUhKgkmTXAgEaNPGLe64Rrv9BizwZZ7uHjMrUHtpK/CJiIic\nXkZmBjNXz2TgkoH8csSt3XzwLw8yqvkoLixxYZ73t3+/G+1LSoLffnNtN93kgl+9enneXYERkMBX\nGCnwiYiI5N6h1EOMXDaSpC+SOJ5xnGLhxYi/Lp7eDXtTNLxonve3b597v2/yZLfQA+Dmm13wq1s3\nz7vL9xT4fKTAJyIicuY2/bqJfgv7MWe9W8d5ccmLGddyHHfWvDPP3+8Dt4XL+PEwZYrb0w/c6R3D\nhkHt2nneXb4VsMDneV4I0A5oCpQBhpnZVs/zYoGNZvazr0UEgwKfiIiI7z7e/DFx8+NYs3sNAE0q\nNyGpTRJXVbjKL/3t2uW2cpk2zZ3iAXDnnS741azply7zlUC9w1cad6LGNcBvQDHgajP72vO8fwK/\nmtkTvhYRDAp8IiIiZycjM4OXvnmJwUsGs+foHjw8Hqr7ECObj6R88fJ+6XPnThgzBqZPh9RUt2Fz\nhw4wdCjUqOGXLvOFsw18uV1oMR6oBFyHG93L3uEioIWvBYiIiEjBFBoSSud6ndnYYyN9GvYhLCSM\nl1e/TMzkGMYsH8Ox9GN53meFCm4blx9/hMcfh7AweO01qFULHngANm7M8y4LhdyO8O0B+pnZTM/z\nwoDjQP2sEb7mwP+ZWQk/15qnNMInIiKStzbu20ifBX34YMMHAFQtVZXxLcfT/vL2fnm/D2DbNhg1\nCl56CdLTITTUBb8hQ+CSS/zSZVAEaoSvOLD9JJ8V4Y8jfiIiInIOiomO4f173mfB/QuoVbYWmw9s\n5o637qDZ35v9vpdfXqtcGZ5/3o3s/e1vrm3mTKheHR55BLZs8Uu3BU5uA98GoPVJPmsKrM2bckRE\nRKSga1mtJasfW83UtlOJLhrN0q1LuWr6VXT+oPPve/nltSpVYMYM+OEH6NQJzNyoX0wMdOniRgLP\nZbmd0u0MTAGGA7OBH4GWwMVZ7Z3N7J9+rDPPaUpXRETE//an7CdxaSJTv5xKemY6UZFRDGk6hCca\nPEFEaITf+t2wAUaMgNmzITMTwsPh0UchPh4qVfJbt34TyG1ZxgB9+eOoYCYw1swG+VpAsCjwiYiI\nBM76vevps6APczfOBeDSMpcyodUE2lVv57f3+wDWr3fn8r7+uhv1i4yEzp1d8KtQwW/d5rmAbrzs\neV4V3MjeBcA+YIGZ/cfXzoNJgU9ERCTw5v04j97ze/Pvvf8GoMUlLZjYaiK1y/l3F+V169wpHW++\n6a6LFIGuXaF/fyhXzq9d54lA7cMXamYZvnaSHynwiYiIBEdaRhrPf/U8CckJ7D+2nxAvhC71upAY\nm0jZYmX92vfatS74vfOOuy5aFLp3h379oKx/uz4rgQp8vwCvAf8ws6987Sw/UeATEREJrn1H9zEs\neRjTvppGhmVQMrIkCdcn0O2abn59vw9g9Wp3Ssd777nrYsWgRw/o2xeio/3atU8CFfieA+4GSgPr\ngX8A/zSzn3ztONgU+ERERPKH73/5nt4LerNg0wIAqkdXZ0KrCdwUc5Nf3+8DWLXKBb8PP3TXxYtD\nz57QuzeUKePXrs9IIBdtRABtgQeyfkYAy4BXgbfN7LCvReSi74rACOBGXOjcAbxuZgNz3Pco0Aeo\nAmwBJpnZ9JM8U4FPREQknzAzPtr4Eb3n92bjr+64jFbVWjGp9SRqlvX/YbkrV0JCAsyb566joiAu\nDnr1glKl/N79aQV00Ua2TksBd+HCXyPgmJkV87WI0/RVBfgU2AQ8C+wGqgLVzCwh232PAs8Do/jf\ncW8DgW5m9vwJnqvAJyIiks8czzjO1JVTSVyayMHUg4R6oXSt35VhscOIPs//c60rVrjgt3Chuy5Z\n0o329ezp/h4sQQl8WR1XBu4DugMVzCy3mzifaT/zgFJA45MtHMk67m0H8JGZPZSt/SXg5qz60nN8\nR4FPREQkn9pzZA8JyQlMXzWdTMukdJHSDIsdRtf6XQkPDfd7/8uXu+C3ZIm7Ll3avd/XoweUCMJh\nsoHeliUKuBM3stcESAXewy3mmOtrEaforxqwEXjAzGad4r4mwFKgpZktztYeCywBbjCz5BzfUeAT\nERHJ59buXkuv+b1YvNn97/3y8y9nYuuJtLm0TUD6X7oUhg6FZcvcdXQ0PPkkdOvmFnoESkDO0vU8\nr53neW8Au4AXAAMeBcqb2T3+CHtZGmf9POZ53kLP8455nver53l/9zwv+6uUtbJ+fpfj++uyfl7u\np/pERETEj2qXq83CBxbyf3f/H9VKV+Pfe//NjbNu5KbZN/HD3h/83v/110NyMixeDI0bw759bu++\nqlVhwgQ4etTvJeSJ3E7DvgfUAZ4CqppZMzN72cwO+a80AC7M+vkybnVwG6A/cBMw3/vf0p3/hr/9\nOb7/a47PRUREpIDxPI9batzC949/z/iW44mKjGLuxrlcMe0Kes3rxf6UnP/7z+v+4YYb4JNPYP58\naNAA9uxxU7yXXAJJSZCS4tcSzlpuA18DM7vczEaZ2R+OH/Y8L9bzvJdz8xDP81p4npeZiz9LctT3\nsZn1MLNkM3sReByoB7TOZf0iIiJSwEWGRdK3UV82dN/Ao1c9SkZmBklfJBEzOYbnvnyO9Mz00z/k\nLHgetGrlFnbMnQv168Pu3W4lb7VqMHkyHDvm1xJ85usq3RigI+5dvspASm5W6XqeVxS4KBddHDWz\n7Z7ndQGmAT3MbGq255TGHe0Wb2ZjPc/rCkzFLc7Yne2+C3DT0N3MbFqOWiwh4fdFvsTGxhIbG5uL\n0kRERCQ/WL1rNb3m9yJ5SzIAV1xwBZNaT6LFJS0C0r8ZfPSRe8fvm29cW8WKMGgQPPywO7fXV8nJ\nySQnJ/9+nZiYGLB9+ErhNl9+ELg2q3k1MB14zR/Tu57nNQY+4fSBrymQzMkXbTQzs6U5nq1FGyIi\nIgWcmTFn/Rz6LujL5gObAWhXvR0TWk0gJjomQDW4EzsSEmDNGtdWuTIMHgydOkF4Hiwq9uuiDc/z\nQj3Pu8nzvDeBnbjRtnK4/fAAepnZdD++y/c5boQu51Kc/15/mfVzBbAXt01MdvfjguGnfqpPRERE\ngsjzPNpf3p513dYxuvloikcU54MNH1DruVr0XdCXg8cOBqAGuPVWN8r31ltQqxZs2wadO0P16vDy\ny5CW5vcyTl3jyUa5PM+bCNwLXIDb4+5N3OkWK7NG+34FYs1smV8L9LyOwEzcSOIc4FLc4pFvzKx5\ntvu6AM/hNl5eDNwADAK655zOzbpfI3wiIiKFzK7fdjFo8SBeWf0KhlH2vLKMaDaCR656hNCQ0IDU\nkJnpgt+wYbB+vWurVs1N/d57L4SFnfkz/bYPn+d5mUAKblXs1OzpKJCBL6u/+7PqiMGN2L2Nm849\nmuO+zrij1S4GtuKOVvvTKRtZ9yrwiYiIFFJf7/yanvN6snzbcgDqlKtDUuskmlVtFrAaMjLgjTcg\nMRE2bHBtMTFu6rdDBwg9g/zpz8D3Iu74tBLAv4E3cCN8GwId+PxBgU9ERKRwMzPeXvc2/Rb2Y+vB\nrQDcVuM2xrccT7Uy1QJWR3o6zJ4Nw4fDpk2urUYNF/zuugtCcrFnil9P2shaVXsrbqFGcyAU+Bb4\nEDddqsAnIiIi+VpKWgoTV0xk9PLRHEk7QkRoBHEN4hjUdBBRkVEBqyMtDf75Txf8tmxxbbVquanf\n9u1PHfwCdrSa53kX4hZBPMj/Tq74HLeQ4y0zy6c7z5yYAp+IiMi5ZcfhHcQvjufVb18F4IJiFzDy\nhpE8VPehgL3fB3D8OMycCSNHusUdAHXquOB3661uEUhOAT1LN1un9XHB7x7cKRaHzKyUr0UEgwKf\niIjIuenLn78kbn4cn/30GQBXlr+SpDZJNL24aUDrSE2FV15xwW/7dtdWt657569duz8Gv6AEvmyd\nR+COOetoZrf5/KAgUOATERE5d5kZr3/3Ov0X9eenQz8BcEfNOxjXYhxVS1cNaC3HjsGMGTBqFOzc\n6drq13cjfm3buuAX1MBXkCnwiYiIyNG0ozz92dOMWT6GlPQUIkMj6d2wN/HXxVMiskR10P0vAAAW\nP0lEQVRAa0lJgRdegNGj3ZFtANdc4975a9NGgc8nCnwiIiLyX9sPbWfAogHMWjsLgPLFyzO6+Wg6\n/qUjIV4ultHmoaNH4fnnYcwY2LPnv60KfD5R4BMREZGcPt/+OT3n9WTlzysBqH9hfZJaJ9G4cuOA\n13LkCEydCuPGwb59Cnw+UeATERGRE8m0TGatmcWAxQPYcXgHAB2u6MDYFmOpXLJywOs5fBiiohT4\nfKLAJyIiIqdy5PgRxn46lvGfjedY+jGKhBWhX6N+9G/cn2IRxQJaixZt+EiBT0T+v737DrOqOvc4\n/v1ZwRos2LBfe6LRWDBRY8FIjDXRa41oLJhExW64xnItUaNRVC4q0XgTS4yxExVFBYnGQlTsWBC7\niCCCFYV57x9rnethc2YYhmme/fs8z372nLXX2Xvtd87MvLPW3mubmTXHm1Pe5KT7TuKG524AYIVF\nV+DcXuey73f2bbfr+5zwtZATPjMzM5sTD7/5MP2G9uOJ954AYLMVNmNA7wH07NGzzY/thK+FnPCZ\nmZnZnGqIBv7y9F/of39/xn8yHoD9vrMf5/Y6lx6L9Wiz4zrhayEnfGZmZtZSH0/7mHMfOpc/PPIH\nps2YRtf5unLSD07ihB+cwELzL9Tqx3PC10JO+MzMzGxujZs8jhPvO5GbXrgJgB6L9eC8Xuexz7f3\nQbUeittCTvhayAmfmZmZtZYHX3+Qo+85mtHjRwOweY/Nubj3xWyywiatsn8nfC3khM/MzMxa04yG\nGVw9+mpOfuBkJnw6AYADNjiAc7Y7h+UXXX6u9u2Er4Wc8JmZmVlbmDptKmePPJsBjw3gyxlfsvD8\nC9N/i/4cu/mxdJ2/a4v26YSvhZzwmZmZWVsa++FYThh2AreOuRWAlRdfmfO3P5891t1jjq/vc8LX\nQk74zMzMrD0MHzeco+85mmfefwaALVfakgG9B7DRchs1ex9O+FrICZ+ZmZm1lxkNM7jyySv57fDf\nMvGziQhx0HcP4uztzmbZRZad7fud8LWQEz4zMzNrbx998RFnPngmlzx+CdMbprPoAoty8pYn069n\nP7rM16XR9znhayEnfGZmZtZRXp70MsffezxDXh4CwKrfWpULfnQBu6+9e83r+5zwtZATPjMzM+to\nw8YO45h7juH5D54HYOtVtuaiHS7iu8t+d6Z6TvhayAmfmZmZdQbTG6Yz+InBnDr8VCZ9PgkhDtno\nEM7a9iy6L9wdcMLXYk74zMzMrDOZ/PlkznjwDAaOGsj0hukstuBinLLVKRy56ZF0mb+LE76WcMJn\nZmZmndGYiWM47t7juOuVuwBYvdvqjO031glfSzjhMzMzs85s6KtDOeaeYxgzcQycjhO+lnDCZ2Zm\nZp3dVzO+4rJ/X0a/nv2c8LWEEz4zMzP7ppjbmzbmac3GmJmZmVnn44TPzMzMrM454TMzMzOrc074\nzMzMzOqcEz4zMzOzOueEz8zMzKzOOeEzMzMzq3NO+MzMzMzqnBM+MzMzszrnhM/MzMyszjnhMzMz\nM6tzTvjMzMzM6pwTPjMzM7M61+kTPklLSxoo6TVJn+X1pZKWqlH3UEljJH2R1307os1mZmZmncl8\nHd2ApkgSMARYDTgFeBFYDzgD2BjYvKruocDlwO+A+4BewCBJiojL27npZmZmZp2GIqKj29AoSWuR\nkry+EfHHqvK+wGXAWhHxiqT5gHeBOyPioKp6VwG7AMtFxPTCvqMzn7uZmZlZhSQiQi19f2cf0p03\nr6cUyiuvK+3fHFgKuLZQ7xpgSWCLNmmdmZmZ2TdAp074IuIF4F7gFEnfk7SIpE2BU4G7IuKlXHW9\nvH6usIsX8nqdtm9tfRgxYkRHN6HTcUxqc1xqc1xqc1xm5ZjU5ri0jU6d8GW7A28Ao4CpwKPAq8Ae\nVXWWyOvJhfd+WNhus+EftFk5JrU5LrU5LrU5LrNyTGpzXNpGuyZ8knpJamjG8kCuPw9wE7Ah0BfY\nCjgc2AS4Kd/UYWZmZmZNaO+7dB8G1m5Gvc/yemfgx8B2ETE8lz0k6TXSUO/OwB183bPXDXi/aj+V\nnr0PMTMzMyupzn6Xbn/gbGDxiPi4qnxxUpL3m4j4vaStgBHA9hFxf1W9rYEHgG0i4sHCvjvviZuZ\nmZkVzM1dup16Hj7g7bzehJS4VWyW1+/k9SPARGA/4P6qevsDk0g9izOZm6CZmZmZfZN09h6+RYDn\nSYnpmcBLpCHh04AvgHUj4rNcty8wiDTx8v3AtsDJwBERcVn7t97MzMysc+jUCR+ApOWB04HtgOWA\n8cAw4PSIeK9Q9zDgOGBl0p29F/kpG2ZmZlZ2nX5aloh4NyIOi4jVI2KhiFgtIvoWk71cd3BErBUR\nXfJ6pmRP0oqSbpL0kaQpkm6WtGL7nU3HktQjP4f4kfxc4gZJK9Wo103SlZI+kPSJpGGSvt0RbW5r\nkvaQdJukN3NMxkj6Xe5drq5XmpgASNpB0gOS3svPpn5L0t8krVOoV6q41CJpaP5ZOrNQXprYSNq6\nkRkXPizUK01MKiTtKGmkpI/z351Rkrap2l6qmEga0cQMHXdX1StVXAAkbZnPc4KkqZKekHRQoU6L\n49LpE77WImkh0nWAawIHAD8H1gCG521l8B/AnqTrGkfWqpCnuhkC/Ag4AvgZMD8pTiu0Uzvb03HA\nV8BvgN6kR/b9EhhWmfanhDGBdMf7KODXwPZAf9IE549W/kkqaVxmImkfYP38MqrKyxqbI4GeVUuv\nyoYyxiRfanQb6WdpN9Lv3xuBhfL20sWE9Pu1Z2E5Nm+7HcoZF0kbkkYv5wEOJs1BPAq4StLhuc7c\nxSUiSrEA/YDpwGpVZauQ/tgf09Hta6cYqOrrQ4AGYKVCnV1z+Q+ryhYjJYkXd/Q5tEFMlqxR9vMc\ng23KGJMmYrVmjsPRjktASorfA/bKcTijalupYgNsnc932ybqlC0mqwCfA0c5JrON1VU5Vt8qa1yA\nc0j3JixUKP8X8K/WiEtpeviAXYBHIuK1SkFEvE66g3fXjmpUe4r86ZiNXYB3omoam4iYSvqvou7i\nFBGTahT/O6+Xz+tSxaQJleG5hrwue1zOA56NiL/V2FbW2DQ1+0HZYvILUidDU9eRly0ms8gjbHsC\nQyLio1xcxrjMS+qA+rxQPpWvf67mKi5lSvjWY9Zn7UJ63u667dyWzqypOK1UkuHvH+b1i3ld2phI\nmlfSApLWAK4gTWx+Q95c5rhsQeoJ/nUjVcoam+skTZc0UdJ1mvka6bLFZAvSzBL7Shor6StJr0j6\nVVWdssWklt2BRYA/V5WVMS5XATOASyQtJ+lbkg4lzThyUa4zV3EpU8LXjVmftQup16JbO7elM1uC\nxuMEdR6rfB3EGcCwiHgyF5c5Jo+RhhleAjYCekXEhLytlHGRtAAp+T0/Il5ppFrZYvMRcAHp2qNt\nSNNo9QIekbR0rlO2mCxPuk7896TpwrYnXaM1UNJRuU7ZYlLLAaR/JO+uKitdXCLiJWAHUm/nO6Rz\nHQj0jYgbc7W5iktnn3jZ2l/nnqenDSndmXs78CVQfWdUaWNCmrx8UWB14HhgqKQtIuINyhuXE4EF\nSU8BakypYhMRo4HRVUX/lDQSeJx0I8epHdKwjjUP6WenT0TclstGSFqFdBPUJR3Urk5Dadq17YAB\nEdFQtalUPz8A+U7bf5AuKbqUNLS7G3CFpGkRcf3cHqNMCd9kame/S+Bn7VabzNfPIK62RNX2uiOp\nK+k6iFVIF8S+W7W5lDEBiIgx+ctRecqE10l3NP+S1KtTqrgoTWN0Mqknq2v+3FR0UXrs4yeU+DNT\nERFPSXqZ9KQkKF9MJpH+URpWKB8G9Ja0LOWLSdH+pMT4z4XyMsblTNLv1J0jYnouGy5pSeBiSX9l\nLuNSpiHd54Fac9WsSxr/tuR50nUCResCb0R+skk9kTQ/cBNpyHLHiHi+UKV0MaklIqYAY0l/xKCc\ncVmN1Lt3LekfxcoCqQd0Mun3TBljU0v1TRxli8nzNH0TS6VOmWJS1AcYHRHPFsrLGJd1gWeqkr2K\nUcCSQHfmMi5lSvjuAHpKWrVSkLvWv5+3WXIHsIKkrSoFkhYDdqYO4yRpHuA60rQSu0XE4zWqlSom\njZG0DOnRhmNz0e2ULy5PkT4r1UtlEt1r8utX8WcGSRuTpvJ5LBeV7fNyS173LpT3Bt6KiPGU+HOS\nPx/rMGvvHpQzLm8DG+QOiGqbkYZ3JzGXP0Od/tFqrSXfvfI0KXC/zcVnAgsD69fpfwyzkLRH/nI7\noC/wK2AiMCEiRuaJHR8CVgROIHUx9yf1WmwQEe+0f6vbjqTLSHE4G7izsPmtiHinbDEBkHQr8ATw\nLGlagDWBY0j/ZW4aEa+WMS6NkdQAnBURp+bXpYqNpGtJie5o0udlQ9L5fgJsFBEfli0mAJLuBzYg\nXQYwjnRB/sHAgRHxlzLGpELSJcDhwPIRMbGwrXRxkbQrcCtwLzCIdLPcLqS/0RdGxPFzHZeOnmyw\nPZccpJuAKaRfSrdQmHi43hfSHGqVZUbV1w9U1elGukV8EvAp6ZqT73R029soHuMKcaheTi1jTPL5\nnki6eHhyPt8xpKeQFCfqLlVcmojXTBMvly02pOs6n85/gL4kPcv8cmCZssYkn++ipDstxwPTSAnx\n3mWOST7n+YEJwO1N1CljXLYnPRFsQs5RniQlxfO0RlxK08NnZmZmVlZluobPzMzMrJSc8JmZmZnV\nOSd8ZmZmZnXOCZ+ZmZlZnXPCZ2ZmZlbnnPCZmZmZ1TknfGZmZmZ1zgmfmc1E0ghJw9tw/w1NLLu0\n1XE7kqQDC+f5iaRxkm6RtGeN+qvkegcUyv9L0puSvpL0VC5bVtIdkibl9xzVXudlZt8c83V0A8ys\n04m8tKWrgStqlL/cxsftaHuQnpm5ILAy8BPgr5IOA3aOiC9yvXeBnsBrlTdK2hQ4C/g9cBvwcd50\nKrAV6UH075GecmFmNhMnfGZWJNo+4XsnIh5v42M0KT+Xcr6I+KodDzs6IipJ3D+BayX9Hfg7KZE7\nCiAivgSK8Vknr6+IiHGF8tERcXtrNFDSghExrTX2ZWadh4d0zUpM0t6Sxkj6QtJzknafg/euKelW\nSe9L+lzSG5JulDRvK7WtQdKZko7Kw59T83DzujXq/lTSo5I+lTQ5t2PFQp3XJV0j6ReSxpCebbpj\n3tZL0lP5PF6RdLCk/5U0Lm9fUNIHki6scezKcO2aLTnPiLgFuB04VFLXvM+ZhnQljSD1igKMzduu\nltQA/BDYqmq4eKX8nlUlXSdpQv7+PiVpt0LbT8/vWU/SPZI+Bv6Wty0k6bwc+2mSXstDyqp6/9b5\n/TtLGphj9EGO8+KFY80n6SRJL+Q4T5B0t6S1quosLelySW/nNr8o6dCWxNXMZuYePrOSktQLuB4Y\nAhwDdAcGkB5sPqYZu7iT9ADvw4GJQA/gx6R/JGfM5r3zSJrl909ETC8U7Z/bciRpGPR84HZJa0fE\njHwehwODgD8BpwOL5fWDktaPiE8quwe2ATYATiM9oPyNnEDeCTwK7JWPcwqweOU8ImKapD8BB0vq\nX+gB6wuMiIi5GY6+G9gN+B7wUI3tv8yx6A/sThq6fQ+4nDQ0Ph34Va47Pie7jwHjgaOBD4C9gZsl\n7RYRQwr7vx24EjgHaMhJ+z2k3sMzgGeBzUlxWQI4vvD+i0mfo32AtUm9lTOAA6vq3ADsClwE3Ad0\nBbYElgNekrRYPvcFSd+fcUBv4LLc6ziwseCZWTNEhBcvXkq4AA8DzxXKNgMagAdm896lcr2dWnDc\nhiaWJQr1XgLmrSr7WS7fPL9eBJgCXFk4xiqkHrx+VWWvA58A3Qt1rwfeB7pUlS0LfAG8VlW2Kimx\n2r+qbP3cnv+czTkfmOut1sj2HfL2Pava3wAcUFXnkFy2UuG9DxW/X8BV+Zy6FcrvBZ6qen163ueR\nhXo/z+VbFMr/K8d1qfx661zv6kK9S4HPq15vm+sd0USMTgE+B1YvlA8mJazzdPTPjBcv3+TFQ7pm\nJZR7cDYGbqouj4jHSIlRdd35qpdcbyLphoLzJB0iaY05bMJV+fjFZUqh3rDIPXnZc3ldGa7dHFgU\nuL7QxrdJyeJWhf09GhETCmU9gbvi6xsmiIjxpISYqrJxpF6vvlXFfUk9hbc0fbqzVRkmba1rJ3sD\ndwFTC3G5F9hA0iKF+rfWeP8bwCOF9w8j9QD3LNS/s/D6OWBBSd3z6x+Rzu2Ps2nzo8DrNdq8JDDL\nUL6ZNZ+HdM3KaSnSH+73a2wrJkRfVr+QtHVEjAS2J/UQnQMsma93Oz8iLm/G8d+LiCebUe/DwuvK\nUGqXvK4kFPc18v5JVV8HaRi0aFlmPWdy2aqFskHAkDwM/AZpmHVQzDoUPacqCWyt9rVEd9Jdu31q\nbAtSAvVJVVnxuN1JdxHXuqGl8v5qs/s+LQl8GE3fDNIdWH0Ojmlmc8AJn1k5TST9YV2mxrZlSNdP\nVWxc2P4y/H+PVx8ASRsARwCDJL0eEUNbvcW1VRK6PsDzNbZ/XHhdqwftPRqPQ9FdpETvcOAZ0pDy\n4Ga1tGk/IQ1nPtEK+4L0/R0JnNfI9mKCV4zLRNJnYJY5ArM5nfplIrCEpC7VPak16owH+jWyvd6n\n7DFrU074zEooImZIGgXsKem/IyIAJG1G6tkZV1V3tj1xEfG0pOOAg4H1gPZK+P5FSurWiIhrWriP\nR4EdJXWNiM8BJC0H/AB4p7piRISkK4CTSEPfw2LmKVLmmKSfATsDA5pIhubUUNJw9wst3OdQ0vWS\nn0bES63QnntIMTsEaOzmi6Gkm3PeiogPWuGYZlbFCZ9ZeZ1Guj7qNkmDgaVJQ7Tj+fqaspokrU+6\nM/MGYCwwL+nGhK+AB5px7B6SiteBAbyer59rloiYKukE4H8kLU1KGqYAK5CmKxkeEX+tNLuR3ZxF\nmhD5HkkXkIYhTyHFoaFG/StJcdoA+Glz25ptmK9rWwBYCdgpH/te0h24LVU8t1NJ8/iNlDSQ1CPX\nDfg2sGpEHDyb/V0HHATcL+kPpN7MBUhDrjsDu1WS4+aIiBGSbgYuzHcQDyddUrAV8I+IeJB09+5e\nwD8lXUTq0VuYdNfvFhGxW+29m1lzOOEzK6mIuF/SfqTk5WbgFdJw2tHM/uaByhMdjiVNx/IFKSnY\nKSKeasbh+zDzlB0VxwOzzHXXlIgYLOkt4ARgX9LvtXdIQ5rVbal5ThHxoqSfkKZ8uZF0w8d5pClm\nVq5Rf6KkkaSbCO5objPz+u95/QXpGsEngL0i4uY53E+xbKbyiHhL0sak7+3vSMn8JNL0Kn9u6r35\n/dMl7QD8BjiMdC3jp8CrpBs0vizsozlt3ZvUy9eH9BmbQkpKB+djTpX0fVKyehIpaf+INC1Pc+Nj\nZo1QHskxM7Ms38X6KjAkIg4tbOsGvAlcGBGndUT7zMzmlHv4zKz0JF1Kuh7wXWB5Uk/n4qRh60qd\npUjDi5WbCga1czPNzFrMCZ+ZWXq6w7mkO3O/JD2loldEPFdVZyfS0zzeAPpERK0pbczMOiUP6ZqZ\nmZnVOT9pw8zMzKzOOeEzMzMzq3NO+MzMzMzqnBM+MzMzszrnhM/MzMyszjnhMzMzM6tz/wcgev8H\nyo2ELQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x73fc870>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(10,6))\n",
"plt.plot(ediff, s2d4_ave, lw=2, label='[Ar]4$s^2$3$d^4$')\n",
"plt.plot(ediff, s1d5_ave, lw=2, label='[Ar]4$s^1$3$d^5$')\n",
"plt.ylabel('Average Configuration Energy')\n",
"plt.xlabel('d-s Energy Difference')\n",
"plt.title('Energy Comparison for\\nTwo Cr e- Configurations')\n",
"plt.legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.4.3"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment