Time scales and mechanisms of relaxation in the energy landscape of polymer glass under deformation: direct atomistic modeling

High-temperature molecular screening of hybrid polyOAPS-imide networks based on octa(aminophenyl)silsesquioxane for increased thermomechanical resistance

A new family of hybrid hyper-cross-linked thin films based on inorganic polyhedral oligomeric silsesquioxane (POSS) cages covalently bound with short organic imides has recently been developed using interfacial polycondensation followed by high-temperature imidization. These polyPOSS-imide networks were aimed at gas separations under harsh conditions, but the aliphatic arms of the initial POSS precursor, octa(aminopropyl)silsesquioxane, were found to be a weak link. This work investigates the replacement of the aliphatic arm by a phenyl derivative, octa(aminophenyl)silsesquioxane (OAPS). Although this new precursor is expected to be more thermoresistant, it introduces extra degrees of complexity since the functional -NH2 group on the phenyl ring can either be attached at a meta, a para or an ortho position. In order to avoid a costly programme of synthesis and testing, molecular dynamics (MD) simulations have been used to efficiently screen a large number of candidate structures based on mixtures of the three OAPS isomers, the initial POSS and three organic precursors, the PMDA, 6FDA and ODPA dianhydrides. Following cross-linking at room temperature, twenty-two model networks were further relaxed at the imidization temperature and directly tested under harsh conditions at 300 °C. The screening stage included the characterization of their intercage single-links and double-links, which reinforce the structures, and intracage links, which have the opposite effect. Carrying out the cross-linking reactions to completion significantly improved the resistance to isotropic dilation. The initial POSS as well as the flexible 6FDA and ODPA linkers were found to be prone to large deformations, whereas the orthoOAPS, metaOAPS, paraOAPS and the PMDA linker prevented volume dilations. Upon uniaxial tension, the Young's moduli varied in the order paraOAPS < POSS ≈ metaOAPS < orthoOAPS for the inorganic precursors and in the order 6FDA < ODPA < PMDA for the organic precursors. In all cases, the networks based on either orthoOAPS and/or PMDA displayed superior resistance. Nine polyOAPS-imides were further heated up to 400 °C, i.e. closer to the expected degradation, and re-submitted to isotropic dilations and uniaxial tensions. They confirmed the trends found at 300 °C with no signs of structural collapse. Using OAPS as the inorganic precursor thus significantly reinforces the thermoresistance of these hybrid hyper-cross-linked networks.

Thermometer Effect: Origin of the Mixed Alkali Effect in Glass Relaxation

Despite the dramatic increase of viscosity as temperature decreases, some glasses are known to feature room-temperature relaxation. However, the structural origin of this phenomenon-known as the "thermometer effect"-remains unclear. Here, based on accelerated molecular dynamics simulations of alkali silicate glasses, we show that both enthalpy and volume follow stretched exponential decay functions upon relaxation. However, we observe a bifurcation of their stretching exponents, with β=3/5 and 3/7 for enthalpy and volume relaxation, respectively, in agreement with Phillips's topological diffusion-trap model. Based on these results, we demonstrate that the thermometer effect is a manifestation of the mixed alkali effect. We show that relaxation is driven by the existence of stressed local structural instabilities in mixed alkali glasses. This driving force is found to be at a maximum when the concentrations of each alkali atom equal each other, which arises from a balance between the concentration of each alkali atom and the magnitude of the local stress that they experience.

The influence of the dianhydride precursor in hyper-cross-linked hybrid polyPOSS-imide networks

Hybrid organic/inorganic hyper-cross-linked membranes based on imides covalently bonded with polyhedral oligomeric silsesquioxanes (POSS) have recently been developed for gas-separation applications under high pressure and/or temperature conditions. Their molecular sieving capabilities have been shown to depend on the nature of the organic dianhydride precursor. In the present work, realistic molecular models of such polyPOSS-imide films based on the flexible 6FDA dianhydride are compared to those based on the shorter and more rigid PMDA dianhydride. The models creation procedure closely mimicks the mixing, polycondensation and imidization steps of the experimental scheme. The resulting networks are found to be highly heterogeneous in terms of both the number of links (from zero to the maximum possible of eight per POSS cage with an average of four) and their structure (interPOSS, intraPOSS, single-links, double-links) because of the eight-equivalent-arms nature of the POSS precursor. For both dianhydride precursors, crosslinking with POSS and the subsequent imidization step decrease the density, create additional void-space and increase the solubility of the resulting membranes. However, when compared to PMDA, the added flexibility of the central 6FDA bridge leads to a larger thermally-induced dilation of the networks and a larger volume loss per H₂O over the imidization step. With their better ability to redensify and to adapt to the added constraints, the cagecage distances and cage(organic bridge)cage angles in the 6FDA polyPOSS-imides span a larger range than in their PMDA counterparts. In addition, the stiffness of the PMDA moiety results in more unrelaxed free volume remaining trapped in the PMDA polyPOSS-imides upon imidization, and as such, to significantly more open structures with less favourable interactions. As expected from their enhanced flexibility, the thermomechanical properties of the 6FDA networks are slightly lower than those based on PMDA. However, the better mechanical resistance of PMDA over 6FDA does not really become significant before very large volume dilations.

Molecular dynamics simulations of uniaxial deformation of thermoplastic polyimides

The results of atomistic molecular-dynamics simulations of mechanical properties of heterocyclic polymer subjected to uniaxial deformation are reported. A new amorphous thermoplastic polyimide R-BAPO with a repeat unit consisting of dianhydride 1,3-bis-(3',4,-dicarboxyphenoxy)diphenyl (dianhydride R) and diamine 4,4'-bis-(4''-aminophenoxy)diphenyloxide (diamine BAPO) was chosen for the simulations. Our primary goal was to establish the impact of various factors (sample preparation method, molecular mass, and cooling and deformation rates) on the elasticity modulus. In particular, we found that the elasticity modulus was only slightly affected by the degree of equilibration, the molecular mass and the size of the simulation box. This is most likely due to the fact that the main contribution to the elasticity modulus is from processes on scales smaller than the entanglement length. Essentially, our simulations reproduce the logarithmic dependence of the elasticity modulus on cooling and deformation rates, which is normally observed in experiments. With the use of the temperature dependence analysis of the elasticity modulus we determined the flow temperature of R-BAPO to be 580 K in line with the experimental data available. Furthermore, we found that the application of high external pressure to the polymer sample during uniaxial deformation can improve the mechanical properties of the polyimide. Overall, the results of our simulations clearly demonstrate that atomistic molecular-dynamics simulations represent a powerful and accurate tool for studying the mechanical properties of heterocyclic polymers and can therefore be useful for the virtual design of new materials, thereby supporting cost-effective synthesis and experimental research.

Mapping Brittle and Ductile Behaviors of Polymeric Glasses under Large Extension

We have carried out a series of tensile extension tests on the two most common polymer glasses to describe their generic mechanical responses as a function of deformation rate at different temperatures. The essentially defect-free polystyrene and poly(methyl methacrylate) both show remarkable re-entrant failure: being ductile at intermediate rates and showing diminishing toughness at both higher and lower rates. We draw phase diagrams to map out the relationship between brittle-like and yield-like states in terms of temperature, rate, and stress. A coherent understanding of the rich phenomenology requires us to describe in more detail the interplay between the chain network and the primary structure bonded by intersegmental van der Waals forces.

Polystyrene Glasses under Compression: Ductile and Brittle Responses

Polystyrene of different molecular weights and their binary mixtures are studied in terms of their various mechanical responses to uniaxial compression at different temperatures. PS of M_w = 25 kg/mol is completely brittle until it is above its glass transition temperature T_g. In contrast, upon incorporation of a high molecular weight component, PS mixtures turn from barely ductile a few degrees below its T_g to ductile over 40° below T_g. In the upper limit, a PS of M_w = 319 kg/mol yields and undergoes plastic flow, even at T = -70 °C. The observed dependence of mechanical responses on molecular weight and molecular weight distribution can be adequately rationalized by the idea that yielding and plastic compression are caused by chain networking.

Stretched Exponential Relaxation of Glasses at Low Temperature

The question of whether glass continues to relax at low temperature is of fundamental and practical interest. Here, we report a novel atomistic simulation method allowing us to directly access the long-term dynamics of glass relaxation at room temperature. We find that the potential energy relaxation follows a stretched exponential decay, with a stretching exponent β=3/5, as predicted by Phillips's diffusion-trap model. Interestingly, volume relaxation is also found. However, it is not correlated to the energy relaxation, but it is rather a manifestation of the mixed alkali effect.

Roughness and ordering at the interface of oxidized polystyrene and water

For the first time, atomistically detailed molecular dynamics calculations revealed molecular ordering of the water-oxidized atactic polystyrene (aPS) interface. Both ordering of the water molecules and the phenyl rings occur. In addition, the natural roughness of the surface has been simulated and compared to experimental values. The composition of the simulated aPS films is based on spin-coated aPS films that have been oxidized and characterized experimentally. The aPS surfaces are oxidized with ultraviolet-ozone radiation and have been characterized by XPS, AFM, and water contact angle measurements. XPS measurements show that the oxygen content in the sample increases rapidly with exposure and reaches saturation near 24 at. % of oxygen. The molecular dynamics simulations show smoothening of an hydrophobic aPS surface upon transition from vacuum to water. The smoothening decreases with increasing hydrophilicity. The calculations reveal ordering of oxidized phenyl rings for aPS surfaces in water. The order increases with increasing hydrophilicity. Additionally, we investigated the water structure near the aPS-water interface as a function of the surface hydrophilicity. With increasing hydrophilicity, the density of water at the aPS-water interface increases. The water density profile is steeper in the presence of hydrophobic aPS. The water shows an ordered layer near both the hydrophobic and hydrophilic surfaces; the position of this layer shifts toward the interface with increasing hydrophilicity.

Microscopic view of accelerated dynamics in deformed polymer glasses

A molecular level analysis of segmental trajectories obtained from molecular dynamics simulations is used to obtain the full relaxation time spectrum in aging polymer glasses subject to three different deformation protocols. As in experiments, dynamics can be accelerated by several orders of magnitude, and a narrowing of the distribution of relaxation times during creep is directly observed. Additionally, the acceleration factor describing the transformation of the relaxation time distributions is computed and found to obey a universal dependence on the strain, independent of age and deformation protocol.

Fracture in glassy polymers: a molecular modeling perspective

Over the past 25 years, molecular modeling and simulations have provided important insights into the physics of deformation and fracture of glassy polymers. This review presents an overview of key results discussed in the context of experimentally observed polymer behavior. Both atomistic and coarse-grained polymer models have been used in different deformation protocols to study elastic properties, shear yielding, creep, physical aging, strain hardening and crazing. Simulations reproduce most of the macroscopic features of plasticity in polymer glasses such as stress-strain relations and creep response, and reveal information about the underlying atomistic processes. Trends of the shear yield stress with loading conditions, temperature and strain rate, and the atomistic dynamics under load have been systematically explored. Most polymers undergo physical aging, which leads to a history-dependent mechanical response. Simulations of strain hardening and crazing demonstrate the nature of polymer entanglements in the glassy state and the role of local plasticity and provide insight into the origin of fracture toughness of amorphous polymers.

Mechanical rejuvenation and overaging in the soft glassy rheology model

Mechanical rejuvenation and overaging of glasses is investigated through stochastic simulations of the soft glassy rheology (SGR) model. Strain- and stress-controlled deformation cycles for a wide range of loading conditions are analyzed and compared to molecular dynamics simulations of a model polymer glass. Results indicate that deformation causes predominantly rejuvenation, whereas overaging occurs only at very low temperatures, small strains, and for high initial energy states. Although the creep compliance in the SGR model exhibits full aging independent of applied load, large stresses in the nonlinear creep regime cause configurational changes leading to rejuvenation of the relaxation time spectrum probed after a stress cycle. During recovery, however, the rejuvenated state rapidly returns to the original aging trajectory due to collective relaxations of the internal strain.