What you know you can't explain, but you feel it. You've felt it your entire life, that there's something wrong with the world. You don't know what it is, but it's there, like a splinter in your mind, driving you mad.
What you know you can't explain, but you feel it. You've felt it your entire life, that there's something wrong with the world. You don't know what it is, but it's there, like a splinter in your mind, driving you mad.
 |
NCTS /
BoseHubbardIn this tutorial, we will examine a more realistic example, of the Bose-Hubbard model in a 1D harmonic trap. The Bose-Hubbard model is  Note that the Hamiltonian preserves particle-number symmetry, and as we are working in a finite 1D system it is not spontaneously broken. We define the trapping potential as  The homogeneous system with no trap has a phase diagram (from Freericks and Monien 1994),  This model can be constructed using the program Next, make a random wavefunction with N=60 particles. The hopping, coulomb and trap terms in the Hamiltonian are implemented as the operators To see some paramters in the superfluid region, try the Hamiltonian "H_J + 1*H_U + 50*H_trap". To see some fringes of Mott insulator, try increasing U, to 10 or 20. m=50 states is probably sufficient. It will take a long time to converge because of the non-uniform groundstate! You will probably need up to a few hundred sweeps. You could try starting off with a smaller number of states so that the program runs faster intially, and then increase the number of states towards the end of the calculation. What is the difference in the local density profile as U is increased? Try the Try Is a maximum of 5 bosons sufficient for small U? The lattice file defines projectors onto the particle number states, P_1, P_2, ...., P_N which can be used to calculate the occupations. |