Glass transition in 1,4-polybutadiene: Mode-coupling theory analysis of molecular dynamics simulations using a chemically realistic model

Structure and dynamics of a 1,4-polybutadiene melt in an alumina nanopore: A molecular dynamics simulation

We present results of Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene (PBD)confined in a cylindrical alumina nanopore of diameter 10 nm. The simulations are done at three different temperaturesabove the glass transition temperature T g . We investigate the density layering across the nanopore as well as theorientational ordering in the polymer melt, brought about by the confinement, on both the segmental and chain scales.For the chain scale ordering, the magnitude and orientation of the axes of the gyration tensor ellipsoid of single chainsare studied and are found to prefer to align parallel to the pore axis. Even though double bonds near the wall arepreferentially oriented along the pore walls, studying the nematic order parameter indicates that there is no nematicordering at the melt-wall interface. As for the dynamics in the melt, we focus here on the mean-square-displacement ofthe monomers for several layers across the nanopore as well as the movement of the chain center of mass which bothdisplay a slowing down of the dynamics in the layer at the wall. We also show the strong adsorption of the monomersto the pore wall at lower temperatures.

Role of torsional potential in chain conformation, thermodynamics, and glass formation of simulated polybutadiene melts

For polymer chains, the torsional potential is an important intramolecular energy influencing chain flexibility and segmental dynamics. Through molecular dynamics simulations of an atomistic model for melts of cis-trans-1,4-polybutadiene (PBD), we explore the effect of the torsions on conformational properties (bond vector correlations and mean-square internal distances), fundamental thermodynamic quantities (density, compressibility, internal energy, and specific heat), and glass transition temperature T_g. This is achieved by systematically reducing the strength of the torsional potential, starting from the chemically realistic chain (CRC) model with the full potential toward the freely rotating chain (FRC) model without the torsional potential. For the equilibrium liquid, we find that the effect of the torsions on polymer conformations is very weak. Still weaker is the influence on the monomer density ρ and isothermal compressibility κ_T of the polymer liquid, both of which can be considered as independent of the torsional potential. We show that a van der Waals-like model proposed by Long and Lequeux [Eur. Phys. J. E 4, 371 (2001)] allows us to describe very well the temperature (T) dependence of ρ and κ_T. We also find that our data obey the linear relation between 1/k_BTρκ_T and 1/T (with the Boltzmann constant k_B) that has recently been predicted and verified on the experiment by Mirigian and Schweizer [J. Chem. Phys. 140, 194507 (2014)]. For the equilibrium liquid, simulations result in a specific heat, at constant pressure and at constant volume, which increases on cooling. This T dependence is opposite to the one found experimentally for many polymer liquids, including PBD. We suggest that this difference between simulation and experiment may be attributed to quantum effects due to hydrogen atoms and backbone vibrations, which, by construction, are not included in the classical united-atom model employed here. Finally, we also determine T_g from the density-temperature curve monitored in a finite-rate cooling process. While the influence of the torsional potential on ρ(T) is vanishingly small in the equilibrium liquid, the effect of the torsions on T_g is large. We find that T_g decreases by about 150 K when going from the CRC to the FRC model.

Diffusivity and glass transition of polymer chains in polymer nanocomposites

The diffusivity and glass transition of polymer chains in polymer nanocomposites are studied by using dynamic Monte Carlo simulation. Nanoparticles are modeled as immobile and distributed in a cubic lattice in the system. The diffusion coefficient D of polymer chains is reduced, while the glass transition temperature Tg is increased by nanoparticles. Our results show that the effect of nanoparticles can be summarized as D = D0[1 - exp(-α·ID/2Rg)] and Tg = Tg,0[1 - exp(-α·ID/2Rg)]-1, with D0 and Tg,0 being the diffusion coefficient and the glass transition temperature in the absence of nanoparticles, Rg the radius of gyration of polymer chains, and ID the surface spacing between nearest-neighbor nanoparticles. The parameter α that governs the dynamics of polymer chains decreases with increasing nanoparticles' size or decreasing the temperature. Our results also show that smaller nanoparticles exert a stronger influence on the polymer dynamics at the same concentration of nanoparticles, whereas larger nanoparticles show a stronger effect at the same ID.

Chain relaxation in thin polymer films: turning a dielectric type-B polymer into a type-A' one

A molecular dynamics simulation study of chain relaxation in a thin polymer film is presented, studying the dielectric response of a random copolymer of cis and trans 1,4-polybutadiene, a type B polymer without net chain dipole moment, confined between graphite walls. We stress the orientational effect of the attractive walls, inducing polarization in the vicinity of the walls, while the center of the film stays bulk-like. This polarization leads to a net dipole moment of the adsorbed chains, which is perpendicular to their end-to-end vector, which we termed as type A' behavior. In this situation, the dipole moment relaxes only upon desorption of the chains from the wall, a dynamic process which occurs on timescales much longer than the bulk relaxation time of the polymer.

Glassy dynamics of model colloidal polymers: The effect of "monomer" size

In recent years, attempts have been made to assemble colloidal particles into chains, which are termed "colloidal polymers." An apparent difference between molecular and colloidal polymers is the "monomer" size. Here, we propose a model to represent the variation from molecular polymer to colloidal polymer and study the quantitative differences in their glassy dynamics. For chains, two incompatible local length scales, i.e., monomer size and bond length, are manifested in the radial distribution function and intramolecular correlation function. The mean square displacement of monomers exhibits Rouse-like sub-diffusion at intermediate time/length scale and the corresponding exponent depends on the volume fraction and the monomer size. We find that the threshold volume fraction at which the caging regime emerges can be used as a rescaling unit so that the data of localization length versus volume fraction for different monomer sizes can gather close to an exponential curve. The increase of monomer size effectively increases the hardness of monomers and thus makes the colloidal polymers vitrify at lower volume fraction. Static and dynamic equivalences between colloidal polymers of different monomer sizes have been discussed. In the case of having the same peak time of the non-Gaussian parameter, the motion of monomers of larger size is much less non-Gaussian. The mode-coupling critical exponents for colloidal polymers are in agreement with that of flexible bead-spring chains.

Are polymers standard glass-forming systems? The role of intramolecular barriers on the glass-transition phenomena of glass-forming polymers

Traditionally, polymer melts have been considered archetypal glass-formers. This has been mainly due to the fact that these systems can easily be obtained as glasses by cooling from the melt, even at low cooling rates. However, the macromolecules, i.e. the structural units of polymer systems in general, are rather different from the standard molecules. They are long objects ('chains') made by repetition of a given chemical motif (monomer) and have intra-macromolecular barriers that limit their flexibility. The influence of these properties on, for instance, the glass-transition temperature of polymers, is a topic that has been widely studied by the polymer community almost from the early times of polymer science. However, in the framework of the glass-community, the relevant influence of intra-macromolecular barriers and chain connectivity on glass-transition phenomena of polymers has started to be recognized only recently. The aim of this review is to give an overview and to critically revise the results reported on this topic over the last years. From these results, it seems to be evident that there are two different mechanisms involved in the dynamic arrest in glass-forming polymers: (i) the intermolecular packing effects, which dominate the dynamic arrest of low molecular weight glass-forming systems; and (ii) the effect of intra-macromolecular barriers combined with chain connectivity. It has also been shown that the mode coupling theory (MCT) is a suitable theoretical framework to discuss these questions. The values found for polymers for the central MCT parameter--the so-called λ-exponent--are of the order of 0.9, clearly higher than the standard values (λ ≈ 0.7) found in systems where the dynamic arrest is mainly driven by packing effects ('standard' glass-formers). Within the MCT, this is a signature of the presence of two competing mechanisms of dynamic arrest, as it has been observed in short-ranged attractive colloids or two component mixtures with dynamic asymmetry. Moreover, recent MD-simulations of a 'bead-spring' polymer model, but including intra-macromolecular potential of different strengths, confirm that the high λ-values found in polymers are due to the effect of intra-macromolecular barriers. Although there are still open questions, these results allow to conclude that there is a fundamental difference between the nature of the glass transition in polymers and in simple (standard) glass-formers.

Simulated glass-forming polymer melts: dynamic scattering functions, chain length effects, and mode-coupling theory analysis

We present molecular-dynamics simulations for a fully flexible model of polymer melts with different chain length N ranging from short oligomers (N = 4) to values near the entanglement length (N = 64). For these systems we explore the structural relaxation of the supercooled melt near the critical temperature T c of mode-coupling theory (MCT). Coherent and incoherent scattering functions are analyzed in terms of the idealized MCT. For temperatures T > T c we provide evidence for the space-time factorization property of the β relaxation and for the time-temperature superposition principle (TTSP) of the α relaxation, and we also discuss deviations from these predictions for T ≈ T c. For T larger than the smallest temperature where the TTSP holds we perform a quantitative analysis of the dynamics with the asymptotic MCT predictions for the late β regime. Within MCT a key quantity, in addition to T c, is the exponent parameter λ. For the fully flexible polymer models studied we find that λ is independent of N and has a value (λ = 0.735 ) typical of simple glass-forming liquids. On the other hand, the critical temperature increases with chain length toward an asymptotic value T c (∞) . This increase can be described by T c (∞) - T c(N) ∼ 1/N and may be interpreted in terms of the N dependence of the monomer density ρ, if we assume that the MCT glass transition is ruled by a soft-sphere-like constant coupling parameter Γ c = ρ c T c (-1/4), where ρ c is the monomer density at T c. In addition, we also estimate T c from a Hansen-Verlet-like criterion and MCT calculations based on structural input from the simulation. For our polymer model both the Hansen-Verlet criterion and the MCT calculations suggest T c to decrease with increasing chain length, in contrast to the direct analysis of the simulation data.

Molecular dynamics simulation of amorphous hydroxypropyl-methylcellulose acetate succinate (HPMCAS): polymer model development, water distribution, and plasticization

Molecular models for HPMCAS polymer have been developed for molecular dynamics (MD) simulation that attempt to mimic the complex substitution patterns in HPMCAS observed experimentally. These molecular models were utilized to create amorphous HPMCAS solids by cooling of the polymeric melts at different water contents to explore the influence of water on molecular mobility, which plays a critical role in stability and drug release from HPMCAS-based solid matrices. The densities found for the simulated amorphous HPMCAS were 1.295, 1.287, and 1.276 g/cm(3) at 0.7, 5.7, and 13.2% w/w water, indicating swelling of the polymer with increasing water content. These densities compare favorably with the experimental density of 1.285 g/cm(3) for commercial HPMCAS-(AQOAT AS-MF) supporting the present HPMCAS models as a realistic representation of amorphous HPMCAS solids. Water molecules were observed to be mostly isolated from each other at a low water content (0.7% w/w), while clusters or strands of water were pervasive and broadly distributed in size at 13.2% w/w water. The average number of first-shell water molecules (n(w)) increased from 0.17 to 3.5, though the latter is still far below that (8.9) expected for the onset of a separate water phase. Increasing water content from 0.7 to 13.2% w/w was found to reduce the T(g) by ~81 K, similar to experimental observations. Plasticization with increasing water content resulted in increasing polymer mobility and water diffusivity. From 0.7 to 13.2% w/w water, the apparent water diffusivity increased from 1.1 × 10(-9) to 7.0 × 10(-8) cm(2)/s, though non-Einsteinian behavior persisted at all water contents explored. This and the water trajectories in the polymers suggest that water diffusion at 0.7% w/w water follows a "hopping" mechanism. At a higher water content (13.2% w/w) water diffusion follows dual diffusive processes: (1) fast water motions within water clusters; and (2) slower diffusion through the more rigid polymer matrix.

Applicability of mode-coupling theory to polyisobutylene: a molecular dynamics simulation study

The applicability of Mode Coupling Theory (MCT) to the glass-forming polymer polyisobutylene (PIB) has been explored by using fully atomistic molecular dynamics simulations. MCT predictions for the so-called asymptotic regime have been successfully tested on the dynamic structure factor and the self-correlation function of PIB main-chain carbons calculated from the simulated cell. The factorization theorem and the time-temperature superposition principle are satisfied. A consistent fitting procedure of the simulation data to the MCT asymptotic power-laws predicted for the α-relaxation regime has delivered the dynamic exponents of the theory-in particular, the exponent parameter λ-the critical non-ergodicity parameters, and the critical temperature T(c). The obtained values of λ and T(c) agree, within the uncertainties involved in both studies, with those deduced from depolarized light scattering experiments [A. Kisliuk et al., J. Polym. Sci. Part B: Polym. Phys. 38, 2785 (2000)]. Both, λ and T(c)/T(g) values found for PIB are unusually large with respect to those commonly obtained in low molecular weight systems. Moreover, the high T(c)/T(g) value is compatible with a certain correlation of this parameter with the fragility in Angell's classification. Conversely, the value of λ is close to that reported for real polymers, simulated "realistic" polymers and simple polymer models with intramolecular barriers. In the framework of the MCT, such finding should be the signature of two different mechanisms for the glass-transition in real polymers: intermolecular packing and intramolecular barriers combined with chain connectivity.

Molecular dynamics simulation of amorphous indomethacin

Molecular dynamics (MD) simulations have been conducted using an assembly consisting of 105 indomethacin (IMC) molecules and 12 water molecules to investigate the underlying dynamic (e.g., rotational and translational diffusivities and conformation relaxation rates) and structural properties (e.g., conformation, hydrogen-bonding distributions, and interactions of water with IMC) of amorphous IMC. These properties may be important in predicting physical stability of this metastable material. The IMC model was constructed using X-ray diffraction data with the force-field parameters mostly assigned by analogy with similar groups in Amber-ff03 and atomic charges calculated with the B3LYP/ccpVTZ30, IEFPCM, and RESP models. The assemblies were initially equilibrated in their molten state and cooled through the glass transition temperature to form amorphous solids. Constant temperature dynamic runs were then carried out above and below the T(g) (i.e., at 600 K (10 ns), 400 K (350 ns), and 298 K (240 ns)). The density (1.312 ± 0.003 g/cm(3)) of the simulated amorphous solid at 298 K was close to the experimental value (1.32 g/cm(3)) while the estimated T(g) (384 K) was ~64 degrees higher than the experimental value (320 K) due to the faster cooling rate. Due to the hindered rotation of its amide bond, IMC can exist in different diastereomeric states. Different IMC conformations were sufficiently sampled in the IMC melt or vapor, but transitions occurred rarely in the glass. The hydrogen-bonding patterns in amorphous IMC are more complex in the amorphous state than in the crystalline polymorphs. Carboxylic dimers that are dominant in α- and γ-crystals were found to occur at a much lower probability in the simulated IMC glasses while hydrogen-bonded IMC chains were more easily identified patterns in the simulated amorphous solids. To determine molecular diffusivity, a novel analytical method is proposed to deal with the non-Einsteinian behavior, in which the temporal evolution of the apparent diffusivity D(t) is described by a relaxation model such as the KWW function and extrapolated to infinite time. The diffusion coefficient found for water diffusing in amorphous indomethacin at 298 K (2.7 × 10(-9) cm(2)/s) compares favorably to results obtained in experimental IMC glasses (0.9-2.0 × 10(-9) cm(2)/s) and is mechanistically associated with β-relaxation processes that are dominant in sub-T(g) glasses.

Slow process in confined polymer melts: layer exchange dynamics at a polymer solid interface

Employing Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene between graphite walls we show that the mass exchange between layers close to the walls is a slow process already in the melt state. For the glass transition of confined polymers this process competes with the slowing down due to packing effects and intramolecular rotation barriers.

Dynamic arrest in polymer melts: competition between packing and intramolecular barriers

We present molecular dynamics simulations of a simple model for polymer melts with intramolecular barriers. We investigate structural relaxation as a function of the barrier strength. Dynamic correlators can be consistently analyzed within the framework of the mode coupling theory of the glass transition. Control parameters are tuned in order to induce a competition between general packing effects and polymer-specific intramolecular barriers as mechanisms for dynamic arrest. This competition yields unusually large values of the so-called mode coupling theory exponent parameter and rationalizes qualitatively different observations for simple bead-spring and realistic polymers. The systematic study of the effect of intramolecular barriers presented here also establishes a fundamental difference between the nature of the glass transition in polymers and in simple glass formers.

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.

Structural and conformational dynamics of supercooled polymer melts: insights from first-principles theory and simulations

We report on quantitative comparisons between simulation results of a bead-spring model and mode-coupling theory calculations for the structural and conformational dynamics of a supercooled, unentangled polymer melt. We find semiquantitative agreement between simulation and theory, except for processes that occur on intermediate length scales between the compressibility plateau and the amorphous halo of the static structure factor. Our results suggest that the onset of slow relaxation in a glass-forming melt can be described in terms of monomer caging supplemented by chain connectivity. Furthermore, a unified atomistic description of glassy arrest and of conformational fluctuations that (asymptotically) follow the Rouse model emerges from our theory.