Skip to content

Instantly share code, notes, and snippets.

@iwatobipen
Created February 7, 2025 13:22
Show Gist options
  • Select an option

  • Save iwatobipen/90e07e0b798ad7e03d09df37048fce96 to your computer and use it in GitHub Desktop.

Select an option

Save iwatobipen/90e07e0b798ad7e03d09df37048fce96 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "ed35f386-45f5-482d-ad50-0e4912392353",
"metadata": {
"execution": {
"iopub.execute_input": "2025-02-07T13:11:03.805574Z",
"iopub.status.busy": "2025-02-07T13:11:03.805487Z",
"iopub.status.idle": "2025-02-07T13:11:03.896937Z",
"shell.execute_reply": "2025-02-07T13:11:03.896659Z",
"shell.execute_reply.started": "2025-02-07T13:11:03.805565Z"
}
},
"outputs": [],
"source": [
"from rdkit import Chem\n",
"from rdkit.Chem import AllChem\n",
"from rdkit.Chem.Draw import IPythonConsole\n",
"from rdkit.Chem import rdDetermineBonds"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "4a2fed26-afb9-445e-bcd9-37fd2f0d17da",
"metadata": {
"execution": {
"iopub.execute_input": "2025-02-07T13:11:03.900967Z",
"iopub.status.busy": "2025-02-07T13:11:03.899790Z",
"iopub.status.idle": "2025-02-07T13:11:03.925206Z",
"shell.execute_reply": "2025-02-07T13:11:03.924982Z",
"shell.execute_reply.started": "2025-02-07T13:11:03.900952Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAY10lEQVR4nO3de1STR94H8AkIBFQoKIGAoqLFrS0qysWui6hUKCrlSBc94gUVVjkoKrYqx27XS/W0ikBRq3VbUVBZuqJ2rQpWlArrimBXUZBKBeQWkkBBbgmXXPaP533ZLFaseZIMid/PHxyeIcnzTY/n18nMZIajVCoJAACoy4h2AAAA/YYyCgDACsooAAArKKMAAKygjAIAsIIyCgDACsoovBLKy8t37NhRVlbWp/348eNffvkllUhgMFBG4ZXw+PHjnTt3PltGk5OTUUaBJZRRAABWUEYBAFhBGQUAYGUQ7QAAuvPVV19lZ2ertlRUVNja2tLKA4YBZRT6U1tbO2jQIHt7e9XGjo4OkUjk4ODA5XJpBVNPRUXF06dPVVs6OjpQRoEllFHoj6+v76hRo77//nvVxsuXLy9cuPDGjRszZsygFUw9n3766fz581VbvL2929raaOUBw4CxUQAAVlBGAQBYQRkFAGAFY6PwAm1tbXfv3lVtefLkCaUsAAMRyii8QH5+/pQpU2inYMvHx6e0tHT06NF92s+ePatQKGgkAsOBD/XwAt7e3k/+1+HDh2mHeml1dXWenp5//vOf+7TzeLw+y7kAXhZ6o/ACXC531KhRqi3Dhw+nFUZt8fHxbW1tfRaNAmgEeqNg+BoaGk6cOMHhcDZt2kQ7CxgglFEwfIcOHZJKpfPnz58wYQLtLGCAUEbBwEkkEmYwd/PmzbSzgGHC2Cj0Z+PGjVZWVn0a33rrrU8++eTZWe+BKTk5ubGx0cPDw9vbm3YWMEwcpVJJOwPopaysrK+//vqbb74xNjamneW55HL5+PHjy8vLz549GxwcTDsOGCZ8qAd1SCSSiIiIs2fPfvbZZ7Sz9CcjI6O8vNzZ2TkoKIh2FjBYKKOgDgsLi9OnTxsZGe3YsSM/P592nOdKTEwkhGzZsmUgd5lB36GMgpp8fHxiYmJkMtnSpUsH5l5z169fv337No/HW758Oe0sYMhQRkF9e/bsmTx5cnl5eUxMDO0svyIuLo4Qsm7dOnNzc9pZwJBhiglYKS0tdXd3l0gk6enpixYtoh3nv4qLiydOnGhubl5VVaWPX7sCPYLeKLDyxhtv7Nu3jxASGRlZXV1NO85/xcXFKZXK8PBw1FDQNvRGgS2lUhkUFPTdd9/NmDHj+vXrA2Eyp66uztnZWS6XP3r0aOzYsbTjgIFDbxTY4nA4x44ds7e3z83NZWbGqUtMTOzu7g4JCUENBR1AbxQ048qVKwEBAYMGDbp586aHhwfFJK2trU5OTi0tLXfu3Jk6dSrFJPCKQG8UNMPf33/t2rU9PT1hYWESiYRikiNHjrS0tPj6+qKGgm6gNwoa09nZ6eXldf/+/aioqC+++IJKhp6eHmdn59ra2qysLH9/fyoZ4FWD3ihoDJfLTUtL43K5hw8fvnDhApUMJ0+erK2tdXV19fPzoxIAXkHojYKGJSQkfPDBB7a2tkVFRXw+X6v3amhoEIlEQqGwvr5eLBbX1dWdOnWqoaEhNTV12bJlWr01QC+UUdAwpVI5b968zMxMPz+/rKwsDoej9kt1dnY2NTXV19cLBALmZ3Nzc+9lbW1td3d3n6cMHTq0q6vr8uXLvr6+7N4HwG+FMgqaJxAIJk2a1NjYmJSUtH79+uc9TCaTicVikUgkEAgaGhrq6urEYnF9fb1QKGS6lu3t7f3fyMbGxt7e3s7OzsHBgcfj8fn8/Pz8c+fOjRgxoqioyMbGRtPvDOBXoIyCVnz77bcLFiwwMzNLT0+3trZ+tkSKxWKxWNz/Pz8ul2tvb8/n83k8noODg52dXe+lo6Mjj8czMzPr8xSZTDZjxoxbt269//77GRkZ2nyLAP8HZRS0ZcWKFadOnZLL5f08xtrams/nW1tbOzg48Pl85mfvJZ/PV2NMoKKiws3NrbW1NSUlBXs7gQ7gEBHQlrFjx8rlcisrq4kTJ/L5fNVP30zXksfjsf/maEtLCzO/NHr0aCcnJ0KIs7NzUlLSypUr165d+/bbb7/++uuaeDcAz4XeKGhFV1eXs7OzQCC4cuUKm7VH3d3djY2NqjNLqhNNdXV1LS0tzCMPHDgQHR3d+8TFixenp6d7eHjcvHnTxMSE7fsBeD70RkErTp8+LRAIXF1d58yZ0/8jmUFSZtiUmW4Si8UCgUAkEolEol9++aX/pw8ePJgZJ7W1tVVtP3LkyL/+9a/CwsI9e/bs2LGD5dsB6Ad6o6B5SqXS1dW1pKTk5MmTS5cuZRqvXbtWWlraO7/Uu9izp6enn5cyMTFhpuDt7e17Z5ZUhwgsLCye99y8vLxZs2ZxOJzc3Ny3335bw28S4P+hjILmfffdd++9996IESMqKip6P1DPnz//0qVLzz6Yy+X+6vwS89POzo7N+GlsbOzevXvHjBlz7949S0tLtV8HoB/4UA+ax5zesWnTJtVByYCAgFGjRqkuWmKmm55dtKRBn3zySU5OTkFBwfr160+cOKG9G8GrDL1R0LDCwkJPT09LS8vq6morKyvaccjjx4/d3Nza29vT0tIWL15MOw4YIGxNAhq2d+9eQsjatWsHQg0lhIwbN27//v2EkKioqKqqKtpxwAChNwqaVFFR4eLiMmjQoMrKSm3vS/JSQkJCMjIyvL29c3JyBsIxJ2BI0BsFTYqLi5PL5UuXLh1QNZQQcuTIET6fn5eXxxzAB6BB6I2CxjBfJers7CwuLp4wYQLtOH1dvXrV39/f2Nj4n//8p5eXF+04YDjQGwWNOXTokFQqDQwMHIA1lBAyZ86cDRs2yGSypUuXtrW10Y4DhgO9UdAMiUQyatSoxsbGvLy8P/zhD7Tj/Lquri4vL6+ioqLVq1cfPXqUdhwwEOiNgmYkJyc3NjZ6eHgM2BpKCDEzM0tLSzM3N//rX/965swZ2nHAQKCMggbI5fLPP/+cEBIbG0s7ywtMmDDhs88+I4SsWbOmpqaGdhwwBCijoAEZGRnl5eXOzs5BQUG0s7xYdHT0vHnzmpubw8PDFQoF7Tig91BGQQMSExMJIVu2bNGLJZkcDic5OdnOzu7q1atJSUm044DewxQTsHX9+nVfX18ej/fkyRNzc3PacX6rzMzMefPmmZqa3r59e9KkSbTjgB5DbxTYYjYiiY6O1qMaSggJCAiIjIzs6uoKDQ2VSqVauktJSYlAIOjT2NjY+ODBg/53CAQ9gjIKrBQXF1+5csXCwiIyMpJ2lpeWkJDg6ur68OHDrVu3aukWbm5uO3fu7NOYnJw8ceJEoVCopZuCjqGMAitxcXFKpTI8PHz48OG0s7w0LpebkpJiamp66NChixcv0o4D+gplFNRXV1eXnp5ubGy8YcMG2lnU5Obmtnv3bqVSuWrVKnQPQT0oo6C+xMTE7u7ukJCQsWPH0s6ivg8++MDX17ehoWHlypWYcQU1YPd7UFNra+vXX39NCImJiaGdhRUjI6OUlJRJkyZlZWUdOXIkKipKs68vFArz8vJUW548eaLZWwBdWPAEatq7d29sbKyvr292djbtLBpw/vz54OBgLpdbUFDg6ur625/Y0dHBnGPae0gfc2afubn5mTNnTE1NnzcjX11dPXLkSA3FB5pQRkEdPT09zs7OtbW1WVlZ/v7+tONoRnh4eHJy8ltvvVVYWMjlcnvbm5ubBQJBfX29QCBobm5mfun92dzc/KuvZm1t3dTUZGpqGhoaGh8fr/qngwcP7ty5E2XUYOBDPajj1KlTtbW1rq6ufn5+tLNoTFJSUm5ubnFxsaenJ5/PZ7qWYrG4/66Gubm56oHPzFF9zMl9zAPMzMyGDRum+pR+DoUGfYQyCi9NqVQy3avNmzdzOBzacVj54YcfZsyYYWRkRAgZMmTI7t27V69eXVxc/ODBg97HWFtbP3vyc+9x0Hw+X9//IwBLKKPw0i5dulRSUjJixIhFixbRzsJKfn7+rFmzvLy88vPzmZasrKzW1tY5c+bExMQwXUs7OzumyAI8D8oovDTm258xMTGmpqa0s7DCnMvk6+vLXAoEgrS0NGNj48OHD48bN45qNNAnKKPwcgoLC3Nzcy0tLcPDw2lnYaWiouLChQtmZmbr1q1jWpKSkrq7uxcuXKjBGrp///4333yzT+Ps2bMTEhKsra01dRegTAkGoa2tLSEhobCwsE97dnZ2UlKSBm/0xz/+kRASGxurwdekYs2aNYSQiIgI5rK1tfW1114jhOTn59MNBnoHZdRAVFdXE0L27t3bp3316tUmJiaaukt5ebmxsbGZmVldXZ2mXpMKkUhkbm7O4XBKSkqYFmakYubMmXSDgT7C2Dm8gEKhqK+vLyoqyszMXLNmjVwuX7JkiYODA/NXpVL57bffKvVt9fHBgwelUul7773HHGLa09Nz8OBBQsjmzZtpRwP9g7FRIFKp9NlV5b1Lzaurq2UyGfNIDodjaWn5zjvv9D535cqVKSkpiYmJGzdupBT/pUkkki+//JKoFM20tLTq6uo33ngjICCAajTQSyijr4rCwkKhUCgSiQQCgVgsVv3+Yv+bFnM4HDs7Ox6Px6xIf/DgwV/+8pfAwMAhQ4YQQhYtWpSamhobGztr1ix92UP+2LFjjY2Nnp6e06dPZ1qYQ1C2bt2KFaCgBnwZ1EDU1NQ4OTm98847vaWBcfHixfv373d3d9vY2Dzva4tcLrfP2nLVy5EjR5qYmDCP7OrqmjZt2r1798LDw5l9SQghkZGRR48enTBhwp07dwb+BvhyuXz8+PHl5eXnzp1bsGABISQzM3Pu3LmOjo4VFRX6voQLqEAZNRBMGeXz+TweT7W9tra2tbWV2c5OKpXyeDxHR0ema9n7/cXBgwf/9huVlpa6u7tLJJL09HRm+b1EIpk6depPP/20YcMG5pjlgSw9PX3x4sUuLi6lpaXMuvrZs2fn5OTs27cPA6OgJrozXKApupmpZxw6dIgQ8tprr1VVVTEtP/74o6mpKYfDuXjxombvpXGenp6EkKNHjzKXhYWFhBBLS8unT5/SDQb6CzP18NKioqICAwOfPn26bNkyuVxOCJkyZcquXbuUSuWqVatEIhHtgM917dq1goICHo+3bNkypoVZ57RmzRorKyuq0UCPoYzCS+NwOMeOHbO3t8/NzWUmZwghmzdvnj17tlgsHsh7yDNFc/369cwYbmVl5blz50xMTKKjo2lHAz2GMgrqsLW1PXHiBIfD2bZtG/O52MjIKDU11cbGJjMz8+jRo7QD/ooHDx58//33qoeYxsfHy2SyJUuWYN9PYANlFNTk7++/du3anp6esLAwiURCCHF0dPzqq68IITExMcXFxbQD9sUcYhoREcHs/tnU1MT8n0CPVrzCAEV5bBY0pLu7+/79+w0NDX3a6+rqiouLtXRTqVQ6ceJEQkhUVFRv44oVKwghbm5uXV1dWrqvGmpqakxNTY2NjR8/fsy0MMfHz507l24wMAAoo8BKcXExc97GP/7xD6alvb3dxcWFELJlyxa62VRt2rSJELJ48WLmUiqVMrvT5+TkUM0FhgBlFNhidsK3tbUVCARMS2FhoYmJiZGRUXZ2Nt1sjJaWFmYi/scff2RaDh8+TAhxd3enGwwMA8ZGga2YmJiAgICGhoYVK1YolUpCiLu7+0cffaRQKMLCwn755RfaAcnhw4dbWlp8fX2nTJlCCFEoFMwCgy1bttCOBgaBdh0HQ1BXVzd8+HBCyIEDB5gWuVzu4+NDCAkODqabrbOzk9mPKisri2k5c+YMIWTMmDE9PT10s4FhQBkFzTh//jwhxMzMrKioiGmpqKiwtLQkhCQnJ1MMduzYMUKIq6urQqFgWqZNm0YI+eKLLyimAkOCMgoa86c//YkQ8uabb0okEqYlNTWVEDJ48OBHjx5RiaRQKJgzPE6ePMm03LhxgxBiY2PT3t5OJRIYHpRR0Jj29vbx48cTQmJiYnobQ0NDCSHu7u7d3d26j3ThwgVCyIgRI3pXX82fP58Qsn37dt2HAUOFMgqa1LtHyaVLl5iWp0+fjho1ihDy8ccf6z6Pt7c3ISQhIYG5ZHZ1srCwEIvFug8DhgplFDRsz549hBAejycUCpmWvLw8Y2NjIyOjH374QZdJCgoKyP/u3rRy5Uryv18WAGAPZRQ0TC6Xz5o1ixASFBTU27ht2zZCyMiRI5uamnSW5P333ycqh5gKhUIul2tsbPzzzz/rLAO8ClBGQfNqampsbGyIyraePT09Xl5ehJClS5fqJsOzh5hu3bqVEBISEqKbAPDqQBkFrcjIyCCEWFhYlJaWMi2PHz8eOnQoIeT06dM6CMBs4xQeHs5c9h5Df+vWLR3cHV4pKKOgLcuXLyeETJkypXeWnNn/ycrKqrKyUqu3bmpq6nMM/f79+wkhPj4+Wr0vvJpQRkFb2traXn/9ddXRSaVSuXDhQkLI9OnTZTKZVu9+69atTz/9lPm9u7vbycmJEDLwzzgBfYQj7UCLCgsLp0+fLpfLr169Onv2bEJIc3PzpEmTampqdu/e/dFHH+kmRmpqalhY2O9+97uSkhLmGDsADcI/KdAiDw+Pbdu2KRSK5cuXNzU1EUKsra2Tk5ONjIx27txZVlammxi9x9CjhoI2oDcK2iWTyWbMmHHr1q3g4OCzZ88yjdu3b3d2dg4LC9NBgKysrICAAAcHh8rKShxDD9qAMgpaV1lZOXny5NbW1hMnTmipdMrlcrFYLBKJBAKBWCwWCAQikUgoFNbX19+/f7+lpWXXrl0ff/yxNm4NMIh2ADB8Y8aMOXDgwIoVK9atW/f73/+emXd6WVKptL6+XiAQNDc3M7+oXlZXV8tksuc919jYuLy8nMU7AOgPeqOgI6GhoX/72988PDxu3rxpYmLS569dXV1isbiurk4sFtfX1wuFwt6uJXPZ2dnZz4tzOBwej8fj8RwdHXk8nr29PZ/PZy5bW1sXLVoklUr//ve/h4SEaPMtwisKZRR05OnTp5MnT66qqgoODp42bVp9fb3qp29mAqofQ4YMebZE2traOjg42NnZ8Xi8QYOe+9EqKSlp48aN1tbWRUVFOEsZNA5lFHTn2rVrc+fONTMza2tre/av1tbWfD7fwcGh96e1tTXzi6OjI/MdJPUolcrAwMBLly7NmTPnypUrHA6HxZsA6Atjo6A7AoGgu7vbysoqIiKC6Vry+Xx7e3s7OztbW1vt3ZfD4SQnJ0+cOPHq1auff/55TEyM9u4FryD0RkF33Nzc7t27d/z4ceYsex3LzMycN2+eqanp7du3J02apPsAYKiwGhl0JDMz8969e46Ojsx++LoXEBAQGRnZ1dUVGhoqlUqpZACDhDIKOhIXF0cI2bBhA8U18AkJCa6urg8fPmQ2zQPQCHyoB124c+eOh4eHpaVlVVUVm8ki9u7evTtt2rSenp4LFy4w5zIBsITeKOgC0xVds2YN3RpKCHFzc9u9e7dSqVy1apVQKKQbBgwDeqOgdZWVlS4uLhwOp7y8fCAs21QoFH5+fteuXXv33XcvX76M9U/AEnqjoHXx8fEymWzJkiUDoYYSQoyMjFJSUoYNG5aVlXXkyBHacUDvoTcK2tXU1OTk5CSRSO7evTuglhmdP38+ODiYy+UWFBS4urrSjgN6DL1R0K6DBw92dHQEBAQMqBpKCFmwYMGqVas6OztDQ0P7/8I+QP/QGwUt6uzsHD16tEgkysnJmTlzJu04fXV0dEyZMqWsrOzDDz9kJsEA1IDeKGjR8ePHRSKRu7v7AKyhhJDBgwenpaWZmJgkJCRkZ2fTjgP6CmUUtEUulzOnd2zZsoV2lueaOnXq9u3bFQrFihUrJBIJ7Tigl/ChHrQlIyMjJCRkzJgxZWVl/exiR51cLg8JCQkLCwsKCqKdBfQSeqOgLfHx8YSQDz/8cCDXUEKIsbHxuXPnmGX5paWlff6alpaGRVHQP/RGQStu3Lgxc+bMYcOGVVVVDR48mHacF8vNzfXx8Xl2h3w/P7+qqqpHjx7RCgYDH3qjoBXMxHd0dLRe1FAANlBGQfN++umnzMxMCwuLqKgo2lkAtA5lFDRv3759CoVi5cqVWt3THmCAGNBj/6CnbGxshg4dqo9ndZw8efL27duqLWVlZWZmZrTygF5AGQVWhEKhkZERj8dTbdy9ezdz2hKtVGr7+eef+5xR2tLS0ufdAfSBMgqsvPvuu0OHDs3Ly1NtzM7ODgwMzMrK8vf3pxVMPbt27frVmXpaeUAvYGwUAIAVlFEAAFZQRgEAWMHYKLAlkUju37+v2vLkyRNKWQAoQBkFtv79738PtC2Z1eDp6VleXm5nZ9enPTU1VSaTUYkE+gJlFNhyd3f/5ptvVFtycnIiIiJo5VEPl8t1dnZ+tt3e3l73YUC/oIwCW88WoIcPH9IKA6B7mGICAGAFZRQAgBWUUQAAVjA2CqxERkY+u3OHi4vL9u3bx40bRyUSgI5h93sAAFbwoR4AgBWUUQAAVlBGAQBYQRkFAGAFZRQAgBWUUQAAVv4DQ9EKA/Ra1MsAAAEWelRYdHJka2l0UEtMIHJka2l0IDIwMjQuMDkuNQAAeJx7v2/tPQYgEABiJgYI4IHiBkY2Bw0gzcxCLM3IAKIZGImluRkYM5gYmRKYmDOYmFkSWFgzmFgZEhjYGBjZGZg4GJg5GVi4GFi5GZwYgarZGFhZmJkYxWeBdEIxA4/tcT97u86J+zMkptuc2b/R3o+VzU708bTdPjtE7N96LrJ/cnXPHpfjfvtdOifa6/JP37Nv/8b9Qqxs+7SCp1u77BDZ/9Jz0f6bT/bYvP+w3b6mjPFAfgSbrbWcnMPC/Dw7zjK2PVO9W+xLYgUcpnNp7BX8uH1/SRmjg1se215jObkDp/Pz9oXVs9lM9G7ZXxgrcOA/s4atGACXHUp10nZqYAAAAXd6VFh0TU9MIHJka2l0IDIwMjQuMDkuNQAAeJx9k0luwzAMRfc+BS9gg6QmapnEQVsUSYA27R267/1RUlZMZVPLQmz56X8OygR2fazvP7+wX2GdJgD85661wndAxOkC9gDH88vbFU73w/Gxcrp9Xe+fQGw32nhmD/fb5bFCcAJcBEthgZkWCgWRdAlJYoy+lRWkJYjknPUzBRRBmBXEmDE5GJpiEqyh6A7ORYQayAETORgVnHfv3Xru3kPYychuHg0w91JhMydxMm+a3X3e7bv7QBZ4Vc+onuo0s35nScXAwEVrvIOiIC+xYKqa5cKBJcqWuDAWB2tTxBCS6vDC+oMtWoyEnB0kVHLu5hG6d2qoqvOQOZGh3T4Y0f27/RCn9rqptgCybWoRBNgCoAE9X9enQ7Adi+PtuvqxiDrZm69hQvAWR5veSBvJu6UvkL0lUSd53UlH8eqqDYjXkFS5DpVqkYz1IBUnGrImk+cxtzETe3/8OfR5+gMvDaYvCUvH9QAAAQJ6VFh0U01JTEVTIHJka2l0IDIwMjQuMDkuNQAAeJxtjztqQzEQRbeS0gZJzE+akbUBp0tvXIRXPYixiV168dHn2SBII3R1RsO5p+N5wWV3Op73/5xYz4/nDgMwx+QoUIzqPAQQBIpaIESDzOowUFKzzoghIpV6NVAlqxB5fEMTzrl4DMIKUjcCkcXOmIwaYbOU2gsymGZXEUjCuqt4CqIQc6NU58U6NYLc6NvGv3Q2G+O2eJTwW4tRArOUt6cfosNThMtL0w/P5Lqm5kxleLbhpjmKV00QxdIt2/Dm6DfJZHHvvh/Xy9fv9XaAsN4/L7efdVkfAQ88xzhHnWOeI+KUn38JyG6JkayncQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x75ee4ff1ff40>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"smiles = 'c1ccccc1'\n",
"mol = Chem.MolFromSmiles(smiles)\n",
"mol = Chem.AddHs(mol)\n",
"\n",
"AllChem.EmbedMolecule(mol,randomSeed=42)\n",
"AllChem.MMFFOptimizeMolecule(mol)\n",
"xyz = Chem.MolToXYZBlock(mol)\n",
"mol"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "40aa3313-0853-4d29-b7ad-011ab1797c64",
"metadata": {
"execution": {
"iopub.execute_input": "2025-02-07T13:11:03.925656Z",
"iopub.status.busy": "2025-02-07T13:11:03.925494Z",
"iopub.status.idle": "2025-02-07T13:11:03.929817Z",
"shell.execute_reply": "2025-02-07T13:11:03.929485Z",
"shell.execute_reply.started": "2025-02-07T13:11:03.925647Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAVk0lEQVR4nO3dfVRT9/0H8O8FAiokSIxioQSQB0kV2qo9WGXyVMUJc/bBbbWznq5CLRzbVbDoERCK2zkgciob87SzMrXD08ksnYB2k2pV1FYsRXnGYSDlQYQEwkMSIrm/P3LG7+5CM+DmgYf36y/y4ZN8P5y2796H3HspmqYJAABMlpWlBwAAmN4QowAAnCBGAQA4QYwCAHCCGAUA4AQxCgDACWIUZqb09PSioiJWsaysLDU1ValUEkJKSkrS0tJYDWq1OjU19fLly2aaEmYExCjMTOnp6efPn2cVy8rK0tLSent7CSHFxcWjY1SlUqWlpSFGYUIQowAAnCBGAQA4QYwCAHBiY+kBAEzl5s2bv/3tb5mVO3fuMF/SNM1q0Gg05pgMZhbEKPw/uVze09Pj6elJURSz3tTUNH/+fKFQqNVqZTKZSCQSCATMhkePHg0MDHh4eJh13P+lq6uroqKCWWltbWX1sBoeP35s8rFg5qEB/iMhIYEQolKpmEW1Wk0IiY+Pp2m6vr6eEJKTk8N64/bt2wUCgfkGHQcejxcTE8MqZmRkEEJaWlpomo6NjaUoitUgl8sJIcnJyWaaEmYEHBsFAOAEMQoAwAliFACAE5xiArbKykpbW9uRl1qtltXwww8/sM7MKBQKc0wGMCVRNB4iAv+xd+/erKysMX8VHx+flZXV0NCwdOnSMRsEAoH+IsspQiaTOTg4ODk5MYt9fX0KhcLV1dXa2lqhUPT397u5uTEbaJqWyWSOjo6Ojo7mnRemMWyNAlt9fb2dnd3Iy6GhIV9fX2ZDWlrajh07mJU9e/ZcunTJTPONDysf9fh8Pp/P1//s5OTECllCCEVRYrHY5MPBzIIYBTaxWDxnzpyRl6O/ke7k5OTu7s6s2Nvbm2MygCkJp5gAADhBjAIAcIIYBQDgxDo1NdXSM8BUYW1t7ePjExwcbGX1//9/pSiKx+OFhoYuWbKEoig+nx8aGuri4sJ64/Lly9esWWP2kQEsD194AgDgBDv1AACcIEYBADhBjAIAcIIYBQDgBDEKAMAJYhQAgBPEKAAAJ7g1CUxXnZ2dJSUlTU1NVlZW/v7+UVFR+htTnTt3bnh4eOvWrczm1tbWwsLCTZs2eXp6WmhemLHw9XuYlo4dO5aQkKDRaHx8fHQ6XWNjo7e3d1FRka+v75o1azQaDetZyl999VV4ePi5c+defPFFS80MMxV26mH6+eKLL2JjYwMDA6VSaW1tbX19fU1NjUgkam5utvRoMBthpx6mn4MHDy5YsODvf//7yH2X/fz8bty4YdmpYNbC1igYgUqlysrKCgwMFIlEixYtCgsLKyws1P9qx44dP//5z1n9d+7ckUgkFy5cmMRa7e3tlZWVmzdvHn3vegCLwNYocDU4OBgeHv7NN99s3rx5+/btNE0XFxe/+OKLf/7zn3fu3Nnc3Nzd3T36LXV1dUqlchLL3b9/nxDi5eVloKejo+N3v/sds9LU1DSJtQDGAzEKXKWlpd26deuPf/xjXFycvrJ79+4zZ85ERUWZYjmVSkUIMbwpKpfL8/PzmZWBgQFTDANAEKPAkU6nO378+DPPPDOSoXqvvvqqiVbUP7Ozvb3dQM9TTz015pl6E40EsxyOjQInDQ0Ncrn8Jz/5idlWlEgkVlZWNTU1ZlsRwDBsjQInnZ2d5EeeZjyirq6Odbf8oaGhSa8oEAg2bNhQXFzc0NDAevIzgEUgRoET/eNGhoeHDfQsWrTowIEDzMq///3v7OzsSS+akZGxevXqDRs2HD16NCIiQr9xeurUqeTk5Amdvq+trT148OD69eujo6MnPQwAYhQ4Wbx4MSFEJpMZ6BEKhbGxsczKtWvXuMRoQEDApUuXoqOjt2zZQgihKIqmaW9v79dff338MarT6aKjox0dHb/77rtJTwJAEKPAkZeXl7Oz8+XLl3U6HfNBeKa2Zs2a6urqe/fuPXjwwNbW1t3d3c/Pj6IoQkheXt7oS5yfe+658vJy5tek/vSnPz355JNr167FYVbgCKeYgBOKomJjY2trazMzM5l1/deSJoSmaYVCMTg4OP63+Pv7b968eePGjRKJRJ+hhJClS5f6+fmxOvl8/sqVK+fPn69/2dbWlpmZeeTIkYkOCTAaYhS42rdv3wsvvLB///6wsLDDhw8fPXo0Ojra1dX1m2++mdDn/OUvfxGLxbt37zbRnExvv/32/v37XV1dzbAWzHjYqQeubG1tS0pKTpw48dlnn+Xk5NjZ2YnF4gMHDuhPoy9btqy3t5f1FoFAsHr1apFINFLp6OhIS0t77733WltbTT3whQsXSkpKFi1a9NZbb1VXV8vl8i+++GL0FasA44Qb5cGUsHXr1pCQEBsbm2+//faTTz4x6VptbW0jx0MLCwubm5s//PBDw1eXAhiArVGwvKKioubm5l27dh0/ftwMy7m4uIx8j7W2tnZ4eBgZClwgRsHCent733777VOnTimVysHBQY1Go1Kp5s6da57VzXMoFmY2xChY2I0bN5ydnffu3UsIefTo0eDg4PHjx5FuMI3g2ChMIR999JEZjo0CGBe2RmEKkUgkZtudBzAWbI0CAHCCr98DAHCCGAUA4AQxCgDACWIUAIATxCgAACeIUQCjaWlpqaurYxU1Gs29e/fkcjkhpL+//969e319fayepqYm/YOjYTpCjAIYzXvvvRcaGsoqNjU1BQQEnDlzhhBy/fr1gICAr7/+mtXz+uuvv/LKK2aaEowNMQoAwAmuYoLZa2hoqKysrK6uTqfTLV++PCgoyNramhBy69YthULx05/+lNnc19f3z3/+c8WKFZ6enhaaF6YoxCjMUqWlpb/5zW9aWlr4fD6Px5PL5f7+/gUFBb6+vhkZGbdv3/7hhx+Y/S0tLa+88sqxY8d27dplqZlhakKMwmx09+7dyMhINze3a9euBQUFEULKy8vfeOONq1ev6m/aP2larfbatWvMSnNzM6unpqbG0dGRWVEqleZ8ICAYF2IUZqOUlBSdTldUVLR06VJ9ZdWqVRUVFTY2XP+L6O7uXrduneGexMTE0cWnn36a49JgKYhRmHW0Wu2XX34ZHBw8kqF63DOUEOLs7FxdXc2sNDY2Pv/888xKfn7+hg0bmJXIyEi1Ws19dbAIxCjMOs3NzWq12sfHx0DP4ODgiRMnmJW2trbxfDhFUQsWLGBWOjs7WT18Pp/VY5QEB0vBPzyYdQYGBgghCxcuNNDT09OTkJDArAwPD5t2LJi2EKPwozo6Ourr63U6na+v78gj3dva2vr6+li7w1qttrGxcfHixUKh0BKTTgyfzyeEPHz40ECPi4sL60x9dXX18uXLTTsZTE84OQhjaGpqioiIcHFxCQkJCQsLc3Nzi4qK6ujoIIS8//77q1evZvW3trYuW7YsLy/PEsNOmFgsnjdvXkNDg6UHgRkCMQpsnZ2dQUFBt27dys3NbWlp6ejoyM/Pr6ioOHbsmKVHMw4bG5vIyMirV6+OPK3eWN54440PPviAVXR2ds7OztZ/rUoikWRnZy9btozV88477+gf6gfTEg3T06effvrZZ5+xik1NTdnZ2VKplKbp77//Pjs7u7u7m9Xz8ccfnz9/3sAnx8XFEUIuXrzILHZ2dup/eO211+bPn896y4MHDwghWVlZk/hDLKKmpmbu3Lmenp5fffWVTqejabqrqysnJ6e1tZWm6S1btri6urLeUlVVRQg5duwYs9jd3Z2YmHjq1CmzTQ5TELZGp6vDhw/n5OSwitXV1Xv27KmtrSWEXL9+fc+ePaOPAKakpBje+y4oKHj66acjIiKYRcMnZKYdiURy8eJFQkhYWJijo+PixYtFItGBAwfu3r07oc+Jj48vKyv7/PPPTTMmTA84xQT/5dGjRw8fPgwPDzfQMzw8fOPGDWZFf9h0elm3bl19ff3NmzcbGhqsrKzc3NxWr16tP/uUlJTU09PD6nd3dz9//nxAQMBI5fLly1VVVfHx8X/729/MOjpMMYhR+C+9vb2EkCeeeMJAT19f39q1a801kQnxeLx169aNvuho5cqVo5sdHByioqJGXg4ODr711lunT59uaWkx7ZQw5SFG4b/MmzePEKK/x/CPcXR0rK+vZ1ZkMtlzzz1n2smmmOTk5MjIyMDAQMQoIEansZaWltTUVGZl9Jd4cnNzRSIRs9Lf32/gM52dnR0cHJqamgz0UBTl7OzMrKhUqvEMPGNUV1fn5OTExMTs27evrq6utrb23LlzL730kqXnAstAjE5jcrm8sLCQWVEqlaye0tJSOzs7ZsXwtdvW1tbh4eElJSVSqdTDw8NIk840QqEwPz9f/7NWq21vb3d3d7fsSGBBiNFp7Jlnnrl+/TqzUlRU9LOf/YxZOXfunEQiYVYMH/ckhCQlJRUXF7/88ssFBQX6WxTTNH316tWnnnpqQufrtVrtyZMnXVxcNm3aNP53TQtPPPHE1q1bR14+ePBgzMOpMEvgC0/AtmrVqr/+9a/379/38fFZtWpVSEjIk08+GRIS8uWXX07oczIzM48cOXL06FETzTlFbNq0acZcmACTg61RGMMvfvGL0NDQwsLC2tpaiqJeeuml0NBQf39/Qsi2bdtGXwwqFAoPHz4cHBw8UmloaPjkk08OHTp08uRJs45udvb29vb29paeAiwJMQpjW7hwYXR09Oj6mHvoAoGAeT8knU63c+fOzMxM/Xl/gJkNO/VgfB999JG9vT2eGAyzBEXTtKVngMm4f/++lZXVkiVLmMW+vj6pVOrp6eng4CCXy1tbW319fVln6uvq6ubNmycWi000WHt7u0QiiY2N9fDwqKqqKi0t/cMf/hAWFmai5QAsDjv105W3t/foIp/P1x/BJIQIhcIx7/7p5+dn0sG0Wu3+/fsJIQqFor+/X6vV6m+TDDBTYWsUTKikpOTo0aMTPcUPML0gRsGEurq6ZDLZs88+a+lBAEwIMQoAwAnO1AMAcIIYBQDgBDEKAMAJYhQAgBPEKAAAJ4hRAABOEKMAAJwgRgEAOEGMwsxRWVkZHx+/cePG0NDQXbt2jTwF+tSpU9u2bWM1q1Sqbdu24dnIwB1iFGaIrKysFStWfPrppwsXLvT29i4rKwsKCjp8+DAhpLy8/MyZM6wL9rRa7ZkzZyorKy00L8wcuMMTzARXrlx5//33w8PDP//8cwcHB0IITdOZmZmjb9QPYHSIUTAhrVYrlUo1Go1YLBYIBPqiSqVSq9VOTk6sZoVCMWfOnLlz505ioezsbBsbm5MnT+ozlBBCUVRiYiKX4QHGCTv1YBJKpTIuLk4oFPr6+vr7+4tEol/96lednZ2EkKSkJKFQqNFomP0qlUooFKakpExiLZqmr127tnLlShcXF+NMDzAR2BoF49NoNC+88MJ333337rvvvvzyy46OjmVlZSkpKfv27Ttx4oTRl+vv7+/p6fHw8DDcFhUVxXz5+PFjo08CsxNiFIwvNzf39u3bH3744bvvvquvLFu2LCIiwtnZ2RTLDQ8PE0JsbW0Nt7G2VbVarSmGgVkIMQrGl5+fLxKJ4uLimEV3d3cTLScQCGxtbWUymeG2jz/+mKKokZdKpXLGP/wZzAMxCkZG0/Tdu3fXrl1rY2Po3y6pVMrcfmQdKp0QKyurZ599trKycmho6H9ukwIYHU4xgZENDAxotVpXV1fDbX5+fksYJBIJl0V37Nghl8szMjK4fAjA5GBrFIxM/zxnpVJpuO327ds8Hm/kpUajCQwMnPSiMTExBQUFBw8ebGxs3Lp1q0gkamhoOHv27M6dO7ds2TL+z5FKpZGRkT4+PoWFhZMeBmYbxCgYGY/Hc3V1lUqlhtv8/f31gaunUqm4LGptbV1cXPz73/8+Ly/v9OnThJB58+YFBwfrT2rxeDw7OzvmgVFCCEVRdnZ2zCgnhOzevXv9+vVlZWVchoHZBjv1YHxhYWFVVVV1dXXmXHTOnDkffPCBTCbr6enp6uoaGBgoKSl5/vnnCSFHjhxRq9Wsfj6fr1arU1NTRyqnT5/WarW//OUvzTk2zACIUTC+hIQEHo/32muvtbW1jRSlUunDhw8n+lFFRUX/+te/JvQWR0fHBQsWTHShrq6upKSk3Nzcib4RADEKxhcQEJCXl1dbW+vt7b1x48Zt27YFBgZ6eXnl5eVN6HMKCgri4uL27t1rojmZdu/evWvXLi8vLzOsBTMMnlMPpiKTyfLz86uqqnQ6naura1BQUGRkpLW19fnz569cuZKZmWltbT3S/Pjx48TExNDQ0JFrjXp7e1esWJGenp6Zmfn999+bdNT6+no/Pz+xWGxjY6NWq7u7u1999dWJhj7MWohRmKLefPNNLy+viIiIN99809QxynTz5s133nnn9u3bZlsRpjvs1MNUdOXKlRs3bsTHx5t/aQcHBz8/P/OvC9MXtkZhylGpVP7+/r/+9a+DgoIaGxuzsrLOnj27YsUKS88FMDZsjcKU09PTEx4e3t7efvbs2UuXLvX29n799deWHgrgR2FrFKa0O3fumPnYKMBEIUZhSlMoFOXl5evXr7f0IAA/CjEKAMAJjo0CAHCCGIVZKisr6x//+AerWF5efujQIblcTggpLS09dOiQ/tb6I4aHhw8dOlRaWmq+QWHKQ4zCLJWenl5QUMAq3rx5Mzk5uaurixBy4cKF5ORk1iObtFptcnLyxYsXzTcoTHmIUQAAThCjAACcIEYBADjB3e9h9rpz505CQgKzUlFRwerZt28f805UrDNOAAQxChPS29vb29vr5ubGeiCHTCZzcHBwcnIaHh5ubW11cnLi8/nMBoVC0d/f7+bmZt55/4eurq5bt24xK+3t7ayeb7/9lvnH4nvWMAYaYNySk5MJIXK5nFWnKCo2Npam6ZaWFkJIRkYGqyEmJobH45lpyvERCATbt29nFXNycggh9fX1NE3r7y+lVquZDfpnRiUkJJhvUJjycGwUAIATxCgAACeIUQAATnCKCSaspqaGdQaJpaOj4+7du8yK/vJKgBkJd3iCCUhJSUlPTx/zV7Gxsbm5uTKZTCwWj9nA4/GGhoZMOd3ESKVSe3v7hQsXMotKpbKrq8vNzY3H48nl8p6eniVLljAbaJp+8ODB/PnzhUKheeeFqQtbozBhFRUVAoGAWfH29ma+TExMjImJYVaSkpJGX8BuWR4eHqOLAoFg5E8TCoWjs5KiKFawAiBGYcLc3d2dnJwMNAiFQlbWGD4IADCt4RQTAAAniFEAAE4QowAAnFinpqZaegaYTjw8PEJDQ21s2EfVQ0JCfHx8KIqys7MLCQlhna+nKMrX1zc4ONiMkwKYCb7wBADACXbqAQA4QYwCAHCCGAUA4AQxCgDACWIUAIATxCgAACf/B5gHSbliYJDiAAAA3npUWHRyZGtpdFBLTCByZGtpdCAyMDI0LjA5LjUAAHice79v7T0GIBAAYiYGCOCB0g2MbAweIAYLsTQjhGYkluZ2AlHis0BsKGbgcT7uZ+/QOXH/ffHpNuf3b7TPYGWz63o6bbfXDhH7t56L7Iuv7tkDVLMfqMa+hn/6nn37N+5XZWXbNz1oujVQzX6gmv2KT/bYvP2w3b62jPGAfQSbrZWcnMPW/Dy75WVse2Z6t9iXxAo4POTS2Mv/cfv+0jJGB5c8tr1GcnIHgGr2cdez2UzwbtlfECtwoItZw1YMAE76SC9JNje4AAABNHpUWHRNT0wgcmRraXQgMjAyNC4wOS41AAB4nH2SS26EMAxA95zCF0jkTz7OcgZGbVUNI7W0d+i+91edED6zaUAiOM9+DmSAOj6m959f2IdMwwCA/9ylFPgWRBzuUCdwvb28zTAul+sWGR9f8/IJxJbQrmf0sjzuW4RgBPSKObOCI0+SEclCSBpCOFLZQPKimlKyZRJURXAGYkgYD1BaxahYJFsGp6xKDWTBSAcYDHS7e1e77j61HSvZ5aEC1Z4LrHLSg0xrzW53u77bT2SGV3MGc5rJsa2zxlxB4WyfeAfVQPYhYyy2S8/CGnTduDLmAyytIopEq8Oe7YGtWwyEnA6Q0EjX5QG6OzbUqvNp50QV7XqpRPd3/alP+9utamsg1aTWgcDaAJ3Q2zw9HYL6vp0fmw9/fcR9zM8cv1sAAADGelRYdFNNSUxFUyByZGtpdCAyMDI0LjA5LjUAAHicdc47DgIxDATQq1CCFFv+JXaUkoY7IG6yhyfB2e3QSlt4Mk/zxOf1vV+vz9/f7bgTBrlLFGBk9UJIHGY2GDWirQMrxcxnQtaoDsIa1OdbRmm+I1GqPODiUoPNDUiurcvyvJf0OFZpg3CKG4y5wtTJCsi8SLSVqHjvQ9Cc6mJEJWwPDCGfJVKtrQhKrTnCmKStFcmlVn/RbMuAk4PL21xfrfRgg+lxfxxf3yhKo+cEQekAAAAASUVORK5CYII=\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x75ee4ffd40b0>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"mol2 = Chem.MolFromXYZBlock(xyz)\n",
"mol2"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "9a78be61-795c-4bce-af86-a68881b6c1ae",
"metadata": {
"execution": {
"iopub.execute_input": "2025-02-07T13:11:03.930292Z",
"iopub.status.busy": "2025-02-07T13:11:03.930151Z",
"iopub.status.idle": "2025-02-07T13:11:03.933898Z",
"shell.execute_reply": "2025-02-07T13:11:03.933649Z",
"shell.execute_reply.started": "2025-02-07T13:11:03.930283Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAY10lEQVR4nO3de1STR94H8AkIBFQoKIGAoqLFrS0qysWui6hUKCrlSBc94gUVVjkoKrYqx27XS/W0ikBRq3VbUVBZuqJ2rQpWlArrimBXUZBKBeQWkkBBbgmXXPaP533ZLFaseZIMid/PHxyeIcnzTY/n18nMZIajVCoJAACoy4h2AAAA/YYyCgDACsooAAArKKMAAKygjAIAsIIyCgDACsoovBLKy8t37NhRVlbWp/348eNffvkllUhgMFBG4ZXw+PHjnTt3PltGk5OTUUaBJZRRAABWUEYBAFhBGQUAYGUQ7QAAuvPVV19lZ2ertlRUVNja2tLKA4YBZRT6U1tbO2jQIHt7e9XGjo4OkUjk4ODA5XJpBVNPRUXF06dPVVs6OjpQRoEllFHoj6+v76hRo77//nvVxsuXLy9cuPDGjRszZsygFUw9n3766fz581VbvL2929raaOUBw4CxUQAAVlBGAQBYQRkFAGAFY6PwAm1tbXfv3lVtefLkCaUsAAMRyii8QH5+/pQpU2inYMvHx6e0tHT06NF92s+ePatQKGgkAsOBD/XwAt7e3k/+1+HDh2mHeml1dXWenp5//vOf+7TzeLw+y7kAXhZ6o/ACXC531KhRqi3Dhw+nFUZt8fHxbW1tfRaNAmgEeqNg+BoaGk6cOMHhcDZt2kQ7CxgglFEwfIcOHZJKpfPnz58wYQLtLGCAUEbBwEkkEmYwd/PmzbSzgGHC2Cj0Z+PGjVZWVn0a33rrrU8++eTZWe+BKTk5ubGx0cPDw9vbm3YWMEwcpVJJOwPopaysrK+//vqbb74xNjamneW55HL5+PHjy8vLz549GxwcTDsOGCZ8qAd1SCSSiIiIs2fPfvbZZ7Sz9CcjI6O8vNzZ2TkoKIh2FjBYKKOgDgsLi9OnTxsZGe3YsSM/P592nOdKTEwkhGzZsmUgd5lB36GMgpp8fHxiYmJkMtnSpUsH5l5z169fv337No/HW758Oe0sYMhQRkF9e/bsmTx5cnl5eUxMDO0svyIuLo4Qsm7dOnNzc9pZwJBhiglYKS0tdXd3l0gk6enpixYtoh3nv4qLiydOnGhubl5VVaWPX7sCPYLeKLDyxhtv7Nu3jxASGRlZXV1NO85/xcXFKZXK8PBw1FDQNvRGgS2lUhkUFPTdd9/NmDHj+vXrA2Eyp66uztnZWS6XP3r0aOzYsbTjgIFDbxTY4nA4x44ds7e3z83NZWbGqUtMTOzu7g4JCUENBR1AbxQ048qVKwEBAYMGDbp586aHhwfFJK2trU5OTi0tLXfu3Jk6dSrFJPCKQG8UNMPf33/t2rU9PT1hYWESiYRikiNHjrS0tPj6+qKGgm6gNwoa09nZ6eXldf/+/aioqC+++IJKhp6eHmdn59ra2qysLH9/fyoZ4FWD3ihoDJfLTUtL43K5hw8fvnDhApUMJ0+erK2tdXV19fPzoxIAXkHojYKGJSQkfPDBB7a2tkVFRXw+X6v3amhoEIlEQqGwvr5eLBbX1dWdOnWqoaEhNTV12bJlWr01QC+UUdAwpVI5b968zMxMPz+/rKwsDoej9kt1dnY2NTXV19cLBALmZ3Nzc+9lbW1td3d3n6cMHTq0q6vr8uXLvr6+7N4HwG+FMgqaJxAIJk2a1NjYmJSUtH79+uc9TCaTicVikUgkEAgaGhrq6urEYnF9fb1QKGS6lu3t7f3fyMbGxt7e3s7OzsHBgcfj8fn8/Pz8c+fOjRgxoqioyMbGRtPvDOBXoIyCVnz77bcLFiwwMzNLT0+3trZ+tkSKxWKxWNz/Pz8ul2tvb8/n83k8noODg52dXe+lo6Mjj8czMzPr8xSZTDZjxoxbt269//77GRkZ2nyLAP8HZRS0ZcWKFadOnZLL5f08xtrams/nW1tbOzg48Pl85mfvJZ/PV2NMoKKiws3NrbW1NSUlBXs7gQ7gEBHQlrFjx8rlcisrq4kTJ/L5fNVP30zXksfjsf/maEtLCzO/NHr0aCcnJ0KIs7NzUlLSypUr165d+/bbb7/++uuaeDcAz4XeKGhFV1eXs7OzQCC4cuUKm7VH3d3djY2NqjNLqhNNdXV1LS0tzCMPHDgQHR3d+8TFixenp6d7eHjcvHnTxMSE7fsBeD70RkErTp8+LRAIXF1d58yZ0/8jmUFSZtiUmW4Si8UCgUAkEolEol9++aX/pw8ePJgZJ7W1tVVtP3LkyL/+9a/CwsI9e/bs2LGD5dsB6Ad6o6B5SqXS1dW1pKTk5MmTS5cuZRqvXbtWWlraO7/Uu9izp6enn5cyMTFhpuDt7e17Z5ZUhwgsLCye99y8vLxZs2ZxOJzc3Ny3335bw28S4P+hjILmfffdd++9996IESMqKip6P1DPnz//0qVLzz6Yy+X+6vwS89POzo7N+GlsbOzevXvHjBlz7949S0tLtV8HoB/4UA+ax5zesWnTJtVByYCAgFGjRqkuWmKmm55dtKRBn3zySU5OTkFBwfr160+cOKG9G8GrDL1R0LDCwkJPT09LS8vq6morKyvaccjjx4/d3Nza29vT0tIWL15MOw4YIGxNAhq2d+9eQsjatWsHQg0lhIwbN27//v2EkKioqKqqKtpxwAChNwqaVFFR4eLiMmjQoMrKSm3vS/JSQkJCMjIyvL29c3JyBsIxJ2BI0BsFTYqLi5PL5UuXLh1QNZQQcuTIET6fn5eXxxzAB6BB6I2CxjBfJers7CwuLp4wYQLtOH1dvXrV39/f2Nj4n//8p5eXF+04YDjQGwWNOXTokFQqDQwMHIA1lBAyZ86cDRs2yGSypUuXtrW10Y4DhgO9UdAMiUQyatSoxsbGvLy8P/zhD7Tj/Lquri4vL6+ioqLVq1cfPXqUdhwwEOiNgmYkJyc3NjZ6eHgM2BpKCDEzM0tLSzM3N//rX/965swZ2nHAQKCMggbI5fLPP/+cEBIbG0s7ywtMmDDhs88+I4SsWbOmpqaGdhwwBCijoAEZGRnl5eXOzs5BQUG0s7xYdHT0vHnzmpubw8PDFQoF7Tig91BGQQMSExMJIVu2bNGLJZkcDic5OdnOzu7q1atJSUm044DewxQTsHX9+nVfX18ej/fkyRNzc3PacX6rzMzMefPmmZqa3r59e9KkSbTjgB5DbxTYYjYiiY6O1qMaSggJCAiIjIzs6uoKDQ2VSqVauktJSYlAIOjT2NjY+ODBg/53CAQ9gjIKrBQXF1+5csXCwiIyMpJ2lpeWkJDg6ur68OHDrVu3aukWbm5uO3fu7NOYnJw8ceJEoVCopZuCjqGMAitxcXFKpTI8PHz48OG0s7w0LpebkpJiamp66NChixcv0o4D+gplFNRXV1eXnp5ubGy8YcMG2lnU5Obmtnv3bqVSuWrVKnQPQT0oo6C+xMTE7u7ukJCQsWPH0s6ivg8++MDX17ehoWHlypWYcQU1YPd7UFNra+vXX39NCImJiaGdhRUjI6OUlJRJkyZlZWUdOXIkKipKs68vFArz8vJUW548eaLZWwBdWPAEatq7d29sbKyvr292djbtLBpw/vz54OBgLpdbUFDg6ur625/Y0dHBnGPae0gfc2afubn5mTNnTE1NnzcjX11dPXLkSA3FB5pQRkEdPT09zs7OtbW1WVlZ/v7+tONoRnh4eHJy8ltvvVVYWMjlcnvbm5ubBQJBfX29QCBobm5mfun92dzc/KuvZm1t3dTUZGpqGhoaGh8fr/qngwcP7ty5E2XUYOBDPajj1KlTtbW1rq6ufn5+tLNoTFJSUm5ubnFxsaenJ5/PZ7qWYrG4/66Gubm56oHPzFF9zMl9zAPMzMyGDRum+pR+DoUGfYQyCi9NqVQy3avNmzdzOBzacVj54YcfZsyYYWRkRAgZMmTI7t27V69eXVxc/ODBg97HWFtbP3vyc+9x0Hw+X9//IwBLKKPw0i5dulRSUjJixIhFixbRzsJKfn7+rFmzvLy88vPzmZasrKzW1tY5c+bExMQwXUs7OzumyAI8D8oovDTm258xMTGmpqa0s7DCnMvk6+vLXAoEgrS0NGNj48OHD48bN45qNNAnKKPwcgoLC3Nzcy0tLcPDw2lnYaWiouLChQtmZmbr1q1jWpKSkrq7uxcuXKjBGrp///4333yzT+Ps2bMTEhKsra01dRegTAkGoa2tLSEhobCwsE97dnZ2UlKSBm/0xz/+kRASGxurwdekYs2aNYSQiIgI5rK1tfW1114jhOTn59MNBnoHZdRAVFdXE0L27t3bp3316tUmJiaaukt5ebmxsbGZmVldXZ2mXpMKkUhkbm7O4XBKSkqYFmakYubMmXSDgT7C2Dm8gEKhqK+vLyoqyszMXLNmjVwuX7JkiYODA/NXpVL57bffKvVt9fHBgwelUul7773HHGLa09Nz8OBBQsjmzZtpRwP9g7FRIFKp9NlV5b1Lzaurq2UyGfNIDodjaWn5zjvv9D535cqVKSkpiYmJGzdupBT/pUkkki+//JKoFM20tLTq6uo33ngjICCAajTQSyijr4rCwkKhUCgSiQQCgVgsVv3+Yv+bFnM4HDs7Ox6Px6xIf/DgwV/+8pfAwMAhQ4YQQhYtWpSamhobGztr1ix92UP+2LFjjY2Nnp6e06dPZ1qYQ1C2bt2KFaCgBnwZ1EDU1NQ4OTm98847vaWBcfHixfv373d3d9vY2Dzva4tcLrfP2nLVy5EjR5qYmDCP7OrqmjZt2r1798LDw5l9SQghkZGRR48enTBhwp07dwb+BvhyuXz8+PHl5eXnzp1bsGABISQzM3Pu3LmOjo4VFRX6voQLqEAZNRBMGeXz+TweT7W9tra2tbWV2c5OKpXyeDxHR0ema9n7/cXBgwf/9huVlpa6u7tLJJL09HRm+b1EIpk6depPP/20YcMG5pjlgSw9PX3x4sUuLi6lpaXMuvrZs2fn5OTs27cPA6OgJrozXKApupmpZxw6dIgQ8tprr1VVVTEtP/74o6mpKYfDuXjxombvpXGenp6EkKNHjzKXhYWFhBBLS8unT5/SDQb6CzP18NKioqICAwOfPn26bNkyuVxOCJkyZcquXbuUSuWqVatEIhHtgM917dq1goICHo+3bNkypoVZ57RmzRorKyuq0UCPoYzCS+NwOMeOHbO3t8/NzWUmZwghmzdvnj17tlgsHsh7yDNFc/369cwYbmVl5blz50xMTKKjo2lHAz2GMgrqsLW1PXHiBIfD2bZtG/O52MjIKDU11cbGJjMz8+jRo7QD/ooHDx58//33qoeYxsfHy2SyJUuWYN9PYANlFNTk7++/du3anp6esLAwiURCCHF0dPzqq68IITExMcXFxbQD9sUcYhoREcHs/tnU1MT8n0CPVrzCAEV5bBY0pLu7+/79+w0NDX3a6+rqiouLtXRTqVQ6ceJEQkhUVFRv44oVKwghbm5uXV1dWrqvGmpqakxNTY2NjR8/fsy0MMfHz507l24wMAAoo8BKcXExc97GP/7xD6alvb3dxcWFELJlyxa62VRt2rSJELJ48WLmUiqVMrvT5+TkUM0FhgBlFNhidsK3tbUVCARMS2FhoYmJiZGRUXZ2Nt1sjJaWFmYi/scff2RaDh8+TAhxd3enGwwMA8ZGga2YmJiAgICGhoYVK1YolUpCiLu7+0cffaRQKMLCwn755RfaAcnhw4dbWlp8fX2nTJlCCFEoFMwCgy1bttCOBgaBdh0HQ1BXVzd8+HBCyIEDB5gWuVzu4+NDCAkODqabrbOzk9mPKisri2k5c+YMIWTMmDE9PT10s4FhQBkFzTh//jwhxMzMrKioiGmpqKiwtLQkhCQnJ1MMduzYMUKIq6urQqFgWqZNm0YI+eKLLyimAkOCMgoa86c//YkQ8uabb0okEqYlNTWVEDJ48OBHjx5RiaRQKJgzPE6ePMm03LhxgxBiY2PT3t5OJRIYHpRR0Jj29vbx48cTQmJiYnobQ0NDCSHu7u7d3d26j3ThwgVCyIgRI3pXX82fP58Qsn37dt2HAUOFMgqa1LtHyaVLl5iWp0+fjho1ihDy8ccf6z6Pt7c3ISQhIYG5ZHZ1srCwEIvFug8DhgplFDRsz549hBAejycUCpmWvLw8Y2NjIyOjH374QZdJCgoKyP/u3rRy5Uryv18WAGAPZRQ0TC6Xz5o1ixASFBTU27ht2zZCyMiRI5uamnSW5P333ycqh5gKhUIul2tsbPzzzz/rLAO8ClBGQfNqampsbGyIyraePT09Xl5ehJClS5fqJsOzh5hu3bqVEBISEqKbAPDqQBkFrcjIyCCEWFhYlJaWMi2PHz8eOnQoIeT06dM6CMBs4xQeHs5c9h5Df+vWLR3cHV4pKKOgLcuXLyeETJkypXeWnNn/ycrKqrKyUqu3bmpq6nMM/f79+wkhPj4+Wr0vvJpQRkFb2traXn/9ddXRSaVSuXDhQkLI9OnTZTKZVu9+69atTz/9lPm9u7vbycmJEDLwzzgBfYQj7UCLCgsLp0+fLpfLr169Onv2bEJIc3PzpEmTampqdu/e/dFHH+kmRmpqalhY2O9+97uSkhLmGDsADcI/KdAiDw+Pbdu2KRSK5cuXNzU1EUKsra2Tk5ONjIx27txZVlammxi9x9CjhoI2oDcK2iWTyWbMmHHr1q3g4OCzZ88yjdu3b3d2dg4LC9NBgKysrICAAAcHh8rKShxDD9qAMgpaV1lZOXny5NbW1hMnTmipdMrlcrFYLBKJBAKBWCwWCAQikUgoFNbX19+/f7+lpWXXrl0ff/yxNm4NMIh2ADB8Y8aMOXDgwIoVK9atW/f73/+emXd6WVKptL6+XiAQNDc3M7+oXlZXV8tksuc919jYuLy8nMU7AOgPeqOgI6GhoX/72988PDxu3rxpYmLS569dXV1isbiurk4sFtfX1wuFwt6uJXPZ2dnZz4tzOBwej8fj8RwdHXk8nr29PZ/PZy5bW1sXLVoklUr//ve/h4SEaPMtwisKZRR05OnTp5MnT66qqgoODp42bVp9fb3qp29mAqofQ4YMebZE2traOjg42NnZ8Xi8QYOe+9EqKSlp48aN1tbWRUVFOEsZNA5lFHTn2rVrc+fONTMza2tre/av1tbWfD7fwcGh96e1tTXzi6OjI/MdJPUolcrAwMBLly7NmTPnypUrHA6HxZsA6Atjo6A7AoGgu7vbysoqIiKC6Vry+Xx7e3s7OztbW1vt3ZfD4SQnJ0+cOPHq1auff/55TEyM9u4FryD0RkF33Nzc7t27d/z4ceYsex3LzMycN2+eqanp7du3J02apPsAYKiwGhl0JDMz8969e46Ojsx++LoXEBAQGRnZ1dUVGhoqlUqpZACDhDIKOhIXF0cI2bBhA8U18AkJCa6urg8fPmQ2zQPQCHyoB124c+eOh4eHpaVlVVUVm8ki9u7evTtt2rSenp4LFy4w5zIBsITeKOgC0xVds2YN3RpKCHFzc9u9e7dSqVy1apVQKKQbBgwDeqOgdZWVlS4uLhwOp7y8fCAs21QoFH5+fteuXXv33XcvX76M9U/AEnqjoHXx8fEymWzJkiUDoYYSQoyMjFJSUoYNG5aVlXXkyBHacUDvoTcK2tXU1OTk5CSRSO7evTuglhmdP38+ODiYy+UWFBS4urrSjgN6DL1R0K6DBw92dHQEBAQMqBpKCFmwYMGqVas6OztDQ0P7/8I+QP/QGwUt6uzsHD16tEgkysnJmTlzJu04fXV0dEyZMqWsrOzDDz9kJsEA1IDeKGjR8ePHRSKRu7v7AKyhhJDBgwenpaWZmJgkJCRkZ2fTjgP6CmUUtEUulzOnd2zZsoV2lueaOnXq9u3bFQrFihUrJBIJ7Tigl/ChHrQlIyMjJCRkzJgxZWVl/exiR51cLg8JCQkLCwsKCqKdBfQSeqOgLfHx8YSQDz/8cCDXUEKIsbHxuXPnmGX5paWlff6alpaGRVHQP/RGQStu3Lgxc+bMYcOGVVVVDR48mHacF8vNzfXx8Xl2h3w/P7+qqqpHjx7RCgYDH3qjoBXMxHd0dLRe1FAANlBGQfN++umnzMxMCwuLqKgo2lkAtA5lFDRv3759CoVi5cqVWt3THmCAGNBj/6CnbGxshg4dqo9ndZw8efL27duqLWVlZWZmZrTygF5AGQVWhEKhkZERj8dTbdy9ezdz2hKtVGr7+eef+5xR2tLS0ufdAfSBMgqsvPvuu0OHDs3Ly1NtzM7ODgwMzMrK8vf3pxVMPbt27frVmXpaeUAvYGwUAIAVlFEAAFZQRgEAWMHYKLAlkUju37+v2vLkyRNKWQAoQBkFtv79738PtC2Z1eDp6VleXm5nZ9enPTU1VSaTUYkE+gJlFNhyd3f/5ptvVFtycnIiIiJo5VEPl8t1dnZ+tt3e3l73YUC/oIwCW88WoIcPH9IKA6B7mGICAGAFZRQAgBWUUQAAVjA2CqxERkY+u3OHi4vL9u3bx40bRyUSgI5h93sAAFbwoR4AgBWUUQAAVlBGAQBYQRkFAGAFZRQAgBWUUQAAVv4DQ9EKA/Ra1MsAAAEQelRYdHJka2l0UEtMIHJka2l0IDIwMjQuMDkuNQAAeJx7v2/tPQYgEABiJgYI4IHiBka2BA0gzcxCLM2oAKIZGYmluRkZMpiYGBOYmTKYWJgTWFkymFgZEtgYGNgZGTiYGDiZGbhYGLhZGZwYgcrZGFhZmJkYxWeBtEIxA4/zcT97h86J+++LT7c5v3+jfQYrm13X02m7vXaI2L/1XGRffHXPHqCa/UA19jX80/fs279xvyor277pQdOtgWr2A9XsV3yyx+bth+32tWWMB+wj2Gyt5OQctubn2S0vY9sz07vFviRWwOEhl8Ze/o/b95eWMTq45LHtNZKTOwBUs4+7ns1mgnfL/oJYgQNdzBq2YgCx2EynvfREEgAAAXh6VFh0TU9MIHJka2l0IDIwMjQuMDkuNQAAeJx9k0tygzAMhvecQheAkWTZFssQMm2nEzLTpr1D973/VDIOdjYF2xjzWb8eZgC/Ptb3n184rrAOAwD+0+Z5hu+AiMMVfALL5eVtg/P9tDxWzrev7f4JxN7Q72f2dL9dHysEZ8BJMWdWGGmikBHJlpBURNpWNpCmoJpSss8UUBVhNBAlYWxgKBaj4hyy7eCUVamAHDBSA8XA8dA+pMeq3bkdnazi4oCr5xl2cdJGpt1mVR8P+arekRleTVNM05RGtu+sMTsYOFuOD1AN5EkyxtminDiwiu6BK2Nu4FwsYgjR7PDE9sDiLQohpwYSGjlWcYGqHQtq1rmLnMjRKh+cqPpVvvPTal2sFgeSbyoeBNgdoA69bOvTIdiPxXLb1nYszC+gVnyxunIrsfhCK6RYiaRVS6wOsZXEX6nlnby37JJ1bjkk66HLFPkgXT7Ih9hFTT6kPrY+En9//Bw2H/4AL0qmL0LmsRoAAADNelRYdFNNSUxFUyByZGtpdCAyMDI0LjA5LjUAAHicbc/JDcQgDAXQVuaYkQB5Axu5gfQQ5ZQ2UvywJae5WIhv7Mexnxde27Gf3z8VW/3cGyZgziVQopw1REggCFQcUjaorAETFTUbETFk9HYyUCVrGfJ8hCbiEZOwgrRpQGR5JExGPWGzUvoNMpjW0CKQguaRkijk2jNq3WIjM4LaVy1GfByLYX3mtMeFn3as/vriBC6f+OOLE1jC8GmtPn29tfPmbxsPsj+6FxeXTr/3D/msSV97bbvnAAAAAElFTkSuQmCC\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x75ee4ffd40b0>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"rdDetermineBonds.DetermineConnectivity(mol2)\n",
"rdDetermineBonds.DetermineBondOrders(mol2)\n",
"mol2"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "7a583769-d831-4cef-8741-98745b832c48",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 5,
"id": "b0244c7d-7d1a-4bc3-83a7-f8f9cf3c673d",
"metadata": {
"execution": {
"iopub.execute_input": "2025-02-07T13:11:03.934326Z",
"iopub.status.busy": "2025-02-07T13:11:03.934230Z",
"iopub.status.idle": "2025-02-07T13:11:03.959338Z",
"shell.execute_reply": "2025-02-07T13:11:03.958926Z",
"shell.execute_reply.started": "2025-02-07T13:11:03.934317Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVhU1fsA8Hdm2BcBWUQUSgXF3YRMQ0UNUwzLDc09cS23r6ZilrnkApqJ4IbggigqKVlmP/eFJTNHXBBBZRBlX2dgYBiWuef3x9GJSFHh3rkDvp/Hp+fecTjnnYSXc+499z0CQggghBCqLyHfASCEUOOGaRQhhBoE0yhCCDUIplGEEGoQTKMIIdQgmEYRQqhBMI2iF8vJAX9/yM19dpqUBFFRvAaEkLbCNIpeLDMTVqyAxYufnd65A+HhvAaEkLbCNIpeqksXSEiACxf4jgMh7YZpFL2USARbt8LcuaBU8h0KQloM0yiqy0cfQY8e4O/PdxwIaTFMo+gVfvoJdu6Ep0/5jgMhbYVpFL1Cq1bg6wubN/MdB0LaCtMoerUFC6BVq2fHjx7BkiVQUsJrQAhpEx2+A0BayskJdu16dqyjA7/8AnfuAADMmQOXLoGDAyxYwGN0CGkRAdYbRa8kl4OrK2RlQUYGXLoEo0aBoyM8eABCnMwghJN6VIf0dJg3D5YuBVNTsLeH0lLYvx8+/RTatIGUFDh3ju/4ENIOmEbRSymVsHMn7NwJRUUwfz4AwPbtIBDA7NkAAEFB/EaHkLbANIpeyskJPDxAoYADB2D4cGjTBiQSOHMGZs0CIyP4v/+Dhw/5DhEhLYBpFNWFDkKDgoAQ+PLLZ8cWFvD550DIP/egEHqb4S0mVBeGAScnSE2FU6fAzQ3s7UGhgKQkUCigZ09o1051+3aFiYkR32EixCccjaK6CIXw1VcAzweh48cDIbBzJ7z3Hkyd+nNGRvPw8DC+Y0SIZ5hG0StMnw7GxnD+PCQnw7x54OaW9tdf8+Ry+bBhpKKiJCgoCCc06C2HaRS9grk5TJwINjaqyMi/u3cHkWjq33/vOHjw4KhRo1q3bp2UlHTp0iW+Y0SIT5hG0astWJAllZr9+KNHSUnJ/PnzAWDbtm0ikWj27NkAEIRLn9DbDW8xodcycODAK1euBAYGfvnll23btk1PTz937lyPHj0cHBwqKytTUlLatGnDd4wI8QNHo+i10EFoUFCQUCicM2cOPba2tvb29mYYZhcufUJvMRyNoteiUqkcHR3T0tLOnDnj4uJib29fWVn56NEjmUzm4uJiYWGRnp5ubGzMd5gI8UC0evVqvmNAjYBQKFQqlRcvXpRKpTNmzEhJSbl9+7auru7UqVMfPnw4YcIEV1dXXV1dvsNEiAc4GkWvq6CgwMHBQalUJicnl5SU9OrVy9LS8t69ey1atOA7NIT4hNdG0euysrL6/PPPCSEhISGurq6WlpYFBQXffPMN33EhxDMcjaI3kJCQEB0dPWXKFFNT08DAwIULF9ra2j5+/NjAwIDv0BDiDaZRVH8uLi7x8fHBwcGzZs3iOxaEeIOTelR/y5YtAwB/f//q6mq+Y0GIN5hGUf2NGTOmffv2qampP//8M18xTJ069Y8//lCfzpgx49dff+UrmKZn3bp1169fV5/6+fnFxcXxGI92wjSK6k8kEtEB6bp16xiG4SWG3NzcsrIy9WleXl5paSkvkTRJ0dHRGRkZ6tPY2NinT5/yGI92wjSKGmTKlCkODg73798/ffo037EgxA/cYBk1iK6u7qJFixYtWrR+/frhw4fzEsO+ffuuXr1KjxMSEsaNG8dLGI1AVhaEh0NSEtjawrBh0L9/He9lGKaqqgoArly5oh7g1xyZIjUcjaKGmjVrlrW19fXr169cucJLAN27dx/ynI2NDS8xNAJPn8L774NMBlOmgLMzTJgAe/f+6w0ZGY8vXfrxxx+nTZvm6upqamoaHBwMALm5uZLnFAoFP8FrNxyNooYyMjJasGDBypUrN27cOGDAAM0H4OLioh4Ih4SEaD6AxsHfHz77DDZufHbarh2MHAmjRsGKFZCYCPfugVT6yMVl6c2b6q+gY89x48aNHj2avhIfH6/xuBsBHI0iFsyfP9/MzOzcuXM3btzgOxb0EvHx4On5z2m/flBZCQUFEBYGMTEglYKVVQ97+/nz5wcHB8fFxclksk2bNvEXbmOCo1HEAjMzszlz5vj7+/v5+Z04cUKTXVtbWxsZGdU8xUJTLyaTQbNm/3rF3ByKi2HPHrC1ha5doUULG4BAnqJr1PApJsSO3NzcNm3aKJXKu3fvdunSRTOd/vDDD87Ozt7e3prprnHz9IRx4+CLL56dKhRgZgY5OWBpWccXSaVSIyMjfX19eiqTyQwMDPDZ31pwUo/Y0aJFCx8fH0LI5s2bNdPjrVu31qxZM2HCBIlEopkeG7fhw2H7dlDfI/rpJ+jdu+4cCgAWFhbqHAoA5ubmmEP/C0ejiDXp6emOjo6EkIcPH7777ruc9sUwTN++fa9du7Z48eItW7Zw2lcToVLBzJlw9Sr06wdpaZCXB6dOQbt2fIfVFGAaRWyaOnXqwYMH586du337dk47Cg0NnTlzpq2tbXJyspmZGad9NSk5OZCcDLa24OQEIhHf0TQRmEYRm5KTkzt37qynp5eamtqyZUuOeikqKnJ2ds7Pzz969Cgutn8zU6bAuXOwbx8MG8Z3KE0HXhtFbHJ2dh4xYoRSqdy2bRt3vSxfvjw/P9/d3X3s2LHc9dI05eRAbi7gdi+swjSqUTExMQUFBerTuLi4/Px8HuPhwnfffScQCHbu3CmVSrloXywW7927V09Pb9euXQKBgIsumjL6j2JhwXccTQqmUY1avHjxnTt31KfLli0Ti8U8xsOF9957b/DgwXK5fMeOHaw3zjDM3LlzGYb5+uuvO3bsyHr7TR+mUQ5gGtWczZs3V1RUqE9/+umn8vJyHuPhDt2gKSAggPWadbt27fr777/t7e2//fZbdlt+W8hkAADm5nzH0aRgGtWcw4cPuxcWOoaGwvLlsHz5X/v3V1ZW8h0UJwYMGODm5lZYWBgaGspis3l5eStXrgSAwMBALX9UiZCKiorH1dWFfAfyb4RAcTEIBIBrG1iFD4NqzkfFxTMLC985epSeBgkEjjpN9v+/r6/vp59+umHDhhYtWlhZWZmamjZr1szU1NTU1NS8vkOhpUuXSqXSIUOGjBgxgt1o2VVScjYtzUdf31GlKmrefJKtrS/fET0nl0N1NTRrBk33G48X+H9Tcz6Uy50rKh526SLs3bttWFiLqqrOjo58B8UVT09PKysrAwODCRMm/PdvDQwMLCwsDA0N6QFV67TWKw8ePAgPD9fX1w8M1PbHvjMyfB0cdpibjwAAQrRplyq8MMoNTKOadhbgDiFBAIZ8R8KprVu3FhQUmJubjxw5Ui6Xl5SU1PyvUqnMzs6uR7PLli1r374969Gyi5DK8vJEmkYFAm36EcM0yg1t+jd+O8xJT9dVKECl4jsQDj148GDVqlUAEBER4VmzONtzMplMnVXlcrlMJispKVGflpSU0Ffkz8lksuLiYoZhrKysNP5p3piDw/a0tKlFRYesrHxsbBYKBHp8R/RMdXGxjp4eXhhlH0GaMuPdd6t0dAgA/XPJwuIbX9/U1FS+42KZSqXq27cvAMyYMYMQkpWVxUqztP6ejY1NSUkJKw1yTFVScjEp6YO0tOl8R/KPX375BQAmenvzHUhTg3fqNedoQcHANm2eXLoEEglIJCvatRvk4dGmTRu+42LZli1bYmNj7ezsNm3adP369TZt2tDb6w00atSoPn365OXlBQUFNbw17glNTQe1auUnl18hpFoqPcZ3PAAARUVFAKBvasp3IE0NplEN2bBhQ2lpaWJBgUm3btC2LbRte+DQoQ8++IDvuFimns6HhoYaGRlNnz69oqKiurr+t1kIIXv37u3WrZtUKvXz8wOATZs2FRZq2UKif8vN3VpWdkOpTCoo2GNs/EFmpm9q6ucSyUiVqpjfwOhzZfVeKYFeBtOoJkilUppc5s6da/m8wmOHDh1Mm9a4gGGYGTNmlJeX+/j4eHp6fv/994mJic7OzvSz149AIDh27FhCQoKfn1///v2HDBlSXFzs7+/PYthsYRhFZuY3KlWxQCDKydn49Ol8PT0HB4cdJib9RSJzmexkcnKv8vIEHiOUyWQAYIG3mFjH91WFt8KgQYMAwMDAoKKigu9YOEQLNtvZ2RUVFf31118ikUgoFMbGxjaw2du3bwuFQgMDg6dPn965c0d9zErMLEpNnSgWQ0rKqP/+lVL5MDGxm1gM8fEG+fmhmo+NmjdvHgAEBQXxFUBThaPRF5DJZGlpaerT4uLix48f17s1sVh8+fJlAFi/fr2enrbctGXdw4cPv//+e6gxnVepVEuXLnVzc2tgy927dx8zZoxSqfzhhx+6des2duxYpVK5du1aNqJmTV5eQFHRYaHQxM7uBYHp6zs5O1+3sprJMMonT2akpU1hGB52KqaTehyNso/vPK6NwsPDhw8frj49duzYkCFD6t0arQNvb2/PRmhaSqVS9evXDwCmTZtGCPH19QUAZ2dnhULBSvsPHz7U1dUViUT3799PTU3V09MTiUSJiYmsNN5wcnnczZt6YrGgqCiy7ncWFITFxxuJxXD//ntKZYpmwlMbNmwYAJw+fVrD/TZ5OBqtLTExseYuwUlJSdeuXat3ayEhIXRgGx4e3vDYtFZAQEBMTIydnd2WLVvi4+N/+uknoVAYGhpqaMjOQwZOTk4+Pj4qler7779v06YNHerSwS/vqqpyUlO9Cam0tfW1sHjF5nqWllM6dIjT13dUKG4lJfWUSjW6iyq9Nnrq1KmG3PRDL8B3Htc6Bw8edO/efXu3biQ4mAQHX/H17dGjR/1Go9XV1bSCRr9+/ViPU3tIJBL6MU+ePKlUKum2oMuWLWO3l6ysLCMjI4FAcO3ataysLNrjn3/+yW4vb4phKpOT+4nF8ODBRwxT/ZpfVV1dLJF4i8UgFguePl3AMJWcBqnm+Pzh4549e96+fVsznb4NMI3WdnT37l0tW6oXyTMA43V0hgwZ8vDhw//7v/8rLCx8/aamTZsGACKRKDs7m7uA+aWezn/xxReEEFoir0OHDmxN52tavnw5ALi7uxNCVqxYAQD9+/dnvZc38uTJl2Ix3L3rUFWV94ZfyuTk/Hjzpo5YDMnJg3Jy2HlI4WXEYrGXl5dAINDT06NXmXR0dHx9fZVKJaf9viUwjdb2q78/AZDr6BSOGaPo1IkA7La1HTJkiHrVTsuWLb29vek0try8/GXtZGdni0QiAPjqq680Gb+G/fTTT/T/SVFRUXx8vK6urlAojImJ4aIvqVTavHlzADh37pxMJqNLx86ePctFX6+jsDCc3nwvK7tRvxbk8pg7d+xCQtytra3PnTvHbniEkISEhOXLlzs4OKhnnwMHDlQoFL6+vvSb09HR8erVq6z3+7bBNFobTaOPdHQ++eSTw05ONI26uLgMGjSoQ4cOtTbpNjAw+PDDDxctWnTkyJHHjx/XbKd3794AYGxsXF39unM9DlVVkbw3HS69mkQiMTExAYBffvlFPZ1funQp6x2p0RX4Li4uDMNs2rQJALp3765Sqbjr8WXKym7Rm0X5+Xsb0k55ea6Hx0d01rJx40aGYRoeW0ZGRkBAQM01ElZWVkKhEADatWt36dIlQkhcXJyzszPtd9GiReVlZQ3v962FabQ2mkbLRSLi4UHatqVptObvcyMjo3bt2nXs2LFVq1a19gJq0aKFz+DBe7/66tCiRS4AxgD79u3j+fMUFpIRI4iDA3FxIa1bkw0b2GpYpVL1798fAKZOnUqez7I5ms6rlZeX29vbA0BkZKT6+OjRo9z1+EJVVYUJCW3EYnjy5MuGt8YwjJ+fH01zn3zyyRtdOKpJJpOFhYV5eXmJnu+cbG5uPnny5PPnzzMMc+fOHVdXVwAQCASTJ08uLCysrKz08/PT09Pzbt+ecXIi5883/LO8nTCN1nZkz54UExP1tdFqPb25Xbv6+vr6+/uPGjWqVatWte7R2dnZ9ezZs0ePHhYWFvYAx55/IQG4rq9PKjV09+ClRowgPj6kqooQQtLTiaMjOXyYlYa3bt1Kp/OFhYVcT+drCg4OBgAnJ6fKysqQkBD1Mdf9qlVXV//88wyxGJKSejMMa89TXLhwwcbGBgDs7e2vXbv2+l+oVCp/++23yZMnGxkZqSdJXl5ekZGRtR73qKqq8vPz09fXp/9wUVFRhJD4+Hjp8OEEgAgEZNYsIpOx9YneHphGazt48ODEYcPI+fP0z+kff/T09Kz5hoyMjKioKF9f3wEDBtAprdoOoZAAZOnoiAGqaDLNz+frgxBCiFxORCJSUPDPK8HB5N8fp35SU1PpZ4+KiqqoqOjatSsALFmypOEtv1J1dTXdzC44OFh9vGvXLg10TS1fvlwgEKxY4VVZmcFuy0+fPu3Tpw8A6OvrBwQE1P1mlUoVExOzYMECdfFAoVDo5uYWHBxcXFxcxxcmJibSXgBg0qRJsvx8UlVF/PyIvj4BIC1bkqgoVj9W04dptLbMzMybN2+qT7Oysm7ceOkNhOrq6nv37oWFhS1YsMDNzW2nSEQA5gEYGBjkAxCA1L//1kjUL3H3LrGx+dcr0dGkfXvy4AFZtIhs20Z++43cvUvk8jdqVaVSubu7A8DkyZMJIXR3ufbt23M6na8pMjKSjqfKysqOHz9Oj0tLSzXQ9cmTJwUCgY6OzpUrV7hov6qqij68AAATJkyQv+if5t69e6tWrapZG6xTp05+fn6vX5OQYZjg4GATE5Ohbdow9vYkLIwQQh49Iu7uz+ZS3t5cXExvqjCNsqly9mwCcNPQMBigHIAAjPvoIz4DSkkhpqb/euXMGdKzJzl6lNS4+EAAiJUVcXUlY8aQJUvIjh1/nD6dmJj4srQYEBCgns7funWLTuejo6M18YkIIYQwDEPv4Pn5+ZHnd/M2sHfZ92WSk5PNzMwA4JVDxQaKioqiHXXo0CEhIYG+mJ6eHhAQ0LNnT3X2dHBw8PX1TUpKql8vEokkb/LkZ98AI0aQrCzCMCQ4mNCLWjY2z9IrehVMo6zasaNmbkoDMBYKMzJYnvq9gaoqYmpKai60/vZbMmUKSUoimzaRr74inp6kY0diYFAzbKZGITVbW9vevXt//vnny5cv371799mzZy9evPjf6fzixYs1/MmuXLlCb6EUFhbWPOauR7lc3qlTJwAYP348d72oJSUlde7cGQBMTExmzJjRr18/9f1MKyuruXPnxsXFsXBbn2FIWBhp3pwAEDMzEhxMGIakphIPDwJApk1j46M0fZhG2XbiBAkO/tbaehbAOwDAe0GdH34gPXqQ69dJbi4JDydWVuTOnRe8LTubXLtGIiLIhg2yJUs+/vjj9u3b03sR/yUSiUaPHk0IOXbsGL3DU8bHchkPDw8A8PX1JYQMHjwYOHh0So1hmDFjxgBA165dNXP1gBCiUCimT58OANbW1jVvHLF/Py0ri4wc+ez3qLs7efSIMAzZv58UFZGUFDJvHnF3JyNGkK1bn92rRP+GaZQTO3bsUCedLVu20BdVKtWdF6YwroWEkOHDiasrmTyZiMWv/3UMw2RkZMTExISHh69du9bHx2fQoEG0PtCsWbPoe06cOBEXF8dN3K8gFosFAoGBgUF6enrNYnpc9LVhwwYAsLCwePToERft1/Lw4UN6IJPJ6DWTHTt2vPA6KZsOHyZWVgSAmJgQ+rs/M5PY2pJNm8jjx+TGDeLuTnx8uI2hccI0ygn63U/TaOfOndUvgiaLGERHk88/J0VF7LaakJAgFAr19PS0YRcpb29vAJgzZw4hZOzYsQAwbNiw2NjYx48fs/iY44ULF2jt1N9//52tNuuQk5MjEom6dOmiUqkOHjwIAIMGDdJAv4QQUlREZs0iAGTsWEIIWbmSTJxYMzJiZETS0zUUTOOBaZQro0ePpguqAeD69etEw2lUKiXvvEMAyMaNrLc9ceJEAJg+nf/N2h48eDBlyhSa0KOjow0NDWnJEsrAwKBt27Zubm7e3t6+vr4BAQGRkZExMTESieT158VPnjyhK4rWrFnD5Uf5B53KfPbZZ4SQzz77DAB27Nihma6fOXmS5OYSQsjIkWT37n/9VY8ehL+nb7UWplGunDp1iq4EAoDZs2f/+eefKSkpmkujEyYQAOLiwsX6f3XFz/v377PeeP0wDEMvj7Zu3bp3796tW7dWzwZeSCgU2travvfee15eXjNnzly1atWePXtOnToVHx+fnZ2tfrq0vLycPvnj5eWlsUdO6V4J4eHhcrnc0NBQKBRmZmZqpuvaPD3J/v3/eqVPH3LyJD/BaDFMo1zp1avXkSNH6A+tQCDQ19engxpN9B0eTgCIsTF58ICjHmbPng0A48aN46j9NxUYGEhvYefk5NBXGIbJzs6+devW77//HhISsnr16tmzZw8fPtzFxaVly5bqicIL6ejotGrV6oMPPmjXrh0AODo6SqVSzXyQ/Px8HR0dPT29oqKiiIgIeF7Rih/z5hFf339OKyuJlRXRmmrZ2gPTKFd69er1888/6+jo0J/M9u3b0/vanHf89CmxsCAAZG+DSmbULTMz09DQUCAQxMfHc9fLa0pJSaFrsH7++efX/6qioqJ79+6dP38+LCzMz89vwYIF3t7etZKsvr6+oaHhiRMnCCEVFRVBQUH37t3j7HMQ8vxRVy8vL0LIqFGjACAwMJDTHusSG0usrUly8rPTH34gfFcm1E6YRrlC06ienh6t7WZsbKyJNKpSPXsQZeRIbjsiZPHixQDwySefcN1R3aqqqug+1bTgKSuUSmVaWlpcXBy9NEkvU9JLNCNGjGCrlxeilyb2799fVlZmbGwsEAh43rwvNJTY25PevUmHDuTDD4kW3FfUQphGudKrVy9ra2uBQKDeQWzbtm2cp9E1awgAadXqX8/RcyM/P79Zs2YAwNFjka9pzZo19JJoEdtrEggheXl5dBPsP//8s+Yx6x1RBQUFOjo6urq6hYWF9Jeum5sbR329maysN31i+K2CaZQrvXr1+vXXXxcuXFhRUUFvH3/44Ydz587lsMsbN4iuLhEKyYULHPZSA90NqW/fvprp7r/i4+P19PSEQuHFixc56oIOQmmZfVrbn7uS+6GhoQAwdOhQ8nz91tatWznqC7EI0yhXevXqpa5nvnDhQnp3ODY2lqMqziUlJatHjiw3MCBcVk3+b6f0ARteStArlUr6KOr//vc/7nqRyWTqkvvq8vsXuPlFNXToUAAIDQ0tKyszMTERCARPnjzhoiPELkyjXKmZRlNTUwUCAd0m09LS0tvbOzg4OJ3VZcxffPEFAEzw8iIVrFXAfB3+/v7wvBy9Jvslzy/OsriN88vQkvuurq4Mw2zcuFF9zG4vUqlUT09PR0cnLy+Plq3q3bs3u10gjmAa5cqBAwfU24rQOxWtW7d2cnJSr6oRCAQ9evRYvnz55cuXG/iU9IkTJ+hq87t377IQ+psoLy9v3bo1ANDb2RoTExMjEol0dHToow2cUigUtFz3iRMn1MdRbBfl3L9/PwAMHjyYEDJ+/HgA2Lx5M7tdII5gGuUc/fGwtLSk+4NKJJLg4GBvb296f4YyMjLy8PAICAiotaHT68jIyKAzzZ07d7If/WugT9106NChSlN1K+RyOV3RuWrVKs30WPMzbt++HQCGDl3C7uUZLy8vAAgODlYqlfR7QyKRsNkB4gymUW5lZWXRO/WHDh2q9VdVVVUxMTG+vr4uLi4193Rq27btrFmzIiMjS0pKXtm+SqWiD714enpqflpNVVZW0qQWpqnylD4+PgDQs2dPje0dUllZ6eTk3LdvyKFDZRUVFZ9++kAoJAcOsNZ+SUmJgYGBSCTKzc09efIkALz//vustY44hmmUWyNGjFCvpq5Dbm5uZGTk5MmT1auj6CTdw8PDz89P/PKyTPRSnY2NjfrpHV6EhYUBwLvvvlvB/ZXZ3377jS6MV9cz1owjRyoByLvvEqWShIURAPLOO4St+ifh4eHwvATJpEmT4HlFatQoYBrlEC3PY2Zm9vp3k6qrq8VisZ+fn4eHh/oJKACwtbWdPHlyZGRkzacSb968qaenJxAINFN5qA7V1dW0pDHXRTTy8/NtbW2B+/rz/6VSke7dCQAJDPznmK1asuoSJBUVFfRXqWYq8iFWYBrlSnZ2Nr1kWe+pbkFBQWRk5KxZs2puRyoSiVxcXFatWhUbG0v3GV+4cCG7kdePek8kTks408p4/fr142Vv+pMnCQCxtiYlJeSXX/45bqCaJUhOnToFAO+99x4b8SINwTTKFfpA9LBhwxreFMMwYrF4/fr1/fr1qzlEFQgEzs7O5eXlDe+i4RiGoQ9l+vv7c9QFvXRgZmaWlpbGURev1KcPASDr1v1z3PAtoGgJErqqny5cW79+fYMjRZojIIQAYltERMTEiRPNzMwSEhLs7e1ZbLm4uPjixYtnzpw5dOhQeXn5oEGDLl68yGL7DXH27NmhQ4daWFikpqaam5sDgKen5+PHj01MTMzNzY2MjIyNjZs1a2ZqampsbGxsbGxubk4Par2hWbNmIpGoVuOZmZldu3aVSqUHDhyYOnUqH58PACA6GtzdwcwMUlPh3j1wdwdzc5BIoHnz+rc5ZsyYEydOBAYGzpkzx9bWtqioKCkpiU41UOPAdx5vgvLz821sbABg37593PXyv//9DwCMjIx4md6+zMCBA6HGOiR6B78e9PX1mzdv7u3tTdthGIY+4UOrhPBr8GAC8KyAHD3+5pv6t1azBMkff/wBAN26dWMrVKQZOBpln7e39/Hjxz08PM6dO1dzJRO7MjMz6br3M2fODBkyhKNe3lRcXFzfvn1NTEwkEomNjU1GRoZcLlcoFFKptKysTKFQyOXy4uJielxcXCyXy+kxfUNZWVlpaWlxcTHDMADw8ccfnz17FgACAwMXLlxobW2dkJDQokULfj/jzZvw/vtgYAAPH0J2NuzfD999B3Z29WxNLpeHhIQ8eiDz13wAAApCSURBVPRo165dM2bM2Lt375o1a2ixAtRYYBpl2fHjx+nS+nv37rE7nf+vli1b5uTkaNW8HgC8vLxOnz7t7u4+YcKEZs2a0am6ubm5iYmJkZERnb+/8rdLeXm5QqGoqqqytbV98OBBz549FQrF8ePHR48erZlPUTdvb4iKgr174YsvWGuzurrazs4uPz8/MTGRLntAjYXOq9+C3kS/fv1GjRo1dOhQrnMoAHh7ewcFBV29erWgoICW1tcGkyZNio2NvXr16tWrV1/2HiMjIyMjI3oZ1NjY2MjIiObZWtdM6VXUFStWKBSKadOmaUkOBQB/f/jiC+jR459X7t6FDh3gJTtSv1plZeWRI0fy8/M7d+6MObTRwdFoQ9FcRhczAsD48eM3bdqkgRwKALGxsf369TM0NBSLxVrys0cIGTBgQHR0tKura8+ePen8vaysrKSkhM7fy8rK6NZ+r8/W1ra6ujolJcXMzIyjsOth4EB4/BgSE4Huode6NVy9Cq+8FCyVSrOysrKzs1NTU1NTU9XHT548UalUNjY2bm5uUVFRGogfsQhHow11+fJlhUKhPr169WpZWZlmuu7duzcdtamTeFRUVLdu3RwdHTUTwH9FRERER0e3aNHiwoULdWQ99WVQmUymUChoni0pKaHHMpmstLRUoVCUlpZmZGRcu3bN1NRUC3/fm5vDxo2wbl3t1ysqIDMTsrIgOxtSUyE19dlx8+bbz5+f/8KmdHV1ra2tc3Jyzp8/L5FI6n1rDvEC02gjlpOT4+jouHbtWrrOPzMzc8uWLYsXL+YrjSoUihUrVgCAn59f3SNHmv3peoZX8vDwuHjx4vbt27/77jt2AmXJypWwYAFMnAgdOwIAZGeDtzekp0NBwYvf36vXO5aWlvb29g4ODu+88469vb362NbWViQSTZw4MSIiYvLkybR+lSY/C2oITKMsOHjwoKWlJT3W2FAUACIiIhS5ueYSCdy8CQB/nD6dlpamsd7/a926dU+fPnVxcZkyZQqLzX777bcXL14MCAhYuHAh3cZDS1hYwPffw5dfwuXLAACmpnDrFgCAnh60bg0tW4KdHbRtC23bPjt2dBxuZvaSFAsAADt37oyNjb127Zqfn596d26k/TCNsqCysrKiokLz/epVVIwoL++/eDE9nWBkFGZoqPkwKIlEsnXrVqFQuGPHjrq3L35TAwcO7Nu3b2xs7J49e77++msWW264mTPhwAE4ehQAwMQExGJo3RrqvSLLzMxs3759H3/88erVqz/++OP333+fxVARh3hdtdoUWFpa1qwL2bJly6SkJM10fWT2bAJQKhQ+MDWV6uoSgHk6OsePH9dM77UMHz4cAHx8fLho/PTp0wBga2vLdaH71zdgAKH7P4nFpE0bYmVFUlLYaZluOdOxY0ft+bCobmyOGhAvbotEgy0sftfXBwC6d57mnT9//tSpU6ampuv+e8OFDcOGDXN1dc3Jydm3bx8X7TeEiwt4er70emg9+Pv7d+3aNSkpiV5oRtoP02ij19Pc/Mn48ZP4uztfVVU1f/58AFi9enXLli056mX58uUA4O/vX1lZyVEXb8TLC9SFt9avhzlzoMZuBg2ir68fFhamp6e3bdu2M2fOsNMo4hTfw+FGT6FQ1Cw7X15errEq9CFLlyoEAgJA/zACwXhLS81P6jdv3gwATk5OSraKGL8IwzBdunQBgL1793LXi/ag4/pWrVoVFhbyHQt6BVx+34ht2rSp8NSp79zd6f3rX1NS9mRlhYSE2NX7Ae83l5ub26FDh+Li4j/++MPT05PTvg4fPjxp0qR27dolJyfXLBjYJDEMM3DgwOjo6IkTJx46dIjvcFCd+M7jqP6io6NPnjypPo2Li9Pw9pzkeX1MzRReqq6uplurRkREaKA73kkkEvoL8siRI3zHguqCo1FUf2VlZS4uLmlpaYmJiZp58CY0NHTmzJmdOnVKSEhgd1mVdtqzZ8/s2bPNzc3v3r2rmSeMUT1gGkVvoKqq6vr163379qWn1dXV0dHROjo6/fv311gATk5OT548iYqKGjlypGY65denn3566tQprusuooZo+r/PEYsKCwtr1jYtKSkZMmSIxnIoAOjq6i5ZsgQAfv/9d411yq+9e/fSGgU7duzgOxb0YphGUSMzffr0AwcObN++Xf1KampqzeowTYy1tXVwcDAALF26NDExke9w0AtgGkWNjKGh4ebNm2/fvq1+ZcyYMTdu3OAxJK599tln06ZNUyqVU6ZMqaqq4jscVBumUfRmysvL33mue/fufIfztggMDHR0dIyPj+foOTHUEJhG0ZsxNDR88tydO3f4DudtYWJicuDAAZFItGHDhr/++ovvcNC/YBpFjdKYMWPaPveWXDF0c3P7+uuvq6urJ02aVFpaync46B+YRlGjFBkZ+ei5zp078x2Ohqxdu7Z79+4SiWTZsmV1vC05Obnm8gmJRNKnTx/uo3t7YRpFb8DExKTmD7ChoeE333zDSyRCoVD0HC8B8EJfXz8iIsLAwGD37t20eOALfffdd48fP1afrly5suYpYh2mUfQGTExMVq1apT41NDRcu3Ytj/G8hTp16vTDDz8QQnx8fPLy8l74nri4uImlpTB7Nv1TdPkywzAajvOt0sTrO6AmaezYsS1qlJgfPXq0Jqux8G7x4sW//vprbGysnZ0dHYwLhUKhUCgQCOggfSbDfFlaCnv20PcfFAi6aM3+200SPgyKUOMzadKko0ePqlSqF/7tKX19r4qKuC5dsqytvS9frgJoZW39sqErajgcjSLUyJw/fz4iIkJPTy8oKMjCwqKoqCi3htWrV1d99hkAZHbq9KBDh2f77SEu4WgUocZEJpN169YtPT1906ZNS5cufeF7ThsafqJUqk+rBQI7KyscjXIHbzEh1Jh89dVX6enpH3744eLnO8L+1x4jI0WNKoKhpqZWeG2US5hGEWo0Tp48eeTIEWNjY/pE08vedsvY+PHdu0AI/fOjtfXhw4c1GefbBif1CDUO+fn5Xbp0ycvL27179+zZs+t4Z2Fhobm5uTrPFhUVNWvWrMlvu8IjTKMINQ6jR4+OiooaPHjw2bNnsX6zVsE0ilAjsG/fvunTp+NuItoJ0yhC2i4jI6Nr164ymezw4cMTJkzgOxxUG95iQkirMQwzZcoUmUw2YsQIzKHaCdMoQlotMDDw8uXL6q1EkBbCNIqQVjt79iwAhIaG2tjY8B0LejG8NoqQFnn69GlZWVnHjh3paUZGhlQqzcjI8PT05DcwVAdcSoaQFjl+/Pj9+/dDQ0Pp6cmTJ//++++DBw/yGxWqG07qEUKoQTCNIoRQg+CkHiHtcuzYsfPnz9NjuVzu5eXFbzzolXA0ipB2GTdunHoLa9yjpVHANIoQQg2CaRQhhBoEr40ipEX69Onj6OioPu3Vq1erVq14jAe9Dlx+jxBCDYKTeoQQahBMowgh1CCYRhFCqEEwjSKEUINgGkUIoQbBNIoQQg3y/ysSgXTRUQA9AAAC3npUWHRyZGtpdFBLTCByZGtpdCAyMDI0LjA5LjUAAHicxdFrSJNRGAfw53131W1uzpxmS+d1RmJIJijbzsFCIiPBSKks3lLLGxQtkj4UJilJQvUhTVKnkNch3snbds6XNNMUNdMwbEZmpSKl0AW7zKYFEfSxB/78eOA55zlwlqwN02AvhT0sOCrIHq09OYwQNHb5Qhxsh/dLhUN2oxc55P059y/FsH7P+pq/szbD8kUOeY7+99kNmZ8C87+UAAMMm8ayPI7H5/iCNJYv1AhFaaxIzImd0lgnZ85ZYkeqkcpA5gIucnBRgNhV46oM5rFuGrdNwaybu0bAcBIhByoADwBPEG0GZy+QbAGZGmRbwcUb5D4g14DcFxR+oPAHRQC4B0L02tcJGJbHFwhFYidnidBjlLE/ChyBIJ8nbfi8TUl9S0t1reYSTPYVEPn8+25Jbg3eNpJt0IryLEFfsnFj7D10x6PYMvZqJw6byLdExdZ2hYYfxTPDZWTq8pxOkWNGmUVm/cTcoO6QNhSp1XVImtXf3XElkiQutqCPexp1cjUlUXlZhvwXVfrJoSNkujiHjAW+1l9NlSBjTAHhDjbo699FUA4PGYz8cf2FgUy662Y0ScoJNjzGVTRF3WQFU6ThTXIz5Z6PGXKHr5Pq2kY6dY5BMWVtKLHDTJaNIZjtje8xLs8Sz+Gz+NlqiOXbwij2botE7Smzlu3ty3jqdoZ+9yWbJd0yjnvv8/FqvNYqVVbiupr9dDYuX9954yEWNKnoW9lJ0vSoH8/kKmlZXCqShplQkJcCP51LsCQXtJLq1E8k4bifQZWVi65d5NFBn5d6kwdHR063ENveCBJafow+gHpyd6AAhVd0UBtUkqGkdMOBujEabgpAzZOF5MyiiSZ8b0B9C0Wk4kQ97RFj64fCr+Tz7Ai13dKhw+mdaMa/nGb4daG6oVYUeqqG9q0u9RTvcMcLJSt4hZPjyfloq+oHA7AE1Pz6beEAAAPqelRYdE1PTCByZGtpdCAyMDI0LjA5LjUAAHicfVbLbiQ3DLz7K/QDbvAt6ri2F0kQrA0kTv4h9/w/UpTGau0lM9PtaU01VSwW2X5q9frj7fd//m37pW9PT63R/3zGGO1vJaKnH62+tJfvv/z23l4/v718rbx+/PX++WfTbDpwD94/Y799fvz4WuH22vwKpxzWnuWSYXg3uogHc79vFQDtShKh0Z75oi5etwBJXVNvpM6QohbiCCQmVLcAGDbcb6A1kLyyW5g2viTVgyawa+8H0BFRrmRK4xXI0xaSWNJuZLR3hDRxqkiIqTIoKh1J8L2BHSH5Uh1p8+fUXjkA58RHMgkcXYG8M+sGgSiPnYPt0GcAiOUuFhqFHD0jZ0TWcSTDVEi+nCi99hal0fskgeoeezOvmOnWeYlOiF6iMov2Y3eWydPNnPNRHug/g+IWOZDaPqrQ6QmZATAq4SsmGdQ8kFbbK5I3Aa1KLmiKc3GY+wmtCj17rfMqoZqOqKQ4WTof0FhQNwXZ2p+1D56ch0g/3Mp9QS2HasVy8kARGfcw+YnMpakaeyi8EsMc+lSdYJA45R+wCKB9wI0Cs/SQeJjJRPTQX6pScflgw1XV1o28kJAf1TuQDFH7Fexzr4u8S5eJpAGfHEgBMq7oQwGQaXWeRDnU6Nxd26/oI1I116oDwyU6PRrueagvBmRcEh31r9oq/q6KZSrJkT2+F1TJB5Kvfs8hWjbN4b1aYSMDSDSQKeoPouLsI6ZMDmudk6EDCUVhpNDZdcPQfmU+6FH1uKE5gxLXrrU9+fRWKUp2mk9GBVWUnM1k1haKPWwI9e3oZXwmFobvffaJSrJVv1IPlyMr5YI6io5JFHNCOBq/Bh68l2f3oWkARalAoAbhhSkFunUPE+mZllapnu0aVSGr5vehmHQFhbfpsL/aIqAyo1bdszsClGxJGYdX1BcDRzQpkbpGIBiktnQ5SqDxIGCcSQVwRPf6QtbtjNnX/qIsMtUkGYL5BNcokx5NpVWsXmlDmHIACrwKwIPUj/y/v7/99FBZj5mXj/e3+zFTb7kfJrBU0/uJYXXcz4V6+z39cdHiHvGMy34PcsOR97w2HOOeyoaDz+FrdWI+hqzViTc52JAnZvPDzON52hTR0jxPmyVGIM/TJsrxWNlcUXAuOpst5pbWiTdhjCfGmpxDSOZpM5YpZY2SY6gYdLmvS1CksTcWW/rvjMTXwk5IYi3sfKS4JtaP3uUSSu6oxRXbyM5nVRq13fkorxU9+2rKojshnfIGfjm6ZK3cFvHHymask3HHL4e318rmrPlY2Zy1OEMv3ZzLvKdV6/rrvyl8f/oPxSrLp3M2vJEAAAKnelRYdFNNSUxFUyByZGtpdCAyMDI0LjA5LjUAAHicXVO7rhRZDPwVpE3mSjMtvx+0iEjYhEuOCNBEV1oWtBDy8Vs+fQOmg+lpdx27ynadzx++PL+/vHt+uvP98v7y+cOXp/WY37/3y+0u9xUD/Osjns9/nHl9PWWdwtd0medPfvP7kluYZl9lEyPp6402blLXPbbI1gRCqSkLCTUaxJtNr7SluJEviIWodpRjrwREnpJcC6MOtt03UQsBJobDi4vC2n23rUiErzcGmfgBgbYkA3nhVPgGjY3jU5ybOackqRoEKFRnjiQK9xKGSIl0myx1s1WySkkckJL3UbBaDFnVnha7biZOMTpEpXnqSQk62mUrprJDtJcRHyJZSrt23hR/S1tp1lKCUiy90xboruqKMwLZRx4mEs3IE9Oomb9zLMjFImW/zXgtNCav8xiqsLYDwydjD6RFm68pj5JAD8CyMdarbhkSh2ITUdXBHGuK2YESJrb0UrfkjmmUF80nOKF0emFCn6BTtGCHOzLQ87THYe4LQ5tmMzPtsdIcghusB8PCsxYmxX6d/QbMdfPJPxpW046h5WJYZrCM9mMLnjX7EVTpleY2C4dGBtfhmJaRD79CxtgLjiZfgphIMRHberkDKrzR2CBInnIqkzMmrXTVgbowocEMlEPu5AFPMjiZfJgcZYYpNQJGAGI1XYHIeG3aHC/zTwaf3sb7LIsIt0zaFatTZdym8aVbcq9eCBu/ruukmbMv7CIOWUazN1wvD4WpYGrner0yFmuTSOaEp4jh6/E35oC1TT0yLFK3gs0NumbBh91wudP96fr11/dvn/77/uMtbS8///7245+X+8uvjd/mY1iPYT+GbI+x8CmWU3w+H6f4xC4nej2J1ROfnvj0xKePfL//B8t9Qh4cK0A+AAAAAElFTkSuQmCC\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x75ee4ffd4040>"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"smiles = 'CC1=C(SC(=N1)C2=CC(=C(C=C2)OCC(C)C)C#N)C(=O)O'\n",
"mol = Chem.MolFromSmiles(smiles)\n",
"mol = Chem.AddHs(mol)\n",
"\n",
"AllChem.EmbedMolecule(mol,randomSeed=42)\n",
"AllChem.MMFFOptimizeMolecule(mol)\n",
"xyz = Chem.MolToXYZBlock(mol)\n",
"mol"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "4b431122-f0f1-4f3f-83d5-ba3efe41842c",
"metadata": {
"execution": {
"iopub.execute_input": "2025-02-07T13:11:03.960540Z",
"iopub.status.busy": "2025-02-07T13:11:03.960403Z",
"iopub.status.idle": "2025-02-07T13:11:03.965870Z",
"shell.execute_reply": "2025-02-07T13:11:03.965537Z",
"shell.execute_reply.started": "2025-02-07T13:11:03.960529Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxN6f8A8M/tar9apEXRnuQijWSNyDZKMxFm7DGEMbYMMZSxxlCa7L6W+GUJk0nJVxhfRohSUpZK+0rLbb11u/f8/jh1HHUn6W7VfN4vf5zz9HTu0xifznnO83w+DIIgACGEUFvJyXoACCHUsWEYRQghkWAYRQghkWAYRQghkWAYRQghkWAYRQghkWAYRZ/YswcSEhqOS0vh8GEAgGvXICbmY5/r1+HRIxmMDaH2CcMo+oSXF8ybBzweAEBxMezeDQAQGgpPnnzsExYG0dGyGR5C7RCGUdRU//7g5yfrQSDUcXSR9QBQu7N7N9jawsyZnzRmZ0NsbMPxhw/SHxRC7ReGUdSUgQF4esJPP4G//8fGW7cgNbXhOD4eRo6UydAQao8wjCIh1qyBoCC4efNjy6JF8NNPDceLF8tkUAi1UxhGkRDy8nDoEEyfDkpKn+lJEMBgSGVMCLVX+IoJCefgABMnfqaPtzcsXQpffw3FxVIZE0LtEgMT5SG627fB0bHhBrO0FF6+BHt7SE4GNTXo2bOhz6tXoKoKhoZQUgLdusHGjeDmBoMGyXDUCMkShlEkkp9/hhcvIDwc5OVlPRSEZATDKBLV8eOgpwcuLrIeB0IygnOjqI0IAlauhEOH4Pp1sLKS9WgQkh28G0Vtx+NBSgr07AlqarIeCkKyg2EUIYREgg/1CCEkEgyjCCEkEgyjCCEkEgyjCCEkEgyjCCEkEgyjCCEkEgyjqLVu3LhhZ2fXo0cPHR2d9evX8/n877777urVq1SHuXPnXrx4UYYjREgmMIyiVklISJg5c6aPj09+fn5qaiqLxaqpqamoqKirq6P6VFZW1tbWynCQCMkE5htFrXL+/Hk3NzcnJycAUFNT8/b2lvWIEGovMIyiVklLS7O3t2/efuTIkZuNWfLj4uK+/fZb6Y4LIdnDMIpahclk8vn85u0ODg6Ojo7kcXp6unQHhVC7gGEUtUrv3r1jYmKat1tZWY0ePZo81tLSku6gEGoX8BUTapUFCxbcuHHj5MmTXC63urr61q1bLffn8/mnTp0qKCiQzvAQkiEMo6hVzMzMoqKirl271rdvX1tb2z/++IPP5xsZGamrq1N9DA0NNTQ0yOODBw+GhISkUkWZEerECITELSMjY/Hixfv373/w4IEkrv/bb7/5+vpSpwcOHNi+fTtBELa2tu/evaPahw8f/ubNG0kMoHMLDg4ODg6mTi9duhQUFEQQxKJFi/Lz86l2Dw+P7OxsGYyv/cG7USR+a9euXbNmTU1NTXV1tSSuX1paWlpaSp2WlZWVlJQAQE5ODo/Ho9pzcnLoy1pRKyUmJiYmJlKnSUlJCQkJAHD9+vXKykqqPSIiory8XAbja38wjCIxq6iosLCwCAoKunv3bmRkpKyHg5DE4Zt6JGZdu3b19fUFgKtXr1paWkroU6Kioqg7o6dPn44YMYI83rZtGzU/S79j7YR8fcHNDczNG0737QNnZ+jTp2m3qioICoL4eFBWhpEjwc2toXx2i16+fBkUFEQeJyQkmJmZkcehoaE6OjrksYQeNToiDKNIUqZNmya5i5uYmEycOJE8podLe3t7fX198jgkJERyA5C9S5dg2LCPYfTyZbC2bhpGuVxwcAALC1iwACoqYNcuuHcPDh367LVLS0vT0tKoY6o9KyuroqKCPK6vrxfHj9EZYBhFHZK5ufmUKVPI47i4OA6HQx6PGTOmd+/e5LGysrJsBtd+BAcDkwnnzzecjhoFpqawejVYWLT8ffb29tu2bSOPfXx8qBv/VatWmTcG7tOnT0tkzB0QhlGEOqyYGKipaThu/EXyiefPYdKkj6fa2mBrC/Hxnw2j6ItgGEUdj4aGBn1nqrq6upycHAAYGBjIy8tT7fr6+vTTTujBA6BW5gqdCC4rA1PTT1rU1aGkROID+5fBAssIdUw2NnDgADTuxIUhQ2DHDhg//pM+mzcDhwOBgR9bBgyAfftgwoQWLkw+wrNYLOqUIIiuXbu+f/++W7duTCaTbG9y+m+GC54Q6lwEAigqAoEAAMDZGS5cgNzchi9FRgKHA6NGtXwBFotFxVDytGvXrgCgra1ND5pNTv/N8KEeoU6kpgZmzYJ+/SAuDq5dg6FDwcsL7Oxg9GioroaYGAgOBiUlWY+ys8GHeoQ6prdvwcAAVFUbTlNTQU8PWCyorgYVFfDwgA0bGiZGORx4+RKUlYHNBkVFGQ65s8IwilCnw+HAzJkQHg5d8HFTGnBuFKHOpbgYli4Ff3+MoVKDYVTiEhMTqQ0hAJCUlJSSkgIAERER9MQZkZGRWA8OiaqwEGxtQV4egoIgJUXWo/m3wId6iVu2bJmJicn69evJ01WrVmlra2/evFlBQSEvL6979+5kO4vFevv2LbWRsYOqqKhISkoqLy8fOnSomppaQUFBZWUlte+lsLCwvLzcAtd+Sw6fD1TWJRYLOvey2XYD70Yl6+LFi3l5edTp5cuXs7OzAWD16tX0bp6entIemQSEhIQYGxv7+vqeOHGib9++jx49CgkJ2bFjB9UhNDTUx8dHhiPs/JhM0NRs+IMxVFowjEpWVFSUTnr66MhI8PICL6+8s2eLiooAICAggN7t4MGDHf2x4P379+7u7hcvXrx27drly5ejoqKsra1FueDDhw9dXFzMzc0HDhx49OhRAFi4cOG9e/eoDh4eHlFRUSIOuzW43NcCAZd2+kYgEJrcSMDj5fF4hVIYEmpXcBJaskzfvx+XnT0kMRHu3QOApQCnFRSobc42NjbkLkYej0fPN9wRPXr0qHfv3uMbd9FYWVmJcrWkpKTJkyf//vvvFy5cyMrK8vf3r6qqSk9Pp+cJzsjI4AjdSC5uqanOpqaXVVRsyNO0NFdj49OqqkPofWprU1JTXZjMrgQhUFKyNDEJlsLAUDuBYVSyehUXDykrS9fR4bu4GNy4oZyXN6ox/5BAIDhw4ABZy2jy5Mkdffd3VlaWkZFR8/YrV65Qt5AVFRVUdruWnTx5ctq0afPnzwcAKyur48ePi2+kElFQ4KupOU1ffwcAEETH/o2IvhSGUWl4ThA3CGIdQdCTQRIEERoaqqSkBAB8Pr+jh1EtLa0SYTkv3Nzczpw5Qx4fPXr0/v37rbna27dvqTx4dOfOnXv8+DF5nNKe3kQLBHVcbipB8BkMJoPRsf8q0ZfCMCoNzhUVU//6C8rK6I1MJvPAgQPkm/pz587JaGhiY2Nj8+zZs/z8/B49eoh+NYFAQE53NGFkZNSvXz/ymNzoLR2pqU5UcOTx8pt30Nf3TktzS0w00tKaq6u7rksXLamNDckcvmKSrFcGBmUqKgpcLrx7BzU1OVpaDBsbR0dHAPD29laldvIBbNy4UU1NTXYjFVWfPn3mzp3r6Oh49uzZiIiIDRs2fPZucfXq1dSNahPm5uZkGbUmRo0aNauRnp6e6MNuJXPziP79M8k/iormzTsoKlr07Ztgbv5nbW1KSsoEgI79whB9EQyjklXUtev5X36BtDTyzz5XV302e9iwYQCwcuVKenr2pUuX0tPqdERHjx7duXPny5cvb968aWZm1qNHD1tbWycnJ6qDjY2Ni4sLeXz+/HklJaV3794JvdScOXPOnj179+5dACAIgl6oUiiBQPDHH3+UfXq/L30qKoOMjE5VV8fx+dJ494XaCXyolywfHx8lJSVorALm6e2toKAAAHFxcfSg+eTJE01NTdkMUXwYDIarq6urqyvVMnz4cHqHIUOGDBkyBACKi4uvXbvm7e39T+WS7OzsTp8+vWbNmuLi4i5dugwZMuTSpUssFos+g8xiscj/mABw/PjxoKAgExMTGxsb8f9gwhAEv6DAt64uS01tvKamW0nJ/ykoGMnL65eUXFRS6s1kakhnGKg9wF1MSAbc3d3Nzc2VlZWjo6ODgoLokxtNEATB+Fwly7y8PG9vbyMjI2dnZ0mE0by8rd27/6Cg0JM8zc/foaU1R05Otbo6QU3N8e1bRwuLW2Vlf5SWXuXxCpSULPT0Nioqmol9GKjdwrtRJAMeHh6VlZWZmZlKSkpCXyVRPhtDAcDLy2vXrl2Sq7Cmr7+Vftqjx2byQE1tHAABIGAwmJqaMzQ1Z0hoAKidwzDaoLa2Ni0trW/fvuRpXV1dSkoKm83Oz89nMBjU24yCggKCIMTyMvrfbOjQoQBQWlpqbW0tYv3ODx8+8Pn87du3x8XFlZSUSO2hnpSbu0lb+yeAz8d61InhK6YG6enp42l1bLKzsx0cHADAz8/vEK2ud2Bg4IEDB6Q/vE5JU1PT1tZWxIt07949ODj42LFjHh4eS5cuFcvAWkMgqMrImNeli46Kik27XW8fGRlpZ2enr6+vra3t6elZX18/e/Zs+nz0ggULgoNxw5WoMIwCABQXF5+nankDlJaWBgUFAcDdu3dzqTo2AP/73/+ysrJkMD7UCj/88IOlpaXUPo7Hy1NSYhNEXWnpZYGgUmqf23ovX76cPn36li1b8vLy0tPTNTQ0qqurKysr6ekZq6qqMD2j6PChHgCgsLDwxIkTs7hcOH4cAOoqKo4ePgwMRkREREZGhplZw+uCmzdvpqWldfRcdkgsFBUt9PQ2yHoULTl//vzUqVPJzWAsFmvLli2yHlGnhWG0wfc1Nbs4HPDwAABdAE8GY59Ww0aUwMDAs2fPAkBZWVltba29vb0sB4pQ66SlpdnZ2TVvP3bs2K1bt8jjp0+ffv3119IdVyeEYbTBjvJyRYBSV1f5oiLWw4cDG9cn1tTUfPfddytWrACAgICApKQkod/+4sWLdevWJSUlcblcJyenQ4cO/f777wDwyy+/kB327t1bVVX166+/SuWnQQi6dOlSX1/fvN3BwYF6DUBmv0UiwjDagFw/u7imZkBtrTetPTMzs7KyMicnBwBev36t07iQnq6srGzcuHFubm5r1qyxs7M7ceJEYmJibm4umb0JALKysnJycsSQfIQggMsF6tV2k9PmCgtBWxtaXFGEOqvevXs/ffq0ebulpeWoxlL1VPEFJAr8B9aAXJ94pbbW+9Nf4H369Pnuu+/Cw8PDw8OnT58u9Hu3LF9uUlKinpr61759jgMGTJky5caNG3FxcVSHM2fOxMTEiGGUb982lMwlZWRAz55CutXXw6pV0LMnTJkChobwww/A5Qrphjq1+fPnR0VFHTt2rLq6urq6+ubNmy33FwgEp0+fphdrQK2Ed6MN7igqTuFy4a+/AAAYjNsKCosWLeLxeL169SL3L5IMDAyabHBMi4nRuXTpiUAAUVEAsFNfn9GYeejGjRuFhYUAEB8fn5OTM2LECCn9ML6+EBsLSUmgrg5VVTB1KvzyC+zfL6VPR+2DsbHx7du3d+zYsXfvXiUlJXt7+3HjxhkZGWlofNyoamhoSO1CPnToUEREhImJCb5E/VIYRhus1NKaEhQEBAEAOeXlZ5ctK/T19fT07NOnj7OzM9XN3Nz8m2++oX9jzuHDWwSCSkVFVr9+kJIin5cHr16RX2Kz2bNnzwYAgUBQU1MjvR8mOBgCA4GcUlBVhV27YPJkDKMdS25uLpfLJVeJVFZWEgRBJQasqqri8/mtyQc2ePDgP//8k95CTtlT9jf+X5GVlZWQkICvm9oGwygAgKGh4alz52DMGPK0W3U1uYx08eLFTCaT6rZgwYLm30vmKo5SVv726VPG11/Df/9LtsfHx79+/To6OhoAOByO2JY0lpSAo2PDsdBHdYKA9HToQ8sQ3bcvFBVBeTm0KRFfTEzM5s2bk5KSCIKYOXPm3r17vby8+vbtu2jRIrKDl5eXkZHRsmXL2nBx1FxkZKSHh4e2trZAIODz+aGhocHBwaWlpf7+/mSHwMDA7Oxs+q4Q0a1bt27r1q1hYWHV1ULLTKGW4NwoAACLxRrTGEMBQEVFhUwJ2qdPH3o1YEtLy+bRkLxHsK6qejZoEPHyJQDEx8cLBIKBAwcuW7YsIyMjIyNjzZo1Tb4rNDQ0OTm5LWNVV4djxxr+7NwppAODAcrKQL/5ra4GJhOUlNrwabm5uRMmTJg5c2ZmZubr169VVFRKS0tLS0srKz8uOC8rK6OfIlEUFhbOmDEjMDAwNjb2+fPne/fu1dKSeAboqqoqU1PTs2fP3r59OzIyUtIf1/ng3aiotIcP550/b8rjmT5/DgA1cnKn796V69tXRUVFW1ub6qasrEy95b93797Ro0e/+eYbagv/F2AyobHsO9DulD/BZsODB0D9AnjwAPr1g8accl8kJCRk6NCh5I2nmpraTqGBW6ZKSkqSk5MVFRW/+uorJpOZmpqqqalJhZ60tDR1dfUO9D769u3bbDabmjia1Fj9UKJUVVV9fX0BIDQ01JT+DhO1DoZRUQ1Yvtzn8mXF9HRHR0clJaVbxcWbDx4MCAiwtbVdu3Yt1c3a2nrDhg0AUF1dvXDhwjlz5ohzEHV1cOQIvH8PixeDkRFs3AhLloCuLtjZwYsXsGYNND4PfqnU1NRBgwY1bw8JCaHuph88eEBt9JIyLy+v06dPjx49uqCgoKKi4u7du56enm5ubnPnzqU6TJo0iZp/aP8yMjJ6NxY9pDt//vyDBw/I44KCgiYT9K31668wcuTHSaFdu2DQIKAVGWzIFVtaCn5+QK4tsbMDT0/QwPSpLcEwKioGg+Fz587Vq1evPn3K4/EGu7jk5eUNGzaMmlR9/fr14MGDq6urU1JSLCwsXFxchg4dWlhY2HKCOOFUVcHB4eOpikrDfO7Tp9C7N3z7LaxcCX/+CU5OcPo0nDwJ3t5gYgIHDkDb/tUByMnJ8fn85u0DBgyg/iWTi2qlLzw8/OTJk/Hx8QYGBgDw+PHjTpD6umvXrkJnSKZOnbp9+3by+Pfffy8uLm7L1Z8//2TBXHw86Oo27cPjwfjxMGgQBAQAAPz2G4wbB48eQQcvuShRGEbFgMlkzpgxY8aMhnSTRkZGERERVOW1KVOmXLp0qaSkZPbs2ZcvX+ZwONra2vfu3cvPz9+8eXOXLl/yV9CzJ1y48PFUVxeuXAEAIJdSVVYCFfImTIAJE0T9wQAsLCzCwsKat/fp04d63rx27ZroH9QGUVFRc+fOJWMoNCbfE1FVVdW1a9fevHmjrq4+e/ZsPT29M2fOjB071tDQkOxw9uxZe3t7ExMT0T9LqH79+u3evbumpqZJ/kAlJSVqakJFRaWNYbQ1IiOhpgaOHGnYtXH8OLDZcPMmCCvUikj4ikmqli9fPn/+/N69e+vo6GhpaX1ZDG0ZQcC6deDlJbYLAgDAjBkz4uLiDh06VFdXJxAIqIrzLbhw4UJ6erp4hyFUdna2sbFx8/Y1a9YYNbpx40brL1hSUjJo0KDLly/36tXr/fv3Y8aMqaurCwwMTEtLo/ocOXLk7du3og/+nzg4OFhYWEyfPj06OjopKSkgIOCzBaa2b99++/bt1n7A7dtw8GDDn9RUIR1evoSRIz/ufGMyYeRIePmy9T/CvxDejYpTUlKSEu2FeHJyMlksKDo62tzcHAB2797N4/HOnj2rr69fUVEhtg+uqIC1a2H6dBg5UmzXBAAAPT29O3fubN68ed++fXJycra2toMHD9bW1qb2uQKAtrY2taL76tWr//nPfxQVFSV3v0bp1q1baWlp83Z/f39qbvSfNp4J5e/vb2hoSN1cb926VaFN7+VEIScnFxkZGRgYuH///rq6ukGDBhEEwWazq6qqqD5WVlbU68qQkJCMjAxVVdVx48a16gPq6j6u4hA2XQNcbtO3kcrKgKugWoRhVJwOHTpEPv1NmDBBXl6+tLRUTk6uqqrK09MzODg4ISGBzFqWnZ2toaFBPSeKwV9/gZwcXL0KN2+Cn5/YLgsAADY2NhEREfSWPXv20E+pOTsOh3PlyhWhq2slYeDAgWfPnvXx8RHXBePi4uhhl/qNmJGR8bLxdkwK2yhUVVW9Pn2qmDZtGv2UmpXmcDihoaHLly8XundeuMmTofF3DAj9LlNTaJLIOTkZ5s1r7fX/lfChXiJu3bqVmZm5cOHCtWvXUg+e1tbWmZmZmZmZLBbrxIkT4vw8F5eGlaTijqFfxMvLa8uWLa2pniQWCxYsqK6unj59enh4eFhY2Lp161quz8jhcObPn/9POboA4P3790KXRh07dmxDo8zMTDEMXUw8PT2dnJzevn2bk5Mj9E1ga3G58OgRvH8PAODkBLGxQG3Av3ED4uOBtpEPNYd3oxJUUFDw7t07DodD/vOurKx88eLFgAEDyK9WV1fHx8cPHDhQpmMUG4IgCIIICAhISUlRU1ObOnWqpD+RxWJFR0efO3cuLCxMWVnZ3t6eIIgmq3GdnZ37NG7o2rJli6qq6nsyWAhjYGAgNHHc7t27qd0Zw4YNE+sPIZKxY8dWVlampqaSYZT5T+uIW0YQ4O4O48fD1q1w+jTo68Mff4CHB5DlWmtqIDQUJL8FoEPDMCpme/bsqaysnDhxYnl5eXV19ahRo3g8Xn19vYeHB4fDGT16tLq6OrmipaysLCwsrNOEUQaDcfToUQC4c+eOoqKidD5UTU3txx9/pLcsXLiQfjp//nzy4P79+02mdJtzcHA4derUkiVLlNq040v6Zs2aBQAZGRnZ2dmtmsa9dAnobzXPnQMmExiMhuUfOTmQkgL6+uDgAG/eALkYAANoK+BDvZht2LBh48aN/v7+CgoKixcvHjNmzOrVq7t06XL9+nU7O7sNGzZkZGSwWCxdXV2hq6w7AUdHx5HiftMlIh6Pt3r1ahMTk/j4+GfPnv1Tt2XLlhkYGHz11Vc///zzsmXLJtLWpQv14MEDCwsL6SxLaIGxsXFrKzIoKn6y801R8WNU5XDg8WOgJTMDLS2Moa2EYVT8tm3blpKSMmXKFPKmxsfHR05Oztzc3MbGpqqqqqioiCAIFotFLSxFUrB37149PT0VFRVV8llVGAUFhYiIiAsXLgwaNMjZ2ZmsoBkQEEDNwwDA/v37yW1dtbW1gYGBzs7OIk1KthP5+bB0Kfj7ty3xAsKHeokoKCjQ1dVtXnPxzJkz//3vf2tqahYsWFBQULB8+XKZDO/fRl5enlwPpKmpSU90IJS1tbW1tTV12uTOmso2u3379lWrVgndm9DB1NTAuHFgaQl+fjBvHkgtK24ngmFUnL799lsyC6S+vv7r16/nzp1L7kXZvn07i8UCgB9//NHLy8vPz2/x4sXe3t6fuRwSN6H5AdogKyvrzJkz79+/f/LkiaKi4rZt28RyWdlQVoZ/Xr2AWgPDqDhNaNx/OWzYsN9++83Pz49MNfTzzz/Tu9FTlqCOyNDQkMwksGHDhnm4pvJfD+dGJcLBwcHNzW3w4MFbt27ds2ePo6Njy1v6ysrKRo8e3SRROWr/Vq5c2atXL1mPAskYo+UVy0gUMTExcXFxDAZjyJAhAwcOfPz4saqqav/+/amvKioqktNwK1as0NfXNzQ0bCGBXm5u7o4dO549e1ZfX+/q6rpp06bjx4/LycktXbqU7HDixAkul/vTTz9J4UdDCFHwoV6C7Ozs7OzsqNMmKYioL/3vf//r0aOHpaVlCxsNa2pqxowZM3HixKtXr/J4PF9f37S0tKysLHq2vezsbKwAgZD04UO97K1YsYLP51+7dq2F+g0RERFycnIBAQGGhoZmZmYnTpwQW3EnhJBo8G5U9u7fvw8A169f5/5zNfk3b94MGzaseabn6OhoqrDH33///dVXX0lunAghoTCMyh6Zs33SpEkCgeCf+vB4PKHJSVVUVKicaS0sLJc0giCeP3+enJyspaXl6OiooKBw7969/v37UzWR7t+/b2Vl9dk1mwh1RPhQ317o6Ojo6en901dNTU1fCkudO3DgwMWNbGxsJDnAf8TlcidOnDh//vy4uLgjR44MGTKEx+OtXbuWPuD169fHx8e3/prV1dXp6enkBqH6+vry8nLqS3w+n8PhiHH8ssXhcKh8KLW1tfS8onV1dVhytUPAMNoxODk5vXnz5tSpU+Tpu3fvPrvE4tmzZ9KpkhQYGFhSUvLs2TM/P7+wsLDQ0FB5Eer2ZGdnjxo1qk+fPm5ubiYmJpGRkTExMSNoW2ueP39Of3HXcT179qxfv362trZOTk5sNjsxMTEoKIhefe/8+fO4KLVDwDDaMWhpaYWHhx8+fNjQ0NDS0nL69OmFhYVMJpOeG41+mpCQsGLFiuvXr0thbA8ePHB3d6eyOlH5VQmCEDRq/dVmzZrFZrMzMjJiY2PDwsJ0m9dc6xS4XK6Li8uSJUtSUlJevHjh7++vgdU3OyycG+0whg8f/uzZs7q6OgAgs6I1qRpP5YHn8/nbtm3z8vLKz8+XwsByc3OpunJ0s2bNojLO5eXlteZSRUVF0dHRYWFh5Ms0MotgdHS0+AYrkpwcqKkBC4uG07Q06NYNlJUhPR2srBoaa2shLQ1oKU+Fe/TokYKCwooVK8jTCeKoP4hkBcNoW2zatOnrr7+mspN5e3s7ODhYW1svWrSIquTD4XBmz54dHh4u3o9uTVpJf3//BQsWSC1ppq6urtBcyBcuXBg9ejR53MqynRkZGbq6us3rJL979456rqfPHkrZgQNw+DAkJoKZGQCApyd8/z1YWYGLC2RkNPRJT4dx4+CzvzXIevTNl17cuXOHSoZSVFSEacA6BHyob4vY2NiCggLq9Pnz5/n5+Vwu986dO1RjXV1dVFSUJD79/v37jx8/pk7//vvvhw8fAsCRI0fIMnlpaWkPHz7ctGnTlStXpJDGbfDgwSEhIWLZDqeqqir0pYqhoeGlRrt27RL9g9rM0hIa7yBFoqqqKvT3wfDhwy82WrZsmRg+CUkehtEOpqys7Pjx4zcbS+WUl5efPHkyIiLizZs3Xl5eZKXMI0eOLFy48O3bt05OTm0sLPEl1qxZk5ub6+LicvLkyYCAgCVLlrTcPzk52d7ePjc3t80sIkUAAAnSSURBVPmXTE1N6+vr4+LimrR36dKlZyNqgZdMfP89lJRASIio12Gz2S9evCgpKWnSrqysTP2kze/KUfuED/VtdPPmTeqGNKPxiY7H4wUGBpLHElqqEhsb+9ft26NUVSE2FgDepKRERET88MMP27Zt4/F4VLfdu3cLBAKqhKREaWhoPHv27M8//3z16pW6ujq5qX/NmjVm5KMvAACsXLmSzPYvEAi2bt1qZWXVPBkrACgrK69evXrOnDkBAQEWFhbx8fFCZ13pwsPDP3z4ILVypAwGBAbCtGkwadLHxpwcaHyvBrS/hJaw2eyxY8dOnz59165dOjo6f//9t4ODQ8vfcurUKT09vcmTJ7dl3EiSMIy2EY/HowIB/U001Ui+CxI7hkBgW1e35PhxOH4cAAYxmeNUVKivZmZmkp9bXl7+Re/HRaSiovL999/TW6hK8SSyahAAHDp0aNasWeQshFA7d+40Nzc/ceJEYWGhhYXFqlWrunbtSlWUAwANDQ1HR0fyOD8///r169IMowBgZweTJwM9xWjPnh/nRl+/hrFjW3WdS5cuBQQE7Ny5s7Kysn///mPGjDE2Nra1taU6GBoaUku7njx58uDBA21tbQyj7RGBvtyECRPI2UCSs7Pz//3f/+Xk5LBYLKqxqKhIQUFB7B/9V1gYAVAP8KZr1wIlJQLgP3JyGzdunDVrlrKy8ogRI8aOHTt27NgePXowmczU1FSxD0AUBQUFVlZWx44dGz9+/KlTp0S/4Lx583JycqZOnSr6pVrD05PYu5cgCKK4mNDXJ/r3Jy5eJBISCCOjj31evSJ69BDz59bW1rq6uiYmJv78889ivjQSB5wb7ZBqGIzxmpr+KipAe3fP5XLT09NTU1NTU1PLy8vJxPvtioaGRnh4+Lhx44yMjPp+dk3Q51y+fJnFYhUUFJSWlpKTwlLTrRts3w6JiS31KS+H1ath0SJISBD143bt2jVixIjU1NSCgoIWEi8gWcGH+g5JWU4u8/vvISEBGt81AYCSktKjR48MDQ0BwN3d/cqVK7IboHCKioqmpqYA8NNPP9FnTttGR0fH2to6Nja2rKysrKxMCi9khg4FNbWGY3d3ePsWzMyge3egbzXS1AR3dwCAlBRYvx6YTPjxRxDxr2LgwIFFRUUvX74sLCzkcrkdpf7zvweG0bYICwujJwq5evUqk8mUk5MrKiqiGrW1tVvOeN82hIJCHpOpz+fDnj1kS4bk38WLHb3WZpuNHj2aXJcqEAhMTExEv+Bnubl9PGYwwNe34Zg+T6qrC+SuiEGD4MkT2L8ffvhB1M/99ttvAaCoqMjKygo3O7VDmP2+g7lz585vy5adnTKFXPeTUlo6JTT06h9/7Nq1S11dff/+/eSzvLu7u4KCgp+fnwzTPiEuF5KT4fBh+M9/ZD0UJEl4N9rBGBsbu6xerdNYmVkhM/NHAwM2m/3NN98YGxtT86HOzs49evTAGCpDd++Cmhp06QKdKB0VEg7vRhGSiJoaOHUKqqvB3R26d5f1aJAkYRhFn/f48eMBAwaoNC5QjYmJYbPZdXV1mZmZZPYQAOBwOGlpaZh+H/0L4YIn9HnTpk3LzMykTmfOnJmamhoTE7O8cW4BAOLi4j67DRShTgnDKOoYcnNz6TvQ8/LyiouLASApKYn+RJWcnCzN7VsIAYZR1FFs2rQpODiYOvXx8Tlz5gyfz+/Xrx+9MHX//v3pFUcQkgJ8U49aZdy4cdRSWSo5U1xcnJGREXlcW1vbs2dP2QwOIZnCu1HUKrdv385s1KtXL7Lxq6++ohovXLgg2xEiJCt4N4o6jF9//dXf3588/vDhA1U0hc1mMxgM8hgnRpH04d0o6jC8vb1TGs2YMYNqT0xMpNqbl+Vo58LCwvbt20ed3rhxw5faZPqp6OjoBQsWDB8+fOrUqUIzDfbu3ZteeprNZn/48CE0NNTDw4NqDA8PX7hwofiGjwAwjKLWWLVqVXfaCvIVK1bo6uqamZnRs3waGxvTiwNLAoPBYDaibj8BQE5OjmqX6AAkITMzM5GWKio7OztBWEqohw8fTp48ecSIEefOnZs7d25gYGCTFd979+599+4dVTPGz88vLS2Nw+Hs3r27sLCQbORyudu2baPXv0FigQ/16PPWr19PP/X09CQPzM3NqUYTExOsHSQ5Bw4cWL58+eLFiwHAzMzM1dW1SYdHjx5NJAjVtWtBUREA6pOTCYLg8XixsbH6+vpkn/r6+qdPn06iJ+5H4oBhFHUMY8aMoVYFAMDo0aP19PTk5OSWLFkiLy9PtS9ZskRRUVEWA2y7169fHzx4kDz+p7oAr169+qHFVFE2hYVfE4RiUBB5upbBON54w56YmLhy5UoA4LWywgn6QhhGUcfQpEzInDlzyINjx47R248cOSK1IYkLn8+nkjHzeDz6bwVKTU1Ny78eTDmcwQSR1KvXSwuLSTEx6pWVanJyAEA++5NzIPSZECRGGEYRkjE2m71u3Try+NixY/fu3Wvex9TU9NWrV58te1dkYPDK3t4hKUm9spJ62yYvL0/mtJZQfTCEYRSh9isjI0NBQUFfX3/GjBl79+51dnYmF+0WFBTo6ek17z/m8eMxjx/TWxgMRp8+fbZu3QoAlZWVu3fvlsrA/13wTT1C7VRhYeG8efPIaYpFixbNmDFj8ODBdnZ25ubmv/zyS5PO93v1yqQtVHjeo0ehhgZ5Q6qlpUXv2R3T9okbJspDSJbq6+sFAgFVl7C+vp7P55PToO7u7jNmzIiOjt6+fTv5VYIg8vLytLS0mpdjcnV1dXV1nddYFmrmzJnjx4+fOnWqjo5OfX092VhZWdm1a1c+n9/hVte2c/hQj5As0Yt6kadkS0RExMCBA42MjKKjo6mvMhgMAwMDodc5c+YM/R3UiRMnFBQUFBUV8/PzqUZVVdXCwkKMoWKHd6MItUerVq2qqqoqLy9PTU29desWPom3ZxhGEWq/3rx5c+XKleYzoahdwTCKEEIiwVkShBASCYZRhBASCYZRhBASCYZRhGSsrKwsJiaGOi0vL3/8+DEAJCcn5+TkUO2vXr3KysqSwfjQ52AYRUjGXrx4Qc+snJycTGZu3bNnT2hoKNXu5+d3+fJlGYwPfQ6GUYQQEgmGUYQQEgluBkVI9pKSkqik1HV1dd26dSOPfXx8qEpNJSUlZKIm1N7g3ShCssdms6lS1fT50F9//ZVq/+6772Q4QtQCDKMIISQSDKMIISQSJs62ICRzGhoatra21Km6uvrgwYMZDIaVlRWVGY/BYFhaWpLZ71G7gqlJEEJIJPhQjxBCIsEwihBCIsEwihBCIsEwihBCIsEwihBCIsEwihBCIvl/EgaoEcsfA6gAAAJLelRYdHJka2l0UEtMIHJka2l0IDIwMjQuMDkuNQAAeJx7v2/tPQYgEABiJgYIUIPSDYxsDB4gBgs6LQChmWB8dgjNjEs9LpoDzRxCNLo9MP0wmhFCMw4Uze0EosQvMQLZDBDMoKZwZYtD0QOhA46z59hsWj3T4bBX536Jtx938zQtdzC9WG6nyN66V/1XucMm7yX2keLT9156Yuwgf71t7xPvFbv0TaId7p6fuz+7/rmNcMNq+5xpq20FXpyxiVbXs5eUWWl/LOfk7v2NlvvD326yn+iy3oZX5sB+m9YcuxX3l9rePRe1/+H0hv2vVJ/ZdqRy2+e7de7vDl5ru/qV2YEUh3N2dixXbUtPZx8wm+C0P6hBw+6Mw9IDhTIb9knPt7R7nrzxQNady3aN57v3L1ux/sC9IkZ7z7lb7MN2rN7/uVjHYefxoD3Fn5/ulzhf4ND8R2fv/zeXHFS2WNoXpzzdq731swP3tCxbv9oHe7P2XnU4sZ3FQS5YfR+P0EKHFct9Dtz2b7Pd0X/CgXmD2IG3vEn7N5w66fCgSejAHP9Uey7D+fZqUgIOi5+H7U3s3Lx/aeqP/eFxSnZCOU32nZXMBzbIP7JdIJ5w4Hz6pv03PMz268yLPXCQYdX+2ac77U0W7Dhwj2Hh/hMxmXZ+Ky8fMJivYr/pRs/+grfzD0T8X2t/4s20/QsSVx1YzuGw70vP3/0/n1488GiijX1o5k77R8rzDmQo7bJffm6zvV7a8gOCf97vmaEv6vB65heHbwn8Ds9eO+0TAwBTa/owBdbRtQAAAw16VFh0TU9MIHJka2l0IDIwMjQuMDkuNQAAeJx9VcuOGzcQvOsr+ANL9LubR3vXSAzDu4C9yT/knv9Hqjlai75EI0CjUbHYXVVN3Ua/frx8++ff8eulL7fbGPQ/77XW+FuJ6PZ99M34/OWPr6/j+f3T548nz29/vb7/HFpYsK/foZ/e375/POHxPHyGUy0bTzJlGa5Bk3gx52OpAGizSITWeOJJKd5LgKTU0gdSN6WohTiIxIR6CYBhy/0BtIEaZ6WF6eAppR60gamZB9DBKLOYyvgi8rILSSxlD2SMV1CaODUTOFUWRbcjhXofwAQlT9VVtn8uze4BOCc+mingaAb6ruoFAlHuOwfboc8CEI9TLDQaubKiNiPrOpphaiRPJyrvvUVpZe4iYO6xN/PFWW7Jl+gE9haVWTSP3Vl2nW7mXHd7oP8mxRI5kDre2ujygswAGLXwzUkGNQ+k9faK5k1QVjcXtMWZHOZ+QtuhJ+/nfFmopiu6KS6W5AMaF9RNUWzvz5qLd81LJI+0cl5Qq6XaXE4eMJGxhslPZF2aqrGHIiuxzKFP+4SAxCn/QkQAzYU0CsKSIXEPk4noob+0UzF9seFbe+tG3kjID/cOJEPUnMG+95rkKSkbSQs5OZACZMzIpQDIjjrvQjnU6Nxdx5+YI1I11/aBkRLdGQ33OtQXAzKmRML/9lbxeTlWpSRH97hvqJIvNN/zXkv6nJi1PHsUfiEDSAyQKfxHoeLsK7ZMjmidJ0MCCUURpNA9dcswfh0+6NF+PKC1SYl7196efGerFSU7wyerSRWWs5lsb6HYPYZQ345ZxntjEfjMPScqxdbzShkuR1fKDXWYjpMo9gnhGPw+8JC9OqcPQwMorEIBfRBOnFIot9cwkZ5taVv1ZHO1Q9bD70tx0jUU2aYj/mpXASqbtX2vdBC0bEUVR1bUrwocbNIipUaADFJbuRwWaNwLMK6iBjjYvW/I0k7OvPYXZZGtJskSnE9IjTLpMVTaZmW3DWE6ATD4MoAXqR/9f3l9+e1Ppb9//B/h/vYfmmZNminUwdAAAAIvelRYdFNNSUxFUyByZGtpdCAyMDI0LjA5LjUAAHicrVM5bhxBDPyKQxuYbfA+0KESRXbg0PBP9HgXuxUrMgZYDJZksQ7O23r78vmJ59d5fq8/7+9//+PPt4/vvsKp+nnJkrb2hxZxM+e2VSTCz4sXpThaUKLU0u1L1ELRKyZEtxSY9i2rmMruP/75Qixlm5dq19lQmjXjJE6sm1ZgV9WDFsHuOxRsufGWgl0xtc7KM8XavkHMiSqGhlLnKRV160yVW3If9nS5MotmYpmbOdenMIszhhbZLwURk6MnA7TPWJg7aj5vHFNT045Zx8WSPDU3tZFmrNnXtBbJKVm1GkpOHrCR0cN06KuxxyMr2ryOaBgVvmN5s81yGG/kUwJ7iN26DJ7FLBAVbBo7CvEMXjZCeBTcJa7zJqKwA6l40TA2dOsQZTJs27mC/brqKSlnFTW8B4vI1gQ/pP5ZCTVSkKg0bAAHgB1LpiUdl0E6RsBJRiIn43AvAZpEus2hKQI4HlcpLgAlJe97g9UyXlU7NuBkxDQKFMSv9eRILEcslCBWvLTRZAgZHooRYoAMGPlkOdaR8UkX12Y2MwphN2bYaj01C8lzFCrF/sxxhOsEmNF+HfesYYdga64MDgFP8E8YJBxkJtLYL1t99APOW8804zAGTmVmxs1K1yuhcMeD52gevNQI4rkVqyEBOOPzgZgbPpCRb2kDJ8py4Eha8AHDX2WElEMKIuEdrDpiuUldf3z8A5ad6dPZXd/5AAAAAElFTkSuQmCC\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x75ee4ffd44a0>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"mol2 = Chem.MolFromXYZBlock(xyz)\n",
"mol2"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "e5208c12-b13e-4c63-8668-a85899eb6122",
"metadata": {
"execution": {
"iopub.execute_input": "2025-02-07T13:11:03.966327Z",
"iopub.status.busy": "2025-02-07T13:11:03.966213Z",
"iopub.status.idle": "2025-02-07T13:11:03.972046Z",
"shell.execute_reply": "2025-02-07T13:11:03.971715Z",
"shell.execute_reply.started": "2025-02-07T13:11:03.966314Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3dd1xTV/sA8CchICSRIaCCbHGAigIqpWhdiK3V2iq43mK1Vtxg3aOKu4gDqYOiWMD3bVXUnxUcWGodSHFEHIgDHGwUAoSVACE5vz+uRkodCPfmBn2+Hz98csPNOU8wPJxzz7nncAghgBBCqKm4bAeAEEItG6ZRhBBqFkyjCCHULJhGEUKoWTCNIoRQs2AaRQihZsE0ihqSSmHTJsjIeH6Ynw9RUWzGg5CGwzSKGqqqguXLYepUoKYUZ2XBTz+xHRNCGgzTKHoFIyMQCCAyku04EGoJMI2iV9u6FVasALGY7TgQ0niYRtGrOTrC11/D4sVsx4GQxuOxHQDSXIGB0K0b9OrFdhwIaTZsjaLXEgph61ZYt47tOBDSbJhG0Zt4e0Pfvs8f5+XBypXw+DGrASGkebBTjxoyMICjR18e/vIL/PUXAMAPP0BUFMhksGULW6EhpIk4uN4oegOlEvr1gytXID0dysvBxQUMDSE3FwQCtiNDSGNgpx69QnExLFsGkyYBlwtdu4JSCbt2gbMzuLuDRAK//sp2fAhpEmyNolcoLQULC5DJ4N49kMnA2RkMDSEnB+LiYOJE6NYNUlOBw2E7SoQ0A7ZG0SsYGcGECUAIhIVBr17g4QESCfzvf+DtDR06QFoaXLjAdogIaQxMo+jV5s4FAPjlFygvf/54xw7g8WDatOePEUIU7NSj1/rkE0hMhF27wM8PbG0hNxfOngUHB7CxAYGg7NYtpaWlEdsxIsQ+bI2i16IaoTt3gpYW+PkBAOzYAWZm4O8fXVtrERa2md3wENIQ2BpFr1VX97wRmpAAPXrAV1+dkstDjx7dW1BQ8NFHH5mYmOTk5Ojq6rIdJkIsw9Yoei0eD2bMACOjktjYo+3agb39QZHoj927d7u5ufXt21csFh88eJDtGBFiH7ZG0ZsUFck6d+5QUVGekZFRWlrq6upqZGSUm5t79OjRSZMmOTs7p6SksB0jQizD1ih6E1NTvS++GKlQKHbv3u3i4uLm5lZaWnrw4MFx48a1a9fuxo0bSUlJbMeIEMswjaK3CAgIAICIiIiqqqq5c+cCwI4dO3R0dKZNm0Y9Zjk+hNiGnXr0du7u7pcvX96zZ8/kyZNtbGzy8/MvXrzYsWNHGxsbQsiTJ08sLCzYjhEh1mBrFL0d1QgNDQ3l8XjfffcdAOzYscPc3HzmzJmrVq0S4Dol6MOGrVH0dnK53NbWNi8v79y5c126dLGyshIKhadPn/7oo4/YDg0h9mFrFL2dtrY21QgNCwszMzOzt7eXSCRLly5lOy6ENAK2RlGjFBQUREVFfffdd6ampvHx8Z9//rmuru6TJ0/atm3LdmgIsQzTKGqKUaNGxcbGrlixYv369WzHghDLMI2iprh69aqbm5u+vn5WVpahoSHb4SDEJrw2ipqib9++gwcPLi8v//nnn9mKYdmyZXv37lUdrlq1ateuXWwF8/75+eefT5w4oTqMiIj4/fffWYxHk2EaRU20bNkyAAgJCZHJZKwEUFxcXFZWpjosKSmpf4iaSSQSpaenqw5TUlLu37/PYjyaDNMoaiJPT093d/fCwsJ9+/axHQtCbMINllHTLVq0aPTo0cHBwX5+fjo6OuoP4Pfff8/MzKQeX7x4cfz48eqPoWUoLYXoaEhNBSMjGDAARo588+lUD0MkEkVHR1PPPHjwwMrKivE4WyZMo6jpvvzyy+7du9+5c+fAgQPffPON+gPo3LnzsGHDqMeqfIoaKisDd3cYOBAmTACxGBYvBpEI1qx5eUJhYfG9e/+7eTMtLS01NTUtLc3Pzw8AxGLxo0ePqFPKy8tZib1FwDSKmo7D4SxevHjSpElBQUG+vr5crrqvETk6Oo580bA6c+aMmmtvMcLCwN4eVIOBffqAkxPMmQMbN8Lt25CaCkVFrW1t5z15onpFbm6uUCj89NNP58+fTz0jFovVH3hLgddGUbNMnDixY8eO9+/fP3bsGNuxoNe4cQM+++zlYceOYGcHd+7AyZPw119QVASGhtqWlrNmzNi5c+e5c+dwQe53ha1R1CxaWlrz58+fPXv2+vXrR48ezVHj7vVt2rTR19dXHRoZGdU/RC9JJNDgJ2NgAKWlsGkT6OpCt25gZcUBwMliTYbT71Fz1dTU2NnZ5efnx8fHq65UMi0sLIzL5U6bNk39VxJanunToUMHWLXq+SEhYGEBJ06As/MbXlReXs7j8fh8vupQS0sLV/N6JfwIouZq1aoVtbTzjz/+qJ4ac3NzFy9ePGPGDFx7v1FGjoR9+6C4+Pnhf/8Lxsbg5PTmF+nr66tyKHWIOfR1sDWKaFBVVWVtbV1cXHzp0iUPDw+mq/Px8Tly5Ii3t/fhw4eZrus9sXw5/Pe/MGAAlJRAWhocPQq9e7Md0/sD0yiix6pVq9atWzdixIi4uDhGK0pISPDy8uLz+Xfv3rW2tma0rvdKSQncuQOGhuDgANrabEfzXsFOPaJHQECAUCg8ceIEo3uF1tbWUkvxr169GnPou9myBby94epVzKG0wzSK6GFsbExtcrd582bmagkODn7w4IGjo+O8efOYq+X9VFgIRUWgxqkUHw5Mo2py/fr1nJwc1eGNGzeys7NZjIcJCxYsaNWq1eHDhzMyMpgoPzs7OygoCAB27NihjU2qd1VaCgCAqxoyANOomvz4449//vmn6nDz5s2nT59mMR4mdOjQYdKkSQqFgqEG6Zw5c6qqqv7zn/8MHjyYifLfcxIJAICREdtxvIcwjapDREREKdUWAACAyMjIYtXsk/fL0qVLeTxeVFQU7W3t+Pj4uLg4fX394OBgekv+UFCfQEyjDMA0qg6xsbFd8vJ6HzkCS5fC0qUZ+/e/rytj2tnZeXt7y+Xy7du301isTCabPXs2AKxbt87c3JzGkmlHSF1tbZZc/kz1jFR6A0AD5sNgp54xeDOoOjiJxd75+T0ePIBTpwBgOYfzG48HU6awHRcjli1bdujQobCwsG7dullYWLRu3VpfX5/6amho2LS7RTdu3Pj48eMePXrMmjWL9oBpJJVef/TIW0eng1IpEwo/trTcAQBi8R5Lyx0cDtu/a9ipZwzb/7UfBofS0l4VFQ+trMDLy+roUWFpqYeDA9tBMcXJycnGxkapVFJ7Mjegq6urp6enq6trVE+DZxoclpaWbtmyhcPh7Ny5k8fT6E9sfv6adu2+b9vWHwAIqQOAgoI1VVVX8/NXGRp+IRB8xFpkSiWUlwOX2/DmekQHjf5Qvmf+BrhIyBqADmxHwqiYmJgnT57o6emNHDmyurq6oqKioqKivLy8oqJCIpFUV1dXV1cDQEFBQePL5HA4EydO/OSTTxiLmh6E1FZX3wVQAnCp5qeZWaBc/tTcfC3LrdGyMlAqwcgIcAkCBmAaVZ/xYvGkc+egqortQBgkFoup6fGhoaHUNNIGqHyq+iqRSFSHVJ4tLy9XHVInSCQSADA2Nlb3m3l3FhabHj8en5pqa2w8uV27hVparQHAwmIz6z16hUQi19XVNTBgN4z3FkHMm+XhUaWrSwCof/dMTBbOmnXz5k2246Kft7c3AAwePFipVObn59NS5rVr1zgcjp6eXk5ODi0FMq2yMjkj47P0dE+2A3mJurXsIzc3tgN5P2ELXx3isrP7Ghqm/v47PHoEjx6t6dfPwdW1Z8+ebMdFs5iYmCNHjggEgr1792ZmZnbu3NnPz0+hUDSz2N69e48ZM0Ymk61fv56WOJkmEHxkabmzouICACkpOUBIc38CzVdSUgIAAqGQ7UDeT5hGGXfo0KG8vLyMkhJdR0ewswM7u43bto0aNYrtuGgmFov9/f0BYNu2bba2ttOnT6+srKyqqtLS0mpymceOHXNycnr48OGGDRt4PF5ERISG7/ErFu+prLxUXZ1eVLRTIHArLPzpyZOJGRlD5PJ3uBDMBGrashEO0zMD0yizFArFd999p1Qqx4wZ06lTJ+pJW1vbFnGl753MmjXr2bNngwcPnjZt2u7duxMSEkxNTZs5e/TkyZOpqakrV67s3LnzlClTFApFYGAgXQHTiBBFXt6KuroiDqdVYeFP2dkzCKmzszvE5ztra5tVVFy4e7dXRcVZFiOkri8b4qRRhrB9VeE9R0364XK5BQUFbMfCoN9//x0ABALBo0ePnjx50rp1awA4cuRIM4vNy8vj8/kcDuf69euqx5cvX6YlZhrl5i4TieDu3d7//pZcXpie7ikSwfXrvLy8QEIUao+OEEKoW78WLVrESu3vPWyN/kNVVVV6errqUCqVPnjwoMmlFRYWRkZGAoCfn1/79u1piE8jicXi6dOnA8DWrVttbW39/PwqKiomTJgwZsyYZpZsbm4+a9YsQsjy5cvNzc3nzJlDCFm5ciUdUdNGIjn+9GkQh8OzsNhCSI1E8n/1v8vjmdrbx5uZBRKiLChY8/DhlwpF6euKYg526pnFchrXMBcvXuzZs6fq8PLlyw4ODk0uzd3dHQAEAkFdXR0d0WmosWPHwovR+V27dgGAqanps2fPaCm8pKSE+uU/e/ZsaWlpmzZtAODPP/+kpfDmk8nu37hhIBLBs2fbCSGZmX4iEeTlrfr3mRJJ3M2bbUQiuH3burLyiprjnDlzJgDs3r1bzfV+ILA1+lJOTs6pU6dUh/n5+bGxsU0u7dKlS8nJyQAQEhLSnGEWDXf8+PGYmBiBQLBnz56srKylS5cCQFhYWNu2bWkp38jIaOHChQCwdOlSAwODRYsWUY+JBuzaoFRWPno0WqEoa9NmQtu2AcXF/xWL93C5uoaGX/z7ZAODEQ4ONwQCt9rarPT0AYWFoeoMFVujzGI7j2uQ8+fPd7W3X2thQcLDSXj4g/Xr7ezsmtwapVbQ6NixI71BahSxWNyuXTsA2L17t1Kp9PLyAoDx48fTW4tUKu3QoQMAHD16VCqVWlhYAMDhw4fpreXdKR898hGJIC2th0JRWVWVkpKiJxKBWPzLm16jrM7O9heJQCSCJ098FYpK9cRK7dgaHx+vnuo+NJhGX7pw9uxqU1PVJHkCsIjLdXBwyM/PP378+DuNEYWGhgIAh8O5du0acwGzbty4cQAwaNAgpVK5e/duADAxMaGrO18fVXiXLl3kcvnPP/8MAJ07d5bL5bRX1HgFBRtFIrhxw6i6OkMuF6em2opEkJ09tzGvFYv3p6QIRCJIS3PKzk5nNM47d+4sWbKEz+fr6enNnz9foWBnjOv9hmn0pYtnzhCAOg6n2Nu7sndvAnC4dWsHB4d9+/ZRLXczM7MRI0YEBgYmJCRUVVW9rpza2lo9PT0A8PLyUmf8anb8+HEAEAgEDx8+zMzMpEbnGWok1tbW2tvbA0BERERdXV3Xrl0BYO/evUzU1Rjl5X9ev64lEnElkhNKZV16updIBPfuuSuVNY0sQSa7l5bWbc2aAa1btz548CDtET569Gjt2rXUD4rC5XIBwN3d/e7du7RX94HDNPoSlUalAJ9//nmIkxOVRm1tbb28vBwcHIT/vANEW1u7d+/es2fPjo6Ovn//vlKpVJVDtdF4PF5JSQmLb+c5hYI8fUp7qaru/K5du1Td+XHjxtFekcpvv/0GAB06dJBKpQcPHgQAc3PzN/wxY05NTdbNmyYiEeTnryGE5OQsFIng1q12tbW571ROXV25r+9/qI/TvHnzamtrmx9bSUlJdHS0p6enakHCNm3a+Pn5JSYmnjx50tLSkvroLlmypKamsRkfvRWm0ZeoNKoAIJ6exNGRSqP1h0r09PQ6derUvXt3KyurBqNGRkZGEwYN+nnatLjVq/twuUYAixcvZvn9SKVkyhTSoQPp04d06EDmzSP09YLHjx+v6s6HhYUx151XUSqVLi4uALBlyxalUunm5gYAwcHBzNX4SgqF7O7d3iIRZGSMIERRWnpMJOJcv65dUXGhaQWGh4fr6OgAQO/evZ88edK0QmQyWWxsrI+Pj2qLKj09PR8fn9jYWCo7U1/Lysr8/f2pZqmTk9P7fcVJnTCNvnT+r7/+FghUF0aVXG6AmZm/v39ISMiECRPs7OwajM6Zmpo6OTm5uLhYW1u3Bwitd1H1HpdLSktZfj/+/mT4cEK110pKiJsb2biRloJf2Z2PiYmhpfA3oHavMjIyKikpOXPmDAAYGhoWFxczXW99588vFYkgNdW+rq5UJrt/44a+SATPnu1oTpkikcjW1pb6U3T69OnGv1ChUCQmJvr5+VH/BQCgpaXl6ekZHR1dUVFBCBGLxeHh4R4eHlOnTlW96uzZs9SHmcfjrV25kshkzQkeEUyj9Z0/f/6jnj1JQgL171ZUVLdu3eqfUFhYeOLEiVWrVg0bNqzB3JGlXC4BKOLxRABVVDK9fp2tN/KcsTGpv4jU6dPE0bH5pYrFYupWgp07dyqVSmoIeOzYsc0vuTGozexWrFhBCBkyZAgALFu2TD1VkxcjXVu2DJZKb1dWlt2540ANuDe/ZLFY/Nlnn1HDkkuWLHnrRGNq4Kj+PR2urq7bt29/+vQpIaSsrCwqKmro0KGqPpO1tXX9wSWpVLpkyRItLa2QAQOIvT05d675b+FDhmn0peLi4kuXLqkOS0tLL168+Ibz8/LyYmJi/P39PTw8lvN4BGATgFAovA5AAC799BPzIb9eaSkB+EdDIyuL6OqSoiLi70+2bSP/938kJaUJTeYJEyYAwMCBA5VKJTVobmJi8pSBy6+vdPXqVQ6HIxAICgoK1LyAXnJycqtWrQAgOjqaWiRh+nSn1NT+CoWUlvKVSmVQUBCV+AYNGvTKH2lmZmZQUFDnzp1V2bNr166BgYHp6emEkLq6uoSEBF9fX4FAQH23VatWI0aMiI6OfuVF5JTLl+u6dycAhMslAQGkUk2zr94/mEbpId+wgQA8FArDAYoACIAPuzNGq6tJgwsL9+4RY2OSnFx/RhcBIIaGpFcv8uWXZN48Ehp69uTJ27dvl5eXv7JU6n4EPp+fkZGRm5tLLXVx6NAhNb0pQgghX331FQBQN4lS95vOmDGD6UqfPn1KzV2dP38+IWTdunXU5YVHjx7RW9G5c+eoNqaFhUVSUhL1ZHFxMdU3Vw0cmZub+/v7JyYmUieIRCJ/f3/VdXwul+vh4bF9+/aioqK31FdbS4KCiI4OASA2NuSPP+h9Ox8ITKM0iY0lXK4qNxUDmAOwfAm/c2cSG/vycM8eMmAAyckh27YRf38yciTp3p0IhQ2yausXExJMTEx69+7t7e29cOHCnTt3njx5MikpiRqd37Fjh1Kp/PTTTwFg1KhRan5b9+/f5/F42traGRkZDx484PF4PB7v3r17zNUol8sHDBgAAB4eHrW1tQkJCVpaWlwu99SpU0xUl5ub269fP2pIfdKkSZ9++qlqByp9ff3JkycnJCRQvf60tLTAwMCOHTuqGqeOjo6BgYHvnNxv3yZ9+jz/DPj4EPVebn4PYBqlz5kzJDx8S+fOfgDdAACAWkqDNdHRxNaWnD1LCgtJXBwxMyOvvImlqIhcu0ZiYkhwsDQgYPjw4Y6OjtS813/T1tbu27evQqFISkricDhMj86/zrfffgsAEyZMIC/W0GJ0rlVAQAAAmJmZ5efnZ2ZmmpiYAMD69euZq1EulwcGBnI4HOpeONXAUWVlJSEkNzd3+/btHh4eqv8XCwsLf3//6825HC+Xkx9/JNQeDZ06EbmcFBSQRYvIkCFk+HCyfj1hY25ZS4FplGbx8fGqD7e/v7/q+dTUVFomBr6bo0fJmDHExYWMHUvecTmPgoKC5OTkAwcObNy40c/Pz8vLq3379lwuV3VPwfnz5+Pi4hgI+u1yc3P19PQaLKAnEomYqIuar6qtrX3p0iWZTObq6goAX3zxRf2ZwkxQKBTUOiwrVqygZiNIJJLo6OgRI0aoGqeGhoa+vr4JCQm0BfPwIRk0iGzfTioqSKdOZNEikpFBUlPJmDHE05MoleTvv8nChWThQkL31YwWDdMozRQKhWpJZkNDw+rqaup5AwOD7OxsNQVx/z4ZMYLQPfBSWFhITaxJTk6mt+QmWLBgAQAMHz6cEEKtV9K3b9+LFy9mZGRIpfSM+RBCbt++TQ3XhIeHE0ImT54MAJ06dZJIJHRV8Trnz58HADs7O0JITU3N6NGjqQEuANDT0xs7duzx48cZmUKvUBCFgkREkH79Xj4pk5H27UlyMqmrI1IpkUoJ3lRaD6ZR+i1YsID7Yhtb1c2R6kujNTXE1ZUAkJkzaS972bJlADB06FDaS35XYrF43LhxqamphJDMzEw+n09dt6Xo6ura2dl5eHj4+Pj4+/sHBQVFR0cnJCTcuXOHmlDZGCUlJdQdqL6+voSQn376CQCEQuGdO3cYfGMvUBusUmtZEUL69OlDDRyFh4eXlZUxXv28eWT58n8888UXhL27bzUcplH6nTt3jpoJBADDhg27ceNGVlaW+tLokiUEgNjaEgZ+2VQrfp49e5b2wpuMulTavn17Dw8PGxsbVavtddq2bdujR4/hw4dPmTJl5cqVu3btOn78+NWrV/Py8lTLnSgUis8//xwAevXqJZVK//77bx0dHQ6Ho55pCUqlkrpxU3WlIiUlJS8vTw1VPzdjBlm9+h/PjBtHQkPVF0CLgmmUfj4+PhEREaoGKZ/Pb9OmDYfDUUcavXiRaGkRHo+8mCtDO2p7zr59+zJ9cbCRqFuqdHV16zcSi4qKUlNTT506FRkZuXbt2tmzZ48aNcrNza1Dhw6qC4uvxOVyzczMXF1du3XrBgAmJiaZmZmq2U5q24Tj0qVLAGBjY8PaDzk4mEyY8I9nevQgJ0+yE4zGe9NHCjUHj8erra0FAHt7+xkzZlDdYWZJJODrCwoFBAbCxx8zVMm8efN27Nhx9erVEydOjBw5kqFaGkm1f0lQUBCV+CgmJiYmJibdu3d/5atKS0vz8/MLCgr+/TUrK6ugoKCgoAAAhELh999/b21tfe/ePT09vcGDB2/cuFE97+vIkSMA4OPjo5ooqm5jxsD69XD5Mnz0EQDA/v1QUgKDBrETjOZjO4+/h6jWqJ6eHjWhr1WrVrt27VJHp378eAJA+vQhDE8JCAkJAYAePXqwvnilt7c3AAwePJiuSORyeW5u7uXLl6dMmQIATk5OVMmlpaVvn8pOE6VSaWVlBQAsb94XG0vs7Ejv3qRHD9KzJ0lJYTMYzYZplH4+Pj7GxsbUtErqb9XcuXMZT6ORkQSACIUkndllgAkhNTU11FIav/32G9N1vUFUVBQAGBgYZGZm0l64TCajrk6q+R4tQgi194ylpaVGXDYpLGR/kR2Nh3sxMWLTpk0BAQGPHz+mpk/HxcVNmzZNtQwP/R4/hoAAAIBdu6BTJ6ZqeUFHR2f58uUAEBgYKJfLma7ulfLy8r7//nsACA0Ntba2pr18XV3dH374AQB++OGHuro62st/g6NHjwK7Pfr6TE0Bd7d/K7bz+HuI6tRTjyMiIgCAw+EcO3aMoen3crn8v19/rTQwIEzeydMAu0vQq1aWYvRWVLbeI7WKXRJjg4SIdphG6Vc/jVZWVnK5XOreSj6f7+npuX379iavzvtKq1atAgDvPn2IehfbP3DgALxYjl6d9RJCtm/fDgCmpqZMryyl/vd49epVqkbWrzujxsM0Sr9jx47dunWLekyNI7dt27bBqHGXLl0CAgJOnz7dzN/PpKQkapkM9U/kVCqVzs7OABASEqLOeqlxcwA4cuQI03Wp3uO2bduYrouyZMkS+OdtxEjzYRpl0F9//cXhcFq1akVNaXz27FlMTIyvr2/9JZ91dXU9PT2DgoKacEt4RUVFp06dQL1LF9cXFxcHACYmJq9bWI92crmc2j7k22+/VU+NJ06coN6jOu4dIoS6b+rNC90iTYNplCmVlZXUhKeN/9q6o66uTiQSBQUFeXp61p8N3r59e19f35iYmNLGjY36+voCgIuLC4vbk33yyScAsG7dOvVUt3r1agCwsLBo5I+IFv369XNzG7V587vtWNcEKSkp1MfgravfI42CaZQps2fPBgBnZ+c3jyyJxeKYmBg/Pz/qPhmKlpaWq6trYGCgSCR63ayXw4cPA4BAILh//z4z76BRLl68CAAGBgZq2BMpJSVFW1uby+X+9ddfTNdVX1KSDIDo6xOmZ45S8x9mz57NbDWIbphGGXHp0iUul6ujo3P79u3Gv+rOnTtUE5XaKpJiamrq4+MTHh6en5+vOjMnJ4e6t33Pnj0MhP9uhg4dCvUW0WBIdXU1dX2ZWn9ezT77jACQBQuYrYWaGKDmPxKo+TCN0q+qqoq6wtXkrm5lZWVCQoK/v3/9GZFcLtfV1XXJkiV//PHHwIEDAeDLL7+kN/KmUe2JlJvLYLd33rx5AODg4KD+iQGEkFu3CJdLdHUJc7dQ3Lp1i/qrKadvH2ykHphG6efv7w8APXv2pGWiaFpa2tatW728vHR1dVUplcfjtW3bVm23J74VtT8Sc6v9X7x4kcvl8ni8q1evMlTFW40dSwCInx9T5VMT19SwrxSiHYcQ0ohJ+qixkpOT+/fvz+Fwrly54uLiQmPJUqn0woUL8fHx0dHRZWVl9vb2GRkZNJbfHGlpaT179uRyuffv36dmj0+dOjUxMVFfX19fX5/P5wsEAkNDQ4FAIBAIhEKhoaEh9WSDE4RCoba2doPCy8vLnZycsrKy1qxZQ+UaVmRkgKMjEAJ37kDXrvSX361bt7t37yYkJHh6etJfOmIU23n8vSKVSqkZSIGBgczV8ssvvwAAl8tV2zSjxqCmDUyZMoU6HNTU1YC0tbWNjIw+/vhjVcnUmvMuLi4s7MLyT1OnEoCGC8i9K6VSmZ+ff/ny5cOHDx89epR68s6dOwBgbGyMPfqWCFujdFqwYMG2bducnJyuXbtWf5iIXrW1tXw+X6FQbN26df78+QzV8q4yMzO7dOmiUChSU1MdHGEnYeEAAAs1SURBVByePn1aVlZWWVlZVlZWVVUllUrLysoqKiqox6WlpVVVVVVVVaoTqqqqysvLKyoqqBvYu3bteu/ePQCIjY0dNWqUrq6uSCSqvxQeK/LyoFMnqK6G69fB2fktJ1dUQHY2ZGdDfv6vmZn3s7KysrOzc3JycnNzqRUUAaB79+6pqakAsGbNmtWrV3/33Xd79+5l+l0g2uF6o7S5fPlyaGgoj8eLjIxkLocCgI6OTo8ePW7evBkWFqY5adTGxmbq1KlhYWFff/319OnThUIhn8+n+u9mZmZ8Pp/qv2tpab25nJqaGqlUWlVVBQBisdjPzw8AgoODWc+hANChA8ycCdu2wcmT/0ijT57A339DTg7k5EB2NmRlQU4OSCTPv9ut2460tCv1y2nXrp2lpaWlpaWDgwP1DLUcCbXuH2pxsDVKm/Ly8kWLFrVr127t2rVM1xUaGkqNXKelpTk6OjJdXSOlpKQMHDhQqVRSSfCVdHV1+Xw+dRmUulSqumbaunVrfX196rGBgYFQKAwJCUlISBgyZEhCQoJGLHcEIBbDhQvg4gK2ts+fefAATp+G779veKaeHtjYgKUluLj8rKv7zNramkqdVlZW9UcL6+rqLl26NGjQIENDw2fPnjH6BxgxBNNo082ZM8fPz8/JyYk6DAgI+Oabb+gdVnqdnJwcKysrLS2thISEJl+FpN3EiRMPHDjg4ODQv39/qv9eVVUlkUhU/XdqQ83GF9iuXbvy8vIHDx5QS39qiG++gbg4uH8f2rYFAPDwgAkTICkJLC3B0hKsrZ8/eLHY7HMymaygoODx48ePHz+mVtqnHmdnZ9fV1bVv397a2vry5cusvCPUTNipb7rk5OT6vbArV65Qm6CpgaWlpbW1tUQiUW2HmZCQYGJi4vzWK3aMSUpKOnjwoJ6e3qlTp2xsbF53mkwmoy6DlpeXU513iURSWVkplUqpPEtdPC0vLxeLxYmJiQCgupKoOaytYeFC2L//+aGjI8yZAwBQWwu5uZCfD+fOwePHkJ8PBQXw+DG0bh13/vwXryxKS0urffv2RUVFYrE4OTnZ3d1dXW8C0QbTaItUVFRkb28/atQoqkdfXFwcFhbWv39/ttKoUqmcN28eIWTp0qVvyKEAoKenp6enZ9KgqfYa3377bWRkZHBwcHh4OD2B0mTWLNi6Fc6ehSFDAADkcujbF3Jy4OnTV5/fubO1gYGBpaWljY2N5QvUY3Nzc21t7WXLlgUFBX399dc3b95kcHlvxAxMo81y7NgxaqQVAJ49e6a2es+ePfvo7t12rq5w/ToAXBOJrly50r9/f7UF0EB4eLhIJLKyslq4cCGNxa5YsWL//v2RkZHLly9nYon7JtPWhpAQmD0bbt16fpieDmVloK0NJiZgbg52dmBnB2ZmqsdORkaSNxS4du3aP//8UyQSLVq0SLU7N2opMI02S21tbU1NDfVYqVSqrV6uXD6wpmZscDAEBwOAp46Oq0CgttobKC0tpWbFh4SE8Pl8Gkvu2LGjj4/PwYMHt23bFhoaSmPJzffZZ7B3L4SEPD+8cAFMTcHMDJo2EqatrR0dHd27d+/w8PARI0aMGDGCxlAR49ictNrCubi4nDt3TnXo5uZ25swZ9VQdFxpKAGo4nAetWxe1akUAgrS01La0cAPUWlaDBw9movC0tDQul6urq1t/ZRZ2TZpEIiMJISQri5iZkU6dCF1LZm/duhUA2rZty/Sq/oheuKVdC5bP4Qw1MorU0wMAakF49UtLSwsPD+fxeNTGHrRzdHQcOXJkdXW1prVGAcDKCgICgMY7cufNmzd48ODCwkJq0wTUUmAabcHMWrXKmjBh0YsZV6yYM2dOXV3d7Nmze/TowVAVq1at4nA4O3fuFIvFDFXxTgYOhC5dnj/+/nsICIB6S8U2C5fLjY6ONjIyOn78OLV9NGoZ2G4Ot2Aymaz+vmPV1dVq24bs6N69xVwuAVD9m29oqP5O/aFDhwDA2NhYLBYzWhG1DyijKxVojv379wOAUCh8+PAh27GgRsHp9y3SwYMHT27aFOrlRS3efO3p0zl///3rr79S65yqh0wmc3R0zMzM3LNnz7Rp0xitKzk5+eOPPzYwMMjKyjIwMGC0Lk0wbty4mJiYfv36nT9//q23zyL2sZ3HUVPcunUrKipKdZiWlrZv3z41x0CNzjs7O6tn4yBq06egoCA11MW6kpISCwsLAAgODmY7FvR22BpFTaFUKgcMGJCUlJSYmOjh4aGGGv/4449hw4aZmJhkZmYK2JvdpTYJCQnDhg3T0dG5cuVKz5492Q4HvQkOMaHGSkxMpFaxAwAul7thw4b4+Hj15FAA8PLycnd3F4vFERER6qmRXUOHDp05c2ZNTc3EiROrq6vZDge9CaZR1FhDhgyRSF7eijN8+HBqjzm1Wbp0KQBQG8d/CDZv3tylS5e7d++yuOY/agxMo6jFGDlyZGRkJLWzNCUnJ6esrIzFkBjF5/N//fVXbW3trVu3nj9/nu1w0GthGkUtBofDiYmJiY+PVz0zd+7c48ePsxgS01xdXZcvX65UKn19fUtLS9kOB70aplH0Dpydna1fkEqlbIfzQfjhhx/c3Nxyc3M1Z6cD1ACmUfQObty4kfUCvauQoNfh8Xi//vqrUCiMioqqf0EDaQ5Mo6iFCQgIsHshISGB7XDUoWPHjkFBQQAwa9asp69b0xSxB9MoamG2b9+e8cLQoUPZDkdNZs2aNXz4cLFYPHny5DfM9a6oqKh/J1t1dbWtra06l3D8MGEaRY21bNmy+h35xYsXs7JOO4fD0XpB/bWzhcPh7N2719jY+MyZM3v27HndaatXr37y5InqcO3atZmZmXiLDdNw2WbUWGvWrKl/iJMZ1czc3Hzv3r2jR49esGDBoEGDOnfu/O9zkpKSvAmBF+vsSXGPPLXANIpakuHDh9fvtA4bNuyV2eR99dVXX02cOPG3335zdHSkGuNcLpfD4XBfmGZoOIkQeNFcDeJw/stqwB8ITKOoJZlD7cD5wsyZM9mKhC329vYcDkehUCgUin9/10lHpwfAzU6dHlhafpGUpFdTo6v+ED88mEYRajHS0tKCg4MBYMOGDfb29iUlJUVFRQUFBdRXPz8/rSVLAOCZre09d/fh167Bi43CEKNwhSeEWoba2to+ffrcvn177ty5P/300yvPOWBuPqGgoP4zHQCy6+o+qOE49cOReoRahsDAwNu3b3fs2HHjxo2vO+d427bPuC9/qePbt69t04bTtN1KUaNhGkWoBUhOTt68eTOXy42KihIKha87LVcoTDlxAgih/gU7OET88guXi7/mzMJOPUKaTiqVOjs7p6enL1++fMOGDW84UyKR8Pl8HR2dVx4ihmAaRUjTzZkzZ9euXT179rx69SrmRA2EaRQhjXb27NmhQ4dqa2tfu3bNidXNtNHr4EUThDRXWVnZt99+SwhZs2YN5lCNha1RhDTXN998s3//fnd398TERJy0pLEwjSKkoWJjY0eNGsXn82/cuPFB3fPa4mCnHiENdfLkSQDYtGkT5lANh61RhDRCYWFhXl6es7MzdSgWi7OysoqKioYNG4bz5zUctkYR0giJiYlLlixRHSYnJy9YsODTTz/FHKr5MI0ihFCzYBpFCKFmwYXyENIUiYmJ1tbW1GOZTObo6MhuPKiRsDWKkKbo37+/av/qffv2sR0OaixMowgh1CyYRhFCqFm0Vq9ezXYMCCHgcDimpqY9e/ZUHRobG/fq1YvdqFBj4PR7hBBqFuzUI4RQs2AaRQihZsE0ihBCzYJpFCGEmgXTKEIINQumUYQQapb/B5Z8J9lYwHMwAAAC8XpUWHRyZGtpdFBLTCByZGtpdCAyMDI0LjA5LjUAAHicxZFrSJNhFMef9+I2N+fmzLvpq84LJlpkgrC9z5NEVIaCklZWzNLyBpkS0odIDaWhUH7RxDtlmw7RppG37XkiKE1R1CShNBMv4QUJJSoMe71+6kPfOvDnB4dzzv8PZ9XaMgWEUgqiwU4FCAoUlE+JOE4gK+KCBDD7VHLBAimG2WuId8iIdLuD/0jJzh697/NXbM3Qux70vqdkO8QmQ+0dobZJUf+LMkqIR6kZhmYZlpHQNixnQwXRIpYTi9JoiVhnK0mjpbY6mTSNlol0drac3A7Yy4HCHijtgYOEUzkEMY50ME3TTswBRzXr5KhmnAFwAcAVADcxcJcCDxnwlIODcuBlD7wVgFMAHwXwVQI/JVArgb8TiNp6og1FM6yNSCyxlcpELqNCLgrsCARw7ztQzrSKHK+q1phNlej1GT12XfnWY1doQMdG8ngfcZEl8FceMkc/hRdcKiyjs0eR94diy2y0sTssPAlNDtfgzHsLGsd8E8wqN2mVXwc1SYGh0M2zCb7J6u/BBZE4ccUMy060auSeBGuKsnjj50bt5NBF/KUiHy/6z2sfpMrgzZN6XBLfojUtRpAUNMTz7Lj29kAmiXgUhePyg/hB1EhuebZZPeoi+YVrz0nGpzG+YLgEPzO2kqkcCp6u6YAJnSa8lhuCut7G9eauzWHX4Wx0fyPEsrk8itQdkTA3Zc5y6MUakpVnaGPuTlsyLOOo7yWLvOIDrXaqBmQ0nCUfY4u1nQ/7ENPmTFbkV3Hbu340Xagi1bGpUHqkDga4K9GThQRLsr4dN6b+wIlXfHlVViHU32FIm/eMtt5FR4ZvmPHEqQgcUnuZvALNuGpAD8PrO8kUaMB9l9L5mKYxcrhODc0TpTh7pY6c32yBfcvluD65mRgkyLpe+hv/nBshM2UaeC69C8741ZI0325oGGqHodcNxGFjtfdxmBNaqlxH33UKNL8UZXX+A5pxBZX7OI1+AAAD8npUWHRNT0wgcmRraXQgMjAyNC4wOS41AAB4nH1Wy3IcNwy86yvmBzSFJwkecrAkl51KWapKFP9D7vn/SjdHGtKX7O6sZnqbIIAGQD0cfP358sc//x73y18eHo5D/uczxjh+uog8/Dh4czx9/fb76/H8/uXpE3l++/v1/a/D6/CBNXj/yv3y/vbjE9Hj+cizpdSI49FOG4H3IafoUO1rqYEYZ4mZjONRT+mWXAKmdC9fTJ8mzaNZwpCFCZeA2GJkLmIccPKsHi380NPKs8kkdu8dxOP5+7ff9JOesGtnqVToZS4rLr6oVSzD7XiF4bAU2oNltyGNQVnB60XsMKmn+6iYP5d3RgJeim4hFXhyNkRfxQWG1Hzs3DS2LA0QAXeL5o3M0avVtKg+tthVyNQzRSq5t7mM3qcT0HjbW/WyWRldr9QLrDO1quZ9211t+pkRqfUhElSYRrHENqYfb5S7spBsEEKYftqUQDY3ZnB7R/BhcIvBNZnJObVF5k6lQo9JXC8hPXw0BqWl1nWjtoua4XCW+6v3odPnYda3mtV+UaOGO22lZIOIijUquTPryqmHZnPUShuRyA91QoG0Pf0DJQJqH6hJQ7H0Zu2jmMLMt/wblWpnDg08UdsMSTKRfqi3MRVJ7WfTnHudkt26TaYM1MlHQT/eFW2GBe1sfTh4Nutep7/aPGR3wo/vaCpxj3TKoSgWn6XaMmsTwQLMdlrrKANK7Ph7CVflYlsScE+qSw7kgM1fw5zVWiM7O+JmNjDRR+EoAzhqqTnazFaiwvYx0cFEYlFPzWfzjUAXsgaRFsqyqDWNinJXbi85S4yJldhr0AaNOpTXCJsSI2Mf1QgRYmtpfCYXdd/7bBe30mDbSm9pW1SupCa0x1hqc1Ak+p/TDyVYexOid0CFVHCAU/HEyIK7XKMivofllOoxzkGFgjMgh2PskYoSl60LPC4H3KZV6l49YYBpK6m2zTTPy4OENWOSurcGY0h1VNomgbcPB0KrhISE9eSNRI/dZr/2N1ezmU2xYRhTqBpX8a23nGJ1ho3EsAIg8CWADvHc4v/6+vLLCXOdOU9vry/rzFEcJ6sL8IhuXceH4mTwdUgYJn+sQ0Ax3nNNeuO15rniun8trK6jr+EcuGobwcGvsU3a4JfuAzX4pbezGHwTuY2irZWo7mNQJ3zHgGmnE77DUHpK+PYVw0snvJBBpJZ/NpMD9A7BmEjC20TRie6TQwnbPiF0wtsc0Ilu7a4T3bpaia6N6S3RRaG3hJeQVJrwHbVT+Zmzra90wnfUTn8Jt61LdMILocOEF0KPCd+6OF0mvBD6THgh9Bmw3T6zePdS5fPnv1a4f/gPwdrOh5XFJhUAAAJdelRYdFNNSUxFUyByZGtpdCAyMDI0LjA5LjUAAHicVVO7bttAEPyVAGlkRGT2/QDhSo0ru0ipqPJv5OMze1QAp6HEG+7szs7c/e3x8fN2uX9sj5fXG/++//rxeL1dPuXzcn97vHxebt/f8fy4rdf1eP79cvLl9N/rs1zm+Q7i/+Fvfy65h2n2VXYxkr5utHOTuh6xR7YmEEpNWUio0SDebHqlPcWNfEEsRHWAjr0SEHnKs4o62A7fRS2mSgzfrlYU1u6H7pUWduVdSj0Wgp7ph+zFVHZ+688/xFJ28K7a5dOqNGuISZxYD9qjSKrAp8KcZxFmSBSJadTIdT4buVjksY0YjBdT1HlKENb2Y8OJMcaSPdq8VhVGiAVlQ8NV9wyJczwTUR3IsZEYvUqdi7CoGxBUedGcYOdQfJ0lGTiPTTG7nT5kQOnI4jD3hUGeocmMOKbNR9i79WAWkrUwKXbIoAzYuPnUn1rVtGPacrEkD5bRs3DsxRPrxbBg6VXmpjajGaMXzzfUQI8NycAYAgjZIV8DMZHGsdnekKqzRG9dxIzioVOZmslDpasO1IUFDWZoOc2dPLxRbczk08lBM51SI4gHsRpVaGS8LDY3WIxfsrQhE2VZjZBnaYdtqozY4qTckntJoTOLyC2CNm7BijinMhrXkGOPCZObOa/N0kRk+YhiQZiIq9EKV8THtKEjg41wA1E896qC5c0wJd6H7Yim8JPNz1uAsJfigoRT9dBJ415METfyC4R0vIDLvNa7boOX4CpK5DgAkZhz0VUpLgIgeNMnHaac9VY77tnLn78Bq+hwSBPsMQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x75ee4ffd44a0>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"rdDetermineBonds.DetermineConnectivity(mol2)\n",
"rdDetermineBonds.DetermineBondOrders(mol2)\n",
"mol2"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "a2aafb35-70a4-48a1-b0fe-535cc64544f0",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "29f473cf-4bf0-427e-8fe4-0c9adde409a0",
"metadata": {},
"outputs": [],
"source": []
}
],
"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.10"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment