Abstract:High Harmonic Generation (HHG), which has been widely used in atomic gas, has recently expanded to solids as a means to study highly nonlinear electronic response in condensed matter and produce coherent high frequency radiation with new properties. Most recently, attention has turned to Topological Materials (TMs) and the use of HHG to characterize topological bands and invariants. Theoretical interpretation of nonlinear electronic response in TMs, however, presents many challenges. In particular, the Bloch wavefunction phase of TMs has undefined points in the Brillouin Zone. This leads to singularities in calculating the inter-band and intra-band transition dipole matrix elements of Semiconductor Bloch Equations (SBEs). Here, we use the laser-electromagnetic velocity gauge ${\boldsymbol p}\cdot {\bf A}(t)$ to numerically integrate the SBEs and treat the singularity in the production of the electrical currents and HHG spectra. We use a prototype of Chern Insulators (CIs), the Haldane model, to demonstrate our approach. We find good qualitative agreement of the velocity gauge compared to the length gauge and the Time-Dependent Density Functional theory in the case of topologically trivial materials such as MoS$_2$. For velocity gauge and length gauge, our two-band Haldane model reproduces key HHG spectra features: ($\textit i$) The selection rules for linear and circular light drivers, ($\textit ii$) The linear cut-off law scaling and ($\textit iii$) The anomalous circular dichroism. We conclude that the velocity-gauge approach captures experimental observations and provides theoretical tools to investigate topological materials.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in topological materials (TMs), the computational difficulties in the semiconductor Bloch equations (SBEs) due to the singularity of the wave - function phase, especially when dealing with high - harmonic generation (HHG). Specifically, the article focuses on how to use the laser - electromagnetic velocity gauge (VG) to numerically integrate SBEs and handle these singularities to accurately simulate the charge current and HHG spectra. ### Problem Background 1. **Singularities in Topological Materials** Topological materials have special electronic structures, and their wave - function phases have undefined points in the Brillouin zone (BZ), which will lead to singular points. These singular points make the calculation of inter - band and intra - band transition dipole moment matrix elements complicated, thus affecting the accuracy of the HHG spectra. 2. **Nonlinear Optical Response** High - harmonic generation (HHG), as a nonlinear optical phenomenon, has been widely used in atomic gases. In recent years, researchers have begun to extend it to solid materials, especially for studying the highly nonlinear electronic response in condensed matter and generating coherent high - frequency radiation with new properties. However, for topological materials, theoretical explanations of the nonlinear electronic response still face many challenges. ### Core Problems of the Paper To meet the above challenges, this paper proposes a new method, that is, using the laser - electromagnetic velocity gauge $ \mathbf{p} \cdot \mathbf{A}(t) $ to numerically integrate SBEs, thereby handling singularities and accurately simulating the charge current and HHG spectra. Specifically: - **The Role of the Velocity Gauge** Using the velocity gauge can avoid the singularity problems that occur in the length gauge, because the velocity gauge decouples adjacent $ k $-points, so that numerical calculations can be carried out on a discrete $ k $-grid without encountering singular points. - **Model Verification** The article uses the Haldane model as a prototype of Chern insulators (CIs) to verify the effectiveness of the velocity gauge method. By comparing the results of the velocity gauge with those of the length gauge and the time - dependent density functional theory (TDDFT), it proves the advantages of the velocity gauge method in dealing with the nonlinear optical response in topological materials. ### Main Conclusions 1. **Qualitative Consistency** In topologically trivial materials such as MoS₂, the HHG spectra calculated by the velocity gauge show good qualitative consistency with those calculated by the length gauge and TDDFT, especially in terms of the plateau structure and the cutoff frequency. 2. **Selection Rules and Circular Dichroism** For topological materials, the velocity gauge can reproduce key HHG features, including selection rules, linear cutoff laws, and anomalous circular dichroism. For example, the calculation results of the selection rules driven by linear and circularly polarized light and circular dichroism (CD) are consistent with those of the length gauge. 3. **Introduction of Theoretical Tools** The velocity gauge method provides a new theoretical tool for studying highly nonlinear optical emissions in topological materials and solves the singularity problems existing in the length gauge. In summary, this paper aims to solve the computational problems in topological materials caused by the singularity of the wave - function phase, proposes the method of using the laser - electromagnetic velocity gauge, and verifies its effectiveness through numerical simulations.

Theory for all-optical responses in topological materials: the velocity gauge picture

Measuring Electromagnetic and Gravitational Responses of Photonic Landau Levels

"Shaking" Photons out of a Topological Material

Probing Chern number by opacity and topological phase transition by a nonlocal Chern marker

Giant room-temperature nonlinearities from a monolayer Janus topological semiconductor

High-order harmonics from elliptically polarized laser for classification of topological states in 2D Chern insulators and 2D topological insulators

Theory of topological exciton insulators and condensates in flat Chern bands

A modal decomposition approach to topological wave propagation

Van Roosbroeck's equations with topological terms: the case of Weyl semimetals

Structure- and laser-gauges for the semiconductor Bloch equations in high-harmonic generation in solids

Giant room-temperature nonlinearities in a monolayer Janus topological semiconductor

Linear and nonlinear optical response based on many-body GW-Bethe-Salpeter and Kadanoff-Baym approaches for two-dimensional layered semiconductors

Are there universal signatures of topological phases in high harmonic generation? Probably not

Role of topological charges in the nonlinear-optical response from Weyl semimetals

Optical absorption measurement of spin Berry curvature and spin Chern marker

Linear and nonlinear optical response based on the GW-Bethe-Salpeter and Kadanoff-Baym approaches for two-dimensional layered semiconductor materials

Ultrafast and Strong-Field Physics in Graphene-Like Crystals: Bloch Band Topology and High-Harmonic Generation

First Principles Calculation of the Topological Phases of the Photonic Haldane Model

Artificial gauge fields in the t-z mapping for optical pulses: spatio-temporal wavepacket control and quantum Hall physics

Topological high-harmonic spectroscopy

Second Harmonic Generation in Topological Insulators under Quantizing Magnetic Fields