Abstract:With the recent advancements in laser technology, there has been increasing interest in nonlinear and nonperturbative phenomena such as nonreciprocal transport, the nonlinear Hall effect, and nonlinear optical responses. When analyzing the nonequilibrium steady state, the relaxation time approximation (RTA) in the quantum kinetic equation has been widely used. However, recent studies have highlighted problems with the use of RTA that require careful consideration. In a study published in Phys. Rev. B, $\textbf{109}$, L180302 (2024), we revealed that the RTA has a flaw in predicting finite linear conductivity even for insulators under weak electric fields, and improved the RTA based on the Redfield equation. In this paper, we further extend our approach to nonlinear responses. This approach provides a simple alternative to RTA and is expected to be useful for the study of nonlinear and nonequilibrium phenomena.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to solve the inaccuracies and problems that occur when using the Relaxation Time Approximation (RTA) method to calculate the nonlinear conductivity in non - centrosymmetric insulators. Specifically:
1. **Limitations of RTA**:
- The paper points out that under a weak electric field, RTA has flaws in predicting the finite linear conductivity of insulators. Even under weak electric field conditions, the linear conductivity predicted by RTA is not zero for insulators, which is inconsistent with the actual situation.
- These problems are not limited to the linear response but also extend to the nonlinear response domain.
2. **Improving the RTA method**:
- The author proposes the Dynamical Phase Approximation (DPA) method based on the Redfield equation to improve RTA. The DPA method can describe the nonlinear response more accurately and eliminate the unphysical contributions in RTA.
- Especially in the nonlinear response, DPA can correct the abnormal second - order nonlinear conductivity term caused by RTA, ensuring that the behavior of insulators under a weak DC electric field is as expected.
3. **Research background and motivation**:
- With the development of laser technology, people's interest in nonlinear and non - perturbative phenomena (such as non - reciprocal transport, the nonlinear Hall effect, and the nonlinear optical response) is increasing. These phenomena cannot be explained by traditional perturbation theory.
- Therefore, identifying and correcting the problems in RTA is crucial for accurately understanding nonlinear and non - equilibrium phenomena.
### Main conclusions
- By introducing the DPA method, the paper solves the key defects of RTA in the nonlinear response, especially in non - centrosymmetric insulators.
- The results show that the DPA method can restore the insulating properties of insulators under a weak electric field, eliminate the unphysical contributions in RTA, and provide a more accurate alternative for studying nonlinear and non - equilibrium phenomena.
### Formula summary
- **Quantum kinetic equation**:
\[
\frac{d\rho_k(t)}{dt}=-i[H_k(t),\rho_k(t)] + D_k(\rho_k(t))
\]
where \(\rho_k(t)\) is the density matrix, \(H_k(t)\) is the Hamiltonian with an external electric field, and \(D_k(\rho_k(t))\) is the dissipation term.
- **Dissipation term in RTA**:
\[
D(\rho_k(t)) = -\frac{\rho_k(t)-\rho_{eq,k}}{\tau}
\]
where \(\rho_{eq,k}\) is the density matrix in the equilibrium state and \(\tau\) is the relaxation time.
- **Dissipation term in DPA**:
\[
D_k = D_k^{(0)}+D_k^{(1)}+D_k^{(2)}
\]
The specific forms are as follows:
\[
[D_k^{(0)}]_{\alpha\beta}=-\frac{[\rho_k(t)]_{\alpha\beta}-f_D(\varepsilon_{k\alpha}(t))\delta_{\alpha\beta}}{\tau}
\]
\[
[D_k^{(1)}]_{\alpha\beta}=-\frac{[W_k(t)]_{\alpha\beta}\delta f_k(t)}{\tau\Delta_k(t)(1 - E/\Delta_k(t)z_{k\alpha}(t))}
\]
\[
[D_k^{(2)}]_{\alpha\beta}=\alpha\left|[W_k(t)]_{+,-}\right|^2\delta f_k(t)\frac{(1 - E/\Delta_k(t)z_{k+}(t))(1 - E/