{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## PbSO4 tutorial: indexing, spacegroup determination & structure solution\n", "In this notebook, you can:\n", "* Load the powder diffraction data and create the PowderPattern object, for two diffrection patterns (X-ray and neutron)\n", "* Find the diffraction peaks and index them (determine the unit cell)\n", "* Perform a profile fit to optimise the background and reflection profiles\n", "* Determine the spacegroup\n", "* Add a molecule to describe the contents of the Crystal structure\n", "* Solve the Crystal structure using a Monte-Carlo/Parallel tempering algorithm, using both X-ray and neutron diffraction patterns\n", "* Save the best result to a CIF file and to Fox .xmlgz format\n", "\n", "Notes:\n", "* This is a simple case, which illustrates the possibility of joint X-ray/neutron optimisation.\n", "* It is important to follow the steps relatively linearly and avoid going back to previous cells until you know better. For example to avoid adding multiple times Scatterer/Molecule objects in the crystal structure, or multiple crystalline phases to the powder pattern with the same crystal, etc..." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Imports" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# 'widget' allows live update and works in both classical notebooks and jupyter-lab.\n", "# Otherwise 'notebook', 'ipympl', 'inline' can be used\n", "%matplotlib widget\n", "\n", "import os\n", "import pyobjcryst\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from pyobjcryst.crystal import *\n", "from pyobjcryst.scatteringpower import ScatteringPowerAtom\n", "from pyobjcryst.atom import Atom\n", "from pyobjcryst.polyhedron import MakeTetrahedron\n", "from pyobjcryst.powderpattern import *\n", "from pyobjcryst.radiation import RadiationType\n", "from pyobjcryst.indexing import *\n", "from pyobjcryst.molecule import *\n", "from pyobjcryst.globaloptim import MonteCarlo\n", "from pyobjcryst.io import xml_cryst_file_save_global\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Create powder pattern object, download data if necessary\n", "We start with the X-ray data which will be used to determine the unit cell and spacegroup." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ " % Total % Received % Xferd Average Speed Time Time Time Current\n", " Dload Upload Total Spent Left Speed\n", " 0 0 0 0 0 0 0 0 --:--:-- --:--:-- --:--:-- 0" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Radiation: X-Ray,Wavelength= tube: Cu, Alpha1/Alpha2= 0.5Imported powder pattern: 6001 points, 2theta= 10.000 -> 160.000, step= 0.025\n", "\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "100 48688 100 48688 0 0 163k 0 --:--:-- --:--:-- --:--:-- 163k\n" ] }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "64f3e219741847358162b2cc58b8d956", "version_major": 2, "version_minor": 0 }, "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4QAAAGQCAYAAAD2lq6fAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjMsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvZiW1igAAAAlwSFlzAAAPYQAAD2EBqD+naQAAch9JREFUeJzt3Xl4FFXePvy7sy8knc0kRMM2IkSCwqCyPoIDskiIir9BjEZ9RRYVMBJAGVxQBxBkHXhUwAVkGZwZARXGCIigyB6ILAYQDRAkISxJZ6XT6T7vHzxVdvWSdCed7k7q/lxXXyRdp6qrTiek7v6eOqURQggQERERERGR6vh4egeIiIiIiIjIMxgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEilGAiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEilGAiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEilGAiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEilGAiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEilGAiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEilGAiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEilGAiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFTKz9M7QA1jMplw8eJFhIWFQaPReHp3iIiIiIgaTAiBsrIyJCQkwMeHNazGxEDYxF28eBGJiYme3g0iIiIiIpfLz8/HLbfc4undaNYYCJu4sLAwADd+WcLDwz28N0REREREDVdaWorExET5XJcaDwNhEycNEw0PD2cgJCIiIqJmhZdENT4OyCUiIiIiIlIp1QXC77//HsOGDUNCQgI0Gg02bdpk1SY3NxepqanQarUICwtDjx49cP78eXm5Xq/HhAkTEBMTg9DQUKSmpuLChQuKbRQXFyM9PR1arRZarRbp6ekoKSlRtDl//jyGDRuG0NBQxMTEYOLEiaiurm6MwyYiIiIiIrKiukBYUVGBO++8E0uXLrW5/Ndff0WfPn3QsWNH7Ny5Ez/99BNee+01BAUFyW0yMjKwceNGrF+/Hrt370Z5eTlSUlJgNBrlNmlpacjJyUFWVhaysrKQk5OD9PR0ebnRaMTQoUNRUVGB3bt3Y/369fj888+RmZnZeAdPRERERERkRiOEEJ7eCU/RaDTYuHEjHnroIfm5kSNHwt/fH6tXr7a5jk6nw0033YTVq1fj0UcfBfDHTJ///e9/MWjQIOTm5uL222/Hvn370L17dwDAvn370LNnT5w8eRIdOnTA119/jZSUFOTn5yMhIQEAsH79ejz99NMoKipy+HrA0tJSaLVa6HQ6XkNIRERE5GFGoxEGg8HTu+H1/P394evra3c5z3Hdh5PKmDGZTNiyZQumTp2KQYMG4ciRI2jbti2mTZsmh8bs7GwYDAYMHDhQXi8hIQHJycnYs2cPBg0ahL1790Kr1cphEAB69OgBrVaLPXv2oEOHDti7dy+Sk5PlMAgAgwYNgl6vR3Z2Nu677z63HTcRERERNYwQAoWFhVaXCJF9ERERiI+P58QxHsZAaKaoqAjl5eV455138Pe//x1z5sxBVlYWhg8fju+++w59+/ZFYWEhAgICEBkZqVg3Li4OhYWFAIDCwkLExsZabT82NlbRJi4uTrE8MjISAQEBchtb9Ho99Hq9/H1paWm9j5eIiIiIXEMKg7GxsQgJCWHIqYUQApWVlSgqKgIAtGzZ0sN7pG4MhGZMJhMA4MEHH8RLL70EAOjSpQv27NmDDz74AH379rW7rhBC8Ytv6z+B+rSxNHv2bLz55pt1HwwRERERuYXRaJTDYHR0tKd3p0kIDg4GcKMgExsbW+vwUWpcqptUpjYxMTHw8/PD7bffrng+KSlJnmU0Pj4e1dXVKC4uVrQpKiqSK37x8fG4dOmS1fYvX76saGNZCSwuLobBYLCqHJqbNm0adDqd/MjPz3f+QImIiIjIZaRrBkNCQjy8J02L1F+85tKzGAjNBAQE4O6778apU6cUz58+fRqtW7cGAHTr1g3+/v7Ytm2bvLygoADHjx9Hr169AAA9e/aETqfDgQMH5Db79++HTqdTtDl+/DgKCgrkNlu3bkVgYCC6detmdx8DAwPlm9DzZvRERERE3oPDRJ3D/vIOqhsyWl5ejjNnzsjf5+XlIScnB1FRUWjVqhWmTJmCRx99FPfeey/uu+8+ZGVl4auvvsLOnTsBAFqtFqNGjUJmZiaio6MRFRWFyZMno3PnzhgwYACAGxXFwYMHY/To0Vi2bBkAYMyYMUhJSUGHDh0AAAMHDsTtt9+O9PR0vPvuu7h27RomT56M0aNHM+QREREREZFbqK5CeOjQIXTt2hVdu3YFAEyaNAldu3bF66+/DgB4+OGH8cEHH2Du3Lno3LkzPvzwQ3z++efo06ePvI2FCxfioYcewogRI9C7d2+EhITgq6++Uox9Xrt2LTp37oyBAwdi4MCBuOOOOxS3svD19cWWLVsQFBSE3r17Y8SIEXjooYcwb948N/UEeZNTp05h1qxZnt4NIiIiItnOnTuh0Wg4c2ozp+r7EDYHvEdL89C+fXucOXMG/HUkIiJqeq5fv468vDy0bdsWQUFBnt4dl9m5cyfuu+8+FBcXIyIiwuXbr63feI7rPqqrEBIREREREdENDIREXoAXVRMREZEn6PV6TJw4EbGxsQgKCkKfPn1w8OBBRZsff/wRd955J4KCgtC9e3ccO3ZMXnbu3DkMGzYMkZGRCA0NRadOnfDf//7X3YdBDcBASERERESkUlOnTsXnn3+OVatW4fDhw7j11lsxaNAgXLt2TW4zZcoUzJs3DwcPHkRsbCxSU1PlW0W88MIL0Ov1+P7773Hs2DHMmTMHLVq08NThUD2obpZRIm/ECiEREVHzU1lZiZMnT7r9dTt27OjQPRErKirw/vvvY+XKlRgyZAgAYMWKFdi2bRs++ugj3H333QCAN954A/fffz8AYNWqVbjllluwceNGjBgxAufPn8cjjzyCzp07AwDatWvXSEdFjYWBkIiIiIioEZw8ebLW+0s3luzsbPz5z3+us92vv/4Kg8GA3r17y8/5+/vjnnvuQW5urhwIe/bsKS+PiopChw4dkJubCwCYOHEinnvuOWzduhUDBgzAI488gjvuuMPFR0SNiYGQyAuwQkhERNT8dOzYEdnZ2R55XUdIs5tbnocIIeo8N5GWP/vssxg0aBC2bNmCrVu3Yvbs2Zg/fz4mTJhQjz0nT2AgJCIiIiJqBCEhIQ5V6jzl1ltvRUBAAHbv3o20tDQAgMFgwKFDh5CRkSG327dvH1q1agUAKC4uxunTpxWhMzExEePGjcO4ceMwbdo0rFixgoGwCWEgJCIiIiJSodDQUDz33HOYMmUKoqKi0KpVK8ydOxeVlZUYNWoUfvrpJwDAW2+9hejoaMTFxWH69OmIiYnBQw89BADIyMjAkCFDcNttt6G4uBg7duxAUlKSB4+KnMVASOQFOGSUiIiIPOGdd96ByWRCeno6ysrKcNddd+Gbb75BZGSkos2LL76IX375BXfeeSe+/PJLBAQEAACMRiNeeOEFXLhwAeHh4Rg8eDAWLlzoqcOhetAIafAwNUmlpaXQarXQ6XQIDw/39O5QPSUlJeHkyZPgryMREVHTc/36deTl5aFt27YICgry9O40GbX1G89x3Yf3ISTyAlKFkIGQiIiIiNyJgZDIC3DIKBERERF5AgMhkRdhhZCIiIiI3ImBkMiLMBASERERkTsxEBJ5EQZCIiIiInInBkIiL8BJZYiIiIjIExgIibwAAyEREREReQIDIRERERERkUoxEBJ5EVYIiYiIqKmYMWMGunTp4undoAZiICTyIgyERERERORODIREXoDXEBIRERGRJzAQEnkBBkIiIiLyBJPJhDlz5uDWW29FYGAgWrVqhZkzZwIAXn75Zdx2220ICQlBu3bt8Nprr8FgMNS6vY8//hidOnVCYGAgWrZsifHjx7vjMKgB/Dy9A0T0BwZCIiIicqdp06ZhxYoVWLhwIfr06YOCggKcPHkSABAWFoaVK1ciISEBx44dw+jRoxEWFoapU6fa3Nb777+PSZMm4Z133sGQIUOg0+nw448/uvNwqB4YCIm8gFQhJCIiInKXsrIyLF68GEuXLsVTTz0FAPjTn/6EPn36AABeffVVuW2bNm2QmZmJzz77zG4g/Pvf/47MzEy8+OKL8nN33313Ix4BuQIDIZEXYYWQiIioGamsBP6v2uZWHTsCISF1NsvNzYVer0f//v1tLv/Pf/6DRYsW4cyZMygvL0dNTQ3Cw8Ntti0qKsLFixftbou8l+quIfz+++8xbNgwJCQkQKPRYNOmTXbbjh07FhqNBosWLVI8r9frMWHCBMTExCA0NBSpqam4cOGCok1xcTHS09Oh1Wqh1WqRnp6OkpISRZvz589j2LBhCA0NRUxMDCZOnIjq6moXHSk1JbyGkIiIqBk6eRLo1s39DwdDaHBwsN1l+/btw8iRIzFkyBBs3rwZR44cwfTp0+2eq9a2LfJuqqsQVlRU4M4778T/9//9f3jkkUfsttu0aRP279+PhIQEq2UZGRn46quvsH79ekRHRyMzMxMpKSnIzs6Gr68vACAtLQ0XLlxAVlYWAGDMmDFIT0/HV199BQAwGo0YOnQobrrpJuzevRtXr17FU089BSEElixZ0ghHTk0BAyEREVEz0rEjkJ3tmdd1QPv27REcHIxvv/0Wzz77rGLZjz/+iNatW2P69Onyc+fOnbO7rbCwMLRp0wbffvst7rvvvvrtN3mE6gLhkCFDMGTIkFrb/P777xg/fjy++eYbDB06VLFMp9Pho48+wurVqzFgwAAAwJo1a5CYmIjt27dj0KBByM3NRVZWFvbt24fu3bsDAFasWIGePXvi1KlT6NChA7Zu3Yqff/4Z+fn5cuicP38+nn76acycOdNuOZ6aNwZCIiKiZiQkBPjznz29F3YFBQXh5ZdfxtSpUxEQEIDevXvj8uXLOHHiBG699VacP38e69evx913340tW7Zg48aNtW5vxowZGDduHGJjYzFkyBCUlZXhxx9/xIQJE9x0RFQfqhsyWheTyYT09HRMmTIFnTp1slqenZ0Ng8GAgQMHys8lJCQgOTkZe/bsAQDs3bsXWq1WDoMA0KNHD2i1WkWb5ORkRQVy0KBB0Ov1yK7lkyS9Xo/S0lLFg5o+DhklIiIiT3jttdeQmZmJ119/HUlJSXj00UdRVFSEBx98EC+99BLGjx+PLl26YM+ePXjttddq3dZTTz2FRYsW4b333kOnTp2QkpKCX375xU1HQvWlugphXebMmQM/Pz9MnDjR5vLCwkIEBAQgMjJS8XxcXBwKCwvlNrGxsVbrxsbGKtrExcUplkdGRiIgIEBuY8vs2bPx5ptvOnVMRERERES2+Pj4YPr06YqhoZK5c+di7ty5iucyMjLkr2fMmIEZM2Yolo8dOxZjx45tjF2lRsIKoZns7GwsXrwYK1eudPo2AEIIxTq21q9PG0vTpk2DTqeTH/n5+U7tJ3knVgiJiIiIyBMYCM388MMPKCoqQqtWreDn5wc/Pz+cO3cOmZmZaNOmDQAgPj4e1dXVKC4uVqxbVFQkV/zi4+Nx6dIlq+1fvnxZ0cayElhcXAyDwWBVOTQXGBiI8PBwxYOaPgZCIiIiIvIEBkIz6enpOHr0KHJycuRHQkICpkyZgm+++QYA0K1bN/j7+2Pbtm3yegUFBTh+/Dh69eoFAOjZsyd0Oh0OHDggt9m/fz90Op2izfHjx1FQUCC32bp1KwIDA9GtWzd3HC55IQZCIiIiInIn1V1DWF5ejjNnzsjf5+XlIScnB1FRUWjVqhWio6MV7f39/REfH48OHToAALRaLUaNGoXMzExER0cjKioKkydPRufOneVZR5OSkjB48GCMHj0ay5YtA3DjthMpKSnydgYOHIjbb78d6enpePfdd3Ht2jVMnjwZo0ePZtVPxRgIiYiIiMidVFchPHToELp27YquXbsCACZNmoSuXbvi9ddfd3gbCxcuxEMPPYQRI0agd+/eCAkJwVdffSXfgxAA1q5di86dO2PgwIEYOHAg7rjjDqxevVpe7uvriy1btiAoKAi9e/fGiBEj8NBDD2HevHmuO1hqchgIiYiIiMidNIJnoE1aaWkptFotdDodK4tN2N13341Dhw6hqKgIN910k6d3h4iIiJxw/fp15OXloU2bNggODvb07jQZVVVVOHv2LNq2bYugoCDFMp7juo/qKoRE3oiTyhARETVd/v7+AIDKykoP70nTIvWX1H/kGaq7hpDIHU6cOIFbb70VgYGBTq3HQEhERNT0+Pr6IiIiAkVFRQCAkJAQp29hpiZCCFRWVqKoqAgRERGKy67I/RgIiVxMCIHk5GSMGTNGnlTImXWJiIio6YmPjwcAORRS3SIiIuR+I89hICRqJIcPH3Z6HQZCIiKipkmj0aBly5aIjY2FwWDw9O54PX9/f1YGvQQDIZGLSaGupqbG4XV4DSEREVHz4Ovry6BDTQonlSFqJEaj0eG2vM6AiIiIiDyBgZDIxepTIbRcl4iIiIjIHRgIiRpJfSqEDIRERERE5E4MhEQuJoU6ZwKh5bpERERERO7AQEjUSFghJCIiIiJvx0BI5GINCXUMhERERETkTgyEpDqlpaX1Gs5JRERERNTcMBCS6mi1WowfP77Rtl+fKh+HjBIRERGRJzAQkipt3Lix0V+jPuGOgZCIiIiI3ImBkFSpMYeMskJIRERERE0FAyGpkrddQ8hASERERESewEBIquRtFUJXrEtERERE5CwGQlIlb6sQEhERERF5AgMhqRIrhEREREREDISkUjU1NY3+Gs6EO15DSERERESewEBIqtSYwYuzjNZt7dq10Ol0nt4NIiIiItVjICTyImoIhOXl5XjiiScwZswYT+8KERERkeoxEBJ5ETUEQh+fG//tFBUVeXhPiIiIiIiBkMjF1BDqXKG6utrTu0BERESkegyEpEreFtrUdA2hdIwGg8HDe0JEREREqguE33//PYYNG4aEhARoNBps2rRJXmYwGPDyyy+jc+fOCA0NRUJCAp588klcvHhRsQ29Xo8JEyYgJiYGoaGhSE1NxYULFxRtiouLkZ6eDq1WC61Wi/T0dJSUlCjanD9/HsOGDUNoaChiYmIwceJEVk2aAU4q4xgGQiIiIiLPU10grKiowJ133omlS5daLausrMThw4fx2muv4fDhw9iwYQNOnz6N1NRURbuMjAxs3LgR69evx+7du1FeXo6UlBTFve3S0tKQk5ODrKwsZGVlIScnB+np6fJyo9GIoUOHoqKiArt378b69evx+eefIzMzs/EOnmTeGry8db9cSTpGNRwrERERkbfz8/QOuNuQIUMwZMgQm8u0Wi22bdumeG7JkiW45557cP78ebRq1Qo6nQ4fffQRVq9ejQEDBgAA1qxZg8TERGzfvh2DBg1Cbm4usrKysG/fPnTv3h0AsGLFCvTs2ROnTp1Chw4dsHXrVvz888/Iz89HQkICAGD+/Pl4+umnMXPmTISHhzdiL1Bjqk/gUWOFkIiIiIg8T3UVQmfpdDpoNBpEREQAALKzs2EwGDBw4EC5TUJCApKTk7Fnzx4AwN69e6HVauUwCAA9evSAVqtVtElOTpbDIAAMGjQIer0e2dnZdvdHr9ejtLRU8SBqShh6iYiIiLwHA2Etrl+/jldeeQVpaWlyxa6wsBABAQGIjIxUtI2Li0NhYaHcJjY21mp7sbGxijZxcXGK5ZGRkQgICJDb2DJ79mz5ukStVovExMQGHaNaeduN6V2xblOjpmMlIiIi8lYMhHYYDAaMHDkSJpMJ7733Xp3thRDysD8Aiq8b0sbStGnToNPp5Ed+fn6d+0beT01DRtVwjERERERNBQOhDQaDASNGjEBeXh62bdumuJ4vPj4e1dXVKC4uVqxTVFQkV/zi4+Nx6dIlq+1evnxZ0cayElhcXAyDwWBVOTQXGBiI8PBwxYO8CwMPERERETUVDIQWpDD4yy+/YPv27YiOjlYs79atG/z9/RWTzxQUFOD48ePo1asXAKBnz57Q6XQ4cOCA3Gb//v3Q6XSKNsePH0dBQYHcZuvWrQgMDES3bt0a8xDJi6khTKrhGImIiIiaCtXNMlpeXo4zZ87I3+fl5SEnJwdRUVFISEjA//t//w+HDx/G5s2bYTQa5SpeVFQUAgICoNVqMWrUKGRmZiI6OhpRUVGYPHkyOnfuLM86mpSUhMGDB2P06NFYtmwZAGDMmDFISUlBhw4dAAADBw7E7bffjvT0dLz77ru4du0aJk+ejNGjR7Pq18Q15LYKagpLajpWIiIiIm+lukB46NAh3HffffL3kyZNAgA89dRTmDFjBr788ksAQJcuXRTrfffdd+jXrx8AYOHChfDz88OIESNQVVWF/v37Y+XKlfD19ZXbr127FhMnTpRnI01NTVXc+9DX1xdbtmzB888/j969eyM4OBhpaWmYN29eYxw2eTleQ0hEREREnqC6QNivX79aT0gdOVkNCgrCkiVLsGTJErttoqKisGbNmlq306pVK2zevLnO16OmhYGndlL/1DZ5EhERERG5B68hJPIiagqTajpWIiIiIm/FQEiq4o4QUp/X4JBRIiIiIvIEBkJSFYYRIiIiIqI/MBASuRhnGa2dGo6RiIiIqKlgICRV8dYwoqYhoxI1HSsRERGRt2IgJHKxhgQdNYQkNRwjERERUVPBQEiq4s4wwuBDRERERN6OgZBUxVtnGXXFuk2FGo6RiIiIqKlgICTyAryGkIiIiIg8gYGQVMXbK4RqwP4hIiIi8h4MhERehGGJiIiIiNyJgZBUxVsrhGoaMqqGYyQiIiJqKhgISVU4y6jnsV+IiIiIvAcDIZGLcZZRIiIiImoqGAhJVbw1cHHIKBERERF5AgMhqYq3hxFv3z9XUtOxEhEREXkrBkIiF2PQqR37h4iIiMh7MBCSqnh7GPH2/SMiIiKi5oWBkMjFpFDnTLjjNYRERERE5AkMhKQqDCPeg+8FERERkecxEJKqeOuN6V2xblOhhmMkIiIiaioYCIncoKCgAFVVVXaXq2nIKBERERF5DwZCUhVPVQgTEhLwwAMP1Gvd5kYNx0hERETUVDAQErnJzp077S6TKoRqwEBIRERE5D1UFwi///57DBs2DAkJCdBoNNi0aZNiuRACM2bMQEJCAoKDg9GvXz+cOHFC0Uav12PChAmIiYlBaGgoUlNTceHCBUWb4uJipKenQ6vVQqvVIj09HSUlJYo258+fx7BhwxAaGoqYmBhMnDgR1dXVjXHY9H94DSERERER0R9UFwgrKipw5513YunSpTaXz507FwsWLMDSpUtx8OBBxMfH4/7770dZWZncJiMjAxs3bsT69euxe/dulJeXIyUlBUajUW6TlpaGnJwcZGVlISsrCzk5OUhPT5eXG41GDB06FBUVFdi9ezfWr1+Pzz//HJmZmY138OT1gcvb988V1HCMRERERE2Fn6d3wN2GDBmCIUOG2FwmhMCiRYswffp0DB8+HACwatUqxMXFYd26dRg7dix0Oh0++ugjrF69GgMGDAAArFmzBomJidi+fTsGDRqE3NxcZGVlYd++fejevTsAYMWKFejZsydOnTqFDh06YOvWrfj555+Rn5+PhIQEAMD8+fPx9NNPY+bMmQgPD3dDb1BjqM99CNWI/UNERETkeaqrENYmLy8PhYWFGDhwoPxcYGAg+vbtiz179gAAsrOzYTAYFG0SEhKQnJwst9m7dy+0Wq0cBgGgR48e0Gq1ijbJyclyGASAQYMGQa/XIzs7u1GPU828PYR4+/65ghqOkYiIiKipUF2FsDaFhYUAgLi4OMXzcXFxOHfunNwmICAAkZGRVm2k9QsLCxEbG2u1/djYWEUby9eJjIxEQECA3MYWvV4PvV4vf19aWuro4ZGb1Cfw8LYTREREROQJrBDaYDnjoxCizlkgLdvYal+fNpZmz54tT1Sj1WqRmJhY636REgOX5/E9ICIiIvIeDIRm4uPjAcCqQldUVCRX8+Lj41FdXY3i4uJa21y6dMlq+5cvX1a0sXyd4uJiGAwGq8qhuWnTpkGn08mP/Px8J49S3TjLKBERERHRHxgIzbRt2xbx8fHYtm2b/Fx1dTV27dqFXr16AQC6desGf39/RZuCggIcP35cbtOzZ0/odDocOHBAbrN//37odDpFm+PHj6OgoEBus3XrVgQGBqJbt2529zEwMBDh4eGKBzV9ahoyqoZjJCIiImoqVHcNYXl5Oc6cOSN/n5eXh5ycHERFRaFVq1bIyMjArFmz0L59e7Rv3x6zZs1CSEgI0tLSAABarRajRo1CZmYmoqOjERUVhcmTJ6Nz587yrKNJSUkYPHgwRo8ejWXLlgEAxowZg5SUFHTo0AEAMHDgQNx+++1IT0/Hu+++i2vXrmHy5MkYPXo0Q14jcmeFsD6vpaawpKZjJSIiIvJWqguEhw4dwn333Sd/P2nSJADAU089hZUrV2Lq1KmoqqrC888/j+LiYnTv3h1bt25FWFiYvM7ChQvh5+eHESNGoKqqCv3798fKlSvh6+srt1m7di0mTpwoz0aampqquPehr68vtmzZgueffx69e/dGcHAw0tLSMG/evMbuAlXz1hBS1zWqzYm3vgdEREREaqQRPDtr0kpLS6HVaqHT6VhZdMDly5flGWAb60f/t99+w5/+9CfExsbK15LWNST0wQcfxJdffonNmzdj6NChjbJf3uLXX3/Frbfeig4dOuDkyZOe3h0iIiLyQjzHdR9eQ0iq4u2ff3j7/rmCGo6RiIiIqKlgICRyMQae2rF/iIiIiLwHAyGpireHEW/fP1dS07ESEREReSsGQlIVb51llLedICIiIiJPYCAk8gJqCoRERERE5D0YCElV3FkhJNvYP0RERETeg4GQyIswLBERERGROzEQkqp4e4VQDYFQDcdIRERE1FQwEJKqeGsYka4hVBNvfS+IiIiI1ISBkMjF6jPLqOW6zZkajpGIiIioqWAgJFXx9jDi7ftHRERERM0LAyERuRVDLxEREZH3YCAkVeGkMp6nhmMkIiIiaioYCElVvDWM8Mb0REREROQJDIRELlafUKemQKiGYyQiIiJqKhgISVXcGUYYfGrH/iEiIiLyPAZCIhfjNYS1U8MxEhERETUVDISkKt4eRrx9/4iIiIioeWEgJFXx1llGpWsI1YChl4iIiMh7MBASeRE1hSU1HSsRERGRt2IgJFXx1gqhK9ZtKtRwjERERERNBQMhqYq3hhE13XaCiIiIiLwHAyGRizHU1Y79Q0REROQ9GAhJVbz9PoRqCktqOlYiIiIib8VASORiDZllVA0hSQ3HSERERNRUMBBaqKmpwauvvoq2bdsiODgY7dq1w1tvvQWTySS3EUJgxowZSEhIQHBwMPr164cTJ04otqPX6zFhwgTExMQgNDQUqampuHDhgqJNcXEx0tPTodVqodVqkZ6ejpKSEnccpmo5E0ZKSkqwcePGRtwbdZLeAzXdaoOIiIjIWzEQWpgzZw4++OADLF26FLm5uZg7dy7effddLFmyRG4zd+5cLFiwAEuXLsXBgwcRHx+P+++/H2VlZXKbjIwMbNy4EevXr8fu3btRXl6OlJQUGI1GuU1aWhpycnKQlZWFrKws5OTkID093a3HqzbOBMIpU6Zg+PDhTod0zjLqGDUdKxEREZG38vP0DnibvXv34sEHH8TQoUMBAG3atME///lPHDp0CMCNk9hFixZh+vTpGD58OABg1apViIuLw7p16zB27FjodDp89NFHWL16NQYMGAAAWLNmDRITE7F9+3YMGjQIubm5yMrKwr59+9C9e3cAwIoVK9CzZ0+cOnUKHTp08MDRk7mKigoAwLVr1xAREeGW11RDSFLDMRIRERE1FawQWujTpw++/fZbnD59GgDw008/Yffu3XjggQcAAHl5eSgsLMTAgQPldQIDA9G3b1/s2bMHAJCdnQ2DwaBok5CQgOTkZLnN3r17odVq5TAIAD169IBWq5Xb2KLX61FaWqp4kOOcCSP1HdLIwENERERETQUrhBZefvll6HQ6dOzYEb6+vjAajZg5cyYee+wxAEBhYSEAIC4uTrFeXFwczp07J7cJCAhAZGSkVRtp/cLCQsTGxlq9fmxsrNzGltmzZ+PNN9+s/wGSwxp6jRtnGbVNOkbz63KJiIiIyDNYIbTw2WefYc2aNVi3bh0OHz6MVatWYd68eVi1apWinWVYEELUGSAs29hqX9d2pk2bBp1OJz/y8/MdOSz6P+4IabyG0DEGg8HTu0BERESkeqwQWpgyZQpeeeUVjBw5EgDQuXNnnDt3DrNnz8ZTTz2F+Ph4ADcqfC1btpTXKyoqkquG8fHxqK6uRnFxsaJKWFRUhF69esltLl26ZPX6ly9ftqo+mgsMDERgYGDDD1SlvLVq15DbTuTn58Pf31/+2fR20jHW1NR4eE+IiIiIiBVCC5WVlfDxUXaLr6+vPLytbdu2iI+Px7Zt2+Tl1dXV2LVrlxz2unXrBn9/f0WbgoICHD9+XG7Ts2dP6HQ6HDhwQG6zf/9+6HQ6uQ15B3dWCOujVatWig8nmgpWCImIiIg8jxVCC8OGDcPMmTPRqlUrdOrUCUeOHMGCBQvwzDPPALhRycnIyMCsWbPQvn17tG/fHrNmzUJISAjS0tIAAFqtFqNGjUJmZiaio6MRFRWFyZMno3PnzvKso0lJSRg8eDBGjx6NZcuWAQDGjBmDlJQUzjDaiOozqYw7r3VTw5BR6RilQFhZWYlhw4Zh5cqVSExM9OSuEREREakOA6GFJUuW4LXXXsPzzz+PoqIiJCQkYOzYsXj99dflNlOnTkVVVRWef/55FBcXo3v37ti6dSvCwsLkNgsXLoSfnx9GjBiBqqoq9O/fHytXroSvr6/cZu3atZg4caI8G2lqaiqWLl3qvoMlh7ijQiito4ZAKJEC4cGDB7Fjxw4sXboUc+bM8fBeEREREakLA6GFsLAwLFq0CIsWLbLbRqPRYMaMGZgxY4bdNkFBQViyZInihvaWoqKisGbNmgbsLTnLndcQqincOcPyGkL2ExEREZHn8BpCUhVPzDLqyPpqrBBaTirT0Nt8EBEREZHzGAiJ7HDVNYQMhEr1CcxERERE1DgYCElV6jOpTH0rhAw6ttnrF1YIiYiIiNyPgZCoDg0NdqwQ1k6Nx0xERETkLRgISVU8cWN6BkIle0NGWSEkIiIicj8GQlIVT0wq4651mzoGQiIiIiL3YyAkssMTk8qoASeVISIiIvIeDISkKu6sEDqznpqGjNrDCiERERGR+zEQkqq4MxDWZ5ZSNQRCNRwjERERUVPBQEhUB3femJ6IiIiIyJ0YCElV6nMfwoZeQ+gMNQRDXkNIRERE5D0YCInq4M4KoZrDEa8hJCIiInI/BkJSFU/cdoKBUIkVQiIiIiLvwUBIquKJSWVIiTemJyIiIvIeDIREdajvNYTOVP3UVCG0h4GQiIiIyP0YCElV6jOpDIeMuhaHjBIRERF5DwZCojpwyKh7sL+IiIiI3I+BkFTFHRXC+rwmK4RERERE5AkMhKQqnphl1Jl11ByOWCEkIiIicj8GQiI76ntjelbAasf+ISIiIvIeDISkKrwPIRERERHRHxgIiergjpCmpkCohmMkIiIiaioYCElVvLVCWN/Xag44KysRERGR5zAQkqrUJ3C54xpCNQVBe8fKQEhERETkfgyERHY09Mb0zgQcNQ0ZtaTGYyYiIiLyFgyENvz+++944oknEB0djZCQEHTp0gXZ2dnyciEEZsyYgYSEBAQHB6Nfv344ceKEYht6vR4TJkxATEwMQkNDkZqaigsXLijaFBcXIz09HVqtFlqtFunp6SgpKXHHIaqWO4aMWq7HIaNKrBASEREReQ8GQgvFxcXo3bs3/P398fXXX+Pnn3/G/PnzERERIbeZO3cuFixYgKVLl+LgwYOIj4/H/fffj7KyMrlNRkYGNm7ciPXr12P37t0oLy9HSkoKjEaj3CYtLQ05OTnIyspCVlYWcnJykJ6e7s7DJQe48z6EasDbThARERF5Dz9P74C3mTNnDhITE/HJJ5/Iz7Vp00b+WgiBRYsWYfr06Rg+fDgAYNWqVYiLi8O6deswduxY6HQ6fPTRR1i9ejUGDBgAAFizZg0SExOxfft2DBo0CLm5ucjKysK+ffvQvXt3AMCKFSvQs2dPnDp1Ch06dHDfQauIOyeVqc+QUzWGI04qQ0REROQ5rBBa+PLLL3HXXXfhr3/9K2JjY9G1a1esWLFCXp6Xl4fCwkIMHDhQfi4wMBB9+/bFnj17AADZ2dkwGAyKNgkJCUhOTpbb7N27F1qtVg6DANCjRw9otVq5jS16vR6lpaWKBznOmcDlzhvTq+kaQvNjFEKo4piJiIiIvBUDoYXffvsN77//Ptq3b49vvvkG48aNw8SJE/Hpp58CAAoLCwEAcXFxivXi4uLkZYWFhQgICEBkZGStbWJjY61ePzY2Vm5jy+zZs+VrDrVaLRITE+t/sOQQ3oew8Qgh5MDNCiERERGR+zEQWjCZTPjzn/+MWbNmoWvXrhg7dixGjx6N999/X9HO8uRVCFHnCa1lG1vt69rOtGnToNPp5Ed+fr4jh0X/x9vvQ6gGlhVCBkIiIiIiz2EgtNCyZUvcfvvtiueSkpJw/vx5AEB8fDwAWFXxioqK5KphfHw8qqurUVxcXGubS5cuWb3+5cuXraqP5gIDAxEeHq54kOPqM2TUHYFQzRVCtR0zERERkTdhILTQu3dvnDp1SvHc6dOn0bp1awBA27ZtER8fj23btsnLq6ursWvXLvTq1QsA0K1bN/j7+yvaFBQU4Pjx43Kbnj17QqfT4cCBA3Kb/fv3Q6fTyW3Is+p7DWF9qCkQ2qsQEhEREZH7cZZRCy+99BJ69eqFWbNmYcSIEThw4ACWL1+O5cuXA7gREjIyMjBr1iy0b98e7du3x6xZsxASEoK0tDQAgFarxahRo5CZmYno6GhERUVh8uTJ6Ny5szzraFJSEgYPHozRo0dj2bJlAIAxY8YgJSWFM4w2Ig4Z9S4mk0m+FYuPDz+fIiIiInI3BkILd999NzZu3Ihp06bhrbfeQtu2bbFo0SI8/vjjcpupU6eiqqoKzz//PIqLi9G9e3ds3boVYWFhcpuFCxfCz88PI0aMQFVVFfr374+VK1fC19dXbrN27VpMnDhRno00NTUVS5cudd/Bqpwj131K7ZzdrrPrq7lCOG7cOABQ/G4QERERkXswENqQkpKClJQUu8s1Gg1mzJiBGTNm2G0TFBSEJUuWYMmSJXbbREVFYc2aNQ3ZVXKSZRhxRyB0Zh01BEJzQgj51imsEBIRERG5H8/ASFXMA5ej1641NKQ1Rsj77LPPUF1d7fLtuoNlKJcwEBIRERG5H8/ASLXqCmquujG9M+s4sm52djZGjhyJefPmOf063sBeICQiIiIi92MgJFVxR4Wwsa8hrKysBACr25o0RebvAcMhERERkfsxEJKquDMQ1uc+ho3V1pvYqxA21eMhIiIiasoYCElV3HkNYWOFO0cmwmkqGAiJiIiIPIuBkFTFmQDiqmsI1VD1c4ble/Dwww9bPU9ERERE7sFASKrirUNG1XzbCbUeOxEREZE3YCAkVfHEpDKNtU5TZfkeMBASEREReQ4DIamKt9+HUG2hiBVCIiIiIs9iICRV8dZrCOsTippqgLJ8D6Tvne1nIiIiImo4BkJSFWcqhPW5bUR92puv48i69d0vb8QKIREREZFnMRCSqpiHwLoCoauCiquDTlO/7YS9CiEDIREREZH7MRCSqjgzZLS+QaWxK4TNCSeVISIiIvIsBkJSFWeGjLoqEPK2E0qsEBIRERF5DwZCUpX6BML6TirTnK71cyUGQiIiIiLvwUBIquKOCqGt13K0rRpmGTXHQEhERETkWQyEpCqNcQ1hVlYW9u3bZ3e7rh4y2tRvz8AKIREREZH38PP0DhC5U2NUCIcMGWKzfWMNGW3qgdAcJ5UhIiIi8ixWCElV3HkNYX3WaUiFsKkEKssKoXQ8TWX/iYiIiJoTBkJSFXfcdqK213RFW3uBsClWDjlklIiIiMizGAhJVZrDbSeaW4WQgZCIiIjIcxgISVU8EQhdvY69IZZNvULYFPefiIiIqKljICRV8dYKoTNtWSEkIiIiIldhICRVceYaQok7KleuGDLaFCtsnGWUiIiIyLMYCOswe/ZsaDQaZGRkyM8JITBjxgwkJCQgODgY/fr1w4kTJxTr6fV6TJgwATExMQgNDUVqaiouXLigaFNcXIz09HRotVpotVqkp6ejpKTEDUelXt5aIazPkNGGbMOTWCEkIiIi8h4MhLU4ePAgli9fjjvuuEPx/Ny5c7FgwQIsXboUBw8eRHx8PO6//36UlZXJbTIyMrBx40asX78eu3fvRnl5OVJSUmA0GuU2aWlpyMnJQVZWFrKyspCTk4P09HS3HZ8aeWsgdKatt19DePHiRfz88892lzMQEhEREXkPBkI7ysvL8fjjj2PFihWIjIyUnxdCYNGiRZg+fTqGDx+O5ORkrFq1CpWVlVi3bh0AQKfT4aOPPsL8+fMxYMAAdO3aFWvWrMGxY8ewfft2AEBubi6ysrLw4YcfomfPnujZsydWrFiBzZs349SpUx45ZjUwD03uuu2EI5rTLKNt2rRBp06dHGrLQEhERETkWQyEdrzwwgsYOnQoBgwYoHg+Ly8PhYWFGDhwoPxcYGAg+vbtiz179gAAsrOzYTAYFG0SEhKQnJwst9m7dy+0Wi26d+8ut+nRowe0Wq3chlzPEzemd/VtJ+y19ZYKocFgqHU5K4RERERE3sPP0zvgjdavX4/Dhw/j4MGDVssKCwsBAHFxcYrn4+LicO7cOblNQECAorIotZHWLywsRGxsrNX2Y2Nj5Ta26PV66PV6+fvS0lIHj4oA7x8y6ghvrxDWxfI9YCAkIiIi8hxWCC3k5+fjxRdfxJo1axAUFGS3nUajUXwvhLB6zpJlG1vt69rO7Nmz5UlotFotEhMTa31NUnJnIKzr56G+r+Xt1xDWhRVCIiIiIu/BQGghOzsbRUVF6NatG/z8/ODn54ddu3bhH//4B/z8/OTKoGUVr6ioSF4WHx+P6upqFBcX19rm0qVLVq9/+fJlq+qjuWnTpkGn08mP/Pz8Bh2v2jhz24mGBhVXDAO1pTlVCBkIiYiIiDyLgdBC//79cezYMeTk5MiPu+66C48//jhycnLQrl07xMfHY9u2bfI61dXV2LVrF3r16gUA6NatG/z9/RVtCgoKcPz4cblNz549odPpcODAAbnN/v37odPp5Da2BAYGIjw8XPEgxzWHIaPefg1hXSwn9mEgJCIiIvIcXkNoISwsDMnJyYrnQkNDER0dLT+fkZGBWbNmoX379mjfvj1mzZqFkJAQpKWlAQC0Wi1GjRqFzMxMREdHIyoqCpMnT0bnzp3lSWqSkpIwePBgjB49GsuWLQMAjBkzBikpKejQoYMbj1hdPDGpjDPrNOSehU0lEFpWCO0NgSUiIiKixsdAWA9Tp05FVVUVnn/+eRQXF6N79+7YunUrwsLC5DYLFy6En58fRowYgaqqKvTv3x8rV66Er6+v3Gbt2rWYOHGiPBtpamoqli5d6vbjURN3DBm1XM9d9yFsKoHK3qQyTSXQEhERETUnDIQO2Llzp+J7jUaDGTNmYMaMGXbXCQoKwpIlS7BkyRK7baKiorBmzRoX7SU5wh1DRuujPqHRPExJ3zcFHDJKRERE5D14DSGpijuvIWysSWXsVTmbSqDipDJERERE3oOBkFTFndcQunvIaFOpENoLhE1l/4mIiIiaEwZCUpX6VNQaOsuoM+s4W01sihVCe0NGKysrPbVLRERERKrFQEiq4s5rCG2tb29b9akiWgbCplJhszepTGlpqad2iYiIiEi1GAhJVTx9DWFdgdDZIaNNsUJob8iowWDw1C4RERERqRYDIamKp+9DWNc9BBsyZLSpVAjtDRk1Go2e2iUiIiIi1WIgJFVxx30Ia1u/MYeMeluF0JFjZSAkIiIi8iwGQlIV8+pUfYaMVldX4/7778epU6ccXs+R0FafCqG334fQ3v44GghPnTqFvLy8xttBIiIiImIgJHVxJkDZCmkXL17E9u3bMX36dIdew9b3ta3jzDWEloHQ2yqENTU1Np+3DOX2AmHHjh3Rrl27xttBIiIiImIgJHVp6JBRf39/AIBer3d4PWcqhI5oKhVCe0NAOWSUiIiIyHswEJKq1KdC6MwwU8vXsFzHlUNGjUZjk6wQMhASEREReQ8GQlKVhg4ZdSS4WbZxJBC6YsioOyqEOp0OJSUlDrW1F/As+0P63l6AJCIiIqLG4+fpHSByJ3cEQndOKmNrG40pIiJCsQ+14ZBRIiIiIu/HCiGpihACPj4+8td1tbVs50glr7agpqZrCNU4ZLSsrAwajQbfffedp3eFiIiIyCEMhKQqQgj4+voCaNg1hN44ZNTbriF0ZMhobbOMNkUXLlwAAHz66ace3hMiIiIixzAQkqqYVwgbMmTUmfUaa1IZk8nk9IQ37qTGCiERERFRU8NASKpiXiGsK3w1dFIZ89e09bW9NnUxD1COhE1PcfQaQukY3BEI8/PzUVVV1eiv423vBREREZE9DISkKs4MGbU1jNOZ207Ym/zFFmcqhOZDRl1VIRRC4Ouvv3ZpkHHkxvTuDoStWrXCX//610bbvkajabRtExERETUGBkJSFVcNGXX1NYSeHjL61Vdf4YEHHsCGDRvqvQ1L3lghBIBt27a55XWIiIiImgIGQlIVo9EIf39/AI13Y/qGDBl1NhC6alKZ4uJiAMClS5fqvQ1LjgRC81DrrkBYXV3tltchIiIiagoYCElVzAOhK64htLWNhlQIHdEYQ0alqqkrQ1ltQ0bNb/1hMpng7+/PSWXqsGbNGnTq1MnTu0FERETNDG9MT6piMplceg2hqwOhp4aM+vnd+K/AlaGstgqhr6+vXOGUAqG9ANmUNOZMry+++CKuXbvWaNsnIiIidWKFkFTFaDTC19cXGo3GJdcQ2gpwls815pBRVwVCKSS7o0Jofh1nc6sQOhPsncUJa4iIiKgxMBCSqkgVQh8fH4eHjNZ2Y3pbIayxh4w2xm0nGiMQ1nZj+uYaCBvzGBgIiYiIqDEwEJKqGI1G+Pj4wMfHxyU3pnfVkFFnKoSNcQ2hJ4aMAn8cAwOh43iPQyIiInIlBkILs2fPxt13342wsDDExsbioYcewqlTpxRthBCYMWMGEhISEBwcjH79+uHEiROKNnq9HhMmTEBMTAxCQ0ORmpqKCxcuKNoUFxcjPT0dWq0WWq0W6enpKCkpaexDVDWpQuiqIaO1VQhtfe+t1xC6e8io9HrNrULYmNcQShXC5tBPRERE5D0YCC3s2rULL7zwAvbt24dt27ahpqYGAwcOREVFhdxm7ty5WLBgAZYuXYqDBw8iPj4e999/P8rKyuQ2GRkZ2LhxI9avX4/du3ejvLwcKSkpipO5tLQ05OTkICsrC1lZWcjJyUF6erpbj1dtnKkQOjKpjK1tWK7nTNXPEfZuO+GKWUZdMbFLXcHFcsioEKLZBEJ3DBk1GAyN9hpERESkPpxl1EJWVpbi+08++QSxsbHIzs7GvffeCyEEFi1ahOnTp2P48OEAgFWrViEuLg7r1q3D2LFjodPp8NFHH2H16tUYMGAAgBtTxicmJmL79u0YNGgQcnNzkZWVhX379qF79+4AgBUrVqBnz544deoUOnTo4N4DV4n6XENY39tO2Pq+rkloGjJk1BVDCV1R4dJoNBBCuKxCKIRoMtfPuSMQNofZWImIiMh7sEJYB51OBwCIiooCAOTl5aGwsBADBw6U2wQGBqJv377Ys2cPACA7OxsGg0HRJiEhAcnJyXKbvXv3QqvVymEQAHr06AGtViu3sUWv16O0tFTxIMdJFUJnhozWNizTFUNGnQ2EjTFkVAoyjgYavV5vd1ld9zR0dpbRr7/+2qF9agxHjx7FjBkzHG7fmLOMSlghJCIiIldiIKyFEAKTJk1Cnz59kJycDAAoLCwEAMTFxSnaxsXFycsKCwsREBCAyMjIWtvExsZavWZsbKzcxpbZs2fL1xxqtVokJibW/wBVyLxC6OprCO09V1fgq28gtJxl1BWB0NHqU23t6qpkWd4LUgqEUji0dPXqVYf2qTHceeedePPNN3H9+nWH2nPIKBERETU1DIS1GD9+PI4ePYp//vOfVsssh7A5MqzNso2t9nVtZ9q0adDpdPIjPz+/rsMgM+bXENYVvmq7hrC2ZZaBsK4KobNBrjFmGZXWdTTQ1BYIHakQ2hoyam8daXsN1ZCqXXl5uUPtOGSUiIiImhoGQjsmTJiAL7/8Et999x1uueUW+fn4+HgAsKriFRUVyVXD+Ph4VFdXo7i4uNY2ly5dsnrdy5cvW1UfzQUGBiI8PFzxIMc1ZoXQ3nBBZwJhQ4aMNiTwODtktLYqVV2TytgbMmpvHVddP9iQwFzbEFlzrBASERFRU8NAaEEIgfHjx2PDhg3YsWMH2rZtq1jetm1bxMfHY9u2bfJz1dXV2LVrF3r16gUA6NatG/z9/RVtCgoKcPz4cblNz549odPpcODAAbnN/v37odPp5Dbkeg29htAyCDoyZLQpXUPoiiGjdc1Yaj5k1JFA2KJFC4f2qS4NCcyODhk1/1DA1ZU8BkIiIiJqDJxl1MILL7yAdevW4YsvvkBYWJhcCdRqtQgODoZGo0FGRgZmzZqF9u3bo3379pg1axZCQkKQlpYmtx01ahQyMzMRHR2NqKgoTJ48GZ07d5ZnHU1KSsLgwYMxevRoLFu2DAAwZswYpKSkcIbRRuTqCqF5gLEXEusKfA0ZMuqq2064csioMxVC6XVrC4RSeGyohlTvnK0QrlmzBmvWrHHp5DIcMkpERESNgYHQwvvvvw8A6Nevn+L5Tz75BE8//TQAYOrUqaiqqsLzzz+P4uJidO/eHVu3bkVYWJjcfuHChfDz88OIESNQVVWF/v37Y+XKlYqT27Vr12LixInybKSpqalYunRp4x6gyrn6GsK6AqHl9405ZNSds4w6EkrstampqUFAQICijfS9rdd3VUWsIYHQGyaVkbBCSERERK7EQGjBkRNyjUaDGTNm1DodfVBQEJYsWYIlS5bYbRMVFYU1a9bUZzepnhpaIZRO+J2pENb1OvW9D6HlLKPuvIawIRXCmpoaq4qgrQqhr68vjEajywKQO64hdMV9HO3hkFEiIiJqDLyGkFSlodcQSoHFVuixFSCdHTLqSKiTXtNbZxmtq43BYLCqCNoLhFJ7V3DnkNHGwCGjRERE1BgYCElVXFUhtPW9vZlHXX3biboCYW5uLrp06YKqqiqnt+mKSWXqCpfOVAgB1wVC87765ZdfnFq3qQ0ZtfzZICIiIrKHgZBUxZlrCG0FQikI1VYhtKwo1hUIzU/wHakQSvtgLxAuWrQIP/30E06dOlXntsz3E7gxY64j3BEIpYlnHN2nuphv29H7CkrqO2TUlQHRmSGj/v7+8gRWthw6dAh/+tOfHD4uIiIiar4YCElVTCaTw0NGpRPv2gKheTCyFQhramrqDITSNvz8/JweMmprOKq0b85cUyjtoysDoSOTylgGQvN1pP1vjAqh9PqOqm+F0JXX+0nvqyPXHZtMJnz33Xd2ly9evBi//fYbzp07Z3P5uXPnsHnz5vrtKBERETUpDISkKgaDAf7+/g4NGbUVCGsbMmqrMmY0Gm0OKzUnhSApFNWlMa4hdKRCaCsY2+JMhdDy2C37Dmicawidvdl9fa8hbIwJYD7++OMGbyM0NBQAUFFRYXP5pEmTMGzYMLcMgSUiIiLPYiAkVTEPhHVV0MyHZlo+Z2vIqK2KVk1NTZ1DQs1vveDICbjU3nI4amMHQvN9qy3oOFIh9PPzU2zTViCUttMYFUJb/bxs2TJcvXpV/t78varvkFFXDXcF/hhCa/l6H374odPBTep/e+/R0aNHAQDXrl1zci/dR6/XY8OGDZ7eDSIioiaPgZBUxVUVQltVMFtByGg02hwGaU5aHhoa6lDwcLRCWJ8ho7VV/iyPq65t2WtjMBjg5+cHjUYjt7F1H0JXB0LzbVseZ1lZGcaNG4f09HSb7es7ZNSVgdDW+/nFF19g9OjRWL16tVPbkgKhvf3TarUAgKKiIif30n3+/ve/45FHHnHqWlkiIiKyxkBIqlJdXY2AgIA6ryEUQljdcxBQVufM/wWUAcb8xut1zdzZkEBoq0JZn/vVOTI80/w47B2TI8NKpQqhRqORjzc4OFixH4B7A6EU+AoLC+XnzMNSfYeMujIQ2upPaXKcyspKp7YlzeBqbybaiIgIAMDly5ed2q47SWHVmdl0iYiIyBoDIamKoxVCe1U9yyBoL8CYB5y6gpT0XEhIiEMBwnwoq+XwVOCPQOhMGHEkfDkSCOsalimtK70H9gKh+f0b6zoOk8nk0NDG2vZN2g97w2JtVQjvu+8+dO3aVfGcZf+58hpCqY/uuOMO+TnLDwEcJQVCe5XP6OhoAEB+fr7Vsv/85z/Ytm2bU69HRERE3ouBkFTF0WsIbc0eav68rfv2SdurqalBUFCQ3M68ja2AID3XokULpyuE5tuzXNeZQOjIfQhrC4SlpaW4/fbbFcP3aguEUoVQCiSWgdC8z+sKVe+88w6io6PrbFdbhVDqK3uB0Nb7snPnTuTk5NjcjuT333+vdZ+cIW1b6ivAuWHB5qQho1VVVZg8eTI++OADxXKpwm1+TaXkr3/9KwYOHFiv13Wl+symS0RERNYYCElVHK0Q2psIxtEKoRQILYeM2gotDRkyaiu0eGLI6MGDB5Gbm6sIFvbCpXQNYWBgIMrKygA0LBDu3LkTAORt2WMr2EukvjN/3vy9cHTIqGW7v/zlLw6t5wipH8xfw5EKoa1gbl4hnD9/Pp577jnFcil8lpaWNmyn3YAzoRIRETUMAyGpihQI67qG0N69AR29htCZCqGzgdB8H2xVF6WTfWduvt7QIaO21rd3ol5VVYXg4GAEBwdDp9MBsA6Ejs5oCvwxQ2ldgbC2CqHUd+YVTvP3wtFJZSwrhLZ+xj744AP87W9/c2h7trZt/hp1BcLffvsNfn5+cmiWSDOW2rv+Ljc3F4AyEF68eNHufQs9QTpmR8M6ERER2cZASKri6JBRKYQEBAQ4HQjNh4xKt52orWrXkAqhFA7CwsLkdaUZIi9evFjntiy36cwso5s3b5avJbM1S2lNTQ3effddnDhxQrGd69ev1xkInakQ1hVuJLVVCG0Nr61PhdCyovbEE09YtXnuuecwe/Zsh7Znrrq6Gj4+PjYD4bhx4+TnzI/zrbfeAgCrm9RLfWor6B4/fhzHjh0DoDyem2++GW3atHF6vxsbAyFQUFDg1TPCEhGRd2MgJFWRZhl1dFIZy0Do7JBRqUIoBZ7aAmFISEi9h4yah0npOXs3HbdFWseZCuGwYcPka8nsVQinTp2Ke++9V7GdqqoqBAUF2QyEtu79WFcglMJ2XVU8RyqEwB8hS3ouKCjI4QrhlStXcPPNN9t9nYaorq5WBH/Adsg1X75q1Sqb+yH1qa2Kn/l1g8XFxS7Y88Z18uRJl2znypUr2Lx5s0u21VieeuopTJkyxer5O+64A3FxcR7YIyIiag4YCElVHL3thHRSHRgYWOukMraqNZYVwpqaGoSEhFi1t3yt8PBwhyaCsQyE0vV4loHQmSGj9QmE5qT9tnVNoxT6JLaGjEr9Y6tf6+oTV1QIzUOUtB3ziqujVSidToebbrrJaluuYDAY0KJFC7u3w5D60pHXlNZbsGCB1TLz4afefGN66X1/4YUXUFJS0uDtPfDAAxg2bFi9J6k5deqUS28zYsunn36KefPmWT1/5cqVRn1dIiJq3hgISVUqKysRHBxcZ4VQOqkODQ1VnCBa3p7AvHJkPjGLMxVC6R5ykZGRTl1DKAVCf39/BAUFyetKJ6XOVAildeo7y6jUX4cOHZKfk67ps+zn69evyxVCqRrVokULAH/0oXk/2OoznU6HTz75BEIIl1QIzU/kpSBtHtRre18+/vhj+evS0lJFIHT2/oD2SD9Hlvti/rVUzbPVD1J4kthqc+TIEQDKQGkeNMyPq76uXLmCLVu2ONz+mWeewahRo2wuMz8mVwSiX3/9FUD9hqDW1NSgY8eOmDx5coP3g4iIyN0YCElVqqqqEBISUuc1hNJJcVhYmCKQSM/bCoRSO51Oh6ioKAA3gobBYJArYHUFQoPBUGtQlfZB2v/q6mr4+/s3uEJoq8JnqbZAKPXD6dOnAdwYZildf2bez0IIRYVQOpGXwoZlIPTz87O5TzNnzsQzzzyDM2fOyIHwySefrPUYHR0yaisQ1hY2N2zYIH/dWBVC6fW1Wq3dCqEUwM1fU/q5M6/66fV6rFixwuo1hg8fDuCPDxKefPJJFBYWysulW1FI/v3vfzt9HGlpaUhJSanzZ1zyySefKAK3OfNjciQQCiHw2Wef2f3QQ7oVR6dOnawmKBJC4KOPPrJbAZQqlOYfiBARETUVDISkGgaDQa7W+fj41DpdvXRSHRsbqwhWloHQ1kyUxcXFaNu2LYAbAaGkpEQOCfYCoa+vrzwZTF0hoqqqCmFhYfLX7gqEtc0gallVCQkJsXnLgsrKStTU1ECr1cpVU61Wq6iomm+vRYsWNvdJWn716lU5GFy4cKGWI3R8yKhlIKxryKh5paq0tBSxsbHy97W9l5cvX651f82ZB0J7s59K+23+nBTG3377bfl5ywlmJFIgkgJh9+7dUVRUJP9sWFYZR4wYgZqaGtx77704fPiwQ8chvUeWw4jrw7wfHOnLnTt3YuTIkfjwww9tLpdmq/3tt9+wd+9exbIff/wRzz77LP7xj3/YXFeqzkoz/Nbl119/xR133OGSfjDnyiHK1PxUVFRgx44dnt4NIvJCDISkGtLJUnBwMFq0aFFrYJLaxsXFKYJNbRXCqqoqVFVVQa/X409/+hOAG4Hl6tWraNmypVV7SWVlJUJCQhAeHg6g7tsnVFVVITo6GgBQWFiIsLAwRSCUXqM+k8rUNmTU/GTTvN2FCxdsBkJbJ7vSiXNkZCRCQ0MB3KgOSifSjgZCqfJVXFwsb2fo0KG1HaKiumOvwgn80W91VQilAP/VV1/J+63T6RAdHY0LFy5g4sSJtZ6gP/LII7XurzlpO+Hh4TAajTY/kJB+nqW2Wq1W8fpStc9eaDEPhD4+PkhKSgJwIyABtoe/Xr58GT/88ANeeeUVh45Dem0pwGVlZUGj0dRr8hrz98SRQCi9r/baSscPWN/s3vwDCFukay0dDYSrVq3CsWPH8OOPP9pcvnHjRrzwwgvy9zU1NTZvLbJnzx7Fz3XLli2duv8oqcuUKVPQv39/xe11iIgABkJSEfNAGB0dXeswM3uBUDoJtRUIr1+/Lp/YxsXFoUWLFrh27RquXbuGli1bIjg42ObkF1IglKp+jgRCqeL4+++/Q6vVIjAwUN4Xad9dXSE0D5jm7UpLS60CYVRUFAoKCqzamwdCKTSfP39ePhm3nCnTXiCUqotXr161uv+iPeb7aBkIzfuqvLwc3333HebOnQsAiIiIsBnszLdnft++8PBw3HzzzQgJCak1EDo6bBJQVggB2zept6wQRkZGKrYh/VzZ+7mQfqYrKysRGhoq32IiPz9fXs/yOkLpd+PcuXMOXS8p7XdOTg4A4J///KfiNczVNbmLs4FQam+v36UKIWDdR9K69j4wkX6ud+3a5dCsp4GBgQDsV/SGDx+O9957D6+88gpMJpPid0la7+TJk+jduzfefvtt+XmdTqcY5ttYhBB4+OGH5Z97ahwLFy7Erl277C6vqKjA+PHjYTKZcO7cOeTl5dW6vePHjwMAvv/+e4def8eOHdBoNI12/9FBgwbh9ddfb9A2rl+/3uiTOdlTXl6OHj161NnvRE0BAyE1CcXFxXj11Vcb9B+/dJ+umJgYxMTEOBUIpZNTW0NGpaF05oEwMjISUVFRcoUwOjoaUVFRNmdtdDYQVlZWWgXCiIgI+bWlE/P6TCpjNBrtnohL2wsMDFRsW6fTWQVCy/6V+l66v11kZCRuueUW+bW1Wi0CAgLkE9+6AqFk2bJl8sl6Xf1mvo+WQ14tA+Ff/vIXHDhwAAAQHR1t85oy8+1VVlbCYDCgqqpKrvQGBwdbnfCbH0u7du1q3V9z5lU/82PR6/W44447FMcgtZWuY5VIxyCFOPMKFHAjmAM3rofTarXyz2N5eTmMRiP0er18bBLpA47Tp09j0KBB8vPV1dVYvny5VfiSfm4effRRAMCXX34JwHblzTxYSe3NORsIpeO39/NtXiG0HO4sHefatWtthjjz3+ukpKRaZz3NycnBq6++CgB13jtwzpw5+PTTT62O7/Lly1i8eDEA4O9//7tiWW3bLCoqcsmtUM6ePYtNmzbh9ttvr9f6Op2uzg9EfvnlF/z88892l3vz/SeLiorkDz3qq6KiApMmTUK/fv1w/fp1mx+4PPbYY/jf//1fJCUloU2bNmjXrl2tl0JI21i3bh0OHz5ssw8rKyuxbt061NTUoH///gDQKPcf3bVrF7Zu3Yq3334b3333HTQaDVavXu3UNubPn4/g4GD5AxZ3GzNmDPbv34927dpBCIGdO3dCo9HI12h/8803GD16tMteb8KECVa3caqvf/7zn/jqq69qbaPT6XDu3Dmv/l0j12EgpCZh7ty5mDlzZoP+yEqfcrZu3RoxMTG1nkSaB8Kamhr55FM68auoqIAQQnHfOctA2LZtW5w5cwbXrl1DVFSUHBAtSYFQqujYG5YG3DiZraiokO85duHCBWi1Wtx0003y8dRVCbLFPKjU1NSgtLQUly5dUrSRTubj4+MVw0FtBULz6+iAP07YpYpQXFwcxowZg4cffhi7du2Cj48PWrdubTU8MSwszOaHANL7c/HiRfk9sXXNornaKoRlZWWIj48HYN1vERERVs+VlJQogsXZs2fl15dCm61AaD400pnrx6RtSyFP+nksLy9HeHg4goKC5PdHWmYvEF6+fBlhYWGYM2eOvKxNmzbQ6/VYtmwZrly5gpiYGEUglLb9t7/9Da1bt5bXMz+h2L17t/z1mjVrMHbsWKsTDvN+LykpkYOTZTu9Xq8IG//6179gNBoRGhqKN954A4ByCGtdgfDy5ct45plnAABffPGFzTBivm+WP0vS+1ZQUIC33nrLal3LIa9Hjx61uy/SBw0AMH78eJsn+ubhtKKiwur/hNatW+ODDz6wuX3L31uJTqdDXFwc0tPTYTQasXjx4npfc+jIh00ZGRnYv3+/1fNGoxERERF47rnn0LlzZ+Tm5uL333+3CrK33XYbOnXqZHPbJSUlCAoKwpo1a+q1/7ZUV1cjPDzcqVlw7YmLi0PXrl2dWsdkMmHFihV49NFHcerUKZw5c0ZeFhwcLA+NN2f+/6DkX//6F77++musWLECGo1GMemV9Huyc+dOdOvWDUFBQZgzZ46iqvzKK6/g8ccfV1TMAWDfvn0YNWqUUyMbavP555/LX//lL38BcGMiK0dGBnz99dcQQihm9XX19biWLl68aDWRlvn3TzzxBO677z4AN4Li1atXMXjwYHz44YfIzMysddtXrlzBpEmTav3w8+rVq1i6dCl++OEHeUZk4MalAGVlZU69L6tWrUJaWhpSU1Ot/haWlJTIQ4ojIiLQpk0bBAUF1fpBAzUPDITUJEjVJmeqXpbOnTsHf39/tGzZEq1bt8bFixftfpIvTdYiVeKkE7LCwkK0bdsW1dXVqKysREFBgVzpKS8vl7cnBcIvvvgC1dXV0Gg0dgNhYWEhoqOj5esMLYeHmbt27RoMBgOSk5PldSMiIhSBUPrD6EwgND8pra6uxt133y0HJEl5eTn8/PwQFRWF33//XX6+pKREEbaioqIUN2cH/nj/br75ZiQlJSEiIgJ+fn7YsGGD/Ilnu3bt8MsvvwC4UXHy9/dHu3bt7N52ArjxB8vyuO0pLi6Gj48PIiIirE74y8vLER0dDX9/f5SWlqJVq1bysrCwMKsKoXTSLQXfRx99VH59qYoWGRmJ0tJSxf5LFY/bbrtNsQ/ffPON4gTJkvRhRvfu3QEA7733HkwmE4qLi+WJZqRqm3RiLX2qn5iYCOCPQHjp0iXExcXJ12ECwGeffQYAGDduHC5duoTo6GgEBATAz88PZWVlcuC5+eabFfs5c+ZM+Wvzk3cpxJv/vAshFMFJ+nAAuDE0rlWrVpg1axbWrl2Lnj17WvXBrbfeisrKSjmQXbt2Db1790bLli0VJ73S67/33nty3589e1ZelpOTY7MSYX6SZfnzYf47WVeFELhRbb1+/TrGjx9v9Ttv/nNqMpkQGhqKO++8U37Oskp/5swZeRv33HOP1WtbunTpEvbs2YNTp04p+lu6XnH9+vUYN24cMjIyEBISUufsvBLzfTKfwEmj0WDBggX4xz/+ge+//x5GoxHXrl3D4sWLkZ6ebrUd6UOf5cuX4/jx4xg5ciRuueUWJCQkoH379njxxRcVfTZx4kRFdXfZsmXykEdb27fn4sWLWLZsGd577z3897//tVp+7NgxlJWVYerUqXjmmWdqDctXrlyR96miosJqEiLJyZMnMWvWLIdO1nfv3o0xY8bgX//6Fzp27IguXbpYtbHcJ+mDQfP/69PS0vDAAw9gzJgxAG5cq/zGG2/g559/lkcBmHvllVfQsmVL/P7777h06RKWLFlitV8A0LNnT3z88cfw9fXFvn37bB7DtWvX0LVrV6shlG+++SZycnJQUVGBixcvwmAwWL2OZNGiRVi6dCk0Gg00Gg1++uknZGdny9tcuHAhHnjgAXlkjvR3c/jw4fI6toagW+6n9PN89uxZVFdX48KFCzhy5Ai++eYbADeCp/TBL3Dj/74RI0bI/zfMnTtXnlQLuFF1NRcTEyN/vWDBAnnfNBoN+vbti5qaGnnm4qSkJCxcuBABAQEoKiqSP9Q7deoUlixZAiEEnn76aXl7t956K0pLS2E0GtGyZUuEh4fD19e3zvvG/vTTT+jatatiW/7+/tizZ498qUBkZCQ6duxodc1yfe/PSk2IoCZNp9MJAEKn03l6VxrVyJEjBQCxYcOGem9jypQpol27dkIIIfbs2SMAiP79+9ts+8orr4jExERx9OhRAUD8+OOPwmQyicDAQDF06FABQJw7d0706dNHPPHEEyIqKkq89dZbYtKkSQKAuH79upg5c6YAIACIqVOninHjxon27dtbvVaHDh3Ec889J4QQcnt7jh07JgCI3bt3y21feuklsXLlSgFAlJaWisTERBEQECCioqIc7pvk5GRx2223CQAiLy/P5n68+uqrIiEhQdx7772iR48ecpv33ntPjBs3Tv4egPjPf/6j+F46voCAAPHAAw/Y3IcZM2YIAGL16tUCgIiKihKjRo0SAITJZFK0feSRRwQAccstt4ioqCjRokULcdNNN9V6jDNnzhQ33XST6NSpk5gwYYJi2TPPPCN69OghEhISxGuvvSYAiH79+okdO3aIVatWCQCivLxcbr9z504BQJw8eVI+xsOHDwsA4uDBg0IIIbZs2SIAiPPnz8vrSW0ffvhh0aVLF6vnS0tLhclkEtXV1Yr9+/vf/y6ioqLElStX5Lbvv/++6Natm3j22WdFZGSkACAMBoOYOHGiaN++vRgzZowAIOLi4gQA8fHHH4vdu3eLwMBA0bdvX6vXNX+/Ro4cKS/v37+/2L59uwAgcnNzhRBCvP3224r20uPcuXNCCCHee+89+bmKigohhJBfIz4+3ua69h6333671XPl5eXyz1VSUpLVz8iaNWsEADF79mwhhBAbNmxQrD9nzhxF/+7bt8/qNaRjFUKIxx9/XPH7ZjKZxIwZM0R2draoqakRL730ktX68+bNEwBETEyM+OWXX4QQQkyfPl0AEJ06dRIJCQmK9pLz588LAIrlixcvFgEBAcJkMonu3bvLz+/cuVOMGDFCsZ2uXbvKX997773ydpcuXWq3j7Ozs61+x4qKiuSvjUajiIuLE/PnzxelpaXi/ffft7utTz/9VOzYsUP+PiYmRt5ObeuZPyZMmKD4vl+/fkIIIQoKCqzabtmyRZSWloq63HHHHTb7WzJgwADF8qVLl9rdlvk2/va3vwkAYtmyZfJ7YrmPLVu2FDU1NbXu36BBg2z2hZ+fn2jVqpX8/Y4dO4QQQpSVlTn1e1Sfh9QHtpb98MMP4uOPPxaXLl2Sj8HW77b0/4+th4+Pj/z1+PHj69yfN954w6HnpMf8+fPlr0eNGiX++c9/yt8nJiaKwMDABvfRiRMnxOeffy5/X11drVi+YsWKRn+f6nrYel+eeeYZh9Z96623xDPPPFPn71djUcs5rjdgIGzi1PLLkpKSIgCITz75pN7buPfee8V9990nhBDCZDIp/gO3NGzYMDF48GBRUVEhAgMDxcKFC0VRUZHiD9DatWsFAJGZmSmSkpIUIUkIIb799lv5e51OJ5+ULl++XH6dH3/8UQAQc+fOFUL88Yd38+bNNo9h8+bNArhx4i213bRpkzhy5IgAIJ588kkBQHTo0EEEBAQ41C81NTUiMDBQPPvsswJQhjkp3AghxP/8z/+I++67T6Snpyv+YDz++OPiwQcfVDx36dIlu39gnn32WZv7UVhYaNVWOsmSToIkffr0UbS75557hK+vr7h+/brd4xw9erTo3LmzGDBggEhMTFQs8/PzE126dFGcTE+dOlUIIcSBAwcEAHHgwAG5/fr16wUAUVxcLPz8/AQAsWjRIgFAXLhwQQghxNmzZ+VjMH8d4EZQASCfyJofy/Lly+X3QSItM/+5lR6vv/66+Ne//iUAiHXr1gkA4q9//ascCGNiYqzWkUL5n//8ZxERESGEEOKDDz6Ql/fq1ctqvwAIo9EohBBiyZIlNt/bFi1aCKPRqAhIWVlZQgghTp06pTg+6WErjJk/PvvsM6vnbr31VgFATJo0Sbz44ouKfRbij+AFQHz11Vfy9/3795ef37Ztm3w80j5Jv0cAxF133SVKSkrEzp07Rffu3UXnzp0FANG3b1/x22+/ye1Gjx4tf52TkyN/Lb3XAERkZKQwGAzy9wEBAYrQJD2qqqrE3LlzBXAjkFoGbyGU/3dVVVXJx7x3716b/Xfo0CEhhBBPPPGESEpKkj9ksXxs27ZN3lZ2drYAIJ544gmxbt06+YMxRx6jRo2yCgEFBQU2f54cfdx2222K/bL1eP3118WxY8fE5cuXhRBCGAwGsWvXLiGEkD9AMH/07t1bDB8+3G5QTk5OFjfddJPo3r27KCoqErm5ufKHRdLDPKg5+tixY4cwGAzCaDQKnU4ndu3apdiHTp06yV/r9Xr5PXn66adtbm/BggUCuPH3qF+/foplV65cUXygEh0dLf8dy87Otgre0kP6P0xSXl4u5s6dKw4dOlTv99DWo6SkRNTU1MgfRjj7YZH0OzBx4kSX7pejjz//+c9yH5WWlsofGkp/nw8dOiRMJpP8/1VmZqb84Zqth/kHT7YeR48etfr/E7jxu+vMfnft2lUIIcTXX39tc/m1a9fE+fPnRX5+vvA0tZzjegMGQi/wv//7v6JNmzYiMDBQ/PnPfxbff/+9w+t68pfF/KTL19dXLF68WLz22mti+/btIi8vT/6jYjQaRWVlZYNeS3qdhQsX1mt9qer15JNPys+1adNG3u7ly5flk0PpU+hp06YJIYQIDw9X/Ge5detWxffTp08XDz30kOI5yZEjR4TBYBBCCKsqjPlj9+7dQgihONGU/qOX1jcajfLzJpNJ/Oc//5GPp6KiQrGe9Ideqs7U5uDBgwKA+Pe//221X2PHjhXnzp2TT1JHjhxp949vZmamuP/+++WTv4yMDLFr1y6r8Cb1sy2LFy9W9IllsGzRooX8dZcuXeSvp02bJn995swZcfjwYVFUVCR27twpXnjhBbmimJmZKZ9YJSUlifHjx8snIf/zP/+jCLvSz69l35o/TCaTHMKAGyHAXGxsrNU6q1evVvT1Y489Znf7UnA0/7k6ffq0iIqKkp+bM2eOVZju37+/MBqNYv78+fIJoPkjLy/PZv8HBATIvw9CCDF48GCbP9fmvwNGo9FmVUR6xMXFiXvuuUfRr9LXHTp0EEIIcfz4cav1WrVqJfr37y9KS0vl971jx46KNkePHhWVlZV2X9v8IYUKy+fT0tIUx2dv/bffflt89NFHtb6GEEKuDNb2kD4UOnPmjHyyaPmQTpRt9f/GjRvFZ599ZvX+1VXxGDNmjBBCiFtuucWhPqvrYTkSwJlHfHy8OHPmjDh9+rSoqakRgLLaM3v2bHH27FmRmppqdxvm/4e74tGyZUuXbUsatXH16lWH19myZYvd/xuFEDZ//upiMpmE0WgUhYWFdtu89957YtGiReLq1atWYdBSTU2N/P+Ercfvv/8ufy39bUlISBARERHi6aefFocPHxa7du1y6JwlNzdX/lmRPuACbgTlI0eOyFVsk8kkNm/eLGpqaoTBYBAFBQXyqAnLR//+/RX/X129elVs375dZGZmip9++kkYjUaxbNkyMW/ePFFcXCwmTJggnnjiCZGTkyOqq6vF888/LwYOHCh/0FIfRqNR8XfQvDpfUlIiLl68KM6ePSuEEKKqqsrqQ2uDwSDmzJkj/v3vfyuel/4ujB49Wq44Dx06VMTExIi2bduKkydP2qxUnzx5UkyfPl18++239T6mxsJA6D4MhB62fv164e/vL1asWCF+/vln8eKLL4rQ0FB56FVdPPnLUtfJkSse0jBGew/p08/27ds7tL0jR45Y9Z29h/QeWFZDrl+/rhgiWV1drfgEvbYqlWVwlB7mfxDM//DZegwePNjmts1P/M0rFY4+SkpKFN+HhoZatTl79qzIz88XAMTEiRPFXXfdJS/76aefbO6XyWQSTzzxhFi8eLFDQ7sMBoMckIUQdk8If/jhB/lr86GUtT1OnjwpLl++bHPZb7/9Jn799VcRFBQkPv30U8U+Pfroo1btx44dK4QQiiFClsPMpOqi+ePSpUuKapGjj1OnTsnbNRqN4u677xbPP/+8/Km0eVvLT3Z/+eUXeZn50Fdb75X0AYT0XrRu3drmegsXLhTZ2dny99IJsPTo2LGj1e9vfHy8EEKIl156SUyaNEnxu1JZWSm++eabWj8wEEIofral35v9+/db9VdERITi+7vuuksIoRzOavmQSD/j5g+9Xi+Ki4utjkf6+vjx4/L63333nQBuhEjzCiIAuWplTqpySo/Ro0fLy/R6vcjPz7ca0mlLTU2N6NChg9i5c6fNIYXSSWZNTY0oKysTq1evFps2bbJbNbR8SB+OPfjgg/JrlpaWipqaGnH48GGxcOFCuW1ubq7NDyPOnj1b6/+RRqNRcdJqOUxUqvJKAUev1zv8O3Tp0iVx7do14efnJzp27Cj/3Xj11VfFxo0bhdFoFCaTSZSVlcmjHhITExXbePbZZ8WHH34o5s6dKyIjI8Xq1auFXq8XK1asEEePHhUpKSlyP0uys7PF/fffL490sfUwH4FQmz179gh/f3/RqVOnOoehNrbS0lKxY8cOce7cOZGXl6f4v4PIFRgI3YeB0MPuueceMW7cOMVzHTt2FK+88opD63v6l6WgoEDs3btXfPHFF+K2224TSUlJNq/58YbH+PHjrfbf3qe3U6ZMsTpOvV6vOJG5cOFCvf4gFxUVibVr14p58+aJkpISm21+/fVXu8dha4irOenEUbru0pHHmjVrhBA3ToQ1Go1YsWKFfJ2l9Pjiiy9svt6VK1ccqkQ2RFVVlZg9e7Z8TB9//LEQ4kaIkKrPdX0Sb15RuXjxokhLSxM+Pj7iySefrPXTc0leXp7o06eP6Nu3r9izZ4/D+15aWip++uknMWvWLPH555/LzxcWFgpfX18BQPz3v/8VlZWVIioqSnTv3l0UFhaKv/zlL2LKlCnis88+Uwwds8doNIpr167ZPSnT6/WN/j5dunRJzJw5U+h0OmEymcT169fF3/72N/H000+LnJycRn3t77//XgA3hvquXLlSVFdXi4qKCvHFF1+In3/+2SqEmEwm+bqxdevWiWvXrimWl5eXi4sXL4qqqipFGNuxY4d47bXX5OcuXLigGLppi8FgEPn5+bV+0Gc0GsXZs2dFWVmZs4deK51OJ44dO1bn/xtC3BjWu2DBAnHlyhVRU1MjLl++LM6dOyeefvppUVxcXO99+Pe//y3efvvtWj+MqE1VVZXIz8+vcx+uX78uzp07J86ePSt+++23er9eY9Lr9Q6FeyLy/DmummiE4NRBnlJdXY2QkBD8+9//xsMPPyw//+KLLyInJ8fmDWn1er1iRsfS0lIkJiZCp9NZ3SOs0X32GbB2rcPNTULcmLlKCBQUFKBFWBgM1dUICg5GdXU1TEYjysrLERYWhkuFhYiPj0fhpUu4KSYGN910E0wmE3Q6HSoqKnBTbCyuXrkCXWkp/Hx9cfnyZbT7058QHBQEHx8fhISGosZguHGzco0Gvr6+0NS9i16pxmhEjcGAwKCgeh9DSUkJBAD//5vOXqPRQOPjg8DAQPj6cLJhIiIissHHB9i0ySMvXVpaCq1W65lzXJXxq7sJNZYrV67AaDTKU0dL4uLiFPcFMjd79my8+eab7tg9l/ORpjHWaJCQkKBYFvp/U+BHR0ff+Pf/7qFmfi81Hx8fREZGyvfrS0hIkLfTvn17q9fz9dDNal3Nz9cXfr6+DdpGRESEa3aGiIiIiJoVBkIvYOt+L5bPSaZNm4ZJkybJ30sVQo949NEbDyIiIiIiapIYCD0oJiYGvr6+VtXAoqIiq6qhJDAwEIHNpPJFRERERESexYuHPCggIADdunXDtm3bFM9v27YNvXr18tBeERERERGRWrBC6GGTJk1Ceno67rrrLvTs2RPLly/H+fPnMW7cOE/vGhERERERNXMMhB726KOP4urVq3jrrbdQUFCA5ORk/Pe//0Xr1q09vWtERERERNTM8bYTTRyn5CUiIiKi5obnuO7DawiJiIiIiIhUioGQiIiIiIhIpRgIiYiIiIiIVIqBkIiIiIiISKUYCImIiIiIiFSKgZCIiIiIiEileB/CJk66a0hpaamH94SIiIiIyDWkc1veIa/xMRA2cWVlZQCAxMRED+8JEREREZFrlZWVQavVeno3mjXemL6JM5lMuHjxIsLCwqDRaDy9O16ptLQUiYmJyM/P541NXYR96nrsU9dif7oe+9S12J+uxz51LU/3pxACZWVlSEhIgI8Pr3JrTKwQNnE+Pj645ZZbPL0bTUJ4eDj/QLgY+9T12Keuxf50Pfapa7E/XY996lqe7E9WBt2DcZuIiIiIiEilGAiJiIiIiIhUioGQmr3AwEC88cYbCAwM9PSuNBvsU9djn7oW+9P12Keuxf50Pfapa7E/1YOTyhAREREREakUK4REREREREQqxUBIRERERESkUgyEREREREREKsVASEREREREpFIMhNRszJ49G3fffTfCwsIQGxuLhx56CKdOnVK0EUJgxowZSEhIQHBwMPr164cTJ054aI+bltmzZ0Oj0SAjI0N+jv3pvN9//x1PPPEEoqOjERISgi5duiA7O1tezj51XE1NDV599VW0bdsWwcHBaNeuHd566y2YTCa5Dfuzdt9//z2GDRuGhIQEaDQabNq0SbHckf7T6/WYMGECYmJiEBoaitTUVFy4cMGNR+E9autPg8GAl19+GZ07d0ZoaCgSEhLw5JNP4uLFi4ptsD+V6voZNTd27FhoNBosWrRI8Tz7VMmRPs3NzUVqaiq0Wi3CwsLQo0cPnD9/Xl7OPm1eGAip2di1axdeeOEF7Nu3D9u2bUNNTQ0GDhyIiooKuc3cuXOxYMECLF26FAcPHkR8fDzuv/9+lJWVeXDPvd/BgwexfPly3HHHHYrn2Z/OKS4uRu/eveHv74+vv/4aP//8M+bPn4+IiAi5DfvUcXPmzMEHH3yApUuXIjc3F3PnzsW7776LJUuWyG3Yn7WrqKjAnXfeiaVLl9pc7kj/ZWRkYOPGjVi/fj12796N8vJypKSkwGg0uuswvEZt/VlZWYnDhw/jtddew+HDh7FhwwacPn0aqampinbsT6W6fkYlmzZtwv79+5GQkGC1jH2qVFef/vrrr+jTpw86duyInTt34qeffsJrr72GoKAguQ37tJkRRM1UUVGRACB27dolhBDCZDKJ+Ph48c4778htrl+/LrRarfjggw88tZter6ysTLRv315s27ZN9O3bV7z44otCCPZnfbz88suiT58+dpezT50zdOhQ8cwzzyieGz58uHjiiSeEEOxPZwEQGzdulL93pP9KSkqEv7+/WL9+vdzm999/Fz4+PiIrK8tt++6NLPvTlgMHDggA4ty5c0II9mdd7PXphQsXxM033yyOHz8uWrduLRYuXCgvY5/WzlafPvroo/L/o7awT5sfVgip2dLpdACAqKgoAEBeXh4KCwsxcOBAuU1gYCD69u2LPXv2eGQfm4IXXngBQ4cOxYABAxTPsz+d9+WXX+Kuu+7CX//6V8TGxqJr165YsWKFvJx96pw+ffrg22+/xenTpwEAP/30E3bv3o0HHngAAPuzoRzpv+zsbBgMBkWbhIQEJCcns48doNPpoNFo5FEC7E/nmUwmpKenY8qUKejUqZPVcvapc0wmE7Zs2YLbbrsNgwYNQmxsLLp3764YVso+bX4YCKlZEkJg0qRJ6NOnD5KTkwEAhYWFAIC4uDhF27i4OHkZKa1fvx6HDx/G7NmzrZaxP53322+/4f3330f79u3xzTffYNy4cZg4cSI+/fRTAOxTZ7388st47LHH0LFjR/j7+6Nr167IyMjAY489BoD92VCO9F9hYSECAgIQGRlptw3Zdv36dbzyyitIS0tDeHg4APZnfcyZMwd+fn6YOHGizeXsU+cUFRWhvLwc77zzDgYPHoytW7fi4YcfxvDhw7Fr1y4A7NPmyM/TO0DUGMaPH4+jR49i9+7dVss0Go3ieyGE1XME5Ofn48UXX8TWrVsV1w1YYn86zmQy4a677sKsWbMAAF27dsWJEyfw/vvv48knn5TbsU8d89lnn2HNmjVYt24dOnXqhJycHGRkZCAhIQFPPfWU3I792TD16T/2ce0MBgNGjhwJk8mE9957r8727E/bsrOzsXjxYhw+fNjp/mGf2iZNyvXggw/ipZdeAgB06dIFe/bswQcffIC+ffvaXZd92nSxQkjNzoQJE/Dll1/iu+++wy233CI/Hx8fDwBWn14VFRVZfQJON/7QFhUVoVu3bvDz84Ofnx927dqFf/zjH/Dz85P7jP3puJYtW+L2229XPJeUlCTP3MafUedMmTIFr7zyCkaOHInOnTsjPT0dL730klzRZn82jCP9Fx8fj+rqahQXF9ttQ0oGgwEjRoxAXl4etm3bJlcHAfans3744QcUFRWhVatW8t+pc+fOITMzE23atAHAPnVWTEwM/Pz86vxbxT5tXhgIqdkQQmD8+PHYsGEDduzYgbZt2yqWt23bFvHx8di2bZv8XHV1NXbt2oVevXq5e3e9Xv/+/XHs2DHk5OTIj7vuuguPP/44cnJy0K5dO/ank3r37m11K5TTp0+jdevWAPgz6qzKykr4+Cj/jPn6+sqfcLM/G8aR/uvWrRv8/f0VbQoKCnD8+HH2sQ1SGPzll1+wfft2REdHK5azP52Tnp6Oo0ePKv5OJSQkYMqUKfjmm28AsE+dFRAQgLvvvrvWv1Xs02bIQ5PZELncc889J7Rardi5c6coKCiQH5WVlXKbd955R2i1WrFhwwZx7Ngx8dhjj4mWLVuK0tJSD+5502E+y6gQ7E9nHThwQPj5+YmZM2eKX375Raxdu1aEhISINWvWyG3Yp4576qmnxM033yw2b94s8vLyxIYNG0RMTIyYOnWq3Ib9WbuysjJx5MgRceTIEQFALFiwQBw5ckSe9dKR/hs3bpy45ZZbxPbt28Xhw4fFX/7yF3HnnXeKmpoaTx2Wx9TWnwaDQaSmpopbbrlF5OTkKP5O6fV6eRvsT6W6fkYtWc4yKgT71FJdfbphwwbh7+8vli9fLn755RexZMkS4evrK3744Qd5G+zT5oWBkJoNADYfn3zyidzGZDKJN954Q8THx4vAwEBx7733imPHjnlup5sYy0DI/nTeV199JZKTk0VgYKDo2LGjWL58uWI5+9RxpaWl4sUXXxStWrUSQUFBol27dmL69OmKk2v2Z+2+++47m/9vPvXUU0IIx/qvqqpKjB8/XkRFRYng4GCRkpIizp8/74Gj8bza+jMvL8/u36nvvvtO3gb7U6mun1FLtgIh+1TJkT796KOPxK233iqCgoLEnXfeKTZt2qTYBvu0edEIIUTj1iCJiIiIiIjIG/EaQiIiIiIiIpViICQiIiIiIlIpBkIiIiIiIiKVYiAkIiIiIiJSKQZCIiIiIiIilWIgJCIiIiIiUikGQiIiIiIiIpViICQiIiIiIlIpBkIiIiIiIiKVYiAkIiIiIiJSKQZCIiIiIiIilWIgJCIiIiIiUikGQiIiIiIiIpViICQiIiIiIlIpBkIiIiIiIiKVYiAkIiIiIiJSKQZCIiIiIiIilWIgJCIiIiIiUikGQiIiIiIiIpViICQiIiIiIlIpBkIiIiIiIiKVYiAkIiIiIiJSqf8fctMnLai9xXUAAAAASUVORK5CYII=", "text/html": [ "\n", "
3Dmol.js failed to load for some reason. Please check your browser console for error messages.
3Dmol.js failed to load for some reason. Please check your browser console for error messages.
3Dmol.js failed to load for some reason. Please check your browser console for error messages.
3Dmol.js failed to load for some reason. Please check your browser console for error messages.