CalculateRho(H) computes the density matrix for a given Tight Binding model.
Input
H : Tight Binding Object
Possible options are:
“Nk” List of integers the format {Nx,Ny,Nz} giving the number of Crystal momentum states used to solve the model. If the reciprocal lattice is zero, then Nk is always set to {0,0,0}, independent of the option given (Standard value {40,40,40})
“mu” Real number giving the chemical potential μ. (Standard value 0)
“kbT” or “kBT” Positive real number indicating the value of β=kBT for the calculation. (Standard value DBL_EPSILON)
Output
ρ : Matrix object with size Nx*Ny*Nz*Norbs, where Norbs is the number of defined orbitals per unit cell in the Tight Binding input.