The mechanical properties of the systems computed by MD simulations with the COMPASS are represented in Table 2, respectively. The shear modulus for each model is calculated using Reuss (GR) and Voigt (GV) methods and then averaged with that obtained by Hill method (GH) to achieve shear modulus averaging GVRH [15].