Electronic response of graphene to an ultrashort intense terahertz radiation pulse

We have recently reported a study (Ishikawa 2010 Phys. Rev. B 82 201402) on a nonlinear optical response of graphene to a normally incident terahertz radiation pulse within the massless Dirac fermion (MDF) picture, where we have derived physically transparent graphene Bloch equations (GBE). Here we extend it to the tight-binding (TB) model and oblique incidence. The derived equations indicate that interband transitions are governed by the temporal variation of the spinor phase along the electron path in the momentum space and predominantly take place when the electron passes near the Dirac point. At normal incidence, the equations for electron dynamics within the TB model can be cast into the same form of GBE as for the MDF model. At oblique incidence, the equations automatically incorporate photon drag and satisfy the continuity equation for electron density. Single-electron dynamics strongly depend on the model and pulse parameters, but the rapid variations are averaged out after momentum-space integration. Direct current remaining after the pulse is generated in graphene irradiated by an intense monocycle terahertz pulse, even if it is linearly polarized and normally incident. The generated current depends on the carrier-envelope phase, pulse intensity and Fermi energy in a complex manner.