{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "0052cdd2",
   "metadata": {},
   "source": [
    "## M3GNet-UP-2022 cohesive energy (naphthalene)\n",
    "\n",
    "This PLAMS workflow mirrors the GUI tutorial and uses the same structures, settings, and steps:\n",
    "\n",
    "- Gas-phase naphthalene geometry optimization\n",
    "- Crystal naphthalene lattice + geometry optimization\n",
    "- D3(BJ) dispersion add-on\n",
    "- Cohesive energy calculation from the two energies"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5fc5ed7b",
   "metadata": {},
   "source": [
    "### 1) Gas-phase naphthalene structure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "bc98730c-5a94-4ff5-aef3-47b0ec4fb538",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:03.859137Z",
     "iopub.status.busy": "2026-02-11T10:35:03.858989Z",
     "iopub.status.idle": "2026-02-11T10:35:06.416989Z",
     "shell.execute_reply": "2026-02-11T10:35:06.416481Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADICAYAAACtWK6eAABN/UlEQVR4nO29eXxdxXn//5k559x9369sybIsb5JtwMY2AYJMbAgEmhWZLaSBNhBISAhNnLYkyGJpA+RH0m8IBZKGFAoxFkkAA04xLXbYYsBgvAjbspE3bXff13Nmfn/onMu1MSSlsiXR83699HpZ915dz3NmnplnZp4F0NHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR0dHR+ZhDxrsBExzS1dVFAKC7u5sD4OPcnrGiJtfq1as5IeTjIpfOCYKsXbtWOPpFzrnAOZ/Mk8rHVS6dE0VXVxfV/r1o0SLLj370o/YHHnhgOoDawJqkg6m+zcJdd901/Z577mm/+uqrLdqLnHN6jL/7P81k7OjjCQXAVqxY0XD22Wd3t7W1nRsKhab6fL6KzWZ7t1Qq9Tz11FN3f/vb385wzskkMk0IAL5kyRLHpz/96Rvnz5/fOWXKlBaPx2NwOp2HS6XSc9u2bev6/Oc/P8g5p4QQNt4N1plgaCvH4sWLT/27v/u7wWeffZb39fXxYrHI64lEIptfeOEFH+ecTpKVhHR1ddFZs2b5vv71r29+/PHH+c6dO3k2mz1CrmQyObh58+ZTAX0l0Xk/hHNOnE6n+4ILLnj33nvv5S+++GI5Fosxzrn2o3DOy5xzfvDgwf8EJsdA0hS/o6PjP++66y6+YcOG8sDAgKLJpSgK0+SKxWLvPv30027OOZkkyn/cmfAdfCLo6OgQCCHc6XRe7fF4phsMhqokSQaLxUIwap4QAJRzbgAgh8Phc1988cVzCCGMc/6+Te8EQuju7mYej+cch8NxrsFgkAVBMFgsFgp1UqCUEgAGAFWv1zt9ypQpV6um40SW64ShKwiAZcuWMQAghJwLgFcqFVqtVqEoyhGfI2R0UhVFkRuNxhXayyeyrf8TOjo6CAAYDIYVgiBwRVEgyzJkWT7WxykALoriuerv+j4EuoIAALq7uwEAiqIgl8uRbDaLdDqNZDJ5rI+TarVKksnkKervE3Ygbdq0iQEA5/yUXC5HMpkMSafTSCQSYIzVFL4OUqlUTnxDJzC6gtRRqVTEaDTKo9EohoeHcfDgQUQikSNWklKpxEdGRngymdwLABs3bpzIz5ACgKIoe+PxOI9EInxkZASHDh3C4ODgEStJpVJBJBLhhw8fFsettRMQ/WGMQgEwQshbQ0NDZ+3fv59brVZIkoRKpYJQKASbzQbGGDKZDEZGRkhfX9/LALBx48bxbflfAKX05Xg8fm1/fz8cDgeMRiMURUEymYTT6QQA5HI5nkqlyJ49e94Caoo/YVfHE4WuIKMwAEQUxZ/ncrmv7dq1yyRJkqIoipDL5RCJRGA2m0EprZZKJWn79u19jz766OPqXYjyZ799/FAAEIfD8Xg2m+3au3fvTKPRWGWMSYVCAbFYDFarFYIgKJVKhe7Zs6ewbdu2nwMgGzdu/D+vHICuIBoMAB0YGOhraGi4Lh6P/3rr1q1IpVJ8ypQpxO12w2QycUqpNDQ0lDx48ODle/fuLa9evXqiz7IcAN27d2956tSpl2ez2f/csWOHO5vN8pGREeL1emG1WrkgCEIsFsOBAweuW7duXV9XVxft7u6eyHKdMCbsCcw4IQBQwuHwpxlj95jN5hk+nw8ulwuCIPBsNvvHgYGBGw8cOPAWVLNsnNv7l0IBsIaGhlMYY3dLknSW1+slbrcbBoMB+Xx+XzQa/ebu3bv/s7OzU+jp6ZnIq6LOOGMAAIPB8I+f//zn+UMPPVTdvHmz8swzz/Curq4lAHD//fdL49vEj4TW5iVnnHEGv+eee5SXX365+vLLL/Nbb731HwFg7dq1hnFs34RkIp/AjBcMgKAoilEQBFitVjgcDthsNlBK7QCo2+2eLCtHPUy9VbcTQmA2m2G322G322E2m42dnZ2C3++fjHIdV3QFOTYKAC7LMorFIvL5PPL5PEqlEsPkMavexy233MIAMEVRUCqVUCgUkM/nUS6XuW5WHRt9k/4hMMYgyzKq1Sqq1eoH3UBPOhhjNZk45zW5JsOR9YlGV5APQVEUVKtVVCoVVCqVj5WCyLKMcrkMAB8buY4HuoJ8CNpMW6lUjlhBenp6xrll/zvq5apfQXTej64gH4K2gpTL5Y/dCqIpCKCvIB+GriAfQr2tXqlU3ufdO1nhnB+xgugOih+MriDHhlBKiaYg2gpSLpcn/amfIAi0fgVhjEFRFC3mRecoJn2HjzWtra0CAF6tVsuyLENRlNqmljFW4JwTt9s9GZ8bZYwRRVEKWqwLY0z7KQPgiURCD5I6isnY0ccTunfv3jIAqbGxcYYkSUgmk2RwcBCJRALVanUeIYQ/8MAD1c7OzkkzmNS2VgkhvLm5eZ7FYkE2m8XQ0BAZHh5GsVicAUD62c9+Vq7P6qKjL6s1NAe9U0899RKXy/UDp9PZbjKZuCRJRBRFCILACSGEMbYhk8n8cM2aNZsng9+S1sYFCxYs9Xg8tzocjnOsVis3GAw1uQRBIIqi7CyXy7f9+7//+xrdWfE9dAUBSGdnJ1UH0X2BQOAap9MJs9nMrVYrsVqtMBqNEAQBiqJwACSbzcrpdPq6hx9++BcTWUm0ts2ZM+drXq/33kAgIJrNZm6xWIjVaoXZbIYoimCMcc45qVQqSCQS9z/44INfV/+W4eOTTfIjcaIUhHDOMRHzSGmDqKWl5e6GhobvOJ1O2WazUYfDQR0OB5xOJ7TgKVmWUSgUlEKhQPP5PIlGo+c8+uijz09QJREAKE1NTSv8fv8Gn8/HbTYbczgcgiaXzWaD0WgE5xylUonlcjkmy7I4MDDwk4cffvjGCSoX1Dgc4AQo73G1Nzs7OwU11SUnhHDOOXnhhRfECZRShvb09CihUOhUk8n0LVEUZUEQBIPBQE0mE+x2O9xuN/x+P4LBIILBIHw+n2Cz2ZjVauVms/nnHR0dprVr1zJMrNWYcM4ZAKPBYPi52WzmgiAwg8EgGI1GWK1WuFyuo+WiTqdTEARBdjgc37rkkktO7enpUSbQnoR0dXWJdQn7OOdcOFYq1bHkuAmvzT4rV65UANCf/OQnLkIIP/vss2VVWcb9wS9atEjo6uqigiDcYDQaBUopKKWEUgpJkmAymWCz2eByueDxeODxeOByuWCz2QRBEJjD4ZjV2Nj4pZUrV9Krr756Ih2ZC4QQHgqFVhoMhlmEEEYpFerlslqtcDqd8Hg88Hq9cLvdcDgcRJVZMBgMN3R1ddGhoaFxP4xQlZR3d3fLhBDe1dXlAkAJIcrKlSuV45l66bh0qhaKeu2116449dRTL29qalrs9XpD+/bt2yvL8h82btz4E0JIehw3g6Sjo0PYtGlTdcuWLWhsbFwGdUbifHTVJoRAEASIoghJkmr7EKPRCIPBAFEUQSnlnPPLe3p6HgGgdHR0iJs2bVIw/nY7BwBK6ZcopdpsW3uTUlqTy2AwgBCCSqWiySUwxrggCMvUvmHjKZc2RlasWOE8++yzvzN79uzzwuFw62WXXTZsNBpfj0QijxBCnj9eqWDHXEG6urooIYR96Utf+kl7e/sNJ510Epqbm+HxeEAI8QJYGggELp0/f/4Vp59++mvjkAtWAKBs2rRJDgaDAb/ff3M2mw3KsswVRaGKokDLH6VdppXLZYiiWHtdvUMQKpUKFwThUxdffPGvBgYGfrRp06Y9QG31HJcNrvp/cwBQFGWWevdB1QvB98mlKYjma8YYI7IsM0pp8KKLLrpnYGDglk2bNkXqvvuE7Uk05TjttNOWLFq06OFFixbNmjVrFsLhMEwmkxdAezgc/mpfX99PCSHfOR5jaUzNnM7OTqG7u5stXLjwG7NmzbqhqalJcTqdssvlYqp2MwBVl8s1a9asWes2bNjgxejsdiLsd9rR0SFiNNZDamtru37mzJnbm5ubvwFArFartH7wlEol5HI5ZDIZJJNJJBIJpFIppNNpFAoFbYARg8FgbG5uvnLmzJmvXXTRRV1Lly4NqoOIj8ddifp/s0WLFl1gMplaK5UKUxSF1itFoVBAJpNBKpVCIpFAMplEOp3WYl6gPguxqanpG7Nnz97+uc997noAUk9Pj9LR0SGeoH0JWb16Nbfb7d6GhoZ1s2fPnuV0Oqtut5uZTCZtLMkGg0FpbW29obe39xvHI9PlWApK1q5dy6xWa8Bms93qdDoVURSJ1WoVBUGgUNN3YjT0s+r1egPTpk375glIc0lUxWCbNm2SZ82atXzJkiUvNjc3/7/GxsZAKBRSLBaLNjBqypHP55HJZBCPxxGJRDA8PIzh4WHEYrHaYCqXyyiXy7xSqcjBYNA5bdq01bNnz96mDiiDtsk93orS1dWlJdKmCxYsuPATn/jEtMWLFx/Q5NLc9bUgKU3pI5EIRkZGMDw8jGg0imQyiXw+j2KxqPmdKaFQKNDc3Pz/Lr/88hcvvPDC5Zs2bZK7u7uZ+kyP28SmpYMNh8PfDAaDAUmSqmazWbJarfVjSWSMEQCK3++/9aGHHgoAYGM54Y6lglBCCDcajWcZDAa3aqrQetu3DoExxgkhV9x1111WAMrxWEU6OzsFQgjftGmT3NDQMGvRokW/amxsfL65uXmp3+9XTCYTr1QqEEWRlctlXiqVUC6Xa1GE6XS6piBDQ0MYGRlBLBZDKpVCLpdDuVzmhw8fJn19fSQWizEAckNDQ6C5ufn/XXHFFVsuuOCCS7q7u1lPT4/S2dkpjPHMSzTF6+7uZoQQvnjx4kBbW9u6adOm/fK+++7bkU6newDQQqEgl0qlmuJns1kkEomaXPUKks1mUSqVeDweZ/v27UMqleJGo1FpbGxcOmXKlOcvueSSX5155pmzNm3aJOP4rZJk48aNSjAYtFJKryCEcEVRBMYYjh5PhBAKAG63293c3HyWOuGO2XMeyz0IAQBCyEnlcpnn83mey+WQSqUQDAYhCEc8R0IpJcVi0fu9731P/u53v4uenh6KUfNnLKBdXV3o7u5WABhOOumkVQ6H47t+v99ptVqZKIqQZZlkMhmWTqeFUqkExhjy+TwEQQAhBISQmjdvoVCAJEkghBwRRFUsFsnu3bthMBiEVCqFQCDAfT4fdzqdbMqUKfMcDsdvLrvssqsSicQdPT09/wWM6f5EC5Mly5cvv83hcLQ+8cQTF4fD4TtGRkYunzZtmqtarT7IGLusVCpRg8EASmkt3agsyyiVSpAkCZTSo2NfSH9/P6lWq8hms0ogECA+n4+5XC64XK4rzWbzF0Oh0I8ff/zxO3t6eiqa4o/VgUtnZycFwEZGRmSLxeItFAq1dLD5fB4Gw/tyS3BZlnmhUDgJwOMbN24cs8l2zDfphJBqMplELBZDNBrF4OAgHA4HwuHwEUoSj8d5Op02PP744+cSQtYBUNauXSt0dnay/8VpRO1WvLu7G3PmzLnE4/Hc5PV659ntdphMJllRFCGfzyupVEqMRqOIRqP/VSwWf2U0Gn9RKpXM6iAiAGobWkmSIIpiTWkYY3KlUhHT6fSLg4OD3xEE4RupVOqKUCgkBoNBFgwGidfrZTabDXa7/RybzXbOJZdc8mA8Hr+7p6dnBwB8xJMh0tXVRVavXs2bm5ud4XD4gj/96U+PhEIhoyAIK0855ZR7nnrqqVsA3A4g63Q637BarYPFYrGBUlq7XdOUoVQqaadx4JxDURQmyzLN5XKpkZGRr2ez2b/NZrMrwuEwwuGwHAgEBK/XKweDQafb7b71K1/5ysXxePz27u7uNcCYKD9Zu3YtXblypUIIwZe//OVzX331VcPw8DCPxWLQ0qYajUZYLKOFsQghKJfLGBkZwcjISPUj/r8f3KAx/C4Bo7bg6Uaj8eW2tjbllFNOEVpbWzF9+nRMnToVLpcLlFIUi0UcOHAAXq8XNpsNu3fvXr93795V11133Q4AWLt2rbBy5cpjPuiuri6q5cNdtmwZent7eVtbG9+4cSNVl320tLTM83q9d7rd7vPdbjfMZrNMCBHK5TJLp9OCajbticfjd6TT6V8BQCgU+mtK6a8551WLxUKtVis1mUzEYDBAkiRtZeEAlGq1KsqyrHDOL3v55ZfXEkLAOT9dFMXVLpfrnFAohFAopAQCAer1epnFYqGccxKPxyupVOr+wcHB21599dUIMKoogUCARyIRomWZ37hxIw0EArxusGkrYm2G/sxnPnNdY2Pjz1Op1IqtW7e+BuD6TCbzH0NDQwcBYObMmZf6/f5/cDgccw8dOiREIhFitVqhXnASo9FYk0s9CmbqhCAUCoXvvfnmmz9W++FKs9n8PZ/PNyscDtfkcrlciiRJYiaTQSKRWB+LxVZt2LChXvlZZ2cnUeXCxo0bsWzZMvYBq0xtYgOAc845Z96sWbPuXLx48fm9vb3o6enBihUrMGvWLLS0tGDatGkIBoMwmUyoVqtIp9NKf3+/sGXLljN++MMfvqK2e0yskbG2++miRYuEQ4cO/cnpdC6cP39+ddasWdLUqVPh9/vhcDggiiKvVCpKOp2m7e3t1OfzMUopfffddysHDx68f9euXbd1d3dHCCF47LHH6gWlfPSw+wNnp2AwGAiFQj9wuVzX+P1+g8ViYYIgQDUVkEwm6dDQUCWVSt0ZjUb/PwCpOXPmeF0u16Wvvvrqv4bD4S8bjcZfV6tViKIIk8mkGI1Grg4kQggRACCRSBxYvHhxKRgMTu3t7Z373HPPHaprxoVGo/FHbre7XR1Qst/vFzweDzOZTIIsy4jFYpFMJvOz3/3ud3cC+NBopfqj1TPPPLPF5/P9wuv13nzKKae8vW7durcPHz784M6dO/8JaraVmTNnrnC73d93u90r1MkBgiDgzTffVA4fPkyNRiMxGo2K0WjkBoMBgiAQSqkAAJlMBu3t7WhqanrroYce6ojH41m1GS5CyHdtNtv3gsGgIRwOs1AoBL/fD5vNBkIIjcfjlWQyeX9fX99t27Zti3yQPJxzsnr1aqIpSmdnp/D4448rnHMsWLAgsHDhwh9Mnz79mrlz5xp8Ph8zGo101apVrFQqsaVLlwpNTU0kFArB4/HAZDIBQLVYLEpbt259s6ur67Suri5lLO/WxlxBALApU6YskGX5BafT6ZkxY4Y8bdo0+P1+2O12CIIgJpNJRKNRtLW1KS0tLbShoYFZrVZBURTs27cvMjg4eNvXv/71+wFUtBVDWx3OOuusuTabbanBYJimmgc7ksnkmwcPHuzw+Xy3+3y+BtXHSFEUheZyOZZOp4WRkREkEonHIpHIj6rV6la1vdKSJUva2tratubz+Ud7enounzFjxjfz+fznCCFnU0oFQRBq+xLGWJZz/ujAwMAvW1paftfa2vrf69ev/+q3vvWtn6dSqZEdO3b885YtW6oYNV2/YjKZbvX5fA3hcBjBYFDx+/3U7XYzg8EglMtlJJPJndls9u5oNBp1OBzLDQbDfFEUIYriNkLIC2+99dYfent7K6eddtoZlNJYMBhMWq3WkeHh4T3PP/98OwAzgCwANDc3n+Tz+W5wOBxf9fl8MJvNTBAEVCoV7USO7t27F8ViMS9JkpVSqq0eUKtnbUin0w82Njb+6+zZs3/0+9///q5vfetbvxwcHPz9mjVrnlGf1zxBEG5yuVyXaKtJKBSiXq+Xmc1mQZZlRKPRwVKpdFN/f/8ml8u10GQyzTOZTCCEHACw+bHHHnsHGF1l6lYUw0UXXXRNY2PjD2bPnh0IBAKw2WyKLMuUc85efPFF4Ve/+hVaWlowY8YMuaGhAW63W0vCLR4+fDjR399/9rp167aN9eXz8TimowDY1KlTF8uy/KAkSe2am4bZbIYsy9VYLPbYnj175rW1tZ28cOFCzJ07V2lqaqKBQIAZjUYhl8uhv79/x/79+2v27cKFC09zOBy32O32s51Op2g2m2sb51KpVEqlUiZJkmCxWBTOOS0Wi0o6nRZV23VnMpn8XjabXa+20bh8+fL/8Hq9Qk9PzxdXrFhx78jIyLRt27Z9CUAJAEKh0FxK6VLG2DQAYIxtk2X5tUQiMaB+xzQA0YaGhtM/9alPbXC5XP9yzz333PClL33pcwMDA89t3ry5aDabw4VC4RqbzXat3+8PhEKhIxRFURRh7969oJTC6XTWbulVz2FUq9XtyWTy1nK5vNpsNrsHBwdbBgcHz69Wq46//uu/fri7u5sFAoHg1KlTv2m321f5fD6DxWLhBoOBVatVmsvlWCaTESKRCOLx+IZcLneHzWbbrSjKEgDzAYBzfqBarW5OJBLvqHK1Atg7d+7c7y5fvvwuSZKuEEXxyd7e3mnPPPPMDtWcvNBoNP6zx+OZ19DQUNufeDweJkmSsH//fiSTyZLH4zFZLJb6y0iZMfZCJpO5ee3atX8CgI6OjktaWlpuam1tndfQ0ACHw6EQQmixWGTxeFzYt28fdu7cuXXnzp07GGMXO51Oyev1wul0glKKQqGwc3Bw8Mo9e/a8fjw8M47XOTbF6GWVdOjQofMZY6eodm4mlUr9oVKp7MJois9VXq/3+vb29sCCBQvQ2tqqNDQ0ULfbzURRFKLRKAYGBp757W9/+xrn/OZgMCiYzWZYLBbFbDZzo9EIURQJAKFSqSj5fJ6m02mSyWSQTqcxPDwcSSQStycSiQcAlFpaWprC4fCKl19++Vdf/vKXb5dl+R/ffPPNM/fs2fOy1vC/YPMsYLRUAgcAURTb7Xb7tYlE4vtXXHHFqU1NTRvj8fhX77vvvofqvqNREIRbbDbbZYFAwBAKheB2u1EqlZgkSdzj8cBisXCLxUJMJhMEQeBqnUAhm80ilUqN9PX1xVKp1Bd37dq1R/1Ow7x5877vcrm+6ff7A1arFQaDQVEU5YjJIRqN7kgmk6vqJocPggAg2k20yWT6nMPhWDYyMvKdq6+++m6bzXb1oUOHTu7p6dmr/f8ArrVYLP+omrYIBoPgnPNyucw8Ho9Q30+SJBHtEi+bzSqxWOyWcrm8ZMaMGRc0NjZCvTej1WqVJZNJYWBgAL29vZE9e/b87N13370TQCUcDs9RFOU8xpiDUkoURXkrHo+vB1A9Xm5Lx/MG+wOTO2t2JwBwzgMAftDQ0HBNe3u7Yc6cOay5uRl+vx9erxeDg4P0V7/6FVwuFxwOh2Kz2QSHwwG73V5zQ1cUhReLRVIoFHg0GsWWLVvyw8PDPalUanWpVDqotWfFihUXt7a2PhqPx6/u6en55YwZM1ZGo9EN6XQ6CeBod3yKI8/TufqjyUTqXgcAhEKhe6ZPn37FJZdc0phKpU5KJBLf7u/v/9ZTTz01qH5kHqX09kAgcE5jY6PJ4/EQl8sFu92Oehd0g8EAxhgKhQLL5/OsUCiIBw8e/Mbjjz9+LwA6Z86cq7xe77c9Hk/tdI4xJpRKJZbJZLRDiIF0On1zIpH4D6imand3N1HbT3Fk39dnjDyi36xWa2DKlCm7Wltb/2v9+vWd3/rWt25PJBKHHn744fvUjzQQQm71eDwrp02bZnW73XC73cRut3OHw0HsdvsRMpVKJaVYLArZbBZ2ux3BYJBp9Uqy2SxGRkbo7t27K/v3779/586dtwGIqKuWgGNcAxBCcPPNNx83n77j7eKh3XgSAOjo6NDKgmnu4QIhRFYvf+ZJknTn1KlTz587dy5mzJihNDc3Cy+99JJy+PBh3tDQINhsNnK0Z63BYNBKpyGVSrFsNksHBgb6HnrooVkA0NTUFG5ra/u91+t99je/+c0tZ5999nPDw8M7du7ceeMYyajJpw28FgB7zzzzzO3z5s1rb2xsDPT29vpTqRR99tlndwQCgelOp/PNUCjkdDqdcDgcxOVywev1wuv1wuVywWw212RKJBJcvdHPb9269QZFUS7x+/0rnE4nTCaTAoCWy2WWyWQE9XY8nUgkfhqPx+8FoG2Wjzm4/gxCR0cHUVfTqQBy06dPP3XFihXPmUym/+9nP/vZdz/zmc/MXb9+/Tucc0yfPn3PlClTZjqdTuZwOGi9p7AagAbOOfL5PJLJJI/H4wqllAQCAaFQKCjxeFzYv38/+vv71+/du3dVPp/fQQjBWWedVb+i1yYtdSxpE9Zx83mbCDEMBKrrsqoo55vN5rumT5/e7vV6EYvFEAqFajOsx+NBIBBAMBiEappAlmVkMhlEo1FEIhGWz+fpu+++e1MkEtn26quvPn3ppZf+dyaTOfvFF1/0ZjKZJNQHejycCtXZTjKbzd/0eDzxgYGBNddcc80bgiC477333mnBYPCuUCh0o9vtlh0Oh1gvUygUgs/ng8ViqVWzUmVCKpXCwYMHoSbUVgghpFKpaBMDHR4eriSTyUcikcgdAHarzREwBgNIu1wUBGGpz+e70WazfeWcc875rN1u/w9Zlk974IEHWqZPn/642+1WHA6H4HK5joij8fl8sFqt4Jwjm83W7shisRgqlQqi0Sj6+/t3Dg0NfW94eHg9ML4On/VMhBgGjlFXEwoAhJD1xWLx+d7e3m8GAoE7GhoahPrYEe30RXNBN5lMtYsvdZNLs9kst9lstzscDkydOnXBq6++eiWAszKZTF6NbCR47yZ6bIUZbWu1WCz+ZGBgAFardd4rr7wypaWlZZUkScxut39BURTOOa+54RBCanKJolgzR7S7ClEUAYC7XC4GANVqleTzeZ5KpbQN+PpYLLaqWq3uUJuhOWWOiXyaG5Asy5uHh4cvBkAVRblh5syZqY0bN77l8Xh+qH1U+xstXKC+n9TvgLp3hNFoZMPDw8rmzZu/PzAwcA/e20vgRHoNfxjjHrRUB8Ooo5kIoHreeefdK0lSjjFGGWNcS1OjuaNrCd00N23tfXXQEUppeWBg4MC2bdsW7d+//8D+/fsfBlDWotGOsxwAIHR2dgr5fH7H9u3b5z755JO/mjlz5lcVRZnOGGOaa71WllnLv1Uul49wMKyTiyiKImSzWQwPD9P+/n6hr69vc39//7lDQ0OfqVarO1S/KApAHmMZtWemeUTz/v7+zueee+6MYDDYYDAYzlQUpeYvpfVRvXv9MfqJc86p1WrNrVq16l5CSBWAqMWgjGHb/1dMhBXkaDgA8uabb54pSZKrWq0yWZapLMtHeKRms1lIklTLUK55oZbLZa4oCqlUKqWXX355Wblc3o/3vIVP5KykaDl8CSERADh8+PAOi8WSNRgMNlmWebVaJUd7D4uiWEsqXSgUkMvlUCqVoJVi2L9/vxCLxXYkEomfpNPph3FiZ12mlZYGMKgGWk0nhFi1wa8pg+ZWn81mYTAYIMsyCCEoFAooFAool8ukUqkwxphr3bp1Z3LO/7uzs5NPtLzHE1VBOKV0L2OsUK1WLdVqlZfLZVLvrk0phSzLtaW7VCohm81qDx/VarX66U9/Ov7UU09ppzLjZsuqZhdva2vbeeDAgVKlUrFXKhVeH5+hOQ1qMmk+RtlstlbDo1wul/fv339HNBr9JwAV7XRHdco80VDOOVwuVyyTyeSr1apFk6l+EtMCzSwWyxEyFQoFrU8LiqLsxajJO96RmO9joioICCFRznmmUqlYS6USlySpZo9TSqEVgakPGVU7hpXLZVosFneuX78+r/owjfeDZwDon/70p3IwGOytVqtnFYtFJkmSoMmjeQrXy6R53BaLRaVSqdBKpfKnaDS6Wg0vFTnnY7bP+Kgy7d69O6/JVCqVmMFgEOr7ST3arSm95iRZKBRQKpVIPp/P7Nq1K6p+53j30/uYqAoiDA0NFcLh8HpZlq/K5/OKKIqiNpC05MvFYvEIL1v14XN1E/ssALZx40YRE8OmpQBkSumznPOOQqHAtVUDwBEyaa9rNrw6M5NsNvsHAMKpp55KAYy55+pHoCYTY6wmk+aao3oI19zq6/upUqko1WpVzOfz64eGhgoTNcXQRDjmPRYUAAuHw3MA7CCEwGazUZvNRsxmM8xm8xHJE7TOkGWZKYpCUqlULJPJtG/dujWmft9EmJkIADQ1Nbmq1eobhJDpJpOJ2e12wWKxwGQyHXFqpbmgy7KsMMZoKpXqT6fTp27bti11onJC/QUQAAiFQj5CyE5CiM9sNnO73U4tFsv7+qnOrZ7Lssyy2Swymcy8LVu27MIErRo8kU6x6mEYXUV2AfguIUTI5XIsl8vVIuLUh4t0Oo1MJoNMJqNks1mSSCRIMpm8buvWrVE18GYiDCRAPXw4ePBgklJ6mZqsTUin03Imk+GaTHVy8UwmI2ezWe0o97Lt27cnV69eTTCxZKLDw8NRzvl1AEixWCSZTEbJ5XLI5XI1ebQfNc6fpdNpIZFIfHfLli271NO3CaccwMRdQTQEAEogEPglpfRvGGPMaDRSu90Og8HAJUlSqArnHMViUS4UCtf19fVN5JSg2up4LoDfEEI8Wjqho2I0oMqUyOfzl+7fv/85TNBZFu/Vl/8agHsppaJ6/8GMRiMzGAyCKIqEMYZKpcIopbRQKPzbrl27/nYC9xOAibkHeR+RSOQX4XD4bxoaGpDL5RCLxaDGNYjqxi/HGPuvcrn8zyMjIxM9qbS2Oj7X0NDwCc7531er1c/LsuyqVCpEtdW5LMspxtgThULhR7FYbA8+mrvIiULBqEy/CAaD2zjn/1AqlZaXy2WbwWCgoiiCcw6LxYJwOAyLxYJ4PP6L8W70X8JENbGOxkYIQWNjI1asWMGXLl0Kg8FwIJPJ3JVIJP7Y1tZm+/SnP90/MjKyuaOjQ5zAyqGhABAGBwf3DA0NXZVIJL5iMpnI/PnzleXLlytLliwh4XD4KwcPHrxqEiiHhtLR0SGOjIxsdrvd/eeee65typQpf8zn83elUqkDlFLMnj2bz5s3DzabDfl83jbeDf5LmCwKorlYgBCiNDU1obW19d8ikciqUqn0ry6X60VCiBEA1cJWJwEKRlMgCZVKJUMIqaUFNRgMsNlsGQDCokWLJEx85QAAqM+ecs6Nfr//xTlz5vxrLBZblc1m/03dqCtaVStK6aTop0lhYmloTnNqQU0RgFgoFNY8+uija7TPdHd3j1fzPgoMo4N/1NmK85q7jCqf0tLSgi1btoxrI/9SNJfzXbt2Xbdr1y7tZZFzLnLOofrBjV8DPwKTZQU5AvVYlwKQHQ7H31588cUvXX311XcDtUTHk433nUpNxFIRfw7t2c+ePfvu66677qUvfOELf4vRe5LJ2CcAJqmCAO95mDLGmgghZ3DO5413m3RqzCOEnEEIaVJ/n1zLRh2TysSqh1IqA0Amk+kihNysva6XDhs/tGe/e/fuc3fv3g2oisEYm7SF2CetgmgxIi6X61MXX3zxJx0Ox75f/OIXDx+vNPg6fx7t2c+dO/eKZcuWzRgcHHzxySef/C/dxBoHNAVhjH1KEIQuAFcBgHrTrDMOaM+ec36VKIpdgiB8Sn1r0o6zSbuCaCZWLpe76dFHH71Je103scaPulOss+tOsXQTazzQVhCHw7H8vPPO002sCYBuYk0g6kyssymlXYSQKwHdxBpP6kysK0VR7BJF8Wz1rUk7zibtCqKZWOedd17Xb37zmx9or+sm1vhRZ2J9SjWxBEA3scYLAgDPPfdc08qVKxscDkfql7/85c7xbpQOMHfu3PazzjrLNTg4OLhu3bp+TOJ7kEm79HHOJTXDxg0Wi+UlSulDwKS9SX/fAJpAteT/YrRnzxh7yOFwvCRJ0g0dHR2iKIrSZHMx0ZiMg0kL3cxv2rRJTqfTXdlsdiEh5JJxLCv9UdHqF2rZJeuROzs7BbfbPWn6SH32lFJ6iaIoCxVF6dq0aZNcLBbzsjw5razJYmLVkqwxxoRCoQCbzfY3F198sY9SOmg2m1+KRqPp+++/n03wWJB6BADVnp4eWCwWh1Z9S3NWpJQ6VDmUySKT1s58Pq+8/fbbJkrp1770pS81xGKxz5XLZS2/rlbDfVIo/mRRkJwWdBOPx4koigiHw9NNJtMN2tJNKY1ccskln12zZs1ED5iqDaQ5c+bMCgQCf2+32z9vMpm4xWIRKpUKjEYjnz59+kNXXXXVE4lE4kc9PT17JotMCxYsWOpwOJ5yu90BLX9AOBzWklIQRVE4IQS5XC433m3+S5jQWtzZ2QkAmDNnztdaWlpqYalOpxNWq5VZLBbZbrfLVqtVcbvdAb/f/9IVV1zxNa2q7Dg3/5h0dXVRVTnO9Xq9r3q93iudTqdby/SuJucmXq/XHQgErgyFQq9efvnl52olpce7/cdCU462travORyOl7xeb8BgMCgmk0k2m82yxWJhLperlqPX5/Ph1FNP/dp4t/svYSLvnASMxkPcEAqFfuJyuRS73V4rfaAqiVb+AMVikeXzeVIoFEg0Gu185JFHHp+Asy4FwJqbm5d6PJ5XA4EAsdlssioXqS/pwBjjpVJJKZVKYjqd5vF4/BNr1qzZPNH2WdozbmpquigYDPZ4PB6ulj6gWsJxq9UKo9FYy5GVz+eVSqUiDA0Nfefhhx/+6QTspxoTVUEoAOb1eucEAoEdfr8fTqeTOp3OWvkDrVSAyWSCLMtaqQCWyWRILBaLHTp0qP0Pf/jDhEr7o2YidAWDwTeCweB0h8PBnE6n4HQ64Xa7a6UCTCYTFEXRSgUo+XyeRqPR/sHBwVOffvrpCZf2x2az+cLh8M5gMOhzOBzc6XRSrZ98Pl+tpINWajuZTPJ0Os0SiQRGRkbm/e53v9s10RRfY0Iu2VAfvMFg+J7BYBAEQeCUUkIphSRJMJvNsNvtOFpZ7HY7lSRJcTgcfo/HcyUA3tXVNVFMLYEQwi0Wy9eMRmMLIUQRRoHBYIDFYoHdbq8pis/ng8fjgcPhEERRVJxOZ4vT6bx65cqV9Oqrr54oe0cBALdarVcajUY/IUSho4UdtSKosFqtqK8V4na7YbfbidFo5DabTbBYLN8DgPb29gk5WU+UB10PwWgYqoVSej4AcM6Fo0sFiKIISZKOVSqAEEI4pfQzAH6MiZMmh2G0DspnCCGjTkt1MmmKYjKZavVBtDy9kiQRQggXBOG8Rx999A7OOdu9e/dHqbN+XGSilNZk0t6o7yeDwQCDwQBgNFxarcMoYPSPz1+0aJFl5cqVBYz2/URYGWtMxBWEAIDL5fITQhyMMc4Yq5U+0NLpa+UBtB9ZlqF+jsqyTBhjbbNnz7aq3znecmoJ7IwA2hhjRC3rcERJB60UwjHKOQiyLEMQhNM+97nPrSaESGrVXz6OhxGkq6sLXq/XSghpVxSFaCUd6mXREnRXKpX3lUBQFIVTSh3BYNCvfec4yfKBjPfAORYEALFYLK2cc4tWJkBTivoM78lkEolEAslkEul0WsuCjmq1yjnn1lmzZvkmSL0JBoDPmzevnRBikmWZ1+XdRblcRrFYRC6XQzqdRiqVQjqdrmWrVz9DRFE0zZgxo+srX/nKls9+9rNXAZC0060TecKl1SHp7u5muVzOyxiT6uucaDJpWTBTqRQSiQQSiUR9P5FKpcIZYxbGWCsA0tnZOeEUZMKZWB0dHeSPf/wjHxwcfCkQCKSq1apbS958dPmDarV6RPmDXC6nZQ0H5zzf1NRU/pu/+ZvWAwcOFJ5//vkhdT/CTsRmUJvZe3p6lGAwGBgZGYlOmTJlXrFYtJdKJaVSqQj1ypHJZGrln3O5HAghKBaLR5R0qFQqoJQqU6ZMmedyuf7tsssuuzqbzf6wu7t7g/Z/trW18eMtn3bi9MlPfvKrgUBg9RtvvGFSi/6Q+klMK2kty/IR5Q/UtKS8UqnQXC6X/MMf/vASRstDTDgFmTAN0mZAtXOlCy+88Jt9fX13ZDIZwev1Uq1aqtVqhfbv+hrcWscUi0W5Wq2K0Wj09+vXr//i9ddfP6Ioyt577733k/iQCrVjyBHpQT/xiU/Ma2pq2mQymVY99NBD/9bU1PSuoijNLpeLOxwOqslUX7nXaDQeXdIB+XyeJ5NJZjQaEQwGid1u55RSIZ1OI5lMro/H46s2bNiwA/iLSln/j9F8w2bOnGkIhUIrHQ7HHwOBgKNSqWx7/fXXkc1mucfjqVW1re8ri8UCo9FYq3+i9hMymQwzGo3K7Nmzv3/TTTfdA6DKOaerV6+eMF7ZE2EFIWvXrqUrV65UAOCcc8453+/339XS0tI+ZcoUPP744ygUCrXs4MB7pQK0wjNaWn3VziWlUgmKojwiSdIpr7zyimv+/Pk/5ZxLN95440NDQ0N3r1mzZvNxlIcBkM4888xvejye+CmnnLJmx44dA3v37l3NOX+wUqn8XhCEGwuFgiJJEtVKOgCoVZHSZNIqNqml2cjevXsFzjnS6bTi9/uJ1+tlagnp82022/LOzs5HBgYG7ti0adNuYEwKYZLOzk7a09OjqEFopK2tzTp16tSHFEV54de//vWnWltb/yqfzy+glN5eKBTYaJWK9/pJvaOqTWTa4UOlUkEmk6HnnXcebW5uvvvnP//532Sz2e8RQtYDwNq1a4XOzs5aPfrx4riXgV67di31+/1k2bJl2LhxI5YtW8bUYvWkq6tL6O7ulgGgo6NjXnNz853hcPj8QCAANTG1sGbNGqW3t5d7PB5BK39gMpm0050jyh8oisLK5TJNp9N9r7zyyiy1DXMBvLNkyZIfn3766X9HKT23WCy+MTg4aHvyyScPd3V1Cb29vf+bgp60s7OTqAORzp07t+Wdd97Zu3Llyu2SJLXv3LkzsHfvXr8kSTSZTO5wOp3TTSbTm5RSpzrbEq38gZZV8aiSDly958lv2bLlhkKhcInH41kRCoUQDAaVQCBAvV4vM5lMAucckUgknc1mf9rX13fvtm3bIsB7l3n/m45ctmzZ3X6/f95vf/vbc88999ybo9HoZw4cOPCFWCw2BADBYHAPIWSm1WplNpuN1pd0qC9/AACKovBSqaQ4HA5y9tlnC6FQSHE4HEI+n0csFls/ODi46qabbtoBAC+88IJ49tlnKxg9sqfLli2jansAtW798VSi46YgH3bxwzkXKKUK5xwLFiwILFy48AfhcPiaxsZGg9lsZlq54GQyiV27dtGNGzeCUgqTyaRYrVbBbDbDYDDUZ0LnnHNSLBa51WrF1KlT89lsdu3hw4d/+MorrwwSQmC1Wr9us9kah4eHb7r22mt7JEla/s4778zZsGFD5Fht7urqohs3bqxtfAOBAD/Kvn+fmbZixYp7AoHAFUuXLm188sknT0okEt+Ox+PfOnTo0CAAtLS0zPN6vbebzeZz3n33XVO1WiVOpxNH1wfRCtAAYIqisEqlIsZisW/s3LnzXs45JYRcZTKZvu31eueFQiGEQiHZ7/cLbrebGY1GoVqtIhKJDGSz2ZufeOKJ/wBQ6erqor29vaSnp4d1dnbSSCRS6/tly5YdsS/TzKnGxkb39OnTz/njH/+49itf+crfAnhg+/btl23duvUxqGZkOBxuCgaDq91ud+fhw4etyWQSdrudWK1WbjabidZPdfVBFMaYkMvlIIoi5s6dy2bOnInm5mYEg0HY7XYaj8cr8Xj8/l27dt3W3d0dAUZXFM3KOMZ4ouqkO+YcFwWpG2jSHXfccf6MGTNOmTp1KhdFMROLxf5w3nnn7QJguOyyy1ZNnz79+unTpwdsNhsEQVDy+TyNx+Ps8OHDQl9fH/bt2/dMMpl8jRByMyFEUC+gFKPRqFUzIpRSoVqtKjabjc6fP584nU5Uq1XEYrFIMpn8pyeeeOJfAVQAwGg0tra2tm6dOXPmA08++eSN119//U+SyeTGhx9++Mm6tn9gJdzOzk5h7dq1tVnrpJNOag+Hw9cuX778+3v27Dk1kUhs3Ldv31e3bt36kPYdDQ0NjYFA4Ba3232Z1+s1OBwOpNNpvPnmm6xYLHJ1VeRGo5GoMnEyipDJZDB16tSRefPmxQ4dOvTFxx57bI/aFAOA71sslm/6/f6AtqL4/X7qcrkUSZLEYrGIWCy2I51Or3r66afXf1ifcc7JsmXLBPX4GACwbNmyM1pbW1/KZDL/tHbt2pumT59+VSqVej6ZTB4EYJo3b97VLpfrJp/PF3A4HKCUYtu2bfzw4cPMZDIJRqNRMRqNXN2sE0KIoNZFUdLp9C2Dg4NLpk6dekF7ezva29uVWbNm0ebmZhYIBASTyYTh4eHI8PDwz1avXn1nb29v5cYbb5zT3t5+3vTp0x02m40AeGvJkiXroe5djoeSjLmCaMpx4YUXLj7ttNMeXLx4cXtbWxvC4TAEQUA2m63u2bPnsccee2ze/PnzT3a5XJAkSSkWizSRSLDBwUGhr68Pu3bt2tHf3397LBZbAwBTpkw5jXN+C2PsbEEQREEQanU0FEUplctlU0NDA0KhkOL3+6nH41EMBoNYKpUQjUZ3JBKJf1i/fv3TajPnATgwd+7cz33qU596mFL6vXvuuefHF154Yeu6dev2AsBZZ50112azLTUYDNNEUYQoituy2exrzzzzzAAAzJ07d9o777wTveCCC05vamracOjQoX95+umnb/B6vZ+Lx+PPASj6fL7wlClTrnE6ndf6fL6A6melyLJMq9Uqi0Qiwp49eyDLMkRRrNVJFwRBy9O7PZ1O3xoKhVbPnTvXfe2117Y8/vjj56dSKUdra+vD6iQUJIR802azrQoGg4ZQKMQDgQDz+/3U4XAwQRCETCaDVCq1IZ/P3zE0NLTb7XYvMRqN89WBewDA5scee+wdAJg2bZqptbX1t16v90BPT891y5cv/10kElEKhcKX9+7dWwaAWbNmne92u+/y+XztqmuMzBgTyuUyy2Qywp49e5BOp0uSJJnq+4kxJjPGXigWizdHIpE/qX1xicViuWnGjBnzTjrpJMyfP1+ZOXMmnT59OvP5fIIgCNi/f//WJ554YsesWbMuXrp0qdTc3AybzQbOOaLR6M6BgYErFy5c+PrxUJIxVZCuri56yy23sPb29gUdHR0vLF++3DN79my5paWldhwL9WDgtddeQ6FQ0MqLsYGBAWHfvn3Ytm1bpL+//7b9+/ffD9U06O7u1mp/IxQKzaWULmWMTQMARVF2VCqVN9PpdIfFYrnd7/c3hMNhhEIhJRAIUKfTyQRB0E571gwNDd2+adMmrazbBQ6H4wsDAwN/e/XVV3/P4XD8/datW6+tVCpXOp3Oc5xOp2A2myFJEgCgWq1mFUV5dPv27b+cMmXK7xwOx38/9NBDX128ePHP0+n0iN1u/+ctW7ZUAYjz5s37it1uvzUQCDTYbDYYjUaFMUaLxSJLp9NCPB5HLBbbmc1m7xYEIUoIWQ5gvvqMtjHGXohEIn/A6Mp3BoCY1WpNfuELXxhxuVx7tm/f3r5jxw5zPB7Pqn9zkiiKN9jt9q8Gg0GEQiEWDAbh9Xpht9tBCKF79uxBPp/Pe71eq2amqsfliizLL4yMjDz5zDPP3HP55Zc/UqlULnvllVdOHhgY6IVaD3HatGkn+3y+v3e73Rd7PB5YLBaFEEI1xUgkEhgeHh7M5/M3AdhECFmI0ckInPMD1Wp1cyKReKduHDBCCOOcGwBc4/F4ftDe3h5YuHAhFixYoLS2ttK5c+eywcFB4b//+7+xePFizJ49W/b7tXvF0e8pFAqJffv2nX3SSSdtY4yNqZKMuYJ0d3cLp59++p8uvPDChSeffHJ1/vz50tSpU2ufYYxxSqmSz+fpyy+/TKPRKNu3bx996623Kn19fffv3LnzNgARrcQx3kv9r91Gf9iGLEAI+YHdbr8mFAoZwuEwC4VCNTdrADQajVYymcxde/fu/fHbb7+dAgCv12v3er2bmpubTymXy/B6vbBYLLBYLIrZbK6ZCFCTECSTyQPpdLp06NChqb29vXOTyeQhrQFz5sy50Ol0/sjn87U7HA6YTCYZgFAqlWqDKBKJRJLJ5M8SicSdUE2/D6H2DIxGY0soFPqFz+e7+fvf//7br7/++tvJZPLBX/7yl/+E946WV0iS9H2Xy7VCVRQEg0EUCgWlUChQj8dDNLlMJhNEUSTqc0a5XEYqlfrqr3/96/+YMWPGtcVi8TeDg4Nxp9PpmjZt2t85nc5Vfr/fYLVamVajPpfLIZlM0uHh4Uoikbg/EoncBiBybFEAqBfBde2t7+MAIeQHDQ0N1yxatMiwZMkS1t7eTt98803W2NjIFi9eLMydO5cYjcb676sCkA4fPvzmrbfeetr999+vTFQFEQAodrv99CVLlry8YsUKZcGCBcLChQsRCoWO+KBaHwJPPfUUenp6sHfv3vXbt29flc/ndxBCcNFFF33Y8STFkR4AmtLUVhkA80RRvNPlcp3f0NCAcDgsBwIBQTvtkWUZ0Wh0Tzwev+vpp59+EIDi9/u/29jYeFcwGFTUOwlqt9uJw+GAxWLR6n3zUqmkFItFMRaLKa+++uplO3bsWAsAs2fPPt3pdK52uVzneDwemM1mBQCtVqssk8nQZDJJRkZGKolE4v5YLFY/iES1/fWDRpsMWJ1sqHsf06ZNu+6kk076eWtr64ovfvGLrz366KPXHzx48D+efvrpg+pHLjUajf/g9XrnTpkyRbDb7UR1FOQOh4M4HA5YrVYYDAZwznmpVGLFYpEVCgUpFot99eGHH/53Va6rPB7P9/1+/yybzQaTyaQoikILhYKSSqXEaDSKWCy2fmRkZFW1Wt1RJxPDe8qAuvYfa/ASjPqpKaprzTxJku6cPXv2+c3NzQiFQjjrrLOwYMECtLW11VZ0AFpJBaVUKgmvvPLKGcuXL3/lwzb0/1PG8h6EAIDRaFyhKAovl8u8WCyiUCgc67McAKZOnVp8/fXXL9m9e/c6Qkjt3P7PHEl+0EPWOoQC2CHL8mdisdgluVzupng8Pk9VFO73+7nb7VbC4fAsp9P5i8suu+ziffv2/TKdTt/kdDohiiKRJImaTCbYbDa43W44nU4YDAbIskxyuZyYSqVkQoh4xhlnfFOW5X1Wq/UbTqfzCp/PJ5pMJiaKIiqVCsnlcjyVSgmRSATJZPLBoaGhuwHUDyIF7yn10bIc63fS1dVFVq9ezQkhjx46dCjNGPuv3t7eHy9cuPDvKKUvhsPhmMPhEHbv3v2bcrm8nlK6U5KkBpPJxCVJIkajkWgetpobOuecFAoFIZlMUs45t1qt95522mnVSqVyVSAQWK5+TgYgFItFnk6nSTweFyORyI54PH57JpNZAxxx7/I/DUDnABT19IwSQnZUq9XP7Nix46/y+fyacDhsLhQKKBQKpFKpQJKko2uNcIPBwDnnKwC84vf7x2ziH/OLQs65lMlkaj5Sw8PD8Hg8cLlctc9QSlEul0m1Wq3s2rXrOQBk5cqVdAyCZjhGB502464plUq/GxgYWJVOp7+bSCScqtlF/X4/s9lsvKWlZUUikViRTqc1r1pa74LucDjgdrthNpshyzIkSYIsy6J6offJpqamN7QbY0mSlGq1StPpNEun00I0GkU8Ht8QjUbvKJfL/6W2Savo+lGyGPDu7m6uFglKMcYe4ZwTn89X7u/vX7t79+4Xb7jhhh/Jsnz5aaedNn/dunWnOp3OBkEQmOaGrlWx0sIFbLbRSmi5XA6MMVIul2E2my3hcPgRQgjsdrtCKSXVapVmMhmeTCbFkZGRdCKR+HEsFrsT7+0TMVb9p5p87Lzzzntu586dlUQiYUkmkzwajWJkZATNzc21OxVNSdTxJn3wV380xlJBtGD8txOJBDl8+DDx+Xyw2WyQJAlNTU3QjgILhQJPJBI4cOBA/Ic//KH43HPPVdSZZ6yot28rAG7L5XJri8Xi3yeTySsTiQTC4bASCATolClTlGKxSCilBGpQEzD64DWXeqPRWFOQcrlcu68ghPBQKEREUVQYYySbzfJ0Ok1isZgQi8V2JJPJo2dXTYHHgppZAuCmeDwOAKGXXnrp+62trc+vXbs2FQ6Hr9TuU7RMlJprvaYoZrMZAKAoSu1Sj1LKXS4XB8AVRakp/MjICOLx+INDQ0M/AqAdNwvd3d1jHQ3IAOCtt94SAcQHBwedQ0ND3Ol0EnXFQyAQ0FZ1pFIpMjAwQAYHB98GgGg0OmYXh2OpIAwAoZT+sVgsJvft2+dQjxlptVpFMpmEasKgVCop2WxW6uvre3jDhg351atXi/hoM+qfQ8HoQBIA7FEU5apEIvFILpe7PZFILE0kEshms4hGo1Rz7Tjatb7eS7X+fVWRCKWUq2aHkEgk6PDwcCSVSt2WSqXqT+HIcQgprVe20ZtSQiJvvPHGX73xxhvbzznnnHn79+/vrFarjDEmajIdHTZQrVa1o/KaXKqpQwqFgpJOp0k0GhWi0ejmeDx+Uz6f11ZCzUQ8HqGyHIA4MjKSD4VCDw8NDa3et2+fbDKZKOcchUIBXq9Xm7RYqVTC9u3bk6+//vof1bZPyE06oG7UGxoavkEIuScYDCptbW28sbGR+v1+YrPZuMFgUCqVitTb2xt5+eWX215//fWE+rfH2+dG29zLGC2e+XWLxfIDn88X4JzXovm0eHctXNTv98PtdteCmNTQXsTjcaTTaUSjUQwMDCASiaTT6fRPIpHIfQBG6p/HcZbrfRBC4HK5LjCZTL+32WyCx+OhDoejZi56vV74/X54vV7YbDZtVUcymUQ8Hq+FEfT392N4eDgSj8dvS6VS9wGoqo6QJyKEgABAQ0ODhzHWa7PZArNmzarOmDFDCIVCxO12c5PJxBhjJJlMCr29vd988MEHfz6WG/RaI8YYCoCFQqGfEEJucLlcaGhogN/vr13uJJPJPfv27bvi7bfffm0cYpGPOFaUJOlmr9d7jd1upy6Xi9rt9pqSuN1uuFwuOJ1OaEeLmmt6Op1GoVDghw8frmzZsuXRSCRyhNmB906gTjRa7ikWDAZ7LRbLXJfLxRwOB1UdG+FyuWqHDxaLRdsTIpvNIpPJIJvNsmQyyV599dX7Dx06dAveO3E70QpPAbCpU6cukWX5YaPROCsYDCIcDsPtdkOSJORyOQwPD//0hRde+M7xGEvHyxeLAOANDQ0rZFm+nFK62Gg0hiRJ2qsoyh/6+/t/AiCNo1zDTyCa2SVTShEIBA4bDIYGNa6daO7a9S729SGjdZk56MjIyB/+8z//8zPqvY1mdox32KgAQAmFQk+IovhZu92uOJ1O0Wq11uSqd62vd0PP5/O8XC4jmUwOrlu3biqA8ZaLAmAtLS3OfD7/Hc75eQaDodVgMAxTSl/PZDKPRCKR549X2Yvj5e7OAQiDg4PPA3i+q6uL/vSnP3Wk0+kUMGoC3HzzzeOZxYJDNbUYYwqAjYqiXJ7P52Vx1F8bwOh9TbVaRbFYrHnY1rugl0olUi6XHwEgLFy4kKq36BMBAgCCIPyWMfa5QqFAtFvzOo9alEolSJKkuYFocimyLIvZbHYj55wuWrRIGGe5GAD67rvvpgGsBrB62rRprnfffTejvofOzk5BPawYc463u7sWL60ARzjETYRZFlBnp8bGxlNlWf4TIYTb7fb3udVrkXHaQJJlWVEUhaZSqb4DBw6ctH///vIESsUDqP3a2tpqyOVy2yilM81mM7Pb7YLFYoHZbK6dWGkKornWK4qiZLNZkkwmT3vrrbfewPit8kejrfq1saOeDALH0ew7URGFxzOC73+LAEAJh8N3U0q/A0C2WCzUarXSY7mgc84VRVFosVgkmUzmnO3btz8/QROfaQcmKwBsAMBNJhOz2WxCfZyGpvicc8Y5Z9VqVYzH4z/Ztm3bjRNULuAEjqcJE3I7jmi370o4HL6PEHINABgMBm6xWEjd3QD46NkuKRaLcqFQuK6vr+8XE3gQAe8p/9cA3EsIEUVR5GazmWjOimr2Hc45J9VqFfl8/v5du3Z9fQyiET8W6AryHhQAa2houIRz/gMA7QA4IYRQSiEIAlcd+zbk8/kfjoyMTPgk2SoCACUYDC4VBOFWzvk5UH2/tGAzSZII53xnpVK57fDhw2swccyqcUdXkCOhANiiRYuk3bt33zd79uyrTj75ZCUYDJJMJkM3b9589euvv/4L9bPjcsfxEam11el0fi0YDD6wcOFCNm3aNM45F956661fbdiw4esY9YzVlaOOiZgXazxhAIxbtmypFgqFfXa7HTNnzuQnn3wy5s6diyVLluwAQBYtWiRh8igHMNpWiXNO0un0DoPBgObmZixYsIAvWLAA7e3t+wBUr7/+eiN05TgCXUHejwKAEEKMAFDvopHL5SwAeEtLy2QcRIxSyjGa0lWLWNSywRgBEI/HM5mU/oQwEdL+TEQ4AF6fL1etTTIZFeNoGKUUBoMBRqMRnHNIkvTnAtH+z6IryIegZZPXki9rZdImO4SQmlyc85qXgM770RXkQ9DcwrXBJIofj8elyaUpyCQpFzgu6E/mQ9BWkPrkZ8B7peEmK9oKUr866hybj8eUeJyon2k/TiuIVrfDZDKBc/6xket4oK8gH4K2gmhKUpe6aFLzQaajms5Tpw5dQY6BWrqAaGG3mj+WGtE2aZ8ZY4wKgkABHLEySpJE1q5d+/E4gRhjJm1nHy/a2tpoT0+PoihKWUuQLQgCjEYjhoaGsoQQ1tLSMume26JFiyghhCmKkq1P/G0wGFCtVssrV65U/H7/pJNL5wSiFb259NJLP33nnXf2Pfvss+ydd95hkUiEpVIpZXBw8IXXXnvtFGA0YfL4tvYvR6u98tnPfvaU22+//YXf/e53yttvv82Gh4eZmtWyb/v27Z8GRpNEj29rdSYk2iD65Cc/+dc333wz37BhA+/v72eyLHMVxjnn2Ww28eabby4GJoeSaHKdfPLJi2+44YbE008/zXft2sVKpdIRclUqFd7b2/vXwOSQS+fEQgEQn8838/Of/3z+vvvuU1566SU5lUrxo6hwzvnw8PCeZ5991qhl/xjvxn8IRG2jcfny5Xvuvvtu/vzzz1eGh4ePEEpRFJlzrmSz2fzGjRtncs6JriSj6A8BQEdHBwXALRbLN5xOp0UQBCaKonCM+wEJgBwMBmd6vd6LCCH8hRdemLAmSUdHh0AI4V6v9yK73T5TEASZUippCSi0HGCUUgEAs9lsFo/H8w01tlsfG9AfAoDRAjLqP0+RZZmXy2WiFqUE8N5AqoNzzs9Q//aEtfOjYjAYzmCMca1oaLlcPtbHtBrlp6i/fxz8zv7X6AoCQE3liWq1KmcyGaKVLY5Go2CM1eeABUYjCkksFmtVf5+wA0nNXwVZlltVuUgymUQkEkGlUjlaLsiyTPL5/PFI4Ddp0a9Q65BlGUNDQ3x4eLiWRE4URQSDwZrfUjab5fF4nEcikbcAQC3TNlGVhGI07ehb0Wh0xfDwMPd6vXA4HDAYDJgyZYqWvBrFYhGJRIIfOHBgvNs8odAVZBSKUTfw56LR6KfeffddZrFYBEopSqUSIpEIbDYbGGMoFAo4ePAg6evrex4Y2zywxwEOAJIkPZ9KpVb19fVBLXWHSqWCeDwOp9MJACgUCiyZTApqMnGsXr16Iiv+CUNXkFEUAMRoND5QLBav2bNnz3RKaaVSqUjZbBbDw8OwWq2cECIXi0XD22+//dy//Mu/bOjq6qJjmebyOKAAoENDQxuCweBzBw8ePFeSpEq1WhVzuRyJRCJQy8JVy+Wyobe3t3/Hjh0PqEnYJrJcJ4yJfER5oqnlyKpWq0+ZzeZwQ0MDpk6dCrUgjmaCvfb2229fcOmllyY+rNjnBIIAIOFw2MM5f8ZgMCwJBoNobGyE3++v5RweGRkZOnjw4GdfeeWVN8YhHeyERVeQI9GUpKFarXZzzs81mUxTzWZzRZKkd2VZ7nnnnXfuBpDBe6c+kwECgLe2tjqy2eyNADoNBkOLxWIxGAyGw4qiPNfb29sFYFBXDp0/R+1kLxwOW9xud7vRaJwONUukevIzGSeWWps7OzuFYDA43Wq1tgOwaK9rt+46On+OWsHO2gtqiThMTuXQeJ9cAD4Och039Ify4dS7kkyG/cZfysdVLh0dHR0dnY8zXV1dVN/w6egcDQcB1/c1OjofCAcn537m3LPvuuuuZWosgo7OpGVszCB11bDdbvOTO4itwdXwtdlzZ59JKeV64I3OZGasfLEoAF5oL/wGDNlLyaX/6lvq26HGUehHiDr/h+kCBQfBeZiBV7EVf8LbOA9Ta+/p6Pyfh0NQPZPOAcc5XV1dFFz3FNbRGYWr7gscnwXHZ494TUdnEjPWs7y+39D5WKHvEXR0PoTjpSD6/YfOx4KxVZCCQAGBq3mVPhCuK5DOJGFsFcSi5A1QRM65u2s1BI5j/xBC+AtdXSIAg37brvPxZy0EAgA/XvG3wW+veOQf/m71g7BZjv1ZQnAaMGXekiU3/+z++2/nnIvQlURngjIGp1hrBaxcqQhi1xLlu+b7Lc0H6adXWR78p6L8V7LBQMVKpRbfLAMElGIj5//wktu9mEUif08Ikbs6OsRuQE9YpjPhGEMTSwInCs2VGvHaC+nXZaUCpfL+Y19RUTil1Hna0NDvmw8cWM8BsnrTJj3FjM7/BbqWjP58OIuAphPRGh2dCQchAOdd9IM26FxNGsD1oCqd/3t0CqM/H06XfkGpo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6NzHPj/AY309sREO31yAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<PIL.Image.Image image mode=RGBA size=200x200>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from scm.base import ChemicalSystem\n",
    "from scm.plams import view\n",
    "\n",
    "naph_gas = ChemicalSystem(\n",
    "    \"\"\"\n",
    "System\n",
    "Atoms\n",
    "C      -2.43705000       0.70853000       0.00000000\n",
    "C      -1.24580000       1.40435000       0.00000000\n",
    "C      -1.24580000      -1.40435000       0.00000000\n",
    "C       0.00000000       0.71830000       0.00000000\n",
    "C      -2.43705000      -0.70853000       0.00000000\n",
    "C       0.00000000      -0.71830000       0.00000000\n",
    "C       1.24580000       1.40435000       0.00000000\n",
    "C       2.43705000       0.70853000       0.00000000\n",
    "C       2.43705000      -0.70853000       0.00000000\n",
    "C       1.24580000      -1.40435000       0.00000000\n",
    "H      -1.24407000       2.49615000       0.00000000\n",
    "H      -3.38504000       1.24787000       0.00000000\n",
    "H      -3.38504000      -1.24787000       0.00000000\n",
    "H      -1.24407000      -2.49615000       0.00000000\n",
    "H       1.24407000      -2.49615000       0.00000000\n",
    "H       3.38504000      -1.24787000       0.00000000\n",
    "H       3.38504000       1.24787000       0.00000000\n",
    "H       1.24407000       2.49615000       0.00000000\n",
    "End\n",
    "End\n",
    "\"\"\"\n",
    ")\n",
    "view(naph_gas, direction=\"along_pca3\", guess_bonds=True, width=200, height=200)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2a84ff0a",
   "metadata": {},
   "source": [
    "### 2) Shared M3GNet settings (maps to GUI panels)\n",
    "\n",
    "These settings correspond to:\n",
    "\n",
    "- **Task → Geometry Optimization**\n",
    "- **Model → M3GNet-UP-2022**\n",
    "- **Details → Engine Add-ons → D3 Dispersion → Yes**\n",
    "- Optional: **Details → Technical → Device → cuda:0**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "ce54d728",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:06.419453Z",
     "iopub.status.busy": "2026-02-11T10:35:06.419317Z",
     "iopub.status.idle": "2026-02-11T10:35:06.422292Z",
     "shell.execute_reply": "2026-02-11T10:35:06.421894Z"
    }
   },
   "outputs": [],
   "source": [
    "from scm.plams import Settings\n",
    "\n",
    "m3gnet_settings = Settings()\n",
    "m3gnet_settings.runscript.nproc = 1\n",
    "m3gnet_settings.input.MLPotential.Model = \"M3GNet-UP-2022\"\n",
    "\n",
    "# D3(BJ) dispersion engine add-on\n",
    "m3gnet_settings.input.ams.EngineAddons.D3Dispersion.Enabled = \"Yes\"\n",
    "\n",
    "# Optional device setting for GPU or CPU\n",
    "# m3gnet_settings.input.MLPotential.Device = \"cuda:0\"   # is determined automatically by default"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9f5557ce",
   "metadata": {},
   "source": [
    "### 3) Gas-phase naphthalene (geometry optimization)\n",
    "\n",
    "Matches the GUI steps for the gas-phase molecule."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "f20b5f0f",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:06.424112Z",
     "iopub.status.busy": "2026-02-11T10:35:06.424015Z",
     "iopub.status.idle": "2026-02-11T10:35:15.017452Z",
     "shell.execute_reply": "2026-02-11T10:35:15.016867Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:06] JOB naphthalene_gas STARTED\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:06] JOB naphthalene_gas RUNNING\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:14] JOB naphthalene_gas FINISHED\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:15] JOB naphthalene_gas SUCCESSFUL\n"
     ]
    }
   ],
   "source": [
    "from scm.plams import AMSJob\n",
    "\n",
    "s_gas = m3gnet_settings.copy()\n",
    "s_gas.input.ams.Task = \"GeometryOptimization\"\n",
    "\n",
    "job_gas = AMSJob(settings=s_gas, molecule=naph_gas, name=\"naphthalene_gas\")\n",
    "job_gas.run();"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "97cb1514",
   "metadata": {},
   "source": [
    "### 4) Crystalline naphthalene (lattice + geometry optimization)\n",
    "\n",
    "Matches the GUI steps for the crystal, including **Optimize lattice → Yes**."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "51032a7b-5a41-4cf0-acae-346e8c403645",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:15.020060Z",
     "iopub.status.busy": "2026-02-11T10:35:15.019892Z",
     "iopub.status.idle": "2026-02-11T10:35:16.066204Z",
     "shell.execute_reply": "2026-02-11T10:35:16.065621Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADICAYAAACtWK6eAAB7IElEQVR4nO2dd1xTZ/fAz703kyySQAYrCSPsDYKChuHCWQdatda2WmpbR6u1amsNuPdoHVXrHlWwatVat+JeOEAQEJUhe0NC5r3394eJL6/F0f7e0b7N9/PpR5rc3PHc5zzPOec55zwANmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsGHDhg0bNmzYsPF3RqVSUQAAsfxnw4YNG//72Ea7/ycIggBJkpSEhIRxcrn8MIVCqd24caMZAMj/9r3Z+P+D/rdv4C+OdYDhtbW1fVdWVna0oqJiBwCAWq22te3/ALaX+P+DJAgCAYAWOp1+3Ww2hwKABADItLQ0AgCAJEmEJMlX2SZocnIy9h+6Xxu/E5uA/DGsHR756KOPKAiCmEwmUytJkqROp3MdO3asYPXq1XQAAARBSARBSGincr0gMERGRgb+H38CG2+EzQb5F9G1a9dbXl5eYTwez2A2m8coFIpbOp2OQaFQTM/kRqdpaGjQfffdd4b2v4uOju7k4OAgPHLkyPHU1FQkLS3tn4TpXwDyLz7f3wqbgPwBRo0axXVxcUFkMhl669Ytj0ePHvVnMpmzY2JizCEhIZTMzMxry5Yti1Gr1RQAoNFoNBZBEAwajUYzGAwUOp1ubmlpMUql0vrNmzevcHFxSThy5Ijni9exzDSAIAiArZP/V6D8t2/gLwQFAMwcDueL2traRSRJGgsKChAKhcJgMBjQ2toK169fxyoqKoBOp/sDAD0tLU0HAGYA0FhPkp6ejtXU1DDb2tq4PB6Po9PpGK2trZqtW7fKi4uLCQBosc40FtXsn/i9QjNmzBjG9u3b9f+iNvjbYZtB3hwUAAgOhxPt7OysQFG0WigU+lVUVGQ9evSouUuXLpkkSTpwOBwjhmE0qVR6WSqVjpk/f/6jV51PIpF8HhwcPJfL5UpCQ0PpCIJwUBSlWmeaV6lnVjoQGqtaZderV6979vb2Y/bt23dFrVajVueBjTfDNoO8OQQAQGtr67X8/PxrR44ccY6JickSCoXNlo55i8Fg9MBxvO7kyZPDwsPDl7S1tRVMnDhxTnFx8U4EQYZyudzaXbt27QQAHAAIJyenkM6dO3/g4+Njh+P43vv374/es2dPMQBAcnIyplKpmM3NzSwKhcJzcXFxnDNnTodC87KZxt7enm42mz1pNJoAAJC8vLw3GhCTk5OxjIwMAmxqnU1A/gDovHnzxNnZ2eb+/fs3t/vcUaFQIFwul9+5c2fNvHnzYsrKyqYEBASk0mi06eHh4Xaenp4gFovHlpaW5jEYDEVzc7M3QRBuTk5ORjc3t75VVVWz1Wr1NLlcTnnvvfdMCIJooJ169jqhIQhCFx4eru/Zs6cOQRACAFq1Wq0BRdE6ACD9/Pzae89e2vnbedX+9ga+TcX6nSxdulT05MkT1NHRscVoNLKKior8mpqaErhc7uzw8HA8NDQUO3HiRNGdO3dCMzMzNQDQMy4u7kRSUpIxJCQEu3HjBvbw4UOgUCiAoijk5+eT7u7uRGhoKJqfn395w4YNXV+4ZPt39JvOmpycjPn5+TFpNBoLQRDzokWLyCFDhtgPGjTIvGLFCi87O7uTbm5u2+vr67/IyMhoePH3Hahn1KSkpE4nTpy4TBCE9fp/WyGxCcibg5IkiYnF4ukhISHfCIVCQqfTUVEUxQiCALPZDB4eHiYfHx8sNze3bc2aNVIA0Mhksmg2m301MDCQlMlkyOPHj8nKyspzOI7XIghSVVNTMzEgIABpamoyR0dHI3369Nml1WpP+vj40BUKxXaSJKEDFQqxfA7QrvMuWbIkJiUlJYHH4/2EIEhT9+7dC2JiYlhBQUHIzZs3r1EolHfpdLqRQqG0VFRUtL3EphEkJiaWkSSZbW9v/8OBAwe2JCcno3/XtRqbivWGqNVqEQDUNjY2Hi4vL79XWVmpq6qqQpubm58YjUaau7t7Zn19vfDBgwcAALQPPvjgGJ/PH798+fJrTCZzcVRU1NTIyEj07NmzE27evLneel6xWNxsMBhmnz9/PiAnJ8elpqZm/9SpU9+Ty+XQ1tYmQxAkjSRJDEGQ9h2UtAgHwDPBJREEkdnZ2Z3l8Xg0giCmpKenN54+fZp99epV0tPT00Cj0SKzs7MdDh48mLVs2TJniUQiVqvVFBqNZsYwTNPU1KQNCwvTDhs2jDSbzU0IggSaTKYxALDZz8+PfEPvGapWq+F/yRFgE5A3YOnSpSKNRtNi6aTZubm52S8eU19f34nBYKQCQHZVVdW9xsbGRZ6enrmTJ0/+qrS0NN3V1XWKQqFArl275oZhzyJL3NzcwsLDwwcolUo8MjJyWVpa2kAajTY9NDR0o4uLC9y5c+dFdatDUBQlAYA8d+4colKpQKvVwsmTJ3+5evXqBAcHB+T8+fN0FEVBIBDsGTNmzDc7duzYZf0tSZLI0qVLxSwWy76mpkY4ZMgQYWVlJUmj0YrMZjOQJInCM4F8lcvZ+h2Rlpb2O1v3z41NxXoFKSkpVJlM5kWj0UqmTZumtXyMJicnIwAAGRkZoFarybS0NGLHjh2K0aNH1wOA8c6dO5zz589Ttm7d+oW9vf04KpXK8PT0pPr4+CAajaYsLy/vVyaTqairq/MWCoVuAwcONGEYRj18+PDyK1eupNJotO9DQkIE5eXlX546dSpXrVYjrxmVEQAgBQLB4K5du4ZRKJSffvrppzsYhg2JiopKN5vNx/R6vSgkJKQTm80GnU53vK2tbW5bWxsDQZDGw4cP32l3Ln5sbGyRl5cXg0KhtPD5/IE8Hq8FQRACRdG2V7mcu3Xr5isWi00ZGRmP3uCe/xLYBOQlpKSkUBUKhZ9ery9JS0trgg6M1WeaDULOnz/fx2g0tqWlpZW2+9p6fM8uXbqc6NatG+7t7Y3dvHkTKisrgU6nA4ZhQJIkERgYaKLT6bSHDx9eWb9+fewLt/KmRvI/HUeSJJaQkODDYDByXF1dd+3bt28Sj8eboFQqZ7m6utLr6uogJCQEhEIhFBQUTH/48OGG8PBwOH369FiZTLZs0KBBZhaLRT19+vQPZrP5E3d3dzscxxkvRgSQJKmjUqltbDa7dd26ddk+Pj7n9+3b9/GL7QTw/44IQP4btpBNQDpArVZTGAyGP4IgxTNmzGi2CkJHx86fP98HwzDNjBkznlqOA4sBjQEAKZFIQoVC4c2wsDBwc3ODwsJCKC4uPqvT6WoJgqg0Go2fiUQiIAjC2KVLFyQuLm53XV1dVkJCgkEmk/3w4nU76mzW+5s7d2738ePHxzk4OPyIIEh5QkLCk8TERK6vry+alZV1ef78+bFUKjVUoVBsVCqVIf379ycVCgXl1KlTSFFREU6j0Yi2tjaqTqcDPp9vDg0NxWpqai6vWrXqN6peSkoKNSgoiN7c3MwyGo0sNpuN79u3byuPx+P07NlzYGNjo0Gr1Wp+x+LmnxKbDfICHQgHBZ69wBdHLmT+/PneLwgHCfD8pZsBAO3ateudY8eOLQoMDJzq7e2NHTly5NPs7OzvrScRCAR7DQZDfHV19c937txxffLkye5PPvnkPR6PBzdu3PD65JNPVvj6+poKCws1L1kURAGApNPpcj6ff9TBwYGO4/iEH3/8sfnSpUv2Op3OTKPRTCRJRsfFxUWdP3/+ekNDw66ampqIW7duGUtLS5HS0tKa4uLiSQRBNJaWlnYJDAxUl5WVUTQaDbBYrNgJEybM++67776xRCYDSZKwceNGEwCY4Nk6DQYAuL+//yWBQPDWl19+WbVo0SI3LpcrnDNnDvZ7FjfbtR/Z7vkIJycn76CgoGWJiYlvW9Td/4j72TaDwD9G4BeF41W/UavVTlQqlZw1a1bly2YYy4o0HhgY2G3lypUnEhMTKQCwOjU19Ut/f39KRkYGblUZrB2PzWZ/snjx4rUjRoyAq1evnj927NhQDw8Pnl6vp5jN5o46mxGedRTpoEGDHqWmpjJ1Ol3Txo0bd2dlZX0ql8vB0dERAAAIgniK43jq9u3b93G53B2TJk16y9XVVTN37tyPnz59utty28zOnTuX0Gg0BoPBuHfr1q1dYWFhq1xdXR85OjpOOHLkyBVvb+9Yk8lUcfTo0Xzrs3p6evoFBATslcvlcrPZvGDNmjWLrG1gXdzsSD17XRiNSqWinD9/HudwOANjYmIO9ujRw/GLL76oA5uA/MdAAICcOHEi3dHRMaCtra140aJFTWq1GikpKfGdPHnyXA8Pj9ulpaWL/P39zRZBcgQATVpaWhu84kWpVCpKZmamWSKRfDRt2rTvx48fD3l5eTcjIyM7kSSJWla7MZVKhWRmZlpjqHh+fn5LIiIiHBsbG2ceOXIk13q+l3U2CoWCMxiM1jVr1vT08/PzNpvNe48fP57NZDIHOzk5fUSSpM7e3l7h7+8fxOPxQKfT/djQ0PD4/fff/yoiIsIolUqHIQhyGMMwwHGc36VLl0IfHx8ahUKpeuutt4L69OnDDwwM3CiXy5MIgigKDg72sbe3NxYUFGxobW2tZ7PZeElJydsuLi7+vXr1InAcR48fP/7euHHjdp8/fx6dM2eOkST/0URvIjRcLtcwZcoUYzv3tltiYuJDLpcbe/DgwZvWtv039Ynn2AQEABYuXMj38PAQjhw5sshs/kebf/jhh0c2btzYDwCgsrIyycnJ6fh3333nVFdX1/Q64bBgbV+Kj4/Pt507d/amUCjTN23adPPFwMFX2DlvvJI+evRoU15eHjx58oTb3NzMnD59ejkA6FAUNZEkSRGLxdP8/f1nymQyTmNjI7i7u5s9PDwoLS0tmzMzM6ckJSXBrl27PvTw8FjWu3dvM4/Ho5w6dWrH999/P4YgCHB0dFwfHR09vm/fvkZ3d3daZmYmlJeXA51OB61WC/fv38f79OlDuLm5IRcvXly0e/fub9o9n9Vd/NLnsAoNiqJMvV5vPH78ON6tWzfu8OHDG6dMmSLlcrl3uVzu/D179sx+8ff/Lpvm7yYgKABAcnIykpubi+Xl5ZkHDx6s5PF4mQ8fPiw2m81XaTSakcViXSAIor6qqurb8ePHh0ZERMC5c+eybt++PSgkJERrUb/+ZVO8VVgsswqZnJyMJicnQ25uLvkSV2mHK+nW71JSUihBQUH0kpISrp2dHQMAaEKhEJydnZvnz5/vq9PpVsvlcuWgQYMocrkcOXr0KJSWlppoNBppNpvpCIIQXbt2NTs4OKBXrly5um7dum4AAM7Ozp+7uLisCAsLM4tEIkphYSHU1tbmUKnUyrq6OgeDwRCmVCrBYDCASqVq6d2796T8/Pyi0NBQk4eHx43f8xxLly6NHDhwYBc7O7tdLi4uLklJSZe7devGFIlEaHZ29hqxWJxKkiTKZrNbJk+e/DpHwB9+T383AfkNKSkpEYcPH76IIAjDsoCXS6fTlSiKUoVCIdDpdOjZsyeBoih6/Pjx+MzMzPMpKSnUgoICUiQSkcnJyTBs2LDXRr4mJydj6enpZGpq6stWmhkA0D5vo8OgQrVajaampqIZGRnksGHD2jsOXiU0kJycjHG5XN4PP/ygd3Fxmczn8xf4+/ubxGIxtby8vKWurm5qW1tbVUlJSbhQKEz18PAAKpUKCIKARCJZWFtb+01GRgaNy+X+MnLkSJWPj495zpw54xsaGvYAgAEAQKlUZri7u8ceP358QI8ePfrFxMTMfuedd0ChUEBjY+NQBweHnzqICnjx2UihUOi7YcOGu0OHDqVqNJqsPXv22F+9etWDJEli6NChxMWLF6Gtra1XVFTUrbKyMonZbKbS6XSz0WhsM5vNra+yaX4vfxcBQQGA8PPzW8xms3sjCFJApVIfy2SyyoKCAvucnJzpXC4XEwqFbSNHjnx/9uzZmRiGhXp5eR1TKpU0DocDFAoFQVE038nJafL8+fNPdnQRixqBAgD+JqOW1Q5ZtGjRyHHjxs3n8/kHHzx4MEOv15ORkZEmkiQhPT0de5kAvsr9bKFDoXFzc+NrNJqd77//fh8nJyfNkiVLJlVVVe09fvw4v7W1lbpkyZL3zGZzTwA4WV9f36pUKufJ5fJHHh4eYw8ePCieMGHCz927dyecnJz6IwhyHADAz8/PT6lUbndzc/NkMplzFi9evNLd3X33pk2bRiYkJMC9e/cWh4SEzCBJkoIgSIe2g9WpQafTE6dOnXp62rRpkJOTY9i9e/eNrKysrhwOB8LCwqCsrKxNr9e3uLm5fbdmzZoFlrbA1qxZY9/U1MSl0WgUq21mMBg0xcXFDX80aexvJSD29vbvUCiUEAaDEcZgMAJxHOdRqVSqh4cHoChqbmtro0RGRmYuWbIkDgBAqVSqhw8fniqVSg3ffffdTQCIDQwMBAcHh5MajWavXq8fLBAIai2Rsk1gyRl5Q6ydlz1+/PjK9evXswAAnj59mujq6no2JSWFevXqVXZOTk4jWNQmqVSK4zge/MknnwxgMBhnBALBpTcQkn9CrVZT0tLSzAMGDJj07bffrpDJZHoAGI0gyMH09HQsNzcXWbBggdloNDKLi4vpNBqN8u6773ppNJoVfD4/nCTJWm9vbydvb28wmUzlN2/e3MhkMrGHDx+O8ff3l/Xq1Ytsbm5GLl269PaRI0fu+fn5HY6MjGwtKysbuXv37sI3jAqgubu7L4iJiYnkcrmfrV279jabzV4aHR09WaPRNDo6Oop0Ot1jHMflMpnsNp1O//TBgwe6CxcuPAYALcCzdZq33noLzcrK6qrVarWLFi26+kcSxv4uAvJPHD9+3PXEiRN1K1eupEokkhUBAQEfKJVKnCAIpLm5ucbBwWHWt99+uzU5OXnftm3bklksFg4AXREEaXN2dl4uFosTlUolBAcHQ0BAANy8efPKnDlzYjZu3JgQFxfX1cvL61sEQRpf13ktLwzz8/M7+O677/YVi8VaOp1O4Di+Iz09vUtgYKCcyWRmqNXqj6yh52lpaaWzZ892AYAmAFCgKNpEEMQbC4k1GSo0NHRbWlrau0lJSXDz5s3FXbp0mXnr1i1KRESECf7ZvkIwDCNxHAeBQPB9aGjoR3FxcWaJREK5cuUKaLVaYLFYoNPpgMVimTp37mzW6/WU69evL925c+fXL1z+D9ltVCoVYmNjZVwutwgAVl65cqUsMjJyNoVCYTQ0NLQ6ODg4ent7Y3Z2duWlpaVDNm/efAMAYNq0ad48Ho8DADBr1qxbv3cwAfj7LBQiAED269fPLigoKOCXX355kp2dbUJRVBcUFLReKBS+e/v2bQwACARBpCiKbh46dOgHZWVlXjt27IDw8HDs3LlzawEgrLy8vHtzc/MndDp9TUBAgMlkMqEtLS1dBgwYkMflcpVeXl6YwWDompyc3AeeLRa+9IVYKpiY8vLyhqampnbS6/VPEhISBjo7O3/n5uYGYWFhcOXKlZTExESenZ3d9+Xl5ZwrV64Irl27ZuJyufb79u37aOjQocsyMjIAfruQ2SEZGRkkAJB37txRL1682P7UqVOmtra2lQBAhoeHW1UfEtqpZziOUwHAxGQy87RaLVRXV4PZbAatVgtlZWU5BEFUajQaAZVKjSgrK6NSKBTo2rXrxAsXLhS1tbUVdO3a1cxisW680BYdqn/WTrxo0aLYoUOHhnt4eGxCEIRPoVDuREVFoRKJZFpAQMDX8+fPD46Ojt4FAJEAgAQFBSECgcClurr6/Ntvv/0AQRBGaWmpL0EQhurq6ksA0AN5djGbgHQAmZKSQh0wYIBk8ODBN4xG4/MvTp48mSsSiTRUKpVPEEQThUIpLC4uVtrZ2cWIRCLYvXs3WVtbSwgEgtDk5OR1GRkZ0+Ry+U91dXULr1y5wikqKkJycnKKSktLL7m7u/vGxcVBWVmZBzxbXTa9ZtSyrn3ojUbjBYt+vkYqlc5UqVSSx48fw/nz58mqqqrhLBZrOJVKBRqNRqrVaujTpw+JYdj0jIyM5QBAWvV3axG6F2OWrMY9/MPVWnz58uWBly9ffn7MC/fZPqTeDM/C6nfn5ub29/b2TpDL5cZz586l1NbW7gKLcAqFwk8bGxsd29rartXX13cuLy/fkpycDE1NTXDlypX3du3a9VNcXBw+bNgw3QvnB4BnsxsAEAKBwNfT0/OMh4cHrbW19b21a9cKCgoK+Ewm0+zk5ETm5ubO/+STT35et25dglgsPsBmswfdu3fPyOfzaQ0NDfV1dXX36+rqmg0Gw3JHR8c6DMNYAACWtIA37jQAfxMVS61W0zw8PNxWr149BcOwaBRFH1Gp1FI7O7vLJSUlbqWlpUsEAgFCp9PbevfuPWv//v37qqurfXx8fM6EhYWhcrkcqaqqInEcR+vr60tKSkoemUymCK1WyxWJRMSQIUM0TCYzeNmyZW9369atm0AgmLV27drbr9N5rUb64sWLRw8fPnwJj8cDe3v7T/h8vlQmk60ViUQmNptNzc3NPVBSUvKYwWD0VCqVQQqFgvT390dqa2uhra2t35YtW34B+MdqvOXcLxVMa8yY1Z2cnJxMvE71sApg165dhyxbtmxvREQEgaJoMoqih2/evEm1OhXa34dSqdyzbt26EV27doWrV69+d/jw4TlisZiP4zhJpVJxg8Ggp1AobTqdTpeWlmZsdznVzJkzz0+bNg1yc3MN6enpp/Ly8vpJJBIQi8VQVlamNRqNRnd39xW//PLL/pqamv2TJ0/2qa+vN1y/fv2b7OzsFQaDwfqs2N69e5UjRox4YFOxOkCtVlMwDPMfPXp0XkpKSiaHw9FTKJQgOp3ui6LoJ1QqlREXF0cCgLm2tpaLomj/6urq7xAEqWlqappx8+bNZQUFBXhLSwtGEASEhobKlEqlrKWlBTAMA5lMRtjb23PPnj37QXl5+ewff/xxkfXabxKiDgCcJ0+erJfJZCwAgL179253dXXNYTAYcPfuXXTKlCkQHh7OnTVr1jQWizWXJMlzXbt2DX3y5Inp8uXLD93d3Y9OmDAhg8PhzF64cGF+YmJiZxaLRSAIct1SHxhNS0szjx071nP69OkjHR0dLyEIctbSWXCLevZaampqEABAmpubIxsaGigoisL169cnkiR5ODw8nLCsOWAqlQo0Gg3CZrPJmzdvLp43b15wVFRUW11d3bebN2+uA4A6tVqN0ul0uslkYtrZ2dlhGCZYsGABhuM4TqPRDLdu3co7cODA3JKSki40Gu3zbdu25QgEguUikegdJpNJc3Nzs9doNFX37t2bGx0d3Q0ArsfExPiHh4ebtVrtxdGjR1Pi4uIoaWlp+hkzZrA5HA79j/We//EZ5A1iq9gikeg7X1/f96RSKXC5XNBqtUZHR8dPV61a9QMAgJub2098Pn/wgAEDgEKhHF2+fLnZ09NzgLu7O4SHh6Oenp7khQsXSAcHhwpfX9/vIiMjOQwG49H333+/IzU1tcNEo3b399xIHz16dF+hUAgYhkFOTg48evQIvLy8QCaTkQ8ePCARBBmyfv36Q4sXL77/5Zdf+gOA5sSJE7KBAweGRkVFLXF2dg7V6/W3BQJBuK+vL9TV1a1ftGjRJ5ZLUadNm5a9ZMkSHwDA9Xq9J5PJLG4X7vImoPBMoJWxsbGH4uLiwGQyjVy8ePHdl3imrAMABv+wj15aMEKtVqNcLpfe0tLCdHR0pMXFxRnXrVtHk8lkKJfLRcaPH18VHBwszM7OJqRSaUpQUNAkGo3GqK2tNbi6uoqcnZ1NISEhVBRFU9999900AIDevXtzw8PDeTQazUWtVl/9IzPI/6yAvCRkHYVneQVkTU0NkpmZaQ4JCXGsra39icPhBHA4nHMajabCy8trgkQiuS0UCqecPn06KDo6+ttevXqZ/fz8+rm7u5+wt7ef7OvruyopKYnw9/dHz549C83NzTBy5EhISkoCAIDy8vIPXFxctr5qYQz+0YloLBari1arJeLi4gJdXFyWcDgcZmJiInLixAkoLi4mJBIJ2tLSkqnT6aLUajWNx+Oh+/fv/zI1NXUpAIBIJFoQHBw8MyYmxtytWzf0xo0b6NmzZyfyeLxTjx8/1nG53LxvvvmG7ubmRtmzZ8/6M2fOTEpNTYX4+Pg/Es9k1Tx+89v09HQsOTn5eVQzgiCExdDvqHO+cnHz3LlzlFu3btE1Gg0dAOxkMhmCIAjm5+fX+s477/iJRKLFBEEEiUQi6ogRI2gMBoM8duxYU0tLSykAcEmSdCEIAqmurj5/4cKFHmBx9/+eB/2rCwhiGb3+qZ6tWq2mUKnUABaLVfz55583dTByIPDMsKUlJCR4RkZGFkdERNABoNHyvU9ERMQcV1fXwRQKRRcSEsJks9mQnZ292sXF5ZsdO3ZwUBQt9Pb25ohEIuT+/fuP8vPzD33wwQdTp0yZAkajEbKzs2cOHTp00blz5yiv6YTPPSsUCgXMZjM4OTlVd+vWTRQaGkpmZGQgpaWlwGQygcVigb29PTAYDLJPnz6kXq9H7927l5SRkXEcAETe3t5lUVFRVH9/f+TkyZPmwsJCA5fLZREEoWOz2TRnZ2e0V69eUFNT06JWqwUAQKSkpFClUikOljCctLS01wkMYg17/+abb95kXeGVs0ZqaioCz7yHz93Kr4sI4HA4dps3b8a//fZbyvLly392dnaOCwsLMwmFQmpBQUHj06dPd7W0tBhRFL3EYDCMbDabc+bMmX22GQQAVq9eTW9ra/NhMBglFuHoUI3oIGMQSJJEhg0bhv700084QRDg4OCwPjo6evygQYPMd+7coWRmZgKVSq1CEOSB0WgMbW5uthcKhcSQIUM0KIqGf//994NCQ0MH2NnZ3T558uSs+vp6a02rl70UhCRJcHJyEq5atWp3jx49gvl8/od2dnZOXl5e3zs6Opo4HA41Pz//UElJyWMajdbN19c3IigoiPTw8CArKyvRurq6wvv373/KZDKby8vLv2axWAPZbDYuFovNAND96NGjJh6PF+fp6blIqVQSHh4eaH19PeA4vsXBwWHWggULKju6MWvnfcns95tOb40jW7JkyZgPPvjgXSqVuuPEiRM7AQDLyMjA09PTyWHDhiF/sOZWh0Jz6dIlzqRJk/yLi4t3vvfee/Kmpib9yZMnPywrK9tfWlrKMZvNxP379xkVFRWSjz/++N7fyUhHAIAMDQ11cnV1Fet0usJTp061qdVqrKWlxYnL5T6ZPHlyy+7du/kIgjS++GOL+uWHIEhxWlqaVf2yujlxAKAlJyfjly5dulhWVjb+1KlTiNFohIaGhraWlhYulUqNVygUEBYWZg4ODsadnZ25R44cGVVaWppWWlq69E0fIjk5GUUQBHd2do5ksVg9+Xw+rFmz5nBISMhdBEHIhw8fohMnToSIiAjWrFmzpioUClpZWdkls9kcWVRURLa0tJBSqVQZHBx8qrCwEGg0WpNQKISamhpCIBDQu3XrFnX06NHlw4cPv33x4sW3QkNDo2tqajSnT5+udHJy+kCv1w/64osvlpeUlJyiUCh9GAzG461bt+6whLfgryjA8BtVCcMwAgD4OI6vFwqFDACITE5O3o0giBHgeZQtqNVqNC8vD8nIyMAXLlzoPmrUqCBXV9czCIJoXlLiCOC3LmFk/vz5olOnThlu3759bfTo0XtnzJjxNZ/P11IolJsIghBqtVpXWVmJisXioObm5pL29/B7+EvOINZcAIVCkebh4THNaDSWSiSSLfv27Vu1bds29P3339cvXrx46eTJk8dSKJQL+fn57/j7+2s7Sop6xaiCAQDi7u5+iM/n9+nXrx94e3t/O3LkSDWHw5FwOJz1fn5+8e7u7lBbWwtBQUG6t95667PW1la2UqlsnThx4pb09PTXuU+thq9HUlLS5SFDhoi2bdu2++HDh4qgoKAuoaGhBplMRr1//z5BEMSwDRs2HAQAJzc3t/t2dnb2UVFRSHNzc2pNTU24m5tb/8bGRlAoFKBQKIiamhqktrb2rEAgmLxq1arcr7/++sa8efMiAaD13r17wVFRUT1cXFwWKxQKe7FYDCEhIeDk5AT37t1buWTJkimbN2/2S0hIECsUinNv8Eqeh83ExcXlTZ482YVGowFJkrfKy8vX//rrr30CAgIEJEnunD9//lbLbySbNm26PW7cOClBEBdRFFXBM/XttTbC0qVLRQCgtWQWor169bq8evXqaG9vb7h48eL0bt26LQUAcsGCBUoKhdLy5ZdfVr3BM3TIX3IGEYlEJAAAjUa7hWGYFgCYJpNJiSCIEUVRAADB48ePvwAAwDBsYGlpaXhAQECmWq2mMRgMX6twqNXql3pxkpOTISMjwwwA6wMDA/t26tQJOBzOeQBo1mq1zQCQnJeXt6isrMzdbDZvb21tlXM4nA3jx48HFosF69evNyMI8jojnYBng1TRr7/+2uXSpUtijUZzhSRJ50ePHt0jCEKYnZ0NDQ0NqJub24FPP/30tKur62c//vjjCZVK9Xb37t21/fv334kgSFpxcfF8X1/fmY6OjmRkZCR68eJFyMnJSaTRaFkjR45ce/nyZd+zZ8+SMpmMc+zYsakGg2FCZWXl4fr6+i2RkZFJbDZbz2KxqDqd7vPBgwd7yOXyeLlczmlpaVnG5XKnveY5SMsqdeutW7d6zZo16y2DwZAdGxs7hsvlbvby8gJ/f3+4dOlS/NChQzEqlZp55swZ2e3bt6W1tbXm2tra2LVr18asW7fu0us8a5YSTJq0tLQ260B58+bNpV999dWakJCQourq6s0AQC5atMgXAFq//PLLqj+iWln5q80g1vvF1Go1sWnTpgRfX999Go1GJxAINhw7dmze6NGj7QwGgz4rK+vbfv36faxUKgtWrVoVsmnTJvzmzZsBBEE8eV067QvXIwUCQSedTgc6ne4G/GN3qd+8RDc3t++XLFnyUVxcHNy+fXtmnz593sRIf34dy98UADBzuVxPLpc7CsOwzOrqapLL5X4lk8l6urm5GVksVoO3t7cIAFofPXo0ZMuWLWfgmZFeERUVhfr7+yOnTp0icnJyahwcHCQsFgtIkiQkEgnas2dPvLW1Fbt27do7hw8f3u3i4tJbIBAcCQoKAjs7O8rly5erDAaDZtq0aZ5jx46F3Nzcm8HBwZ1+p0v4OR4eHtqkpCS6q6sruXHjRgzDMIRCoQCCIE0MBoMjl8uRXr16IUVFRVlLliyJ3LBhA/Wjjz6yhuhYZyaygxJML9owbADQAQDeQSGNv0U+yG9cdA4ODn1jY2OPenl5kTiOH1ixYsVQy1e0Hj16/BISEpLo4uJiKC0tXcdms9fT6fTGr776ql6tVqNXr16NEIvFj3bu3FkPrzYYn3/3QmMj8KySIHn06FGMzWaTmZmZ9qGhocsDAgLKd+7cuRD+UXj6N7FGHT1fcnKy1Yjt8H7odHpPR0fHjNjYWO5bb72FMxgMbNu2bSU6ne5TOp1efefOna/t7OzeYrFYuFQqNWMY1vPw4cM6JpPZJSAgYFVYWBjh6uqKlJeXIzqdrrm0tDSroqIilyCI4QRBOEokEqRz585PDh482I3L5a5OSEhwZTKZE+fOnXv9DaMCyKVLl34watSoMUwms9He3n66WCye6e7u/q6bmxve2tpqvnHjxvW2tjaEzWaH+Pv7c+RyOREWFoYWFBToKysrfQ4cOFAC8I+VewAAtVrNYDAY3q8owYSiKEqQJAkLFy70BYDWf4VwAPy1VCwiPDycGhAQwJLL5WAymVgXL16cyeFwwMnJiaiurh4yevToPTQaraa1tTXEZDKpXFxccKVSycjOzp54+vTpuffu3WuyjOis2NjY61QqdSAAHH5NvSUSALDk5GR4QcUgAQBPS0uzxviQAFB3586dMXfu3OnoPM9Hw5eMxkS7VW0SAFCVSoVactXBkrd+kiTJ9MrKynE3b94kW1payFu3bvHpdPoBGo1GY7PZLc7OzkRdXR1BkiQ1Kioq/PDhwyu//PLLe5s2bRphMBii8vLyiPr6eqBSqbyQkJAEGo2WYDabwWAwgL+/P6FQKBTOzs4umZmZQ27fvv3sZp51tFdGBbQz0tdKpVI6AMD69evFAoHgkU6ng4sXL+JffPEFLTQ09M6CBQs+Y7FYytbW1ky5XC4+f/68saqqChQKRe60adOm2NnZ/ZCWlobHxsbyAQB4PB6/ubm5JC0trSOXPQAAQRAEzJkz518qHAB/gRnE8qB0pVK5LywsLJEkSdRsNqMkSdJramoAAEgcx8HFxQX0en1DS0tLXktLizgqKsozLS0NvXr1KnnkyBHCZDL9SqFQvtiyZUsBAFBiYmJKJRLJ5oCAAPXRo0ex6dOnE+0WuN5YlbC+iBfSZbEtW7Ygx48f/01W29atW+3ff//9pj/QFCgAIBwOx10sFl9nsVh2oaGhdDc3t2Fz5sw5DgC+Li4u2yIjI309PT1NRqOR2tDQcE0qlY5bsmRJLgBwnJ2d77m4uMjDwsKgvr7+0sWLF138/PwUjo6OZEhICCISicgbN26Qbm5upwIDA1d27doVr66uLvf29n5dHNNzI12lUuV+9tlnriiKQnNzMxw9ehR0Oh2IxWIIDg4mc3Jy9FQqddqpU6f2TZs27fa4ceNcS0pKSuRyuUqpVKZ4eXlNEYvFpQ0NDdcdHBz6uLq6kk1NTeqVK1eus7ZDcnIyYl3otV7/ZSWY/r/86QUELDq/k5PTFC6Xizk5OT0CAHFWVlYBm80eFRIS8p5WqyU5HA5CpVKbFArFpOXLl+9MSUlp/Pjjj+09PT3NX3zxBW4ymegUCqXJYDCsFAgES2/cuHFNKBSePXz48Od/VL+2olar7dPS0pqsakFQUNAMuVze7fDhw32sZWsQBLHbtm3b7qFDh3YjSXItm82eDW/otbFiPb+3t/e+Xr16DevZsyewWKzQhISEuwAAAoEgks1mZ4jFYhmKoucIgmB5eXl1cnBwWCmXy+ds3759xcCBA99PSEgwdu7cuSuNRitycHDYFxoa2v2tt97CpVIpdvr0aUBRFEaMGAFdunQBs9ncqtPpOnG53PzXtJNV7VH6+PiM1uv1T7t37x7GZDJHYxjGiIiIQE6dOgVlZWXA5/NBr9fXOzs7C9PS0giNRmPavHlzzKJFi7IAQOrv7/+Dp6dnn9jYWIiIiIDjx4/r7969GxIdHV2elpamefG6arVa+roSTH+Uv4KAPI8OpVAosH37dvfRo0c/JggC5HJ5lEKhOKfX680cDueH3Nzc6rCwsEWOjo7XL126RFOpVMGRkZFlK1asmEej0fr5+voOZLPZUFlZWebu7i6xt7dvrqurm/L999/v++WXX3bHxsYKcByfIhAI7r1OaCzfIxs3blyVnJw8ksFgrFm0aFGav78/MnPmzLV+fn5dw8LCAo8ePYplZWWZ6HR69xUrVpz65JNPoKmpSbdq1SqnF1QG5CUxTc+bod2/9lFRUe9RKJScy5cvnwLL7ALP1nC4PB7Ps7m5+TYAgIODw9iQkJBUBwcHutFoRMLDw/kUCgW/d+/eO3v27MmQSCS9HR0djwQFBQGLxaJcuXKlqqGhoVStVncaPXo0NDc3w+PHj7vHxMSceY0nq0M8PDz0SUlJdGdnZ3Lz5s0IQRCAoiiwWCyg0+mEm5sbJCQkoCUlJU/y8vI6HzlypBoAIiIiIm7GxMSYvby80AMHDqClpaXAYrEMGIYVUanUGi6Xe8Xe3n5+VFQUq7W1te0Nq8z8bv70AmLtQBMnTnR0dHR0ycnJybbaC0KhcEBsbOwBpVJJ4Dh+aMWKFcMAwF6pVB5QKBQxwcHBmFAoRJqbmxsqKysHXb58eYBIJJosFAppb731FkgkEjhz5gzp6elZ2rt3b5lMJoPq6upzEokk4TWdweowiFiyZMnNadOmAUEQxOHDh4WDBg1qcnJy+iQgIGDByZMn7dsZm9zAwMDjH3zwQTSNRtvz6aefvmNdkGt/Yosx/E+hM9bjOhJapF2Ie7v7er4oZvkOVSgUOdHR0X6DBw/GT548id26dauNyWReM5vNOc3NzW+bTCaxSCQiY2NjS06ePJlApVIXJiQkSHU6Xd6FCxe+yM7ObrO+kpe0CUKSJHTu3JmvVqtXRkdHi+zt7Sc5ODh87eXlNcbV1RXXarXmGzdu3GxrayOYTKZfaGiog7+/P+nu7k6UlpZiFRUV2y5fvjynd+/e8tOnT89xcHCIEQqFSFNT0+GrV69+x2KxRHQ6PYxKpXrI5XJxv379RnO53KeWqib/lkJyf2oBsQrH6tWruUqlUrBly5ZqiURCdXJywmpqauzu3bu3y9XVNS48PJyoqalBi4qK9vJ4vBKNRhOo0+mSunXrRvr7+6OXLl2Cx48fm3Acb3v06BHPksIJQqEQLy0txbRa7a8mkylk6NChDs3Nzd+MHj16SUZGBvpi522HVefmdOrU6ef+/ftH2tvbN4pEIk1paan6zJkz8WKxeDxJkiN37Njxo0qlsjcYDK5cLvfn0NBQNx6P11xWVjZ+3bp1+1NTU7G0tDRzly5dnCQSCffAgQP5L7lmey8eolKpMJFIRHbgXLC6oUkAID09PemPHj0ySKXSjQqF4sPAwEDz06dPsczMTCOTyaRjGAaenp7A4XCIqKgoQiKRULZt29bp6tWrN3/Pu7KuSYSGhn69d+/eeUqlEpYtW5a/f//+bL1eP6yurs44depUal1d3bcLFiz4TCwWB7DZ7MsymYzDYDCgoaGBlMvlqF6vJ8vLyysqKytr2traglAUxTQazTS9Xr+s/fWePHnCOHDgADJ16lQd/JuEA+BPLCBW4fjhhx849+/f9/r1119nh4aGJpAkSTebzWA2m2l0Oh1KS0txgiAIqVSK6vX65paWlgc4jktQFPVQKpUgFArJ1tZWUqfTNVVWVs6+fPnyQ7FYvFuhUDgwmUwgSZIcNWrUg5ycnIN0Op0cMGDA1k6dOj1+k4ohFuOcPX369FURERGarVu3OtPp9KEKhQKio6OhsLAQcnNzL5MkqSQIgms2m+n9+/fHJRIJduLEiZurVq3qlJyczMzIyNAFBQUtdnZ2TsZx/LpAIJi3d+/eXMvMAWvXrp33zjvvDGCxWEsxDNsOv892QQAAYbPZQpFIdFMqlUojIyOpEolk7owZM47R6XS+VCrdFxERwXV1dYWGhgaIiIi4plKpJrS1tfGcnJxq3Nzc7r9Be2AAgDs4OAwbMWLEvvj4eFi4cOG9+vp6f19fX8Tb2xsPDAyES5cuGREEUW/dunUFjuM9XV1dj3O5XMLPzw/T6XQXjUajJ5/Pl5aWloK3t7fZ0dGRkp+fv/bIkSMTAADz8/PDevfuLdbpdJr169c3wr+5BOmf0s3bLi+ZV1tbK79z5052Y2NjZk5Ozg0Gg3GzsbGRUl9fb6LRaG/LZLKxAoEAa2pqAqFQiHh5eW1ds2bNZj6fv0qr1Y6SyWQCAED4fL5AqVR+5eXlteSXX36JN5lMR6hUatHp06fH0en01M6dO8/44osvSAqF8nFzc3NnACh6lR2iVqutI/+I9957b6yPjw9IpdLs8ePH4wqFgsAwjEqj0XAEQWJwHAeDwQAxMTHkwIED0Xv37kF9fb3yo48+ct6wYUM5ACBarZZuNBopADCYRqP9CAC5w4cPxykUSqTJZJrB5XJBr9f/sHv37iNjxoypb99hX7NGQVrc2LVOTk7H+/Tp85FKpQI2m505Y8aM60ajEdra2t67devWvIKCAqrBYNgJAN0FAsGVQYMG0ZhMpqa2tjYSQZDXGek4ACB1dXXpmzdvbrl8+bLbnTt3NnI4nFFPnz7dxWQyseLiYtDr9TQqlbp89OjRY1gs1uIbN24UJCUl+SQkJDSqVKrhCIKYgoKC1isUircCAgKQiIgIkkajjW5ra/u2Z8+eleHh4dxevXqV4fjzifPfWp/3Tycg7YWDJEn5kydP8izuvJUWt257TlOp1EsGgyGgurr6OJ1OjwoICPjh448/HmdnZ/fhw4cPmZ9//vmHWq2WWLNmTWFjY6PSyclpVUhISJG7u7tEIBDQwsPDIxYvXvx+9+7dhxqNRpZWqxXeuHHDrWfPng8tda46xN/fnwQAbPDgwUkikYg0Go14TU1NU1NTU/GdO3c86urqSIIgMK1WixsMBsOTJ09+6dat22A+n48FBAQ0FxcX17u4uBROmzYtZenSpbtNJtNlgiAmGo1Gws7OriY5ORmTSCSU9PT0vN27d1/DMKyTvb391TFjxjS8GCLzupBzy5bOaGFh4Yxdu3YV/fzzz8U3btw4CwAUkiShpqbmIAAcgWejsamwsHBuaWnptaCgoCgvLy92YWGhMwDkw+s1DhIAkLa2tuOWNRS0V69ee8+fP4+1tLS4G41GzNnZ+S0XF5eAurq6oMbGxt3e3t56CoVCnjlzpmHq1KkmAKjLzs5eSqPRhpaWlhI0Gg25cuUKt7S0tGDbtm3m/fv3l4aHhz/h8/nXZDLZgo0bN/5bjHMrfyoV60Xh0Ov1ue3yEzCVSoWIRCJrshO+aNEi7ogRI1CZTNY+YtcjPDx8h4ODQ7SdnR0RHh6OOjo6Io8ePTqRmZkpRRDEVywWtzfSwcvLqywxMdHVw8MDCgsL73t7e4eTJGl+xWj53EhftGjRzenTp4PJZDLl5eXxQ0JCeFKpdGNgYGCvR48enfT09Owjl8tBq9W2PXz4kD506FCMwWCcmThx4nC5XD4lMDDwK4lEcuTBgwfXJRLJPDabDTiOT925c+cKy7VY3bt3PxscHBzB5/PLKyoqJsbFxR3Nzc1FLG2D9u3b1+OXX355+Efa2tq2YMn68/Pzo+Xm5pqEQmFiTEzMVz4+PvknTpyY9gZGenusBbk7CrGhyWSyd5RK5ReNjY1Kf39/bOjQobhGo8EOHz68S6vVLmxtbcWLi4s3Ojo6dhEIBGhzc/NxnU63+9GjRyiGYUF0Ot1LoVBIBg0a1HP69Omt8HcQkNcIR4fHmc1mRWlpaeHGjRsNycnJUFNTg1y8eNFMEATIZLJ73bp1C+zfvz9JEASSmZmJ4DgOOTk5gGHYb4x0g8EQPnz4cGFra+uX77zzzop9+/b9xsPUjvZG+uEBAwaEc7ncZrFY3NrQ0DDj6NGjCY6OjhNoNNqAjRs3NnTq1Gmzm5ubr0ajIQMCAoDNZkNzc/PdFStWhCMIEhoeHr5TIpH4qVQqwtPTE7ly5Qre3Nz8gdFolDQ1NQ2hUChRPXr0MMtkMsqJEycerlq1SkmSJJKamoqkpaWxe/bs+dTR0bHn7t27r71G3XqdcQ/wZh0NSU9PR3/HwipFpVIBAID1/QAA18vL65a3t7dnbGws+ejRIzh+/DjK4XCATqfrm5ubNfX19Q4YhoFQKPy+rKzsY73+DxVH/H/xpxCQPyIcLzsuPDycevv2bZOLi8sP7u7uY0NDQ000Go2an5//RKfTPa6urkYbGhrCPDw8eHZ2dgAA+IgRI/Lv3bt3gEajcQYMGLCuc+fOD3+HkW73+eefb+/SpcvTrVu3ujIYjCFyuRw6d+4MhYWFUFhYeEKr1YYTBCHEMIzo27cvJhKJ4NKlS/D06dMqBEF0Dx48cEQQhNm3b1/Mx8cHbty4QWRnZ2tJktSYzWYkODhYPHHiRHj06BGxe/duwmQyLSwvL59rGaH58fHxtW5ubgPlcvmv58+fR8+fP2+979dWK2lPu3RZ67OBNcbpX0lycjJz9uzZDoMHD05EUXQrl8vVe3t7MwDg/V27dj3kcrlyk8nkz+FwRBQKpVtNTc1qs9m8HgCoAIAnJyeTfn5+Lyvs/S/lv26DdGRzbNiwgUhNTaXCsxq3REfHvUyIsrKycABACYJYnJeXFzpw4MDgmpqa4sOHD8cBQCkAAJ1Od3dzcztBkmTx8ePHJ6MoOqNHjx6zJk+ebAKAkTU1NbHwZkY63qNHj3GTJ08eKpPJwMnJ6d5HH32Ey+VyAkVRKo1GM5tMpl4IggCVSgWSJLGsrCxgMplEXV0d1NXVPcnOzj7AZrOb29raNmRlZUFeXh6h0+lQFov1lE6nj9q7d6+mT58+Be7u7qRAIMC2bduG2dvbz2YymT39/PwmrV+//qbBYCBMJlOTpUL87wqVaU/7GdNq54SFha2QyWR1Bw8eXGApwE3I5XLeokWLNvTv39+BRqPNpNFoN940GsGSDu39yy+/PHn48OG20NDQmcnJycqQkBAyKCjowu7dux+3trZeBgC4ePGi2927dxvHjRunhWcqrQHg2eap/yleaoT+J+io02/cuNGEPEv0N7Vr8DcSDgtEeno6Ul5e/nDEiBFnP/74Y1i4cCG1tbWVm56ejq1evZpuMBgem0ymFkdHx8sAkPfLL7+8e+jQobL8/HxaW1ub6O7du27W1e2X3fv58+cBAMgnT57YHTt2DGpqaqCysrK1sbHxUVZWFvXQoUPk3bt3KRqNBm9ubjY+ffp0WmZm5v5z587hjx49IisqKoDJZLqPHTv2cX5+/iYcxz9qbm6G0tJSpLq6enpxcXExhmG3P/roo09++umn0tWrV6P79+8vzszMnJKTk/MUw7BorVZ7qV+/fhcCAwNRJyen+cOHD/dFEAQ7e/bs9vLy8m3Z2dl8kiQR0rINwKveAwDQ9uzZs6ShoeFsbW3tgNTUVFCr1SiKoq44jr9nKSGEIc/y0ZMdHR2HsVisBK1Wu+Z157diTVajUCjWlAPk6dOnw+7evZu6Y8eOwRKJ5DFJkpRu3boxZs6cGbF161btuHHjWuF3Flr4V/JfnUEQBCFnzJgh1Ol0iuLi4vtyuRwQBIE1a9a8NWbMmGAAOMRms7PhWQLMmwgHAABY1YQHDx6oCgsLsYCAAOcLFy50HjZs2P309HQMAFCtVsupqanxh2cVA8HBwWHil19++ZmHh8etVatWXbQI70vDKjIzM3EAQIqKijasXr3a59y5c0hGRsY4APAEgPVMJjP28ePHVxQKRZS7uzsNRdEUkUi0HsfxiEmTJsllMllJr169DhqNxp/GjRv3a35+/nlHR0eSw+EgDAajcd26dX3u3LkzODY2dp1AIHB8/PgxzuFwXEaNGqX66aefhrW2tn4pEAj6SSSSriqVimQwGKrLly8f2rJlS014eHgsl8sFJpOZjyDIIvLZPosvraiOIAjO4/G60un0aZY4KQ8EQQ4DAAQGBubiOO6elpZG9O7dG0lOTsaOHTt2ddWqVeVtbW0SvV6fOWzYMDI9Pf2Vg+3LMjlra2vv7dmz556lP8Ds2bOBwWB4IwjycMGCBa/cQPU/wX/NBlGr1airqytLJBLxBwwY0H77ZPctW7YUvf/++4jZbC6kUqnearWay2Kx3LVa7f03qLoBYPEyMZnMTj169Fjo6elZd+bMmY/u3bvXbA1t9/Ly+opGoznk5uZOgT9QDuZV11UqlTEeHh7nTp48ScNxPDI0NHRNQEBAJ5PJBAiCgK+vL+Hk5IR07tzZ29/fnx0REbFbLBb7Wox09PLly2Rzc/NHBEHwysrKBvP5/Kg+ffqAWCxGL1y4AJWVlU06nY6Rl5fHcHNzM8fHx1PodLqpoKAAbWtru69QKIJ79eoFZWVlk4cMGfLduXPnsFckbaEAQNJoNM9+/fpdGjp0qD2CIFQul7vrxo0bizIzM1dLpVK5wWDodeDAgccAAFKp1C0wMDAzPDzclSTJeyUlJW8rlcpHAB27nF+T5oympKRgUqkUr6ysxBQKhd9Ljvuv8J8SkOchykwmEzt+/Lhp4sSJH2s0msn379+/wGKxajAMy7Wzsyu+ceOGV3x8/NYpU6botVot/cCBAylSqfQ8i8Uqe9lOQv8frGsdVqs0IyMDeYX3qiMQawnPYcOGIQBA8Pn8pMjIyJ/j4+MdZ86c2QgAFKlU+o2rq+vM0NBQtH///mhVVRWcO3euGkVRsqCggI0gCKtv376ot7c3ZGVl4Xfu3NGhKFoLAEyj0SgJDg4GPp9PPn36lKivr297/Pjx+oqKCtze3n6mQqEABoMBLBYL+vXr99PGjRvLlEqlS1ZW1pisrKy2N3E4wDMh8Ro7duwnDAajtaGhYbBOp/Nxc3PDAgICICcnp7m4uDiLzWZDfX29u52dnXzAgAEmLpdLPX78+OZ169aNS0lJoVp2v33Om9YA+D1bb/8n+a/MIFQqFRQKxXoAGG8wGMBsNtfRaDQ2g8FgIAgCXC6XkEqlZHx8PHrnzp2rW7dujQF41oh5eXmkn58fAvBGNZxQkiTB0ulfuwuUhQ7dnRYPDwKv3hwHBQBCJBJFOzk5nZRIJC7Hjx/XWBOq3NzcCgMDA71CQ0PNJElSrl+/fr2ysvIngiCqNBrNlrCwMIzBYJB6vR5FUTSPSqWOTk9Pv+Pq6prl4OAQwuPxgMlkgpOTE0GlUk/weLxZK1eudA4PD/8ZQZA1eXl5Z0eOHLkoODiYk5KSwgOArOvXrw+IiopqBYCXbkVmDYYcPHjwlLVr1y6RSCT1S5YsydqzZ09ScnIyHh4ejmZlZSEFBQWAoii0tLRAly5dyOTkZPLKlSvIsWPHcnbt2hUM8Cwmy7rzVm5uLp1KpSqtNsdfTTgA/r02CAIApLe3N4fNZp+kUqlVGIY9tLOzy/b19SV2797tptPpCJIkid69e+cVFRWl5OTkYA4ODjOFQuE7LBbLfOfOHQTH8YDPP/98FoqiG9PS0n6zlA7wrAJfXFwc2YHNQCC/LfWCkCQJEonEbs2aNUt79uzpYzKZFgiFwjMZGRnUYcOGGREEgaFDh2Ltazi94axCAADU1tZeY7PZ4Xfv3tVY7oGSnJxMXrp0aVphYeF3AwcOdM7MzLxx+vTpRLBs+OLh4eHY2tq6vKGhAcFxfFZbW1uUj49P1tSpUxfU1NT0aGpq2jdixIhErVaL79692yiXy/vodLo+ffr0OSYWiwmBQODRo0ePxWlpaRGbNm1qNBgMVJPJpCooKFBGR0ffssyUHXY6i81GCwsLe5/JZGJtbW0iJpMZ0NDQAFlZWVhVVRU0NTXBvXv3WuRyOa2xsbFeqVQ6y2QyKC8vx588eRLw0UcfXXZwcHh3/vz5jwCeeZp2795NLysrezJjxowWtVpNGTZsGKlWq5G8vDyyXTrtn1Y4AP69M4h1ddNOKBR+RaFQotlstgJBEDcURSkMBgO8vb3h7t27RHx8POrm5hY1a9asG0ql0kEkEl1NSkryePz4ceXZs2cp3t7eIolE0sBms5cXFxeXubq6euj1+qKtW7f+qFarn/vD36RxreHnbm5uiRs2bDjdu3dvMJvN16lUajTAs3WUrKwsEgDMKpWKEhcXB2lpaeadO3cO7dq1q7yxsfGHQ4cOtbyu7u6LWDMOFyxYcGHmzJldAOAGAHRPSkoyczgcc35+fqRIJLpoNptRT0/Ptzdv3pxhb28/Kjw8fCODwajicDhCX19fnq+vL37mzJkDN2/e7MPlchlOTk5Yv379SCaTiVy4cOFRaGhoeb9+/brw+Xzs3r17h0JCQpItrt8OBdzaHjweL3HTpk2nk5OTyaqqKq1UKuXQaLQhcrn8WxcXF3udTkeIRCK8rq6uhk6nO7NYLPoHH3yAtrW1bRo1atR30dHRm+RyeTSfz591586d60FBQb0Jgrj/ww8/bHtZm1iL9/1ZhQPg3zuDWB+0rb6+fhaKorB48WLn9957rwUApF5eXjcxDOMkJSWRGo3GXFBQMG337t0fjxo1ivrhhx8KvvjiCwQAahAEiW9ubp5ob2//uaen53yZTAaRkZGAYRjweDwyLS1t99GjR1OCgoJKEQQ5/rpGtmwgA1VVVeWbNm3StLa2sk0mU+jp06czDx06dLqpqWl03759TWazedaCBQsOZmZmQlhYWIKHh0eGTCYDiUTiGRoaOj41NbV9UeaO+KfwdBRFCQAQPnjwoHNDQwNKo9Gid+zY4Xf8+PGbCIKAVCoV8Xg8gk6nVzY0NMSQJLkfQZDdZ86cOe/n55fbtWtXXlhYGFlZWYmaTKZkPz8/uHfvHtjZ2eHl5eUYg8EwGY1G+dmzZ7WlpaWU3r17Q2lp6a8AgJ8/f/6l79nSHohOpyv98ccf63Ec5yEIwjx+/PjWurq61E2bNj2QSCQu169fn0GS5BgXF5culZWVGgzD7K5cuYIjCDJ46tSpuuXLl3cuLCwc6+3tvVQkEvH9/f3BwcEB6HR6bG1t7V6BQECl0WjmpqamSIIgdB4eHt9Zi/f9WYUD4N/v5kUAgBw5ciRfIpGwf/nll2qSJHEURQvDwsJ2mkymT2/dukUiCEKxt7cfevr06Zjo6OiiW7duCa5fv47rdLqQqVOnJixfvnyuwWDYqNPpTtvb2/swGAxcp9NR6urqtqWkpEzr1KlTsKOjI9TU1AxHECT9NclOBAAgRqMx/8SJE53v378vam1tbU5MTFxDpVLnBAYGglgshiNHjmzr37+/DkGQxpycHNX58+dJhUJhunv3bu/k5GRHAHhdNZT23hzSouI0nDlzZtznn38+TSqVHl68ePHtdiP4o4aGhgNUKjVKo9EYEAQho6OjmdevXy9va2srrKqqiszKyiJwHMfq6uoMOp2uhUqlFubn58fgOA40Go3KZrOhX79+hevXrz9XXV3Nrqqq+hEAIC4u7lWCbG2PhwcPHuzk5OT0PpVKRRoaGgbpdLonkZGRSGBgICgUijU5OTn3SJIUOTk5AYqiZFRUFIaiqMPVq1cnjxo1qh+NRqvKyckhDQYDTqFQCDabjdJotLEGg2FsVVWV0WAwAEEQBgaDcZZKpWaiKPpw+vTprX9W4QD4D6hYy5cvZ06ZMoVEEOSfAmkcHR2/EAgES1tbW3GSJLNpNJoIQRAnLpeLMJlMwsHBAUlKSoKSkpLmy5cvx1+5cuWum5vbSolEMtnf39+k0Who169fT3d0dAyaM2eOT1xcHGRlZc3q1q3b/DesRfUikaGhoTf69u1rbm5uRg8dOoQyGAzrVsgGe3t7mlwux1UqFeX27dubvv/++5SJEyfSLdsNd1hEG+Clxv1L3coBAQE+zc3NLWVlZRXwbADDJRLJUAzDvp8wYQLvxo0bpSdPnhyo1WrLAKCJxWKNCAkJ2UOlUr+/f//+ueTk5LTY2Fj2yJEj+QBw8/r16wPf0Egnxo0b9+XSpUvn2tvb1y1ZsuTu7t27k4YNG/bcSC8sLARLLg7QaDTw9PQEKpVKPHnyBKmoqPjVZDLllJSUyOh0+tvBwcEEk8lEq6urzx8+fPgdAGiDZzOuaf369WKj0dgwefLklj+zcAD8+1bSEQAgjx07Ri8oKHAKDw8/Gx0dfbBr167L4uPjR/bt27e72Wzu+vTpU2hqakIjIyO1MTExo4uLiz3Ky8u3MZlM1MvLi6ioqIC6ujp7uVx+vFu3bmeYTGaERqMxHj9+HONyufDpp5/qsrKyBixevPj4zJkzd8ybN289ACDx8fGvGjERkiQRsVjM2r9//3q9Xp/Z0tKy9+eff66rqakpPnv2LFZRUYF6enpWl5WVrS4uLj7b0tKCs1gshMfjIU1NTaDT6WJIkkTb7cVttYNIy4rzc4YNG4YjCGJuLxwoihKWBcsX2wy9f/9+vkU4AADM6enpaFVVVcaoUaN+/fzzz9GMjIwqjUZTgmFYE0mSmFar/ZEkyWZnZ+ecurq69PXr13c+ffq0k16vZ5nN5riCggLl66ICrM+Qk5PTKzs7m4rjuJTBYPg1NDSQt27dwo4ePYo8ePCAvHfvXotOpzPX1dWdvXTpUuWvv/5KnDlzhiwqKjKJxWLmkiVLFhUWFo5AECSjoqICf/DgAY5h2EoAKEeebWqqXb16tU9paan2ryAcAP8eFQuBZx3FrqyszO7YsWMVOp3uLJVK7UKn03ujKDqBRqPR3dzcwMfHh8jOziZoNFqsXC5vBYAnUql0RlFRkaqxsVHBZDKJyspKUiaTiRUKhZhOpwONRgNnZ2fS29sbGhoa3gaA9y5cuJB04cKFN7o5a8FoNze3aA6HM55OpwNBEPDrr78qBAIBo6amBsFx3Dx27FhhQEDA8e++++4zBoORzOFwdvn5+WFHjx6t0+v1jh988EHJ1KlTU5cvX74ZAKBfv35hV65cKUpLS2tJTk7GrMF0W7Zs6aFSqQJ1Ot22jIyMJotx35FXzDr7oO3+Rt5++20cAOyrqqoGG41GhE6ndz59+nQnHMdPZ2VloQBAYhhW3traGgQAGEmSLZ6envM+//zzCa6urnu+/vrrO68qsQoAMGzYMBIA4Pbt25+npaWtCA4Obli5cuUoe3v74QUFBSs0Gg1Hq9Wi7u7upqqqqodOTk6evr6+l93d3aPff/99p/r6+typU6faLV26tHHixIkzrl69+pTP5+tJkqTQaLRqgGe58YsWLZLNnDnzedGwP7twAPyLVSzriDBv3jypyWRC0tLSKuG3OjobAJw9PT3vhoSEMDgcjnVfjGMeHh7jZs2aVcnj8cL5fP5lV1dXanJyMnrs2LGf79y5ExYUFOTi7OwMXbt2RbRaLZmXl6fv3LnzTj6ff9vT09OFTqfv9PDwKHzNyIQCAEGj0Xz69et34+233+ZotVrAcRwOHz4MLBYLvL29cTabjd2/f/80n88ffuLECc8tW7ZcjIqKojU0NGwVCoVTgoKCZisUis+FQuGNhoYGrVKpjGez2SV1dXXJ33777U0AAJlMFrJ+/fo7SUlJoNfr05lM5vA/UBkEBQDS1dV1TlJS0qfu7u4Xzpw5M+rUqVNtKpUKy8zMNAcEBBxDUfRCdnb2onZFIuzh2VbRfwhr2Lynp+dqX1/fLkeOHPmwc+fOK6VSaXRzc3MtADj4+PgwFQoFyGSyJotKNzoyMvJbgUBgHx8fT4pEIuTGjRvbSkpKdkskEm97e/vqqqoqD4IgcA8PjxUdqaR/Nv5lM0i7wEMXHMfZaWlpBdCu3qw1RHnOnDkakiQLcByfevPmzTShUIjgOH6MxWJFGY3Gx5MnT15RVla2sb6+vmrgwIEyPz+/sokTJw5FEMQ+Nzf3or29vQ+dTicIgkANBgOzqKgoZcKECSASiUCj0Yy5dOmSPwBoXiEkz430AwcOROXk5DibTCasW7dun9vb2/fw8/NDhEIhduDAAZJCoXSn0Wi3WCzW05MnT1JdXFzIvLy8uMWLF+PTp0+fkp2dvTkkJOQHPz+/+ICAAJNYLJbt379/e9++fd/r1q3bg5UrV3pdvXoVwsPDDTk5OT369u3rlZGR8fh3bmhPAACUlZV9s3HjxuXQrtNbE5L0ev2gUaNGkffu3cMAgEhNTUXnzp3b9JJdnV7F84XVxYsXowAAJpOphSRJFxRF7169ejXeyclpiq+v7xSdTkdQKBTSz88PycvL4w4ZMqSAxWJV3717FzeZTCSdTgcOh0NSqdT3Wltb38Nx3FBWVkaQJElQqdQz58+fXwOW6Nw/M/+SGaS9cBiNRs7s2bMfvOJ6pFqttg8MDBSsWrWq8dKlSwDPdnZC3dzcximVyhU8Hg9jsVhoREQEtbKysjo/P7/HwYMH77u4uKyUSqWTAwICTFqtlnb58uU99vb2HsuWLYuKjo6GnJwcQ1lZmXTUqFGNbxpe0e7/o8PCwq4mJSURGo0GPXDgAAAAsNlsYDAYQKfTSZlMRsbFxaF379693tDQ0C0jI8MoEAg+DwoKWhYbG4tTqVRs165dGgzDuFQq1UiSpJ7NZnM8PT3xxMREyo0bN35Yv379h9aIAAAAS1QA8TqBUavV6Ny5cwlLpwd4/cj7Ug+b1f551eJnuyJ1A7hc7s92dnb88+fPN1va1M/X1zc3PDwcgoKC8La2NuzKlSuZBoPhbnl5uYjJZI4IDg7GaTQaVlpaevn06dPDwbIgavnX9LLr/tn4fwvIy4Sjg/yA58LBYDBk+fn5Bdu3b9cDPHsZ+/fvxy2JObKAgID8uLg4RteuXYktW7aglZWVLQwG447RaERMJlNUXV0dJSkpCfPx8dkxe/bshVFRUYuVSiVfp9Nt37Nnz5bXFGB7rj588sknUZ9++ul3fn5+ZefPn588atSoCzKZTCGVSsnGxsba69ev7yNJ0l8gEISHh4fz/Pz8SC6XSz558gStr6//7P79+1fodDqzpaXlpIODA4XH42EYhh06derUfBqNJhaLxWu9vb1lCoUCd3V1xZ48edIiFouTFy1adPIl74K0Ji29xGZ4sdMjJEmCj48Pe+XKlTNDQ0OdWSzW1KysrKYff/wRkUqluL+//4thNi8rxt3htSIjI4VNTU0eDx8+vAUAxJQpUxxGjhxp17dv30kuLi6fjhgxgn78+PEzp0+f7gsARgCA4ODg/SKRaLBOp0O4XG6/Y8eO/dI+omH27Nm/Zwb962LNA1i0aJHL/PnzfSwfY2q1mv3CoQgAwNKlS1nz58/3S0lJobb7/HnFwPDwcCqCIODm5naze/fuREpKCh4TE2Nms9mNDg4OpLOzM9m7d29iyJAhpu3bt5tnzZqlhT+ANUBx8ODBN0iSJOvr68n33nvvYVhYWKWPjw8ZExNj2rhxo2nixIm9AQBcXFxGeXt7k/Hx8ebu3buTnTt3xseMGUMOHDiQ9Pf3r/T29sa7detGCoVCYuzYseaffvpJBADA4/ESk5KSmlauXIknJSVVBgQEGIYOHUpOmDDhxOeff67q1atX3Lvvvqvu37+/J8CzsIsX2/ZVqFQqCgBAZGTkjNzcXJIkSVKj0Xzf0bFqtRq1njM9Pb1TaWlppze9jpVFixbx5syZE3rr1i0qANCXL1/+wGAwEHq9fi1Jktg777zDAgAkJCTkm549ez6Nj4+v//DDD7vDs0LcFPjn9/2X4A/bIB3NHJGRkcK5c+cejImJ8fzss88m8Pn8A9Y1CbVabQcArK+//jqv/Wna/23JBgSdTvf506dPLxiNRmL06NGUiRMnJr/99tslPB5PVVxcvCk8PJxy+fJl0Gq1dnv37s3gcDjbXVxcgng83iG5XJ73OvXKEnULWVlZmfPnz4+kUCitx48fzxMKhQNCQ0Px0NBQs1arZeh0ukmffPLJ3XXr1u0WiUSdDQbDp87OzmRQUFDb9evXNzo5OX3q7e0t0ev1EBwcDJ06dSI1Gg166NChyQDwNZ1OfzJ16lRKYmIiOn78+PNMJnNeY2PjKolE0lMmk/W0RgU0NTVN5PP5wWlpaeWZmZldBAJBNYIgj173HNbK76WlpdU7d+6EgQMHQm1t7Ud3796tSEtLK7a3tx/j6Oj42N7efsZXX31Vn5aWBgMHDnwrICDgoKurK5Ak2Q9BkF9e5TiwzrbWfBwcx3MjIiJMFApFpdfrfWg0GtTW1n6ydu3ab3bt2tWAIAiYzWaZwWDATSYTQRAEE54FcP5b61f9u/jDAoI827BRYjab7SxqFdLc3DyWIIiubDYb6urqpqnV6kMJCQnmpUuXsgCANW3atA6DDdtBWF7IZZlMdmnEiBFdPT09K8vKyu4DQFVzc/NDBoOBX716dbSdnV2DyWQ6xuPxPn7rrbeOBAUFAY7jHx86dMgPXm2kgzVQrqSkZNrcuXOPBQcHV1dVVeXhOL7i4cOHn+v1eqylpQVwHE+Sy+UPvvrqq+lHjx5dLhQK3x06dCjH19f3wYYNG6bevXt3S3h4+GpPT88Eb29vUiQSoTt27ACTyfTVlClTAs+ePZt17tw5lo+PD1lQUNDr008/nbp27doetbW1g5qamtITExNRsVhswjBMiCBI+vTp05+EhYWNYrPZNU1NTTEA8Og1qaw4ACDV1dU79u7dS7lw4QKHy+VqAgICVgkEAmZgYCA4OjomnDp1yiM5Ofljd3f3qtOnT8cVFhaSCoUCzp49+zEAnICXdFxLGxIvJqtZPs/68ccffzYYDN1YLNbmtLS0hjFjxjC2bdtGeHt7F9bX18ciCOJ4//79CAA4UlNT85eaOaz87pu2drxJkya5CYVCGUEQWcXFxcT27duNdDpdERkZeXnMmDEirVY7/7PPPvtm0aJFPBzHGV9//XUN/L5wc3ZAQEDi/fv3swHgMTwrTUPCb1egqcOHD29YtmwZW6vV6g8ePOg0c+bMNzHSX6xriwEA7ujo+K6dnZ07juPH6urqejg7O0/18/Pj29vbP6TT6VI/Pz/248ePr129erVbVlaWic/nTw0ODl4WGxtrolKp1J07dzYRBGEvEonAZDIBiqKkh4cH2b17d/T27ds3161bFwsARplMdtLT07NHUFAQ/vTpUzQrK+uev79/yPLly8HLywtu3rzZp1OnTr/+AbcwoCiqVqlUqQMHDjQ8fvwYO3ToEIXNZgOKom0kSYJYLLYLDw/HXV1dsXv37o1qaWnZB892ozVam8YS4PgmmZx8AGhcvXo1t6WlRWg2m8tSU1Px4OBgO5PJpCQIorqgoKAC/s0VEP9d/GEBWb9+fbDBYKhiMpl6NpuNjhw5sjU4OHiop6fnptDQUCYAPKmoqBgjEonyrNssW69JvmIf7Fdd0/K/mFqtJvPy8pDc3FwsNTUVnzhx4ujOnTu/6+DgsOOHH37Y3oGR/vyFv3ju9ot6L7m8lMvlzvPy8vqgV69ekJCQQGRnZ6PHjh2bYWdnd764uJja1NR0WigUUng8Hkaj0X4+e/bsfLPZzHN2dl4bEhKi9Pb2Jng8HvL06VOksbHx1uPHjy+0tLS4Yxj2lslkMqlUKqqrq+v4pUuXVvTv33+Tv7//vfLy8lEbN26stzbBS+7tuZE+f/78ebGxsWECgeCwr6/vQQDI8vb2ZvN4PDQnJ6epsLDwLoPB8HN2dhZ5enqCv7+/2cXFhXLt2rXvd+3a9bG1LdpvPGp1qDx58iTvxUQoy/doWloasXXrVvuqqiqZnZ1d/r8jqe2/yR8WkH79+g338vJq1ev1tNbW1lgmk0mrqqp6h8/n85OTk001NTXUM2fOHNu9e/fABQsW8Gk0WptGo9G9pCN2JDRIcnIy+orO+0Z1nNqvObxhaANFpVJBZmYmAc9mFRMAoEqlsiEkJIQbHBwM586dI/Ly8ox8Pp9JEIQWQRC6r68v5erVq2RSUpJJpVK5vvvuuzUuLi69GQzGr1KplKBQKEhZWRkSEBAAHA4HioqKgMViQXNzMzl27FiksLCwZsWKFeIOnvGl92stGB0ZGTljx44dC318fKCxsRE+++yz5bdv335Lq9V6SCQS8q233mq9f/+++86dOxscHR2XxsTEfObt7Q1Xr17Vsdlspru7+wEOhzN74cKF+QAAsbGxfHt7e4NKpZJrNJqitLQ044v3Ym1LtVptT6VS5SaT6cV06H8qov2aNv/T8rsFxNrhpFLpqdDQ0Diz2UxSKJSnBEG01NXVuU2aNMk+JiYGfvzxRyI/P7/Jzs5ujkwm28dkMlGz2WzHYrGMra2tJJVKbSMIwvDFF18YMAyzFhP7p3t7xUyDIJa6TRYdHdplDnakjlhtrd8bwPj8t1Kp9FuhUPixvb29OTg42HDt2jW/rKwsJpPJ7OHj4/NdVFQUgWEYYTabKSaTaeGWLVu+AgAQi8XTWCzWEmdnZyIyMtJ46NCh63Z2dipPT0/C0dERjYiIgLq6OrKyslLXpUuXL4VCYY6fnx/6ww8/XElLSzO+RqgxAMDFYvH777///paBAwdCRUUF5Ofnw9mzZ4HNZoO3tzchFovRvLy8X+Ry+adXr171Xb58+REfHx/K2bNnJw0cODA7PDx8hYuLSxiLxVpTVVXl5OXlFU+hUO7X1NQM2bp1a61FEElrZUvrLPN7Cmn8VfnDM8jkyZM7Hz169OmjR4+aAaAFAMDNzU2dlpaWOmTIEHNubi5l4cKF4OrqCjqd7meCIKZt27atRKVSCUQiUW10dLTEbDaT+/fvJ7RarTE5OZkGAEYul6traWkxvGamebFk5m9GWrVajaamppIjRowI+fLLL/d4eHiA0Wh8y8HBoRB+X3V0FEEQQiQSyb28vIrGjh2LOTs75wsEgqBOnTqZCIIAqVS6TigUfuzo6AiNjY1ZYrE4TCaT3REKhV/u3Lnzjq+vb8nw4cPZQUFBuZ07dw4AgLfkcvlPiYmJSJ8+fZDKykq4ffs2+Pr6wpgxY0AgEEBjY+MvGRkZg1JSUvDX3Kt1B653XVxcPNhsdmt4ePjE6upqp9DQUITD4SCnT58mGhoaUIlEoqmqqtJ8+OGH4n79+iGXL18+0L179yEAACKRaJxSqVyrVCppvXr1AgaDAb/++us2sVg80RL98E8X/TsIB8D/Q0CWLFkSIpfLn7z99tvNo0ePZhQXF5tJkgxtbm4+PXnyZG5TU9P82bNnV0RHR6/w9vam19fX60wmU6mXl5cjQRBblyxZ8sV7770XOGfOnJV8Pv9Kenr6d3V1dVQURekmkwmjUqk4hmH65uZmw4tCo1KpPGfPnv1jaGionclkGpmXl5dbWFiIVFRUkAAAaWlpeHp6Ojps2DA8Pj7+1A8//NDd3d0dysrKFru5uc20GL6/56UiAEC6u7sn+Pn5xT548GDPo0ePiiwjKwEAhEAgiDKZTGhra+tVAPAJDw+fo1AoklEUzbOzs/MKDw+nVFRU3CktLY3fvXt3i5ub2ykvL6/uQUFBeHl5OXblypUrPj4+gevXr+e4u7tDdna29uTJk9Lp06drXmY/vQwMw+aoVKpvBgwYYH78+DHlwIEDwGQygcFgAI1GAy6XC6GhobhMJsOys7O/2Lx583IAAIlEcrJTp0494uLiTI2NjZSDBw8iKIrqKRTKIwqFUmVnZ5eNYdj19evXH9q7d68fjuM5/8vCAfD/EJDt27f7MpnMQqtKQ6FQoEuXLlNEItHcoKAgOkmSjwmC6J2WlsYODQ1dYm9v39Pf3x/p378/lJaWQmlp6dGgoKAeQ4cOpQMAPHz4ME6pVGZa1SW1Wk1jMplMs9lsR6fTGSaTCTObzYhYLG788ccfZ69du3ain58f1NTU/CgWi0e+eJ+W1FmTUCgc0Lt374zExES0urr67ZkzZ/7U0c5OHbWNpf4swD9y21+56gwA/7SzE4ZhfSMjI492794dPDw8yP379yONjY3ZKIqebGpqcjcajYP1er0pLi6OKpPJPv32228rExMT13l4eIBWq/1y3bp1O18Xt2X9ft68eWEjRoyY5+7ufkEikexls9n3lEol297eHs3JyWkpKirKolAoSicnJ+egoCAiICAAYbPZZH5+PlFQUJDIYDDyW1tbhzc3N89zc3OzQxCEvH379pdNTU2lbDY7hkKhuLHZbI+AgIBr/fv3X/H06dNSi23yP80fFpA+ffqM8PHx0TY1NTnU1dWFiMViXmVlZYJEInEZMGCAkSRJ2k8//bRvx44dbwMAx9PTsyY6OpoSEREBGo2Gkp+fX4qiaEPXrl1DZDKZsaSkJGTcuHH5r9rZ6fPPP2d27twZmzFjxpCuXbtuS0hIACqVCo6OjnPmzZt31d3dPRrDsNoffvhhLcAzr0xFRcU7Li4uG5VKJcVgMJwtLCx8Ozg4uPEPRpJiKpUKiYuL6yh2yprfgYNl2+iMjAzEw8OjJjQ01F4ul8PPP/8M9fX1CIPBACaTCV5eXjidTkdGjBgB169fr1q5cqUzPLN5nu/e+pp7tKqcjIkTJ+Z9++238traWlCr1ZuuXr0a39LS4ikSichBgwa1lpaWuq9du7be2dl5vVQqHe/i4oLr9XqUIAhEJBJBeXl5Q1NTU2FVVZXSbDYLjEZjo729vVNJScnzRDcMw2Djxo3uZrO57KOPPvrLxFP9f/j/GumxJpOJbjAYHtDpdH1jY6Mvn89nBgcHExQKBa2qqqrFcfzrnTt3bhKJRAsdHR1nKBQKHABarly5EtjQ0NDo6+v7oZ2dXVZWVtYleEmsUUejN5PJ7GQymShKpdI7MjJydWtrK6dr167g6OgIV65cudXQ0NCAIIjUYDAESiQS6Nevn7m2tpby66+/Tt67d++3lhpO5hfPbx0A+vfvL/7yyy/XhoSEUM1m82Q+n1/8mkW7F8EAAKRS6RcCgWCRUCg09e3b13T48OG3L1++3CoQCMZ7eXkN79SpEyAIAjiOw8CBAydKpdKfhUKh4sSJEzfef/99/es8b5b3QQkPDz/2zTffJD548ABfv359oYuLi6+bm5s5PDycxDAMu3v37i/e3t4zvv766zyRSHTZzs6ui0QiAaFQWNPS0vLQxcUlpra2FnQ6nTkkJARtaGjACwoKOt++fTsLACgTJ07kCwQCRWVl5f1/954cfyb+8Azy2WefRZSVlUFJSUnrrVu3Cixfd3J3d9/P5/NdURRt9PT05FtKWR5TKBSfLF++/N2IiIi03NxceVVVVemLht/vuX77z5hM5vz4+Php/fr1Mzs7O1MuXbpEffz4MRAEAXq9Hj755BPo3Lmzefv27ZSbN29u2Lt37/h2z/9PQXxWwVGpVAu3bNky3d3dHYqLi3coFIr3LALypot2CPKsFpa9q6vrkwkTJtg7OTk97dq1q7XuLzg5OW3gcrnJDAajmMFgXOnTp8+7iYmJaJcuXVgajeaXqVOnDtqwYcObGOkkAPCVSuUkb2/vi0eOHDnr7Ox82dfXt4tUKgW9Xg8mkwkwDDNKJJKNN2/evGFnZ7d5zJgxVD8/v4yoqKjhDg4O7/v7+8/z8PAQJSUlkXQ6HTt69OgWOp0+5dtvvzVduHDBQaVSlb35m/rf4HeHmqSmpiIAQHI4HGpUVBRr//79ORs3bqRaptwbjx8/9uHxeMrm5ubSioqKST4+PtO9vb375Obm5vbo0aNILpdDVFTUhPnz53/53nvvhanV6lU8Hu/mtWvXvkpKSjK+brQEADI+Pt571qxZe4KCgjgODg4DnJ2ds6uqqqinTp2iSKVS5MmTJ6U6nW57Y2MjYW9v/6G3t7eTUCiksNnsivLy8g8mT54cLBaLp8hksodOTk7m27dvG6ZOnapTq9UUa+psYWHh9Tlz5mgHDx7MaGhoyAIA8vz5879nQCFJkkRKSkqaCILov2vXrsSqqqrDt27dQsLDwylZWVnmioqKjyoqKmYAgAYATNeuXTsAAGeioqKgqqoqzmw2MxAEeWXYDPzD9mksLCxMKywsBLVaja5YsWLAgwcPlj169IhLp9Ob3d3dR1AoFEZJSckEmUzW7OjoaG5sbKTu3r07FwDIurq6LYWFhe/yeDxpaWmpuampCbl48eJYFEXfDQkJaWYymTkqleo2i8W6e+zYsV2/M6flL8vvFhBro5SWlt6Ry+WChQsXysxmM7Z48WKz0Wg0u7m5tbz77rvZGIYR5eXlqeXl5Xtqa2uXCIXCAQEBAcEJCQlQVlY27auvvooOCwvrJJfL6QDQ1d3d/RiCIK/bixtFEMQcGxv7vlQqDXNwcIDa2tpvKioqRhmNRieSJGe3tbUVZWVljayrq3sIAECn0y+npaUdS0xMRJuamj66ePFiNUEQm+Vy+emSkpJZ5eXlHDc3t56TJ08+yWQyt61YscI0aNCgxt69e/M0Gg1x8+ZNVK/X9+zbt+/O+Pj4xnb38iYRASQAIGVlZZfKysouWT/Mysoi4NnIj8GzXBhISUmhXrp06dK+fftW1tXVDXJwcPh6y5Ytra9Ll31+IZJEVCoV1s4+qm9tbX3f+n1+fv7yLl26zLKzs0tubm7mBQcHQ0BAANnU1DSzR48eFz09Pe8cOXJkVVFRUbDBYLAjSVJnMpkWYxjWVFpaKmUwGO48Hq+HUqn8XWEvf3X+ZQFkEydOpAsEAiaFQuGYzWYWAKBGoxHGjBlT89FHHzGePn1aGB0dTY2MjEQ0Gg324MGDRyiKNnXt2jVcJpOZHj9+HDJ+/PgH8Oo1CgwAcIFAMGTAgAH74+LigEajgaOj45xFixYdZbPZh+3t7S9VV1eP9fLy4mo0GmF9fb2KTqcv9fHxoeh0uszBgwf37dKli9Hd3X2us7PzTIVCAYmJiWAwGODatWupTCYzi0qlSgoKCpa7ublx+/TpY2ppaaGeOXPmq8GDB6/Oyspi+Pv7N7/EkfAyocHUavXLEqM6ighgAMBvtlOyrO2ApY1e64WDZzWRAeAfAZp8Pv8jd3f37yMjI4mAgAD0p59+grKyMrCzszMxGIzCiooKf51OBwiC6D/99NOu8+bNu9Vuw8y/Hf8vAbHmErxs+p84cSKdSqVynJ2d6cuWLZvJYrE+dXV1JWg0mo5CocT88ssvT729vd/RarVXnj59+uK+3K8coSkUSiez2Uz18vLy7NKly7ctLS3c2NhYcHBwgMuXL7fU1ta2mkwmLoqiHHd3d+jbt6+5pqaGcujQoSk//fTTSoIgQC6XX4iJiencu3dvgiAIytmzZ9G6ujqD0WhEcBynjhkzBmJjY8mDBw+ab9++XSiVSic7Ozs/aG5u5tFoNEBR1KjX6/WOjo4tVCrVkJKSYu6gLV75HJZOb30PRGpqKmLJHPw9DoHXgQAALF68mE2hUKgrVqz41cHBIVwgEICDg8OdjIyMSTQaTcjlcru1trbSJBJJ5+bmZjmTyVRVVlYWgqX8kFqttq4z/c+rVlb+1SHIiNXw7kho+Hz+nJCQkK91Ol3IV199VVleXs6PiYlpefz4MfbgwQMaQRA6Ho+nzczM1HWwh97z81v25iOsayZMJnNuQkLCjKSkJMLV1RU7d+4clp+f/wAAzCiKBrJYLAgKCjIajUZaWVnZ99u2bfsYAMDJyWm7QqF4V6FQmCkUCuXp06dfnD59ehMAQGho6Ioff/xxrLe3t/mnn36ipKenA5/Pb9XpdAvff//9pfHx8eaUlBRq165dmeXl5c50Ol2nVqtb/Pz8oH///lSj0ajNy8t76XO0+7ujweWfPERWff/jjz+OnDhx4mK5XF7AZDIngcUd/CapxRs2bLCrqKiwT0tLq5DJZL3d3d1/HT9+PIjF4qNxcXH9MQwDHMeBTqfDzp07PefOnVufk5PT+Irz/i34j8ToW4PqACCyZ8+elzp37uyVlpZWCvBsrcLFxYXG5XLpTCaTjSAI02AwUOh0utloNBpRFNW/QmhoKpWKePz48SihULhNoVCQYrEYKS0tbSFJ8q7BYEB1Ol1jc3NzD6FQyDCZTPkCgcDT1dX1V4lEkpKWltbYqVOnbCqVKqHT6UaJRKLcu3dvI4Ig4OrqOiA2NjZ90KBBtPv37y/Ytm2be1RU1AhHR0eorKzMRhDkiVwuZ2o0mg3r168/MGnSpCC1Wv12XV1dQUZGxkkURTlmsxl73XMMGTLEf8aMGQu9vLxqCwoKvmxra2tet26dtRgCmZGRgVvdyyNHjryxa9euSARB4Pbt2+PDw8M3kCRJeUVUwPMSTADATktLq7V8RomIiJgRFBQUVl9f/01TU9ODzMxM8tNPP7UXCATypqam+9aCePA3cOW+iv9UEou1oR3i4+NrWSzWoCNHjvycmpqKvSxUwWrTvKnQcLncGUql8isHB4fie/fuvVdZWXnbei4Oh+NtZ2fnXF1dfQkAunTu3HmNu7u7B5vN/iwvL+8dhUIRbDQa9ZcuXZI9ffpUBwCgUChGhYWFrVcqlWwcxw9dv3793fz8/OG+vr6pNBrNJTQ0FBITE+H27duAIMhaX1/fUf3797cHAF1WVpZHREREJUmSyKRJk2gvPgeFQsERBDGTJNl848aN/Tt37oyzJJl96ujouA5FUbAGbyYnJ2MAgPn5+Zn37du3ZNiwYZNjYmLMZWVlvcaNG3fhZVEBHZRgqnjxmPZ0EFv1txcOgP/wFmxBQUG6qqqqoxQK5YFFLfhNzgbAM5XBMoIZoF2Zm/ZCo9VqpYGBgZTQ0FAcQRDjjBkzlnTq1Gmkvb39IYtwPH/Bra2tBa2trQUWO+D81atXAysrK+e6urp+r1AoICEhAQwGA4fFYk2n0WhXMAzjFxUVrXB1deUkJCQY6+rqBpWVlY3NzMxcXV1dfT8oKOiyVqs1V1ZWYiaTidLU1NS/pqYG4fF4wGAw0LNnz9JJyzbNL3sOFovFFggE1OLi4ptr166N69SpE2i12rRr164VTZ48uczd3X0gSZLn9u7dex0sRbIDAgKKCgoKULPZjNXV1SUBwOXc3NzfdOKXlGD6J5KTk7FPPvkEiY+PN78k8PBvLxwAf/4E+lfaNAD/6Gx0Op1+6NChn318fFojIiLe12g0GI1G0/J4PG12drahfcKP1XaRy+UXu3TpEtO7d28cx3HK+fPnoba21mQ2m1GTyYQJBAIyIiLCbDQakcLCwrs4jqfs3bv3jqur6yWJRBLj7u5O6HS6lsOHDysAwLlXr14jWCzW+QMHDpyGFxYhX/EciJeX1zgAoPv7+4e7uLiMKikpIRMTE2kEQUBBQcGqxsZGKZVKddRoNAmenp4QFxcHjx8/hjNnzkTS6fQ7NTU1yPnz53Hrc7UXjq+//jr/ZW0Lv3Pvx78j/w0BweDZzPFHR6iXdjZ/f//1LBYrbuzYsUEVFRUcJpNp15F61tLSom9sbGw7ePDgOi8vr3FyudxEpVIpFRUVX5w6dWozAACbzf5IJpMttre3ByqVCl5eXkAQhAHH8Xl37txZU1ZW9m14ePg7t27dCmxoaMhv73Z9w8SsDo8Ti8WX4uLiugwYMMDIYrGop06dQisqKsCigsGECRPAy8sL37VrF1ZcXNxrx44dJ+EFdWjRokVuOI7bff311/kvuZd/KsH0soxBG3/+GeSNsHh5SD8/Pw+tVkstKSn5TeG69uoZQRB2NBoNvXTpEre2tvZnDMO4NBpNLxKJPPfs2fPcc2NnZxfK4XD49fX1el9f3298fHx6SyQSqK6uvgEANQ4ODn3s7OyGLVu27KeJEyeGpqamDjEYDLednJwOvGEMFTlo0CD/GTNmLJHL5Q0ikWi8i4tLqkKhmOrl5WUSCATUhw8fltbV1RWYzeZipVLZe9OmTa4MBoOcPn266f79+/rg4OAf5HL5Mr1eb540aVLrt99+S0MQRDhp0qQnarUazcvLe/6O22dnLl26lGU0GmUlJSUPbcLxcv4nBORlvGqdxvIdEh0dfV4gEHBJknSIjIwcwmQy67VaLU6hUHSpqanV7X/L4XDeCQkJmWtnZycPCQkBlUoF9+7dAxRF1/j6+r5jMdJNZ8+e9UpMTCx5VTiG1bgeMGDA6d27dyey2WzQ6/UpTCbzB5FI9LNMJuvD5XJzSkpKPigqKroDAODu7r5g9OjRM7p06WK+e/du9OLFi3tFR0cvcHFxuU2lUmfU1NQEc7ncoMbGxlTrjrQdoVar7dhsNvsNqsz87fmv7pP+b8C69QAB0LG+b1XPrGsp/v7+QKPRPJlMZrOdnV2hXq/X4jjO4nA4rNTUVG9LMTczm802fPHFF+lcLveRXC6/6O3tba6traWZzWakoaGhb01NDfB4PKDRaOSNGzdeWyAtIyMDAADy8/Mvrlu3LjEiIgJaW1sX5+fnPxo2bNgXjo6O/RkMxpGioqI7KSkpdnQ6nXXr1q3y3NxcUq/XIy0tLcMaGhpmHDt27Ofw8PDveTzeUTc3N1rXrl2hqqqqj0QiUeE4bjSZTAo6nU5wuVxmZWVlyfbt23NTU1O506ZNq/pXNvz/Kv/TM8irsNSeJfz9/T9AEGQclUoNcnR0dD158mTDi8daHAF2GIZx6XQ67fvvv98hFAqjZTIZCQDae/fueZeVldmrVKohNBrt7KlTpy63+/krc+sBAGQy2RgAMAQHB/d0cXF55+nTp6RKpaICAJmbm3uhpaVFjKIoz2g0SgMCAkClUpEPHjxATp061fmXX365RhAE28fHpzI+Pt6ud+/eeF1dHfXUqVOg0+kAwzDcbDabzGYzSqFQ5v7888+LU1NT8b/Tavj/h/+1GeSNsayfILm5uZsBYLNSqYxBUbTV8vXLXM6NAABcLrcvSZLLRCLRmIqKireWL19OQVFU7+TktLaqqgqJi4sTt19Jt2YZtuO50CAIQpaUlGwDACgpKflRJBJ5xsfHx4jFYjOXy0Xz8/NjGxsbzzMYDIrZbIbKykpzYWEhXldXR+fz+VyCIEClUsHDhw+rHj586IkgCGkymaC8vHzOzZs3z+r1+jx4VjCa+8477xAoipr+SJrB35W/7QzSjjddEENIkoS4uDgsMzPT7OPjE+3p6Xn+yJEjbMtKNjp16lSmUChkoCjK/J0RAZhKpaLGxcWZt27dutTFxeUzpVJp4vP51EePHhmNRmM9nU4va2xsZGu1Wj9HR0doa2trE4vFuEQiSf3uu+9WAEBwVFTUVSqVSmez2XeOHz8eYXX77tmzR4zjOPrOO+9U/wvju/4W2ATkGQg8s1/eJGwVAwDCwcGhb0hIyJGIiAjBwoULm14WjPh7IwIAgC8QCHYqlcq+bDY7p6ioaEJxcXEeANQBAEUsFr9Dp9MdS0tL97u4uIwOCAiYLZPJruv1+i8LCwvTpVKpGABOHThwIAkAID09nVdWVsaeOnVq+b+iof5u2ATkd2L1THl4eHRhMpnb6XR6cFZWVvsU1Dda3OxIaNra2kzu7u6tR48eNVRUVFQJhcJvDh8+vNL6O+uM8MLfHhEREdt5PF6MTCYzxsbG0mpqaurKy8tVZWVlPBcXl3AAyCMIQtzc3Pxk9+7d1/4uyU7/CmwC8v9AJpMx2hc1eAVvLDQAwGUwGOiJEycOubi4ZEdFRc2orKxEnJycmjtQz6y75TJ9fX2r4+LiOD179iSbmpqQn3/+GVpaWoDP5+uNRiOJ4zjNbDbPOHny5LJ2waM2XoNNQP57vFJogoKCihwcHM4OHjx4SmNjo5BGo9E6Us/WrVtnzMvLM0skksv+/v7RPj4+5paWFsqNGzfmlJaWHtHpdI/gH5ELLf/ZR/zrYxOQ/x//6ojX5/ut+/r69tbr9fVPnjy50f6AF9Uzs9lM5fP5+u+++y6Cy+VuxzCMTpLk7atXr0b804nb1euyYeN/DeRVO0EdO3aMnpycLIiJiSlSqVTEgAEDfgV4VjwP/lFE2jYY/gFsjfbn5WX7oQC8RD3r2rXrOQzD4lgs1mVnZ+f4jRs3msEWtv7/An39ITb+S+DQsXAAAJAIgpBW4bCUQ0IaGxur6uvr7zQ1NUU9evSICxav2n/ofv8nsTXe/w5IeHg4RaPROKIoOuDBgwc/wB/b7sGGDRs2bPwdQQAAsW4PbcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNiwYeN/FbVajVp2P7Jhw0Y7/lCSlVqtRlNSUqgAQH1VfrUNG/8DkFiXLt3i33777RGWD16ZupucnIydA6BEdOqkXrNu3SKSJCk2IbHxV+B35qSTmKOjig2A4MHBvgPeeustJQBAenr6qzo7kpGRgcdTqeaoyEiVsbm5GkEQ87Bhw2z58Db+V7CO9kf5CHItB0F+/Pnjj6cdyMrK7gwAQJJkh52dBEBIADQZwPsnOv1YtodH9fEJE3ys3/2n7t6Gjf8AJAoACMDB6wBPdUwm6YSiAADkbwx1a88nn5WugbsAygoAbQbAJcvnttnDxv8a1llkrC9AnR9A2waARoXlu+dCkpycTAMA2CqTMQAAGgEUJoB1dQB+H9FofiRJIrbZw8bfANIdwLjFKiR+fuk0IAELjAmMnj7l86OPSNKtHsBVjyDbSAD3//bd2rDxH0KN/mPG0CkA8C3P/gUAQOBr+EC2cuuu70iVSqJHkK0kgALgmbqltqlWNv4+WISETsoBITdD96u+FBK+QUjewpylwyJJgA0knf5cOP6rt2rDxn8HEgMUAJikExKSdxoqPEhFOZC5AXCCZIIzoDbhsPF3J5nEEASA/nm1pyBmR9WOGEFV2TSOFyAInFP9fTcJtWGjHSQKKECK3U7pTrsUqWXmsNkbNmz8g3+EjtjcuDZsdIhtjcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNiwYcOGDRs2bNj4U/J/1z3dnmEYGJAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<PIL.Image.Image image mode=RGBA size=200x200>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "naph_crystal = ChemicalSystem(\n",
    "    \"\"\"\n",
    "System\n",
    "Atoms\n",
    "C      -0.87388368       0.11073012       2.36589887\n",
    "C       0.15273024       3.05568012       4.76889870\n",
    "C       4.17013024       5.77916988       4.76889870\n",
    "C       3.14351632       2.83421988       2.36589887\n",
    "H      -1.09410897       0.35103804       3.23705766\n",
    "H       0.37295554       3.29598804       3.89773991\n",
    "H       4.39035554       5.53886196       3.89773991\n",
    "H       2.92329103       2.59391196       3.23705766\n",
    "C      -0.12579372       0.96653259       1.59034638\n",
    "C      -0.59535971       3.91148259       5.54445119\n",
    "C       3.42204029       4.92336741       5.54445119\n",
    "C       3.89160628       1.97841741       1.59034638\n",
    "H       0.15997330       1.77874980       1.94280538\n",
    "H      -0.88112673       4.72369980       5.19199219\n",
    "H       3.13627327       4.11115020       5.19199219\n",
    "H       4.17737330       1.16620020       1.94280538\n",
    "C       0.21357858       0.61726152       0.25185835\n",
    "C      -0.93473202       3.56221152       6.88293921\n",
    "C       3.08266798       5.27263848       6.88293921\n",
    "C       4.23097858       2.32768848       0.25185835\n",
    "C      -3.75675609       1.48307682       6.55901941\n",
    "C       3.03560266       4.42802682       0.57577816\n",
    "C       7.05300266       4.40682318       0.57577816\n",
    "C       0.26064391       1.46187318       6.55901941\n",
    "H      -3.47488266       2.30353989       6.89435489\n",
    "H       2.75372923       5.24848989       0.24044268\n",
    "H       6.77112923       3.58636011       0.24044268\n",
    "H       0.54251734       0.64141011       6.89435489\n",
    "C       6.72508585       4.75845021       1.84862605\n",
    "C       0.58856072       1.81350021       5.28617152\n",
    "C      -3.42883928       1.13144979       5.28617152\n",
    "C       2.70768585       4.07639979       1.84862605\n",
    "H       6.21662119       4.18477395       2.37588759\n",
    "H       1.09702538       1.23982395       4.75890998\n",
    "H      -2.92037462       1.70512605       4.75890998\n",
    "H       2.19922119       4.65007605       2.37588759\n",
    "End\n",
    "Lattice\n",
    "       8.03480000       0.00000000       0.00000000 \n",
    "       0.00000000       5.88990000       0.00000000 \n",
    "      -4.73855343       0.00000000       7.13479757 \n",
    "End\n",
    "End\n",
    "\"\"\"\n",
    ")\n",
    "naph_crystal_supercell = naph_crystal.make_supercell((2, 2, 2))  # for visualization\n",
    "view(naph_crystal_supercell, direction=\"tilt_y\", guess_bonds=True, width=200, height=200)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "c5c4e994",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:16.068826Z",
     "iopub.status.busy": "2026-02-11T10:35:16.068566Z",
     "iopub.status.idle": "2026-02-11T10:35:36.296336Z",
     "shell.execute_reply": "2026-02-11T10:35:36.295872Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:16] JOB naphthalene_crystal STARTED\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:16] JOB naphthalene_crystal RUNNING\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:36] JOB naphthalene_crystal FINISHED\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[11.02|11:35:36] JOB naphthalene_crystal SUCCESSFUL\n"
     ]
    }
   ],
   "source": [
    "s_crystal = m3gnet_settings.copy()\n",
    "\n",
    "# Geometry optimization with lattice optimization enabled\n",
    "s_crystal.input.ams.Task = \"GeometryOptimization\"\n",
    "s_crystal.input.ams.GeometryOptimization.OptimizeLattice = \"Yes\"\n",
    "\n",
    "\n",
    "job_crystal = AMSJob(settings=s_crystal, molecule=naph_crystal, name=\"naphthalene_crystal\")\n",
    "job_crystal.run();"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d13e749e",
   "metadata": {},
   "source": [
    "### 5) Lattice optimization results\n",
    "\n",
    "After running the crystal job, you can extract the optimized cell lengths and angles."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "30a9ca2d-753c-4bfe-a997-b3d451a22f3b",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:36.298681Z",
     "iopub.status.busy": "2026-02-11T10:35:36.298563Z",
     "iopub.status.idle": "2026-02-11T10:35:36.302439Z",
     "shell.execute_reply": "2026-02-11T10:35:36.301950Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Initial lengths: [8.0348 5.8899 8.565 ]\n",
      "Initial angles: [ 90.         123.58999997  90.        ]\n"
     ]
    }
   ],
   "source": [
    "lengths = naph_crystal.lattice.get_lengths(unit=\"angstrom\")\n",
    "angles = naph_crystal.lattice.get_angles(unit=\"degree\")\n",
    "print(\"Initial lengths:\", lengths)\n",
    "print(\"Initial angles:\", angles)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "b5e1f4a7-748b-436c-a641-e9bf2da9ffd6",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:36.304352Z",
     "iopub.status.busy": "2026-02-11T10:35:36.304246Z",
     "iopub.status.idle": "2026-02-11T10:35:36.310303Z",
     "shell.execute_reply": "2026-02-11T10:35:36.309889Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Optimized lengths: [7.98578655 6.02392682 8.59150217]\n",
      "Optimized angles: [ 90.         122.86278172  90.        ]\n"
     ]
    }
   ],
   "source": [
    "opt_naph_crystal = job_crystal.results.get_main_system()\n",
    "opt_naph_crystal.symmetrize()  # clean up some numerical noise\n",
    "lengths = opt_naph_crystal.lattice.get_lengths(unit=\"angstrom\")\n",
    "angles = opt_naph_crystal.lattice.get_angles(unit=\"degree\")\n",
    "print(\"Optimized lengths:\", lengths)\n",
    "print(\"Optimized angles:\", angles)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f4313100",
   "metadata": {},
   "source": [
    "### 6) Cohesive energy calculation\n",
    "\n",
    "Use the same formula as the GUI tutorial:\n",
    "\n",
    "E_coh = (2 * E_gas - E_crystal) / 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "a1fbada2",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-02-11T10:35:36.312333Z",
     "iopub.status.busy": "2026-02-11T10:35:36.312223Z",
     "iopub.status.idle": "2026-02-11T10:35:36.317005Z",
     "shell.execute_reply": "2026-02-11T10:35:36.316593Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "E_gas = -11584.06 kJ/mol\n",
      "E_crystal = -23346.80 kJ/mol\n",
      "E_coh = 89.34 kJ/mol\n"
     ]
    }
   ],
   "source": [
    "E_gas = job_gas.results.get_energy(unit=\"kJ/mol\")\n",
    "E_crystal = job_crystal.results.get_energy(unit=\"kJ/mol\")\n",
    "E_coh = (2 * E_gas - E_crystal) / 2\n",
    "print(f\"E_gas = {E_gas:.2f} kJ/mol\")\n",
    "print(f\"E_crystal = {E_crystal:.2f} kJ/mol\")\n",
    "print(f\"E_coh = {E_coh:.2f} kJ/mol\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}