I need to determine numerically the ground state energy of the Fermi-Hubbard model when total filling is n=1, but the number of fermions with spin up: 'n_up', and spin down: 'n_down' can take any value (but of course n_up + n_down = 1).
Do you know of any C++/Fortran/Python etc. ready to use package that would allow me to determine ground state energy for different hopings 'J', interaction energies 'U' and compositions n_up, n_down?
Did you try the following script?
https://github.com/georglind/fermihubbard
Thank you @Ramin! Sorry for such a late response, but I was very busy lately. The script is based on exact diagonalization. Any idea for approximate techniques that could be applied to large lattices with conserved number of particles? Also, I would be satisfied with "mean-field" methods for 3D lattices, but where two chemical potentials can be specified.
In the case of bosons, there is TEBD or DMRG for 1d lattices. 3d lattices can be attacked with quantum Monte Carlo.
- Badges
- Science topic
- Similar topics
- Computer Science and Engineering
- Programming Languages
- Python