Documentation
The Hubbard-I solver approximates the solution of the impurity model by neglecting any hybridization and solving an atomic problem. The atomic problem is defined by the local interaction Hamiltonian and the orbital dependent high-frequency behavior of the non-interacting bath Green’s function.
C++ reference manual
Python reference manual
- class triqs_hubbardI.Solver(beta, gf_struct, n_iw=1025, n_tau=10001, n_l=30, n_w=500, w_min=-15, w_max=15, idelta=0.01)[source]
Class providing initialization and solve function. Contains all relevant Greensfunctions and self energy.
Methods
solve
(**params_kw)Solve the impurity problem: calculate G(iw) and Sigma(iw)
- __init__(beta, gf_struct, n_iw=1025, n_tau=10001, n_l=30, n_w=500, w_min=-15, w_max=15, idelta=0.01)[source]
Initialise the solver.
- Parameters:
- betascalar
Inverse temperature.
- gf_structlist of pairs [ (str,int), …]
Structure of the Green’s functions. It must be a list of pairs, each containing the name of the Green’s function block and its linear size. For example:
[ ('up', 3), ('down', 3) ]
.- n_iwinteger, optional
Number of Matsubara frequencies used for the Green’s functions.
- n_tauinteger, optional
Number of imaginary time points used for the Green’s functions.
- n_linteger, optional
Number of legendre polynomials used for the Green’s functions.
- n_winteger, optional
Number of real frequency points used for the Green’s functions.
- w_mininteger, optional
Lower limit of the range of real frequencies
- w_maxinteger, optional
Upper limit of the range of real frequencies
- ideltafloat, optional
Broadening of Green’s function on real frequencies
- solve(**params_kw)[source]
Solve the impurity problem: calculate G(iw) and Sigma(iw)
- Parameters:
- params_kwdict {‘param’:value} that is passed to the core solver.
- Only required parameter is
h_int (Operator object): the local Hamiltonian of the impurity problem to be solved,
- Other parameters are
calc_gtau (bool): calculate G(tau)
calc_gw (bool): calculate G(w) and Sigma(w)
calc_gl (bool): calculate G(legendre)
calc_dm (bool): calculate density matrix