Generalization of the Fortuin-Kasteleyn transformation and its application to quantum spin simulations,

We generalize the Fortuin-Kasteleyn (FK) cluster representation of the partition function of the Ising model to represent the partition function of quantum spin models with an arbitrary spin magnitude in arbitrary dimensions. This generalized representation enables us to develop a new cluster algorithm for the simulation of quantum spin systems by the worldline Monte Carlo method. Because the Swendsen-Wang algorithm is based on the FK representation, the new cluster algorithm naturally includes it as a special case. As well as the general description of the new representation, we present an illustration of our new algorithm for some special interesting cases: the Ising model, the antiferromagnetic Heisenberg model with $S=1$, and a general Heisenberg model. The new algorithm is applicable to models with any range of the exchange interaction, any lattice geometry, and any dimensions.Comment: 46 pages, 10 figures, to appear in J.Stat.Phy

Strong-coupling expansion for the momentum distribution of the Bose Hubbard model with benchmarking against exact numerical results

A strong-coupling expansion for the Green's functions, self-energies and correlation functions of the Bose Hubbard model is developed. We illustrate the general formalism, which includes all possible inhomogeneous effects in the formalism, such as disorder, or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third-order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator to superfluid transition along with a generalization of the RPA-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions; the accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling theory results are benchmarked against numerically exact QMC simulations in two and three dimensions and against DMRG calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.Comment: 48 pages 14 figures RevTe

Dual Monte Carlo and Cluster Algorithms

We discuss the development of cluster algorithms from the viewpoint of probability theory and not from the usual viewpoint of a particular model. By using the perspective of probability theory, we detail the nature of a cluster algorithm, make explicit the assumptions embodied in all clusters of which we are aware, and define the construction of free cluster algorithms. We also illustrate these procedures by rederiving the Swendsen-Wang algorithm, presenting the details of the loop algorithm for a worldline simulation of a quantum $S=$ 1/2 model, and proposing a free cluster version of the Swendsen-Wang replica method for the random Ising model. How the principle of maximum entropy might be used to aid the construction of cluster algorithms is also discussed.Comment: 25 pages, 4 figures, to appear in Phys.Rev.

Magneto-hydrodynamic Simulations of a Jet Drilling an HI Cloud: Shock Induced Formation of Molecular Clouds and Jet Breakup

The formation mechanism of the jet-aligned CO clouds found by NANTEN CO observations is studied by magnetohydrodynamical (MHD) simulations taking into account the cooling of the interstellar medium. Motivated by the association of the CO clouds with the enhancement of HI gas density, we carried out MHD simulations of the propagation of a supersonic jet injected into the dense HI gas. We found that the HI gas compressed by the bow shock ahead of the jet is cooled down by growth of the cooling instability triggered by the density enhancement. As a result, cold dense sheath is formed around the interface between the jet and the HI gas. The radial speed of the cold, dense gas in the sheath is a few km/s almost independent of the jet speed. Molecular clouds can be formed in this region. Since the dense sheath wrapping the jet reflects waves generated in the cocoon, the jet is strongly perturbed by the vortices of the warm gas in the cocoon, which breaks up the jet and forms a secondary shock in the HI-cavity drilled by the jet. The particle acceleration at the shock can be the origin of radio and X-ray filaments observed near the eastern edge of W50 nebula surrounding the galactic jet source SS433.Comment: 30 pages, 16 figure