{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Energy volume curve\n", "## Theory \n", "By fitting the energy-volume curve, we can determine several equilibrium properties, including the equilibrium energy $E_0$, equilibrium volume $V_0$, equilibrium bulk modulus $B_0$, and its derivative $B_0'$. These properties are then utilized in the Einstein model to provide an initial estimation for the thermodynamic properties, namely the heat capacity $C_v$ and the free energy $F$." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Initialisation" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We start by importing matplotlib, numpy and the pyiron project class. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2019-09-04T12:58:08.037567Z", "start_time": "2019-09-04T12:58:06.149845Z" }, "tags": [] }, "outputs": [ { "data": { "text/html": [ "
\n" ], "text/plain": [] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "\n" ], "text/plain": [] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "dd260ca2d4ea4a9296c26b7e12a1ce1a", "version_major": 2, "version_minor": 0 }, "text/plain": [] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from pyiron import Project\n", "%matplotlib inline " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In the next step we create a project, by specifying its relative path/name." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "ExecuteTime": { "end_time": "2019-09-04T12:58:09.049344Z", "start_time": "2019-09-04T12:58:08.040033Z" }, "tags": [] }, "outputs": [], "source": [ "pr = Project('thermo')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Atomistic structure" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To analyse the energy volume dependence a single super cell is sufficient, so we create an iron super cell as an example." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "ExecuteTime": { "end_time": "2019-09-04T12:58:09.222395Z", "start_time": "2019-09-04T12:58:09.051675Z" }, "tags": [] }, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "122085497932408fbab5cda44864e959", "version_major": 2, "version_minor": 0 }, "text/plain": [ "NGLWidget()" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "basis = pr.create.structure.bulk(name='Fe', cubic=True, a=2.75)\n", "basis.plot3d()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Measure energies by looping over different volumes\n", "\n", "Let's first see how we can get the energy as a function of the volume mainly based on the knowledge you acquired in the [first tutorial notebook](https://pyiron.readthedocs.io/en/latest/source/notebooks/first_steps.html). For this, we are going to vary the volume of our initial structure by up to 5 %.\n", "\n", "Energy volume curves are commonly calculated with ab-initio simulation codes. In this example, we use [GPAW](https://wiki.fysik.dtu.dk/gpaw/), but the functionality below would work just as well with any DFT code, just by changing the job name prefix `Gpaw` to the DFT code of your choice (e.g. `Vasp`, `Sphinx`...). In this tutorial, we select an energy cutoff of 320 eV, but the simulation would run without changing any parameter. For the rest of the DFT parameters, you can look up functions starting with `job.set_*`." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Launch calculations" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "tags": [] }, "outputs": [], "source": [ "# Exponent of 1/3 to translate volume strain of 5 % to line strain\n", "line_strain_list = np.linspace(0.95, 1.05, 7)**(1 / 3) - 1" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The job sample_m0d01695243 was saved and received the ID: 19694133\n", "The job sample_m0d01123691 was saved and received the ID: 19694134\n", "The job sample_m0d00558671 was saved and received the ID: 19694135\n", "The job sample_0d0 was saved and received the ID: 19694136\n", "The job sample_0d00552497 was saved and received the ID: 19694137\n", "The job sample_0d01098989 was saved and received the ID: 19694138\n", "The job sample_0d01639636 was saved and received the ID: 19694139\n" ] } ], "source": [ "for strain in line_strain_list:\n", " dft = pr.create.job.Gpaw(job_name=('sample', strain))\n", " dft.set_encut(320.0) # Optional, among other plane wave dft parameters\n", " dft.structure = basis.apply_strain(strain, return_box=True)\n", " dft.run()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If your environment allows you to launch calculations to a queueing system, it might be useful to set\n", "\n", "```python\n", "dft.server.queue = \"name_of_the_queue\"\n", "dft.server.cores = number_of_cores\n", "dft.server.run_time = run_time_in_seconds\n", "```\n", "\n", "Queueing system setup can be found in the pyiron docs, the homepage for which is [here](https://pyiron.readthedocs.io/en/latest/index.html)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Analyse results" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "tags": [] }, "outputs": [], "source": [ "results = {'energy': [], 'volume': []}\n", "for strain in line_strain_list:\n", " dft = pr.load(('sample', strain))\n", " results['volume'].append(dft.structure.get_volume())\n", " results['energy'].append(dft.output.energy_pot[-1])" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "tags": [] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlQAAAG1CAYAAAA/aGqTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABpi0lEQVR4nO3dd3gU1f7H8femJ4SEEiDBRFpooQYpUqSICaggWBAUkCZXuSIgP+DKla4C0hTsXpGAevGqYEExEAEVpZeAFEM3CEF6ChCyZOf3x5iFmAAJKZtNPq/n2YfdmTOz39nDJt+cc+Yci2EYBiIiIiJyy1wcHYCIiIiIs1NCJSIiIpJHSqhERERE8kgJlYiIiEgeKaESERERySMlVCIiIiJ5pIRKREREJI+UUImIiIjkkZujAygpbDYbx48fp3Tp0lgsFkeHIyIiIjlgGAbJyclUrlwZF5frt0MpoSokx48fJyQkxNFhiIiIyC04evQowcHB192vhKqQlC5dGjArxM/Pz8HRFG1Wq5WVK1cSGRmJu7u7o8ORm1B9OR/VmXNRfTlWUlISISEh9t/j16OEqpBkdPP5+fkpoboJq9WKj48Pfn5++uHhBFRfzkd15lxUX0XDzYbraFC6iIiISB4poRIRERHJIyVUIiIiInmkMVQiIk4oPT0dq9V6S8darVbc3NxITU0lPT09nyOT/Kb6Klju7u64urrm+TxKqEREnIhhGJw4cYLz58/n6RyBgYEcPXpU8+I5AdVXwStTpgyBgYF5+nyVUImIOJGMZKpixYr4+Pjc0i8Am81GSkoKvr6+N5yoUIoG1VfBMQyDixcvcvLkSQCCgoJu+VxKqEREnER6ero9mSpfvvwtn8dms5GWloaXl5d+QTsB1VfB8vb2BuDkyZNUrFjxlrv/VDMiIk4iY8yUj4+PgyMRKV4yvlO3Oi4RlFCJiDgdjaMRyV/58Z1Sl58TS7cZbDp8lpPJqVQs7UXzauVwddEPWhERkcKmhMpJRe9KYPKyPSQkptq3Bfl7MbFrGJ3r3/qgOhERyR9RUVGMGDEiT3dkivNQl58Tit6VwJCPtmVKpgBOJKYy5KNtRO9KcFBkIiLZ69+/PxaLBYvFgru7O5UqVSIiIoIPPvgAm82WqWzVqlXtZTMewcHBTJo0Kcv2vz+OHDnimAvMgUmTJtG4cWNHhyEFRAmVk0m3GUxetgcjm30Z2yYv20O6LbsSIiKmdJvBhkNn+Cr2GOsPnimUnxmdO3cmISGBI0eO8N1339GhQweGDx9Oly5duHLlSqayU6ZMISEhwf7Yvn07o0aNyrQtODg4S7mQkJCbxmEYRpb3E8krJVROZtPhs5laptzTrXTfvYb2BzcDZlKVkJjKpsNnHRShiBR10btOcN/bW3j8/U0M/ySWx/6zgTavrC7w1m1PT08CAwO57bbbaNKkCf/+97/56quv+O6774iKispUtnTp0gQGBtofFSpUwNfXN9M2V1fXLOWyu+X9hx9+wGKxsGLFCpo2bYqnpydr167FMAxmzJhB9erV8fb2plGjRnz++ef2486dO0fv3r2pUKEC3t7e1KxZkwULFmQ657XdebGxsddtJYuKimLy5Mns2LHD3pqWcc2TJk3i9ttvx9PTk8qVKzNs2LA8f9ZS+DSGysmcTM7czffE1m8Yv2Y+eypW44fqTeGvOxX+Xk5EBMwhA8/8d3uWVu6MIQNv92lSqOMw7777bho1asTSpUt58sknC/S9xowZw6xZs6hevTplypRh3LhxLF26lLfffpuaNWvy008/0adPHypUqEC7du0YP348e/bs4bvvviMgIIADBw5w6dKlW3rvnj17smvXLqKjo/n+++8B8Pf35/PPP+fVV1/lk08+oV69epw4cYIdO3bk52VLIVFC5WQqlvbK9PrzBvfwfz9/RNjJw7SM38n6Ko2yLScicrMhAxbMIQMRYYGFesdwnTp12LlzZ6Zt//rXvxg3bpz99dSpU/PccjNlyhQiIiIAuHDhAnPmzGH16tW0bNkSgOrVq/Pzzz/z7rvv0q5dO+Lj4wkPD6dp06aAObbrVnl7e+Pr64ubmxuBgYH27fHx8QQGBnLPPffg7u7O7bffTvPmzW/9IsVh1OXnZJpXK0eQvxcZP+oSvUvzWYN7ABi86QssmHf7Na9WzmExikjR9PchA3/nqCEDhmFkmQdo9OjRxMbG2h9PPPFEnt8nIzEC2LNnD6mpqURERODr62t/LFq0iIMHDwIwZMgQPvnkExo3bsyYMWNYt25dnmP4ux49enDp0iWqV6/O4MGD+eKLLzS+y0kpoXIyri4WJnYNA7AnVR807YYNC3cf2kKN00eZ2DVM81GJSBY5HQpQ2EMG9u7dS7Vq1TJtCwgIIDQ01P4oU6ZMnt+nVKlS9ucZdxZ+++23mRK3PXv22MdR3Xvvvfz++++MGDGC48eP07FjR0aNGgVgXwLGMK62993KLNshISHExcXx5ptv4u3tzT//+U/atm2bpxm7xTGUUDmhzvWDeLtPEwL9zW6938tWJqZmCwAWnP9Z81CJSLZyOhSgMIcMrF69ml9//ZWHH3640N4TICwsDE9PT+Lj4zMlbqGhoZnuFKxQoQL9+/fno48+4rXXXuO9996zbwdISLg6kD82NvaG7+nh4UF6enqW7d7e3jzwwAPMmzePH374gfXr1/Prr7/mw1VKYdIYKifVuX4QEWGB9pnSb2v8AjzWlZBvPodTr8JfX3YRkQwZQwZOJKZmO47KAgQW4JCBy5cvc+LECdLT0/nzzz+Jjo5m2rRpdOnSJV+69HKjdOnSjBo1iueeew6bzUabNm1ISkpi3bp1+Pr60q9fPyZMmMAdd9xBvXr1uHz5Mt988w1169YFsCdekyZN4qWXXmL//v3Mnj37hu9ZtWpVDh8+TGxsLMHBwZQuXZrFixeTnp5OixYt8PHx4cMPP8Tb25sqVaoUxscg+UgJlRNzdbHQssZfK843qgxzmkGpUnDmjBIqEckiY8jAkI+2YYFMSVXGIIGCHDIQHR1NUFAQbm5ulC1blkaNGjFv3jz69etn70IrTC+++CIVK1Zk2rRpHDp0iDJlytincwCzRWns2LEcOXIEb29v7rrrLj755BMA3N3dWbx4MUOGDKFRo0Y0a9aMl156iR49elz3/R5++GGWLl1Khw4dOH/+PAsWLKBMmTJMnz6dkSNHkp6eToMGDVi2bBnly5cvlM9A8o/FuLYDWApMUlIS/v7+JCYm4ufnVzBvkpICvr4Fc+5CZLVaWb58Offddx/u7u6ODkduQvVVeFJTUzl8+DDVqlXDy+vWu+WW7zzO5GW7+TM5zb5NS1cVXTabjaSkJPz8/BySeJYEN/pu5fT3t1qoipNikEyJSMHrXD+QFsHexJ29wqmUNC2uLpIPlFAVRydPwuLFMGyYfaJPEZFrubpYuLN6ebV4iOQTJVTFTVoa1K8Pp05B7drQubOjIxIRESn29KdJcePhAX36mM/nzHFsLCIiIiWEEqriaNgwcHGBmBj423IOIiIikv+UUBVHVavCI4+Yz1991aGhiIiIlARKqIqrkSPNfz/+GK6ZyVdERETynxKq4qpFC2jVCqxWePNNR0cjIiJSrCmhKs7+7//A09NMqkREiomoqKh8WSw5J44cOYLFYrnpOn2O8sMPP2CxWDh//ryjQ7mpW6m3qlWr8tprr92wzKRJk6hUqRIWi4Uvv/yS/v37071791uO81YpoSrOunWD+Hh45RVHRyIiJdzbb79Nw4YN8fPzw8/Pj5YtW/Ldd985OiwpRD179mTfvn35es69e/cyefJk3n33XRISErj33nuZO3cuUVFR9jLt27dnxIgR+fq+2dE8VMWZqytUrOjoKERECA4OZvr06YSGhgKwcOFCunXrxvbt26lXr162x6SlpeHh4VEo8RXkexXmdRRl3t7eeHt75+s5Dx48CEC3bt2w/DWRtaenZ76+R06phaqk2LoV1q1zdBQiUkJ17dqV++67j1q1alGrVi1efvllfH192bBhg71M1apVeemll+jfvz/+/v4MHjwYMLuKbr/9dnx8fHjwwQc5c+bMTd/v119/5e6778bb25vy5cvzj3/8g5SUFPv+jG6hadOmUblyZWrVqgXApk2bCA8Px8vLi6ZNm7J9+/Ys596zZw/33Xcfvr6+VKpUib59+3L69Gn7/vbt2zN06FBGjhxJQEAAERER2caYEcPkyZOpWLEifn5+PPXUU6SlXV1j8fLlywwfPpyaNWvi4+NDmzZt2Lx5c7bnu3DhAn5+fnz++eeZti9btoxSpUqRnJxs78LMWKTZx8eHRo0asX79+kzHLFmyhHr16uHp6UnVqlWZPXt2pv0ZdfXEE0/g6+tLlSpV+Oqrrzh16hTdunXD19eXBg0asGXLFvsxf+/yO3jwIN26daNSpUr4+vrSrFkzvv/++2yvLTuTJk2ia9euALi4uNgTqmu7/Pr378+PP/7I3LlzsVgsWCwWjhw5kuP3yA0lVCXBggXQtKk5P5XWwhYRB0tPT+eTTz7hwoULtGzZMtO+mTNnUr9+fbZu3cr48ePZuHEjAwcO5J///CexsbF06NCBl1566Ybnv3jxIp07d6Zs2bJs3ryZzz77jO+//56hQ4dmKrdq1Sr27t1LTEwM33zzDRcuXKBLly7Url2brVu3MmnSJEaNGpXpmISEBNq1a0fjxo3ZsmUL0dHR/Pnnnzz66KOZyi1cuBA3Nzd++eUX3n333evGmhHDmjVrWLx4MV988QWTJ0+27x8zZgxLly7lrbfeYsuWLYSGhtKpUyfOnj2b5VylSpWiV69eLFiwINP2BQsW8Mgjj1C6dGn7thdeeIFRo0YRGxtLrVq1eOyxx7hy5QoAW7du5dFHH6VXr178+uuvTJo0ifHjx2fqRgN49dVXad26Ndu3b+f++++nb9++PPHEE/Tp04dt27YRGhrKE088gXGd3zspKSncd999fP/992zfvp1OnTrRtWtX4uPjr/t5XWvUqFH2a01ISCAhmzva586dS8uWLRk8eLC9TEhISI7On2uGFIrExEQDMBITEwv/zU+eNAxvb8MAw1i9uvDfP5fS0tKML7/80khLS3N0KJIDqq/Cc+nSJWPPnj3GpUuXsu5MSbn+42/l05OSjHN//GGkJyVlLXvxYs7Oewt27txplCpVynB1dTX8/f2Nb7/9NtP+KlWqGN27d8+07bHHHjM6d+6caVvPnj0Nf3//677Pe++9Z5QtW9ZIuSbOb7/91nBxcTFOnDhhGIZh9OvXz6hUqZJx+fJle5l3333XKFeunHHhwgX7trffftsAjO3btxuGYRjjx483IiMjM73f0aNHDcCIi4szDMMw2rVrZzRu3Pgmn4YZQ3bv5+vra6SnpxspKSmGu7u78eGHHxrnzp0z0tPTjbS0NKNy5crGjBkzDMMwjDVr1hiAce7cOcMwDGPjxo2Gq6urcezYMcMwDOPUqVOGu7u78cMPPxiGYRiHDx82AOP999+3v+fu3bsNwNi7d69hGIbx+OOPGxEREZliHT16tBEWFmZ/XaVKFaNPnz721wkJCQZgjB8/3r5t/fr1BmAkJCQYhmEYCxYsuGG9GYZhhIWFGa+//nqm93n11VevW/6LL74w/p7K9OvXz+jWrZv9dbt27Yzhw4ff8H1v9N3K6e9vtVCVBBUqwMCB5vMZMxwbi4gUDF/f6z8efjhTUUtgIGWCg3Hx88ta9t57M5+3atXsz3kLateuTWxsLBs2bGDIkCH069ePPXv2ZCrTtGnTTK/37t2bpRXr76//bu/evTRq1IhSpUrZt7Vu3RqbzUZcXJx9W4MGDTKNbco4zsfH57rvtXXrVtasWYOvr6/9UadOHeDqeJ7sruN6snu/lJQUjh49ysGDB7FarbRu3dq+393dnebNm7N3795sz9e8eXPq1avHokWLAPjwww+5/fbbadu2baZyDRs2tD8PCgoC4OTJk/bP4dr3BPPz279/P+np6dmeo1KlSoD5mf59W8Z5/+7ChQuMGTOGsLAwypQpg6+vL7/99luOW6iKGiVUJcXIkeZyNNHRWo5GRBzCw8OD0NBQmjZtyrRp02jUqBFz587NVObaJAi4bnfRjRiGYR9P83fXbr+V97LZbHTt2pXY2NhMj/3792dKWv5+7tyyWCz2eP5+LTe6PoAnn3zS3hW2YMECBgwYkKW8u7t7pvfKuLbrnT+7zya7c9zovH83evRolixZwssvv8zatWuJjY2lQYMGmcaQORMlVCVF9erQo4f5fOZMx8YiIvkvJeX6jyVLMhU1Tpzg/B9/YEtKylr271MZHDmS/TnzgWEYXL58+YZlwsLCMg1cB7K8zu6Y2NhYLly4YN/2yy+/4OLiYh98fr3jduzYwaVLl677Xk2aNGH37t1UrVqV0NDQTI9bSaKyez9fX1+Cg4MJDQ3Fw8ODn3/+2b7farWyZcsW6tate91z9unTh/j4eObNm8fu3bvp169frmIKCwvL9J4A69ato1atWri6uubqXDeydu1a+vfvz4MPPkiDBg0IDAwskAHjHh4emVrWCooSqpJk9Gjz38WL4fffHRuLiOSvUqWu//DyynnZv9/Wfr1yufTvf/+btWvXcuTIEX799VdeeOEFfvjhB3r37n3D44YNG0Z0dDQzZsxg3759vPHGG0RHR9/wmN69e+Pl5UW/fv3YtWsXa9as4dlnn6Vv3772bqjsPP7447i4uDBo0CD27NnD8uXLmTVrVqYyzzzzDGfPnuWxxx5j06ZNHDp0iJUrVzJw4MBb+qWdlpZmf7/vvvuOiRMnMnToUFxcXChVqhRDhgzhX//6F99//z179uxh8ODBXLx4kUGDBl33nGXLluWhhx5i9OjRREZGEhwcnKuY/u///o9Vq1bx4osvsm/fPhYuXMgbb7yRZYB+XoWGhrJ06VJiY2PZsWMHjz/++HVbs/KiatWqbNy4kSNHjnD69OkCeQ9QQlWy3HEHdOwIgYFw6JCjoxGREuTPP/+kb9++1K5dm44dO7Jx40aio6OvO6VAhjvvvJP333+f119/ncaNG7Ny5UrGjRt3w2N8fHxYsWIFZ8+epVmzZjzyyCN07NiRN95444bH+fr6smzZMvbs2UN4eDgvvPACr/xtYuTKlSvzyy+/kJ6eTqdOnahfvz7Dhw/H398fF5fc/0rt2LEjNWvWpG3btjz66KN07dqVSZMm2fdPnz6dhx56iKeffpqmTZty4MABVqxYQdmyZW943kGDBpGWlsbAjPGzudCkSRM+/fRTPvnkE+rXr8+ECROYMmUK/fv3z/W5buTVV1+lbNmytGrViq5du9KpUyeaNGmSr+8B5t2Arq6uhIWFUaFChQIbo2UxbqWDWnItKSkJf39/EhMT8fPzc1wgx46Zg9SL8CRzVquV5cuXc99992Xqj5eiSfVVeFJTUzl8+DDVqlXD6++tTrlgs9lISkrCz8/vlpIAyR/9+/fn/PnzfPnllzcsdyv19fHHHzN8+HCOHz+uSUVz4EbfrZz+/tZM6SXNbbc5OgIRESkgFy9e5PDhw0ybNo2nnnpKyVQh0p8mJdWVK/DZZ3DNYEgREXFuM2bMoHHjxlSqVImxY8c6OpwSRQlVSRUZCY8+CgsXOjoSEZESKSoq6qbdfbk1adIkrFYrq1atwvcW5wuTW6OEqqT6a50jZs2CQridVEREpDhTQlVSDRoE5crBwYPwxReOjkZERMSpKaEqqUqVgmeeMZ/PmKFFk0WciG7OFslf+fGdUkJVkg0dak74t3kz/Pijo6MRkZvImJbi4sWLDo5EpHjJ+E7lZeoXp5k24eWXX+bbb78lNjYWDw8Pzp8/n225qKgo5syZw759+yhTpgyPPPLIDSdza9++PT/+LZno2bMnn3zySZayly9fpkWLFuzYsYPt27fTuHHjvFyS41WsCAMGwNtvw/Tp0L69oyMSkRtwdXWlTJky9sVmfXx8brim2/XYbDbS0tJITU3VPFROQPVVcAzD4OLFi5w8eZIyZcrkaWkdp0mo0tLS6NGjBy1btmT+/PnZlpkzZw6zZ89m5syZtGjRgtTUVA7lYEbwwYMHM2XKFPtr778vvfCXMWPGULlyZXbs2HFrF1EUjRoF770HFy+aj2tWPReRoicwMBDAnlTdCsMwuHTpEt7e3reUkEnhUn0VvDJlyti/W7fKaRKqyZMnA2YLVHbOnTvHuHHjWLZsGR07drRvr1ev3k3P7ePjc9MP8rvvvmPlypUsWbKE7/6+eKgzq14d9uyBGywYKiJFh8ViISgoiIoVK2K1Wm/pHFarlZ9++om2bdtqdnsnoPoqWO7u7vmy6LPTJFQ3ExMTg81m49ixY9StW5fk5GRatWrF7NmzCQkJueGxH3/8MR999BGVKlXi3nvvZeLEiZQuXdq+/88//2Tw4MF8+eWX+OSwBefy5cuZVlFPSkoCzC/Grf4QLDDVqkERiinj8ylyn5NkS/XlOLf6S8Bms3HlyhVcXV3z5ReJFCzVV8Gy2Ww3XDA5pz/bik1CdejQIWw2G1OnTmXu3Ln4+/szbtw4IiIi2Llz53Wn3+/duzfVqlUjMDCQXbt2MXbsWHbs2EFMTAxgNrX279/fvjDlkSNHchTPtGnT7K1q11q5cmWOk7LC5p6URMCePSTceaejQwGw14E4B9WX81GdORfVl2Pk9CYQhyZUkyZNyjbpuNbmzZtp2rTpTc9ls9mwWq3MmzePyMhIABYvXkxgYCBr1qyhU6dO2R43ePBg+/P69etTs2ZNmjZtyrZt22jSpAmvv/46SUlJuZ7Cf+zYsYwcOdL+OikpiZCQECIjIx27OPL1HDuGW/36cPkyV/buhSpVHBaK1WolJiaGiIgINW87AdWX81GdORfVl2Nl9DDdjEMTqqFDh9KrV68blqlatWqOzhUUFARAWFiYfVuFChUICAggPj4+xzE1adIEd3d39u/fT5MmTVi9ejUbNmzA09MzU7mmTZvSu3dvFl5n6RZPT88sx4DZV1skvxBVq8Kdd8KqVbjPnQuvv+7oiIruZyXZUn05H9WZc1F9OUZOP3OHJlQBAQEEBATky7lat24NQFxcHMHBwQCcPXuW06dPUyUXrS27d+/GarXaE7R58+bx0ksv2fcfP36cTp068b///Y8WLVrkS+xFxr//DatWwfvvw7hxUKmSoyMSERFxCk4zoUV8fDyxsbHEx8eTnp5ObGwssbGxpKSkAFCrVi26devG8OHDWbduHbt27aJfv37UqVOHDh06AHDs2DHq1KnDpk2bADh48CBTpkxhy5YtHDlyhOXLl9OjRw/Cw8PtCdrtt99O/fr17Y9af90NV6NGDXviVmx06AAtWkBqKrz2mqOjERERcRpOk1BNmDCB8PBwJk6cSEpKCuHh4YSHh7NlyxZ7mUWLFtGiRQvuv/9+2rVrh7u7O9HR0fbmOqvVSlxcnH2AmYeHB6tWraJTp07Url2bYcOGERkZyffff18y76SwWMxWKoC33oLrTJ4qIiIimTnNXX5RUVHXnYMqg5+fH/Pnz7/uxJ9Vq1bNtF5PSEhIllnSb+bv5yh2unSB+vVh1y4zqcpIsEREROS6nKaFSgqJiwuMHQvu7mqhEhERySGnaaGSQvToo9C2LRS3MWIiIiIFRC1UkpWbm5IpERGRXFBCJTe2YwesWePoKERERIo0JVRyfUuWQOPG8NRTkJ7u6GhERESKLCVUcn2RkVC2LOzfbyZXIiIiki0lVHJ9pUvDsGHm85dfhhusxi0iIlKSKaGSGxs2DHx9YedOWLbM0dGIiIgUSUqo5MbKlYNnnzWfT5kCxXlSUxERkVukhEpubuRIKFUKtm2D775zdDQiIiJFjhIqubmAAPjnP6FSJUhOdnQ0IiIiRY4SKsmZcePg8GHo2dPRkYiIiBQ5WnpGcsbPz9ERiIiIFFlqoZLcSU+HTz+F9esdHYmIiEiRoYRKcufll81uv7FjHR2JiIhIkaGESnJn4EDw8IAffzQfIiIiooRKcik42EyqAF580bGxiIiIFBFKqCT3nn8e3Nxg1SpYt87R0YiIiDicEirJvSpVoF8/87laqURERJRQyS0aOxZcXSE6GjZtcnQ0IiIiDqV5qOTW1KgBffrAb79pfT8RESnxlFDJrXvrLfD2BovF0ZGIiIg4lBIquXU+Po6OQEREpEjQGCrJu8REmDIFtm51dCQiIiIOoRYqybsxY+C992DjRvj2W0dHIyIiUujUQiV5N3q0ecff8uWwYYOjoxERESl0Sqgk70JD4YknzOeTJjk0FBEREUdQQiX5Y9w4c/b0FSvgl18cHY2IiEihUkIl+aN6dRgwwHw+caJjYxERESlkSqgk/7zwAri7m2v8/fSTo6MREREpNLrLT/JPlSrw1FOQlAQhIY6ORkREpNAooZL8NW+eZk4XEZESR11+kr+UTImISAmkhEoKxm+/weOPQ0yMoyMREREpcOryk4Lx7ruweDEcPgz33KOWKxERKdbUQiUF41//Am9vc+b06GhHRyMiIlKglFBJwQgMhGeeMZ9PmACG4dh4RERECpASKik4Y8ZAqVKwZQt89ZWjoxERESkwSqik4FSoACNGmM/HjYP0dIeGIyIiUlCUUEnBGjUKypSB3bvhf/9zdDQiIiIFQnf5ScEqUwZefBGSk+GBBxwdjYiISIFQQiUFb+hQR0cgIiJSoNTlJ4XLZgOr1dFRiIiI5CslVFJ41qyBpk1h7lxHRyIiIpKvlFBJ4fn9d9i+HaZNg8RER0cjIiKSb5RQSeHp0wfq1IGzZ2HOHEdHIyIikm+UUEnhcXMz7/gDM6E6dcqx8YiIiOQTJVRSuB56CJo0gZQUmD7d0dGIiIjkCyVUUrhcXGDqVPP5m2/CH384Nh4REZF8oIRKCl9kJLRtC5cvw4IFjo5GREQkzzSxpxQ+iwVmzYKjR+HBBx0djYiISJ4poRLHaNbMfIiIiBQD6vITx0tKgvh4R0chIiJyy5RQiWOtXAk1asCgQY6ORERE5JblOKG6cuVKQcYhJVWtWmYL1fffm8mViIiIE8pxQhUUFMSoUaPYu3dvQcYjJU3VqvDPf5rP//Uvc/FkERERJ5PjhGrkyJEsW7aM+vXr07JlS+bPn09KSkpBxiYlxQsvgJ8fxMbC4sWOjkZERCTXcpxQjR07lri4OH744Qfq1KnDiBEjCAoKYsCAAfzyyy8FGaMUdwEB8Pzz5vNx48z5qURERJxIrgel33XXXSxYsIATJ07w2muvceDAAe666y5q167NjBkzCiJGAF5++WVatWqFj48PZcqUuW65qKgoGjZsiJeXF4GBgQwdOvSG523fvj0WiyXTo1evXlnKffvtt7Ro0QJvb28CAgJ46KGH8npJcq3hw6FyZThyBJd33nF0NCIiIrlyy3f5lSpVikGDBrF27VqWLVvG6dOnGTt2bH7GlklaWho9evRgyJAh1y0zZ84cXnjhBZ5//nl2797NqlWr6NSp003PPXjwYBISEuyPd999N9P+JUuW0LdvXwYMGMCOHTv45ZdfePzxx/N8TXINHx+YPNmc9FNTKIiIiJO55Yk9L168yP/+9z8WLFjAL7/8Qo0aNRg9enR+xpbJ5MmTAbMFKjvnzp1j3LhxLFu2jI4dO9q316tX76bn9vHxITAwMNt9V65cYfjw4cycOZNB19zaX7t27VxELznSvz+0aIGtTh1YvtzR0YiIiORYrluo1q5dy8CBA+3dadWqVWPNmjXs27eP5zPGwThATEwMNpuNY8eOUbduXYKDg3n00Uc5evToTY/9+OOPCQgIoF69eowaNYrk5GT7vm3btnHs2DFcXFwIDw8nKCiIe++9l927dxfk5ZRMbm7QoIGjoxAREcm1HLdQTZ06laioKA4ePEjTpk2ZOXMmjz32GH5+fgUZX44dOnQIm83G1KlTmTt3Lv7+/owbN46IiAh27tyJh4dHtsf17t2batWqERgYyK5duxg7diw7duwgJibGfl6ASZMmMWfOHKpWrcrs2bNp164d+/bto1y5ctme9/Lly1y+ZnB1UlISAFarFavVmp+XXuxYrVZKJSRg+89/sA4e7Ohw5CYy/j/r/7XzUJ05F9WXY+X0c89xQvXqq6/Sp08fBg0aRP369W85sGtNmjTJ3pV3PZs3b6Zp06Y3PZfNZsNqtTJv3jwiIyMBWLx4MYGBgaxZs+a6Y6kGX/MLu379+tSsWZOmTZuybds2mjRpgu2veZFeeOEFHn74YQAWLFhAcHAwn332GU899VS25502bVq217Zy5Up8fHxuej0lmffJk9zz7LNYbDZW22ykhIQ4OiTJgYw/QsR5qM6ci+rLMS5evJijcjlOqI4fP467u3umbampqXh5eeUusmsMHTo02zvqrlW1atUcnSsoKAiAsLAw+7YKFSoQEBBAfC4GOTdp0gR3d3f2799PkyZNsj2vp6cn1atXv+F5x44dy8iRI+2vk5KSCAkJITIyssi06hVVVquVP+fPJ2jjRjqsXEn6kiWODkluwGq1EhMTQ0RERJafEVI0qc6ci+rLsTJ6mG4mxwlVRiXabDZefvll3nnnHf7880/27dtH9erVGT9+PFWrVs00cPtmAgICCAgIyHH5G2ndujUAcXFxBAcHA3D27FlOnz5NlSpVcnye3bt3Y7Va7YnUHXfcgaenJ3FxcbRp0wYw/3MfOXLkhuf19PTE09Mzy3Z3d3d9IXJgT58+BG7ejMuyZbhs3Ah/ffZSdOn/tvNRnTkX1Zdj5PQzz/Wg9JdeeomoqChmzJiRaVxSgwYNeP/993N7uhyLj48nNjaW+Ph40tPTiY2NJTY21j5be61atejWrRvDhw9n3bp17Nq1i379+lGnTh06dOgAwLFjx6hTpw6bNm0C4ODBg0yZMoUtW7Zw5MgRli9fTo8ePQgPD7cnaH5+fjz99NNMnDiRlStXEhcXZ5+6oUePHgV2vSVdSkgItoEDzRf/939akkZERIq0XCdUixYt4r333qN37964urratzds2JDffvstX4O71oQJEwgPD2fixImkpKQQHh5OeHg4W7ZsyRRbixYtuP/++2nXrh3u7u5ER0fbs0ur1UpcXJy9P9TDw8M+V1Xt2rUZNmwYkZGRfP/995mubebMmfTq1Yu+ffvSrFkzfv/9d1avXk3ZsmUL7HoFbBMmgK8vbNoE//ufo8MRERG5rlzPQ3Xs2DFCQ0OzbM8YFF5QoqKirjsHVQY/Pz/mz5/P/Pnzs91ftWpVDMOwvw4JCeHHH3+86Xu7u7sza9YsZs2alauYJY8CA80lacaNM//t3h28vR0dlYiISBa5bqGqV68ea9euzbL9s88+Izw8PF+CErEbORLCwuAf/zBnURcRESmCct1CNXHiRPr27cuxY8ew2WwsXbqUuLg4Fi1axDfffFMQMUpJ5u0NO3fCNV2wIiIiRU2uW6i6du3K//73P5YvX47FYmHChAns3buXZcuWERERURAxSkl3bTJ1TZetiIhIUXFLa/l16tQpR4sOi+SrVavgX/+ChQshB2s0ioiIFJZct1CJOMybb8LWrVCAi3CLiIjcCiVU4jxeeQXc3eG772DlSkdHIyIiYqeESpxHzZrwzDPm81GjID3dsfGIiIj8RQmVOJfx46FsWfj1V1iwwNHRiIiIALeQUP3www8FEIZIDpUrBxMmmM/HjYPkZMfGIyIiwi0kVJ07d6ZGjRq89NJLHD16tCBiErmxf/4TQkPhzz/hww8dHY2IiEjuE6rjx48zfPhwli5dSrVq1ejUqROffvopaWlpBRGfSFYeHuYdf//7H/y1ULWIiIgj5TqhKleuHMOGDWPbtm1s2bKF2rVr88wzzxAUFMSwYcPYsWNHQcQpkllkJDz6qJajERGRIiFPg9IbN27M888/zzPPPMOFCxf44IMPuOOOO7jrrrvYvXt3fsUocmOJibBvn6OjEBGREuyWEiqr1crnn3/OfffdR5UqVVixYgVvvPEGf/75J4cPHyYkJIQePXrkd6wiWf38M9SqZbZWaRoFERFxkFwnVM8++yxBQUE8/fTT1KpVi+3bt7N+/XqefPJJSpUqRUhICNOnT+e3334riHhFMqtTB9LSYMcOmD/f0dGIiEgJleuEas+ePbz++uscP36c1157jfr162cpU7lyZdasWZMvAYrcUEAATJ5sPn/hBTh3zrHxiIhIiZTrhGrVqlU89thjeHh4XLeMm5sb7dq1y1NgIjk2ZAiEhcHp0zBliqOjERGREsgttwd8/fXX2W63WCx4eXkRGhpKtWrV8hyYSI65u8Nrr5l3/r3xBvzjH1C3rqOjEhGREiTXCVX37t2xWCwYhpFpe8Y2i8VCmzZt+PLLLylbtmy+BSpyQxER8MAD8PXX8Nxz5gLKmlJBREQKSa67/GJiYmjWrBkxMTEkJiaSmJhITEwMzZs355tvvuGnn37izJkzjBo1qiDiFbm+2bPB0xMCA82B6iIiIoUk1y1Uw4cP57333qNVq1b2bR07dsTLy4t//OMf7N69m9dee42BAwfma6AiNxUaCocOQeXKjo5ERERKmFy3UB08eBA/P78s2/38/Dh06BAANWvW5PTp03mPTiS3lEyJiIgD5DqhuuOOOxg9ejSnTp2ybzt16hRjxoyhWbNmAOzfv5/g4OD8i1Iktw4dgscfh4QER0ciIiIlQK67/N5//326d+9OcHAwISEhWCwW4uPjqV69Ol999RUAKSkpjB8/Pt+DFcmxgQPhxx/NOwAXLnR0NCIiUszlOqGqU6cOe/bsYeXKlezbtw/DMKhTpw4RERG4uJgNXt27d8/vOEVyZ8YMaNECFi0yp1Fo3drREYmISDGWq4TqypUreHl5ERsbS+fOnencuXNBxSWSN82bw5NPwvvvwzPPwJYt4Jbrvx9ERERyJFdjqNzc3KhSpQrpWoRWnMG0aVC2LOzYge2tt1h/8AxfxR5j/cEzpNuMmx8vIiKSQ7n+k33cuHGMHTuWjz76iHLlyhVETCL5IyAApk6FIUO4OGYsz+4vx+lS5mSzQf5eTOwaRuf6QQ4OUkREioNc3+U3b9481q5dS+XKlalduzZNmjTJ9BApSqJbdmFnYCi+ly/y1MYl9u0nElMZ8tE2onfpLkAREcm7W1p6RsQZpNsMJi+Po2LEENof2so7LR627zMACzB52R4iwgJxddEyNSIicutynVBNnDixIOIQyXebDp8lITGVhMq12VG5dpb9BpCQmMqmw2dpWaN84QcoIiLFRq67/ADOnz/P+++/z9ixYzl79iwA27Zt49ixY/kanEhenExOzbLNxZZOnZOHb1pOREQkN3LdQrVz507uuece/P39OXLkCIMHD6ZcuXJ88cUX/P777yxatKgg4hTJtYqlvTK9rpByjqjPJhJy/gQdB7/LKd+y2ZYTERHJrVy3UI0cOZL+/fuzf/9+vLyu/iK69957+emnn/I1OJG8aF6tHEH+XmSMjjrj44fV1RW/tIs8/8MHWDDv9mteTXeriohI3uQ6odq8eTNPPfVUlu233XYbJ06cyJegRPKDq4uFiV3DAHMAus3FlQkRQ7Bh4eHda2h+dBcTu4ZpQLqIiORZrhMqLy8vkpKSsmyPi4ujQoUK+RKUSH7pXD+It/s0IdDfbE3dGVSLxY07AfD+pig61w5wZHgiIlJM5Dqh6tatG1OmTMFqtQLYF0d+/vnnefjhh29ytEjh61w/iJ//dTeLB9/J3F6NqfWfeRjly1P6wG/w6quODk9ERIqBXCdUs2bN4tSpU1SsWJFLly7Rrl07QkNDKV26NC+//HJBxCiSZ64uFlrWKE+3xrfRrGlNLLNmmTsmTYIjRxwZmoiIFAO5vsvPz8+Pn3/+mdWrV7Nt2zZsNhtNmjThnnvuKYj4RApGv34QFQVWK6SlOToaERFxcrlOqDLcfffd3H333fkZi0jhsVhg6VIoUwZcbmk6NhEREbtbSqhWrVrFqlWrOHnyJDabLdO+Dz74IF8CEylwf1/c2zDMREtERCSXcp1QTZ48mSlTptC0aVOCgoKw6BeQOLuUFHMslc0Gc+Y4OhoREXFCuU6o3nnnHaKioujbt29BxCNS+DZvhtmzzdapxx6DZs0cHZGIiDiZXA8eSUtLo1WrVgURi4hjdOgAvXubXX5PPQVXrjg6IhERcTK5TqiefPJJ/vvf/xZELCKOM3u2OUB9+3Z44w1HRyMiIk4m111+qampvPfee3z//fc0bNgQd3f3TPvnaAyKOKNKlWDGDPjHP2DcOHj4YQgJcXRUIiLiJHKdUO3cuZPGjRsDsGvXrkz7NEBdnNqgQbBwIfzyCwwbBl984eiIRETESeQ6oVqzZk1BxCHieC4u8M47EB4Oq1dDfDzcfrujoxIRESdwyzMaHjhwgBUrVnDp0iUADMPIt6BEHKZ+fbOVas8eJVMiIpJjuU6ozpw5Q8eOHalVqxb33XcfCQkJgDlY/f/+7//yPUCRQvf443DbbY6OQkREnEiuE6rnnnsOd3d34uPj8fHxsW/v2bMn0dHR+RqciMNFR8PatY6OQkREirhcj6FauXIlK1asIDg4ONP2mjVr8vvvv+dbYCIOt3Ah9O8PNWrAzp1wzR8QIiIi18p1C9WFCxcytUxlOH36NJ6envkSlEiR8OCDZtffwYPm0jQiIiLXkeuEqm3btixatMj+2mKxYLPZmDlzJh06dMjX4EQcys8P3n7bfD57NmzZ4th4RESkyMp1l9/MmTNp3749W7ZsIS0tjTFjxrB7927Onj3LL7/8UhAxijhO167Qqxd88ok5T9WWLfC3yWxFRERy3UIVFhbGzp07ad68OREREVy4cIGHHnqI7du3U6NGjYKIUcSx5s6FcuXMcVQzZzo6GhERKYJy3UIFEBgYyOTJk/M7FpGiqWJFeO01eOIJmDIF+vXTtAoiIpLJLSVUIiVOnz6wZg088ICSKRERyUIJlUhOWCzwwQeOjkJERIqoW156prC9/PLLtGrVCh8fH8qUKXPdclFRUTRs2BAvLy8CAwMZOnToDc/bvn17LBZLpkevXr0yldm3bx/dunUjICAAPz8/WrdurTUNS7oTJ8y1/kRERHCihCotLY0ePXowZMiQ65aZM2cOL7zwAs8//zy7d+9m1apVdOrU6abnHjx4MAkJCfbHu+++m2n//fffz5UrV1i9ejVbt26lcePGdOnShRMnTuT5usQJff891Ktnjqmy2RwdjYiIFAFO0+WXMQg+Kioq2/3nzp1j3LhxLFu2jI4dO9q316tX76bn9vHxITAwMNt9p0+f5sCBA3zwwQc0bNgQgOnTp/PWW2+xe/fu6x4nxVj16nD5Mvz4I7z5Jjz7rKMjEhERB8tRQhUeHo7FYsnRCbdt25angG5VTEwMNpuNY8eOUbduXZKTk2nVqhWzZ88mJCTkhsd+/PHHfPTRR1SqVIl7772XiRMnUrp0aQDKly9P3bp1WbRoEU2aNMHT05N3332XSpUqcccdd1z3nJcvX+by5cv210lJSQBYrVasVms+XHHxlfH5FNnPKSQEl2nTcB02DOP557lyzz0QGuroqBymyNeXZKE6cy6qL8fK6eeeo4Sqe/fueYmlUBw6dAibzcbUqVOZO3cu/v7+jBs3joiICHbu3ImHh0e2x/Xu3Ztq1aoRGBjIrl27GDt2LDt27CAmJgYwZ4KPiYmhW7dulC5dGhcXFypVqkR0dPQNx3JNmzYt26klVq5cme3SPZJVRh0UScHBtGrQgAq//kriI4/wy0svgYvT9KAXiCJdX5It1ZlzUX05xsWLF3NUzmIYhlHAsVzXpEmTbjqf1ebNm2natKn9dVRUFCNGjOD8+fOZyk2dOpUXXniBFStWEBkZCcCpU6cIDAxk+fLlORpLBbB161aaNm3K1q1badKkCYZh0L17d6xWKy+88ALe3t68//77fP3112zevJmgoKBsz5NdC1VISAinT5/Gz88vR7GUVFarlZiYGCIiInAvyrOSHzmCW5MmWFJSSJ89G1sJ7fpzmvoSO9WZc1F9OVZSUhIBAQEkJibe8Pe3Q8dQDR06NMsddX9XtWrVHJ0rI7EJCwuzb6tQoQIBAQHE5+JurCZNmuDu7s7+/ftp0qQJq1ev5ptvvuHcuXP2D/Ktt94iJiaGhQsX8vzzz2d7Hk9Pz2wXi3Z3d9cXIoeK/GdVs6Y5c/qQIbiOG4dr167mthKqyNeXZKE6cy6qL8fI6Wee64QqPT2dV199lU8//ZT4+HjS0tIy7T979myOzxUQEEBAQEBuQ8hW69atAYiLiyM4ONgey+nTp6lSpUqOz7N7926sVqs9Qcto6nP5W3eOi4sLNt3hJU89BUuXmonUdVorRUSk+Mv1oI/JkyczZ84cHn30URITExk5ciQPPfQQLi4uTJo0qQBCNMXHxxMbG0t8fDzp6enExsYSGxtLSkoKALVq1aJbt24MHz6cdevWsWvXLvr160edOnXo0KEDAMeOHaNOnTps2rQJgIMHDzJlyhS2bNnCkSNHWL58OT169CA8PNyeoLVs2ZKyZcvSr18/duzYwb59+xg9ejSHDx/m/vvvL7DrFSdhscC335p3+/n6OjoaERFxkFwnVB9//DH/+c9/GDVqFG5ubjz22GO8//77TJgwgQ0bNhREjABMmDCB8PBwJk6cSEpKCuHh4YSHh7NlyxZ7mUWLFtGiRQvuv/9+2rVrh7u7O9HR0fbmOqvVSlxcnL3VycPDwz5XVe3atRk2bBiRkZF8//33uLq6AmYrWnR0NCkpKdx99900bdqUn3/+ma+++opGjRoV2PWKE7m2Odhmg1y00oqISPGQ6y6/EydO0KBBAwB8fX1JTEwEoEuXLowfPz5/o7tGVFTUdeegyuDn58f8+fOZP39+tvurVq3KtWPwQ0JC+PHHH2/63k2bNmXFihW5ildKoPh46N8fLlyAX34BN6eZ5k1ERPIo1y1UwcHBJCQkABAaGsrKlSsB82687AZhi5QYFgts2wabNsHUqY6ORkREClGuE6oHH3yQVatWATB8+HDGjx9PzZo1eeKJJxg4cGC+ByjiNEJCzLFUAFOmwObNjo1HREQKTa77JKZPn25//sgjjxAcHMy6desIDQ3lgQceyNfgRJzO44/D11/Dp59Cnz6wfTtoIlcRkWIvz4M87rzzTu688878iEXE+Vks8Pbb8PPPsG8fjBkDb7zh6KhERKSA5Sih+vrrr7n33ntxd3fn66+/vmFZtVJJiVeuHERFQWSk2QXYpQt07uzoqEREpADleC2/EydOULFixRuu62exWEhPT8+v2EScV0QEPPssrFoFgYGOjkZERApYjhKqa2cE1+zgIjn0yitmF6CXl6MjERGRApbru/wWLVqUadHfDGlpaSxatChfghIpFry9MydTf83qLyIixU+uE6oBAwbYJ/O8VnJyMgMGDMiXoESKlfR0mDYNqlc3J/8UEZFiJ9cJlWEYWCyWLNv/+OMP/P398yUokWLFMMypFE6dgieeMBMsEREpVnI8bUJ4eDgWiwWLxULHjh1xu2ZZjfT0dA4fPkxn3ckkkpWbG3z4ITRuDD/+aI6t+ve/HR2ViIjkoxwnVBl398XGxtKpUyd8fX3t+zw8PKhatSoPP/xwvgcoUiyEhprzUQ0YABMmwN13g+ZvExEpNnKcUE2cOBEwFxju2bMnXrpzSSR3+vWDFSvgk0/MGdW3bwd1k4uIFAu5nim9X79+gHlX38mTJ7NMo3D77bfnT2QixY3FAu+8Axs2wOHDMGQIfPyxuV1ERJxarhOq/fv3M3DgQNatW5dpe8ZgdU3sKXID/v6weDHcc4+6/EREipFcJ1T9+/fHzc2Nb775hqCgoGzv+BORG7jzTvj9dyhf3tGRiIhIPsl1QhUbG8vWrVupU6dOQcQjUjJcm0wlJ4OnJ3h4OC4eERHJk1zPQxUWFsbp06cLIhaRkmfLFmjSBMaNc3QkIiKSB7lOqF555RXGjBnDDz/8wJkzZ0hKSsr0EJFc+OMPOHAAZs6ElSsdHY2IiNyiXHf53XPPPQB07Ngx03YNShe5Bd27w9NPm3f/PfEE7NgBlSo5OioREcmlXCdUa9asKYg4REquOXNg7VrYvducn2rlSnB1dXRUIiKSC7lOqNq1a1cQcYiUXN7e8Nln0KwZrF4NL74IkyY5OioREcmFXI+hAli7di19+vShVatWHDt2DIAPP/yQn3/+OV+DEykx6tY1u/0ApkwxEysREXEauU6olixZQqdOnfD29mbbtm1cvnwZgOTkZKZOnZrvAYqUGH36wJNPQmQkNGjg6GhERCQXcp1QvfTSS7zzzjv85z//wd3d3b69VatWbNu2LV+DEylx3ngDli+HChUcHYmIiORCrhOquLg42rZtm2W7n58f58+fz4+YREouT09wueZruWOH42IREZEcy3VCFRQUxIEDB7Js//nnn6levXq+BCVS4qWnw6BBEB6u+alERJxArhOqp556iuHDh7Nx40YsFgvHjx/n448/ZtSoUfzzn/8siBhFSh5XV3B3B8OA3r3hr5s/RESkaMr1tAljxowhMTGRDh06kJqaStu2bfH09GTUqFEMHTq0IGIUKZleew02boTYWHjsMfPOP7dcf2VFRKQQ3NK0CS+//DKnT59m06ZNbNiwgVOnTvHiiy/md2wiJZuXlzk/VenS5sSfY8c6OiIREbmOXCdUAwcOJDk5GR8fH5o2bUrz5s3x9fXlwoULDBw4sCBiFCm5QkNhwQLz+axZZoIlIiJFTq4TqoULF3Lp0qUs2y9dusSiRYvyJSgRucbDD8Po0ebzJ5+Ec+ccG4+IiGSR4wEZSUlJGIaBYRgkJyfj5eVl35eens7y5cupWLFigQQpUuJNnQqHDpkJVdmyjo5GRET+JscJVZkyZbBYLFgsFmrVqpVlv8ViYfLkyfkanIj8xc0NPv/c0VGIiMh15DihWrNmDYZhcPfdd7NkyRLKlStn3+fh4UGVKlWoXLlygQQpIn9z6JA5UL1fP0dHIiIi5CKhateuHQCHDx8mJCQEF5dbukFQRPLqjz+gaVNITITgYOjY0dERiYiUeLme1KZKlSqcP3+e+fPns3fvXiwWC2FhYQwcOBB/f/+CiFFErnXbbfDgg/DBB9CrF2zdCrff7uioRERKtFw3M23ZsoUaNWrw6quvcvbsWU6fPs2cOXOoUaOGFkcWKQwWC7z5JtxxB5w+bd4FmM2dtyIiUnhynVA999xzPPDAAxw5coSlS5fyxRdfcPjwYbp06cKIESMKIEQRycLLC5YsgfLlYcsW+Mc/zGVqRETEIW6phepf//oXbtcsgeHm5saYMWPYsmVLvgYnIjdQpQp8+qm57t9HH8Hs2Y6OSESkxMp1QuXn50d8fHyW7UePHqV06dL5EpSI5NDdd5tr/gF8+SVcueLIaERESqxcD0rv2bMngwYNYtasWbRq1QqLxcLPP//M6NGjeeyxxwoiRhG5kWeeMdf769lTiyeLiDhIrn/6zpo1C4vFwhNPPMGVv/4adnd3Z8iQIUyfPj3fAxSRm7BYss5HdeWKkisRkUKU65+4Hh4ezJ07l2nTpnHw4EEMwyA0NBQfH5+CiE9EcsNmg3HjYPt2+OYbc3yViIgUuByPobp48SLPPPMMt912GxUrVuTJJ58kKCiIhg0bKpkSKSoOHDDHVEVHw/PPOzoaEZESI8cJ1cSJE4mKiuL++++nV69exMTEMGTIkIKMTURyq1YtWLDAfD5rFnz4oWPjEREpIXLc5bd06VLmz59Pr169AOjTpw+tW7cmPT0dV3UriBQdPXvCzp0wdSo8+SRUqwZt2jg6KhGRYi3HLVRHjx7lrrvusr9u3rw5bm5uHD9+vEACE5E8ePFFc3matDTo3h0OHnR0RCIixVqOE6r09HQ8PDwybXNzc7Pf6SciRYiLi9ndd8cdcOYMdOkCly87OioRkWIrx11+hmHQv39/PD097dtSU1N5+umnKVWqlH3b0qVL8zdCEbk1pUrBsmXQrh1MmADXfHdFRCR/5Tih6vf3eW4wx1GJSBEWFAS7d4O7u6MjEREp1nKcUC3IuHNIRJzLtclUQgLExMATTzguHhGRYkhTKYuUFGfOQPPm8Mcf4OMDjzzi6IhERIqNXC+OLCJOqnx5884/gL59YdMmx8YjIlKMKKESKUlefRXuvx9SU807/w4ccHREIiLFghIqkZLE1RUWL4YmTeDUKejcGU6edHRUIiJOTwmVSElTujR8+605g/rBg2aLVUqKo6MSEXFqSqhESqLAQHMB5fLlze6/5GRHRyQi4tR0l59ISVWrFqxaBVWqQJkyjo5GRMSpOU0L1csvv0yrVq3w8fGhzA1++EdFRdGwYUO8vLwIDAxk6NChNz33+vXrufvuuylVqhRlypShffv2XLp0yb7/3Llz9O3bF39/f/z9/enbty/nz5/Ph6sScbBGjTInU7t3OywUERFn5jQJVVpaGj169GDIkCHXLTNnzhxeeOEFnn/+eXbv3s2qVavo1KnTDc+7fv16OnfuTGRkJJs2bWLz5s0MHToUF5erH83jjz9ObGws0dHRREdHExsbS9++ffPt2kSKhLlzoUEDePddR0ciIuJ0nKbLb/LkyYDZApWdc+fOMW7cOJYtW0bHjh3t2+vVq3fD8z733HMMGzaM559/3r6tZs2a9ud79+4lOjqaDRs20KJFCwD+85//0LJlS+Li4qhdu/atXpJI0XL2LBgG/POfpAdUYFPjdpxMTqViaS+aVyuHq4vF0RGKiBRZTpNQ3UxMTAw2m41jx45Rt25dkpOTadWqFbNnzyYkJCTbY06ePMnGjRvp3bs3rVq14uDBg9SpU4eXX36ZNm3aAGYLlr+/vz2ZArjzzjvx9/dn3bp1102oLl++zOXLl+2vk5KSALBarVit1vy67GIp4/PR51TIXngB16NHcVmwgPSevXjn0YlsqNYIgEA/L56/tw731K2U5TDVl/NRnTkX1Zdj5fRzLzYJ1aFDh7DZbEydOpW5c+fi7+/PuHHjiIiIYOfOnXh4eGR7DMCkSZOYNWsWjRs3ZtGiRXTs2JFdu3ZRs2ZNTpw4QcWKFbMcW7FiRU6cOHHdeKZNm2ZvVbvWypUr8fHxycOVlhwxMTGODqHEsXTpQtO9e6m8YQPzv3iJdVOmcK5WLeACaYe3svzw9Y9VfTkf1ZlzUX05xsWLF3NUzqEJ1aRJk7JNOq61efNmmjZtetNz2Ww2rFYr8+bNIzIyEoDFixcTGBjImjVrsh1LZbPZAHjqqacYMGAAAOHh4axatYoPPviAadOmAWCxZO3qMAwj2+0Zxo4dy8iRI+2vk5KSCAkJITIyEj8/v5teT0lmtVqJiYkhIiIC92sX9pUClW4z6PTaT5xpO4b3Tkym1ZEdNJj4In17T2N/xSpYgEp+XqwY0TZT95/qy/mozpyL6suxMnqYbsahCdXQoUPp1avXDctUrVo1R+cKCgoCICwszL6tQoUKBAQEEB8fn+NjAOrWrWs/JjAwkD///DPLsadOnaJSpazdHxk8PT3x9PTMst3d3V1fiBzSZ1W4thw8w+/nLoPFgycfHMfHn4wjPCGONvu3sKt8VQB+P3eZ7X8k07JG+SzHq76cj+rMuai+HCOnn7lDE6qAgAACAgLy5VytW7cGIC4ujuDgYADOnj3L6dOnqVKlSrbHVK1alcqVKxMXF5dp+759+7j33nsBaNmyJYmJiWzatInmzZsDsHHjRhITE2nVqlW+xC5SFJxMTrU/v+jhTf8ek+i0bz2fNoq8bjkRETE5zbQJ8fHxxMbGEh8fT3p6OrGxscTGxpLy15IZtWrVolu3bgwfPpx169axa9cu+vXrR506dejQoQMAx44do06dOmzatAkwu/JGjx7NvHnz+Pzzzzlw4ADjx4/nt99+Y9CgQYDZWtW5c2cGDx7Mhg0b2LBhA4MHD6ZLly66w0+KlYqlvTK9TvQunSmZ8rRexi81JUs5ERFxokHpEyZMYOHChfbX4eHhAKxZs4b27dsDsGjRIp577jnuv/9+XFxcaNeuHdHR0fbmOqvVSlxcXKYBZiNGjCA1NZXnnnuOs2fP0qhRI2JiYqhRo4a9zMcff8ywYcPsY7MeeOAB3njjjYK+ZJFC1bxaOYL8vTiRmIrxt32+ly8yf8kUSnGFupPvdUh8IiJFmcUwjL//7JQCkJSUhL+/P4mJiRqUfhNWq5Xly5dz3333abxAIYvelcCQj7YBZEqqapw5yucfjaFsajLcdRd89x2UKgWovpyR6sy5qL4cK6e/v52my09ECl7n+kG83acJgf6Zu/UuVq/J3qjPwN8f1q6Frl0hh7cSi4iUBE7T5ScihaNz/SAiwgLZdPhs1pnSq66AiAhYswa6d4evvwZXV0eHLCLicGqhEpEsXF0stKxRnm6Nb6NljfJX551q0QKWLze7+2Ji4OGH4ZoVAURESiolVCKSO23awDffgLc3bNgAR444OiIREYdTl5+I5F779rBsGVSoALVrw8GDjo5IRMSh1EIlIremY0do2PDq619/hStXHBePiIgDKaESkTwrt3cvbm3bQq9ekJbm6HBERAqdEioRyTP3lBSwWmHJEnjkEQ1UF5ESRwmViOTZn82akb5kCXh5mWOruneHS5ccHZaISKFRQiUi+cLo1Onq3X/R0fDAA5r8U0RKDCVUIpJ/OnY0kylfX/j+e7j3XvhrAXMRkeJMCZWI5K+2bWHFCvDzM5eq8fBwdEQiIgVO81CJSP5r1QrWrYPQUCVUIlIiqIVKRApGvXrg6Wk+NwyYPBmOHnVsTCIiBUQJlYgUvFdegUmToHVriItzdDQiIvlOCZWIFLzevaFOHbOFqk0b2LrV0RGJiOQrJVQiUvBCQuCnn+COO+D0aejQAX74wdFRiYjkGyVUIlI4KlSA1avNhZWTk6FzZ/jqK0dHJSJOLt1msP7gGb6KPcb6g2dItxkOiUN3+YlI4fHzg+++M9f8++or89/DhyEw0NGRiYgTit6VwORle0hITLVvC/L3YmLXMDrXDyrUWNRCJSKFy8sLPv8cBg2C995TMiUityR6VwJDPtqWKZkCOJGYypCPthG9K6FQ41FCJSKFz80N3n8f+va9ui0+HtLSHBeTiDiNdJvB5GV7sHfuGVe7+TKeTV62p1C7/5RQiYjjnTgB7drBffdBYqKjoxGRIm7T4bP2lqmmf+xmedQwApNO2/cbQEJiKpsOny20mJRQiYjj7dsHp07BqlVw113wxx+OjkhEirCTyalgGAza/CWf/HcsYScPM2rth9mXKyRKqETE8dq2NadVCAyEX3+Fli1hxw5HRyUiRVTF0l5EHNjI+NXv42bY+KpuOyZEPJ1tucKihEpEioYmTWD9eqhb12yhat0avv7a0VGJSBHUvFo5dt3Rjq/qtmN8xNMM7zqKix7e9v0WzLv9mlcrV2gxKaESkaKjalX45Rfo2BEuXIDu3eHTTx0dlYgUBYYBixdDUhKuLhYmPlCPEV1H8VGTLmCx2ItlPJvYNQxXF0v25yoASqhEpGgpW9acq2rIEKhVCyIjHR2RiDhacjL06QOPPw5PPw2GQef6Qbzd9w4C/TN36wX6e/F2nyaFPg+VJvYUkaLH3R3efNO8469MGXObYUBKCpQu7dDQRKSQbdsGPXvCgQPg6goNGth3da4fRERYIJsOn+VkcioVS5vdfIXZMpVBLVQiUjRZLFeTKTATrIYNzUHrIlL8GQbMm2fepHLggLkm6I8/wtixmbr4XF0stKxRnm6Nb6NljfIOSaZACZWIOIPLl+GNN+DIEbjzTo2rEinuzp2Dhx6C4cPNCX+7dYPYWPNmlSJKCZWIFH2enlcHq1+8aDb/jx4NV644OjIRKQhXrsDGjeDhYbZSffEFlCu8O/ZuhRIqEXEO5ctDdDSMGWO+njULOnUyJwQVEedntV5dQqZCBbMlev16ePbZTF18RZUSKhFxHm5u8Mor5g/aUqVg9WqzC/DSJUdHJiJ5sWcPNG8OH310dVubNub8dE5CCZWIOJ8ePWDDBggNhWeeAW/vmx8jIkWPYcDrr8Mdd5hjpCZONFuqnJCmTRAR51S/PmzfbrZUZTh8GCpWzLxNRIqm48dh4EBYscJ83akTLFhgTpvihNRCJSLOy9f36tiKlBS47z6z22D3bsfGJSI39tln5jQoK1aAl5fZSvXddxBUuJNx5iclVCJSPBw5Yk4EumcPNGtm/qWbMcBVRIqOvXvNO3XPnIHwcHPizqFDnWLg+Y0ooRKR4qF+fXMMRmSkOUh94EDo189suRKRoqNuXRg1CsaPN8dC1q3r6IjyhRIqESk+KlY0uw2mTgUXF/jwQ7O1avt2R0cmUnL9+ae5Bt/+/Ve3zZgBU6aY80wVE0qoRKR4cXExl6b44Qe47Tb47TeYMMHRUYmUPIYBixdDWJj579NPOzqiAqWESkSKp7vuMrsA+/eHd991dDQiJcvRo9C9u9kydfYsNG4Mc+Y4OqoCpYRKRIqvgABzcHrlyle3jRoF//2v42ISKc7S02HuXLNV6uuvzSkQpkyBTZugUSNHR1egNA+ViJQcq1fD7Nnm82++gTffhLJlHRuTSHGyYAGMGGE+b9UK3nsP6tVzaEiFRS1UIlJytG0LkyeDq6s5pqNePVi2zNFRiRQf/fqZ3e3vvgtr15aYZAqUUIlISeLmZg5Q/+UXqFULEhLggQegTx9zThwRyTnDgC+/hM6dIS3N3ObuDj/+CP/4h3mDSAlSsq5WRASgRQtzwPro0eYP/Y8/ho4dNRGoSE7FxZkrEzz4oDnb+VtvXd3n5BN03iolVCJSMnl7m3PhrFtnDqCdNKnE/iIQybHkZBgzBho0gOhocx6psWPNFqkSToPSRaRky2itunZB1k8+gQsXYMCAEtdtIZItwzBbcseMMbvKAe6/H159FWrWdGxsRYR+UoiIXJtMnTljriv25JPm4Npff3VcXCJFyYcfmslUaKh5l+w33yiZuoYSKhGRa/n7w7//DaVKmd2B4eHmWKvrrAmYbjNYf/AMX8UeY/3BM6TbNA5Lion4eDh1ynxusZjzS02bBrt2ma1TkokSKhGRa7m5wciRsHcvPPSQOVHhrFnmOKsvv8w0cD16VwJtXlnNY//ZwPBPYnnsPxto88pqonclOC5+kbxKTITnnzfvhL122aY6dcztnp6Oi60IU0IlIpKdkBBYssScp6pKFXMpjYcegn37ADOZGvLRNhISUzMddiIxlSEfbVNSJc4nLQ3mzYMaNeCVV+DyZThwwPyjQm5KCZWIyI106QJ79ph3Mj39NNSuTbrNYPKyPbhfsWYpntF+NXnZHnX/iXOw2eDTT81JOIcPN8cR1q1r/jGxcqU5Ea7clBIqEZGb8fGBqVPNpWqATYfP4nX4IL+8M4D+W77GLf1KpuIGkJCYyqbDZx0QrEguzZkDPXuarVGVKsE778DOneYfE5pKJMeUUImI5NRfv1xOJqfSd9u3VLhwnkmr3iP6g6Hcs39jlolBTyanZncWEccyDHOcVIYBAyA42FyWaf9+eOopcyyh5Io+MRGRXKpY2ovn7h7E/oDb+b+1HxJ69g/eX/oiWyvX4ZX2/dkUUt9eTqRI+eEHGD/e7MZbs8b8I6F8eTh0KPP0IZJraqESEcml5tXKUalsKT5p3JkO/3iPt+58hEtuntxx/Dc+/e/zzP16JkH+XjSvVs7RoYqYLVIrVpiLg3foAD//DBs2wMGDV8somcozJVQiIrnk6mJhYtcwAFI8SzGjXX/aPvUfPgy/D6uLK3EVqjCxaxiuLhp/Ig5ks8HSpdCsmbmA8dq1ZuI0ZIiZTIWGOjrCYsVpEqqXX36ZVq1a4ePjQ5kyZa5bLioqioYNG+Ll5UVgYCBDhw696bnXr1/P3XffTalSpShTpgzt27fn0qVLABw5coRBgwZRrVo1vL29qVGjBhMnTiQtY2VtESmROtcP4u0+TQj0N7v1TvmWY3zkP3nsuQWEz5xA5/pBZsGvv4bevc07BUUK05dfwsMPw9at5o0Vzz0Hhw+bCxnfdpujoyt2nGYMVVpaGj169KBly5bMnz8/2zJz5sxh9uzZzJw5kxYtWpCamsqhQ4dueN7169fTuXNnxo4dy+uvv46Hhwc7duzA5a/1u3777TdsNhvvvvsuoaGh7Nq1i8GDB3PhwgVmzZqV79cpIs6jc/0gIsIC2XT4LCeTU6lY2uzms7dMGQa8+CJs2QKLF5vzWI0bB40bOzRuKfrSbYb9/1WATw5/VSclQVyc2SIF8MAD0KQJ3HcfDBsGFSoUXMCCxTAMp5ooJSoqihEjRnD+/PlM28+dO8dtt93GsmXL6NixY47Pd+eddxIREcGLL76Y42NmzpzJ22+/fdNk7VpJSUn4+/uTmJiIn59fjo8riaxWK8uXL+e+++7DXf36RZ7q6ya2b4eXXjK7XjJ07WomVs2bOyQk1VnRFr0rgcnL9tgnjfV0NZjRPB2Pandwb8PgrAccOmROyPnBB+DnZ7ZCZdSrYWjqgzzK6e9vp+nyu5mYmBhsNhvHjh2jbt26BAcH8+ijj3L06NHrHnPy5Ek2btxIxYoVadWqFZUqVaJdu3b8/PPPN3yvxMREypXTYFMRyYHwcHPG9V9/hV69zF9uy5ZBixZmq4HINa43Az/Ac/+LvToDv2GYY6IeeshcoHjuXEhONhOq+PirBymZKjRO0+V3M4cOHcJmszF16lTmzp2Lv78/48aNIyIigp07d+Lh4ZHtMQCTJk1i1qxZNG7cmEWLFtGxY0d27dpFzWxW0T548CCvv/46s2fPvmE8ly9f5vLly/bXSUlJgPmXodWadXZluSrj89Hn5BxUXzlUuzYsWgQvvIDrjBlYFi8mvXVrjIzPLSXF/OVXqlSBh6I6K5rSbQbTvt2Nh2vmjiNPF8P+77Rvd9Phz714jB6FZft2exlbZCS2YcMwIiLM/0eq23yT0++JQxOqSZMmMXny5BuW2bx5M02bNr3puWw2G1arlXnz5hEZGQnA4sWLCQwMZM2aNXTq1CnbYwCeeuopBgwYAEB4eDirVq3igw8+YNq0aZnKHz9+nM6dO9OjRw+efPLJG8Yzbdq0bK9t5cqV+Pj43PR6xGx1FOeh+sqFhx/Gq317Ut3dYflyAGp+9hmhX3/NkchIjkRGcqlSpQIPQ3VW9Iysk/12F6uVKU3dgQus33mI9tu3k+7hwdF27Tj0wAMkh4TAlSvw3XeFGm9JcPHixRyVc2hCNXToUHr16nXDMlWrVs3RuYKCzDtqwsLC7NsqVKhAQEAA8dc2f97kGIC6detmOeb48eN06NCBli1b8t577900nrFjxzJy5Ej766SkJEJCQoiMjNQYqpuwWq3ExMQQERGh8R1OQPWVP1xnzsQlOZlaS5ZQc+lSjHvvxTZ4MEbnzvm+lprqrGha/msCY5bstL/2SbtE572/8Fjsd5SpHUyXViO4bLMw4+GHuTLfH+O++7itfHl0v17ByuhhuhmHJlQBAQEEBATky7lat24NQFxcHMHB5qC9s2fPcvr0aapUqZLtMVWrVqVy5crExcVl2r5v3z7uvfde++tjx47RoUMH7rjjDhYsWGC/A/BGPD098fT0zLLd3d1dP8BySJ+Vc1F95dGPP5pTLLz1Fpbvv8eyfDkuy5fD7bebY63+7//y/S1VZ0VLRf9SWK02Wsb/ykO7VnHvvnX4WM2hI9Zzf8AdaVx28aKifyncBg50cLQlR06/I04zhio+Pp6zZ88SHx9Peno6sbGxAISGhuLr60utWrXo1q0bw4cP57333sPPz4+xY8dSp04dOnToAJiJUceOHVm0aBHNmzfHYrEwevRoJk6cSKNGjWjcuDELFy7kt99+4/PPPwfMlqn27dtz++23M2vWLE6dOmWPKTAwsNA/BxEpplxd4cEHzcf+/fDuu7BggTnAeNeuq+UMA1JTwdvbcbFKgWjx7cesf+8lKiVe/T1zqGxllja6h/r97yZtnydBvpqBv6hymoRqwoQJLFy40P46PDwcgDVr1tC+fXsAFi1axHPPPcf999+Pi4sL7dq1Izo62p5dWq1W4uLiMvWHjhgxgtTUVJ577jnOnj1Lo0aNiImJoUaNGoA55unAgQMcOHDA3vKVwclmnBARZ1GzJsyaZU638Nln0KjR1X1btsA998Cjj0K/ftC6te7kclb79pkTbP51I4JLSgqVEk9x3suXZXXbsrTe3WyvXBtPN5hRJh1AM/AXYU43D5Wz0jxUOac5cpyL6quQvfACTJ169XWVKvDII+ajRYscJVeqMwfat89Mkj/7DHbsgP/+Fx57zNx37Bhs2MCKancwacXBnM9DJQUqp7+/naaFSkREMGdej4gwp2D47DP4/XeYPdt8hITATz9BDm/mkUKQnm4uRPzNN+bj2u5bV1ezezfDbbfBww/TCbincZVMM6Wf3ruBe+oW/J2fcuuUUImIOBMXF2jf3ny8+SZER5uJ1bJl5txDISFXy37wAZQuDZGR4O/vqIhLnvT0q3dm/vkntGlzdZ+bG3TsCD16QPfuUL58tqdwdbHQsoa5z2q1snxvAccseaaESkTEWXl7Xx3InppqtnZk/CK32eDf/zZ/obu5mb/U77vPbN2qW9excRc3Fy/CL7/A6tWwapWZJGXMB1W5svmZV6gAXbpAp06glTaKJSVUIiLFgZcXNGhw9fWlS9C7N3z7rblg7g8/mA/ArVw5wtq3NxMsuTXr15vJ06pVsG4dpKVd3efubs587+trvl650jExSqFSQiUiUhyVKnV1bNXBg2ZiFR0Na9diOXvWbMHKkJgITzwBd95pPpo2NbsKxZym4vBh8+7KHj2uDvqfPt2cNyzDbbeZXXkdO8Ldd19NpqTEUEIlIlLc1ahhTg46bBhYrVzZuJHfd+2iasb+TZvM5CAjQXBxgXr1oHlzc8qGTp2gVi0HBV+I0tJg716IjTXvwMv49+xZc3+LFuZdlQD33w+entChg5lE1ayp6StKOCVUIiIlibs7RosWXDhz5uq22rVhzhzzbrSNG807B3/91XwAvPfe1YRq715zwtGaNc1HaKg5TigHK0gUGWfPmtMXxMWZ45oyBoZPmACvvJK1vIcHNGliHpeRUP3jH+ZD5C9KqERESrrbb4fnnrv6OiHBTKy2bjWTqmbNru5buxZmzsx8vLe32QpWrRqMG2e2bAGcPm0Oiq9YEcqWJd3F1T4VQMXS5ozf+TpJpdVqJj1+fldnkl+3Dj7+GI4ehT/+MJPFjBYnMMc3RUSYzxs3Nu+GbNTIfJ7xCAszW6NEbkAJlYiIZBYUZN7S37171n1168Kzz5p3FB44YI4vunTJnF9p1y64ZlF4li6Fp56yv7zk6UOwV2lKe/ly3suXQV2eotfT3elcP8hMfD77zExcLBZz7FLGA2DAgKt3J8bEwNy55t11KSlw5oz5SEw093/zjdklB2acb72V9TqCg81Wt2snNn3kEejZU113ckuUUImISM7ddZf5yGC1mq0++/eb6w7Wq3d1X1oalC0L584B4Hv5Ir6XLxKS+CcAb545z5CPtvF2nyZ03rEDXnvtxu+bkVAlJJiD7K8nKenq86ZNzdnlQ0LMJCokxGxN+2u5l0zc9CtRbp3+94iIyK1zdzfHUYWGZt03dCjp/3yGdlNjuHDyNGVSU/BPTaHMpWT8LqewP8CchHTysj1EtGuM6/PPm/NpgdlKlNFSZLFA9epXz9uqFcyfDz4+ZmJUrpw5DiogwEzgMubiAjPBe+mlArp4kauUUImISIHZdPgsf6RYwcefcz7Zz9aekJjKpkqNaTltWs5Oer0ETsSBnOi2DBERcTYnk1PztZxIUaWESkRECkzF0l75Wk6kqFJCJSIiBaZ5tXIE+XtxvfvmLECQvzmFgogzU0IlIiIFxtXFwsSuYQBZkqqM1xO7huXvfFQiDqCESkREClTn+kG83acJgf6Zu/UC/b3MKRPqBzkoMpH8o7v8RESkwHWuH0REWGDBzpQu4kBKqEREpFC4ulhoWaO8o8MQKRDq8hMRERHJIyVUIiIiInmkhEpEREQkj5RQiYiIiOSREioRERGRPFJCJSIiIpJHSqhERERE8kgJlYiIiEgeKaESERERySPNlF5IDMMAICkpycGRFH1Wq5WLFy+SlJSEu7u7o8ORm1B9OR/VmXNRfTlWxu/tjN/j16OEqpAkJycDEBIS4uBIREREJLeSk5Px9/e/7n6LcbOUS/KFzWbj+PHjlC5dGotFi4HeSFJSEiEhIRw9ehQ/Pz9HhyM3ofpyPqoz56L6cizDMEhOTqZy5cq4uFx/pJRaqAqJi4sLwcHBjg7Dqfj5+emHhxNRfTkf1ZlzUX05zo1apjJoULqIiIhIHimhEhEREckjJVRS5Hh6ejJx4kQ8PT0dHYrkgOrL+ajOnIvqyzloULqIiIhIHqmFSkRERCSPlFCJiIiI5JESKhEREZE8UkIlBWratGk0a9aM0qVLU7FiRbp3705cXFymMoZhMGnSJCpXroy3tzft27dn9+7dNz33kiVLCAsLw9PTk7CwML744ouCuowSo6DqKyoqCovFkuWRmppakJdT7OWkvpYuXUqnTp0ICAjAYrEQGxubo3Pr+1UwCqrO9B1zPCVUUqB+/PFHnnnmGTZs2EBMTAxXrlwhMjKSCxcu2MvMmDGDOXPm8MYbb7B582YCAwOJiIiwL9eTnfXr19OzZ0/69u3Ljh076Nu3L48++igbN24sjMsqtgqqvsCclDAhISHTw8vLq6AvqVjLSX1duHCB1q1bM3369ByfV9+vglNQdQb6jjmcIVKITp48aQDGjz/+aBiGYdhsNiMwMNCYPn26vUxqaqrh7+9vvPPOO9c9z6OPPmp07tw507ZOnToZvXr1KpjAS6j8qq8FCxYY/v7+BR1uiff3+rrW4cOHDcDYvn37Tc+j71fhya8603fM8dRCJYUqMTERgHLlygFw+PBhTpw4QWRkpL2Mp6cn7dq1Y926ddc9z/r16zMdA9CpU6cbHiO5l1/1BZCSkkKVKlUIDg6mS5cubN++veACL6H+Xl+3St+vwpNfdQb6jjmaEiopNIZhMHLkSNq0aUP9+vUBOHHiBACVKlXKVLZSpUr2fdk5ceJEro+R3MnP+qpTpw5RUVF8/fXXLF68GC8vL1q3bs3+/fsL7gJKmOzq61bp+1U48rPO9B1zPC2OLIVm6NCh7Ny5k59//jnLPovFkum1YRhZtuXHMZJz+Vlfd955J3feeaf9devWrWnSpAmvv/468+bNy7+gS7Ab1det0Per4OVnnek75nhqoZJC8eyzz/L111+zZs0agoOD7dsDAwMBsvzle/LkySx/IV8rMDAw18dIzuV3ff2di4sLzZo101/P+eR69XWr9P0qePldZ3+n71jhU0IlBcowDIYOHcrSpUtZvXo11apVy7S/WrVqBAYGEhMTY9+WlpbGjz/+SKtWra573pYtW2Y6BmDlypU3PEZurqDqK7v3iY2NJSgoKN9iL4luVl+3St+vglNQdZbd++g7VsgcMRJeSo4hQ4YY/v7+xg8//GAkJCTYHxcvXrSXmT59uuHv728sXbrU+PXXX43HHnvMCAoKMpKSkuxl+vbtazz//PP217/88ovh6upqTJ8+3di7d68xffp0w83NzdiwYUOhXl9xU1D1NWnSJCM6Oto4ePCgsX37dmPAgAGGm5ubsXHjxkK9vuImJ/V15swZY/v27ca3335rAMYnn3xibN++3UhISLCX0fer8BRUnek75nhKqKRAAdk+FixYYC9js9mMiRMnGoGBgYanp6fRtm1b49dff810nnbt2hn9+vXLtO2zzz4zateubbi7uxt16tQxlixZUghXVLwVVH2NGDHCuP322w0PDw+jQoUKRmRkpLFu3bpCuqriKyf1tWDBgmzLTJw40V5G36/CU1B1pu+Y41kMwzAKtg1MREREpHjTGCoRERGRPFJCJSIiIpJHSqhERERE8kgJlYiIiEgeKaESERERySMlVCIiIiJ5pIRKREREJI+UUImIiIjkkRIqERERkTxSQiUi4gTeffddgoOD6dixI3/++ad9e3JyMs2aNaNx48Y0aNCA//znPw6MUqTk0tIzIiJFXHJyMrVr12bJkiUsXrwYb29vXnnlFQDS09O5fPkyPj4+XLx4kfr167N582bKly/v4KhFSha1UIlIida+fXtGjBjh6DBuyNPTkzJlylCzZk2Cg4MpV66cfZ+rqys+Pj4ApKamkp6ejv5OFil8bo4OQETkVnTt2pVLly7x/fffZ9m3fv16WrVqxdatW2nSpIkDostfHh4eDBgwgEqVKlG2bFmOHTuWaf/58+dp164d+/fvZ+bMmQQEBDgoUpGSSy1UIuKUBg0axOrVq/n999+z7Pvggw9o3LhxsUimMqxbt45nn32WixcvEhcXl2lfmTJl2LFjB4cPH+a///1vpjFWIlI4lFCJiFPq0qULFStWJCoqKtP2ixcv8r///Y9BgwYBcPnyZYYNG0bFihXx8vKiTZs2bN68+brnrVq1Kq+99lqmbY0bN2bSpEn21+3bt+fZZ59lxIgRlC1blkqVKvHee+9x4cIFBgwYQOnSpalRowbfffddpvMYhsGMGTOoXr063t7eNGrUiM8///ym13rq1Cm+/fZbhgwZwgMPPMCCBQuyLVepUiUaNmzITz/9dNNzikj+UkIlIk7Jzc2NJ554gqioqExjhj777DPS0tLo3bs3AGPGjGHJkiUsXLiQbdu2ERoaSqdOnTh79mye3n/hwoUEBASwadMmnn32WYYMGUKPHj1o1aoV27Zto1OnTvTt25eLFy/ajxk3bhwLFizg7bffZvfu3Tz33HP06dOHH3/88Ybv9dFHH9GoUSNq165Nnz59+Pjjj7FarQD8+eefJCUlAZCUlMRPP/1E7dq183RtInILDBERJ7V3714DMFavXm3f1rZtW+Oxxx4zDMMwUlJSDHd3d+Pjjz+2709LSzMqV65szJgxwzAMw2jXrp0xfPhw+/4qVaoYr776aqb3adSokTFx4kT763bt2hlt2rSxv75y5YpRqlQpo2/fvvZtCQkJBmCsX7/eHouXl5exbt26TOceNGiQPd7radCggfHaa68ZhmEYVqvVCAgIMJYuXWoYhmFs2bLFaNSokdGwYUOjQYMGxltvvXXDc4lIwdCgdBFxWnXq1KFVq1Z88MEHdOjQgYMHD7J27VpWrlwJwMGDB7FarbRu3dp+jLu7O82bN2fv3r15eu+GDRvan7u6ulK+fHkaNGhg31apUiUATp48CcCePXtITU0lIiIi03nS0tIIDw+/7vts3bqVPXv20KtXL8BsmevZsycLFizgwQcf5I477iA2NjZP1yIieaeESkSc2qBBgxg6dChvvvkmCxYsoEqVKnTs2BHA3hVosVgyHWMYRpZtGVxcXLJMO5DRvXYtd3f3TK8tFkumbRnnt9lsmf799ttvue222zId6+nped3rW7BgAenp6ZmOMQwDFxcXTpw4QWBg4HWPFZHCozFUIuLUHn30UVxdXfnvf//LwoULGTBggD2ZCQ0NxcPDg59//tle3mq1smXLFurWrZvt+SpUqEBCQoL9dVJSEocPH85znGFhYXh6ehIfH09oaGimR0hISLbHXL58mcWLFzN79mxiY2Ptjx07dlC9enU++uijPMclIvlDLVQi4tR8fX3p2bMn//73v0lMTKR///72faVKlWLIkCGMHj2acuXKcfvttzNjxgwuXrxovwvw7+6++26ioqLo2rUrZcuWZfz48bi6uuY5ztKlSzNq1Ciee+45bDYbbdq0ISkpiXXr1uHr60u/fv2yHPPVV1+RkpLCoEGD8Pf3z7TvkUceYcGCBYwaNSrPsYlI3imhEhGnN2jQIObPn09kZCS33357pn3Tp0/HZrPRt29fkpOTadq0KStWrKBs2bLZnmvs2LEcOnSILl264O/vz4svvpgvLVQAL774IhUrVmTatGkcOnSIMmXK0KRJE/79739nW37BggXcc889WZIpgIcffpipU6eyceNGWrRokS/xicit01p+IiIiInmkMVQiIiIieaSESkRERCSPlFCJiIiI5JESKhEREZE8UkIlIiIikkdKqERERETySAmViIiISB4poRIRERHJIyVUIiIiInmkhEpEREQkj5RQiYiIiOSREioRERGRPFJCJSIiIpJH/w9hVgAZwcPxGQAAAABJRU5ErkJggg==", "text/plain": [ "\n", " | Parameter | \n", "Value | \n", "Comment | \n", "
---|---|---|---|
0 | \n", "num_points | \n", "11 | \n", "number of sample points | \n", "
1 | \n", "fit_type | \n", "polynomial | \n", "['polynomial', 'birch', 'birchmurnaghan', 'mur... | \n", "
2 | \n", "fit_order | \n", "3 | \n", "order of the fit polynom | \n", "
3 | \n", "vol_range | \n", "0.1 | \n", "relative volume variation around volume define... | \n", "
4 | \n", "axes | \n", "[x, y, z] | \n", "Axes along which the strain will be applied | \n", "
5 | \n", "strains | \n", "None | \n", "List of strains that should be calculated. If... | \n", "