-
Cumulant methods for electron-phonon problems. II. The self-consistent cumulant expansion
Authors:
Paul J. Robinson,
Ian S. Dunn,
David R. Reichman
Abstract:
In this work we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattices sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demons…
▽ More
In this work we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattices sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demonstration that the SC-CE corrects many deficits of the standard perturbative CE in the (non-dispersive) Holstein model. The natural phonon damping that arises in the more realistic dispersed Holstein model renders the SC-CE stable, allowing for a complete assessment of the method. Here we find that self-consistency dramatically corrects many of the failures found in the perturbative CE, but also introduces some unphysical features. Finally, we comment on the potential use of SC-CE as a tool for calculating Green's functions in generic many-body problems.
△ Less
Submitted 8 June, 2022; v1 submitted 4 March, 2022;
originally announced March 2022.
-
Cumulant methods for electron-phonon problems. I. Perturbative expansions
Authors:
Paul J. Robinson,
Ian S. Dunn,
David R. Reichman
Abstract:
In this work we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CE for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as w…
▽ More
In this work we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CE for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for non-zero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant, and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.
△ Less
Submitted 8 June, 2022; v1 submitted 4 March, 2022;
originally announced March 2022.
-
Removing instabilities in the hierarchical equations of motion: exact and approximate projection approaches
Authors:
Ian S. Dunn,
Roel Tempelaar,
David R. Reichman
Abstract:
The hierarchical equations of motion (HEOM) provide a numerically exact approach for computing the reduced dynamics of a quantum system linearly coupled to a bath. We have found that HEOM contains temperature-dependent instabilities that grow exponentially in time. In the case of continuous-bath models, these instabilities may be delayed to later times by increasing the hierarchy dimension; howeve…
▽ More
The hierarchical equations of motion (HEOM) provide a numerically exact approach for computing the reduced dynamics of a quantum system linearly coupled to a bath. We have found that HEOM contains temperature-dependent instabilities that grow exponentially in time. In the case of continuous-bath models, these instabilities may be delayed to later times by increasing the hierarchy dimension; however, for systems coupled to discrete, non-dispersive modes, increasing the hierarchy dimension does little to alleviate the problem. We show that these instabilities can also be removed completely at a potentially much lower cost via projection onto the space of stable eigenmodes; furthermore, we find that for discrete-bath models at zero temperature, the remaining projected dynamics computed with few hierarchy levels are essentially identical to the exact dynamics that otherwise might require an intractably large number of hierarchy levels for convergence. Recognizing that computation of the eigenmodes might be prohibitive, e.g. for large or strongly-coupled models, we present a Prony filtration algorithm that may be useful as an alternative for accomplishing this projection when diagonalization is too costly. We present results demonstrating the efficacy of HEOM projected via diagonalization and Prony filtration. We also discuss issues associated with the nonnormality of HEOM.
△ Less
Submitted 10 May, 2019; v1 submitted 12 February, 2019;
originally announced February 2019.