Validation of Force Fields of Rubber through Glass-Transition Temperature Calculation by Microsecond Atomic-Scale Molecular Dynamics Simulation

Glass Transition Temperatures and Thermal Conductivities of Polybutadiene Crosslinked with Randomly Distributed Sulfur Chains Using Molecular Dynamic Simulation

The thermal conductivities and glass transition temperatures of polybutadiene crosslinked with randomly distributed sulfur chains having different lengths from mono-sulfur (S1) to octa-sulfur (S8) were investigated. The thermal conductivities of the related models as a function of the heat flux autocorrelation function, applying an equilibrium molecular dynamic (EMD) simulation and the Green-Kubo method, were studied for a wide range of temperatures. The influence of the length of sulfur chains, degree of crosslinking, and molar mass of the crosslinker on the glass transition temperature and final values of thermal conductivities were studied. First, the degree of crosslinking is considered constant for the eight simulation models, from mono-sulfur (S1) to octa-sulfur (S8), while the molar mass of the sulfur is increases. The results show that the thermal conductivities of the crosslinked structure decrease with increasing temperature for each model. Moreover, by increasing the lengths of the sulfur chains and the molar weight of the crosslinker, thermal conductivity increases at a constant temperature. The MD simulation demonstrates that the glass transition temperature and density of the crosslinked structure enhance as the length of the sulfur chains and molar mass of the sulfur increase. Second, the molar weight of sulfur is considered constant in these eight models; therefore, the degree of crosslinking decreases with the increase in the lengths of the sulfur chains. The results show that the thermal conductivities of the crosslinked structure decrease with the increase in the temperature for each model. Moreover, by increasing the lengths of sulfur chains and thus decreasing the degree of crosslinking, the trend in changes in thermal conductivities are almost the same for all of these models, so thermal conductivity is constant for a specific temperature. In addition, the glass transition temperature and density of the crosslinked structure decrease.

Structural and dynamical properties of thermoplastic polyurethane/fullerene nanocomposites: a molecular dynamics simulations study

In this work, the structural and dynamical properties of thermoplastic polyurethane (TPU)/fullerene (C60) nanocomposites are investigated using atomistic molecular dynamics simulations, focusing on the glass transition, thermal expansion, polymer mobility, polymer-C60 interactions, and diffusion behavior of C60. The results show a slight increase in the glass transition temperature (Tg) with increasing C60 weight fraction (wt%), attributed to hindered polymer dynamics, and a remarkable reduction in the coefficient of thermal expansion above Tg. Results of the mean squared displacement and the time decay of bond-reorientation autocorrelation indicate that the mobility of TPU hard segments is more restricted than that of soft segments, owing to the electrostatic attractions and the π-π stacking between isocyanate groups and C60 molecules. Analysis of TPU-C60 interaction energy reveals that the electrostatic interactions are weakened with an increase in the C60 wt%, while the van der Waals contributions become more significant due to the TPU-C60 interfacial characteristics. Further analysis shows that the translational and rotational diffusion of C60 are both increasingly suppressed with the increase of C60 wt%, indicating a violation of Stokes-Einstein (SE) and Stokes-Einstein-Debye (SED) relations, presumably due to the polymer chain-mediated hydrodynamic interactions arising from chain bridges between neighboring C60 particles. This is highlighted by a stronger decoupling of translational-rotational diffusion and a lower ratio of translational-rotational diffusion coefficient (DT/DR) with increasing C60 wt%. This work elucidates an atomistic understanding of the structure and properties of polymer/C60 nanocomposites.

The effect of chemical constitution on polyisoprene dynamics

Polyisoprene (PI) melts have been studied, with most reports focusing on systems with high 1,4-cis content. In contrast, 1,4-trans PI homopolymers or random copolymers have seldom been examined, despite a handful of investigations suggesting a distinct dynamic behavior. Herein, we employ all-atom simulations to investigate the effect of chemical architecture on the dynamics of cis and trans-PI homopolymers, as well as copolymers. We examine the thermodynamic, conformational, and structural properties of the polymers and validate the performance of the models. We probe chain dynamics, revealing that cis-PI presents accelerated translation and reorientation modes relative to trans as recorded by the mean square displacement of the chain center-of-mass as well as by the characteristic times of the lower modes in a Rouse analysis. Interestingly, progressing to higher modes, we observe a reversal with trans units exhibiting faster dynamics. This was further confirmed by calculations of local carbon-hydrogen vector reorientation dynamics, which offer a microscopic view of segmental mobility. To obtain insight into the simulation trajectories, we evaluate the intermediate incoherent scattering function that supports a temperature-dependent crossover in relative mobility that extends over separations beyond the Kuhn-length level. Finally, we analyzed the role of non-Gaussian displacements, which demonstrate that cis-PI exhibits increased heterogeneity in dynamics over short-timescales in contrast to trans-PI, where deviations persist over times extending to terminal dynamics. Our all-atom simulations provide a fundamental understanding of PI dynamics and the impact of microstructure while providing important data for the design and optimization of PI-based materials.

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.

Specific Ion and Electric Field Controlled Diverse Ion Distribution and Electroosmotic Transport in a Polyelectrolyte Brush Grafted Nanochannel

Controlling ion distribution inside a charged nanochannel is central to using such channels in diverse applications. Here, we show the possibility of using a charged polyelectrolyte (PE) brush-grafted nanochannel for triggering diverse nanoscopic ion distribution and nanofluidic electroosmotic transport by controlling the valence and size of the counterions (that screen the charges of the PE brushes) and the strength of an externally applied axial electric field. We atomistically simulate separate cases of fully charged polyacrylic acid (PAA) brush functionalized nanochannels with Na⁺, Cs⁺, Ca²⁺, Ba²⁺, and Y³⁺ counterions screening the PE charges. Four key findings emerge from our simulations. First, we find that the counterions with a greater valence and a smaller size prefer to remain localized inside the brush layer. Second, for the case where there is an added chloride salt with the same cation (as the screening counterions), there are more coions (Cl^- ions) in the brush-free bulk than counterions (for counterions Na⁺, Ca²⁺, Ba²⁺, Y³⁺): this is a manifestation of the overscreening (OS) of the PE brush layer. Contrastingly, the number of Cs⁺ ions remain higher than the Cl^- ions inside the brush-free bulk, ensuring that there is no OS effect for this case. Third, large applied electric field enables a few Na⁺, Cs⁺, and Ba²⁺ counterions to leave the brush layer and to go to the bulk: this makes the OS of the PE brush layer disappear for the cases of PE brushes being screened by the Na⁺ and Ba²⁺ ions. On the other hand, no such electric-field-mediated disappearance of OS is observed for the cases of Ca²⁺ and Y³⁺ screening counterions; we attribute this to the firm attachment of these counterions to the negatively charged monomers. Free energy associated with a counterion binding to a PE chain corroborates this diversity in the counterion-specific response to the applied electric field. Finally, we demonstrate that such diverse ion distributions, along with specific electric-field-strength-dependent ion properties, lead to (1) electroosmotic (EOS) transport in nanochannels grafted with PAA brushes screened with Cs⁺ ions to be always counterion dominated, (2) EOS transport in nanochannels grafted with PAA brushes screened with Ca²⁺ and Y³⁺ ions to be always coion-dominated, and (3) EOS transport in nanochannels grafted with PAA brushes screened with Na⁺ and Ba²⁺ ions to be coion dominated for smaller electric fields and counterion dominated for larger electric fields.

Nanocomposite of Fullerenes and Natural Rubbers: MARTINI Force Field Molecular Dynamics Simulations

The mechanical properties of natural rubber (NR) composites depend on many factors, including the filler loading, filler size, filler dispersion, and filler-rubber interfacial interactions. Thus, NR composites with nano-sized fillers have attracted a great deal of attention for improving properties such as stiffness, chemical resistance, and high wear resistance. Here, a coarse-grained (CG) model based on the MARTINI force field version 2.1 has been developed and deployed for simulations of cis-1,4-polyisoprene (cis-PI). The model shows qualitative and quantitative agreement with the experiments and atomistic simulations. Interestingly, only a 0.5% difference with respect to the experimental result of the glass transition temperature (T_g) of the cis-PI in the melts was observed. In addition, the mechanical and thermodynamical properties of the cis-PI-fullerene(C₆₀) composites were investigated. Coarse-grained molecular dynamics (MD) simulations of cis-PI-C₆₀ composites with varying fullerene concentrations (0-32 parts per hundred of rubber; phr) were performed over 200 microseconds. The structural, mechanical, and thermal properties of the composites were determined. The density, bulk modulus, thermal expansion, heat capacity, and T_g of the NR composites were found to increase with increasing C₆₀ concentration. The presence of C₆₀ resulted in a slight increasing of the end-to-end distance and radius of the gyration of the cis-PI chains. The contribution of C₆₀ and cis-PI interfacial interactions led to an enhancement of the bulk moduli of the composites. This model should be helpful in the investigations and design of effective fillers of NR-C₆₀ composites for improving their properties.

The structure of CO2 and CH4 at the interface of a poly(urethane urea) oligomer model from the microscopic point of view

The world desperately needs new technologies and solutions for gas capture and separation. To make this possible, molecular modeling is applied here to investigate the structural, thermodynamic, and dynamical properties of a model for the poly(urethane urea) (PUU) oligomer model to selectively capture CO₂ in the presence of CH₄. In this work, we applied a well-known approach to derive atomic partial charges for atoms in a polymer chain based on self-consistent sampling using quantum chemistry and stochastic dynamics. The interactions of the gases with the PUU model were studied in a pure gas based system as well as in a gas mixture. A detailed structure characterization revealed high interaction of CO₂ molecules with the hard segments of the PUU. Therefore, the structural and energy properties explain the reasons for the greater CO₂ sorption than CH₄. We find that the CO₂ sorption is higher than the CH₄ with a selectivity of 7.5 at 298 K for the gas mixture. We characterized the Gibbs dividing surface for each system, and the CO₂ is confined for a long time at the gas-oligomer model interface. The simulated oligomer model showed performance above the 2008 Robeson's upper bound and may be a potential material for CO₂/CH₄ separation. Further computational and experimental studies are needed to evaluate the material.

Using the group contribution method and molecular dynamics to predict the glass transition temperatures and mechanical properties of poly-(p-phenylenediamine-alt-2, 6-diformyl multiphenyl)

This article reports the glass transition temperatures of poly-( p-phenylenediamine-alt-2,6-diformyl multiphenyl) predicted by both the group contribution method and the molecular dynamics simulations. The related modeling method and the degree of polymerization, density, specific volume, radius of volume, radius of rotation, and non-bonding energy terms with temperature are analyzed in depth. The bulk modulus, shear modulus, compressibility, Young’s modulus, and Poisson’s ratio of poly-( p-phenylenediamine-alt-2,6-diformyl multiphenyl) at room temperature are simulated by molecular dynamics. The results show that the simulated glass transition temperatures of poly-( p-phenylenediamine-alt-2,6-diformyl multiphenyl) are greater than 480 K, which indicates that poly-( p-phenylenediamine-alt-2,6-diformyl multiphenyl) can be expected to be used as a high-temperature-resistant material. As the number of rigid benzene rings on the molecular side chain increases, the glass transition temperature decreases, with an average decrease of 10 K for each additional benzene ring. The free volume theory can explain the glass transitions of poly-( p-phenylenediamine-alt-2,6-diformyl multiphenyl). The modulus and density of poly-( p-phenylenediamine-alt-2,6-diformyl multiphenyl) change accordingly with an increase of rigid benzene rings on the side chain, probably due to the fact that the flexibility of the polymers changes accordingly as the number of benzene rings on the side chain increases.

Overscreening, Co-Ion-Dominated Electroosmosis, and Electric Field Strength Mediated Flow Reversal in Polyelectrolyte Brush Functionalized Nanochannels

Controlling the direction and strength of nanofluidic electrohydrodyanmic transport in the presence of an externally applied electric field is extremely important in a number of nanotechnological applications. Here, we employ all-atom molecular dynamics simulations to discover the possibility of changing the direction of electroosmotic (EOS) liquid flows by merely changing the electric field strength in a nanochannel functionalized with polyelectrolyte (PE) brushes. In exploring this, we have uncovered three facets of nanoconfined PE brush behavior and resulting EOS transport. First, we identify the onset of an overscreening effect: such overscreening refers to the presence of more counterions (Na⁺) within the brush layer than needed to neutralize the negative brush charges. Accordingly, as a consequence of the overscreening, in the bulk liquid outside the brush layer, there is a greater number of co-ions (Cl^-) than counterions in the presence of an added salt (NaCl). Second, this specific ion distribution ensures that the overall EOS flow is along the direction of motion of the co-ions. Such co-ion-dictated EOS transport directly contradicts the notion that EOS flow is always dictated by the motion of the counterions. Finally, for large-enough electric fields, the brush height reduces significantly, causing some of the excess overscreening-inducing counterions to squeeze out of the PE brush layer into the brush-free bulk. As a result, the overscreening effect disappears and the number of co-ions and counterions outside the PE brush layer become similar. Despite that there is an EOS transport, this EOS transport, unlike the standard EOS transport that occurs due to the imbalance of the co-ions and counterions, occurs since a larger residence time of the water molecules in the first solvation shell of the counterions (Na⁺) ensures a water transport in the direction of motion of the counterions. The net effect is the reversal of the direction of the EOS transport by merely changing the strength of the electric field.

Machine learning discovery of high-temperature polymers

To formulate a machine learning (ML) model to establish the polymer's structure-property correlation for glass transition temperatureT g , we collect a diverse set of nearly 13,000 real homopolymers from the largest polymer database, PoLyInfo. We train the deep neural network (DNN) model with 6,923 experimentalT g values using Morgan fingerprint representations of chemical structures for these polymers. Interestingly, the trained DNN model can reasonably predict the unknownT g values of polymers with distinct molecular structures, in comparison with molecular dynamics simulations and experimental results. With the validated transferability and generalization ability, the ML model is utilized for high-throughput screening of nearly one million hypothetical polymers. We identify more than 65,000 promising candidates withT g > 200°C, which is 30 times more than existing known high-temperature polymers (∼2,000 from PoLyInfo). The discovery of this large number of promising candidates will be of significant interest in the development and design of high-temperature polymers.

Towards realistic simulations of polymer networks: tuning vulcanisation and mechanical properties

Simulations of coarse-grained network models have long been used to test theoretical predictions about rubber elasticity, while atomistic models are still largely unexplored. Here we devise a novel algorithm for the vulcanisation of united-atom poly(cis-1,4-butadiene), characterize the topology of the resulting networks and test their mechanical properties. We observe clear changes in the network structure when using slower vulcanisation, contrary to the traditional view that cross-linking simply freezes the melt configuration. Non-ideality of our networks reverberates on the distribution of strand length and on the strands deformation, which is highly non-affine, especially for short strands. Nevertheless, we do recover some of the trends observed on ideal bead-and-spring networks and controlled laboratory experiments, such as the linear relationships linking the degree of cross-linking and the density. We also compare different deformation methods and find step-equilibrium protocols to be more reliable. Regardless of the adopted method, it is advisable to precede the deformation by a pre-stretching cycle in order to release internal stresses accumulated during the vulcanisation.

Coarse-graining of polyisoprene melts using inverse Monte Carlo and local density potentials

Bottom-up coarse-graining of polymers is commonly performed by matching structural order parameters such as distribution of bond lengths, bending and dihedral angles, and pair distribution functions. In this study, we introduce the distribution of nearest-neighbors as an additional order parameter in the concept of local density potentials. We describe how the inverse-Monte Carlo method provides a framework for forcefield development that is capable of overcoming challenges associated with the parameterization of interaction terms in polymer systems. The technique is applied on polyisoprene melts as a prototype system. We demonstrate that while different forcefields can be developed that perform equally in terms of matching target distributions, the inclusion of nearest-neighbors provides a straightforward route to match both thermodynamic and conformational properties. We find that several temperature state points can also be addressed, provided that the forcefield is refined accordingly. Finally, we examine both the single-particle and the collective dynamics of the coarse-grain models, demonstrating that all forcefields present a similar acceleration relative to the atomistic systems.

Molecular dynamics study of natural rubber-fullerene composites: connecting microscopic properties to macroscopic behavior

Macroscopic and microscopic properties of fullerene (C₆₀)-cis-polyisoprene (cis-PI) composites at varying fullerene concentrations were investigated using atomistic molecular dynamics (MD) simulations over microsecond time scales. Results show that the introduction of fullerenes into a polymer matrix increases density, bulk modulus and heat capacity while thermal expansivity decreases. The presence of fullerenes slowed the diffusion of both C₆₀ and cis-PI. Moreover, increasing fullerene concentration results in ordering of the cis-PI chains at the cis-PI-fullerene interfaces and shrinking of bulk PI regions. Free energy calculations of fullerene dimerization suggest that fullerenes disperse at low and aggregate at high fullerene concentrations. Our multi-scaled analysis approach demonstrates the role of 'ordered' regions adjacent to the interface between cis-PI and fullerene in controlling the level of order and mobility of the cis-PI chains. The relationship between the microscopic behavior and the changes in mechanical and thermal properties are discussed. Our study is beneficial for further studies and development of advanced rubber technology for novel, cost-effective, material with very high stiffness and thermal endurance with optimizing conditions of filler contents.

Hydration of diblock copolymer micelles: Effects of hydrophobicity and co-solvent

Diblock polymer micelles dispersed in an aqueous environment are being actively investigated for various applications, but there is only a qualitative understanding of the effect of the chemical structure on the micelle hydration and water dynamics as these properties are difficult to assess experimentally. Using all-atom molecular dynamics simulations, we investigate aqueous solutions of three comparable in size diblock copolymer micelles with core-forming blocks of different hydrophobicity: polybutadiene (PB), polycaprolactone (PCL), and polytetrahydrofuran (pTHF) with the same hydrophilic block, polyethylene oxide (PEO). We found that core-block hydrophobicity and ability to form hydrogen bonds with water strongly affect the water dynamics near the core: water molecules spend considerably less time in contact with the PB block than with PCL and pTHF blocks. We obtained polymer and solvent volume fraction profiles and determined that the interfacial width systematically increases with a decrease of core block hydrophobicity with water penetration into the core being negligible for PB-PEO and PCL-PEO micelles, while for pTHF-PEO micelles the interface is more diffuse and there is a noticeable penetration of water (17% by volume). For PCL-PEO micelles, which are commonly used in biomedical applications, we also investigated tetrahydrofuran (THF) penetration into the micelles from mixed THF:water solution at early stages of micelle dissolution. We found an inhomogeneous solvent distribution with a maximum of THF volume fraction in the interfacial core-corona region and partial exclusion from the PEO corona, which slows down micelle dissolution. These results can have important implications for micelle stability and use in biomedical applications.

Molecular View of CO2 Capture by Polyethylenimine: Role of Structural and Dynamical Heterogeneity

The molecular thermodynamics and kinetics of CO₂ sorption in Polyethylenimine (PEI) melt have been investigated systematically using GCMC and MD simulations. We elucidate presence of significant structural and dynamic heterogeneity associated with the overall absorption process. CO₂ adsorption in a PEI membrane shows a distinct two-stage process of a rapid CO₂ adsorption at the interfaces (hundreds of picoseconds) followed by a significantly slower diffusion limited release toward the interior bulk regions of PEI melt (hundreds of nanoseconds to microseconds). The spatial heterogeneity of local structural features of the PEI chains lead to significantly heterogeneous absorption characterized by clustering and trapping of CO₂ molecules that then lead to subdiffusive motion of CO₂. In the complex interplay of interaction and entropy, the latter emerges out to be the major determining factor with significantly higher solubility of CO₂ near the interfaces despite having lower density of binding amine groups. Regions having higher free-volume (entropically favorable) viz. interfaces, pores and loops demonstrate higher CO₂ capture ability. Various local structural features of PEI conformations, for example, inter- and intrachain loops, pores of different radii, and di- or tricoordinated pores are explored for their effects on the varying CO₂ adsorption abilities.

Molecular dynamics simulations of oligoester brushes: the origin of unusual conformations

We present results from all-atom molecular dynamics simulations for the structural properties of oligomeric lactic acid chains (OLA) grafted to the surface of cellulose nanocrystals (CNCs) and immersed in the melt of polylactic acid (PLA). Earlier, we have found that the distribution of free ends of OLA molecules is bimodal [Glova et al., Polym. Int., 2016, 65(8), 892]. The results cannot be explained within the standard picture of uncharged polymer brushes exposed to the melt of a chemically identical polymer. Although the oligomeric brushes of the OLA chains are uncharged, they have partial polarization charges producing a non-zero dipole moment of the monomeric chain unit. We study the influence of partial charges on the structure of the layer of OLA chains grafted to the CNC surface. A detailed analysis of the conformations of the grafted chains shows that interaction of partial charges in the models causes bending of the OLA molecules toward the cellulose surface, forming a hairpin structure. The observed separation of the grafted chains into two populations increases with grafting density. We demonstrate that hydrogen bonds can be formed between the free ends of the grafted chains and the CNC surface, but they do not affect the brush structure significantly. Thus, dipole-dipole interactions turn out to be the key factor governing the unusual conformations of grafts.

Filler reinforcement in cross-linked elastomer nanocomposites: insights from fully atomistic molecular dynamics simulation

Using a fully atomistic model, we perform large-scale molecular dynamics simulations of sulfur-cured polybutadiene (PB) and nanosilica-filled PB composites. A well-integrated network without sol fraction is built dynamically by cross-linking the coarse-grained precursor chains in the presence of embedded silica nanoparticles. Initial configurations for subsequent atomistic simulations are obtained by reverse mapping of the well-equilibrated coarse-grained systems. Based on the concept of "maximally inflated knot" introduced by Grosberg et al., we show that the networks simulated in this study behave as mechanically isotropic systems. Analysis of the network topology in terms of graph theory reveals that mechanically inactive tree-like structures are the dominant structural components of the weakly cross-linked elastomer, while cycles are mainly responsible for the transmission of mechanical forces through the network. We demonstrate that quantities such as the system density, thermal expansion coefficient, glass transition temperature and initial Young's modulus can be predicted in qualitative and sometimes even in quantitative agreement with experiments. The nano-filled system demonstrates a notable increase in the glass transition temperature and an approximately two-fold increase in the nearly equilibrium value of elastic modulus relative to the unfilled elastomer even at relatively small amounts of filler particles. We also examine the structural rearrangement of the nanocomposite subjected to tensile deformation. Under high strain-rate loading, the formation of structural defects (microcavities) within the polymer bulk is observed. The nucleation and growth of cavities in the post-yielding strain hardening regime mainly take place at the elastomer/nanoparticle interfaces. As a result, the cavities are concentrated just near the embedded nanoparticles. Therefore, while the silica nanofiller increases the elastic modulus of the elastomer, it also creates a more defective structure of higher energy in comparison with the unfilled network.