Arne Schobert, Jan Berges, Erik G. C. P. van Loon, Michael A. Sentef, Sergey Brener, Mariana Rossi, Tim O. Wehling
SciPost Phys. 16, 046 (2024) ·
published 9 February 2024
|
· pdf
The interplay of electronic and nuclear degrees of freedom presents an outstanding problem in condensed matter physics and chemistry. Computational challenges arise especially for large systems, long time scales, in nonequilibrium, or in systems with strong correlations. In this work, we show how downfolding approaches facilitate complexity reduction on the electronic side and thereby boost the simulation of electronic properties and nuclear motion-in particular molecular dynamics (MD) simulations. Three different downfolding strategies based on constraining, unscreening, and combinations thereof are benchmarked against full density functional calculations for selected charge density wave (CDW) systems, namely 1H-TaS$_2$, 1T-TiSe$_2$, 1H-NbS$_2$, and a one-dimensional carbon chain. We find that the downfolded models can reproduce potential energy surfaces on supercells accurately and facilitate computational speedup in MD simulations by about five orders of magnitude in comparison to purely ab initio calculations. For monolayer 1H-TaS$_2$ we report classical and path integral replica exchange MD simulations, revealing the impact of thermal and quantum fluctuations on the CDW transition.
Arne Schobert, Jan Berges, Tim Wehling, Erik van Loon
SciPost Phys. 11, 079 (2021) ·
published 20 October 2021
|
· pdf
Charge-density waves are responsible for symmetry-breaking displacements of atoms and concomitant changes in the electronic structure. Linear response theories, in particular density-functional perturbation theory, provide a way to study the effect of displacements on both the total energy and the electronic structure based on a single ab initio calculation. In downfolding approaches, the electronic system is reduced to a smaller number of bands, allowing for the incorporation of additional correlation and environmental effects on these bands. However, the physical contents of this downfolded model and its potential limitations are not always obvious. Here, we study the potential-energy landscape and electronic structure of the Su-Schrieffer-Heeger (SSH) model, where all relevant quantities can be evaluated analytically. We compare the exact results at arbitrary displacement with diagrammatic perturbation theory both in the full model and in a downfolded effective single-band model, which gives an instructive insight into the properties of downfolding. An exact reconstruction of the potential-energy landscape is possible in a downfolded model, which requires a dynamical electron-biphonon interaction. The dispersion of the bands upon atomic displacement is also found correctly, where the downfolded model by construction only captures spectral weight in the target space. In the SSH model, the electron-phonon coupling mechanism involves exclusively hybridization between the low- and high-energy bands and this limits the computational efficiency gain of downfolded models.
