Modeling entangled dynamics: comparison between stochastic single-chain and multichain models

To test the effectiveness of stochastic single-chain models in describing the dynamics of entangled polymers, we systematically compare one such model; the slip-spring model; to a multichain model solved using stochastic molecular dynamics(MD) simulations (the Kremer-Grest model). The comparison involves investigating if the single-chain model can adequately describe both a microscopic dynamical and a macroscopic rheological quantity for a range of chain lengths. Choosing a particular chain length in the slip-spring model, the parameter values that best reproduce the mean-square displacement of a group of monomers is determined by fitting toMDdata. Using the same set of parameters we then test if the predictions of the mean-square displacements for other chain lengths agree with the MD calculations. We followed this by a comparison of the time dependent stress relaxation moduli obtained from the two models for a range of chain lengths. After identifying a limitation of the original slip-spring model in describing the static structure of the polymer chain as seen in MD, we remedy this by introducing a pairwise repulsive potential between the monomers in the chains. Poor agreement of the mean-square monomer displacements at short times can be rectified by the use of generalized Langevin equations for the dynamics and resulted in significantly improved agreement.

Large deviations of Rouse polymer chain: first passage problem

The purpose of this paper is to investigate several analytical methods of solving first passage (FP) problem for the Rouse model, a simplest model of a polymer chain. We show that this problem has to be treated as a multi-dimensional Kramers' problem, which presents rich and unexpected behavior. We first perform direct and forward-flux sampling (FFS) simulations, and measure the mean first-passage time $\tau(z)$ for the free end to reach a certain distance $z$ away from the origin. The results show that the mean FP time is getting faster if the Rouse chain is represented by more beads. Two scaling regimes of $\tau(z)$ are observed, with transition between them varying as a function of chain length. We use these simulations results to test two theoretical approaches. One is a well known asymptotic theory valid in the limit of zero temperature. We show that this limit corresponds to fully extended chain when each chain segment is stretched, which is not particularly realistic. A new theory based on the well known Freidlin-Wentzell theory is proposed, where dynamics is projected onto the minimal action path. The new theory predicts both scaling regimes correctly, but fails to get the correct numerical prefactor in the first regime. Combining our theory with the FFS simulations lead us to a simple analytical expression valid for all extensions and chain lengths. One of the applications of polymer FP problem occurs in the context of branched polymer rheology. In this paper, we consider the arm-retraction mechanism in the tube model, which maps exactly on the model we have solved. The results are compared to the Milner-McLeish theory without constraint release, which is found to overestimate FP time by a factor of 10 or more.

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.

Oestrogenic activity of p-hydroxybenzoic acid (common metabolite of paraben esters) and methylparaben in human breast cancer cell lines

This paper addresses the question of whether p-hydroxybenzoic acid, the common metabolite of parabens, possesses oestrogenic activity in human breast cancer cell lines. The alkyl esters of p-hydroxybenzoic acid (parabens) are used widely as preservatives in consumer products to which the human population is exposed and have been shown previously to possess oestrogenic activity and to be present in human breast tumour tissue, which is an oestrogen-responsive tissue. Recent work has shown p-hydroxybenzoic acid to give an oestrogenic response in the rodent uterotrophic assay. We report here that p-hydroxybenzoic acid possesses oestrogenic activity in a panel of assays in human breast cancer cell lines. p-Hydroxybenzoic acid was able to displace [H-3]oestradiol from cytosolic oestrogen receptor of MCF7 human breast cancer cells by 54% at 5 x 10(6)-fold molar excess and by 99% at 10(7)-fold molar excess. It was able to increase the expression of a stably integrated oestrogen responsive reporter gene (ERE-CAT) at a concentration of 5 x 10(-4) M in MCF7 cells after 24 h and 7 days, which could be inhibited by the anti-oestrogen ICI 182 780 (Faslodex, fulvestrant). Proliferation of two human breast cancer cell lines (MCF7, ZR-75-1) could be increased by 10(-5) M p-hydroxybenzoic acid. Following on from previous studies showing a decrease in oestrogenic activity of parabens with shortening of the linear alkyl chain length, this study has compared the oestrogenic activity of p-hydroxybenzoic acid where the alkyl grouping is no longer present with methylparaben, which has the shortest alkyl group. Intrinsic oestrogenic activity of p-hydroxybenzoic acid was similar to that of methylparaben in terms of relative binding to the oestrogen receptor but its oestrogenic activity on gene expression and cell proliferation was lower than that of methylparaben. It can be concluded that removal of the ester group from parabens does not abrogate its oestrogenic activity and that p-hydroxybenzoic acid can give oestrogenic responses in human breast cancer cells. Copyright (C) 2005 John Wiley & Sons, Ltd.

Oestrogenic activity of benzylparaben

Previous work has demonstrated that the alkyl esters of p-hydroxybenzoic acid (parabens) possess oestrogenic activity, which increases with length of alkyl chain from methylparaben to n-butylparaben and with branching in the alkyl chain from n-butylparaben to isobutylparaben. This study reports on the oestrogenic activity of benzylparaben in a variety of assays in vitro and in vivo. Benzylparaben was able to displace [H-3]oestradiol from cytosolic oestrogen receptor (ER) of MCF7 human breast cancer cells by 22% at 1000-fold molar excess, by 40% at 10000-fold molar excess, by 57% at 100000-fold molar excess and by 100% at 1000000-fold molar excess. It was able to increase expression of a stably transfected oestrogen responsive reporter gene (ERE-CAT) in MCF7 cells after 24 h at 10(-5)M/10(-4)M and after 7 days at 10(-6)M/10(-5)M/10(-4)M. Proliferation of MCF7 cells could be increased by 10(-6)M/10(-5)M benzylparaben and this could be inhibited by 10(-7)M pure anti-oestrogen ICI 182,780, indicating that growth effects were ER mediated. Further evidence for ER-mediation was provided from the ability of benzylparaben to increase the growth of a second oestrogen-dependent human breast cancer cell line ZR-75-1, but not the oestrogen-insensitive NIDA-MB-231 cell line. When tested in the presence of 10(-10)M 17beta-oestradiol, benzylparaben gave no antagonist response on the growth of either MCF7 or ZR-75-1 cells. Finally, benzylparaben could increase uterine weight in the immature mouse following topical application of three daily doses of 33 mg to dorsal skin. These results demonstrate that the oestrogenicity of methylparaben can be increased by the addition of an aryl group as well as by lengthening or branching the alkyl grouping. Copyright (C) 2003 John Wiley Sons, Ltd.

A mechanistic study of cyclic siloxane pyrolyses at low pressures

Matrix isolation IR spectroscopy has been used to study the vacuum pyrolysis of hexamethylcyclotrisiloxane (D3), octamethylcyclotetrasiloxane (D4) and decamethyl cyclopentasiloxane (D5), and the results interpreted in the context of various kinetic models. In particular, it is shown that the significant pyrolysis products - which include CH3, CH4, C2H2, C2H4, C2H6 and SiO - may be satisfactorily accounted for by radical reactions involving dimethylsiloxane (D1), and estimates are made of the various chain lengths for the proposed reactions based on a range of ambient conditions.

Self-assembled columns of fullerene

Columnar mesophases based on alternating triphenylene and hexaphenyltriphenylene moieties are exceptionally stable and able to accommodate bulky side-chain substituents within the alkyl chain continuum between the columns. This paper presents a system in which the triphenylene bears a fullerene on its side-chain and the hexaphenyltriphenylene equivalent is the aza-derivative hexakis(4-nonylphenyl)dipyrazino[2,3-f : 2'3'-h] quinoxalene, PDQ9. The mesophase formed was identified as hexagonal columnar (Col(h)) by X-ray diffraction (a = 25.2 angstrom and c = 3.5 angstrom) but, in addition to the expected peaks, there is indication of a two-dimensional hexagonal superlattice with d-spacing 59 angstrom. This superlattice is believed to arise from ordering of the fullerenes within the liquid crystal matrix. It can be explained on the assumption that, to maximise fullerene-fullerene contact, the fullerenes form chains which wrap around the central column in every group of seven columns of the triphenylene : PDQ9 Col(h) array.

Dicopper(II) complexes of a new di-para-xylyldioxatetraazamacrocycle and cascade species with dicarboxylate anions: thermodynamics and structural properties

The new dioxatetraazamacrocycle (L-1) was synthesized by a 2 + 2 condensation and characterized. Stability constants of its copper(II) complexes were determined by spectrophotometry in DMSO at 298.2 K in 0. 10 mol dm(-3) KClO4. Mainly dinuclear complexes are formed and the presence of mononuclear species is dependent on the counterion (Cl- or ClO4-). The association constants of the dinuclear copper(II) complexes with dicarboxylate anions [oxalate (oxa(2-)), malonate (mal(2-)), succinate (suc(2-)), and glutarate (glu(2-))] were also determined by spectrophotometry at 298.2 K in DMSO, and it was found that values decrease with an increase of the alkyl chain between the carboxylate groups. X-Band EPR spectra of the dicopper(II) complexes and of their cascade species in frozen DMSO exhibit dipole-dipole coupling, and their simulation, together with their UV-vis spectra, showed that the copper centres of the complexes in solution had square pyramidal geometries though with different distortions. From the experimental data, it was also possible to predict the Cu...Cu distances, the minimum being found at 6.4 angstrom for the (Cu2LCl4)-Cl-1 complex and then successively this distance slightly increases when the chloride anions are replaced by dicarboxylate anions, from 6.6 angstrom for oxa(2-) to 7.8 for glu(2-). The crystal structures of the dinuclear copper cascade species with oxa(2-) and suc(2-) were determined and showed one anion bridging both copper centres and Cu...Cu distances of 5.485(7) angstrom and 6.442(8) angstrom, respectively.

Cascade dicopper architectures of a dibenzodioxatetraazamacrocycle

The dibenzodioxatetraazamacrocycle [26]pbz(2)N(4)O(2) was characterised by single crystal X-ray diffraction and the protonation constants of this compound and the stability constants of its copper(II) and lead(II) complexes were determined by potentiometry in water at 298.2 K in 0.10 mol dm(-3) in KNO3. Mono- and dinuclear complexes were found for both metal ions, the dinuclear complexes being the main species in the 5-7.5 pH range for copper(II) and 7.5-8.5 for lead(II). As expected the values of the stability constants for the copper(II) complexes are lower than those for related macrocycles containing only nitrogen atoms. The presence of mono- and dinuclear copper complexes was also confirmed by electrospray ionization mass spectrometry. These results suggest that the symmetric macrocyclic cavity of [26]pbZ(2)N(4)O(2) has enough space for the coordination of two metal ions. Additionally, NMR spectroscopy showed that the dinuclear complex of lead(II) has high symmetry. The equilibrium constants of the dinuclear copper(II) complexes and dicarboxylate anions (oxalate, malonate and succinate) were also determined in 0.10 mol dm-3 aqueous KNO3 solution. Only species containing one anion, Cu(2)H(h)LA((2+h)), were found, strongly suggesting that the anion bridges the two copper(II) ions. The binding constants of the cascade species formed by [Cu-2[26]pbZ(2)N(4)O(2)(H2O)(4+) with dicarboxylate anions decrease with the increase in length of the alkyl chain of the anion, a fact which was attributed to a higher conformational energy necessary for the rearrangement of the macrocycle to accommodate the larger anions bridging the two copper(II) centres. The variation of the magnetic susceptibility with temperature Of [Cu-2(H-2[26]pbz(2)N(4)O(2))(oxa)(3)]-4H(2)O and [Cu-2([26]pbz(2)N(4)O(2))(suc)Cl-2] were measured and the two complexes showed different behaviour. (c) 2007 Elsevier Ltd. All rights reserved.

Ultrasonic chemical oxidative degradations of 1,3-dialkylimidazolium ionic liquids and their mechanistic elucidations

A highly efficient process for oxidative degradation of 1,3-dialkylimidazolium ionic liquids in hydrogen peroxide/acetic acid aqueous medium assisted by ultrasonic chemical irradiation is, for the first time, described. It is shown that more than 93% of the 1,3-dialkylimidazolium cation with the corresponding Cl-, Br-, BF4- and PF6- counter-anions at a concentration of 2.5 mM can be degraded at 50 degrees C within 12 h while at 72 h the conversions approach 99%. A tentative mechanism for the degradation of these ILs is for the first time proposed through a detailed kinetic analysis of several characteristic transients and/or immediate products, which are identified during the ILs degradation using GC-MS. The results clearly indicate that three hydrogen atoms in the imidazolium ring are the first sites preferably oxidized, followed by cleavage of the alkyl groups attached to the N atoms from the ring. The nature of the alkyl chain length on the imidazolium ring and the type of counter anion do not seem to affect the degradation process. Further, selective fragmentations of C-N bonds of the imidazolium or derived ring lead to ring opening, forming degraded intermediates. It is also shown that acetoxyacetic acid and biurea are the final kinetically stable degraded products from the ILs under the degradation conditions.

Structure, rheology and shear alignment of Pluronic block copolymer mixtures

The structure and flow behaviour of binary mixtures of Pluronic block copolymers P85 and P123 is investigated by small-angle scattering, rheometry and mobility tests. Micelle dimensions are probed by dynamic light scattering. The micelle hydrodynamic radius for the 50/50 mixture is larger than that for either P85 or P123 alone, Clue to the formation of mixed micelles with a higher association number. The phase diagram for 50/50 mixtures contains regions Of Cubic and hexagonal phases similar to those for the parent homopolymers, however the region of stability of the cubic phase is enhanced at low temperature and concentrations above 40 wt%. This is ascribed to favourable packing of the mixed micelles containing core blocks with two different chain lengths, but similar corona chain lengths. The shear flow alignment of face-centred cubic and hexagonal phases is probed by in situ small-angle X-ray or neutron scattering with simultaneous rheology. The hexagonal phase can be aligned using steady shear in a Couette geometry, however the high modulus Cubic phase cannot be aligned well in this way. This requires the application of oscillatory shear or compression. (C) 2008 Elsevier Inc. All rights reserved.

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.

Segmental dynamics in entangled linear polymer melts

We present molecular dynamics (MD) and slip-springs model simulations of the chain segmental dynamics in entangled linear polymer melts. The time-dependent behavior of the segmental orientation autocorrelation functions and mean-square segmental displacements are analyzed for both flexible and semiflexible chains, with particular attention paid to the scaling relations among these dynamic quantities. Effective combination of the two simulation methods at different coarse-graining levels allows us to explore the chain dynamics for chain lengths ranging from Z ≈ 2 to 90 entanglements. For a given chain length of Z ≈ 15, the time scales accessed span for more than 10 decades, covering all of the interesting relaxation regimes. The obtained time dependence of the monomer mean square displacements, g1(t), is in good agreement with the tube theory predictions. Results on the first- and second-order segmental orientation autocorrelation functions, C1(t) and C2(t), demonstrate a clear power law relationship of C2(t) C1(t)m with m = 3, 2, and 1 in the initial, free Rouse, and entangled (constrained Rouse) regimes, respectively. The return-to-origin hypothesis, which leads to inverse proportionality between the segmental orientation autocorrelation functions and g1(t) in the entangled regime, is convincingly verified by the simulation result of C1(t) g1(t)−1 t–1/4 in the constrained Rouse regime, where for well-entangled chains both C1(t) and g1(t) are rather insensitive to the constraint release effects. However, the second-order correlation function, C2(t), shows much stronger sensitivity to the constraint release effects and experiences a protracted crossover from the free Rouse to entangled regime. This crossover region extends for at least one decade in time longer than that of C1(t). The predicted time scaling behavior of C2(t) t–1/4 is observed in slip-springs simulations only at chain length of 90 entanglements, whereas shorter chains show higher scaling exponents. The reported simulation work can be applied to understand the observations of the NMR experiments.

In vitro antioxidant activity and antigenotoxicity of 5-n-alkylresorcinols

Incubation with 5-n-alkylresorcinols (chain lengths C15:0, C17:0, C19:0, C21:0, and C23:0) increased the self-protection capacity of HT29 human colon cancer cells against DNA damage induced by hydrogen peroxide and genotoxic fecal water samples using comet assay (single-cell gel electrophoresis assay). The alkylresorcinols did not exert potent antioxidant activity in the FRAP (ferric reduction ability of plasma) and DPPH (2,2-diphenyl-1-picrylhydrazyl) radical assays. However they were able to significantly inhibit copper-mediated oxidation of human LDL (low-density lipoprotein) in vitro, and pentadecylresorcinol at 25 micromol/L increased lag time by 65 min. The results show that alkylresorcinols have antigenotoxic and antioxidant activity under in vitro conditions.