Abstract

Recent experiments on graphene multilayers under displacement field have demonstrated a wide variety of electronically ordered phases, including valley and/or spin polarized phases as well as potentially unconventional superconducting phases. In addition, quantum oscillation measurements in Bernal bilayer graphene and rhombohedral trilayer graphene showed the presence of electronic nematic states. Here, we investigate the emergence of nematic order in these systems, with emphasis on Bernal bilayer graphene, within a self-consistent Hartree–Fock framework, using a realistic band structure with trigonal warping and a dual-gate screened Coulomb interaction. We compute the phase diagram as carrier density and displacement field are varied, and find a sequence of isospin (spin and valley) polarized states consistent with experiment, including partially polarized phases with large and small Fermi pockets. Within these partially isospin polarized states, we identify regions with three C₃-symmetric pockets for the minority carriers and regions with only one or two such pockets, implying that the system develops a spontaneous nematic order that breaks C₃ symmetry. We find numerically that nematicity emerges near the boundary between fully and partially polarized phases and is controlled by the strength of the screened interaction. We analytically derive a criterion for the nematic order, which agrees with our numerical results.

Abstract

Recent experiments on graphene multilayers under displacement field have demonstrated a wide variety of electronically ordered phases, including valley and/or spin polarized phases as well as potentially unconventional superconducting phases. In addition, quantum oscillation measurements in Bernal bilayer graphene and rhombohedral trilayer graphene showed the presence of electronic nematic states. Here, we investigate the emergence of nematic order in these systems, with emphasis on Bernal bilayer graphene, within a self-consistent Hartree–Fock framework, using a realistic band structure with trigonal warping and a dual-gate screened Coulomb interaction. We compute the phase diagram as carrier density and displacement field are varied, and find a sequence of isospin (spin and valley) polarized states consistent with experiment, including partially polarized phases with large and small Fermi pockets. Within these partially isospin polarized states, we identify regions with three C₃-symmetric pockets for the minority carriers and regions with only one or two such pockets, implying that the system develops a spontaneous nematic order that breaks C₃ symmetry. We find numerically that nematicity emerges near the boundary between fully and partially polarized phases and is controlled by the strength of the screened interaction. We analytically derive a criterion for the nematic order, which agrees with our numerical results.

Abstract

Recent experiments on graphene multilayers under displacement field have demonstrated a wide variety of electronically ordered phases, including valley and/or spin polarized phases as well as potentially unconventional superconducting phases. In addition, quantum oscillation measurements in Bernal bilayer graphene and rhombohedral trilayer graphene showed the presence of electronic nematic states. Here, we investigate the emergence of nematic order in these systems, with emphasis on Bernal bilayer graphene, within a self-consistent Hartree–Fock framework, using a realistic band structure with trigonal warping and a dual-gate screened Coulomb interaction. We compute the phase diagram as carrier density and displacement field are varied, and find a sequence of isospin (spin and valley) polarized states consistent with experiment, including partially polarized phases with large and small Fermi pockets. Within these partially isospin polarized states, we identify regions with three C₃-symmetric pockets for the minority carriers and regions with only one or two such pockets, implying that the system develops a spontaneous nematic order that breaks C₃ symmetry. We find numerically that nematicity emerges near the boundary between fully and partially polarized phases and is controlled by the strength of the screened interaction. We analytically derive a criterion for the nematic order, which agrees with our numerical results.

Abstract

Recent experiments on graphene multilayers under displacement field have demonstrated a wide variety of electronically ordered phases, including valley and/or spin polarized phases as well as potentially unconventional superconducting phases. In addition, quantum oscillation measurements in Bernal bilayer graphene and rhombohedral trilayer graphene showed the presence of electronic nematic states. Here, we investigate the emergence of nematic order in these systems, with emphasis on Bernal bilayer graphene, within a self-consistent Hartree–Fock framework, using a realistic band structure with trigonal warping and a dual-gate screened Coulomb interaction. We compute the phase diagram as carrier density and displacement field are varied, and find a sequence of isospin (spin and valley) polarized states consistent with experiment, including partially polarized phases with large and small Fermi pockets. Within these partially isospin polarized states, we identify regions with three C₃-symmetric pockets for the minority carriers and regions with only one or two such pockets, implying that the system develops a spontaneous nematic order that breaks C₃ symmetry. We find numerically that nematicity emerges near the boundary between fully and partially polarized phases and is controlled by the strength of the screened interaction. We analytically derive a criterion for the nematic order, which agrees with our numerical results.