Enhancement of high-order harmonic generation in graphene by mid-infrared and terahertz fields

Wenwen Mao,Angel Rubio,Shunsuke A. Sato

2023-06-30

Abstract:We theoretically investigate high-order harmonic generation (HHG) in graphene under mid-infrared (MIR) and terahertz (THz) fields based on a quantum master equation. Numerical simulations show that MIR-induced HHG in graphene can be enhanced by a factor of 10 for fifth harmonic and a factor of 25 for seventh harmonic under a THz field with a peak strength of 0.5 MV/cm by optimizing the relative angle between the MIR and THz fields. To identify the origin of this enhancement, we compare the fully dynamical calculations with a simple thermodynamic model and a nonequilibrium population model. The analysis shows that the enhancement of the high-order harmonics mainly results from a coherent coupling between MIR- and THz-induced transitions that goes beyond a simple THz-induced population contribution.

Materials Science,Mesoscale and Nanoscale Physics,Optics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to enhance the high - order harmonic generation (HHG) induced by mid - infrared (MIR) light pulses in graphene by introducing terahertz (THz) fields. Specifically, the author explored the enhancement mechanism of MIR - induced HHG in graphene in the presence of THz fields and analyzed the microscopic physical mechanism behind this enhancement.

### Research Background and Motivation

High - order harmonic generation (HHG) is an extremely nonlinear optical phenomenon involving efficient up - conversion of photons. This phenomenon was initially observed in atomic gases and has been widely studied in solid materials. Graphene is considered an efficient HHG material due to its unique electronic structure (such as Dirac cones). However, in order to further improve the HHG efficiency and develop new HHG light sources and spectroscopic techniques, researchers need to explore new methods to enhance HHG.

### Research Questions

1. **How to enhance MIR - induced HHG by introducing THz fields?**
- Through theoretical calculations and numerical simulations, the author studied the influence of THz fields on MIR - induced HHG. The results show that at the optimized relative angle, the THz field can significantly enhance the intensities of the fifth - order and seventh - order harmonics.

2. **What is the enhancement mechanism?**
- The author compared the fully dynamic calculation results with the thermodynamic model and the non - equilibrium population model and found that the enhancement of HHG mainly stems from the coherent coupling between MIR and THz fields, not just the population effect caused by the THz field.

### Main Findings

- **Enhancement Effect**:
- When the peak intensity of the THz field is 0.5 MV/cm, the intensity of the fifth - order harmonic is enhanced by 10 times, and the intensity of the seventh - order harmonic is enhanced by 25 times.

- **Enhancement Mechanism**:
- The enhancement of HHG is mainly due to the coherent coupling between MIR and THz fields, which goes beyond the simple population contribution caused by the THz field. Specifically, the THz field not only changes the population distribution of the system but also enhances HHG through coherent driving.

### Conclusion

By introducing THz fields, the high - order harmonic generation induced by MIR light pulses in graphene can be significantly enhanced. This enhancement is mainly attributed to the coherent coupling between MIR and THz fields, not just the population effect caused by the THz field. This finding provides new ideas for developing efficient HHG light sources and spectroscopic techniques.

### Formula Summary

1. **Quantum Master Equation**:
\[
\frac{d}{dt}\rho_k(t)=\frac{1}{i\hbar}\left[H_k + \frac{eA(t)}{\hbar},\rho_k(t)\right]+\hat{D}[\rho_k(t)]
\]
where $\rho_k(t)$ is the one - body reduced density matrix of Bloch wave vector $k$, $H_k+\frac{eA(t)}{\hbar}$ is the one - body Hamiltonian, and $\hat{D}[\rho_k(t)]$ is the relaxation operator.

2. **Current Density**:
\[
J(t)=\frac{2}{(2\pi)^2}\int dk\text{Tr}\left[\hat{J}_k(t)\rho_k(t)\right]
\]
where $\hat{J}_k(t)=-\frac{e}{m_e\hbar}\frac{\partial(H_k + \frac{eA(t)}{\hbar})}{\partial k}$ is the current operator.

3. **Harmonic Power Spectrum**:
\[
I_{\text{HHG}}(\omega)\sim\omega^2\left|\int_{-\infty}^{\infty}dtJ(t)e^{i\omega t}\right|^2
\]

4. **Harmonic Integral Intensity**:
\[
I_{\text{nth}}^{\text{total}}=\int_{\frac{n - 1}{2}\om