{
  "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# Hund Coupling as L. de'Medici PRB 83, 205112 (2011)\ndef hund_coup(bands, dop, u_lim):\n    hund_cu = [0, 0.05, 0.1, 0.15, 0.2, 0.25]\n    plt.figure()\n    for int_max, hund in zip(u_lim, hund_cu):\n        name = 'Z_'+str(bands)+'bands_U'+str(int_max)+'J'+str(hund)+'n'+str(dop)\n        data = ssplt.calc_z(2, [dop], np.arange(0, int_max, 0.1), hund, name)\n        plt.plot(data['u_int'], data['zeta'][0, :, 0], label='$J/U={}$'.format(str(hund)))\n    plt.title('Quasiparticle weight for {} electron(s) in {} bands'.format(str(dop*2*bands), str(bands)))\n    ssplt.label_saves('Z_{}bands_{}n_Hund.png'.format(str(bands), str(dop)))\n\nif __name__ == \"__main__\":\n    hund_coup(2, 0.5, [5.2, 3.8, 3.4, 3.1, 3.0, 3.0])\n    hund_coup(2, 0.25, [5.0, 5.5, 6.0, 7.5, 10.0, 15.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"
      }
    }
  }
}