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

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.

Assessment of the partial saddle point approximation in field-theoretic polymer simulations

Field-theoretic simulations are numerical treatments of polymer field theory models that go beyond the mean-field self-consistent field theory level and have successfully captured a range of mesoscopic phenomena. Inherent in molecularly-based field theories is a "sign problem" associated with complex-valued Hamiltonian functionals. One route to field-theoretic simulations utilizes the complex Langevin (CL) method to importance sample complex-valued field configurations to bypass the sign problem. Although CL is exact in principle, it can be difficult to stabilize in strongly fluctuating systems. An alternate approach for blends or block copolymers with two segment species is to make a "partial saddle point approximation" (PSPA) in which the stiff pressure-like field is constrained to its mean-field value, eliminating the sign problem in the remaining field theory, allowing for traditional (real) sampling methods. The consequences of the PSPA are relatively unknown, and direct comparisons between the two methods are limited. Here, we quantitatively compare thermodynamic observables, order-disorder transitions, and periodic domain sizes predicted by the two approaches for a weakly compressible model of AB diblock copolymers. Using Gaussian fluctuation analysis, we validate our simulation observations, finding that the PSPA incorrectly captures trends in fluctuation corrections to certain thermodynamic observables, microdomain spacing, and location of order-disorder transitions. For incompressible models with contact interactions, we find similar discrepancies between the predictions of CL and PSPA, but these can be minimized by regularization procedures such as Morse calibration. These findings mandate caution in applying the PSPA to broader classes of soft-matter models and systems.

Fluctuation-Corrected Phase Diagrams for Diblock Copolymer Melts

New developments in field-theoretic simulations (FTSs) are used to evaluate fluctuation corrections to the self-consistent field theory of diblock copolymer melts. Conventional simulations have been limited to the order-disorder transition (ODT), whereas FTSs allow us to evaluate complete phase diagrams for a series of invariant polymerization indices. The fluctuations stabilize the disordered phase, which shifts the ODT to higher segregation. Furthermore, they stabilize the network phases at the expense of the lamellar phase, which accounts for the presence of the Fddd phase in experiments. We hypothesize that this is due to an undulation entropy that favors curved interfaces.

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+c₂χ_b ²+c₃χ_b ³+⋯, 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 χ.

Well-tempered metadynamics applied to field-theoretic simulations of diblock copolymer melts

Well-tempered metadynamics (WTMD) is applied to field-theoretic simulations (FTS) to locate the order-disorder transition (ODT) in incompressible melts of diblock copolymer with an invariant polymerization index of N̄=10⁴. The polymers are modeled as discrete Gaussian chains with N = 90 monomers, and the incompressibility is treated by a partial saddle-point approximation. Our implementation of WTMD proves effective at locating the ODT of the lamellar and cylindrical regions, but it has difficulty with that of the spherical and gyroid regions. In the latter two cases, our choice of order parameter cannot sufficiently distinguish the ordered and disordered states because of the similarity in microstructures. The gyroid phase has the added complication that it competes with a number of other morphologies, and thus, it might be beneficial to extend the WTMD to multiple order parameters. Nevertheless, when the method works, the ODT can be located with impressive accuracy (e.g., ΔχN ∼ 0.01).

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).

Study of Segregation in Non-Dilute Solutions of Linear Diblock Copolymers and Symmetric Miktoarm (or Janus Star) Polymers Using Monte Carlo Simulations with the Bond Fluctuation Model

The bond fluctuation model was employed to characterize the approach to the mesophase separation transition of pure linear AB copolymers and symmetric miktoarms, also called Janus, star polymers, A_f/2B_f/2, where f = 6 or 12 is the total number of arms, in a common good solvent. We consider a concentration sufficiently high to mimic the melting behavior and also a lower concentration. The segregation between A and B units is represented by a repulsive interaction parameter, ε. Different total numbers of units are also considered. Results for different properties, such as the molecular size, the asphericity and orientational correlation of blocks, or arms, of different compositions are obtained as a function of the segregation parameter. We also calculate scattering structure factors. The initial effect of segregation on the scattering with opposite contrast factors between the A and B blocks can be explained with a common description based on the random phase approximation for both the linear copolymers and the f = 6 miktoarms, once the numerical form factors of the different molecules in their particular systems are considered. However, the results for f = 12 clearly deviate from this description probably due to some degree of ordering in the position of highly armed molecules.