{ "cells": [ { "cell_type": "markdown", "id": "b440535c", "metadata": {}, "source": [ "## Initial Imports \n", "Import PLAMS components and set up to run jobs in tandem." ] }, { "cell_type": "code", "execution_count": 1, "id": "3dfcc7e8", "metadata": {}, "outputs": [], "source": [ "from scm.plams import Settings, AMSJob, Molecule, Atom, FCFJob, config, JobRunner, init, view\n", "from scm.plams.recipes.fcf_dos import FCFDOS\n", "from scm.base import ChemicalSystem" ] }, { "cell_type": "code", "execution_count": 2, "id": "eea44527", "metadata": {}, "outputs": [], "source": [ "config.default_jobrunner.parallel = JobRunner(parallel=True, maxjobs=2)" ] }, { "cell_type": "markdown", "id": "cf343e36", "metadata": {}, "source": [ "## Setup Molecules\n", "Create the NO2 molecules using pre-optimized geometries (usually the geometry optimization step would come first)." ] }, { "cell_type": "code", "execution_count": 3, "id": "d2c09231-a063-4de5-b391-8a383dda24d9", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[30.03|16:39:56] Starting Xvfb...\n", "[30.03|16:39:57] Xvfb started\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADICAYAAACtWK6eAABGHUlEQVR4Ae19CZgdR3Xuqe6+y9zZRxpJM6PFFjb4ScY2CJywShgwBExYgmzi4JAYkNkch/eR5IEDsoAATkiw0XMA8YCPxxJbSsILmCXGi8TiOBAZL1hekIW1zWgkS7PftZf3/9W37vSMRkKWZq5MbtV851Z1d3V1zenz91mqulrEJssBywHLAcsBywHLAcsBywHLAcsBywHLAcsBywHLAcsBywHLAcsBy4HfYg6o3+K+N2zXIxF1HWglqBs0EyMOiURrRUIcRHWbTpYD3smeaM+rLwcIii0iDq8KoQ+QnZDgoxLPIREsIXKbngIHLECeArNOR1UDjCooCAxZv2xZdmWlsqxZqeVZkcVOFLlAwHPSInuySh0KgyCPao9UyuXH1JEjoyhrYKAtl+dX22LRpt/AAQuQ38Cg03hYrYdAQ5h99CEAKDoWVyqvykXRm1orlQubRJZ04XirUpIBARySQs4bGiEvRZH42eyB/b2994ZR9KN8EGxRBw/u4v9D0CFDLatRyI/jJTLKpqcZB+A7uDCntLa4dtGiZfNd9+rOKLqs23EWt6GvGQh/ijkEHEAJ0wQIiPsIEEcpBY3iADQEARARyWgY5oMo+t5wFH1x+cDAbdx9F6q/LAYgN22agQMWIDMw5XTugtbwNkBoL1+6tPPsILh2XhS9o9dx2rMQ8qYoCnLoXA4AaALBnFJZx9EahFrEAwEYQEX1tkY4KYpCCUMoDWCHm2EoI1F065FK5UPPOHjwQRzgKRHOYB2bpnHAAmQaQ07jJk0qBXCE6/r6Xrsoim5Y4jhnNUGgmwGYFmgV+BwKJACHNAEYJIBEMsg9kNYXBAfJJIACACFFEgQEiwJQnGIQFKFVNiwcGPgUqxIoOEv7KuZUm4OVlglPCw4omFUOzaor+/o+3idy7QJ0qzmK/NbYz1AAiTQDBNAekkNeA4jriktwIJeZQEJwECQBLDZDvh+grB32kSD49n/6/uWvGhycsCA5WhYsQI7mSb331MDx1p6em5a57nvawjBox8ML/oZDJzwJDg0SgIEAySH3PHgdBIehJEj4nxgNMgkO6CP4/ZUKvHgAJYq8Md+/vz8I3njO4OCvN6OlS6v+T70Z8XS8ngXI6b0rajUEchtE9k0AxxkAR3sYVgAOjwCB9pAWAqSqNQiOZgCBwCBlCQ4DEFMmQAzxf0sChMCIwSFSLhMk3PZhdnljQfDrx5V6yXP27dsfoQUIhjW3wD4GPWw6TRww4Hh5b+/G+Y7znhTAgchTCqRvjAtgGKIDTj+D5pQHcKQMIJinEL9KbicBwv+NZhY1iAEHj6O9KnkAS6XVdc88w/dvvnvx4otl374yQMIayBo7WYCcpvsPn4OhXH/VokXv6nKc96UBDog5RX0SHCgjZAuF4NRIA4TbAIkGBc8wZEDCYwYE/P+SADH7Kf5M1DBRlCJIOjzvxWdUKjfhyJXYSx8FqGrsZAFyeu4/HfJw+eLF58Lf+CwiVQFuBMVbSyWe77VwrQYIQZIECgCgCAIDEqNFjCYxxwwIDEC4n/vM/hgcxgxLweSqLPK8P93V2/sD1d+/mSBB7YYGCe+FTXXmALQHH99RSxBsbFPKaA3FRzYPJEl3DQLNsQ1N0ACKWoBkgECAzAQWo1nSMNpmIgMoc77rEhARggI33LdwITKtX9idhk0WIHW+9QCHHiXvW7ToTXC617gY/KOoTwGFecqbnH00ZQOO6XkSLEbgmRsyYDG52Z/MXUzpUipY4Lo9rY5zDfpER72hZaSh//k6Y0NfjqaVrFiRxsj3xzBVRPsbvAnTH9P0jmsEcIQgzrHivhpYpoMmuW0Ac6x8KjBijRSDzsXcrRADkB+6p69vMfoWrG9gkFiAUODql+hmRAuGhl6bcd0V1B4ABi0rnQwg+Ng25WSu9wMEzKeARJ+d2Ge2mRvQmHy65pm+DTsOQAwXOE5zaxS9jddfYwFCTtpUBw5o2QYi1nmYJ4WnE+VPJwMEAw7mAYRaE8vTCE95nH0c0q3ih3VmSgYwJidQqmVgxMH1IgD4irtWr/bWNLCjbjXITMIzN/vI67C3t3cJtMZqFUUKoutSfAmGJBkwEASYgTsFHD62STMChNGq4xHOqwGGZZMIDCYDFgCE0+UxpeWZ7s6dq3A04gh7XKmxfi1A6ne/Na8h/K9AsAiz1LXcq5mAYQCSzPlSiAFHBcJL4szcGiDMVBKTm1Fzs21ycw4BQmIyebwV/8JZxyi+yobhy7mjG/BJHm6UMm1im+rIAUjZC3g5iqcBRxIIpqwBgXrMCQaCwwCjAiEvgzDyLp4RfONLGG1AoTdmEy9ogGHqMzf7pgGE/SJ0dK7Ui3n6mniTxYZKFiD1u92UfQem1fl8TQPCpwwYmCcBMaWMuklw1IBRBYceZae2MOYR/x8KPIV/JoBM1ywGJMyrQKn2S5tZ0FIrv4JXfNXu3UW0yqsQOw2TLEDqc6tpnkTPetazWsbGxpZr3wLiOwUIqGC2KyxDWJlXIJI6x3YZBOdePAizCw3AaScOcr4TwkFEnXC8phm4b/p+ao6ZQMLzQOxbtR8K0xkJsq5MGHah1K/bb7AfC5D63HANkGKxOA+Xy/GSBgzT8xo4INjTgeECGHrSIo6xTJDQsWHjHFOBdoq1ALUBgXAsgPBYPJM3rsf6VTJmHPqlAMgIbbTAlFuC5vu3AC7IcXLjJAuQ+txr/XjHVKcePOlzeEpHFEANAFx/es4ndwrEnDcI4VatOZgTGJyXpSDkulEcZ6Ldw3fVHQo6BweNeUWQMOGYAUFtZm9Sk+A8goNayoCE/SIaMC7C5hsyWYDU8bbDX4DlFMsaNQcBMJ0IjBo4KLDYZnyV4KDw60mLBiTYb0ACWED+Q4E5pDWLnq81XYPw2tQeqKPNLJZBEagCon9DgJg+JYGL5hsyWYDU8bYDHJR1LdQU6KQgmrIBCH0N/S4I6vEkrH2l7RuCxIACu2vRJoIjgObwkaehPVIgnl+rS3CQcFxTFRhBFRwmMkaAcAyE/SlVCft0iJrXa7RkAVLfO56Prat4gNCAgoLIl6S4bQCitQa2teZATkHX/gYFPJHoVFMrhQAEHXsfoCgTIMhJBAnPI8AMQDh+ogEFcPgkbJO09kA9AsSABB32R5TKJy7ZUEULkPrcbi3VGCDcCz9kHGZSBy4bwYRR5ilNgBAcvCGkmllFAYfAEiBTtAEEGiKvo046KoY6BiB6ATmARLeTAIhx4g1ANEgIjipIjA/CPhXQPwBWjUfR+AHP24tdshaXY95IyQKkPnebsoxXMtJ5AOQQnN4OPvaxUxWxnxGomQBinvwM4VK4ayAhOLCP0spZvjWAYD/DwhogyBnx0uYZ6mntwyuCCBCtdZAHVdJhZZ6P41qDoGYJF56Ioicn0O/qtfX/geYaJlmA1O9WY1rTzhLmYv0SGuRsCDXl26FZhaf1UdpDCzT2a8EkOCDsWrixj1KqCQJNgHBCY1J7GIAYM00DjW1VAaLbQZnA0loEuT4fddgfahDYVCHA5kyI3LcR/d4MpXZpg4V4wQathZnbNPccoKzzyX0Psjeay1HQqUWoQYxpZQQ6CRIKNcc1+OTnOVp7oFyb8UuQgDQ42BbKbI9tECg1oLEdJK1BUDYAYTjXhHfZHwBE92tM5EcoykNxEyw2VLIapH63W0smNMFtEM7rcVnKrU4Mp1IgzROfuQEHK2hkIdfgmAYSCjZn/BpwaD+mCg7THtsy7VGLsB0mrUEIEp6PbWoRo9HQH/dwGPpD6fTtunKMyWqxcTILkPrda8qy6ujo2DE0NPQgxkTOhaBynwYKn9pJDaKf+DwBxESR1kSBBgAo1Boc2E/hpgaigNMH0dNRsD0TQIAOHEFiO8g0SJDrc5FXARLA/3AAkHtv3L37EZyBVhvPQQc7rIlFJtQxuTt27Cj39fV9GZrkMwBIVVpj4acWqT3pE53SwMA286RQG99DgwPHtPZAkwRa0kEnUDTgcKwGOJRNe+adE+N/MIKVRwdHRL6EU6PrYjkhDhsuWQ1S31tO+caKPe4tCK1eiyXV5zu0eKqKghoBTnEszIl+GYBocGA/c9blU1+DA8/3FMoaINivAYKcwDDm1fRQsW4T57Ad3RbyKkBCAMQdCIKBXzU13cK+AVSs0pDJAqS+t52y7ezZs2dgYVP7z5raml+T9yRMYdCcAsvExzRBwmSAwfwocGAf65IMODRAaGJhn563xRyU1Eo6koV9pm0DEPpBVfMqnHAc71AQXH/7rl0j1zWw9gBLrIlFJtQxUVb52D+7tRRdtPKwF+3qiNz92YrkQqdm5FNYCRIjxEZjGGE2T/waQFCX4NAAgVbQAMG2AYfRIsa8wiHdtmmX7fCa8IOColIeFrJ+8JZnPvOm9QMDzoYG1h7kk9Ug5EL9EmXVbwsy6+arXLYtTPuvGM5693Tm5dF0XnKY8kShZUqCxAiyAYgRaJPTzOKN1ACplgmOJEAIjukASbaHsY+IgYIjyJ90nHWybZu/I27CdAlHGy9ZgNTvnnPOOJf5aW+X9BWdmJyeUimnSWXk8rFOub31iGxNH5JWLJMbTyKJzadx9I9ASAqzAQZBRHDQNKppEJQJDN5YopFlAw4DEKOZjCbi+QCHn3ecFLTHtT8aGLhnLU7d0uDaA2yxGoRMqFPi6ol+zkldOk9yWNcz4wMcnuekMAvXkz8sn4FZuJ78uzcgbQmQ8PGdByVBQWCY6SkECCmpQZIAOZZ5ZQDHtqA9KgWAYzAINv/HwMAnVqM5gIOXbPhkNUidRKCqPZxWyby3U5rwl1b4w28GnzVIYwauK3/sn42nvQOQ7J+iSdhFPuUpscyzIAq2AUdSe/CGkpLmFUFiUlJ7sA20p8FxKIruuC8Mr1wPxbPBag7DLqtBapyY2wK1B17XcF/WIU3nN0sWX6ZNuwRICuAgQATaI48+/FF5OaaeB3Jn+oC0JDQJu0dtwnlb1af+FC0yXYMkAWJMK9MG2wHYOFs3oFl1JAju+JXrvk4GBwsAB/FEHNkEDlgNUgcxgD0vMFlgOmWv7sLXBnOSDjP49GbaofbIYPGFNDSHB6Hl55or8qbxpVJo8uWe3GFpDhEHniav1CQkTio02sOYWUZ7JAGCajpR6qt+h19GtAqDgR5m635uX3Pz+wUTEnGY4Ghop1wzKvFjAZJgxhwV9bdA8Lg/q83PvqYFH3NOS8aLwRFrD3wBARBABAuzc0sQzxJmoKwZXSCDTlH2ZgpBGiOKEG4Kb1IZ1Bx3I9nG92BuAGJOQLOcPhwCWG6AcY5SGB7BW4TvOzgw8E/V/5tVLTiqzDCZBYjhxNzlFGy/xW+6qkta0zk45xk45ymYV9QeDs0rAISSyXlRFRQKEOUiys8fapeh7rJb8gAP/alzrTgoyBosLFArJDSD9lGwSyMJx9FKTehdvKfu4v0PfpPwmzjyYYBjD6qyLdMMT7UpwQELkAQz5qBIGaZV094m2SvaMNIB7eHQMU9VzSu4ADjs6hm5BEcJQCiCxnHaSFQpZQvhe/Kt8mZP1KuUozAPkbJMqY8nG3pau0SEDz6AXksMKTt4qYrTRPDWLSY3huEwXpT6NrZv6O/v/0W1JhUN+2fTMThgAXIMxszSbgqgn3Eyl3ZI68KcZKE9st4kOKg90lpCIzjm5So4sIShX8Di6kektPXhkf4vy4h8ecmSJSvxQdpLIOyvwfStczIV1e24jhpyK/BelJNRHr9Mq1UB+04ggQ4AKPdi89s477v79+/fx2NI7BeVlgUHuXGcZAFyHObMwiEKIEK72fe2SQvCshmVYlgX5NK8AjiM9gjgd8TaI5Q8Pjo1JCU1pCqfYYVVssrZvnf7QyiTrm9P5dY9O1jwhXO9Bf54q+PtSA//1yPO8H9gtkqxKXQx1hc8iXqPtLS0PIa3GEdRNonAAGYsMAxDflNuAfKbOHTyx/VTGqHdNW3Sej5GPhDazboaIDpylcFzn76Hgyc916WKzSv4H2ERa5QcjgoP/Tx84i5IM0yn7XzaO6sAlu0ozA+a3r5AtQZdTkv0/HBe9PrS8jMnJsYef6I09H/+Un5y+7Qu08cgWY0xjTEnsknG2TQ3HOCTOmqNcn/VLtAhGPtIqSx8jyyc82zsnMM44uKhSd8DCyaEYxjpGJHyJpxfXiOrzVNfAI5KKpVasVC1PKfH6XRbUi0pjICHg0Gp64j4lz3ijnx+mSzLrpNVsWMT++oEBqPCzG16ihywAHmKDDvB6vETOyUrm1Xzmpw0Y8gPvgcAEkeuoD105ApjH7CLYoAwchVEBYwnHpKJ4f1h6Wu4ltom22imMVHA1XmV83b0BZ0XpJzUp7Nubl+banYDX9Ref0IORoWP7JbdxUcxxIi6PI+5TafAAQuQU2DecU7VfG0JW9a1S0c6K00BzCssqEDfo6o9NECwEiJEmL6H9j9ga+Uxn3FMyl/eI3uG1spatmOEXK2FqXSJXBJskft2fKJ0x19kosyzxyqV94/6lYcOSH7gV+HQd1A/CarjdNEeOhEOWB/kRLj01OowmMSnd2dOclc06wkjWddTTRoc2jnH2AfekdJvBPrQCzp6pcc+AvewTJSHvJHPsYUVssWAgz2IMBqPvRtYVjfKq9N/OvLVYZRveLWc9bn2THPPo8HhMR4DJc/Dpk0nywGrQU6Wc8c+T/sMaSd9Rau0d2YkF6RUE4KwxvfAPFyjPTQ4MNWcI+hR6BepPaLi93aU9u2E9nABBW1WVS/VeY7b940r3NWXYNu5Rn7AAXS5Wl6d+YHsLN1Suv8JbFpwkCmzmKwGmUVmVpui9ki3SOu6FsGbH9LkaN8DJpaD6JVCaJfv/PF98gCDfFxRHdM+MHruO8OYrjiiJjZO6xIB57c6uTcukK7L29zWyz/R/Cf3daSbv1op5r9xzfhXDhEVHwFoqoCadrrdPBUOWA1yKtw7+lz9wEFo96Jm1b4yg6mG0B5O0ryKQ7sEiCPlIPY9oD2CEkK7w1H+/geCx7ZiiE9tkS0Emn7JCnlTl2r7X53SXkYUrJj1mi9IRZnPTPjhjne6F//tUkSurrNm1dF3Yxb2WIDMAhMTTdAkinJR6zUt0oFpJblwKjgy0B54Qaoa2qXvQe3B6NU45uaOSf4mnB+ukTXUGjpROyB1dkUdE/OdeemOVFc2VE44BPVRCP35Q9HEGxC5Apis2xGzanZ/LUBmj5/kJT49Kytzqu2iDMY+PGnyDEAcNztpXk0N7YYlhHaHZHRwLJzYjDaSUSgt9QBJ/73hR5/neem1rpv+ocJrVmHg5A4HBUFIeD3OqSRBNXv/km3JAmT2ZEDzMhe2wffoSqelOdDggO+hQ7sI8UZwzvHpTR3aNb4HprZjHaqympDi13bJrpHViYHBateil8pqT8na8AuFLf/84eH/fXHkR6uGK4UvHYhGHvhluPdW1kuMl8zef2Rb0lEPy4ZT50DVEpKOLrf38R45u6Nd5kub06lavXbBiLfkvCaMg8BFgXPu+yUpVfJS8GFUhaPRQemvDHi7Vj5eevxxdIVt0VQ7Kq2XtZi8tSLcIBs4Mi6Yo5XaLtv5gqFNc8QBG8WaHcbqSFPaacK4x7zOlLT4nsp5rhOPfTB6FYd2+VKUj5FzLBIN36McBYCKj7f68t8DOHRo1zjn6BbNq44znTNvX6lW/nhJ05JNG8Y/9zC7yxDwCjmoNsg2goOAsg4IGTMHyQJkVpjKRagVRszb12Wlk74HIlc5mFZNCO1iiQWYV3ROOAU9gG4o4aeE6e1FhHbH8ObHhBrToV0MBJqkAZdDaLdLda9qd+etyji59143b8O/lSrFTZ8c/eQPqxUtOAzH5ii3PsipMxYPGXzJCaHdJtWxMi0tiFw1Ow5GzkmKAMHIOb5wrkO7nHel/Q98WrYsAEg0dv9jwQNbGdrFm+vV0C4BJ02dav76DllAoOU9r8lzJPvmsh/e9gfu2rvPkXN6UIeagyCxaY44YAFyyoxdr0O72ajjmiaZF3nSGrpV7eEyckXzKhHaZVi3jGklZXja8EKoP2YI7WqZT7dK149yTvtIR3pBLu20qNFKvjBaKfij0UTWEefwelnP+2fNq1O+h8duwALk2Lw5kSPg3wYObTw7ozouSmPk3EVolwBxYF5p7YGRc609dGg39j1KoY/Qru+OytDgk+GhmUK7QIgaeSS4749b3fbzxUldi1m/u/zAbZqIit4RNfqJHbKjvFW22vt3InfpFOpYBp8C83Aq+Rdlw46rm6Q77Upr4GJZH21eESAI8dI5nwztEiB4tRaT3EtSxLeV818bmjm0C60UYfbu2vQ3i/+4+7PDH/5E5Yh77nh54m0Ho0O3bA9+8l1224Z2yYW5TdZ+PXn+AhzwukUtaHOXPjpPVra1yCJpduarZq9Tcim8JOXlENpNoRq+JuuXEdrFaIc/ismJQ9GQ7PYPe4+u2FfacVRot10WLB+Rg7tM166WGzMb5Ro9OdHss3l9OGCjWCfPZwBE+Wmv9YpstLDdlXbfUS2ecnIwrSbNK04r4ZojWINKa49K5AMqFa8o4/8KcEwL7WqfuyPndvz8fHnVzkXZpZ/taPO+tXHgmjy7uVY2I7qFQHDszJ98z+2ZJ8wBC5ATZlWtIrQuX2R6D3a8LJ0KW65KSyec8yxCu3wpKh77UAjt0jnXH8iE3Ot5V9AkGDl38px1pUbonOsVF5kjQfiV3+S0Xd4uvZ1NbteFOa/r6245t+vtLZ/58uj48Nc2y9q9GFFHVav4Ncfq8EMb2qYT4wCkkk9wRo0Yjn2Zj9V1L8m6i89y3YVhmOlxKuluKaY74H43I4CbwfRDD2MdjuBza1iUoSxBkA/8cNTJh/0/3uv/7Mer5L9gf21hFAxt69Buc6ta+FcdshQfRmguR44HXHnLw9D7+GF34IGzMy9ajqq4vo5enVivba1T4oDVICfEvrU0bQCKSzk+oXp7P3iRSNfaKFJvUNAeFWlxKk4b9AK0BubwuvhzoFO4nCif9al0GfvL4vrDWO3tyagYdPxUxm7HIgzPq04TIfAUgYKFSbs/knJb1uXS81+QdedJoTISjFaG/PFoeMfO0t17GNrFVBNdl/VtmlsOWF19fP6CPzSn9ACe29f3gT+Mos73K9X1XKW64aJnQTSl+EEC5lyEmi9EcUE4KhsSvxpVVdR8+MPocqJyINHYA45T+od9+37nZuwE8CJUdphjU+Qt7f94UdZtvrJQHH/9UHl/y97wZ5c9HN6+GZMZPUSv9FwsXdH+zCkHLECOyV6aMfGTuq/vba8VWfRxpXouAEAIDEhxC4S5CUKdw8KeDOcSIPFqJfGa61TOMUAUYKPFHkEvLhgaRViGBBMXlSqhfOQ+kUN/vX//6u+yK2fJ1ZnnyEv8LbG2kjdk/vYZ+4OH3vKE//MbD8mOcVThPYtRxBNsmlMOWIDMyN7YpFq48Lxmz1v1aZHed4kshDA3AxRt4FkzwNMMAW/C2fF4R+yUQ4Po1RJTyOFeVAFCmWZAmCleEpQgwZTFqMLlc4GiERx58vPth2/+5I7i9XtYb7Ws9xZApWyRDbDPbDpdHLAAOYrzMTiWLHnZM4Jg6b84zhnnh2EbbP521GwDMFqQI5SL0fKY+CIUidojfuec4FAKVROJAJlOIcK/UQQ/PCohNOw6bunnI2fnH7h5kde16RsjV9xrTidYtlWnuJt9Nq8PByxApvB5NeyibX5f3wsuEln+TZGlC6OoA/Z+B/a3omYrBL8ZeQwOPZ2EExIBEM7ajRdloOkEnVEl03wSHBgSwbT3ScJiihgrKQSB8t2uYJ88q/JwFJZKtx8uD2w6UL7vtp3y9bHYsrKmleFnvfKpj7l6XfVpeR1qjm3+4sXnnRtFC28VWQBwNMGkygEcdMIntQSBYMh1sRi1hyV9NHmSgmWF1wmPogxcFO7ncUMeWiY5DqJdTrPrShqfYO7170s/S41H3ivLwfCWYurIMoAD+me9fZidBrnB7bEJHMCDYnO4ePEz+8Kw59+UWtRUBQdAQ3DEppOOUmnzKQ2hjokAIVgw3R1EYYfnUc2pRUwyGoSaA8MimliXZSbWDcMmlQ75abQl/uPNvsoW7rh1X+X7D8I4Q804DBzXtr/14gAY3/AJorkWpFQQdH3DcRYsBzhgVjFCxciUcbiZ03xKAQQxYQEFgAEf4kzFmmO69qDWyAJfhrhtyGgTo0UMqATmmhdlVCX9DHd44bt5f9C3LSCbTgcHLECq4xw9Pef8het2r8ZTHIN3GWhW2EMYCTchW4IDn/arkYtvm7tuCiYSKTabjGllQJDMjwcSAssAhFoF75C4nN/ouWe+rqfnp38WD1DeZbX9aUBIoz+Z+ICIliyZ1xMEyx+DY4647XzwpAvUgUNt0Bh0zFs0OQ6msoNcTCUheR4J868g4AYk0zUCBZ6JJlbSvKoAhobKCOSy7ENvkWITrIyIcB5nDYylUoPn7N69ZlDkOvTLjqLHHK3Pb4NrEJpWgvm1Hdc7TgfCUxkIZBr7+LCGhVUbx4ATjbAtCd/B1BT7HJ4GhgGHcb6TmoTlpCZh2Rxnbs5hG0aL0B+BGYdfN3ScRe3lcmYD9Ar6thL7bKonBxqZ4Xw4wDHvfHYY9v5C5EzwYgH2dWF3ByjWHgztOk4LiNqjpao5WrT2SGE5n1TK0UKeFHajUag9tMmElpMahFqCGoOag1TCmx6mbLRKrEU40j4Gt/4A+vOr/9Hf/yZoOeuwgxl1SxSSRk36f4fPcTWEH+oiBUGk1uBuPjdimoxEYcIINgwxNIuVGTQA+OQnGfPK5EY7MJ9JaxggGe1hNIgBlR5k11qkG53qfFt8o7Y28j2LWVDH30ZlNqXf7+npyUWR93thCAnWoj4VHPF9ICgAlwQZAWZOMgBJ5gYkJ5Kb80x7pv3YrHMdjrjDUHtrT88XciJrEBjmCig21YMDjQqQ6v/tvwS2/mI8+6E9KJbJRKFkSuaRBgr3GcBMCnMMoqSQs5wUfpan7zPbyfNM27G24pEgVKp9KQIIL8JV0KEt0/oa99T+zj4HGpXR+gkMAXxNbMY4AEgyERQxqXiKOgBBUHBfmChze6p2mRTuqfsp5jxmclM228xNOdkGzsIVHAAkh4mN2ZfrC0q31SAxI+b8t1EBogERRc55eOkJTDYmy9HAICAmKYCQm20OgU/D1XFuF510kwiAZEpuTwUHa+kTFWK+IPd34vO2nviFkxey5afMgUYECMUzXLhwIcK6ck5VcMEHCmJS+AkAQ9wfA4LfNGfZ5NjQESq2kySOeXDbTEo05WQdls35cWn6LytoAjh4XXVmT896+CEcCzGgnn6O3Z5NDjQqQDA2kVkARrZPBcXRQIiB4WtA8EWnGBzMuc8HAPAKFE4zZOZZcduUTZ7cZ+oznwk08U1mf2jShQA1r6k6o6iLwLapThxoxOkL2sAJgmAhQrWYlMgQUYB9Rlswj4FgwIBZWtV9BEVFUxhWAIoyCCsnIgpmQJD0JSj4xq/g/eS2qcf8eCCJtUutH+hfJUJYucV1c0vR1KGqo86O2TSHHGhEgGh2wkEP+HZfDIxYI8TAmKkcg8KAI84JDizrA5AEAUfXsXoJTjU+BC9CADByZXwMAwiCw0wrmZxaEtdnl2Jw0Iyb3pcQIV89WKP/B/sz9xxoWIBMspZmTFIQK9Vt5nE5fod8EiSx5iA4qDmKAAcBwhemJsFAIec2wWAAwn0ESRIgLJMMeJjHr+Wa6/O67J/ZZtmmenGgYQECIeSzHXyGdGK13EkBNIIY53glFgLOb3uQYq3Bmb1YqwqCzenv8dys+BVbrGhS1QAUdKM9ZgIIQWEmKB4NEppy5RpN9q+EuVkT7LBNdeJAIwIEIkzhdQcrlUoBphZm8Jaxrwy0GKDEOU0ppWLNwWMYda+SC4C4GhxBwEmMBBoT3zEnSJTWCr8JIEkzy5ha9G2iqAjCBC3dH/aF/fNVGE6MBcHhPbwS3mGh6rNpjjnQsAAplUoHYRqNgL8ACIXRkBZIbJuc01DIpqPBASMK++OpKChAqBnRCkF8kcrVICF2SNQsTMaUMlrDgCTAN9P5SZ3YfOMkxZiq/aoCuDCk1OGJuCX7Ww8ONCpAnMHBwYment5H8PRfBE2Bp3ERJhdBQaBgTrrOaVbF4IhB4mKb3xmMp74bcFDIY7+By/kEOO5D2PlKbmyCJW+kAQjzGCQhHPYAZRMVM6YVtUhM6AtG0bl/4tcDA5vy8dKj9hXcJF/nqtyIACEv+ejH9A15AABZA6HGE7qIXSSCAy9q4FXbmIz2AH5wGgHCPOTQhE40q6geCA5qDx/AIMW+S/yaLs8luAgkhTrUNNQYBBPBQUARIIaM9jB9ygMgpIn/jK8Z979attkccqBRAUKJpmB/D/RnKELqqTnwcJ4CkElwxE64AUhsKqEy26gSFhjFyqExOOi7cBE5ahC2EYOLGicJkNgci0ESgyPpnBMkBAj7lHei6DBG04fuwAbSDt3/uGx/55IDjQoQ4+D+GEK6D8K/GE4w9uUBlCz4bTRIDJBJzWG0RnxLqAXMBEaaVsa84qIONInM6oqMdBnNE2sRAozAosYxADEj83TSjalXwIXy6FcBADm0R+Sxn2IHOqFXhI87YX/nlAONChA+gb2BgYF8b2/v9zFC/U761tAiAAh9YJpX5snPpz+BERM1QJxizeFgIjA1AUfbzXgJx0fi1RUnAaaVlDbRqHUmNQ/PjcHFaSvxSD19oVh7ECDjDO1igHDk6wMDt0KdrEejG+xgSPUuzHXWqAAhX7UWwQj4RgjpnwAAQAKHGAgQapAkQICbqnkUR6SIr5hicBAgXJA6JgJsqnkVm1gxuOJ1epOmmdE+DCtzPCbWIMa8Ggc4+tGxg1/FQSbd77hof+eaA40OEHffvn0P4s3Cf0FY9i0QdjyZK+DJGPgeC3X85DdaI37yG3DEskpQGe3BBzv9j1gDTfU/CDL6MMg02JLOfWyexeAgSOgPESATvuOMe0FwENrja4/FH/DR3yjBMZvqwYFGBgj5S3GF8xt9FOB4I8rkB/YVgYhxFGPNgQJSrDFMTv8hNqsMOOIBRTPqTnDEIeKpQDMmWgyUOPJFoE3Ou6L/obUHLpin9hjxvMG/wnXRp+vYCZvqyIFGBwjNFffAgQMPQ4ushxb5FICCRzjmlWhTi5qDxGQAwlNIMTA4j4um1eR8KWoPM35CcBhKgo1m1mR7NLHi+WDUHgQIfY8J+B5jbhD0r9+799v9Imtwr7ZZ3wOcqWcyd7+e13y6XYs8oPRGAMmdAAlWV6SpRW3C3S2gVlAzKFelJuSMdmUSxLBuDI5YEfHZQ5qqQQzgYoBMgi0GSKw9lBSDyBlxw2Dvdw4MfOv1ia9coT2b6skBSkCjJ/MojwCOPwI4dmHwkJKNxzoFeAw0Ws3HEznLUymKxqEZSBPV3Gwn8zEcY5tJYvsjoGHoq1HMCtsd5Srl4KxDr2QHovWygn206TRwwGqQSabzUR8sXrz4XIDkZwAJXqbStg/2k03UGNQk1CLUIEktwqiXodhBn9QebJbPIVKS3dO1Bx3zMia2Hwpales9v/QOKYyUZF/lrvO2y/UP2o93gj2nIfFJaVPMAToCHqJav+zr67sE5W8CJPhGCF/G4IQs+gWsQj/BEE0iY2YRIMbEMiAx4GBOcBiAUCEYgNTawto+R9wuv927QN4ctXmdt5e8vZu6K727cX21QS/3g9NsqisHLECmspu+h7t///47lyxZ8iLMlcIn2Bx8gg3D3Vq6y1ADrMKnPSNNM2kRAw6ylsAwIDHgwC7t8LNJtsWJkiVMURl3VWl0pGvi3JvdlLdp08gF97Im02o5gk+w6crxDvtbNw4k71rdLvpbcCFKdcCVT+CX3AiQvJ1OdWxyMdzKDxCyCrUHQTJdixAcBiDTzSujPTi1hR/xLKGhEdQf/Xx42PvkQHHHHmwAFHfhI56H8BHPS6mmbDpNHLAAOTbjKdm0gwTTUV4Hc+ujoAuqQIGUG//EBQ+pNWhiESjG1JoECFbwxX6yOgZHpAcT8eFOxensE/fB2cdnoAe/iwr4DPSr8RnoK2ufgV7X9r1nHMg/+pad/t037pAt46hiGmJ1m+aYAxYgx2cw+UMiUFz4Jn8IgLwfQHkuiBpFUyz4rEYgYGJizSmPgYFv4aAV1NXAQAQ3cgJsPoAp8f+wb9/Bm3ESKsRf10UZScna9r+/qKf5rCsdP/P68ZGxlicqP7vs9vBvN9sv3sYcqtevBciJcVqbXNWqChrlIpTXYjT9DQDFAq0hIPMECt70EDfKASLUIBR1V5qibkAMn+j0u3wvaHZLwcj1D41954O6QvzD9gnCXKvTvfYM9YJ1y7KrXnBm67mScTJBcbTk9+cfvnd3uGf1JdITbLAf0Umwbm6L8V2c22v8d2i96qRr1RD09/ffgX/qDvgnt7lO6l+aZRnEvtdtUt3S5HRJh7tc2tKLJJdqlQy+p54O28Uv+1IuFFWpPK6ORI+8SOQ7skrWpbbLJnrq2pQjo4pq9KN5eXLJWGWgPFxsk97mM9NdTS3ueKltRXcoSwGOx23Il5yqT6INYNOJcYAOhAEKn/iYRBjc6gelndAITlCOQlXC59jKHeL4WNDBj7BmfAoeCb+CWxLPLYmbmnCVNx7mnOxLXui97yUAR2WtbOY9YNtsc6ISla4fkt1RMTycLvtj6XIwuit0C3/d5i467/vyhV3QUQj52s+wgVd1SVaDPHU2G6CQd2Vx1Bcq0cjfVeQwXpptcSoARAXA8PWaWZiTBfH3uHAKxR8f/oycVOiGbU4xan0v9vw4cfkYfKF8c8I98tHRcHDnuL/os786sO9btwrfQxd5sWx2t8gWFlnXpjpwwPogJ89k8+Rf4LjOo83S29aBGFSr6lVtzkJpSy2AmdUtrekOfFAaH/3EMlwBzKxicVyK5bFoNOz3n/SeWHFP6Us0maZohYxklpektMt07eqzbsxs3HlNyWzbvH4csCbWyfNaR7Zw+iAc881lGcEKWiOBj3lWPuZi+eEEFmIYx1uCfPjDvMKM35Qb4QPTDhxvN2hXnalWv+UqHIy2xuEvFOMUgyNSK2QtY8ZCcPz94s1NH+3Y9MdXeR+8+Xfldzn4wmQfcDEf5uzXAuTUWEuQKMj8Rl8K5ZIMuT4mMPqYsOiHJIIkL1GI5UkR4k3h9dwMhk0yruOklYNYV9MVqxDQ3SYbjG9jeoP7oiKMe5Rfk7182bvarv1QMRj+ZWs699WF7sLLniHnvpYVoXlouNk0hxywADk15mqAYAnRBzGv8c4SZuOWZdTXWoQAgQYJgolJLYJVT9IQaWgRJ+2ooEM6Fy5wui5FFyKMbxhhp1agn9Pe7nb936HKwfsrQeFvKuHE8sDNF3Ip129zWj60NtYuvL5Nc8gBC5BTZy55iHHD6MayTKiSjMBRh4lVNbN8AERrEayUqLWIG0qWWgQrRWTx/kiTZOisO1vluumOdzkv4y8dloPtT5b350dKB6Mm12lqzqa8NjdX9CU9rxrNsmbWqd/DY7ZgAXJM1pzwAY5jAARyZxj5DxVlGEvGjYUVmlk1XwRmVlCAw1BGRCvWIllPufjEetimWs//Pfeta7CAaQStQC1iQr6FSlTecFgGJR+N5PKVEX+iMvzPkoouvr7yiRd+S74xgLpG27ALNs0BByxAZoepFGxMKlSbSvBBSjKK93YnEO6NSWsROutYCM7Rzjp9EUeyDgK+qllaotzVcTfWmt7E2iSUb41HY9ufjAZvyFfGzv/LwQ+svXrg3T/cDCDB/2CYmWCyGsRwbQ5yy9zZYarhYwdCvo+3S3dHtyyTTrVQdTrzpSM1TzrTXdKebpecmxMPgyNBuSKl4oTky+PRUHiwst/rX3lbactRId9k9+h3rJAVIUwrrbXWYSR+EwYbk3VseXY5YAcKZ4effJKTl0MY6f5aQcb/rChjfjlq8cpRE6JZoCCLr50jaouP7bjwPVyEfAXhr9B3gq6oLV3whxny/QsQtdEU53u1rPa2ylZM9sJyjUjf7P30c/L5/HufGNvz/BcGL3zx3XL3GHZbc4vMmeVkTazZY6gWas+J8Egvlgsy6sJpxxsf+aqplQdIChLCzFJ49wn1dMg3q0O+KmqR9BWvkFe0X3d0yFdtw2omGEF31rqvfPNft6+7rSD57e3p7Nt7vfbznis9l/BfWC+rCSybZpkDFiCzx1ACxEHI96EwCu8sYFEG+CPQIhNShv9RQTTLx5hIgDERuiuuDvly4FA5GVhcXdKKkG/npVADWKRhSsiXPex9l7Puv8ZlYkuhMvHKCX8kqniFfEc6I/PczIZ1sip1HTQMK9o0uxywAJl9fipHOTeWpKDyMoZpigCH1iKxJgmgRczAoceRdYR8swj55rCQSrOkdcj3uqNDvkNjarx5QAbLhypPFkuVCadJqVwqrZ7sdDL/ryQtVnvM7n2stWYBUmPFrBSqId/gTh8h3wkNkImwBICUCRJoEFKsRSq1kC+0CEK+Udilms9/q/u6NTOHfP1PHZThNGZxZUfKw/eNVkbfH6XCFX9e+cpfflW2FQEqEyiYlX/ENhJzwDrpsy8JfJrjg4dqU1HyN2KwLyxJs1MOszC1sgBIBk47v2OIj3/iGyLaWYcWwUfQQxVmnZEww5DvnQz4bon7Vgv5Drv5lx+JRv/p4xM3fx/ujd5/o1yduUY2lhKDhgwY2DRLHLBPnVliZKIZw1Md8u2Sro4+6ZH5ap6a73TK/FSHzMMM3450KwyjJoz74YtTCPkWSwj5liaiQ+GRSr93aOXXS7fthC/iVAU/0bwuqhvl1elr5Ad6hu/VeI+9VzI9H5R/ewJHeX0LEs2mU/+xJtap83B6CxROahEd8s1j+gkoiM2sArRIAVqkoM0sfqgTU7L0LF9EszgFJehWuXSbn3o3G10pOwzYuMkyBwh5zyKC41vLrum4re+aP39t+xnbW9zw7t+XZ7Xy4qDkedhj08lywALkZDl3/PN0yDdyPIysl8vjMu4WsfBcKYoBEke1GPLFyiaY5ethli8nMQIgbhrWVqukr3ytPLvz0vhLUkbYKfvhrbLBvUCWrni7+zt/t3do/4P4EslnmjLeyh4v1/M8Wfg6VI6usyHf49+dp3DUAuQpMOspVNUhX6lUGPLdOi55hbELnwApQYOUoUGMFuEqjW4NIJjACJWyUHIdy5wFV+B6CPnWxje05tiekhX97tB9Y1L4wLA/svhA5VAQeJVoaaZF5rvpj66XZdleGVeboW2sJnkKd+wYVS1AjsGYWdjNJ7+KlHN9kSMiMuEUqlpkEiRFTIWnFsHAYS3kK04O69K1S3rdWlmRro5vqPVobB0W4FpVkR2Ho4lf7MYrvocro5WJUsFJR86RKJW6pSvoetcG2V28SrZXMIdeL7MCkDh3Ycwe59t7fRI31UaxToJpJ3gKo0wK03y3+k54/5jKnw+ABMUo45aiNHyRNLQI3l13YVtBdPGKCCJaUBko+EEUdkfZlaPuvJddGqjbMRDobIDQoz1tuqVc+dITwdiFZ8todNgN1H+qh39dyBw+6C2Ri28M+86NguBJVHxkpKXlMbVz56g5j1rlIVxiw7SpLCf4/zRkNQuQub3tdNZ9hHxvmpDiJphaMLKyMLMyIICEAAk9+CKO/rYItQgmakmIWb4qTKmO0Hv/TSL/LrI9eN2SJSu7ff+SNqVek5XonEIliPrcEbwbcgTrzavntan08xgR41d3I8/DuEskfj5/4Jbe3nuxvum3B0W+e+n+/fv47xIoCCPzS++4oE3H4wDNAJvmjgOGv21c2KFb2hacIfOjHtXuLHTaZEGqTboR7p2XbpZWNyNpHfItI+RbkInSuAwGo6X7W/e+p9R65M3dofeqRRiib4Xg59DfFhDWTNFh4gzG5lMQeJAOn0Eh8bqOh9e4XKwAydHLkTAcxiSXbw+J3PC2/v5fYJcGCk0xlm2amQPmBs581O6dDQ5QS/uQ1L/LSeYDZ8g8f7F0eIsAkIVuqywAQLozrdjTJE340C7sKykBIMViAV8KGZEH5z8i7V5ZWmEzARR+M4S/CcKfRQ6a8avuVFvVGwvIwc3BNnUTFuGWI/h61hgmBO+Jog9fNTCwB/vpm7AeijZN54AFyHSOzP62FkAsl/UM13cf6sWCQGdKl+pRbbIIAFnoxQCZl85hVa00NAJMpFJZCgDIT5t/IUF6OJiHUXcAwwFAsE7j5EcXCBBMoK8tl00NQjQSILwoby6Jkg81wdW2+W05DZTBMDwyFEXvu2Rg4J9YByDSX7hG0aYEB8hHm+aWA3yAOwhk7UTI93ujWFwU7xz6BYyBFDGzl8SwL1ZohCBjlq8K8VKVKw+3/ErK6SMwyDwXwACUBO8eHv2lxBzEegpBS+Sq1Iy8GcdbQPBdAElxQZKBFulVquts1/3m1r6+f3zfWWdhQjFDyhpXc8uN37LWLUDqeMPgQmzEdEV8DaTg5IEYgqQEgJQCgAQ5P7BLjXCP96Ds9/ZIN6JdTYCN0RowrbRJpXMIPWb0SrZKLDcBEDpHuQYaAxjsa0G5FTnA4uWwNH0OQDlLqXdfOjHx3U8vXty0AYrEgmSqQFC72lQfDpDXSjnOvfNU0/lnS2ewGLOxeqEfFnrNsijdIj3pNnnQ2yl75AnM3krD54i0M57UHBoUaIjAyCDPIE+TUE4hJxkfpHZzYbZhdiRqMMMyd6AKqAQqcKTScVJPRNEdP/H9139gcBC7rE+imYUfq0EMJ+Y+p9yGENObxmFKjUgRE+DxtgjeMCxhsDDAyr53qx2yE+CYB3FHZEprC/oZBIL2NapA0IBAGUsHTRLMsiwoBXJACqQHVkyO0C+AgG9jYbwFhHdQtPkFjZJyobqgSV7+Asf5MkAVbsWZuKRN4EDtIWO5MeccMLzWIV98HGHB2dIeLcVU+MVwv8fa8HpVdkh68L3QVoCDvgJ8D/11dmoQmkw0rYwJRQHXWqQKFIKGH/UhCHRimQlaQifGslhO5ihTmxRB+Siq+NAkDwfBtRcNDHwCNTEiY0PAVoPE4lOPX0oqn8wjIRZ24NfQx6IS1l8sy6Pth+VQdhhwgQYAOLS5hIqMSmmCsNdyACAFSlcpAw2RpsaghkihlslZJqXRmimb3GgV5Az9VkHnQZMEfY7zsdsXL76Q4OCAIrrQ0MkCpL63H49xaG032pSP/OIRp+we6ixLOYtxDoymazMKFTQYTA5wMHSLQT8NDO1nVEFCcyoFQGhwEBgGHAYUydyUTR3mBAraotkFkOhxlfnY0xYEG9cjx2UbPlkm1FcEYoCU5VfpTOZONa9JpTISNGPShwGFBgP6xJyPb00AB0fESZyzRSI4qDnob9SAkdQYGXguBIXJWTZkwMRzSWjPQdsACT6wGAVLXPfC1X1977gUWoSmVn1Z9PS6Gu+DTfXjAB9I0fylS3syQXAh/YxMhE/qIJ8ODFY0oDA5tYgBCHPXaAEDkoTA62E/1NfJ+B5YH5VrpGo/xceQoUk8DmIfsiGmgYVh1BZFGzYvX36L7No1iqNsCVnjJQuQ+t5zyr2fCoLLIOzzOWAHcFDMaxpDAwPbBAXLNaDwKU9QgDwAgaQjVUYbJMFCoKCeFmu0oR1zOuckAoPibsDD41WAMEfY2MFER7/XdReNFgqXoeamu+LuJRDFkxojkf821Y8DwYoVK9JRGF7pYUACRg8W+4ltGAMGk/PZT7OHxOgUcw0a5gSK0RYGIMa8MmZUMqeZNd3UYn1zrgEX2uV14OeoNPqHs96Obqg1DRzNsgCBBNQp0ZaPhoeHV0ATPNuNTRqKeU1TEBRHEcHB/cjpTFOLkKaMcyQF/XhASYKG9QxIjMZhu7gOtJsboH8YI3nuvy5bdg76RB3TkLJiTSwIX50SZV+cKLoYS/xgAesogMRhdd5Y8pLAoPmjKyPXidtVqgEEbWgHezo4zDaPk5go3knzymxzH30So424jetQU0EwgnmO4x0ql1+BFh7eGncTFRorVTnYWP/0afpvtZMbKvW7FHsynlSFwNFdgpDqxBxkAKLLRvjNk5+5AYbJjYY4Vs56BhimPV4LZfarChKGsF7KfhwCzHR/GuzHAqR+Nzw4C7NmYdqfyyc6GH8U7ymB08l0D28FaqBorVAVZIJFbxsBNwKfzA1gTG6OJQHCfaYNtEkwIjk0s9DJC25EvxnyNX1ppByPEZvqwAFKXFQulzkxt5sTBpG0FE4HxIzbEFi9PxbcGCgsG6FO5gYARuh5JZO4j2aUOWbOY1tHk2I/ISDzU3G/S+gDa7ErDZOOeoo1zH9e339U8zkIgiV4OnOKFYUMb49Pagwa99yekkMcYZJN2ae3KaYmsXw8Igh43OTJujPtQ7u6czgLGiSCw97SgX7zclvQirlso+Sz/g+vX7/eITUKA5/i/4k5h8ozQCAYZiLaMtxv8ullbrONU0oEiknJMvaZ/vH6CPsykEDN15Bpdk2sSNQGtYH3z6YZOABw8KGsj1D4kjQdBBocqEs/QJdN/eo+1ofBBGlGe8ej6XW4zWTOMWW9M/5h24Z4bYCEmw2ZZu9Jvx7qF6byG1/zxldv/NzGV2pucp9NNQ5gCod+IBEiSXCwzGFqkyfLph73+RBqkxM4WsjpUxyPGMZNHk9umzaSeaIfpj8V7GvUNDsCHGFwdwPWz7xL/mTw4sFNVy27aowMja6Ln5aNytzE/60ZkUqlBiDUeXoNEPRICzsqMTfCmASE2aeBQXCA+CagJgg9RuSPFn4CgNNJDCW3WTaUPDcBEPO2IftEYJRBeF8kYY9hRwOlUzexYFaBX2H4nLBXjshN9y2/L7PkfUt2ca+6rmH5Ol2ENECy2ezhsUolj4M5I4AUQlNmXiMCAb5BDRAJYHBRuDIEHHO69IxeLfTJOVZG4OmEm0RAJMHC+gYszHkc52kQ4hxcI8IeVYii8Sdddy+bWYtaprlGyU8dIMQAb/8IPuniyQfylXyYb82XwUptcjUKI3/D/0kOqUcffXSit6dnF6aMzIeJFOLp7BIgRxHBYMDBnHW4D0RwpCDMJZAHciHcnJtVi2SxIxR2kgEIzqntw0cUa9rFgKRaP0DO99R5DdSihlN5LA+UT6WOsNlGTKcOkBgejuwCRH5ffu5hJo//B/4wmMnHl35yNiJjZ/if6VP7EOT74axfSG+9hB00YUgGJMnyUcCA4BIU1Byc+q4nLSLnq7c1XWHAQK2QBAj3c58xvZI59tNcIzgMASQhQOoUHOehDbt3F3G2g2eh1SC4VyebuMZli6+NBG12nWw7/93P+w/8g+8E6Y+eJ0FiwMEcUwn1InJ8gxDztmKCEHMKCMlB2YFgc/YttQcXd8A8LzySqkAgOGYCCEGS1CIEB4gaqUjC+UW0VSCh3UIY/gRF2YrWkFmAkBknmTjsahiIu2RBMo2Phje3Y8AQsqeaYMJgTRNRFEgKOMGCKfDxO+kQVP16LXKqeYKEM4Dx3ngMEIID+42Xx/AxXnbSxzUwCA6CxyQc02YWAZLQHiHKRezTVAUI/A58m1fcQ2h0xHHuYBOHGtQamA0Ty9wCmx+fAwSI09/fv7enp2cbAPIqRIywpqJ4BAgJ71/UAKJBgW0NDOQGHAQFtQZFn8QnEcGBEDI+6xZy1UT97rrRLDjMCjHhmHHMCQyfmgNUJkBwjNpDg4P+ETp4OAwfG3/mM7dH+/cTakBW46WnABCG+pSsXbvW2bJly4zMygzHY1e42clnV+Nx9dj/sTFTNkH+Xg3B1WYWw1rUIASINq2Q17QGOEnnheaT1hjI8U4s9lTBwRz7OC7iY38Fws5VT7RpZoCEYxokOK6BhDoEB4HB+jSvjImF5X8QbZEQKs4bUeprG7ZtY2CN3WHecOkEAbIZ9whfm0SKwaHj4nyAxemqq3D8i2GpIxrBDrNS+ORxU8/mFDLV2dn53aGhoYfwbscKahGaMwQIzaskQAgMTRBaMlMDpAoObMbaA7kGBwECYODJrwFCJ543V5/DyjiPQKIzTk1DgBAcFZQZMta+B46PoxL64+wPw7G9Sn2VpyLFiIzLDfV7ggC5lODoXrZMKi960Uufs2WLugt+HrV7nL6oP8Yi654rvz+Kd0Gf+cXP794g7xzCQd7XyXrV6g2euTt27Cj39vZ+DFrkZk6ZpVOOp3YNILwpBhzUHKiHmY0xSMg7CjmZykE9khZ4gKMC4gJyGhzI2UZN61TPowbhOdQ2WuMgNwCheTUOvCFy5T0ZBJs29vfv24yuXNqg5hVYph8yzI+ROOmwF3xecA3W9PtfxeIPv7B8+bOe3929zT14sFDx/Sa8tzyMT4p1BO9Sz3vF8l9s//Dj510gC372syLiNH+zfvVqr6qij9F+Q+6mFnHhi2yGL/JuhGpXQwNoLUINMgUc2OYTRoMDAq/NJGwTHAQJZ/Zq7UFwQLi1NjAAQR0d7aqer89FHQ0qgqMKEjO2Qg2iTStgam8YDuxKpz++HuW1Daw9wLrfBBBWYXJwT7A8TaptxPeXfqm/v4KH3jTlg3FXlfLknKwrEz/96cH4PPt7DA5oGYeyfR+e6PdC9B2gJhoDFgxAjGlEgOgEAeZLUzSTaO+QCI4A+/ToN3JtXhEgOEYtwjamaBC2wXNIVZAQIIyeVSNXAaJWqcEw/POv7t49TO2B6zeseQW26AcU8xNJC1DpCOiYztrHRN7NGwL0fGFD9UF3Ig03aB1aQMGiRYve5Xne5wAUjgummvDTAWoDtVQph5zEY1wm1Dj0WKInXt0d+2laacccZT2IiJwXMABBUYPDAExrEICDo+YIquhV3osAx+NB8JUbBgauhPbwcA+Pea/ZXiOk2gPq+P8snXT6IUpWr/6It23bdTM8VTBZkcofN0qr8+M3aI/GHODD3oc/shEOO7WJBgnB0AFqBREkzaAaQFDWACE4UCZIaJpp5x5lrT2wbcwr7cNgm4ngoOrSGgR5TXugWAI4YFr95B7Hufh/7ttXht/Be8zqDZ1OECDkURzmZeFYHFtftbvsk+dYHDpqP/nPB70Pf+QmmFzvqYLEAyBUOw4QJAYgBiQEiNEi+lMI2NYAQZ6MXrFhXkDf5Co4KPXab0FeBkF7+NAcXn8Y/voRx3nJnfv27cd9dHAPZ3gI4oQGS5p3DfY/P93+Xd4DWkIBQYIVE98D/4BRQy4mjU99TjW1aGZNAQi2p2sQAoOkzSvkBiCUeBLtJoCD0bOA4DgQhvfvDYI3/nBw8Ndwyt0t2I8qNoEDFiBPDzGogaSvr+/j6NK17BYiTj6En98oVEaTaD8Ex4wGMYOLBEnNvEI5CRBsarVPcFDyCQyQi3CuHAyCb28Lw8sHBwcnrOYgp6YmC5Cp/DidW7wXpBAgeS38hRvgl5wFk4s7fQCj9iHPpBYhMGoaBGWCJAkOc4MBDv2V2wqugfES50gYFsewQPUPBgY+hVPEgoNcODoZ/h19xO45XRygjPtLly7txKTGawGUd0Ce2zl+gflYAcEBn4Sfg6YJpsz8LeODTDOv6DgCYhocHkfaAQoOBt6KaSQf2r5v34NVYERolmTTNA5YgExjyNNkk3Ku/QCEgZfBeb8aQLkML1otZv8YmcqAAJAQACHV5nFVtQdjiZivqPTSQnz5Ci8+4S02+R6mknzxoYGB29jOaiicbbFLwk2bZuCABcgMTHma7OK9obzTp5Zly5Z1VCqVV2FU/U0YMLwQj/slQIDLcK7xPYx5xROpcQCqA0DZvRjz+NFoEGw5ePDgLraFcwkgXY3bNh2bA2SSTU9vDvAe6SiX6SbAklWVyjJohuUwnxZzVByJ9TDlFy/+BUEewHjkcEvLY7Jz56g5jxEqlm2UynDE5v+dOEAAUMC1kJ/oP7Ye4KIphfoEmU1PkQNWgzxFhj1NqvO+JWmmbsGSsqPhMzHG7rMcsBywHLAcsBywHLAcsBywHLAcsBywHLAcsBywHLAcsBywHLAcsBxoGA78f/cYQKpsN3cIAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "no2_radical = ChemicalSystem()\n", "no2_radical.add_atom(\"N\", coords=(0.0, 0.0, -0.01857566))\n", "no2_radical.add_atom(\"O\", coords=(0.0, 1.09915770, -0.49171967))\n", "no2_radical.add_atom(\"O\", coords=(0.0, -1.09915770, -0.49171967))\n", "view(no2_radical, width=200, height=200, direction=\"along_x\", guess_bonds=True)" ] }, { "cell_type": "code", "execution_count": 4, "id": "b070d1e0-7381-4f69-b773-2450f47dba63", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADICAYAAACtWK6eAABIz0lEQVR4Ae19CZgdR3Xu6eq737kzo9k3aWRJtrCNsfGSOBiQl2B4GMILIH0scRII2BDIi4PBIQlkPBgIGIwBJwHbWQGHROIlgQdfzCoUMBi8xMa2jMCLZGkWbbPfvZf3/9W37vSMJWNby8hzu6QzVd1dXd331Pn7nFOrSBQiDkQciDgQcSDiQMSBiAMRByIORByIOBBxIOJAxIGIAxEHIg5EHIg4EHEg4kDEgYgDEQciDkQciDgQcWB5ccAXsfCLSFGIONC4HCAQNovYW0VijHlcA4dmSu1YIY5tBSFWuBABp3FFpiF+uQZFTdgP+YP/rr09txGAOeRFnAyB5XBZGvp89AV5blY/gaE2ibjm9T/T09PZYdvnZURelPb99SmRlXHfF9uyTkr4/lTKsiY93y+iwrfHlLrfUeqO03fvfhj3eywDICOIfFzXxzwXhUjFPudkgNpgSw0Yr+rtzVys1CtXiPxuzvfPb7OszhbLkjR+VQKUBNGOYpzA+RjiOIhhyvNcAObnVZGvIn3b88fHt/N8DSgegIJkFMCHKDxHOGABHIrgOH9gIP2/PO9dbSLv7FJqzQpoCmgJAsHLWJYH7WFBY0gSFLcsBXD4AIZPFUHBt5AXccxUft51q2Xf3zyl1MfX7tnzAPkBdChcb3htYnj0HJGRhn1NKgItrO/u77+s2/f/qk+pM5og6Bnfd7O4mAYQspZlIWZaUkppgBAkMRAuAR2h6vZxs+eRPPH9GEhKrlss+P6N31Xqw5v27CluBYguEnEaluv44SGONTIbTtzfbkyqdevWJS8tFD7ea1l/3InXBTCcZnzlAQpFgEBzSAag0ABhTICAbJAY4s80IAEgAA5NvuuK5eKP59kEypzr3rdf5PfXjIzc3+ggiQBCoTlBwxBEexgifHFXV/eaePy2Qcu6pAm+A4BhERxNEHaCIwsA1AHCtG1LCqTBgVhIBAnBcQiACLAhDhSF4/gglxql5Dj5g77/poHR0a8BSjburDcInKDsOiavFQHkmLD1yAs14HhJd/dJK23726uVWtvseVUAIw6SXAgc0CJaexAYhuIERQxuuQFIGCTm9ahBCA5StRqAJIhdAMWm9TUGXwcg+VuAhD5Lw5lbEUCMsJxYsfY5Tunt7TjVsn66Gk21LTCpWiGkGhx4V609CBJoBq1BAIAsiABJAxgWwREmAoREDUJiMAAJtIdIpRIAJYg9ahQfZtcT1eqrV+/d+/VG1CSsiCicWBxga5Ulp52W6LeszR1KnRQDONBsG0NLlG6qRd8G+zfmCSBB34Ymag4NjjhykxK4k5REGxcphTauMKXRKEwy55iH+eNxFAj3Hq1iXbHYvz4yMPB8wMoFSBpKZtg0HoUTiAMAh27KPWtq6sZ2pS6Kw6yiuLKi8P3X0kkJ1QSQKEMAiAI4SFpTUHsYkDA22gT5tD/C31xz0rV5xfMoSwc68AEp5HHh+GfbHOfLWwcHz5Ndu6q4ypyIln8gn6Nw4nBAdwKuHhjYgA6/P0p4nkMxJzDq4IBossm2ThBsq0aKsQEIYwOKMFCMRjFaxcRGw8xrEHM/VJLlrLDt56+tVj9ILYLXaRi5aZgfeuJg4PBvMhR8lVXadT+LISPonIC816Sx9m3XN+s0v+Eh0oCpAWWBY06QHA4sYRMsMKuerHUCwPE13FaR9zzQ03Mq0uhAaQyQRADRIndC/LGHIXi9/f1vg7P9AvgdHiqHzatP6qyibWNIvzmA4ofAUgcOAKPNqbBWMWmjXcJx2BQz55lfoVNdKT9n2ym0mP0l3omP56st+xAB5MSpYo+OOXyOP8WwER8VQ5GvS6EGBIGAcwSDib3Fad51OCJgeM0Ax8QGNDw26XAc5IvBtvKaLet1dw8MnIFSGsJhjwByYgAEdpD4nRMTr0rY9hoF7UEx5qtpIIRiAgK9F1KPkXbNOZ4HPe1ggMQbDgUcAyDGuI538uCLxDOue3ntGctefpb9D3zawrK0GYkDjpl6IwwaDjnXxwYcGhC4vjimt6wJLU6M0WsRdHfPt0KhJJTC1qpDxSYf7tXXGZtgwLMwVhXcg3d8/T8ODqK9GNg0+ZdpHAFk6SuWMub29fW1Q5QvRsccjzkjUKMkDIo6IHDNgMGc4zG7uRlzHGIdEKYpl73lpmPQxObcYgChnMMARqEr0cc4r9VrHedMAvnf8K7MvlxDBJClr1ldB8rzzrNtu43fe5BFYBjhXxDXgMBzGhQ8BlVJEPQKyKHAk8wwEhOzx5xp03NuYnPd3McY5T0JJAGvXPTiW2iCvoSHnXjX4PTy/EvbNwpLywEtYL5tv5ivQYBAhFUYFFoz4Fo9Jihg+vBYA8MApAaSMgTchtCzX0SDwZhJ+gF4gvEtzLEBiAHQYqDU8hG0JL4bfJJzGWHULwpcviECyNLXrRYwmEWnsy/DCKABSB0UeE+tMQ4BDPoFMQg1WsB0THBw6AnnhGjnmwBhQD6tWcz58DmCxGiUxUDhfQh8F+SyMKwRxXircZYF85WXbYhMrKWv2kDAfH+AAkwBNKCgIJq0iXmuTshPDUKAkKg5NEHYMflJKkbgOUKXxEGIi8lcM7EBBwET0iRu7VkoRWH2obhKrfz04GALEIKj5WtmRRpkaQHCL7C/fv363MzMTB9FDWhRBANBYGKmFxCFFecq0AzoUIRHH2gQPYARgs3xWRYFHYHffozlwhpASFHgKfiLNQiv8TzJgMoc4x7MXdcA1O+ANGMAJpb02eC2vEMEkBOgfvP5PEYT6rlP2r7Hd34hIHDMc8xEQHAKLSuOwNCggBCzKckMXMTswPonnQChgCeRR1HoCY5DAcSAxwCkpkk4J4SOP7WG1lQojzFAwqKXfYgAcgJUMZxpzHp1dfsu7S2CYTFRCxiQMNaaA7EGCWJ0LqLpy6sDA6e0BLPJl0LughIABn0TAkn7JjoT5Bx56uZUTYtwGi5n4RIcBiAaJLin9m4RQMi/KBx7DsA5x+gRC5/3IGjzCckwSMo45tcsrDnoQBIYjA04wjYPNYcHQNB/YNNvFWmscqKJmkffh3vDADHAcAAOUpUgwf3ax0FWgoTvUgI29wZ9Nixh2YZIgyxt1fIrbO3atWu6p6dnN/pBOuArsJnXhgAKu6qpOWhaGdIAwbE2qSCsGiQ4JjCMn8FCSRogBAcJ4KjWwKFNtBpAtObBdYIEnZRa22iNQ4AQVDWiFikhD8Dhw/exkN47kUrl9XODx+GJyy9EAFn6OqWM4UOv9ppXoXDzK00KA8SAw2gONgtrATc3Qoh9AIGCrk0rCLTRHgSI0R4ECAGmQYZYA6sGkLBJRrOs3s+CfAQtyKdmy4uM3PTII2W867JePysCCCp7iQPlna7HwxD4V0BAddMQ26CKoMWaw2iMMEi0iYS8DBocAAAHM9Jm0wBBmtqAADHawwCE6CTIdBmIqT0IEoKjDi4c08SiBsE7+dREecv6GZ93LW5HxPdfliECyIlTrXfWXoUyqwO/2EaDGIGuf/VxTWeE0GqHGzE1DyXVkBZyCDNBxp53Aw6jicJlGYAQHAYgNNHMMBbjgxTwtDmUl/e8HyBa9iECyNJXsXbO4X/8sOq5U/gct+KVKOvaa4dABo454rBAExx1JNVAooUb57XvAUBwGDyBQZDVAYK0AQg//SRdFssgOHDMcliG1iCICRKaVdQeIHuf583BHvxvnJJroaSGmVimIQLIiVGxavfu3QfW5NonSs2pVi2hNfmnYBIkBhx1gV703lqoAQYt2LWYZhrX7CU4uNI7Kzvsf7BMgoNlEhw6INZl4cAAxLSqQaO5Zcuyp133jpvGx/cP4Vbcv2zNK/IjAoiWiqX7swGyvw0f+D47+6HT5hJrppXtPpyr2in0ilDyKLb0RQxAKNAMjGsirWNtVkG4te+BmGpJaw7EBIf2P5Cu958gTWAsABzBgXNhgLAcAoTmHoBqFeAozVjWzThk4O0RQDQroj/HggP2fwMc2Xj8jF4v/b4WlfZOn0spJ16QXyRLksbEWy34eDK1yAJhxrEWZsQmjwGJds6hNWgaESQUcPohWoMgNmAzsQYd8prAFLUHy+O9NYC4JawDPOa6D42tWPENf2wMT9CjYcxtyzKONMjSVSsXiJMtkPseN3PLgNWkcpJwU1ZCvTbfJP8YG5f9dkWSNZDQXMqH3tWAw8QUZg0MxMxLcGjnHGmaVcbEIigMEXALQId7THmmLILDmHnUHhOW9eEt27dXrg3wxkct6xABZImqdwPkFOBwelTTlQNW7vxWSblZK2XH7aS0WWm5vLxSPpt+TAu9MasorAYkRpDDwDBCrQGCvAQI7zHgOJQGYdmGkKyba6YsgsP4HiOu+8DuU075ytDoqBpuAO1BfkQAIReOf1Dfxwc/JbK6w8pc3yEZ7PGRVBmVxCLUSXGw5M7pfpNsdFbKF+K7pAkrZKF3Qr8lBX4OpM0pxBRkAxIDDOYh0bwiSFjJJl6sPYwGQRYd+BRTHstAZ6VP/2MCpw8odcW2bducrkAJ6fzL/U8EkCWoYWgPtv44q1TrDf2Sa26G9siolJUCOBLQIEmVkDmsiPUyp1d+ac3Kj2MTkgNIsC+BflsCgZrEfOUZG3AwNgBhbIDBmOBghTN+knmFcyydZMqrAcTJKxWH7/EX3xkbu3NjoPmYpSFCBJDjX8261arNTr+2R3KvbZW0k7ZSsZRFYKSwzmgC+3oE1VL0PfntYr88msnLtKpCuIP1sPjKlFA67kaYw8BI4jxNI6NBDEgMQAgOpmlaMW0CwbFIe1RLATi+/oOxsY82GjjIlzB/DJ+i+NhxQDvmKL65Q3I3dEIvZLAXFACCDW9SWnvY0B5Ynh1NrYq91ZKq2nLZXBc6/QLtEX41nqH5MwuaqcVzoZhpc8w8Jm1icx/vJU2DpkCTIJhUVZhU8V2+/91feN4bhiArWwL84GrjhFjj/NSl/6X4AuuV23vtlut7pHV1zmgPaA6aVjFoDwVw0N/QAwWxduGU70h/MSVnxJvkZ9k5N4VhjQAGP/71QO1BTUKtwQGO1CAkag5DrGhqDUP8MrIQEoFGYhkgDkZ0igDHhOd97yGlXi179xaHg48pszVUiABy/Kpbr9yei6XP7/Bb3t4qWTdjpe00wJEEOOLwPyCTkFgMEcGA9yoW/ym7GGKO9By8ijUzKX9nqmSXFCZFBV0UlPF6oOTSZyBYqFVoXhliJRsiQHjjgptxXDPVXIDDJjhmXPffdom8VbCZZy07ra+GCxFAjk+Vmz4Pu81rvqXLasXGm2n0ecDz0H5HCkviJoCNYHITAVIBOAqI86Bp9GpMeWXLLjg/kJz9Egi3wpASyjQVwAJZXwwUAoJahBV9KO1BqQe5AJaFIfF20fNK6EN5/9jY2GdwiYHlNyQ4zI9nHIVjywFtWnWqFe/ttlrPyEkGjnnaTqo0TCs45tQeVhxmjo1h6TCT4HsUDUAwxW8Wo6lG7OK3H5kefylWdLgMILgP80cwKRCICqwjKo66EBM1DDxBs6sAmgFNgyZBE7jnIC4fgMJBLJMcX4WyZn3/62Wlfm1sZITgYNksql4u0g0XDC8b7ocfxx+sv8CJRGL9arf7/tXSE2+3mqhDrBY7Lbl4StJ2XM8XFygF16lIySlJ0SnKrFfwp/2Cv1MOVh9x9535uMz+AhVGJWH39/e/EfHVoLMUcOKhxYtjqBCoWXSCB+FQO8khItwUB//RKhbc9B3MI/kktMa3avmpbFhOw4fIxDrGIrARX2G0/kin0/rZLqsjmUGnYMpCd6AxrdBqZaFZl2OfONlW+x2Y7pr3HMn7jjsp5dh+vzC0U2Z3XBhYSpRzd2Rk5EuI/7W/o+ut00lZm/Vjr4dKWQOpp3A/ZfBcDKsS2Y6ZhN/Gog+bnxgbu7d2A8HMEIEj4IM2TWvJKDoGHNCOeatqfnO71XZpTrIwrTKxwLRKo9UqcMyxHgm8DAeLJGDWHswr+h4FgGNWqrFRf/aB+72xTwJouiy8o78BQNkG86hPZV93/kTLBy+Qvs98NbfznG258Y6VTvsFrpJTkK+99nuAhUCjQGM8gfQOtJDdu3fv3sdr1xkxD8ERAYPcCAUyJgrHhgPma9y+2l71wKD0d66QZoFppZpjGWmKpSUTS2jTyoJ55LpVKVeLMK0K0B4Ff8af80Zkwtpp7b/gEWffnXhFY/awmZfTctvPtfse+A2rv+d5yQ5pTaTv7Xey775w9os/fgY/hxYEfYyG9jOeil+mEp8qT3Tt2XGAvPW67e6PdEhnd1qyfgqjrRJ2Br3lgfaw4JhzybcqxNM06Rbhh9C0mpGKvVdmbyU4NspGAw6hyQZw+M+zu64/XXp6VthNJSuWqKIH/uwnZOZHH4pdcMtrkutXnyanJc6Rc0wDFoFgiGWZel/g3D+7n7m87zKMWt6/8vj/Ogqhk7GbL22VFW9vkhz9DjuJUboJtFzF4H8o+B4+3AW2WlXQ58G1dIv0O7yqNwdwjMv0zhF37prAtNpS/8KzNxsaxOqT3BfbVPYHvYn2VCaWic9Uq9X9laIz41ffPuOUzt4u2ytrZA3vIwjCRDOqXh7SUXgKDkQAeQrmPMtLxmxNrfBbP7tCOvyUZK2klYXmyMDvYJ8H+rmhPVwMJwm0RwCOAgBSAGQmsGbIAZm9ekImZmrvAEzogP6Ujepa3P0999Gtn6r+9wZRsff4vj3pOH58wi3FHvInvrHV3fXvQzKEpuUtkU8R8O1Z/zWV+awLiG58EgdoymBf8c6PDMjgn8O8cpqt5liz3QS/IyPZGHvOY3r7Wh/dc1WnLGX4HUU3L3PenDstM/ZOf+TfH/Aeed0G2QBnfBu//k8KV8B8ukXu4TX/r5s3rpsozt70c29sw/+4u059WKZ24Tw/fpGmeBLnntmJCCDPjF+/KrcWyng8e/qAt/LuHhmIt6DHI6dyVi6WBTjS6PNI6PVxuVSo45axRQEd8zxarub8WX/WG5Xx/E5355njEHL2VOCBFHJdbkc8c/bJ3qrXtrjOdbfLI+UhAKhP1ltXyi0cZSIXxFeefkd190NMR+HocICMj8LR4QCEeSMFWrW5bX/XIl2ppDRJwspiEi1MK/geNnrMBeOtPG1a1cZaQYsUvSoc84o3LXl7wp+9Zi/8jwvlQvoxWgPAXNJviHJvWGP1/MVZqbN+9MmWt14yDO1CcNCJp19SAwffIQpHiQMRM48SI1GMNq2aVfuVPdZJn2+XHrfZarWb7BxMq0B7JNFjzvnh7BCsore8QtMq0B7o85iyd/t77nzQe+gCCDw6F7Vj7gdOuri9quOqU61VN56s+surUp3JlkRGypXi52dmZ4aHZcs4IKSGg98SmVVHr07rzX1HsciGLAqaWA8e7G2yOq5vkg4vKTkVx5DEOLUHmnUVtQcccwfagwMRy2i1KlF7+BW/gPG3B2TSnbKnrwD3PICDTNSO+WkyxNjulo6LBq0BaY23JlUsWcUoX6/q+u8Ys6d/dn7s5F8DOHyjaXhzFI4OByIT66jwcQP4aPnt9sCnW6W3OSHNflw1WQRIDC1XNpp1aVqx1cpBI63ebwMtVjStCj6bdQv2lEx/ck91zwN0zPFKpvXJGpZhPQjlRu9zr0uo5HuyieaprN0Uz1er7lh1qrJHJisHnMIofwbzHpWfExVS50AEkDornnUCvsI2J223vLZJujelZIWTsJrsOJp1Y9AeNpp1LfZ5wAJz0f9N7cE+j5r28EpSsvf7B3f80t0xxKZZtFoZcPCFKPD+JtmkLpKLnJsrW26MxRPnzVUr/69YdeP7vNnEiEy+9REZ2cPm31p+3heFo8SByAc5MkaCfxTMLU2d9in3dcm6wZy0+1nVqppiGNQO3wOzzdE5iMGz0AMupiNV4XeU4XcU4XVgQIk7IaPoMd/18t3uTo6kpWNOgOh6GZTB5C7ZxXN5apb1aLG6pdZi9Rb1mj943Bo55/vu3X9IYEF7RL4HGHW0QwSQI+IozaFtTrPdd2OnrL+qWXodDGSPZe0WgKNJMhjKnoBpxd2cPPgbDpp06ZiXCA5vGj0fE/aYPH7bo+4DvxPu82CrFDv5+u2BG9bJusvWxta88x/K/7CVr3qFXBG/WW520ARszCnWoUkf0a+Jbn4yByKAPJknT/eM/tqnY1j0zV/3gzZZgwbddgxFbLWy0B5prT3Q58E5TejzcNnnUYX2cOcwYncawxGnZL/s2n/AffTMffLO/fQgEKgFdLnxePyc53nP//Hz1Knx3mS3NCdyn/eU86GPHvjoGDMSRBhv5Ueag9w4diHyQZ4db/FhQQMsXO+0131LVvpjcexaELdyVgyOuQ3HPGi1SsDv4GBE9nlgrgec8pJXxkY0ZW8WC/nM+BNX75N9ezdiagbKMn0e1AbJld7aL/Zbq+Npu7kUi6e8hJV+R7FYvudN6vI3QdukauCINAdr4RiGSIM8K+ZydO0WLPbW/aft1qkfy8kqJ2O11Uwrao8MptJifSvd50HTCn0e1TxAMkvt4RbkIBzzR7/1mHf3ywE0XVboNVgnsXPVhqtOip38Z/3pgRXxmHKrVQxGKc0knnAfl52xR0/aUd6xE/n4gYt8DzDhWIVIgzxzzoJnX4EjnVydsXo/mJRuL2bl7JhqQnNuVhS0h4U+Dx/L97DVqopWqwr8j7JXwQ6xZSzjWcTc8P2l6djo/wkevSWsBRRbo+BfVO/2tn2iOdXxa5aKf93yknax6iQOYH7hmIx/mOCgiYX7I3A88/p7RndEAHn67MKXfQj8upmCiW7AgZszMpiNSasfs5phWgEgNgCi+zwS9eEkFT2NFpOhAJCSX4L2mLFn/f3XHayM7pCgz8MIuTaz6Jy/Xl5vnwNn/O9nPvbIjZN/8ep8pfDGaWdu35g//ujP3Hs/wlarWk/703/7KOez4kBkYj0l2wgImjGn4yu/qd4/ERd5c5u64EtNss7JqM5Y2m5Dc26zdsyT2rTCLegIpGlVRZNuxZmBaYWFe2RCHfR/8dBO70fnwrTCAMNgOEnwDPHa7ZWv6fJXJR/27tjM12KLVa/0YouzYe+F8qK+SrzQ+VD1vvtNfuaJwrHlAHtto/AkDhAYp+PjoUGhv/CDg7+XqlTOWI1ZgCeXZ++5Ie40e8prhjkUdAiilxvDSmJYvsdGoxVXRvThnMO8AlAqGE5C13xW9nkzauJtMIy4ths1Ec0rPMfCM/xcVtpv6VHP6zoj9ZsXtmbi191y4ANj1BYbZAjD3odHsTIce8z5UTNaB8koHEsOkNlRqHOAU703ARzBRKP+/j8aEOm/zPczLxPJniOSWWlZK9Drl8ZyOciKFiqF+eRxgCGFoVhZTNzLwpRKoofDqsxItToN/wPaQ2bdEvo8DsiDt467d2O8VdB/EjxWO+nSaa+5eaWc9Zae2FqnN70qEVdqrFCeu/oL+T/9MvPVOgOZjMBBLhynEAGkzuj51qS+vne8UKTnKmw4+1tKdbVirWkAAo63n0DuJAQ0Br6lQDC2uNA0hpFwZRKtDgCYGFa3TWLqU7I6Lnb5Cc+r7lAzzk9+vtv/ydnwYTDJ6Uo90QmFAYziYZZI34B9/vbV6tyWXLLVbU40e7ZnxfOlaRmp/Pz2R7yfXD4qOw7CeUd2/RjGUTgOHIhMLG2yDEHyht3e3o2rLGvldVjF6k1K9cY8iK7nZeB7ZFAVaeRJYuWcFISaQOFSoWQfidaSgvgGIMH+m1JQnZJPnIRbz/Zj3gGsj3jJY11A2r59r9mrVU+wyzMwNaRK8qHRPvfUc2OJpo/kEt2bYLfZM5WDpRnnYHJa9jVj/fbCtXItSo3AAR4c10CmN3LQX3AyoK/v8ndAa3xEqYE2z8viTDOAkcX1LEBBgGDQIbYpAEgQk6A9aqveEija5NJfeL1aIUSZKx1yyxsHVPWxsidW9xkfFzlw1ejoi/8taBH7EIwz7ktL4Q/GUv1ux9+/sVopfyxfmlo14t6Xf8y959RJeXQ3HlZ/V6SjcJw40MgA0QLX29ubUeriv7aswbf4fhuEugXAaIZKaEIVZAEE9GtoCgASgIMAoRaJI2Yx8wFjElFGmLjvOIDiFbEsu4fMRRyP3Dg6eu57cFcaVOTd58jN8TXyHc4FcV8jQ60H7L037fN3bP+lt/WvzNgs5ovC8eVAowJEg2NgYCDteS/6ulJrL/a8FjS7tsBeagFPchB8ahEDEE54SuMcl+sJNIhCixU+/ZrCVRYGBxqxOAxLE897Hhb48Uu+QlOX4+/+D3/v7+c6/NP3ZpKnvfeuwjXQLkPqHOmz75Er8S71wDrC3VFYCg40IkDwm4es3t5bUkr9+tcs66RLPG8FBHIFbCaY+xocTRogBAm1B8EREAFCQusVIEaAMA4HAxADDAzB0gAxMbWJ5+WxXhz6yPO3S//sA9KSGhzLxpvf/5WJTV8IyhpSQ0ig/4PAiMARZvBxTjcaQPB79fwNr6/vom8qdcrLPK+9Bo5WsL4ZQk9wkGhasWecAOHEJyy6oKfO2ogDYBAcBiim3gJNMQ8KAiNMAXDgn3gF17EKVlfxTn+dP2srq0VmCnu+MVMaed9PxEYv+7UARn1Iuyk+io8zB9gE00BBg8Pt6zsPzngvwJEBONLQHHS+2SoVON50upnmjk82FlpgHIslkQ7AQYAYMgChNmEgAAgSAwoHDbokY44xTX/bwwJANoahHEifKW3ln/np0gEHswwvm1aPrxLvtnORiTkj7UF2LWFoJIBArLegKXfty9HZ9+cABwQwCUQEzna4RcqAg35GAA6MzMVib9i+fAGFNYmpQ6NBCBCCIaxlDIh4Hqshotc9jUUcXHk8fpLVPzeiZr1fypS16zqUVcHuBHjfKCw1BxZZ0Ev9Osfs+fi+D/kbNnDVnZYb0QGIB7Gjz4CDGoN7ypiYafoZWAGRqyByuR70DcYBJ1ICtxlKJgEzUApKKEzmPGNzHwG2QPNQQ6HZuKK63fHsGnvOe/xfxt0fbjnEEPhjxpio4KfmQIMAZAO+xsPeL36x+l1KtZ6KXnF8w+M4RwXKD7UhdPZpoCiAg+YUiQCx65rDCHsYKGEwLE4bIDG/0UAESWCaoWsRrWLKR997+nQpdV95e1Bdf1gz2J668qKrx54DjQAQCNs2d3CwpRX+xp+idxx2fRy/24CDoCBBWEMx9gCEEAdk2+jaJpxCJlYYIAYEh4rDgAprkHmQaM2lBENZLKvzmtNOG4J+uhAGWhROBA40AkCoHvxKJfl6TIXthXEE4bPxu/nTDfGDTYDQmbYADKYDwAQgCQASBonRCAYAYXAY8JhrJq8ByEK/hM+JKzT/urbd8fyJifNeFbRebSWCo7DEHGgEgKBdiUhI/g5MK2gPSL8GBkERDgFIeCYACnLhLoLFCLSJCZQwWIzghzXModLmPsamLD4rMPHieLc0hqRk3sgzIvtxHIWl5sByBwh/n9fbu2IAInm+j7l/iCGe82AwFWChyyGg4IwBCWMjzE8Vh4V/cdrctxgY88/Q2sv2PKzmI4mL+/puag/monBMfRSWkgONABAIfuwS206gPUnBvDIyF3ygCYqgu2E+tjB/iQRs1SioonmBntcyhzpHQBwKWObc4nuC0qnZPM+2mzFYsu+84NyW5V4/wc88gf82SAWoF9W0BlAwDwQDgIVgMKAIx7znyAOBYYIBCY/nzys8lB2SsRcH+TpDd5g7o/h4cmC5A8S0BqFpl2zV0wAZg8IAYLYwBUPVg3NoEa5N4mMZT0WmF93EJi8K0CF4B3MUjqmx9DthhVKO/vVrGuRCvmQUlpADy7mlhF9fn8PZfd+DDxIIYCDsYTBgKDrmbliWOYf5G9iznOBgTMKUQBo/ehhJMJYqGEpiTCnWX1j4w3lMOgwa5l8YiAM+n+/I51pttes8EYUl5MBy1yDs5EMft9VaAwZAY4AQjucBwSFQATCqoZgA8TRAzBgrE3PYiCEzvCR8bPIZoBitYuLgvYLni1TR3IvhYeL2rV//vhzeGwCJHPUlxIfuLVvK5x/zZ6MfAz0MRIUBBIVxnggGyyIoqrWY1wgOUgWgwBpXWIjB87C/h8uh7gEgjN9AQafwG21ijnluMVAMSBgbgBgtxTh4Rw2QTD7voO89CkvNgeVsYoV4S2uLZgzGAIbAEaQpkKQAJAEwAoAEwCgDGBy0SGKvN2cRIjuCEXKCg2QCAWA0hwEJY3POAIVzQ4J3Cp4fvId5H1NaFC8lB5Y9QGAaUZzR90FhNMJnYgImnI5D6APyPA5cjEE7BIR+bhwTBeyzCNhmtMWhAGJAQmBU8QgDFIIkAAgbAip1WgSUQjYb44tFYYk5sOwB4rpYkkqsKXzvWwSr40IQARgCYyE4giV9aGZRaLkIgwGHjS8/wUGc1VQHppITOOwtp8Ab84pZjFYxAOH1MDgCLUJwlFFGGbF5F/0+WMGhauPa6I4dn5hFaShRN28tsZg07uOXM0AABrHGxsYKfX39e9CMOghhxLkyTpOMYFJIaTbxeB4YTFNjuG6gNQJwECBshqXDngIlAI5gKIoxu5AB1wMtEYAhANE8SDjlluAgEEsgPj8g0QsulvHcwgTLEb3aiW7aCg6jv8edA8sZIGQmTCs6Hf7DEPYL4F9AdKFQ6iAhUDgnhOCgkDI7WWIjjVl/XIQE47bmtUcAjgAgFHQHAOFMRGN+ITuCAci8FuH2a9xEh61hhvi8ACTBO/G9ShiLVURcvIvlIPAF6DxFYYk4sNwBYtj6IyTeBlmDCqAAamFETHCwsSjQHBT0AFMYfQ7rJgAITiGwA4+yasBh20H/SOC0UwMFwKKm4b0BSLDqFVACM08TweFiK4TA+WdDgNEefB++VwGAmEKeuR/iAGE7HxqFJeTAcgeI/vpCuL8LIS1DE2A8FrUGhTEFMgAxIJkHiDGpqAWCQFk1GsRoD/osATgWgysMkKAPJbiHWodAMQ76/PsUqN2U6+6fUOqJmgbZjKfTrIvCUnGgEQCi4IfsQY/6nQDIS9HfAKHLAwnAitYeATiM5jAawACEFcPV2o32CJbtoWlFU6lSB0hgmgUmmdEgAUjmtY6HRRrY3xEAxICkrj1cpYq26858b3T0iwdFNuMddfc+XyEKS8SB5Q4QspVSy164LwEgGyC0kHZqkTzIaJB5zUHfgz4HhTsIATioBTiokQKOliYt5E82r+YBwnupfWia0Swzphnvn+9r4XsQIAXQLFTFFMZiHfgyDhAeilRHwIgl/dsIAEFDq1iJROIrlUplGM56L9qVqEUgzQQIwREIdqA1Aic7aLKd1xzBEHgX4Ah63gkSaqDA9wgWewice8p14IcgsQAg1D4BQAiSoGEgAEgezbsFaI/RB9vafvj18XHdvEtQR2GJOdAIAKGUx3bt2jXV19d3A5plb0DfIUBTBSpmcWleewQACT7cgVMeAITmFUAFomAHGiQAR2CezTvoCzVIoIWoQYwWMWUQXCRjXmGlResgxmEduH779oeBnGtZLxFAwISlDqzRRgiUTAig/3k46ztgQlEAcY7OOkEyB6LJZYgmDykP4Q7Svp+HyVTAMWOeD2KmPW+uRgvTJo/vzyE/n8HYPIPPnYWumRXHGrFaymk5ef8lr2BDwEa95RsuR2HJOdBIdi5VBRaO670UAxi/CbDwCw2gkAU5UDOoCZQFZUBpUCpEZkV3mmU0rag9SMRaWAvxm2O0EJK6n48ahBgl8bGBeWVBg7gygZ1HynJ2eZNbminY+927N/3U/bMtG+Gkbwnti4iborAEHGDtNkqgdMbQovUtmFofBUiwuiLaWwXd6FqL0JwyJpURZhNTqNnqxJavoGmXve8BOA4PEGTAPaZMmmlBOQFA0IuuNdZeOcW92G9NtHjVdM6uzPV98DTZ+NXT5CFmjsISc6CRNAhZzd/LTzwWr+7D4tUK6/MSJFQFvETtkQNRgxgtQk1CYBgKNMi89lgMEJYTZqsBCMFmtAd9jxK2UT9gNedT/sn+xXZP7hSsTSTfODg1/r4zK3ftGI4WrwaPlj40kgYht420cmrr/wZ9DSDB9gdGk8whC4U4TMQPTS2aRaRnAhA+jmTAwXLZtFvwlJqzK8VHxZtZLTPJXWOZcub9m/f/7y/gopyJfUKGgLJhfW80WJE8WaoQ/tQt1TssxXO1Fgk20PGwgY7CBjpak/CDAZ4wojYhGV8krEEIEuahmWX8D8bz/gcOEAweCRACrQI7ixvozGHE7r7/2Dd2d67Fet7eLvfs926Xfx7H+sHqCmygc8uCDXSiEb3k5FKFRgUI+a1BEmzBprAFm4Ut2HSTLKSZoxZ5maAwpha1CI8JjsMBhOwMs9Q0DxMgFfR1FFHoFPySSWzB9uB7cJLoK4IWbMH2x4NbW6fGx2/aV310+395H8AWbJHDTh4tRQjX5lI8f6mfqUHCl4BP8g5E2DdEYV0qCrZudsJ1LMyrAUGAhEGy2ElnUfMACRhLDYLpupL3LVXE6Me5cZFpbOK5+9+oLbigNp6DrNeCmMaIyq5/eWNrovtjGA2zas/sQ/kR58FTfyCf2z2fn7micLw4ENTj8Xraifkc8oCEFRh7V0GTXIc0toFW2AZam0j8/CNwYbcYhNpoEKNFaGrRvMKQd2gd6CD8w1AU9KN4Vglnqm7Ka1aW6/5XQe5/6759sje0vYF+NuYoemvkpes6Yms/sib7wk092VUS92Ol4kwhOV785Y+L3s6XnSO9pWhLNtbD8Q2soCgEHKCUazBAm7wQ6atAvwWgtCLWzbUcT4WAPxyIQs0SaA3qAMz4AKFDHNZZHH0qOVklzc46SZcGvWSpS5WqMz8v+aWzq7LDuUduhreunW8WwEL7Mnbb9pPk/Jae2MnuquzJXm9mMB6rxGV08rHbx6v3Xf7/5OaDABKfT9RG4ThxQHP8OD3rufAY8oNCq4HS398/gPRl8E1+E/G5EOrVhAa0DADDAYhIY7UT22uWlN8DWKyTNvtUWWGfIk1+H/xyW8qlWclXDrqOO2NP+A/deof74Ss2yFBsmwyzSYuBwBRMRbm5XQbfcpL1686a1PMTfdnBsZxqvfpDY6/5Mq8PaZMMhljNFOO5KBx7DkQAOTSPA9VQAwqzDA4OphzHWQ+wrJqdnblVnExn2huUjLdaNcmg5OwBaU70SC7Zil0PMxLnLEPoBrdSlUqp4Jcqc9gpfa81Yz12wR3OjXeG9j5nHVAr5OJ24pFVcmbXGvvszzXF2q/7j8JHx2rAUACGAZTJz9eKwjHmAA3oKDyZA9qWwmkCheRjsGMJ8f01utpW5S9Z0urErRlVsfZLCZMVY25FYk4RnkqbxGM5bNOTwMIO2IQnjlXlvZgVq7Yqx2v5u0HZAG2k232NsPMZs5h7eMW4vyP5aOWuzexyuUJuRkvAqEtwvEHe1efFrc7N1b++n6CJNAk4dhxCBJCnZjKBYsBCYa4BxrrNF+d3y9bYpQnJunE/aSewiIPjpDBGOCkO9h70uIYWcsew54HChlbKjamSp5w26T19rZzzni3epz4aMrX4DPjx7lfzMoPBihvtx+QxZfpDbui77Q3F2cJn9hQend0gG56PvOyxNOBCMgrHigNkchSeGQcIEgh0crWyqw9mpDvdLGutnPRZTapbcrFOmFodoBWSjTdLErN8ObjerZSlVJ7DTldz3ow3Xt1nP3rWjyr/sGORNmDZqBPegVG9yT9Y15dZdeNgds2rrKorM9OTsrP8iw//o/+5D25G38imaDDjM6u5Z5GbFRKFZ8YBfu3hWJd3QoavK8tBVZJ9bhmjciv+FLoDQc60VJ0ZLr6ArGVoEVcScUvSsbiVxKagrVZbaoXT/lk+drucHv5IsWyQH1+h2t834jzx05nyxKsOlkdciVUqrZmsdMfaP/B78nurCY6af8JionCMOBAB5NkxVoMEHYqfcP3qA0XZGwNA3Io/CYBMS8WdkjIAUgVAOIfEgkUUt31JwdQCSOyEUk6H1XXpK9QfvplD2tlTHnoNAkbNWjNv2StPrBiv7CpNlMbh8nsJbAE01pZofjOuj9c0T9TkG2LcsUhGAHl2XDWCiUGH9hUVmXMK6P8ryaRfhgYpa5AAKDWQcIE4G1N1E/D4Uth9GiBRaZX0WqzcDRfLG7oxtB3lsWddBwKk7Kjq5WOyp7rfH0tNVw6qA8Xxzxfjc+dck7/mX/5Z/rl0umy3AJKw9qndHkVHkwORk/7sucm+khiafu+0bOuvi3LwqqQ0OUlJxcpYcbHsxiHlMTjtNpx2JUloEO65rmBqKddGz7pyO/yO7qo9c8OwO/w7G3TfiG4QYLk22rjuKdqFm/bIE5d1Oe3v/NDEF7fCipOb5Yo4WrewmWGw4gmQpafPP/ufEd35VByIvkBPxZ1ffU2bQ8jWpGzrvrS0DrbLSVgEuFu1qE5piXVIa7xNWhKtkkvkMHMwhU7FoG+kXM5LqTLrzrmT9n7Z/fL/dL/wrUV9I3w6R0fS/Mqj9Sq2XtZbt8gtHBYsn8382R+MlsfO+Zj7T39YM7do9kXhKHMg0iBHxlBjak1DL1xdltn/Oyf7XWgRlcTswxS1CKa/V2pahJMQFZqAoUjw/Yfcu3FLeVmYZLnPvkJeceZpchqFn6Az5bLvRTf7bpEtzjbZJjd0vn9dd2LFp+Jl69XKqcgfua/7T/SJhMHFW6JwlDhg7N6jVFxDFqNNLfRh/DsWGN2cl4OxAhz2sj+NzsMZzBsEObNw3OfEcTmyHZ2JUCNJOuzxmEraltul2te32d3D7PyDqbXYYbcADg9aIvbK2Iv/ZLwwfpfjlV/tJ6rVllSq0mvn/uFt8ooBNPtyaGVkERxlEYwAcnQYSvPG8t30VVUpzczKAasgU37JBzgAknINJLpVCw4719bSDjt62emwpyzbXSHN790Y33jGNvSa09SqvRY1iTbjPq4+9n8n/KlP7avua32i+ES1apXsjmxToifekmiTFAZ+WXJt5LQfndoMlRIBJMSMI0gSIBDqwhiMo2uKMq3m5KCHGPMHAy1ShhYpO3No+sXSQegbsbFle5J9IxiGkkJve5uVs5vd9C0oR20MXsRoA8bY5KS89RHZLWPuvvJUeTJerMyqiip/vikXf8H18p8/jcARMO1o/zWVcLTLbcTy9JcePxxTaq07Ulb2/C4ZcNuly0anH0b4rpC2+ApZkWiB056TDIakxDBO3uNgxnJBiuVZt+gW7BF/9zv+yfvKzUNwyodlmxmgyA8ZmpTV1vXWwIXnWGvuPTnZf81Q/rbvktGboXE2wgzDC+gBxox5PgpHzoHIST9yHpoSKJQkD825byt7hbtn5GA8jYHwKYwyTMNZL2Gvw5KDfhBbSQrTQRTW1uIuVQpbTAvGaik/5rZ5Tde/Rl7+zWvl9l0wm9Sw7lkPHuHZ3tV7nP2vfdh/4jqsK1EmiPrQsrWp1rL1gfjLT7eq33zIvFAUHzkHIhPryHkYLoGmVqxarT6E2SKfmpNpe1YmYVTNStEDubMASc3UggHmw2G36bBjsG8GplYSwOmyWpr7Ve4G9HMAbBtM/QRNuFW5d9YvfgDPKF8h58ShYdwrAY6tA3+y7put7/yvQT9719XyG4N8oSHgjnEUjowDEROPjH+Hult39Pmuf13Vr+6Ylkk7L9Newa+BBL5ICb5I2cmL48FhR48gh6Ek4bBnYjE7qcTpsVpe+wf2K19LE4vmU+ghsJ7EpvDfIvdguyxfPpH77T8puJWftmQyrxhMtqQHrKa/CfIjVxSOmAMRQI6YhU8qwK+dKSlL/Z+C5LGD6CRW5J2VAtbnLWId32INJHTYXawWp+Cw68GMHIZi26rZSvidkrnhzfJrzVuCwggMBpbtXYs4aycverlav22umv/UwcrkilmrUE2lY84ZiY7LPmG/FOAa9haBi/dH4RlyIPJBniHDnmZ20zfyLcu2bp2RmbdnJe2m/bidxu64aXQSpuGLlDAEJY21IHw7iclVFjoQ9eQqbI6onF63eXXBbr/+s+6Wd1DQN8kWlsmAnIISnMt3WpMv6aymSkkM9eqIJ+PtybTMVNStTWLdu9k9LTEpj/nQIwZcciFu3g9wwUnxh0O+DQuNwqE5UGfeoS9HZ4+AA0Y7tytbPYBlHDoHpFs6pU11qBbpiLVIO+aLtCebpQVTdDN02NHVx3kj5VLeL1QK3qQ3Ze2y9l/weWdbeIouy6UmaU/YsQfOk96eF9v9cmqm+96zM6vf/YLRz/z46bwzK/578JcuBFCQDnycp3Njg+WJNMixq3AKHf2H/RDnq+es/JemZNpJS0KlsQd7BjoggyEoJYdaBCvIxdIwtRRatTBoJY4VU9CyZVfTas7TfSMv3IiCYG5RrlkucsmBiu9ctcOf+GSf2/SZ29MP/12rmu64qqfv91JKTmkCgOLIiMWJuNeD3r9XWdYTuHlH3nHufd/evY9fFKyxSrSht17UJmgm3BKFEAfI8CgcWw4QJK6lrG9mrNSlA9Ll9qBXpNNqkc5Ys3RAi3SgX2RFIosRjwmshwVprlaw0ENRCpW8g76R2CPe/vd/wvv+x9msu122+RBmI8ix3+hY8dYz0sm1Ta56fYul1rSiRgGO+oKpZvUuvgRVD2/Muy73m97uiHx7zvc3v3Vs7F6cJlDUtYiHI40CLgQhAojhxLGLKZeY7ZRY77jO/e3SHF8pnVa3tFidqlmDpBMAaU80wdRKw9SKYfEHT5xKSUpYDaVYzfvj3lT1UfvgmZ+r3P0LVBjNK/uV/f1vbBe5ulOsszpxMoMuwhQoAwykkIfDgEk1LaJVGV8EWTlV3raD9YtkAsu0YIL7d+Y875MbsTUEsvABXPTLgJCnGjZEJtaxr3ptamF/xB0woYamrfzHmiTpZGBMZWD4ZGBKZWBmZeB609RKxZICn0U3+9ZMLa/PzyaLXuEmCO2lv97T/8peWz7Sa1lntQMQWd/zm+G6ABiKhMV+4fIH4DAAIUioQYwWQTl6EWJgw0OeGMDysjlsBfHd/v6vj1nWn1t79jxAbYJb0NevAYlkYwb8/igcBw6QzxQ4fL7V/+Ss1BmD0uH2SovdrXIYjJKTznhOOpJNGI6ShqkVh6mFdRorFTjsBZmr5P1Zp2B9O/f4D/K56ZecBB8m53tuC4qDOYV1uRYusQ0NUl9mm+CIQ1vwS2gAol8E5TMQvVAVbhVl2UDwfs8rTVrW+y8dGfkMX/ovUf5wA5tckQahlBz7EEgjBBFCeEXeL99xEP0iWcAAGsTKYMZ5xlWSxTzBErUIOtETOJcIHHYuGWRV8M1PpssvacU3vcn3sJaj2ASGAQe0R93v0BoEoKD/oX0QpJ+kRXCOI7f4YgCJ7SBdwUbuXXh8h2V9+of9/Wd+2/PePTw2VkAeNi0TSw0XIoAcvyqnTW/XpujeOiH5K3OScLLo/cjCzMoCIHmYWlmIYoUgsROSgCdgYcRvyknIHZlxLPtbcFf4cWgMX2VRmAEF91Cg1tAEwV8AkBo4YojDWoTaAbsIaYDA1NIOh4P3K9Mn8TznFKXeAs0z2Dww8CqYXMWhBtUk5FMUjh8HyG9aOFnYM/evkPTgGlmBVXxzqgdi3x3LSle8STrRotWWSGGKri0t6Fjc7u6Rr8XukV4YSS24GdpDa44F2gMA0EAhQEjIwzgBoomlNQjiBRUOYOjxv7WYQKEmKePeou9XYQ7Gd/r+dx92nNdcvndvEacbzifR5ih+eBSODwcgkTrMKLGuxlQqLFqaxxooRcEgRJnDjMM81kcpgIqYAaLwXZ9WefkOGnebIdppkPEvCABjQhEEOs0YpAHChSJImGvC+SYxDBu2EOvhwxxCHE7Xjnk9jnsyKCOHbXwVdt1aa1mXnKzUvwJYNLEaTl4a7gcH8rmkf2lqxfQUXd/99wNSiE1IwZ0GQGYxwncO47PyoCLW+PXgOn/DfxB74eYBEBtaAYMacbMBBrWC1hCItZYgQCDgmiDsBEcCVAdGHHcQDKRw2pxDXgKHHZYplJUFSODuVNfa9qu+39v75wCJuxk58LiGCQs0bsP86qX/ofwwsa9isApTq1PS2bXSogYgkr1orO2KZWRNvEUeTU/KXfZOWQnvoRngoHmVA9G0og8CAdZ+CL/4aQh1uhanKOAgahCYSfjug3CtHuNeHWhacTctxi5wy/SiuAS7q4Dm4Ans+fCYyItePjLyE4JkE3KaYpZzHGmQpaldmis2bP2dkPtrJqWMpX8KWHKuJNOYs17BKvE75ID8yN4trfhgx5GJWoNOdp0g8DrNGACgE07ziGQ0iGU0A7VFAiUwNmkeGzLnGPOemiYhoGCuWTTZugCvDt//BF6hoT6qEUBQ40sU+AW2sXzprVXfu3OfFO2DnE7ll2U/Bsej1QpmU2BS1UGBG5imjaMJgssecZIGCYEC4Y5DyFUYHGFgJKG3DkUmzyKQUANBM2FXXt/tV+olP+zr27QJ2mNr8Cp4k+UdIoAsXf3CrtHBg1BfMYvpU3ulYE1iLZT7muekEvMwfATCjyxhULDCDCgw+LCe1ucIEJBtNABBYgTfaIvFMcHCcyY2+Y0mQXl8Dn0SDGXx077/wc2nnZa4MDKxdOVFf44tB6hFOEX3AU+8Txz0SuqxnOPl0xhTBSNMgwEZGB+WILh0qtEBqUnvR0KAGA1iQGJAcCjtcShwGICwLJSNRgBqEa/Htk/vmpy8DHaW3whahHyPwtJygCBBi67/l3ZT5jG3KWknPZ/LOGjtEQYGjX/uj2hiftnrRIBAmBdoDwMOoxXC2sMAJXzOpMNmFsolQPgc+CLUIH7C968kyy5sgN51au8oLD0HvJUrV3a4ntcWQ4sSwMERt3WtoQGB43ocAgkBo0FCIa5R3ck2WsQIvNEKuCdY8hpWnmm5IhB4noHnzHmTxnWoENvzfQstZhds7unptMbH96ME3mXMxeD+ZfSX9RCFpeUAlQVaV90XQ9hbUSFctVcDpA6IQ72fEfJarDUL0vVm3bAfcihNQtCENYbRMgZUxkwjcGrggZ9jKTjrbUo1rRB5KV/r+4GiO9QbLotzEUBOnGo8n6+CCuFwjqcMCz7XIaDUwWGE2oDExEb4DWDMcVjDmHOMTTkhEEJb+ezNR1PyS57yJZfJxcjEWvqKZJ8Iw6nssAM4FuCDYKgTBFWnF8U0cjBDMDB2GBvBNrEBCI+ZZsx8DDShmDadhuG85v5aft4B4ibxHOj4At7+fZTAeLmGSIMsfc1S5i0M6eg2r1IHBE4wTQkk6TSEORwTGDo/Yo9pkBb4xXFY2J9JenE5MLP4LhCc/s+sW5f8UAQQcCMKx4oD/Cj7g4ODmE6uVqKBiMKuNABwYQEwQsfmfD2GEJs0762HsHDzJI9NfKhr5tziPLXj2nthbXoN1u5MqZTlOVCtYGZcXiHSICdAfQIYNFtg+8yDwgg824BJPD5UrK8DWOHryBoEmk2GeIZpE5vz4XPB1cP+5d18DmMgwq1iNeHDZl4mFyKAnAAVieEm/ALrujCAOGwMgTagCOfhPA4XxLgOinCavgaPTWzSphnXHIfvYToU+DwDRGiRZas1Qj85qJTwiSNNDw0NKdKRltNI92ezWayzizHtCBRAmjAURhOH0+YcY4JhMVVxDoALgGDi8AhdpsN0qDzm3CKw6GfW3gvphgDI0W3FQuPGsDXMOo7C0+MAP9HWjh07Zvt6e0fh//ZBC3hYQMEGaYAsEEoCAn6CBgVj5sE5gkITBLtaoySF3ACBTrnD3AhG6OlvmLTJxzxMG4CYGPmonfQzUATfrYJtsvaiyZdFLudw9L70Q9BGsKR/+5W//YqbPnfTyzTTeC4Kv4oDAY8saw87+yDWXKdKCyEF0RDF24BGnzMCG4orSGM+OTYF9djxGAg7hX4xVVHCoSicz9xPkLBcPgfvAEDyDJ+xe3jXrmkgREMNp5ZlODoCjNERalh5slV+f++le2+5cvDKWXLLv3bZf2COhlBoUwXgeMgURoAY0mDA8SFjCi7zhoBBcJQh3CWQTyE3Qm8AgaWExBDPhdM8NvkZ836UR62ExRw04Vk+tIeUldqJR3OVx6MjQyjsRAxHbmLpvi2w8YVeHza6/5v71tyXXPnulY/RQrWubQgz9UjrVX9F8MX/IeeNIyiCgQsnGJCYmOeZ5ppZGCiip9ly7BaP9RguCDLm8uoh8HpcFgSZvd66FpBHm06mD0R/93HOnDfACIMEZVETlUIAKaE8EmYZ3o0Ii3EHxTO9HMORA4TcZxVPw8mMyXsL1YJXyBUqsBW0ybUcmXaUfxMtFnRuq7sgjBMYytFGowZCaBEMi4GiwQHh1gCpAcOAoz5PBEJNc81UDZckxQIM873oBIkJBAjyawqDBOecGjgIEBJWOpECStkPLTKt1HdZxP6g9k1pyy4+coAEDFLyGCDyW3IXVpa1nNc5U+AUa4HQicJTc4A8skdHRw/29vZ+Dx2Gr0MrlAtgxPilToEIEkxp0qSBgbTWGogxwpbTEjkCOBjVi1hrDGoIBIy+heLANm8gahc9XovXSAQHCdfCIPEBlCqIZpoGB/IQHCAPoFUTnvf4fYnE/biTgyqBruUbjgZADHfYrNLk6LaV5a12zQ8+irGWZnyYv4wyX0/xpTmFr3V9mR8DEFaY1hgQcIKCAOGXCKNsxYKg64Io/BBwlKcBgmH0mK/o6bnqXP2Emob5dF7kIUDoensABE0qgqOCmL4MqYhrMKn4PpjSpdSM637ln3ftKq3mq+BJKGrZhqMJEBjCuhmfzALfIpA8A6mhkFltbW1fn5ycfBQyuAZffiwCJIogMeCg9tAAQaw1B2KCg4JOE0oLPNIMdXBAsB2YVHS09dJASLMMgoSgogZhXjr0GkiIqzXSDj+uYWUT3UmTx2PGseriftv+on4IAFOLl210NAGybJl0nH6YvX379grMrJvgP3wao2W5oDSWjQsWojbg0FoD57RJhVhrDsQaHAABv0wLwAHh1gABILjiiQEH7+c9evlR3od8BAg1CLWNbrlCTJDQvMr7vgPtETvgups/jNXfN+MVNi1z8wrs0R8gxk8j6GkK1saNC3ZdXXBfcgpL74NQSeGP2YI80cFhOUBbXp1yyil/A5/hYYCEfMRai0EXOzUJicd1guDqNGMQzSBtDkGoCyQIOzbLkTmYTJrQQmXiWaZBOsb1WZOH+Um4Pw+aQ5mzQM8c6nQPQHsgFvsw63cL3qMRwtPUIJvxwbG0M7ZlCzeTDMBSZ9CVV+L6rV651Z/GOV7khywCSZ1BTytBnqlt27Y50CJ/AoDcTj6S6RDOedOKmUBkrmYwBNg43CyAxHu0uQStoXu/EbPvglpIt3QxBukycL82sRDTodfjuQAM3se+j5p55VJ77PW8j1+/e/dD5+J2AETLA4pZ1uFpAmQTmdE5OCjVCy546Qu3bLG24uPDugjCrXojSLnibPmtGawacMqtn981LG+fxEXWwXy+WvYoOiwHyGd7bGzsm319fTfBF/kjaBNsyCbxWV4ALQCHKQaCzHkgusUK5yjkLo45LKVKwjm9sBxiA5B6OcgLhCy4l8NX2PlIgEBruWXLij3heQ/ek8kMDeEVYFote9/DsPZXAISDDvtQL11/jN0n3l8qffvmNWvWn9fZuc3et6+Ixo40ls+cwpL5re47rHN/c83/3PPBR19wlnT99Kclebt8ZGjDhtgwvojmYVH8tDhA4bNbW1vfOzU19XyA5CKABCv1YskrXNBao1aM+fLwBv31R0yEkeGkuvZAmg46K5sgCwMNh08CCAFFgKCZmQ0F1j5YeaNKvfH2Rx4p/3qAUfNoffty/vMrAGJ+Osf9+yoeb552nFV/PzrKj9qiW320osRj8ryULfk77gBPo/AsOUDh82sO+yaA46cwt06CGeTAYV/AdJ0RmQkKQ2YgI8GRwHks0VPvcT8cQHQ5yKfBhZgAQd8LweFPwhfa7ftv+OLIyIMbga3h4FHI0Rgh/EH6Vb+4CxkmQIfVCNeJvJNfOaDnZjDS1N+vKje6fmgOkJVed3f3SbBavw1NshZ+RRUVhl3Vn7xHSBrn2KlIStZIA6SmObT/gfNGg7DiTeXTB6EWIkC09kAS2sMGOGTE99+1eXT0b4cATtTpYesety7LYHj0K34cnXT6IZZs2PCXsW3briU/FwUMViQmwFTatFE4KhzQIOnq6uqOx+O3QZNcAo1CObYAAuxuGKz2TsCQDEjCACEwSNQeRoMYU42Vb75iNXD4AIgLhz62H2bVft9/09dGR79GzbGlwTQHWKPD0wQI87JxSmc/rPQP1UyA4Qb80gTsPCZ/+dF312GBhEKh8HGA5I/5FPgcDgSeW7FpoBiQPEmDIK8ByCG1B8sCQTWw38XGKF056Pv3ASC/v3Vk5P4NqNNtDVyfWuLBnyic2BzQmoSv2N/ffxlMor+CyXUGHXNoa5cag0ABWWEtQhMrrEEIEBZEQgjuBj4Ajhg7E6c8r5gXuRH23IcF+xI2OjjIpAgg5MJzI7CuKNvuwMBAGqbWu5B+JwYfrtGvj74LgMEDQDxoE8toksUAAUiAKv0nxqZhNgPPwrcpi2yeUerjD6CXnNeH8KzhQLnwsGFDBJDnXtVTEdBlEHQoZqBJXonk70KrnA/zqxOkHXEOaNROOi4SJEZ78DoRwi2fUcjPKyJfRa/5bU+Mj2/HadkY+Bu0ujSQeK6RQwSQ52bt17WJef0eLCaN1q7zcPwimF7rEa+sNZYAEmij13/9AtIPe0rdX1Dqjsndux9GPt3gUgMGQaGPEUch4sBzngMEilEOz+rHbAgaVmpuybMqYlnfRAZHYXlwgHVJQWd8OE1gARAaDNsCTcF8pChEHIg4EHEg4kDEgYgDEQciDkQciDgQcSDiQMSBiAMRByIORByIOBBxIOJAxIGIAxEHIg5EHIg4EHEg4kDEgYgDEQeeBgf+PzAHGiWSVB5fAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "no2_anion = ChemicalSystem()\n", "no2_anion.add_atom(\"N\", coords=(0.0, 0.0, 0.12041))\n", "no2_anion.add_atom(\"O\", coords=(0.0, 1.070642, -0.555172))\n", "no2_anion.add_atom(\"O\", coords=(0.0, -1.070642, -0.555172))\n", "no2_anion.charge = -1\n", "view(no2_anion, width=200, height=200, direction=\"along_x\", guess_bonds=True)" ] }, { "cell_type": "markdown", "id": "75a6c70d", "metadata": {}, "source": [ "## Calculate Vibrational Frequencies\n", "Create the settings objects for the ADF donor/acceptor vibrational frequencies calculations. Run the calculations in parallel." ] }, { "cell_type": "code", "execution_count": 5, "id": "3413ccb6", "metadata": {}, "outputs": [], "source": [ "settings_freq = Settings()\n", "settings_freq.input.adf.symmetry = \"NoSym\"\n", "settings_freq.input.adf.basis.type = \"DZP\"\n", "settings_freq.input.adf.basis.core = \"None\"\n", "settings_freq.input.adf.xc.lda = \"SCF VWN\"\n", "settings_freq.input.ams.Task = \"SinglePoint\"\n", "settings_freq.input.ams.Properties.NormalModes = \"Yes\"\n", "settings_freq.input.adf.title = \"Vibrational frequencies\"\n", "\n", "settings_freq_radical = settings_freq.copy()\n", "settings_freq_radical.input.adf.spinpolarization = 1\n", "settings_freq_radical.input.adf.unrestricted = \"Yes\"\n", "\n", "settings_freq_anion = settings_freq.copy()" ] }, { "cell_type": "code", "execution_count": 6, "id": "42f1e4c0", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[30.03|16:40:16] JOB fsradical STARTED\n", "[30.03|16:40:16] JOB fsanion STARTED\n", "[30.03|16:40:16] JOB fsradical RUNNING\n", "[30.03|16:40:16] JOB fsanion RUNNING\n", "[30.03|16:40:24] JOB fsanion FINISHED\n", "[30.03|16:40:24] JOB fsanion SUCCESSFUL\n", "[30.03|16:40:25] JOB fsradical FINISHED\n", "[30.03|16:40:25] JOB fsradical SUCCESSFUL\n" ] } ], "source": [ "freq_job_radical = AMSJob(molecule=no2_radical, settings=settings_freq_radical, name=\"fsradical\")\n", "freq_job_anion = AMSJob(molecule=no2_anion, settings=settings_freq_anion, name=\"fsanion\")\n", "\n", "freq_results = (freq_job_radical.run(), freq_job_anion.run())" ] }, { "cell_type": "code", "execution_count": 7, "id": "21b88bb8-e040-4b1c-802d-1db291563550", "metadata": {}, "outputs": [], "source": [ "rkf = freq_results[0].rkfpath(file=\"adf\")\n", "!amsspectra {rkf}" ] }, { "cell_type": "markdown", "id": "7ef3f382", "metadata": {}, "source": [ "## Calculate Vibronic Spectra\n", "Use Frank-Condon jobs to calculate the vibronic spectra." ] }, { "cell_type": "code", "execution_count": 8, "id": "943500b6", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[30.03|16:42:14] JOB fcfabs STARTED\n", "[30.03|16:42:14] JOB fcfemi STARTED\n", "[30.03|16:42:14] JOB fcfabs RUNNING\n", "[30.03|16:42:14] JOB fcfemi RUNNING\n", "[30.03|16:42:16] JOB fcfemi FINISHED\n", "[30.03|16:42:16] JOB fcfemi SUCCESSFUL\n", "[30.03|16:42:16] JOB fcfabs FINISHED\n", "[30.03|16:42:16] JOB fcfabs SUCCESSFUL\n" ] } ], "source": [ "def fcf_job(state1, state2, spctype, name):\n", " settings_fcf = Settings()\n", " settings_fcf.input.spectrum.type = spctype\n", " settings_fcf.input.state1 = state1\n", " settings_fcf.input.state2 = state2\n", " return FCFJob(inputjob1=state1, inputjob2=state2, settings=settings_fcf, name=name)\n", "\n", "\n", "freq_radical = freq_results[0].rkfpath(file=\"adf\")\n", "freq_anion = freq_results[1].rkfpath(file=\"adf\")\n", "\n", "fc_abs = fcf_job(freq_radical, freq_anion, \"absorption\", \"fcfabs\")\n", "fc_emi = fcf_job(freq_anion, freq_radical, \"emission\", \"fcfemi\")\n", "\n", "fc_results = (fc_abs.run(), fc_emi.run())" ] }, { "cell_type": "markdown", "id": "062357d2", "metadata": {}, "source": [ "## Calculate Density of States\n", "\n", "Calculate the DOS by computing the overlap of the absorption and emission FCF spectra." ] }, { "cell_type": "code", "execution_count": 9, "id": "586448c2", "metadata": {}, "outputs": [], "source": [ "job = FCFDOS(fc_results[0].kfpath(), fc_results[1].kfpath(), 10000.0, 10000.0)\n", "dos = job.dos()" ] }, { "cell_type": "code", "execution_count": 10, "id": "4ce80a1a", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The density of states is 1.30090295e-08\n" ] } ], "source": [ "print(f\"The density of states is {dos:.8e}\")" ] }, { "cell_type": "code", "execution_count": null, "id": "c8215e21-7e50-493f-b4eb-ecd14efad989", "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 }