# Collective excitations and the nature of Mott transition in undoped gapped graphene

###### Abstract

Particle-hole continuum (PHC) for massive Dirac fermions in presence of short range interactions, provides an unprecedented opportunity for formation of two collective split-off states, one in the singlet and the other in the triplet (spin-1) channel in undoped system. Both poles are close in energy and are separated from the continuum of free particle-hole excitations by an energy scale of the order of gap parameter . They both disperse linearly with two different velocities reminiscent of spin-charge separation in Luttinger liquids. When the strength of Hubbard interactions is stronger than a critical value, the velocity of singlet excitation which we interpret as a charge boson composite becomes zero, and renders the system a Mott insulator. Beyond this critical point, the low-energy sector is left with a linearly dispersing triplet mode – a characteristic of a Mott insulator. The velocity of triplet mode at the Mott criticality is twice the velocity of underlying Dirac fermions. The phase transition line in the space of and is in qualitative agreement with a more involved dynamical mean field theory (DMFT) calculation.

###### pacs:

73.22.Pr, 73.20.Mf, 75.10.Kt## I Introduction

Graphene is a huge network of bonds formed between carbon atoms on a two-dimensional (2D) sheet Novoselov (). The low-energy effective model for electrons in this system is the 2+1 dimensional Dirac theory which contains an energy scale proportional to – the Fermi velocity of electrons NetoRMP (). Effects of interactions and/or disorder can be taken into account on top of this non-interacting fixed point NetoRMP (). Breaking the sub-lattice symmetry gives rise to a gap in the single-particle excitations spectrum Lanzara (), which introduces the gap parameter as another energy scale. The gapped graphene can be modeled with 2+1 dimensional massive Dirac fermions Alireza (). Gapped graphene supports interesting new class of excitations in the spectrum such as, domain walls Semenoff (). Moreover the single-particle gap (mass) gives rise to suppression of Coulomb fields at nano-meter length scales Pereira (). The gap could also give rise to universal nonlinear optical response nloJafari (). When the many-body interactions of various form are added to the massive Dirac theory, the nature of many-body excitations becomes even more interesting. Recent ab-initio estimates of the strength of the Coulomb repulsion in various forms of graphene suggests that the Hubbard parameter in graphene can be quite remarkable Wehling (), which introduces a third energy scale, . Therefore, an interesting question here would be the effects of local Hubbard type interactions on the excitation spectrum of gapped graphene, and the nature of transition to Mott insulating phase in gapped graphene. Our recent full-fledged DMFT investigation of the so called ionic-Hubbard model on the honeycomb lattice suggested three phases moradDMFT (): (i) Band insulating phase for . (ii) Mott insulating phase for , and (iii) a semi-metallic phase for . The above study is about the nature of ground state. In order to study the connection of excited states and the ground state, one notes that; deep in the Mott insulating phase, the energy scale is expected to be dominant over , such that the low-energy excitations are expected to be spin fluctuations arising from super-exchange interaction induced by the large . We would like to focus on, the evolution of the collective spin and charge dynamics in this system as a function of the Hubbard for a system with a non-zero gap parameter . It turns out that the presence of a non-zero gap parameter facilitates the separation of two collective states each of which has a different velocity, reminiscent of the spin-charge separation in Luttinger liquids.

## Ii Formulation of the problem

The Hamiltonian we consider here is the so called ionic-Hubbard model, which is defined by

(1) |

where with corresponding to number operator at the ’th unit cell at sub-lattices, A, B, respectively. The spinor notation has been used. Here denotes the component of the physical spin, is a two-dimensional momentum vector, and stands for Pauli matrices. is the strength of the on-site Coulomb interaction known as Hubbard parameter, and is the mass parameter. The quadratic part of this Hamiltonian can be diagonalized by a simple unitary transformation to the basis of conduction () and valence () states, and the corresponding eigen-values are given by,

(2) |

where and being the C-C bond length. The hopping amplitude between the nearest neighboring carbon atoms is eV. Due to paramagnetic nature of the non-interacting system (i.e. ), the occupation numbers and energies of the conduction and valence bands are independent of spin orientation. We consider the undoped system with precisely one electron per carbon atom, and assume the temperature to be zero. Then we extend the collective mode analysis of Ref. JafariBaskaranJPCM () to the case with , whereby the single-particle spectrum changes from to the one given by Eq. (2). This change in the single-particle spectrum, changes the borders of the particle-hole continuum. However, as will be show here, at the border corresponding to massless spectrum, interesting collective quanta can be made in presence of short range Coulomb interactions.

Using equation of motion it can be shown that a simple RPA-like expression governs the eigen-value equation for the collective excitations in the triplet and singlet channels of the two-band arising from Dirac fermions in presence of the short range Coulomb interactions are given by,

(3) |

where the () sign corresponds to triplet (singlet) channel BickersScalapino (), and the sign difference can be traced back to fermion anti-commutation relation. Note that although in the above expression, is the particle-hole (polarization) bubble, the (approximate) bosonic operators which satisfy the above equation are of a peculiar form, which does not necessarily coincide with the precise particle-hole fluctuation form JafariBaskaranJPCM (). The non-interacting susceptibility employed in the above equation, is defined by,

(4) |

For the dispersion relation (2), the susceptibility has been calculated by many authors, the analytic form of which reads Pyatkovskiy (); Kotov ():

(5) |

To understand the structure of this function, in Fig. 1 we have plotted the real and imaginary part of this function versus for some representative value of . As can be seen, the imaginary part is non-zero for , where

(6) |

corresponds to the lower boundary of PHC of massive Dirac fermions. The singularity of the above inter-band response at is of a weak logarithmic form which could give rise to a solution for Eq. (3) only in the triplet channel, and for extremely large values of which is unphysical. The however possesses another much more interesting simple pole structure at a lower energy scale, given by,

(7) |

Note that this equation defines a line in the plane of and , which corresponds to the lower boundary of the PHC of massless () Dirac fermions and is determined by the characteristic velocity, of the underlying Dirac fermions. This could be considered as a property of a ”parent” massless Dirac system which manifests itself as a simple-pole structure when a mass parameter, is turned on in the non-interacting part of the Hamiltonian. At the gapless point, where , the two energy scales at and merge and the resulting singularity of will be of inverse-square-root form BaskaranJafari ().

As can be seen, to search for solutions of Eq. (3), one has to look for the intersection of the real part (green line) with the constant horizontal line , where the upper (lower) sign corresponds to singlet (triplet) channel. As can be very clearly seen in the figure, near , where the imaginary part of is zero, there can be poles both in singlet and triplet channel. The simple pole structure of the particle-hole propagator near this energy scale implies that the poles exist for every value of the Hubbard , no matter how small or how large it is. To see this more clearly, note that the real part of the susceptibility near behaves as,

(8) |

The solutions of the eigenvalue equations (3) for arbitrary value of define two collective branches in singlet and triplet channel as,

(9) |

where upper sign corresponds to the singlet channel, and the lower sign corresponds to the triplet channel. These solutions are available even for arbitrarily small of . This formula can be interpreted as the spin charge separation in the sense that the singlet and the spin-1 modes move at different velocities given by,

(10) |

The mere existence and simple pole structure of the at the energy scale is a unique consequence of the gap opening in the single-particle spectrum of excitations. The presence of gap in the single particle spectrum, pushes the continuum of free particle-hole pairs to higher energies, and hence the above collective excitations which appear around line are well separated from the boundary of PHC and therefore will be protected from Landau damping to the incoherent background of free particle-hole excitations.

## Iii Discussion

The opening of a single particle gap in the spectrum of excitations in graphene, although pushes the lower edge of the particle-hole continuum from the energy scale to energy scale , but still leaves behind a signature of underlying massless Dirac fermions in the form of a singular behavior for on the line defined by corresponding to the PHC boundary of the underlying massless Dirac fermions. The simple pole left on this line has two consequences: (i) Due to the sign change of the susceptibility across , there will be collective mode solutions at both singlet and triplet channels, for arbitrary values of the Hubbard . The explicit form of the solutions are given by Eq. (9). (ii) The different velocities for the two modes is reminiscent of the situation one encounters in one-dimension, where a faster ”spin mode” overtakes the slower ”charge mode”. Albeit the difference is that in 1D the divergences in particle-hole bubbles are of the characteristic inverse square root type. In the case of 2D massive Dirac fermions, where is finite, the divergence in the particle-hole propagator is of a simple pole form. When the limit is taken, the simple pole merges with the logarithmic singularity at the boundary of PHC, , and gives rise again to a inverse-square-root behavior for gapless Dirac fermions in 2D BaskaranJafari (). It is therefore the non-zero value of which provides a solution in the singlet channel in addition to the known solution of the triplet channel at point BaskaranJafari (); JafariBaskaranJPCM (). In this sense, it appears that the non-zero value of enhances the separation of the triplet and the singlet modes. Therefore the smallest value of the gap parameter, , totally changes the nature of collective excitations in graphene. This observation may have implications for novel approaches to the bosonization of the 2+1 dimensional massive Dirac fermions. Note that if the interaction employed was of a long-range Coulomb form, the solution in the spin-1 channel would be of a linearly dispersing gapped form, and the solution in the singlet channel would not exist at all. Therefore the two modes discussed here are peculiar feature of short range interactions, which due to remarkable value of the Hubbard in graphene, may have relevance to physically fabricated samples of gapped graphene.

Now let us discuss the nature of the phase transition marked by vanishing of the velocity of the singlet mode. As can be seen in Eq. (10), when the Hubbard is large-enough to satisfy,

(11) |

the velocity of singlet excitations becomes zero, and the low-energy sector is exhausted by only triplet excitations of the form

(12) |

Thinking in the spirit of slave-boson approach Vaezi (), an electron (hole) can be assumed to be composed of its spin part, the spinon, and the charge part the doublon (holon). Therefore a composite object constructed from an electron and a hole can be thought of a separate spin-1 composite of two spinons, and a charge composite of a doublon and a holon giving rise to a spin zero boson. Hence the triplet mode can be interpreted as a triplet bound state of two spinons BaskaranJafari (); JafariBaskaranJPCM (), and the new singlet mode the emergence of which is being facilitated by a non-zero gap in the single-particle spectrum, can be interpreted as a composite boson constructed from a doublon and a holon the total charge of which is zero, as it should be. In this sense, the vanishing of the velocity of singlet mode can be associated with infinite enhancement of the effective mass of charge bosons. This gives rise to localization of charge carriers and hence is a Mott transition Jafari2009 (). Indeed the decreasing trend as a function of the gap parameter in Eq. (11) is in qualitative agreement with our previous DMFT result for the phase boundary of the Mott insulating phase in the ionic Hubbard model moradDMFT (). Note that this agreement holds for very small to moderate values of . For large values of , the phase boundary in DMFT will be given by moradDMFT (). However, as long as experimental realization of gap parameter in graphene is concerned, will be on the scale of few tens of meV, which is quite small in the scale of the hopping term eV. Moreover, the fact that the only low-energy excitations left in the system when goes beyond is a another evidence that at the systems becomes a Mott insulator. Because in the Mott insulator, the only possible low-energy modes are spin excitations.

When is further increased, it can be seen in Fig. 1 that the real part of and (triplet channel) will intersect in another point which is immediately below the edge of the PHC, , and therefore a second triplet mode at slightly higher energy than the first one is expected to form, which is anticipated due to weak logarithmic singularity at . However, this mode is not likely to be realized as it requires quite a large value of which falls already in the Mott insulating phase, where the ground state maybe totally deformed with respect to the initial starting point of massive Dirac fermions employed in this work.

It is interesting to note that, although the velocity of the underlying Dirac fermions is , the velocity of the spin-1 mode is always more than . Specially at the transition to Mott insulating phase, the velocity of the spin-1 mode will be , i.e. the triplet mode (if we call it triplon) moves twice faster than the underlying Dirac electrons. This prediction can be tested not only in gapped graphene system, but also on a platform based on cold atoms engineered to mimic a massive Dirac theory at nano-Kelvin energy scales. The existence of a mode whose velocity can be twice the velocity of underlying Dirac fermions may have interesting implications in teleportation of quantum information, as well as in the spin-only forms of transport.

The above two modes discussed here exist for small values of too. Therefore it should be possible to check for their experimental consequences. The well isolation of the energy scale from the PHC edge makes the gapped graphene more interesting for performing neutron scattering experiments. In this work we find that the triplet excitation exists below the Mott transition. Moreover, the Mott insulating phase has its own magnetic excitations. Therefore in the gapped graphene we anticipate a triplet excitation over a large range of values of the Hubbard parameter . Such low-energy spin-1 branch of excitations will have a characteristic contribution in the specific heat at constant volume, where is the velocity of spin modes. Below the Mott transition, there will be another contribution coming from singlet charge modes. The influence of these modes on various properties of gapped graphene remains to be investigated.

## Iv acknowledgments

This work was supported by the National Elite Foundation (NEF) of Iran. We are grateful to K. Yamada for hospitality at the Institute for Materials Research, Tohoku University, and the suggestion of the present calculation. We thank S. Rouhani and A. Vaezi for useful discussions.

## References

- (1) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306 (2004) 666; K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, A. A. Firsov, Nature 438 (2005) 197.
- (2) For a review see: A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, A. K. Geim, Rev. Mod. Phys. 81 (2009) 109.
- (3) S. Y. Zhou, G.-H. Gweon, A. V. Fedorov, P. N. First, W. A. de Heer, D.-H. Lee, F. Guinea, A. H. Castro Neto, and A. Lanzara Nature Materials 6 (2007) 770.
- (4) Alireza Qaiumzadeh, and Reza Asgari, Phys. Rev. B 79 (2009) 075414.
- (5) G. W. Semenoff, V. Semenoff, and Fei Zhou, Phys. Rev. Lett. 101 (2008) 087204.
- (6) Vitor M. Pereira, Valeri N. Kotov, and A. H. Castro Neto, Phys. Rev. B 78 (2008) 085101. Valeri N. Kotov, Vitor M. Pereira, and Bruno Uchoa, Phys. Rev. B 78 (2008) 075433 .
- (7) S. A. Jafari, arxiv:1202.1070 (2012).
- (8) T. O. Wehling, E. Sasioglu, C. Friedrich, A. I. Lichtenstein, M. I. Katsnelson, and S. Blügel, Phys. Rev. Lett. 106 (2011) 236805.
- (9) M. Ebrahimkhas, S. A. Jafari, arxiv:1112.3616 (2011).
- (10) G. Baskaran, S.A. Jafari, Phys. Rev. Lett. 89 (2002) 016402.
- (11) S. A. Jafari, G. Baskaran, J. Phys. Condens. Matter 24 (2012) 09560.
- (12) N. E. Bickers, D. J. Scalapino, Ann. Phys. 193 (1989) 206.
- (13) P. K. Pyatkovskiy, J. Phys. Codens. Matter 21 (2009) 025506.
- (14) V. N. Kotov, V. M. Pereira, B. Uchoa, Phys. Rev. B 78 (2008) 075433.
- (15) A. Vaezi, X. G. Wen, arxiv:1010.5744 (2010).
- (16) S. A. Jafari, Eur. Phys. Jour. B 68, 537 (2009).