No Arabic abstract
Obtaining a rigorous and reliable method for linking computer simulations of polymer blends and composites at different length scales of interest is a highly desirable goal in soft matter physics. In this paper a multiscale modeling procedure is presented for the efficient calculation of the static structural properties of binary homopolymer blends. The procedure combines computer simulations of polymer chains on two different length scales, using a united atom representation for the finer structure and a highly coarse-grained approach on the meso-scale, where chains are represented as soft colloidal particles interacting through an effective potential. A method for combining the structural information by inverse mapping is discussed, allowing for the efficient calculation of partial correlation functions, which are compared with results from full united atom simulations. The structure of several polymer mixtures is obtained in an efficient manner for several mixtures in the homogeneous region of the phase diagram. The method is then extended to incorporate thermal fluctuations through an effective chi parameter. Since the approach is analytical, it is fully transferable to numerous systems.
We study the thermodynamic stability of fluid-fluid phase separation in binary nonadditive mixtures of hard-spheres for moderate size ratios. We are interested in elucidating the role played by small amounts of nonadditivity in determining the stability of fluid-fluid phase separation with respect to the fluid-solid phase transition. The demixing curves are built in the framework of the modified-hypernetted chain and of the Rogers-Young integral equation theories through the calculation of the Gibbs free energy. We also evaluate fluid-fluid phase equilibria within a first-order thermodynamic perturbation theory applied to an effective one-component potential obtained by integrating out the degrees of freedom of the small spheres. A qualitative agreement emerges between the two different approaches. We also address the determination of the freezing line by applying the first-order thermodynamic perturbation theory to the effective interaction between large spheres. Our results suggest that for intermediate size ratios a modest amount of nonadditivity, smaller than earlier thought, can be sufficient to drive the fluid-fluid critical point into the thermodinamically stable region of the phase diagram. These findings could be significant for rare-gas mixtures in extreme pressure and temperature conditions, where nonadditivity is expected to be rather small.
We compute the fourth virial coefficient of a binary nonadditive hard-sphere mixture over a wide range of deviations from diameter additivity and size ratios. Hinging on this knowledge, we build up a $y$ expansion [B. Barboy and W. N. Gelbart, J. Chem. Phys. {bf 71}, 3053 (1979)] in order to trace the fluid-fluid coexistence lines which we then compare with the available Gibbs-ensemble Monte Carlo data and with the estimates obtained through two refined integral-equation theories of the fluid state. We find that in a regime of moderately negative nonadditivity and largely asymmetric diameters, relevant to the modelling of sterically and electrostatically stabilized colloidal mixtures, the fluid-fluid critical point is unstable with respect to crystallization.
Mesoscopic molecular dynamics simulations are used to determine the large scale structure of several binary polymer mixtures of various chemical architecture, concentration, and thermodynamic conditions. By implementing an analytical formalism, which is based on the solution to the Ornstein-Zernike equation, each polymer chain is mapped onto the level of a single soft colloid. From the appropriate closure relation, the effective, soft-core potential between coarse-grained units is obtained and used as input to our mesoscale simulations. The potential derived in this manner is analytical and explicitly parameter dependent, making it general and transferable to numerous systems of interest. From computer simulations performed under various thermodynamic conditions the structure of the polymer mixture, through pair correlation functions, is determined over the entire miscible region of the phase diagram. In the athermal regime mesoscale simulations exhibit quantitative agreement with united atom simulations. Furthermore, they also provide information at larger scales than can be attained by united atom simulations and in the thermal regime approaching the phase transition.
We extensively investigated the critical behavior of mixtures of colloids and polymers via the two-component Asakura-Oosawa model and its reduction to a one-component colloidal fluid using accurate theoretical and simulation techniques. In particular the theoretical approach, hierarchical reference theory [Adv. Phys. 44, 211 (1995)], incorporates realistically the effects of long-range fluctuations on phase separation giving exponents which differ strongly from their mean-field values, and are in good agreement with those of the three-dimensional Ising model. Computer simulations combined with finite-size scaling analysis confirm the Ising universality and the accuracy of the theory, although some discrepancy in the location of the critical point between one-component and full-mixture description remains. To assess the limit of the pair-interaction description, we compare one-component and two-component results.
There are many proteins or protein complexes which have multiple DNA binding domains. This allows them to bind to multiple points on a DNA molecule (or chromatin fibre) at the same time. There are also many proteins which have been found to be able to compact DNA in vitro, and many others have been observed in foci or puncta when fluorescently labelled and imaged in vivo. In this work we study, using coarse-grained Langevin dynamics simulations, the compaction of polymers by simple model proteins and a phenomenon known as the bridging-induced attraction. The latter is a mechanism observed in previous simulations [Brackley et al., Proc. Natl. Acad. Sci. USA 110 (2013)], where proteins modelled as spheres form clusters via their multivalent interactions with a polymer, even in the absence of any explicit protein-protein attractive interactions. Here we extend this concept to consider more detailed model proteins, represented as simple patchy particles interacting with a semi-flexible bead-and-spring polymer. We find that both the compacting ability and the effect of the bridging-induced attraction depend on the valence of the model proteins. These effects also depend on the shape of the protein, which determines its ability to form bridges.