The Boson-Hubbard Model on a Kagome Lattice with Sextic Ring-Exchange Terms

High order ring-exchange interactions are crucial for the study of quantum fluctuations on highly frustrated systems. We present the first exact quantum Monte Carlo study of a model of hard-core bosons with sixth order ring-exchange interactions on a two-dimensional kagome lattice. By using the Stochastic Green Function algorithm, we show that the system becomes unstable in the limit of large ring-exchange interactions. It undergoes a phase separation at all fillings, except at 1/3 and 2/3 fillings for which the superfluid density vanishes and an unusual mixed valence bond and charge density ordered solid is formed.Comment: 4 pages, 7 figure

Parallel Tempering Simulation of the three-dimensional Edwards-Anderson Model with Compact Asynchronous Multispin Coding on GPU

Monte Carlo simulations of the Ising model play an important role in the field of computational statistical physics, and they have revealed many properties of the model over the past few decades. However, the effect of frustration due to random disorder, in particular the possible spin glass phase, remains a crucial but poorly understood problem. One of the obstacles in the Monte Carlo simulation of random frustrated systems is their long relaxation time making an efficient parallel implementation on state-of-the-art computation platforms highly desirable. The Graphics Processing Unit (GPU) is such a platform that provides an opportunity to significantly enhance the computational performance and thus gain new insight into this problem. In this paper, we present optimization and tuning approaches for the CUDA implementation of the spin glass simulation on GPUs. We discuss the integration of various design alternatives, such as GPU kernel construction with minimal communication, memory tiling, and look-up tables. We present a binary data format, Compact Asynchronous Multispin Coding (CAMSC), which provides an additional $28.4\%$ speedup compared with the traditionally used Asynchronous Multispin Coding (AMSC). Our overall design sustains a performance of 33.5 picoseconds per spin flip attempt for simulating the three-dimensional Edwards-Anderson model with parallel tempering, which significantly improves the performance over existing GPU implementations.Comment: 15 pages, 18 figure

Complex phases in the doped two-species bosonic Hubbard Model

We study a two-dimensional bosonic Hubbard model with two hard-core species away from half filling using Quantum Monte Carlo simulations. The model includes a repulsive interspecies interaction and different nearest-neighbor hopping terms for the two species. By varying the filling we find a total of five distinct phases, including a normal liquid phase at higher temperature, and four different phases at lower temperature. We find an anti-ferromagnetically ordered Mott insulator and a region of coexistent anti-ferromagnetic and superfluid phases near half filling. Further away from half filling the phase diagram displays a superfluid phase and a novel phase inside the superfluid region at even lower temperatures. In this novel phase separated region, the heavy species has a Mott behavior with integer filling, while the lighter species shows phase separated Mott and superfluid behaviors.Comment: 5 pages, 4 figure

Periodic Anderson model with Holstein phonons for the description of the Cerium volume collapse

Recent experiments have suggested that the electron-phonon coupling may play an important role in the $\gamma \rightarrow \alpha$ volume collapse transition in Cerium. A minimal model for the description of such transition is the periodic Anderson model. In order to better understand the effect of the electron-phonon interaction on the volume collapse transition, we study the periodic Anderson model with coupling between Holstein phonons and electrons in the conduction band. We find that the electron-phonon coupling enhances the volume collapse, which is consistent with experiments in Cerium. While we start with the Kondo Volume Collapse scenario in mind, our results capture some interesting features of the Mott scenario, such as a gap in the conduction electron spectra which grows with the effective electron-phonon coupling.Comment: 8 pages, 6 figure

Phase diagram of the Bose-Hubbard model on a ring-shaped lattice with tunable weak links

Motivated by recent experiments on toroidal Bose-Einstein condensates in all-optical traps with tunable weak links, we study the one-dimensional Bose-Hubbard model on a ring-shaped lattice with a small region of weak hopping integrals using quantum Monte Carlo simulations. Besides the usual Mott insulating and superfluid phases, we find a phase which is compressible but non superfluid with a local Mott region. This `local Mott' phase extends in a large region of the phase diagram. These results suggest that the insulating and conducting phases can be tuned by a local parameter which may provide a new insight to the design of atomtronic devices.Comment: 5 pages, 5 figure

Periodic Anderson model with electron-phonon correlated conduction band

This paper reports dynamical mean field calculations for the periodic Anderson model in which the conduction band is coupled to phonons. Motivated in part by recent attention to the role of phonons in the $\gamma$-$\alpha$ transition in Ce, this model yields a rich and unexpected phase diagram which is of intrinsic interest. Specifically, above a critical value of the electron-phonon interaction, a first order transition with two coexisting phases develops in the temperature-hybridization plane, which terminates at a second order critical point. The coexisting phases display the familiar Kondo screened and local moment character, yet they also exhibit pronounced polaronic and bipolaronic properties, respectively.Comment: 4 pages, 6 figure

Study of multiband disordered systems using the typical medium dynamical cluster approximation

We generalize the typical medium dynamical cluster approximation to multiband disordered systems. Using our extended formalism, we perform a systematic study of the non-local correlation effects induced by disorder on the density of states and the mobility edge of the three-dimensional two-band Anderson model. We include inter-band and intra-band hopping and an intra-band disorder potential. Our results are consistent with the ones obtained by the transfer matrix and the kernel polynomial methods. We apply the method to K$_x$Fe$_{2-y}$Se$_2$ with Fe vacancies. Despite the strong vacancy disorder and anisotropy, we find the material is not an Anderson insulator. Our results demonstrate the application of the typical medium dynamical cluster approximation method to study Anderson localization in real materials.Comment: 10 pages, 8 figure

Locally self-consistent embedding approach for disordered electronic systems

We present a new embedding scheme for the locally self-consistent method to study disordered electron systems. We test this method in a tight-binding basis and apply it to the single band Anderson model. The local interaction zone is used to efficiently compute the local Green's function of a supercell embeded into a local typical medium. We find a quick convergence as the size of the local interaction zone which reduces the computational costs as expected. This method captures the Anderson localization transition and accurately predicts the critical disorder strength. The present work opens the path towards the development of a typical medium embedding scheme for the $O(N)$ multiple scattering methods.Comment: 7 pages, 5 figure

Experimental investigation of flow-induced vibration of a rotating circular cylinder

While flow-induced vibration of bluff bodies has been extensively studied over the last half-century, only limited attention has been given to flow-induced vibration of elastically mounted rotating cylinders. Since recent low-Reynolds-number numerical work suggests that rotation can enhance or suppress the natural oscillatory response, the former could find applications in energy harvesting and the latter in vibration control. The present experimental investigation characterises the dynamic response and wake structure of a rotating circular cylinder undergoing vortex-induced vibration at a low mass ratio (m∗ = 5.78) over the reduced velocity range leading to strong oscillations. The experiments were conducted in a free-surface water channel with the cylinder vertically mounted and attached to a motor that provided constant rotation. Springs and an air-bearing system allow the cylinder to undertake low-damped transverse oscillations. Under cylinder rotation, the normalised frequency response was found to be comparable to that of a freely vibrating non-rotating cylinder. At reduced velocities consistent with the upper branch of a non-rotating transversely oscillating cylinder, the maximum oscillation amplitude increased with non-dimensional rotation rate up to α ≈ 2. Beyond this, there was a sharp decrease in amplitude. Notably, this critical value corresponds approximately to the rotation rate at which vortex shedding ceases for a non-oscillating rotating cylinder. Remarkably, at α = 2 there was approximately an 80% increase in the peak amplitude response compared to that of a non-rotating cylinder. The observed amplitude response measured over the Reynolds-number range of (1100 Re 6300) is significantly different from numerical predictions and other experimental results recorded at significantly lower Reynolds numbers