quant-met 0.0.23__tar.gz → 0.0.24__tar.gz

{quant_met-0.0.23 → quant_met-0.0.24}/CHANGELOG.rst

@@ -2,6 +2,25 @@
 ..
 .. SPDX-License-Identifier: MIT
 
+.. _changelog-0.0.24:
+
+0.0.24 — 2025-03-13
+-------------------
+
+Changed
+^^^^^^^
+
+- Add orbital positions into the Fourier transform
+
+- Implemented a more efficient version of the gap equation using direkt matrix multiplication
+
+Fixed
+^^^^^
+
+- Fix orbital positions in TBLattice for decorated Graphene
+
+- Fix definition of TBLattice for dressed Graphene
+
 .. _changelog-0.0.23:
 
 0.0.23 — 2025-03-10

{quant_met-0.0.23 → quant_met-0.0.24}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: quant-met
-Version: 0.0.23
+Version: 0.0.24
 Summary: Calculate superconductivity in flat-band systems.
 Author-email: Tjark Sievers <tsievers@physnet.uni-hamburg.de>
 License-File: LICENSE.txt

{quant_met-0.0.23 → quant_met-0.0.24}/docs/source/example_input_dmft.yaml

@@ -8,8 +8,8 @@ model:
   hopping_x: 0
   hopping_x_gr_a: 1
   chemical_potential: 0.0
-  hubbard_int_orbital_basis: [1.0, 1.0, 1.0]
-  lattice_constant: 3
+  hubbard_int_orbital_basis: [-1.0, -1.0, -1.0]
+  lattice_constant: 1.73205080757
   q: [0.0, 0.0]
   beta: 100
 
@@ -17,12 +17,12 @@ control:
   calculation: dmft-scf
   prefix: test
   outdir: test/
-  conv_treshold: 1e-2
+  conv_treshold: 1e-4
   max_iter: 100
   n_spin: 1
   n_success: 1
   wmixing: 0.5
-  n_bath: 3
+  n_bath: 2
   n_iw: 1024
   n_w: 4000
   broadening: 0.005

{quant_met-0.0.23 → quant_met-0.0.24}/docs/source/example_input_scf.yaml

@@ -5,10 +5,10 @@
 model:
   name: DressedGraphene
   hopping_gr: 1
-  hopping_x: 0.01
+  hopping_x: 0
   hopping_x_gr_a: 1
   chemical_potential: 0.0
-  hubbard_int_orbital_basis: [0.001, 0.001, 0.001]
+  hubbard_int_orbital_basis: [0.1 , 0.1, 0.1]
   lattice_constant: 3
   q: [0.0, 0.0]
   beta: 100
@@ -17,7 +17,7 @@ control:
   calculation: scf
   prefix: test
   outdir: test/
-  conv_treshold: 1e-6
+  conv_treshold: 1e-1
   max_iter: 1000
   calculate_additional: true
 

quant_met-0.0.24/docs/source/examples/EG-X bands.pdf.license

@@ -0,0 +1,3 @@
+SPDX-FileCopyrightText: 2025 Tjark Sievers
+
+SPDX-License-Identifier: MIT

quant_met-0.0.24/docs/source/examples/band_structures.ipynb

@@ -0,0 +1,165 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "5a227fd4395b199e",
+   "metadata": {},
+   "source": "# Band structure calculations"
+  },
+  {
+   "metadata": {},
+   "cell_type": "code",
+   "outputs": [],
+   "execution_count": null,
+   "source": [
+    "import matplotlib as mpl\n",
+    "import matplotlib.pyplot as plt\n",
+    "import numpy as np\n",
+    "from quant_met import geometry, mean_field, parameters, plotting"
+   ],
+   "id": "b22b860fe454346"
+  },
+  {
+   "metadata": {},
+   "cell_type": "code",
+   "outputs": [],
+   "execution_count": null,
+   "source": [
+    "mpl.style.use(\"default\")\n",
+    "\n",
+    "graphene_lattice = geometry.GrapheneLattice(lattice_constant=np.sqrt(3))\n",
+    "n = 1000  # number of points across the whole k space path\n",
+    "t_gr = 1\n",
+    "t_x = 0.01\n",
+    "band_path, band_path_plot, ticks, labels = graphene_lattice.generate_high_symmetry_path(\n",
+    "    number_of_points=n\n",
+    ")"
+   ],
+   "id": "f85888e34d3c736f"
+  },
+  {
+   "metadata": {},
+   "cell_type": "code",
+   "outputs": [],
+   "execution_count": null,
+   "source": [
+    "V_list = [0, 0.5, 1.4, 3]\n",
+    "\n",
+    "fig, axs = plt.subplots(nrows=1, ncols=len(V_list), figsize=(len(V_list) * 10, 10))\n",
+    "\n",
+    "for V, ax in zip(V_list, axs, strict=False):\n",
+    "    egx_h = mean_field.hamiltonians.DressedGraphene(\n",
+    "        parameters.hamiltonians.DressedGrapheneParameters(\n",
+    "            hopping_gr=t_gr,\n",
+    "            hopping_x=t_x,\n",
+    "            hopping_x_gr_a=V,\n",
+    "            lattice_constant=graphene_lattice.lattice_constant,\n",
+    "            chemical_potential=0,\n",
+    "            hubbard_int_orbital_basis=[0.0, 0.0, 0.0],\n",
+    "        )\n",
+    "    )\n",
+    "\n",
+    "    band_structure = egx_h.calculate_bandstructure(\n",
+    "        band_path, overlaps=np.array([[0, 0, 1], [1, 0, 0]])\n",
+    "    )\n",
+    "\n",
+    "    bands = band_structure[[\"band_0\", \"band_1\", \"band_2\"]]\n",
+    "    overlaps = band_structure[[\"wx_0\", \"wx_1\", \"wx_2\"]]\n",
+    "\n",
+    "    plotting.plot_bandstructure(\n",
+    "        bands=bands.to_numpy().T,\n",
+    "        overlaps=overlaps.to_numpy().T,\n",
+    "        k_point_list=band_path_plot,\n",
+    "        ticks=ticks,\n",
+    "        labels=labels,\n",
+    "        ax_in=ax,\n",
+    "        fig_in=fig,\n",
+    "        overlap_labels=[r\"$w_{\\mathrm{Gr}_1}$\", r\"$w_X$\"],\n",
+    "    )\n",
+    "\n",
+    "    ax.set_title(f\"V = {V:.2f}\", fontsize=30)\n",
+    "    ax.tick_params(axis=\"x\", labelsize=20)\n",
+    "\n",
+    "fig.savefig(\"EG-X bands.pdf\", bbox_inches=\"tight\")"
+   ],
+   "id": "be9baaf277e035e0"
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 4,
+   "id": "8de39516c2786a21",
+   "metadata": {
+    "ExecuteTime": {
+     "end_time": "2025-03-11T17:32:31.243574Z",
+     "start_time": "2025-03-11T17:32:31.075630Z"
+    }
+   },
+   "outputs": [
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAACVMAAAMzCAYAAACh4xBCAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAPwhJREFUeJzs3X+o3Xd9+PHnTbLcuB/3iu28thqzuOkIC3P0Brumyx+6eUeVQmHQiGDUVViYW2kzZcsKbi3CZYPJfrhkiu2KUCW4qfhHUO9fa7UO1pCMsRY21rIb3Y0hEe6tbqRrevaH316435ukOUlumlMfDzh/3DefzznvC+ePzwue5/MZGwwGgwAAAAAAAAAAAH7ErXu5NwAAAAAAAAAAAHAtEFMBAAAAAAAAAAAkpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEB1CTHVo48+2u23396NN97Y2NhYX/7yl1/ynH/4h39oenq6TZs29aY3vam/+Zu/uZS9AgAAXBbzDAAAMMrMNAAAsPaGjql+8IMf9Na3vrVPfvKTF3X8M88807ve9a527drV0aNH+8M//MPuvvvu/v7v/37ozQIAAFwO8wwAADDKzDQAALD2xgaDweCSTx4b60tf+lJ33HHHeY/5/d///b7yla/01FNPLa/t3bu3f/7nf+5b3/rWpX40AADAZTHPAAAAo8xMAwAAa2PDWn/At771rWZmZlas/fqv/3oPPvhg//u//9uP/diPrTrnzJkznTlzZvnvF154oe9973tdd911jY2NrfWWAQBgKIPBoGeffbYbb7yxdeuGvvkr1zDzDAAAr3TmmVc2Mw0AAK90azHTrHlMdeLEiaamplasTU1N9fzzz3fq1KluuOGGVefMzs52//33r/XWAADgijp+/HhveMMbXu5tcAWZZwAA+FFhnnllMtMAAPCj4krONGseU1Wrfqnw4pMFz/cLhv3797dv377lvxcXF3vjG9/Y8ePHm5iYWLuNAgDAJVhaWmrz5s391E/91Mu9FdaAeQYAgFcy88wrn5kGAIBXsrWYadY8pnrd617XiRMnVqydPHmyDRs2dN11153znPHx8cbHx1etT0xMuFAHAOCa5XEHrzzmGQAAflSYZ16ZzDQAAPyouJIzzZo/AP2WW25pbm5uxdrXv/71duzYcc5ncQMAAFwrzDMAAMAoM9MAAMDwho6pvv/973fs2LGOHTtW1TPPPNOxY8ean5+vfnj71z179iwfv3fv3v7zP/+zffv29dRTT/XQQw/14IMP9pGPfOTK/AcAAAAXyTwDAACMMjMNAACsvaEf8/fEE0/09re/ffnvF5+b/f73v7+HH364hYWF5Yv2qq1bt3b48OHuvffe/vqv/7obb7yxv/zLv+w3fuM3rsD2AQAALp55BgAAGGVmGgAAWHtjg8Fg8HJv4qUsLS01OTnZ4uKi53EDAHDNcb3Khfh+AABwLXO9ykvxHQEA4Fq2FterQz/mDwAAAAAAAAAA4JVITAUAAAAAAAAAAJCYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAACqS4ypDhw40NatW9u0aVPT09M99thjFzz+kUce6a1vfWs//uM/3g033NAHP/jBTp8+fUkbBgAAuBzmGQAAYJSZaQAAYG0NHVMdOnSoe+65p/vuu6+jR4+2a9eubrvttubn5895/De+8Y327NnTXXfd1b/+67/2hS98oX/6p3/qQx/60GVvHgAAYBjmGQAAYJSZaQAAYO0NHVN94hOf6K677upDH/pQ27Zt68///M/bvHlzBw8ePOfx//iP/9jP/MzPdPfdd7d169Z+5Vd+pd/6rd/qiSeeuOzNAwAADMM8AwAAjDIzDQAArL2hYqrnnnuuI0eONDMzs2J9Zmamxx9//Jzn7Ny5s29/+9sdPny4wWDQd7/73f7u7/6ud7/73ef9nDNnzrS0tLTiBQAAcDnMMwAAwCgz0wAAwNUxVEx16tSpzp4929TU1Ir1qampTpw4cc5zdu7c2SOPPNLu3bvbuHFjr3vd63r1q1/dX/3VX533c2ZnZ5ucnFx+bd68eZhtAgAArGKeAQAARpmZBgAAro6hH/NXNTY2tuLvwWCwau1FTz75ZHfffXcf+9jHOnLkSF/96ld75pln2rt373nff//+/S0uLi6/jh8/finbBAAAWMU8AwAAjDIzDQAArK0Nwxx8/fXXt379+lW/cDh58uSqX0K8aHZ2tltvvbWPfvSjVf3iL/5iP/ETP9GuXbv6+Mc/3g033LDqnPHx8cbHx4fZGgAAwAWZZwAAgFFmpgEAgKtjqDtTbdy4senp6ebm5lasz83NtXPnznOe89///d+tW7fyY9avX1/98NcSAAAAV4N5BgAAGGVmGgAAuDqGfszfvn37+sxnPtNDDz3UU0891b333tv8/PzyLWH379/fnj17lo+//fbb++IXv9jBgwd7+umn++Y3v9ndd9/d2972tm688cYr958AAAC8BPMMAAAwysw0AACw9oZ6zF/V7t27O336dA888EALCwtt3769w4cPt2XLlqoWFhaan59fPv4DH/hAzz77bJ/85Cf7vd/7vV796lf3jne8oz/5kz+5cv8FAADARTDPAAAAo8xMAwAAa29sMAL3cV1aWmpycrLFxcUmJiZe7u0AAMAKrle5EN8PAACuZa5XeSm+IwAAXMvW4np16Mf8AQAAAAAAAAAAvBKJqQAAAAAAAAAAABJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEB1iTHVgQMH2rp1a5s2bWp6errHHnvsgsefOXOm++67ry1btjQ+Pt7P/uzP9tBDD13ShgEAAC6HeQYAABhlZhoAAFhbG4Y94dChQ91zzz0dOHCgW2+9tU996lPddtttPfnkk73xjW885zl33nln3/3ud3vwwQf7uZ/7uU6ePNnzzz9/2ZsHAAAYhnkGAAAYZWYaAABYe2ODwWAwzAk333xzN910UwcPHlxe27ZtW3fccUezs7Orjv/qV7/ae97znp5++ule85rXXNIml5aWmpycbHFxsYmJiUt6DwAAWCuuV0eHeQYAAFZyvTpazDQAALDSWlyvDvWYv+eee64jR440MzOzYn1mZqbHH3/8nOd85StfaceOHf3pn/5pr3/963vLW97SRz7ykf7nf/7nvJ9z5syZlpaWVrwAAAAuh3kGAAAYZWYaAAC4OoZ6zN+pU6c6e/ZsU1NTK9anpqY6ceLEOc95+umn+8Y3vtGmTZv60pe+1KlTp/rt3/7tvve97533mdyzs7Pdf//9w2wNAADggswzAADAKDPTAADA1THUnaleNDY2tuLvwWCwau1FL7zwQmNjYz3yyCO97W1v613velef+MQnevjhh8/7y4f9+/e3uLi4/Dp+/PilbBMAAGAV8wwAADDKzDQAALC2hroz1fXXX9/69etX/cLh5MmTq34J8aIbbrih17/+9U1OTi6vbdu2rcFg0Le//e3e/OY3rzpnfHy88fHxYbYGAABwQeYZAABglJlpAADg6hjqzlQbN25senq6ubm5Fetzc3Pt3LnznOfceuut/dd//Vff//73l9f+7d/+rXXr1vWGN7zhErYMAAAwPPMMAAAwysw0AABwdQz9mL99+/b1mc98poceeqinnnqqe++9t/n5+fbu3Vv98Pave/bsWT7+ve99b9ddd10f/OAHe/LJJ3v00Uf76Ec/2m/+5m/2qle96sr9JwAAAC/BPAMAAIwyMw0AAKy9oR7zV7V79+5Onz7dAw880MLCQtu3b+/w4cNt2bKlqoWFhebn55eP/8mf/Mnm5ub63d/93Xbs2NF1113XnXfe2cc//vEr918AAABcBPMMAAAwysw0AACw9sYGg8Hg5d7ES1laWmpycrLFxcUmJiZe7u0AAMAKrle5EN8PAACuZa5XeSm+IwAAXMvW4np16Mf8AQAAAAAAAAAAvBKJqQAAAAAAAAAAABJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEB1iTHVgQMH2rp1a5s2bWp6errHHnvsos775je/2YYNG/qlX/qlS/lYAACAy2aeAQAARpmZBgAA1tbQMdWhQ4e65557uu+++zp69Gi7du3qtttua35+/oLnLS4utmfPnn71V3/1kjcLAABwOcwzAADAKDPTAADA2hsbDAaDYU64+eabu+mmmzp48ODy2rZt27rjjjuanZ0973nvec97evOb39z69ev78pe/3LFjxy76M5eWlpqcnGxxcbGJiYlhtgsAAGvO9eroMM8AAMBKrldHi5kGAABWWovr1aHuTPXcc8915MiRZmZmVqzPzMz0+OOPn/e8v/3bv+0//uM/+qM/+qOL+pwzZ860tLS04gUAAHA5zDMAAMAoM9MAAMDVMVRMderUqc6ePdvU1NSK9ampqU6cOHHOc/793/+9P/iDP+iRRx5pw4YNF/U5s7OzTU5OLr82b948zDYBAABWMc8AAACjzEwDAABXx1Ax1YvGxsZW/D0YDFatVZ09e7b3vve93X///b3lLW+56Pffv39/i4uLy6/jx49fyjYBAABWMc8AAACjzEwDAABr6+J+hvD/XH/99a1fv37VLxxOnjy56pcQVc8++2xPPPFER48e7Xd+53eqeuGFFxoMBm3YsKGvf/3rveMd71h13vj4eOPj48NsDQAA4ILMMwAAwCgz0wAAwNUx1J2pNm7c2PT0dHNzcyvW5+bm2rlz56rjJyYm+pd/+ZeOHTu2/Nq7d28///M/37Fjx7r55psvb/cAAAAXyTwDAACMMjMNAABcHUPdmapq3759ve9972vHjh3dcsstffrTn25+fr69e/dWP7z963e+850++9nPtm7durZv377i/Ne+9rVt2rRp1ToAAMBaM88AAACjzEwDAABrb+iYavfu3Z0+fboHHnighYWFtm/f3uHDh9uyZUtVCwsLzc/PX/GNAgAAXC7zDAAAMMrMNAAAsPbGBoPB4OXexEtZWlpqcnKyxcXFJiYmXu7tAADACq5XuRDfDwAArmWuV3kpviMAAFzL1uJ6dd0VeRcAAAAAAAAAAIARJ6YCAAAAAAAAAABITAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAAlZgKAAAAAAAAAACgElMBAAAAAAAAAABUYioAAAAAAAAAAIBKTAUAAAAAAAAAAFCJqQAAAAAAAAAAACoxFQAAAAAAAAAAQCWmAgAAAAAAAAAAqMRUAAAAAAAAAAAA1SXGVAcOHGjr1q1t2rSp6enpHnvssfMe+8UvfrF3vvOd/fRP/3QTExPdcsstfe1rX7vkDQMAAFwO8wwAADDKzDQAALC2ho6pDh061D333NN9993X0aNH27VrV7fddlvz8/PnPP7RRx/tne98Z4cPH+7IkSO9/e1v7/bbb+/o0aOXvXkAAIBhmGcAAIBRZqYBAIC1NzYYDAbDnHDzzTd30003dfDgweW1bdu2dccddzQ7O3tR7/ELv/AL7d69u4997GMXdfzS0lKTk5MtLi42MTExzHYBAGDNuV4dHeYZAABYyfXqaDHTAADASmtxvTrUnamee+65jhw50szMzIr1mZmZHn/88Yt6jxdeeKFnn32217zmNec95syZMy0tLa14AQAAXA7zDAAAMMrMNAAAcHUMFVOdOnWqs2fPNjU1tWJ9amqqEydOXNR7/Nmf/Vk/+MEPuvPOO897zOzsbJOTk8uvzZs3D7NNAACAVcwzAADAKDPTAADA1TFUTPWisbGxFX8PBoNVa+fy+c9/vj/+4z/u0KFDvfa1rz3vcfv3729xcXH5dfz48UvZJgAAwCrmGQAAYJSZaQAAYG1tGObg66+/vvXr16/6hcPJkydX/RLi/3fo0KHuuuuuvvCFL/Rrv/ZrFzx2fHy88fHxYbYGAABwQeYZAABglJlpAADg6hjqzlQbN25senq6ubm5Fetzc3Pt3LnzvOd9/vOf7wMf+ECf+9zneve7331pOwUAALgM5hkAAGCUmWkAAODqGOrOVFX79u3rfe97Xzt27OiWW27p05/+dPPz8+3du7f64e1fv/Od7/TZz362+uFF+p49e/qLv/iLfvmXf3n5FxOvetWrmpycvIL/CgAAwIWZZwAAgFFmpgEAgLU3dEy1e/fuTp8+3QMPPNDCwkLbt2/v8OHDbdmypaqFhYXm5+eXj//Upz7V888/34c//OE+/OEPL6+///3v7+GHH778/wAAAOAimWcAAIBRZqYBAIC1NzYYDAYv9yZeytLSUpOTky0uLjYxMfFybwcAAFZwvcqF+H4AAHAtc73KS/EdAQDgWrYW16vrrsi7AAAAAAAAAAAAjDgxFQAAAAAAAAAAQGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAAAAAEAlpgIAAAAAAAAAAKjEVAAAAAAAAAAAAJWYCgAAAAAAAAAAoBJTAQAAAAAAAAAAVGIqAAAAAAAAAACASkwFAAAAAAAAAABQiakAAAAAAAAAAAAqMRUAAAAAAPB/7d1vaNXl+wfwa7rpUnBgli4t09DsD1lONBURohb5pQiKhKJ/FDQiNKU/KyMzAqkoyFADU3tiJWmFD6zcg7L5hyJbEU0o1DRJkxmaZVnq/X0Q+vutzfPx7LtzznSvF5wH59792a4jl+fcb7jOOQAAAESEYSoAAAAAAAAAAICIMEwFAAAAAAAAAAAQEYapAAAAAAAAAAAAIsIwFQAAAAAAAAAAQEQYpgIAAAAAAAAAAIgIw1QAAAAAAAAAAAARYZgKAAAAAAAAAAAgIgxTAQAAAAAAAAAARIRhKgAAAAAAAAAAgIgwTAUAAAAAAAAAABARhqkAAAAAAAAAAAAiwjAVAAAAAAAAAABARBimAgAAAAAAAAAAiAjDVAAAAAAAAAAAABFhmAoAAAAAAAAAACAiDFMBAAAAAAAAAABEhGEqAAAAAAAAAACAiDBMBQAAAAAAAAAAEBGGqQAAAAAAAAAAACLCMBUAAAAAAAAAAEBEGKYCAAAAAAAAAACICMNUAAAAAAAAAAAAEWGYCgAAAAAAAAAAICIMUwEAAAAAAAAAAESEYSoAAAAAAAAAAICI6OAw1cKFC2PYsGFRWVkZNTU10djYmHP/unXroqamJiorK2P48OHx2muvdahYAACA/5U8AwAAnM5kGgAAKKy8h6lWrFgRDz/8cMyePTuamppi8uTJccMNN8TOnTvb3b99+/aYOnVqTJ48OZqamuLJJ5+M6dOnx6pVq/7n4gEAAPIhzwAAAKczmQYAAAqvLKWU8rlg/PjxMWbMmFi0aNGJtUsuuSRuvvnmmDdvXpv9jz/+eKxevTq2bNlyYq2uri6+/vrr2LRp0yn9zV9//TWqqqriwIED0a9fv3zKBQCAgnNePX3IMwAA0Jrz6ulFpgEAgNYKcV4tz2fzX3/9FZs3b476+vpW67W1tbFx48Z2r9m0aVPU1ta2Wrv++utjyZIl8ffff0dFRUWbaw4fPhyHDx8+cf/AgQMR8c8/AAAAdDXHz6l5vk+BIpNnAACgLXnm9CHTAABAW4XINHkNU7W0tMTRo0dj4MCBrdYHDhwYe/bsafeaPXv2tLv/yJEj0dLSEtXV1W2umTdvXsydO7fN+vnnn59PuQAAUFT79u2LqqqqUpfBScgzAABwcvJM1yfTAADAyXVmpslrmOq4srKyVvdTSm3Wsva3t37cE088EbNmzTpxf//+/TF06NDYuXOnMEe7fv311zj//PPjxx9/9DHDtKE/yKJHyKJHyHLgwIG44IILon///qUuhVMgz9DVeJ0hix4hix4hF/1BFnnm9CPT0NV4rSGLHiEX/UEWPUKWQmSavIapBgwYED179mzzDoe9e/e2eWfDcYMGDWp3f3l5eZx99tntXtO7d+/o3bt3m/Wqqir/OcipX79+eoST0h9k0SNk0SNk6dGjR6lLIAd5hq7O6wxZ9AhZ9Ai56A+yyDNdn0xDV+e1hix6hFz0B1n0CFk6M9Pk9Zt69eoVNTU10dDQ0Gq9oaEhJk6c2O41EyZMaLN/7dq1MXbs2Ha/ixsAAKAQ5BkAAOB0JtMAAEBx5D2WNWvWrHj99ddj6dKlsWXLlpg5c2bs3Lkz6urqIuKfj3+96667Tuyvq6uLHTt2xKxZs2LLli2xdOnSWLJkSTzyyCOd9ygAAABOgTwDAACczmQaAAAovLy+5i8iYtq0abFv37549tlnY/fu3XH55ZfHmjVrYujQoRERsXv37ti5c+eJ/cOGDYs1a9bEzJkzY8GCBXHeeefF/Pnz45Zbbjnlv9m7d++YM2dOux8rCxF6hNz0B1n0CFn0CFn0yOlDnqEr0iNk0SNk0SPkoj/IokdOLzINXZEeIYseIRf9QRY9QpZC9EhZSil12m8DAAAAAAAAAAA4TeX9NX8AAAAAAAAAAABnIsNUAAAAAAAAAAAAYZgKAAAAAAAAAAAgIgxTAQAAAAAAAAAAREQXGqZauHBhDBs2LCorK6OmpiYaGxtz7l+3bl3U1NREZWVlDB8+PF577bUiVUop5NMf7777blx33XVxzjnnRL9+/WLChAnx0UcfFbFaSiHf55DjNmzYEOXl5XHllVcWtkBKLt8eOXz4cMyePTuGDh0avXv3josuuiiWLl1apGophXx7ZPny5TF69Ojo06dPVFdXx7333hv79u0rUrUU06effho33nhjnHfeeVFWVhbvv/9+5jXOqt2PPEMWmYYsMg1ZZBpykWfIRabhVMg0ZJFpyEWeIYs8QxaZhpMpWZ5JXcDbb7+dKioq0uLFi1Nzc3OaMWNG6tu3b9qxY0e7+7dt25b69OmTZsyYkZqbm9PixYtTRUVFWrlyZZErpxjy7Y8ZM2ak559/Pn3++efpu+++S0888USqqKhIX375ZZErp1jy7ZHj9u/fn4YPH55qa2vT6NGji1MsJdGRHrnpppvS+PHjU0NDQ9q+fXv67LPP0oYNG4pYNcWUb480NjamHj16pFdeeSVt27YtNTY2pssuuyzdfPPNRa6cYlizZk2aPXt2WrVqVYqI9N577+Xc76za/cgzZJFpyCLTkEWmIRd5hiwyDVlkGrLINOQiz5BFniGLTEMupcozXWKYaty4camurq7V2qhRo1J9fX27+x977LE0atSoVmsPPPBAuvrqqwtWI6WTb3+059JLL01z587t7NLoIjraI9OmTUtPPfVUmjNnjoP6GS7fHvnggw9SVVVV2rdvXzHKowvIt0defPHFNHz48FZr8+fPT0OGDClYjXQNp3JQd1btfuQZssg0ZJFpyCLTkIs8Qz5kGtoj05BFpiEXeYYs8gxZZBpOVTHzTMm/5u+vv/6KzZs3R21tbav12tra2LhxY7vXbNq0qc3+66+/Pr744ov4+++/C1YrxdeR/vi3Y8eOxcGDB6N///6FKJES62iPLFu2LLZu3Rpz5swpdImUWEd6ZPXq1TF27Nh44YUXYvDgwTFy5Mh45JFH4o8//ihGyRRZR3pk4sSJsWvXrlizZk2klOLnn3+OlStXxn/+859ilEwX56zavcgzZJFpyCLTkEWmIRd5hkJwXu1eZBqyyDTkIs+QRZ4hi0xDZ+uss2p5ZxeWr5aWljh69GgMHDiw1frAgQNjz5497V6zZ8+edvcfOXIkWlpaorq6umD1Ulwd6Y9/e+mll+L333+P2267rRAlUmId6ZHvv/8+6uvro7GxMcrLS/40SIF1pEe2bdsW69evj8rKynjvvfeipaUlHnzwwfjll198J/cZqCM9MnHixFi+fHlMmzYt/vzzzzhy5EjcdNNN8eqrrxajZLo4Z9XuRZ4hi0xDFpmGLDINucgzFILzavci05BFpiEXeYYs8gxZZBo6W2edVUv+yVTHlZWVtbqfUmqzlrW/vXXODPn2x3FvvfVWPPPMM7FixYo499xzC1UeXcCp9sjRo0fj9ttvj7lz58bIkSOLVR5dQD7PI8eOHYuysrJYvnx5jBs3LqZOnRovv/xyvPHGG975cAbLp0eam5tj+vTp8fTTT8fmzZvjww8/jO3bt0ddXV0xSuU04Kza/cgzZJFpyCLTkEWmIRd5hs7mvNr9yDRkkWnIRZ4hizxDFpmGztQZZ9WSj/sOGDAgevbs2WaqcO/evW2mxY4bNGhQu/vLy8vj7LPPLlitFF9H+uO4FStWxH333RfvvPNOXHvttYUskxLKt0cOHjwYX3zxRTQ1NcVDDz0UEf8cylJKUV5eHmvXro1rrrmmKLVTHB15Hqmuro7BgwdHVVXVibVLLrkkUkqxa9euGDFiREFrprg60iPz5s2LSZMmxaOPPhoREVdccUX07ds3Jk+eHM8995x3YHZzzqrdizxDFpmGLDINWWQacpFnKATn1e5FpiGLTEMu8gxZ5BmyyDR0ts46q5b8k6l69eoVNTU10dDQ0Gq9oaEhJk6c2O41EyZMaLN/7dq1MXbs2KioqChYrRRfR/oj4p93Otxzzz3x5ptv+m7UM1y+PdKvX7/45ptv4quvvjpxq6uri4svvji++uqrGD9+fLFKp0g68jwyadKk+Omnn+K33347sfbdd99Fjx49YsiQIQWtl+LrSI8cOnQoevRofYzq2bNnRPzfdDvdl7Nq9yLPkEWmIYtMQxaZhlzkGQrBebV7kWnIItOQizxDFnmGLDINna3TzqqpC3j77bdTRUVFWrJkSWpubk4PP/xw6tu3b/rhhx9SSinV19enO++888T+bdu2pT59+qSZM2em5ubmtGTJklRRUZFWrlxZqodAAeXbH2+++WYqLy9PCxYsSLt37z5x279/f6keAgWWb4/825w5c9Lo0aOLVC2lkG+PHDx4MA0ZMiTdeuut6dtvv03r1q1LI0aMSPfff3+pHgIFlm+PLFu2LJWXl6eFCxemrVu3pvXr16exY8emcePGleohUEAHDx5MTU1NqampKUVEevnll1NTU1PasWNHSslZFXmGbDINWWQassg05CLPkEWmIYtMQxaZhlzkGbLIM2SRacilVHmmSwxTpZTSggUL0tChQ1OvXr3SmDFj0rp160787O67705Tpkxptf+TTz5JV111VerVq1e68MIL06JFi4pcMcWUT39MmTIlRUSb29133138wimafJ9D/j8H9e4h3x7ZsmVLuvbaa9NZZ52VhgwZkmbNmpUOHTpU5Koppnx7ZP78+enSSy9NZ511Vqqurk533HFH2rVrV5Grphg+/vjjnGcLZ1VSkmfIJtOQRaYhi0xDLvIMucg0nAqZhiwyDbnIM2SRZ8gi03AypcozZSn5nDMAAAAAAAAAAIAe2VsAAAAAAAAAAADOfIapAAAAAAAAAAAAwjAVAAAAAAAAAABARBimAgAAAAAAAAAAiAjDVAAAAAAAAAAAABFhmAoAAAAAAAAAACAiDFMBAAAAAAAAAABEhGEqAAAAAAAAAACAiDBMBQAAAAAAAAAAEBGGqQAAAAAAAAAAACLCMBUAAAAAAAAAAEBEGKYCAAAAAAAAAACIiIj/AjIRJ3FSMqpQAAAAAElFTkSuQmCC",
+      "text/plain": [
+       "<Figure size 3000x1000 with 3 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "import matplotlib as mpl\n",
+    "\n",
+    "mpl.style.use(\"default\")\n",
+    "\n",
+    "V_list = [0.01, 0.5, 2]\n",
+    "\n",
+    "fig, axs = plt.subplots(nrows=1, ncols=len(V_list), figsize=(len(V_list) * 10, 10))\n",
+    "\n",
+    "egx_h = mean_field.hamiltonians.OneBand(\n",
+    "    parameters.hamiltonians.OneBandParameters(\n",
+    "        hopping=t_gr,\n",
+    "        lattice_constant=graphene_lattice.lattice_constant,\n",
+    "        chemical_potential=0,\n",
+    "        hubbard_int_orbital_basis=[0.0],\n",
+    "    )\n",
+    ")\n",
+    "\n",
+    "band_structure = egx_h.calculate_bandstructure(band_path, overlaps=np.array([[0, 0, 1], [1, 0, 0]]))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "d98e492c13d6bf",
+   "metadata": {
+    "ExecuteTime": {
+     "end_time": "2025-03-11T17:32:31.254281Z",
+     "start_time": "2025-03-11T17:32:31.252711Z"
+    }
+   },
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.11.3"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

{quant_met-0.0.23 → quant_met-0.0.24}/pyproject.toml

@@ -4,7 +4,7 @@
 
 [project]
 name = "quant-met"
-version = "0.0.23"
+version = "0.0.24"
 description = "Calculate superconductivity in flat-band systems."
 authors = [{name="Tjark Sievers", email="tsievers@physnet.uni-hamburg.de"}]
 readme = "README.md"
@@ -50,8 +50,11 @@ ignore = [
     "S101", # This checks against using assert
     "PLR0913", # This checks against using more than 5 function arguments
     "PLR0915", # This checks against too long functions
+    "PLR0912", # Check against too many branches
+    "C901", # Check against too complex functions
     "PLR2004", # Check against "magic values"
     "T201", # Warn about print statements
+    "PLW0603" # Warn about global updates
 ]
 
 [tool.ruff.lint.pydocstyle]

{quant_met-0.0.23 → quant_met-0.0.24}/src/quant_met/cli/_utils.py

@@ -41,7 +41,11 @@ def _tbl_factory(h: BaseHamiltonian[GenericParameters]) -> TBLattice:
         0.5 * lattice_constant * np.array([1, np.sqrt(3), 0]),
         0.5 * lattice_constant * np.array([1, -np.sqrt(3), 0]),
     ]
-    orbital_positions = [(0, 0.5, 0), (0, -0.5, 0), (0, -0.5, 0)]
+    orbital_positions = [
+        (0.5 * (np.sqrt(3) - 1), 0, 0),
+        (0.5 * (np.sqrt(3) + 1), 0, 0),
+        (0.5 * (np.sqrt(3) - 1), 0, 0),
+    ]
     hoppings = {
         (0, 0): [
             [0, h.hopping_gr, h.hopping_x_gr_a],
@@ -50,13 +54,13 @@ def _tbl_factory(h: BaseHamiltonian[GenericParameters]) -> TBLattice:
         ],
         (1, 0): [[0, 0, 0], [h.hopping_gr, 0, 0], [0, 0, 0]],
         (-1, 0): [[0, h.hopping_gr, 0], [0, 0, 0], [0, 0, 0]],
-        (0, 1): [[0, 0, 0], [h.hopping_gr, 0, 0], [0, 0, 0]],
-        (0, -1): [[0, h.hopping_gr, 0], [0, 0, 0], [0, 0, 0]],
+        (0, 1): [[0, h.hopping_gr, 0], [0, 0, 0], [0, 0, 0]],
+        (0, -1): [[0, 0, 0], [h.hopping_gr, 0, 0], [0, 0, 0]],
     }
 
     return TBLattice(
         units=basis_vectors,
         hoppings=hoppings,
         orbital_positions=orbital_positions,
-        orbital_names=["C1", "C2", "X"],
+        orbital_names=["A", "B", "X"],
     )

{quant_met-0.0.23 → quant_met-0.0.24}/src/quant_met/cli/dmft.py

@@ -8,7 +8,8 @@ import logging
 from pathlib import Path
 
 from h5 import HDFArchive
-from triqs.gf import Gf, Idx
+from mpi4py import MPI
+from triqs.gf import Gf
 
 from quant_met.cli._utils import _hamiltonian_factory, _tbl_factory
 from quant_met.dmft.dmft_loop import dmft_loop
@@ -43,7 +44,13 @@ def dmft_scf(parameters: Parameters) -> None:
         h0_nambu_k[k][:n_orbitals, :n_orbitals] = enk(k)
         h0_nambu_k[k][n_orbitals:, n_orbitals:] = -enk(-k)
 
-    xmu = -h.hubbard_int_orbital_basis[0] / 2
+    ust = 0
+    jh = 0
+    xmu = (
+        h.hubbard_int_orbital_basis[0] / 2
+        + (tbl.n_orbitals - 1) * ust / 2
+        + (tbl.n_orbitals - 1) * (ust - jh) / 2
+    )
 
     solver = dmft_loop(
         tbl=tbl,
@@ -69,23 +76,21 @@ def dmft_scf(parameters: Parameters) -> None:
     g_iw, g_an_iw = get_gloc(s_iw, s_an_iw, h0_nambu_k, xmu, parameters.control.broadening, kmesh)
     g_w, g_an_w = get_gloc(s_w, s_an_w, h0_nambu_k, xmu, parameters.control.broadening, kmesh)
 
-    n_iw0 = int(0.5 * len(s_iw.mesh))
-    iw_0 = s_iw.mesh[n_iw0].value.imag
-
-    gap = s_an_iw[Idx(0)][0, 0].real / (1 - (s_iw[Idx(0)][0, 0].imag / iw_0))
+    comm = MPI.COMM_WORLD
+    rank = comm.Get_rank()
 
-    data_dir = Path("data/DressedGraphene/dmft/sweep_V/")
-    data_dir.mkdir(parents=True, exist_ok=True)
+    if rank == 0:
+        data_dir = Path("data/DressedGraphene/dmft/sweep_V/")
+        data_dir.mkdir(parents=True, exist_ok=True)
 
-    # Save calculation results
-    result_file = result_path / f"{parameters.control.prefix}.hdf5"
-    with HDFArchive(f"{result_file}", "w") as ar:
-        ar["s_iw"] = s_iw
-        ar["s_an_iw"] = s_an_iw
-        ar["g_iw"] = g_iw
-        ar["g_an_iw"] = g_an_iw
-        ar["g_w"] = g_w
-        ar["g_an_w"] = g_an_w
-        ar["gap"] = gap
+        # Save calculation results
+        result_file = result_path / f"{parameters.control.prefix}.hdf5"
+        with HDFArchive(f"{result_file}", "w") as ar:
+            ar["s_iw"] = s_iw
+            ar["s_an_iw"] = s_an_iw
+            ar["g_iw"] = g_iw
+            ar["g_an_iw"] = g_an_iw
+            ar["g_w"] = g_w
+            ar["g_an_w"] = g_an_w
 
-    logger.info("Results saved to %s", result_file)
+        logger.info("Results saved to %s", result_file)

{quant_met-0.0.23 → quant_met-0.0.24}/src/quant_met/dmft/dmft_loop.py

@@ -18,7 +18,7 @@ from triqs.operators import c, c_dag, dagger, n
 from quant_met.mean_field.hamiltonians import BaseHamiltonian
 from quant_met.parameters import GenericParameters
 
-from .utils import _dmft_weiss_field, get_gloc
+from .utils import _check_convergence, _dmft_weiss_field, get_gloc
 
 logger = logging.getLogger(__name__)
 
@@ -80,8 +80,27 @@ def dmft_loop(
         h_loc[o1, o2] * c_dag(spin, o1) * c(spin, o2) for spin, o1, o2 in product(spins, orbs, orbs)
     )
 
+    ust = 0
+    jh = 0
+    jx = 0
+    jp = 0
+
     # Interaction part
-    hamiltonian += -h.hubbard_int_orbital_basis[0] * sum(n("up", o) * n("dn", o) for o in orbs)
+    hamiltonian += h.hubbard_int_orbital_basis[0] * sum(n("up", o) * n("dn", o) for o in orbs)
+    hamiltonian += ust * sum(
+        int(o1 != o2) * n("up", o1) * n("dn", o2) for o1, o2 in product(orbs, orbs)
+    )
+    hamiltonian += (ust - jh) * sum(
+        int(o1 < o2) * n(s, o1) * n(s, o2) for s, o1, o2 in product(spins, orbs, orbs)
+    )
+    hamiltonian -= jx * sum(
+        int(o1 != o2) * c_dag("up", o1) * c("dn", o1) * c_dag("dn", o2) * c("up", o2)
+        for o1, o2 in product(orbs, orbs)
+    )
+    hamiltonian += jp * sum(
+        int(o1 != o2) * c_dag("up", o1) * c_dag("dn", o1) * c("dn", o2) * c("up", o2)
+        for o1, o2 in product(orbs, orbs)
+    )
 
     # Matrix dimensions of eps and V: 3 orbitals x 2 bath states
     eps = np.array([[-1.0, -0.5, 0.5, 1.0] for _ in range(tbl.n_orbitals)])
@@ -119,8 +138,6 @@ def dmft_loop(
         bath_fitting_params=fit_params,
     )
 
-    gooditer = 0
-    g0_prev = np.zeros((2, 2 * n_iw, tbl.n_orbitals, tbl.n_orbitals), dtype=complex)
     for iloop in range(max_iter):
         print(f"\nLoop {iloop + 1} of {max_iter}")
 
@@ -151,39 +168,11 @@ def dmft_loop(
 
         # Check convergence of the Weiss field
         g0 = np.asarray([g0_iw.data, g0_an_iw.data])
-        errvec = np.real(np.sum(abs(g0 - g0_prev), axis=1) / np.sum(abs(g0), axis=1))
-        # First iteration
-        if iloop == 0:
-            errvec = np.ones_like(errvec)
-        errmin, err, errmax = np.min(errvec), np.average(errvec), np.max(errvec)
-
-        g0_prev = np.copy(g0)
-
-        if err < epsilon:
-            gooditer += 1  # Increase good iterations count
-        else:
-            gooditer = 0  # Reset good iterations count
-
-        conv_bool = ((err < epsilon) and (gooditer > n_success) and (iloop < max_iter)) or (
-            iloop >= max_iter
-        )
+        # Check convergence of the Weiss field
+        g0 = np.asarray([g0_iw.data, g0_an_iw.data])
+        err, converged = _check_convergence(g0, epsilon, n_success, max_iter)
 
-        # Print convergence message
-        if iloop < max_iter:
-            if errvec.size > 1:
-                print(f"max error={errmax:.6e}")
-            print("    " * (errvec.size > 1) + f"error={err:.6e}")
-            if errvec.size > 1:
-                print(f"min error={errmin:.6e}")
-        else:
-            if errvec.size > 1:
-                print(f"max error={errmax:.6e}")
-            print("    " * (errvec.size > 1) + f"error={err:.6e}")
-            if errvec.size > 1:
-                print(f"min error={errmin:.6e}")
-            print(f"Not converged after {max_iter} iterations.")
-
-        if conv_bool:
+        if converged:
             break
 
     return solver