{
"cells": [
{
"cell_type": "markdown",
"source": [
"# Anyonic models"
],
"metadata": {}
},
{
"cell_type": "markdown",
"source": [
"We solve the almost-bosonic anyon model of https://arxiv.org/pdf/1901.10739.pdf"
],
"metadata": {}
},
{
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"n Energy log10(ΔE) log10(Δρ) Δtime\n",
"--- --------------- --------- --------- ------\n",
" 1 +64.36949742807 -2.49 6.13s\n",
" 2 +58.27460242935 0.78 -1.39 72.7ms\n",
" 3 +42.75547208388 1.19 -1.01 14.3ms\n",
" 4 +32.72182473424 1.00 -0.87 40.0ms\n",
" 5 +28.65270254117 0.61 -0.76 12.3ms\n",
" 6 +13.62042284706 1.18 -0.71 16.1ms\n",
" 7 +12.88111147864 -0.13 -0.91 8.53ms\n",
" 8 +10.92717077925 0.29 -0.63 14.1ms\n",
" 9 +9.334307576085 0.20 -0.67 13.8ms\n",
" 10 +7.680212136437 0.22 -0.73 13.8ms\n",
" 11 +6.919436428620 -0.12 -0.66 8.65ms\n",
" 12 +6.389827978826 -0.28 -0.68 12.0ms\n",
" 13 +6.207808846761 -0.74 -0.76 11.9ms\n",
" 14 +6.090523098309 -0.93 -0.77 8.83ms\n",
" 15 +5.954415594367 -0.87 -0.83 12.0ms\n",
" 16 +5.944013546273 -1.98 -0.88 6.68ms\n",
" 17 +5.892941365545 -1.29 -0.90 9.89ms\n",
" 18 +5.826726255386 -1.18 -0.95 6.83ms\n",
" 19 +5.745943672255 -1.09 -0.99 9.91ms\n",
" 20 +5.707653718832 -1.42 -1.01 8.54ms\n",
" 21 +5.659841885200 -1.32 -0.85 11.8ms\n",
" 22 +5.632943051951 -1.57 -1.00 11.7ms\n",
" 23 +5.607223982781 -1.59 -1.25 6.80ms\n",
" 24 +5.595945709876 -1.95 -1.22 11.8ms\n",
" 25 +5.588827308814 -2.15 -1.27 8.57ms\n",
" 26 +5.580525271263 -2.08 -1.36 11.8ms\n",
" 27 +5.571857154876 -2.06 -1.50 6.70ms\n",
" 28 +5.568465971100 -2.47 -1.24 9.93ms\n",
" 29 +5.567657317603 -3.09 -1.57 6.63ms\n",
" 30 +5.565756121210 -2.72 -1.82 9.79ms\n",
" 31 +5.565038182986 -3.14 -1.81 6.63ms\n",
" 32 +5.563665466976 -2.86 -1.60 11.8ms\n",
" 33 +5.562528860565 -2.94 -1.72 11.8ms\n",
" 34 +5.560377725829 -2.67 -1.76 6.83ms\n",
" 35 +5.559276863161 -2.96 -1.70 11.8ms\n",
" 36 +5.558900469181 -3.42 -1.86 6.72ms\n",
" 37 +5.558524933287 -3.43 -1.85 11.8ms\n",
" 38 +5.557951739196 -3.24 -1.72 6.66ms\n",
" 39 +5.557832976259 -3.93 -1.82 9.92ms\n",
" 40 +5.556855126483 -3.01 -1.83 6.73ms\n",
" 41 +5.555854816213 -3.00 -1.63 9.84ms\n",
" 42 +5.554582532236 -2.90 -1.19 6.69ms\n",
" 43 +5.549979991108 -2.34 -1.34 9.92ms\n",
" 44 +5.549594975354 -3.41 -0.83 6.75ms\n",
" 45 +5.547443001357 -2.67 -0.74 9.78ms\n",
" 46 +5.537982479542 -2.02 -0.68 8.64ms\n",
" 47 +5.528001898549 -2.00 -0.68 12.0ms\n",
" 48 +5.500301552405 -1.56 -0.58 8.74ms\n",
" 49 +5.480504836116 -1.70 -0.64 12.0ms\n",
" 50 +5.462094303668 -1.73 -0.67 11.8ms\n",
" 51 +5.433076189719 -1.54 -0.78 8.52ms\n",
" 52 +5.375713620894 -1.24 -0.69 11.8ms\n",
" 53 +5.339099116652 -1.44 -0.75 11.9ms\n",
" 54 +5.299007602051 -1.40 -0.85 8.60ms\n",
" 55 +5.268589761222 -1.52 -0.71 11.8ms\n",
" 56 +5.230538071352 -1.42 -0.66 8.63ms\n",
" 57 +5.200569989734 -1.52 -0.69 11.8ms\n",
" 58 +5.186201854578 -1.84 -0.67 11.6ms\n",
" 59 +5.178026122302 -2.09 -0.60 8.62ms\n",
" 60 +5.166705937078 -1.95 -0.69 11.7ms\n",
" 61 +5.159235757399 -2.13 -0.92 6.68ms\n",
" 62 +5.150561361500 -2.06 -0.73 11.8ms\n",
" 63 +5.142022328686 -2.07 -0.81 11.7ms\n",
" 64 +5.131878680754 -1.99 -0.78 8.70ms\n",
" 65 +5.125306712830 -2.18 -0.80 11.8ms\n",
" 66 +5.117754240311 -2.12 -0.89 11.6ms\n",
" 67 +5.108754060644 -2.05 -0.80 8.76ms\n",
" 68 +5.101050315837 -2.11 -0.76 11.9ms\n",
" 69 +5.093934098088 -2.15 -0.78 8.55ms\n",
" 70 +5.086316484253 -2.12 -0.74 12.0ms\n",
" 71 +5.075851161539 -1.98 -0.82 11.6ms\n",
" 72 +5.062946042065 -1.89 -0.70 8.67ms\n",
" 73 +5.049776569998 -1.88 -0.70 11.9ms\n",
" 74 +5.036114281197 -1.86 -0.70 11.6ms\n",
" 75 +5.021672160622 -1.84 -0.67 8.74ms\n",
" 76 +5.008058393962 -1.87 -0.66 12.0ms\n",
" 77 +4.996245868626 -1.93 -0.65 8.66ms\n",
" 78 +4.986162664338 -2.00 -0.67 11.9ms\n",
" 79 +4.977846698493 -2.08 -0.71 11.8ms\n",
" 80 +4.968023636162 -2.01 -0.69 8.68ms\n",
" 81 +4.957119226039 -1.96 -0.68 11.8ms\n",
" 82 +4.945097415211 -1.92 -0.65 8.52ms\n",
" 83 +4.922012929261 -1.64 -0.57 11.8ms\n",
" 84 +4.903536885185 -1.73 -0.61 11.7ms\n",
" 85 +4.872822998541 -1.51 -0.68 8.60ms\n",
" 86 +4.855520432095 -1.76 -0.74 11.8ms\n",
" 87 +4.834607441567 -1.68 -0.78 11.8ms\n",
" 88 +4.819736764542 -1.83 -0.75 8.62ms\n",
" 89 +4.809283086474 -1.98 -0.83 11.9ms\n",
" 90 +4.800138687829 -2.04 -0.84 8.54ms\n",
" 91 +4.792276413720 -2.10 -0.86 11.6ms\n",
" 92 +4.784736826280 -2.12 -0.86 11.7ms\n",
" 93 +4.779688310993 -2.30 -0.93 6.67ms\n",
" 94 +4.770099838134 -2.02 -0.81 11.9ms\n",
" 95 +4.763014950498 -2.15 -0.68 8.62ms\n",
" 96 +4.748364517153 -1.83 -0.63 12.1ms\n",
" 97 +4.741930961418 -2.19 -0.70 29.8ms\n",
" 98 +4.736145228996 -2.24 -1.17 6.93ms\n",
" 99 +4.728828738111 -2.14 -0.99 12.0ms\n",
" 100 +4.722858671743 -2.22 -1.14 8.51ms\n",
" 101 +4.716337310643 -2.19 -1.06 11.7ms\n",
" 102 +4.712169918229 -2.38 -0.94 11.7ms\n",
" 103 +4.709036010736 -2.50 -0.92 8.67ms\n",
" 104 +4.705871059790 -2.50 -0.95 11.6ms\n",
" 105 +4.703398638534 -2.61 -0.96 8.56ms\n",
" 106 +4.700141117206 -2.49 -0.96 11.8ms\n",
" 107 +4.695721177212 -2.35 -0.92 11.7ms\n",
" 108 +4.692035292275 -2.43 -0.84 8.61ms\n",
" 109 +4.687819958886 -2.38 -0.91 11.8ms\n",
" 110 +4.683794280151 -2.40 -0.81 11.7ms\n",
" 111 +4.679577650317 -2.38 -0.77 8.62ms\n",
" 112 +4.675646356510 -2.41 -0.77 11.9ms\n",
" 113 +4.673253557760 -2.62 -0.93 8.52ms\n",
" 114 +4.670612244034 -2.58 -0.99 11.7ms\n",
" 115 +4.667159396724 -2.46 -0.85 11.8ms\n",
" 116 +4.665245443918 -2.72 -0.89 8.55ms\n",
" 117 +4.663234598496 -2.70 -1.04 11.9ms\n",
" 118 +4.661719188026 -2.82 -0.81 11.8ms\n",
" 119 +4.660864167439 -3.07 -0.95 8.61ms\n",
" 120 +4.660136295826 -3.14 -1.20 11.8ms\n",
" 121 +4.659229474186 -3.04 -1.22 8.62ms\n",
" 122 +4.658174908941 -2.98 -0.94 11.9ms\n",
" 123 +4.657206678834 -3.01 -0.91 11.7ms\n",
" 124 +4.656539691829 -3.18 -1.06 8.75ms\n",
" 125 +4.655718288245 -3.09 -1.19 11.8ms\n",
" 126 +4.655074618551 -3.19 -1.18 11.6ms\n",
" 127 +4.654475324300 -3.22 -1.17 8.58ms\n",
" 128 +4.653904395163 -3.24 -1.14 11.8ms\n",
" 129 +4.653181522739 -3.14 -1.13 8.66ms\n",
" 130 +4.652548595928 -3.20 -0.91 11.8ms\n",
" 131 +4.651798539700 -3.12 -0.90 11.8ms\n",
" 132 +4.650818410843 -3.01 -0.90 8.64ms\n",
" 133 +4.649985963347 -3.08 -0.96 11.8ms\n",
" 134 +4.649183807112 -3.10 -1.17 11.8ms\n",
" 135 +4.648380073304 -3.09 -1.18 8.78ms\n",
" 136 +4.647525884258 -3.07 -1.09 11.8ms\n",
" 137 +4.646816309371 -3.15 -1.08 8.69ms\n",
" 138 +4.646347574989 -3.33 -1.07 12.0ms\n",
" 139 +4.645883982502 -3.33 -1.15 11.7ms\n",
" 140 +4.645320775206 -3.25 -1.18 8.70ms\n",
" 141 +4.644922992984 -3.40 -1.28 11.9ms\n",
" 142 +4.644646376723 -3.56 -1.40 11.6ms\n",
" 143 +4.644327081774 -3.50 -1.32 8.67ms\n",
" 144 +4.643995855926 -3.48 -1.16 11.7ms\n",
" 145 +4.643660605706 -3.47 -1.19 8.52ms\n",
" 146 +4.643423233897 -3.62 -1.25 11.7ms\n",
" 147 +4.643160769994 -3.58 -1.39 11.6ms\n",
" 148 +4.642804953965 -3.45 -1.34 8.65ms\n",
" 149 +4.642562214683 -3.61 -1.47 11.8ms\n",
" 150 +4.642334627774 -3.64 -1.43 8.73ms\n",
" 151 +4.642131262001 -3.69 -1.34 12.0ms\n",
" 152 +4.641971796184 -3.80 -1.41 11.6ms\n",
" 153 +4.641823561638 -3.83 -1.33 8.63ms\n",
" 154 +4.641680379947 -3.84 -1.38 11.8ms\n",
" 155 +4.641564877971 -3.94 -1.50 11.8ms\n",
" 156 +4.641474820629 -4.05 -1.47 8.59ms\n",
" 157 +4.641390287007 -4.07 -1.56 11.8ms\n",
" 158 +4.641330561129 -4.22 -1.58 8.49ms\n",
" 159 +4.641263834241 -4.18 -1.61 11.8ms\n",
" 160 +4.641203411815 -4.22 -1.73 11.7ms\n",
" 161 +4.641162942194 -4.39 -1.84 8.60ms\n",
" 162 +4.641125533739 -4.43 -1.98 11.9ms\n",
" 163 +4.641096411580 -4.54 -1.97 11.6ms\n",
" 164 +4.641059231808 -4.43 -1.86 8.74ms\n",
" 165 +4.641036365711 -4.64 -1.86 11.9ms\n",
" 166 +4.641012712625 -4.63 -2.00 8.58ms\n",
" 167 +4.640989362761 -4.63 -2.11 11.8ms\n",
" 168 +4.640960472489 -4.54 -2.06 11.6ms\n",
" 169 +4.640941703346 -4.73 -2.02 8.79ms\n",
" 170 +4.640922063548 -4.71 -1.97 11.9ms\n",
" 171 +4.640900945749 -4.68 -1.95 11.9ms\n",
" 172 +4.640878312622 -4.65 -1.95 8.68ms\n",
" 173 +4.640849761811 -4.54 -1.94 11.7ms\n",
" 174 +4.640820488212 -4.53 -1.95 8.49ms\n",
" 175 +4.640796527453 -4.62 -1.83 11.7ms\n",
" 176 +4.640773749136 -4.64 -1.90 11.7ms\n",
" 177 +4.640755346773 -4.74 -1.85 8.70ms\n",
" 178 +4.640734824353 -4.69 -1.90 11.8ms\n",
" 179 +4.640720032238 -4.83 -2.00 11.4ms\n",
" 180 +4.640705276384 -4.83 -2.07 8.75ms\n",
" 181 +4.640693964022 -4.95 -2.07 11.9ms\n",
" 182 +4.640685989244 -5.10 -2.01 8.59ms\n",
" 183 +4.640678048050 -5.10 -2.10 11.7ms\n",
" 184 +4.640668246646 -5.01 -2.07 11.6ms\n",
" 185 +4.640660946961 -5.14 -2.10 8.57ms\n",
" 186 +4.640654574525 -5.20 -2.20 11.6ms\n",
" 187 +4.640647310593 -5.14 -2.29 8.59ms\n",
" 188 +4.640641768344 -5.26 -2.14 11.9ms\n",
" 189 +4.640635041669 -5.17 -2.06 11.7ms\n",
" 190 +4.640628757081 -5.20 -2.16 8.53ms\n",
" 191 +4.640620486203 -5.08 -2.05 11.7ms\n",
" 192 +4.640610543550 -5.00 -2.05 11.8ms\n",
" 193 +4.640601292369 -5.03 -1.99 8.67ms\n",
" 194 +4.640593807345 -5.13 -1.90 32.6ms\n",
" 195 +4.640584153370 -5.02 -1.89 8.69ms\n",
" 196 +4.640576025855 -5.09 -1.96 11.9ms\n",
" 197 +4.640567866610 -5.09 -2.01 11.7ms\n",
" 198 +4.640560856329 -5.15 -2.06 8.75ms\n",
" 199 +4.640553183725 -5.12 -2.08 11.8ms\n",
" 200 +4.640546837934 -5.20 -1.99 11.7ms\n",
" 201 +4.640541665762 -5.29 -2.01 8.70ms\n",
" 202 +4.640535835510 -5.23 -2.15 11.7ms\n",
" 203 +4.640530783046 -5.30 -2.19 8.56ms\n",
" 204 +4.640526287877 -5.35 -2.23 11.8ms\n",
" 205 +4.640522778315 -5.45 -2.25 11.8ms\n",
" 206 +4.640520184378 -5.59 -2.33 8.62ms\n",
" 207 +4.640518039249 -5.67 -2.31 11.8ms\n",
" 208 +4.640515350266 -5.57 -2.41 11.6ms\n",
" 209 +4.640513025613 -5.63 -2.49 8.69ms\n",
" 210 +4.640511199883 -5.74 -2.47 11.9ms\n",
" 211 +4.640508984829 -5.65 -2.48 8.58ms\n",
" 212 +4.640507044400 -5.71 -2.47 12.0ms\n",
" 213 +4.640505353849 -5.77 -2.40 11.6ms\n",
" 214 +4.640503174248 -5.66 -2.46 8.67ms\n",
" 215 +4.640500803781 -5.63 -2.49 11.8ms\n",
" 216 +4.640498708346 -5.68 -2.46 8.54ms\n",
" 217 +4.640496745789 -5.71 -2.42 12.0ms\n",
" 218 +4.640494839801 -5.72 -2.40 11.8ms\n",
" 219 +4.640493113532 -5.76 -2.31 8.56ms\n",
" 220 +4.640491301012 -5.74 -2.40 11.9ms\n",
" 221 +4.640489842272 -5.84 -2.50 11.7ms\n",
" 222 +4.640488641514 -5.92 -2.55 8.64ms\n",
" 223 +4.640487716055 -6.03 -2.69 11.9ms\n",
" 224 +4.640487076062 -6.19 -2.63 8.57ms\n",
" 225 +4.640486542766 -6.27 -2.74 11.8ms\n",
" 226 +4.640486152829 -6.41 -2.79 11.7ms\n",
" 227 +4.640485706201 -6.35 -2.94 8.53ms\n",
" 228 +4.640485378963 -6.49 -3.00 12.0ms\n",
" 229 +4.640485102995 -6.56 -3.03 11.9ms\n",
" 230 +4.640484874210 -6.64 -2.96 8.79ms\n",
" 231 +4.640484639727 -6.63 -2.98 11.8ms\n",
" 232 +4.640484384228 -6.59 -2.96 8.51ms\n",
" 233 +4.640484098832 -6.54 -2.95 11.9ms\n",
" 234 +4.640483766620 -6.48 -2.98 11.8ms\n",
" 235 +4.640483451967 -6.50 -2.95 8.64ms\n",
" 236 +4.640483155727 -6.53 -2.78 11.7ms\n",
" 237 +4.640482872362 -6.55 -2.79 11.7ms\n",
" 238 +4.640482544118 -6.48 -2.78 8.64ms\n",
" 239 +4.640482176640 -6.43 -2.81 11.9ms\n",
" 240 +4.640481831183 -6.46 -2.83 8.56ms\n",
" 241 +4.640481462143 -6.43 -2.81 11.8ms\n",
" 242 +4.640481195534 -6.57 -2.83 11.6ms\n",
" 243 +4.640480934298 -6.58 -2.75 8.79ms\n",
" 244 +4.640480638412 -6.53 -2.78 11.9ms\n",
" 245 +4.640480439309 -6.70 -2.91 11.7ms\n",
" 246 +4.640480234806 -6.69 -2.88 8.69ms\n",
" 247 +4.640480073269 -6.79 -2.99 11.8ms\n",
" 248 +4.640479942719 -6.88 -3.02 8.53ms\n",
" 249 +4.640479855452 -7.06 -2.96 11.8ms\n",
" 250 +4.640479745736 -6.96 -3.08 11.7ms\n",
" 251 +4.640479650558 -7.02 -3.10 8.65ms\n",
" 252 +4.640479558344 -7.04 -3.11 11.8ms\n",
" 253 +4.640479457582 -7.00 -3.15 8.62ms\n",
" 254 +4.640479350209 -6.97 -3.07 12.1ms\n",
" 255 +4.640479257560 -7.03 -2.92 12.0ms\n",
" 256 +4.640479118823 -6.86 -2.94 8.58ms\n",
" 257 +4.640478986126 -6.88 -2.89 12.0ms\n",
" 258 +4.640478823391 -6.79 -2.89 11.8ms\n",
" 259 +4.640478664329 -6.80 -2.81 8.99ms\n",
" 260 +4.640478523651 -6.85 -2.81 11.8ms\n",
" 261 +4.640478402063 -6.92 -2.85 8.50ms\n",
" 262 +4.640478286020 -6.94 -2.77 11.9ms\n",
" 263 +4.640478162543 -6.91 -3.00 11.7ms\n",
" 264 +4.640478066594 -7.02 -3.04 8.65ms\n",
" 265 +4.640477977485 -7.05 -3.07 11.9ms\n",
" 266 +4.640477921523 -7.25 -3.13 11.8ms\n",
" 267 +4.640477872842 -7.31 -3.23 8.65ms\n",
" 268 +4.640477836013 -7.43 -3.25 11.7ms\n",
" 269 +4.640477798360 -7.42 -3.45 8.54ms\n",
" 270 +4.640477769127 -7.53 -3.60 11.7ms\n",
" 271 +4.640477744929 -7.62 -3.50 11.6ms\n",
" 272 +4.640477723655 -7.67 -3.58 8.57ms\n",
" 273 +4.640477705576 -7.74 -3.55 11.7ms\n",
" 274 +4.640477683623 -7.66 -3.43 11.6ms\n",
" 275 +4.640477663889 -7.70 -3.54 8.66ms\n",
" 276 +4.640477644288 -7.71 -3.59 11.9ms\n",
" 277 +4.640477623690 -7.69 -3.53 8.56ms\n",
" 278 +4.640477608209 -7.81 -3.54 11.9ms\n",
" 279 +4.640477590737 -7.76 -3.41 11.7ms\n",
" 280 +4.640477574881 -7.80 -3.43 8.71ms\n",
" 281 +4.640477556835 -7.74 -3.47 11.8ms\n",
" 282 +4.640477538253 -7.73 -3.47 11.8ms\n",
" 283 +4.640477520027 -7.74 -3.53 8.71ms\n",
" 284 +4.640477506438 -7.87 -3.53 11.9ms\n",
" 285 +4.640477494334 -7.92 -3.48 8.53ms\n",
" 286 +4.640477483386 -7.96 -3.55 11.9ms\n",
" 287 +4.640477473697 -8.01 -3.61 11.8ms\n",
" 288 +4.640477463703 -8.00 -3.65 8.73ms\n",
" 289 +4.640477456292 -8.13 -3.83 11.9ms\n",
" 290 +4.640477450260 -8.22 -3.76 9.05ms\n",
" 291 +4.640477445832 -8.35 -3.74 29.4ms\n",
" 292 +4.640477441374 -8.35 -3.82 12.5ms\n",
" 293 +4.640477437358 -8.40 -3.88 8.62ms\n",
" 294 +4.640477433769 -8.44 -3.95 11.9ms\n",
" 295 +4.640477430538 -8.49 -4.02 12.0ms\n",
" 296 +4.640477427843 -8.57 -4.01 8.88ms\n",
" 297 +4.640477425674 -8.66 -4.04 12.0ms\n",
" 298 +4.640477422721 -8.53 -3.86 8.56ms\n",
" 299 +4.640477419853 -8.54 -3.93 12.5ms\n",
" 300 +4.640477416271 -8.45 -3.88 11.8ms\n",
" 301 +4.640477412303 -8.40 -3.79 8.62ms\n",
" 302 +4.640477407819 -8.35 -3.73 11.7ms\n",
" 303 +4.640477404150 -8.44 -3.71 11.7ms\n",
" 304 +4.640477400405 -8.43 -3.58 8.64ms\n",
" 305 +4.640477395338 -8.30 -3.74 11.8ms\n",
" 306 +4.640477391559 -8.42 -3.69 8.71ms\n",
" 307 +4.640477387899 -8.44 -3.72 12.1ms\n",
" 308 +4.640477383837 -8.39 -3.74 11.7ms\n",
" 309 +4.640477381421 -8.62 -3.74 8.65ms\n",
" 310 +4.640477378699 -8.57 -3.71 11.8ms\n",
" 311 +4.640477375769 -8.53 -3.77 11.8ms\n",
" 312 +4.640477373500 -8.64 -3.81 8.62ms\n",
" 313 +4.640477371264 -8.65 -3.92 11.7ms\n",
" 314 +4.640477369651 -8.79 -3.88 8.65ms\n",
" 315 +4.640477368291 -8.87 -4.02 11.9ms\n",
" 316 +4.640477367160 -8.95 -3.93 11.7ms\n",
" 317 +4.640477366034 -8.95 -4.11 8.63ms\n",
" 318 +4.640477364796 -8.91 -4.10 11.8ms\n",
" 319 +4.640477363615 -8.93 -4.10 11.6ms\n",
" 320 +4.640477362604 -9.00 -4.06 8.80ms\n",
" 321 +4.640477361376 -8.91 -3.94 11.8ms\n",
" 322 +4.640477360344 -8.99 -3.92 8.60ms\n",
" 323 +4.640477359000 -8.87 -3.86 11.9ms\n",
" 324 +4.640477357368 -8.79 -3.96 11.9ms\n",
" 325 +4.640477355816 -8.81 -3.91 8.70ms\n",
" 326 +4.640477354358 -8.84 -3.84 11.8ms\n",
" 327 +4.640477353075 -8.89 -3.85 8.57ms\n",
" 328 +4.640477351886 -8.92 -3.84 11.9ms\n",
" 329 +4.640477350769 -8.95 -3.90 11.7ms\n",
" 330 +4.640477349811 -9.02 -4.08 8.56ms\n",
" 331 +4.640477349077 -9.13 -4.22 12.0ms\n",
" 332 +4.640477348543 -9.27 -4.24 11.9ms\n",
" 333 +4.640477348190 -9.45 -4.40 8.68ms\n",
" 334 +4.640477347877 -9.50 -4.42 11.9ms\n",
" 335 +4.640477347615 -9.58 -4.58 8.69ms\n",
" 336 +4.640477347392 -9.65 -4.70 11.9ms\n",
" 337 +4.640477347186 -9.69 -4.84 11.7ms\n",
" 338 +4.640477346979 -9.68 -4.81 8.58ms\n",
" 339 +4.640477346842 -9.86 -4.78 11.8ms\n",
" 340 +4.640477346687 -9.81 -4.63 11.6ms\n",
" 341 +4.640477346545 -9.85 -4.76 8.78ms\n",
" 342 +4.640477346388 -9.80 -4.68 12.3ms\n",
" 343 +4.640477346229 -9.80 -4.67 8.61ms\n",
" 344 +4.640477346061 -9.77 -4.69 11.8ms\n",
" 345 +4.640477345911 -9.82 -4.48 11.7ms\n",
" 346 +4.640477345753 -9.80 -4.43 8.61ms\n",
" 347 +4.640477345576 -9.75 -4.43 11.8ms\n",
" 348 +4.640477345403 -9.76 -4.44 11.7ms\n",
" 349 +4.640477345239 -9.78 -4.45 8.65ms\n",
" 350 +4.640477345085 -9.81 -4.53 11.7ms\n",
" 351 +4.640477344973 -9.95 -4.50 8.67ms\n",
" 352 +4.640477344882 -10.04 -4.48 11.9ms\n",
" 353 +4.640477344789 -10.03 -4.58 12.0ms\n",
" 354 +4.640477344708 -10.09 -4.66 8.83ms\n",
" 355 +4.640477344638 -10.15 -4.76 12.1ms\n",
" 356 +4.640477344586 -10.29 -4.74 8.68ms\n",
" 357 +4.640477344538 -10.32 -4.83 11.9ms\n",
" 358 +4.640477344507 -10.50 -4.82 11.7ms\n",
" 359 +4.640477344467 -10.40 -4.83 8.55ms\n",
" 360 +4.640477344431 -10.44 -4.98 11.8ms\n",
" 361 +4.640477344400 -10.51 -4.90 11.6ms\n",
" 362 +4.640477344361 -10.41 -4.90 8.65ms\n",
" 363 +4.640477344322 -10.41 -4.81 11.8ms\n",
" 364 +4.640477344286 -10.45 -4.72 8.56ms\n",
" 365 +4.640477344228 -10.24 -4.60 12.5ms\n",
" 366 +4.640477344170 -10.23 -4.65 11.8ms\n",
" 367 +4.640477344101 -10.16 -4.62 8.73ms\n",
" 368 +4.640477344042 -10.23 -4.51 11.9ms\n",
" 369 +4.640477343987 -10.25 -4.59 11.6ms\n",
" 370 +4.640477343940 -10.33 -4.51 8.61ms\n",
" 371 +4.640477343891 -10.31 -4.58 11.8ms\n",
" 372 +4.640477343847 -10.35 -4.65 8.56ms\n",
" 373 +4.640477343796 -10.30 -4.72 12.0ms\n",
" 374 +4.640477343758 -10.42 -4.73 11.6ms\n",
" 375 +4.640477343726 -10.49 -4.72 8.65ms\n",
" 376 +4.640477343700 -10.59 -4.74 11.8ms\n",
" 377 +4.640477343680 -10.70 -4.81 11.7ms\n",
" 378 +4.640477343660 -10.70 -4.91 8.73ms\n",
" 379 +4.640477343643 -10.76 -5.06 11.9ms\n",
" 380 +4.640477343628 -10.84 -5.08 8.69ms\n",
" 381 +4.640477343616 -10.92 -5.06 12.0ms\n",
" 382 +4.640477343604 -10.90 -5.00 11.8ms\n",
" 383 +4.640477343591 -10.90 -5.05 8.71ms\n",
" 384 +4.640477343578 -10.87 -5.00 12.0ms\n",
" 385 +4.640477343562 -10.81 -5.10 8.54ms\n",
" 386 +4.640477343548 -10.84 -5.01 12.0ms\n",
" 387 +4.640477343535 -10.92 -4.99 26.1ms\n",
" 388 +4.640477343523 -10.91 -4.92 8.83ms\n",
" 389 +4.640477343510 -10.89 -4.87 12.0ms\n",
" 390 +4.640477343498 -10.91 -4.94 12.0ms\n",
" 391 +4.640477343486 -10.93 -5.04 8.65ms\n",
" 392 +4.640477343475 -10.94 -4.98 11.9ms\n",
" 393 +4.640477343466 -11.04 -5.09 8.56ms\n",
" 394 +4.640477343459 -11.18 -5.05 11.9ms\n",
" 395 +4.640477343450 -11.05 -5.04 12.0ms\n",
" 396 +4.640477343445 -11.25 -5.21 8.61ms\n",
" 397 +4.640477343439 -11.23 -5.22 11.8ms\n",
" 398 +4.640477343434 -11.31 -5.34 11.7ms\n",
" 399 +4.640477343430 -11.45 -5.44 8.66ms\n",
" 400 +4.640477343428 -11.65 -5.52 11.8ms\n",
" 401 +4.640477343426 -11.66 -5.56 8.51ms\n",
" 402 +4.640477343424 -11.68 -5.72 11.8ms\n",
" 403 +4.640477343422 -11.82 -5.74 11.7ms\n",
" 404 +4.640477343420 -11.77 -5.79 8.69ms\n",
" 405 +4.640477343419 -11.85 -5.83 11.9ms\n",
" 406 +4.640477343418 -11.89 -5.66 12.1ms\n",
" 407 +4.640477343416 -11.86 -5.66 8.65ms\n",
" 408 +4.640477343415 -11.80 -5.64 11.9ms\n",
" 409 +4.640477343413 -11.79 -5.63 8.62ms\n",
" 410 +4.640477343412 -11.84 -5.60 11.8ms\n",
" 411 +4.640477343410 -11.79 -5.57 12.3ms\n",
" 412 +4.640477343409 -11.87 -5.55 8.78ms\n",
" 413 +4.640477343407 -11.85 -5.43 11.9ms\n",
" 414 +4.640477343405 -11.74 -5.44 11.7ms\n",
" 415 +4.640477343404 -11.75 -5.48 8.57ms\n",
" 416 +4.640477343402 -11.80 -5.49 11.8ms\n",
" 417 +4.640477343401 -11.87 -5.41 8.56ms\n",
" 418 +4.640477343400 -11.96 -5.50 11.8ms\n",
" 419 +4.640477343399 -12.01 -5.45 11.7ms\n",
" 420 +4.640477343398 -11.98 -5.51 8.66ms\n",
" 421 +4.640477343397 -12.06 -5.69 11.8ms\n",
" 422 +4.640477343396 -12.24 -5.69 8.45ms\n",
" 423 +4.640477343396 -12.25 -5.78 11.7ms\n",
" 424 +4.640477343395 -12.36 -5.77 11.7ms\n",
" 425 +4.640477343395 -12.44 -5.75 8.64ms\n",
" 426 +4.640477343394 -12.34 -5.80 11.7ms\n",
" 427 +4.640477343394 -12.29 -5.79 11.6ms\n",
" 428 +4.640477343393 -12.25 -5.78 8.67ms\n",
" 429 +4.640477343393 -12.26 -5.70 11.8ms\n",
" 430 +4.640477343392 -12.31 -5.62 8.67ms\n",
" 431 +4.640477343392 -12.22 -5.57 11.9ms\n",
" 432 +4.640477343391 -12.20 -5.56 11.7ms\n",
" 433 +4.640477343390 -12.11 -5.53 8.61ms\n",
" 434 +4.640477343390 -12.11 -5.59 11.9ms\n",
" 435 +4.640477343389 -12.07 -5.47 11.6ms\n",
" 436 +4.640477343388 -12.25 -5.48 8.59ms\n",
" 437 +4.640477343387 -12.19 -5.46 11.7ms\n",
" 438 +4.640477343387 -12.28 -5.53 8.66ms\n",
" 439 +4.640477343386 -12.29 -5.60 12.0ms\n",
" 440 +4.640477343386 -12.38 -5.80 11.8ms\n",
" 441 +4.640477343386 -12.54 -5.85 8.64ms\n",
" 442 +4.640477343385 -12.67 -5.84 11.8ms\n",
" 443 +4.640477343385 -12.74 -5.89 11.6ms\n",
" 444 +4.640477343385 -12.78 -5.99 8.65ms\n",
" 445 +4.640477343385 -12.74 -6.05 11.7ms\n",
" 446 +4.640477343385 -12.91 -6.09 8.62ms\n",
" 447 +4.640477343385 -12.86 -6.13 12.2ms\n",
" 448 +4.640477343385 -12.98 -6.13 11.8ms\n",
" 449 +4.640477343384 -12.95 -6.02 8.81ms\n",
" 450 +4.640477343384 -13.01 -6.15 12.1ms\n",
" 451 +4.640477343384 -12.97 -6.15 11.7ms\n",
" 452 +4.640477343384 -13.10 -6.20 8.65ms\n",
" 453 +4.640477343384 -12.88 -6.16 11.7ms\n",
" 454 +4.640477343384 -13.01 -6.14 8.54ms\n",
" 455 +4.640477343384 -13.03 -5.97 11.8ms\n",
" 456 +4.640477343384 -12.84 -6.02 11.8ms\n",
" 457 +4.640477343384 -12.88 -6.05 8.62ms\n",
" 458 +4.640477343383 -12.98 -5.99 11.6ms\n",
" 459 +4.640477343383 -12.97 -6.14 8.55ms\n",
" 460 +4.640477343383 -13.30 -6.18 11.8ms\n",
" 461 +4.640477343383 -13.25 -6.25 11.9ms\n",
" 462 +4.640477343383 -13.51 -6.32 8.93ms\n",
" 463 +4.640477343383 -13.32 -6.47 11.9ms\n",
" 464 +4.640477343383 -13.85 -6.56 12.3ms\n",
" 465 +4.640477343383 -13.46 -6.66 8.98ms\n",
" 466 +4.640477343383 -14.35 -6.65 11.8ms\n",
" 467 +4.640477343383 -13.82 -6.68 8.57ms\n",
" 468 +4.640477343383 -13.85 -6.63 11.8ms\n",
" 469 +4.640477343383 -13.62 -6.66 11.7ms\n",
" 470 +4.640477343383 -14.45 -6.65 8.50ms\n",
" 471 +4.640477343383 -13.64 -6.57 11.9ms\n",
" 472 +4.640477343383 -13.82 -6.45 11.6ms\n",
" 473 +4.640477343383 -13.65 -6.48 8.55ms\n",
" 474 +4.640477343383 -13.91 -6.36 11.7ms\n",
" 475 +4.640477343383 -14.27 -6.49 6.59ms\n",
" 476 +4.640477343383 -13.91 -6.42 11.7ms\n",
" 477 +4.640477343383 -13.64 -6.50 11.6ms\n",
" 478 +4.640477343383 -13.50 -6.42 8.60ms\n",
" 479 +4.640477343383 -13.67 -6.38 11.7ms\n",
" 480 +4.640477343383 -13.80 -6.44 8.47ms\n",
" 481 +4.640477343383 -13.80 -6.52 11.9ms\n",
" 482 +4.640477343383 -13.47 -6.54 11.6ms\n",
" 483 +4.640477343383 -13.91 -6.64 8.68ms\n",
" 484 +4.640477343383 + -Inf -6.61 29.5ms\n",
" 485 +4.640477343383 -13.85 -6.66 8.66ms\n",
" 486 +4.640477343383 -13.80 -6.44 12.0ms\n",
" 487 +4.640477343383 -14.75 -6.60 11.6ms\n",
" 488 +4.640477343383 -13.85 -6.67 8.49ms\n",
" 489 +4.640477343383 -14.35 -6.75 15.5ms\n",
" 490 +4.640477343383 -15.05 -6.73 16.2ms\n",
" 491 +4.640477343383 -15.05 -6.68 9.85ms\n",
" 492 +4.640477343383 -14.21 -6.49 8.67ms\n",
" 493 +4.640477343383 -14.45 -6.46 11.7ms\n",
" 494 +4.640477343383 -13.80 -6.41 8.60ms\n",
" 495 +4.640477343383 -14.35 -6.49 11.8ms\n",
" 496 +4.640477343383 -14.10 -6.60 11.7ms\n",
" 497 +4.640477343383 + -Inf -6.66 6.65ms\n",
" 498 +4.640477343383 -13.97 -6.74 11.6ms\n",
" 499 +4.640477343383 -13.67 -6.74 8.60ms\n",
" 500 +4.640477343383 -14.75 -6.48 11.6ms\n",
" 501 +4.640477343383 -13.91 -6.73 11.7ms\n",
" 502 +4.640477343383 -14.15 -6.43 8.58ms\n",
" 503 +4.640477343383 -13.91 -6.56 11.5ms\n",
" 504 +4.640477343383 -15.05 -7.02 11.8ms\n",
" 505 +4.640477343383 -14.35 -6.96 8.72ms\n",
" 506 +4.640477343383 -14.15 -7.12 11.7ms\n",
"┌ Warning: DM not converged.\n",
"└ @ DFTK ~/work/DFTK.jl/DFTK.jl/src/scf/scf_callbacks.jl:60\n",
"e(1,1) / (2π) = 1.2112891860627029\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": "Plot{Plots.GRBackend() n=1}",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deXgUVbr48VOd7oSEQEIwSNgJIAQQdNiioCjbVVmMiAgugDvx8XoRr1cdZhguzKPiMiKKXnEBFVCMXmeugj6jjoKIMDAgsggiYVEQAoQEQrZe6vdHz/Svz6lYvVBQ6dT38/BHn6rq6krTydvnvPWeo+m6LgAAcCqX3RcAAICdCIQAAEcjEAIAHI1ACABwNAIhAMDRCIQAAEcjEAIAHI1ACABwNAIhAMDRCIQAAEdzn5uXCQQCv5s1++EZM87NyyU6n8/ndp+j/5oGw+/3JyUl2X0VicTv97tcLk3T7L6QROLw382MFE/EY3btKq2q8sV02n37fjxyZOs999wT73WdqXP0P1pTU3PYq32+79i5eTkAgOXGds2JeMyNN364ZcvRGE+8raCg3MZAyNAoAMDRnNvHBwBYTtNEwg230yMEADgaPUIAgGU0TYv1Dixdt7kLSSAEAFiKoVEAABIIPUIAgHUS8GYZAiEAwDJx5AiF3TlChkYBAI5GIAQAOBpDowAAy8RRUG/3yCiBEABgnThyhLbfXMPQKADA0egRAgCsoyVeQT2BEABgGS32oU7bh0YJhAAAy8RTR2h3JCRHCABwNHqEAADrxJEjZGgUANCQJFocZGgUAOBs9AgBAJbRNBFzQb3dfUICIQDAOuQIAQBOpok4plijfAIAAPvQIwQAWCaO1SdsRyAEAFgnAXOEDI0CAByNHiEAwDKa0GIth7C7Q0ggBABYKPYcoe05RQIhAMA6CbgeITlCAICj0SMEAFgmjoJ628dGCYQAAMvEUUdodxxkaBQA4Gz0CAEA1knAm2UIhAAAy2ha4k26TSAEAFgmATuE5AgBAM5GjxAAYJ0EnHSbQAgAsE4cOcKzdCVRY2gUAOBo9AgBAJbRYi+Qt/umUQIhAMA6mhZ7OYTdkZBACACwTgLeLEOOEADgaPQIAQCWScAOIYEQAGCdeKZYszsUMjQKAHA0eoQAAOsk4NgogRAAYJk46ghtD4QMjQIAHI1ACACwjha7KM5aVVV19OjROC7n5MmTZWVlyka/33/kyBGv1xtsEggBANbR4vpnasaMGTk5Ob169brsssuOHTtmPGDy5Ml5eXlZWVlr1qwJbfT5fLfcckuHDh06d+58ww031NTUBLd/8803HTp0GDBgQOvWrT/88ENBIAQAWEjT4vln4quvvnrttde+//77Q4cOtWnTZtasWcZj8vPzX3/9dbfbHerkCSEWL168devWn3/++eDBgz/99NPLL78shNB1/bbbbps5c+a+ffuWLl06efLkyspKAiEAoP5asmTJ+PHjc3JyNE27//77lyxZouu6ckxhYeEll1ySlJQUvnHp0qV33313WlpaSkrK1KlTlyxZIoTYsGHD4cOHp0yZIoQYPnx4Tk7OihUrCIQAAMvE1SE06xLu3bv3ggsuCD7u2rVreXn5iRMnormS4uLi8Cfu3bs3eLaOHTt6PJ7w7dGWT5w+ffrxxx//8ssv3W73yJEjp0+fHoy9Bw8efPTRR3fs2NG7d+/HH3+8RYsWUZ4QANAAxV5HqAv/sWPHPvvss/CN3bp1a9OmjRDi5MmTjRs3Dm4MPigrK8vKyop42lOnTqWlpYWeWF5eHjxbaKMQIj09vaysLNpA+Mgjj2zZsmXBggXV1dVTpkxJSkqaPn26EOK6667r37//4sWLn3322YkTJ37++edRnhAAACGE31e+a9euuXPnhm+85pprHnjgASFEixYtQrd9BvuCUfa4srOzg8Ev+MTgs8LPFtzep0+faAPhxo0b77rrrt69ewshxo8fv3HjRiHE+vXrd+/e/fXXX3s8nvnz52dnZ+/YsaN79+5RnhMA0MBEvPnFyO3JGpg/8IMPPqhz74UXXrhhw4bg4w0bNuTm5qanp0dz2uATr7nmmuATe/XqJYTo2bPnjz/+WF5enpGREQgENm7cOH369GhzhAUFBW+//fa+ffu+//77v/zlLwUFBUKI77777uKLLw4OtjZu3LhHjx5btmyJ8oQAgIZHi6eQ0Cxy3nHHHR999FFRUdHWrVtnzpx57733Brfffvvt77zzTvDxp59+WlRUVF1dvWrVqqKiopMnTwohCgsLX3jhhVWrVn3zzTd/+tOfgk/s1KnTlVdeOW3atB9++OF3v/tdVlbWFVdcEW2P8J577vnoo4/69Onj9XpHjBhx7bXXCiFKSkoyMzNDxzRr1qykpKTOp+u67vf7o3wtAEA95PV6Q7eZ/Ko45ho1lZub+957782dO/fEiRNjx46dNm1acHuLFi1CXcO//vWv+/fvHz58+I4dO3bs2HHJJZc0bdp06NChzzzzzKOPPhoIBGbPnj1q1KjgwW+99dZ//dd/jR8/vkuXLh999JGmadEGwnHjxuXn569atcrn802aNOm+++575ZVXmjZtWllZGTqmoqIiIyOjzqdrmqbc2AoASCyRo+DZMWLEiBEjRigbn3jiidDjp556qs4nTpo0adKkScrG7OzsRYsWhW+JamhU1/VVq1bddNNNLpcrOTl5woQJX375pRCiQ4cOu3fvDh4TCASKi4s7dOgQzQkBAA1TPGOjNl9yVIFQ07SuXbuuXLlSCBEIBFauXJmXlyeEGDFixKlTp1asWCGEKCoqSklJueyyy87q5QIA6jOrp1c7F6IdGn311VdvvfXWxYsX19TUtGjRIpiiTElJWbRo0W233ZaZmVlRUbFkyRLGPwEAiSXaQJifn7979+6jR496PJ7wG2RGjhx58ODBw4cP5+TkuN2sbggAztbgF+bNzs42bvR4PG3btrXoegAACSyOpJ/5FGvnAH04AIBl4lihPjFulgEAoKGiRwgAsE6DzxECAGAinsJAu8dGGRoFADgaPUIAgGXiWH3C7pFReoQAAGejRwgAsE4cdYR2dwnpEQIAHI0eIQDAMnHkCG1PEhIIAQCWsjuwxYqhUQCAo9EjBABYJp5Jt+2+W4ZACACwTDx1hHYPpRIIAQCW0YQW+7JKTLEGAIB96BECAKwTx+oTdiMQAgCsk4A5QoZGAQCORo8QAGAZLfZyCNt7hARCAIB1yBECAJwsEesIyRECAByNHiEAwDJxTLFme5eQHiEAwNEIhAAAR2NoFABgmXhuljk7VxI9AiEAwDKaSLxlmBgaBQA4Gj1CIIx+hvutFOFLsu3DSUCd4iiot/vDTCAEAFgn9hyh7QiEAADLxFFHSI4QAAA70SNEQ6ebtM74dDEdG+O3XsMraea7zV4q0YaqgHOJQAgAsExcQ6Nn6VqixdAoAMDR6BEiAZ3RaKfZk2MeOI3pCbGeXf2arJvtVMdNYxhHFfXgKzkajDhmlrF96J5ACACwlN2BLVYEQgCAZTQt5nIIyicAALATPULUP4aEVmwTn5mn/UzPpdexV//VRqSznSnNpCV0+Uu0ppZqxJBQFELouma2P9FGumCnOFafsPsDRiAEAFgnjrlG7cbQKADA0egRAgAsowlNS7TlJwiEsF9dmTn1kKhbhr3y2c1fSzfsVjeoL23pBG5qUtAsSagpWUAlZaj8YYmUhNGUVKhmVoZoe0YH9Vk8K9Tb/YkiEAIALJVgHUJyhAAAZ6NHCACwTDyTbp+lS4kagRA2iJAUNCbqzPebZgH1gOneSBnEM3x6TDQ1C6g05SygS2lKr+1SDzY7c12bdJOdatGh7X/GUJ8k4lyjDI0CAByNHiEAwDpxFNTb3SMkEOLsiFBWYDqiaDzadAQyoIxeyk3D3hieW8cBunJy9UJNrjMidSxUHvx0mTeTTJumzxXGsVPDAfJ1Kj8mI6UIE0+OkDpCAEBDkYAdQnKEAABno0cIALBMHOUTto+nEwhhkTNICpoXPBi3BMybfqnp95vtDfgDJgdHPLl5ijHmHKE8QKNk8pKSpN1KFtDtlvYmueXnupXnqkNBSUnyS6tXJjeVyg1ShlAk2tgoQ6MAAEejRwgAsEw8k26fnSuJHoEQAGAZTVNnQYrmKfYiEMIasSUFTXNpgUiJOr/PLLGn7vVKTZ+y12f23IivZZ4yDBiSnebMSwOVPJ/bIzU9ctOdnCTv1eW9xtWm5BSjUOoO5RSj8r8dS8pQ1IO/eji7ErB+gkAIAKjX1q5d+8wzz5w4ceK666677777jD3OkpKSmTNn7ty5s2fPnrNnz87KyhJCPPnkk3v27Akdk5ub+/DDDwshpk+ffvr06eDGPn363H333QRCAIBlLO8Q/vzzz1dfffVTTz3VrVu3u+66y+12FxYWKsdcf/31Xbp0mTt37vPPP3/TTTd98sknQoi+fft27NgxeMDs2bObNm0afPzmm2/OnDkzJydHCNG2bVtBjxAAYCWr6whfffXVESNG3H333UKIxx9/fMaMGUog3LRp05YtWz777LOUlJQePXq0aNFi165dXbt2HTJkSPCA48eP33rrrcuXLw895eqrr+7SpUuoSSBEnMzXSjJPCgZMJ/wMGBJ1amJPTvt5lSygsrfWL+2tNX+udLDxbH6vUqRolnE0pELlN8Xwy68WDsq1gGpSMEXKAibLzZRGUtPfSPpNTw7IZYNCeJSsoS4f4JZ+TE0puzL9o2f8ExfpbQAkmzdvvvzyy4OP8/Pzd+7cWVlZmZaWFn7AxRdfnJKSIoRIT0/v2bPn5s2bu3btGjrgjTfe6NOnT/fu3UNbHn300dTU1P79+999990pKSkx1BF+++23BQUF3bp1u/LKK7/++uvgxn379l1//fV5eXkTJkw4dOjQmfy0AICGQIvxn6kjR440a9Ys+DiY/Dty5Ej4ASUlJZmZmaFmVlbW4cOHww944403br/99lCzsLDwhhtuGDZs2BtvvDFmzBhd16PtEe7cuXP48OEzZsx44oknfvnll1A0LigoGDFixNNPP/3MM8/ceOONX331VZQnBAA0PHHUEdZUlXz66aedOnUK33jTTTfNmTNHCNGkSZOqqqrgxsrKSiFEKNsXlJ6eXl1dHWqePn06/IB169bt2bNn/PjxoS3B0wohRo0alZOTs3379mgD4ezZs2+++eZp06YJIbp16xbcuHbt2p9++umxxx5zu91PPfXUeeedt3Xr1gsvvDDKcyKB1DVbmPlYqFkzEDAbUVRGI4UQPnl4s7YmIDd9clM+uFpqepW9clMZRxVCeOWhVGVk1VCMITWV4oqIQ4KGoVFlLFRqpsijnSmpUrNRY6mZqgzw+j3KSwfkwVLD/7Y8UioXSMhVHsYlnAwnYzC0QYtjrtHUxi0GDhz40ksvhW9s2bJl8EH79u1DN3/u2bMnPT092C8MCT9A1/W9e/e2b98+tPf111+/8cYbmzRpYnzd5s2bN23atLS0NNqh0Q0bNnTt2vW22267+uqrn3vuOb/fL4TYtm3bRRdd5Ha7hRCpqak9evTYunVrlCcEAEAIIYSWlpaWKwuNO06YMKGoqOj48eNCiJdffnnChAnBQLt06dJNmzYJIUaMGFFWVvbZZ58JIT766KNAIBDKKZ4+fXr58uXh46IlJSXHjh0LPn7ttde8Xm+vXr2i7RH+9NNPTz755AsvvJCRkTF16tQTJ07MmjXr6NGjGRkZoWOaNWtWUlJS59N1XQ/GTgBAgvJ6vR6POpygiqN+wtTQoUPHjBmTl5fXvHnz5OTkjz/+OLh9/vz5EydO/M1vftOoUaNggGzbtu3BgwcXLVoUusiioqKWLVteeumlobPt3Llz5MiRrVq1qqqq0nX97bffzszMjDYQZmRk3HnnnSNHjhRC/Pd///cjjzwya9aszMzMUFmiEOLUqVOhlKZC07SkJPVGNQBAAokcBePKEUa68VhbsGDBrFmzysrKOnfuHBp3Xbt2bWjOo+uvv/6qq676+eef27Vrl5qaGnrulClTpkyZEn62yy+//Pjx4/v27UtNTW3Tpk3wbNEGwk6dOoU6f5mZmRUVFUKIjh077tq1K7gxEAj8+OOPubm5UZ4Q9V3E2cF0s6ZSKmCYI01OCtaaFTwIIWqqpSxgTZV0QHWlstcX/V7lzEpC0XgxMaUMlTWeIq7KlKSsrKTWS8hZwDQ5CygnBdNqksObSsGJcRI7oUt/4JQ/ZJr8h0rdq5lVU7gMf+RYtglxyM7Ozs7ODt+idK4aN24cXjJhIjk5+YILLgjfEm2OcPLkycuWLausrPT7/a+//nqwUHHYsGE1NTXvv/++EOKtt95KT08fOHBglCcEADQ82j/7hDH9s1m0PcI77rhjw4YN7du393g8vXv3Xrx4sRAiOTl5yZIlkyZNuv/++z0ez7Jly5TJeQEAztKAJ912u92vvvrq888/r2lao0aNQtuHDh36008/lZeXZ2ZmxjytDgAAdottirXwJGSIy+X6tXtkkLjMl1USxrWTzJOC/hhmQauW03hCiOrT0paq097wZmWF1KxSmqeVpmnK0PDSaqGh0vTGVEcoNY1fHJV1l5SVlZIbmeUI05pIST7lOtULM64Ppab9NLOmkhNUl2GS9upJhu/HTLHWwMVcR6jZ3SVkrlEAgGUsv2v0HCAQAgCsk4A5Qu5tAQA4Gj1C/ItplZuxBk4tI1SWUvKbFw5KGSylmE/JCAohKk/VhjcrTkrN0yfN9laeMssgKlWGxvRkrXxthjpCORWnJAUDZtkwlzonp7rukidZKpPypEhX0qhS+uVVcplKFraOpKBMuRhXkvQVWSlwVK5T2etyydOWaupLq5ORUlbYsMQx1yg5QgBAw5GAI6MMjQIAnI0eIf4pwtincbey/Lqy0JI8SGiol5CaypRpSsGDMIx2njxRE948JTcryqXm6ZNKrYV0KnVotDK28gllTjXlTVDGh5RlldzJ6tdQj0deVl4dXjYr1VAKVMzHQl2Gkga3cm0es2aS0pTHUV1JcpWIYQRYuTJ16NT23gHOUAJ2CQmEAADLaHXVyEZ+jq0IhAAAy8RRR2h3HCRHCABwNnqE+BURUoSGKdYizKkmpdaU1Y6qK+U0nlwsIYSoKJNzhKU1crM6vHmqTM4RnjKbca3GdJEmIURNtbIMk9nUZcqIkEfOAjZKk2ZBSzIk6jwp8rpLcipOSbYpWVjlypUrMSzwpC4O6kmRtiQ3UprSH4pk+WB/svQmJAXklGGkyhtdmb/N9GD7+w6IiBwhAMDR4qgjtLtclKFRAICj0SMEAFiGSbeRUMznVFPbxg2m6y75zOoIa2vMqveUrJ4QouKkXCkoZwGVpKBSdKjMqabk/JQZ1JS9wrhOk3y8kqgzrxRMSZVSaxlZjYQsPTMlvKnkCJWaRSXPar6ylbL4VHKKmoVVFnVKbSxP56ZUT8ploP4Us9n19CTDJ0dZtsm0rlAXphlEwAoEQgCAZbQEXI+QHCEAwNHoEQIALBNPQb3dQ94EQueKaXLROorBzNddUuoIa5UcoVJHKGWkjHONKvOFnpYLDZXjleSZkrxU1htSKuqU1KYQorZGOl75qZXMnFpHKCf5lKRgTvsmyms1z2ksPV1OMSqrU5UdqwpvKvOv1pqWPyqTuwpDKrS2KoYZVpU5UZU33PjJMU9OUziY6OJZhsnuSMjQKADA0egRAgCsE8fMMnYjEOJfYhqwMgx5qSOlSvmEz2xoVClaMK5Qr8zBpqwjr45PysMcqcrEZh75d1T+KYzLMCmTkylzxRmGRqXnKjOTZZ4nDY227pShvFab3KYmTz95QppJLnmf+ZpN0sHKwLUyUiqMq03Vmv2YysRyAb+yFpUyNGoonxAUSDRo1BECAJxME1qs5RB2x0FyhAAAZ6NHCACwDqtPoMHQ1WakKdbkpt+0zECppvCqKUM1UacUA/hNixaUpGB6RrLclKcxk0sUKivUyo2SgxXSS/uUmcykpvL7rOQX05pIV9K8ZZryWq06SjnC1MbSD1J6RGoq6czy41JS8LQ8z5z56lHC8JaqBRJqUlBO8sknU3OCEcsnlCfYfSc9zpAW+/+h7f/lDI0CAByNHiEAwDKaFnuBvN3DAARCAIB1yBGi4VDrBI0HyC0lURQwm3ZLyTmZF+cJQ8pKuTRltaMm8mJGLVqnhzfPbyNNY6bk7ZTUmjCsT6TUOCorRinZTaW2T8mqKhlEIURyilQamJLqNml6kqWDlRWglK/kyjsWUN9g9T/IME1aDE3lxZhhzWnimWv07FxJ9MgRAgAcjR4hAMAyTLoNAECCoUeIf4pxqtEIqzipE06apqCUujSlKQwJSOXro5ItayIXDrZoLSUF23bJDG9mNJfm/zTmCJVLPX64Mrx59JfT4U0lR6g0K8qlk5cdkyr/hBClJdLKSsoUrMpCS1UVyuJTcnrSZzb/Z0Rn9P3c9oQPECMCIQDASky6DQBwrnhyhGfpUqJGIIRFYht7O6OnavIq80rlgLJ6UWq6NDNZ02ZScUWz7NTwZkojaZRVGNZOUiZsS2lk9huk1Foow6o/7ylXjnfJP1fjptJrVZ6SRlZLDkqjsqfkcddaZaElZRYzw70BLrmWQ23KF6aZNkWkZZXM/+rZ/jcRZyoB6wi5WQYA4Gj0CAEAlknE8gkCIQDAMnHMLGM7AiH+SfnoRpzpSt0if/a1WJpKCsql5pyESx7CV37NlCm+DCslyWs8VflMmrXyglDCME2acuXGadLCVcsnLz0i5QgP7C5Tj5dXVkqTs5vKkk8V5TVyU8oRKtPUKatNud1qTsTtkQ/wmM3fFillKJ864kcn4f5qosEhRwgAcDR6hAAAyyRijpAeIQDA0egR4lco39E0Q7Gfsl9pyl+xlLo0JeekpKySPOr3M/M1hpSkYKU899gJed4yZVklZR0lZaIyYSj+U4r5lAyikln1yenJk2VSVi/pQIXyWkrdoXKpyhvu98qZUXlpJfO3NNlQLqkURHrkBaGUha6UhGKSkuI1zQcLIbRIhYZIdIm2Li+BEABgoQQsqCcQAgAsE0eO0HbkCAEAjkaPEHUzTB9prBs0KyZzJclpJLdZkklZRyk5Rc1gKQe4kqRcmk8usFOWUir5WUrF1cprGzVqLP0KKLV6Qojj8kJLytpJSk5RVzKGctMrFykqs4MKQ8YxJVW6NjUz55becOUtUg5uJJ+qUZpUoSgM74My56ryP6IkIF1u07JCJht1mDgK6m3vQBIIAQDWScAcIUOjAABHo0eIf4lxjjVlNEOZdkuZe0wZTFPuzk9OlZrKkKAQIkU+oLpKairTiSkTlSl7K+SBU+XCvF51aLRCrnkoL5WGRpXjlSoRZXxSGTM0rhqvXGqSWxpKVU7ucZlVRKQ2lgY/GzeRVnRq3FQdGlWOVyo31KFRuZoiSR4GV67TsEiTYRzMWF5hshP1nlZXJiXSU5h0GwDQUGhx1BGenSuJHoEQAGCd+HKEZ7Cy95kjRwgAcDR6hM5lyAmazalmHOsw1Esok6jJ2TI5q6TknJSb+1PT1Y9ldaWU4lIWSzLkCL1yU0oZVp2W9iqJOmOO0LBsk/TSyrui/Fx6slkNiVLwIAxFC8q0Zyly3s6QBZSaTTJTwpvpSrOp1DQ+3ZAjjKGQQ12kqY6PTixZQNtHzRCrOArqNc3eHiGBEABgmXjqCM/OlUSPoVEAgKPRIwQAWCau9QjP0rVEi0DoYKaFg+pO42I6EeoIzeZUU4relIxU43QpIygMk5P5vFJTmZlMl2c2U3KE3hopC6gkBX2GKdaUkyupUKXkMbmR9FxlLSrzYkphTJ2muU2aaelmpYHpGVIWsEkzJWWovsNpcqGhMgeb8v/l9khNtY4w0hRrpnWDgA0IhAAAy8SRI4z4ZaiysvL//u//ysrK/u3f/q1jx451HvO3v/1t586dPXv2vPzyy4Nb9u3b9+OPP4YOGDhwYGpqavDx2rVrv/32286dO48YMUKQIwQA1GfV1dUDBw5ctGjRtm3b+vTp88033xiPeeSRR6ZOnVpcXDxlypTZs2cHNy5fvvyee+5Z+C8VFf+cf//pp5+eMGFCcXHxtGnT7r//fkGPEABgIU3EniM07RIuX77c5XKtXLkyKSmpXbt2c+bMWblyZfgBR44cee65577//vsOHTrceeedffv2vf/++zMzM4UQQ4cOXbhwYfjBp0+f/uMf//j555/36dPnwQcf7NSp03/+538SCPErzKoKhYiYI1QXCZLGHlLkZJhPLonzGYr5fD5pS8BvVnNknq2sdkspwyQ5++hNkprCkN0MNDJbaUkzn3BVLhxU3hMhRLJcOGhYO0mptpTetDS5EDC9qZJBTDbZKwxVicrkrsqlKu+JYXJRITdjSxKSMUx4Vq8+8fHHH48ZMyYpKUkIUVBQ8Nvf/tbn87nd///X4fPPP8/Ly+vQoYMQolu3bm3btl21atW1114rhDh8+PB7773XqlWrAQMGBM+wdu3aJk2a9OnTRwiRk5PTt2/fv/71r7ENjfp8vptuumnq1KmhLT/88MNVV13VsWPH0aNH79+/P6azAQBg7uDBgzk5OcHHrVq18vv9R44cCT/g0KFDrVq1CjVbtWp18OBBIYTH4zl58mRRUdGkSZMGDBhw4sSJXzs4th7h3Llzt27dGgrFuq5fd91148aNe+ONN55++unx48evX78+rp8UANAQxHGzTE1V+f5d2x955JHwjfn5+QUFBf86pxb+IBCQhoiUG8VdLlfwgAceeGD69OlCCK/XO3z48CeeeGLu3LnGg3VdjyEQ7ty5s6ioaPr06fPnzw9uWbNmTUlJycyZM5OSkubMmZOdnf3tt99edNFF0Z8T9Yfy2dX1CGOj6kLkysdLXaFe2qtUDjSShzoDfnVoNGBYsUh6LflK1FIN+bWUagpliXnjFGsB09oMJRdimGdOrpfwxFY+odRmKM3UxmYzrqkDp3JTOVjUsUK9PKeaPKirzqlmuu5SHQvUs0J9wxZ7HWGS25OSktKsWbPwjc2bNw8+yMnJCXUBDx8+7HK5WrZsGX5k+AHBY4J9vtBleDyekSNHfvHFF8GDS0pKwg8ePJyQ15YAABYzSURBVHhwtIEwEAjcddddCxYsCB//3L59+0UXXRQceG3UqFH37t2DW6I8JwCgoYk9R+hJady5c+eHH364zr1Dhw598803Z8yYoWnaxx9/PHjwYI/HI4Q4cuRIWlpakyZNrrjiijvvvPPw4cMtW7bcv3//jz/+eNlllyknWb9+fTCJmJ+ff/To0R07dnTv3r20tHTdunULFy6MNhD+6U9/uuiiiwYOHBgeCI8fP96kSZNQMzMz8+jRo3U+Xdd1n89X5y4AQEJQ7lI5N26++eZnnnlm4sSJXbt2ff7554uKioLbx4wZM3HixGnTprVp02bKlClXX331+PHjly1bdu+992ZnZwshrrnmmtzc3Ozs7LVr13733Xfr1q0TQmRkZDzwwANjx46dPHnyBx98MG7cuM6dO0f1Ix04cGDevHkrVqwoLi4uKSmpqakpLi7u2LFjs2bNQpUZQoiTJ09mZWXVeQZN08792wcAsFA0f8atvmlUpKenr1+//p133iktLf3yyy979eoV3D5nzpz27dsHH7/00kvvv//+999/P2fOnOD9okKI2bNnr127try8/JZbbikoKAh12+bMmTNgwIDNmzc/+OCDN9xwg4iyjvDQoUMpKSljx44VQlRUVJSVlQ0fPnzbtm25ubk7d+7UdV3TNL/fv3v37k6dOsX4DiAx1DHmr0zJpq7KJO1VZhfT5bv3lcSbrqsZLOW1DKUaZrk3ZbYww4xrUtNYueE3LdVQa0iUieXkC1PWojKWTygLMyWbr8pkWlwRaa9xBSjpAOVK1HWXTOdUcymrLNXxR456iYZMq3PtLfOnRDq8WbNmhYWFysbgpDD/OoM2btw45YC+ffv27du3zhOOGjVq1KhRoWZU5RP5+fl7/uXZZ5/t3r37nj17UlNThw0bFggEli5dKoR45ZVXsrKyLr300mhOCABAPRHzcGV6enrr1q3/+WS3+5133pk0adJ9993XsmXLZcuWxbweIwCgIbF8bPTsizkQjhkzZsyYMaHmoEGDiouLq6urGzVqZPIsAIATxDPptt2suYGFKNjwRCgrFEKTKwtdypc6edBdV7JlyjxlQspI6casnHm5nmnhYGpjKQtYU222KpPfF6GEUbk2Q47QLHlpnjIUhqyhmjKUfy4lFaokFJWDk1Pk5aJS1Jd2mycFlcJBdaElpeRUmDWNEu2PJszFtR6hzR8CVp8AADgaJQ0AAOs4IUcIAMCv08yXVarjCXZHQgIh4qUO65suSKSk/dzKc6UElfF3wrCyUgyTi9bWSB9yr7zuklI46Pep+Ul1mlPTuUbVn9qtXLacePOoP6jyg7g9clLQY16VaJbkU56rZFWFIX/pcptVChoWWpKbpqss1XUAYDMCIQDAMloCftchEAIArEOOEA2V8SueeSGByqW0lNE0uWhBUycAU+dvM59iLVledF4e/DSMhUrNgGFCNfMVoBTqZGPqtHNmA7zC8IMYmvKAsNtstFN9rmldR8RLVesj1PIJ5WQRp1hDQ0b5BAAACYYeIQDAMnHMLGP7qAGBEABgHXKEcA7zOdiUCdiEkoJSai3kG/A1oabllBSCIUcoNf1K5YBpFtDvl/caMoK6OueaGfNcmnkezniAkthT1ngyzzgqe82Tl8ZLNc8CRkgKCiDBEAgBAJaJ42YZ2789EQgBAJZJwJFR7hoFADgbPUJY40xShkpO0PjtTDl5wGWaLZMXeXL7pYOVLKC6ypIxRxhDGaEheaamNpWdEXKE5ilGZWKz2J5r+AZuvpRSbElB27/ew3aJ9hkgEAIArJOABfUEQgCAZeKpI7S7B0mOEADgaPQIcVbEljJUnmuoctM1swNcStovSUn7yadSkoK60lQvJqYcocJQfhcpUReheu8M5v+M+NJq2+wruu3f31GvxXHbqN0IhAAAy2ha7Dk/u79bMTQKAHA0eoQ4FyKNlCqHG4YjlZE9ZdF4ZZBQbusuZfDTrHJDN770GQyNqkUIhrZ6uOlYaGwHm79UHSdmLBTWiGNhXts/YARCAIBlNKGp3/kiP8VmBEIAgHUScI41coQAAEejRwgbqCkBNednTGmZZurMM4hCzU/++isb8o3WipAiVLdFSJzElgWM8HMxTRosE3tBve0IhAAAy8SxDJNmXll89jE0CgBwNHqEAADrJODNMgRC1AORkmdq8V+E9Jj+6y01ARmxgtFKZ7hq9xnk8UgB4pyJZ9Lts/ybFxGBEABgmTjqCG3/akaOEADgaPQIAQCWSrAOIYEQiSBCysG8cDBCAlI33x2Bktk4hxNLkfZD/cTCvAAAJBh6hAAAy8RTUH+WLiVqBEIkvshTl4VRb9PWzHbG+NIWP9X2Pw9AHKgjBAA4WhxzjdodCMkRAgAcjR4hAMAymhCx5wht7hISCOEwpr9xdo/QALABQ6MAAEejRwgAsEx86xGepYuJEoEQAGCZOGaWsT0nQSAEAFgnjjpCu5EjBAA4Gj1CAIBl4soRnqVriRaBEABgHVafAAAgsRAIAQCOxtAoAMAymohjGSbqCAEADUUi1hEyNAoAcDR6hAAA67AwLwDAyRJxrlGGRgEAjkaPEABgmThulrF7ZJRACACwVqLlCBkaBQA4Gj1CAICVbC+QjxU9QgCAZYI5wlj/RVRVVVVaWmpygM/nO3bsmN/vV7YfP37cuFFBIAQAWEaLQ6Qe5O9///ucnJy8vLzBgwcfP37ceMCnn37atm3b3/zmN+3bt1+9enVw40MPPdSsWbO8vLzMzMwHHnggEAgEt3fu3DnrX26//XZBIAQA1Gdr1qxZuHDhjh07Dh8+3LJly1mzZikHeL3eSZMmLViw4MCBA3Pnzp08eXKwC9ijR48dO3aUlJRs37793XffXbZsWfD4srKy9evXl5aWlpaWvv7664JACACwkhbXv1+3ZMmS8ePHt2rVStO0//iP/1iyZImu6+EHfPrpp8nJyWPHjhVCTJw4saqqas2aNUKIKVOm5OTkCCHatWt36aWX/vDDD7/2EtEGwscff7xv374tWrTo06fP8uXLQ9u3bds2ePDg5s2bDxs2bPfu3VGeDQDQIMWTIDQNhMXFxV27dg0+7tatW1lZ2YkTJ37tAJfL1aVLl+Li4vADSkpKvvzyy6FDh4a29OnTJzU1dfDgwVu2bBHRB8Ljx48/99xzO3bs+P3vf3/bbbf9/e9/F0IEAoHrr7/+mmuu2b9//8CBA2+88cYozwYAQFAg4C8rK/uH7PDhw8G95eXljRs3Dj4OPigrKwt/+smTJ9PS0kLN9PT08ANqamomTpw4bty4wYMHB7esXLny+PHjx44dGzBgwKhRoyorK6Mtn3j66aeDDwoKCvr167d+/fr+/fuvXr26rKzsoYcecrlcM2bMmDdv3saNG/v27RvPOwEASHxxrEd4qvToli1b7rnnnvCNEydOfPDBB4UQLVq0CAW2YF+wRYsW4UdmZ2eHR74TJ06EDvB6vTfccEPz5s1feOGF0AH9+/cXQng8nieeeGLhwoXffvttzHWEpaWlW7dufeyxx4QQ33//fa9evVwulxAiOTk5Ly9v586dBEIAQPQyz8sZPHjwBx98UOfeCy+8cOPGjcHHGzdu7NixY3p6unLAli1bvF6vx+Opqqratm1br169hBB+v3/SpEk+n++tt95KSkoynrm2ttbr9aakpMQWCH0+3+TJk0ePHj1w4EAhRGlpaZMmTUJ7MzIyjh07VucTdV33+XwxvRYAoF4JBhvzY+KZa9T0+DvuuOPiiy9+//33u3XrNnPmzMLCwtD2ESNG3Hjjjfn5+bm5uY888khhYeG8efMuvvjiCy+8UAhx++23r1u3bsGCBV999ZUQonXr1nl5eVu2bNm4cWO/fv2qqqrmzp3bpUuXXr16xRAIA4HAlClTvF7vwoULg1uysrIqKipCB5SXlzdv3vxXfk7N7WYWGwBIYBGj4NnQqVOnd999d+7cuSdOnLjuuuumT58e3N68efNQavCDDz546KGHrrvuul69eoVu50xPT+/Xr9/ixYuDzSFDhuTl5aWmpn7yySfz589PTU3Nz89/6aWXPB5PtMFJ1/XCwsIDBw58/PHHKSkpwY2dO3fevn27ruuapvl8vh9++KFLly5W/fAAgMQT+3qEEWdku+qqq6666ipl45NPPhl63K5du/ByhqAFCxYYT3XBBRcUFRUpG6O9a7SwsPCrr75asGDBkSNHiouLg1PdDBkyxO12L1y4MBAIzJ8/v2XLlgMGDIjyhACABsjSIsJzI9pAuGbNmpqamoKCguHDhw8fPjxYjZ+UlFRUVPTiiy+mp6cvW7bsnXfesX2hYQAAYhLt0Oi2bdvq3N6/f/9gQSIAAJbfLHMOcAMLAMAyWuw5QtuHEgmEAABL2d3DixWTbgMAHI0eIQDAMvXjPtDYEAgBAJZJxBwhQ6MAAEejRwgAsE4cY6N2j6USCAEA1tFiHuq0fWiUQAgAsEw8BfVn50qiR44QAOBo9AgBAJYiRwgAcCxNaJrtkS1GDI0CAByNHiEAwDKsPgEAcLYEnGONQAgAsAxTrAEAkGDoEQIALBNHjtD2oVQCIQDAOgmYI2RoFADgaPQIAQAWSrybZQiEAADLJOAqTARCAICFEjASkiMEADgaPUIAgGXiKKi3HYEQAGCZRJxrlKFRAICj0SMEAFhG02Iuh7B9KJVACACwlN1DnbFiaBQA4Gj0CAEAVkq0DiGBEABgnURcj5BACACwDjPLAACQWOgRAgAsk4gF9QRCAIBlEnGKNYZGAQCORo8QAGCdBLxZhkAIALCMFnvOz/aRVAIhAMAymtC0GLt4sR5vOXKEAABHo0cIALBOHDlCuxEIAQDWib2O0PbAydAoAMDR6BECACzDpNsAAGcjRwgAcDJNxNEjPEvXEi1yhAAAR6NHCACwTALOsEYgBABYKAEjIUOjAABHo0cIALCMpsVRDkH5BACgwWCFegCAo5EjBAAgsdAjBABYJp6Ceru7hARCAIBlNFafAAAgsRAIAQCOxtAoAMAycS3DdJauJVoEQgCAZeLIEdodBxkaBQA4Gz1CAIB1ErCgnkAIALBO7DlC25OEDI0CAByNHiEAwDJa7B08uzuEBEIAgHXiKZ+wO0l4jgKhy+V66+k/rnr3zXPzcolu//797dq1i31NL+eqqqo6efLk+eefb/eFJJKjR4+mpaU1btzY7gtJJPv27evQoYPdV2GbzTfdNGfOHPNjhnY4L9bTpha33JeSEu9FWUDTdf3cvNL+/fv9fv+5ea1EV1NTk2LrxyLh6Lru9XqTk5PtvpBEUltb6/F4+L4VE4f/bubk5KSmplp+Wl3Xa2trbXxjz10gBACgHuKuUQCAoxEIAQCORiAEADgagRAA4GhJs2bNsvsaHE3X9Y0bN7755pt/+ctfDh482L17d7f7nzUtJ06ceP75599//31d1zt37mzvddZPJSUlS5cuTUtLy87ODm7Zs2fPvHnzPvnkk+zs7JYtW9p7efWNz+dbunTp4sWLN2zYkJ2dHXzTdF1/++23Fy1atHv37l69eoU+fhBCeL3eZcuWvfXWWxs2bGjbtm2zZs2C20tLS4O/m0IIfjcbAHqENjt06NC4ceNKS0vbtm37yiuvDB061OfzCSG8Xu9ll122efPm3Nzcu+6669VXX7X7SuujwsLChx56aNWqVcHmgQMH+vfvX1tbm5WVNXjw4L///e/2Xl694vP5rrrqqv/5n/9p06aN1+tdt25dcPtvf/vbxx9/vEuXLh9++OHYsWPtvcj65tZbb33xxRd79ux56tSp3r177969WwhRW1s7aNCgLVu2dOzY8Y477li0aJHdl4kzpsNWXq/X6/UGH586dapRo0b/+Mc/dF0vKirKy8vz+/26rn/44Ye5ubnBxwhZvnz5hAkTLr/88pdeeim45eGHH7755puDj//whz+MGzfOvqurd1588cXevXuHPmxB5eXl6enp27Zt03W9srIyIyNj06ZNNl1gveP3+z0ez5YtW4LNIUOGzJ8/X9f1t99+u2fPnoFAQNf1P//5z126dAk+RuKiR2gzt9sdGozy+Xw+ny89PV0IsXr16qFDh7pcLiHE8OHD9+7d+9NPP9l5ofXM8ePHZ82aNW/evPCNq1evHjFiRPDx8OHDQz1FCCFWrlx5yy23rFy58tlnn12/fn1w46ZNm5o2bdqjRw8hRGpq6qBBg3jTQlwuV9euXTdv3iyEKC8v37NnT/CNWr169bBhw4ITEQwfPnz37t2HDh2y+VpxZgiE9ci0adNGjx59wQUXCCF++eWXUN4rJSUlIyPjl19+sfXq6pf777//oYceUuZUC3/TWrRocezYsdraWjuurj7au3fvwoULP/jgg1OnThUUFLz44otCiMOHD4feMSHE+eefz9/0cP/7v//7hz/8oWvXru3bty8sLBwyZIiQP2ZpaWnp6en8biY6EuP1xaxZszZt2vTll18Gm263O3xGOuYPC7dixYpffvllypQpyna32x3MsAohfD5fUlJSUlLSub64+srlcuXl5QUTWn369JkyZcq9995r/Jg5ef4whd/vnzx5ckFBwdSpU/fu3Xv33XcPGDDgiiuuCP+YCSF8Ph+/m4mOQFgvzJ0799133/3iiy+ysrKCW1q3bn3w4MHg47KystOnT7dq1cq+C6xf3nrrrX379vXr108IsWvXrn379h06dGj27NmtW7cOdWgOHjzYsmVLAmFImzZtunfvHnzco0ePY8eOBT9Uv/zySyAQCA7CHzx4MD8/39bLrEe+++67f/zjH6tXr3a73d26dZswYcKiRYuuuOKK8I/ZsWPHqqur+d1MdAyN2m/evHmvvfbaZ599Fj7QN3r06I8//riiokII8d577/Xr149igJDHHnusqKjo5Zdffvnll7t06XLzzTfffvvtQojRo0cXFRXpui6EKCoqGj16tN1XWo8UFBSEUoPffPNNhw4dGjdu3K9fP4/H87e//U0IcejQoXXr1o0cOdLWy6xHmjdv7vV6f/7552Bzz5495513nhBi9OjRK1asOH36tBDivffeu+SSS4LbkbiYdNtmP/74Y5cuXTp06NC8efPglqeeeurKK68UQowdO3bPnj0XXXTRihUrli9fPnToUFuvtJ4aPHjwxIkTp06dKoQoLy8fNGjQ+eefn5mZuW7dujVr1jh5xRxFZWXlkCFDUlNTc3NzP/zww9deey34RWHx4sUPP/zwyJEjV69eXVBQ8PTTT9t9pfXIv//7v//5z38eNWpUcXHxjh07vv7663bt2um6fu211x44cKB3794rVqwoKioK/sIicREIbVZdXb19+/bwLZ06dcrMzBRCBAKBVatWHTlyZNCgQW3atLHpAuu7Xbt2ZWVlhW5eqK6u/uyzz2pqaoYNG5aRkWHvtdU3Xq/3iy++qKioyM/PDx/N27Vr16ZNm3JzcwcMGGDj5dVP27Zt27FjR7NmzQYNGhRagSj0u3nZZZe1bt3a3ivEmSMQAgAcjRwhAMDRCIQAAEcjEAIAHI1ACABwNAIhAMDRCIQAAEcjEAIAHI1ACABwNAIhAMDRCIQAAEcjEAIAHO3/AXvEslUr126ZAAAAAElFTkSuQmCC",
"text/html": [
"\n",
"\n"
],
"image/svg+xml": [
"\n",
"\n"
]
},
"metadata": {},
"execution_count": 1
}
],
"cell_type": "code",
"source": [
"using DFTK\n",
"using StaticArrays\n",
"using Plots\n",
"\n",
"# Unit cell. Having one of the lattice vectors as zero means a 2D system\n",
"a = 14\n",
"lattice = a .* [[1 0 0.]; [0 1 0]; [0 0 0]];\n",
"\n",
"# Confining scalar potential\n",
"pot(x, y, z) = ((x - a/2)^2 + (y - a/2)^2);\n",
"\n",
"# Parameters\n",
"Ecut = 50\n",
"n_electrons = 1\n",
"β = 5;\n",
"\n",
"# Collect all the terms, build and run the model\n",
"terms = [Kinetic(; scaling_factor=2),\n",
" ExternalFromReal(X -> pot(X...)),\n",
" Anyonic(1, β)\n",
"]\n",
"model = Model(lattice; n_electrons, terms, spin_polarization=:spinless) # \"spinless electrons\"\n",
"basis = PlaneWaveBasis(model; Ecut, kgrid=(1, 1, 1))\n",
"scfres = direct_minimization(basis, tol=1e-14) # Reduce tol for production\n",
"E = scfres.energies.total\n",
"s = 2\n",
"E11 = π/2 * (2(s+1)/s)^((s+2)/s) * (s/(s+2))^(2(s+1)/s) * E^((s+2)/s) / β\n",
"println(\"e(1,1) / (2π) = \", E11 / (2π))\n",
"heatmap(scfres.ρ[:, :, 1, 1], c=:blues)"
],
"metadata": {},
"execution_count": 1
}
],
"nbformat_minor": 3,
"metadata": {
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.11.4"
},
"kernelspec": {
"name": "julia-1.11",
"display_name": "Julia 1.11.4",
"language": "julia"
}
},
"nbformat": 4
}