Drop of quasiparticle weight by increasing interaction

The quasiparticle weight of the electronic system drops as the local interaction is increased. When studying the two band case, the half-filled scenario and other commensurate fillings undergoe a transition into the Mott insulator. All other doping configurations retain the strongly correlated metal.

Traceback (most recent call last):
  File "/home/docs/checkouts/readthedocs.org/user_builds/slaveparticles/envs/latest/local/lib/python2.7/site-packages/sphinx_gallery/gen_rst.py", line 467, in execute_script
    exec(code_block, example_globals)
  File "<string>", line 19, in <module>
  File "<string>", line 15, in plot_dop_phase
  File "build/bdist.linux-x86_64/egg/slaveparticles/utils/plotter.py", line 41, in calc_z
    hopping=[0.5]*2*bands, populations=[dop]*2*bands)
  File "build/bdist.linux-x86_64/egg/slaveparticles/spins.py", line 118, in __init__
    self.set_filling(self.param['populations'])
  File "build/bdist.linux-x86_64/egg/slaveparticles/spins.py", line 143, in set_filling
    spin_gen_op(self.oper, estimate_gauge(populations))
  File "build/bdist.linux-x86_64/egg/slaveparticles/spins.py", line 79, in spin_gen_op
    oper['Sfliphop'] = spinflipandhop(slaves)
  File "build/bdist.linux-x86_64/egg/slaveparticles/spins.py", line 53, in spinflipandhop
    sfh += Sup[2*n+1] * Sdw[2*n    ] * Sup[2*m    ] * Sdw[2*m+1]
  File "/usr/lib/python2.7/dist-packages/scipy/sparse/base.py", line 387, in __iadd__
    raise NotImplementedError
NotImplementedError

from __future__ import division, absolute_import, print_function
import slaveparticles.utils.plotter as ssplt
import numpy as np
import matplotlib.pyplot as plt

#band dop_phasediag
def plot_dop_phase(bands, int_max, hund_cu):
    """Phase plot of Quasiparticle weight for N degenerate bands
    under doping shows transition only at integer filling
    the rest are metallic states"""
    name = 'Z_dop_phase_'+str(bands)+'bands_U'+str(int_max)+'J'+str(hund_cu)
    dop = np.sort(np.hstack((np.linspace(0.01,0.99,50),
                    np.arange(1./2./bands, 1, 1/2/bands))))
    data = ssplt.calc_z(bands, dop, np.arange(0, int_max, 0.1), hund_cu, name)

    ssplt.surf_z(data, name)

plot_dop_phase(2, 6, 0.)

Total running time of the script: (0 minutes 0.000 seconds)

Download Python source code: plot_z_dop_2orb.py

Download IPython notebook: plot_z_dop_2orb.ipynb