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 spectrum. In comparison to exact vibrational frequencies, PILD faithfully captures most nuclear quantum effects in vibrational dynamics as temperature changes and as the isotopic substitution occurs.