Drop of quasiparticle weight by increasing interaction

The quasiparticle weight of the electronic system drops as the local interaction is increased. Multi orbital degenerate systems are studied

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 21, in <module>
  File "<string>", line 17, in plot_degbandshalffill
  File "/usr/lib/pymodules/python2.7/matplotlib/pyplot.py", line 2980, in plot
    ax = gca()
  File "/usr/lib/pymodules/python2.7/matplotlib/pyplot.py", line 803, in gca
    ax =  gcf().gca(**kwargs)
  File "/usr/lib/pymodules/python2.7/matplotlib/pyplot.py", line 450, in gcf
    return figure()
  File "/usr/lib/pymodules/python2.7/matplotlib/pyplot.py", line 423, in figure
    **kwargs)
  File "/usr/lib/pymodules/python2.7/matplotlib/backends/backend_tkagg.py", line 79, in new_figure_manager
    return new_figure_manager_given_figure(num, figure)
  File "/usr/lib/pymodules/python2.7/matplotlib/backends/backend_tkagg.py", line 87, in new_figure_manager_given_figure
    window = Tk.Tk()
  File "/usr/lib/python2.7/lib-tk/Tkinter.py", line 1767, in __init__
    self.tk = _tkinter.create(screenName, baseName, className, interactive, wantobjects, useTk, sync, use)
TclError: no display name and no $DISPLAY environment variable

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

#Degenerate bands
def plot_degbandshalffill():
    """Plot of Quasiparticle weight for degenerate
    half-filled bands, showing the Mott transition"""
    ulim = [3.45, 5.15, 6.85, 8.55]
    bands = range(1, 5)
    for band, u_int in zip(bands, ulim):
        name = 'Z_half_'+str(band)+'band'
        dop = [0.5]
        data = ssplt.calc_z(band, dop, np.arange(0, u_int, 0.1),0., name)
        plt.plot(data['u_int'], data['zeta'][0, :, 0], label='N={}'.format(str(band)))

    ssplt.label_saves('Z_half_multiorb.png')

plot_degbandshalffill()

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

Download Python source code: plot_z_half_multiorb.py

Download IPython notebook: plot_z_half_multiorb.ipynb