🔬 condensed matter

A flexible implementation of strong segregation theory for two dimensional ABC star terpolymer morphologies

This paper presents a novel computational implementation of strong segregation theory that models two-dimensional ABC star terpolymer morphologies by assembling complex structures from flexible "Strongly Segregated Polygons" to efficiently calculate free energies and construct phase diagrams.

Original authors: Merin Joseph, Daniel J. Read, Alastair M. Rucklidge

Published 2026-05-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Merin Joseph, Daniel J. Read, Alastair M. Rucklidge

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have a bag of three different types of magical, sticky string: Red, Blue, and Yellow. You tie them together at a single knot in the middle to make a three-armed star. If you have a huge pile of these stars and you heat them up, they don't just mix into a messy soup. Instead, because the colors hate touching each other, they try to sort themselves out. The Red strings want to be with other Reds, the Blues with Blues, and the Yellows with Yellows.

But here's the catch: they are all tied together at the knot. They can't run away to their own separate rooms. So, they have to arrange themselves in a very specific, organized pattern where the knots stay in the middle, and the strings stretch out into their own colored zones.

This paper is about a new, clever way to predict exactly what those patterns will look like.

The Problem: Too Many Shapes to Count

Scientists have known for a long time that these "ABC star" molecules can form many different shapes, like honeycombs, squares, or complex mosaics. To figure out which shape is the most stable (the one nature will actually choose), you usually have to do incredibly difficult math. It's like trying to solve a giant 3D puzzle where every piece is constantly changing shape. Doing this for every possible mixture of Red, Blue, and Yellow strings takes a supercomputer a very long time.

The Solution: The "Strongly Segregated Polygon" (SSP)

The authors of this paper created a simpler, faster tool. They realized that no matter how complex the final pattern looks, it can be built by tiling together a few basic building blocks. They call these blocks Strongly Segregated Polygons (SSPs).

Think of an SSP like a slice of pizza, but instead of a triangle, it's a six-sided shape (a hexagon) with a knot right in the center.

The Center: This is where the Red, Blue, and Yellow strings are tied together.
The Slices: The pizza is cut into six triangular slices. Two slices are Red, two are Blue, and two are Yellow.
The Stretch: The strings stretch out from the center knot to the edge of the slice.

The magic of their method is that they treat these shapes like flexible tiles. You can take a bunch of these "pizza slices," glue them together edge-to-edge (making sure Red touches Red, Blue touches Blue, etc.), and then gently push and pull the corners of the tiles until the whole pattern settles into the most comfortable, lowest-energy position.

How It Works in Practice

The Building Blocks: Instead of trying to model the whole messy blob at once, the computer just looks at one of these six-sided tiles. It calculates the energy of the strings stretching and the cost of the different colors touching each other.
Tessellation: They repeat this tile over and over to build a full pattern (like tiling a bathroom floor).
The Optimization: The computer then plays a game of "tug-of-war." It moves the corners of the tiles slightly to see if the whole pattern becomes more stable. It keeps doing this until it finds the perfect shape for that specific mix of Red, Blue, and Yellow.

What They Found

Using this "tile-pushing" method, the authors mapped out a "menu" of all the possible shapes these molecules can make.

Simple Shapes: When you have equal amounts of all three colors, they form a perfect honeycomb pattern (the [6.6.6] shape).
Complex Shapes: When you have more of one color than the others, the pattern shifts. You get squares, triangles, or even more complicated mosaics where some knots are surrounded by 8 neighbors and others by 4.
The "Lamellar" Surprise: They found that when the interactions between the colors are very strong, some patterns turn into layers (like a lasagna) mixed with cylinders. This was a bit of a surprise in their calculations compared to previous, more complex computer models.

Why This Matters

The authors didn't just find new shapes; they found a shortcut.

Speed: Their method is much faster than the traditional "heavy lifting" methods used by other scientists. It allows them to scan thousands of different mixtures in the time it would take others to check just a few.
Flexibility: Because their method is based on simple tiles, they can easily change the rules. They can make the Red strings "hate" the Blue strings more than the Yellow ones, and instantly see how the whole pattern shifts.
Accuracy: Even though it's a simpler method, their results matched up very well with the more complex, expensive computer models and even real-world experiments.

The Bottom Line

Think of this paper as inventing a new, super-efficient way to design a mosaic floor. Instead of trying to calculate the stress on every single tile in a massive, complex floor, the authors realized you only need to understand how one specific type of tile behaves and how it fits with its neighbors. By mastering that one tile, they could quickly predict the design of the entire floor for any combination of colors, saving time and computing power while still getting the right answer.

Technical Summary: A Flexible Implementation of Strong Segregation Theory for 2D ABC Star Terpolymer Morphologies

Problem Statement
The study of microphase separation in ABC star terpolymers is complicated by the vast variety of possible morphologies, ranging from simple periodic tilings to complex multi-core structures and quasicrystals. While Self-Consistent Field Theory (SCFT) offers high accuracy across weak to strong segregation limits, it is computationally intensive, particularly as molecular complexity and simulation box sizes increase. Conversely, Strong Segregation Theory (SST) is computationally efficient but has historically been limited in its application to specific, simple morphologies (often single-core structures) or required complex analytical derivations for each new geometry. There is a need for a flexible, computationally efficient SST framework capable of screening a wide range of compositions and interaction strengths for complex, multi-core ABC star terpolymer morphologies before resorting to more expensive SCFT calculations.

Methodology
The authors present a novel computational implementation of SST based on a tessellation approach using a flexible base motif called "Strongly Segregated Polygons" (SSPs).

The SSP Motif: The central concept is that any 2D morphology can be constructed by assembling SSPs. Each SSP represents a single ABC star core (modeled as a cylinder in 3D, or a point in 2D) surrounded by A, B, and C domains. The SSP is modeled as a polygon (typically six-sided, but extended to eight-sided for lamellar phases) divided into six triangles radiating from the central core node (Node 0) to six outer boundary nodes.
Free Energy Formulation: The total free energy per chain ( $f_c$ $f_{c}$ ) is the sum of stretching energy ( $f_{str}$ $f_{s t r}$ ), interfacial energy ( $f_{int}$ $f_{in t}$ ), and core energy ( $f_{core}$ $f_{cor e}$ ).
- Stretching Energy: Calculated by integrating the stretching free energy of polymer brushes grafted onto a convex surface over the triangular wedges of the SSP. The core radius ( $R_{core}$ ) is not fixed but derived from the area of the SSP ( $A_T$ ) and the degree of polymerization ( $N$ ), introducing a logarithmic dependence on $N$ .
- Interfacial Energy: Determined by the surface tensions between domains (related to Flory interaction parameters $\chi_{IJ}$ ) and the lengths of the interfaces within the SSP.
- Core Energy: Modeled as a mixing energy and stretching penalty within the core region, proportional to the number of monomers in the core.
Optimization Framework: To find the equilibrium morphology, the authors construct a periodic patch of SSPs representing a candidate structure (e.g., [6.6.6], [8.8.4], or multi-core variants). The positions of the "free nodes" (boundary vertices) and the periodicity vectors are varied to minimize the total free energy of the patch. The minimization is performed using the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm.
Phase Diagram Construction: The method systematically scans the ternary composition space ( $\phi_A, \phi_B, \phi_C$ ) for various interaction strengths ( $N\chi$ ). For each point, the free energies of all candidate morphologies (including single-core and multi-core structures) are calculated and compared to identify the global minimum. The authors employ a "growing" strategy, using converged configurations from neighboring compositions as initial guesses to ensure robust minimization.

Key Contributions

Generalized SST Framework: Unlike previous SST implementations limited to single-core structures, this method handles morphologies with multiple non-equivalent core types (e.g., [8.6.4; 8.4.6; 8.6.6]) and complex tilings like the $\Sigma$ -phase (Archimedean tiling 3.3.4.3.4).
Computational Efficiency: The SSP method allows for the rapid calculation of free energies and the construction of entire ternary phase diagrams with significantly less computational cost than SCFT, Monte Carlo, or Dissipative Particle Dynamics.
Multi-Core Analysis: The framework successfully identifies and analyzes the stability of complex multi-core phases that are often obscured in other theoretical treatments or difficult to resolve experimentally.
Extension to Lamellar Phases: The method is extended to include eight-sided SSPs to accurately model lamellar and lamellar+cylinder (L+C) phases, which cannot be represented by simple six-sided polygons.

Results

Symmetric Interactions: For symmetric interaction strengths ( $N\chi_{AB} = N\chi_{BC} = N\chi_{AC} = 60$ ) and $N=1000$ , the phase diagram is dominated by single-core structures ([6.6.6], [8.8.4], [12.6.4], and L+C) in their respective composition regimes. Multi-core structures, such as [8.6.4; 8.4.6; 8.6.6] and the $\Sigma$ -phase, appear at the intersections of these single-core regions.
Effect of $N$ : Varying the degree of polymerization ( $N$ ) alters the balance between stretching and interfacial energies. Increasing $N$ reduces the stability region of the L+C phase, revealing previously masked multi-core structures like [10.6.4; 10.4.6; 10.6.6].
Asymmetric Interactions: When interaction strengths are asymmetric, the phase boundaries shift. The authors observe that higher repulsion between specific monomer pairs favors morphologies where those domains have a higher number of neighbors (approximating a circle to minimize interfacial area). This leads to systematic shifts in the phase diagram, such as the displacement of $\Sigma$ -phases by [12.6.4] phases when B-monomer interactions are strongest.
Structural Deformation: The minimization process reveals that domain interfaces are not straight lines but curve to minimize energy. For instance, in the $\Sigma$ -phase, square tiles contract while triangular tiles expand, and interfaces curve based on local tension differences.

Significance and Claims
The authors position this work as a significant addition to the SST toolbox, bridging the gap between simple analytical models and computationally expensive SCFT. They claim that their method provides a "useful tool" for exploring a wide range of complex morphologies with "relative computational ease."

The paper modestly asserts that while the method is valid in the strong segregation limit, the resulting domain distributions can serve as high-quality initial conditions for more accurate SCFT calculations. The authors emphasize that their approach successfully reproduces the qualitative features of phase diagrams reported in the literature (including SCFT results) while offering a faster route to screening compositions and interaction parameters. They explicitly state that the method is not intended to replace SCFT for high-precision intermediate segregation studies but rather to facilitate the exploration of the vast parameter space of ABC star terpolymers, particularly for complex, multi-core, and potentially quasicrystalline structures.