Spatial correlations in polycarbonates: Neutron scattering and simulation

Molecular basis of secondary relaxation in stiff-chain glassy polymers

Recent progress in establishing local order in polycarbonate-like glasses using rotational echo double resonance and centerband-only detection of exchange solid-state nuclear magnetic resonance (NMR) has stimulated a renewed attempt to connect molecular motion within glassy polymers and the mechanical properties of the glass. We have in fact established a correlation between molecular motion characterized by NMR and the mechanical secondary relaxation (tan δ) for nine polycarbonate-like glasses. All of the NMR and mechanical data are for T ≪ Tg. The resulting structural insights suggest that the chains of these polymers are simultaneously both Flory random coils and Vol'kenstein bundles. The cooperative motions of groups of bundles can be described qualitatively by a variety of constrained-kinetics models of the glass. All of the models share a common trait for large-amplitude motion: an exponential increase in the time required for an inter-bundle dilation event with a linear increase in bundle group size. This dependence and a locally ordered Vol'kenstein bundle lead to an understanding of the surprising 60° (K) shift of tan δ to higher temperature for ring-fluoro-polycarbonate relative to that of polycarbonate by the apparently minor substitution of a fluorine for a hydrogen on every fourth ring.

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.

The Perfect Glass Paradigm: Disordered Hyperuniform Glasses Down to Absolute Zero

Rapid cooling of liquids below a certain temperature range can result in a transition to glassy states. The traditional understanding of glasses includes their thermodynamic metastability with respect to crystals. However, here we present specific examples of interactions that eliminate the possibilities of crystalline and quasicrystalline phases, while creating mechanically stable amorphous glasses down to absolute zero temperature. We show that this can be accomplished by introducing a new ideal state of matter called a "perfect glass". A perfect glass represents a soft-interaction analog of the maximally random jammed (MRJ) packings of hard particles. These latter states can be regarded as the epitome of a glass since they are out of equilibrium, maximally disordered, hyperuniform, mechanically rigid with infinite bulk and shear moduli, and can never crystallize due to configuration-space trapping. Our model perfect glass utilizes two-, three-, and four-body soft interactions while simultaneously retaining the salient attributes of the MRJ state. These models constitute a theoretical proof of concept for perfect glasses and broaden our fundamental understanding of glass physics. A novel feature of equilibrium systems of identical particles interacting with the perfect-glass potential at positive temperature is that they have a non-relativistic speed of sound that is infinite.

Automated parametrization of the coarse-grained Martini force field for small organic molecules

The systematic exploration of chemical compound space holds many promises toward structure-function relationships and material design. In the context of computer simulations, progress is hampered by both the sheer number of compounds and the efforts associated with parametrizing a force field for every new molecule. A coarse-grained (CG) representation provides not only a reduced phase space but also a smaller number of compounds, due to the redundancy of CG representations mapping to the same structure. Though many CG models require the explicit force-field parametrization of a molecule with all others, others assume transferability by means of mixing rules, such as the Martini force field. To alleviate the burden associated with tedious parametrizations for each new compound, the present work aims at automating the mapping and parametrization of common small organic molecules for Martini. We test the method by analyzing the water/octanol partitioning of more than 650 neutral molecules, the hydration free energy of 354 others, and the free energies of hydration and solvation in octanol of another 69 compounds. Last, we compare with all-atom simulations the thermodynamics of insertion of four individual solute molecules in a phospholipid membrane. The protocol demonstrates the feasibility of an automated parametrization scheme for Martini and provides prospects for high-throughput simulation methodologies.

Effect of chemical substituents on the structure of glassy diphenyl polycarbonates

Polycarbonates offer a wide variety of physical property behavior that is difficult to predict due to complexities at the molecular scale. Here, the physical structure of amorphous glassy polycarbonates having aliphatic and cycloaliphatic chemical groups is explored through atomistic simulations. The influence of chemical structure on solubility parameter, torsion distributions, radial distribution function, scattering structure factor, orientation distributions of phenylene rings and carbonate groups, and free volume distributions, leading to interchain packing effects, are shown. The effect of the cyclohexyl ring at the isopropylidene carbon as compared to the effect of the methyl groups positioned on the phenylene rings results in a larger reduction in the solubility parameter (δ). The interchain distance estimated for polycarbonates in this work is in the range of 5-5.8 Å. The o-methyl groups on the phenylene rings, as compared to a cyclohexyl ring, lead to higher interchain distances. The highest interchain distance is observed with a trimethylcyclohexylidene group at the isopropylidene carbon. Atomistic simulations reveal two different types of packing arrangement of nearest-neighbor chains in the glassy state, one type of which agrees with the NMR experimental data. The fundamental insights provided here can be utilized for design of chemical structures for tailored macroscopic properties.

Chain packing in polycarbonate glasses

Chain packing in homogeneous blends of carbonate (13)C-labeled bisphenol A polycarbonate with either (i) CF(3)-labeled bisphenol A polycarbonate or (ii) ring-F-labeled bisphenol A polycarbonate has been characterized using (13)C{(19)F} rotational-echo double-resonance (REDOR) nuclear magnetic resonance. In both blends, the (13)C observed spin was at high concentration, and the (19)F dephasing or probe spin was at low concentration. In this situation, an analysis in terms of a distribution of isolated heteronuclear pairs of spins is valid. Nearest-neighbor separation of (13)C and (19)F labels was determined by accurately mapping the initial dipolar evolution using a shifted-pulse version of REDOR. Based on the results of this experiment, the average distance from a ring-fluorine to the nearest (13)C=O is more than 1.2 A greater than the corresponding CF(3)-(13)C=O distance. Next-nearest and more-distant-neighbor separations of labels were measured in a 416-rotor-cycle constant-time version of REDOR for both blends. Statistically significant local order was established for the nearest-neighbor labels in the methyl-labeled blend. These interchain packing results are in qualitative agreement with predictions based on coarse-grained simulations of a specially adapted model for bisphenol A polycarbonate. The model itself has been previously used to determine static and dynamic properties of polycarbonate with results in good agreement with those from rheological and neutron scattering experiments.

Effect of stretching on the sub-Tg phenylene-ring dynamics of polycarbonate by neutron scattering

We have investigated the effect of cold drawing on the motion of phenylene rings in bisphenol-A polycarbonate by neutron scattering. The intensity scattered by isotropic and stretched polycarbonate is different, the quasielastic broadening being larger for the isotropic sample. This difference can be well accounted for by considering that preferentially oriented rings in the stretched polymer have their motion sterically hindered. The extent of the effect of stretching on the phenylene-ring motion obtained from this neutron-scattering investigation is in good agreement with that obtained when studying the effect of cold drawing on the gamma relaxation by dielectric spectroscopy.

Coarse-grained force field for simulating polymer-tethered silsesquioxane self-assembly in solution

A coarse-grained model has been developed for simulating the self-assembly of nonyl-tethered polyhedral oligomeric silsesquioxane (POSS) nanoparticles in solution. A mapping scheme for groups of atoms in the atomistic molecule onto beads in the coarse-grained model was established. The coarse-grained force field consists of solvent-mediated effective interaction potentials that were derived via a structural-based coarse-graining numerical iteration scheme. The force field was obtained from initial guesses that were refined through two different iteration algorithms. The coarse-graining scheme was validated by comparing the aggregation of POSS molecules observed in simulations of the coarse-grained model to that observed in all-atom simulations containing explicit solvent. At 300 K the effective coarse-grained potentials obtained from different initial guesses are comparable to each other. At 400 K the differences between the force fields obtained from different initial guesses, although small, are noticeable. The use of a different iteration algorithm employing identical initial guesses resulted in the same overall effective potentials for bare cube corner bead sites. In both the coarse-grained and all-atom simulations, small aggregates of POSS molecules were observed with similar local packings of the silsesquioxane cages and tether conformations. The coarse-grained model afforded a savings in computing time of roughly two orders of magnitude. Further comparisons were made between the coarse-grained monotethered POSS model developed here and a minimal model developed in earlier work. The results suggest that the interactions between POSS cages are long ranged and are captured by the coarse-grained model developed here. The minimal model is suitable for capturing the local intermolecular packing of POSS cubes at short separation distances.

Long time atomistic polymer trajectories from coarse grained simulations: bisphenol-A polycarbonate

Based on coarse grained simulations of a specially adapted model for bisphenol-A polycarbonate (BPA-PC) we generate by inverse mapping, the reintroduction of chemical details, well equilibrated all-atom conformations and time trajectories of dense polymeric melts for up to 7.8 µs. This is several orders of magnitude more than any direct all-atom simulations have reached so far. These polymer melts contain up to 68600 atoms in = 100 chains of molecular weight = 5217. By comparison with short all-atom simulations we show that these trajectories are physically meaningful, providing us with a powerful tool to compare long time simulations to experiments, which probe specific local dynamics on long time scales, such as NMR relaxation.

The role of temperature and density on the glass-transition dynamics of glass formers

A correlation between the monomeric volume and the dynamic quantity E*(V)/H*, used to provide a quantitative measure of the role of temperature and density on the dynamics, is demonstrated for a series of polymers and glass-forming liquids. We show that monomeric volume and local packing play a key role in controlling the value of this ratio and thus the dynamics associated with the glass temperature.

Influence of density and temperature on the microscopic structure and the segmental relaxation of polybutadiene

We investigate the influence of temperature and density on the local structure and the dynamics of polybutadiene by controlling both hydrostatic pressure and temperature in polarized neutron diffraction experiments on deuterated polybutadiene and in inelastic incoherent scattering experiments on protonated polybutadiene. We observe that the static structure factor S(Q) does not change along macroscopic isochores. This behavior is contrary to the relaxations observed on the nanosecond and picosecond time scales and viewed by the dynamic incoherent scattering law S(Q,omega), which differ strongly along the same thermodynamic path. We conclude that the static behavior, i.e., S(Q), is dominated by macroscopic density changes, similar to the vibrational excitations in the meV range. However, the relaxation dynamics is more sensitive to thermal energy changes. This is confirmed by the finding that lines of identical relaxation behavior (in time, shape, and Q dependence), isochrones on the 10(-9) sec time scale, clearly cross the constant density lines in the (P,T) plane. Concerning S(Q), we can reasonably relate the variation of the main-peak position to the average neighbor chain distance and deduce crude microscopic thermal expansion and compressibility coefficients. In the low-Q regime, the observed pressure and temperature variation of S(Q) exceeds the compressibility contribution and suggests the existence of additional scattering, which might originate from structural correlations arising at higher temperature and low pressure.