Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userReal-Time Diagrammatic Monte Carlo for Nonequilibrium Quantum Transport
187
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We propose a novel approach to nonequilibrium real-time dynamics of quantum impurities models coupled to biased non-interacting leads, such as those relevant to quantum transport in nanoscale molecular devices. The method is based on a Diagrammatic Monte Carlo sampling of the real-time perturbation theory along the Keldysh contour. We benchmark the method on a non-interacting resonant level model and, as a first non-trivial application, we study zero temperature non-equilibrium transport through a vibrating molecule.
rate research
Read More
Diagrammatic expansions are a central tool for treating correlated electron systems. At thermal equilibrium, they are most naturally defined within the Matsubara formalism. However, extracting any dynamic response function from a Matsubara calculation ultimately requires the ill-defined analytical continuation from the imaginary- to the real-frequency domain. It was recently proposed [Phys. Rev. B 99, 035120 (2019)] that the internal Matsubara summations of any interaction-expansion diagram can be performed analytically by using symbolic algebra algorithms. The result of the summations is then an analytical function of the complex frequency rather than Matsubara frequency. Here we apply this principle and develop a diagrammatic Monte Carlo technique which yields results directly on the real-frequency axis. We present results for the self-energy $Sigma(omega)$ of the doped 32x32 cyclic square-lattice Hubbard model in a non-trivial parameter regime, where signatures of the pseudogap appear close to the antinode. We discuss the behavior of the perturbation series on the real-frequency axis and in particular show that one must be very careful when using the maximum entropy method on truncated perturbation series. Our approach holds great promise for future application in cases when analytical continuation is difficult and moderate-order perturbation theory may be sufficient to converge the result.
We present the first approximation free diagrammatic Monte Carlo study of a lattice polaron interacting with an acoustic phonon branch through the deformation potential. Weak and strong coupling regimes are separated by a self-trapping region where quantum resonance between various possible lattice deformations is seen in the ground state properties, spectral function, and optical conductivity. The unique feature of such polaron is the interplay between long- and short wavelength acoustic vibrations creating a composite phonon cloud and leading to persistent self-trapping due to the existence of multiple quasi-stable states. This results in a spectral response whose structure is much more complex than in any of the previously considered polaron models.
The Quantum Monte Carlo (QMC) method can yield the imaginary-time dependence of a correlation function $C(tau)$ of an operator $hat O$. The analytic continuation to real-time proceeds by means of a numerical inversion of these data to find the response function or spectral density $A(omega)$ corresponding to $hat O$. Such a technique is very sensitive to the statistical errors in $C(tau)$ especially for large values of $tau$, when we are interested in the low-energy excitations. In this paper, we find that if we use the flat histogram technique in the QMC method, in such a way to make the {it histogram of} $C(tau)$ flat, the results of the analytic continuation for low-energy excitations improve using the same amount of computational time. To demonstrate the idea we select an exactly soluble version of the single-hole motion in the $t-J$ model and the diagrammatic Monte Carlo technique.
We present a simple trick that allows to consider the sum of all connected Feynman diagrams at fixed position of interaction vertices for general fermionic models. With our approach one achieves superior performance compared to Diagrammatic Monte Carlo, while rendering the algorithmic part dramatically simpler. As we consider the sum of all connected diagrams at once, we allow for cancellations between diagrams with different signs, alleviating the sign problem. Moreover, the complexity of the calculation grows exponentially with the order of the expansion, which should be constrasted with the factorial growth of the standard diagrammatic technique. We illustrate the efficiency of the technique for the two-dimensional Fermi-Hubbard model.
High order perturbation theory has seen an unexpected recent revival for controlled calculations of quantum many-body systems, even at strong coupling. We adapt integration methods using low-discrepancy sequences to this problem. They greatly outperform state-of-the-art diagrammatic Monte Carlo. In practical applications, we show speed-ups of several orders of magnitude with scaling as fast as $1/N$ in sample number $N$; parametrically faster than $1/sqrt{N}$ in Monte Carlo. We illustrate our technique with a solution of the Kondo ridge in quantum dots, where it allows large parameter sweeps.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
