Prediction of Structure and Energy of Trans-1,4-Polybutadiene Glassy Surface by Atomistic Simulations of Free-Standing Ultrathin Films

Atomistic modeling of amorphous trans-1,4-polybutadiene (TPBD), using molecular mechanics and molecular dynamics (MD) simulations, is performed to generate three-dimensionally periodic bulk and two-dimensionally periodic thin film condensed phases. The condensed structures are constructed using multiple polymer chains. Structural and energetic relaxations and sampling of properties are performed using MD in the canonical ensemble (NVT) by a procedure that relieves local high-energy spots and brings the system to realistic thermodynamic states. The calculated surface energy for TPBD, 30.72 erg/cm², is in excellent agreement with the reported experimental value of 31 erg/cm². The structure of the surface layers is probed in terms of the atomic mass density variations, bond-bond orientation function profiles, and the distribution of the dihedral angles about the rotatable backbone bonds. The thickness of the surface layer over which the density varies smoothly but rapidly is found to be approximately 15 Å. The level of agreement of the calculated surface energy with the experimental value is superior in comparison to previous investigations in the literature using the atomistic approach for flexible polymers.