Terahertz-induced high-order harmonic generation and nonlinear charge transport in graphene

We theoretically study the terahertz-induced high-order harmonic generation (HHG) and nonlinear electric transport in graphene based on the quantum master equation with the relaxation time approximation. To obtain microscopic insight into the phenomena, we compare the results of the fully dynamical calculations with those under a quasistatic approximation, where the electronic system is approximated as a nonequilibrium steady state. As a result, we find that the THz-induced electron dynamics in graphene can be accurately modeled with the nonequilibrium steady state at each instance. The population distribution analysis further clarifies that the THz-induced HHG in graphene originates from the reduction of effective conductivity due to a large displacement of electrons in the Brillouin zone. By comparing the present nonequilibrium picture with a thermodynamic picture, we explore the role of the nonequilibrium nature of electron dynamics on the extremely nonlinear optical and transport phenomena in graphene.