Undergrad How to simulate a lattice with boundary conditions?

Click For Summary
The discussion centers on simulating a lattice with boundary conditions for the Hubbard model, specifically focusing on creating a simple lattice, such as a square lattice, and modeling two sites with two electrons. The user seeks guidance on coding the lattice structure and linking sites with a hopping parameter to implement periodic boundary conditions. There is mention of the need to create a Hamiltonian matrix for the kinetic part of the Hubbard model, followed by diagonalization to obtain eigenvalues and compute the density of states. The conversation also touches on the distinction between mean field approximations and quantum Monte Carlo simulations, clarifying their roles in the modeling process. Overall, the user is looking for practical coding advice and resources to facilitate their simulation efforts.
Amentia
Messages
108
Reaction score
5
Hello,

My question is very simple but I do not have a lot of experience with simulation. I want to write some code to simulate a lattice with boundary conditions and then I will perform calculations with the Hubbard model to find different kinds of properties of interest. I would like to know how to do this first step, for any lattice (you can choose a simple one like a square lattice with few atoms to give an example).

Thank you!
 
Physics news on Phys.org
Maybe, you might use "periodic boundary" conditions, but I am not really sure whether it works when simulating a Hubbard model on a two dimensional lattice. I am sure, however, that there is many stuff around if you search in this direction.
 
Hello,

Thank you for your answer. In fact, I did not find anything, maybe because it is too simple. My question is even more simple than that: assume no boundary conditions, what should I write to model a system like two sites with two electrons that can be on either site with a hopping parameter t=1? Because I think once I am sure how to link the sites numerically, I should just have to link bottom atomic sites with top atomic sites with the correct t or 0 hopping paramater, in order to simulate the periodic boundary conditions. I know it is possible for a Hubbard model on a 2D lattice because there is already literature on that. But they give their results, not their codes.
 
You should search for quantum monte carlo simulations regarding the Hubbard model.
 
Ok thank you. By the way, is the mean field approximation included into the quantum monte carlo simulations?
 
Amentia said:
Ok thank you. By the way, is the mean field approximation included into the quantum monte carlo simulations?

Mean field approximations provide analytical, approximative solutions. Monte Carlo simulations provide numerical solutions.
 
I am still confused. Maybe I will detail more what I want to do.

1) Create a lattice with sites that can have up to two electrons and link them with a hopping parameter t (main question of this thread)
2) Use this lattice to make a hamiltonian matrix for the kinetic part of the Hubbard model
3) Diagonalize and get the eigenvalues
4) Compute density of states
5) Knowing the number of electrons and the number of sites, find the chemical potential for finite temperature
6) Use a mean field approximation for the second part of the Hamiltonian
7) Plot a curve to see how the material will behave, for example if it is supposed to be ferromagnetic, if there are phase transitions...
 
Thank you for the book. I think it is much more involved than what I want to do now but I will keep it for later. It could be of great help.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
13K
  • · Replies 1 ·
Replies
1
Views
6K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
14
Views
4K
Replies
1
Views
4K
  • · Replies 0 ·
Replies
0
Views
3K