Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

Hierarchical modeling of entangled polymers

We demonstrate that it is possible to link multi-chain molecular dynamics simulations with the tube model using a single chain slip-links model as a bridge. This hierarchical approach allows significant speed up of simulations, permitting us to span the time scales relevant for a comparison with the tube theory. Fitting the mean-square displacement of individual monomers in molecular dynamics simulations with the slip-spring model, we show that it is possible to predict the stress relaxation. Then, we analyze the stress relaxation from slip-spring simulations in the framework of the tube theory. In the absence of constraint release, we establish that the relaxation modulus can be decomposed as the sum of contributions from fast and longitudinal Rouse modes, and tube survival. Finally, we discuss some open questions regarding possible future directions that could be profitable in rendering the tube model quantitative, even for mildly entangled polymers

Relative proportions of polycyclic aromatic hydrocarbons differ between accumulation bioassays and chemical methods to predict bioavailability

Chemical methods to predict the bioavailable fraction of organic contaminants are usually validated in the literature by comparison with established bioassays. A soil spiked with polycyclic aromatic hydrocarbons (PAHs) was aged over six months and subjected to butanol, cyclodextrin and tenax extractions as well as an exhaustive extraction to determine total PAH concentrations at several time points. Earthworm (Eisenia fetida) and rye grass root (Lolium multiflorum) accumulation bioassays were conducted in parallel. Butanol extractions gave the best relationship with earthworm accumulation (r2 ≤ 0.54, p ≤ 0.01); cyclodextrin, butanol and acetone–hexane extractions all gave good predictions of accumulation in rye grass roots (r2 ≤ 0.86, p ≤ 0.01). However, the profile of the PAHs extracted by the different chemical methods was significantly different (p < 0.01) to that accumulated in the organisms. Biota accumulated a higher proportion of the heavier 4-ringed PAHs. It is concluded that bioaccumulation is a complex process that cannot be predicted by measuring the bioavailable fraction alone. The ability of chemical methods to predict PAH accumulation in Eisenia fetida and Lolium multiflorum was hindered by the varied metabolic fate of the different PAHs within the organisms.

Microscopic Mechanisms of Strain Hardening in Glassy Polymers

The mechanisms underlying the increase in stress for large mechanical strains of a polymer glass, quantified by the strain-hardening modulus, are still poorly understood. In the present paper we aim to elucidate this matter and present new mechanisms. Molecular-dynamics simulations of two polymers with very different strain-hardening moduli (polycarbonate and polystyrene) have been carried out. Nonaffine displacements occur because of steric hindrances and connectivity constraints. We argue that it is not necessary to introduce the concept of entanglements to understand strain hardening, but that hardening is rather coupled with the increase in the rate of nonaffine particle displacements. This rate increases faster for polycarbonate, which has the higher strain-hardening modulus. Also more nonaffine chain stretching is present for polycarbonate. It is shown that the inner distances of such a nonaffinely deformed chain can be well described by the inner distances of the worm-like chain, but with an effective stiffness length (equal to the Kuhn length for an infinite worm-like chain) that increases during deformation. It originates from the finite extensibility of the chain. In this way the increase in nonaffine particle displacement can be understood as resulting from an increase in the effective stiffness length of the perturbed chain during deformation, so that at larger strains a higher rate of plastic events in terms of nonaffine displacement is necessary, causing in turn the observed strain hardening in polymer glasses.

Uptake Pathways of Polycyclic Aromatic Hydrocarbons in White Clover

An understanding of the primary pathways of plant uptake of organic pollutants is important to enable the risks from crops grown on contaminated soils to be assessed. A series of experiments were undertaken to quantify the importance of the pathways of contamination and the Subsequent transport within the plant using white clover plants grown in solution culture. Root uptake was primarily an absorption process, but a component of the contamination was a result of the transpiration flux to the shoot for higher Solubility compounds. The root contamination can be easily predicted using a simple relationship with K-OW, although if a composition model was used based on lipid content, a significant under prediction of the contamination was observed. Shoot uptake was driven by the transpiration stream flux which was related to the solubility of the individual PAH rather than the K-OW. However, the experiment was over a short duration, 6 days, and models based on K-OW may be better for crops grown in the field where the vegetation will approach equilibrium and transpiration cannot easily be measured, A significant fraction of the shoot contamination resulted from aerial deposition derived from volatilized PAH. This pathway was more significant for compounds approaching log K-OA > 9 and log K-AW < -3. The shoot uptake pathways need further investigation to enable them to be modeled separately, There was no evidence of significant systemic transport of the PAR so transfer outside the transpiration stream is likely to be limited.

Bile metabolites of polycyclic aromatic hydrocarbons (PAHs) in European eels Anguilla anguilla from United Kingdom estuaries

A total of 94 European eels (Anguilla anguilla) were collected from five estuaries in the UK. The deconjugated metabolites of polycyclic aromatic hydrocarbons (PAHs) in the bile of the eels were separated using HPLC. Six PAH metabolites were identified: 1-hydroxy (1-OH) metabolites of phenanthrene, pyrene and chrysene; and the 1-OH, 3-OH and 7,8 dihydrodiol metabolites of benzo[a]pyrene (BaP). The mean concentration of the six metabolites was greatest in eels from the Tyne (49 muM) followed by the Wear (33 muM), Tees (19 muM), Thames (4 muM) and Severn (2 muM) estuaries. Although 1-OH pyrene was always the dominant compound, there were significant differences (P<0.05) between sites and between estuaries for some metabolites. Normalising the molar concentration of the bile metabolites to the bile biliverdin absorbance reduced sample variation. When the metabolites identified were-each expressed as a percentage of the total detected, the metabolite profile was characteristic for each estuary. (C) 2002 Elsevier Science B.V. All rights reserved.

Linear nucleation of polymers

The use of high-melting fibres as linear nuclei for quiescent polymeric melts is instrumental in providing the superior mechanical properties of polymeric self-composites. It also has inherent advantages in the elucidation of fundamental aspects of polymeric crystallization and self-organization, not least in allowing systematic microscopic studies of polymeric crystallization from nucleation through to the growth interface. This has demonstrated explicitly that lamellae develop in two distinct ways, for slower and faster growth, depending on whether fold packing has or has not time to order before the next molecular layer is added with only the former leading to banded growth in linear polyethylene. Other gains in understanding concern cellulation and morphological instability, internuclear interference, isothermal lamellar thickening and banded growth being a consequence of the partial relief of initial surface stress. (c) 2006 Elsevier Ltd. All rights reserved.