A Monte Carlo algorithm for the random Boson Hubbard model is reported. The analytic expression of the matrix elements is presented, and the ergodicity of the Monte Carlo flips is discussed. The results in one dimension supports a previously proposed perturbational scaling argument.
We study the effect of diagonal disorder on a one-dimensional (1-D) Hubbard model by means of a Monte Carlo simulation. Contrary to previous analytical studies of interacting disordered systems in one dimension, the numerical simulation is not restricted to the weak coupling regime in disorder and/or interaction. Quantum fluctuations are introduced using the scheme developed by J. Hirsch [Phys. Rev. B 31, 4403 (1985)] for Grand Canonical ensembles and arbitrary dimensions. A second Monte Carlo loop creates the configurations for the disorder potential. Preliminary results for various correlation functions are presented.
We have studied the effects of random-hopping matrix elements and random potentials on the properties of the one-dimensional Hubbard model. Using a quantum Monte Carlo technique, disorder-averaged static spin- and charge-density susceptiblities have been evaluated for various strengths of the disorder. Results for the spin susceptibility at wave number q=2kF indicate that this quantity, which is the fastest diverging susceptibility of the pure system, diverges as T\ensuremath{\rightarrow}0 also when there is randomness in the hopping matrix elements, but not in the presence of random potentials. Both types of disorder cause a divergence of the uniform magnetic susceptibility. However, for random potentials a finite critical strength of the disorder appears to be required. At half-filling the transition from Mott (gapped) to Anderson (gapless) insulating behavior has been studied. A critical disorder strength is needed to destroy the gap, in agreement with Ma's renormalization group calculations.
We extend a Quantum Monte Carlo (QMC) algorithm to the Hubbard model containing bond-dependent off-diagonal hoppings, which describes the manganites and other similar systems, such as systems with spin-flip hoppings. To show correctness of this extension, we detect two simple examples of one-dimensional chains. The first one is half-filled and the hopping matrixes between arbitrary nearest neighboring pairs are uniform, i.e., ti,i+1=ti+1,i+2(i=1,2,3,…). The other one is quarter-filled and the hopping matrixes alternate every four lattice constants with ${t}_{i,i+1}{=t}_{i+1,i+2}\ensuremath{\ne}{t}_{i+2,i+3}{=t}_{i+3,i+4}$ (i=1,5,9,…). Results of the two examples are reasonable.