Fluids density functional theory and initializing molecular dynamics simulations of block copolymers

Phase Behavior of Gradient Copolymer Melts with Different Gradient Strengths Revealed by Mesoscale Simulations

Design and preparation of functional nanomaterials with specific properties requires precise control over their microscopic structure. A prototypical example is the self-assembly of diblock copolymers, which generate highly ordered structures controlled by three parameters: the chemical incompatibility between blocks, block size ratio and chain length. Recent advances in polymer synthesis have allowed for the preparation of gradient copolymers with controlled sequence chemistry, thus providing additional parameters to tailor their assembly. These are polydisperse monomer sequence, block size distribution and gradient strength. Here, we employ dissipative particle dynamics to describe the self-assembly of gradient copolymer melts with strong, intermediate, and weak gradient strength and compare their phase behavior to that of corresponding diblock copolymers. Gradient melts behave similarly when copolymers with a strong gradient are considered. Decreasing the gradient strength leads to the widening of the gyroid phase window, at the expense of cylindrical domains, and a remarkable extension of the lamellar phase. Finally, we show that weak gradient strength enhances chain packing in gyroid structures much more than in lamellar and cylindrical morphologies. Importantly, this work also provides a link between gradient copolymers morphology and parameters such as chemical incompatibility, chain length and monomer sequence as support for the rational design of these nanomaterials.

Recent progress in the science of complex coacervation

Complex coacervation is an associative, liquid-liquid phase separation that can occur in solutions of oppositely-charged macromolecular species, such as proteins, polymers, and colloids. This process results in a coacervate phase, which is a dense mix of the oppositely-charged components, and a supernatant phase, which is primarily devoid of these same species. First observed almost a century ago, coacervates have since found relevance in a wide range of applications; they are used in personal care and food products, cutting edge biotechnology, and as a motif for materials design and self-assembly. There has recently been a renaissance in our understanding of this important class of material phenomena, bringing the science of coacervation to the forefront of polymer and colloid science, biophysics, and industrial materials design. In this review, we describe the emergence of a number of these new research directions, specifically in the context of polymer-polymer complex coacervates, which are inspired by a number of key physical and chemical insights and driven by a diverse range of experimental, theoretical, and computational approaches.

Tailoring the mechanical properties of polymer nanocomposites via interfacial engineering

The improvement of mechanical properties of polymer nanocomposites (PNCs) has been studied for many years, with the main focus on the structure of the nanofillers. Much less effort has been devoted to unraveling the factors controlling the structure of the grafted chains. Herein, through coarse-grained molecular-dynamics simulations, we have successfully fabricated an ideal, mechanically-interlocked composite structure composed of end-functionalized chains grafted to the nanoparticle surface forming rings and making the matrix chains thread through these rings. Depending on the details of the grafting, the reinforcement effect can be remarkable, improving the tensile stress of the system significantly up to 700%. Meanwhile, anisotropy of the system's mechanical response is also observed. Furthermore, the influence of the grafted chain distribution on the mechanical properties of the system has been investigated as well. We observe that the mechanical properties of the system are closely related to the total number of the beads in the grafted chains or the synergistic effect between the length and density of the grafted chains leads to no significant difference in the performance of systems. At constant grafting density, the mechanical properties of the systems correlate negatively to the grafted chain length. In general, our study should help to design and fabricate high-performance PNCs with excellent mechanical properties.

Ion Correlation Effects in Salt-Doped Block Copolymers

We apply classical density functional theory to study how salt changes the microphase morphology of diblock copolymers. Polymers are freely jointed and one monomer type favorably interacts with ions, to account for the selective solvation that arises from different dielectric constants of the microphases. By including correlations from liquid state theory of an unbound reference fluid, the theory can treat chain behavior, microphase separation, ion correlations, and preferential solvation, at the same coarse-grained level. We show good agreement with molecular dynamics simulations.

Effect of sequence dispersity on morphology of tapered diblock copolymers from molecular dynamics simulations

Tapered diblock copolymers are similar to typical AB diblock copolymers but have an added transition region between the two blocks which changes gradually in composition from pure A to pure B. This tapered region can be varied from 0% (true diblock) to 100% (gradient copolymer) of the polymer length, and this allows some control over the microphase separated domain spacing and other material properties. We perform molecular dynamics simulations of linearly tapered block copolymers with tapers of various lengths, initialized from fluids density functional theory predictions. To investigate the effect of sequence dispersity, we compare systems composed of identical polymers, whose taper has a fixed sequence that most closely approximates a linear gradient, with sequentially disperse polymers, whose sequences are created statistically to yield the appropriate ensemble average linear gradient. Especially at high segregation strength, we find clear differences in polymer conformations and microstructures between these systems. Importantly, the statistical polymers are able to find more favorable conformations given their sequence, for instance, a statistical polymer with a larger fraction of A than the median will tend towards the A lamellae. The conformations of the statistically different polymers can thus be less stretched, and these systems have higher overall density. Consequently, the lamellae formed by statistical polymers have smaller domain spacing with sharper interfaces.

Diffusion of Selective Penetrants in Interfacially Modified Block Copolymers from Molecular Dynamics Simulations

To show the influence of the interface on structure and dynamics of microphase separated polymer systems, we study interfacially modified AB block copolymers with small molecule penetrants. The polymers have a random midblock or tapered midblock whose composition varies from pure A to pure B (or from pure B to pure A for an inverse taper) between two pure blocks of A and B. We perform simple coarse-grained molecular dynamics simulations of symmetric polymers that form lamellae. With normal tapering, both polymer and penetrant diffusion parallel to the lamellae increases as taper length increases. Inverse tapered polymers exist in different conformational states (e.g., stretched vs folded back and forth across the interface) with different dynamic behavior, leading to nonmonotonic trends in their diffusion. However, the local mixing of monomers (rather than polymer conformation) appears to be the most important factor in determining penetrant diffusion.

A Review of Multiscale Computational Methods in Polymeric Materials

Polymeric materials display distinguished characteristics which stem from the interplay of phenomena at various length and time scales. Further development of polymer systems critically relies on a comprehensive understanding of the fundamentals of their hierarchical structure and behaviors. As such, the inherent multiscale nature of polymer systems is only reflected by a multiscale analysis which accounts for all important mechanisms. Since multiscale modelling is a rapidly growing multidisciplinary field, the emerging possibilities and challenges can be of a truly diverse nature. The present review attempts to provide a rather comprehensive overview of the recent developments in the field of multiscale modelling and simulation of polymeric materials. In order to understand the characteristics of the building blocks of multiscale methods, first a brief review of some significant computational methods at individual length and time scales is provided. These methods cover quantum mechanical scale, atomistic domain (Monte Carlo and molecular dynamics), mesoscopic scale (Brownian dynamics, dissipative particle dynamics, and lattice Boltzmann method), and finally macroscopic realm (finite element and volume methods). Afterwards, different prescriptions to envelope these methods in a multiscale strategy are discussed in details. Sequential, concurrent, and adaptive resolution schemes are presented along with the latest updates and ongoing challenges in research. In sequential methods, various systematic coarse-graining and backmapping approaches are addressed. For the concurrent strategy, we aimed to introduce the fundamentals and significant methods including the handshaking concept, energy-based, and force-based coupling approaches. Although such methods are very popular in metals and carbon nanomaterials, their use in polymeric materials is still limited. We have illustrated their applications in polymer science by several examples hoping for raising attention towards the existing possibilities. The relatively new adaptive resolution schemes are then covered including their advantages and shortcomings. Finally, some novel ideas in order to extend the reaches of atomistic techniques are reviewed. We conclude the review by outlining the existing challenges and possibilities for future research.