Subsequent to the equilibration, the Type I production is performed to compute the thermodynamic properties. Specifically, a 10 ns production run is conducted under the NPT ensemble, with the thermodynamic output, such as mass density, volume of simulation box, and potential energy, recorded every 1 ps. As one of the most common thermodynamic properties of PDMS materials, the mass density is calculated from simulations following the equation: