role for studying the vibrational and mechanical properties of glasses. Firstly, MD simulations enable to assess the theoretical predictions. In fact, various MD simulations have been performed on simple atomic glasses, e.g., Lennard-Jones (LJ) systems [@Schober_1996; @Mazzacurati_1996; @Shintani_2008; @Monaco2_2009; @Mizuno2_2013; @Lerner_2016; @Wang_2019]. Concerning the isostaticity and marginal stability [@Wyart_2005; @Wyart_2006; @Wyart_2010; @DeGiuli_2014], the systems with a finite-ranged, purely repulsive potential have also been studied [@Silbert_2005; @Silbert_2009; @Vitelli_2010; @Xu_2010], and are considered as the simplest model of glasses. In particular, it is crucial for MD simulations to solve finite-dimensional effects that are not captured by the mean-field treatments [@Mizuno_2017; @Shimada_2018; @Mizuno_2018]. Secondly, MD simulations perform quasi-experiments on well-defined systems and access data that cannot be examined experimentally. Relevant systems to experiments and applications have been simulated, including covalent-bonding [@Taraskin_1997; @Taraskin_1999; @Horbach_2001; @Leonforte_2006; @Beltukov_2016; @Beltukov_2018], metallic [@Derlet_2012; @Fan_2014; @Crespo_2016; @Brink_2016], polymeric [@Jain_2004; @Schnell:2011he; @Ness:2017cc; @Milkus:2018ha; @Giuntoli:2018fv] glasses. These simulation studies complete theoretical understandings based on simple systems and experimental observations of more complex systems.

The vibrational properties and the BP of polymeric glasses have been studied by both of experiments [@Niss_2007; @Hong_2008; @Caponi_2011; @Perez-Castaneda2_2014; @Terao_2018; @Zorn_2018] and MD simulations [@Jain_2004; @Schnell:2011he; @Ness:2017cc; @Milkus:2018ha; @Giuntoli:2018fv]. The effects caused by non-covalent bonds including bending 