{ "cells": [ { "cell_type": "markdown", "id": "6a0f6d8f-ff7c-479b-893c-64a27bdfc59e", "metadata": {}, "source": [ "## Initialization" ] }, { "cell_type": "code", "execution_count": 1, "id": "40789f7b-9c05-4170-88b2-e3714a317260", "metadata": {}, "outputs": [], "source": [ "from scm.plams import Molecule, Settings, ReorganizationEnergyJob, view, AMSJob" ] }, { "cell_type": "markdown", "id": "28876aeb-a195-451e-be23-e5745174f6f7", "metadata": {}, "source": [ "## Define molecule\n", "\n", "Normally you would read it from an xyz file, but here is for convenience explicit code" ] }, { "cell_type": "code", "execution_count": 2, "id": "8bf306b3-5302-46ce-b391-efccd3e42171", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADICAYAAACtWK6eAAA+iUlEQVR4Ae2dCZxU1Z3v762q3ulmBxuaVQVZYnDBqE8FF9SMRn2JoPjxM24zGmMi8ybjCyYfFZ2ZNyZxXsbERMInijoaXIk6Eo0iIqNPo0YQCZsgCAo2sjTQe1fVfd/fqTpF0dAsWl3dXXXu53PrVt17695zfuf/+y9n9Ty3OQQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOAQcAg4Bh4BDwCHgEHAIOATyBYE77rgj9OSTT4Zfe+21SBAEYfYQeffzJf8ZzKcvLIWhxVK/M/h896hsIpAszDYLUAWdzfR05XdJwbSVfimcpNJp65YufT4nNakKdMqUKTGVzMMPPzx44MCBY3r37l0ViUS2UZjLjz322FVc4mtg8u/7ftClS7H9Eu+DkZfEx589e/bII444YjR778LCwk+59texY8du0Ov5HuY+g3n7Jcc9+SsjYM3+nXfeefpTTz0176233qrdvHlzEIvFKMMgqKuri+7cufMvGzduvNK+jNM5qShs/r7kMYXJrbfeeiXk+Mv8+fOjH3/8cVBfX2/wbGxsrAXPeVu2bDld7wDHNi32l0yD+1smEZg8ebJxBW655ZYfP/roo/EPP/wwoPCCaDQqdrSwR9lTG9eee/7550vRfCrclEBkMk1d9FlA4nsnnHBC6Q033PDcrFmzggULFgR//etfg5qaGmEoLBMahy/gG4c0P1Ze+dmmO9YVscgZxoscWIzY1KlTfzx8+PB/raysjJeUlMR69OgRhMNh5dME6SrDOBvHlr59+140fvz45x988MFifsuXdiShAgMshUVx//79nx80aNBFZWVlLbincY4BuwgQEV7sck1j4Cus/xXy/FhuVi6RJCcIkiRH/MQTTxzbs2fPf6GwohRoqLS0NFxQUJASemlFNj/ExrGAvRl/+mw05W0qWH7nBB7k40tvIocUzYgRI27r1avX2Qh/M3AVCM+ioiKf3+nPFqA6IdyikOdfcGfH8j1+oMCe611myyWBCCDE9IqKigQLkHi2gxWESBLjP9Mo0MHcH7MxzMH+mKPXDTmKi4sHo1imQQpCt1iBDC5Hs7eRbwM09/sE79PBMYBobdzatU7nAkF8abxu3br1pRAvYA+am5vD7B5B5IEKVSWlgg2wOmXsl+jExIkTcwETZeXLbCbvYHkJlqKspaVFWPpNTU0Gy4aGBo9zbT1XliSAHBf88Y9/7Ctlkwsuay4Ig8kDBXIkwWIPFSIBo19bW+tRW+XpeKBNBQ6RJAgTdB8EyecqX5N3yDEBPAJhqZ2aKoOj8Ny1a9d+4bRkoBx6lJeXH6mbUFxdXr4UbOXEJo1FQQYqQO0EjN62bds84hFPfjNacZ98QiiRyNe+e/fuIckbFMDn62byTswxBKXhg6N2o2iEHzGGh/vl4Up53bt3V/tICid9l8UGS1UB50x7SM4QhIKLQBB/69atcplMYRJUeorHRYQ+ffqYcypcbXLBrJWhUL0dO3Y0pUrbfWmSZd2+fbtHoG6UixQN8YXBUzGJLAuWwiNWMWjJDZNSAn//s88+yxm5yoWMGLcAImxCg9VRi1KKdgtU4yJyqDBVeNKE0oIqZGm7pGsl9yEQQT7//PMVKukZM2bIl46aUs+/D5N38BEWp6I0gi+++MJYYRFBeOJKGeykXISnlJDwBGPjkn366ad1NCZuEnRqN+nqEOYCQeQWhBDwT2j7eAe3aiKt5DEshdo9TGESkxg3QS6C1YLyAkQStJ5PoXobNmx4tasXZqbSD27C4jqshF9dXW1IkE4OWQ8pHIunCIKVjrGHIcc7t9122yeqDaQ3Q5d3V3OBIJIL4wxTiD+nkM6EIHtZCWk7uQNUA6cIIvLgQgT8J4zl2URXij8kCzVn/OcvQRhTzY1g/+Goo47ahAIZsGnTJjW0mj5ZclVljaVwRBDrduk9stRSNhDk5/q9fPnyPQGKTnTRLVcIIqEO4/u+SMfEJyjIy9avX9+MpiuAHL5cKBFEBWpdAtXMjBo1yjv//POl+XrSMHbhFVdc8TTEydtOd2okhBwxyHHhGWec0fOYY46J/fd//3fo9ddfN8RQEC7crLsqPGWosSAtxHSFYP4E7UkvJhtuc0LR5ATLk8pJefHpHlGCxnuGQjtPVgKrESVo91XrIgsigkCI8HHHHefR41fBe4Cm9NesWbPr/fffP37atGlr+V+I/3d59+AwlbbJs6poIcf7p556asXRRx8dDB482H/22We93/3ud6rBiqkCRFgq/oAgssARWRaC8z+98cYb38EKN0AyxR5dPv4QfrlEEJsfFUwISzIDQb8ZQe+uC6q9UvwhV4BCDyZNmuTLgiAEHv21VDMZXrp06buPPPLIaWhBab84/82JQlb+D7KZ/le0W4RPO+20N04//fTxRx55ZIyYLmzd0nnz5gUPPfSQkRdhqaBdR/DcifX45SeffDKDd0ip6J6cwS1XXCzKxGwqGBVQHHfrdvpZPYA1uRyinEghDmIPaO+owUU4H8sR51ooGbSH0Y5RhGL8t771rXsgxs0aNcdz8qI2a8KECeroGf3617/+i6997WvjwS2K0ogIGwXkxHTqjBhiTM1L9IDuIYyxuhshx3u4r4/jen0CZrZHdM6Qg3yaXpk65tJmCyismi0y9tPWmSMueeIvf/nLFPpgRXG5ItKGWJwIhIkOGzbsB/fdd98bZ5555pO4CxHchZwmCeSIEGNEhw4dOmXkyJE/oPeuIYcaBFEoCryj4Bh59913n2T4wGWtsdTvtJjDYr+/27rkOWMyu2TKDy3R6uqgXQVnqoOT37tBhvcZWXgkvnacUXEhgnQTj2BtAmpiaqnDP57qylyPR1Jxx0knnfQ+vaG7Kebo16+fOh16WNs4FiP0zjvvrJ0zZ87xKIxaFIax0Elc7fecjddyzcWi3PbaVHDphWdquzi3i5hj6sqVK9/EioTQlgECoYZFn9/xAQMGVNBINgfNeJo6QsIZPIici0d8rIdqqCLUVs3BelTgQsUI0kUa1VQFqgZHUbSsWLFiqjCj6lYNibZ2yh45lbtbl+9M9iWKRgUrxfAu/vXNuA3h1atXx9auXavGQmlNzXwil2M81uUejW1Ptq5/iVd13r8o7pBrRUXFPZBjPFYjinIIK/BWOwcKIka1bXjVqlU3L1my5F25YlIWnTdH7ZOyXLcgbaGmuEJ5n4mWPOmDDz64RtXBWBHkI+JRexNBi0arqqp+cNddd71x++23p+IRNSaqEQzXQ65IMHr06KCTVmtqmh5/P2mN27gDS3kVFRM/UFAOOSKqAlfcgfUwcQf/nU07yEx7f1tg5vL5XI9BDlR2yrssaCHxyDto0rGMLNRIujBC41Hfr6G5AdWXDWjQk2fOnLnsQK6WAnqeFYcs6S7dgd7fXtdCCLRcp7YqF4y7iEIYO27cuHeIw4qooPAZfuyrsoIKjBhtGqryXkb7x0m4mc1YDuUp5wLwQymAfLUgwsYWeIPiERoK36HxCyNSGCAoJh7hdxzBKaNleQ7VvuPHjBkTv+CCC0ahabtBlkHELh8Tt+x87LHHPrK1XWk1OoeCf0bvse+GHBLo0HnnnXc0lqE7DXnDqabdyLGWdK9A4EPU1s1h7H4JikCDzYyrjXsZ4Fp5H330UQMVFYo7GtgVd1is+JpfWz5bEFvSUhLSttdQ9/8gViRK0BqRFcHFUoe8OAIjjfwOFqUnv4+WKwKRzP85J798JcfnEK5Zc+fO/cRcyG6DmcrRCDE9BIbQxnM9abwYUhwDidXGkxq7AaE/ogPiDu4/aciQIXF6HoTUKo6SUIfOKHFYZNmyZdf++c9/np3PrlWyDHOuJd3m63CPliQP4m5cg5BFaWGPoGE9Out5xCgBQuSr96rcEMgR52hmS5G/wmbeRz+lnTSa3f3AAw/8jBNxxStZcLmk/Y3FOOWUU/436dS4/O4ihdIKUVRDFye20oQLIcVY2tTNX5aS/BryYD2ixFURKi1mL1q06FpHDgOTI0gCBoNDKh6BGGMJvuO4IyH1BNagIXx2tSar/1FIR1kQCSAbtwUashundTmCJVHj2p8g1rdfeOGFeq6ntHvyXZk8mGdTqVBKwD2XdJ6ntEGOaDKdPkdfZJHVgyAa0xGXtSCtIXU85LfujxOYh6j2Xvbyyy/nfdyRXkAJ1Zd+Jn+/W008Fi38Cm5WfzSx+mlJyExHR86bMRCyJDongkjAJHDqBk61sWYdjEKSAsalvMrgrYu4r6mdglyVXQhyFEGM53GrzuZdLewRyOyzp9Jq0ytSp6dX1bmKO9RLl8qIaqzHpHXr1qkywmKRv9KQzLmAcFsCAbkpCkiX0d1iPUIlzRuXQMktkQaWJpawqeu8erSqZ6usi44QyYNQcsMKcLmaOXc2DW+/UNsBwXPGcU4+M0aafgEZztY7sRCav8pXepPulSGyiJ0ku0mvTbPywPk4+fFJ63oC82WXXnqpMBAWbgOBfK7F2ksArM9N4HoFQn4yAheDHMhc2BBEQicNnE4Sfdd1WRBdUyySdF8Kcc2iCOb1l1122Rys0SIEWlF9RhraaN0O4wJGCbAn8t7reW+UtBbq/dotqUUSpVHWTiSRVVE+5Aba9KrrP+5WjNq6k6+++uor6LH7e1VZ21q5vUDKwx/OxdpT6MKikPjjL7groyFJHKEKy2JIuKRtZSVkLbRLI0vwJHAihXq94scr/jC7GtuISyKM6X75P//zP8/b85rMfSPA/hNpOhcLEEX45Vql3Cqlz1o5pVvflRcRQwSRW6X0qlqXxlLNCBOi/WM5tVgnvPjii80QLW+rdtNLyFmQBBqmjxH+/BkI/BhOaSyIupyYSQrSAbNa2mpqaWvdp6Osid35rWkJA543kVGL/xOyfMa5EML5ldwXnqkpQONU1Q7U+HvSxuuDVDnaNKcflX6bXh1tGtOOYZ4bx9qMQTmcwT2v2DYV/TeftxSw+QwCeTeWFKGaiHCYya2RY03gbLStZFpWQrtGz+GSpHZ7j72ue3VOG0ce5xfi88/laM591Q89W1ZLMZFmHNFv7XqvTac92vTaNKsbid10Lj2tpE/jY5TIieyv4MJlJsH2hV306AiSKDjrToyWUKcLmxUykUICJtdE7pSE1Aqr7tc11WTpqHv1P/scBcLEA2big68qJ3qnYguerzTLIqWIoXe2Tm+yds1WSZu0KQ1Ko/Ki9Iosyf8q76N1nV4DFhP9zNvNESRR9MbtQYsOSgqLr6PdJUzpgibXRIIq4bJE0T0ijhU6/db/JcB8D+l7prbke03fKb0nnRh6j87pHqVFZFJ6temafsua6R5LduWNaybPPGuQ7sXF+kquoJ6RC5sjSFopIvRNVgNLmNKFTUTQrlhD5NA1CVi68FkSWZLYZxAEG6Jkws3Su5UGvUNCboVdgt86vUqbdt0jolpC67fy2Tq9eh67m2EyTSYcQRJgSMVGEZxVCNKpCJK6ZxgBtEImoZRgaZNwSdh0jyWIzknAdF6Cp13fdY6AOs5v/p4Zt17PURrZQva9eo8IbNNriaF7k1Ys5Rqmn1MaRTbSp5WilLdVyuPEiRMNJvqez5sjSFrpQ4IPrKaV4FiLkU4OXZcg6bq0tr2m8+nCquvaRRRmbgxpbi49Txbgq256Bs8yfcP0TpsWkULpsWnSe3Rv+j1Kgwii89biiFxKp35Dlg++avpy6f+OIInSNP42wvMywmKm0JTApGth3WaFTYIkckjYrDDqmiVJUvA035Y6BTZwnEy3jq08IhNdOMwzaPPow/4UzyxJCriZi9imx6bXpkkkUnptnnTdkj2Z3jBpjBFHvaxrWJA4PZj1Na+3zNj83IBQghfQ6W8BQjaR7zEEMKzGQLVGy52yuyWHhK21wFmBRGiRu2iEhrhX3n777XOt1s4kVKT1ZdI6CQsRJU0Rm1ZVAdu06mjTawktEmlTWrWTVo27D2PlFrKdxSXJhQvSAcFZEEBIbkYoEPjbEThp0UL55lawrTDJXZHGteSwGjv9Pq4bfx6NHNDl5A6epTYV+fSZEjpJuLqt3MFzz9G79U7SaxQe54y1s2RVeq31aJ1e5YfdJ638vfF2nqnl07QUG1/d5izIHhkQFgoQejBGew2C1BsBM7OdSCNbrZyujUUSK3B6jAQzKZTNaOVCup38knHd07gkcmSkH5bek9zMM7Ei95KGm3mvWgE1qN5YPKVXaT1QevUcXC8z9SrdVbZhKY967rnnajhtsdAteb0lbG1eQ5DKvAROrskp9F3qrdF2CJsvK6JdYyfQsmZXe4f9nn7knjh7C1ajkFbuVyHHLeqywWMzZTlSiU0+M0xab8ECvApJ1BlSfp3WLE+lU2ltK73qi8U1n8FhcUZQ9oYsp+gFdNw0WOh7vm/OxdojAaZ6iY6Jl6hTH/2y4kyiFlq8eLHp0IcmjmEVNKFayGpla0H0CDS4GWDFPSGE7nnGV2hMd0vSVTHP3vOqjHzTM+OQME5aL8J6zYEkOsplipOGEGkxFkRxiE0r7pjIqu40msXeo0uJh8WMowS0bPYlXHtRs7VwdBsIOBcrIQbWpShF2FYxSUMVQqIJG0KyHsz6bvo9Scjks0sIJVzaky6WZjlXN40NCOW/MVXnzKR02ee2p7Cl3oG79V3SdCv7YNKotQJN9xalV+lU+kUW/SatHvd7zP/lKa+QPsR8xp8y/n4k07K290jI9sQjo892LlYCToMDAnQqFqQKX1w9W00jHEIW4G6py/uzuCCL5NvT+TCu8epMmaPZ4WOaNgd3bA6u1bFJckhoU4Kb0RLb92HS9uZ9DPOdSfqOZZTgHKqpfdwvrfXhMSG1xzGuSSjYFtEL+Fm1y0B+VSKoFkv5V/f+KnrznqqbFKjrmO+bAyEhARIwjRSchJAoyI1L02pT+wEC5zE46XqOT2ucBf56nEFQOnqMPtTaIxKun1GluxOXRX+U0GbTTbHvK2RUILK/82eKhyC8SOIxc6LHkgZxZor0WN7gafJxvfKlWEXxiawkW1yuJW7WJP3QxHg65vvmCIIEyBvhEIEYFyEgckM0lFBVp5rgwEeIVrM4zDYI9Le6LsGTm4KAKSAOEeyuZeaT5YzE06yLqk3qqE3vVpkuhwBrIYLSFtcgLvJgBB4y/C3Xt3F9NSRSJYTiFeXV5Jm8X6TRlQsXLsx0rVtHYfKV3usIQhUsPrsGNp2A8I+EJFohF5fdTMYQl6aFDI+Bcl+uH09jXMA9cMT06I2rjYF7FnK9GaHqDHgqDc3EGAtlIUQAWQjSafJJ2o+/9dZb+0L8x5IWxNR6yc0SDsIAkpyg78kaOB6Xv1tnKNCORt9oVgT/IiyEJmqIyb2ShZDQoGVFgKe5filulGaCN9eVaM6b69xjWtU6Se2Pce1I/1OyDLJwIoK6ziD0MVzEEDOXXEryn9Y5Xdd9+s4WEwa4Whfph3OzEuZYWOTthiDJldDUPmfLfcJ6mFlB0LCa50oE2MTcvKshxxRdhzxmHAYaWtZDwrcRn/9NAUiVbnu0dxxu2dg0vInQb0wSwFgR8mrW/SDdU3joagi+ifyJIMaKkCWTdzA4m+u+c7McQYxbgTCMRPhPkPuEcIStewVBVK07l+tcKhnPLgKl3C+uBwjhAua/qk+6I9kMzEnWfjelQQ199aR9AQQILAFEeOWNPIz/6U9/WgJR5sq62DgEwhg3jOsnnHvuuSPlZimu2u9b8uRkXmeeMjb5xzrIvVKQnnKv5D5JeBCyh2gj+BvuEUm0ZBtyY/o+KXhXkPtHyUonc0eM24jF+CN5MOmUG0WedF5LHZRQqfA3fH9I8Qn37OVmcV1LIRg3q5PEVYK4Q7a8Jgja0bgjCP6kdPcKV0OdDeVebaeR8AME5urklDnWvVL/pTDXt1FDNF8lR9fwzlTrY9MyHwu3jXSGIYgsiZLqI/yKr67m+weQZjvXRRBz3bphKAtT3TuRbu/6U75u+UwQO7F0FeQ4NUkQ6z6pW4lI8hKCoW7kp7Mb90q1V1gVXVdg+zqDobZ3IvfKyrF1s7aTh9dlCSGAxnqYige5WeTj9NmzZ0fIx0tcFznMdVlOWUhIdCrVvVVMICeC5K2c5G3Gk4Wu4PwC3KdSBELulXGfEBrjlsi9QpBO43qZrrOn3Kuka/IyzzGrTVnJ7ERHuVNqr3kZCyE3SnkybhaZiNG4WbZgwYLTuEdulrq7axfp9T/VdpXidl6gZ0CUvJWTvM04BS/NqHp/0zkRl0LLAcg6GPeJLhg11E4tQkiulXsFkcx1uV8QKMz1elrO5+kZuFed0Q0x+SN98xB6ZL9eIwbNRNVKM/mVVbmW74u4XqPr/DbXIZUqK2RFLtG9+exm5StBpCUlQL0giK2dsrVXxn3C1XiTVae0jIBZUgDrYha4hCCq/pWmfUOd+5K1PJ2VICrfT3Gz3pAVwVLYVnMMo5lz+Dy4ECU/b0IQXY/pKAUgNwwSjT/55JN7Jd0sYZZ3W74SxIx3QPjPxn3qjSCk3CuEyYwBQXAeRRrGYUF6cT1u3SuER1W7Eqb/krR08loeU75w+r8k+OxansG6Weo/1uvv//7vxymvcr+4tpcbRt570/ipNpG8HSOSrwQx7RUQZDLuk9yN1u5TI/2Y5uGnXyX3SgSR+4UgGe2Ke9VC7dWrEpxO6l4padqsZXsVQrdAALlRttXcdE7EjbyK++SGNSavm9osuZJyK7Gck/WgTtJLQEnJ6paPBDFBKCiXIwBnsqt12bgcxB/GfcJKLGJ98HqufcfWbmkoK9c144m07GIGRK2EMNZVy2qhHcbLRBClcSV5WoybpbQbN4qsmM6JBOzfefLJJ+txKRfJDeN6yg3DzRI2Z9ILuJxeAqbHwWG8OyduzUeCmDzjMn0D69AHEqTcJ4TfuE8I07OU7ghckAFJ98pU/9rrHJ/huoJX46p1cklQGjWC8Bm5VxDAWAjyqOpcuVkDXnvttRGQ/Vmupa5DGEN+YcRyD99QHvNxjEjeEYRVbDXDiDrkTUY4Uu4TAiTfPIxrFcNazKWaU+5X6nqaexXQhVzVu3KvjKum7514s2l8WeSQGyWi4HIpycbNYmpUuVFzFaRjRfa5jpKYLMwgSt7JS15lWOMcGE7aguZUzPE3cp84Gutg3SeswxKWQK7GwkwVQXAxzMqwXFfnRLkoqxiWulwCg1DZFmsJW2fdrGu0HKuxSm4WRDBulNwsuY7kRePnq9mXEIv47LY2y7hhwkqYzZo1q0UYdtaMtke68iazqo6lulJroFcyyu6HEKAKjSiihLAcalkOy6IwhHYUQ2iXsczZCAkP2tMICUKjsSGaCOFZCqIZ90rYRdujUNrhmXKzmiHEs1iP6ZBcBFFeQlICZH0Egr8Md3MYCkF5DutIZYTyrmC96sorr7wHcv37H/7wh81JLG0FQDskt/M8Mh8I4st3hhyx8ePHfxfX6XYIUokAmDmvJAgiAoJgJjagaEohzhgNVWVJMnNNxSVhUhCLC/aCfnexWh3rZr0AMaZDEtMRE4voMT5de0h5hhhmUodWmMhSSpH8kP9dcdVVV90FljPVvSbZvd8+W7Dk3KbM5/JmyKEamOOOO+4+JiS4KRlXaHXYsAQCS2KEQiTBjTBYIDiEJPEQVsP46nKvpE3pd/XZ2rVrRzJbex03GsHpIuDZtJaRz1VYyIHMZqL4wxBFmOBumokq0jGBFEZxCBcwiIFJWLEaPQh+/cADD3yf2q/wlClTZElyliS5HoNoCs0Ys3r8CotwE5qxhUJXrZWZc1cxiNwqGsRSi3RqsUusjOaIMqvCSniwLmboLdr3eZFD2rOLCYUEWGmug/zPK0jHEsZpA5E7ZZSEXE1ZDikMXCqz4KfFRZiwa0FTTWbRwnRIN1177bW/ghztssQ16ew0W84SJCnEMabnUWPf99GcIoc2IwgSDJHACoFcKrtLIHRe13WfSIILIuIsouQ6a+fEgwmVrIhq7xZJKchiigxyLbXrN9iYc1ZxCAeLifLP7xD/LeD+Fqzx91k2+iopoCTWB3t/l7yeqzGIL/+YSeD6IAT3UPCaAkebEQSrKa0gSGBEBF2XCyHXSsIi10K/VXslYWGKn+so5ccJ0Lvi0gDGFULQr1N+yWtquQT1u1LeLUGEhe6RYrC44GEZd1P3UeulsTBxLM09U6dOnff73/9+G1hZN65LEqGtROekBaFGxrhAFOY1EKQPZRdXXS5bSluKJCp8kUSCIIGQ1dCu79Z6SMsiOPK9JRDnXHLJJRMIUuNdTGsKjzhTjE4Ak3P0HUx0zigByXY6NsqzsGmNi8UGHFQ9HOd3H3C8hv8H1GyZ5+mZubTlJEFowFPdv1yB69i1FrMZBGQFQcJgXQprTSQQ6buERNd0H8+QAJnJ5CDLlRKATjbE9mAyKe2u7UrlHRz2qaIVNtqVV+2yJunYKC4RPsIlqTR4VEjTJV1H42vBjBkzukKbUAKFw/jMRYIoTwE+8jAKfHiyQU+t5wYWe0zHyAqHhCd9t+d1ZDNLLvN9gmZPXLRoUVdpA1Hao0qz0p7Mv8FD3/e36w92sxhYXJLKQjgJD5/fw5nWdBj3yYrknDzlXIYoWCPNFFwVhVaAAGjWdSMIiifsLp9au9oCaOPYZ9c13Wv/y9Gsn87z+1OLU6Hz9l360ok3gwdVsxWksb/SLcG2OOhosbB4tMbE/rb3WUz4r2K7AnCuUv6ZVdJaqk4Mx+ElLScJovgABTdaBUmhavrQlBDYwlbrufojUXWrcRBmV0OgvuucdgXruk//QRa0jrhYUS6tKZh5T1fAz6RRaUaQy8mHOi4aglhC2KMURWtchInFRdikYyJshQPH0cK8i7mdSvpBt65QwAfNRPIGNQqG8YdVgxWjnn+jCCLhThcACYEK2RKDBkAzObXaBOzOf1NCYQVC/9NzOLZw3KV34rYk/LZDTWHH3GfSCEGU5hbIYTCxuKRjk04OkUI40DFzL1yElyWKnqGdezYKcybJ9tV4yHtyxpLkQkZSreVW/iZNmnQWBXfj6tWrv6OaGGpbfNVK2V3nbA2NjjpvA3JLKksia004BgiQv2NHzY7PP988jE6PO5OC0NlJojIOaA/q3tDQuI5wqidBdkBNnSasMI2C9mjxsVW8+p0MyA20VrlYi8Ix0Hdmvn+GoP5+JqNbYMtAREFhqbass+Njk7zfY2ciiJkxXakEVB0OCKwCwjFjxvhqzdXNbKXf+ta3vk09/zU06J0lwX/66adNbKHGLhW2FQQdrTDou2pnVGujGisRRFpW2tRqVJGFXWPVw9u2Vc+fP3/hJN4fobq3iwTq6oH7erSysv8roVABi356MUhgehPY/AuD1vjonBQH1seUiayFSCLXVLvIIayYhdH8xrIswOrMZp3DuclFeDwRhY6fgarGE8XU5udhlX+bT8nwhQ5vKJRbxO5J0JPEMFlEUNUZzpxPz7PuB3RV3RrATzrppCpqUS6nwG+k5mq46uqltWjYi+MChd98801jLSTkKmjVxtj3WHdD10QQe82elzBoF1G4h27i9QyoGlg2evTkIyjwzydMuCPy+ut3iqAHJHN6+rP8HU2jNqHXo6WlfY8Ih3uVBUENMYivLu2pKl3l2+Zd6ZOSkOulfNtqbmGW7o4JF7lgWiMFpSTloR4KZ2Fpz8Ki3AmB7mfSi8cp10+Tz1S3n3SFZqA43PI3f8riR0daEB8t7KdplvD06dMrlPe7775b/rKxDFYD8Tt01113RVV42s4///zjWELsBqzDZfSs7SGLQGHGRAIE3DRarV+/3nv44YeNFpTbYDWkrcvXURrSWg8rJFZA9ghEjJilBmINpt7/JNysTZt379763cWLf/28SYynvllmSGriZ6f4VJqeBsPA6917JFOrVs7s2XNoZW3tCgZ8rWGGSHWjKTKWVDjYNg7h0RYuUhxWeYg8ikXGjRvnDRs2zDawxoQlZWA6NbKQaQ33PMHxtyy5sFiwiGjMx6X1R7S+ookXk3CFb7zxxgrWhfS4N1X+lEWI61KIHaKEOoogyrSxAD/60Y/OYgzGFLT/GWj//ioo8K1GOBexpNjjaJiFSQB18BmXcA5LoE3DQpxHVxLNIWuIgXCrZsaM7ZDpV8BN4cQZThow3X9IlgWCGL/bCkA6OaQp0wmyhyTqelKPYJTga1+MFi1iJdvm0K5dO9Gg1TM3bpz3T9XVS+uS1qSTuFwJlwq8yvr3P+me8vIBdPM/AjKU0xugKLR27R9wj7agHBSDJBpIhYnFRUIubKz1SMclSRDFY7Ioap33tZYju6edshFZ1CFUjYhm+iBqt6KUxZ8ok3uZRWU+6UoJ+8UXXzyRzqSX08H4DBRdf/0f97iastIY+SeZdmiBCl4udZoy1amsbB1BEL0zoDtID5YF+83YsWOnjhgxwoMgRgtJy8sKSNNAEA8r8Mijjz56Fyb8bArjBgrheN2bFG5ZFGkrM1WPalhEDGYciVHv7zGUNLxhwwZPu55JwcexNFrjw/zfCkFrQRDyEoSEMMTiTU0N/H8sYyZGI1DdcUGKGdcdMHy1JYQ1+WD37k03rlx571tSdJ43g/2g/nZ7FS6KRybWD8rLjzylvHzY/d27D/56aWkP8l2GcijS+idgtNZbvXouCsr0KAgJ89YEESbaRY5WBAGWeEjulXDt3bu3h3DHcKtEkDBCrrEypjxRYpqYT1Y9ItJRLuoq/z5l9Ntnn332VX7fzjCEv0UG9IwUwRQzKk0iIWNV5jCBxvfOPPPMGsoasciuJck2QdSCGwBKnzPOOGMRBDkGkxrDVVLLtwTXpkcaRhYmDEj+I488ohWTClUIFJYGOplCImg002WqsAA7APw4pAhh0n0mdfNYUHMFv39NIa3jf08DbgmFpcngjOVpTRCRUrs2pTMaVRtKc3j79hgEGcJ0nQPCFRVV9NXqj/BUUIAhNGMsUlPzRVNd3Rczli277af6a8dYk5TVQKN/40clJQNnlJcPLiop6Ub8UcbM9Yqx4sRUtaS7PrZlyzvhXbtWkdwwAhwKQRAzqXVS8RhypBNEmLBpoEyEOKMBzC9dsWLFMO65CSKMwpp7LBIaUJZxCBKirPwkaRSjqNZM685rCW0D8Isvvtj83nvvFSILAf+TDJj/yAohB7b8dW8YQq1k7cUzuHfrjBkz0t1yk6j2/MhqkI6ZNAJ42WWXPTRq1KhjAK6ZvVAWQZotbRMwpjsEBRaDSIUbN240c+NKe0EaM7+TrEWSFHFG/0UYqxFmOh4t2fwK+yzK8Tme06K6fDTUhcj8o7gFlXLBeEaUZ2taGy0aY3u2qmDMLkGAIEwxumvpxo1belVU7K7CkjDBcxMzoTf55eUDSXNvJrYupaW+qigS6ftvxx5779l1detvJHBf43lqD5gMydtd44GVGiyfYlmDqqNKSkbfX1Y29JzS0j5o4XIsRzHaWxNuq+1nN5ZvB8Ntv6A3bsOnKOTtyOyx5FXtQuibqI+ltpjI6mp5a+v/S1YiBPebsdRXIrALpEz4zyy0/MWUw/WU0SSsfJhez+r5HOV7KLmLLGFZG56n5bWl5ApxpWIorzB7RLWOshxJOUiVP+9ohoDHUE39EN9Vhh6uFknJzqaEZGXDpaLG5/UoWuByvs+hQa+FTBeoFkSaZn+bwFAhUCABVYW+TK6EHWuRIgakiDDST65YPaSYy2I2v+FZuDuIJv/lGQrYlc8oRNRQ27v5PpWCUhcJ4z601pQSGAq+BmH5Dc9WafTj/pllZf0v6N79aCzIUHqyDgpVVAzCr5dVK6d2y2dGkFhk587NBKab/2nNmukP8D82ESVVFZ04lbHPPZUD3bufdF1p6eB7ysqG9ygp6YGVLKUat4i+UsqLZjHZwSpT2xg2vInYaf28LVs+/O7o0d6WmpoBd5Cc75G/HhZvvu8Tf4CHGkjnYDmmo4w2J8tTysRUpiSzdIqehbB/mzItVcDNyEWqlytDyd2TpZFCvP/++01bjFbgFaF0r5appnySj9rr0MKvAryCqdzzOHFlBJcrK/Fe1ghigyyWIX6VzJ1JG0YMckQI0M1ovr3gaPVDBYc59qg2FDliFJAPEUKyFvjBn0CSRyHHg/ztY/2V+5F931YEqBC12d8eWm0U1y/m3LnsQ9l7sOtPWzgugyDzIeM83DNTRck5syE30woKyu/Gvy8uLx8R7dZtUAQ3BsIMpGB7474UMBuIrI5a4jfPqal5+39t2fJA9eTJqrI21dI2LfaROlIGsgCjOU5MO/8F97ZpgZL/Uc1ZGYHt//hFaemwqSUlA7AaFTGRo7BQCr+ZNO0SOaKNjVsi1F4xe+Kn02tqPr437UWyrlX8vgDczuE4lmM/jsKwhuN6fr/M/pxcVn5rS2HJd9KS+m3zN5xz1+K2XYm1YJn5Iarp0jrtAWUe5jne4sWLTQ2YrokcOoooUlb72UQGXN3tr2GFziYtZk6v/dyX8VPKXDY2vUcB2xFMnLAK7VOB5gggh68AHRN60DSsXLnSe+yxx9TN3MO8ixiLAfq3WJTH+fNOCtO79NJL1XaiQlL8sr9N6dCeuk4bShFkKNXNPHNn+jW+y/pIQ0ogtOl/44g5Z5eWDhpXXn4M47qHQZAh1JINRjhlCctwYzxquuL4+BvWNzV99v0NG/5Rs8CzpVsT9XwdQ1oOZl1kAWeQJxv476m+LS8/4YKioiPvKy09aijuHsQoI44rwWVU/ygtdSCrsRX3aRNB9eol9fVrrqmv37mEOiHefacSpPzYPOq3thACawoEbV6PUmpKnE5c41P4at/fJpyEr2nTQpD1nMsR+huonTouSRQPV8wE5bjZhhiWIFibtiyI3ufjUu9iOe6RF1544ec8WwRuKx3cnpltv3TNzKP3eoqAU23GEMx0BQJpum3YzoIHI4gatag6BJMgSnD3Cj7vr6k21ORtURGDgF/uW1z9gfZ6674/BKh2pUd7PCkA6UIggVEhS3js8yyhhNcSJiA9pbZ23R1o5+nNzTupBm5EQzdEunVroopyID58L55RgNYbMbS2ttcLgwc/fO+GDbdPhwyNnvcaz5jBYxKt8JWVJ5SGw+efGAS9jwlwOlgOEdVdVEs76EqWGHxv82a/npvZ7P+ekjYtrqg49+7i4uHTiouH4rfLpSqj4qGAS1EsRsJqNDVVRyAFBPnk7h07/ipG8H4F83v1ALAKQHgImxiKYgdHu6Ww4oTFwV5rfbTXFT/qfzspn9/inj2AYjsXK38TVe6TiDsjuGBmomxiH8U/qpI3rfH7c7EsGXD96CRRMYTnfk5Z6/m2fFqnI2O/s0WQVIIBI1B1rHbVPhHcpQdnqfvsF8AxMQiWwmdOq1pcs2+/9NJLTfKTcYUiXI8ptrH3H+JRBWkLU2RI3w4Eut6jgmkkWbfi078Rjb7/H83NNUc1N++OEcSHIIqPRaHatD+BZylt1gOJUfpNGzToVxOamlZ+b8uWM018VFl5zuBQaOzNnjfg0iCoHOL7lRCjO48uYddE2c3876hPBg68+ul4/KNfbt585gYuQIaxpxQXj/lNcfHIcUVF/ZiJpALLVYKS1rruxmoQD23DcnwWqatbtaahYe0/1Nd/hgVTNiW0+8XqQHikX1MSDmWz/1H3ESkcKbI/nnfeea+yYvBmFGRPyl5y4Fs5UGyJC2UaLlu7WVKC2vQf3ZfNLVsEEWAUYGwDGawHiDKIEXBEmBLjwYkLVNB75Z37TVDOfXHuDxN7bJIlSQaIEuTDJcZez0/+kNY8nE15UYmp4OdFo41v1tWt+E1Ly46pLS01uDUNMUjCJHRN5GcgLk9PbitGQMaM87yKhf36zb01HJ6NZj36575/TE+5+0HQk5qwHpCpG/fK2ytiRwN4EYjj/dD3+15bWfn/btmx45+6FxeP+7fi4qMLi4p6Jq1GEcRSt5BaLNk2SLol3Ni4llq+1XN27lz8PR5Uk7AaZpSlKQfOHWg7XDwO9iyVkSEKlkm1Ypuo/OhJeWpWlbAqXNSIq1os1WCpvUTtKDraDUWoKVMVd9YRqBtFQQPyoeTFPuJLH7NFEIGuOZiqqWJdQUZPoCZD8zKZGfyUerlRmF5DElkHQDSd4aQxRCYFdtQovYVpbSbgl0uVyYI8XAD1bhW8SFITiwVXNDRs/q9odPdv8Pt7IKxRWt+ZJqeBXsRD8KtlTQohwNjC3bsX/nssVgUpjsIC9uQZ3Vl6uSdZ7sGufmRl7Oo0qRo4L6CxOhYKDYdlfX7HYldUvTYjSLIcJab6Nh5XR8qdshqx5ubPaKNYVtPUtO57dXUfzUnw+FLSaNyyw81jJu83ZUXXkmbaO97CIoxRmaLwDDGkGOVa4YKT58C4WnK71Xip38gNMV1dCHKsuOGGG6pV4cMzcoogAhuNqNrT6C/I6GPUeWu2PlMIlgzEFimCyHqINDKraBkfcAPihV+aP+xxj5I/O+wgKyZrQt78OTTCvRuLLfmPaHTbBVgUtHodrk4D03s2kq9hdObbgcC/hkIcAAm6UcgVEXZsRQW1N90REM02opbrRH6obWaQVoj2GI9Gyx6sUXIxgvQOVqmIKUE126Ma/bYxZ1c1bUMfRxobV87btevdf0DdrKEGTMRAiA4alyVe1v6fRqDRAr9E6K8j/vTpKaFpTdXuQp6Ny0x+mo37Lc/CxiPIgbrV+1igX5BMPUeKPecIIo0boubiSboO3MwE0N9ACzRDmEIFaWr0k6mVmZUmQYoMQdAeLViPwqVLl/6SbtQfqvcnDUUHihN4TVY3aUfSE6jQ1sTj0QvR4NNpRiGI317c0lKHNanFgmjxndfoLxZgMXpxaykkKSWv6mSpSbIlECFcC8yS7BIbOgKBMTu9b9WroA+4jEarfoSWlUu1leXTPoUYyxpbWlbfWV+/7O7EPxWId7jVSCRlz6cEOowX8CFW5FeU8c24zs0oSY2VN+UtRWllgRpPQxDJCMqykLE9f/7JT37yZLK5QLKUlS1bLlYqM7hGMar7LsdavEpXheEEaTG5UbhXikdMXyGCNA0LjUMOxl9sK+S+V+kBeovIgYuVFc2RSvChfzEKQLfD7bvj8R0vxeNvzY5Gt4yTpo9Gl0KEHZDiWO7QnFtFaE11CFS/MPWs1ZqARB+EHyKJ4lIsByTwEJpEIvA2IFkVruc6jh+zb440Ny9Z0tS04pqWlk1L9lTf7jcQTzykYz8NSfAebkEhjoEkZ1O7peW2idsa5YL7kgXFmRBEnR3lYhVyz8cQ63K+x2bQ1SSbWcjqy5IZkwOhSd0Gk+FfY0YvVA9OdUOgCs9YEJlbuVeQI4aW+fkHH3wwg3ubpWnYjD+bfFZnPaB4fMQ7KCbJdxQXHzGtZ88hiP7XyPsQ0qyarSMgR3/I0R/3qx+1NwXs1GFRT6HGZLlZIoiW8aDCz+z0kDGEiUarUR4/A6JF9zY0vHcnD0xW33ZaYqSXkylE2p8KIcQMLtxCuZp5khWDShYkB4o/II5crhfwOG7Cxd6QtB5ZVZAdQRCBZUiiLxDlXIT/WkhxLOa2v9wrAFsHOPOpN38cEJeIGNyjtHYFcihb2sijAsmArjR9Xy8oGH1GPD4Kp2kQLtYAQxCyiyBUQoy+NDaG2WlmLEtYEUsQWQ8RBA+UeMwjYPVYBi4ebml5dVF19bkTEvUEsRSeiVd3+s9UWeJujaPsL8eNOocyHqa4lL0aOViKy/UgXYheVm46ghx6b0cRRO9WoUrgrdCHqbmoUCzCtkMf2tLcKntf4kLn/zTYDh/es6KhoWwF1bqVNBySh0HwvZI44wgE4QgsZiWk6GcIQghmSCI3S3GIYhARRFX/gkU7VoRGVjp+xVZuLil5adTHH0/nqqBp/1blDEMufCQD5DKx0aLeU+XPTp4S56UYtfE7q5YjmSRTG2C/Z/toM5wMSb0YwKSIkWzrOJTW8Wyn+1DfZwq/qSkylOLtT3kjxTEKGr9pr1398GJYTVnOxKPtUb/0vdXOmRakJujf1NR3KLd8QEAuDFOCxveusInVSrNw0h5Nb8GXYuQceaeRZ48S1amsbvjKHb7Zgk2Kh0lP8CVaxzs8I/tPQCFxSEDXi2YEgiqptD0ImnAdGwnqm7AWhSbmoCLHWA+5WLIg+q0jgT/3apdcNSE4ce5gOGDX36QorbJMycAhdBvKSs47A0FsRlXyObopa8TRHv6SOTYi6A3sBBSxOoixm2C0hABcI/gSZLAulmqxtFuiJEil7lnat7Pn1NbpZKAzESSnSprMGK1IsLkhGqUfiFdPCF6HANTR7YJAgj0er4UgpQh/MSRQY1lPzoVN24clitpBFIeIJKrVicd38Yx6nrGjNhLZtCEBWna6XeRaAR1KfhxBDgWlL3ePtKFaf6upqVsbCjV+PQh2Qpoe+NbdIEgpwl4EQYogiJZ/C3FOM0GqVb1QvjfXRQqRI2BPjO1gUkfO1oaDYPvaDRv+TzVOF25JlwvQvxyiHfAvR5D2BV2BJiIeeoVYE4JshzTdOVUKGYrZCyCIikDV2JokogVCNGBJ1MpexG/10NU4cpGjge+1nNtFh78a/lP9Ch+QZYYeoMjfbe2AgCNIO4Ca9kgEWEG3/zAdA/6Rlgy0/VZIUsJRXdrVpUZxqWafl/skItRDkBLO2euafV3zCyuQl4u2m+4q65hVfd3DejabeUfiq/vMNAIqHbe1LwKmCraycsBMiHJDPF5AvdQgOpNU8tY+7L2xFj0gRQXHbuzqzavKKW6h9pO+W5BHSzEosGeeodDuglhs2W83b/6/3012SLS1gNzvtkwj4AiSaUT3fZ4w9ula0a2urv49gu+jg6AEl2gg1rsfl3qz94AUFRxFDvqbmPEgCYIoLpF1IRJhTEl9hMFTH5WV/fnENWu+WcvIQCxKqqGVr27LNALOxco0ovs+T0Ls01V/FwOBzmdNn/kE7MPi8Y1YEqJzrxkCqcqXjlYmeJd7VcSuoqHO18OEEJmEQg0F8fin62Kx9eevWfPOLs/7pi7q2W5rRwScBWlHcFs92rhajMXWZGtzqaQapxgkCCpwkfpwrSfWo5y/7CEIgTuTH0Tpu6Vq4S1LotGN32aa03XOtWqFbDv+dARpR3D382hjEujuzziRltsgyfVU5/ZhfgO69xfRtlECEVS7ZQw7A8o002TjVqzLrIKCpf/MLEe0NpoZSVxgvh9w2+OUI0h7oHrgZxqS6JbKyhF9YrGa7zLU9p+Li3tBiFKm7NHcxBWMhejtb9/+9m0tLYtnbt7sbU0+MvXfA7/CXc0UAo4gmULy8J4j3InCfdrJg1G9ex+//Oijrw4GDTrd79dvjC4xXvszf9u2HaPnz//aitGj7yhcvvxOYhYXcxwezF/9bmkkt2UfAQXX6oIoovTQpNKRiBaraaYLfJM5hsM1dDlZ2IPrvmah5OgCckDI9uZqsbKN+F7vCxmiqN8V82czeMpnfIg4o+9hja50xNgLr+z/cATJPub7vFHjPTQOXUNtNVhKmx2bnvjlPjsKAUeQjkI+7b0JC7KHICKMI0gaQB341RGkA8G3r25tQRxBLDIdf3QE6cAyEBE0QDDdgmjqH52Xu7X/lQA6MMF5+GpHkE5Q6NaCyK3S7ixIJyiUZBIcQTpBWaRbEBFEv10M0gkKhiQ4gnR8OZhJbeRO2VosjUcXQVg+DWfLbR2JgGso7ED0n3jiCXVgjIdCPuRQ24fW/0jsie8FLItg1tfowFTm96sdQTqo/Jn3y58yxSy/xhgPjSaMEnsEBOZaQFOE0XQ/O5n93Y+NHt3XWZIOKicHfPaB95lG02eG+vg3v/lNVpod/aO+ffsP7tdvqF9VdSR9sSohRwGGo5alEmo3bNlS/dNzzjlnJie0JoZa3l2Xk+yXmXtjthCwMwaeeuqp9/3whz8MHnrooYCZ64Ply5exSNBnrIOxk9nOGyFHs2Y91zJlwVtvvXWf0ufcrWyV0p73uCB9Dxbt/i05z3CM4bf/m+WQb2LZ42bmI2Yp6fJQ9+49WQKir6xHejq0Clf0xBNPvGnhwoUbsCA/E0nkdqXf5L63HwIuBmk/bFs/OaTpNFlRqYop/u9gZGGMNTAKWBNFS7eahUxbkUP/V/kUMAJRa6rc8cwzz1SJHJrpvPXD3e/2QcAB3T647vNUgnKDNYvHaM3wUogSsAaGr3UwtEafVtZqY1OcGLCGeCmkukH3TJw40ZVbG2Bl+rQDOtOItvE8JuPWMFlGDEbOghQBU5L6WihIu9bE0PEAmyEJRDpLz4AgbsjtAcDK5KUDlkomX5Tnz5KAx4k3unMcySRyLM4ZDWmhUu38Phg8zA4fMGY9GDlt2rTuuFlqH3E1kAdDLQPXHUEyAOKhPkLBOGSIaHk57azLZ3YtXKmJqQ+0afVX/hPBRXNldiCgMnzNgZ1hQA/0OKxFI/sulrY2y1tThWuWPNaqSjrX1gapAl3XykusMd+IBWnrVnc+wwg4gmQY0DYep8a9MMJdz/FtBF1rv8fZvR07dnhbt271WOp6vySR+6V7Wcw04N63Z82aVc85VfW6BsM2wM7kadcOkkk0D/wsq/YfxxpM3rRpk6nateuBa/FSxSNa7ZeqXxO0y+2S5RCJIJfPgpaP6xUzsrwU8oGzldtXHUGyV74KMkJU8b6ANXgbgpxM+0cLNVgFqsESORSLiAwiDbVdJi7hXAuuWMHSpUvffuyxx15IrvZ64IAle3nK+Tc5gmS5iFnvu5kW9L9D6N9lvt4SXKUWyBEhYPdlLexa8VgUKqqCKAQp4L4G/vd3+q9IluUk5/XrrNnPaxCynHkJeLyqqmoCxPgd7tRRAwcO9NgD3KuYCFJYWAg/wr5quqqrq9esXbv27+bPn/96R60VnmV8OtXrHEE6pjgMSbQuOCT4Ce7UpVQBD8H9MvGH4hEC8U+wKk/TwPivJHGHI0fHFJQjSMfgrrdqsJTpdMgahqUQ5WvEImNFDia2Xkac8iHXVevl2U6O+u42h0A+ISAFJaLss6mtI9k13imxfdDJ3gkHfvawPtCbVA5yu2x5qI1D/U9cWwcguM0h4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah4BBwCDgEHAIOAYeAQ8Ah0CkQ+P8CWgWTJE8sswAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# molecule = Molecule(\"pyrrole.xyz\")\n", "\n", "\n", "def get_molecule(input_string):\n", " job = AMSJob.from_input(input_string)\n", " return job.molecule[\"\"]\n", "\n", "\n", "molecule = get_molecule(\n", " \"\"\"\n", "System\n", " Atoms\n", " C -1.12843000 0.00000000 -0.35463200\n", " C -0.71293000 0.00000000 0.96463800\n", " C 0.71293000 0.00000000 0.96463800\n", " C 1.12843000 0.00000000 -0.35463200\n", " N 0.00000000 0.00000000 -1.14563200\n", " H 0.00000000 0.00000000 -2.15713200\n", " H -2.12074000 0.00000000 -0.79100200\n", " H -1.36515000 0.00000000 1.83237800\n", " H 1.36515000 0.00000000 1.83237800\n", " H 2.12074000 0.00000000 -0.79100200\n", " End\n", "End\n", "\"\"\"\n", ")\n", "\n", "molecule.properties.name = \"pyrrole\" # normally the name of the xyz file\n", "\n", "view(molecule, guess_bonds=True, width=200, height=200, direction=\"along_y\")" ] }, { "cell_type": "markdown", "id": "4993e4c1-5448-4b11-96b1-fcc9ae1b3f5d", "metadata": {}, "source": [ "## Setup and run job" ] }, { "cell_type": "code", "execution_count": 3, "id": "cee28e69-7b62-48be-a617-1422c2b68a67", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[23.03|17:38:48] JOB pyrrole STARTED\n", "[23.03|17:38:48] JOB pyrrole RUNNING\n", "[23.03|17:38:48] JOB pyrrole/go_A STARTED\n", "[23.03|17:38:48] JOB pyrrole/go_A RUNNING\n", "[23.03|17:38:53] JOB pyrrole/go_A FINISHED\n", "[23.03|17:38:53] JOB pyrrole/go_A SUCCESSFUL\n", "[23.03|17:38:53] JOB pyrrole/go_B STARTED\n", "[23.03|17:38:53] JOB pyrrole/go_B RUNNING\n", "[23.03|17:38:58] JOB pyrrole/go_B FINISHED\n", "[23.03|17:38:58] JOB pyrrole/go_B SUCCESSFUL\n", "[23.03|17:38:58] JOB pyrrole/sp_A_geo_B STARTED\n", "[23.03|17:38:58] JOB pyrrole/sp_A_geo_B RUNNING\n", "[23.03|17:39:01] JOB pyrrole/sp_A_geo_B FINISHED\n", "[23.03|17:39:01] JOB pyrrole/sp_A_geo_B SUCCESSFUL\n", "[23.03|17:39:01] JOB pyrrole/sp_B_geo_A STARTED\n", "[23.03|17:39:01] JOB pyrrole/sp_B_geo_A RUNNING\n", "[23.03|17:39:04] JOB pyrrole/sp_B_geo_A FINISHED\n", "[23.03|17:39:04] JOB pyrrole/sp_B_geo_A SUCCESSFUL\n", "[23.03|17:39:04] JOB pyrrole FINISHED\n", "[23.03|17:39:04] JOB pyrrole SUCCESSFUL\n" ] } ], "source": [ "# Generic settings of the calculation\n", "# (for quantitatively better results, use better settings)\n", "common_settings = Settings()\n", "common_settings.input.adf.Basis.Type = \"DZ\"\n", "\n", "# Specific settings for the neutral calculation.\n", "# Nothing special needs to be done for the neutral calculation,\n", "# so we just use an empty settings.\n", "neutral_settings = Settings()\n", "\n", "# Specific settings for the anion calculation:\n", "anion_settings = Settings()\n", "anion_settings.input.ams.System.Charge = -1\n", "anion_settings.input.adf.Unrestricted = \"Yes\"\n", "anion_settings.input.adf.SpinPolarization = 1\n", "\n", "# Create and run the ReorganizationEnergyJob:\n", "job = ReorganizationEnergyJob(\n", " molecule, common_settings, neutral_settings, anion_settings, name=molecule.properties.name\n", ")\n", "job.run();" ] }, { "cell_type": "markdown", "id": "ae2e3890-cadf-4bf9-a735-1c8226a9c366", "metadata": {}, "source": [ "## Fetch and print the results:" ] }, { "cell_type": "code", "execution_count": 4, "id": "43b1bc35-ba6f-4c69-af03-b51508104386", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "== Results ==\n", "\n", "Molecule: pyrrole\n", "State A: neutral\n", "State B: anion\n", "\n", "Reorganization energy: 0.473683 [eV]\n", "\n", "| State | Optim Geo | Energy [eV]\n", "| A | A | -63.801633\n", "| A | B | -63.487503\n", "| B | A | -61.702138\n", "| B | B | -61.861691\n", "\n" ] } ], "source": [ "energy_unit = \"eV\"\n", "energies = job.results.get_all_energies(energy_unit)\n", "reorganization_energy = job.results.reorganization_energy(energy_unit)\n", "\n", "print(\"\")\n", "print(\"== Results ==\")\n", "print(\"\")\n", "print(f\"Molecule: {molecule.properties.name}\")\n", "print(\"State A: neutral\")\n", "print(\"State B: anion\")\n", "print(\"\")\n", "print(f\"Reorganization energy: {reorganization_energy:.6f} [{energy_unit}]\")\n", "print(\"\")\n", "print(f\"| State | Optim Geo | Energy [{energy_unit}]\")\n", "print(f\"| A | A | {energies['state A geo A']:.6f}\")\n", "print(f\"| A | B | {energies['state A geo B']:.6f}\")\n", "print(f\"| B | A | {energies['state B geo A']:.6f}\")\n", "print(f\"| B | B | {energies['state B geo B']:.6f}\")\n", "print(\"\")" ] }, { "cell_type": "code", "execution_count": null, "id": "ebd5620e-77cd-4cec-b27e-138ab94aa69c", "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.8.16" } }, "nbformat": 4, "nbformat_minor": 5 }