Computer simulation studies of the influence of side alkyl chain on glass transition behavior of carbazole trimer

In this work,all-atom molecular dynamics simulations were employed to study the influence of the side alkyl chain on glass transition behavior of several carbazole trimers(CT) in a temperature range from 423 to 183 K.The glass transition temperatures were obtained from the break in the slope of the volume-temperature curves and found to agree with the experimental values.The short time dynamics of four CT molecules were probed by usingvelocity autocorrelation functions and mean-square displacements.The current studies showed that the dynamics of CT systems can be easily interpreted through the cage effect.Furthermore,the investigation of the torsional autocorrelation function and P_(2-state)/P_(3-state) functions showed that the rotational barriers of side chains can slow down the conformational relaxation and lead to stronger temperature dependence of conformational relaxation.The relaxation time,characteristic time of P_(2-state)(t) and P_(3-state)(t) functions were all found to have Arrhenius-type temperature dependence.     

The role of the torsional potential in relaxation dynamics: a molecular dynamics study of polyethylene

The role of the torsional potential in bulk polymer chain dynamics is investigated via molecular dynamics simulation using polyethylene as a model system. A number of three-fold barrier values, both greater and less than the standard one, were invoked. The one-fold potential that determines the gauche vs trans energy difference was also varied. For each of the selected torsional potentials, the MD volumetric glass transition temperature, T_g, was located. It was found that T_g is quite sensitive to the three-fold barrier magnitude, moving from below 100 K to nearly 400 K as the barrier goes from zero to twice the standard value. However T_g was found to be quite insensitive to the gauche trans energy difference. Details of the conformational dynamics were studied for the case of a zero torsional potential. This included the rate and location of conformational transitions, the decay of the torsional angle autocorrelation function (ACF) and the cooperativity of conformational transitions, all as a function of temperature. The temperature dependence of the conformational transition rate remains Arrhenius at all temperatures. The relaxation time characterizing the torsional angle ACF decay exhibits WLF temperature behavior. The conformational transitions are randomly distributed over the bonds at high temperature, but near T_g they become spatially heterogeneous and localized. The transitions show next-neighbor correlation as well as self-correlated forward-backward transitions. All of these features are similar to those found in previous simulations under the standard torsional potential.     

Ultrafast dynamics of polybutadiene probed by optically heterodyne-detected optical-Kerr-effect spectroscopy

The ultrafast dynamics of polybutadiene have been studied using ultrafast optical-Kerr-effect spectroscopy. The data are compared with measurements on 1,3- and 1,4-pentadiene. The two diene derivatives have quite distinct subpicosecond dynamics, with an important contribution from an intramolecular torsional mode in the 1,4-derivative. The main part of the polymer spectral density can be assigned, by analogy with the data for 1,4-pentadiene, to intramolecular torsional motion about carbon–carbon single bonds. The picosecond diffusive orientational relaxation times of the dienes are not well described by simple hydrodynamics.     

Nuclear magnetic resonance investigation of dynamics in poly(ethylene oxide)-based lithium polyether-ester-sulfonate ionomers

Nuclear magnetic resonance spectroscopy has been utilized to investigate the dynamics of poly(ethylene oxide)-based lithium sulfonate ionomer samples that have low glass transition temperatures. (1)H and (7)Li spin-lattice relaxation times (T(1)) of the bulk polymer and lithium ions, respectively, were measured and analyzed in samples with a range of ion contents. The temperature dependence of T(1) values along with the presence of minima in T(1) as a function of temperature enabled correlation times and activation energies to be obtained for both the segmental motion of the polymer backbone and the hopping motion of lithium cations. Similar activation energies for motion of both the polymer and lithium ions in the samples with lower ion content indicate that the polymer segmental motion and lithium ion hopping motion are correlated in these samples, even though lithium hopping is about ten times slower than the segmental motion. A divergent trend is observed for correlation times and activation energies of the highest ion content sample with 100% lithium sulfonation due to the presence of ionic aggregation. Details of the polymer and cation dynamics on the nanosecond timescale are discussed and complement the findings of X-ray scattering and quasi-elastic neutron scattering experiments.     

Investigations of the migration of elements in tree rings by Laser-ICP-MS

The autocorrelation function of the concentration of elements in growth rings has been used to study their migration from ring to ring. The correlation time which can be calculated by the autocorrelation function gives an indication of the transport of elements. However, a large correlation time might only be taken as a necessary criterion for the migration but not for a sufficient one. Laser sampling-ICP-MS has been used for the determination of the contents of elements in the rings of an 150 years old oak - the autocorrelation functions of selected elements have been calculated. Large correlation times have been only obtained for Ca, Fe and Cr. It may be concluded that Mg, Zn, Cu, Sr, Pb, Ni, Ba, Al and Mn hardly migrate radially in the trunk of an oak.     

Investigations of the migration of elements in tree rings by Laser-ICP-MS

The autocorrelation function of the concentration of elements in growth rings has been used to study their migration from ring to ring. The correlation time which can be calculated by the autocorrelation function gives an indication of the transport of elements. However, a large correlation time might only be taken as a necessary criterion for the migration but not for a sufficient one. Laser sampling-ICP-MS has been used for the determination of the contents of elements in the rings of an 150 years old oak – the autocorrelation functions of selected elements have been calculated. Large correlation times have been only obtained for Ca, Fe and Cr. It may be concluded that Mg, Zn, Cu, Sr, Pb, Ni, Ba, Al and Mn hardly migrate radially in the trunk of an oak.     

Molecular dynamics simulation of the polymer electrolyte poly(ethylene oxide)/LiClO(4). II. Dynamical properties

The dynamical properties of the polymer electrolyte poly(ethylene oxide) (PEO)LiClO(4) have been investigated by molecular dynamics simulations. The effect of changing salt concentration and temperature was evaluated on several time correlation functions. Ionic displacements projected on different directions reveal anisotropy in short-time (rattling) and long-time (diffusive) dynamics of Li(+) cations. It is shown that ionic mobility is coupled to the segmental motion of the polymeric chain. Structural relaxation is probed by the intermediate scattering function F(k,t) at several wave vectors. Good agreement was found between calculated and experimental F(k,t) for pure PEO. A remarkable slowing down of polymer relaxation is observed upon addition of the salt. The ionic conductivity estimated by the Nernst-Einstein equation is approximately ten times higher than the actual conductivity calculated by the time correlation function of charge current.     

The dynamics of single chains within a model polymer melt

Discontinuous molecular dynamics simulations are performed on a system containing 32 hard chains of length 192 at a volume fraction of phi = 0.45 to explore the idea that localized entanglements have a significant effect on the dynamics of the individual chains within an entangled polymer melt. Anomalous behavior can still be observed when studying the dynamics of the individual chains, although increased time averaging causes the anomalous relaxation-memory-release behavior that was observed previously in the system to smooth out. First, the individual chain mean squared displacements and apparent diffusion coefficients are calculated, and a wide distribution of diffusive behavior is found. Although the apparent diffusion coefficient curve averaged over all chains displays the predicted long-time diffusive behavior, the curves for the individual chains differ both qualitatively and quantitatively. They display superdiffusive, diffusive, and subdiffusive behavior, with the largest percentage of chains exhibiting superdiffusive behavior and the smallest percentage exhibiting the predicted diffusive behavior. Next, the individual chain end-to-end vector autocorrelation functions and relaxation times are determined, and a wide distribution of stress relaxation behavior is found. The times when the end-to-end vector autocorrelation functions relax completely span almost an order of magnitude in reduced time. For some chains, the end-to-end vector autocorrelation function relaxes smoothly toward zero similar to the system average; however, for other chains the relaxation is slowed greatly, indicating the presence of additional entanglements. Almost half of the chains exhibit the anomalous behavior in the end-to-end vector autocorrelation function. Finally, the dynamic properties are displayed for a single chain exhibiting anomalous relaxation-memory-release behavior, supporting the idea that the relaxation-memory-release behavior is a single-chain property.     

Equilibrium conformational dynamics of a polymer in a solvent

Molecular dynamics simulations were used to study the conformational dynamics of a bead-spring model polymer in an explicit solvent under good solvent conditions. The dynamics of the polymer chain were investigated using an analysis of the time autocorrelation functions of the Rouse coordinates of the polymer chain. We have investigated the variation of the correlation functions with polymer chain length N, solvent density rho, and system size. The measured initial decay rates gamma(p) of the correlation functions were compared with the predictions from a theory of polymer dynamics which uses the Oseen tensor to describe hydrodynamic interactions between monomers. Over the range of chain lengths considered (N = 30-60 monomers), the predicted scaling of gamma(p) proportional to N(-3nu) was observed at high rho, where nu is the polymer scaling exponent. The predicted gamma(p) are generally higher than the measured values. This discrepancy increases with decreasing rho, as a result in the breakdown in the conditions required for the Oseen approximation. The agreement between theory and simulation at high rho improves considerably if the theoretical expression for gamma(p) is modified to avoid sum-to-integral approximations, and if the values of (R(p)2), which are used in the theory, are taken directly from the simulation rather than being calculated using approximate scaling relations. The observed finite-size scaling of gamma(p) is not quantitatively consistent with the theoretical predictions.     

On the short-time limit of ring polymer molecular dynamics

We examine the short-time accuracy of a class of approximate quantum dynamical techniques that includes the centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD) methods. Both of these methods are based on the path integral molecular dynamics (PIMD) technique for calculating the exact static equilibrium properties of quantum mechanical systems. For Kubo-transformed real-time correlation functions involving operators that are linear functions of positions or momenta, the RPMD and (adiabatic) CMD approximations differ only in the choice of the artificial mass matrix of the system of ring polymer beads that is employed in PIMD. The obvious ansatz for a general method of this type is therefore to regard the elements of the PIMD (or Parrinello-Rahman) mass matrix as an adjustable set of parameters that can be chosen to improve the accuracy of the resulting approximation. We show here that this ansatz leads uniquely to the RPMD approximation when the criterion that is used to select the mass matrix is the short-time accuracy of the Kubo-transformed correlation function. In particular, we show that the leading error in the RPMD position autocorrelation function is O(t(8)) and the error in the velocity autocorrelation function is O(t(6)), for a general anharmonic potential. The corresponding errors in the CMD approximation are O(t(6)) and O(t(4)), respectively.     

Ion and polymer dynamics in polymer electrolytes PPO-LiClO4. II. 2H and 7Li NMR stimulated-echo experiments

We use 2H NMR stimulated-echo spectroscopy to measure two-time correlation functions characterizing the polymer segmental motion in polymer electrolytes PPO-LiClO4 near the glass transition temperature Tg. To investigate effects of the salt on the polymer dynamics, we compare results for different ether oxygen to lithium ratios, namely, 6:1, 15:1, 30:1, and infinity. For all compositions, we find nonexponential correlation functions, which can be described by a Kohlrausch function. The mean correlation times show quantitatively that an increase of the salt concentration results in a strong slowing down of the segmental motion. Consistently, for the high 6:1 salt concentration, a high apparent activation energy Ea=4.1 eV characterizes the temperature dependence of the mean correlation times at Tg相似文献   

Dynamics of polymers in a particle-based mesoscopic solvent

We study the dynamics of flexible polymer chains in solution by combining multiparticle-collision dynamics (MPCD), a mesoscale simulation method, and molecular-dynamics simulations. Polymers with and without excluded-volume interactions are considered. With an appropriate choice of the collision time step for the MPCD solvent, hydrodynamic interactions build up properly. For the center-of-mass diffusion coefficient, scaling with respect to polymer length is found to hold already for rather short chains. The center-of-mass velocity autocorrelation function displays a long-time tail which decays algebraically as (Dt)(-3/2) as a function of time t, where D is the diffusion coefficient. The analysis of the intramolecular dynamics in terms of Rouse modes yields excellent agreement between simulation data and results of the Zimm model for the mode-number dependence of the mode-amplitude correlation functions.     

Analysis of uncertainties in polymer viscoelastic properties obtained from equilibrium computer simulations

We critically evaluate the uncertainties in the stress autocorrelation function obtained from equilibrium molecular dynamics simulation of model polymer melts. This quantity is central to evaluating transport properties, e.g., the complex modulus and the viscosity. In contrast to the intuitive expectation that simulations have to be run five to six orders of magnitude longer than the chain relaxation time to reduce uncertainties to acceptable levels, our analysis shows that the majority of the uncertainty is associated with rapidly oscillating bonded interactions. These fluctuations occur on time scales which are approximately 10(4) times shorter than the relaxation time of a chain of length 80. Consequently, the effects of these oscillations on the stress autocorrelation function can be dramatically reduced by (i) conducting long simulations (typically 10(6) times longer than the bond relaxation times or only 10(2) chain relaxation times) and (ii) by performing running averages with time windows whose time scales are much longer than these oscillations. Conducting such long simulations also allows for the accurate determination of the melt viscosity and the low-frequency complex modulus, but performing running averages do not impact these quantities since they are time integrals of the stress autocorrelation function.     

High Pressure Induced Coil-Globule Transitions of Smart Polymers

Summary: In current work time-resolved optical spectroscopy (TROS) has been used to study coil-globule transitions monitored by the local segmental dynamics of anthracene labeled poly (N-isopropymethacrylamide), PNIPMAM as a function of pressure (0.1 MPa–200 MPa) over a temperature range of 283 K to 333 K. The positions of temperature-induced transition were observed to be independent on molecular weight of polymer at low pressures. The positions of pressure-induced transition were observed to be dependent on molecular weight of polymer at temperatures below LCST at atmospheric pressure. Double globule-coil-globule transition was observed to occur with pressure increasing at temperatures nearly above LCST. All these results along with values of intrinsic viscosity evaluated from values of correlation times measured for globules formed at different pressure/temperature conditions suggest the different mechanisms of compactisation governed by pressure and temperature and, correspondently, the different types of final structures. At low pressures with temperature increasing the compact, well-packed globules are forming via initial interactions between neighboring parts of polymer chain and further collapse. Relatively loosened particles are forming with pressurizing at low temperatures. Interaction between remote along the chain units takes part from the first stage of globule formation. And finally, rather solvated and irregularly twisted particles are forming at high pressure and high temperatures, i.e. at conditions, when both processes are involved.     

Viscoelasticity of stretched polymer chains: Analytical theory and computer-aided simulation

The viscoelasticity of stretched polymer chains has been studied by the method of collisional dynamics. To this end, time correlation functions of the fluctuations of the microscopic stress tensor are modeled and relaxation moduli are expressed. Before, for stretched polymer networks, correlation functions used to be calculated in terms of an approximate theory that allowed one to estimate the strain dependences of loss modulus. The calculated dependences are shown to agree qualitatively with the results of measurements performed over a wide interval of strains, including prefracture strains. This theory is verified by comparing the time correlation functions of stress tensor fluctuations for a single stretched chain; these functions are found by computer-aided simulation and calculated on the basis of the existing analytical theory. In this case, a simple theory is adopted according to which a polymer molecule represents a chain composed of N atoms connected by freely jointed elastic bonds. The first and Nth atoms of this chain are attached by harmonic springs to immobile points located at a fixed distance. The decay of time correlation functions under study can be resolved into three stages. After a short initial interval provided by local motions, one can observe a region of power-law decay, which is followed by monoexponential decay at long times. The results of computer-aided simulation generally agree with the predictions of analytical theory. Certain discrepancies primarily concern the dependences of the exponent of power-law relaxation on the degree of chain stretching.     

Communication: a minimal model for the diffusion-relaxation backbone dynamics of proteins

We present a model for the local diffusion-relaxation dynamics of the C(α)-atoms in proteins describing both the diffusive short-time dynamics and the asymptotic long-time relaxation of the position autocorrelation functions. The relaxation rate spectra of the latter are represented by shifted gamma distributions, where the standard gamma distribution describes anomalous slow relaxation in macromolecular systems of infinite size and the shift accounts for a smallest local relaxation rate in macromolecules of finite size. The resulting autocorrelation functions are analytic for any time t ≥ 0. Using results from a molecular dynamics simulation of lysozyme, we demonstrate that the model fits the position autocorrelation functions of the C(α)-atoms exceptionally well and reveals moreover a strong correlation between the residue's solvent-accessible surface and the fitted model parameters.     

Orientational dynamics in the isotropic phase of a calamitic liquid-crystal model

We report a molecular dynamics simulation study on the isotropic phase of an idealized calamitic liquid crystal model with a length-to-width ratio of approximately 5-6. The study focuses on the characterization of single-particle and collective orientational dynamics on approaching the phase transition to the nematic phase. Recent experimental and simulation works have suggested that a power law behavior exists at relatively short times in the decay of the time derivative of the orientational correlation functions. Qualitatively, our simulation data are consistent with these findings. Both single-particle and collective time correlation function derivatives possess, in their respective log-log plots, a linear region at very short times, whose slope is essentially independent from the thermodynamic state. Nevertheless, the single-particle orientational correlation functions are better described by a function which is the sum of a fast exponential, an intermediate stretched-exponential and a slow exponential, while the collective orientational correlation functions are satisfactorily described by a sum of two exponentials, at higher density, or by just one exponential, at lower density.     

Investigation of benzene-hexafluorobenzene dynamics in liquid binary mixtures

The structure and microscopic dynamics of liquid mixtures of benzene and hexafluorobenzene at room temperature and several compositions have been studied by molecular-dynamics simulations. In this implementation we have rescaled the intermolecular H-F cross potential parameters obtained from the Lorentz-Berthelot combining rules, in order to avoid the substantial overestimation of the energy of mixing predicted by the model when the usual rules are employed. We found that a reduction in the strength of cross H-F interactions by 50% relative to the geometric mean is required in order to get a good agreement with experiments. Radial-angular pair-correlation functions between like and unlike species have been computed and analyzed, by comparing them with the correlations in the corresponding neat liquids. We have also studied the microscopic intermolecular momentum transfer, by computing the time correlation function between the initial velocity of a central molecule and later velocities of neighboring molecules. Structural and dynamical information extracted from the mentioned functions seem to be consistent with the picture of relatively long-lived benzene-hexafluorobenzene (Bz-Hf) complexes present in the mixtures, which would be responsible for the considerable perturbation of the structure in the first shell of like species, and would be moving within the liquid in a parallel face-to-face configuration. Using the tools developed originally to estimate hydrogen-bond lifetimes in liquids, we have computed the lifetimes of the Bz-Hf complexes as a function of the mixture composition, by two different methods: the direct time-averaging scheme and from the autocorrelation function of bond occupation numbers. The obtained lifetimes are strongly dependent on the scheme chosen to compute the characteristic times. We have obtained for the Bz-Hf dimer in solution, at room temperature, lifetimes in the range of 30-40 ps from averaging schemes and around 60-120 ps from autocorrelation function methods. In the latter case, the longest times correspond to the equimolar mixture.     

Ab initio molecular-dynamics study of supercritical carbon dioxide

Car-Parrinello molecular-dynamics simulations of supercritical carbon dioxide (scCO(2)) have been performed at the temperature of 318.15 K and at the density of 0.703 g/cc in order to understand its microscopic structure and dynamics. Atomic pair correlation functions and structure factors have been obtained and good agreement has been found with experiments. In the supercritical state the CO(2) molecule is marginally nonlinear, and thus possesses a dipole moment. Analyses of angle distributions between near neighbor molecules reveal the existence of configurations with pairs of molecules in the distorted T-shaped geometry. The reorientational dynamics of carbon dioxide molecules, investigated through first- and second-order time correlation functions, exhibit time constants of 620 and 268 fs, respectively, in good agreement with nuclear magnetic resonance experiments. The intramolecular vibrations of CO(2) have been examined through an analysis of the velocity autocorrelation function of the atoms. These reveal a red shift in the frequency spectrum relative to that of an isolated molecule, consistent with experiments on scCO(2). The results have also been compared to classical molecular-dynamics calculations employing an empirical potential.