Phase separation patterns for diblock copolymers on spherical surfaces: a finite volume method

We explore phase separation on spherical surfaces by solving the Cahn-Hilliard equation modified for diblock copolymers using a finite volume method. The spherical surface is discretized into almost uniform triangles by employing successive dyadic refinements of the spherical icosahedron, a methodology that avoids potential mathematical and numerical problems related to the poles in spherical coordinates. The finite volume method is based on averaging Voronoi cells built from triangular meshes to calculate the Laplace-Beltrami operator on the curved surface, which greatly improves both the accuracy and speed of calculation as compared to the conventional finite difference method. By using this method we simulate the phase separation of diblock copolymers on a spherical surface. It is found that stable and intrinsic defects, which would not occur in a flat space after sufficient annealing, appear in the periodic arrangement of the domains on the curved surface due to the distinct Euler characteristic of the surface.