Spin One.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from KitaevHoneycomb import run_atomic, save_data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "Jx = 1.0\n",
    "Jy = 1.0\n",
    "Jz = 1.0\n",
    "\n",
    "S = 1.0\n",
    "\n",
    "Lx = 2\n",
    "Ly = 4\n",
    "shift = 0\n",
    "bc = 'periodic'\n",
    "bc_MPS = 'infinite'\n",
    "\n",
    "chi_list = [int(2**i) for i in range(2, 9)]\n",
    "\n",
    "result_dict = {}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "chi_list = [int(2**i) for i in range(2, 11)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'shift'=0 for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  512 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -0.6401133650166546, S = nan, norm_err = 2.1e-06\n",
      "Current memory usage 782.9 MB, time elapsed: 790.6 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.4924e-05, max_E_trunc = 5.0301e-05\n",
      "MPS bond dimensions: [512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -0.6401137494566491, S = nan, norm_err = 2.0e-07\n",
      "Current memory usage 782.9 MB, time elapsed: 1118.4 s\n",
      "Delta E = -3.8444e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.4640e-05, max_E_trunc = 4.9520e-05\n",
      "MPS bond dimensions: [512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -0.6401137765582803, S = 2.0908915971324720, norm_err = 7.7e-08\n",
      "Current memory usage 782.9 MB, time elapsed: 1315.7 s\n",
      "Delta E = -2.7102e-09, Delta S = nan (per sweep)\n",
      "max_trunc_err = 7.3096e-06, max_E_trunc = 4.9466e-05\n",
      "MPS bond dimensions: [512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -0.6401137600911682, S = 2.0909575328583321, norm_err = 9.3e-08\n",
      "Current memory usage 782.9 MB, time elapsed: 1500.2 s\n",
      "Delta E = 1.6467e-09, Delta S = 6.5936e-06 (per sweep)\n",
      "max_trunc_err = 7.3111e-06, max_E_trunc = 4.9475e-05\n",
      "MPS bond dimensions: [512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 2256, maximum chi = 512\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000053+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512, 512]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=512 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'shift'=0 for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  1024 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -0.6400875545036975, S = nan, norm_err = 2.2e-06\n",
      "Current memory usage 2624.6 MB, time elapsed: 4692.4 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.4805e-06, max_E_trunc = 8.5601e-06\n",
      "MPS bond dimensions: [1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -0.6401413726991976, S = nan, norm_err = 9.1e-08\n",
      "Current memory usage 2624.6 MB, time elapsed: 6794.0 s\n",
      "Delta E = -5.3818e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.4014e-06, max_E_trunc = 8.3584e-06\n",
      "MPS bond dimensions: [1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -0.6401414096337135, S = nan, norm_err = 2.2e-08\n",
      "Current memory usage 2624.6 MB, time elapsed: 8867.9 s\n",
      "Delta E = -3.6935e-09, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.3956e-06, max_E_trunc = 8.3532e-06\n",
      "MPS bond dimensions: [1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -0.6401414102433378, S = 2.1256223689047840, norm_err = 1.2e-08\n",
      "Current memory usage 2624.6 MB, time elapsed: 10129.2 s\n",
      "Delta E = -6.0962e-11, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.1971e-06, max_E_trunc = 8.3516e-06\n",
      "MPS bond dimensions: [1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -0.6401414067152444, S = 2.1256485271645329, norm_err = 2.7e-08\n",
      "Current memory usage 2624.6 MB, time elapsed: 11458.9 s\n",
      "Delta E = 3.5281e-10, Delta S = 2.6158e-06 (per sweep)\n",
      "max_trunc_err = 1.2104e-06, max_E_trunc = 8.4467e-06\n",
      "MPS bond dimensions: [1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024]\n",
      "================================================================================\n",
      "DMRG finished after 50 sweeps.\n",
      "total size = 2736, maximum chi = 1024\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000042+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024, 1024]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=1024 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    }
   ],
   "source": [
    "initial_psi = None\n",
    "for chi in chi_list:\n",
    "    param = (Jx, Jy, Jz, Lx, Ly, chi)\n",
    "    if param in result_dict:\n",
    "        result = result_dict[param]\n",
    "    else:\n",
    "        result = run_atomic(\n",
    "            Jx=Jx, Jy=Jy, Jz=Jz,\n",
    "            Lx=Lx, Ly=Ly, \n",
    "            S=S,\n",
    "            shift=shift,\n",
    "            chi=chi,\n",
    "            bc=bc, bc_MPS=bc_MPS,\n",
    "            initial_psi=initial_psi,\n",
    "        )\n",
    "    initial_psi = result['psi'].copy()\n",
    "    result_dict[param] = result"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2.70183595385866\n",
      "2.1206040326889926\n",
      "2.4034627364208627\n",
      "3.2606305864847034\n",
      "4.6860632932379875\n",
      "7.244590060174479\n",
      "10.806348139717414\n",
      "14.52775472232833\n",
      "17.925858676965426\n"
     ]
    }
   ],
   "source": [
    "for key, element in result_dict.items():\n",
    "#     Jx, Jy, Jz, chi = key\n",
    "    psi = element['psi']\n",
    "    print(psi.correlation_length())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "# some api for the file operation\n",
    "import h5py\n",
    "from tenpy.tools import hdf5_io\n",
    "import os.path\n",
    "\n",
    "# functools\n",
    "from functools import wraps\n",
    "\n",
    "# path\n",
    "from pathlib import Path\n",
    "\n",
    "def save_data(\n",
    "    result, \n",
    "    S,\n",
    "    Lx, \n",
    "    Ly, \n",
    "    chi,\n",
    "    Jx, \n",
    "    Jy, \n",
    "    Jz, \n",
    "    prefix='data/', \n",
    "):\n",
    "    file_name = full_path(chi, Jx, Jy, Jz, L, prefix=prefix, suffix='.h5')\n",
    "    with h5py.File(file_name, 'w') as f:\n",
    "        hdf5_io.save_to_hdf5(f, result)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "from tenpy.networks.site import SpinSite"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "site = SpinSite(S=1.5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'-1.5': 0, '-0.5': 1, '0.5': 2, '1.5': 3, 'down': 0, 'up': 3}"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "site.state_labels"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "531441"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "3**12"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  4 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "."
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      ".........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.2362683458304451, S = nan, norm_err = 1.3e-05\n",
      "Current memory usage 675.7 MB, time elapsed: 6.7 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 9.8701e-03, max_E_trunc = 2.8746e-02\n",
      "MPS bond dimensions: [4, 4, 4, 4, 4, 4, 4, 4, 3, 2, 4, 4, 4, 4, 4, 4]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.2362688851752530, S = nan, norm_err = 2.0e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 13.4 s\n",
      "Delta E = -5.3934e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 9.8509e-03, max_E_trunc = 2.8741e-02\n",
      "MPS bond dimensions: [4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.2362688619944464, S = 0.3449626971240179, norm_err = 2.2e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 17.0 s\n",
      "Delta E = 2.3181e-09, Delta S = nan (per sweep)\n",
      "max_trunc_err = 4.9235e-03, max_E_trunc = 2.8739e-02\n",
      "MPS bond dimensions: [3, 3, 4, 4, 4, 4, 4, 2, 2, 2, 4, 4, 4, 4, 4, 3]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.2362688619942617, S = 0.3449626970929214, norm_err = 2.4e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 20.7 s\n",
      "Delta E = 1.8474e-14, Delta S = -3.1097e-12 (per sweep)\n",
      "max_trunc_err = 4.9235e-03, max_E_trunc = 2.8739e-02\n",
      "MPS bond dimensions: [3, 3, 4, 4, 4, 4, 4, 2, 2, 2, 4, 4, 4, 4, 4, 3]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 2256, maximum chi = 4\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [3, 3, 4, 4, 4, 4, 4, 2, 2, 2, 4, 4, 4, 4, 4, 3]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000009+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [3, 3, 4, 4, 4, 4, 4, 2, 2, 2, 4, 4, 4, 4, 4, 3]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=4 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  5 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['P_tol'] for config lanczos_params\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.2508656464152956, S = nan, norm_err = 1.3e-05\n",
      "Current memory usage 675.7 MB, time elapsed: 20.2 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.9557e-03, max_E_trunc = 7.3272e-03\n",
      "MPS bond dimensions: [5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.2508738247962157, S = nan, norm_err = 2.1e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 28.4 s\n",
      "Delta E = -8.1784e-07, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.9158e-03, max_E_trunc = 7.2772e-03\n",
      "MPS bond dimensions: [5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.2508662137399256, S = 0.4934437841466527, norm_err = 1.8e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 34.2 s\n",
      "Delta E = 7.6111e-07, Delta S = nan (per sweep)\n",
      "max_trunc_err = 9.6838e-04, max_E_trunc = 7.3191e-03\n",
      "MPS bond dimensions: [4, 4, 5, 5, 5, 5, 5, 4, 4, 4, 5, 5, 5, 5, 5, 4]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.2508662137224462, S = 0.4934437681072582, norm_err = 1.8e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 39.8 s\n",
      "Delta E = 1.7479e-12, Delta S = -1.6039e-09 (per sweep)\n",
      "max_trunc_err = 9.6838e-04, max_E_trunc = 7.3191e-03\n",
      "MPS bond dimensions: [4, 4, 5, 5, 5, 5, 5, 4, 4, 4, 5, 5, 5, 5, 5, 4]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 2256, maximum chi = 5\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [4, 4, 5, 5, 5, 5, 5, 4, 4, 4, 5, 5, 5, 5, 5, 4]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000009+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [4, 4, 5, 5, 5, 5, 5, 4, 4, 4, 5, 5, 5, 5, 5, 4]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=5 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  8 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.2829822388690921, S = nan, norm_err = 1.1e-05\n",
      "Current memory usage 675.7 MB, time elapsed: 23.0 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 8.0677e-03, max_E_trunc = 2.7946e-02\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.2830087512431660, S = nan, norm_err = 1.7e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 32.1 s\n",
      "Delta E = -2.6512e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 8.0103e-03, max_E_trunc = 2.7813e-02\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.2830091021933612, S = nan, norm_err = 2.9e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 41.1 s\n",
      "Delta E = -3.5095e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 8.0095e-03, max_E_trunc = 2.7813e-02\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.2830091104267041, S = 0.6461973929230547, norm_err = 2.5e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 47.3 s\n",
      "Delta E = -8.2333e-10, Delta S = nan (per sweep)\n",
      "max_trunc_err = 4.0047e-03, max_E_trunc = 2.7813e-02\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.2830093722652123, S = 0.6461523728552305, norm_err = 2.8e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 53.5 s\n",
      "Delta E = -2.6184e-08, Delta S = -4.5020e-06 (per sweep)\n",
      "max_trunc_err = 4.0046e-03, max_E_trunc = 2.7813e-02\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "================================================================================\n",
      "DMRG finished after 50 sweeps.\n",
      "total size = 2736, maximum chi = 8\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000009+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=8 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  11 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.2902971133005749, S = nan, norm_err = 1.3e-05\n",
      "Current memory usage 675.7 MB, time elapsed: 25.1 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.0492e-03, max_E_trunc = 6.7538e-03\n",
      "MPS bond dimensions: [11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.2903051912253147, S = nan, norm_err = 2.0e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 33.4 s\n",
      "Delta E = -8.0779e-07, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.0153e-03, max_E_trunc = 6.7081e-03\n",
      "MPS bond dimensions: [11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.2903146805557952, S = 0.7421868511576812, norm_err = 5.0e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 38.7 s\n",
      "Delta E = -9.4893e-07, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.0093e-03, max_E_trunc = 6.7196e-03\n",
      "MPS bond dimensions: [11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.2903058518858899, S = 0.7418958057750772, norm_err = 5.5e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 44.8 s\n",
      "Delta E = 8.8287e-07, Delta S = -2.9105e-05 (per sweep)\n",
      "max_trunc_err = 1.0178e-03, max_E_trunc = 6.7456e-03\n",
      "MPS bond dimensions: [11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 2256, maximum chi = 11\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000004+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=11 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  16 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3083494240340698, S = nan, norm_err = 8.3e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 32.6 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.5674e-03, max_E_trunc = 1.0476e-02\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3083645960874293, S = nan, norm_err = 1.3e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 43.9 s\n",
      "Delta E = -1.5172e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.5526e-03, max_E_trunc = 1.0444e-02\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3097190335347619, S = nan, norm_err = 2.1e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 55.0 s\n",
      "Delta E = -1.3544e-04, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 3.4364e-03, max_E_trunc = 1.2461e-02\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3103355304198487, S = nan, norm_err = 3.4e-08\n",
      "Current memory usage 675.7 MB, time elapsed: 66.3 s\n",
      "Delta E = -6.1650e-05, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 3.5526e-03, max_E_trunc = 1.3269e-02\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "disable mixer after 50 sweeps, final amplitude 1.10e-08\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3103549512800328, S = 0.9021315491875989, norm_err = 7.3e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 75.9 s\n",
      "Delta E = -1.9421e-06, Delta S = nan (per sweep)\n",
      "max_trunc_err = 3.5524e-03, max_E_trunc = 1.3291e-02\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 3216\n",
      "Energy = -1.3103576674399164, S = 0.9021289839635425, norm_err = 6.9e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 81.3 s\n",
      "Delta E = -2.7162e-07, Delta S = -2.5652e-07 (per sweep)\n",
      "max_trunc_err = 1.7759e-03, max_E_trunc = 1.3291e-02\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 3216, maximum chi = 16\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=16 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  22 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3174587311162185, S = nan, norm_err = 8.8e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 33.7 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.1419e-03, max_E_trunc = 5.4696e-03\n",
      "MPS bond dimensions: [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3174771840614738, S = nan, norm_err = 1.4e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 45.2 s\n",
      "Delta E = -1.8453e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.1154e-03, max_E_trunc = 5.4043e-03\n",
      "MPS bond dimensions: [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3174775771734772, S = nan, norm_err = 2.2e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 56.8 s\n",
      "Delta E = -3.9311e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.1153e-03, max_E_trunc = 5.4043e-03\n",
      "MPS bond dimensions: [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3174775776320189, S = 1.0519984591335703, norm_err = 8.5e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 62.3 s\n",
      "Delta E = -4.5854e-11, Delta S = nan (per sweep)\n",
      "max_trunc_err = 5.5763e-04, max_E_trunc = 5.4043e-03\n",
      "MPS bond dimensions: [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3174728205868718, S = 1.0534201413646775, norm_err = 8.7e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 68.2 s\n",
      "Delta E = 4.7570e-07, Delta S = 1.4217e-04 (per sweep)\n",
      "max_trunc_err = 5.6441e-04, max_E_trunc = 5.4307e-03\n",
      "MPS bond dimensions: [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 3216\n",
      "Energy = -1.3174728361417749, S = 1.0534202232396073, norm_err = 8.8e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 74.2 s\n",
      "Delta E = -1.5555e-09, Delta S = 8.1875e-09 (per sweep)\n",
      "max_trunc_err = 5.6441e-04, max_E_trunc = 5.4308e-03\n",
      "MPS bond dimensions: [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 3216, maximum chi = 22\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000036+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=22 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  32 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3226447362726717, S = nan, norm_err = 1.0e-05\n",
      "Current memory usage 675.7 MB, time elapsed: 36.3 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 9.7864e-04, max_E_trunc = 4.4215e-03\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3226625247753034, S = nan, norm_err = 1.4e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 48.4 s\n",
      "Delta E = -1.7789e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 9.5690e-04, max_E_trunc = 4.3719e-03\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3226628286573572, S = nan, norm_err = 2.3e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 60.4 s\n",
      "Delta E = -3.0388e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 9.5652e-04, max_E_trunc = 4.3712e-03\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3226628380611061, S = 1.1778733964099817, norm_err = 7.0e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 68.3 s\n",
      "Delta E = -9.4037e-10, Delta S = nan (per sweep)\n",
      "max_trunc_err = 4.7822e-04, max_E_trunc = 4.3711e-03\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3226553288658920, S = 1.1792159622888492, norm_err = 8.7e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 76.5 s\n",
      "Delta E = 7.5092e-07, Delta S = 1.3426e-04 (per sweep)\n",
      "max_trunc_err = 4.8721e-04, max_E_trunc = 4.4445e-03\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 3216\n",
      "Energy = -1.3226558658129477, S = 1.1790059539783466, norm_err = 8.3e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 84.8 s\n",
      "Delta E = -5.3695e-08, Delta S = -2.1001e-05 (per sweep)\n",
      "max_trunc_err = 4.8598e-04, max_E_trunc = 4.4335e-03\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 3216, maximum chi = 32\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (0.9999999999999964+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=32 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  45 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3292338029937341, S = nan, norm_err = 9.5e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 39.9 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 7.1681e-04, max_E_trunc = 3.4672e-03\n",
      "MPS bond dimensions: [45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3295451476625146, S = nan, norm_err = 1.4e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 53.8 s\n",
      "Delta E = -3.1134e-05, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 8.6849e-04, max_E_trunc = 4.0896e-03\n",
      "MPS bond dimensions: [45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3295464510176913, S = nan, norm_err = 2.3e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 67.3 s\n",
      "Delta E = -1.3034e-07, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 8.5158e-04, max_E_trunc = 4.0204e-03\n",
      "MPS bond dimensions: [45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3295464571494335, S = 1.3162578375292542, norm_err = 9.3e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 76.1 s\n",
      "Delta E = -6.1317e-10, Delta S = nan (per sweep)\n",
      "max_trunc_err = 4.2578e-04, max_E_trunc = 4.0204e-03\n",
      "MPS bond dimensions: [45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3295451626137691, S = 1.3163463297138984, norm_err = 8.8e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 85.1 s\n",
      "Delta E = 1.2945e-07, Delta S = 8.8492e-06 (per sweep)\n",
      "max_trunc_err = 4.3411e-04, max_E_trunc = 4.0892e-03\n",
      "MPS bond dimensions: [45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45]\n",
      "================================================================================\n",
      "DMRG finished after 50 sweeps.\n",
      "total size = 2736, maximum chi = 45\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000056+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=45 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  64 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3325232634072179, S = nan, norm_err = 9.3e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 47.3 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 5.5203e-04, max_E_trunc = 2.5776e-03\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3325385116995747, S = nan, norm_err = 1.3e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 64.2 s\n",
      "Delta E = -1.5248e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 5.4081e-04, max_E_trunc = 2.5306e-03\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3325386650295457, S = nan, norm_err = 2.0e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 81.3 s\n",
      "Delta E = -1.5333e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 5.4072e-04, max_E_trunc = 2.5306e-03\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3325386687315159, S = 1.4118557925175028, norm_err = 2.5e-14\n",
      "Current memory usage 675.7 MB, time elapsed: 92.7 s\n",
      "Delta E = -3.7020e-10, Delta S = nan (per sweep)\n",
      "max_trunc_err = 2.7035e-04, max_E_trunc = 2.5306e-03\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3325369554444535, S = 1.4125064360581987, norm_err = 9.2e-15\n",
      "Current memory usage 675.7 MB, time elapsed: 104.7 s\n",
      "Delta E = 1.7133e-07, Delta S = 6.5064e-05 (per sweep)\n",
      "max_trunc_err = 2.7316e-04, max_E_trunc = 2.5439e-03\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "================================================================================\n",
      "DMRG finished after 50 sweeps.\n",
      "total size = 2736, maximum chi = 64\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000064+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=64 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  90 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3342603268911404, S = nan, norm_err = 6.4e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 65.2 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 4.7337e-04, max_E_trunc = 2.2313e-03\n",
      "MPS bond dimensions: [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3343165091954532, S = nan, norm_err = 9.8e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 88.3 s\n",
      "Delta E = -5.6182e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 4.4988e-04, max_E_trunc = 2.1409e-03\n",
      "MPS bond dimensions: [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3343213051203975, S = nan, norm_err = 1.3e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 111.4 s\n",
      "Delta E = -4.7959e-07, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 4.4885e-04, max_E_trunc = 2.1378e-03\n",
      "MPS bond dimensions: [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3343195561181460, S = 1.4932647740265237, norm_err = 2.0e-13\n",
      "Current memory usage 675.7 MB, time elapsed: 127.1 s\n",
      "Delta E = 1.7490e-07, Delta S = nan (per sweep)\n",
      "max_trunc_err = 2.2332e-04, max_E_trunc = 2.1248e-03\n",
      "MPS bond dimensions: [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3343187730433783, S = 1.4947269576888713, norm_err = 2.8e-13\n",
      "Current memory usage 675.7 MB, time elapsed: 144.4 s\n",
      "Delta E = 7.8307e-08, Delta S = 1.4622e-04 (per sweep)\n",
      "max_trunc_err = 2.2658e-04, max_E_trunc = 2.1491e-03\n",
      "MPS bond dimensions: [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 3216\n",
      "Energy = -1.3343188367896062, S = 1.4947344372235065, norm_err = 2.8e-13\n",
      "Current memory usage 675.7 MB, time elapsed: 161.6 s\n",
      "Delta E = -6.3746e-09, Delta S = 7.4795e-07 (per sweep)\n",
      "max_trunc_err = 2.2661e-04, max_E_trunc = 2.1493e-03\n",
      "MPS bond dimensions: [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 3216, maximum chi = 90\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000024+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=90 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  128 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3340084889375348, S = nan, norm_err = 4.4e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 96.2 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 3.1178e-04, max_E_trunc = 1.3974e-03\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3353294686126773, S = nan, norm_err = 5.8e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 134.1 s\n",
      "Delta E = -1.3210e-04, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.8959e-04, max_E_trunc = 1.3110e-03\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3353282690177224, S = nan, norm_err = 9.3e-08\n",
      "Current memory usage 675.7 MB, time elapsed: 171.5 s\n",
      "Delta E = 1.1996e-07, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.9335e-04, max_E_trunc = 1.3341e-03\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3353284174756794, S = 1.5615444134770424, norm_err = 3.2e-12\n",
      "Current memory usage 675.7 MB, time elapsed: 197.0 s\n",
      "Delta E = -1.4846e-08, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.4664e-04, max_E_trunc = 1.3340e-03\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3353231947791357, S = 1.5623254548506962, norm_err = 5.2e-12\n",
      "Current memory usage 675.7 MB, time elapsed: 225.4 s\n",
      "Delta E = 5.2227e-07, Delta S = 7.8104e-05 (per sweep)\n",
      "max_trunc_err = 1.4748e-04, max_E_trunc = 1.3458e-03\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "================================================================================\n",
      "DMRG finished after 50 sweeps.\n",
      "total size = 2736, maximum chi = 128\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.000000000000007+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=128 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  181 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3345801324892079, S = nan, norm_err = 5.8e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 162.6 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.0889e-04, max_E_trunc = 9.5465e-04\n",
      "MPS bond dimensions: [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3360220841168342, S = nan, norm_err = 6.5e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 227.6 s\n",
      "Delta E = -1.4420e-04, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.0471e-04, max_E_trunc = 9.3596e-04\n",
      "MPS bond dimensions: [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3360221250425326, S = nan, norm_err = 9.9e-08\n",
      "Current memory usage 675.7 MB, time elapsed: 289.8 s\n",
      "Delta E = -4.0926e-09, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.0457e-04, max_E_trunc = 9.3551e-04\n",
      "MPS bond dimensions: [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3360221246160364, S = 1.6212780398863258, norm_err = 1.5e-11\n",
      "Current memory usage 675.7 MB, time elapsed: 332.2 s\n",
      "Delta E = 4.2650e-11, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.0229e-04, max_E_trunc = 9.3554e-04\n",
      "MPS bond dimensions: [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3360192645171480, S = 1.6228814061627057, norm_err = 2.3e-11\n",
      "Current memory usage 675.7 MB, time elapsed: 382.8 s\n",
      "Delta E = 2.8601e-07, Delta S = 1.6034e-04 (per sweep)\n",
      "max_trunc_err = 1.0424e-04, max_E_trunc = 9.5035e-04\n",
      "MPS bond dimensions: [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 3216\n",
      "Energy = -1.3360192828608035, S = 1.6228822599695754, norm_err = 2.3e-11\n",
      "Current memory usage 675.7 MB, time elapsed: 433.1 s\n",
      "Delta E = -1.8344e-09, Delta S = 8.5381e-08 (per sweep)\n",
      "max_trunc_err = 1.0424e-04, max_E_trunc = 9.5036e-04\n",
      "MPS bond dimensions: [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 3216, maximum chi = 181\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000038+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181, 181]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=181 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  256 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3364282736622499, S = nan, norm_err = 4.0e-06\n",
      "Current memory usage 675.7 MB, time elapsed: 326.4 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.4118e-04, max_E_trunc = 6.6904e-04\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3364659374094288, S = nan, norm_err = 6.4e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 456.3 s\n",
      "Delta E = -3.7664e-06, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.4093e-04, max_E_trunc = 6.6976e-04\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3364659652819597, S = nan, norm_err = 1.0e-07\n",
      "Current memory usage 675.7 MB, time elapsed: 585.9 s\n",
      "Delta E = -2.7873e-09, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.4127e-04, max_E_trunc = 6.7130e-04\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3364659535332635, S = 1.6702839474417737, norm_err = 4.5e-11\n",
      "Current memory usage 675.7 MB, time elapsed: 673.4 s\n",
      "Delta E = 1.1749e-09, Delta S = nan (per sweep)\n",
      "max_trunc_err = 7.0666e-05, max_E_trunc = 6.7161e-04\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3364636033303441, S = 1.6717932741628503, norm_err = 7.2e-11\n",
      "Current memory usage 675.7 MB, time elapsed: 776.7 s\n",
      "Delta E = 2.3502e-07, Delta S = 1.5093e-04 (per sweep)\n",
      "max_trunc_err = 7.1992e-05, max_E_trunc = 6.8251e-04\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 3216\n",
      "Energy = -1.3364636039450488, S = 1.6717937944086079, norm_err = 7.2e-11\n",
      "Current memory usage 675.7 MB, time elapsed: 882.0 s\n",
      "Delta E = -6.1470e-11, Delta S = 5.2025e-08 (per sweep)\n",
      "max_trunc_err = 7.1995e-05, max_E_trunc = 6.8254e-04\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 3216, maximum chi = 256\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000009+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=256 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=2 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=1.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=1.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  362 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 816\n",
      "Energy = -1.3338751336590420, S = nan, norm_err = 3.8e-06\n",
      "Current memory usage 846.0 MB, time elapsed: 657.8 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.0769e-04, max_E_trunc = 4.9799e-04\n",
      "MPS bond dimensions: [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 1296\n",
      "Energy = -1.3367606148174005, S = nan, norm_err = 5.4e-07\n",
      "Current memory usage 846.0 MB, time elapsed: 957.6 s\n",
      "Delta E = -2.8855e-04, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.0370e-04, max_E_trunc = 4.8490e-04\n",
      "MPS bond dimensions: [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 1776\n",
      "Energy = -1.3367608749149724, S = nan, norm_err = 8.9e-08\n",
      "Current memory usage 846.0 MB, time elapsed: 1237.8 s\n",
      "Delta E = -2.6010e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.0298e-04, max_E_trunc = 4.8337e-04\n",
      "MPS bond dimensions: [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 2256\n",
      "Energy = -1.3367611670706481, S = 1.7125694921817129, norm_err = 1.1e-10\n",
      "Current memory usage 846.0 MB, time elapsed: 1401.1 s\n",
      "Delta E = -2.9216e-08, Delta S = nan (per sweep)\n",
      "max_trunc_err = 5.1503e-05, max_E_trunc = 4.8358e-04\n",
      "MPS bond dimensions: [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 2736\n",
      "Energy = -1.3367581932288033, S = 1.7140611455091763, norm_err = 8.9e-11\n",
      "Current memory usage 846.0 MB, time elapsed: 1672.1 s\n",
      "Delta E = 2.9738e-07, Delta S = 1.4917e-04 (per sweep)\n",
      "max_trunc_err = 5.4992e-05, max_E_trunc = 5.1216e-04\n",
      "MPS bond dimensions: [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 3216\n",
      "Energy = -1.3367582293241469, S = 1.7141503792149260, norm_err = 1.1e-10\n",
      "Current memory usage 846.0 MB, time elapsed: 1938.5 s\n",
      "Delta E = -3.6095e-09, Delta S = 8.9234e-06 (per sweep)\n",
      "max_trunc_err = 5.5026e-05, max_E_trunc = 5.1242e-04\n",
      "MPS bond dimensions: [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 3216, maximum chi = 362\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.00000000000001+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362, 362]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=362 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    }
   ],
   "source": [
    "Jx = 1.0\n",
    "Jy = 1.0\n",
    "Jz = 1.0\n",
    "\n",
    "S = 1.5\n",
    "\n",
    "Lx = 2\n",
    "Ly = 4\n",
    "shift = 0\n",
    "bc = 'periodic'\n",
    "bc_MPS = 'infinite'\n",
    "\n",
    "chi_list = [int(2**i) for i in np.arange(2, 9, 0.5)]\n",
    "\n",
    "result_dict = {}\n",
    "\n",
    "initial_psi = None\n",
    "for chi in chi_list:\n",
    "    param = (Jx, Jy, Jz, Lx, Ly, chi)\n",
    "    result = run_atomic(\n",
    "        Jx=Jx, Jy=Jy, Jz=Jz,\n",
    "        Lx=Lx, Ly=Ly, \n",
    "        S=S,\n",
    "        shift=shift,\n",
    "        chi=chi,\n",
    "        bc=bc, bc_MPS=bc_MPS,\n",
    "        initial_psi=initial_psi,\n",
    "    )\n",
    "    initial_psi = result['psi'].copy()\n",
    "    result_dict[param] = result\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "S_list = []\n",
    "xi_list = []\n",
    "for key, result in result_dict.items():\n",
    "    psi = result['psi']\n",
    "    S_list.append(np.max(psi.entanglement_entropy()))\n",
    "    xi_list.append(psi.correlation_length())\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEaCAYAAAAR0SDgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAZfklEQVR4nO3df7wddX3n8dfbcNUg1FiTqglguiumSleNjWL9iWs1/jZaXX+C2lqsvypdmiIo6q5bW031oRaVR6oYtRZtJcbfjdRVkHVFww8NEKMoiiQgQYyI3rokfvaPmYsnN/fXCTnn3HvP6/l4nEfOmfnOnM+ZMznvme/MnUlVIUnS7QZdgCRpdjAQJEmAgSBJahkIkiTAQJAktQwESRJgIKgHkjw/yRf6/J6nJXlfP9+zn5K8KMkFU4z/cpKX9LmmRyTZPsX45UkqySH9rEsHzkCYY5L8IMkf3cZ59PTHo6o+UlWP63a6JEckOSfJDUl+lmRrkhfN8D3fXFUH/JmSrE+yPcmvp3vPJBuS/K+O18ckuTbJyRO07fsPdb9U1VeqasXY69u6bo5frr2Q5FFtSPX0feYqA0H7GeAW3YeBHwH3BO4KnAD8uE/v/U3g5cDF3UyU5AHAl4C/qaq39aCunhqmrfckI8A7gQsHXctsZSDMYUnuleS8dmv6hiQf6xj30CTfaMd9I8lD2+F/AzwCOCPJzUnOaIdXklck+S7w3XbYnyW5MsmNST6VZGnH/CvJXyT5fvve65Lcrh23T/dGuwV9bjufHyc5bZKP9CBgQ1X9oqr2VNUlVfX5dh5j3Q8nJtk5fos8yRuT/NO4ti9McnVb32unWpZV9e6q+iLwH10s/wcD/w6cVlVnzHS6jun/Ncl17Xd0fpJjOsbdtV3mNyX5OvCfx0372CTfbqc9A8i48X+SZFuSnybZnOSeHeP2+67HTfvBsWWbZFnb/uXt63u132OSHJfkmnb4h4GjgE+369Vfd8zy+TP9Hiao5bNJXjVu2LeSrOlmPq2TgS8A3z6AaYeCgTC3vYlmBb8LcATwDwBJfhv4LPAumi3ttwOfTXLXqnot8BXglVV1WFW9smN+a4Bjgfsm+a/A3wL/DbgH8EPgo+Pe/+nAKuCBwNOAPxlfYJLDaX40/w1YCtwL+OIkn+drwLuTPCfJUZO0eTRwNPA44DXTdFE8HFgBPAZ4fZL7TNG2Ww+m+Ux/WVUHeuzi8zSf5Xdo9kw+0jHu3TThdA+a5Xrrsk2yGDgHeB2wGPge8LCO8WuA04BnAEtovu+zx733GtrveoK6zgOOa58/Cvh++y/AI4Gv1Lhr3lTV8cDVwFPa9eqtHaNvy/fwQeAFHZ/t/sAy4HPt628l2T3J4z0d092TZhn+zy7ee+gYCHPbLTTdK0ur6j+qamyr/EnAd6vqw+2W9tk0W0VPmWZ+f1tVN1bVKPB84KyquriqfgWcCvxhkuUd7d/Str8aeAfw3Anm+WTguqp6W1vjz6tqsl32Z9H8eJ0OXJXk0iQPGtfmf7R7EFuBD0zynp1tR6vqmzRdQvef8tN35yHAz2h+1A9IVZ3VLo9fAW8E7p/kzkkWAH8MvL79rJfR/DCOeSJwRVV9vKpuoVn213WMfynNd7mtqvYAbwYe0LmXwL7f9XjnAY9o9/geCbyV3wTOo9rx3bgt38MngaOTHN2+Ph74WFX9P4Cqul9VLZrk8fKO+bwLOL2qbu6y9qFiIMxtf03TVfD1JJcnGduKXEqzRd/phzRbVlP5UcfzfebR/kf6ybh5dLb/YTvNeEfSbMFOq6p+WlWvqapjgLsBlwKbknR2h8zkPcd0/kj+EjgMoO3SGHtMticynXcD3wDOTXKXdr5ndsx3sm4x2rYLkvxdku8luQn4QTtqMc1W/SHs/1nHLO0c126td7a9J/DOsS1l4Eaa9WSy724fVfU94GbgATTdi58BdiZZwYEFwoTfw0y0YfkvwAvagHouzbGmGUvyFODwqvrYtI2HnIEwh1XVdVX1Z1W1lGar8D1J7gXspPlR6HQUsGNs0slm2fF8n3kkuRNN99OOjjZHjpv/zgnm+SPG9X/PRFXdAPw9zY/fb3f5ntPN+7COx9XdTt/aS7MXdTWwOclvVdWfd8z3zdNM/zyabrY/Au4MLG+HB9gF7GH/zzrm2s5xbWB2tv0R8NJxW8sLq+qrHW2mu8zxecAzgdtX1Y729Qk03ZOXTjJNry6d/EGaZf0Y4JdV9X/HRrQbQjdP8jizbfYYYFV7vOY64NnASUk+2aN65ywDYQ5L8qwkR7Qvf0rzH3IvTf/qvZM8L8khSZ5N01f8mbbtj4H/NM3s/xl4cZIHJLkDTbfDhVX1g442a5PcJcmRwKuBibbAPgPcPclJSe6Q5PAkx07yed6S5Pfbmg8HXgZcWVU/6Wh2epJD0xyAffEk79m1JLdPckeaH+SRJHdst0gn1XbXPAu4AfhcG5oTOaSd39hjBDgc+BXNXtehNMt3bL57gY3AG9vPel/ghR3z+yxwTJJnpDlL6C+Au3eMPxM4tV1GtN1Qz5rpsmidB7wSOL99/WXgVcAFbX0Tmcl6NZ0F45bV7dsA+DXwNsbtHVTVMeMCvvPx522z04F70+zxPAD4FPCPNOuPOhgIc9uDgAuT3Eyzkr+6qq5qf0CfTHNWxU9oupae3G51Q3Pq3TPTnIHyrolm3J5xczrNwctrabbynzOu2SeBi2i2GD8LvH+C+fwceCzN8YvraM5qefQkn+dQ4BPAbpoDmfcEnjquzXnAlTQHpv++qg7WH8B9ARgFHgqsb58/crqJ2r7sZ9AcAP50koUTNHtvO7+xxweAD9F0A+0ArqA5oN7plTRdK9cBG9ppxt7zBpog+jua7/do4P90jP8E8Bbgo2131GXAE6b7LOOcRxNaY4FwAc33c/6kUzQnIbyu7ar6qy7fb8xr2HdZ/e92+IeA/wL8U7czbI/TXDf2aOf7i6q68QBrnLfiDXJ0IJIUcHRVXdmn91sOXAWMtAdKNUSSnACcWFUPH3Qt85l7CJJmtSSH0vzR4PpB1zLfGQiSZq0kq2kOsv+Y5riWesguI0kS4B6CJKk1Zy9stXjx4lq+fPmgy5CkOeWiiy66oaqWTDRuzgbC8uXL2bJly6DLkKQ5Jcn4qxjcyi4jSRJgIEiSWgaCJAkwECRJLQNBkgTM4bOMJGnYbLpkB+s2b2fn7lGWLlrI2tUrWLNyutuczJyBIElzwKZLdnDqxq2M3tJcfXzH7lFO3bgV4KCFgl1GkjQHrNu8/dYwGDN6y17Wbd5+0N7DQJCkOWDn7olufz358ANhIEjSHLB00UT3Xpp8+IEwECRpDli7egULRxbsM2zhyALWrl5x0N7Dg8qS1KVen+0zkbH5z+mzjNobsH+I5ibgvwbWV9U7J2n7IJp7yz67qj7e69okqVv9ONtnMmtWLuvpe/Sjy2gPcHJV3Qd4CPCKJPcd3yjJApobg2/uQ02SdED6cbbPoPQ8EKrq2qq6uH3+c2AbMFHEvQo4B7i+1zVJ0oHqx9k+g9LXg8pJlgMrgQvHDV8GPB04c5rpT0yyJcmWXbt29axOSZpMP872GZS+BUKSw2j2AE6qqpvGjX4HcEpV7d1vwg5Vtb6qVlXVqiVLJrzhjyT1VD/O9hmUvpxllGSEJgw+UlUbJ2iyCvhoEoDFwBOT7KmqTf2oT5Jmqh9n+wxKP84yCvB+YFtVvX2iNlX1ux3tNwCfMQwkzVa9PttnUPqxh/Aw4Hhga5JL22GnAUcBVNWUxw0kSf3R80CoqguAdNH+Rb2rRpI0GS9dIUkCDARJUstAkCQBBoIkqeXVTiXNGoO4iqh+w0CQNCsM8iqiathlJGlWmM9XEZ0rDARJs8J8voroXGEgSJoV5vNVROcKA0HSrDCfryI6V3hQWdKsMJ+vIjpXGAiSZo35ehXRucIuI0kSYCBIkloGgiQJMBAkSS0DQZIEGAiSpJaBIEkCDARJUstAkCQBBoIkqdXzQEhyZJIvJdmW5PIkr56gzfOTfKt9fDXJ/XtdlyRpX/24ltEe4OSqujjJ4cBFSc6tqis62lwFPKqqfprkCcB64Ng+1CZJavU8EKrqWuDa9vnPk2wDlgFXdLT5asckXwOO6HVdkqR99fUYQpLlwErgwima/Snw+UmmPzHJliRbdu3a1YMKJWl49S0QkhwGnAOcVFU3TdLm0TSBcMpE46tqfVWtqqpVS5Ys6V2xkjSE+nI/hCQjNGHwkaraOEmb+wHvA55QVT/pR12SpN/ox1lGAd4PbKuqt0/S5ihgI3B8VX2n1zVJkvbXjz2EhwHHA1uTXNoOOw04CqCqzgReD9wVeE+TH+ypqlV9qE2S1OrHWUYXAJmmzUuAl/S6Fmku23TJDu83rJ7ynsrSHLDpkh2cunEro7fsBWDH7lFO3bgVwFDQQeOlK6Q5YN3m7beGwZjRW/aybvP2AVWk+chAkOaAnbtHuxouHQgDQZoDli5a2NVw6UAYCNIcsHb1ChaOLNhn2MKRBaxdvWJAFWk+8qCyNAeMHTj2LCP1koEgzRFrVi4zANRTdhlJkgADQZLUMhAkSYCBIElqGQiSJMBAkCS1DARJEmAgSJJaBoIkCTAQJEktA0GSBBgIkqSWgSBJAgwESVKr54GQ5MgkX0qyLcnlSV49QZskeVeSK5N8K8kDe12XJGlf/bgfwh7g5Kq6OMnhwEVJzq2qKzraPAE4un0cC7y3/VeS1Cc930Ooqmur6uL2+c+BbcD4u3w8DfhQNb4GLEpyj17XJkn6jb4eQ0iyHFgJXDhu1DLgRx2vr2H/0CDJiUm2JNmya9euntUpScOob4GQ5DDgHOCkqrpp/OgJJqn9BlStr6pVVbVqyZIlvShTkoZWXwIhyQhNGHykqjZO0OQa4MiO10cAO/tRmySp0Y+zjAK8H9hWVW+fpNmngBPas40eAvysqq7tdW2SpN/ox1lGDwOOB7YmubQddhpwFEBVnQl8DngicCXwS+DFfahLktSh54FQVRcw8TGCzjYFvKLXtUiSJudfKkuSAANBktQyECRJgIEgSWoZCJIkwECQJLUMBEkSYCBIkloGgiQJMBAkSS0DQZIEGAiSpJaBIEkC+nP5a+mg2XTJDtZt3s7O3aMsXbSQtatXsGblfndblXQADATNGZsu2cGpG7cyesteAHbsHuXUjVsBDAXpILDLSHPGus3bbw2DMaO37GXd5u0DqkiaXwwEzRk7d492NVxSdwwEzRlLFy3sarik7hgImjPWrl7BwpEF+wxbOLKAtatXDKgiaX7xoLLmjLEDx55lJPWGgaA5Zc3KZQaA1CM97zJKclaS65NcNsn4Oyf5dJJvJrk8yYt7XZMkaX/9OIawAXj8FONfAVxRVfcHjgPeluT2fahLktSh54FQVecDN07VBDg8SYDD2rZ7el2XJGlf0wZCki8mOabj9VOTvC7JsQephjOA+wA7ga3Aq6vq15PUcmKSLUm27Nq16yC9vSQJZraHcERVXQ6Q5KHAh4GjgA8kefpBqGE1cCmwFHgAcEaS35qoYVWtr6pVVbVqyZIlB+GtJUljZhIIN3U8PwE4s6pOpOnvP+Ug1PBiYGM1rgSuAn7vIMxXktSFmQTClUmemeR3gDXAJwGq6nrgDgehhquBxwAkuRuwAvj+QZivJKkLM/k7hL+k6SY6Gzi3qr4KkGSE5iDwlJKcTbM3sTjJNcAbgBGAqjoTeBOwIclWIMApVXVD9x9FknRbTBsIVXUd8Ngktxt3sPfRwJdmMP1zpxm/E3jcdPORJPXWjP9SefyZP1X1BeALB70iSdJAeHE7SRJgIEiSWgaCJAkwECRJLQNBkgQYCJKkloEgSQIMBElSy0CQJAEGgiSpZSBIkgADQZLUMhAkSYCBIElqGQiSJMBAkCS1DARJEmAgSJJaBoIkCTAQJEktA0GSBPQhEJKcleT6JJdN0ea4JJcmuTzJeb2uSZK0v37sIWwAHj/ZyCSLgPcAT62qY4Bn9aEmSdI4PQ+EqjofuHGKJs8DNlbV1W3763tdkyRpf7PhGMK9gbsk+XKSi5KcMFnDJCcm2ZJky65du/pYoiTNf7MhEA4B/gB4ErAaOD3JvSdqWFXrq2pVVa1asmRJP2uUpHnvkEEXAFwD3FBVvwB+keR84P7AdwZbliQNl9mwh/BJ4BFJDklyKHAssG3ANUnS0On5HkKSs4HjgMVJrgHeAIwAVNWZVbUtyb8B3wJ+DbyvqiY9RVWS1Bs9D4Sqeu4M2qwD1vW6FknS5GZDl5EkaRYwECRJgIEgSWoZCJIkYHb8HYJ6bNMlO1i3eTs7d4+ydNFC1q5ewZqVywZdlqRZxkCY5zZdsoNTN25l9Ja9AOzYPcqpG7cCGAqS9mGX0Ty3bvP2W8NgzOgte1m3efuAKpI0WxkI89zO3aNdDZc0vAyEeW7pooVdDZc0vAyEeW7t6hUsHFmwz7CFIwtYu3rFgCqSNFt5UHmeGztw7FlGkqZjIAyBNSuXGQCSpmWXkSQJMBAkSS0DQZIEGAiSpJaBIEkCDARJUstAkCQBBoIkqWUgSJKAPgRCkrOSXJ/ksmnaPSjJ3iTP7HVNkqT99WMPYQPw+KkaJFkAvAXY3Id6JEkT6HkgVNX5wI3TNHsVcA5wfa/rkSRNbODHEJIsA54OnDmDticm2ZJky65du3pfnCQNkYEHAvAO4JSq2jtdw6paX1WrqmrVkiVLel+ZJA2R2XD561XAR5MALAaemGRPVW0aaFWSNGQGHghV9btjz5NsAD5jGEhS//U8EJKcDRwHLE5yDfAGYASgqqY9biBJ6o+eB0JVPbeLti/qYSmSpCnMhoPKkqRZwECQJAEGgiSpZSBIkgADQZLUMhAkSYCBIElqGQiSJMBAkCS1DARJEmAgSJJaBoIkCTAQJEktA0GSBMyCG+TMRZsu2cG6zdvZuXuUpYsWsnb1CtasXDbosiTpNjEQurTpkh2cunEro7c0t4DesXuUUzduBTAUJM1pdhl1ad3m7beGwZjRW/aybvP2AVUkSQeHgdClnbtHuxouSXOFgdClpYsWdjVckuYKA6FLa1evYOHIgn2GLRxZwNrVKwZUkSQdHB5U7tLYgWPPMpI03/Q8EJKcBTwZuL6qfn+C8c8HTmlf3gy8rKq+2eu6bos1K5cZAJLmnX50GW0AHj/F+KuAR1XV/YA3Aev7UJMkaZye7yFU1flJlk8x/qsdL78GHNHrmiRJ+5ttB5X/FPj8ZCOTnJhkS5Itu3bt6mNZkjT/zZpASPJomkA4ZbI2VbW+qlZV1aolS5b0rzhJGgKz4iyjJPcD3gc8oap+Muh6JGkYDXwPIclRwEbg+Kr6zqDrkaRh1Y/TTs8GjgMWJ7kGeAMwAlBVZwKvB+4KvCcJwJ6qWtXruiRJ++rHWUbPnWb8S4CX9LoOSdLUZsUxhH7xPgaSNLmhCQTvYyBJUxv4QeV+8T4GkjS1oQkE72MgSVMbmkDwPgaSNLWhCQTvYyBJUxuag8rex0CSpjY0gQDex0CSpjI0XUaSpKkZCJIkwECQJLUMBEkSYCBIklqpqkHXcECS7AJ+OOg6DpLFwA2DLmKWctlMzmUzNZfPxO5ZVRPecnLOBsJ8kmSL94CYmMtmci6bqbl8umeXkSQJMBAkSS0DYXZYP+gCZjGXzeRcNlNz+XTJYwiSJMA9BElSy0CQJAEGQt8kOSvJ9Ukum2R8krwryZVJvpXkgf2ucVBmsGyOS/KzJJe2j9f3u8ZBSXJkki8l2Zbk8iSvnqDNUK47M1w2Q7vuHIihuvz1gG0AzgA+NMn4JwBHt49jgfe2/w6DDUy9bAC+UlVP7k85s8oe4OSqujjJ4cBFSc6tqis62gzrujOTZQPDu+50zT2EPqmq84Ebp2jyNOBD1fgasCjJPfpT3WDNYNkMraq6tqoubp//HNgGjL+px1CuOzNcNuqCgTB7LAN+1PH6Gly5O/1hkm8m+XySYwZdzCAkWQ6sBC4cN2ro150plg247syYXUazRyYY5jnBjYtprr9yc5InAptoukeGRpLDgHOAk6rqpvGjJ5hkaNadaZbN0K873XAPYfa4Bjiy4/URwM4B1TKrVNVNVXVz+/xzwEiSxQMuq2+SjND84H2kqjZO0GRo153pls2wrzvdMhBmj08BJ7RnjDwE+FlVXTvoomaDJHdPkvb5g2nW258Mtqr+aD/3+4FtVfX2SZoN5bozk2UzzOvOgbDLqE+SnA0cByxOcg3wBmAEoKrOBD4HPBG4Evgl8OLBVNp/M1g2zwRelmQPMAo8p4bnT+wfBhwPbE1yaTvsNOAoGPp1ZybLZpjXna556QpJEmCXkSSpZSBIkgADQZLUMhAkSYCBIElqGQiSJMBAkCS1DASpC0luvg3TPqq978PeJFclOflg1ibdVv6lstQ/dwP+FXhtVe0ecC3SftxDkA5Akv+e5LL2cVLH8NOTfDvJuUnOTvJXHZO9EPh34Gf9rleaCfcQpC4l+QOa6wUdS3Pp6QuTnAcsAP6Y5rr8h9Bcevmijkn/Afg4sDfJy6rqn/tauDQNA0Hq3sOBT1TVLwCSbAQeQbPH/cmqGm2Hf3psgiS/B7wVeArwZS+wptnIQJC6N9ENaaYaDvBS4O1V9aUe1CMdFB5DkLp3PrAmyaFJ7gQ8HfgKcAHwlCR3bO/i9aSOae4I3L3/pUoz5x6C1KWqujjJBuDr7aD3VdUlAEk+BXwT+CGwhd8cQF4H/GOSFwDXAicMw01sNLd4PwTpIEpyWHv/3kNp9iROrKqLx7X5IPAvVfXZgRQpTcI9BOngWp/kvjRdRB+cIAyeDNyJ5vRTaVZxD0GSBHhQWZLUMhAkSYCBIElqGQiSJMBAkCS1DARJEmAgSJJa/x9P3rM8CeDJlgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.log(xi_list)[1:], S_list[1:], 'o')\n",
    "plt.xlabel(r'$\\log \\xi$')\n",
    "plt.ylabel(r'$S$')\n",
    "plt.title(f'Isotropic Spin-{int(S)} K-Ladder with Ly={Ly}')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "9"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(result_dict)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}