dmft.ipt_real

IPT in real frequencies

Functions

dmft.ipt_real.dimer_dmft(U, tp, nfp, w, dw, gss, gsa, conv=1e-07, t=0.5)

Solve DMFT equations in real frequencies for the dimer

Parameters:
  • U (float) – Couloumb interaction
  • tp (float) – Dimerization
  • npf (1D real ndarray) – Thermal Fermi function
  • w (1D real ndarray) – frequency grid. Has to be equispaced and symmetric
  • dw (float) – frequency separation
  • gss (1D complex ndarray) – Starting guess for the symmetric Green function
  • gsa (1D complex ndarray) – Starting guess for the asymmetric Green function
  • conv (float) – convergence criteria
  • t (float) – hopping
Returns:

  • (gss, gsa) (tuple of 1D complex ndarray, Green Functions)
  • (ss, sa) (tuple of 1D complex ndarray, Self-Energy)

Examples using dmft.ipt_real.dimer_dmft

dmft.ipt_real.dimer_solver(w, dw, tp, U, nfp, gss, gsa, t=0.5, eta=0.003j)
dmft.ipt_real.ph_hf_sigma(Aw, nf, U)

Imaginary part of the second order diagram

because of particle-hole symmetry at half-fill in the Single band one can work with A^+ only

dmft.ipt_real.sigma(Aw, nf, U)
dmft.ipt_real.ss_dmft_loop(gloc, w, u_int, beta, conv)

DMFT Loop for the single band Hubbard Model at Half-Filling

Parameters:
  • gloc (complex 1D ndarray) – local Green’s function to use as seed
  • w (real 1D ndarray) – real frequency points
  • u_int (float) – On site interaction, Hubbard U
  • beta (float) – Inverse temperature
  • conv (float) – convergence criteria
Returns:

  • gloc (complex 1D ndarray) – DMFT iterated local Green’s function
  • sigma (complex 1D ndarray) – DMFT iterated self-energy

Examples using dmft.ipt_real.ss_dmft_loop