Unexpected linear conductivity in Landau-Zener model: limitations and improvements of the relaxation time approximation in the quantum master equation

Ibuki Terada,Sota Kitamura,Hiroshi Watanabe,Hiroaki Ikeda
2024-01-30
Abstract:The nonequilibrium steady states of quantum materials have many challenges. Here, we highlight issues with the relaxation time approximation (RTA) for the DC conductivity in insulating systems. The RTA to the quantum master equation (QME) is frequently employed as a simple method, yet this phenomenological approach is exposed as a fatal approximation, displaying metallic DC conductivity in insulating systems within the linear response regime. We find that the unexpected metallic behavior is caused by the fact that the density matrix in the RTA incompletely incorporates the first order of the external field. To solve this problem, we have derived a new calculation scheme based on the QME that ensure correct behavior in low electric fields. Our method reproduces well the overall features of the exact electric currents in the whole field region. It is not time-consuming, and its application to lattice systems is straightforward. This method will encourage progress in this research area as a simple way to more accurately describe nonequilibrium steady states.
Mesoscale and Nanoscale Physics
What problem does this paper attempt to address?
This paper attempts to solve the inaccuracy problems that occur when using the Relaxation Time Approximation (RTA) method to calculate the DC conductivity in insulating systems. Specifically, when dealing with the Quantum Master Equation (QME), the RTA method exhibits metallic behavior for the DC conductivity in the linear response region, even though these systems are actually insulators. The reason for this inaccuracy lies in the fact that the density matrix in RTA fails to correctly include the first - order contribution of the external electric field. To solve this problem, the author proposes a new calculation scheme based on QME, which ensures that the physical behavior of the system can be correctly described under a low - electric - field condition. This scheme improves the description of the non - equilibrium steady - state by gradually introducing perturbative corrections of the external electric field, and its computational complexity is comparable to that of the traditional RTA method and can be conveniently applied to lattice systems. ### Key Problem Summary: 1. **Problems with the RTA method**: - In insulating systems, the RTA method causes the DC conductivity to exhibit metallic characteristics under a low - electric - field condition. - The reason for this abnormal behavior is that the density matrix fails to fully include the first - order contribution of the external electric field. 2. **Advantages of the new method**: - By more precisely handling the influence of the external electric field, the new method solves the flaws of the RTA method. - This method can correctly describe the current characteristics across the entire electric - field range. - It has high computational efficiency and is suitable for lattice systems. ### Mathematical Formula Representation: - **Density matrix elements under RTA**: \[ [\rho_k(t)]_{\alpha\alpha} \sim f_D(\epsilon_k^\alpha(t)) + E\tau_1 \frac{\partial f_D(\epsilon_k^\alpha(t))}{\partial k} \] \[ [\rho_k(t)]_{+-} \sim -\frac{W_k(t)}{\Delta_k(t) - i\tau_2^{-1}} \left(f_D(\epsilon_k^-(t)) - f_D(\epsilon_k^+(t))\right) \] - **The improved QME method**: \[ D[\rho_k(t)] = D_0 + D_1 + D_2 \] where, \[ [D_0]_{\alpha\beta} = -2\Gamma([\rho_k(t)]_{\alpha\beta} - f_D(\epsilon_k^\alpha(t))\delta_{\alpha\beta}) \] \[ [D_1]_{\alpha\beta} = -2\Gamma \frac{[W_k(t)]_{\alpha\beta}}{\Delta_k(t)} \delta f_k(t) \] \[ [D_2]_{\alpha\beta} = 2\Gamma \frac{|W_k(t)|^2}{\Delta_k^2(t)} \left(\delta f_k(t) + \Delta_k(t) f'_D(\epsilon_k^{-\alpha}(t))\right) \delta_{\alpha\beta} \] Through this method, the author successfully eliminates the unphysical linear terms in the RTA method and obtains the correct exponential behavior in the low - electric - field region.