Bicontinuous microemulsion in binary blends of complementary diblock copolymers.

The phase behavior of binary blends of AB diblock copolymers of compositions f and 1 - f is examined using field-theoretic simulations. Highly asymmetric compositions (i.e., f ≈ 0) behave like homopolymer blends macrophase separating into coexisting A- and B-rich phases as the segregation is increased, whereas more symmetric diblocks (i.e., f ≈ 0.5) microphase separate into an ordered lamellar phase. In self-consistent field theory, these behaviors are separated by a Lifshitz critical point at f = 0.2113. However, its lower critical dimension is believed to be four, which implies that the Lifshitz point should be destroyed by fluctuations. Consistent with this, it is found to transform into a tricritical point. Furthermore, the highly swollen lamellar phase near the mean-field Lifshitz point disorders into a bicontinuous microemulsion (BµE), consisting of large interpenetrating A- and B-rich microdomains. BµE has been previously reported in ternary blends of AB diblock copolymer with its parent A- and B-type homopolymers, but in that system the homopolymers have a tendency to macrophase separate. Our alternative system for creating BµE is free of this macrophase separation.

Accounting for the ultraviolet divergence in field-theoretic simulations of block copolymer melts.

This study examines the ultraviolet (UV) divergence in field-theoretic simulations (FTSs) of block copolymer melts, which causes an unphysical dependence on the grid resolution, Δ, used to represent the fields. Our FTSs use the discrete Gaussian-chain model and a partial saddle-point approximation to enforce incompressibility. Previous work has demonstrated that the UV divergence can be accounted for by defining an effective interaction parameter, χ=z∞χb+c2χb 2+c3χb 3+⋯, in terms of the bare interaction parameter, χb, used in the FTSs, where the coefficients of the expansion are determined by a Morse calibration. However, the need to use different grid resolutions for different ordered phases generally restricts the calibration to the linear approximation, χ ≈ z∞χb, and prevents the calculation of order-order transitions. Here, we resolve these two issues by showing how the nonlinear calibration can be translated between different grids and how the UV divergence can be removed from free energy calculations. By doing so, we confirm previous observations from particle-based simulations. In particular, we show that the free energy closely matches self-consistent field theory (SCFT) predictions, even in the region where fluctuations disorder the periodic morphologies, and similarly, the periods of the ordered phases match SCFT predictions, provided the SCFT is evaluated with the nonlinear χ.

Entropic surface segregation from athermal polymer blends: Polymer flexibility vs bulkiness.

We examine athermal binary blends composed of conformationally asymmetric polymers of equal molecular volume next to a surface of width ξ. The self-consistent field theory (SCFT) of Gaussian chains predicts that the more compact polymer with the shorter average end-to-end length, R0, is entropically favored at the surface. Here, we extend the SCFT to worm-like chains with small persistence lengths, ℓp, relative to their contour lengths, ℓc, for which R0≈2ℓpℓc. In the limit of ℓp ≪ ξ, we recover the Gaussian-chain prediction where the segregation depends only on the product ℓpℓc, but for realistic polymer/air surfaces with ξ ∼ ℓp, the segregation depends separately on the two quantities. Although the surface continues to favor flexible polymers with smaller ℓp and bulky polymers with shorter ℓc, the effect of bulkiness is more pronounced. This imbalance can, under specific conditions, lead to anomalous surface segregation of the more extended polymer. For this to happen, the polymer must be bulkier and stiffer, with a stiffness that is sufficient to produce a larger R0 yet not so rigid as to reverse the surface affinity that favors bulky polymers.

Fluctuation correction for the order-disorder transition of diblock copolymer melts.

The order-disorder transition (ODT) of diblock copolymer melts is evaluated for an invariant polymerization index of N¯=104, using field-theoretic simulations (FTS) supplemented by a partial saddle-point approximation for incompressibility. For computational efficiency, the FTS are performed using the discrete Gaussian-chain model, and results are then mapped onto the continuous model using a linear approximation for the Flory-Huggins χ parameter. Particular attention is paid to the complex phase window. Results are found to be consistent with the well-established understanding that the gyroid phase extends down to the ODT. Furthermore, our simulations are the first to predict that the Fddd phase survives fluctuation effects, consistent with experiments.

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.

Calibration of a lattice model for high-molecular-weight block copolymer melts.

The Morse calibration is applied to a lattice model designed for efficient simulations of two-component polymer melts of high molecular weight. The model allows multiple occupancy per site, which results in high invariant polymerization indices, and interactions are limited to monomers within the same site, which enhances the computational speed. The calibration maps the interaction parameter of the lattice model, α, onto the Flory-Huggins χ parameter of the standard Gaussian-chain model, by matching the disordered-state structure function, S(k), of symmetric diblock copolymers to renormalized one-loop predictions. The quantitative accuracy of the calibration is tested by comparing the order-disorder transition of symmetric diblock copolymer melts to the universal prediction obtained from previous simulations. The model is then used to confirm the universality of fluctuation corrections to the critical point of symmetric binary homopolymer blends.

Calibration of the Flory-Huggins interaction parameter in field-theoretic simulations.

Field-theoretic simulations (FTS) offer a versatile method of dealing with complicated block copolymer systems, but unfortunately they struggle to cope with the level of fluctuations typical of experiments. Although the main obstacle, an ultraviolet divergence, can be removed by renormalizing the Flory-Huggins χ parameter, this only works for unrealistically large invariant polymerization indexes, N¯. Here, we circumvent the problem by applying the Morse calibration, where a nonlinear relationship between the bare χb used in FTS and the effective χ corresponding to the standard Gaussian-chain model is obtained by matching the disordered-state structure function, S(k), of symmetric diblock copolymers to renormalized one-loop predictions. This calibration brings the order-disorder transition obtained from FTS into agreement with the universal results of particle-based simulations for values of N¯ characteristic of the experiment. In the limit of weak interactions, the calibration reduces to a linear approximation, χ ≈ z∞χb, consistent with the previous renormalization of χ for large N¯.

Effect of chain stiffness on the entropic segregation of chain ends to the surface of a polymer melt.

Entropic segregation of chain ends to the surface of a monodisperse polymer melt and its effect on surface tension are examined using self-consistent field theory (SCFT). In order to assess the dependence on chain stiffness, the SCFT is solved for worm-like chains. Our focus is still on relatively flexible polymers, where the persistence length of the polymer, ℓ p , is comparable to the width of the surface profile, ξ, but still much smaller than the total contour length of the polymer, ℓ c . Even this small degree of rigidity causes a substantial increase in the level of segregation, relative to that of totally flexible Gaussian chains. Nevertheless, the long-range depletion that balances the surface excess still exhibits the same universal shape derived for Gaussian chains. Furthermore, the excess continues to reduce the surface tension by one unit of k B T per chain end, which results in the usual N -1 reduction in surface tension observed by experiments. This enhanced segregation will also extend to polydisperse melts, causing the molecular-weight distribution at the surface to shift towards smaller N n relative to the bulk. This provides a partial explanation for recent quantitative differences between experiments and SCFT calculations for flexible polymers.

Entropic segregation of short polymers to the surface of a polydisperse melt.

Chain ends are known to have an entropic preference for the surface of a polymer melt, which in turn is expected to cause the short chains of a polydisperse melt to segregate to the surface. Here, we examine this entropic segregation for a bidisperse melt of short and long polymers, using self-consistent field theory (SCFT). The individual polymers are modeled by discrete monomers connected by freely-jointed bonds of statistical length a , and the field is adjusted so as to produce a specified surface profile of width [Formula: see text]. Semi-analytical expressions for the excess concentration of short polymers, [Formula: see text], the integrated excess, [Formula: see text] , and the entropic effect on the surface tension, [Formula: see text], are derived and tested against the numerical SCFT. The expressions exhibit universal dependences on the molecular-weight distribution with model-dependent coefficients. In general, the coefficients have to be evaluated numerically, but they can be approximated analytically once [Formula: see text]. We illustrate how this can be used to derive a simple expression for the interfacial tension between immiscible A- and B-type polydisperse homopolymers.

Segregation of chain ends to the surface of a polymer melt: Effect of surface profile versus chain discreteness.

Silberberg has argued that the surface of a polymer melt behaves like a reflecting boundary on the random-walk statistics of the polymers. Although this is approximately true, independent studies have shown that violations occur due to the finite width of the surface profile and to the discreteness of the polymer molecule, resulting in an excess of chain ends at the surface and a reduction in surface tension inversely proportional to the chain length, N . Using self-consistent field theory (SCFT), we compare the magnitude of these two effects by examining a melt of discrete polymers modeled as N monomers connected by Hookean springs of average length, a , next to a polymer surface of width [Formula: see text]. The effects of the surface width and the chain discreteness are found to be comparable for realistic profiles of [Formula: see text] ∼ a. A semi-analytical approximation is developed to help explain the behavior. The relative excess of ends at the surface is dependent on the details of the model, but in general it decreases for shorter polymers. The excess is balanced by a long-range depletion that has a universal shape independent of the molecular details. Furthermore, the approximation predicts that the reduction in surface energy equals one unit of kBT for every extra chain end at the surface.

Segregation of chain ends to the surface of a polymer melt.

Self-consistent field theory (SCFT) is used to examine the surface of an incompressible polymer melt of freely jointed chains, each consisting of N discrete monomers connected by bonds of an arbitrary potential. As a result of entropic considerations, the end monomers tend to accumulate in a narrow region next to the surface (on the monomer scale), which causes a slight depletion of ends further into the melt (on the molecular scale). Due to the reduced configurational entropy available to polymers in the vicinity of a surface, we find an entropic contribution to the surface tension that increases with the degree of polymerization, N. While many quantities are dependent on the precise bond potential connecting the monomers, the excess of ends at the surface in the limit of infinite N, the functional form of the long-range depletion of ends, and the N dependence of the surface tension turn out to be universal.

Comparison of A-block polydispersity effects on BAB triblock and AB diblock copolymer melts.

Recent experiments on triblock copolymer melts suggest that polydispersity effects are dramatically enhanced when polydisperse blocks are constrained by both ends to the internal interfaces of an ordered morphology. To quantify the relevance of architecture, we compare BAB triblock and AB diblock copolymer melts with polydisperse A blocks and monodisperse B blocks, using self-consistent field theory (SCFT). We do, in fact, find an enhanced shift in the order-order transitions (OOTs) of the triblock copolymer system in good agreement with the experiments, which we attribute to a reduction of entropy in the A-rich domains due to the absence of chain ends. There is also a slightly enhanced dilation of the domains, but not nearly to the same degree as reported by the experiments. Unlike in the experiments, our calculations indicate that the polydispersity-induced shifts in the order-disorder transition (ODT) should be quantitatively similar for both diblocks and triblocks. It is possible that some of the pronounced effects observed in the experiments have more to do with the detailed shape of the molecular-weight distribution than the triblock architecture.

Effect of salt on the compression of polyelectrolyte brushes in a theta solvent.

Classical strong-stretching theory (SST) predicts that, as opposing polyelectrolyte brushes are compressed together in a salt-free theta solvent, they contract so as to maintain a finite polymer-free gap, which offers a potential explanation for the ultra-low frictional forces observed in experiments despite the application of large normal forces. However, the SST ignores chain fluctuations, which would tend to close the gap resulting in physical contact and in turn significant friction. In a preceding study, we examined the effect of fluctuations using self-consistent field theory (SCFT) and illustrated that high normal forces can still be applied before the gap is destroyed. We now look at the effect of adding salt. It is found to reduce the long-range interaction between the brushes but has little effect on the short-range part, provided the concentration does not enter the salted-brush regime. Consequently, the maximum normal force between two planar brushes at the point of contact is remarkably unaffected by salt. For the crossed-cylinder geometry commonly used in experiments, however, there is a gradual reduction because in this case the long-range part of the interaction contributes to the maximum normal force.

Dynamics of interacting edge defects in copolymer lamellae.

It is known that terraces at the air-polymer interface of lamella-forming diblock copolymers do not make discontinuous jumps in height. Despite the underlying discretized structure, the height profiles are smoothly varying. The width of a transition region of a terrace edge in isolation is typically several hundreds of nanometres, resulting from a balance between surface tension, chain stretching penalties, and the enthalpy of mixing. What is less well known in these systems is what happens when two transition regions interact with one another. In this study, we investigate the dynamics of the interactions between copolymer lamellar edges. We find that the data can be well described by a model that assumes a repulsion between adjacent edges. While the model is simplistic, and does not include molecular level details, its agreement with the data suggests that some of the the underlying assumptions provide insight into the complex interplay between defects.

Efficiency of pseudo-spectral algorithms with Anderson mixing for the SCFT of periodic block-copolymer phases.

This study examines the numerical accuracy, computational cost, and memory requirements of self-consistent field theory (SCFT) calculations when the diffusion equations are solved with various pseudo-spectral methods and the mean-field equations are iterated with Anderson mixing. The different methods are tested on the triply periodic gyroid and spherical phases of a diblock-copolymer melt over a range of intermediate segregations. Anderson mixing is found to be somewhat less effective than when combined with the full-spectral method, but it nevertheless functions admirably well provided that a large number of histories is used. Of the different pseudo-spectral algorithms, the 4th-order one of Ranjan, Qin and Morse performs best, although not quite as efficiently as the full-spectral method.

Reversible sphere-to-lamellar wetting transition at the interface of a diblock copolymer system.

We use ellipsometry to investigate a transition in the morphology of a sphere-forming diblock copolymer thin-film system. At an interface the diblock morphology may differ from the bulk when the interfacial tension favours wetting of the minority domain, thereby inducing a sphere-to-lamella transition. In a small, favourable window in energetics, one may observe this transition simply by adjusting the temperature. Ellipsometry is ideally suited to the study of the transition because the additional interface created by the wetting layer affects the polarisation of light reflected from the sample. Here we study thin films of poly(butadiene-ethylene oxide) (PB-PEO), which order to form PEO minority spheres in a PB matrix. As temperature is varied, the reversible transition from a partially wetting layer of PEO spheres to a full wetting layer at the substrate is investigated.

Compression of polyelectrolyte brushes in a salt-free theta solvent.

This paper examines the normal force between two opposing polyelectrolyte brushes and the interpenetration of their chains that is responsible for sliding friction. It focuses on the special case of semi-dilute brushes in a salt-free theta solvent, for which Zhulina and Borisov (J. Chem. Phys. 107, 5952 (1997)) have derived analytical predictions using the classical strong-stretching theory (SST). Interestingly, SST predicts that the brushes contract as they are compressed together maintaining a polymer-free gap, which provides an explanation for the ultra-low frictional forces observed in experiment. We examine the degree to which the SST predictions are affected by chain fluctuations by employing self-consistent field theory (SCFT). While the normal force is relatively unaffected, fluctuations are found to have a strong impact on brush interpenetration. Even still, the contraction of the brushes does significantly prolong the onset of interpenetration, implying that a sizeable normal force can be achieved before the sliding friction becomes significant.

Strong-segregation limit of the self-consistent field theory for diblock copolymer melts.

The self-consistent field theory (SCFT) introduced by Helfand for diblock copolymer melts is expected to converge to the strong-segregation theory (SST) of Semenov in the asymptotic limit, χN --> ∞ . However, past extrapolations of the lamellar/cylinder and cylinder/sphere phase boundaries, within the standard unit-cell approximation, have cast some doubts on whether or not this is actually true. Here we push the comparison further by extending the SCFT calculations to χN = 512000 , by accounting for exclusion zones in the coronae of the cylindrical and spherical unit cells, and by examining finite-segregation corrections to SST. In doing so, we provide the first compelling evidence that SCFT does indeed reduce to SST.

Monte Carlo phase diagram for diblock copolymer melts.

The phase diagram for diblock copolymer melts is evaluated from lattice-based Monte Carlo simulations using parallel tempering, improving upon earlier simulations that used sequential temperature scans. This new approach locates the order-disorder transition (ODT) far more accurately by the occurrence of a sharp spike in the heat capacity. The present study also performs a more thorough investigation of finite-size effects, which reveals that the gyroid (G) morphology spontaneously forms in place of the perforated-lamellar (PL) phase identified in the earlier study. Nevertheless, there still remains a small region where the PL phase appears to be stable. Interestingly, the lamellar (L) phase next to this region exhibits a small population of transient perforations, which may explain previous scattering experiments suggesting a modulated-lamellar (ML) phase.

Fast and accurate SCFT calculations for periodic block-copolymer morphologies using the spectral method with Anderson mixing.

We study the numerical efficiency of solving the self-consistent field theory (SCFT) for periodic block-copolymer morphologies by combining the spectral method with Anderson mixing. Using AB diblock-copolymer melts as an example, we demonstrate that this approach can be orders of magnitude faster than competing methods, permitting precise calculations with relatively little computational cost. Moreover, our results raise significant doubts that the gyroid (G) phase extends to infinite chi N . With the increased precision, we are also able to resolve subtle free-energy differences, allowing us to investigate the layer stacking in the perforated-lamellar (PL) phase and the lattice arrangement of the close-packed spherical ( S (cp) phase. Furthermore, our study sheds light on the existence of the newly discovered Fddd ( O(70) morphology, showing that conformational asymmetry has a significant effect on its stability.