.. _sphx_glr_auto_examples_twosite_plot_halffill_z.py:


================================================
Following the Metal to Mott insulator Transition
================================================

Plot of the Quasiparticle weigth decay for the Hubbard
Model in the Bethe Lattice as the local interaction is raised.




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


    *

      .. image:: /auto_examples/twosite/images/sphx_glr_plot_halffill_z_001.png
            :scale: 47

    *

      .. image:: /auto_examples/twosite/images/sphx_glr_plot_halffill_z_002.png
            :scale: 47





.. code-block:: python


    # Code source: Óscar Nájera
    # License: BSD 3 clause

    from __future__ import division, absolute_import, print_function
    import matplotlib.pyplot as plt
    import numpy as np
    from dmft.twosite import dmft_loop


    def plot_z(axis='matsubara', du=0.05):
        fig = plt.figure()
        u_int = np.arange(0, 6.2, du)
        for beta in [6, 10, 20, 30, 50, 100, 1e3]:
            out_file = axis + '_halffill_b{}_dU{}'.format(beta, du)
            try:
                res = np.load(out_file + '.npy')
            except IOError:
                res = dmft_loop(u_int, axis, beta=beta, hop=1)
                np.save(out_file, res)

            plt.plot(res[:, 0] / 2, res[:, 1], '+-',
                     label='$\\beta = {}$'.format(beta))

        plt.legend(loc=0)
        plt.title('Quasiparticle weigth, estimated in {} frequencies'.format(axis))
        plt.ylabel('Z')
        plt.xlabel('U/D')
        fig.savefig(out_file + '_Z.png', format='png',
                    transparent=False, bbox_inches='tight', pad_inches=0.05)


    plot_z('real')
    plot_z('matsubara')

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



.. container:: sphx-glr-footer


  .. container:: sphx-glr-download

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



  .. container:: sphx-glr-download

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

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

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