I’ve been using GULP to calculate phonon modes for instantaneous liquid configurations which exhibit imaginary frequencies (I think solids near melting point also do this). It does this perfectly; however, the .dens file created by the “output phonon” command does not include the imaginary frequencies in its calculation. Is there someway to include these frequencies in the PDOS as negative numbers? GULP already uses that convention for the frequency output.
I attached the input and outputs incase that is helpful.
LJ_liq.gin (12.2 KB)
LJ_liq_phon.dens (1.7 KB)
LJ_liq (9 KB)