help for phonon spectral function of single phonon mode

[DFT13]

Dear kbspooner, I wonder how to plot the frequency-depencent phonon spectral function of single phonon mode, like this, can you give me some suggestions?

[kbspooner]

I suppose that depends on the format of your data. If you can put it into a plain text file of the same format as the phonopy projected_dos.dat it could be plotted fairly directly with ThermoParser, or depending on your confidence with python (or desire to learn) you could write a small function to load it into the ThermoParser format and plot it in the same way. In the latter case I would be happy to incorporate it into this code (and credit you, of course). Hope this helps!