A Tale of Two Vectors: A Lanczos Algorithm For Calculating RPA Mean Excitation Energies
Research output: Contribution to journal › Journal article › Research › peer-review
Documents
- Postprint
Other version, 503 KB, PDF document
The experimental and theoretical determination of the mean excitation energy,
I(0), and the stopping power, S(v), of a material is of great interest in particle and
material physics, as well as radiation therapy. For calculations of I(0), the complete set of electronic transitions in a given basis set is required, effectively limiting such calculations to systems with a small number of electrons, even at the random-phase approximation (RPA)/time-dependent Hartree-Fock (TDHF) or time-dependent density functional theory (TDDFT) level. To overcome such limitations, we present here the implementation of a Lanczos algorithm adapted for the paired RPA/TDHF eigenvalue problem in the Dalton program and show that it provides good approximations of the entire RPA eigenspectra in a reduced space. We observe rapid convergence of I(0) with the number of Lanczos vectors as the algorithm favors the transitions with large contributions. In most cases, the algorithm recovers RPA I(0) values of up to 0.5 % accuracy at less than a quarter of the full space size. The algorithm not only exploits the RPA paired structure to save computational resources, but it is also preserves certain sum-over-states properties, as first demonstrated by Johnson et al. [Comput. Phys. Commun. 1999, 120, 155]. The block Lanczos RPA solver, as presented here, thus shows promise for computing mean excitation energies for systems larger than what was computationally feasible before.
I(0), and the stopping power, S(v), of a material is of great interest in particle and
material physics, as well as radiation therapy. For calculations of I(0), the complete set of electronic transitions in a given basis set is required, effectively limiting such calculations to systems with a small number of electrons, even at the random-phase approximation (RPA)/time-dependent Hartree-Fock (TDHF) or time-dependent density functional theory (TDDFT) level. To overcome such limitations, we present here the implementation of a Lanczos algorithm adapted for the paired RPA/TDHF eigenvalue problem in the Dalton program and show that it provides good approximations of the entire RPA eigenspectra in a reduced space. We observe rapid convergence of I(0) with the number of Lanczos vectors as the algorithm favors the transitions with large contributions. In most cases, the algorithm recovers RPA I(0) values of up to 0.5 % accuracy at less than a quarter of the full space size. The algorithm not only exploits the RPA paired structure to save computational resources, but it is also preserves certain sum-over-states properties, as first demonstrated by Johnson et al. [Comput. Phys. Commun. 1999, 120, 155]. The block Lanczos RPA solver, as presented here, thus shows promise for computing mean excitation energies for systems larger than what was computationally feasible before.
Original language | English |
---|---|
Article number | 014102 |
Journal | The Journal of Chemical Physics |
Volume | 156 |
Issue number | 1 |
Number of pages | 16 |
ISSN | 0021-9606 |
DOIs | |
Publication status | Published - 3 Jan 2022 |
- Faculty of Science - Mean excitation energy, stopping power, RPA, Lanczos, random phase approximation, time-dependent Hartree Fock
Research areas
Number of downloads are based on statistics from Google Scholar and www.ku.dk
No data available
ID: 286431719