.. _sphx_glr_auto_examples_twosite_plot_dop_z.py:


=============================================================
Quasiparticle decay for different local interaction strenghts
=============================================================

Ploting the Quasiparticle weight at different local interaction strenghts
at different doping amounts




.. image:: /auto_examples/twosite/images/sphx_glr_plot_dop_z_001.png
    :align: center





.. 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_dop import dmft_loop_dop


    def plot_z(axis='real', beta=1e3, u_int=[2., 4., 5.85, 6., 8., 10., 100.]):
        fig = plt.figure()
        for u in u_int:
            out_file = axis + '_dop_b{}_U{}'.format(beta, u)
            try:
                res = np.load(out_file + '.npy')
            except IOError:
                res = dmft_loop_dop(u)
                np.save(out_file, res)

            zet = [sim.imp_z() for sim in res[:, 1]]
            plt.plot(res[:, 0], zet, '+-', label='$U/t= {}$'.format(u))

        plt.legend(loc=0)
        plt.title(
            'Quasiparticle weigth, estimated in real freq at $\\beta={}$'.format(beta))
        plt.ylabel('Z')
        plt.xlabel('n')
        plt.xlim([0, 1])
        fig.savefig(out_file + '_Z.png', format='png',
                    transparent=False, bbox_inches='tight', pad_inches=0.05)

    plot_z()

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



.. container:: sphx-glr-footer


  .. container:: sphx-glr-download

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



  .. container:: sphx-glr-download

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

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

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