Antiferromagnetism Induced by Bond Su-Schrieffer-Heeger Electron-Phonon Coupling: A Quantum Monte Carlo Study

Non-Hermitian Strongly Interacting Dirac Fermions

Exotic quantum phases and phase transition in the strongly interacting Dirac systems have attracted tremendous interests. On the other hand, non-Hermitian physics, usually associated with dissipation arising from the coupling to environment, emerges as a frontier of modern physics in recent years. In this Letter, we investigate the interplay between non-Hermitian physics and strong correlation in Dirac-fermion systems. We generalize the projector quantum Monte-Carlo (PQMC) algorithm to the non-Hermitian interacting fermionic systems. Employing PQMC simulation, we decipher the ground-state phase diagram of the honeycomb Hubbard model with spin resolved non-Hermitian asymmetric hopping processes. The antiferromagnetic (AFM) ordering induced by Hubbard interaction is enhanced by the non-Hermitian asymmetric hopping. Combining PQMC simulation and renormalization group analysis, we reveal that the quantum phase transition between Dirac semi-metal and AFM phases belongs to Hermitian chiral XY universality class, implying that a Hermitian Gross-Neveu transition is emergent at the quantum critical point although the model is non-Hermitian.

Emergence of Competing Orders and Possible Quantum Spin Liquid in SU(N) Fermions

In the past few decades, tremendous efforts have been made toward understanding the exotic physics emerging from competition between various ordering tendencies in strongly correlated systems. Employing state-of-the-art quantum Monte Carlo simulation, we investigate an interacting SU(N) fermionic model with varying interaction strength and value of N, and we unveil the ground-state phase diagram of the model exhibiting a plethora of exotic phases. For small values of N-namely, N=2, 3-the ground state is an antiferromagnetic (AFM) phase, whereas in the large-N limit, a staggered valence bond solid (VBS) order is dominant. For intermediate values of N such as N=4, 5, remarkably, our study reveals that distinct VBS orders appear in the weak and strong coupling regimes. More fantastically, the competition between staggered and columnar VBS ordering tendencies gives rise to a Mott insulating phase without spontaneous symmetry breaking (SSB), existing in a large interacting parameter regime, which is consistent with a gapped quantum spin liquid. Our study not only provides a platform to investigate the fundamental physics of quantum many-body systems-it also offers a novel route toward searching for exotic states of matter such as quantum spin liquid in realistic quantum materials.

Resonating Valence Bond States in an Electron-Phonon System

We study a simple electron-phonon model on square and triangular versions of the Lieb lattice using an asymptotically exact strong coupling analysis. At zero temperature and electron density n=1 (one electron per unit cell), for various ranges of parameters in the model, we exploit a mapping to the quantum dimer model to establish the existence of a spin-liquid phase with Z_{2} topological order (on the triangular lattice) and a multicritical line corresponding to a quantum critical spin liquid (on the square lattice). In the remaining part of the phase diagram, we find a host of charge-density-wave phases (valence-bond solids), a conventional s-wave superconducting phase, and with the addition of a small Hubbard U to tip the balance, a phonon-induced d-wave superconducting phase. Under a special condition, we find a hidden pseudospin SU(2) symmetry that implies an exact constraint on the superconducting order parameters.

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and we report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but it incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

Real-time evolution of static electron-phonon models in time-dependent electric fields

We present an exact Monte Carlo method to simulate the nonequilibrium dynamics of electron-phonon models in the adiabatic limit of zero phonon frequency. The classical nature of the phonons allows us to sample the equilibrium phonon distribution and efficiently evolve the electronic subsystem in a time-dependent electromagnetic field for each phonon configuration. We demonstrate that our approach is particularly useful for charge-density-wave systems experiencing pulsed electric fields, as they appear in pump-probe experiments. For the half-filled Holstein model in one and two dimensions, we calculate the out-of-equilibrium response of the current and the energy after a pulse is applied as well as the photoemission spectrum before and after the pump. Finite-size effects are under control for chains of 162 sites (in one dimension) or 16×16 square lattices (in two dimensions).