A cluster algorithm for Lattice Gauge Theories

A new algorithm for simulating compact U(1) lattice gauge theory in three dimensions is presented which is based on global changes in the configuration space. We show that this algorithm provides an effective way to extract partition functions at given external flux. As an application, we study numerically the finite temperature deconfinement phase transition.Comment: 4 pages, 2 figures. Talk given at the Conference on Computational Physics, Genova, Italy, Sept. 200

Phase diagram of an extended classical dimer model

We present an extensive numerical study of the critical behavior of dimer models in three dimensions, focusing on the phase transition between Coulomb and crystalline columnar phases. The case of attractive interactions between parallel dimers on a plaquette was shown to undergo a continuous phase transition with critical exponents close to those of the O(N) tricritical universality class, a situation which is not easily captured by conventional field theories. That the dimer model is exactly fine-tuned to a highly symmetric point is a non trivial statement which needs careful numerical investigation. In this paper, we perform an extensive Monte Carlo study of a generalized dimer model with plaquette and cubic interactions and determine its extended phase diagram. We find that when both interactions favor alignment of the dimers, the phase transition is first order, in almost all cases. On the opposite, when interactions compete, the transition becomes continuous, with a critical exponent \eta ~ 0.2. The existence of a tricritical point between the two regimes is confirmed by simulations on very large size systems and a flowgram method. In addition, we find a highly-degenerate crystalline phase at very low temperature in the frustrated regime which is separated from the columnar phase by a first order transition.Comment: 12 pages, 13 figure

Neural network setups for a precise detection of the many-body localization transition: finite-size scaling and limitations

Determining phase diagrams and phase transitions semi-automatically using machine learning has received a lot of attention recently, with results in good agreement with more conventional approaches in most cases. When it comes to more quantitative predictions, such as the identification of universality class or precise determination of critical points, the task is more challenging. As an exacting test-bed, we study the Heisenberg spin-1/2 chain in a random external field that is known to display a transition from a many-body localized to a thermalizing regime, which nature is not entirely characterized. We introduce different neural network structures and dataset setups to achieve a finite-size scaling analysis with the least possible physical bias (no assumed knowledge on the phase transition and directly inputing wave-function coefficients), using state-of-the-art input data simulating chains of sizes up to L=24. In particular, we use domain adversarial techniques to ensure that the network learns scale-invariant features. We find a variability of the output results with respect to network and training parameters, resulting in relatively large uncertainties on final estimates of critical point and correlation length exponent which tend to be larger than the values obtained from conventional approaches. We put the emphasis on interpretability throughout the paper and discuss what the network appears to learn for the various used architectures. Our findings show that a it quantitative analysis of phase transitions of unknown nature remains a difficult task with neural networks when using the minimally engineered physical input.Comment: v2: published versio

The semiflexible fully-packed loop model and interacting rhombus tilings

Motivated by a recent adsorption experiment [M.O. Blunt et al., Science 322, 1077 (2008)], we study tilings of the plane with three different types of rhombi. An interaction disfavors pairs of adjacent rhombi of the same type. This is shown to be a special case of a model of fully-packed loops with interactions between monomers at distance two along a loop. We solve the latter model using Coulomb gas techniques and show that its critical exponents vary continuously with the interaction strenght. At low temperature it undergoes a Kosterlitz-Thouless transition to an ordered phase, which is predicted from numerics to occur at a temperature T \sim 110K in the experiments.Comment: 4 pages, 4 figures, v2: corrected typo, v3: minor modifications, published versio

Critical Correlations for Short-Range Valence-Bond Wave Functions on the Square Lattice

We investigate the arguably simplest $SU(2)$-invariant wave functions capable of accounting for spin-liquid behavior, expressed in terms of nearest-neighbor valence-bond states on the square lattice and characterized by different topological invariants. While such wave-functions are known to exhibit short-range spin correlations, we perform Monte Carlo simulations and show that four-point correlations decay algebraically with an exponent $1.16(4)$. This is reminiscent of the {\it classical} dimer problem, albeit with a slower decay. Furthermore, these correlators are found to be spatially modulated according to a wave-vector related to the topological invariants. We conclude that a recently proposed spin Hamiltonian that stabilizes the here considered wave-function(s) as its (degenerate) ground-state(s) should exhibit gapped spin and gapless non-magnetic excitations.Comment: 4 pages, 5 figures. Updated versio

Phase Diagram of Interacting Bosons on the Honeycomb Lattice

We study the ground state properties of repulsively interacting bosons on the honeycomb lattice using large-scale quantum Monte Carlo simulations. In the hard-core limit the half-filled system develops long ranged diagonal order for sufficiently strong nearest-neighbor repulsion. This staggered solid melts at a first order quantum phase transition into the superfluid phase, without the presence of any intermediate supersolid phase. Within the superfluid phase, both the superfluid density and the compressibility exhibit local minima near particle- (hole-) density one quarter, while the density and the condensate fraction show inflection points in this region. Relaxing the hard-core constraint, supersolid phases emerge for soft-core bosons. The suppression of the superfluid density is found to persist for sufficiently large, finite on-site repulsion.Comment: 4 pages with 5 figure

Quantum Monte Carlo simulations of fidelity at magnetic quantum phase transitions

When a system undergoes a quantum phase transition, the ground-state wave-function shows a change of nature, which can be monitored using the fidelity concept. We introduce two Quantum Monte Carlo schemes that allow the computation of fidelity and its susceptibility for large interacting many-body systems. These methods are illustrated on a two-dimensional Heisenberg model, where fidelity estimators show marked behaviours at two successive quantum phase transitions. We also develop a scaling theory which relates the divergence of the fidelity susceptibility to the critical exponent of the correlation length. A good agreement is found with the numerical results.Comment: 4 pages, 3 figures; v2: added scaling theory; v3: published versio

Universal Reduction of Effective Coordination Number in the Quasi-One-Dimensional Ising Model

Critical temperature of quasi-one-dimensional general-spin Ising ferromagnets is investigated by means of the cluster Monte Carlo method performed on infinite-length strips, L times infty or L times L times infty. We find that in the weak interchain coupling regime the critical temperature as a function of the interchain coupling is well-described by a chain mean-field formula with a reduced effective coordination number, as the quantum Heisenberg antiferromagnets recently reported by Yasuda et al. [Phys. Rev. Lett. 94, 217201 (2005)]. It is also confirmed that the effective coordination number is independent of the spin size. We show that in the weak interchain coupling limit the effective coordination number is, irrespective of the spin size, rigorously given by the quantum critical point of a spin-1/2 transverse-field Ising model.Comment: 12 pages, 6 figures, minor modifications, final version published in Phys. Rev.