{ "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": "2024-03-15T22:41:46.624147Z", "iopub.status.busy": "2024-03-15T22:41:46.623701Z", "iopub.status.idle": "2024-03-15T22:41:46.642129Z", "shell.execute_reply": "2024-03-15T22:41:46.641681Z" } }, "outputs": [ { "data": { "text/plain": [ "'0.19.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": "2024-03-15T22:41:46.644141Z", "iopub.status.busy": "2024-03-15T22:41:46.643832Z", "iopub.status.idle": "2024-03-15T22:41:46.954880Z", "shell.execute_reply": "2024-03-15T22:41:46.954332Z" } }, "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": "2024-03-15T22:41:46.957498Z", "iopub.status.busy": "2024-03-15T22:41:46.957066Z", "iopub.status.idle": "2024-03-15T22:41:47.524735Z", "shell.execute_reply": "2024-03-15T22:41:47.524250Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkcAAAG4CAYAAABPb0OmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAACWLElEQVR4nOzdd3yV9f3//8fZOdl7DxLCEgh7r7AEB4rWWT+uj9pv6/5ZP221n49iW8XRYVsUrW3FUStaxVZAZYa9IeyAYWUvyDxJzrx+f1zJSQ5JIECSk4TX/Xa7buec67rOdd7nQHKeeU+NoigKQgghhBACAK23CyCEEEII0Z1IOBJCCCGEaEbCkRBCCCFEMxKOhBBCCCGakXAkhBBCCNGMhCMhhBBCiGYkHAkhhBBCNCPhSAghhBCiGb23C9ATuVwuCgoKCAgIQKPReLs4QgghhGgHRVGorq4mNjYWrbbt+iEJR5ehoKCAhIQEbxdDCCGEEJchNzeX+Pj4No9LOLoMAQEBgPrhBgYGerk0QgghhGiPqqoqEhIS3N/jbZFwdBkam9ICAwMlHAkhhBA9zMW6xEiHbCGEEEKIZiQcCSGEEEI0I81qncTpdGK3271dDCGuWgaDAZ1O5+1iCCEukaIouCor0QUHe60MEo46mKIoFBUVUVFR4e2iCHHVCw4OJjo6WqbcEKIHUBSFmg0bKFu8GJwu+nz+mdd+diUcdbDGYBQZGYmvr6/8UhbCCxRFoba2lpKSEgBiYmK8XCIhRFsUl4uadesoe3sx9UeOAKAxmbCdOoUpJcUrZZJw1IGcTqc7GIWFhXm7OEJc1cxmMwAlJSVERkZKE5sQ3YziclG9ahVli9/BeuwYABqzmZC77ybswQfQR0R4rWwSjjpQYx8jX19fL5dECAFNP4t2u13CkRDdhOJ0UrXyG8reeQfbiRMAaP38CLnnHkIfuB99aKiXSyjhqFNIU5oQ3YP8LArRfSh2O5XLV3D2nXewnTkDgDYggNB77yX0vnu92gH7fBKOhBBCCNFpXDYblcu+4ux772HPywNAFxRE6IMPEHLPPeguMlu1N0g4EkIIIUSHc9XXU/H5vzj7t7/hKCoCQBcaqoaiu3+Izt/PyyVsm4QjcdXr06cPTz/9NE8//bS3iyKEED2ey2Kh/NOlnH3/fZxlZQDoIyMJe/ghgm+/HW3DYInuTGbIFh6Kiop46qmnSE1NxcfHh6ioKCZNmsTixYupra3t8vLMmzePuXPntnps06ZNaDQaDhw40MWlEkIIcT5ndTVlixeTPWMmJW+8gbOsDENsLNELXqTv6lWE3ndfjwhGIDVHopmTJ08yadIkgoODeeWVVxg6dCgmk4mDBw/yl7/8hbi4OG666aYWz7Pb7RgMhk4p00MPPcQPfvAD8vLyiI+P9zj2/vvvM3r0aNLS0lo8T1EUnE4ner38FxdCiM7kKC/n3IcfUv7xP3BVVwNgSEok/Ef/j6Cb5qHppO+HziQ1R51NUcBm8c6mKJdU1EcffRS9Xs/u3bu54447GDRoECkpKdx8882sWLGCefPmAeoIoMWLF3PTTTfh5+fHyy+/DMDixYvp27cvRqORAQMG8NFHHzX7GBQWLFhAYmIiJpOJ2NhYnnzySffxt99+m379+rlrq2677TYAbrzxRiIiIliyZIlHWWtqavj888956KGHAMjIyECj0fDNN98watQoTCYTmzdv5sSJE9x8881ERUXh7+/PmDFjWLNmTYv3Xl1dzd13342fnx9xcXG89dZbHscrKip4+OGHiYiIIDAwkBkzZrB//36Pc77++mvGjBmDj48P4eHh3HLLLZf0+QshRE9iLymh+LXXyZ45i7OL38FVXY0xtS+xb7xB3xUrCP7BrT0yGIHUHHU+ey28Euud136+AIzt6/B29uxZVq1axSuvvIKfX+vPaT4sesGCBbz66qu8+eab6PV6li1bxlNPPcWbb77JrFmzWL58OQ8++CDx8fFMnz6dL774gj/84Q98+umnDB48mKKiIne42L17N08++SQfffQREydO5Ny5c2zatAkAvV7Pfffdx5IlS/jlL3/pLsPnn3+O0+nk7rvv9ijjL37xC37729+SkpJCSEgIubm5XH/99bz88suYTCY+/PBD5s2bx7Fjx0hMTHQ/74033uD555/npZde4rvvvuOpp56if//+zJ49G4Dbb78ds9nMN998Q1BQEO+++y4zZ87k+PHjhIaGsmLFCm655RZ++ctf8uGHH2Kz2Vi5cmU7/6GEEKLnsOfnc/Zvf6PiX1+g2GwAmK4ZRPiPf0zArFlotD2/3kWjKJdYvSCoqqoiKCiIyspKAgMD3fvr6+s5deoUycnJ+Pj4qDttlh4Rjnbs2MH48eP58ssvPWo8wsPDqa+vB+Cxxx7jtddeQ6PR8PTTT/OHP/zBfd6kSZMYPHgwf/nLX9z77rjjDiwWCytWrOD3v/897777LocOHWrRBPfll1/y4IMPkpeXR0ArQzqzsrIYNGgQ69evJz09HYCpU6eSlJTkrp3KyMhg+vTpfPXVV9x8880XfK9Dhgzhxz/+MY8//jigdsgeNGgQ33zzjfucu+66i6qqKlauXMnmzZu54YYbKCkpwWQyuc9JTU3lZz/7GT/60Y+YOHEiKSkpfPzxxxd8bdG1Wv2ZFEJcFuupU5x9769U/uc/4HAAYB4xgvCf/Bi/KVN6xLxibX1/n09qjjqbwVcNKd567Su0c+dOXC4X99xzD1ar1b1/9OjRHucdPXqUH/3oRx77Jk2axB//+EdArXl58803SUlJYe7cuVx//fXMmzcPvV7P7NmzSUpKch+bO3cut9xyi3t244EDBzJx4kT+/ve/k56eTnZ2Nps2beJXv/pVi/KeX66amhoWLFjAihUrKCwsxOFwUFdXR05Ojsd5EyZMaPH4zTffBGD//v3U1NS0WBKmrq6OEw2zu2ZmZvLII49c8LMUQoieqP7Ycc6++y5V334LLhcAvhPGE/7jn+A7dkyPCEWXSsJRZ9No2l17402pqaloNBqONaxv0yilYdE/83kjDNpqemtLQkICx44dY82aNaxevZpHH32UN954gw0bNhAQEMDevXvJyMhg1apVvPDCCyxYsIBdu3YR3DBj6kMPPcQTTzzBW2+9xfvvv0/fvn2ZNm1ai9c5v1zPPvssq1ev5re//S2pqamYzWZuu+02bA1Vwe1RU1NDTEwMGRkZLY41lu/8z0cIIXq6uv37KXvnXWrWr3fv809PJ/zH/w/z8OHeK1gX6PkNg6JDhIWFMXv2bBYtWoTFYrnk5w8aNIgtW7Z47NuyZQvXXHON+7HZbGbevHn86U9/IiMjg23btnHw4EFA7Vs0a9YsXn/9dQ4cOMDp06dZt26d+7l33HEHWq2WTz75hA8//JD//u//btdfK1u2bOGBBx7glltuYejQoURHR3P69OkW523fvr3F40GDBgEwcuRIioqK0Ov1pKamemzh4eEApKWlsXbt2vZ9WEII0U0pioJl+3bOPPAgp++8Sw1GGg0Bc+eSvOxLEt5Z3OuDEUjNkWjm7bffZtKkSYwePZoFCxaQlpaGVqtl165dZGVlMWrUqDaf+z//8z/ccccdjBgxglmzZvH111/z5ZdfukeGLVmyBKfTybhx4/D19eXjjz/GbDaTlJTE8uXLOXnyJFOnTiUkJISVK1ficrkYMGCA+/r+/v7ceeedPPfcc1RVVfHAAw+06z3169ePL7/8knnz5qHRaPi///s/XA3Vws1t2bKF119/nfnz57N69Wo+//xzVqxYAcCsWbOYMGEC8+fP5/XXX6d///4UFBS4O2GPHj2aF198kZkzZ9K3b1/uuusuHA4HK1eu5Oc///kl/AsIIYR3KIpCzfoMyt59h/r9DXPH6fUE3XQTYQ8/jCkl2bsF7GqKuGSVlZUKoFRWVnrsr6urU44cOaLU1dV5qWRXrqCgQHn88ceV5ORkxWAwKP7+/srYsWOVN954Q7FYLIqiKAqgLFu2rMVz3377bSUlJUUxGAxK//79lQ8//NB9bNmyZcq4ceOUwMBAxc/PTxk/fryyZs0aRVEUZdOmTcq0adOUkJAQxWw2K2lpacrSpUtbXH/r1q0KoFx//fUtjq1fv14BlPLyco/9p06dUqZPn66YzWYlISFBWbRokTJt2jTlqaeecp+TlJSkvPTSS8rtt9+u+Pr6KtHR0cof//hHj+tUVVUpTzzxhBIbG6sYDAYlISFBueeee5ScnBz3OV988YUyfPhwxWg0KuHh4cqtt9560c9bdK7e8DMpRGdyORxKxfLlyol5NylHBgxUjgwYqBxNG6YU/urXii0vz9vF63BtfX+fT0arXYZLGq0mhPAa+ZkUonUum43Kf/+bs3/9K/Yz6gAVrZ8fIT+8m9D770ff0GWgt5HRakIIIYTw4LJYKP/sc869/z6OkhIAdEFBhNx/H6H33IMuKMjLJeweJBwJIYQQvZyzooJzH/+D8o8+wllZCYA+KorQBx8g5Pbb0V7iCOTeTsKREEII0UvZi0s4t2QJ5UuXojQsHm5ISiTs4YcJuvlmtEajl0vYPUk4EkIIIXoZ25kznP3b36lctgzFbgfANHAg4T96hIA5c9DodF4uYfcm4UgIIYToJeqPHOHsX/9K1bffuWezNo8aRfj/+1GPWeKjO5BwJIQQQvRgiqJQu2sXZ9/7K5aGRbsB/KZNJfyRR/A9b1klcXESjoQQQogeSHG5qMnI4Oxf3qMuM1PdqdUSeN11hD3yMD4DB3q1fD2ZhCMhhBCiB1HsdqpWruTsX/+K9ftsADRGI0G33kLYQw9hTEjwcgl7PllbTXS4jIwMNBoNFRUVbZ6zZMkS96KtbVmwYAHDr4I1fDqSfPZC9F6u2lrOffgR2XPmUPDzX2D9Phutnx9hjzxM6to1xCxYIMGog0g4Eh6Kiop46qmnSE1NxcfHh6ioKCZNmsTixYupbRgG2lWeffbZDlvMVaPRXHBbsGBBh7zO5epOnztcXZ+9EN2do7yc0kVvkT1jJsWvvIKjoBBdWBgR/9//R+r6dUT+9KfoIyK8XcxeRZrVhNvJkyeZNGkSwcHBvPLKKwwdOhSTycTBgwf5y1/+QlxcHDfddFOXlcff3x9/f/8OuVZhYaH7/tKlS3nhhRc4duyYx2t1liVLlrBkyRIyMjJaPd7dPnfoPZ+9ED2ZPT+fs0s+oOJf/0KpqwPAkJBA2EP/TdD8+WhlSZxOIzVHwu3RRx9Fr9eze/du7rjjDgYNGkRKSgo333wzK1asYN68eZw+fRqNRkNmY+c/oKKiAo1G0+LLf8uWLaSlpeHj48P48eM5dOhQi9f86quv6NevHz4+PsyZM4fc3Fz3sfObdh544AHmz5/Pb3/7W2JiYggLC+Oxxx7D3jCHB8Dbb7/tvl5UVBS33XYbANHR0e4tKCgIjUZDdHQ0ZrOZuLg4srKyAHC5XISGhjJ+/Hj3NT/++GMSGqqqW2u2yszMRKPRcPr06Uv9yIH2fe6AfPad8NkL0R3VHz9Owc9/Tva1cyj/6COUujpM1wwi7ve/o+83Kwm56y4JRp1Mao46maIo1DnqvPLaZr253XNanD17llWrVvHKK6/g18Y08pc6P8b//M//8Mc//pHo6Gief/555s2bx/HjxzEYDADU1tby8ssv8+GHH2I0Gnn00Ue566672LJlS5vXXL9+PTExMaxfv57s7GzuvPNOhg8fziOPPMLu3bt58skn+eijj5g4cSLnzp1jU7Nhra0JCgpi+PDhZGRkMHr0aA4ePIhGo2Hfvn3U1NTg7+/Phg0bmDZt2iW99/bqjM8d5LMXoqdRFIW6PXs4+9e/UdPsjx3f8eMJe+Rh/CZOlDmKupCEo05W56hj3CfjvPLaO364A1+Db7vOzc7ORlEUBgwY4LE/PDyc+vp6AB577DF+8pOftPv1X3zxRWbPng3ABx98QHx8PMuWLeOOO+4AwG63s2jRIsaNG+c+Z9CgQezcuZOxY8e2es2QkBAWLVqETqdj4MCB3HDDDaxdu5ZHHnmEnJwc/Pz8uPHGGwkICCApKYkRI0ZctJzp6elkZGTw7LPPkpGRwezZs8nKymLz5s3MnTuXjIwMfvazn7X7fV+K9n7ur7322iVdVz57IXoGxeWieu1azv31b9Tt36/u1GgImD2bsEcexjx0qHcLeJWSZjVxQTt37iQzM5PBgwdjtVov6bkTJkxw3w8NDWXAgAEcPXrUvU+v1zNmzBj344EDBxIcHOxxzvkGDx6Mrtm09zExMZQ0rCw9e/ZskpKSSElJ4d577+Uf//hHuzozT5s2jc2bN+N0OtmwYQPp6enuL+2CggKys7NJT09v9/vOyclx99nx9/fnxz/+MZs2bfLY98orr1zwGlfyucPV+9kL0VO4bDbKP/+ck9ffQP4TT1K3fz8ao5HgO+6g7zcrif/THyUYeVGPrjlauHAhX375JVlZWZjNZiZOnMhrr73W4q/w5pYsWcKDDz7osc9kMrn/Su9oZr2ZHT/c0SnXbs9rt1dqaioajcajoyxASkqKei2zei2tVs3TiqK4z2ne76SzNTYLNdJoNLgapsgPCAhg7969ZGRksGrVKl544QUWLFjArl27Ljh0ferUqVRXV7N37142btzIK6+8QnR0NK+++irDhg0jNjaWfv36Ae17/7GxsR79gr788ku++OIL/vGPf7j3hYaGAu3/3Nv72p2pJ3z2QnR3zupqyj/9lHMffoiztAwAbUAAIXffTei9/yWjzrqJHl1ztGHDBh577DG2b9/O6tWrsdvtXHvttVgslgs+LzAwkMLCQvd25syZTiujRqPB1+Drle1S2qfDwsKYPXs2ixYtuuDnF9Hwg9t8BFLzINDc9u3b3ffLy8s5fvw4gwYNcu9zOBzs3r3b/fjYsWNUVFR4nHOp9Ho9s2bN4vXXX+fAgQOcPn2adevWXfA5wcHBpKWlsWjRIgwGAwMHDmTq1Kns27eP5cuXe/R5ac/71+v1pKamurfIyEjMZrPHvsZw1N7Pvb2v3ehq/eyF6K7sRUUUv/4G2enTKf3d73GWlqGPiiLy5z8ndf16Ip/5/yQYdSM9uubo22+/9Xi8ZMkSIiMj2bNnD1OnTm3zeY2jZYSnt99+m0mTJjF69GgWLFhAWloaWq2WXbt2kZWVxahRozCbzYwfP55XX32V5ORkSkpK+N///d9Wr/erX/2KsLAwoqKi+OUvf0l4eDjz5893HzcYDDzxxBP86U9/Qq/X8/jjjzN+/Pg2+7xczPLlyzl58iRTp04lJCSElStX4nK5LliT2Cg9PZ0///nP7hFWoaGhDBo0iKVLl/LWW2+5z0tNTSUhIYEFCxbw8ssvc/z4cX73u99dVnkbtedzB+Sz74TPXojOVn/8OOf+/j6VK1ZAQ02nMbUvYQ89TNAN16MxGr1cQtGaHl1zdL7KykqgqcmiLTU1NSQlJZGQkMDNN9/M4cOHL3i+1WqlqqrKY+uN+vbty759+5g1axbPPfccw4YNY/To0fz5z3/m2Wef5de//jUAf//733E4HIwaNYqnn36a3/zmN61e79VXX+Wpp55i1KhRFBUV8fXXX2Ns9ovA19eXn//85/zwhz9k0qRJ+Pv7s3Tp0ssuf3BwMF9++SUzZsxg0KBBvPPOO/zzn/9k8ODBF33utGnTcDqdHv1b0tPTW+wzGAz885//JCsri7S0NF577bU23397tfdzB/nsO/qzF6IzKIqCZcdOcn70I07ddDOVX30Fdju+o0cTv/htUv7zH4JvmS/BqBvTKM0b8Hswl8vFTTfdREVFBZs3b27zvG3btvH999+TlpZGZWUlv/3tb9m4cSOHDx8mPj6+1ecsWLCAl156qcX+yspKAgMD3Y/r6+s5deoUycnJ+MgcFEJ4nfxMiq6kOBxUr17N2b/9nfrGucUaR5499N+Yhw3zbgF7EKfLiU6ru/iJl6iqqoqgoKAW39/n6zXh6Cc/+QnffPMNmzdvbjPktMZutzNo0CDuvvtuj7/Qm7NarR4jhqqqqkhISJBwJEQ3Jz+Toiu4amupWLaMc0s+wN4wmarGZFIXgn3gAYxJSV4uYc9Qaa1kfe56Vp1eRYW1gk9u+KTDX6O94ahH9zlq9Pjjj7N8+XI2btx4ScEI1Kr6ESNGkJ2d3eY5JpMJk8l0pcUUQgjRizjKyij/5BPK//EJzoZuHbrgYEJ++ENC/use9Bfp4iHUQLQuZx2rzqxie+F2HC6H+1h+TT5x/nFeKVePDkeKovDEE0+wbNkyMjIySE5OvuRrOJ1ODh48yPXXX98JJRRCCNHbWE+e5Nz7S6j8979RbDZAXfMs9IH7Cb71VrTm9k+jcjVqDETfnfmOHQU7cChNgSg1OJVr+1zLtUnXei0YQQ8PR4899hiffPIJ//73vwkICKCoqAhQlyVonB/mvvvuIy4ujoULFwLqKJ7x48eTmppKRUUFb7zxBmfOnOHhhx/22vsQQgjRvbmX9/jb36lZv96932dYGmH//RABs2ai0XV8H5neory+3N1ktqPQMxD1C+nHtUnXcm2fa0kJSvFiKZv06HC0ePFigBYz6L7//vs88MADgDpbcePkcaDO+fLII49QVFRESEgIo0aNYuvWrVxzzTVdVWwhhBA9hOJwUL1mDWf//j71Bw6oOzUa/GfMIOy/H8Q8cqSsedaGc/XnWJuzltWnV7OzaCdOxek+1j+kvzsQJQddeqtPZ+s1HbK7UlsduqTzpxDdi/xMisvlrLFQ+eUXnPvgQ+z5+QBojEaCbrmF0Pvvx5TS/b7Qu4OyujLWnlnL6jOr2VW8C5fich8bFDqI2UmzmZ00mz5BfbxSvquqQ7YQQgjREezFxZR/9BHlSz/DVV0NgC4khJC77ybknh+iDwvzcgm7n5LaEtacWcPqM6vZU7wHhaY6l8Fhg92BKDEw0YulvDQSjoQQQlz16rOyOPf++1SuWAkOtT+MsU8fQh94gKD5N6OVmkcPRZYiVp9Zzeozq8ksyfQIREPDh3Jt0rXMSppFfMCljSDvLiQcCSGEuCopLheWzZs5+/771G5rWo/Qd8wYQh98EP/0aWi0vWohiSuSV53nriE6UHbA49jwiOHMSprF7KTZxPrHeqmEHUfCkehwGRkZTJ8+nfLy8jZXZF+yZAlPP/00FRUVbV5nwYIFfPXVV7K46BXq06cPTz/9NE8//bS3iyJEt+Cqr6fyP//h3AcfYjtxQt2p0xE4dy6hDzyAeegQ7xawGzldeZo1OWogOnL2iHu/Bg0jo0YyO2k2sxJnEeUX5cVSdjwJR8JDUVERCxcuZMWKFeTl5REUFERqair/9V//xf3334+vr2+XleXZZ5/liSee6JBrXWw0yYsvvsiCBQs65LUu1QMPPEBFRQVfffWVV15fiKuFOmnjPyn/5z9xlpcDoPXzI/j22wm9714MsT2/xuNKKYrCiYoTrM5Rm8y+L//efUyr0TImagyzk2YzM2km4eZwL5a0c0k4Em4nT55k0qRJBAcH88orrzB06FBMJhMHDx7kL3/5C3Fxcdx0001dVh5/f3/8/f075FqFhYXu+0uXLuWFF17g2LFjHq/VWZYsWcKSJUvIyMjotNcQQrSt/vhxzn3wAVX/+RrFbgfAEBtLyH33Enzbbeg68ee/J1AUhWPlx1h1ehVrctZwqvKU+5heo2dszFhmJc1iZuJMQn2ujlm/pTFVuD366KPo9Xp2797NHXfcwaBBg0hJSeHmm29mxYoVzJs3j9OnT6PRaDyauioqKtBoNC2+/Lds2UJaWho+Pj6MHz+eQ40LMTbz1Vdf0a9fP3x8fJgzZw65DesSgdqsNnz4cPfjBx54gPnz5/Pb3/6WmJgYwsLCeOyxx7A3/LIDePvtt93Xi4qK4rbbbgMgOjravQUFBaHRaIiOjsZsNhMXF0dWVhagLmAcGhrK+PHj3df8+OOPSUhIANQmQ41G49EcmJmZiUaj4fTp05f6kV/UkiVLWjRNfvXVVy1qwr7++mvGjBmDj48P4eHh3HLLLW1e869//SvBwcGsXbsWgN///vcMHToUPz8/EhISePTRR6mpqXGff+bMGebNm0dISAh+fn4MHjyYlStXdtybFKITKIpCzabN5Dz0MKduupnKL75EsdvxGZZG3Jt/oO+q7wh74IGrNhgpisKB0gP8fvfvuf7L67n969t57+B7nKo8hUFrYFr8NH496ddk3JnBu7Pf5fb+t181wQik5qjTKYqCUlfnldfWmM3tnpzs7NmzrFq1ildeeQU/P7/Wr3eJE539z//8D3/84x+Jjo7m+eefZ968eRw/fhyDwQBAbW0tL7/8Mh9++CFGo5FHH32Uu+66iy1btrR5zfXr1xMTE8P69evJzs7mzjvvZPjw4TzyyCPs3r2bJ598ko8++oiJEydy7tw5Nm3adMEyBgUFMXz4cDIyMhg9ejQHDx5Eo9Gwb98+ampq8Pf3Z8OGDUybNu2S3ntXWrFiBbfccgu//OUv+fDDD7HZbG2Gl9dff53XX3+dVatWMXbsWAC0Wi1/+tOfSE5O5uTJkzz66KP87Gc/4+233wbUmehtNhsbN27Ez8+PI0eOdGpNmxBXwt2f6MMPsWU39CfSagmYPZvQB+7Hd8QI7xbQi5wuJ/tK9rEmZw1rzqyhuLbYfcxH58PkuMnMSprFtPhp+Buv7p9xCUedTKmr49jIUV557QF796BpZx+h7OxsFEVhwIABHvvDw8Opr68H1C/Jn/zkJ+1+/RdffJHZs2cD8MEHHxAfH8+yZcu44447ALDb7SxatIhx48a5zxk0aBA7d+50f3GfLyQkhEWLFqHT6Rg4cCA33HADa9eu5ZFHHiEnJwc/Pz9uvPFGAgICSEpKYkQ7fhGmp6eTkZHBs88+S0ZGBrNnzyYrK4vNmzczd+5cMjIy+NnPftbu993VXn75Ze666y5eeukl975hw4a1OO/nP/85H330ERs2bGDw4MHu/c07avfp04ff/OY3/PjHP3aHo5ycHH7wgx8wdOhQAFJSusf0/kI0Zy8pofyTT6j4dCnOhppdra8vwbffRsi992K8xEXJewu7y86uwl2syVnD2py1nKs/5z7mq/dlWvw0ZiXNYnLcZHwNXdentLuTcCQuaOfOnbhcLu655x6sVuslPXfChAnu+6GhoQwYMICjR4+69+n1esaMGeN+PHDgQIKDgzl69Gib4Wjw4MHomq1fFBMTw8GDBwGYPXs2SUlJpKSkMHfuXObOncstt9xy0U7k06ZN429/+xtOp5MNGzZw7bXXEh0dTUZGBmlpaWRnZ7dYouZCcnJyPJajcTgc2O12j9qW559/nueff77d17yQzMxMHnnkkQue87vf/Q6LxcLu3btbhJs1a9awcOFCsrKyqKqqwuFwUF9fT21tLb6+vjz55JP85Cc/YdWqVcyaNYsf/OAHpKWldUjZhbhS9UeOcO6DD6hc+Q007090770E3/YDdAEBXi5h17M6rWwr2MbqM6vJyM2gylblPhZgDGB6wnRmJ81mQuwETDqT9wrajUk46mQas5kBe/d47bXbKzU1FY1G49FJGZpqCRoX8m1cp675qjPN+/x0tsYmuUYajQaXS52ePiAggL1795KRkcGqVat44YUXWLBgAbt27WpzSgGAqVOnUl1dzd69e9m4cSOvvPIK0dHRvPrqqwwbNozY2Fj69esHtO/9x8bGevTJ+vLLL/niiy/4xz/+4d4XGtq+tnutVsv5K/yc/3rmdvw7T5kyhRUrVvDZZ5/xi1/8wr3/9OnT3HjjjfzkJz/h5ZdfJjQ0lM2bN/PQQw9hs9nw9fXl4YcfZs6cOaxYsYJVq1axcOFCfve733XYSEIhLpXidFKzfj3nlnxA7e7d7v3mkSMJvf9+AmbOQKO/ur7eau21bMzfyNoza9mYt5FaR637WKhPKDMSZzA7cTZjYsZg0BoucCUBEo46nUajaXfTljeFhYUxe/ZsFi1axBNPPNFmv6OIiAhAHf3V2GTV1jxE27dvJzFRnS6+vLyc48ePM2jQIPdxh8PB7t273bVEx44do6KiwuOcS6XX65k1axazZs3ixRdfJDg4mHXr1nHrrbe2+Zzg4GDS0tJYtGgRBoOBgQMHEhkZyZ133sny5cs9+hs1f/8hISGtvn+9Xk9qaqr7cWRkJGaz2WNfe0VERFBdXY3FYnH/m5z/emlpaaxdu5YHH3ywzeuMHTuWxx9/nLlz56LX63n22WcB2LNnDy6Xi9/97nfu4PfZZ5+1eH5CQgI//vGP+fGPf8xzzz3He++9J+FIdDlnTQ2VX37JuY8+xt44eEOvV+cnuv8+zA1Nv1eLSmslG/I2sObMGrYWbMXqbKrdj/SNZFbiLGYlzWJk5Eh0Wt0FriTOJ+FIuL399ttMmjSJ0aNHs2DBAtLS0tBqtezatYusrCxGjRqF2Wxm/PjxvPrqqyQnJ1NSUsL//u//tnq9X/3qV4SFhREVFcUvf/lLwsPDmT9/vvu4wWDgiSee4E9/+hN6vZ7HH3+c8ePHt9mkdjHLly/n5MmTTJ06lZCQEFauXInL5WrRj6o16enp/PnPf3aPbgsNDWXQoEEsXbqUt956y31eamoqCQkJLFiwgJdffpnjx4/zu9/97rLK21xlZWWL0BMWFsa4cePw9fXl+eef58knn2THjh0sWbLE47wXX3yRmTNn0rdvX+666y4cDgcrV67k5z//ucd5EydOZOXKlVx33XXo9XqefvppUlNTsdvt/PnPf2bevHls2bKFd955x+N5Tz/9NNdddx39+/envLyc9evXX1GAFeJS2XJyOPfxx1R+8SUuiwUAbVAQIXfcQcg9P8QQHe3lEnadsroy1uWsY23OWnYW7sShONzHEgISmJU0i1mJsxgSPgStRgakXzZFXLLKykoFUCorKz3219XVKUeOHFHq6uq8VLIrV1BQoDz++ONKcnKyYjAYFH9/f2Xs2LHKG2+8oVgsFkVRFOXIkSPKhAkTFLPZrAwfPlxZtWqVAijr169XFEVR1q9frwDK119/rQwePFgxGo3K2LFjlf3797tf5/3331eCgoKUL774QklJSVFMJpMya9Ys5cyZM+5zXnzxRWXYsGHux/fff79y8803e5T3qaeeUqZNm6YoiqJs2rRJmTZtmhISEqKYzWYlLS1NWbp0aYv32PjazS1btkwBlMWLF3tcG1CysrI8zt28ebMydOhQxcfHR5kyZYry+eefK4By6tSpVj/T999/313G1tx///0K0GJ76KGH3GVLTU1VzGazcuONNyp/+ctflPN/dL/44gtl+PDhitFoVMLDw5Vbb73VfSwpKUn5wx/+4H68YcMGxc/PT/nTn/6kKIqi/P73v1diYmIUs9mszJkzR/nwww8VQCkvL1cURVEef/xxpW/fvorJZFIiIiKUe++9VykrK2vz/XQnveFn8mrlcrmUmm3blZyfPKocGThIOTJgoHJkwEAl+/oblHP//FRx1tZ6u4hdJq86T/ng0AfKfSvvU4YuGaoMWTLEvc3/ar6yaN8iJetsluJyubxd1G6vre/v82kU5bwODeKiqqqqCAoKorKyksDAQPf++vp6Tp06RXJyMj6ySKEQXic/kz2Py2qlavlyzn34EdZmfSD9pk4h9L778Zs08ZKnFemJTlacdA+5P3ruqMexIWFDmJk0k1mJs+gT1Mc7Bews1mo4vRkq82DshQeaXI62vr/PJ81qQgghvM5eXEz5P/9JxdLP3Et7aMxmgm+ZT8h/3YspJdnLJexciqJw5OwR95D75rNUazVaRkWNYmbiTGYmziTarxc1I7qcUJAJJ9bByfWQuwNcDtD7wIh7weCdP2okHAkhhPAKRVGoy8yk/KOPqVq1Chxq/xl9bAyh9/yXOhQ/KMjLpew8DpdDnZTxzBrW5a6jyFLkPqbX6pkQM4FZSbNIT0jvXbNTV+SoYejEejiZAfUVnsdD+kDfGWCzSDgSQghxdXDZbFR/+y3nPvyI+mbLCvmOHk3IffcSMKP3DsW3Oq1sL9jOmpw1ZORmUGGtcB8z681MjpvMzMSZTI2fSoCxl8zRVF+lNpU11g6dzfY8bgqC5ClqIOo7HUK9P9Fs7/zfJ4QQottxlJZS/ulSypcuxVlWBoDGaCTwxhsJvfe/8OmloyBrbDVszNvI2py1bM7f7DEHUZApiPT4dGYlzWJ8zHh89L2gb5zTAQX7mjWV7QTF2XRco4P40Q1haAbEjgRd94oj3as0Qgghep26/fs59/E/qPr2W/cs1vrISEJ+eDfBd9yBvp2TovYkZXVlrM9dz9qctewo3IHD1TTkPtI30t1/aFTUKPTaXvBVfO5kQzPZeji5EayVnsdDUyBluhqG+kwGc7BXitleveBfpPtpnLFZCOFd8rPoPS6bjepvvuHcx/+gvmGJHwDzsGGE3Hcvgddei8bQu2Zqzq3Odc9BlFmSiULTYPA+gX2YmTiTGYkzesccRLXn4NRGNQydWA8VZzyP+wRBSnpDIJqu9iPqQSQcdSCj0YhWq6WgoICIiAiMRuNVMeRUiO5GURRsNhulpaVotVqMRqO3i3TVsBeXULH0U8qXfobz7FkANAYDgddfT8h//RfmoUO8XMKOoygKx8qPuQPR8fLjHscHhw121xClBHu/H80Vcdggb6cahE6sg8JMUJr98aHVQ8K4ptqh2OHQg2fllnDUgbRaLcnJyRQWFlJQUODt4ghx1fP19SUxMdG9NIroHIqiULcvk/KPzxt1FhlJyN13qU1nYWFeLmXHaBxhti5nHetz15Nfk+8+ptPoGB01mhmJM5iROKNnD7lXFCg5qo4mO7keTm8Bu8XznPABaq1QynToMwlMvaQDORKOOpzRaCQxMRGHw4HT6bz4E4QQnUKn06HX66X2thO56uupWrGCc//4B9YjTRMVmkeNIvS/7iFg1qxe0XRW76hnW8E21uWuazHCzEfnw6S4ScxInMHUuKkE+wR7q5hXrrpIDUONQ+xrijyP+0U0NZWlpENQXNeXsYtIOOoEGo0Gg8HQYgV5IYToDWx5eZT/859U/usLnJVqx1uNyUTgDTcQ+l/34HPNNV4u4ZWrqK9gY/5G1uWsY2vBVuocde5jQaYgpsVPY2biTCbETsCsN3uxpFfAWgNntjZ0os6AkiOex/U+kDSxqd9Q5GC4SmphJRwJIYS4KMXlwrJlK+X/+Ac1GzaozS6AIS6OkB/eTdCtt6IPCfFyKa9MQU0B63PXsy5nHXuK9+BsNvw8xi+GGYkzmJk4kxGRI3rmCLPGIfaNYSh3J7jszU7QQMywpqayhHFem4TR2zr1X3fTpk1MmTKFLVu2MGnSpM58KSGEEJ3AWVVF5VdfUf6PT7CdaRqR5DdpEiH33IP/tKlodD2z462iKBwvP866nHWsy11H1rksj+P9Q/qr/YcSZjAwdGDPa6JVFDh7oikMndrUcoh9cGJTzVDyNPDtfdMqXI5ODUfffPMNer2eFStWSDgSQogepD4ri/J/fELl8uUodWqTktbfn6BbbyHk7rsxJffMtc7sLjt7i/eyPnc963PWU2BpGjyj1WgZGTmSGYkzmJ4wnfiAeC+W9DLVlDQNsT+5ASpzPY/7BEPy1IbaofRuMRt1d9Rp4eill17C4XAwY8YMnnzySX71q1/xwgsvdNbLCSGEuEKKzUbVd6so/+c/qdu7173f1C+VkHvuIWjePLR+fl4s4eWx2C1syd/C+tz1bMzbSJWtyn3MR+fD+NjxzEycybT4aYT49LCmQWsN5GxrGFWWAcWHPI/rjGrzWGMYihneo4fYdxWNoijKxU+7PO+99x6VlZUEBwfz8MMPd9bLdLmqqiqCgoKorKwkMDDQ28URQogrYi8ooPyzz6j4/F/uuYnQ6wmYPYvQH/4Q8+jRPa5JqbS2lIy8DNbnrGd74XbszfrWhJhCmJYwjekJ03teh2qnHfL3wqkNbfQbAqLTGkaVpUPiBDD6eqGg3VN7v787tVnN4XDw7LPP8u6773bmywghhLhEisuFZds2yv/5T2rWrYeG2cT1kZEE33kHwbffjiEy0sulbD9FUThRcYL1uevJyM3gQNkBj+MJAQnMSJjB9MTpDI8Yjq6n1J4oCpRmqU1kJzPUBVxt1Z7nBCVC33Q1DCVPA79wLxS0d+nUmqPeSmqOhBA9lbOigoplX1H+6T+xn8lx7/cdN46Qu+8mYOaMHjM3UeOEjI2BKLfas3/NkLAh7v5DfYP7dnjtV35FHeUWW5vHQ/yMxAVfRq1UZX5TzdDJDKgp9jxuDlH7DSVPa1iaIxl6WM2et3SLmqNGa9asYebMmT2uWlYIIXoDRVGoP3CA8n9+StU336BYrQBo/fwIuvlmQn54N6bUVC+Xsn0sdgtbC7aSkZvBhrwNVDYbfWXUGhkXM47pidOZFj+NSN/Oq/nKr6hjxm8zsDraXr/PpNey7tn0iwek2nNqjdDJDDUUnc32PK73UZvHGpvKotOumvmGvKVLwtGcOXMoLCwksgdV0QohRE/nqq2lcsUKKv75KfVHmib4Mw0aRMhddxF04w09ooN1saWYDXkbWJ+7nh2FOzz6DzVOyJiekM6k2En4Grqmf025xXbBYARgdbgot9hahiNbLeRub2oqK9wPzRapRaOF2BFNYSh+7FU735C3dEk4kpY7IYToOtbsbMo/XUrlv/+Nq1rtn6IxGgm87jpC7r4Ln2HDunVNfuOCro3NZUfOes7cnBCQQHpCOtMTpveMCRmdDijYq4ahUxsgdwc4z2uOCx8AKdPUprI+k8Ec7JWiClU3/x8lhBCiPVw2G9XfraJ86afU7d7j3m9ITCTkzjsJuvWWbj2Dtc1pY3fRbtbnrmdD3gYKLYXuYxo0DIsY5g5EyUHJ3TrcgcIATS6TtIdJ/O5vULSzZSfqwDg1CKVMU/sPBcZ6p6iiVV0WjhYvXszkyZMZOXIkId34B1QIIXoS25kzlH/2GZVfLsNZXq7u1OkImDGd4DvuxG/SRDTdtH9KRX0Fm/I3sT53PVvyt1DrqHUf89H5MCF2AtMTpjM1fiph5jAvlvTiEjTFTNIeZpL2EBO0RwjXNMyl1Njn3ScYkqc0BKJ0CEuVTtTdWJeFo0WLFvHSSy+h0WhISEhg5MiRHlt0dHRXFUUIIXo0xW6nev16Kj5dimXrVvd+fXQ0wbffRvBtt2GIivJiCdt2qvIUGbkZZORmkFmaiUtp6rcTYY5gavxU0hPSGR8zHh999+1no68t5ibtFiZpDzNRe5gEbanH8VrFxC7XAAZMuJHo4XMgeqhMvtiDdFk4Onz4MA6Hg3379rF371727t3Le++9R25uLhqNhujoaPLz8y/pmgsXLuTLL78kKysLs9nMxIkTee211xgwYMAFn/f555/zf//3f5w+fZp+/frx2muvcf3111/J2xNCiE5ny8un4l+fU/nFlzhKG76MNRr8pk4h5M678J86BY2+e/WWaBxu3zi67EzVGY/jA0IGkJ6QTnpCOteEXYNW0z1ruagrh9Nb1D5DpzYysDSLPxmbDtsVHfuUVLa6BrPFOYRMJRU7epYPm0x0bJD3yi0uS5f8FDW2DcfGxhIbG8sNN9zgPnb27Fn27NlDZmbmJV93w4YNPPbYY4wZMwaHw8Hzzz/Ptddey5EjR/BrYwTG1q1bufvuu1m4cCE33ngjn3zyCfPnz2fv3r0MGTLkst6fEEJ0FsVupzojg4rPPseyebM6KSCgCw8n+Ac/IPj22zHGx3m5lJ4qrZVsyd9CRl4Gm/M3U92sv41eq2ds9FjSE9KZFj+NWP9u2tfGZlGX5Ti1Ud0K90OzWi4FDYdcSWxxDWGbazC7XAOopfvWdIlL0yWTQGq1WoqKijp9KH9paSmRkZFs2LCBqVOntnrOnXfeicViYfny5e5948ePZ/jw4bzzzjvteh2ZBFII0dlarSUC/CZOJPiOOwiYMR2N0XiBK3StM1Vn3LVDe4v34lSc7mPBpmCmxk9lWvw0JsZOxN/o772CtsVhhbxdTWEob3fLZTnC+6t9hpKnUhA6mumL9nfMPEeiy3SrSSC//fZbgoI6v1qxslKdDCw0NLTNc7Zt28YzzzzjsW/OnDl89dVXbT7HarVibZg0DdQPVwghOlqbtURhYQTfeivBt9+GMTHRy6VU2V12Mksy2ZC7gQ15GzhdddrjeN+gvkxLUOcfSgtP637LdTgdULBPbSY7vQlytoOj3vOcoERImdo0vL7ZiLJYYN2zIZ0zQ7bwui4JR9dee22nv4bL5eLpp59m0qRJF2weKyoqIuq8jopRUVEUFRW1+ZyFCxfy0ksvdVhZhRCiOVtODhWf/4uKr5bhLC1z7+9utUQV9RVsLtjMxtyNbC44r7lMo2dU9CjS49OZljCNhIAEL5a0FS4nFB1Ug9CpjXBmW8vh9f5RDctyNGwhfS54ybhgs4SfXqp79dy7Ao899hiHDh1i8+bNHX7t5557zqO2qaqqioSEbvaDL4ToUVw2G9WrV1Px+b+o3b7dvV8XHk7wLfPVvkReriVSFIWTlSfZkLeBDbkbWowuCzGFMCV+ClPjpzIxdiIBxgAvlvY8igIlR5vC0OnNUF/heY45RK0RamgqI7y/DK8XQC8JR48//jjLly9n48aNxMfHX/Dc6Ohoios9F/ErLi6+4FQCJpMJk8nUIWUVQlzdrNnZVHz+OZVf/RtnQ1cANBr8Jk8m+PbbCJg+3asLv1qdVnYX7WZD3gY25m0kv8ZzFHFqcKq7M/XQ8KHdp7lMUaDsezi9EU5tUsNQbZnnOaZASJqoBqE+UyBqiKxRJlrVo8ORoig88cQTLFu2jIyMDJKTky/6nAkTJrB27Vqefvpp977Vq1czYcKETiypEOJq5qqtperb76j417+o27vXvV8fHa32JfrBrRjivDfirLS2lE35m9iQu4Fthduoc9S5jxm0BsZGj2VawjSmxk8lzr+bjIxTFDh3sqFmqCEM1ZzXPUJvhqQJahBKngYxw0DXo7/2RBfp0f9LHnvsMT755BP+/e9/ExAQ4O43FBQUhNmstgPfd999xMXFsXDhQgCeeuoppk2bxu9+9ztuuOEGPv30U3bv3s1f/vIXr70PIUTvoygK9QcPUvH5v6hauRKXxaIe0OnwT08n+Pbb8J8yBY2u62teXIqLw2WH2Zi/kY15G1usXdY4GePU+KmMjxnfZYu5XlT5aTUENYahqjzP4zoTJIxtqhmKGwV67/fVEj1Pjw5HixcvBiA9Pd1j//vvv88DDzwAQE5ODtpm1aYTJ07kk08+4X//9395/vnn6devH1999ZXMcSSE6BCO8nKqvv6ais//hfX77937DYmJBN96K0G33IIhqnOnNWlNta2abQXb2Ji3kU35mzhXf87j+JCwIUxNUIfbDwwd2D0mY6zIaQpCpzdDZY7nca0B4seoy3L0maLel9XrRQfoknmOehuZ50gI0ZzicmHZto2Kf/2LmjVrUezq/Dgak4mAOdcS/IPb8B0zukvXOFMUhVNVp9iUt4mNeRvZW7wXh+JwH/cz+DExdiJT4qYwJX4K4ebwLitbmypym4LQ6Y1qOGpOq1drg/pMVsNQwjgwdpNaLdEjdKt5joQQojey5eVR+eUyKr5ahqOgaRV5n2uuIei2HxB0443ouvAPKKvTyq6iXe5AlFfj2ezUJ7CPe3TZqMhRGHTe6/gNnBeGNkGF59IiaHQQN9IzDJm64QSSoteRcCSEEJfAVVenDsH/cpnHEHxtQABB8+YRfNsP8Lnmmi4rT2FNIZvyN7EpbxM7ina06Ew9Omq0u/9QYqCXJ5BsTxiKHd4QhqZC4jgwdaPpAcRVQ8KREEJchLtz9RdfUrViBa6aGvWARoPfhAkE3XorAbNmovXp/P4udped/SX72ZSv1g5lV2R7HI/0jWRK3JTu0Zm6/IwahM5saQhD5zWTeYShhpohH+mqILxPwpEQQrTBUVpK5dfLqfjyC2zZJ9z7DXFxBN16C8Hz53fJEPyyujI2529mU94mthVso9reNLOzVqNlWMQwpsZPZUrcFPqH9Hcv9t2lFEUdTXZmS0PN0JaWHag1OnU4fWMHaglDopuScCSEEM24bDZq1mdQuWwZNZs2gVNdQNXdufrWW/EdO7ZTO1c7XU4Olh10N5cdPXfU43iwKZhJcZOYGjeVSXGTCDJ1/tqVLSgKnD0BZxqC0JmtLYfWa/UQO0KtGUqaLM1koseQcCSEuOopikL9kSNUfrmMquXLm2auBszDhhF0yy0E3nA9uoDO+2I/W3eWrQVb2ZS/ia0FW6m0VnocHxw2mCnxU5gSN4XBYYO7fmZqRYHSrKZmsjNbocZztQG0hobRZJPUQBQ/VjpQix5JwpEQ4qrlKCuj8uvlVC5bhvX4cfd+fWQkQTffTNAt8zGlpHTKazfWDm3O38zm/M0cOXsEhaaZVQKMAUyKncSU+ClMjJ3YrqH2+RV1HbdKvMsJxYfVEHRms3pbe9bzHJ1JnVuozyR1WY74sTK0XvQKEo6EEFcVl9VKzbp1VH71b2o2b25qNjMaCZg1i6BbbsFv4oROmbn6YrVDg0IHMTluMlPipzA0fCh6bft/RedX1DHjtxlYHa42zzHptax7Nr31gOS0Q+H+hpqhrZCzHc4rH3qzOgN1n8mQNEmtJZJJF0UvJOFICNHrKYpC3b59VH71b6q++QZXdVOHZp+0NIJvvYXA665DF9SxfXccLgeHyg6xKX8TW/K3tKwdMgQwIXYCk+MmMzluMhG+EZf9WuUW2wWDEYDV4aLcYlPDkb0O8vfAmW1qzVDuTrDXej7BGACJ49X1yZImq/2HZDkOcRWQcCSE6LVseXlU/vvfVP77P9hzmkZO6WNiCLrpJoJuvqnDm81KakvYkr+Fzfmb2Va4jWpbtcfxxtqhyXGTSYtIu6TaoSvhTy2jtN8TtWsLfLdXDUbO85rgzCFqjVDSRHWLGioLtYqrkvyvF0L0Ks6qKqq+/ZbK//yHut173Ps1vr4EXnstQfPn4zt2TIeNNrM77WSWZroD0bHyYx7HA42BTIydyOS4yUyMnXhFtUOXIpQqxmiPMVabxVjtUa7RnEGnUSCz2Un+0Q21QpPULWIgdOESJ0J0VxKOhBA9nmKzUbNpE5X//g8169e71zZzT9I4/2YCZs1C69sxnYXzqvPYWrCVzfmb2VG4g1pHU3OUBg2DwwYzOX4yk2InMTR8aOePLFMUDNW53KLdxBhtFuO0WfTVFrY47YwrksCB0wgZlA6JEyA0BbwxJ5IQ3ZyEIyFEj6T2I8qk8uv/UL3yG4/h96Z+/Qi6+SYCb7wRQ3T0Fb9WnaOOXUW72JK/ha0FWzldddrjeKhPqLvv0MTYiYT6hF7xa16QywWlR5s6TudsY0BVPn84rztQliuBna6B7HINYKdrIMWEsnzaZELivDAvkhA9iIQjIUSPYjt9Wh1+//XXnv2IIiIInDePoJvmYRow4IpmiVYUhe8rvmdr/la2FGxhb/FebK6m/jk6jY5hEcOYHDeZSXGTGBg6EK2mE5ujHFYoyIScxjC0HeorPMus0bPP2ccdhna7BlCJzDEkxOWQcCSE6PbsJSVUf/MNlV8vp/7QIfd+ja8vgbNnE3TzTfiOG3dFw+8r6ivYVriNLflb2FawjZK6Eo/jMX4xTIqbxOTYyYyNGUuAsRNneq6rgLxdkLNNHU2WvwecVs9zDL7qHENJEyFxAke0/bn1nb2dVyYhriISjoQQ3ZKzuprqVaupWrEcy/YdalMSgE6H38SJBN10EwEzZ1x2PyK7y87B0oNsKdjC1vytHD572GOYvY/Oh1HRo5gUO4lJsZNIDkruvDXLKvPczWPk7IDiQ9CsLAD4hqudpxMnqMPro9NAZ3AfDq6ow6TXXnSeoxA/GYovxMVIOBJCdBsum42aDRuoWr5C7Vhta2rKMg8fTuC8GwmcOxd9WNhlXT+3KlcNQwVb2Vm0E4vd4nG8X0g/JsVOYmLsREZGjcSkM13R+2mVywklR5qax3K2t1yTDNTO0okTmrawvhfsPB0XbGbds+kdN0O2EFcxCUdCCK9SHA4s23dQtWIF1WvWeEzQaOzbl6B5NxJ4ww0YExIu+drVtmp2Fu5ka8FWthZsJa/GM4QEm4KZEDOBiXETmRg7kUjfyCt+Py3YLJC3W51kMXe7emut8jxHo4OYNDUEJYxTbwOiLvml4oLNEn6E6AASjoQQXU5xuajbt4+qFSuo+vY7nOfOuY/po6IIvOEGgubdiGngwEtqymqckXprwVa2FWzjYNlBnIqz6doaPcMjhzMxdiIT4yYyKHRQx3ekripQa4Nyd6i3RQehWRkAdebphDFNYSh+NBj9OrYcQojLJuFICNElFEWh/vARqlaupOqbb3AUNs3DowsJIWDuHIKuvx7zqFHtnqBRURRyqnPYVrCNrQVb2VW0ixp7jcc5fQL7qGEodiKjo0fjZ+jAENK4OGvujoYwtAMqc1qeFxgPieMgYbzaXyhqMHT23EdCiMsm4UgI0WkURcF6/DhV33xD9TffYjtzxn1M6+9PwKxZBN5wA37jx6ExGC5wpSYV9RXsKNrBtoJtbC/cTn5NvsfxIFMQ46LHMSF2AhNiJxDnH9dxb6i+Uh1FlrtTDUN5u8HmGcbQaCFqiBqCEsapt0HxHVcGIUSnk3AkhOhw1uxsqlZ+Q9W332I7edK9X2My4T99OoE3XI//1KloTRfv8Gx1WtlXso/tBdvZVriNo2ePeowq02v1jIgcofYdip3IwNCBHTMjtaLAuZNNQSh3p9qR+vxRZKZAtVksYVxTE5mpE4f5CyE6nYQjIUSHsJ48RdW331D9zTdYv89279cYjfhNnULg3OvwT09H53/hZi2X4iLrXBbbC7ezvWA7e0v2Yj1vjp++QX3dNUOjo0bja+iAZUFsFijY1xCEdkHeTqg92/K8kOSGIDRWrRWKGChNZEL0MhKOhBCXzXryJNXffUfVt99hPdZswVWDAf/Jkwm8bi7+M2ag87/wTM251blsL9zOjsId7CzcSbm13ON4hDmCCbETGB8znnEx4658VJmiQMWZphCUuwOKDrXsOK0zQexwdbLFxmYy/04Y0SaE6FYkHAkh2k1RFKzff0/1d6uoXvWdRw0Rej1+EycQeN31BMycgS4wsM3rnKs/x87CnWrtUCv9hnz1voyJHsP4mPFMiJ1ASlDKlU3AaKtVa4XydjYEol1gKWl5XkCsWiOUMFYNQtFDQd8Jcx0JIbo1CUdCiAtSFAXrsWNUffcd1d+t8uhDhMGgBqJr56iBKDi41WvU2mvZU7yHHYU72F64nWPlxzyO6zV60iLS3DVDQyOGYtC2r4N2KwWG8lNNcwvl7Wy9VkhrUOcWih/TFIak47QQAglHQohWKC4X9QcPUr16NVWrV2M/0zQ8XWMw4DdlCoFzrsV/+vRWa4hsThv7S/ezs2gnOwp3cLD0IA7F4XFO/5D+jIsZx/iY8YyKGnX5Q+zrq9S1x/J3q4Eob1frfYUCYpqCUPwYiBkGBpkwUQjRkoQjIQSgzlRdu3s31avXUL1mDY7iYvcxjcmE/9QpBFw7B//p6S36EDldTo6eO8qOwh3sKNzBvpJ91DvrPc6J8IlhSOgohoSOZnDoKIKMIZe+nIXLCaVZTSEob7f6+PwRZDqjuvZYYxCKH6PWCnXW2mhCiF5FwpEQVzGX1Ypl61aqV6+hZt06nBUV7mNaX1/809MJuHY2/lOmoPVrqtlxKS6+L/+eXUW72FG0gz1Fe6i2V3tcO9gUytmyBOw1fXFYUqm2h3IS+A8AhwF1IdR1z6a3HZCqCptqhPL3qP2Gzp9XCCA4qSEEjVZvpa+QEOIKSDgS4irjrKykZuNGqtesxbJpE67aWvcxXXAw/jNnEDB7Nn4TJrjnIVIUhVOVp9hZuJOdRTvZVbSrxYgyf4M/o6NHMz5mPGOjx1JviWDeoi0XLIvV4aLcYlPDkc0ChfsbglBDIKrKb/kkoz/EjmiqEYofLSPIhBAdSsKREFcBe34+1WvXUb1uHbW7d4Ojqf+PPiqKgNmzCZg9G99RI9Ho9SiKQl51HjvP7GRX8S52Fu6ktK7U45pmvZmRkSMZEz2GcTHjGBg6EL226VfKodrKNsujw0k/TT7DtCeI3fgfqDikTrCouDxP1Ggh8hqIG6WGoLjREDFA5hUSQnQqCUdC9EKKomDNylID0dq1WI8e9Thu6peK/4yZBMycgc+QIaDRkF+Tz6pTX7OraBc7i3ZSXFvs8Ryj1sjwyOHuMDQkbAgGXXtGlCnEa8oYpjlBmvYEw7UnGKo5ha+mYWLH5gPXAuMaaoUamsdihoPpwnMkCSFER5NwJEQv4bJaqd2xg+r166nJ2OCxsCtaLeaRIwhoCESGxETya/L5rng3u7YuZVfRLgothR7X02v1pIWnMTp6NOOixzEschgmXTv68dSUQsFeIrO28r5hPWnak4RpqlucVq2YOeBKof+odCIGTFRrhwJjrvRjEEKIKybhSIgezFFaSs2GDVRnZGDZug2lWf8hjY8PfpMmETBjBn7p0yg0WNhQvJvdOe+ya9cuiixFHtfSa/QMjRjK6KjRjI0Zy7CIYZj1FxlJVl+p9hMq2NcwnH6fe1X6SCCyofXLpujIUhI54EohU0kl09WXE0osClqWj5lMRFxQR34sQghxRSQcCdGDNDaX1WRkUL0+g/oDBzyO6yMj8U9Pxz89ndLB0WyvPMSuou3sWfsWJXWeM0LrNXoGhw9mTPQYxkSPYXjE8AuvUWazQNHBhiC0V709+30rJ2ogvB/lIUP5wxF/9rv6kqUkYsXYAZ+AEEJ0PglHQnRzzpoaLFu3UrNxI5YNG3GUenaM9hkyBL/0aZwblcLOoHPsKdnLnuIFnPv2nMd5jc1ko6JGMSZ6DMMihrUdhuz1UHxIDUAFmept6dGWHaYBghIhboTaVyhulNpPyCeQ/PxKPjy4uWM+BCGE6EISjoToZhRFwXbiBDUbNlKzcSO1e/Z4jC7TmM2Yx4+lavQA9vfTsd1+nL0ln1B93LNfj1FrJC1C7TM0Omo0aRFprTeTOaznBaFMNQi5HC3P9Y+GuJFqEIodqS7K6hfe6vsI8TNi0muxOloJVA1Mei0hflKjJIToXiQcCdENOGtqqN2+nZrNm7Fs3IS9oMDjuL5PEpbRA8jq70tGWAl7K/ZS59gCp5rOMevNjIgcwaioUYyKGsWQ8CEtO1Db66H4MBTua+grlKkOoW8tCPmGq+EnZnhDIBp5SR2m44LNrHs2nXKLrc1zLnmGbCGE6AI9Ohxt3LiRN954gz179lBYWMiyZcuYP39+m+dnZGQwffr0FvsLCwuJjo7uxJIK4Ulxuag/fBjL5s3UbN5CXWYmOJstjGo0Up+WyolBgayPr2Kr5iQOJV9dJaNMPSXQGMjIqJGMjhrNqKhRLeYZUlei39kUggozoeRoywVYAcyhahCKbWgeixneIcttxAWbJfwIIXqcHh2OLBYLw4YN47//+7+59dZb2/28Y8eOEdhssczISJldV3Q+e3EJli1bsGzejGXrVo+lOgBssWGcGRTC5oRa1oYXYzMcbzqoQJRvFCOjRjIyciQjo0aSGpyKVqNVj9dXQs52NQgVHlBvy4613kfIN0wNPzHD1C12BAQnyrpjQgjRoEeHo+uuu47rrrvukp8XGRlJcHBwxxdIiGac1dXU7tqFZdt2LNu2Yss+4XHc4Wskp18w2xLq2JpQS2lwJdA4q7SG1OBURkSOcAeiWP9Y9VBNKRTth0NfNwWh8lO0yi+yqWksZph6PzBOgpAQQlxAjw5Hl2v48OFYrVaGDBnCggULmDRp0gXPt1qtWK1W9+OqqqrOLqLogVw2G3WZmdRu345l6zbqDh70aCpTNFCY4MfOBCt7khWyY504deqIMoPWyIjwIQyPHM6IiBGMiBxBsDFQDT1FB2HHX9XbwgNQU9R6AYISISZNDUHRDbcyqaIQQlyyqyocxcTE8M477zB69GisVit//etfSU9PZ8eOHYwcObLN5y1cuJCXXnqpC0sqegLF4aD+6FFqd+zAsmMntbt3o9TVeZxTGqZnb6KTQ300HErSYDGrITvEFMrUyOGMiFSD0DWBKRjPnYCiQ3DoGyh6Tb1vazmzNGggrG9Ts1hjEPIN7YJ3LYQQvZ9GURTF24XoCBqN5qIdslszbdo0EhMT+eijj9o8p7Wao4SEBCorKz36LoneTXE6qT9ylNqdO9Vtzx5cNTUe51T6aTiQBIf6aDjYR0NZkNp8lRqcyrCIYQyLGMYI/ySSas6hKTmsBqCig1B2vPWO0joTRF0D0UPVEBSdBlGDZb0xIYS4DFVVVQQFBV30+/uqqjlqzdixY9m8+cIT1ZlMJkymdqwpJXoVtWYoi9rdu9Xaod27Uc4LQzU+cDRBw+FENQzlRoC/MYC0iDRuDxvCMGMYQ20OAs9mw4kDsPUTqClu/QXNoRA9pCkERQ+F8H7QrsVdhRBCdJSrPhxlZmYSEyP9Mq4W+RV1bc+7U19PwOlj+GYdpHbPXmoz90GtZzOZxdQQhpLUQJQTqSE5pC9Dg1P5b30ww2wO+pbno80+BNu+AJe9lRdqaBaLHgpRQxpqhYZCQIx0lBZCiG6gR4ejmpoasrOz3Y9PnTpFZmYmoaGhJCYm8txzz5Gfn8+HH34IwJtvvklycjKDBw+mvr6ev/71r6xbt45Vq1Z56y2ILpRfUceM32a4Z2wOsFm45uxphpw9yZBz2aSWF1CrKNQ2e06NDxyLawpDVYnBDAnsQ5rOn5utNoZUFBFw/ADUrm/9RU2BajNY1BC1eSxqqHpr9Ov8NyyEEOKy9OhwtHv3bo9JHZ955hkA7r//fpYsWUJhYSE5OTnu4zabjZ/+9Kfk5+fj6+tLWloaa9asaXViSNG7KA4H5fsPMfP7zQyqyGbQudPEVbccdXg2ALLiNRxN0HAiyUBAQgxDdGam1NXx2Lk84o4fRMPBli+g0UJoX7VZzB2GBkNQgtQGCSFED9NrOmR3pfZ26BLe4ygvp25fJmW7t1C5Zxe6rFPorS2buPLCICtBw9F4LcWRPgSbNMwy2BlZUUBfm7X1vx78o9Xan8iGLeoaiBgIBpkJWgghujPpkC2uGs6aGmoPHqRoz2Yq9+9Bm3USc2nTEPjGrvS1RsiO1XA8Fkqj9BjCnKRpqphiq+cRmx2fWgWPNjVTEEQOVINP1OCGIDRYhswLIUQvJ+FI9Ciu2loqjxwgf9cGqg7sQ5d1ioCiKjQN9Z/Ne/Lkh8L3cRoqYjQYQuuJ8a9nsN3GPJsdX0WBhoFnNYoPWUoKx13xHFPiuXfeXPoMGgWBsdIkJoQQVyEJR6LbspYUk7tvE6X7d2I9moXpZAGBJRa0CuiAkGbnlgbCqRiojnChC7MTHmSlv8bGJJsdc0PLsdMWSGZ9El+74vheiSNbiSVbiSdfCUNB677WLQmTISioa9+sEEKIbkPCkfA6h7We/MM7KTq4k+qsw/D9KQLPnCWg2gFA8Hnnl/tBXjRYwp3ow+xEBlvpp7Myye5AB2qfoPBr1DmCIgZCxACIGMjRSh9+sGhLF787IYQQPY2EI9FlbPW15BzdQeHBndQcOwKncvHLPUtomQ2dAoGoWyMXUBQK5yJcOEKd+AbbiQ6sp7/OxkR06lxB4f0gvD+ENdyGp4JPG7U+VZWt7xdCCCGakXAkOpTD6SA/9whFR/dQnn0E28lstDlF+JdUE3bWid4F4ahbcxYTlIYp1IW40IXaCQmy0ce/nkHBsWhD+0JYqhqGwlIhNAWCk0B3af99Q/yMmPRa9zxHrTHptYT4GS/9jQshhOg1JByJS1ZVW07BqYOUnjhM9ckj2HLOoCkowVxiIeysA19ry1qgRnVGOBuqUB/qRBfkICjQSWxMEKnRyRjCUiA0uSEA9YWQPmDw6bByxwWbWfdsetszZKMGqLhgGZIvhBBXMwlHwoPNYaOk7BRnTx2h/PQRavNOYi0ogOJzGM/VElDuILQKdErrNUCgNodVBCnUBCm4gpwYg7UERwYQHR/PgOTBaMNSICRZDUJBCV26dlhcsFnCjxBCiAuScHQVqHfUc674JJW5WVQXnaK2JJ/6cyXYystxnqtAU16LscqOb7WToBrwsYMRiLrANe06qApQqAtQcAVpMIT4EBgTRlSfZGJSR6CPSoHgRDX8+MhEmUIIIXoOCUfdhMvp5PS2L3Da63HarDjtdTjsNlx2G06bFZfDitNux+mw4bJasVuqsddacNbW4aqzotTbUKx2tFYHWqsLndWFqc6FuQ7869SaHoCAhu1iak1g8Vew+WlwBRnQh/niFxVOaEIfolKGEJA8FE1wPAREg9508QsKIYQQPYSEo25CcTmxPvyi+7EWtfamI9UboM4MVjM4zDpcvno0QT4YQwPxjYgkMCaR8D4DCUkehi48AcwhMgmiEEKIq46Eo25CZzBS4wOKBlxaz03Rgkurabqv0+AyanGZtCgmPRqTAa2PAa3ZB53ZF4OfL0b/QHzDogiMSiI4NpWAuAHogiJBq714YYQQQoirmISjbmRM5lFvF0EIIYS46kk1ghBCCCFEMxKOhBBCCCGakXAkhBBCCNGMhCMhhBBCiGYkHAkhhBBCNCOj1S6DoqgzKlZVVXm5JEIIIYRor8bv7cbv8bZIOLoM1dXVACQkJHi5JEIIIYS4VNXV1QQFBbV5XKNcLD6JFlwuFwUFBQQEBKDpwBmkq6qqSEhIIDc3l8BAWY+ss8nn3bXk8+5a8nl3Lfm8u9blft6KolBdXU1sbCzaC0yKLDVHl0Gr1RIfH99p1w8MDJQfri4kn3fXks+7a8nn3bXk8+5al/N5X6jGqJF0yBZCCCGEaEbCkRBCCCFEMxKOuhGTycSLL76IyWTydlGuCvJ5dy35vLuWfN5dSz7vrtXZn7d0yBZCCCGEaEZqjoQQQgghmpFwJIQQQgjRjIQjIYQQQohmJBwJIYQQQjQj4aibeOutt+jTpw8+Pj6MGzeOnTt3ertIvdbGjRuZN28esbGxaDQavvrqK28XqVdbuHAhY8aMISAggMjISObPn8+xY8e8Xaxea/HixaSlpbknx5swYQLffPONt4t1VXj11VfRaDQ8/fTT3i5Kr7VgwQI0Go3HNnDgwA5/HQlH3cDSpUt55plnePHFF9m7dy/Dhg1jzpw5lJSUeLtovZLFYmHYsGG89dZb3i7KVWHDhg089thjbN++ndWrV2O327n22muxWCzeLlqvFB8fz6uvvsqePXvYvXs3M2bM4Oabb+bw4cPeLlqvtmvXLt59913S0tK8XZReb/DgwRQWFrq3zZs3d/hryFD+bmDcuHGMGTOGRYsWAerabQkJCTzxxBP84he/8HLpejeNRsOyZcuYP3++t4ty1SgtLSUyMpINGzYwdepUbxfnqhAaGsobb7zBQw895O2i9Eo1NTWMHDmSt99+m9/85jcMHz6cN99809vF6pUWLFjAV199RWZmZqe+jtQceZnNZmPPnj3MmjXLvU+r1TJr1iy2bdvmxZIJ0TkqKysB9QtbdC6n08mnn36KxWJhwoQJ3i5Or/XYY49xww03ePweF53n+++/JzY2lpSUFO655x5ycnI6/DVk4VkvKysrw+l0EhUV5bE/KiqKrKwsL5VKiM7hcrl4+umnmTRpEkOGDPF2cXqtgwcPMmHCBOrr6/H392fZsmVcc8013i5Wr/Tpp5+yd+9edu3a5e2iXBXGjRvHkiVLGDBgAIWFhbz00ktMmTKFQ4cOERAQ0GGvI+FICNFlHnvsMQ4dOtQpfQREkwEDBpCZmUllZSX/+te/uP/++9mwYYMEpA6Wm5vLU089xerVq/Hx8fF2ca4K1113nft+Wloa48aNIykpic8++6xDm40lHHlZeHg4Op2O4uJij/3FxcVER0d7qVRCdLzHH3+c5cuXs3HjRuLj471dnF7NaDSSmpoKwKhRo9i1axd//OMfeffdd71cst5lz549lJSUMHLkSPc+p9PJxo0bWbRoEVarFZ1O58US9n7BwcH079+f7OzsDr2u9DnyMqPRyKhRo1i7dq17n8vlYu3atdJHQPQKiqLw+OOPs2zZMtatW0dycrK3i3TVcblcWK1Wbxej15k5cyYHDx4kMzPTvY0ePZp77rmHzMxMCUZdoKamhhMnThATE9Oh15Wao27gmWee4f7772f06NGMHTuWN998E4vFwoMPPujtovVKNTU1Hn9lnDp1iszMTEJDQ0lMTPRiyXqnxx57jE8++YR///vfBAQEUFRUBEBQUBBms9nLpet9nnvuOa677joSExOprq7mk08+ISMjg++++87bRet1AgICWvSd8/PzIywsTPrUdZJnn32WefPmkZSUREFBAS+++CI6nY677767Q19HwlE3cOedd1JaWsoLL7xAUVERw4cP59tvv23RSVt0jN27dzN9+nT342eeeQaA+++/nyVLlnipVL3X4sWLAUhPT/fY//777/PAAw90fYF6uZKSEu677z4KCwsJCgoiLS2N7777jtmzZ3u7aEJcsby8PO6++27Onj1LREQEkydPZvv27URERHTo68g8R0IIIYQQzUifIyGEEEKIZiQcCSGEEEI0I+FICCGEEKIZCUdCCCGEEM1IOBJCCCGEaEbCkRBCCCFEMxKOhBBCCCGakXAkhBBCCNGMhCMhhBBCiGYkHAkhhBBCNCNrq10Gl8tFQUEBAQEBaDQabxdHCCGEEO2gKArV1dXExsai1bZdPyTh6DIUFBSQkJDg7WIIIYQQ4jLk5uYSHx/f5nEJR5chICAAUD/cwMBAL5dGCCGEEO1RVVVFQkKC+3u8LRKOLkNjU1pgYKCEIyGEEKKHuViXGOmQLYQQQgjRjIQjIYQQQohmJBwJIYQQQjQjfY6EEEII4XUuiwVHaSn2khJc1dUEzJzptbJIOBJCCCFEp3FZLNhLSnCUlOIoKW64LcFRWqpuDfddtbVNTzIYGHhgv9fmEpRwJIQQQohLptjtak1PcTGO4pKG4FOiBqHiEjX0lJTgsljafU2try/6yEj0EREotbVo/Pw68R20TcKREEIIITy4amvV0FNUhL2oGEdxUVMIKirCXlKC8+xZUJR2XU/r7+8OPfrISPSREer9iAgMjfsjItB6KQydT8KREEIIcRVx1ddjLyxUQ05hEfaiQhxFxdiLi3AUqiHIVVnZvosZDBgiItBHRaGPisIQFYk+MqohAKkhyBAZ2W1CT3tJOBJCCCF6CcXhUJu6CgqwFxQ2hKDChhBUhKOwEGdFRbuupfH1xRAdjSE6Cn1UNPqoSAzR0Wr4abivCwlBc4E1ynoqCUdCCCFED+GyWBqCTwH2wkI1ADXeLyzAUVwCTudFr6P19UUfE6OGnegoDNEx7ltDdBT66Gi0/v5X7eLqEo6EEEKIbkBRFFxVVdjz87Hl5+NoDEEFBerj/AKc7Wnu0uvVGp+YGAyxMeijYzDERKNv3BcdjTYw8KoNPu0h4UgIIYToIs6G8NO42fLyPR67amoueg1tYCCG2NimrSEEGWJi0MfEog8PQ6PTdcG76b0kHAkhhBAdxGW1qkEnLw9bbi72vIb7+XnY8/JxVVVd9Bq6sDAMcXFN4SeuMQjFYYiLRefv3wXv5Oom4UgIIYRoJ0VRcJaVYcvNw56b03Cbiy1PvXWUlFz0Gu7wExeLMS6u4X4chvh4DDExaM3mLngn4kIkHAkhhBDNKHa72s8nJwdbTg72nFy1FqghBCl1dRd8vtbXF0NCAob4eIzxcRji4jEkxGOMj8cQF4fW17eL3om4XBKOhBBCXHVc9fVq2MnJwXYmB1tuDvYzOWoIKii48IgvjQZ9TDTGhMSG0JOAMVENQ4aEBHTBwdLZuYfr9uHorbfe4o033qCoqIhhw4bx5z//mbFjx7Z67pdffskrr7xCdnY2drudfv368dOf/pR7773Xfc4DDzzABx984PG8OXPm8O2333bq+xBCCNG1XHV12HJzsZ05gz0nB9vpM+7aIEdR0QVnd9b4+GBMiMeQmIQxIUENQYmJ6v3YWDRGYxe+E9HVunU4Wrp0Kc888wzvvPMO48aN480332TOnDkcO3aMyMjIFueHhobyy1/+koEDB2I0Glm+fDkPPvggkZGRzJkzx33e3Llzef/9992PTSZTl7wfIYQQHUux2bDl5WM7fRrbmTMet46iogs+V+vvjzEpCWNSIoaERDX8JCZgSExCHxHeKyc3FO2jUZR2LoziBePGjWPMmDEsWrQIAJfLRUJCAk888QS/+MUv2nWNkSNHcsMNN/DrX/8aUGuOKioq+Oqrry67XFVVVQQFBVFZWUlgYOBlX0cIIcTFKS4XjuJibKdPYz11Ctup0+4QZM/LA5erzedqAwPVAJSYiDEpEWNSEobERIx9+kjz11Wovd/f3bbmyGazsWfPHp577jn3Pq1Wy6xZs9i2bdtFn68oCuvWrePYsWO89tprHscyMjKIjIwkJCSEGTNm8Jvf/IawsLA2r2W1WrFare7HVe0YiimEEOLSOGss2E6dwnb6FLZTp9QgdFqtBbpQJ2iNr68agPokNdz2cd9KABKXo9uGo7KyMpxOJ1FRUR77o6KiyMrKavN5lZWVxMXFYbVa0el0vP3228yePdt9fO7cudx6660kJydz4sQJnn/+ea677jq2bduGro1JsxYuXMhLL73UMW9MCCGuYorLhaOoCOtJNQDZTp1U7588eeFh8Ho9xvh4jMnJavjp0wdjch+MSX3QR0ZIABIdqtuGo8sVEBBAZmYmNTU1rF27lmeeeYaUlBTS09MBuOuuu9znDh06lLS0NPr27UtGRgYzZ85s9ZrPPfcczzzzjPtxVVUVCQkJnfo+hBCiJ3PZbNjPnMF64iTWkyewnTyl3p66cC2QLjwcU58+aghqDELJfTDGx6MxGLrwHYirWbcNR+Hh4eh0OoqLiz32FxcXEx0d3ebztFotqampAAwfPpyjR4+ycOFCdzg6X0pKCuHh4WRnZ7cZjkwmk3TaFkKIVqhNYSexZp/AdvIE1hMnsZ04gS0vr+3h8AYDxsRETCnJGJNTMCYnN9xPRif9OEU30G3DkdFoZNSoUaxdu5b58+cDaofstWvX8vjjj7f7Oi6Xy6O/0Pny8vI4e/YsMTExV1pkIYTotZxVVVhPnMB24gTW7BNYT5zAeiIbR0Fhm8/R+vtj7JuCKTlFve3bV60NSkhAo++2Xz9CdN9wBPDMM89w//33M3r0aMaOHcubb76JxWLhwQcfBOC+++4jLi6OhQsXAmrfoNGjR9O3b1+sVisrV67ko48+YvHixQDU1NTw0ksv8YMf/IDo6GhOnDjBz372M1JTUz2G+gshxNXKWV2NNTsba3Y2tuxsrN+r9y/UH0gXHo4pJQVTal+MKX0x9U3BmNJX+gKJHqtbh6M777yT0tJSXnjhBYqKihg+fDjffvutu5N2Tk4O2mbzUFgsFh599FHy8vIwm80MHDiQjz/+mDvvvBMAnU7HgQMH+OCDD6ioqCA2NpZrr72WX//619JsJoS4qrhqa9Xan+PfY/3+e3cgutDcQProaEx9+zaFoNS+GFNS0IeEdGHJheh83Xqeo+5K5jkSQvQUis2G9fRpNQA1BqHvv8eem9vmc/RRUZhSU9WtXyrGvn0xpaaiCwjowpIL0fF6/DxHQggh2k9RFBwFBdQfP4712HGsx49j/f441lOnweFo9Tm6iPCGENSvIQj1w5TaVzpFi6uehCMhhOhhnNXVavg5fpz6Y8fUGqHjx3HV1LR6vtbfXw0+/furt/36YeqXij40tItLLkTPIOFICCG6KcXpxJaTg/XYceqPZak1QllZ6qrxrTEYMCUnYxowAFP/fvj074+pf3/00dHSMVqISyDhSAghugFnTQ3WY8eoP5ql3h47hvX779ucMFEfHY1pQH98+g9Qa4T698eU3EdWixeiA0g4EkKILqQoCo6iIuqPZlGfdRRr1jHqs7Kw5+S0er7GZFKbwQYOwGfAQDUQDRiALiioi0suxNVDwpEQQnQSxeHAduoU9UePUn/kKPVZWViPHsVZWdnq+froaHwGDMA0cCA+AwdgGjAAY1ISmjbWfRRCdA4JR0II0QFcdXVqB+nGIHT0KNbjx1Fam6Ffr8eUkoLPoIGYBg5Sg9DAgTJfkBDdhIQjIYS4RM6aGqxHj1J/5Ih7s544CS5Xi3O1vr5qTdCgQfhcMwifQYMwpqailb5BQnRbEo6EEOICnBUVam3Q4cNqEDp8BNuZM62eqwsLU0NQsyBkSExE02wmfyFE9yfhSAghGjgrKqg/coS6Q4fVMHT4MPa8vFbP1cfE4HPNNWoIuuYafK65Bn1kpAyZF6IXkHAkhLgqOauqqD98mLpDh6g/fIT6Q4faDEKGhAR8Bg92hyCfawbJBIpC9GISjoQQvZ7LYqH+6FHqDh6i/tAh6g4dxH6m9aHzhsREfAZfg3nwYHyGDMHnmmtkOQ0hrjISjoQQvYrLZsOalUXdwYPUHzxE/eFDbXaWNsTH4zNkCOYhg901QzJ/kBBCwpEQosdSnE5sp05Rd+Ag9YcOqrfHjoHd3uJcfXQ0PkMGYx4yBJ8hQ/EZfI0MnRdCtErCkRCiR1AUBUdxMXX7D1B/8EBDIDqEq7a2xbm6kBB8hg7BPGQoPkOH4DN4MIbISC+UWgjRE0k4EkJ0S84ai9o/6MAB6g7sp37/ARylpS3O05jNah+hoWmY04biM3Qohrg4GTUmhLhsEo6EEF6nOJ1Ys7Op27+fuv37qT9wAGv2CVAUzxN1Okz9+2MeOhTzsDR8hgzF1DcFjV5+lQkhOo78RhFCdDlHWZlaI5TZEIYOHmy1eUwfG4M5bVhTGLrmGrS+vl4osRDiaiLhSAjRqRSbjfpjx6jbl+muGWptPiGtry8+aWmY09LUIDR0qPQTEkJ4hYQjIUSHspeUUJeZqdYKZWZSf/hwy8VXNRpMqX3xGTYM87BhmNOGYUrtK6vPCyG6BQlHQojLptjt1Gcdo27fvoZAlIm9oKDFebqgIMzDh2MeroYhn6FD0QUEeKHEQghxcRKOhBDt5igvV5vH9u1Tt0OHUOrrPU/SajH169cQhtRAZOzTR0aPCSF6DAlHQohWKS6XOoKsWRhqbTV6bVAQ5uHD8G0IQz5D09D5+3mhxEII0TEkHAkhAHDV1lJ34CB1+/ZSu28fdZn7cVVVtTjP2Lcv5hHD8R0xAvPw4RiTk9FotV4osRBCdA4JR0JcpexFRdTt20ft3n3U7d1LfVYWOJ0e52h8fdVh9CNHqGEoLQ1dcLB3CiyEEF1EwpEQVwHF6cT6/ffU7t1LXUMYaq3jtD4mRg1BI0ZgHjkCnwEDZIJFIcRVR37rCdELqU1kB6jds0cNQ/v346qp8TxJp8NnwADMI0fiO1INRIaYGO8UWAghuhEJR0L0AvbiErWv0J69bTaRaf381NFjI0fgO3Ik5rQ0tH7ScVoIIc4n4UiIHkZxubB+n62Gob17qduzF3t+fovz9LEx+I4Y6Q5Dpv79ZZJFIYRoBwlHQnRzrrq6plFke/e2PopMq8U0YIBaI9QQhqSJ7OqTX1FHucXW5vEQPyNxweYuLJEQPZOEIyG6GXtJidpPaN9eavfuo/7oUXA4PM7R+PpiHpam1gyNGol52DB0/v5eKrHoDvIr6pjx2wysDleb55j0WtY9my4BSYiLkHAkhBcpTmfDRIv71FqhfZnYc3NbnKePjMQ8amRDM9lIfAb2jFFkUpPRdcottgsGIwCrw0W5xSafuRAX0f1/uwrRizhrLNQf2K9OstjWKDKNBlP//k0dp0eMxBAX2+OW35CaDCFETyXhSIhOoigK9pwc6jIzqW1Ypd567Bi4PMOC1tdXXZB1eMP8QsOH9YpFWaUmQwjRU0k4EqKDuCwW6g4eUlen37+fusxMnOXlLc4zxMZiHjlSXYJj5EhM/fr1iCYyIYS4WshvZCEug+J0Yjt5Uh1FduAAdQcOtForpDEY8Bk8uGmF+hHDMURFeanUQggh2kPCkRDtYC8uoe7AfuobwlD9oUO4LJYW5+ljYppWqB82DNM116A1Gr1QYiGEEJdLwpEQ53GcPUv94cPUHz5M3eHD1B88hKO4uMV5GrMZ8+DB+AxLwzw0TWqFhBCil5BwJK5qLYLQocM4iopanqjVYkpNxTwsDZ+hQ9Vaob59pa+Q6DZC/IyY9NqLjg4M8ZOaTCEuRn6zi6uC4nBgO32a+qxjWI8do/5YFtZjx1utEUKjwdinDz6DBzds12AePFjWIRPdWlywmXXPpsu8UkJ0AAlHoldRFAVHaSm2kyexHj+uhqGsLKzZ2Si2Vr40zgtC5iGDMQ0aJLNNdwCpyeh6ccFmCT9CdAAJR6JHUpxO7Hl5WE+cxHbyBNYTJ7GePIHt5Clc1dWtPkfr64upf39MAwfgM2AApgEDMfXvj85faoQ6g9RkCCF6KglHotty2WzY8/Kx5+Viy83FnpOLLS8Pe04OtjNnWq8JAtBqMSYkYOyXis+Age4wZIiPp6DKytnmX9aVDqisdD+UL+uOJTUZQoieqNuHo7feeos33niDoqIihg0bxp///GfGjh3b6rlffvklr7zyCtnZ2djtdvr168dPf/pT7r33Xvc5iqLw4osv8t5771FRUcGkSZNYvHgx/fr166q3JADF5cJ59iz2khIcxSU4SopxlJRgLyzCnpeHLTdX7Q+kKG1eQ2MyYUxOxpSSgrFvCqa+fTGmpGDs06fV4fOynIUQQoj26NbhaOnSpTzzzDO88847jBs3jjfffJM5c+Zw7NgxIiMjW5wfGhrKL3/5SwYOHIjRaGT58uU8+OCDREZGMmfOHABef/11/vSnP/HBBx+QnJzM//3f/zFnzhyOHDmCj49PV7/Fjudygr2uYatFsdfhsNVgs9VQb63CZq/BaqvBaq/FZq/F5rRic9qxK3bsThs2lx27y9lwa8euOHEoCi6tBgcaXBoNTk3DLRqcGnChQUGLwabFVAdGKxhrFUx1Loy1Toy1dgw19Ziq6vGpqMPnnAVjRS1a54WXlgDA7IMmLgZ9fBzGhARMSUmYE/rgk5KCITYWjU7X7o9GlrMQQgjRHhpFucCf5l42btw4xowZw6JFiwBwuVwkJCTwxBNP8Itf/KJd1xg5ciQ33HADv/71r1EUhdjYWH7605/y7LPPAlBZWUlUVBRLlizhrrvuatc1q6qqCAoKorKyksDAwMt7c+dRnE5W/iQdBRcuRcGlUVBQcOJCQcGFus/VsM+BC7tGabrVKNg1qJu2YdNocADKBRYs1bpA7wKdE3RKw60L9E4FnUu9b7SDjx1MdvCxgcmuqI9t6n4fO/hYQXuJ79kFVPjDOX8oD9Bwzh/OBWgoCYbiYA3FIVBtBlopv5/BDz+DH/4Gf/wN/up9o797n5/Bj0BjIEGmIPdtSYWWn3yYheIyg9L23wXLn5jMkLigS3w3Qgghurv2fn9325ojm83Gnj17eO6559z7tFots2bNYtu2bRd9vqIorFu3jmPHjvHaa68BcOrUKYqKipg1a5b7vKCgIMaNG8e2bdvaDEdWqxWr1ep+XFVVdblvq20aDSkbyzr4osp5t53PYQSbD9hNLmw+CnaTot76QL1ZweKnUOOvUOMP1X6aphCn0WAH7BotVoMJjd5IsFaHWaOhDhf1Lgd2xeF+HYvdgsVuoYSSSyqff3/1VnEZUZxmFKcvisMfxemH4vBDcfqzLr+cEmcsoT6h7s3P4IfmAiFTCCFE79Ftw1FZWRlOp5Oo82YcjoqKIisrq83nVVZWEhcXh9VqRafT8fbbbzN79mwAihom92vtmkWtTfzXYOHChbz00kuX+1baRaPVcmJiCBpAiwYtGjSKVr1FgxYtWlDvKRq0Gh06dA23WrTo1VuNDi069TyNDq1Gi1ajQafRodGo13JTFNDr0bg3XcNjg8djrckHra8Zja8vWrMvWl9ftL5mtGYzWl9fNGYzOn9/dEFBaJr39XE6wFoF1mp1qyuHmmKoKTnvtmGzlHGhIOcArOYg6oITqQ1JoiYkAUtQLDV+EdSY/LA4aqmx12CxW6i2VVNtq6bSVkmVtYpKayXn6iuotlWj0ShotDY0WhsYKlu8zntHv4Ojnvt8dD6EmcOIMEcQbg73uN+4RfpGEuoTik7b/qY+IYQQ3c8lh6NNmzYxZcoUtmzZwqRJkzqjTFckICCAzMxMampqWLt2Lc888wwpKSmkp6df9jWfe+45nnnmGffjqqoqEhISOqC0nm78+9YOv6ZX6fTgG6pu7eF0qCGp4gyUn27Ymu7ra4rQ11XiV3cQCg96PlfvA+H9IXIQRAyEyFEQNQSC4t3NcofyK7nxzxtBW49GV4dGV6tu+ho0OgsanQWtvoZxqSbsVHOu/hzl9eXUOmqpd9aTX5NPfk3+hd+yRucOSs23KN8oInwjiPaNJsovCrNe+jQJIUR3dcnh6JtvvkGv17NixYpODUfh4eHodDqKz5vBuLi4mOjo6Dafp9VqSU1NBWD48OEcPXqUhQsXkp6e7n5ecXExMTExHtccPnx4m9c0mUyYTKYreDeiXXR6CIpTt6SJLY/baqEiB86dhNIsdSs5CmXHwVEPRQfUrTnfMIgZDrHDCTT1JxYrBa4wFJcvij2s1WL8z22efY7qHHWU1ZVxtu4sZXVllNWVUVpX6vm4tpSy+jKcipPi2mKKa1uZebuZYFMw0X7R7rAU7RftfhzjH0OkbyQGreFSP0EhhBAd4JLC0UsvvYTD4WDGjBk8+eST/OpXv+KFF17olIIZjUZGjRrF2rVrmT9/PqB2yF67di2PP/54u6/jcrnc/YWSk5OJjo5m7dq17jBUVVXFjh07+MlPftLRb0F0NKMvRA5Ut4HXN+13OdXapcaw1Py29iycWAsn1pIIbPWBs0oAh1zJHFSSOeBKYa+rP2W03QHbrDeTEJBAQsCFawudLidn689SUltCcW0xJbUlHltxbTHFlmJqHbVUWCuosFaQda71JmKtRkuEOYJY/1hi/GLct83v+xp8L+NDFEIIcTGXPFrtvffeo7KykuDgYB5++OHOKhegDuW///77effddxk7dixvvvkmn332GVlZWURFRXHfffcRFxfHwoULAbVv0OjRo+nbty9Wq5WVK1fyi1/8gsWLF7vL+tprr/Hqq696DOU/cODAJQ3l74zRaqIT2Ouh5DAUZEJhJrbcvWhKjmLQOFucetoVxR6lHwcYwKP33UNU3+HQCX2HFEWh2l5NkaXIYyuuLabIUkShpZAiSxF2l/2i1wr1CSXWL5a4gDhi/WOJ948n1j9W3fxi8dH3gqkphBCiA3XaaDWHw8Gzzz7Lu+++e0UFbI8777yT0tJSXnjhBYqKihg+fDjffvutu0N1Tk4OWm3TAHKLxcKjjz5KXl4eZrOZgQMH8vHHH3PnnXe6z/nZz36GxWLhRz/6ERUVFUyePJlvv/22d8xxJDwZfCBulLoBRiC/rJz6vIP4lB7AXHYQ39JMTOeO0UdbTB+K+QGb4R9/A2MAxI+GxPGQMBYSxoHxypcZ0Wg0BBoDCTQG0j+kf6vnuBQXZ+vOUmgppMBSQGFNIQU1BRRZityPq+1qn6hz9ec4dPZQq9cJN4cT7x9PXECceusfR3xAPPH+8UT6RkrHcSGEaEO3nueou5Kao16mrgLyd0PuTsjdAXm7wVbjeY5Wr4asPpPVrYPC0uWqslVRUFOgdhKvzqfAUuDuMJ5fnU+to/aCz9dr9cT5x6mByT/e3WwYHxBPfEA8fgZZb04I0fu09/tbwtFlkHDUy7mcUHIEcrargenMVqjK8zynm4Wl5hRFodJaSX5NPnk1eeptddNtgaUAh8txwWuE+oS6a5kag1NCQAKJgYmE+YTJnE9CoC5JJAsr9yxdGo5ycnLYuHEjJpOJESNGuEeL9VYSjq4yiqJOL3B6s7qd2tRKWDKozW8p6ZAyHWJHqKPvuiGny0lJbQl5NXnkVeeRW53rvp9XnUe5tfyCzzfrzcQHxJPgr4alxhqnxIBEYvxipLlOXBVkrcaeqcvC0Z/+9CeeeeYZfH190Wg01NTUMGrUKP7617+SlpZ2JZfutiQcXeXaE5ZMQZA8BfpOV8NSaEqry6B0R9W2ajUo1ajBqXHLq86j0FKIS2n7y0Cv1btrmxqDU2JAIomBicT6x8r0BKLXUOdN23zR82Q5ou6ly8JRREQE/+///T9+9atfodVqyc7O5u233+bdd9/lu+++Y/LkyVdy+W5JwpHwoChQfgpOZsCJ9XBqA9SfN/N2UCL0TYe+M9XaJXNw15ezA9iddvJr8smtziWnOoe86jxyqnPIqcohvyb/gqPsdBodMX4xJAYmugNTUmCSWvPkH49BJ8FJ9BwSjnqmLgtHZrOZgwcPtmhK+81vfsPXX3/Njh07ruTy3ZKEI3FBLqc6fcDJdXAiQ+3k3Tw0aHRqE1zqLHWLTgPtpS7b2/04XeoEmI1hKbc6l5yqHHeIqnfWt/lcrUZLjF+MOywlBSaRFJhEYkAicf5xEpxEtyPhqGfqsnA0btw4Hn/8ce69916P/cePH2fYsGHU1dVdyeW7JQlH4pLYLGqn7hPrIHuNOqN3c36RDUFpJvSd0f7lVnoQRVEoqS0hp1oNTWeqznjc1jna/j3RWOPUPDg11jpJU53wFglHPVOnzXN0vt/97nfceuutGI1G7rjjDvcolh07dtCvX78rvbwQPZ/RD/rNVjcWquvFZa9Rt5MbwFIC+z9RN40W4kZD/znqFjWkx/RVuhCNRkOUXxRRflGMiR7jcUxRFMrqyjwCU2PtU051DnWOOrXDeE1ei+vqNDri/ONICEwgKaApNCUFJBHjH4Ne2z07xQshurcOGa22cuVKfvzjH1NfX8/w4cOx2WwcOnSIjz76iOuuu64jytmtSM2R6DAOK+RsU4PS92ug9Kjn8cA46Hct9J8LyVPVJVSuIoqiUFpX6hmcqnI4U32G3KrcCzbVNXYOb+zj1LzGKdo3WkbViSsiNUc9U5fPc2S1Wlm7di0ZGRns27ePffv2UV5eTmhoKEOHDiUtLY0333yzI17K6yQciU5TmQffr4Lj36m1Ss2bm/Q+akDqd61aqxSc6L1ydgMuxUVJbYlnaGpW62RztT3/jFFrdI+oc4emhpqnSN9ItJqe3wdMdC4JRz1Tt5gEMjc3l8zMTHdYWrZsWWe9VJeScCS6hL1OnSbg++/g+CqozPE8HjUEBlwH/a9T51XqBZ26O4pLcVFsKeZMdbPQ1FjjVJ17wUkwfXQ+LZrpGmuews3hMgGmAGSeo56qW4Sj3krCkehyigIlRxuC0nfqCLjm8w35RzX0U7pOnSrgKmt+uxQOl4NCS6FHTVNjeMqvyceptFyYuJGv3rfVZrrEgERCfUIlOF1lZIbsnkfCUSeScCS8rvac2vx27BvIXgu26qZjeh81IA24Tu2rFBDttWL2NHaXncKawhah6UzVGQosBRecANPf4N9iDqfG+yGmEAlOQnQDEo46kYQj0a04bHBmMxz7Vg1L5ze/xY1Wg9LAGyBiYK8Y/eYNdqedvJq8VmucCi2FKLT9qzTAENDq5JdJgUkEm4IlOAnRRSQcdSIJR6LbUhR10dxjK9WglL/H83hIHxhwgxqWEid02/Xfehqr06rOFt4w/UDzPk5FlqILPjfAGNBqbVNiQKIEJyE6mISjTiThSPQY1UVqSDr2jbq8idPadMwnWO2nNOB6dQJKU4C3Stmr1TvqyavO8+gc3jjCrri2+ILPDTAEuDuHJwQ2rVOXEJBAmE+YBCchLpGEo04k4Uj0SNYaOLleDUrHv4Xas03HdEZIngYDr1fDkvRT6hJ1jrqWNU4NUxFcLDg1dg5vvrhv4/0I3wiZjkCIVkg46kQSjkSP53JC7k44tgKyVsC5k57H40Y3BKUbIGKA9FPygsYap+ZLruRU55BblXvRPk4+Oh/iA+JJDEh0z+eUEJBAQkAC0X7RMnO4uGpJOOpEEo5Er6IoUHqsISithPzdnsdDU9TapIE3QMI4kJmlvc7mtJFXk0duVa7nQr/VORTUFFxwOgK9Vq8uudIQlhprmxICEogLiMOkM3XhOxGia0k46kQSjkSv5u6ntLKhn1KzeVx8w9TpAQZcry6SK/MpdTuN0xE01jg1BqfGze6yt/lcDRoifSM9apriA+Ld9wON8vtO9GwSjjqRhCNx1bBWq/MoHVupTj5ZX9F0TO8DKdPV5rf+14F/hNeKKdrH6XK6l1zJqVb7OeVV57lDVK2j9oLPDzIFkeDvGZga78uyK6InkHDUiSQciauS064ukpu1Um2Cq2g+n5IGEsaqNUoDroeI/l4rprg8iqJwrv6cu4apMTQ1bmfrz17w+UatkVj/WHdYivePJz6gYfOPx9cgtYzC+yQcdSIJR+KqpyhQfFitUcpaAYWZnsfDUpuCUsJY6afUC9Taa9V+Ts2CU+NtQU0BDqXt9eoAQn1C3UGp+W2cfxxRvlHo5P+I6AISjjqRhCMhzlOZ1zSf0qmN0Lxfi0c/pelg9PNeOUWncLgcFNcWk1edp27NQlReTR6V1soLPl+v1RPrF0ucf5w7MMUFxKkByj+eIFOQzOkkOoSEo04k4UiIC6ivghNr1ea377+D+mZfjHofdT6lxnXfAmO8V07RZapsVe5apvyafPKr88mryVPv1+TjcF241snP4KcGpobNHaAaNmmyE+0l4agTSTgSop0u2E8JiB2h1ij1nwvRQ2U+pauQ0+WktK7Uo6apMUDl1+RTWld60WuEmEKI848j1j+WuIA44vzUmqdY/1hi/WLx0ft0wTsRPYGEo04k4UiIy6AoUHK0YeTbt5C3G5pPZBgYDwPmqrVKfaaAXubbEepkmAWWAvKqPUNT41Zlq7roNcJ8wtzhKdZfbb6L8YtRb/1jMOvNXfBORHcg4agTSTgSogNUF6vNbse+hRPrwFHXdMzgp/ZP6j9XXf/NP9J75RTdWrWtmoKaAvJq8iioKXCHpsb7FrvlotcI9Qkl1i+WGP8YYv1i3SEqxi+GWP9YAoyy7mBvIeGoE0k4EqKD2evg5Iam+ZRqzlvJPm6UOpdS/znS/CbaTVEUqmxVFNQUuMNSgaXAHZ4Kagqosddc9DoBxgA1PPnFuANUtH+0O0iF+oTKHE89hISjTiThSIhO5HJB0X41JB3/Fgr2eR4PjFNDUr85kDxVZukWV6TSWkmhpZCCmgIKLYXk1+RTWFNIgUUNTxXWiotew6g1Eu0XTYx/DDF+Mer9826l6a57kHDUiSQcCdGFqgrh+1VqWDq5HuzNZnHW+6j9k/rPgX7XQkiS98opeqVae61HeDr/fkltCS7FddHrhJhCiPaL9ghMze+Hm8NlQeAuIOGoE0k4EsJL7HVwerNao3R8FVSeN/otfAD0v1atVUocDzqDd8oprhp2l53S2lJ3YCqyFLmDU5GliIKagosuywKg0+iI8I0g2jfaHZyi/aKJ8o1y35fmuysn4agTSTgSohtQFCjNUmuUvl+tThnQfDV6UxD0TYfU2ZA6S+ZUEl6hKArV9moKa5qCU/Pb4tpiii3FF51hHNTJMqN8o9TNL4po32ii/KI89oX5hMls4xcg4agTSTgSohuqq1BHvX2/Sg1LtWWex6OGQr/Z6hY/FnTShCG6B6fLydn6sx6hqXlwKqotoqyurF3Nd401UI2BKdI3kmi/aCJ9I91blG8URp2xC95Z9yPhqBNJOBKim3O5oGCvGpKyV0P+XjzmVJJaJdHD2F12ztaddQcnd3hq3CzFlNWV4Wxee3oBIaaQprDkF0WkWb3fGKwifCMIMYX0umVbJBx1IglHQvQwlrKmWqXstVB3zvN45GBInQF9Z0LiBDDIjMqi52msgSq2qIGpyFJESW2JO0CV1JZQUluC1Wlt1/UMWgMR5gh3aIr0jfR4HGGOIMI3ggBDQI8JURKOOpGEIyF6MJdTnR6grVolvRn6TIbUmWpYCu8n8yqJXqNx7qfG4OTe6tTb0tpSimuLOVd/7uIXa+Cj83GHpXBzOBG+6m2kb6T62Kwe6w4LCEs46kQSjoToRSxn1SkCTqxTa5XOn4AyKFGtVUqZrs6r5BvqnXIK0YXsTjtldWXu0NQYnErrSj3ut2f5lkYGrcEdlsLMYe4wFe4bTrhPU6gK8wnD0EkjTSUcdSIJR0L0UooCxYfhxFo1KOVsA6et6bhGqy6WmzJdXd4kfizor86OrUKAuvZdaV0ppbWllNSVUFZb5n5cWldKWZ36uNJaeUnXDTYF890PvsPX0LGTvEo46kQSjoS4Stgs6rxKJxpqlsqOeR43+EGfSdC3oWYpYoA0wQnRCpvT5g5KZXVllNWWUVZfRmltKWfrzrqPna07i0Nx4KPzYec9Ozu8GU7CUSeScCTEVaoyH05mqM1wJzPAUup53D9abXpLmQbJ0yA4wRulFKLHcikuKq2VVFgrSA5K7vDrSzjqRBKOhBC4XFByWK1ROrFebYJz1HueE5qihqSUadBnKviFeaesQghAwlGnknAkhGjBXg95O+HkBji1QR0Fd/6cM9FD1bDUZwokTQCfIO+UVYirVHu/v7v9Ii1vvfUWffr0wcfHh3HjxrFz5842z33vvfeYMmUKISEhhISEMGvWrBbnP/DAA2g0Go9t7ty5nf02hBC9ncFHbVKb+X/w8Br4+Sm4658w7scQMUg9p+ggbFsE/7wTXusDf5kOq19QpxWwVnu1+EKIJt265mjp0qXcd999vPPOO4wbN44333yTzz//nGPHjhEZGdni/HvuuYdJkyYxceJEfHx8eO2111i2bBmHDx8mLi4OUMNRcXEx77//vvt5JpOJkJCQdpdLao6EEJesuhhOb4JTG9Xbcyc9j2t0EDdSrVXqMxkSxoHJ3ztlFaKX6hXNauPGjWPMmDEsWrQIAJfLRUJCAk888QS/+MUvLvp8p9NJSEgIixYt4r777gPUcFRRUcFXX3112eWScCSEuGKV+WpIOr0JTm2CijOexzU6iB0OSZOawpI52BslFaLXaO/3d7ddedFms7Fnzx6ee+459z6tVsusWbPYtm1bu65RW1uL3W4nNNRz0raMjAwiIyMJCQlhxowZ/OY3vyEsrO2OklarFau1abr1qqr2T3olhBCtCoqDYXepG0BFjhqSTm+C01ugMgfy96jb1j8BGrXPUtIkdfqAxInSwVuITtJtw1FZWRlOp5OoqCiP/VFRUWRlZbXrGj//+c+JjY1l1qxZ7n1z587l1ltvJTk5mRMnTvD8889z3XXXsW3bNnQ6XavXWbhwIS+99NLlvxkhhLiY4EQYcY+6gRqWzmxV51k6sxXOnYCiA+q2Y7F6TvgASByvrgeXNAGCk2SeJSE6QLcNR1fq1Vdf5dNPPyUjIwMfn6ZFJO+66y73/aFDh5KWlkbfvn3JyMhg5syZrV7rueee45lnnnE/rqqqIiFB5i8RQnSi4ER1a6xZqiqEM1vU7fQWdULKxm3vB+o5ATENYWmiehs1GLSt/9EnhGhbtw1H4eHh6HQ6iouLPfYXFxcTHR19wef+9re/5dVXX2XNmjWkpaVd8NyUlBTCw8PJzs5uMxyZTCZMJtOlvQEhhOhIgTEw9DZ1A3VNuNwdkLMVcrari+lWF8LhZeoGYAyA+NFqf6XEcRA3Gnykn6QQF9Ntw5HRaGTUqFGsXbuW+fPnA2qH7LVr1/L444+3+bzXX3+dl19+me+++47Ro0df9HXy8vI4e/YsMTExHVV0IYTofH5hMPB6dQOw1ar9k3K2qxNS5u4EW3XDbN7rG56kUWuTEsZCwnj1NqSPNMUJcZ5uPVpt6dKl3H///bz77ruMHTuWN998k88++4ysrCyioqK47777+P/bu/eYps4+DuDftiDwQqkitxaQqZsXVKggMqc455yGGDKTLXHOTNyW/TNcZni3xWWJYGaiy7LEeWPGP+b+MeqWwJJt6ozLvCwaEaxSfFW8I7RcFUqVa/v+8ZRyKuAE2556+H6SJ6c9PT399STSr895znMSEhKwZcsWAMA333yDjRs3Yv/+/Zg/f757PxEREYiIiEB7ezs2bdqEt956C/Hx8bhx4wa++OIL2Gw2VFZWPnXvEK9WI6KA19sDNFwWvUs158Ty8SviACA8FkjMFD1MiZnixrqcQoAU6rm/Wg0AVq5cicbGRmzcuBFWqxVGoxFHjhxxD9K+e/cu1Or+eSyLi4vR1dWFt99+22M/hYWFKCoqgkajwaVLl/DTTz/hwYMHMBgMWLp0Kb7++mueNiMiZdEEAfpU0eZ+JNbZrP1BqeYcYDEB9gbg6u+iAYBKDcTOcIUlV2Aa/xKgDvg5g4m8JqB7jgIVe46ISBG6OwDLRaD2PHCvDLh3HmitGbhdiA5ImA0kZIhmSBdjoIieM4qYBDJQMRwRkWK1WTzDUm0F0PNo4HZag5jROyHdFZhm815xFPAYjnyI4YiIRo3eHqChSoSk2nKxbPwf4HQM3Hb8i4DeKIKSwQjo04AQrb8rJhoSw5EPMRwR0ajWZXedjivvD0yDDfaGCoh+yTMwxadywDfJhuHIhxiOiIgeY28C6kyA5YJY1pmAtnuDbKgCxk8WIUmfJgaMx6fxVijkFwxHPsRwRET0FNobxRVxdSYxSWXdBcBWN/i2kQkiLMW7rrCLmylmCOccTORFDEc+xHBERDRC7Y2A9SJgcd0nznJJ3DduMCE6MWll/CwgfqYITLHTgeAw/9ZMisFw5EMMR0REXtTRBtSbRVCyXASslUDjFcDRPXBblVrMuxQ/UwSn2BlAXAqgS2IvE/0rhiMfYjgiIvKxni6g6ZoITdZK0erNwMPmwbcP0YlepbgUIDalv5cpbKxfy6bAxnDkQwxHREQycDrFLN99ganhMlB/GWi6Cjh6Bn+P1iBCUux0IGaaazmVUwyMUgxHPsRwREQUQHq6gOZqEZQaqsSyvmqIq+VcdEn9gSlmmghM0VOAUP5NVzKGIx9iOCIieg50tAINV8SkldJlu3Xo92j1rqA0FYiZ4lpOA8KjOaZJARiOfIjhiIjoOfawRQz4bvifWDZeARqvPTk0hY4VE1qOf0ks+x5HTQKCxvitdHo2DEc+xHBERKRAjx4ATdViDFPjVTEgvPEKcP8OgCF+KlUaYFxyf2gaPxmImixupaLVA2q1P78B/QuGIx9iOCIiGkW6HwEtN0VYarouls3VIkh1tQ/9vqAwV1iaJMKSOzhNBsJjeJpOBk/7+x3kx5qIiIieP8FhYk6luBme6/uunmuudgWmG0DzdbG8fxvoeSSurKs3D7LPcBGaoia6lpLHWgN7nGTGcERERDQSKhUQqRdt4kLP13q7gQd3+8NSS19wugm01gDddqC+UrTHaULEqbpxLwzexoT7+puNegxHRERE3qYJFqfPxk8e+FpPpwhOLTdd7Vb/4wd3gN5O1ym8a4PvOzxWhKexya7lBPF47AQxRQEHiD8zhiMiIiJ/Cgrpv+Ltcb09Yn6m+7cHtpZbQMcDwN4g2r2yQXauAiINksCUJALT2CRANwHQJfDedE+B4YiIiChQaIL6T58N5tF9cfXc/duil+nBXdHuux73PALaakW7e2bwfYTHSAJTX0sAdIlAZCLndALDERER0fMjbJxoBuPA15xOwN7kCk13RGBqrQEe1PQvu+2AvVG0uorBP0MTIsJSpCsw6RLF48gE0SsVaRA1KDhAMRwREREpgUoFRMSIljhn4OtOp+h5ejwwtdaInqbWWqC9Xox56hsDNZSgsP6g5G4JgDZeXG2njQci4kRP2HPo+ayaiIiIhkelAv4TJZo+bfBteroAW50ISq33xPinvse2OqCtDnjYLE7ftbiuwhv6A4GIWDEZplYvAlOkQYSmvvCkjRen+dQan3zlkWI4IiIiIiFozJPHPAFAdwdgs4ig1FYnep1sFleAsorWbgUcPaInqr0esJiG3p9KLQLS46Ep+7+yDR5nOCIiIqKnFxzqmrBy4tDbOBzAwyZXiLKIpc0qep9s9SI82erFVXdOR3+Isl4S71epgUVf+uf7DILhiIiIiLxLrRan1CJihz6FBwCOXjE43GYV4ahv2WmT9VQbwxERERHJQ61xDeKOl7sSD7x5CxEREZEEwxERERGRBMMRERERkQTDEREREZEEB2SPgNPpBAC0tbXJXAkRERE9rb7f7b7f8aEwHI2AzWYDACQlJclcCREREQ2XzWaDTqcb8nWV89/iEw3gcDhQV1cHrVYLlRdvvNfW1oakpCTU1NQgMjLSa/ulwfF4+xePt3/xePsXj7d/jfR4O51O2Gw2GAwGqNVDjyxiz9EIqNVqJCYm+mz/kZGR/MflRzze/sXj7V883v7F4+1fIzneT+ox6sMB2UREREQSDEdEREREEgxHASQkJASFhYUICQmRu5RRgcfbv3i8/YvH2794vP3L18ebA7KJiIiIJNhzRERERCTBcEREREQkwXBEREREJMFwRERERCTBcBQgdu3ahRdeeAGhoaHIysrCuXPn5C5JsU6ePInc3FwYDAaoVCqUlpbKXZKibdmyBZmZmdBqtYiNjcWKFStw9epVuctSrOLiYqSmpronx5s3bx4OHz4sd1mjwtatW6FSqbB+/Xq5S1GsoqIiqFQqjzZt2jSvfw7DUQA4ePAgCgoKUFhYiIqKCqSlpWHZsmVoaGiQuzRFstvtSEtLw65du+QuZVQ4ceIE8vPzcfbsWRw7dgzd3d1YunQp7Ha73KUpUmJiIrZu3Yry8nKcP38eixcvxptvvomqqiq5S1O0srIy7NmzB6mpqXKXongzZsyAxWJxt9OnT3v9M3gpfwDIyspCZmYmdu7cCUDcuy0pKQmffPIJNmzYIHN1yqZSqVBSUoIVK1bIXcqo0djYiNjYWJw4cQILFy6Uu5xRISoqCt9++y0+/PBDuUtRpPb2dqSnp2P37t3YvHkzjEYjtm3bJndZilRUVITS0lKYTCaffg57jmTW1dWF8vJyLFmyxL1OrVZjyZIlOHPmjIyVEflGa2srAPGDTb7V29uLAwcOwG63Y968eXKXo1j5+flYvny5x99x8p3q6moYDAZMmjQJq1evxt27d73+GbzxrMyamprQ29uLuLg4j/VxcXG4cuWKTFUR+YbD4cD69esxf/58zJw5U+5yFKuyshLz5s1DR0cHIiIiUFJSgpSUFLnLUqQDBw6goqICZWVlcpcyKmRlZWHfvn2YOnUqLBYLNm3ahOzsbJjNZmi1Wq99DsMREflNfn4+zGazT8YIUL+pU6fCZDKhtbUVv/zyC/Ly8nDixAkGJC+rqanBp59+imPHjiE0NFTuckaFnJwc9+PU1FRkZWUhOTkZhw4d8uppY4YjmUVHR0Oj0aC+vt5jfX19PeLj42Wqisj71q1bh99++w0nT55EYmKi3OUo2pgxY/Diiy8CADIyMlBWVobvv/8ee/bskbkyZSkvL0dDQwPS09Pd63p7e3Hy5Ens3LkTnZ2d0Gg0MlaofGPHjsWUKVNw/fp1r+6XY45kNmbMGGRkZOD48ePudQ6HA8ePH+cYAVIEp9OJdevWoaSkBH/99RcmTpwod0mjjsPhQGdnp9xlKM7rr7+OyspKmEwmd5szZw5Wr14Nk8nEYOQH7e3tuHHjBvR6vVf3y56jAFBQUIC8vDzMmTMHc+fOxbZt22C32/H+++/LXZoitbe3e/wv49atWzCZTIiKisKECRNkrEyZ8vPzsX//fvz666/QarWwWq0AAJ1Oh7CwMJmrU54vv/wSOTk5mDBhAmw2G/bv34+///4bR48elbs0xdFqtQPGzoWHh2P8+PEcU+cjn332GXJzc5GcnIy6ujoUFhZCo9Fg1apVXv0chqMAsHLlSjQ2NmLjxo2wWq0wGo04cuTIgEHa5B3nz5/Ha6+95n5eUFAAAMjLy8O+fftkqkq5iouLAQCLFi3yWP/jjz9i7dq1/i9I4RoaGrBmzRpYLBbodDqkpqbi6NGjeOONN+QujeiZ3bt3D6tWrUJzczNiYmKwYMECnD17FjExMV79HM5zRERERCTBMUdEREREEgxHRERERBIMR0REREQSDEdEREREEgxHRERERBIMR0REREQSDEdEREREEgxHRERERBIMR0REREQSDEdERC59NwwoKiryeP4kNTU1WLRoEVJSUpCamoqff/7ZlyUSkR/w9iFERC67d+9GUFAQqqurodFokJOTg1dfffWJ77FYLKivr4fRaITVakVGRgauXbuG8PBwP1VNRN7GniMiUpzm5mbExsbi9u3bw3rfxx9/jNbWVmzfvh25ubkDgtE777yD7777zmOdXq+H0WgEAMTHxyM6OhotLS3PUj4RyYw9R0SkOAUFBbDZbNi7d++w3vfDDz9ArVa7e46WL1+O7Oxs9+tmsxkLFy7ErVu3oNPpBry/vLwceXl5MJvNz/wdiEg+DEdEpCgPHz6EXq/H0aNH8fLLLw/rvU6nEyqVCkVFRSgqKnI/l8rMzMTatWuRn5/vsb6lpQXZ2dnYu3cvXnnllWf+HkQkH55WIyJF+eOPPxASEuIRjE6fPo3g4GB0dHS4192+fRsqlQp37txxr+sLQn0Dsh8PRgCQm5uLAwcOeKzr7OzEihUrsGHDBgYjIgVgOCIiRTl16hQyMjI81plMJkyfPh2hoaHudRcuXMC4ceOQnJw8rP3PnTsX586dQ2dnJwDR27R27VosXrwY77333rN/ASKSHcMRESnKnTt3YDAYPNZdvHgRs2fP9lhnMpmQlpY27P0bDAZ0dXXBarUCAP755x8cPHgQpaWlMBqNMBqNqKysHPkXICLZBcldABGRNz169MijhwgQQejdd9/1WHfhwgX3VWbDERYWBkCMbQKABQsWwOFwjKxYIgpI7DkiIkWJjo7G/fv33c97e3thNpsH9BxVVFSMKBz1XaYfExPzTHUSUeBiOCIiRZk9ezYuX77sfn716lV0dHR4nGo7c+YMamtrRxSOzGYzEhMTER0d7Y1yiSgAMRwRkaIsW7YMVVVV7t4jk8kEANixYweqq6tx+PBhrFmzBgDQ1dU17P2fOnUKS5cu9Vq9RBR4GI6ISFFmzZqF9PR0HDp0CIAIR8uWLcPNmzcxa9YsfPXVV9i0aRMiIyOxffv2Ye27o6MDpaWl+Oijj3xROhEFCE4CSUSK8/vvv+Pzzz+H2WxGTk4OMjMzsXnz5mfeb3FxMUpKSvDnn396oUoiClS8Wo2IFGf58uWorq5GbW0tLl68iA8++MAr+w0ODsaOHTu8si8iClzsOSIixbJardDr9aiqqkJKSorc5RDRc4LhiIiIiEiCA7KJiIiIJBiOiIiIiCQYjoiIiIgkGI6IiIiIJBiOiIiIiCQYjoiIiIgkGI6IiIiIJBiOiIiIiCQYjoiIiIgkGI6IiIiIJBiOiIiIiCQYjoiIiIgkGI6IiIiIJP4Pc6kZU/WY9pEAAAAASUVORK5CYII=", "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.0" } }, "nbformat": 4, "nbformat_minor": 5 }