from EasyED.ensembles import GrandcanonicalEnsemble from EasyED.hamiltonians import Hubbard, Hubbard2 from EasyED.util import report from reduceddensitymatrix import ReducedDensityMatrix as RDM beta = 10 t = -1 tp = 0 u = 3 mu = u*.5 hop = [[0,t,t,tp],[t,0,tp,t],[t,tp,0,t],[tp,t,t,0]] hamiltonian = Hubbard2(hop, u) plaq = GrandcanonicalEnsemble(hamiltonian, beta, mu) plaq.calcEigensystem() e0 = plaq.hamiltonian.getGroundStateEnergy() energies = plaq.hamiltonian.eigenEnergies eigenstates = plaq.hamiltonian.eigenStates.toarray() rho_dim = RDM(beta, eigenstates, energies, 4, 4, [], [0]) rho_dim.calculate() print rho_dim.get_entanglement_spectrum() print rho_dim.get_entanglement_entropy() hamiltonian = Hubbard(hop, u) print hamiltonian.singleParticleBasis plaq = GrandcanonicalEnsemble(hamiltonian, beta, mu) plaq.calcEigensystem() e0 = plaq.hamiltonian.getGroundStateEnergy() energies = plaq.hamiltonian.eigenEnergies eigenstates = plaq.hamiltonian.eigenStates.toarray()
