Interfacial and topological effects on the glass transition in free-standing polystyrene films

The glass transition and enthalpy recovery of polystyrene nanorods using Flash differential scanning calorimetry

The glass transition (Tg) behavior and enthalpy recovery of polystyrene nanorods within an anodic aluminum oxide (AAO) template (supported nanorods) and after removal from AAO (unsupported nanorods) is studied using Flash differential scanning calorimetry. Tg is found to be depressed relative to the bulk by 20 ± 2 K for 20 nm-diameter unsupported polystyrene (PS) nanorods at the slowest cooling rate and by 9 ± 1 K for 55 nm-diameter rods. On the other hand, bulk-like behavior is observed in the case of unsupported 350 nm-diameter nanorods and for all supported rods in AAO. The size-dependent Tg behavior of the PS unsupported nanorods compares well with results for ultrathin films when scaled using the volume/surface ratio. Enthalpy recovery was also studied for the 20 and 350 nm unsupported nanorods with evolution toward equilibrium found to be linear with logarithmic time. The rate of enthalpy recovery for the 350 nm rods was similar to that for the bulk, whereas the rate of recovery was enhanced for the 20 nm rods for down-jump sizes larger than 17 K. A relaxation map summarizes the behavior of the nanorods relative to the bulk and relative to that for the 20 nm-thick ultrathin film. Interestingly, the fragility of the 20 nm-diameter nanorod and the 20 nm ultrathin film are identical within the error of measurements, and when plotted vs departure from Tg (i.e., T - Tg), the relaxation maps of the two samples are identical in spite of the fact that the Tg is depressed 8 K more in the nanorod sample.

Surface Characteristics of Poly(alkyl methacrylate)s from Molecular Dynamics Simulations Using All-Atom Force Field

Molecular dynamics (MD) simulations of melt films of poly(alkyl methacrylate)s (PAMAs) with methyl, ethyl, and n-butyl substituents, respectively, have been performed using an all-atom model to investigate their surface and thin film properties. The applied all-atom force fields predict the bulk densities of PAMAs in good agreement with experiments. Moreover, predictions of the surface tensions of PMMA, PEMA, and Pn-BMA melts are in reasonably good agreement with experiments. The density profiles and orientational-order parameters of chain segments show atomic-scale characteristics in the air/polymer interfacial region. In the surface region, the backbone segments of PAMAs form a well-defined layer structure with the chain vectors oriented parallel to the surface, while the ester side-chains strongly segregate to the surface region and show perpendicular orientation to the surface, with the most pronounced surface segregation noted for Pn-BMA. Such surface segregations of chain segments make it difficult to apply a simple relationship between the cohesive energy density and the surface tension of polymers, for example, and should be taken into account in relating the surface/thin film characteristics to the bulk properties of polymers in general.

Conformation and dynamics of ring polymers under symmetric thin film confinement

Understanding the structure and dynamics of polymers under confinement has been of widespread interest, and one class of polymers that have received comparatively little attention under confinement is that of ring polymers. The properties of non-concatenated ring polymers can also be important in biological fields because ring polymers have been proven to be a good model to study DNA organization in the cell nucleus. From our previous study, linear polymers in a cylindrically confined polymer melt were found to segregate from each other as a result of the strong correlation hole effect that is enhanced by the confining surfaces. By comparison, our subsequent study of linear polymers in confined thin films at similar levels of confinements found only the onset of segregation. In this study, we use molecular dynamics simulation to investigate the chain conformations and dynamics of ring polymers under planar (1D) confinement as a function of film thickness. Our results show that conformations of ring polymers are similar to the linear polymers under planar confinement, except that ring polymers are less compressed in the direction normal to the walls. While we find that the correlation hole effect is enhanced under confinement, it is not as pronounced as the linear polymers under 2D confinement. Finally, we show that chain dynamics far above Tg are primarily affected by the friction from walls based on the monomeric friction coefficient we get from the Rouse mode analysis.

Interface Characteristics of Neat Melts and Binary Mixtures of Polyethylenes from Atomistic Molecular Dynamics Simulations

Molecular dynamics simulations of free-standing thin films of neat melts of polyethylene (PE) chains up to C150H302 and their binary mixtures with n-C13H28 are performed employing a united atom model. We estimate the surface tension values of PE melts from the atomic virial tensor over a range of temperatures, which are in good agreement with experimental results. Compared with short n-alkane systems, there is an enhanced surface segregation of methyl chain ends in longer PE chains. Moreover, the methyl groups become more segregated in the surface region with decreasing temperature, leading to the conclusion that the surface-segregation of methyl chain ends mainly arises from the enthalpic origin attributed to the lower cohesive energy density of terminal methyl groups. In the mixtures of two different chain lengths, the shorter chains are more likely to be found in the surface region, and this molecular segregation in moderately asymmetric mixtures in the chain length (C13H28 + C44H90) is dominated by the enthalpic effect of methyl chain ends. Such molecular segregation is further enhanced and dominated by the entropic effect of conformational constraints in the surface for the highly asymmetric mixtures containing long polymer chains (C13H28 + C150H3020). The estimated surface tension values of the mixtures are consistent with the observed molecular segregation characteristics. Despite this molecular segregation, the normalized density of methyl chain ends of the longer chain is more strongly enhanced, as compared with the all-segment density of the longer chain itself, in the surface region of melt mixtures. In addition, the molecular segregation results in higher order parameter of the shorter-chain segments at the surface and deeper persistence of surface-induced segmental order into the film for the longer chains, as compared with those in neat melt films.

Theory of Spatial Gradients of Relaxation, Vitrification Temperature and Fragility of Glass-Forming Polymer Liquids Near Solid Substrates

We employ a new force-level statistical mechanical theory to predict spatial gradients of the structural relaxation time and T_g of polymer liquids near microscopically rough and smooth hard surfaces and contrast the results with vapor interface systems. Repulsive rough (smooth) surfaces induce large slowing down (modest speeding up) of the relaxation time compared to the bulk. Nevertheless, a remarkable degree of universality of distinctive dynamical behaviors is predicted for different polymer chemistries and all interfaces, including a double exponential form of the alpha time gradient, power law decoupling of the relaxation time from its bulk value with exponential spatial variation of the exponent, exponential spatial gradient of T_g, weak dependence of normalized T_g gradients on vitrification criterion, and near linear growth with cooling of the slowed down layer thickness near a rough hard interface. The results appear consistent with simulations and experiments, and multiple testable predictions are made.

Contributions of force field interaction forms to Green-Kubo viscosity integral in n-alkane case

Decades of molecular simulation history proved that the Green-Kubo method for shear viscosity converges without any problems in atomic and simple molecular liquids, unlike liquids with high values of viscosity. In the case of highly viscous liquids, the time decomposition method was developed in 2015 by Maginn and co-authors [Y. Zhang, A. Otani, and E. J. Maginn, J. Chem. Theory Comput. 11, 3537-3546 (2015)] which allows us to improve the convergence of the Green-Kubo integral. In this paper, the contributions of intramolecular and intermolecular force field parts to the viscosity integral are discovered to gain the understanding of the Green-Kubo method. The n-alkanes from n-ethane to n-pentane at 330 K in the optimized potentials for liquid simulations-all atom force field are used as reference models. The dependencies of these contributions and decay times of the corresponding correlation functions on the chain length are observed. The nonequilibrium simulations are carried out to verify the Green-Kubo results. The obtained values of viscosity are compared with experimental data.

A simulation study on the glass transition behavior and relevant segmental dynamics in free-standing polymer nanocomposite films

In polymer/nanoparticle composite (PNC) thin films, polymer chains experience strong confinement effects not only at the free surface area but also from nanoparticles (NPs). In this work, the influence of NP-polymer interaction and NP distribution on the polymer segmental dynamics and the glass transition behavior of PNC free-standing films are investigated through molecular dynamics simulations. We demonstrate that NPs will migrate to the film surface area and form an NP-concentrated layer when NP-polymer interactions are weak, while NPs are well dispersed in the bulk region when NP-polymer interactions are strong. In both cases, we find increases in the glass transition temperature Tg compared with the pure film without NPs, although with a different degree. The weakly interacting system has the same Tg as the pure bulk system without NPs. The NP layer formed at the surface area reduces both the mobility of the surface polymer beads and the mobility gradient in the film normal direction (MGFND), therefore resulting in an increase in the Tg which highlights the vital role of the mobile surface layer. In contrast, the NPs in the bulk region enlarge the MGFND. NPs have opposite influences on the polymer bead dynamic anisotropy when they interact weakly or strongly with polymers, weakened for the former and enhanced for the latter. These findings offer a clear picture of the segmental dynamics and glass transition behavior in free-standing PNC films with different NP-polymer interaction strengths. We hope these results will be helpful for the property design of related materials.

Theory of the spatial transfer of interface-nucleated changes of dynamical constraints and its consequences in glass-forming films

We formulate a new theory for how caging constraints in glass-forming liquids at a surface or interface are modified and then spatially transferred, in a layer-by-layer bootstrapped manner, into the film interior in the context of the dynamic free energy concept of the Nonlinear Langevin Equation (NLE) theory approach. The dynamic free energy at any mean location (cage center) involves contributions from two adjacent layers where confining forces are not the same. At the most fundamental level of the theory, the caging component of the dynamic free energy varies essentially exponentially with distance from the interface, saturating deep enough into the film with a correlation length of modest size and weak sensitivity to the thermodynamic state. This imparts a roughly exponential spatial variation of all the key features of the dynamic free energy required to compute gradients of dynamical quantities including the localization length, jump distance, cage barrier, collective elastic barrier, and alpha relaxation time. The spatial gradients are entirely of dynamical, not structural or thermodynamic, origin. The theory is implemented for the hard sphere fluid and diverse interfaces which can be a vapor, a rough pinned particle solid, a vibrating (softened) pinned particle solid, or a smooth hard wall. Their basic description at the level of the spatially heterogeneous dynamic free energy is identical, with the crucial difference arising from the first layer where dynamical constraints can be weakened, softened, or hardly changed depending on the specific interface. Numerical calculations establish the spatial dependence and fluid volume fraction sensitivity of the key dynamical property gradients for five different model interfaces. A comparison of the theoretical predictions for the dynamic localization length and glassy modulus with simulations and experiments for systems with a vapor interface reveals good agreement. The present advance sets the stage for using the Elastically Collective NLE theory to make quantitative predictions for the alpha relaxation time gradient, decoupling phenomena, T_g gradient, and many film-averaged properties of both model and experimental (colloids, molecules, and polymers) systems with diverse interfaces and chemical makeup.

Influence of wall heterogeneity on nanoscopically confined polymers

We investigate via molecular dynamics simulations the behaviour of a polymer melt confined between surfaces with increasing spatial correlation (patchiness) of weakly and strongly interacting sites. Beyond a critical patchiness, we find a dramatic dynamic decoupling, characterized by a steep growth of the longest relaxation time and a constant diffusion coefficient. This arises from dynamic heterogeneities induced by the walls in the adjacent polymer layers, leading to the coexistence of fast and slow chain populations. Structural variations are also present, but they are not easy to detect. Our work opens the way to a better understanding of adhesion, friction, rubber reinforcement by fillers, and many other open issues involving the dynamics of polymeric materials on rough, chemically heterogeneous and possibly "dirty" surfaces.