{ "cells": [ { "cell_type": "markdown", "id": "8482721f-0615-4b73-83ee-187b2f890521", "metadata": {}, "source": [ "# SAFT-VR-Mie with polar contributions" ] }, { "cell_type": "code", "execution_count": 1, "id": "d1842386", "metadata": { "execution": { "iopub.execute_input": "2025-10-15T23:11:58.922823Z", "iopub.status.busy": "2025-10-15T23:11:58.922654Z", "iopub.status.idle": "2025-10-15T23:11:58.938044Z", "shell.execute_reply": "2025-10-15T23:11:58.937404Z" } }, "outputs": [ { "data": { "text/plain": [ "'0.23.1'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import teqp\n", "teqp.__version__" ] }, { "cell_type": "code", "execution_count": 2, "id": "16f05ec8", "metadata": { "execution": { "iopub.execute_input": "2025-10-15T23:11:58.939256Z", "iopub.status.busy": "2025-10-15T23:11:58.939128Z", "iopub.status.idle": "2025-10-15T23:11:59.199231Z", "shell.execute_reply": "2025-10-15T23:11:59.198743Z" } }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import math" ] }, { "cell_type": "code", "execution_count": 3, "id": "be8268a7", "metadata": { "execution": { "iopub.execute_input": "2025-10-15T23:11:59.201280Z", "iopub.status.busy": "2025-10-15T23:11:59.201009Z", "iopub.status.idle": "2025-10-15T23:11:59.688150Z", "shell.execute_reply": "2025-10-15T23:11:59.687730Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkcAAAG4CAYAAABPb0OmAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjcsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvTLEjVAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAlixJREFUeJzs3Xd8lfX9///H2TnZew8SwhIIe6+wBAeK1lk/ro/ab+v+WT9ttZ+PYlvF0WFbFK1txVErWsVWQGWGvSHsgGFlL8g8Sc68fn9cyUkOSSBAkpOE1/12u27nnOu6znXe50BynnlPjaIoCkIIIYQQAgCttwsghBBCCNGdSDgSQgghhGhGwpEQQgghRDMSjoQQQgghmpFwJIQQQgjRjIQjIYQQQohmJBwJIYQQQjQj4UgIIYQQohm9twvQE7lcLgoKCggICECj0Xi7OEIIIYRoB0VRqK6uJjY2Fq227fohCUeXoaCggISEBG8XQwghhBCXITc3l/j4+DaPSzi6DAEBAYD64QYGBnq5NEIIIYRoj6qqKhISEtzf422RcHQZGpvSAgMDJRwJIYQQPczFusRIh2whhBBCiGYkHAkhhBBCNCPNap3E6XRit9u9XQwhrloGgwGdTuftYgghLpGiKLgqK9EFB3utDBKOOpiiKBQVFVFRUeHtoghx1QsODiY6Olqm3BCiB1AUhZoNGyhbvBicLvp8/pnXfnYlHHWwxmAUGRmJr6+v/FIWwgsURaG2tpaSkhIAYmJivFwiIURbFJeLmnXrKHt7MfVHjgCgMZmwnTqFKSXFK2WScNSBnE6nOxiFhYV5uzhCXNXMZjMAJSUlREZGShObEN2M4nJRvWoVZYvfwXrsGAAas5mQu+8m7MEH0EdEeK1sEo46UGMfI19fXy+XRAgBTT+LdrtdwpEQ3YTidFK18hvK3nkH24kTAGj9/Ai55x5CH7gffWiol0so4ahTSFOaEN2D/CwK0X0odjuVy1dw9p13sJ05A4A2IIDQe+8l9L57vdoB+3wSjoQQQgjRaVw2G5XLvuLse+9hz8sDQBcUROiDDxByzz3oLjJbtTdIOBJCCCFEh3PV11Px+b84+7e/4SgqAkAXGqqGort/iM7fz8slbJuEI3HV69OnD08//TRPP/20t4sihBA9nstiofzTpZx9/32cZWUA6CMjCXv4IYJvvx1tw2CJ7kxmyBYeioqKeOqpp0hNTcXHx4eoqCgmTZrE4sWLqa2t7fLyzJs3j7lz57Z6bNOmTWg0Gg4cONDFpRJCCHE+Z3U1ZYsXkz1jJiVvvIGzrAxDbCzRC16k7+pVhN53X48IRiA1R6KZkydPMmnSJIKDg3nllVcYOnQoJpOJgwcP8pe//IW4uDhuuummFs+z2+0YDIZOKdNDDz3ED37wA/Ly8oiPj/c49v777zN69GjS0tJaPE9RFJxOJ3q9/BcXQojO5Cgv59yHH1L+8T9wVVcDYEhKJPxH/4+gm+ah6aTvh84kNUedTVHAZvHOpiiXVNRHH30UvV7P7t27ueOOOxg0aBApKSncfPPNrFixgnnz5gHqCKDFixdz00034efnx8svvwzA4sWL6du3L0ajkQEDBvDRRx81+xgUFixYQGJiIiaTidjYWJ588kn38bfffpt+/fq5a6tuu+02AG688UYiIiJYsmSJR1lramr4/PPPeeihhwDIyMhAo9HwzTffMGrUKEwmE5s3b+bEiRPcfPPNREVF4e/vz5gxY1izZk2L915dXc3dd9+Nn58fcXFxvPXWWx7HKyoqePjhh4mIiCAwMJAZM2awf/9+j3O+/vprxowZg4+PD+Hh4dxyyy2X9PkLIURPYi8pofi118meOYuzi9/BVV2NMbUvsW+8Qd8VKwj+wa09MhiB1Bx1PnstvBLrndd+vgCM7evwdvbsWVatWsUrr7yCn1/rz2k+LHrBggW8+uqrvPnmm+j1epYtW8ZTTz3Fm2++yaxZs1i+fDkPPvgg8fHxTJ8+nS+++II//OEPfPrppwwePJiioiJ3uNi9ezdPPvkkH330ERMnTuTcuXNs2rQJAL1ez3333ceSJUv45S9/6S7D559/jtPp5O677/Yo4y9+8Qt++9vfkpKSQkhICLm5uVx//fW8/PLLmEwmPvzwQ+bNm8exY8dITEx0P++NN97g+eef56WXXuK7777jqaeeon///syePRuA22+/HbPZzDfffENQUBDvvvsuM2fO5Pjx44SGhrJixQpuueUWfvnLX/Lhhx9is9lYuXJlO/+hhBCi57Dn53P2b3+j4l9foNhsAJiuGUT4j39MwKxZaLQ9v95FoyiXWL0gqKqqIigoiMrKSgIDA9376+vrOXXqFMnJyfj4+Kg7bZYeEY527NjB+PHj+fLLLz1qPMLDw6mvrwfgscce47XXXkOj0fD000/zhz/8wX3epEmTGDx4MH/5y1/c++644w4sFgsrVqzg97//Pe+++y6HDh1q0QT35Zdf8uCDD5KXl0dAK0M6s7KyGDRoEOvXryc9PR2AqVOnkpSU5K6dysjIYPr06Xz11VfcfPPNF3yvQ4YM4cc//jGPP/44oHbIHjRoEN988437nLvuuouqqipWrlzJ5s2bueGGGygpKcFkMrnPSU1N5Wc/+xk/+tGPmDhxIikpKXz88ccXfG3RtVr9mRRCXBbrqVOcfe+vVP7nP+BwAGAeMYLwn/wYvylTesS8Ym19f59Pao46m8FXDSneeu0rtHPnTlwuF/fccw9Wq9W9f/To0R7nHT16lB/96Ece+yZNmsQf//hHQK15efPNN0lJSWHu3Llcf/31zJs3D71ez+zZs0lKSnIfmzt3Lrfccot7duOBAwcyceJE/v73v5Oenk52djabNm3iV7/6VYvynl+umpoaFixYwIoVKygsLMThcFBXV0dOTo7HeRMmTGjx+M033wRg//791NTUtFgSpq6ujhMNs7tmZmbyyCOPXPCzFEKInqj+2HHOvvsuVd9+Cy4XAL4TxhP+45/gO3ZMjwhFl0rCUWfTaNpde+NNqampaDQajjWsb9MopWHRP/N5IwzaanprS0JCAseOHWPNmjWsXr2aRx99lDfeeIMNGzYQEBDA3r17ycjIYNWqVbzwwgssWLCAXbt2EdwwY+pDDz3EE088wVtvvcX7779P3759mTZtWovXOb9czz77LKtXr+a3v/0tqampmM1mbrvtNmwNVcHtUVNTQ0xMDBkZGS2ONZbv/M9HCCF6urr9+yl7511q1q937/NPTyf8x/8P8/Dh3itYF+j5DYOiQ4SFhTF79mwWLVqExWK55OcPGjSILVu2eOzbsmUL11xzjfux2Wxm3rx5/OlPfyIjI4Nt27Zx8OBBQO1bNGvWLF5//XUOHDjA6dOnWbdunfu5d9xxB1qtlk8++YQPP/yQ//7v/27XXytbtmzhgQce4JZbbmHo0KFER0dz+vTpFudt3769xeNBgwYBMHLkSIqKitDr9aSmpnps4eHhAKSlpbF27dr2fVhCCNFNKYqCZft2zjzwIKfvvEsNRhoNAXPnkrzsSxLeWdzrgxFIzZFo5u2332bSpEmMHj2aBQsWkJaWhlarZdeuXWRlZTFq1Kg2n/s///M/3HHHHYwYMYJZs2bx9ddf8+WXX7pHhi1ZsgSn08m4cePw9fXl448/xmw2k5SUxPLlyzl58iRTp04lJCSElStX4nK5GDBggPv6/v7+3HnnnTz33HNUVVXxwAMPtOs99evXjy+//JJ58+ah0Wj4v//7P1wN1cLNbdmyhddff5358+ezevVqPv/8c1asWAHArFmzmDBhAvPnz+f111+nf//+FBQUuDthjx49mhdffJGZM2fSt29f7rrrLhwOBytXruTnP//5JfwLCCGEdyiKQs36DMrefYf6/Q1zx+n1BN10E2EPP4wpJdm7BexqirhklZWVCqBUVlZ67K+rq1OOHDmi1NXVealkV66goEB5/PHHleTkZMVgMCj+/v7K2LFjlTfeeEOxWCyKoigKoCxbtqzFc99++20lJSVFMRgMSv/+/ZUPP/zQfWzZsmXKuHHjlMDAQMXPz08ZP368smbNGkVRFGXTpk3KtGnTlJCQEMVsNitpaWnK0qVLW1x/69atCqBcf/31LY6tX79eAZTy8nKP/adOnVKmT5+umM1mJSEhQVm0aJEybdo05amnnnKfk5SUpLz00kvK7bffrvj6+irR0dHKH//4R4/rVFVVKU888YQSGxurGAwGJSEhQbnnnnuUnJwc9zlffPGFMnz4cMVoNCrh4eHKrbfeetHPW3Su3vAzKURncjkcSsXy5cqJeTcpRwYMVI4MGKgcTRumFP7q14otL8/bxetwbX1/n09Gq12GSxqtJoTwGvmZFKJ1LpuNyn//m7N//Sv2M+oAFa2fHyE/vJvQ++9H39BloLeR0WpCCCGE8OCyWCj/7HPOvf8+jpISAHRBQYTcfx+h99yDLijIyyXsHiQcCSGEEL2cs6KCcx//g/KPPsJZWQmAPiqK0AcfIOT229Fe4gjk3k7CkRBCCNFL2YtLOLdkCeVLl6I0LB5uSEok7OGHCbr5ZrRGo5dL2D1JOBJCCCF6GduZM5z929+pXLYMxW4HwDRwIOE/eoSAOXPQ6HReLmH3JuFICCGE6CXqjxzh7F//StW337lnszaPGkX4//tRj1niozuQcCSEEEL0YIqiULtrF2ff+yuWhkW7AfymTSX8kUfwPW9ZJXFxEo6EEEKIHkhxuajJyODsX96jLjNT3anVEnjddYQ98jA+Awd6tXw9mYQjIYQQogdR7HaqVq7k7F//ivX7bAA0RiNBt95C2EMPYUxI8HIJez5ZW010uIyMDDQaDRUVFW2es2TJEveirW1ZsGABw6+CNXw6knz2QvRertpazn34Edlz5lDw819g/T4brZ8fYY88TOraNcQsWCDBqINIOBIeioqKeOqpp0hNTcXHx4eoqCgmTZrE4sWLqW0YBtpVnn322Q5bzFWj0VxwW7BgQYe8zuXqTp87XF2fvRDdnaO8nNJFb5E9YybFr7yCo6AQXVgYEf/f/0fq+nVE/vSn6CMivF3MXkWa1YTbyZMnmTRpEsHBwbzyyisMHToUk8nEwYMH+ctf/kJcXBw33XRTl5XH398ff3//DrlWYWGh+/7SpUt54YUXOHbsmMdrdZYlS5awZMkSMjIyWj3e3T536D2fvRA9mT0/n7NLPqDiX/9CqasDwJCQQNhD/03Q/PloZUmcTiM1R8Lt0UcfRa/Xs3v3bu644w4GDRpESkoKN998MytWrGDevHmcPn0ajUZDZmPnP6CiogKNRtPiy3/Lli2kpaXh4+PD+PHjOXToUIvX/Oqrr+jXrx8+Pj7MmTOH3Nxc97Hzm3YeeOAB5s+fz29/+1tiYmIICwvjsccew94whwfA22+/7b5eVFQUt912GwDR0dHuLSgoCI1GQ3R0NGazmbi4OLKysgBwuVyEhoYyfvx49zU//vhjEhqqqltrtsrMzESj0XD69OlL/ciB9n3ugHz2nfDZC9Ed1R8/TsHPf072tXMo/+gjlLo6TNcMIu73v6PvNysJuesuCUadTGqOOpmiKNQ56rzy2ma9ud1zWpw9e5ZVq1bxyiuv4NfGNPKXOj/G//zP//DHP/6R6Ohonn/+eebNm8fx48cxGAwA1NbW8vLLL/Phhx9iNBp59NFHueuuu9iyZUub11y/fj0xMTGsX7+e7Oxs7rzzToYPH84jjzzC7t27efLJJ/noo4+YOHEi586dY1OzYa2tCQoKYvjw4WRkZDB69GgOHjyIRqNh37591NTU4O/vz4YNG5g2bdolvff26ozPHeSzF6KnURSFuj17OPvXv1HT7I8d3/HjCXvkYfwmTpQ5irqQhKNOVueoY9wn47zy2jt+uANfg2+7zs3OzkZRFAYMGOCxPzw8nPr6egAee+wxfvKTn7T79V988UVmz54NwAcffEB8fDzLli3jjjvuAMBut7No0SLGjRvnPmfQoEHs3LmTsWPHtnrNkJAQFi1ahE6nY+DAgdxwww2sXbuWRx55hJycHPz8/LjxxhsJCAggKSmJESNGXLSc6enpZGRk8Oyzz5KRkcHs2bPJyspi8+bNzJ07l4yMDH72s5+1+31fivZ+7q+99tolXVc+eyF6BsXlonrtWs799W/U7d+v7tRoCJg9m7BHHsY8dKh3C3iVkmY1cUE7d+4kMzOTwYMHY7VaL+m5EyZMcN8PDQ1lwIABHD161L1Pr9czZswY9+OBAwcSHBzscc75Bg8ejK7ZtPcxMTGUNKwsPXv2bJKSkkhJSeHee+/lH//4R7s6M0+bNo3NmzfjdDrZsGED6enp7i/tgoICsrOzSU9Pb/f7zsnJcffZ8ff358c//jGbNm3y2PfKK69c8BpX8rnD1fvZC9FTuGw2yj//nJPX30D+E09St38/GqOR4DvuoO83K4n/0x8lGHlRj645WrhwIV9++SVZWVmYzWYmTpzIa6+91uKv8OaWLFnCgw8+6LHPZDK5/0rvaGa9mR0/3NEp127Pa7dXamoqGo3Go6MsQEpKinots3otrVbN04qiuM9p3u+kszU2CzXSaDS4GqbIDwgIYO/evWRkZLBq1SpeeOEFFixYwK5duy44dH3q1KlUV1ezd+9eNm7cyCuvvEJ0dDSvvvoqw4YNIzY2ln79+gHte/+xsbEe/YK+/PJLvvjiC/7xj3+494WGhgLt/9zb+9qdqSd89kJ0d87qaso//ZRzH36Is7QMAG1AACF3303ovf8lo866iR5dc7RhwwYee+wxtm/fzurVq7Hb7Vx77bVYLJYLPi8wMJDCwkL3dubMmU4ro0ajwdfg65XtUtqnw8LCmD17NosWLbrg5xfR8IPbfARS8yDQ3Pbt2933y8vLOX78OIMGDXLvczgc7N692/342LFjVFRUeJxzqfR6PbNmzeL111/nwIEDnD59mnXr1l3wOcHBwaSlpbFo0SIMBgMDBw5k6tSp7Nu3j+XLl3v0eWnP+9fr9aSmprq3yMhIzGazx77GcNTez729r93oav3sheiu7EVFFL/+Btnp0yn93e9xlpahj4oi8uc/J3X9eiKf+f8kGHUjPbrm6Ntvv/V4vGTJEiIjI9mzZw9Tp05t83mNo2WEp7fffptJkyYxevRoFixYQFpaGlqtll27dpGVlcWoUaMwm82MHz+eV199leTkZEpKSvjf//3fVq/3q1/9irCwMKKiovjlL39JeHg48+fPdx83GAw88cQT/OlPf0Kv1/P4448zfvz4Nvu8XMzy5cs5efIkU6dOJSQkhJUrV+JyuS5Yk9goPT2dP//5z+4RVqGhoQwaNIilS5fy1ltvuc9LTU0lISGBBQsW8PLLL3P8+HF+97vfXVZ5G7Xncwfks++Ez16IzlZ//Djn/v4+lStWQENNpzG1L2EPPUzQDdejMRq9XELRmh5dc3S+yspKoKnJoi01NTUkJSWRkJDAzTffzOHDhy94vtVqpaqqymPrjfr27cu+ffuYNWsWzz33HMOGDWP06NH8+c9/5tlnn+XXv/41AH//+99xOByMGjWKp59+mt/85jetXu/VV1/lqaeeYtSoURQVFfH1119jbPaLwNfXl5///Of88Ic/ZNKkSfj7+7N06dLLLn9wcDBffvklM2bMYNCgQbzzzjv885//ZPDgwRd97rRp03A6nR79W9LT01vsMxgM/POf/yQrK4u0tDRee+21Nt9/e7X3cwf57Dv6sxeiMyiKgmXHTnJ+9CNO3XQzlV99BXY7vqNHE7/4bVL+8x+Cb5kvwagb0yjNG/B7MJfLxU033URFRQWbN29u87xt27bx/fffk5aWRmVlJb/97W/ZuHEjhw8fJj4+vtXnLFiwgJdeeqnF/srKSgIDA92P6+vrOXXqFMnJyfjIHBRCeJ38TIqupDgcVK9ezdm//Z36xrnFGkeePfTfmIcN824BexCny4lOq7v4iZeoqqqKoKCgFt/f5+s14egnP/kJ33zzDZs3b24z5LTGbrczaNAg7r77bo+/0JuzWq0eI4aqqqpISEiQcCRENyc/k6IruGprqVi2jHNLPsDeMJmqxmRSF4J94AGMSUleLmHPUGmtZH3ueladXkWFtYJPbvikw1+jveGoR/c5avT444+zfPlyNm7ceEnBCNSq+hEjRpCdnd3mOSaTCZPJdKXFFEII0Ys4ysoo/+QTyv/xCc6Gbh264GBCfvhDQv7rHvQX6eIh1EC0Lmcdq86sYnvhdhwuh/tYfk0+cf5xXilXjw5HiqLwxBNPsGzZMjIyMkhOTr7kazidTg4ePMj111/fCSUUQgjR21hPnuTc+0uo/Pe/UWw2QF3zLPSB+wm+9Va05vZPo3I1agxE3535jh0FO3AoTYEoNTiVa/tcy7VJ13otGEEPD0ePPfYYn3zyCf/+978JCAigqKgIUJclaJwf5r777iMuLo6FCxcC6iie8ePHk5qaSkVFBW+88QZnzpzh4Ycf9tr7EEII0b25l/f429+pWb/evd9nWBph//0QAbNmotF1fB+Z3qK8vtzdZLaj0DMQ9Qvpx7VJ13Jtn2tJCUrxYimb9OhwtHjxYoAWM+i+//77PPDAA4A6W3Hj5HGgzvnyyCOPUFRUREhICKNGjWLr1q1cc801XVVsIYQQPYTicFC9Zg1n//4+9QcOqDs1GvxnzCDsvx/EPHKkrHnWhnP151ibs5bVp1ezs2gnTsXpPtY/pL87ECUHXXqrT2frNR2yu1JbHbqk86cQ3Yv8TIrL5ayxUPnlF5z74EPs+fkAaIxGgm65hdD778eU0v2+0LuDsroy1p5Zy+ozq9lVvAuX4nIfGxQ6iNlJs5mdNJs+QX28Ur6rqkO2EEII0RHsxcWUf/QR5Us/w1VdDYAuJISQu+8m5J4fog8L83IJu5+S2hLWnFnD6jOr2VO8B4WmOpfBYYPdgSgxMNGLpbw0Eo6EEEJc9eqzsjj3/vtUrlgJDrU/jLFPH0IfeICg+TejlZpHD0WWIlafWc3qM6vJLMn0CERDw4dybdK1zEqaRXzApY0g7y4kHAkhhLgqKS4Xls2bOfv++9Rua1qP0HfMGEIffBD/9GlotL1qIYkrkled564hOlB2wOPY8IjhzEqaxeyk2cT6x3qphB1HwpHocBkZGUyfPp3y8vI2V2RfsmQJTz/9NBUVFW1eZ8GCBXz11VeyuOgV6tOnD08//TRPP/20t4siRLfgqq+n8j//4dwHH2I7cULdqdMROHcuoQ88gHnoEO8WsBs5XXmaNTlqIDpy9oh7vwYNI6NGMjtpNrMSZxHlF+XFUnY8CUfCQ1FREQsXLmTFihXk5eURFBREamoq//Vf/8X999+Pr69vl5Xl2Wef5YknnuiQa11sNMmLL77IggULOuS1LtUDDzxARUUFX331lVdeX4irhTpp4z8p/+c/cZaXA6D18yP49tsJve9eDLE9v8bjSimKwomKE6zOUZvMvi//3n1Mq9EyJmoMs5NmMzNpJuHmcC+WtHNJOBJuJ0+eZNKkSQQHB/PKK68wdOhQTCYTBw8e5C9/+QtxcXHcdNNNXVYef39//P39O+RahYWF7vtLly7lhRde4NixYx6v1VmWLFnCkiVLyMjI6LTXEEK0rf74cc598AFV//kaxW4HwBAbS8h99xJ8223oOvHnvydQFIVj5cdYdXoVa3LWcKrylPuYXqNnbMxYZiXNYmbiTEJ9ro5Zv6UxVbg9+uij6PV6du/ezR133MGgQYNISUnh5ptvZsWKFcybN4/Tp0+j0Wg8mroqKirQaDQtvvy3bNlCWloaPj4+jB8/nkONCzE289VXX9GvXz98fHyYM2cOuQ3rEoHarDZ8+HD34wceeID58+fz29/+lpiYGMLCwnjsscewN/yyA3j77bfd14uKiuK2224DIDo62r0FBQWh0WiIjo7GbDYTFxdHVlYWoC5gHBoayvjx493X/Pjjj0lISADUJkONRuPRHJiZmYlGo+H06dOX+pFf1JIlS1o0TX711VctasK+/vprxowZg4+PD+Hh4dxyyy1tXvOvf/0rwcHBrF27FoDf//73DB06FD8/PxISEnj00Uepqalxn3/mzBnmzZtHSEgIfn5+DB48mJUrV3bcmxSiEyiKQs2mzeQ89DCnbrqZyi++RLHb8RmWRtybf6Dvqu8Ie+CBqzYYKYrCgdID/H7377n+y+u5/evbee/ge5yqPIVBa2Ba/DR+PenXZNyZwbuz3+X2/rdfNcEIpOao0ymKglJX55XX1pjN7Z6c7OzZs6xatYpXXnkFPz+/1q93iROd/c///A9//OMfiY6O5vnnn2fevHkcP34cg8EAQG1tLS+//DIffvghRqORRx99lLvuuostW7a0ec3169cTExPD+vXryc7O5s4772T48OE88sgj7N69myeffJKPPvqIiRMncu7cOTZt2nTBMgYFBTF8+HAyMjIYPXo0Bw8eRKPRsG/fPmpqavD392fDhg1Mmzbtkt57V1qxYgW33HILv/zlL/nwww+x2WxthpfXX3+d119/nVWrVjF27FgAtFotf/rTn0hOTubkyZM8+uij/OxnP+Ptt98G1JnobTYbGzduxM/PjyNHjnRqTZsQV8Ldn+jDD7FlN/Qn0moJmD2b0Afux3fECO8W0IucLif7SvaxJmcNa86sobi22H3MR+fD5LjJzEqaxbT4afgbr+6fcQlHnUypq+PYyFFeee0Be/egaWcfoezsbBRFYcCAAR77w8PDqa+vB9QvyZ/85Cftfv0XX3yR2bNnA/DBBx8QHx/PsmXLuOOOOwCw2+0sWrSIcePGuc8ZNGgQO3fudH9xny8kJIRFixah0+kYOHAgN9xwA2vXruWRRx4hJycHPz8/brzxRgICAkhKSmJEO34Rpqenk5GRwbPPPktGRgazZ88mKyuLzZs3M3fuXDIyMvjZz37W7vfd1V5++WXuuusuXnrpJfe+YcOGtTjv5z//OR999BEbNmxg8ODB7v3NO2r36dOH3/zmN/z4xz92h6OcnBx+8IMfMHToUABSUrrH9P5CNGcvKaH8k0+o+HQpzoaaXa2vL8G330bIvfdivMRFyXsLu8vOrsJdrMlZw9qctZyrP+c+5qv3ZVr8NGYlzWJy3GR8DV3Xp7S7k3AkLmjnzp24XC7uuecerFbrJT13woQJ7vuhoaEMGDCAo0ePuvfp9XrGjBnjfjxw4ECCg4M5evRom+Fo8ODB6JqtXxQTE8PBgwcBmD17NklJSaSkpDB37lzmzp3LLbfcctFO5NOmTeNvf/sbTqeTDRs2cO211xIdHU1GRgZpaWlkZ2e3WKLmQnJycjyWo3E4HNjtdo/alueff57nn3++3de8kMzMTB555JELnvO73/0Oi8XC7t27W4SbNWvWsHDhQrKysqiqqsLhcFBfX09tbS2+vr48+eST/OQnP2HVqlXMmjWLH/zgB6SlpXVI2YW4UvVHjnDugw+oXPkNNO9PdO+9BN/2A3QBAV4uYdezOq1sK9jG6jOrycjNoMpW5T4WYAxgesJ0ZifNZkLsBEw6k/cK2o1JOOpkGrOZAXv3eO212ys1NRWNRuPRSRmaagkaF/JtXKeu+aozzfv8dLbGJrlGGo0Gl0udnj4gIIC9e/eSkZHBqlWreOGFF1iwYAG7du1qc0oBgKlTp1JdXc3evXvZuHEjr7zyCtHR0bz66qsMGzaM2NhY+vXrB7Tv/cfGxnr0yfryyy/54osv+Mc//uHeFxravrZ7rVbL+Sv8nP965nb8O0+ZMoUVK1bw2Wef8Ytf/MK9//Tp09x444385Cc/4eWXXyY0NJTNmzfz0EMPYbPZ8PX15eGHH2bOnDmsWLGCVatWsXDhQn73u9912EhCIS6V4nRSs34955Z8QO3u3e795pEjCb3/fgJmzkCjv7q+3mrttWzM38jaM2vZmLeRWket+1ioTygzEmcwO3E2Y2LGYNAaLnAlARKOOp1Go2l305Y3hYWFMXv2bBYtWsQTTzzRZr+jiIgIQB391dhk1dY8RNu3bycxUZ0uvry8nOPHjzNo0CD3cYfDwe7du921RMeOHaOiosLjnEul1+uZNWsWs2bN4sUXXyQ4OJh169Zx6623tvmc4OBg0tLSWLRoEQaDgYEDBxIZGcmdd97J8uXLPfobNX//ISEhrb5/vV5Pamqq+3FkZCRms9ljX3tFRERQXV2NxWJx/5uc/3ppaWmsXbuWBx98sM3rjB07lscff5y5c+ei1+t59tlnAdizZw8ul4vf/e537uD32WeftXh+QkICP/7xj/nxj3/Mc889x3vvvSfhSHQ5Z00NlV9+ybmPPsbeOHhDr1fnJ7r/PswNTb9Xi0prJRvyNrDmzBq2FmzF6myq3Y/0jWRW4ixmJc1iZORIdFrdBa4kzifhSLi9/fbbTJo0idGjR7NgwQLS0tLQarXs2rWLrKwsRo0ahdlsZvz48bz66qskJydTUlLC//7v/7Z6vV/96leEhYURFRXFL3/5S8LDw5k/f777uMFg4IknnuBPf/oTer2exx9/nPHjx7fZpHYxy5cv5+TJk0ydOpWQkBBWrlyJy+Vq0Y+qNenp6fz5z392j24LDQ1l0KBBLF26lLfeest9XmpqKgkJCSxYsICXX36Z48eP87vf/e6yyttcZWVli9ATFhbGuHHj8PX15fnnn+fJJ59kx44dLFmyxOO8F198kZkzZ9K3b1/uuusuHA4HK1eu5Oc//7nHeRMnTmTlypVcd9116PV6nn76aVJTU7Hb7fz5z39m3rx5bNmyhXfeecfjeU8//TTXXXcd/fv3p7y8nPXr119RgBXiUtlycjj38cdUfvElLosFAG1QECF33EHIPT/EEB3t5RJ2nbK6MtblrGNtzlp2Fu7EoTjcxxICEpiVNItZibMYEj4ErUYGpF82RVyyyspKBVAqKys99tfV1SlHjhxR6urqvFSyK1dQUKA8/vjjSnJysmIwGBR/f39l7NixyhtvvKFYLBZFURTlyJEjyoQJExSz2awMHz5cWbVqlQIo69evVxRFUdavX68Aytdff60MHjxYMRqNytixY5X9+/e7X+f9999XgoKClC+++EJJSUlRTCaTMmvWLOXMmTPuc1588UVl2LBh7sf333+/cvPNN3uU96mnnlKmTZumKIqibNq0SZk2bZoSEhKimM1mJS0tTVm6dGmL99j42s0tW7ZMAZTFixd7XBtQsrKyPM7dvHmzMnToUMXHx0eZMmWK8vnnnyuAcurUqVY/0/fff99dxtbcf//9CtBie+ihh9xlS01NVcxms3LjjTcqf/nLX5Tzf3S/+OILZfjw4YrRaFTCw8OVW2+91X0sKSlJ+cMf/uB+vGHDBsXPz0/505/+pCiKovz+979XYmJiFLPZrMyZM0f58MMPFUApLy9XFEVRHn/8caVv376KyWRSIiIilHvvvVcpKytr8/10J73hZ/Jq5XK5lJpt25WcnzyqHBk4SDkyYKByZMBAJfv6G5Rz//xUcdbWeruIXSavOk/54NAHyn0r71OGLhmqDFkyxL3N/2q+smjfIiXrbJbicrm8XdRur63v7/NpFOW8Dg3ioqqqqggKCqKyspLAwED3/vr6ek6dOkVycjI+skihEF4nP5M9j8tqpWr5cs59+BHWZn0g/aZOIfS++/GbNPGSpxXpiU5WnHQPuT967qjHsSFhQ5iZNJNZibPoE9THOwXsLNZqOL0ZKvNg7IUHmlyOtr6/zyfNakIIIbzOXlxM+T//ScXSz9xLe2jMZoJvmU/If92LKSXZyyXsXIqicOTsEfeQ++azVGs1WkZFjWJm4kxmJs4k2q8XNSO6nFCQCSfWwcn1kLsDXA7Q+8CIe8HgnT9qJBwJIYTwCkVRqMvMpPyjj6latQocav8ZfWwMoff8lzoUPyjIy6XsPA6XQ52U8cwa1uWuo8hS5D6m1+qZEDOBWUmzSE9I712zU1fkqGHoxHo4mQH1FZ7HQ/pA3xlgs0g4EkIIcXVw2WxUf/st5z78iPpmywr5jh5NyH33EjCj9w7FtzqtbC/YzpqcNWTkZlBhrXAfM+vNTI6bzMzEmUyNn0qAsZfM0VRfpTaVNdYOnc32PG4KguQpaiDqOx1CvT/RbO/83yeEEKLbcZSWUv7pUsqXLsVZVgaAxmgk8MYbCb33v/DppaMga2w1bMzbyNqctWzO3+wxB1GQKYj0+HRmJc1ifMx4fPS9oG+c0wEF+5o1le0Exdl0XKOD+NENYWgGxI4EXfeKI92rNEIIIXqduv37OffxP6j69lv3LNb6yEhCfng3wXfcgb6dk6L2JGV1ZazPXc/anLXsKNyBw9U05D7SN9Ldf2hU1Cj02l7wVXzuZEMz2Xo4uRGslZ7HQ1MgZboahvpMBnOwV4rZXr3gX6T7aZyxWQjhXfKz6D0um43qb77h3Mf/oL5hiR8A87BhhNx3L4HXXovG0Ltmas6tznXPQZRZkolC02DwPoF9mJk4kxmJM3rHHES15+DURjUMnVgPFWc8j/sEQUp6QyCarvYj6kEkHHUgo9GIVquloKCAiIgIjEbjVTHkVIjuRlEUbDYbpaWlaLVajEajt4t01bAXl1Cx9FPKl36G8+xZADQGA4HXX0/If/0X5qFDvFzCjqMoCsfKj7kD0fHy4x7HB4cNdtcQpQR7vx/NFXHYIG+nGoROrIPCTFCa/fGh1UPCuKbaodjh0INn5ZZw1IG0Wi3JyckUFhZSUFDg7eIIcdXz9fUlMTHRvTSK6ByKolC3L5Pyj88bdRYZScjdd6lNZ2FhXi5lx2gcYbYuZx3rc9eTX5PvPqbT6BgdNZoZiTOYkTijZw+5VxQoOaqOJju5Hk5vAbvF85zwAWqtUMp06DMJTL2kAzkSjjqc0WgkMTERh8OB0+m8+BOEEJ1Cp9Oh1+ul9rYTuerrqVqxgnP/+AfWI00TFZpHjSL0v+4hYNasXtF0Vu+oZ1vBNtblrmsxwsxH58OkuEnMSJzB1LipBPsEe6uYV666SA1DjUPsa4o8j/tFNDWVpaRDUFzXl7GLSDjqBBqNBoPB0GIFeSGE6A1seXmU//OfVP7rC5yVasdbjclE4A03EPpf9+BzzTVeLuGVq6ivYGP+RtblrGNrwVbqHHXuY0GmIKbFT2Nm4kwmxE7ArDd7saRXwFoDZ7Y2dKLOgJIjnsf1PpA0sanfUORguEpqYSUcCSGEuCjF5cKyZSvl//gHNRs2qM0ugCEujpAf3k3QrbeiDwnxcimvTEFNAetz17MuZx17ivfgbDb8PMYvhhmJM5iZOJMRkSN65gizxiH2jWEodye47M1O0EDMsKamsoRxXpuE0ds69V9306ZNTJkyhS1btjBp0qTOfCkhhBCdwFlVReVXX1H+j0+wnWkakeQ3aRIh99yD/7SpaHQ9s+OtoigcLz/Oupx1rMtdR9a5LI/j/UP6q/2HEmYwMHRgz2uiVRQ4e6IpDJ3a1HKIfXBiU81Q8jTw7X3TKlyOTg1H33zzDXq9nhUrVkg4EkKIHqQ+K4vyf3xC5fLlKHVqk5LW35+gW28h5O67MSX3zLXO7C47e4v3sj53Petz1lNgaRo8o9VoGRk5khmJM5ieMJ34gHgvlvQy1ZQ0DbE/uQEqcz2P+wRD8tSG2qH0bjEbdXfUaeHopZdewuFwMGPGDJ588kl+9atf8cILL3TWywkhhLhCis1G1XerKP/nP6nbu9e939QvlZB77iFo3jy0fn5eLOHlsdgtbMnfwvrc9WzM20iVrcp9zEfnw/jY8cxMnMm0+GmE+PSwpkFrDeRsaxhVlgHFhzyP64xq81hjGIoZ3qOH2HcVjaIoysVPuzzvvfcelZWVBAcH8/DDD3fWy3S5qqoqgoKCqKysJDAw0NvFEUKIK2IvKKD8s8+o+Pxf7rmJ0OsJmD2L0B/+EPPo0T2uSam0tpSMvAzW56xne+F27M361oSYQpiWMI3pCdN7Xodqpx3y98KpDW30GwKi0xpGlaVD4gQw+nqhoN1Te7+/O7VZzeFw8Oyzz/Luu+925ssIIYS4RIrLhWXbNsr/+U9q1q2HhtnE9ZGRBN95B8G3344hMtLLpWw/RVE4UXGC9bnrycjN4EDZAY/jCQEJzEiYwfTE6QyPGI6up9SeKAqUZqlNZCcz1AVcbdWe5wQlQt90NQwlTwO/cC8UtHfp1Jqj3kpqjoQQPZWzooKKZV9R/uk/sZ/Jce/3HTeOkLvvJmDmjB4zN1HjhIyNgSi32rN/zZCwIe7+Q32D+3Z47Vd+RR3lFlubx0P8jMQFX0atVGV+U83QyQyoKfY8bg5R+w0lT2tYmiMZeljNnrd0i5qjRmvWrGHmzJk9rlpWCCF6A0VRqD9wgPJ/fkrVN9+gWK0AaP38CLr5ZkJ+eDem1FQvl7J9LHYLWwu2kpGbwYa8DVQ2G31l1BoZFzOO6YnTmRY/jUjfzqv5yq+oY8ZvM7A62l6/z6TXsu7Z9IsHpNpzao3QyQw1FJ3N9jyu91GbxxqbyqLTrpr5hrylS8LRnDlzKCwsJLIHVdEKIURP56qtpXLFCir++Sn1R5om+DMNGkTIXXcRdOMNPaKDdbGlmA15G1ifu54dhTs8+g81TsiYnpDOpNhJ+Bq6pn9NucV2wWAEYHW4KLfYWoYjWy3kbm9qKivcD80WqUWjhdgRTWEofuxVO9+Qt3RJOJKWOyGE6DrW7GzKP11K5b//jata7Z+iMRoJvO46Qu6+C59hw7p1TX7jgq6NzWVHznrO3JwQkEB6QjrTE6b3jAkZnQ4o2KuGoVMbIHcHOM9rjgsfACnT1KayPpPBHOyVogpVN/8fJYQQoj1cNhvV362ifOmn1O3e495vSEwk5M47Cbr1lm49g7XNaWN30W7W565nQ94GCi2F7mMaNAyLGOYORMlByd063IHCAE0uk7SHSfzub1C0s2Un6sA4NQilTFP7DwXGeqeoolVdFo4WL17M5MmTGTlyJCHd+AdUCCF6EtuZM5R/9hmVXy7DWV6u7tTpCJgxneA77sRv0kQ03bR/SkV9BZvyN7E+dz1b8rdQ66h1H/PR+TAhdgLTE6YzNX4qYeYwL5b04hI0xUzSHmaS9hATtEcI1zTMpdTY590nGJKnNASidAhLlU7U3ViXhaNFixbx0ksvodFoSEhIYOTIkR5bdHR0VxVFCCF6NMVup3r9eio+XYpl61b3fn10NMG330bwbbdhiIryYgnbdqryFBm5GWTkZpBZmolLaeq3E2GOYGr8VNIT0hkfMx4fffftZ6OvLeYm7RYmaQ8zUXuYBG2px/FaxcQu1wAGTLiR6OFzIHqoTL7Yg3RZODp8+DAOh4N9+/axd+9e9u7dy3vvvUdubi4ajYbo6Gjy8/Mv6ZoLFy7kyy+/JCsrC7PZzMSJE3nttdcYMGDABZ/3+eef83//93+cPn2afv368dprr3H99ddfydsTQohOZ8vLp+Jfn1P5xZc4Shu+jDUa/KZOIeTOu/CfOgWNvnv1lmgcbt84uuxM1RmP4wNCBpCekE56QjrXhF2DVtM9a7moK4fTW9Q+Q6c2MrA0iz8Zmw7bFR37lFS2ugazxTmETCUVO3qWD5tMdGyQ98otLkuX/BQ1tg3HxsYSGxvLDTfc4D529uxZ9uzZQ2Zm5iVfd8OGDTz22GOMGTMGh8PB888/z7XXXsuRI0fwa2MExtatW7n77rtZuHAhN954I5988gnz589n7969DBky5LLenxBCdBbFbqc6I4OKzz7HsnmzOikgoAsPJ/gHPyD49tsxxsd5uZSeKq2VbMnfQkZeBpvzN1PdrL+NXqtnbPRY0hPSmRY/jVj/btrXxmZRl+U4tVHdCvdDs1ouBQ2HXElscQ1hm2swu1wDqKX71nSJS9Mlk0BqtVqKioo6fSh/aWkpkZGRbNiwgalTp7Z6zp133onFYmH58uXufePHj2f48OG888477XodmQRSCNHZWq0lAvwmTiT4jjsImDEdjdF4gSt0rTNVZ9y1Q3uL9+JUnO5jwaZgpsZPZVr8NCbGTsTf6O+9grbFYYW8XU1hKG93y2U5wvurfYaSp1IQOprpi/Z3zDxHost0q0kgv/32W4KCOr9asbJSnQwsNDS0zXO2bdvGM88847Fvzpw5fPXVV20+x2q1Ym2YNA3UD1cIITpam7VEYWEE33orwbffhjEx0culVNlddjJLMtmQu4ENeRs4XXXa43jfoL5MS1DnH0oLT+t+y3U4HVCwT20mO70JcraDo97znKBESJnaNLy+2YiyWGDdsyGdM0O28LouCUfXXnttp7+Gy+Xi6aefZtKkSRdsHisqKiLqvI6KUVFRFBUVtfmchQsX8tJLL3VYWYUQojlbTg4Vn/+Liq+W4Swtc+/vbrVEFfUVbC7YzMbcjWwuOK+5TKNnVPQo0uPTmZYwjYSABC+WtBUuJxQdVIPQqY1wZlvL4fX+UQ3LcjRsIX0ueMm4YLOEn16qe/XcuwKPPfYYhw4dYvPmzR1+7eeee86jtqmqqoqEhG72gy+E6FFcNhvVq1dT8fm/qN2+3b1fFx5O8C3z1b5EXq4lUhSFk5Un2ZC3gQ25G1qMLgsxhTAlfgpT46cyMXYiAcYAL5b2PIoCJUebwtDpzVBf4XmOOUStEWpoKiO8vwyvF0AvCUePP/44y5cvZ+PGjcTHx1/w3OjoaIqLPRfxKy4uvuBUAiaTCZPJ1CFlFUJc3azZ2VR8/jmVX/0bZ0NXADQa/CZPJvj22wiYPt2rC79anVZ2F+1mQ94GNuZtJL/GcxRxanCquzP10PCh3ae5TFGg7Hs4vRFObVLDUG2Z5zmmQEiaqAahPlMgaoisUSZa1aPDkaIoPPHEEyxbtoyMjAySk5Mv+pwJEyawdu1ann76afe+1atXM2HChE4sqRDiauaqraXq2++o+Ne/qNu7171fHx2t9iX6wa0Y4rw34qy0tpRN+ZvYkLuBbYXbqHPUuY8ZtAbGRo9lWsI0psZPJc6/m4yMUxQ4d7KhZqghDNWc1z1Cb4akCWoQSp4GMcNA16O/9kQX6dH/Sx577DE++eQT/v3vfxMQEODuNxQUFITZrLYD33fffcTFxbFw4UIAnnrqKaZNm8bvfvc7brjhBj799FN2797NX/7yF6+9DyFE76MoCvUHD1Lx+b+oWrkSl8WiHtDp8E9PJ/j22/CfMgWNrutrXlyKi8Nlh9mYv5GNeRtbrF3WOBnj1PipjI8Z32WLuV5U+Wk1BDWGoao8z+M6EySMbaoZihsFeu/31RI9T48OR4sXLwYgPT3dY//777/PAw88AEBOTg7aZtWmEydO5JNPPuF///d/ef755+nXrx9fffWVzHEkhOgQjvJyqr7+morP/4X1++/d+w2JiQTfeitBt9yCIapzpzVpTbWtmm0F29iYt5FN+Zs4V3/O4/iQsCFMTVCH2w8MHdg9JmOsyGkKQqc3Q2WO53GtAeLHqMty9Jmi3pfV60UH6JJ5jnobmedICNGc4nJh2baNin/9i5o1a1Hs6vw4GpOJgDnXEvyD2/AdM7pL1zhTFIVTVafYlLeJjXkb2Vu8F4ficB/3M/gxMXYiU+KmMCV+CuHm8C4rW5sqcpuC0OmNajhqTqtXa4P6TFbDUMI4MHaTWi3RI3SreY6EEKI3suXlUfnlMiq+WoajoGkVeZ9rriHoth8QdOON6LrwDyir08quol3uQJRX49ns1Cewj3t02ajIURh03uv4DZwXhjZBhefSImh0EDfSMwyZuuEEkqLXkXAkhBCXwFVXpw7B/3KZxxB8bUAAQfPmEXzbD/C55pouK09hTSGb8jexKW8TO4p2tOhMPTpqtLv/UGKglyeQbE8Yih3eEIamQuI4MHWj6QHEVUPCkRBCXIS7c/UXX1K1YgWumhr1gEaD34QJBN16KwGzZqL16fz+LnaXnf0l+9mUr9YOZVdkexyP9I1kStyU7tGZuvyMGoTObGkIQ+c1k3mEoYaaIR/pqiC8T8KREEK0wVFaSuXXy6n48gts2Sfc+w1xcQTdegvB8+d3yRD8sroyNudvZlPeJrYVbKPa3jSzs1ajZVjEMKbGT2VK3BT6h/R3L/bdpRRFHU12ZktDzdCWlh2oNTp1OH1jB2oJQ6KbknAkhBDNuGw2atZnULlsGTWbNoFTXUDV3bn61lvxHTu2UztXO11ODpYddDeXHT131ON4sCmYSXGTmBo3lUlxkwgydf7alS0oCpw9AWcagtCZrS2H1mv1EDtCrRlKmizNZKLHkHAkhLjqKYpC/ZEjVH65jKrly5tmrgbMw4YRdMstBN5wPbqAzvtiP1t3lq0FW9mUv4mtBVuptFZ6HB8cNpgp8VOYEjeFwWGDu35makWB0qymZrIzW6HGc7UBtIaG0WST1EAUP1Y6UIseScKREOKq5Sgro/Lr5VQuW4b1+HH3fn1kJEE330zQLfMxpaR0yms31g5tzt/M5vzNHDl7BIWmmVUCjAFMip3ElPgpTIyd2K6h9vkVdR23SrzLCcWH1RB0ZrN6W3vW8xydSZ1bqM8kdVmO+LEytF70ChKOhBBXFZfVSs26dVR+9W9qNm9uajYzGgmYNYugW27Bb+KETpm5+mK1Q4NCBzE5bjJT4qcwNHwoem37f0XnV9Qx47cZWB2uNs8x6bWseza99YDktEPh/oaaoa2Qsx3OKx96szoDdZ/JkDRJrSWSSRdFLyThSAjR6ymKQt2+fVR+9W+qvvkGV3VTh2aftDSCb72FwOuuQxfUsX13HC4Hh8oOsSl/E1vyt7SsHTIEMCF2ApPjJjM5bjIRvhGX/VrlFtsFgxGA1eGi3GJTw5G9DvL3wJltas1Q7k6w13o+wRgAiePV9cmSJqv9h2Q5DnEVkHAkhOi1bHl5VP7731T++z/Yc5pGTuljYgi66SaCbr6pw5vNSmpL2JK/hc35m9lWuI1qW7XH8cbaoclxk0mLSLuk2qEr4U8to7TfE7VrC3y3Vw1GzvOa4Mwhao1Q0kR1ixoqC7WKq5L8rxdC9CrOqiqqvv2Wyv/8h7rde9z7Nb6+BF57LUHz5+M7dkyHjTazO+1klma6A9Gx8mMexwONgUyMncjkuMlMjJ14RbVDlyKUKsZojzFWm8VY7VGu0ZxBp1Egs9lJ/tENtUKT1C1iIHThEidCdFcSjoQQPZ5is1GzaROV//4PNevXu9c2c0/SOP9mAmbNQuvbMZ2F86rz2Fqwlc35m9lRuINaR1NzlAYNg8MGMzl+MpNiJzE0fGjnjyxTFAzVudyi3cQYbRbjtFn01Ra2OO2MK5LAgdMIGZQOiRMgNAW8MSeSEN2chCMhRI+k9iPKpPLr/1C98huP4femfv0IuvkmAm+8EUN09BW/Vp2jjl1Fu9iSv4WtBVs5XXXa43ioT6i779DE2ImE+oRe8WtekMsFpUebOk7nbGNAVT5/OK87UJYrgZ2ugexyDWCnayDFhLJ82mRC4rwwL5IQPYiEIyFEj2I7fVodfv/11579iCIiCJw3j6Cb5mEaMOCKZolWFIXvK75na/5WthRsYW/xXmyupv45Oo2OYRHDmBw3mUlxkxgYOhCtphOboxxWKMiEnMYwtB3qKzzLrNGzz9nHHYZ2uwZQicwxJMTlkHAkhOj27CUlVH/zDZVfL6f+0CH3fo2vL4GzZxN08034jht3RcPvK+or2Fa4jS35W9hWsI2SuhKP4zF+MUyKm8Tk2MmMjRlLgLETZ3quq4C8XZCzTR1Nlr8HnFbPcwy+6hxDSRMhcQJHtP259Z29nVcmIa4iEo6EEN2Ss7qa6lWrqVqxHMv2HWpTEoBOh9/EiQTddBMBM2dcdj8iu8vOwdKDbCnYwtb8rRw+e9hjmL2PzodR0aOYFDuJSbGTSA5K7rw1yyrz3M1j5OyA4kPQrCwA+IarnacTJ6jD66PTQGdwHw6uqMOk1150nqMQPxmKL8TFSDgSQnQbLpuNmg0bqFq+Qu1YbWtqyjIPH07gvBsJnDsXfVjYZV0/typXDUMFW9lZtBOL3eJxvF9IPybFTmJi7ERGRo3EpDNd0ftplcsJJUeamsdytrdckwzUztKJE5q2sL4X7DwdF2xm3bPpHTdDthBXMQlHQgivUhwOLNt3ULViBdVr1nhM0Gjs25egeTcSeMMNGBMSLvna1bZqdhbuZGvBVrYWbCWvxjOEBJuCmRAzgYlxE5kYO5FI38grfj8t2CyQt1udZDF3u3prrfI8R6ODmDQ1BCWMU28Doi75peKCzRJ+hOgAEo6EEF1Ocbmo27ePqhUrqPr2O5znzrmP6aOiCLzhBoLm3Yhp4MBLaspqnJF6a8FWthVs42DZQZyKs+naGj3DI4czMXYiE+MmMih0UMd3pK4qUGuDcneot0UHoVkZAHXm6YQxTWEofjQY/Tq2HEKIyybhSAjRJRRFof7wEapWrqTqm29wFDbNw6MLCSFg7hyCrr8e86hR7Z6gUVEUcqpz2Fawja0FW9lVtIsae43HOX0C+6hhKHYio6NH42fowBDSuDhr7o6GMLQDKnNanhcYD4njIGG82l8oajB09txHQojLJuFICNFpFEXBevw4Vd98Q/U332I7c8Z9TOvvT8CsWQTecAN+48ehMRgucKUmFfUV7CjawbaCbWwv3E5+Tb7H8SBTEOOixzEhdgITYicQ5x/XcW+ovlIdRZa7Uw1DebvB5hnG0GghaogaghLGqbdB8R1XBiFEp5NwJITocNbsbKpWfkPVt99iO3nSvV9jMuE/fTqBN1yP/9SpaE0X7/BsdVrZV7KP7QXb2Va4jaNnj3qMKtNr9YyIHKH2HYqdyMDQgR0zI7WiwLmTTUEod6fakfr8UWSmQLVZLGFcUxOZqROH+QshOp2EIyFEh7CePEXVt99Q/c03WL/Pdu/XGI34TZ1C4Nzr8E9PR+d/4WYtl+Ii61wW2wu3s71gO3tL9mI9b46fvkF93TVDo6NG42vogGVBbBYo2NcQhHZB3k6oPdvyvJDkhiA0Vq0VihgoTWRC9DISjoQQl8168iTV331H1bffYT3WbMFVgwH/yZMJvG4u/jNmoPO/8EzNudW5bC/czo7CHews3Em5tdzjeIQ5ggmxExgfM55xMeOufFSZokDFmaYQlLsDig617DitM0HscHWyxcZmMv9OGNEmhOhWJBwJIdpNURSs339P9XerqF71nUcNEXo9fhMnEHjd9QTMnIEuMLDN65yrP8fOwp1q7VAr/YZ89b6MiR7D+JjxTIidQEpQypVNwGirVWuF8nY2BKJdYClpeV5ArFojlDBWDULRQ0HfCXMdCSG6NQlHQogLUhQF67FjVH33HdXfrfLoQ4TBoAaia+eogSg4uNVr1Npr2VO8hx2FO9heuJ1j5cc8jus1etIi0tw1Q0MjhmLQtq+DdisFhvJTTXML5e1svVZIa1DnFoof0xSGpOO0EAIJR0KIViguF/UHD1K9ejVVq1djP9M0PF1jMOA3ZQqBc67Ff/r0VmuIbE4b+0v3s7NoJzsKd3Cw9CAOxeFxTv+Q/oyLGcf4mPGMihp1+UPs66vUtcfyd6uBKG9X632FAmKaglD8GIgZBgaZMFEI0ZKEIyEEoM5UXbt7N9Wr11C9Zg2O4mL3MY3JhP/UKQRcOwf/6ekt+hA5XU6OnjvKjsId7Cjcwb6SfdQ76z3OifCJYUjoKIaEjmZw6CiCjCGXvpyFywmlWU0hKG+3+vj8EWQ6o7r2WGMQih+j1gp11tpoQoheRcKREFcxl9WKZetWqlevoWbdOpwVFe5jWl9f/NPTCbh2Nv5TpqD1a6rZcSkuvi//nl1Fu9hRtIM9RXuotld7XDvYFMrZsgTsNX1xWFKptodyEvgPAIcBdSHUdc+mtx2QqgqbaoTy96j9hs6fVwggOKkhBI1Wb6WvkBDiCkg4EuIq46yspGbjRqrXrMWyaROu2lr3MV1wMP4zZxAwezZ+Eya45yFSFIVTlafYWbiTnUU72VW0q8WIMn+DP6OjRzM+Zjxjo8dSb4lg3qItFyyL1eGi3GJTw5HNAoX7G4JQQyCqym/5JKM/xI5oqhGKHy0jyIQQHUrCkRBXAXt+PtVr11G9bh21u3eDo6n/jz4qioDZswmYPRvfUSPR6PUoikJedR47z+xkV/EudhbupLSu1OOaZr2ZkZEjGRM9hnEx4xgYOhC9tulXyqHayjbLo8NJP00+w7QniN34H6g4pE6wqLg8T9RoIfIaiBulhqC40RAxQOYVEkJ0KglHQvRCiqJgzcpSA9HatViPHvU4buqXiv+MmQTMnIHPkCGg0ZBfk8+qU1+zq2gXO4t2Ulxb7PEco9bI8Mjh7jA0JGwIBl17RpQpxGvKGKY5QZr2BMO1JxiqOYWvpmFix+YD1wLjGmqFGprHYoaD6cJzJAkhREeTcCREL+GyWqndsYPq9eupydjgsbArWi3mkSMIaAhEhsRE8mvy+a54N7u2LmVX0S4KLYUe19Nr9aSFpzE6ejTjoscxLHIYJl07+vHUlELBXiKztvK+YT1p2pOEaapbnFatmDngSqH/qHQiBkxUa4cCY670YxBCiCsm4UiIHsxRWkrNhg1UZ2Rg2boNpVn/IY2PD36TJhEwYwZ+6dMoNFjYULyb3TnvsmvXLoosRR7X0mv0DI0Yyuio0YyNGcuwiGGY9RcZSVZfqfYTKtjXMJx+n3tV+kggsqH1y6boyFISOeBKIVNJJdPVlxNKLApalo+ZTERcUEd+LEIIcUUkHAnRgzQ2l9VkZFC9PoP6Awc8jusjI/FPT8c/PZ3SwdFsrzzErqLt7Fn7FiV1njNC6zV6BocPZkz0GMZEj2F4xPALr1Fms0DRwYYgtFe9Pft9KydqILwf5SFD+cMRf/a7+pKlJGLF2AGfgBBCdD4JR0J0c86aGixbt1KzcSOWDRtxlHp2jPYZMgS/9GmcG5XCzqBz7CnZy57iBZz79pzHeY3NZKOiRjEmegzDIoa1HYbs9VB8SA1ABZnqbenRlh2mAYISIW6E2lcobpTaT8gnkPz8Sj48uLljPgQhhOhCEo6E6GYURcF24gQ1GzZSs3EjtXv2eIwu05jNmMePpWr0APb307Hdfpy9JZ9QfdyzX49RayQtQu0zNDpqNGkRaa03kzms5wWhTDUIuRwtz/WPhriRahCKHakuyuoX3ur7CPEzYtJrsTpaCVQNTHotIX5SoySE6F4kHAnRDThraqjdvp2azZuxbNyEvaDA47i+TxKW0QPI6u9LRlgJeyv2UufYAqeazjHrzYyIHMGoqFGMihrFkPAhLTtQ2+uh+DAU7mvoK5SpDqFvLQj5hqvhJ2Z4QyAaeUkdpuOCzax7Np1yi63Ncy55hmwhhOgCPTocbdy4kTfeeIM9e/ZQWFjIsmXLmD9/fpvnZ2RkMH369Bb7CwsLiY6O7sSSCuFJcbmoP3wYy+bN1GzeQl1mJjibLYxqNFKflsqJQYGsj69iq+YkDiVfXSWjTD0l0BjIyKiRjI4azaioUS3mGVJXot/ZFIIKM6HkaMsFWAHMoWoQim1oHosZ3iHLbcQFmyX8CCF6nB4djiwWC8OGDeO///u/ufXWW9v9vGPHjhHYbLHMyEiZXVd0PntxCZYtW7Bs3oxl61aPpToAbLFhnBkUwuaEWtaGF2MzHG86qECUbxQjo0YyMnIkI6NGkhqcilajVY/XV0LOdjUIFR5Qb8uOtd5HyDdMDT8xw9QtdgQEJ8q6Y0II0aBHh6PrrruO66677pKfFxkZSXBwcMcXSIhmnNXV1O7ahWXbdizbtmLLPuFx3OFrJKdfMNsS6tiaUEtpcCXQOKu0htTgVEZEjnAHolj/WPVQTSkU7YdDXzcFofJTtMovsqlpLGaYej8wToKQEEJcQI8OR5dr+PDhWK1WhgwZwoIFC5g0adIFz7darVitVvfjqqqqzi6i6IFcNht1mZnUbt+OZes26g4e9GgqUzRQmODHzgQre5IVsmOdOHXqiDKD1siI8CEMjxzOiIgRjIgcQbAxUA09RQdhx1/V28IDUFPUegGCEiEmTQ1B0Q23MqmiEEJcsqsqHMXExPDOO+8wevRorFYrf/3rX0lPT2fHjh2MHDmyzectXLiQl156qQtLKnoCxeGg/uhRanfswLJjJ7W7d6PU1XmcUxqmZ2+ik0N9NBxK0mAxqyE7xBTK1MjhjIhUg9A1gSkYz52AokNw6Bsoek29b2s5szRoIKxvU7NYYxDyDe2Cdy2EEL2fRlEUxduF6AgajeaiHbJbM23aNBITE/noo4/aPKe1mqOEhAQqKys9+i6J3k1xOqk/cpTanTvVbc8eXDU1HudU+mk4kASH+mg42EdDWZDafJUanMqwiGEMixjGCP8kkmrOoSk5rAagooNQdrz1jtI6E0RdA9FD1RAUnQZRg2W9MSGEuAxVVVUEBQVd9Pv7qqo5as3YsWPZvPnCE9WZTCZMpnasKSV6FbVmKIva3bvV2qHdu1HOC0M1PnA0QcPhRDUM5UaAvzGAtIg0bg8bwjBjGENtDgLPZsOJA7D1E6gpbv0FzaEQPaQpBEUPhfB+0K7FXYUQQnSUqz4cZWZmEhMj/TKuFvkVdW3Pu1NfT8DpY/hmHaR2z15qM/dBrWczmcXUEIaS1ECUE6khOaQvQ4NT+W99MMNsDvqW56PNPgTbvgCXvZUXamgWix4KUUMaaoWGQkCMdJQWQohuoEeHo5qaGrKzs92PT506RWZmJqGhoSQmJvLcc8+Rn5/Phx9+CMCbb75JcnIygwcPpr6+nr/+9a+sW7eOVatWeestiC6UX1HHjN9muGdsDrBZuObsaYacPcmQc9mklhdQqyjUNntOjQ8ci2sKQ1WJwQwJ7EOazp+brTaGVBQRcPwA1K5v/UVNgWozWNQQtXksaqh6a/Tr/DcshBDisvTocLR7926PSR2feeYZAO6//36WLFlCYWEhOTk57uM2m42f/vSn5Ofn4+vrS1paGmvWrGl1YkjRuygOB+X7DzHz+80Mqshm0LnTxFW3HHV4NgCy4jUcTdBwIslAQEIMQ3RmptTV8di5POKOH0TDwZYvoNFCaF+1WcwdhgZDUILUBgkhRA/Tazpkd6X2dugS3uMoL6duXyZlu7dQuWcXuqxT6K0tm7jywiArQcPReC3FkT4EmzTMMtgZWVFAX5u19b8e/KPV2p/Ihi3qGogYCAaZCVoIIboz6ZAtrhrOmhpqDx6kaM9mKvfvQZt1EnNp0xD4xq70tUbIjtVwPBZKo/QYwpykaaqYYqvnEZsdn1oFjzY1UxBEDlSDT9TghiA0WIbMCyFELyfhSPQortpaKo8cIH/XBqoO7EOXdYqAoio0DfWfzXvy5IfC93EaKmI0GELrifGvZ7DdxjybHV9FgYaBZzWKD1lKCsdd8RxT4rl33lz6DBoFgbHSJCaEEFchCUei27KWFJO7bxOl+3diPZqF6WQBgSUWtArogJBm55YGwqkYqI5woQuzEx5kpb/GxiSbHXNDy7HTFkhmfRJfu+L4XokjW4klW4knXwlDQeu+1i0JkyEoqGvfrBBCiG5DwpHwOoe1nvzDOyk6uJPqrMPw/SkCz5wloNoBQPB555f7QV40WMKd6MPsRAZb6aezMsnuQAdqn6Dwa9Q5giIGQsQAiBjI0UoffrBoSxe/OyGEED2NhCPRZWz1teQc3UHhwZ3UHDsCp3Lxyz1LaJkNnQKBqFsjF1AUCuciXDhCnfgG24kOrKe/zsZEdOpcQeH9ILw/hDXchqeCTxu1PlWVre8XQgghmpFwJDqUw+kgP/cIRUf3UJ59BNvJbLQ5RfiXVBN21oneBeGoW3MWE5SGKdSFuNCF2gkJstHHv55BwbFoQ/tCWKoahsJSITQFgpNAd2n/fUP8jJj0Wvc8R60x6bWE+Bkv/Y0LIYToNSQciUtWVVtOwamDlJ44TPXJI9hyzqApKMFcYiHsrANfa8taoEZ1RjgbqlAf6kQX5CAo0ElsTBCp0ckYwlIgNLkhAPWFkD5g8OmwcscFm1n3bHrbM2SjBqi4YBmSL4QQVzMJR8KDzWGjpOwUZ08dofz0EWrzTmItKIDicxjP1RJQ7iC0CnRK6zVAoDaHVQQp1AQpuIKcGIO1BEcGEB0fz4DkwWjDUiAkWQ1CQQldunZYXLBZwo8QQogLknB0Fah31HOu+CSVuVlUF52itiSf+nMl2MrLcZ6rQFNei7HKjm+1k6Aa8LGDEYi6wDXtOqgKUKgLUHAFaTCE+BAYE0ZUn2RiUkegj0qB4EQ1/PjIRJlCCCF6DglH3YTL6eT0ti9w2utx2qw47XU47DZcdhtOmxWXw4rTbsfpsOGyWrFbqrHXWnDW1uGqs6LU21CsdrRWB1qrC53VhanOhbkO/OvUmh6AgIbtYmpNYPFXsPlpcAUZ0If54hcVTmhCH6JShhCQPBRNcDwERIPedPELCiGEED2EhKNuQnE5sT78ovuxFrX2piPVG6DODFYzOMw6XL56NEE+GEMD8Y2IJDAmkfA+AwlJHoYuPAHMITIJohBCiKuOhKNuQmcwUuMDigZcWs9N0YJLq2m6r9PgMmpxmbQoJj0akwGtjwGt2Qed2ReDny9G/0B8w6IIjEoiODaVgLgB6IIiQau9eGGEEEKIq5iEo25kTOZRbxdBCCGEuOpJNYIQQgghRDMSjoQQQgghmpFwJIQQQgjRjIQjIYQQQohmJBwJIYQQQjQjo9Uug6KoMypWVVV5uSRCCCGEaK/G7+3G7/G2SDi6DNXV1QAkJCR4uSRCCCGEuFTV1dUEBQW1eVyjXCw+iRZcLhcFBQUEBASg6cAZpKuqqkhISCA3N5fAQFmPrLPJ59215PPuWvJ5dy35vLvW5X7eiqJQXV1NbGws2gtMiiw1R5dBq9USHx/fadcPDAyUH64uJJ9315LPu2vJ59215PPuWpfzeV+oxqiRdMgWQgghhGhGwpEQQgghRDMSjroRk8nEiy++iMlk8nZRrgryeXct+by7lnzeXUs+767V2Z+3dMgWQgghhGhGao6EEEIIIZqRcCSEEEII0YyEIyGEEEKIZiQcCSGEEEI0I+Gom3jrrbfo06cPPj4+jBs3jp07d3q7SL3Wxo0bmTdvHrGxsWg0Gr766itvF6lXW7hwIWPGjCEgIIDIyEjmz5/PsWPHvF2sXmvx4sWkpaW5J8ebMGEC33zzjbeLdVV49dVX0Wg0PP30094uSq+1YMECNBqNxzZw4MAOfx0JR93A0qVLeeaZZ3jxxRfZu3cvw4YNY86cOZSUlHi7aL2SxWJh2LBhvPXWW94uylVhw4YNPPbYY2zfvp3Vq1djt9u59tprsVgs3i5arxQfH8+rr77Knj172L17NzNmzODmm2/m8OHD3i5ar7Zr1y7effdd0tLSvF2UXm/w4MEUFha6t82bN3f4a8hQ/m5g3LhxjBkzhkWLFgHq2m0JCQk88cQT/OIXv/By6Xo3jUbDsmXLmD9/vreLctUoLS0lMjKSDRs2MHXqVG8X56oQGhrKG2+8wUMPPeTtovRKNTU1jBw5krfffpvf/OY3DB8+nDfffNPbxeqVFixYwFdffUVmZmanvo7UHHmZzWZjz549zJo1y71Pq9Uya9Ystm3b5sWSCdE5KisrAfULW3Qup9PJp59+isViYcKECd4uTq/12GOPccMNN3j8Hhed5/vvvyc2NpaUlBTuuececnJyOvw1ZOFZLysrK8PpdBIVFeWxPyoqiqysLC+VSojO4XK5ePrpp5k0aRJDhgzxdnF6rYMHDzJhwgTq6+vx9/dn2bJlXHPNNd4uVq/06aefsnfvXnbt2uXtolwVxo0bx5IlSxgwYACFhYW89NJLTJkyhUOHDhEQENBhryPhSAjRZR577DEOHTrUKX0ERJMBAwaQmZlJZWUl//rXv7j//vvZsGGDBKQOlpuby1NPPcXq1avx8fHxdnGuCtddd537flpaGuPGjSMpKYnPPvusQ5uNJRx5WXh4ODqdjuLiYo/9xcXFREdHe6lUQnS8xx9/nOXLl7Nx40bi4+O9XZxezWg0kpqaCsCoUaPYtWsXf/zjH3n33Xe9XLLeZc+ePZSUlDBy5Ej3PqfTycaNG1m0aBFWqxWdTufFEvZ+wcHB9O/fn+zs7A69rvQ58jKj0cioUaNYu3ate5/L5WLt2rXSR0D0Coqi8Pjjj7Ns2TLWrVtHcnKyt4t01XG5XFitVm8Xo9eZOXMmBw8eJDMz072NHj2ae+65h8zMTAlGXaCmpoYTJ04QExPTodeVmqNu4JlnnuH+++9n9OjRjB07ljfffBOLxcKDDz7o7aL1SjU1NR5/ZZw6dYrMzExCQ0NJTEz0Ysl6p8cee4xPPvmEf//73wQEBFBUVARAUFAQZrPZy6XrfZ577jmuu+46EhMTqa6u5pNPPiEjI4PvvvvO20XrdQICAlr0nfPz8yMsLEz61HWSZ599lnnz5pGUlERBQQEvvvgiOp2Ou+++u0NfR8JRN3DnnXdSWlrKCy+8QFFREcOHD+fbb79t0UlbdIzdu3czffp09+NnnnkGgPvvv58lS5Z4qVS91+LFiwFIT0/32P/+++/zwAMPdH2BermSkhLuu+8+CgsLCQoKIi0tje+++47Zs2d7u2hCXLG8vDzuvvtuzp49S0REBJMnT2b79u1ERER06OvIPEdCCCGEEM1InyMhhBBCiGYkHAkhhBBCNCPhSAghhBCiGQlHQgghhBDNSDgSQgghhGhGwpEQQgghRDMSjoQQQgghmpFwJIQQQgjRjIQjIYQQQohmJBwJIYQQQjQja6tdBpfLRUFBAQEBAWg0Gm8XRwghhBDtoCgK1dXVxMbGotW2XT8k4egyFBQUkJCQ4O1iCCGEEOIy5ObmEh8f3+ZxCUeXISAgAFA/3MDAQC+XRgghhBDtUVVVRUJCgvt7vC0Sji5DY1NaYGCghCMhhBCih7lYlxjpkC2EEEII0YyEIyGEEEKIZiQcCSGEEEI0I32OhBBCCOF1LosFR2kp9pISXNXVBMyc6bWySDgSQgghRKdxWSzYS0pwlJTiKCluuC3BUVqqbg33XbW1TU8yGBh4YL/X5hKUcCSEEEKIS6bY7WpNT3ExjuKShuBTogah4hI19JSU4LJY2n1Nra8v+shI9BERKLW1aPz8OvEdtE3CkRBCCCE8uGpr1dBTVIS9qBhHcVFTCCoqwl5SgvPsWVCUdl1P6+/vDj36yEj0kRHq/YgIDI37IyLQeikMnU/CkRBCCHEVcdXXYy8sVENOYRH2okIcRcXYi4twFKohyFVZ2b6LGQwYIiLQR0Whj4rCEBWJPjKqIQCpIcgQGdltQk97STgSQgghegnF4VCbugoKsBcUNoSgwoYQVISjsBBnRUW7rqXx9cUQHY0hOgp9VDT6qEgM0dFq+Gm4rwsJQXOBNcp6KglHQgghRA/hslgagk8B9sJCNQA13i8swFFcAk7nRa+j9fVFHxOjhp3oKAzRMe5bQ3QU+uhotP7+V+3i6hKOhBBCiG5AURRcVVXY8/Ox5efjaAxBBQXq4/wCnO1p7tLr1RqfmBgMsTHoo2MwxESjb9wXHY02MPCqDT7tIeFICCGE6CLOhvDTuNny8j0eu2pqLnoNbWAghtjYpq0hBBliYtDHxKIPD0Oj03XBu+m9JBwJIYQQHcRltapBJy8PW24u9ryG+/l52PPycVVVXfQaurAwDHFxTeEnrjEIxWGIi0Xn798F7+TqJuFICCGEaCdFUXCWlWHLzcOem9Nwm4stT711lJRc9Bru8BMXizEuruF+HIb4eAwxMWjN5i54J+JCJBwJIYQQzSh2u9rPJycHW04O9pxctRaoIQQpdXUXfL7W1xdDQgKG+HiM8XEY4uIxJMRjjI/HEBeH1te3i96JuFwSjoQQQlx1XPX1atjJycF2Jgdbbg72MzlqCCoouPCIL40GfUw0xoTEhtCTgDFRDUOGhAR0wcHS2bmH6/bh6K233uKNN96gqKiIYcOG8ec//5mxY8e2eu6XX37JK6+8QnZ2Nna7nX79+vHTn/6Ue++9133OAw88wAcffODxvDlz5vDtt9926vsQQgjRtVx1ddhyc7GdOYM9Jwfb6TPu2iBHUdEFZ3fW+PhgTIjHkJiEMSFBDUGJier92Fg0RmMXvhPR1bp1OFq6dCnPPPMM77zzDuPGjePNN99kzpw5HDt2jMjIyBbnh4aG8stf/pKBAwdiNBpZvnw5Dz74IJGRkcyZM8d93ty5c3n//ffdj00mU5e8HyGEEB1Lsdmw5eVjO30a25kzHreOoqILPlfr748xKQljUiKGhEQ1/CQmYEhMQh8R3isnNxTto1GUdi6M4gXjxo1jzJgxLFq0CACXy0VCQgJPPPEEv/jFL9p1jZEjR3LDDTfw61//GlBrjioqKvjqq68uu1xVVVUEBQVRWVlJYGDgZV9HCCHExSkuF47iYmynT2M9dQrbqdPuEGTPywOXq83nagMD1QCUmIgxKRFjUhKGxESMffpI89dVqL3f39225shms7Fnzx6ee+459z6tVsusWbPYtm3bRZ+vKArr1q3j2LFjvPbaax7HMjIyiIyMJCQkhBkzZvCb3/yGsLCwNq9ltVqxWq3ux1XtGIophBDi0jhrLNhOncJ2+hS2U6fUIHRarQW6UCdoja+vGoD6JDXc9nHfSgASl6PbhqOysjKcTidRUVEe+6OiosjKymrzeZWVlcTFxWG1WtHpdLz99tvMnj3bfXzu3LnceuutJCcnc+LECZ5//nmuu+46tm3bhq6NSbMWLlzISy+91DFvTAghrmKKy4WjqAjrSTUA2U6dVO+fPHnhYfB6Pcb4eIzJyWr46dMHY3IfjEl90EdGSAASHarbhqPLFRAQQGZmJjU1Naxdu5ZnnnmGlJQU0tPTAbjrrrvc5w4dOpS0tDT69u1LRkYGM2fObPWazz33HM8884z7cVVVFQkJCZ36PoQQoidz2WzYz5zBeuIk1pMnsJ08pd6eunAtkC48HFOfPmoIagxCyX0wxsejMRi68B2Iq1m3DUfh4eHodDqKi4s99hcXFxMdHd3m87RaLampqQAMHz6co0ePsnDhQnc4Ol9KSgrh4eFkZ2e3GY5MJpN02hZCiFaoTWEnsWafwHbyBNYTJ7GdOIEtL6/t4fAGA8bEREwpyRiTUzAmJzfcT0Yn/ThFN9Btw5HRaGTUqFGsXbuW+fPnA2qH7LVr1/L444+3+zoul8ujv9D58vLyOHv2LDExMVdaZCGE6LWcVVVYT5zAduIE1uwTWE+cwHoiG0dBYZvP0fr7Y+ybgik5Rb3t21etDUpIQKPvtl8/QnTfcATwzDPPcP/99zN69GjGjh3Lm2++icVi4cEHHwTgvvvuIy4ujoULFwJq36DRo0fTt29frFYrK1eu5KOPPmLx4sUA1NTU8NJLL/GDH/yA6OhoTpw4wc9+9jNSU1M9hvoLIcTVylldjTU7G2t2NrbsbKzfq/cv1B9IFx6OKSUFU2pfjCl9MfVNwZjSV/oCiR6rW4ejO++8k9LSUl544QWKiooYPnw43377rbuTdk5ODtpm81BYLBYeffRR8vLyMJvNDBw4kI8//pg777wTAJ1Ox4EDB/jggw+oqKggNjaWa6+9ll//+tfSbCaEuKq4amvV2p/j32P9/nt3ILrQ3ED66GhMffs2haDUvhhTUtCHhHRhyYXofN16nqPuSuY5EkL0FIrNhvX0aTUANQah77/Hnpvb5nP0UVGYUlPVrV8qxr59MaWmogsI6MKSC9Hxevw8R0IIIdpPURQcBQXUHz+O9dhxrMePY/3+ONZTp8HhaPU5uojwhhDUryEI9cOU2lc6RYurnoQjIYToYZzV1Wr4OX6c+mPH1Bqh48dx1dS0er7W318NPv37q7f9+mHql4o+NLSLSy5EzyDhSAghuinF6cSWk4P12HHqj2WpNUJZWeqq8a0xGDAlJ2MaMABT/3749O+PqX9/9NHR0jFaiEsg4UgIIboBZ00N1mPHqD+apd4eO4b1++/bnDBRHx2NaUB/fPoPUGuE+vfHlNxHVosXogNIOBJCiC6kKAqOoiLqj2ZRn3UUa9Yx6rOysOfktHq+xmRSm8EGDsBnwEA1EA0YgC4oqItLLsTVQ8KREEJ0EsXhwHbqFPVHj1J/5Cj1WVlYjx7FWVnZ6vn66Gh8BgzANHAgPgMHYBowAGNSEpo21n0UQnQOCUdCCNEBXHV1agfpxiB09CjW48dRWpuhX6/HlJKCz6CBmAYOUoPQwIEyX5AQ3YSEIyGEuETOmhqsR49Sf+SIe7OeOAkuV4tztb6+ak3QoEH4XDMIn0GDMKamopW+QUJ0WxKOhBDiApwVFWpt0OHDahA6fATbmTOtnqsLC1NDULMgZEhMRNNsJn8hRPcn4UgIIRo4KyqoP3KEukOH1TB0+DD2vLxWz9XHxOBzzTVqCLrmGnyuuQZ9ZKQMmReiF5BwJIS4Kjmrqqg/fJi6Q4eoP3yE+kOH2gxChoQEfAYPdocgn2sGyQSKQvRiEo6EEL2ey2Kh/uhR6g4eov7QIeoOHcR+pvWh84bERHwGX4N58GB8hgzB55prZDkNIa4yEo6EEL2Ky2bDmpVF3cGD1B88RP3hQ212ljbEx+MzZAjmIYPdNUMyf5AQQsKREKLHUpxObKdOUXfgIPWHDqq3x46B3d7iXH10ND5DBmMeMgSfIUPxGXyNDJ0XQrRKwpEQokdQFAVHcTF1+w9Qf/BAQyA6hKu2tsW5upAQfIYOwTxkKD5Dh+AzeDCGyEgvlFoI0RNJOBJCdEvOGovaP+jAAeoO7Kd+/wEcpaUtztOYzWofoaFpmNOG4jN0KIa4OBk1JoS4bBKOhBBepzidWLOzqdu/n7r9+6k/cABr9glQFM8TdTpM/ftjHjoU87A0fIYMxdQ3BY1efpUJITqO/EYRQnQ5R1mZWiOU2RCGDh5stXlMHxuDOW1YUxi65hq0vr5eKLEQ4moi4UgI0akUm436Y8eo25fprhlqbT4hra8vPmlpmNPS1CA0dKj0ExJCeIWEIyFEh7KXlFCXmanWCmVmUn/4cMvFVzUaTKl98Rk2DPOwYZjThmFK7SurzwshugUJR0KIy6bY7dRnHaNu376GQJSJvaCgxXm6oCDMw4djHq6GIZ+hQ9EFBHihxEIIcXESjoQQ7eYoL1ebx/btU7dDh1Dq6z1P0mox9evXEIbUQGTs00dGjwkhegwJR0KIVikulzqCrFkYam01em1QEObhw/BtCEM+Q9PQ+ft5ocRCCNExJBwJIQBw1dZSd+Agdfv2UrtvH3WZ+3FVVbU4z9i3L+YRw/EdMQLz8OEYk5PRaLVeKLEQQnQOCUdCXKXsRUXU7dtH7d591O3dS31WFjidHudofH3VYfQjR6hhKC0NXXCwdwoshBBdRMKREFcBxenE+v331O7dS11DGGqt47Q+JkYNQSNGYB45Ap8BA2SCRSHEVUd+6wnRC6lNZAeo3bNHDUP79+OqqfE8SafDZ8AAzCNH4jtSDUSGmBjvFFgIIboRCUdC9AL24hK1r9CevW02kWn9/NTRYyNH4DtyJOa0NLR+0nFaCCHOJ+FIiB5Gcbmwfp+thqG9e6nbsxd7fn6L8/SxMfiOGOkOQ6b+/WWSRSGEaAcJR0J0c666uqZRZHv3tj6KTKvFNGCAWiPUEIakiezqk19RR7nF1ubxED8jccHmLiyRED2ThCMhuhl7SYnaT2jfXmr37qP+6FFwODzO0fj6Yh6WptYMjRqJedgwdP7+Xiqx6A7yK+qY8dsMrA5Xm+eY9FrWPZsuAUmIi5BwJIQXKU5nw0SL+9RaoX2Z2HNzW5ynj4zEPGpkQzPZSHwG9oxRZFKT0XXKLbYLBiMAq8NFucUmn7kQF9H9f7sK0Ys4ayzUH9ivTrLY1igyjQZT//5NHadHjMQQF9vjlt+QmgwhRE8l4UiITqIoCvacHOoyM6ltWKXeeuwYuDzDgtbXV12QdXjD/ELDh/WKRVmlJkMI0VNJOBKig7gsFuoOHlJXp9+/n7rMTJzl5S3OM8TGYh45Ul2CY+RITP369YgmMiGEuFrIb2QhLoPidGI7eVIdRXbgAHUHDrRaK6QxGPAZPLhphfoRwzFERXmp1EIIIdpDwpEQ7WAvLqHuwH7qG8JQ/aFDuCyWFufpY2KaVqgfNgzTNdegNRq9UGIhhBCXS8KREOdxnD1L/eHD1B8+TN3hw9QfPISjuLjFeRqzGfPgwfgMS8M8NE1qhYQQopeQcCSuai2C0KHDOIqKWp6o1WJKTcU8LA2foUPVWqG+faWvkOg2QvyMmPTai44ODPGTmkwhLkZ+s4urguJwYDt9mvqsY1iPHaP+WBbWY8dbrRFCo8HYpw8+gwc3bNdgHjxY1iET3VpcsJl1z6bLvFJCdAAJR6JXURQFR2kptpMnsR4/roahrCys2dkotla+NM4LQuYhgzENGiSzTXcAqcnoenHBZgk/QnQACUeiR1KcTux5eVhPnMR28gTWEyexnjyB7eQpXNXVrT5H6+uLqX9/TAMH4DNgAKYBAzH174/OX2qEOoPUZAgheioJR6Lbctls2PPyseflYsvNxZ6Tiy0vD3tODrYzZ1qvCQLQajEmJGDsl4rPgIHuMGSIj6egysrZ5l/WlQ6orHQ/lC/rjiU1GUKInqjbh6O33nqLN954g6KiIoYNG8af//xnxo4d2+q5X375Ja+88grZ2dnY7Xb69evHT3/6U+699173OYqi8OKLL/Lee+9RUVHBpEmTWLx4Mf369euqtyQAxeXCefYs9pISHMUlOEqKcZSUYC8swp6Xhy03V+0PpChtXkNjMmFMTsaUkoKxbwqmvn0xpqRg7NOn1eHzspyFEEKI9ujW4Wjp0qU888wzvPPOO4wbN44333yTOXPmcOzYMSIjI1ucHxoayi9/+UsGDhyI0Whk+fLlPPjgg0RGRjJnzhwAXn/9df70pz/xwQcfkJyczP/93/8xZ84cjhw5go+PT1e/xY7ncoK9rmGrRbHX4bDVYLPVUG+twmavwWqrwWqvxWavxea0YnPasSt27E4bNpcdu8vZcGvHrjhxKAourQYHGlwaDU5Nwy0anBpwoUFBi8GmxVQHRisYaxVMdS6MtU6MtXYMNfWYqurxqajD55wFY0UtWueFl5YAwOyDJi4GfXwcxoQETElJmBP64JOSgiE2Fo1O1+6PRpazEEII0R4aRbnAn+ZeNm7cOMaMGcOiRYsAcLlcJCQk8MQTT/CLX/yiXdcYOXIkN9xwA7/+9a9RFIXY2Fh++tOf8uyzzwJQWVlJVFQUS5Ys4a677mrXNauqqggKCqKyspLAwMDLe3PnUZxOVv4kHQUXLkXBpVFQUHDiQkHBhbrP1bDPgQu7Rmm61SjYNaibtmHTaHAAygUWLNW6QO8CnRN0SsOtC/ROBZ1LvW+0g48dTHbwsYHJrqiPbep+Hzv4WEF7ie/ZBVT4wzl/KA/QcM4fzgVoKAmG4mANxSFQbQZaKb+fwQ8/gx/+Bn/8Df7qfaO/e5+fwY9AYyBBpiD3bUmFlp98mIXiMoPS9t8Fy5+YzJC4oEt8N0IIIbq79n5/d9uaI5vNxp49e3juuefc+7RaLbNmzWLbtm0Xfb6iKKxbt45jx47x2muvAXDq1CmKioqYNWuW+7ygoCDGjRvHtm3b2gxHVqsVq9XqflxVVXW5b6ttGg0pG8s6+KLKebedz2EEmw/YTS5sPgp2k6Le+kC9WcHip1Djr1DjD9V+mqYQp9FgB+waLVaDCY3eSLBWh1mjoQ4X9S4HdsXhfh2L3YLFbqGEkksqn39/9VZxGVGcZhSnL4rDH8Xph+LwQ3H6sy6/nBJnLKE+oe7Nz+CH5gIhUwghRO/RbcNRWVkZTqeTqPNmHI6KiiIrK6vN51VWVhIXF4fVakWn0/H2228ze/ZsAIoaJvdr7ZpFrU3812DhwoW89NJLl/tW2kWj1XJiYggaQIsGLRo0ila9RYMWLVpQ7ykatBodOnQNt1q06NVbjQ4tOvU8jQ6tRotWo0Gn0aHRqNdyUxTQ69G4N13DY4PHY63JB62vGY2vL1qzL1pfX7S+ZrRmM1pfXzRmMzp/f3RBQWia9/VxOsBaBdZqdasrh5piqCk577Zhs5RxoSDnAKzmIOqCE6kNSaImJAFLUCw1fhHUmPywOGqpsddgsVuotlVTbaum0lZJlbWKSmsl5+orqLZVo9EoaLQ2NFobGCpbvM57R7+Do577fHQ+hJnDiDBHEG4O97jfuEX6RhLqE4pO2/6mPiGEEN3PJYejTZs2MWXKFLZs2cKkSZM6o0xXJCAggMzMTGpqali7di3PPPMMKSkppKenX/Y1n3vuOZ555hn346qqKhISEjqgtJ5u/PvWDr+mV+n04Buqbu3hdKghqeIMlJ9u2Jru62uK0NdV4ld3EAoPej5X7wPh/SFyEEQMhMhREDUEguLdzXKH8iu58c8bQVuPRleHRlerbvoaNDoLGp0Frb6Gcakm7FRzrv4c5fXl1DpqqXfWk1+TT35N/oXfskbnDkrNtyjfKCJ8I4j2jSbKLwqzXvo0CSFEd3XJ4eibb75Br9ezYsWKTg1H4eHh6HQ6is+bwbi4uJjo6Og2n6fVaklNTQVg+PDhHD16lIULF5Kenu5+XnFxMTExMR7XHD58eJvXNJlMmEymK3g3ol10egiKU7ekiS2P22qhIgfOnYTSLHUrOQplx8FRD0UH1K053zCIGQ6xwwk09ScWKwWuMBSXL4o9rNVi/M9tnn2O6hx1lNWVcbbuLGV1ZZTVlVFaV+r5uLaUsvoynIqT4tpiimtbmXm7mWBTMNF+0e6wFO0X7X4c4x9DpG8kBq3hUj9BIYQQHeCSwtFLL72Ew+FgxowZPPnkk/zqV7/ihRde6JSCGY1GRo0axdq1a5k/fz6gdsheu3Ytjz/+eLuv43K53P2FkpOTiY6OZu3ate4wVFVVxY4dO/jJT37S0W9BdDSjL0QOVLeB1zftdznV2qXGsNT8tvYsnFgLJ9aSCGz1gbNKAIdcyRxUkjngSmGvqz9ltN0B26w3kxCQQELAhWsLnS4nZ+vPUlJbQnFtMSW1JR5bcW0xxZZiah21VFgrqLBWkHWu9SZirUZLhDmCWP9YYvxi3LfN7/safC/jQxRCCHExlzxa7b333qOyspLg4GAefvjhzioXoA7lv//++3n33XcZO3Ysb775Jp999hlZWVlERUVx3333ERcXx8KFCwG1b9Do0aPp27cvVquVlStX8otf/ILFixe7y/raa6/x6quvegzlP3DgwCUN5e+M0WqiE9jroeQwFGRCYSa23L1oSo5i0DhbnHraFcUepR8HGMCj991DVN/h0Al9hxRFodpeTZGlyGMrri2myFJEoaWQIksRdpf9otcK9Qkl1i+WuIA4Yv1jifePJ9Y/Vt38YvHR94KpKYQQogN12mg1h8PBs88+y7vvvntFBWyPO++8k9LSUl544QWKiooYPnw43377rbtDdU5ODlpt0wByi8XCo48+Sl5eHmazmYEDB/Lxxx9z5513us/52c9+hsVi4Uc/+hEVFRVMnjyZb7/9tnfMcSQ8GXwgbpS6AUYgv6yc+ryD+JQewFx2EN/STEznjtFHW0wfivkBm+EffwNjAMSPhsTxkDAWEsaB8cqXGdFoNAQaAwk0BtI/pH+r57gUF2frzlJoKaTAUkBhTSEFNQUUWYrcj6vtap+oc/XnOHT2UKvXCTeHE+8fT1xAnHrrH0d8QDzx/vFE+kZKx3EhhGhDt57nqLuSmqNepq4C8ndD7k7I3QF5u8FW43mOVq+GrD6T1a2DwtLlqrJVUVBToHYSr86nwFLg7jCeX51PraP2gs/Xa/XE+cepgck/3t1sGB8QT3xAPH4GWW9OCNH7tPf7W8LRZZBw1Mu5nFByBHK2q4HpzFaoyvM8p5uFpeYURaHSWkl+TT55NXnqbXXTbYGlAIfLccFrhPqEumuZGoNTQkACiYGJhPmEyZxPQqAuSSQLK/csXRqOcnJy2LhxIyaTiREjRrhHi/VWEo6uMoqiTi9werO6ndrUSlgyqM1vKemQMh1iR6ij77ohp8tJSW0JeTV55FXnkVud676fV51HubX8gs83683EB8ST4K+GpcYap8SARGL8YqS5TlwVZK3GnqnLwtGf/vQnnnnmGXx9fdFoNNTU1DBq1Cj++te/kpaWdiWX7rYkHF3l2hOWTEGQPAX6TlfDUmhKq8ugdEfVtmo1KNWowalxy6vOo9BSiEtp+8tAr9W7a5sag1NiQCKJgYnE+sfK9ASi11DnTdt80fNkOaLupcvCUUREBP/v//0/fvWrX6HVasnOzubtt9/m3Xff5bvvvmPy5MlXcvluScKR8KAoUH4KTmbAifVwagPUnzfzdlAi9E2HvjPV2iVzcNeXswPYnXbya/LJrc4lpzqHvOo8cqpzyKnKIb8m/4Kj7HQaHTF+MSQGJroDU1Jgklrz5B+PQSfBSfQcEo56pi4LR2azmYMHD7ZoSvvNb37D119/zY4dO67k8t2ShCNxQS6nOn3AyXVwIkPt5N08NGh0ahNc6ix1i04D7aUu29v9OF3qBJiNYSm3Opecqhx3iKp31rf5XK1GS4xfjDssJQUmkRSYRGJAInH+cRKcRLcj4ahn6rJwNG7cOB5//HHuvfdej/3Hjx9n2LBh1NXVXcnluyUJR+KS2Cxqp+4T6yB7jTqjd3N+kQ1BaSb0ndH+5VZ6EEVRKKktIadaDU1nqs543NY52v490Vjj1Dw4NdY6SVOd8BYJRz1Tp81zdL7f/e533HrrrRiNRu644w73KJYdO3bQr1+/K728ED2f0Q/6zVY3FqrrxWWvUbeTG8BSAvs/UTeNFuJGQ/856hY1pMf0VboQjUZDlF8UUX5RjIke43FMURTK6so8AlNj7VNOdQ51jjq1w3hNXovr6jQ64vzjSAhMICmgKTQlBSQR4x+DXts9O8ULIbq3DhmttnLlSn784x9TX1/P8OHDsdlsHDp0iI8++ojrrruuI8rZrUjNkegwDivkbFOD0vdroPSo5/HAOOh3LfSfC8lT1SVUriKKolBaV+oZnKpyOFN9htyq3As21TV2Dm/s49S8xinaN1pG1YkrIjVHPVOXz3NktVpZu3YtGRkZ7Nu3j3379lFeXk5oaChDhw4lLS2NN998syNeyuskHIlOU5kH36+C49+ptUrNm5v0PmpA6netWqsUnOi9cnYDLsVFSW2JZ2hqVutkc7U9/4xRa3SPqHOHpoaap0jfSLSant8HTHQuCUc9U7eYBDI3N5fMzEx3WFq2bFlnvVSXknAkuoS9Tp0m4Pvv4PgqqMzxPB41BAZcB/2vU+dV6gWdujuKS3FRbCnmTHWz0NRY41Sde8FJMH10Pi2a6RprnsLN4TIBpgBknqOeqluEo95KwpHocooCJUcbgtJ36gi45vMN+Uc19FO6Tp0q4CprfrsUDpeDQkuhR01TY3jKr8nHqbRcmLiRr9631Wa6xIBEQn1CJThdZWSG7J5HwlEnknAkvK72nNr8duwbyF4LtuqmY3ofNSANuE7tqxQQ7bVi9jR2l53CmsIWoelM1RkKLAUXnADT3+DfYg6nxvshphAJTkJ0AxKOOpGEI9GtOGxwZjMc+1YNS+c3v8WNVoPSwBsgYmCvGP3mDXannbyavFZrnAothSi0/as0wBDQ6uSXSYFJBJuCJTgJ0UUkHHUiCUei21IUddHcYyvVoJS/x/N4SB8YcIMalhIndNv133oaq9OqzhbeMP1A8z5ORZaiCz43wBjQam1TYkCiBCchOpiEo04k4Uj0GNVFakg69o26vInT2nTMJ1jtpzTgenUCSlOAt0rZq9U76smrzvPoHN44wq64tviCzw0wBLg7hycENq1TlxCQQJhPmAQnIS6RhKNOJOFI9EjWGji5Xg1Kx7+F2rNNx3RGSJ4GA69Xw5L0U+oSdY66ljVODVMRXCw4NXYOb764b+P9CN8ImY5AiFZIOOpEEo5Ej+dyQu5OOLYCslbAuZOex+NGNwSlGyBigPRT8oLGGqfmS67kVOeQW5V70T5OPjof4gPiSQxIdM/nlBCQQEJAAtF+0TJzuLhqSTjqRBKORK+iKFB6rCEorYT83Z7HQ1PU2qSBN0DCOJCZpb3O5rSRV5NHblWu50K/1TkU1BRccDoCvVavLrnSEJYaa5sSAhKIC4jDpDN14TsRomtJOOpEEo5Er+bup7SyoZ9Ss3lcfMPU6QEGXK8ukivzKXU7jdMRNNY4NQanxs3usrf5XA0aIn0jPWqa4gPi3fcDjfL7TvRsEo46kYQjcdWwVqvzKB1bqU4+WV/RdEzvAynT1ea3/teBf4TXiinax+lyupdcyalW+znlVee5Q1Sto/aCzw8yBZHg7xmYGu/LsiuiJ5Bw1IkkHImrktOuLpKbtVJtgqtoPp+SBhLGqjVKA66HiP5eK6a4PIqicK7+nLuGqTE0NW5n689e8PlGrZFY/1h3WIr3jyc+oGHzj8fXILWMwvskHHUiCUfiqqcoUHxYrVHKWgGFmZ7Hw1KbglLCWOmn1AvU2mvVfk7NglPjbUFNAQ6l7fXqAEJ9Qt1BqfltnH8cUb5R6OT/iOgCEo46kYQjIc5Tmdc0n9KpjdC8X4tHP6XpYPTzXjlFp3C4HBTXFpNXnaduzUJUXk0eldbKCz5fr9UT6xdLnH+cOzDFBcSpAco/niBTkMzpJDqEhKNOJOFIiAuor4ITa9Xmt++/g/pmX4x6H3U+pcZ13wJjvFdO0WWqbFXuWqb8mnzyq/PJq8lT79fk43BduNbJz+CnBqaGzR2gGjZpshPtJeGoE0k4EqKdLthPCYgdodYo9Z8L0UNlPqWrkNPlpLSu1KOmqTFA5dfkU1pXetFrhJhCiPOPI9Y/lriAOOL81JqnWP9YYv1i8dH7dME7ET2BhKNOJOFIiMugKFBytGHk27eQtxuaT2QYGA8D5qq1Sn2mgF7m2xHqZJgFlgLyqj1DU+NWZau66DXCfMLc4SnWX22+i/GLUW/9YzDrzV3wTkR3IOGoE0k4EqIDVBerzW7HvoUT68BR13TM4Kf2T+o/V13/zT/Se+UU3Vq1rZqCmgLyavIoqClwh6bG+xa75aLXCPUJJdYvlhj/GGL9Yt0hKsYvhlj/WAKMsu5gbyHhqBNJOBKig9nr4OSGpvmUas5byT5ulDqXUv850vwm2k1RFKpsVRTUFLjDUoGlwB2eCmoKqLHXXPQ6AcYANTz5xbgDVLR/tDtIhfqEyhxPPYSEo04k4UiITuRyQdF+NSQd/xYK9nkeD4xTQ1K/OZA8VWbpFlek0lpJoaWQgpoCCi2F5NfkU1hTSIFFDU8V1oqLXsOoNRLtF02MfwwxfjHq/fNupemue5Bw1IkkHAnRhaoK4ftValg6uR7szWZx1vuo/ZP6z4F+10JIkvfKKXqlWnutR3g6/35JbQkuxXXR64SYQoj2i/YITM3vh5vDZUHgLiDhqBNJOBLCS+x1cHqzWqN0fBVUnjf6LXwA9L9WrVVKHA86g3fKKa4adped0tpSd2AqshS5g1ORpYiCmoKLLssCoNPoiPCNINo32h2cov2iifKNct+X5rsrJ+GoE0k4EqIbUBQozVJrlL5frU4Z0Hw1elMQ9E2H1NmQOkvmVBJeoSgK1fZqCmuaglPz2+LaYootxRedYRzUyTKjfKPUzS+KaN9oovyiPPaF+YTJbOMXIOGoE0k4EqIbqqtQR719v0oNS7VlnsejhkK/2eoWPxZ00oQhugeny8nZ+rMeoal5cCqqLaKsrqxdzXeNNVCNgSnSN5Jov2gifSPdW5RvFEadsQveWfcj4agTSTgSoptzuaBgrxqSsldD/l485lSSWiXRw9hdds7WnXUHJ3d4atwsxZTVleFsXnt6ASGmkKaw5BdFpFm93xisInwjCDGF9LplWyQcdSIJR0L0MJayplql7LVQd87zeORgSJ0BfWdC4gQwyIzKoudprIEqtqiBqchSREltiTtAldSWUFJbgtVpbdf1DFoDEeYId2iK9I30eBxhjiDCN4IAQ0CPCVESjjqRhCMhejCXU50eoK1aJb0Z+kyG1JlqWArvJ/MqiV6jce6nxuDk3urU29LaUoprizlXf+7iF2vgo/Nxh6VwczgRvuptpG+k+tisHusOCwhLOOpEEo6E6EUsZ9UpAk6sU2uVzp+AMihRrVVKma7Oq+Qb6p1yCtGF7E47ZXVl7tDUGJxK60o97rdn+ZZGBq3BHZbCzGHuMBXuG064T1OoCvMJw9BJI00lHHUiCUdC9FKKAsWH4cRaNSjlbAOnrem4RqsulpsyXV3eJH4s6K/Ojq1CgLr2XWldKaW1pZTUlVBWW+Z+XFpXSlmd+rjSWnlJ1w02BfPdD77D19Cxk7xKOOpEEo6EuErYLOq8SicaapbKjnkeN/hBn0nQt6FmKWKANMEJ0Qqb0+YOSmV1ZZTVllFWX0ZpbSln6866j52tO4tDceCj82HnPTs7vBlOwlEnknAkxFWqMh9OZqjNcCczwFLqedw/Wm16S5kGydMgOMEbpRSix3IpLiqtlVRYK0gOSu7w60s46kQSjoQQuFxQclitUTqxXm2Cc9R7nhOaooaklGnQZyr4hXmnrEIIQMJRp5JwJIRowV4PeTvh5AY4tUEdBXf+nDPRQ9Ww1GcKJE0AnyDvlFWIq1R7v7+7/SItb731Fn369MHHx4dx48axc+fONs997733mDJlCiEhIYSEhDBr1qwW5z/wwANoNBqPbe7cuZ39NoQQvZ3BR21Sm/l/8PAa+PkpuOufMO7HEDFIPafoIGxbBP+8E17rA3+ZDqtfUKcVsFZ7tfhCiCbduuZo6dKl3HfffbzzzjuMGzeON998k88//5xjx44RGRnZ4vx77rmHSZMmMXHiRHx8fHjttddYtmwZhw8fJi4uDlDDUXFxMe+//777eSaTiZCQkHaXS2qOhBCXrLoYTm+CUxvV23MnPY9rdBA3Uq1V6jMZEsaByd87ZRWil+oVzWrjxo1jzJgxLFq0CACXy0VCQgJPPPEEv/jFLy76fKfTSUhICIsWLeK+++4D1HBUUVHBV199ddnlknAkhLhilflqSDq9CU5tgooznsc1OogdDkmTmsKSOdgbJRWi12jv93e3XXnRZrOxZ88ennvuOfc+rVbLrFmz2LZtW7uuUVtbi91uJzTUc9K2jIwMIiMjCQkJYcaMGfzmN78hLKztjpJWqxWrtWm69aqq9k96JYQQrQqKg2F3qRtARY4akk5vgtNboDIH8veo29Y/ARq1z1LSJHX6gMSJ0sFbiE7SbcNRWVkZTqeTqKgoj/1RUVFkZWW16xo///nPiY2NZdasWe59c+fO5dZbbyU5OZkTJ07w/PPPc91117Ft2zZ0Ol2r11m4cCEvvfTS5b8ZIYS4mOBEGHGPuoEals5sVedZOrMVzp2AogPqtmOxek74AEgcr64HlzQBgpNkniUhOkC3DUdX6tVXX+XTTz8lIyMDH5+mRSTvuusu9/2hQ4eSlpZG3759ycjIYObMma1e67nnnuOZZ55xP66qqiIhQeYvEUJ0ouBEdWusWaoqhDNb1O30FnVCysZt7wfqOQExDWFponobNRi0rf/RJ4RoW7cNR+Hh4eh0OoqLiz32FxcXEx0dfcHn/va3v+XVV19lzZo1pKWlXfDclJQUwsPDyc7ObjMcmUwmTCbTpb0BIYToSIExMPQ2dQN1TbjcHZCzFXK2q4vpVhfC4WXqBmAMgPjRan+lxHEQNxp8pJ+kEBfTbcOR0Whk1KhRrF27lvnz5wNqh+y1a9fy+OOPt/m8119/nZdffpnvvvuO0aNHX/R18vLyOHv2LDExMR1VdCGE6Hx+YTDwenUDsNWq/ZNytqsTUubuBFt1w2ze6xuepFFrkxLGQsJ49TakjzTFCXGebj1abenSpdx///28++67jB07ljfffJPPPvuMrKwsoqKiuO+++/j/27v3mKbOPg7g37Yg8EKpIrcWkKmbF1SoIDKnOOechhgyky1xzkzclv0zXGZ4t8VliWBmosuyxHljxj/m/jHqlsCSbeqMy7wsGhGsUnxVvCO0XBVKlWv7/vGUcirgBNueevh+kienPT09/fUk0q/Pec5zEhISsGXLFgDAN998g40bN2L//v2YP3++ez8RERGIiIhAe3s7Nm3ahLfeegvx8fG4ceMGvvjiC9hsNlRWVj517xCvViOigNfbAzRcFr1LNefE8vEr4gAgPBZIzBQ9TImZ4sa6nEKAFOq5v1oNAFauXInGxkZs3LgRVqsVRqMRR44ccQ/Svnv3LtTq/nksi4uL0dXVhbfffttjP4WFhSgqKoJGo8GlS5fw008/4cGDBzAYDFi6dCm+/vprnjYjImXRBAH6VNHmfiTW2az9QanmHGAxAfYG4OrvogGASg3EznCFJVdgGv8SoA74OYOJvCage44CFXuOiEgRujsAy0Wg9jxwrwy4dx5orRm4XYgOSJgNJGSIZkgXY6CInjOKmAQyUDEcEZFitVk8w1JtBdDzaOB2WoOY0Tsh3RWYZvNecRTwGI58iOGIiEaN3h6goUqEpNpysWz8H+B0DNx2/IuA3iiCksEI6NOAEK2/KyYaEsORDzEcEdGo1mV3nY4r7w9Mgw32hgqIfskzMMWncsA3yYbhyIcYjoiIHmNvAupMgOWCWNaZgLZ7g2yoAsZPFiFJnyYGjMen8VYo5BcMRz7EcERE9BTaG8UVcXUmMUll3QXAVjf4tpEJIizFu66wi5spZgjnHEzkRQxHPsRwREQ0Qu2NgPUiYHHdJ85ySdw3bjAhOjFpZfwsIH6mCEyx04HgMP/WTIrBcORDDEdERF7U0QbUm0VQslwErJVA4xXA0T1wW5VazLsUP1MEp9gZQFwKoEtiLxP9K4YjH2I4IiLysZ4uoOmaCE3WStHqzcDD5sG3D9GJXqW4FCA2pb+XKWysX8umwMZw5EMMR0REMnA6xSzffYGp4TJQfxlougo4egZ/j9YgQlLsdCBmmms5lVMMjFIMRz7EcEREFEB6uoDmahGUGqrEsr5qiKvlXHRJ/YEpZpoITNFTgFD+TVcyhiMfYjgiInoOdLQCDVfEpJXSZbt16Pdo9a6gNBWImeJaTgPCozmmSQEYjnyI4YiI6Dn2sEUM+G74n1g2XgEarz05NIWOFRNajn9JLPseR00Cgsb4rXR6NgxHPsRwRESkQI8eAE3VYgxT41UxILzxCnD/DoAhfipVGmBccn9oGj8ZiJosbqWi1QNqtT+/Af0LhiMfYjgiIhpFuh8BLTdFWGq6LpbN1SJIdbUP/b6gMFdYmiTCkjs4TQbCY3iaTgZP+/sd5MeaiIiInj/BYWJOpbgZnuv7rp5rrnYFphtA83WxvH8b6HkkrqyrNw+yz3ARmqImupaSx1oDe5xkxnBEREQ0EioVEKkXbeJCz9d6u4EHd/vDUktfcLoJtNYA3XagvlK0x2lCxKm6cS8M3saE+/qbjXoMR0RERN6mCRanz8ZPHvhaT6cITi03Xe1W/+MHd4DeTtcpvGuD7zs8VoSnscmu5QTxeOwEMUUBB4g/M4YjIiIifwoK6b/i7XG9PWJ+pvu3B7aWW0DHA8DeINq9skF2rgIiDZLAlCQC09gkQDcB0CXw3nRPgeGIiIgoUGiC+k+fDebRfXH13P3bopfpwV3R7rse9zwC2mpFu3tm8H2Ex0gCU19LAHSJQGQi53QCwxEREdHzI2ycaAbjwNecTsDe5ApNd0Rgaq0BHtT0L7vtgL1RtLqKwT9DEyLCUqQrMOkSxePIBNErFWkQNSg4QDEcERERKYFKBUTEiJY4Z+DrTqfoeXo8MLXWiJ6m1lqgvV6MeeobAzWUoLD+oORuCYA2Xlxtp40HIuJET9hz6PmsmoiIiIZHpQL+EyWaPm3wbXq6AFudCEqt98T4p77HtjqgrQ542CxO37W4rsIb+gOBiFgxGaZWLwJTpEGEpr7wpI0Xp/nUGp985ZFiOCIiIiIhaMyTxzwBQHcHYLOIoNRWJ3qdbBZXgLKK1m4FHD2iJ6q9HrCYht6fSi0C0uOhKfu/sg0eZzgiIiKipxcc6pqwcuLQ2zgcwMMmV4iyiKXNKnqfbPUiPNnqxVV3Tkd/iLJeEu9XqYFFX/rn+wyC4YiIiIi8S60Wp9QiYoc+hQcAjl4xONxmFeGob9lpk/VUG8MRERERyUOtcQ3ijpe7Eg+8eQsRERGRBMMRERERkQTDEREREZEEwxERERGRBAdkj4DT6QQAtLW1yVwJERERPa2+3+2+3/GhMByNgM1mAwAkJSXJXAkRERENl81mg06nG/J1lfPf4hMN4HA4UFdXB61WC5UXb7zX1taGpKQk1NTUIDIy0mv7pcHxePsXj7d/8Xj7F4+3f430eDudTthsNhgMBqjVQ48sYs/RCKjVaiQmJvps/5GRkfzH5Uc83v7F4+1fPN7+xePtXyM53k/qMerDAdlEREREEgxHRERERBIMRwEkJCQEhYWFCAkJkbuUUYHH2794vP2Lx9u/eLz9y9fHmwOyiYiIiCTYc0REREQkwXBEREREJMFwRERERCTBcEREREQkwXAUIHbt2oUXXngBoaGhyMrKwrlz5+QuSbFOnjyJ3NxcGAwGqFQqlJaWyl2Som3ZsgWZmZnQarWIjY3FihUrcPXqVbnLUqzi4mKkpqa6J8ebN28eDh8+LHdZo8LWrVuhUqmwfv16uUtRrKKiIqhUKo82bdo0r38Ow1EAOHjwIAoKClBYWIiKigqkpaVh2bJlaGhokLs0RbLb7UhLS8OuXbvkLmVUOHHiBPLz83H27FkcO3YM3d3dWLp0Kex2u9ylKVJiYiK2bt2K8vJynD9/HosXL8abb76JqqoquUtTtLKyMuzZswepqalyl6J4M2bMgMVicbfTp097/TN4KX8AyMrKQmZmJnbu3AlA3LstKSkJn3zyCTZs2CBzdcqmUqlQUlKCFStWyF3KqNHY2IjY2FicOHECCxculLucUSEqKgrffvstPvzwQ7lLUaT29nakp6dj9+7d2Lx5M4xGI7Zt2yZ3WYpUVFSE0tJSmEwmn34Oe45k1tXVhfLycixZssS9Tq1WY8mSJThz5oyMlRH5RmtrKwDxg02+1dvbiwMHDsBut2PevHlyl6NY+fn5WL58ucffcfKd6upqGAwGTJo0CatXr8bdu3e9/hm88azMmpqa0Nvbi7i4OI/1cXFxuHLlikxVEfmGw+HA+vXrMX/+fMycOVPuchSrsrIS8+bNQ0dHByIiIlBSUoKUlBS5y1KkAwcOoKKiAmVlZXKXMipkZWVh3759mDp1KiwWCzZt2oTs7GyYzWZotVqvfQ7DERH5TX5+Psxms0/GCFC/qVOnwmQyobW1Fb/88gvy8vJw4sQJBiQvq6mpwaeffopjx44hNDRU7nJGhZycHPfj1NRUZGVlITk5GYcOHfLqaWOGI5lFR0dDo9Ggvr7eY319fT3i4+NlqorI+9atW4fffvsNJ0+eRGJiotzlKNqYMWPw4osvAgAyMjJQVlaG77//Hnv27JG5MmUpLy9HQ0MD0tPT3et6e3tx8uRJ7Ny5E52dndBoNDJWqHxjx47FlClTcP36da/ul2OOZDZmzBhkZGTg+PHj7nUOhwPHjx/nGAFSBKfTiXXr1qGkpAR//fUXJk6cKHdJo47D4UBnZ6fcZSjO66+/jsrKSphMJnebM2cOVq9eDZPJxGDkB+3t7bhx4wb0er1X98ueowBQUFCAvLw8zJkzB3PnzsW2bdtgt9vx/vvvy12aIrW3t3v8L+PWrVswmUyIiorChAkTZKxMmfLz87F//378+uuv0Gq1sFqtAACdToewsDCZq1OeL7/8Ejk5OZgwYQJsNhv279+Pv//+G0ePHpW7NMXRarUDxs6Fh4dj/PjxHFPnI5999hlyc3ORnJyMuro6FBYWQqPRYNWqVV79HIajALBy5Uo0NjZi48aNsFqtMBqNOHLkyIBB2uQd58+fx2uvveZ+XlBQAADIy8vDvn37ZKpKuYqLiwEAixYt8lj/448/Yu3atf4vSOEaGhqwZs0aWCwW6HQ6pKam4ujRo3jjjTfkLo3omd27dw+rVq1Cc3MzYmJisGDBApw9exYxMTFe/RzOc0REREQkwTFHRERERBIMR0REREQSDEdEREREEgxHRERERBIMR0REREQSDEdEREREEgxHRERERBIMR0REREQSDEdEREREEgxHREQufTcMKCoq8nj+JDU1NVi0aBFSUlKQmpqKn3/+2ZclEpEf8PYhREQuu3fvRlBQEKqrq6HRaJCTk4NXX331ie+xWCyor6+H0WiE1WpFRkYGrl27hvDwcD9VTUTexp4jIlKc5uZmxMbG4vbt28N638cff4zW1lZs374dubm5A4LRO++8g++++85jnV6vh9FoBADEx8cjOjoaLS0tz1I+EcmMPUdEpDgFBQWw2WzYu3fvsN73ww8/QK1Wu3uOli9fjuzsbPfrZrMZCxcuxK1bt6DT6Qa8v7y8HHl5eTCbzc/8HYhIPgxHRKQoDx8+hF6vx9GjR/Hyyy8P671OpxMqlQpFRUUoKipyP5fKzMzE2rVrkZ+f77G+paUF2dnZ2Lt3L1555ZVn/h5EJB+eViMiRfnjjz8QEhLiEYxOnz6N4OBgdHR0uNfdvn0bKpUKd+7cca/rC0J9A7IfD0YAkJubiwMHDnis6+zsxIoVK7BhwwYGIyIFYDgiIkU5deoUMjIyPNaZTCZMnz4doaGh7nUXLlzAuHHjkJycPKz9z507F+fOnUNnZycA0du0du1aLF68GO+9996zfwEikh3DEREpyp07d2AwGDzWXbx4EbNnz/ZYZzKZkJaWNuz9GwwGdHV1wWq1AgD++ecfHDx4EKWlpTAajTAajaisrBz5FyAi2QXJXQARkTc9evTIo4cIEEHo3Xff9Vh34cIF91VmwxEWFgZAjG0CgAULFsDhcIysWCIKSOw5IiJFiY6Oxv37993Pe3t7YTabB/QcVVRUjCgc9V2mHxMT80x1ElHgYjgiIkWZPXs2Ll++7H5+9epVdHR0eJxqO3PmDGpra0cUjsxmMxITExEdHe2NcokoADEcEZGiLFu2DFVVVe7eI5PJBADYsWMHqqurcfjwYaxZswYA0NXVNez9nzp1CkuXLvVavUQUeBiOiEhRZs2ahfT0dBw6dAiACEfLli3DzZs3MWvWLHz11VfYtGkTIiMjsX379mHtu6OjA6Wlpfjoo498UToRBQhOAklEivP777/j888/h9lsRk5ODjIzM7F58+Zn3m9xcTFKSkrw559/eqFKIgpUvFqNiBRn+fLlqK6uRm1tLS5evIgPPvjAK/sNDg7Gjh07vLIvIgpc7DkiIsWyWq3Q6/WoqqpCSkqK3OUQ0XOC4YiIiIhIggOyiYiIiCQYjoiIiIgkGI6IiIiIJBiOiIiIiCQYjoiIiIgkGI6IiIiIJBiOiIiIiCQYjoiIiIgkGI6IiIiIJBiOiIiIiCQYjoiIiIgkGI6IiIiIJBiOiIiIiCT+D3OpGVP1mPaRAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ek = 100 # [K]\n", "sigma_m = 3e-10\n", " \n", "N_A = 6.022e23\n", "fig, (ax1, ax2) = plt.subplots(2, 1)\n", "\n", "# # From https://arxiv.org/pdf/mtrl-th/9501001.pdf which pulled from M. van Leeuwen and B. Smit, Phys. Rev. Lett. 71, 3991 (1993)\n", "# These data need to be rescaled according to Hentschke et al. (DOI: https://doi.org/10.1103/physreve.75.011506)\n", "# mustar2 = [2.5, 3.0, 3.5, 4.0]\n", "# T = [2.63, 3.35, 4.20, 5.07]\n", "# rho = [0.29, 0.25, 0.24, 0.24]\n", "# ax1.plot(mustar2, T, 'd')\n", "# ax2.plot(mustar2, rho, 'd')\n", "\n", "# Comparing with Hentschke, DOI: https://doi.org/10.1103/physreve.75.011506\n", "mustar2 = [1, 2, 3, 4]\n", "T = [1.41, 1.60, 1.82, 2.06]\n", "rho = [0.30, 0.31, 0.312, 0.289]\n", "ax1.plot(mustar2, T, 's')\n", "ax2.plot(mustar2, rho, 's')\n", "\n", "kB = 1.380649e-23 # Boltzmann's constant, J/K\n", "epsilon_0 = 8.8541878128e-12 # Vacuum permittivity\n", "\n", "for polar_model in ['GrossVrabec','GubbinsTwu+GubbinsTwu','GubbinsTwu+Luckas']:\n", " \n", " x = []; y = []; TT = []; DD = []\n", " rhostar_guess = 0.27\n", " Tstar_guess = 1.5\n", " for mustar2 in np.arange(0.001, 5, 0.1):\n", " z = np.array([1.0])\n", " mu2_C2m2 = 4.0*np.pi*epsilon_0*sigma_m**3*ek*kB*mustar2\n", " mu_Cm = mu2_C2m2**0.5\n", " model = teqp.make_model({\n", " \"kind\": 'SAFT-VR-Mie',\n", " \"model\": {\n", " \"polar_model\": polar_model,\n", " \"coeffs\": [{\n", " \"name\": \"Stockmayer\",\n", " \"BibTeXKey\": \"me\",\n", " \"m\": 1.0,\n", " \"epsilon_over_k\": ek, # [K]\n", " \"sigma_m\": sigma_m,\n", " \"lambda_r\": 12.0,\n", " \"lambda_a\": 6.0,\n", " \"mu_Cm\": mu_Cm,\n", " \"nmu\": 1.0\n", " }]\n", " }\n", " })\n", "\n", " T, rho = model.solve_pure_critical(Tstar_guess*ek, rhostar_guess/(N_A*sigma_m**3))\n", " # Store the values\n", " x.append(mustar2)\n", " TT.append(T/ek)\n", " DD.append(rho*N_A*sigma_m**3)\n", " # Update the guess for the next calculation\n", " Tstar_guess = TT[-1]\n", " rhostar_guess = DD[-1]\n", "\n", " ax1.plot(x, TT, label=polar_model)\n", " ax2.plot(x, DD)\n", " \n", "ax1.legend(loc='best')\n", "ax1.set(ylabel=r'$T^*$')\n", "ax2.set(ylabel=r'$\\rho^*$', xlabel=r'$(\\mu^*)^2$')\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.14" } }, "nbformat": 4, "nbformat_minor": 5 }