1. Efficient update of determinants of rmany-electron wave function overlaps
- Author
-
Universitat Politècnica de València. Departamento de Sistemas Informáticos y Computación - Departament de Sistemes Informàtics i Computació, Universitat Politècnica de València. Departamento de Informática de Sistemas y Computadores - Departament d'Informàtica de Sistemes i Computadors, Generalitat de Catalunya, AGENCIA ESTATAL DE INVESTIGACION, Agencia Estatal de Investigación, Croatian Science Foundation (CSF), European Regional Development Fund, Ministerio de Economía y Competitividad, Alonso-Jordá, Pedro, Davidovic, Davor, Sapunar, Marin, Herrero, Jose R., Quintana-Ortí, Enrique S., Universitat Politècnica de València. Departamento de Sistemas Informáticos y Computación - Departament de Sistemes Informàtics i Computació, Universitat Politècnica de València. Departamento de Informática de Sistemas y Computadores - Departament d'Informàtica de Sistemes i Computadors, Generalitat de Catalunya, AGENCIA ESTATAL DE INVESTIGACION, Agencia Estatal de Investigación, Croatian Science Foundation (CSF), European Regional Development Fund, Ministerio de Economía y Competitividad, Alonso-Jordá, Pedro, Davidovic, Davor, Sapunar, Marin, Herrero, Jose R., and Quintana-Ortí, Enrique S.
- Abstract
[EN] The calculation of overlaps between many-electron wave functions at different nuclear geometries during nonadiabatic dynamics simulations requires the evaluation of a large number of determinants of matrices that differ only in a few rows/columns. While this calculation is fast for small systems, its cost grows faster than the alternative electronic structure calculation used to obtain the wave functions. For wave functions that can be written as a CIS expansion, all determinants can be computed using the set of level-2 minors of the reference matrix. However, this is still a costly computation for large systems. In this paper, we provide an algorithm for efficiently calculating all level-2 minors of a matrix by re-utilizing and updating the LU factorization for the determinants of the minors. This approach results in a parallel version of the algorithm that is up to an order of magnitude faster then the current best parallel implementation. The algorithm thus allows the computation of exact wave function overlaps for relatively large systems, with a high density of states, at virtually no cost compared with the electronic structure calculations. Furthermore, the new algorithm opens the path to further investigations in efficient computing of the exact wave function overlaps for complex wave functions such as MR-CIS and MR-CISD. Program summary Program Title: CIS Overlap Licensing provisions: MIT Programming language: FORTRAN 2008, C Nature of problem: Calculation of overlaps between CIS type wave functions at different nuclear geometries during nonadiabatic dynamics simulation requires calculating a large number of connected determinants and scales with the seventh power of the size of the system being studied. Without additional approximations, for large systems this computation becomes more costly than the electronic structure calculation used to obtain the wave functions. Solution method: All of the determinants required for a CIS wave function overlap calc
- Published
- 2021