Fourier acceleration of iterative processes in disordered systems

Technical details are given on how to use Fourier acceleration with iterative processes such as relaxation and conjugate gradient methods. These methods are often used to solve large linear systems of equations, but become hopelessly slow very rapidly as the size of the set of equations to be solved increases. Fourier acceleration is a method designed to alleviate these problems and result in a very fast algorithm. The method is explained for the Jacobi relaxation and conjugate gradient methods and is applied to two models: the random resistor network and the random central-force network. In the first model, acceleration works very well; in the second, little is gained. We discuss reasons for this. We also include a discussion of stopping criteria.     

Phase separation in the two-dimensional bosonic Hubbard model with ring exchange

We show that soft-core bosons in two dimensions with a ring exchange term exhibit a tendency for phase separation. This observation suggests that the thermodynamic stability of normal Bose liquid phases driven by ring exchange should be carefully examined.     

Disorder averaging and finite-size scaling

We propose a new picture of the renormalization group (RG) approach in the presence of disorder, which considers the RG trajectories of each random sample (realization) separately instead of the usual renormalization of the averaged free energy. The main consequence of the theory is that the average over randomness has to be taken after finding the critical point of each realization. To demonstrate these concepts, we study the finite-size scaling properties of the two-dimensional random-bond Ising model. We find that most of the previously observed finite-size corrections are due to the sample-to-sample fluctuation of the critical temperature and scaling predictions are fulfilled only by the new average.     

Supersolid phases in the one-dimensional extended soft-core bosonic Hubbard model

We present results of quantum Monte Carlo simulations for the soft-core extended bosonic Hubbard model in one dimension exhibiting the presence of supersolid phases similar to those recently found in two dimensions. We find that in one and two dimensions, the insulator-supersolid transition has dynamic critical exponent z = 2 whereas the first order insulator-superfluid transition in two dimensions is replaced by a continuous transition with z = 1 in one dimension. We present evidence that this transition is in the Kosterlitz-Thouless universality class and discuss the mechanism behind this difference. The simultaneous presence of two types of quasi-long-range order results in two solitonlike dips in the excitation spectrum.     

Exact numerical study of pair formation with imbalanced fermion populations

We present an exact quantum Monte Carlo study of the attractive one-dimensional Hubbard model with imbalanced fermion population. The pair-pair correlation function, which decays monotonically in the absence of polarization P, develops oscillations when P is nonzero, characteristic of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase. The pair momentum distribution peaks at a momentum equal to the difference in the Fermi momenta. At strong coupling, the minority and majority momentum distributions are shown to be deformed, reflecting the presence of the other species and its Fermi surface. The FFLO oscillations survive the presence of a confining potential, and the local polarization at the trap center exhibits a marked dip, similar to that observed experimentally.