This is an old revision of the document!
Calcualtion of e-ph coupling Here we describe how to calculate the vibrational modes of molecules using dynamical matrix. Note, the calculation of the e-ph coupling is meant only for non periodic systems or with only Gamma k-point.
iephc = 1
The calculation e-ph coupling is toggled by the keywords iephc
in section &OPTION&
of the input file fireball in
&OPTION basisfile = answer.bas icluster = 1 sigmatol = 0.00000001 max_scf_iterations = 100 iephc = 1 &END
Control parameters of the e-ph coupling calculation are set in file ephc.optional
, which must be present in the work directory.
ephc.dat ! output file eigens.dat ! input file with reference eigenvalues (for zero displacement) 20 ! effective temperature [K] 28.08 ! effective mass [a.u.] 4 ! number of eigenstates which will be analyzed 40,41,42,43 ! list of eigenstates (those listed in eigens.dat) 4 ! number of modes to be analyzed (0...means all) 35,36,37,38 ! list of vibrational modes
We have to define which molecular orbitals will be analyzed. Eigenvalues of the selected molecular orbitals must be listed in the file eigens.dat