中国物理B ›› 2020, Vol. 29 ›› Issue (10): 104206-.doi: 10.1088/1674-1056/abab76

• • 上一篇    下一篇

Zhong Guan(管仲)1, Lu Liu(刘璐)2, Guo-Li Wang(王国利)1,†(), Song-Feng Zhao(赵松峰)1, Zhi-Hong Jiao(焦志宏)1, Xiao-Xin Zhou(周效信)1,‡()   

  • 收稿日期:2020-04-26 修回日期:2020-06-15 接受日期:2020-08-01 出版日期:2020-10-05 发布日期:2020-10-05
  • 通讯作者: Guo-Li Wang(王国利), Xiao-Xin Zhou(周效信)

An improved method for the investigation of high-order harmonic generation from graphene

Zhong Guan(管仲)1, Lu Liu(刘璐)2, Guo-Li Wang(王国利)1,†, Song-Feng Zhao(赵松峰)1, Zhi-Hong Jiao(焦志宏)1, and Xiao-Xin Zhou(周效信)1,‡   

  1. 1 College of Physics and Electronic Engineering, Northwest Normal University, Lanzhou 730070, China
    2 Department of Physics, College of Science, National University of Defense Technology, Changsha 410073, China
  • Received:2020-04-26 Revised:2020-06-15 Accepted:2020-08-01 Online:2020-10-05 Published:2020-10-05
  • Contact: Corresponding author. E-mail: wanggl@nwnu.edu.cn Corresponding author. E-mail: zhouxx@nwnu.edu.cn
  • About author:
    †Corresponding author. E-mail: wanggl@nwnu.edu.cn
    ‡Corresponding author. E-mail: zhouxx@nwnu.edu.cn
    * Project supported by the National Natural Science Foundation of China (Grant Nos. 11764038, 11864037, 11765018, and 11664035) and the Science Foundation of Northwest Normal University, China (Grant No. NWNU-LKQN-17-1).

Abstract:

High-order harmonic generation (HHG) of bulk crystals in strong laser field is typically investigated with semiconductor Bloch equations (SBEs). However, in the length gauge, it suffers from the divergence for the crystals with a zero band gap, such as graphene, using both Bloch- and Houston-states expansion methods. Here, we present a method of solving the SBEs based on time-dependent Bloch basis, which is equivalent to semiconductor Bloch equations in the velocity gauge. Using this method, we investigate the HHG of a single-layer graphene. It is found that our results for population are in good agreement with the other results. For a initial condition py = 0, we find the electrons just move in single valence band or conduction band, which are in accord with classical results. Our simulations on the HHG dependence of polarization of driving laser pulse confirm that 5th, 7th, and 9th harmonic yields increase to the maximal value when laser ellipticity ε ≈ 0.3. What is more, similar to the case of atoms in the laser field, the total strength of 3rd harmonic decrease monotonically with the increase of ε. In addition, we simulate the dependence of HHG on crystallographic orientation with respect to the polarization direction of linear mid-infrared laser pulse, and the results reveal that for higher harmonics, their radiation along with the change of rotation angle θ reflects exactly the sixfold symmetry of graphene. Our method can be further used to investigate the behaviors of other materials having Dirac points (i.e., surface states of topological insulators) in the strong laser fields.

Key words: high-order harmonic generation, graphene, velocity gauge, divergence

中图分类号:  (Frequency conversion; harmonic generation, including higher-order harmonic generation)

  • 42.65.Ky
42.50.Hz (Strong-field excitation of optical transitions in quantum systems; multiphoton processes; dynamic Stark shift) 42.65.Re (Ultrafast processes; optical pulse generation and pulse compression) 72.20.Ht (High-field and nonlinear effects)