The emergence of superconductivity in doped Mott insulators has been debated for decades. In this paper, we report the theoretical discovery of a time-reversal symmetry breaking superconducting ground state in the doped Mott insulator (described by the well known $t$-$J$ model) on honeycomb lattice, based on a recently developed variational method: the Grassmann tensor product state approach. As a benchmark, we use exact diagonalization and density-matrix renormalization methods to check our results on small clusters. We find systematic consistency for the ground-state energy as well as other physical quantities, such as the staggered magnetization. At low doping, the superconductivity coexists with antiferromagnetic ordering.

• Received 9 February 2012

DOI:https://doi.org/10.1103/PhysRevB.88.155112