Open-source code for self-consistent field theory calculations of block polymer phase behavior on graphics processing units

Equilibrium phase behavior of gyroid-forming diblock polymer thin films

Thin-film confinement of self-assembling block polymers results in materials with myriad potential applications-including membranes and optical devices-and provides design parameters for altering phase behavior that are not available in the bulk, namely, film thickness and preferential wetting. However, most research has been limited to lamella- and cylinder-forming polymers; three-dimensional phases, such as double gyroid (DG), have been observed in thin films, but their phase behavior under confinement is not yet well understood. We use self-consistent field theory to predict the equilibrium morphology of bulk-gyroid-forming AB diblock polymers under thin-film confinement. Phase diagrams reveal that the (211) orientation of DG, often observed in experiments, is stable between nonpreferential boundaries at thicknesses as small as 1.2 times the bulk DG lattice parameter. The (001) orientation is stable between modestly B-preferential boundaries, where B is the majority block, while a different (211)-oriented termination plane is stabilized by strongly B-preferential boundaries, neither of which has been observed experimentally. We then describe two particularly important phenomena for explaining the phase behavior of DG thin films at low film thicknesses. The first is "constructive interference," which arises when distortions due to the top and bottom boundaries overlap and is significant for certain DG orientations. The second is a symmetry-dependent, in-plane unit-cell distortion that arises because the distorted morphology near the boundary has a different preferred unit-cell size and shape than the bulk. These results provide a thermodynamic portrait of the phase behavior of DG thin films.

Lattice-Cluster-Theory-Informed Cross-Fractionation Chromatography Revealing Degree of Crystallinity of Single Macromolecular Species

The relationship between macromolecular architecture and crystallization properties is a relevant research topic in polymer science and technology. The average degree of crystallinity of disperse polymers is a well-studied quantity and is accessible by various experimental methods. However, how the different macromolecular species contribute to the degree of crystallinity and, in particular, the relationship between a certain macromolecular architecture and the degree of crystallinity are not accessible today, neither experimentally nor theoretically. Therefore, in this work, a lattice cluster theory (LCT)-informed cross-fractionation chromatography (CFC) approach is developed to access the degree of crystallinity of single and nonlinear macromolecular species crystallizing from solution. The method entangles high-throughput experimental data from CFC with the LCT for semicrystalline polymers to predict the degree of crystallinity of polymer species with different molecular weights and branching. The approach is applied to a linear low-density polyethylene (ethylene/1-octene copolymer) and a high-density polyethylene, which have specific and different bivariate distributions. The degree of crystallinity of individual macromolecular species of these polymer samples is calculated, and the predicted average degree of crystallinity is compared with experimental measurements, thus successfully validating the approach. Furthermore, the average segment length between branches is introduced as a characteristic molecular feature of branched polyethylene, and its relationship with the degree of crystallinity of certain species is established.

Predicting compatibilized polymer blend toughness

Polymer blends can yield superior materials by merging the unique properties of their components. However, these mixtures often phase separate, leading to brittleness. While compatibilizers can toughen these blends, their vast design space makes optimization difficult. Here, we develop a model to predict the toughness of compatibilized glassy polymer mixtures. This theory reveals that compatibilizers increase blend toughness by creating molecular bridges that stitch the interface together. We validate this theory by directly comparing its predictions to extensive molecular dynamics simulations in which we vary polymer incompatibility, chain stiffness, compatibilizer areal density, and blockiness of copolymer compatibilizers. We then parameterize the model using self-consistent field theory and confirm its ability to make predictions for practical applications through comparison with simulations and experiments. These results suggest that the theory can optimize compatibilizer design for industrial glassy polymer blends in silico while providing microscopic insight, allowing for the development of next-generation mixtures.

Boundary Frustration in Double-Gyroid Thin Films

Self-consistent field theory for thin films of AB diblock polymers in the double-gyroid phase reveals that in the absence of preferential wetting of monomer species at the film boundaries, films with the (211) plane oriented parallel to the boundaries are more stable than other orientations, consistent with experimental results. This preferred orientation is explained in the context of boundary frustration. Specifically, the angle of intersection between the A/B interface and the film boundary, the wetting angle, is thermodynamically restricted to a narrow range of values. Most termination planes in the double gyroid cannot accommodate this narrow range of wetting angles without significant local distortion relative to the bulk morphology; the (211)-oriented termination plane with the "double-wave" pattern produces relatively minimal distortion, making it the least frustrated boundary. The principle of boundary frustration provides a framework to understand the relative stability of termination planes for complex ordered block polymer phases confined between flat, nonpreferential boundaries.

Liquid-Like States in Micelle-Forming Diblock Copolymer Melts

Large cell self-consistent field theory (SCFT) solutions for a neat, micelle-forming diblock copolymer melt, initialized using the structure of a Lennard-Jones fluid, reveal the existence of a vast number of liquid-like states, with free energies of order 10^-3 k_BT per chain higher than the body-centered cubic (bcc) state near the order-disorder transition (ODT). Computation of the structure factor for these liquids at temperatures below the ODT indicates that their intermicellar distance is slightly swollen compared to bcc. In addition to providing a mean-field picture of the disordered micellar state, the number of liquid-like states and their near-degeneracy with the equilibrium bcc morphology suggest that self-assembly of micelle-forming diblock copolymers navigates a rugged free energy landscape with many local minima. This picture provides a basis for the anomalously slow ordering kinetics of particle-forming diblock copolymer melts observed in experiments.

Combining particle and field-theoretic polymer models with multi-representation simulations

Particle-based and field-theoretic simulations are both widely used methods to predict the properties of polymeric materials. In general, the advantages of each method are complementary. Field-theoretic simulations are preferred for polymers with high molecular weights and can provide direct access to chemical potentials and free energies, which makes them the method-of-choice for calculating phase diagrams. The trade-off is that field-theoretic simulations sacrifice the molecular details present in particle-based simulations, such as the configurations of individual molecules and their dynamics. In this work, we describe a new approach to conduct "multi-representation" simulations that efficiently map between particle-based and field-theoretic simulations. Our approach involves the construction of formally equivalent particle-based and field-based models, which are then simulated subject to the constraint that their spatial density profiles are equal. This constraint provides the ability to directly link particle-based and field-based simulations and enables calculations that can switch between one representation to the other. By switching between particle/field representations during a simulation, we demonstrate that our approach can leverage many of the advantages of each representation while avoiding their respective limitations. Although our method is illustrated in the context of complex sphere phases in linear diblock copolymers, we anticipate that it will be useful whenever free energies, rapid equilibration, molecular configurations, and dynamic information are all simultaneously desired.

Tuning conformational asymmetry in particle-forming diblock copolymer alloys

Self-consistent field theory is employed to compute the phase behavior of binary blends of conformationally asymmetric, micelle-forming diblock copolymers with miscible corona blocks and immiscible core blocks (a diblock copolymer "alloy"). The calculations focus on establishing conditions that promote the formation of Laves phases by tuning the relative softness of the cores of the two different Laves phase particles via independent control of their conformational asymmetries. Increasing the conformational asymmetry of the more spherical particles of the Laves structure has a stabilizing effect, consistent with the expectations of increased imprinting of the Wigner-Seitz cells on the core/corona interface as conformational asymmetry increases. The resulting phase diagram in the temperature-blend composition space features a more stable Laves phase field than that predicted for conformationally symmetric systems. The phase field closes at low temperatures in favor of macrophase separation between a hexagonally-packed cylinder (hex) phase and a body-centered cubic phase. Companion calculations, using an alloy whose components do not produce a hex phase in the neat melt state, suggest that the Laves phase field in such a blend will persist at strong segregation.

Stabilizing a Double Gyroid Network Phase with 2 nm Feature Size by Blending of Lamellar and Cylindrical Forming Block Oligomers

Molecular dynamics simulations are used to study binary blends of an AB-type diblock and an AB₂-type miktoarm triblock amphiphiles (also known as high-χ block oligomers) consisting of sugar-based (A) and hydrocarbon (B) blocks. In their pure form, the AB diblock and AB₂ triblock amphiphiles self-assemble into ordered lamellar (LAM) and cylindrical (CYL) structures, respectively. At intermediate compositions, however, the AB₂-rich blend (0.2 ≤ x _AB ≤ 0.4) forms a double gyroid (DG) network, whereas perforated lamellae (PL) are observed in the AB-rich blend (0.5 ≤ x _AB ≤ 0.8). All of the ordered mesophases present domain pitches under 3 nm, with 1 nm feature sizes for the polar domains. Structural analyses reveal that the nonuniform interfacial curvatures of DG and PL structures are supported by local composition variations of the LAM- and CYL-forming amphiphiles. Self-consistent mean field theory calculations for blends of related AB and AB₂ block polymers also show the DG network at intermediate compositions, when A is the minority block, but PL is not stable. This work provides molecular-level insights into how blending of shape-filling molecular architectures enables network phase formation with extremely small feature sizes over a wide composition range.

Open-source platform for block polymer formulation design using particle swarm optimization

Facile exploration of large design spaces is critical to the development of new functional soft materials, including self-assembling block polymers, and computational inverse design methodologies are a promising route to initialize this task. We present here an open-source software package coupling particle swarm optimization (PSO) with an existing open-source self-consistent field theory (SCFT) software for the inverse design of self-assembling block polymers to target bulk morphologies. To lower the barrier to use of the software and facilitate exploration of novel design spaces, the underlying SCFT calculations are seeded with algorithmically generated initial fields for four typical morphologies: lamellae, network phases, cylindrical phases, and spherical phases. In addition to its utility within PSO, the initial guess tool also finds generic applicability for stand-alone SCFT calculations. The robustness of the software is demonstrated with two searches for classical phases in the conformationally symmetric diblock system, as well as one search for the Frank-Kasper [Formula: see text] phase in conformationally asymmetric diblocks. The source code for both the initial guess generation and the PSO wrapper is publicly available.

Field-Theoretic Simulations for Block Copolymer Melts Using the Partial Saddle-Point Approximation

Field-theoretic simulations (FTS) provide an efficient technique for investigating fluctuation effects in block copolymer melts with numerous advantages over traditional particle-based simulations. For systems involving two components (i.e., A and B), the field-based Hamiltonian, Hf[W-,W+], depends on a composition field, W-(r), that controls the segregation of the unlike components and a pressure field, W+(r), that enforces incompressibility. This review introduces researchers to a promising variant of FTS, in which W-(r) fluctuates while W+(r) tracks its mean-field value. The method is described in detail for melts of AB diblock copolymer, covering its theoretical foundation through to its numerical implementation. We then illustrate its application for neat AB diblock copolymer melts, as well as ternary blends of AB diblock copolymer with its A- and B-type parent homopolymers. The review concludes by discussing the future outlook. To help researchers adopt the method, open-source code is provided that can be run on either central processing units (CPUs) or graphics processing units (GPUs).

RuSseL: A Self-Consistent Field Theory Code for Inhomogeneous Polymer Interphases

In this article, we publish the one-dimensional version of our in-house code, RuSseL, which has been developed to address polymeric interfaces through Self-Consistent Field calculations. RuSseL can be used for a wide variety of systems in planar and spherical geometries, such as free films, cavities, adsorbed polymer films, polymer-grafted surfaces, and nanoparticles in melt and vacuum phases. The code includes a wide variety of functional potentials for the description of solid–polymer interactions, allowing the user to tune the density profiles and the degree of wetting by the polymer melt. Based on the solution of the Edwards diffusion equation, the equilibrium structural properties and thermodynamics of polymer melts in contact with solid or gas surfaces can be described. We have extended the formulation of Schmid to investigate systems comprising polymer chains, which are chemically grafted on the solid surfaces. We present important details concerning the iterative scheme required to equilibrate the self-consistent field and provide a thorough description of the code. This article will serve as a technical reference for our works addressing one-dimensional polymer interphases with Self-Consistent Field theory. It has been prepared as a guide to anyone who wishes to reproduce our calculations. To this end, we discuss the current possibilities of the code, its performance, and some thoughts for future extensions.

Field theoretic approach for block polymer melts: SCFT and FTS

This perspective addresses the development of polymer field theory for predicting the equilibrium phase behavior of block polymer melts. The approach is tailored to the high-molecular-weight limit, where universality reduces all systems to the standard Gaussian chain model, an incompressible melt of elastic threads interacting by contact forces. Using mathematical identities, this particle-based version of the model is converted to an equivalent field-based version that depends on fields rather than particle coordinates. The statistical mechanics of the field-based model is typically solved using the saddle-point approximation of self-consistent field theory (SCFT), which equates to mean field theory, but it can also be evaluated using field theoretic simulations (FTS). While SCFT has matured into one of the most successful theories in soft condensed matter, FTS are still in its infancy. The two main obstacles of FTS are the high computational cost and the occurrence of an ultraviolet divergence, but fortunately there has been recent groundbreaking progress on both fronts. As such, FTS are now well poised to become the method of choice for predicting fluctuation corrections to mean field theory.