The discovery of quantum spin-Hall (QSH) insulators has brought topology to the forefront of condensed matter physics. While a QSH state from spin-orbit coupling can be fully understood in terms of band theory, fascinating many-body effects are expected if it instead results from spontaneous symmetry breaking. Here, we introduce a model of interacting Dirac fermions where a QSH state is dynamically generated. Our tuning parameter further allows us to destabilize the QSH state in favour of a superconducting state through proliferation of charge-2e topological defects. This route to superconductivity put forward by Grover and Senthil is an instance of a deconfined quantum critical point (DQCP). Our model offers the possibility to study DQCPs without a second length scale associated with the reduced symmetry between field theory and lattice realization and, by construction, is amenable to large-scale fermion quantum Monte Carlo simulations.

No abstract

We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.

No abstract

No abstract