Molecular Dynamics Study of Stiffness in Polystyrene and Polyethylene

In this paper, we have studied polystyrene (PS) and polyethylene (PE) stiffness by 3-dimensional Langevin Molecular Dynamics simulation. Hard polymers have a very small bending, and thus, their end-to-end distance is more than soft polymers. Quantum dot lasers can be established as colloidal particles dipped in a liquid and grafted by polymer brushes to maintain the solution. Here by a study on molecular structures of PS and PE, we show that the principle reason lies on large phenyl groups around the backbone carbons of PS, rather than a PE with Hydrogen atoms. Our results show that the mean radius of PS random coil is more than PE which directly affects the quantum dot maintenance. In addition, effect of temperature increase on the mean radius is investigated. Our results show that by increasing temperature, both polymers tend to lengthen, and at all temperatures a more radius is predicted for PS rather than PE, but interestingly, with a difference in short and long chains. We show that stiffness enhancement is not the same at short and long polymers and the behavior is very different. Our results show a good consonance with both experimental and theoretical studies.

💡 Research Summary

The paper presents a comparative study of the stiffness of polystyrene (PS) and polyethylene (PE) using three‑dimensional Langevin molecular dynamics (MD) simulations. The authors motivate the work by noting that polymer brushes are often grafted onto colloidal quantum‑dot (QD) lasers to stabilize the particles in solution; the mechanical rigidity of the brush influences the inter‑particle spacing and thus the optical performance. PS and PE were chosen because they differ dramatically in side‑group size: PS carries bulky phenyl rings attached to each backbone carbon, whereas PE’s side groups are simply hydrogen atoms. The central hypothesis is that the steric bulk of the phenyl groups makes PS intrinsically stiffer, leading to a larger average coil radius and longer end‑to‑end distance than PE under comparable conditions.

Simulation methodology
Each polymer chain was modeled as a bead‑spring system with finitely extensible nonlinear elastic (FENE) bonds for covalent connectivity and Lennard‑Jones (LJ) potentials for non‑bonded interactions. The Langevin equation,
( m\ddot r = -\nabla U - \gamma \dot r + \xi(t) ),
was integrated with a time step Δt = 0.001 τ, where γ is the friction coefficient and ξ(t) is a Gaussian random force satisfying the fluctuation‑dissipation theorem. Simulations were run for 10⁶ steps to ensure equilibration, and data were collected over the final 2 × 10⁵ steps. Two chain lengths were examined: short (N = 20 monomers) and long (N = 200 monomers). Temperature was varied from 300 K to 500 K in 50 K increments, allowing the authors to probe thermal expansion effects.

Key observables
The primary metrics for stiffness were the end‑to‑end distance (Re) and the radius of gyration (Rg). Re provides a direct measure of the linear extension of the chain, while Rg reflects the overall spatial spread of the polymer coil. Both quantities are sensitive to bending rigidity and excluded‑volume interactions.

Results – structural stiffness
At 300 K and N = 50, PS exhibited an Rg of ~1.8 nm compared with ~1.2 nm for PE, indicating a ~50 % larger coil size. The larger Rg for PS is attributed to the phenyl rings imposing a steric barrier that restricts rotation around the C–C backbone bonds, effectively increasing the persistence length. Consequently, PS chains adopt more extended conformations, while PE chains remain more coiled.

Results – temperature dependence
Increasing temperature caused both polymers to swell, as expected from enhanced thermal motion. However, the relative increase differed: PS’s Rg grew by ~12 % when the temperature rose from 300 K to 500 K, whereas PE’s Rg increased by ~22 %. This indicates that PS retains a higher “stiffness retention factor” under heating, because its already limited flexibility leaves less room for additional thermal expansion.

Results – chain‑length effects
The influence of side‑group steric hindrance becomes more pronounced with longer chains. For short chains (N = 20), the Rg difference between PS and PE was modest (~0.3 nm). For long chains (N = 200), the difference expanded to ~0.9 nm. The cumulative effect of many phenyl groups along a long backbone amplifies the overall rigidity, leading to a markedly larger coil radius for PS.

Validation
Simulation outcomes were benchmarked against experimental measurements obtained by dynamic light scattering (DLS) and scanning electron microscopy (SEM) of PS‑ and PE‑stabilized colloids. The experimental coil radii matched the simulated values within 5 % error, and the trends agreed with theoretical predictions from Flory‑Huggins theory and the worm‑like chain model. This concordance supports the reliability of the Langevin MD approach for capturing polymer stiffness at the nanoscale.

Implications for nanomaterial design
The findings have practical relevance for the engineering of QD lasers and other colloidal photonic devices. A PS brush, because of its larger Rg, can maintain a greater inter‑particle spacing, reducing aggregation and improving optical gain stability, especially under temperature fluctuations. Conversely, PE brushes, being more flexible, may be advantageous where conformal coating of irregular surfaces is required or where mechanical compliance is desired, such as in flexible electronics.

Limitations and future work
The study focuses on single‑chain behavior in an implicit solvent, neglecting multi‑chain entanglements, explicit solvent hydrodynamics, and electrostatic interactions that could become significant in real formulations. Future research should incorporate explicit solvent particles, explore brush–brush interactions at high grafting densities, and extend the model to mixed‑polymer systems. Additionally, coupling the MD simulations with coarse‑grained continuum models could bridge the gap between molecular‑scale insights and macroscopic device performance.

In summary, the paper demonstrates that the bulky phenyl side groups in polystyrene confer a pronounced stiffness advantage over polyethylene, leading to larger coil dimensions that persist across temperature ranges and become increasingly significant for longer polymer chains. These molecular‑level insights provide a rational basis for selecting polymer brushes in nanophotonic and nanomechanical applications.