observed in multiple glasses [@Schober_1996; @Mazzacurati_1996; @Taraskin_1997; @Taraskin_1999]. Importantly, $P^k$ below $\omega_\mathrm{BP}$ collapses between different values of $\varepsilon_\mathrm{bend}$. This result indicates that the variations of not only the vDOS and the vibrational states due to $\varepsilon_\mathrm{bend}$ can be characterized by the macroscopic shear modulus changes. However, $P^k$ does not collapse above $\omega_\mathrm{BP}$, as also shown in Fig. \[fig7\]. This result is attributed to the fact that the high-frequency modes above $\omega_\mathrm{BP}$ reflect microscopic vibrations that cannot be captured by the macroscopic elasticity.

### Comparison with LJ glasses

The low-frequency spectra are comparable to that of atomic LJ glasses reported in Ref. [@Shimada:2018fp]. As observed in Fig. \[fig4\], the height of $g(\omega_\mathrm{BP})/\omega_\mathrm{BP}^2$ of LJ glasses is higher than that of polymer glasses, and $\omega_\mathrm{BP}$ is lower than that of polymer glasses. These observations are different from the study reported in Ref. [@Giuntoli:2018fv], which demonstrated that the low-frequency spectra of polymer glasses correspond to those atomic LJ glasses. In Ref. [@Giuntoli:2018fv], the bonded monomers interact via a harmonic potential with a large bonding energy scale of $k = 2500$. This value is two orders of magnitude larger than $\varepsilon_\mathrm{FENE} = 30$, investigated in this study. With respect to the large bonding energy, 