Within the hard-core boson model and the 1/S spin wave approximation, we have constructed the linear response theory for disordered boson system. Our numerical results shows that the instability of the superfluidity occurs when the paramagnetic part cancels out the diamagnetic part with a critical finite amout of disorder.
The ground-state of hard-core bosons in a random one-body potential is studied numerically to observe the transition from superfluid to insulating ``Bose glass'' phase. Finite-size scaling analyses are performed to estimate the transition point and critical exponents. The dynamical critical exponent is z=1.95\ifmmode\pm\else\textpm\fi{}0.25 and is consistent with the prediction of Fisher et al. [Phys. Rev. B 40, 546 (1989)]. Values for other exponents and preliminary results regarding the universal conductivity are presented.
Excitation spectra of hard-core bosons with nearest- and next-nearest-neighbor interactions on square and triangular lattices in the superfluid phase are studied using the spin-wave theory. We go beyond the standard linear spin-wave approximation calculating corrections due to spin-wave interaction. Account for the spin-wave interaction considerably improves the quantitative description of the excitation spectra. The behavior of the spectra is analyzed in detail at different particle densities and relations between the interaction parameters. Particular attention is given to the spectrum minima with the aim to find out the conditions for instability of the superfluid phase. We show that for a triangular lattice at high enough next-nearest-neighbor repulsion there appears a true roton minimum inside the Brillouin zone located on the line connecting the zone center with the midpoint of the face of the zone boundary. The obtained analytical results are in excellent quantitative agreement with the numerical simulations known from literature.
In this dissertation, we study the superfluid-insulator quantum phase transition in disordered boson systems. Recently, there has been considerable controversy over the validity of the scaling relations of the superfluid--Bose-glass quantum phase transition in three dimensions. Results from experimental and numerical studies on disordered quantum magnets contradict the scaling relations and the associated conventional scaling hypothesis for the singular part of the free energy. We determine various critical exponents of the superfluid--Bose-glass quantum phase transition in three-dimensional disordered Bose-Hubbard model through extensive Monte Carlo simulations. Our numerical study shows the previous studies on disordered quantum magnets were performed outside the quantum critical region, and our results for the critical exponents are in perfect agreement with previous theoretical predictions. We next move to study theoretically and numerically the superfluid-insulator quantum phase transition in one-dimensional disordered boson systems. While the superfluid-insulator transition in the weak disorder limit is well understood through the perturbative renormalization group study by Giamarchi and Schulz, transitions in the strong disorder regime are beyond the reach of this method. This problem was recently attacked by Altman et al. with the real space renormalization group method. They reached the conclusion that the superfluid-insulator transition in the strong disorder regime can be explained by the Coulomb blockade physics of weak links. However, their method is actually uncontrolled. Taking account of the crucial role of the hydrodynamic renormalization of weak links, finally, Pollet et al. put forward an asymptotically exact renormalization group theory of the superfluid-insulator transition. Based on this theory, we are able to provide an accurate description of the interplay between the well-known Giamarchi-Schulz criticality and the new weak-link criticality. A significant part of the ground-state phase diagram of one-dimensional disordered Bose-Hubbard model at unit filling is also determined numerically. In particular, we established the position of the multicritical point beyond which the new weak-link criticality takes over of the Giamarchi–Schulz criticality.
We study the T=0 critical properties of the superfluid-insulator transition in 2D hard core Bose systems with disorder. Using quantum Monte Carlo simulations and finite size scaling on 64\ifmmode\times\else\texttimes\fi{}64 lattices we find the dynamical exponent z=0.5\ifmmode\pm\else\textpm\fi{}0.1 and the correlation length exponent \ensuremath{\nu}=2.2\ifmmode\pm\else\textpm\fi{}0.2. At the transition, the system is metallic with a conductivity ${\mathrm{\ensuremath{\sigma}}}_{\mathit{c}}$=(1.2\ifmmode\pm\else\textpm\fi{}0.2)(e∗)2/h and a nonzero compressibility \ensuremath{\kappa}\ensuremath{\ne}0. These conclusions differ from the existing scaling theory as well as from simulations on simplified models argued to be in the same universality class. Our results are suggestive of new low lying collective excitations (modified from usual phonons) in the disordered system.
We investigate the critical phenomena of hard-core bosons in a disordered medium. Such a system is mapped onto a quantum spin-1/2 XY model with transverse random field. The system then is studied through a quantum real-space renormalization-group method. We find that randomness is always relevant in a one-dimensional (1D) system, in agreement with exact results. In two and three dimensions, there is a critical amount of disorder, below which the superfluid phase is stable. In 2D, the dynamic exponent z=1.7 for compressible states, and is close to the value of z=d as predicted by Fisher et al. $z--- is smaller for incompressible states. The correlation length exponent \ensuremath{\nu} is insensitive to z, and roughly equals 1.4. Unlike the superfluid--Mott-insulator transition without disorder, which has two distinct universality classes, we find there is only one universality class for the superfluid--Bose-glass transition.