Elastoviscous properties of polyisobutylene. IV. Relaxation time spectrum and calculation of bulk viscosity

The elastoviscous behavior of polyisobutylene may be interpreted in terms of a mechanical model consisting of a distribution of Maxwell elements connected in parallel. The structure of this “generalized Maxwell model” is specified by the distribution of relaxation times of the component elements. The relaxation of stress curve of the material is directly related to the distribution of relaxation times, and general expressions for the bulk viscosities (tensile and shear) of such a system in terms of the distribution of relaxation times are readily obtained. A simple “box distribution” of relaxation times is described which can be used to approximate the relaxation behavior of polyisobutylene at the long-time end of the relaxation time spectrum, and in terms of which the expressions for bulk viscosity reduce to very simple form. The parameters specifying this distribution may be determined from experimental relaxation curves by a simple graphical method. Values of these parameters as a functions of molecular weight and temperature are computed, by use of these data. It is shown that bulk viscosity values calculated from relaxation data by this method are in good agreement with experimental values for both tensile and shear deformations, and for both unfractionated and fractionated polymers. Measurements of viscosity and of relaxation of stress can thus be directly correlated, and could be used in combination to characterize elastoviscous properties over wide ranges of molecular weight and temperature.     

An accelerated algorithm for discrete stochastic simulation of reaction-diffusion systems using gradient-based diffusion and tau-leaping

Stochastic simulation of reaction-diffusion systems enables the investigation of stochastic events arising from the small numbers and heterogeneous distribution of molecular species in biological cells. Stochastic variations in intracellular microdomains and in diffusional gradients play a significant part in the spatiotemporal activity and behavior of cells. Although an exact stochastic simulation that simulates every individual reaction and diffusion event gives a most accurate trajectory of the system's state over time, it can be too slow for many practical applications. We present an accelerated algorithm for discrete stochastic simulation of reaction-diffusion systems designed to improve the speed of simulation by reducing the number of time-steps required to complete a simulation run. This method is unique in that it employs two strategies that have not been incorporated in existing spatial stochastic simulation algorithms. First, diffusive transfers between neighboring subvolumes are based on concentration gradients. This treatment necessitates sampling of only the net or observed diffusion events from higher to lower concentration gradients rather than sampling all diffusion events regardless of local concentration gradients. Second, we extend the non-negative Poisson tau-leaping method that was originally developed for speeding up nonspatial or homogeneous stochastic simulation algorithms. This method calculates each leap time in a unified step for both reaction and diffusion processes while satisfying the leap condition that the propensities do not change appreciably during the leap and ensuring that leaping does not cause molecular populations to become negative. Numerical results are presented that illustrate the improvement in simulation speed achieved by incorporating these two new strategies.     

Understanding diffusion in confined systems: methane in a ZK4 molecular sieve. A molecular dynamics simulation study

The equilibrium probability distribution of N methane molecules adsorbed in the interior of n alpha cages of the ZK4 zeolite, the all-silica analogue of zeolite A, is modeled by a modified hypergeometric distribution where the effects of mutual exclusion between particles are extracted from long molecular dynamics simulations. The trajectories are then analyzed in terms of time-correlation functions for the fluctuations in the occupation number of the alpha cages. The analysis digs out the correlations induced by the spatial distribution of the adsorbed molecules coupled with a migration mechanism where a molecule can pass from one alpha cage to another, one-by-one. These correlations lead to cooperative motion, which manifests itself as a nonexponential decay of the correlators. Our results suggest ways of developing improved lattice approaches that may be useful for studying diffusion in much larger systems and for a much longer observation time.     

Relaxation characteristics and intrinsic birefringence and viscosity of polystyrene solutions for a wide range of molecular weights

A comparison is made between measurements on polystyrene solutions and the relaxation characteristics and intrinsic birefringence and viscosity given by the theory for the flexible Gaussian chain of variable number of segments and with internal viscosity and internal hydrodynamic interaction. This is done in order to determine the applicability of the theory to polymers over a wide range of molecular weights, including the low molecular weight range in which there may be conflict with the theoretical assumption of chains having a large number of segments. The longest, terminal relaxation time and the number of chain segments are determined from measurements of the frequency dependence of oscillatory flow birefringence while the intrinsic birefringence and viscosity are determined from steady flow measurements. The range of molecular weights studied is from approximately 900 at 10⁶. It is found that the segment weight is approximately 1000 and the number of segments is in direct proportion to the molecular weight for the range from 1 to 1000 segments. The terminal relaxation time has a molecular weight dependence of the type given by the theory but with better agreement for higher molecular weights. While the measured dependences of the intrinsic viscosity and birefringence are in agreement with theory for molecular weights greater than 5 × 10⁴, they deviate significantly for molecular weights below 1 × 10⁴. The ratio of the intrinsic birefringence to intrinsic viscosity, which in theory is a constant independent of molecular weight, is found to change at the lower molecular weights.     

Prediction of viscosity for molecular fluids at experimentally accessible shear rates using the transient time correlation function formalism

Nonequilibrium molecular dynamics (NEMD) simulations were performed and the transient time correlation function (TTCF) method applied to calculate the shear viscosity of n-decane. Using the TTCF method we were able to calculate the viscosity at shear rate orders of magnitude lower than is possible by direct NEMD simulation alone. For the first time for a molecular fluid, we were able to simulate shear rates accessible by experimental measurements, which are typically performed at shear rates well below those accessible by NEMD simulation. The TTCF method allows us to close the gap between the lowest shear rates accessible by MD simulation and the highest shear rates possible in experimental studies. Additionally a multiple time step method for Gaussian thermostatted SLLOD equations of motion was developed following earlier work [G. A. Pan et al., J. Chem. Phys. 122, 4114 (2005)] for atomic fluids.     

Relaxation time measurements and molecular dynamics of dipalmitoyl-α-lecithin in the solid phase

The relaxation times T₁ and T_D were measured for dipalmitoyl-L-α-lecithin and dipalmitoyl-DL-α-lecithin at various temperatures. It was possible to separate the contribution of the slow motions to the relaxation rates from the fast ones. For some types of motions the activation energies could be estimated.     

Structure,thermodynamics and diffusion in asymmetric binary mixtures: a molecular dynamics simulation study

Structural and thermodynamic properties as well as diffusion coefficients of binary fluid mixtures with asymmetry in mass, size, charge and their combinations have been studied using classical molecular dynamics simulations. The fluid mixture is modelled as spherical particles interacting via the Weeks–Chandler–Andersen and Coulomb potential. The diameter, charge and mass of the fluid particles are in the range 6–60 Å, 1–10e and 1—500 amu, respectively. Systematic variations in pair-correlation functions, thermodynamic properties as well as the self-diffusion coefficient are found with the size, charge and mass ratio of the particles. The self-diffusion coefficient for systems having more than one type of asymmetry is calculated and expressed in terms of diffusion coefficients of systems with only one type of asymmetry.     

Diffusion and viscosity of a calamitic liquid crystal model studied by computer simulation

We report a molecular dynamics simulation study on an ensemble of rod-like particles, each composed of nine soft spheres held rigidly along a line. We have calculated translational mean square displacements and velocity autocorrelation functions in the fluid phases exhibited by the model, i.e., smectic A, nematic and isotropic. These quantities have then been used to compute diffusion coefficients. In addition, we have calculated viscosities in the nematic and isotropic phases. Despite its crude nature, the model is capable of providing a faithful reproduction of many features of the transport behavior observed in real liquid-crystalline materials. The simulation results have been compared with the predictions of the modified affine transformation theory, finding only a fair agreement.     

A systematic molecular simulation study of ionic liquid surfaces using intrinsic analysis methods

In this paper, we apply novel intrinsic analysis methods, coupled with bivariate orientation analysis, to obtain a detailed picture of the molecular-level structure of ionic liquid surfaces. We observe pronounced layering at the interface, alternating non-polar with ionic regions. The outermost regions of the surface are populated by alkyl chains, which are followed by a dense and tightly packed layer formed of oppositely charged ionic moieties. We then systematically change the cation chain length, the anion size, the temperature and the molecular model, to examine the effect of each of these parameters on the interfacial structure. Increasing the cation chain length promotes orientations in which the chain is pointing into the vapor, thus increasing the coverage of the surface with alkyl groups. Larger anions promote a disruption of the dense ionic layer, increasing the orientational freedom of cations and increasing the amount of free space. The temperature had a relatively small effect on the surface structure, while the effect of the choice of molecular model was clearly significant, particularly on the orientational preferences at the interface. Our study demonstrates the usefulness of molecular simulation methods in the design of ionic liquids to suit particular applications.     

Studies of translational diffusion in the smectic A phase of a Gay-Berne mesogen using molecular dynamics computer simulation

Molecular dynamics computer simulations are used to determine the self-diffusion coefficients for a Gay-Berne model mesogen GB (4.4,20,1,1) in the isotropic, nematic and smectic A phases along two isobars. The values of the parallel and perpendicular diffusion coefficients, D(parallel) and D(perpendicular), are calculated and compared in the different phases. For the phase sequence isotropic-smectic A, D(perpendicular)*> or =D(parallel)* over the whole smectic A range with the ratio D(parallel)*/D(perpendicular)* decreasing with decreasing temperature. At a higher pressure, a nematic phase is observed between these two phases and we find that D(parallel)*>D(perpendicular)* throughout the nematic region and the inequality D(parallel)*>D(perpendicular)* remains on entering the smectic A phase. However, the ratio D(parallel)*/D(perpendicular)* decreases with decreasing temperature within the smectic A range and eventually this ratio inverts such that D(perpendicular)*>D(parallel)* at low temperatures. The temperature dependence of the parallel diffusion coefficient in the smectic A phase for this model mesogen is compared to that predicted by a theoretical model for diffusion subject to a cosine potential.     

Finite time thermodynamics: limiting performance of diffusion engines and membrane systems

In this paper, the limiting performance of membrane systems with inhomogeneous composition is studied within the class of fixed rate processes. The problem of maintaining a nonequilibrium state in such a system using minimal power (separation problem) and the problem of extracting maximal power from such a system (diffusion engine problem) are formulated and solved. Results are obtained for diffusion engines with constant and periodic contact between the working body and the reservoirs.     

Prediction of CO2 and H2 solubility,diffusion, and permeability in MFI zeolite by molecular dynamics simulation

Structural Chemistry - Molecular dynamics simulation has been employed to calculate the amounts of solubility, diffusion coefficient, and permeability for the pure and volumetric binary mixture of...     

Solvent and hydrogen confinement in molecular capsules-Hirshfeld surface and molecular simulation analysis

Hirshfeld surface analysis of the 'ordered' inner phase of the molecular capsule complex, [(chloroform)(6)@C-n-butylpyrogallol[4]arene)(6)], provides insight into the intermolecular contacts and orientation of the solvent molecules. Molecular simulations show that adding two or three hydrogen molecules to the six solvent molecules is energetically favoured, and this correlates with NMR studies.     

Electrolysis of water in the diffusion layer: first-principles molecular dynamics simulation

With Car-Parrinello molecular dynamics simulations the elementary reaction steps of the electrolysis of bulk water are investigated. To simulate the reactions occurring near the anode and near the cathode, electrons are removed or added, respectively. The study focuses on the reactions in pure water. Effects depending on a particular electrode surface or a particular electrolyte are ignored. Under anodic conditions, the reaction continues till molecular oxygen is formed, under cathodic conditions the formation of molecular hydrogen is observed. In addition the formation of hydrogen peroxide is observed as an intermediate of the anodic reaction. The simulations demonstrate that the electrochemistry of oxygen formation without direct electrode contact can be explained by radical reactions in a solvent. These reactions may involve the intermediate formation of ions. The hydrogen formation is governed by rapid proton transfers between water molecules.     

Structure and diffusion in amorphous aluminum silicate: a molecular dynamics computer simulation

The amorphous aluminum silicate (Al2O3)2(SiO2) [AS2] is investigated by means of large scale molecular dynamics computer simulations. We consider fully equilibrated melts in the temperature range 6100 K> or =T> or =2300 K as well as glass configurations that were obtained from cooling runs from T=2300 to 300 K with a cooling rate of about 10(12) K/s. Already at temperatures as high as 4000 K, most of the Al and Si atoms are fourfold coordinated by oxygen atoms. Thus, the structure of AS2 is that of a disordered tetrahedral network. The packing of AlO4 tetrahedra is very different from that of SiO4 tetrahedra in that Al is involved with a relatively high probability in small-membered rings and in triclusters in which an O atom is surrounded by four cations. We find as typical configurations two-membered rings with two Al atoms in which the shared O atoms form a tricluster. On larger length scales, the system shows a microphase separation in which the Al-rich network structure percolates through the SiO2 network. The latter structure gives rise to a prepeak in the static structure factor at a wave number q=0.5 A(-1). A comparison of experimental x-ray data with the results from the simulation shows good agreement for the structure function. The diffusion dynamics in AS2 is found to be much faster than in SiO2. We show that the self-diffusion constants for O and Al are very similar and that they are by a factor of 2-3 larger than the one for Si.     

Relaxation of chemical potential and a generalized diffusion equation

The effect of slow structural relaxation in a solvent of high viscosity on the chemical potential driving the diffusion of penetrant molecules is described by a generalized diffusion equation with a memory term. The linearized version of this equation is solved for some special cases, and the correlation function of concentration fluctuations in thermodynamic equilibrium is calculated. As a result of the memory term, for very slow relaxation two different stages of the diffusion process can be distinguished.     

Relation between amorphous structure of polymers and penetrant diffusion: A molecular dynamics simulation

Molecular dynamics simulation has been performed for studying the relation between amorphous structure of polymers and penetrant diffusion. The self-diffusion coefficients of O₂ and He in various polymer models, which differ from each other in view of the amorphous structure, were calculated above their glass transition temperatures. The amorphous structure was characterized by considering the percolation of the unoccupied volume. A good correlation was found between the self-diffusion coefficients and the number of clusters in the unoccupied volume at the critical point of the percolation. Based on the simulated cluster size distribution at the critical point, we defined a parameter into which effects of both the amorphous structure and the penetrant size are well incorporated. It was confirmed that the penetrant diffusion is intimately associated with the amorphous structure of polymers.     

Study on laminar viscosity and zero shear viscosity of latex systems

The flowing nature and rheological properties of polymethyl methacrylate latex systems in a coaxial cylinder viscometer were studied on the basis of laminar shear flow model and rheological experimental data. The physical meaning of laminar viscosity (eta(i,j)) and zero shear viscosity (eta(0)) were described. We assumed that laminar shear flows depended on position and shear time, so microrheological parameters were the function of position and shear time. eta(i,j) was the viscosity of any shear sheet i between two neighboring laminar shear flows at time t; j was denoted as j=t/Deltat; and Deltat was the interacting time of two particles or two laminar shear flows. tau(i,j) and gamma(i,j) were shear stress and shear rate of any shear sheet i at j moment. According to Newton regulation tau(i,j)=eta(i,j)gamma(i,j), apparent viscosity eta(a) should be a statistically mean value of j shear sheets laminar viscosity at j moment, i.e., eta(a)= summation operator(i=j)eta(i,j)gamma(i,j)/ summation operator(i=j)gamma(i,j). eta(0) was defined as shear viscosity between a laminar shear flow and a still fluid surface, i.e., eta(0)=(tau(i,j)/gamma(i,j))(j-i-->0). These new ideas described above may be helpful in the study of the micromechanisms of latex particle systems and worthy of more research.     

Comparison of DNA hydration patterns obtained using two distinct computational methods, molecular dynamics simulation and three-dimensional reference interaction site model theory

Because proteins and DNA interact with each other and with various small molecules in the presence of water molecules, we cannot ignore their hydration when discussing their structural and energetic properties. Although high-resolution crystal structure analyses have given us a view of tightly bound water molecules on their surface, the structural data are still insufficient to capture the detailed configurations of water molecules around the surface of these biomolecules. Thanks to the invention of various computational algorithms, computer simulations can now provide an atomic view of hydration. Here, we describe the apparent patterns of DNA hydration calculated by using two different computational methods: Molecular dynamics (MD) simulation and three-dimensional reference interaction site model (3D-RISM) theory. Both methods are promising for obtaining hydration properties, but until now there have been no thorough comparisons of the calculated three-dimensional distributions of hydrating water. This rigorous comparison showed that MD and 3D-RISM provide essentially similar hydration patterns when there is sufficient sampling time for MD and a sufficient number of conformations to describe molecular flexibility for 3D-RISM. This suggests that these two computational methods can be used to complement one another when evaluating the reliability of the calculated hydration patterns.