CrossingBarriers
This page gives hints on how to calculate barriers for crossings with the ABINIT package, as well as ensemble DFT, or the pSIC approach to polaron formation.
Introduction¶
The knowledge of geometries at which crossings between two electronic states happen, with minimal energy, or geometries at which the energy difference between the ground state and the excited state is small, and the energy is still low, plays an important role in the study of non-radiative transitions.
It is possible to formulate the search for such geometries in terms of minimisation of a functional that is the linear combination of the energy of the two states at the same geometry, with Lagrange multipliers [Jia2019]. This is also related with a simple approach to Ensemble DFT: just make a linear combination of the DFT energies, the XC correlation energy being not computed with a single common density, but from each density separately. Also, the pSIC, polaron self-interaction corrected method [Sadigh2015], [Sadigh2015a], can be formulated in the same terms.
In ABINIT, with imgmov==6, it is possible to deal with such linear combination of systems with the same geometry, but differing occupation factors occ, and even with different cellcharge.. It is possible to find the geometry at which the resulting energy is minimal, for a given value of the mixing factors mixesimgf. Set nimage=2, and set the occupation numbers for image 1 to the ground-state occupations, and for image 2 to the excited-state occupations.
Related Input Variables¶
useful:
- imgmov IMaGe MOVEs
- imgwfstor IMaGe WaveFunction STORage
- mixesimgf MIXing Electronic Structure IMAGE Factors
Selected Input Files¶
tutoparal:
- tests/tutoparal/Input/tpsic_01.abi
- tests/tutoparal/Input/tpsic_02.abi
- tests/tutoparal/Input/tpsic_03.abi
v9: