Partial Structure Factors in 1,4-Polybutadiene. A Combined Neutron Scattering and Molecular Dynamics Simulations Study

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.

Strategies for the Development of Conjugated Polymer Molecular Dynamics Force Fields Validated with Neutron and X-ray Scattering

Conjugated polymers (CPs) enable a wide range of lightweight, lower cost, and flexible organic electronic devices, but a thorough understanding of relationships between molecular structure and dynamics and electronic performance is critical for improved device efficiencies and for new technologies. Molecular dynamics (MD) simulations offer in silico insight into this relationship, but their accuracy relies on the approach used to develop the model's parameters or force field (FF). In this Perspective, we first review current FFs for CPs and find that most of the models implement an arduous reparameterization of inter-ring torsion potentials and partial charges of classical FFs. However, there are few FFs outside of simple CP molecules, e.g., polythiophenes, that have been developed over the last two decades. There is also limited reparameterization of other parameters, such as nonbonded Lennard-Jones interactions, which we find to be directly influenced by conjugation in these materials. We further provide a discussion on experimental validation of MD FFs, with emphasis on neutron and X-ray scattering. We define multiple ways in which various scattering methods can be directly compared to results of MD simulations, providing a powerful experimental validation metric of local structure and dynamics at relevant length and time scales to charge transport mechanisms in CPs. Finally, we offer a perspective on the use of neutron scattering with machine learning to enable high-throughput parametrization of accurate and experimentally validated CP FFs enabled not only by the ongoing advancements in computational chemistry, data science, and high-performance computing but also using oligomers as proxies for longer polymer chains during FF development.

Insight into the Structure and Dynamics of Polymers by Neutron Scattering Combined with Atomistic Molecular Dynamics Simulations

Combining neutron scattering and fully atomistic molecular dynamics simulations allows unraveling structural and dynamical features of polymer melts at different length scales, mainly in the intermolecular and monomeric range. Here we present the methodology developed by us and the results of its application during the last years in a variety of polymers. This methodology is based on two pillars: (i) both techniques cover approximately the same length and time scales and (ii) the classical van Hove formalism allows easily calculating the magnitudes measured by neutron scattering from the simulated atomic trajectories. By direct comparison with experimental results, the simulated cell is validated. Thereafter, the information of the simulations can be exploited, calculating magnitudes that are experimentally inaccessible or extending the parameters range beyond the experimental capabilities. We show how detailed microscopic insight on structural features and dynamical processes of various kinds has been gained in polymeric systems with different degrees of complexity, and how intriguing questions as the collective behavior at intermediate length scales have been faced.

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.

Atomic motions in the alphabeta-merging region of 1,4-polybutadiene: a molecular dynamics simulation study

We present fully atomistic molecular dynamics simulations on 1,4-polybutadiene in a wide temperature range from 200 to 280 K, i.e., in the region where the alpha- and beta-relaxations merge and above. A big computational effort has been performed-especially for the lowest temperatures investigated-to extend the simulation runs to very long times (up to 1 mus for 200 K). The simulated sample has been carefully validated by using previous neutron scattering data on the real sample with similar microstructure. Inspecting the trajectories of the different hydrogens in real space, we have observed a heterogeneous dynamical behavior (each kind of hydrogen moves in a different way) with signatures of combined hopping and diffusive motions in the whole range investigated. The application of a previously proposed model [Colmenero et al., Europhys. Lett. 71, 262 (2005)] is successful and a characterization of the local motions and diffusion is possible. The comparison of our results to those reported in the literature provides a consistent scenario for polybutadiene dynamics and puts into a context the different experimental observations. We also discuss the impact of the hopping processes on the observation and interpretation of experimentally accessible magnitudes and the origin of the deviations from Gaussian behavior in this system.

Positron-annihilation-lifetime response and broadband dielectric relaxation spectroscopy: diethyl phthalate

We report the results of a combined phenomenological analysis of the data from positron-annihilation-lifetime spectroscopy (PALS) and the relaxation data from broadband dielectric spectroscopy (BDS) on diethyl phthalate (DEP). The ortho-positronium (o-Ps) lifetime, tau3, as a function of temperature over a temperature range from 67 K up to 300 K is compared with the spectral features and the relaxation parameters of the BDS spectra decomposed into the primary alpha and the secondary beta processes in the temperature range from 140 K up to 380 K by using the Williams-Watts scheme. Phenomenological model-free analysis of the tau3-T plot provides the three characteristic PALS temperatures, where the two most pronounced ones at TgPALS=185 K and Tb2=245 K=1.32TgPALS are related to the glass-liquid transition and the onset of a quasiplateau region, respectively. In the case of a weaker bend effect at Tb1=210 K=1.14TgPALS, a number of new coincidences with changes in the dielectric beta process have been found. They concern the changes in width parameter of the distribution function for the beta relaxation time and the activation energy of the betaeff process, a crossover from the Arrhenius to the non-Arrhenius type of temperature dependence as well as with the onset of a short-time tail of the beta relaxation time distribution and finally, with changes in the relaxation strength of the alpha and betaeff processes. All these findings indicate a close connection of the o -Ps annihilation parameters and relaxation characteristics of BDS response for the DEP matrix.

"Self-concentration" effects on the dynamics of a polychlorinated biphenyl diluted in 1,4-polybutadiene

The mobility of isolated polychlorinated biphenyl (PCB54) in 1,4-polybutadiene (PB) has been investigated by means of broadband dielectric spectroscopy. The aim was to provide new insights about the effect of the environment on the dynamics of PCB54. The authors' results indicate that PCB54 structural dynamics is neither independent of the PB matrix nor slaved to the matrix itself. The authors interpret these results as a consequence of the limited size of cooperatively rearranging regions (CRRs) involved in PCB54 structural relaxation possessing an effective concentration different from the macroscopic one. This implies a non-negligible influence of "self-concentration," already proven for the component segmental dynamics in polymer blends, also in the relaxation of binary mixtures involving low molecular weight glass formers. This allowed the evaluation of the size of CRR, which was about 1 nm for PCB54 in PB. This means that the cooperativity extends over the first shell around PCB54 molecules.