Search

We present a computer simulation model of polymer melts representing each chain as one single particle. Besides the position coordinate of each particle, we introduce a parameter nij for each pair of particles i and j within a specified distance from each other. These numbers, called entanglement numbers, describe the deviation of the system of ignored coordinates from its equilibrium state for the given configuration of the centers of mass of the polymers. The deviations of the entanglement numbers from their equilibrium values give rise to transient forces, which, together with the conservative forces derived from the potential of mean force, govern the displacements of the particles. We have applied our model to a melt of C800H1602 chains at 450 K and have found good agreement with experiments and more detailed simulations. Properties addressed in this paper are radial distribution functions, dynamic structure factors, and linear as well as nonlinear rheological properties. [ABSTRACT FROM AUTHOR]

It has been suggested that, due to topological constraints, rings in the melt may assume a more compact shape than Gaussian chains. In this paper, we exploit the availability of narrow fractions of perdeuterated linear and cyclic polydimethylsiloxane (PDMS) and, through the analysis of the small angle neutron scattering (SANS) profiles, demonstrate the difference in scattering properties of linear and cyclic PDMS molecules. As expected for Gaussian chains, for the H/D linear PDMS samples, log-log plots of the scattered intensity versus scattering vector Q display a Q(-2) dependence. However, for H/D cyclic blends, the scaling exponent is higher than 2, as predicted by computer simulations reported in the literature. We show that cyclic molecules in bulk display the characteristic maximum in plots of scattered intensity versus Q(-2) that is expected on the basis of Monte Carlo calculations and from the Casassa equation [E. F. Casassa, J. Polym. Sci. A 3, 605 (1965)]. It is also shown that, for rings, the Debye equation [P. Debye, J. Appl. Phys. 15, 338 (1944)] is no longer appropriate to describe the SANS profiles of H/D cyclic blends, at least up to Mw≈10 000. For these samples, the Casassa form factor gives a better representation of the SANS data and we show that this function which was developed for monodisperse cyclics is still adequate to describe our slightly polydisperse samples. Deviations from all above observations are noted for Mw>11 000 and are attributed to partial contamination of cyclic samples with linear chains. The failure of both the Debye and the Casassa form factors could be due to contamination of the cyclic fractions by linear polymers or to a real conformational change. [ABSTRACT FROM AUTHOR]

We focus on highly grafted binary polymer brushes with compatible components in the cases of different chain lengths. Layered structures parallel to the surface that indicating “phase separation” are observed in a series of dissipative particle dynamics simulations. The stretch parameters indicate that the short chains are suppressed in the lower layer of the film, whereas the longer chains are much stretched in the region dominated by the short chains (lower layer) but possess relaxed conformations in the upper layer. By slightly changing the solvent selectivity to prefer the short chains, we find a reversion of the layered structure. Such a sensitive switch of film property implies its potential application as tuning the wettability and adhesion of the surface in industry. [ABSTRACT FROM AUTHOR]

The depletion interaction between proteins caused by addition of either uncharged or partially charged oligomers was studied using the canonical Monte Carlo simulation technique and the integral equation theory. A protein molecule was modeled in two different ways: either as (i) a hard sphere of diameter 30.0 Å with net charge 0, or +5, or (ii) as a hard sphere with discrete charges (depending on the pH of solution) of diameter 45.4 Å. The oligomers were pictured as tangentially jointed, uncharged, or partially charged, hard spheres. The ions of a simple electrolyte present in solution were represented by charged hard spheres distributed in the dielectric continuum. In this study we were particularly interested in changes of the protein-protein pair-distribution function, caused by addition of the oligomer component. In agreement with previous studies we found that addition of a nonadsorbing oligomer reduces the phase stability of solution, which is reflected in the shape of the protein-protein pair-distribution function. The value of this function in protein-protein contact increases with increasing oligomer concentration, and is larger for charged oligomers. The range of the depletion interaction and its strength also depend on the length (number of monomer units) of the oligomer chain. The integral equation theory, based on the Wertheim Ornstein-Zernike approach applied in this study, was found to be in fair agreement with Monte Carlo results only for very short oligomers. The computer simulations for a model mimicking the lysozyme molecule (ii) are in qualitative agreement with small-angle neutron experiments for lysozyme-dextran mixtures. [ABSTRACT FROM AUTHOR]

The authors measure 2H NMR multitime correlation functions to investigate the segmental motion of poly(propylene glycol) containing various amounts of the salt LiClO4. 2H NMR two-time correlation functions indicate that addition of salt affects not only the time scale of the segmental motion, but also the degree of the nonexponential relaxation behavior. To quantify the origin of the nonexponential segmental motion, the authors analyze 2H NMR three-time correlation functions. In general, nonexponential relaxation can result from homogeneous dynamics, i.e., intrinsic nonexponentiality, and from heterogeneous dynamics, i.e., existence of a distribution of correlation times G(ln τ). For the studied high and low salt concentrations, including neat poly(propylene glycol), the analysis shows that both homogeneous and heterogeneous contributions are important. 2H NMR four-time correlation functions allow the authors to measure the lifetime of the dynamical heterogeneities. For the studied salt concentrations, the rate exchange occurs on the same time scale as the segmental motion, indicating short-lived dynamical heterogeneities. To arrive at these results, the authors reconsider the interpretation of 2H NMR three-time correlation functions. Results of analytical calculations and computer simulations show that it is necessary to extend the previous way of analysis so as to include effects due to correlated back-and-forth jumps. [ABSTRACT FROM AUTHOR]

In this work we investigate how a pattern imposed in a copolymer film at a certain distance from the surface propagates through the film onto an adsorbing heterogeneous surface. We bias the copolymer film to adopt a specified target pattern and then use simulation to design a surface pattern that helps the adsorbed film to maintain that target pattern. We examine the effect of varying the copolymer chain length, the size of the target pattern, and the distance from the surface where the target pattern is applied, z′, on the extent of pattern transfer. For each chain length, target pattern, and z′ we compare the energy of the system when a pattern is applied in the bulk to the energy when no pattern is applied in order to understand why a certain pattern size is transferred to the surface with higher fidelity than the others. At constant chain length, pattern transfer is best when the pattern size brings the energy of the system close to the energy when no pattern is applied. At constant pattern size, pattern transfer is best in the systems with longer chains. This is because longer chains are more likely to adsorb as brushes and loops which then helps transfer the pattern through the adsorbed film down to the surface. [ABSTRACT FROM AUTHOR]

Large-scale computer simulations show that two orthogonal external fields can control the orientation of lamellar microdomains in diblock copolymers in three dimensions and lead to an enhanced long-range ordering. [ABSTRACT FROM AUTHOR]

There is a rapidly growing interest in performing computer simulations in a closed space, avoiding periodic boundary conditions. To extend the range of potential systems to include also macromolecules, we describe an algorithm for computer simulations of polymer chain molecules on S3, a hypersphere in four dimensions. In particular, we show how to generate initial conformations with a bond angle distribution given by the persistence length of the chain and how to calculate the bending forces for a molecule moving on S3. Furthermore, we discuss how to describe the shape of a macromolecule on S3, by deriving the radius of gyration tensor in this non-Euclidean space. The results from both Monte Carlo and Brownian dynamics simulations in the infinite dilution limit show that the results on S3 and in R3 coincide, both with respect to the size and shape as well as for the diffusion coefficient. All data on S3 can also be described by master curves by suitable scaling by the corresponding values in R3. We thus show how to extend the use of spherical boundary conditions, which are most effective for calculating electrostatic forces, to polymer chain molecules, making it possible to perform simulations on S3 also for polyelectrolyte systems. [ABSTRACT FROM AUTHOR]

Fluorescence lifetime measurements in a polymer chain are modeled using a memory function expansion, computer simulations, and simple scaling arguments. Unless the quenching rate is localized and infinitely fast, the fluorescence lifetime is generally not equivalent to the first passage time. The fluorescence lifetime distribution is decomposed into memory functions that can be measured separately in single-molecule experiments. The leading order of the expansion gives the Wilemski–Fixman (WF) approximation, and the convergence of higher order terms determines its validity. Simulations of the fluorescence quenching on a Rouse chain verify the accuracy of the WF approximation at small contact radii, short contour lengths, and small quenching rates. Detailed investigation of the average fluorescence lifetime reveals two competing mechanisms: the independent motion of end-to-end vector, which dominates at small contact radius, and the slowest relaxation of polymer, which dominates at large contact radius. The Wilemski–Fixman rate is used in combination with scaling arguments to predict the dependence of fluorescence lifetime on the contour length. Our predictions for the scaling of the average lifetime with the contour length are in good agreement with both simulations and recent experiments by Eaton and his group [L. J. Lapidus, W. A. Eaton, and J. Hofrichter, Proc. Natl. Acad. Sci. U.S.A. 97, 7220 (2000)]. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The mean first passage time for the close approach of the ends of a polymer is investigated by computer simulation of the Brownian dynamics of a Rouse model. There are N bonds, of length b, and the contact distance is a. The results for N=50, 75, and 100 are compared with theoretical predictions by Wilemski and Fixman (WF), and by Szabo, Schulten, and Schulten (SSS). Theory differs from simulation by a factor of 1.3 to 3 when a is comparable to b but improves considerably when a is reduced. Closest agreement is obtained at a/b=0.1 for WF with a δ function sink, where theory is within the statistical error (≊5%) of the simulation for each case. [ABSTRACT FROM AUTHOR]