{
  "nbformat_minor": 0, 
  "nbformat": 4, 
  "cells": [
    {
      "execution_count": null, 
      "cell_type": "code", 
      "source": [
        "%matplotlib inline"
      ], 
      "outputs": [], 
      "metadata": {
        "collapsed": false
      }
    }, 
    {
      "source": [
        "\nCreated on Wed Jul 16 22:37:48 2014\n\n@author: oscar\nplots of report\n"
      ], 
      "cell_type": "markdown", 
      "metadata": {}
    }, 
    {
      "execution_count": null, 
      "cell_type": "code", 
      "source": [
        "from __future__ import division, absolute_import, print_function\nimport slaveparticles.utils.plotter as ssplt\nimport numpy as np\nimport matplotlib.pyplot as plt\n\n#Degenerate bands\ndef plot_degbandshalffill():\n    \"\"\"Plot of Quasiparticle weight for degenerate\n       half-filled bands, showing the Mott transition\"\"\"\n    ulim = [3.45, 5.15, 6.85, 8.55]\n    bands = range(1, 5)\n    for band, u_int in zip(bands, ulim):\n        name = 'Z_half_'+str(band)+'band'\n        dop = [0.5]\n        data = ssplt.calc_z(band, dop, np.arange(0, u_int, 0.1),0., name)\n        plt.plot(data['u_int'], data['zeta'][0, :, 0], label='$N={}$'.format(str(band)))\n\n    ssplt.label_saves('Z_half_multiorb.png')\n\n#band dop\ndef plot_dop(bands, int_max, dop, hund_cu, name):\n    \"\"\"Plot of Quasiparticle weight for N degenerate bands\n       under selected doping shows transition only at half-fill\n       the rest are metallic states\"\"\"\n    data = ssplt.calc_z(bands, dop, np.arange(0, int_max, 0.1), hund_cu, name)\n    ssplt.plot_curves_z(data, name)\n\n#band dop_phasediag\ndef plot_dop_phase(bands, int_max, hund_cu):\n    \"\"\"Phase plot of Quasiparticle weight for N degenerate bands\n       under doping shows transition only at interger filling\n       the rest are metallic states\"\"\"\n    name = 'Z_dop_phase_'+str(bands)+'bands_U'+str(int_max)+'J'+str(hund_cu)\n    dop = np.sort(np.hstack((np.linspace(0.01,0.99,50),\n                    np.arange(1./2./bands, 1, 1/2/bands))))\n    data = ssplt.calc_z(bands, dop, np.arange(0, int_max, 0.1), hund_cu, name)\n\n    ssplt.imshow_z(data, name)\n    ssplt.surf_z(data, name)\n\nif __name__ == \"__main__\":\n    plot_degbandshalffill()\n    plot_dop(1, 4.6, [0.5, 0.499, 0.495, 0.49, 0.45, 0.4, 0.2, 0.1], 0., 'Z_dop_1orb')\n    plot_dop(3, 7, np.asarray([3, 4, 5])/6, 0., 'Z_3b_int_fil')\n    for U, N in zip([4.6, 6, 8], range(1, 4)):\n        plot_dop_phase(N, U, 0.)"
      ], 
      "outputs": [], 
      "metadata": {
        "collapsed": false
      }
    }
  ], 
  "metadata": {
    "kernelspec": {
      "display_name": "Python 2", 
      "name": "python2", 
      "language": "python"
    }, 
    "language_info": {
      "mimetype": "text/x-python", 
      "nbconvert_exporter": "python", 
      "name": "python", 
      "file_extension": ".py", 
      "version": "2.7.6", 
      "pygments_lexer": "ipython2", 
      "codemirror_mode": {
        "version": 2, 
        "name": "ipython"
      }
    }
  }
}