.. _sphx_glr_dimer_lattice_plot_phase_bandwidth.py:


===========================
Phase diagrams on Bandwidth
===========================





.. code-block:: python

    # Author: Óscar Nájera

    from __future__ import division, absolute_import, print_function
    import os
    from math import log, ceil
    import numpy as np
    import matplotlib.pyplot as plt
    import dmft.common as gf
    import dmft.dimer as dimer
    import dmft.ipt_imag as ipt


    def dos_at_fermi_level_fixbeta(xlist, ylist, beta, phasename, datapath):
        dos_fl = []
        w_n = gf.matsubara_freq(beta, 3)
        save_file = 'phase_dimer_ipt_{}.npy'.format(phasename)
        if os.path.exists(save_file):
            dos_fl = np.load(save_file)
            return dos_fl

        for xdat in xlist:
            filestr = datapath.format(phasename, xdat)
            gfs = np.load(filestr)
            dos_fl.append(np.array([gf.fit_gf(w_n, gfs[i][0][:3])(0.)
                                    for i in ylist]))

        dos_fl = -np.array(dos_fl).T
        np.save(save_file, dos_fl)

        return dos_fl


    def plot_phase_diagram_fixbeta(xlist, ylist, phasename, beta, ax, alpha=1):

        x, y = np.meshgrid(xlist, ylist)
        yrange = list(range(len(ylist)))
        if alpha < 1:
            yrange = yrange[::-1]

        ax.pcolormesh(x, y,
                      dos_at_fermi_level_fixbeta(xlist, yrange, beta, phasename,
                                                 'disk/phase_Dimer_ipt_{}/tp{:.3}/giw.npy'),
                      cmap=plt.get_cmap(r'viridis'), alpha=alpha, vmin=0, vmax=3)
        ax.axis([x.min(), x.max(), y.min(), y.max()])


    TPR = np.arange(0, 1.1, 0.02)
    Drange = np.linspace(0.05, 1.1, 81)

    for beta in [1000., 100., 30.]:
        f, ax = plt.subplots()
        plot_phase_diagram_fixbeta(
            TPR, Drange, 'D_met_B{}'.format(beta), beta, ax, 1)
        plot_phase_diagram_fixbeta(
            TPR, Drange, 'D_ins_B{}'.format(beta), beta, ax, .21)
        ax.set_xlabel(r'$t_\perp/U$')
        ax.set_ylabel(r'$D/U$')
        ax.set_title(
            'Phase diagram $\\beta={}$,\n color represents $-\\Im G_{{AA}}(0)$'.format(beta))
    plt.show()




.. rst-class:: sphx-glr-horizontal


    *

      .. image:: /dimer_lattice/images/sphx_glr_plot_phase_bandwidth_001.png
            :scale: 47

    *

      .. image:: /dimer_lattice/images/sphx_glr_plot_phase_bandwidth_002.png
            :scale: 47

    *

      .. image:: /dimer_lattice/images/sphx_glr_plot_phase_bandwidth_003.png
            :scale: 47





.. code-block:: python



    def dos_at_fermi_level_temp(xlist, temp, phasename, datapath):
        dos_fl = []
        save_file = 'phase_dimer_ipt_{}.npy'.format(phasename)
        if os.path.exists(save_file):
            dos_fl = np.load(save_file)
            return dos_fl

        for T in temp:
            w_n = gf.matsubara_freq(1 / T, 3)
            filestr = datapath.format(phasename, 1 / T)
            gfs = np.load(filestr)
            dos_fl.append(np.array([gf.fit_gf(w_n, gfs[i][0][:3])(0.)
                                    for i in xlist]))

        dos_fl = -np.array(dos_fl)
        np.save(save_file, dos_fl)

        return dos_fl


    def plot_phase_diagram_temp(xlist, temp, phasename, ax, alpha=1):

        x, y = np.meshgrid(xlist, temp)
        x_range = list(range(len(xlist)))
        if alpha < 1:
            x_range = x_range[::-1]

        ax.pcolormesh(x, y, dos_at_fermi_level_temp(x_range, temp, phasename,
                                                    'disk/phase_Dimer_ipt_{}/B{:.5}/giw.npy'),
                      cmap=plt.get_cmap(r'viridis'), alpha=alpha, vmin=0, vmax=3)
        ax.axis([x.min(), x.max(), y.min(), y.max()])

    TEMP = np.arange(1 / 512., .08, 1 / 400)
    for tp in [.15, .3, .5]:
        f, ax = plt.subplots()
        plot_phase_diagram_temp(Drange, TEMP, 'D_met_tp{}'.format(tp), ax, 1)
        plot_phase_diagram_temp(Drange, TEMP, 'D_ins_tp{}'.format(tp), ax, 0.2)
        ax.set_ylabel(r'$T/U$')
        ax.set_xlabel(r'$D/U$')
        ax.set_title(
            'Phase diagram $t_\\perp={}$,\n color represents $-\\Im G_{{AA}}(0)$'.format(tp))

    plt.show()



.. rst-class:: sphx-glr-horizontal


    *

      .. image:: /dimer_lattice/images/sphx_glr_plot_phase_bandwidth_004.png
            :scale: 47

    *

      .. image:: /dimer_lattice/images/sphx_glr_plot_phase_bandwidth_005.png
            :scale: 47

    *

      .. image:: /dimer_lattice/images/sphx_glr_plot_phase_bandwidth_006.png
            :scale: 47




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



.. container:: sphx-glr-footer


  .. container:: sphx-glr-download

     :download:`Download Python source code: plot_phase_bandwidth.py <plot_phase_bandwidth.py>`



  .. container:: sphx-glr-download

     :download:`Download Jupyter notebook: plot_phase_bandwidth.ipynb <plot_phase_bandwidth.ipynb>`

.. rst-class:: sphx-glr-signature

    `Generated by Sphinx-Gallery <https://sphinx-gallery.readthedocs.io>`_