Comment on "Spin-orbit coupling induced ultrahigh-harmonic generation from magnetic dynamics" with prescriptions on how to validate scientific software for computational quantum transport

Branislav K. Nikolic,Jalil Varela Manjarres
DOI: https://doi.org/10.48550/arXiv.2212.06895
2023-03-21
Abstract:In a recent paper [Phys. Rev. B {\bf 105}, L180415 (2022)], Ly and Manchon used open source code {\tt TKWANT} for time-dependent computational quantum transport to predict surprising features in the mature field of current pumping by magnetization dynamics in spintronics -- in the presence of spin-orbit (SO) coupling, the pumped charge current oscillates at both the frequency $\omega_0$ of magnetization precession and high harmonics $N=\omega/\omega_0$, reaching {\em astonishingly high} cutoff \mbox{$N_\mathrm{max} \simeq 1000$} by increasing the SO coupling. However, results in the paper violate two basic theorems of time-dependent quantum transport: ({\em i}) current response to time-periodic external field {\em must be perfectly periodic} itself in the long time limit for a two-terminal device because its active region is attached to two semi-infinite leads bringing continuous energy spectrum; and ({\em ii}) no DC component of charge current is allowed in the left-right symmetric two-terminal devices, or in asymmetric devices its value cannot be changed by simply increasing the SO coupling. We illustrate these two theorems by using completely different calculations applied to one-dimensional two-terminal devices with either ferromagnetic (for which the device is left-right symmetric) and antiferromagnetic (for the device is left-right asymmetric) active region hosting the Rashba SO coupling. We conclude that harmonics in pumped current in the presence of SO coupling do exist, but their ``ultrahigh'' cutoff is an artifact of either ``bugs'' or inadequate algorithms selected within {\tt TKWANT}. Finally, we suggest strategies for {\em validating} time-dependent quantum transport codes, or selection of algorithms by a user within putatively validated (by developers) code, prior to deploying them to produce research papers.
Mesoscale and Nanoscale Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is related to the inconsistencies and potential errors in the research carried out using the time - dependent quantum transport calculation software TKWANT. Specifically, the paper comments on and analyzes an article published by Ly and Manchon in 2022 (Phys. Rev. B 105, L180415). Ly and Manchon's research predicted that in the presence of spin - orbit coupling (SO coupling), the charge current pumped by magnetization dynamics would oscillate not only at the magnetization precession frequency \(\omega_0\) but also at high harmonics \(N\omega_0\), and the highest harmonic order \(N_{\text{max}}\) could reach approximately 1000. However, this paper points out that Ly and Manchon's results violate two fundamental theorems in time - dependent quantum transport: 1. **Periodicity Theorem**: For a two - terminal device, in the long - time limit, the current in response to a time - periodic external field must be completely periodic. This is because the active region of the device is connected to two semi - infinite leads with a continuous energy spectrum. 2. **DC Component Theorem**: For a left - right symmetric two - terminal device, a DC component is not allowed; for an asymmetric device, the value of the DC component cannot be changed only by increasing the SO coupling. To verify these theorems, the authors used a completely different calculation method, applied to a one - dimensional two - terminal device, where the active region contains ferromagnetic or antiferromagnetic materials and Rashba SO coupling is introduced. The results show that although high - order harmonics do exist in the presence of SO coupling, the "ultra - high" cutoff frequency may be due to a "bug" in the TKWANT code or the selection of an inappropriate algorithm. In addition, the authors propose several strategies to verify the validity of time - dependent quantum transport codes and methods for selecting appropriate algorithms to ensure the accuracy and reliability of research results. These include: - Using simple examples to verify whether the code satisfies the above two fundamental theorems. - Comparing the code results with known analytical solutions. - Using Floquet formalism for numerical calculations to obtain a more accurate number of high - order harmonics. In summary, this paper aims to reveal the potential problems in Ly and Manchon's research and provide suggestions for improving and verifying calculation methods to ensure the reliability and accuracy of future research.