Improved convergence in block copolymer self-consistent field theory by Anderson mixing

Temperature dependence of micelle shape transitions in copolymer solutions: the role of inter-block incompatibility

The nature of the transition between worm-like and spherical micelles in block copolymer dispersions varies between systems. In some formulations, heating drives a transition from worms to spheres, while in other systems the same transition is induced by cooling. In addition, a sphere-worm interconversion can be accompanied either by an increase or a decrease in the core solvation, even if the direction of the temperature dependence is the same. Here, self-consistent field theory is used to provide a potential explanation of this range of behaviour. Specifically, we show that, within this model, the dependence of the transition on the incompatibility χBS of the solvophobic block B and the solvent S (the parameter most closely related to the temperature) is strongly influenced by the incompatibility χAB between B and the solvophilic block A. When χAB is small (χAB < 0.1), it is found that increasing χBS produces a transition from worm-like micelles to spheres (or, more generally, from less curved to more curved structures). When χAB is above 0.1, increasing χBS drives the system from spheres to worm-like micelles. Whether a transition is observed within a realistic range of χBS is also found to depend on the fraction of solvophilic material in the copolymer. The relevance of our calculations to experiments is discussed, and we suggest that the direction of the temperature dependence may be controlled not only by the solution behaviour of the solvophobic block (upper critical solution temperature-like versus lower critical solution temperature-like) but also by χAB.

The impact of intramolecular polydispersity on the self-assembly of ABn miktoarm star copolymers

The self-assembly behaviors of AB_n miktoarm star copolymers as one typical type of asymmetric architecture have been studied well in the past few decades due to their deflected phase boundaries. In particular, recently, they have attracted renewed theoretical interest due to their expanded spherical phase region that stabilizes complex Frank-Kasper spherical phases. However, previous theoretical studies have never considered AB_n copolymers with unequal arm lengths, which is more or less the case for synthesized copolymers. In this work, we investigate the self-assembly behaviors of AB_n miktoarm star copolymers with unequal B-arms using self-consistent field theory. We propose an intramolecular polydispersity index (iĐ) to quantify the distribution of unequal B-blocks. Accordingly, we further propose a simple quantity of an effective arm number n_equ = n/iĐ for quantitatively comparing the phase boundaries between various AB_n copolymer samples with different arm numbers or different distributions of B-blocks. Our results indicate that different AB_n copolymers with equal n_equ exhibit similar phase diagrams. On the other hand, we also found that the phase boundaries of two different samples with same n_equ are not exactly overlapped. We speculate that the effect of spontaneous curvature may be mainly controlled by n_equ, but the packing frustration of B-blocks may also be dependent on the other quantities that are closely related to the shape of the distribution of B-arms, such as higher order polydispersity indexes.

Theoretical Study of Phase Behaviors of Symmetric Linear B1A1B2A2B3 Pentablock Copolymer

The nanostructures that are self-assembled from block copolymer systems have attracted interest. Generally, it is believed that the dominating stable spherical phase is body-centered cubic (BCC) in linear AB-type block copolymer systems. The question of how to obtain spherical phases with other arrangements, such as the face-centered cubic (FCC) phase, has become a very interesting scientific problem. In this work, the phase behaviors of a symmetric linear B₁A₁B₂A₂B₃ (f_A1 = f_A2, f_B1 = f_B3) pentablock copolymer are studied using the self-consistent field theory (SCFT), from which the influence of the relative length of the bridging B₂-block on the formation of ordered nanostructures is revealed. By calculating the free energy of the candidate ordered phases, we determine that the stability regime of the BCC phase can be replaced by the FCC phase completely by tuning the length ratio of the middle bridging B₂-block, demonstrating the key role of B₂-block in stabilizing the spherical packing phase. More interestingly, the unusual phase transitions between the BCC and FCC spherical phases, i.e., BCC → FCC → BCC → FCC → BCC, are observed as the length of the bridging B₂-block increases. Even though the topology of the phase diagrams is less affected, the phase windows of the several ordered nanostructures are dramatically changed. Specifically, the changing of the bridging B₂-block can significantly adjust the asymmetrical phase regime of the Fddd network phase.

Efficient Langevin and Monte Carlo sampling algorithms: The case of field-theoretic simulations

We introduce Langevin sampling algorithms to field-theoretic simulations (FTSs) of polymers that, for the same accuracy, are ∼10× more efficient than a previously used Brownian dynamics algorithm that used predictor corrector for such simulations, over 10× more efficient than the smart Monte Carlo (SMC) algorithm, and typically over 1000× more efficient than a simple Monte Carlo (MC) algorithm. These algorithms are known as the Leimkuhler-Matthews (the BAOAB-limited) method and the BAOAB method. Furthermore, the FTS allows for an improved MC algorithm based on the Ornstein-Uhlenbeck process (OU MC), which is 2× more efficient than SMC. The system-size dependence of the efficiency for the sampling algorithms is presented, and it is shown that the aforementioned MC algorithms do not scale well with system sizes. Hence, for larger sizes, the efficiency difference between the Langevin and MC algorithms is even greater, although, for SMC and OU MC, the scaling is less unfavorable than for the simple MC.

Deep learning and inverse discovery of polymer self-consistent field theory inspired by physics-informed neural networks

We devise a deep learning solver inspired by physics-informed neural networks (PINNs) to tackle the polymer self-consistent field theory (SCFT) equations for one-dimensional AB-diblock copolymers. The PINNs framework comprises two parallel feedforward neural networks that separately represent the segmental partition functions and self-consistent chemical potential fields. The two networks are coupled through a loss function incorporating the governing equation, initial and boundary conditions, and the incompressibility constraint. To avoid the metastable homogeneous solution, the network parameters are initialized based on known self-consistent fields obtained from the numerical pseudospectral method. For copolymers of length N at a given volume fraction of A block (f) and the reduced Flory-Huggins interaction parameter (χN), the minimization of the loss function leads to the converged network parameters that successfully capture the stable lamellar phase. The periodicity of the lamellar structure is correctly reproduced for the explored sets of [f,χN], irrespective of the presumed computational domain size for initialization. Moreover, the proposed PINNs are applicable to the inverse discovery of the interaction parameter and the embedded chemical potential fields for an observed structure. This capability of solving the inverse SCFT problem demonstrates the potential of using PINNs to accelerate the exploration of new polymeric materials.

An experimental and theoretical investigation into the self-assembly of a chemically modified high-χ coil-rod diblock copolymer

A precursor diblock copolymer with a silicon backbone, polystyrene-block-poly(methylvinylsiloxane), was synthesized, and 1H,1H,2H,2H-perfluorodecanethiol was quantitatively introduced into the backbone via a thiol-ene reaction to yield a novel coil–rod diblock copolymer, poly(styrene-block-poly(2-((3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,10-heptadecafluorodecyl)thio)ethyl)methylsiloxane). The ultra-hydrophobicity of the introduced perfluoroalkyl side chain enhanced the segregation between counter-blocks and significantly increased the χ value, which is essential for minimizing the size of self-assembled domains for lithographic applications. Thus, self-assembled domains with a minimal spacing of approximately 10 nm were formed. A hexagonally packed array with significant anisotropy was observed in the self-assembled morphology by small-angle X-ray scattering and transmission electron microscopy. Such an array was precisely reproduced by modified self-consistent field theory (SCFT) calculation developed for the coil–rod structure. Furthermore, the phase diagram was estimated, and the morphological dependence on the relative scale of the rod unit was investigated by SCFT prediction.

A high-χ coil-rod diblock copolymer capable of forming unique anisotropic hexagonally packed cylindrical domains evidenced by experimental and theoretical study.

Tilt Modulus of Bilayer Membranes Self-Assembled from Rod-Coil Diblock Copolymers

Quantitatively understanding membrane fission and fusion requires a mathematical model taking their underlying elastic degrees of freedom, such as the molecule's tilt, into account. Hamm-Kozlov's model is such a framework that includes a tilt modulus along with the bending modulus and Gaussian modulus. This paper investigates the tilt modulus of liquid-crystalline bilayer membranes by applying self-consistent field theory. Unlike the widely used method in molecular dynamics simulation which extracts the tilt modulus by simulating bilayer buckles with various single modes, we introduce a tilt constrain term in the free energy to stabilize bilayers with various tilt angles. Fitting the energy curve as a function of the tilt angle to Hamm-Kozlov's elastic energy allows us to extract the tilt modulus directly. Based on this novel scheme and focused on the bilayers self-assembled from rod-coil diblock copolymers, we carry out a systematic study of the dependence of the tensionless A-phase bilayer's tilt modulus on the microscopic parameters.

Interplay between nematic and cholesteric interactions in self-consistent field theory

Chirality is a design feature of a number of biomolecules (e.g., collagen). In these molecules, cholesteric (chiral-nematic) behavior emerges from a combination of the tendency for the biopolymers to align (nematic interactions) and for the alignment direction to change with position, rotating around an axis normal to the alignment direction. This paper presents self-consistent field theory (SCFT) of chiral-nematic polymers, which takes into account polymer flexibility and the orientational degrees of freedom of polymer segments. Using the resulting SCFT, we construct a phase diagram showing regions of stability for isotropic, nematic, and cholesteric phases. Furthermore, we find that nematic interactions can stabilize the cholesteric phase, pushing the isotropic-cholesteric phase transition to lower cholesteric interaction strength, until the isotropic-nematic-cholesteric triple point is reached.

Periodic Patchy Spheres Self-Assembled by AmBCAn' Multiblock Terpolymers

We have designed A_mBCA_n^' multiblock terpolymers and studied their self-assembly using self-consistent field theory, aiming to generate the periodically arranged patchy spheres and thus to clarify the regulation mechanism of the number of patches. A number of two-dimensional phase diagrams are constructed for three typical architectures A₂BCA₂^', A₂BCA₃^', and A₃BCA₂^'. Four kinds of stable patchy spheres with the number of patches as 2 (S₂), 4 (S₄), 5 (S₅), and 6 (S₆) are obtained. These phases follow a common transition sequence of S₂ → S₄ → S₅ → S₆ along with the increasing of the volume fraction of C-block (f_C), which forms the core sphere patched with B-domains. Moreover, the S₆ phase exhibits the widest stability window, while S₅ has the narrowest one. The increased arms of A'-blocks in A₂BCA₃^' architecture deflect the phase boundaries toward large f_C and accordingly expand the regions of these patchy spheres due to the amplified effect of spontaneous curvature. In contrast, the increased arms of A-blocks in A₃BCA₂^' remarkably expands the window of S₆ but narrows those of the other patchy spheres, which is mainly caused by increased packing frustration resulting from the reduced extension of the more divided A-blocks. The widest window of the S₆ phase reaches Δf_C ∼ 0.13, which is readily accessed by experiment. Our work not only demonstrates a self-assembly strategy to engineer the patchy spheres, but also sheds light on the regulation mechanism of the patchy number.

Self-consistent field theory of chiral nematic worm-like chains

Many macromolecules of biological and technological interest are both chiral and semi-flexible. DNA and collagen are good examples. Such molecules often form chiral nematic (or cholesteric) phases, as is well-documented in collagen and chitin. This work presents a method for studying cholesteric phases in the highly successful self-consistent field theory of worm-like chains, offering a new way of studying many biologically relevant molecules. The method involves an effective Hamiltonian with a chiral term inspired by the Oseen-Frank (OF) model of liquid crystals. This method is then used to examine the formation of cholesteric phases in chiral-nematic worm-like chains as a function of polymer flexibility, as well as the optimal cholesteric pitch and distribution of polymer segment orientations. Our approach not only allows for the determination of the isotropic-cholesteric transition and segment distributions, beyond what the OF model promises, but also explicitly incorporates polymer flexibility into the study of the cholesteric phase, offering a more complete understanding of the behavior of semiflexible chiral-nematic polymers.

Molecular Bonding in an Orbital-Free-Related Density Functional Theory

A density functional theory based on polymer self-consistent field theory is applied to systems of two atoms in order to show that this approach is capable of predicted molecular bonding. Periodic table elements from hydrogen up to neon are examined and homonuclear diatomic molecules are found to form for H₂, N₂, O₂, and F₂, in agreement with known results. The heteronuclear molecules CO and HF, which are known to exist under ambient conditions, are also found to be stable. Bond lengths for most of these molecules agree with experimental results to within less than 8%, with the exception of O₂ and F₂ which deviate more significantly. The bonding energy for H₂ is given and is within 16% of the known value, but fundamental vibrational frequencies do not agree well with experiment. The main approximations of the theory are very simple and include a Fermi-Amaldi correction to the electron-electron interaction to account for self-interactions and a basic expression for the Pauli potential to account for the exclusion principle. The self-consistent equations are solved in terms of basis functions that encode the cylindrical symmetry of diatomic molecules. Since orbitals are not used, the approach is related to orbital-free density functional theory.

The stress deformation response influenced by the chain rigidity for mesostructures in diblock copolymers

A self-consistent field theory formalism based on the wormlike chain model is developed to investigate the stress-strain relation for mesostructures in diblock copolymers under the influence of chain rigidity, involving the adjustable simulation cell in the non-orthogonal coordinates by means of optimization of free energy. We elucidate the effect of the chain persistency broadly spanning from the Gaussian chain to the rigid rodlike chain on the elastic response of mesophases that deviate from the initial equilibrium structures. We analytically and numerically demonstrate that our current approach in the long chain limit recovers to the Gaussian-chain-based theory. Being ascribed to the distinct conformational behaviors for flexible chains and rigid rodlike chains, the tensile and compressive stresses applied to lamellae exhibit asymmetric deformation behaviors and the shear stress applied to the initial equilibrium hexagonal cylinders results in noticeable deviations in the shape and spatial arrangement of cylindroids for various chain rigidity values. For the zero stress, in addition, our approach can be straightforwardly utilized to explore the optimal size and shape of the simulation cell in order to achieve a stress free configuration of systems.

The C36 Laves phase in diblock polymer melts

The C14 and C15 Laves phases form as micelle packing structures in many types of soft matter, but the related C36 phase, which consists of alternating C14-type and C15-type layers, has not been observed in any such system. To understand this absence in the context of diblock polymers, we used self-consistent field theory to relate the morphology and energetics of C36 to other known mesophases. Two case studies were conducted: blends of AB diblock polymers with A homopolymers (where A forms the micelle core), in which C14 and C15 have stability windows, and neat AB diblock melts, in which Laves phases are metastable. Laves phases exhibit nearly identical micelle morphologies and nearly degenerate free energies, with the free energy of C36 being a near-perfect bisector of the C14 and C15 free energies in all cases, revealing an intrinsic symmetry in free energy that is attributed solely to the structural relationship between the phases in which the packing of C36 is intermediate between C14 and C15. Based on this connection between structure and free energy, C36 is thus not expected to form in flexible diblock polymers, since C14 and C15 can always form instead via facile mass transfer.

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

Field-theoretic simulations beyond δ-interactions: Overcoming the inverse potential problem in auxiliary field models

Modern field-theoretic simulations of complex fluids and polymers are constructed around a particle-to-field transformation that brings an inverse potential u^-1 in the model equations. This has restricted the application of the framework to systems characterized by relatively simple pairwise interatomic interactions; for example, excluded volume effects are treated through the use of δ-function interactions. In this study, we first review available nonbonded pair interactions in field-theoretic models and propose a classification. Then, we outline the inverse potential problem and present an alternative approach on the basis of a saddle-point approximation, enabling the use of a richer set of pair interaction functions. We test our approach by using as an example the Morse potential, which finds extensive applications in particle-based simulations, and we calibrate u^-1 with results from a molecular dynamics simulation. The u^-1 thus obtained is consistent with the field-theoretic model equations, and when used in stand-alone self-consistent field simulations, it produces the correct fluid structure starting from a random initial state of the density field.

Analytical Calculation of the Elastic Moduli of Self-Assembled Liquid-Crystalline Bilayer Membranes

Liquid-crystalline orders are ubiquitous in membranes and could significantly affect the elastic properties of the self-assembled bilayers. Calculating the free energy of bilayer membranes with different geometries and fitting them to their theoretical expressions allow us to extract the elastic moduli, such as the bending modulus and Gaussian modulus. However, this procedure is time-consuming for liquid-crystalline bilayers. In paper reports a novel method to calculate the elastic moduli of the self-assembled liquid-crystalline bilayers within the self-consistent field theory framework. Based on the asymptotic expansion method, we derive the analytical expression of the elastic moduli, which reduces the computational cost significantly. Numerical simulations illustrate the validity and efficiency of the proposed method.

Frank-Kasper Phases Self-Assembled from a Linear A1B1A2B2 Tetrablock Copolymer

The Frank-Kasper (FK) phases self-assembled from block copolymer systems have attracted abiding interest. In this work, the formation mechanism of the complex FK phases from the self-assembly of simple A₁B₁A₂B₂ tetrablock copolymers is investigated using self-consistent field theory (SCFT). For a typical set of parameter spaces, we utilize SCFT to construct a number of phase diagrams. In these phase diagrams, the FK phases exhibit a notable stability region. The stable region of the FK phases reveals that the distribution of A₁ and A₂ blocks can be precisely regulated by tuning the ratio of the A₁/A₂ block, wherein the long A₁ blocks can aggregate within the "core" while the short A₂ blocks can form the "shell" of a spherical domain in the FK phases, respectively, to accommodate the sizes and shapes of the spherical domains in the complex spherical packing phases. Besides, we also demonstrate that the existence of the B₂ block plays a crucial factor to stabilize the FK phases.

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¯=10⁴, 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.

Elastic properties of self-assembled bilayer membranes: Analytic expressions via asymptotic expansion

Bilayer membranes self-assembled from amphiphilic molecules are ubiquitous in biological and soft matter systems. The elastic properties of bilayer membranes are essential in determining the shape and structure of bilayers. A novel method to calculate the elastic moduli of the self-assembled bilayers within the framework of the self-consistent field theory is developed based on an asymptotic expansion of the order parameters in terms of the bilayer curvature. In particular, the asymptotic expansion method is used to derive analytic expressions of the elastic moduli, which allows us to design more efficient numerical schemes. The efficiency of the proposed method is illustrated by a model system composed of flexible amphiphilic chains dissolved in hydrophilic polymeric solvents.

The size and affinity effect of counterions on self-assembly of charged block copolymers

The effect of counterions' size and affinity on the microphase separated morphologies of neutral-charged diblock copolymers is investigated systematically using a random phase approximation (RPA) and self-consistent field theory (SCFT). The phase diagrams as a function of χ_AB and f_A at different counterion sizes and different affinities to neutral blocks are constructed, respectively. Stability limits calculated using the RPA are in good agreement with the disorder-body-centered cubic phase boundaries from SCFT calculations. It was found that increasing the size of counterions causes the phase diagram to shift upward and leftward, which is attributed to electrostatic interactions and the intrinsic volume of counterions. The domain size of the ordered phase shows an unexpected tendency that it decreases with increasing counterions' size. The counterions' distributions in H and G phases demonstrate that it is electrostatic interaction, instead of packing frustration, that plays a leading role in such systems. For finite size counterions, with the increase in affinity between counterions and neutral blocks, the phase diagram shifts upward, indicating the improved compatibility between different blocks. Furthermore, the affinity effect between counterions and neutral blocks can be mapped into an effective Flory parameter χ_AB ^' = χ_AB + 0.27χ_BC.

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.

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

Self-consistent field theory (SCFT) is a powerful approach for computing the phase behavior of block polymers. We describe a fast version of the open-source Polymer Self-Consistent Field (PSCF) code that takes advantage of the massive parallelization provided by a graphical processing unit (GPU). Benchmarking double-precision calculations indicate up to 30× reduction in time to converge SCFT calculations of various diblock copolymer phases when compared to the Fortran CPU version of PSCF using the same algorithms, with the speed-up increasing with increasing unit cell size for the diblock polymer problems examined here. Where double-precision accuracy is not needed, single-precision calculations can provide speed-up of up to 60× in convergence time. These improvements in speed within an open-source format open up new vistas for SCFT-driven block polymer materials discovery by the community at large.

Elastic properties of liquid-crystalline bilayers self-assembled from semiflexible-flexible diblock copolymers

The mechanical response and shape of self-assembled bilayer membranes depend crucially on their elastic properties. Most of the studies focused on the elastic properties of fluid membranes, despite the ubiquitous presence of membranes with liquid-crystalline order. Here the elastic properties of liquid-crystalline bilayers self-assembled from diblock copolymers composed of a semiflexible block are studied theoretically. Specifically, the self-consistent field theory (SCFT) is applied to a model system composed of semiflexible-flexible diblock copolymers dissolved in flexible homopolymers that act as solvents. The free energy of self-assembled tensionless bilayer membranes in three different geometries, i.e. planar, cylindrical and spherical, is obtained by solving the SCFT equations using a hybrid method, in which the orientation-dependent functions are treated using the spherical harmonics, whereas the position-dependent operators are treated using the compact difference schemes. The bending modulus κ_M and Gaussian modulus κ_G of the bilayer are extracted from the free energies. The effects of the molecular parameters of the system, such as the chain rigidity and the orientational interaction, are systematically examined.