Stability conditions of armchair graphene nanoribbon bipolarons

Graphene nanoribbons are 2D hexagonal lattices with semiconducting band gaps. Below a critical electric field strength, the charge transport in these materials is governed by the quasiparticle mechanism. The quasiparticles involved in the process, known as polarons and bipolarons, are self-interacting states between the system charges and local lattice distortions. To deeply understand the charge transport mechanism in graphene nanoribbons, the study of the stability conditions for quasiparticles in these materials is crucial and may guide new investigations to improve the efficiency for a next generation of graphene-based optoelectronic devices. Here, we use a two-dimensional version of the Su-Schrieffer-Heeger model to investigate the stability of bipolarons in armchair graphene nanoribbons (AGNRs). Our findings show how bipolaron stability is dependent on the strength of the electron-phonon interactions. Moreover, the results show that bipolarons are dynamically stable in AGNRs for electric field strengths lower than 3.0 mV/angstrom. Remarkably, the system's binding energy for a lattice containing a bipolaron is smaller than the formation energy of two isolated polarons, which suggests that bipolarons can be natural quasiparticle solutions in AGNRs.