skelton-group / Phonopy-Spectroscopy

A collection of tools for simulating vibrational spectra, which interfaces with the Phonopy package.
MIT License
148 stars 56 forks source link

a warning is displayed:invalid value encountered in divide #14

Open studyvasp opened 1 year ago

studyvasp commented 1 year ago

When I calculate the Raman spectrum of glycine, a warning is displayed:invalid value encountered in divide intensities = np.divide(intensities, norm) The I(v) [Ang^4 amu^-1] in Raman-Spectrum.dat is 0.00000,but the data in Raman is normalRaman-PeakTable.dat is right, Could you give me some help?

JMSkelton commented 1 year ago

When you say the data in Raman-PeakTable.dat is "right", do you mean the band intensities are non-zero?

Without any additional information, this error means that, for whatever reason, the calculated band intensities passed to the plotting code are all zero.

In order to determine whether this is a bug in the code or an issue with the calculation, can you please post your Raman.yaml file and the command you used to generate the plot?

studyvasp commented 1 year ago

Command according to the example:phonopy-raman -p --ir-reps --linewidth-hdf5="kappa-m484848-g0.hdf5" --linewidth-temperature=300 Raman.yaml is: distance_units: Ang frequency_units: THz step_units: sqrt(amu) * Ang

lattice_vectors:

displacement_sets: