Molecular dynamics simulation of polystyrene copolymer with octyl short-chain branches in toluene

In this study, dimensional, conformational and dynamic behaviors of a short-chain branched styrene/1-octene copolymer chain with different 1-octene percentages, i.e., 0, 2, 4 and 6%, in toluene are investigated at the temperature of 298.15 K via molecular dynamics simulation. The chain dimensions and flexibility in the solvent are evaluated by calculating the radius of gyration (R_g), end-to-end distance (<r₀>), surface area (A_ch), and volume (V_ch) of the copolymer chain. The mean square displacement (MSD) and diffusivity coefficient for each copolymer chain are measured to determine its dynamic behavior and mobility in aromatic media. To consider the effect of increasing the 1-octene co-monomer percentage on the copolymer chain affinity to the solvent molecules, the interaction energy (E_int) and Flory-Huggins (FH) interaction parameter are calculated for each equilibrated solution model. The simulation results indicate that the co-monomer level increment in the copolymer structure reduces the chain R_g amount and its interaction with the solvent. The <r₀> of the chain increases up to 4% co-monomer content, while further co-monomer content decreases the <r₀> value. Additionally, the viscosity of the equilibrated dilute solutions is calculated via non-equilibrium molecular dynamics simulation (NEMD). Moreover, the steric hindrance of the copolymers and the solvent molecules capturing in the dilute solution is determined via radial distribution function (RDF) analysis. Helmholtz free energy and the system entropy changes are calculated to evaluate the tendency of the copolymer to the solvent molecules and its dilute solution irregularity, respectively. [Figure not available: see fulltext.].