Example: Franck-Condon Factors: NO2

Download FranckCondon_NO2.run

Note: this is a revised example.

As an example of a Franck-Condon calculation, lets look at the transition of NO2 to NO2 - . NO2 is a small molecule with only three vibrational modes. Putting an extra electron on the molecule will cause a big displacement, resulting in large electron-phonon couplings.

In general, the larger the molecule, the smaller the displacement and hence the electron-phonon couplings and Franck-Condon factors. Moreover, larger molecules have more vibrational modes, meaning that the already smaller displacement will generally be smeared out over more modes, resulting in an additional decrease in electron-phonon couplings. This is fortunate, since the number of Franck-Condon factors increases factorially with the number of vibrational modes, making it prohibitively expensive to take more than a few vibrational quanta into account for most molecules.

In order to calculate the Franck-Condon factors for Nitrite and Nitrogen dioxide, the equilibrium positions of the nuclei and the vibrational modes have to be obtained:

$ADFBIN/adf << eor
TITLE Nitrogen dioxide
ATOMS
    N         0.000000    0.000000   -0.016179
    O         0.000000    1.098646   -0.492918
    O         0.000000   -1.098646   -0.492918
END
BASIS
    CORE NONE
    TYPE DZP
END
XC
    LDA SCF VWN
END
ANALYTICALFREQ
END
UNRESTRICTED
CHARGE 0 1
GEOMETRY
END
eor

mv TAPE21 NO2.t21
rm logfile

$ADFBIN/adf << eor
TITLE Nitrite
ATOMS
    N         0.000000    0.000000    0.093662
    O         0.000000    1.120366   -0.540999
    O         0.000000   -1.120366   -0.540999
END
CHARGE -1.0
BASIS
    CORE NONE
    TYPE DZP
END
XC
    LDA SCF VWN
END
ANALYTICALFREQ
END
GEOMETRY
END
eor

mv TAPE21 NO2-.t21
rm logfile

This runscript produces two TAPE21 files containing the frequencies and the normal modes for both charge states. Lets first look at the ground state to ground state overlap:

$ADFBIN/fcf << eor
STATES NO2.t21 NO2-.t21
QUANTA 0 0
TRANSLATE
ROTATE
eor

Here, zero vibrational quanta are specified for both charge states, which corresponds to the vibrational ground state. Looking at the standard output, we see for NO2 :

Frequency (cm-1 ) \(\lambda\) (dimensionless)
756 1.979
1380 1.489
1716 0.000

And for NO2 - :

Frequency (cm-1 ) \(\lambda\) (dimensionless)
785 1.552
1265 0.000
1338 2.231

Both states have two vibrational modes with a significant electron-phonon coupling. The ground state to ground state Franck-Condon factor is therefore expected to be quite small. And indeed, looking at the output, we see that it is 0.349*10-1 , less than four percent of the total.

Since NO2 has only three vibrational modes, many quanta can be included, and this indeed turns out to be necessary. Setting the maximum number of quanta at 20 results in 1771 permutations for both states and a total of 3136441 Franck-Condon factors. Even with so many factors, the average sum is still only 0.463. Including one extra vibrational quanta results in an additional 960135 Franck-Condon factors, but an average sum of only 0.473, i.e. a percent more. This one percent is smeared out over so many factors that their individual contributions become negligible.