Dynamical properties of quasiparticles in a tunable Kekulé graphene superlattice

Xiao-Yu Xiong,Xi-Dan Hu,Qizhong Zhu,Zhi Li

DOI: https://doi.org/10.1016/j.physleta.2023.128957

2023-03-20

Abstract:We investigate the dynamical properties of quasiparticles in graphene superlattices with three typical Kekulé distortions (i.e., Kekulé-O, Kekulé-Y and Kekulé-M). On the one hand, we numerically show the visualized evolution process of Kekulé quasiparticles; while on the other hand, we analytically obtain the centroid trajectory of the quasiparticles, and both of them agree well with each other. The results reveal that the relativistic Zitterbewegung (ZB) phenomenon occurs in the Kekulé systems. Furthermore, through analyzing the frequency of ZB, we unveil the one-to-one relationship between ZB and Kekulé textures, i.e., the ZB frequenies of Kekulé-O, Kekulé-Y and Kekulé-M quasiparticles feature single, double and six frequencies, respectively. Finally, we propose a scheme to distinguish among different Kekulé textures from the dynamical perspective. The predictions in this paper are expected to be experimentally verified in the near future, so as to facilitate further research of Kekulé structures in solid materials or artificial systems.

Mesoscale and Nanoscale Physics,Quantum Physics

What problem does this paper attempt to address?

This paper aims to explore the dynamic characteristics of quasiparticles in graphene superlattices with three typical Kekulé distortions, namely Kekulé - O, Kekulé - Y and Kekulé - M. Specifically, the researchers demonstrated the evolution process of Kekulé quasiparticles through numerical simulation and obtained the expression of the centroid trajectory of quasiparticles by analytical methods. The results of the two are in good agreement. The study reveals the existence of the relativistic Zitterbewegung (ZB, zero - point oscillation) phenomenon in these systems. By analyzing the ZB frequency, a one - to - one correspondence between the ZB frequency and the Kekulé texture was discovered. For example, the ZB frequencies of Kekulé - O, Kekulé - Y and Kekulé - M quasiparticles are respectively manifested as single - frequency, double - frequency and six - frequency.

In addition, the study also proposed a scheme to distinguish different Kekulé textures from a dynamic perspective. The predictions in the paper are expected to be verified experimentally in the future, thereby providing support for further research on Kekulé structures in solid materials or artificial systems.

### Main problems:
1. **Research background**: Kekulé graphene is a superlattice material formed by periodically manipulating the carbon - carbon bond density wave. Due to its Brillouin zone folding, the Kekulé system can induce quasiparticles with chiral symmetry breaking.
2. **Research objective**: To explore the dynamic characteristics of quasiparticles in Kekulé graphene superlattices with different structures, especially the Zitterbewegung phenomenon and its frequency characteristics.
3. **Research method**: By combining numerical simulation and analytical derivation, the researchers obtained the expression of the centroid trajectory of quasiparticles and further explored the dynamic characteristics of different Kekulé systems through Fourier analysis.

### Key formulas:
1. **Hamiltonian of the Kekulé system**:
\[
H = \begin{pmatrix}
m_0 v^2 & 0 & v_0 k_- & -\eta \Delta Q^* \\
\nu, + 2\Delta & m_0 v^2 & 0 & v_0 k_+ \\
-m_0 v^2 & 0 & 0 & \eta \Delta Q^* \\
\nu, - & 2\Delta & m_0 v^2 & v_0 k_- \\
\end{pmatrix}
\]
where \( k_\pm = k_x \pm i k_y \), \( v_0 \) and \( m_0 \) are the Fermi velocity and effective mass in the original graphene system respectively, \( J \) is the hopping intensity, and \( \Delta \) is the Kekulé coupling strength.

2. **Energy bands of different Kekulé structures**:
- Kekulé - O:
\[
E_{\text{Kek - O}} = \pm \sqrt{9 \Delta^2 J^2 + v_0^2 k^2}
\]
- Kekulé - Y:
\[
E_{\text{Kek - Y}} = \alpha v_0 (1 + \beta \Delta) k
\]
- Kekulé - M:
\[
E_{\text{Kek - M}} = \beta \Delta m_0 v_0^2 + \alpha \sqrt{v_0^2 k^2 + m_0^2 v_0^4 (1 + \beta \Delta)^2}
\]

3. **Time evolution of quasiparticle coordinates**:
\[
\hat{r}(t) = \hat{r}_0 + \hat{\alpha} t + \hat{\xi}(t)
\]
where \( \hat{r}_0 \) is the initial position of the quasiparticle centroid, \( \hat{\alpha} \) is the drift velocity term, and \( \hat{\xi}(t) \) corresponds to Zitterbewegung.