gem.solvers.simple_ed.SimpleED.solve_Hemb

SimpleED.solve_Hemb(num_eig=None, which=None, tol=None, dense_cutoff=None, bw_cutoff=None, verbose=0, T=0.0)[source]

Diagonalise each sector and build the global (thermal) partition function.

Parameters:

  • num_eig – int, optional. Number of eigenvalues to compute. Read from solver_params when not given; if it is None there too, only the ground state is computed at T=0 while at T>0 the full Hamiltonian of each sector is diagonalised.

  • dense_cutoff – int, optional. Sectors of dimension smaller than this are diagonalised fully (eigh), the larger ones partially (eigsh). Read from solver_params when not given, default 4000.

  • bw_cutoff – float, optional. At T>0, states whose Boltzmann weight exp(-beta*(E-gs_ene)) falls below this are dropped from the partition function. Read from solver_params when not given, default 1e-8.