Competing pairing channels in the doped honeycomb lattice Hubbard model

Proposals for superconductivity emerging from correlated electrons in the doped Hubbard model on the honeycomb lattice range from chiral $d+id$ singlet to $p+ip$ triplet pairing, depending on the considered range of doping and interaction strength, as well as the approach used to analyze the pairing instabilities. Here, we consider these scenarios using large-scale dynamic cluster approximation (DCA) calculations to examine the evolution in the leading pairing symmetry from weak to intermediate coupling strength. These calculations focus on doping levels around the van Hove singularity (VHS) and are performed using DCA simulations with an interaction-expansion continuous-time quantum Monte Carlo cluster solver. We calculated explicitly the temperature dependence of different uniform superconducting pairing susceptibilities and found a consistent picture emerging upon gradually increasing the cluster size: while at weak coupling the $d+id$ singlet pairing dominates close to the VHS filling, an enhanced tendency towards $p$-wave triplet pairing upon further increasing the interaction strength is observed. The relevance of these systematic results for existing proposals and ongoing pursuits of odd-parity topological superconductivity are also discussed.Comment: 7 pages, 5 figure

Self-Learning Determinantal Quantum Monte Carlo Method

Self-learning Monte Carlo method [arXiv:1610.03137, 1611.09364] is a powerful general-purpose numerical method recently introduced to simulate many-body systems. In this work, we implement this method in the framework of determinantal quantum Monte Carlo simulation of interacting fermion systems. Guided by a self-learned bosonic effective action, our method uses a cumulative update [arXiv:1611.09364] algorithm to sample auxiliary field configurations quickly and efficiently. We demonstrate that self-learning determinantal Monte Carlo method can reduce the auto-correlation time to as short as one near a critical point, leading to $\mathcal{O}(N)$-fold speedup. This enables to simulate interacting fermion system on a $100\times 100$ lattice for the first time, and obtain critical exponents with high accuracy.Comment: 5 pages, 4 figure

Itinerant quantum critical point with frustration and non-Fermi-liquid

Employing the self-learning quantum Monte Carlo algorithm, we investigate the frustrated transverse-field triangle-lattice Ising model coupled to a Fermi surface. Without fermions, the spin degrees of freedom undergoes a second-order quantum phase transition between paramagnetic and clock-ordered phases. This quantum critical point (QCP) has an emergent U(1) symmetry and thus belongs to the (2+1)D XY universality class. In the presence of fermions, spin fluctuations introduce effective interactions among fermions and distort the bare Fermi surface towards an interacting one with hot spots and Fermi pockets. Near the QCP, non-Fermi-liquid behavior are observed at the hot spots, and the QCP is rendered into a different universality with Hertz-Millis type exponents. The detailed properties of this QCP and possibly related experimental systems are also discussed.Comment: 9 pages, 8 figure