We predict the existence of spontaneous spin and valley symmetry-broken states of interacting massive Dirac Fermions in a gated bilayer graphene quantum dot based on the exact diagonalization of the many-body Hamiltonian. The dot is defined by a vertical electric field and lateral gates, and its single-particle (SP) energies, wave functions, and Coulomb matrix elements are computed by using the atomistic tight-binding model. The effect of the Coulomb interaction is measured by the ratio of Coulomb elements to the SP level spacing. As we increase the interaction strength, we find the electrons in a series of spin and valley symmetry-broken phases with increasing valley and spin polarizations. The phase transitions result from the competition of the SP, exchange, and correlation energy scales. A phase diagram for N = 1-6 electrons is mapped out as a function of the Coulomb interaction strength.
Keywords: atomistic tight-binding model; bilayer graphene; configuration interaction approach; gated quantum dot; phase diagram.