In this paper, we study the massive Dirac equation with the presence of the Morse potential in polar coordinate. The Dirac Hamiltonian is written as two second-order differential equations in terms of two spinor wavefunctions. Since the motion of electrons in graphene is propagated like relativistic fermionic quasi-particles, then one is considered only with pseudospin symmetry for aligned spin and unaligned spin by arbitrary 𝑘. Next, we use the confluent Heun’s function for calculating the wavefunctions and the eigenvalues. Then, the corresponding energy spectrum obtains in terms of 𝑁 and 𝑘. Afterward, we plot the graphs of the energy spectrum and the wavefunctions in terms of 𝑘 and 𝑟, respectively. Moreover, we investigate the graphene band structure by a linear dispersion relation which creates an energy gap in the Dirac points called gapped graphene. Finally, we plot the graph of the valence and conduction bands in terms of wavevectors.
