Numerical simulations of continuum field theories and exotic quantum phase transitions

Abstract

In this thesis, we investigate several topics pertaining to emergent collective quantum phenomena in the domain of correlated fermions, using the quantum Monte Carlo method. They display exotic low temperature phases as well as phase transitions which are beyond the Landau–Ginzburg theory. The interplay between three key points is crucial for us: fermion statistics, many body effects and topology. We highlight the following several achievements: 1. Successful modeling of continuum field theories with lattice Hamiltonians, 2. their sign-problem-free Monte Carlo simulations of these models, 3. and numerical results beyond mean field descriptions. First, we consider a model of Dirac fermions with a spin rotational invariant inter- action term that dynamically generates a quantum spin Hall insulator. Surprisingly, an s-wave superconducting phase emerges due to the condensation of topological de- fects of the spin Hall order parameter. When particle-hole symmetry is present, the phase transition between the topological insulator and the superconducting phase is an example of a deconfined quantum critical point(DQCP). Although its low energy effec- tive field theory is purely bosonic, the exact conservation law of the skyrmion number operator rules out the possibility of realizing this critical point in lattice boson models. This work is published in Ref. [1]. Second, we dope the dynamically generated quantum spin Hall insulator mentioned above. Hence it is described by a field theory without Lorentz invariance due to the lack of particle-hole symmetry. This sheds light on the extremely hot topic of twisted bilayergraphene: Why is superconductivity generated when the repulsive Coulomb interaction is much stronger than the electron-phonon coupling energy scale? In our case, Cooper pairs come from the topological skyrmion defects of the spin current order parameter, which are charged. Remarkably, the nature of the phase transition is highly non-mean-field-like: one is not allowed to simply view pairs of electrons as single bosons in a superfluid-Mott insulator transition, since the spin-current order parameter can not be ignored. Again, due to the aforementioned skyrmions, the two order parameters are intertwined: One phase transition occurs between the two symmetry breaking states. This work is summarized in Ref. [2]. Third, we investigate the 2 + 1 dimensional O(5) nonlinear sigma model with a topological Wess-Zumino-Witten term. Remarkably, we are able to perform Monte Carlo calculations with a UV cutoff given by the Dirac Landau level quantization. It is a successful example of simulating a continuous field theory without lattice regularization which leads to an additional symmetry breaking. The Dirac background and the five anti-commuting Dirac mass terms naturally introduce the picture of a non-trivial Berry phase contribution in the parameter space of the five component order parameter. Using the finite size scaling method given by the flux quantization, we find a stable critical phase in the low stiffness region of the sigma model. This is a candidate ground state of DQCP when the O(5) symmetry breaking terms are irrelevant at the critical point. Again, it has a bosonic low energy field theory which is seemingly unable to be realized in pure boson Hamiltonians. This work is summarized in Ref. [3].