Glass transition of polymers: atomistic simulation versus experiments

All-atom molecular dynamics simulation of solvent diffusion in an unentangled polystyrene film

The diffusion behavior of low molecular weight solvents within an unentangled polystyrene film below and above its glass transition temperature is investigated. The diffusion behavior in the glassy state exhibits a distinct behavior known as case II or class II diffusion, noticeably diverging from conventional Fickian diffusion observed above the glass transition temperature of the polymer film. In the context of case II diffusion, the primary experimental observation entails the emergence of a well-defined concentration front moving at a constant speed, delineating a swollen, rubbery region from a glassy region within the polymer system. Despite the prevalence of this phenomenon in experimental settings, simulating case II diffusion has posed a significant challenge, primarily due to the computationally intensive nature of the diffusion process. To address this, we have developed an all-atom molecular dynamics simulation approach for the observation of case II diffusion in glassy polymers. This method aims to unravel the intricacies of the diffusion process, providing valuable insights into the dynamic interactions between solvents and the polymer matrix.

Nanoscale Simulation of Plastic Contaminants Migration in Packaging Materials and Potential Leaching into Model Food Fluids

Polymers are the most commonly used packaging materials for nutrition and consumer products. The ever-growing concern over pollution and potential environmental contamination generated from single-use packaging materials has raised safety questions. Polymers used in these materials often contain impurities, including unreacted monomers and small oligomers. The characterization of transport properties, including diffusion and leaching of these molecules, is largely hampered by the long timescales involved in shelf life experiments. In this work, we employ atomistic molecular simulation techniques to explore the main mechanisms involved in the bulk and interfacial transport of monomer molecules from three polymers commonly employed as packaging materials: polyamide-6, polycarbonate, and poly(methyl methacrylate). Our simulations showed that both hopping and continuous diffusion play important roles in inbound monomer diffusion and that solvent-polymer compatibility significantly affects monomer leaching. These results provide rationalization for monomer leaching in model food formulations as well as bulky industry-relevant molecules. Through this molecular-scale characterization, we offer insights to aid in the design of polymer/consumer product interfaces with reduced risk of contamination and longer shelf life.

Toward a Better Understanding of the Poly(glycerol sebacate)-Water Interface

Molecular simulations were conducted to provide a better description of the poly(glycerol sebacate) (PGS)-water interface. The density and the glass-transition temperature as well as their dependencies on the degree of esterification were examined in close connection with the available experimental data. The work of adhesion and water contact angle were calculated as a function of the degree of esterification. A direct correlation was established between the strength of the hydrogen bond network in the interfacial region and the change in the water contact angle with respect to the degree of esterification. The interfacial region was described by local density profiles and orientations of the water molecules.

Dynamics and Internal Structure Evolution during the Glass Transition of the Ethylene-Cyclic Olefin Copolymers: A Molecular Dynamics Simulation

Amorphous ethylene-cyclic olefin copolymers (COCs) which can be used in cell phone lenses and prefilled syringes have attracted increasing attention due to their excellent and tunable thermal properties. In order to better explain the influence of COC microstructure (cyclic olefin types and content) on the glass transition mechanism, we used molecular dynamics (MD) simulations to track the evolution of free volume, diffusion coefficients, atomic mobility, trans conformation probabilities, and characteristic parameters of α-relaxation kinetics during the quenching process. MD results show that for the classic COC E-co-NB (ethylene-norbornene copolymer), an increase in cyclic olefin content from 25 to 50 mol % reduces atomic mobility, limiting the molecular chain movement at higher temperatures and improving Tg. Compared to NB, the more rigid rings in tricyclopentadiene (TCPD) and exo-1,4,4a,9,9a,10-hexahydro-9,10(1',2')-bridged phenylidene-1,4-bridged methylideneanthracene (HBM) have the following effects: (1) reducing the thermal expansion coefficient and overall chain mobility; (2) enhancing the diffusion energy barrier; (3) promoting the formation of local ordered structures; (4) accelerating α-relaxation dynamics at high temperatures and improving the dynamic fragility m. These lead to an upward shift in the temperature region where chain movement is limited and thus improve Tg and high-temperature dimensional stability. In this simulation, the correlation equation between Tg, m, and the microstructural parameters of COCs is established, which is of great significance for the development of COCs with high performance.

Determining glass transition in all-atom acrylic polymeric melt simulations using machine learning

The functionality of many polymeric materials depends on their glass transition temperatures (Tg). In computer simulations, Tg is often calculated from the gradual change in macroscopic properties. Precise determination of this change depends on the fitting protocols. We previously proposed a robust data-driven approach to determine Tg from the molecular dynamics simulation data of a coarse-grained semiflexible polymer model. In contrast to the global macroscopic properties, our method relies on high-resolution microscopic details. Here, we demonstrate the generality of our approach by using various dimensionality reduction and clustering methods and apply it to an atomistic model of acrylic polymers. Our study reveals the explicit contribution of the side chain and backbone residues in influencing the determination of the glass transition temperature.

Amorphous cis-1,4-polybutadiene P-V-T properties from atomistic simulations

CONTEXT

As a result of the diversity of microstructures encountered in cis-1,4-polybutadiene and the variety of measurement methods used, experimental values of variation of glass transition temperature (Tg) with pressure are relatively dispersed. However, atomistic simulations enable access to valuable information for very well-controlled chemistry and structures with a well-defined and systematic acquisition protocol. By varying the temperature and pressure, the specific volume of the melt was computed, yielding results that deviated by only 2% from experimental data. A linear relationship between Tg and pressure was observed, with Tg predicted to be 162 K at zero pressure and a rate of change of Tg with respect to pressure (dTg/dP) of 0.24 K/MPa.

METHOD

The atomistic dilatometry experiments were conducted on a model of amorphous cis-1,4 polybutadiene with an approximate molecular weight of 5400 g/mol using the LAMMPS code and the all-atom forcefield pcff + . The dilatometry process involved cooling and heating at a rate of 9 × 10¹² K/min. The specific volume was calculated by averaging over seven independent configurations for each temperature. The Tait equation was employed to fit the specific volume evolution within the temperature range of 10 to 700 K under different pressures of 0, 60, and 100 MPa.

Data-Driven Identification and Analysis of the Glass Transition in Polymer Melts

Understanding the nature of glass transition, as well as the precise estimation of the glass transition temperature for polymeric materials, remains open questions in both experimental and theoretical polymer sciences. We propose a data-driven approach, which utilizes the high-resolution details accessible through the molecular dynamics simulation and considers the structural information on individual chains. It clearly identifies the glass transition temperature of polymer melts of weakly semiflexible chains. By combining principal component analysis and clustering, we identify the glass transition temperature in the asymptotic limit even from relatively short time trajectories, which just reach into the Rouse-like monomer displacement regime. We demonstrate that fluctuations captured by the principal component analysis reflect the change in a chain's behavior: from conformational rearrangement above to small fluctuations below the glass transition temperature. Our approach is straightforward to apply and should be applicable to other polymeric glass-forming liquids.

Molecular Simulation Strategies for Understanding the Degradation Mechanisms of Acrylic Polymers

Acrylic polymers, commonly used in paints, can degrade over time by several different chemical and physical mechanisms, depending on structure and exposure conditions. While exposure to UV light and temperature results in irreversible chemical damage, acrylic paint surfaces in museums can also accumulate pollutants, such as volatile organic compounds (VOCs) and moisture, that affect their material properties and stability. In this work, we studied the effects of different degradation mechanisms and agents on properties of acrylic polymers found in artists' acrylic paints for the first time using atomistic molecular dynamics simulations. Through the use of enhanced sampling methods, we investigated how pollutants are absorbed into thin acrylic polymer films from the environment around the glass transition temperature. Our simulations suggest that the absorption of VOCs is favorable (-4 to -7 kJ/mol depending on VOCs), and the pollutants can easily diffuse and be emitted back into the environment slightly above glass transition temperature when the polymer is soft. However, typical environmental fluctuations in temperature (<16 °C) can lead for these acrylic polymers to transition to glassy state, in which case the trapped pollutants act as plasticizers and cause a loss of mechanical stability in the material. This type of degradation results in disruption of polymer morphology, which we investigate through calculation of structural and mechanical properties. In addition, we also investigate the effects of chemical damage, such as backbone bond scission and side-chain cross-linking reactions on polymer properties.

A Molecular Dynamics Approach to the Impacts of Oxidative Aging on the Engineering Characteristics of Asphalt

Oxidative aging is an inevitable environmental factor that accelerates asphalt pavement deterioration. This study employed a molecular dynamics simulation to investigate the impact of aging on asphalt cement from the perspectives of thermodynamic properties, and diffusion and adhesion characteristics. Results indicate that aging increased bulk density from 1.008 to 1.081 g/cm3 and cohesive energy density by 15.6%, which was attributed to the promoted molecular polarity and intermolecular attractiveness. The enhanced molecular interactions also reduced molecular mobility, which led to an increase in the glass transition temperature by 30 K, suggesting that aging diminished the resistance of asphalt to thermal cracking. Simulations of the diffusion behaviors across different temperatures demonstrated that the Arrhenius relationship described well the temperature dependence of the diffusion coefficient, and that aging considerably slowed down the diffusion process as represented by Arrhenius prefactor D0, which dropped by 38.2%. The asphalt–aggregate adhesion was assessed using layered models with and without a water interlayer of different thicknesses. The adhesion was enhanced upon aging due to the significantly improved electrostatic interactions at the interface. Evaluation of the residual adhesion with the presence of interfacial water suggested that aging would raise the moisture susceptibility of asphalt pavement. The increase in molecular polarity was considered to be highly responsible for these aging consequences, and was thus further investigated via the electrostatic potential surface and dipole moment.

Curing and Molecular Dynamics Simulation of MXene/Phenolic Epoxy Composites with Different Amine Curing Agent Systems

Herein, the curing kinetics and the glass transition temperature (T_g) of MXene/phenolic epoxy composites with two curing agents, i.e., 4,4-diaminodiphenyl sulfone (DDS) and dicyandiamine (DICY), are systematically investigated using experimental characterization, mathematical modeling and molecular dynamics simulations. The effect of MXene content on an epoxy resin/amine curing agent system is also studied. These results reveal that the MXene/epoxy composites with both curing agent systems conform to the SB(m,n) two-parameter autocatalytic model. The addition of MXene accelerated the curing of the epoxy composite and increased the T_g by about 20 K. In addition, molecular dynamics were used to simulate the T_g of the cross-linked MXene/epoxy composites and to analyze microstructural features such as the free volume fraction (FFV). The simulation results show that the introduction of MXene improves the T_g and FFV of the simulated system. This is because the introduction of MXene restricts the movement of the epoxy/curing agent system. The conclusions are in good agreement with the experimental results.

The Unusual Conductivity of Na+ in PEO-Based Statistical Copolymer Solid Electrolytes: When Less Means More

The low conductivity of Na⁺ electrolytes in solid polymer electrolytes (SPEs) curtails the development of Na polymer batteries. In this study, NaClO₄ (3-24 wt %, 90-9:1 O:Na) is dissolved in statistical copolymers of ethylene oxide (EO) and propylene oxide (PO) (0-20 mol %). Remarkably, the conductivity of these SPEs increases as the concentration of Na⁺ decreases, thus departing from the usual Nernstian behavior. Using a combination of calorimetric measurements and molecular dynamic simulations, this unusual phenomenon is attributed to the presence of physical cross-links generated by Na⁺ . As a result, polymers containing a low salt concentration (3 wt %) display a drastically enhanced ionic conductivity (up to 0.2 10^-4  S cm^-1 at 25 °C), thus paving the way for the design of all-solid-state PEO-based sodium batteries operational at room temperature.

Atomistic simulation of volumetric properties of epoxy networks: effect of monomer length

Properties of epoxy thermosets can be varied broadly to suit design requirements by altering the chemistry of the component agents. Atomistically-detailed molecular dynamics simulations are well-suited for molecular insight into the structure-property relationship for a rational tailoring of the chemistry. Since the macroscopic properties of interest for applications emerge hierarchically from molecular-scale chemical interactions, seamless integration of experiment, computation, and theory is of great interest. Recently, a Specific Volume-Cooling Rate analysis protocol was successfully developed to quantitatively compare the volumetric properties of an epoxy network model with experimental results in the literature, in spite of the nine orders of magnitude mismatch in the accessible time-scales. Here, we extend the application of the method for two epoxy networks in the same class of chemistry but whose monomers have a higher number of repeating units compared to the previous one for validating the generality of our approach. We observed that atomistic simulations are able to predict the experimental temperature trend of the specific volume within 0.4% for both these networks. Using the William-Landel-Ferry equation to account for rate differences, we also see good agreement between the computational and experimental values of the glass transition temperature.

Acrylic Paints: An Atomistic View of Polymer Structure and Effects of Environmental Pollutants

Most of the artwork and cultural heritage objects are stored in museums under conditions that are difficult to monitor. While advanced technologies aim to control and prevent the degradation of cultural heritage objects in line with preventive conservation measures, there is much to be learned in terms of the physical processes that lead to the degradation of the synthetic polymers that form the basis of acrylic paints largely used in contemporary art. In museums, stored objects are often exposed to temperature and relative humidity fluctuations as well as airborne pollutants such as volatile organic compounds (VOCs). The glass transition of acrylic paints is below room temperature; while low temperatures may cause cracking, at high temperatures the sticky surface of the paint becomes vulnerable to pollutants. Here we develop fully atomistic models to understand the structure of two types of acrylic copolymers and their interactions with VOCs and water. The structure and properties of acrylic copolymers are slighlty modified by incorporation of a monomer with a longer side chain. With favorable solvation free energies, once absorbed, VOCs and water interact with the polymer side chains to form hydrogen bonds. The cagelike structure of the polymers prevents the VOCs and water to diffuse freely below the glass transition temperature. In addition, our model forms the foundation for developing mesoscopic and continuum models that will allow us to access longer time and length scales to further our understanding of the degradation of artwork.

Extending the timescale of molecular simulations by using time-temperature superposition: rheology of ionic liquids

Molecular dynamics simulations are used to determine the temperature dependence of the dynamic and rheological properties of a model imidazolium-based ionic liquid (IL). The simulation results for the volumetric properties of the IL are in good agreement with the experimental results. The temperature dependence of the diffusion coefficient of anions and cations follows the Vogel-Fulcher-Tammann equation over the range of the temperatures studied. The shear viscosity of the IL shows a Newtonian plateau at low shear rates and shear-thinning behavior at high shear rates. The dynamic modulus values indicate that the IL behaves like a viscous liquid at high temperatures and low frequencies, while its viscoelastic response becomes similar to that of an elastic solid at low temperatures and high frequencies. Using the time-temperature superposition (TTS) principle, the dynamic moduli, shear viscosity, and mean squared displacement of cations and anions in the diffusive regime can be collapsed onto master curves by applying a single set of shift factors. Due to the large mismatch in the timescale investigated by the atomistically detailed simulations and experiments, the glass transition temperature predicted in simulations shifts to higher values. When this timescale mismatch is accounted for by using appropriate shift factors, the master curves of the dynamic moduli obtained in simulations closely match those obtained in experiments. This result demonstrates the exciting ability of TTS to overcome the large timescale disparity between simulations and experiments which will enable the use of molecular simulations for quantitatively predicting the rheological property values at frequencies of practical interest.

Impact of Amorphization Methods on the Physicochemical Properties of Amorphous Lactulose

The influence of the amorphization technique on the physicochemical properties of amorphous lactulose was investigated. Four different amorphization techniques were used: quenching of the melt, milling, spray-drying, and freeze-drying, and amorphous samples were analyzed by differential scanning calorimetry, NMR spectroscopy, and powder X-ray diffraction analysis. Special attention was paid to the tautomeric composition and to the glass transition of amorphized materials. It was found that the tautomeric composition of the starting physical state (crystal, liquid, or solution) is preserved during the amorphization process and has a strong repercussion on the glass transition of the material. The correlation between these two properties as well as the plasticizing effect of the different tautomers was clarified by molecular dynamics simulations.

Local dynamics within the glass transition domain

The glass transition of an amorphous material is a fundamental property characterized by an abrupt change in viscosity. Its very knowledge was a conundrum as no satisfying theory existed at the molecular level. We herein relate this complex phenomenon to events occurring at the molecular scale. By studying conformational transitions in the carbon-chain polymer of polyethylene, we clearly establish a relation between local dynamics and the classical dihedral potential energy diagram of a carbon-carbon bond. This methodology is applied to a carbon-chain polymer with a side-group, polystyrene. A direct link is proved between activation energy and glass transition temperature. This work thus provides the cornerstone for linking molecular structure to macroscopic polymer properties, and in particular, the glass transition temperature.

Model of the DGEBA-EDA Epoxy Polymer: Experiments and Simulation Using Classical Molecular Dynamics

Polyepoxy samples are synthesized from diglycidylether of bisphenol A (DGEBA) and ethylene diamine (EDA) monomers at a stoichiometric ratio of 2 DGEBA : 1 EDA in model conditions in order to promote a high degree of polymerization and a low density of defects and to try to approach the ideal models obtained by simulation. A slow polymerization (>24 h at ambient temperature) and a postcuring achieved in an inert atmosphere lead to a conversion degree of 92±2% and a midpoint glass transition temperature of 391±1 K. In parallel, a model is created with a multistep cross-linking procedure. In this work, all-atom molecular dynamics (MD) simulations are performed with LAMMPS and the GAFF 1.8 force field. In the initial liquid mixture of reactants (600 molecules), proper mixing is demonstrated by the calculation of the partial radial distribution functions (RDF), which show a minimum intermolecular distance of 2.8 Å and similar distributions for EDA-EDA, DGEBA-DGEBA, and DGEBA-EDA molecules in the simulation boxes. Then, in alternation with MD equilibrations, cross-linking is performed on frozen configurations by creating covalent bonds between reactive pairs within a reaction radius of 3 Å. The resulting boxes show conversion rates of 90-93% and densities close to the experimental value. Finally, a cooling ramp from 700 K to 25 K is applied in order to monitor the specific volume and the coefficient of volumetric thermal expansion (CVTE) of the polymer and to derive the glass transition temperature. Experimental thermomechanical analyses (TMA) compares well with simulations for both the specific volume and the CVTE evolutions with temperature. Whereas the uncertainty remains high with the fitting procedure used, we calculate a glass transition temperature of 390±8 K which compares very well with the experimental values (391±1 K from DSC and 380 K from TMA).

Computing thermomechanical properties of dry homopolymers used as raw materials for formulation of biomedical hydrogels

Different static properties have been calculated with COMPASS force field for polyacrylamide, poly(2-hydroxyethylacrylate) (HEA), poly(2-hydroxyethylmethacrylate) (HEMA), poly(glycidylmethacrylate) (GMA), polyethylene glycol (PEG), and poly(2,2,2-trifluoroethylmethacrylate) (TFEM). For each polymers, the calculated values were averaged on five equilibrated configurations of amorphous cell composed of one atactic chain containing 100 repeat units. The ranking obtained from the densities calculated at 300 K is TFEM > HEA ≈ xpolycrylamide > HEMA ≈ GMA > PEG. Concerning the glass transition temperature we have obtained polyacrylamide > HEMA ≈ GMA ≈ HEA > PEG, and polyacrylamide > HEMA ≈ HEA > GMA ≈ PEG > TFEM for the bulk modulus. The calculated results, when available, have been compared with experimental data coming from literature.

Effect of Interfacial Bonding on Interphase Properties in SiO2/Epoxy Nanocomposite: A Molecular Dynamics Simulation Study

Atomistic molecular dynamics simulations have been performed to explore the effect of interfacial bonding on the interphase properties of a nanocomposite system that consists of a silica nanoparticle and the highly cross-linked epoxy matrix. For the structural properties, results show that interfacial covalent bonding can broaden the interphase region by increasing the radial effect range of fluctuated mass density and oriented chains, as well as strengthen the interphase region by improving the thermal stability of interfacial van der Waals excluded volume and reducing the proportion of cis conformers of epoxy segments. The improved thermal stability of the interphase region in the covalently bonded model results in an increase of ∼21 K in the glass transition temperature (Tg) compared to that of the pure epoxy. It is also found that interfacial covalent bonding mainly restricts the volume thermal expansion of the model at temperatures near or larger than Tg. Furthermore, investigations from mean-square displacement and fraction of immobile atoms point out that interfacial covalent and noncovalent bonding induces lower and higher mobility of interphase atoms than that of the pure epoxy, respectively. The obtained critical interfacial bonding ratio when the interphase and matrix atoms have the same mobility is 5.8%. These results demonstrate that the glass transitions of the interphase and matrix will be asynchronous when the interfacial bonding ratio is not 5.8%. Specifically, the interphase region will trigger the glass transition of the matrix when the ratio is larger than 5.8%, whereas it restrains the glass transition of the matrix when the ratio is smaller than 5.8%.