Formaldehyde and hydrogen peroxide are two important realistic molecules in atmospheric chemistry. We implement path integral Liouville dynamics (PILD) to calculate the dipole-derivative autocorrelation function for obtaining the infrared (IR) spectrum. In comparison to exact vibrational frequencies, PILD faithfully capture most nuclear quantum effects in vibrational dynamics as temperature changes and as the isotopic substitution occurs.