The optimized Rouse-Zimm theory of excluded volume effects on chain dynamics

Based on the optimized Rouse-Zimm (ORZ) approximation to the Kirkwood diffusion equation, we investigate the effects of excluded volume interactions on the single chain dynamics. By incorporating the nonuniformly expanded moments of interbead distances into the expressions for the diffusion and structure matrices appearing in the ORZ diffusion equation, we obtain the general relaxation spectrum for flexible chains that is valid over the range from theta; solvents to good solvents. The present theory involves four parameters: the Kuhn statistical length b(0), the bead number N, the excluded volume parameter z, and the hydrodynamic interaction parameter h(*). These model parameters are determined from structural data of polymers with the aid of the quasi-two-parameter theory. The set of relaxation times of ORZ normal modes calculated with these bead-and-spring model parameters enables the theoretical prediction of various frictional and dynamical properties of polymers within a unified framework. The present ORZ theory generalizes the Ptitsyn-Eizner-type approaches by incorporating the nonuniform chain expansion effect into the structure matrix as well as the diffusion matrix.     

Response of the mean-square optical anistropy of chain molecules to the imposition of excluded volume

The influence of the chain expansion produced by excluded volume on the mean-square optical anisotropy has been studied in six types of polymers. The mean-square optical anistropy for a specified configuration is calculated using the valence optical scheme. Realistic rotational isomeric state models are used for the configurational statistics of the unperturbed chains. Excluded volume is introduced by hard sphere interactions. Results obtained with chains of 100, 200, 300, and 400 bonds permit extrapolation to the behavior expected for much longer chains. The mean-square optical anisotropy of polyethylene is insensitive to excluded volume. A similar conclusion was obtained several years ago in a study of chains confined to a tetrahedral lattice and weighted in a manner appropriate for the short-range interactions in polyethylene.² Different behavior is seen in poly(vinyl chloride), poly(vinyl bromide), polystyrene, poly(p-chlorostyrene), and poly(p-bromostyrene). The mean-square optical anisotropy of these five vinyl polymers is sensitive to the imposition of excluded volume if the stereochemical composition is exclusively racemic. Much smaller effects are seen in meso chains and in chains with Bernoullian statistics and an equal probability for meso and racemic dyads.     

Concentration dependence of excluded volume effects

The concentration dependence of the excluded volume effects in polymer solutions is investigated. Through thermodynamic arguments for the interpenetration of polymer segments and the free energy change, we show that the disappearance of the excluded volume effects should occur at medium concentration. The result is in accord with the recent experimental observations.     

Influence of excluded volume on the conformational property of tail-like chain on the simple cubic lattice

The influence of excluded volume on the conformational property of linear tail-like chain with one end attached to a flat surface is investigated by means of dynamic Monte Carlo method. Conformational properties such as mean-square end-to-end distance 〈R²〉, mean-square radius of gyration 〈S²〉 and mean asphericity parameter 〈A〉 are calculated for random walking (RW) and self-avoiding walking (SAW) tail-like chains on the simple cubic lattice. We find that the EV has nearly the same effect on 〈R²〉 as on 〈S²〉: (1) 〈R²〉_SAW/〈R²〉_RW≈〈S²〉_SAW/〈S²〉_RW∝n^0.204±0.05, where n is the chain length, and (2) the limiting value of 〈R²〉/〈S²〉≈7.7 for both chains. The distribution P(R) of the SAW tail-like chain can be expressed as a R⁴ correction of that of the RW one. We find that the value 〈A〉 of the SAW tail-like chain is bigger than that of the RW tail-like chain for all chain lengths, and the limiting values are 0.446±0.006 and 0.403±0.005 respectively.     

The importance of excluded solvent volume effects in computing hydration free energies

Continuum dielectric methods such as the Born equation have been widely used to compute the electrostatic component of the solvation free energy, DeltaG(solv)(elec), because they do not need to include solvent molecules explicitly and are thus far less costly compared to molecular simulations. All of these methods can be derived from Gauss Law of Maxwell's equations, which yields an analytical solution for the solvation free energy, DeltaG(Born), when the solute is spherical. However, in Maxwell's equations, the solvent is assumed to be a structureless continuum, whereas in reality, the near-solute solvent molecules are highly structured unlike far-solute bulk solvent. Since we have recently reformulated Gauss Law of Maxwell's equations to incorporate the near-solute solvent structure by considering excluded solvent volume effects, we have used it in this work to derive an analytical solution for the hydration free energy of an ion. In contrast to continuum solvent models, which assume that the normalized induced solvent electric dipole density P(n) is constant, P(n) mimics that observed from simulations. The analytical formula for the ionic hydration free energy shows that the Born radius, which has been used as an adjustable parameter to fit experimental hydration free energies, is no longer ill defined but is related to the radius and polarizability of the water molecule, the hydration number, and the first peak position of the solute-solvent radial distribution function. The resulting DeltaG(solv)(elec) values are shown to be close to the respective experimental numbers.     

Solution of the excluded volume problem for biaxial particles

We report the solution of the excluded volume problem for a pair of biaxial hard molecules; namely, sphero-platelets. As an application of this result we study the isotropic to nematic liquid crystal transition for a fluid composed of these particles in the Onsager limit (length δ breadth or width). We show that the range of stability of the isotropic phase decreases with increasing particle biaxiality.     

An analytical approximation for the orientation-dependent excluded volume of tangent hard sphere chains of arbitrary chain length and flexibility

Onsager-like theories are commonly used to describe the phase behavior of nematic (only orientationally ordered) liquid crystals. A key ingredient in such theories is the orientation-dependent excluded volume of two molecules. Although for hard convex molecular models this is generally known in analytical form, for more realistic molecular models that incorporate intramolecular flexibility, one has to rely on approximations or on computationally expensive Monte Carlo techniques. In this work, we provide a general correlation for the excluded volume of tangent hard-sphere chains of arbitrary chain length and flexibility. The flexibility is introduced by means of the rod-coil model. The resulting correlation is of simple analytical form and accurately covers a wide range of pure component excluded volume data obtained from Monte Carlo simulations of two-chain molecules. The extension to mixtures follows naturally by applying simple combining rules for the parameters involved. The results for mixtures are also in good agreement with data from Monte Carlo simulations. We have expressed the excluded volume as a second order power series in sin?(γ), where γ is the angle between the molecular axes. Such a representation is appealing since the solution of the Onsager Helmholtz energy functional usually involves an expansion of the excluded volume in Legendre coefficients. Both for pure components and mixtures, the correlation reduces to an exact expression in the limit of completely linear chains. The expression for mixtures, as derived in this work, is thereby an exact extension of the pure component result of Williamson and Jackson [Mol. Phys. 86, 819-836 (1995)].     

Solution of the excluded volume problem for biaxial particles

Abstract

We report the solution of the excluded volume problem for a pair of biaxial hard molecules; namely, sphero-platelets. As an application of this result we study the isotropic to nematic liquid crystal transition for a fluid composed of these particles in the Onsager limit (length δ breadth or width). We show that the range of stability of the isotropic phase decreases with increasing particle biaxiality.     

Equilibrium and dynamical properties of gaussian stiff chain molecules

We present a novel analytically tractable model for stiff chain molecules. The equilibrium distribution function of the chain is derived using the maximum-entropy principle. For that purpose, we first formulate a discrete chain model, where the connections of the points and the restriction on bending are taken into account via constraints. We then perform the limit to a continuous chain and show that the mean-square end-to-end distance and the radius of gyration of the continuous chain are identical with the same quantities of the Kratky-Porod wormlike chain. The dynamics of our chain is investigated in dilute solution without hydrodynamic interactions. The linear dynamical equation is solved by a normal mode analysis. We discuss the dependence of the relaxation times on the single parameter of the model, the persistence length. For small persistence lengths we obtain the well known relaxation times of the Rouse model. In the stiff-chain limit, we find the pure bending relaxation times and, in addition, the rotational relaxation time.     

Effects of excluded volume interaction and dimensionality on diffusion-mediated reactions

The kinetic problem of a diffusion-mediated reaction, in which minority reactants are immobile and majority reactants are mobile, is known as the target problem. The standard theory of the target problem ignores the excluded volume interaction between the mobile reactants. Recently, a new theory of the target problem was proposed where the effect of excluded volume interaction was analytically investigated using a lattice model with prohibited double occupancy of the lattice sites. The results of that theory are approximate and need verification. In this work, we perform Monte Carlo simulations on lattices and use their results to assess the accuracy of the analytical theory. We also generalize our theory to the case of different dimensionality and perform calculations for lattices in one- and two-dimensional systems. The analytical results accurately reproduce the simulation results except in the dilute limit in one dimension. For any dimensions, the decay of the target survival probability is accelerated by the presence of excluded volume interaction.     

Algorithm for rapid calculation of excluded volume of large molecules

An algorithm is described for rapid calculation of excluded volume of large molecules. The excluded volume is defined based on coordinates of constituent atoms as the volume of overlapping spheres, each standing for a space around an atom inaccessible for a solvent molecule. A computer program based on the algorithm has been tested on a protein, ovomucoid. The accuracy of the numerical calculation is discussed.     

Kinetic effects of excluded volume and selective adsorption in macromolecule solutions

Kinetic effects of excluded volume and selective adsorption in macromolecule solutions have been studied for an arbitrary position of reacting groups along the chain. It has been shown that the reactivities of similar functional groups can depend upon their location along the chain and also upon the nature of the solvent.     

Macroion adsorption: the crucial role of excluded volume and co-ions

The adsorption of charged colloids (macroions) onto an oppositely charged planar substrate is investigated theoretically. Taking properly into account the finite size of the macroions, unusual behaviors are reported. It is found that the role of the co-ions (the little salt-ions carrying the same sign of charge as that of the substrate) is crucial in understanding the mechanisms involved in the process of macroion adsorption. In particular, the co-ions can accumulate near the substrate's surface and lead to a counterintuitive surface charge amplification.     

Effect of excluded volume interactions on the interfacial properties of colloid-polymer mixtures

We report a numerical study of equilibrium phase diagrams and interfacial properties of bulk and confined colloid-polymer mixtures using grand canonical Monte Carlo simulations. Colloidal particles are treated as hard spheres, while the polymer chains are described as soft repulsive spheres. The polymer-polymer, colloid-polymer, and wall-polymer interactions are described by density-dependent potentials derived by Bolhuis and Louis [Macromolecules 35, 1860 (2002)]. We compared our results with those of the Asakura-Oosawa-Vrij model [J. Chem. Phys. 22, 1255 (1954); J. Polym Sci 33, 183 (1958); Pure Appl. Chem. 48, 471 (1976)] that treats the polymers as ideal particles. We find that the number of polymers needed to drive the demixing transition is larger for the interacting polymers, and that the gas-liquid interfacial tension is smaller. When the system is confined between two parallel hard plates, we find capillary condensation. Compared with the Asakura-Oosawa-Vrij model, we find that the excluded volume interactions between the polymers suppress the capillary condensation. In order to induce capillary condensation, smaller undersaturations and smaller plate separations are needed in comparison with ideal polymers.     

Molecular silverware. I. General solutions to excluded volume constrained problems

General mathematical solutions to excluded volume constrained problems in computational chemistry are reported. The solutions have been used to create a new family of molecular modeling algorithms to facilitate the study of molecular interactions in condensed phases. The new algorithms, collectively known as Molecular Silverware, are for the most part interactive and designed for packing, solvating, and sampling molecules embedded in simple or complex topological environments. Multifolded, disconnected, or porous molecular structures are permitted. Molecular Silverware assists the preparation of Monte Carlo and molecular dynamics simulations at a small fraction of the total simulation time. Primary targets for applications include the study of molecular recognition mechanisms and the selective binding of DNA, RNA, peptides, saccharides and other biopolymers in solution as well as the prediction of phase separation behavior and physical properties of non-crystalline condensed phases such as bulk polymers, polymer blends, organic liquids, membranes, micelles, gels, crosslinked networks, glasses, and amorphous heterogeneous catalysts. As a result of this new approach to excluded volume constraints, the computer simulation of noncrystalline condensed phases is no longer hampered by the lack of a general and efficient method for the creation and configurational sampling of small and large molecular assemblies at high densities.