Dear users,
I have a question about el-ph calculation in magnetic systems (nspin). In
PH/phq_readin.f90, there is
the line
IF (elph.AND.lsda) CALL errore ('phq_readin', 'El-ph and spin not &
implemented', 1)
that stops a el-ph calculation when nspin=2.. But having a look at the
file PH/elphon.
Andrea Floris wrote:
> IF (elph.AND.lsda) CALL errore ('phq_readin', 'El-ph and spin not &
> implemented', 1)
I think spin-polarized electron-phonon calculations are implemented,
but untested. Since nobody has ever reported if it works, the check
is still there.
P.
--
Paolo Giannozzi, Democrito