Phase Behavior of Colloidal Superballs: Shape Interpolation from Spheres to Cubes

TL;DR

This study explores how the shape of colloidal superballs influences their phase behavior, revealing shape-dependent ordering, phase transitions, and the impact of simulation artifacts, with implications for experimental synthesis.

Contribution

It provides a detailed analysis of the phase behavior of convex superballs with shape interpolation from spheres to cubes, highlighting the role of asphericity and identifying phase transition characteristics.

Findings

Cubatic order depends on shape parameter q.

Phase transition from crystal to reduced order occurs for 1<q<3.

Long-range order persists for q≥3 until melting.

Abstract

The phase behavior of hard superballs is examined using molecular dynamics within a deformable periodic simulation box. A superball's interior is defined by the inequality $|x|^{2q} + |y|^{2q} + |z|^{2q} \leq 1$, which provides a versatile family of convex particles ($q \geq 0.5$) with cube-like and octahedron-like shapes as well as concave particles ($q < 0.5$) with octahedron-like shapes. Here, we consider the convex case with a deformation parameter q between the sphere point (q = 1) and the cube (q = 1). We find that the asphericity plays a significant role in the extent of cubatic ordering of both the liquid and crystal phases. Calculation of the first few virial coefficients shows that superballs that are visually similar to cubes can have low-density equations of state closer to spheres than to cubes. Dense liquids of superballs display cubatic orientational order that extends over several particle lengths only for large q. Along the ordered, high-density equation of state, superballs with 1 < q < 3 exhibit clear evidence of a phase transition from a crystal state to a state with reduced long-ranged orientational order upon the reduction of density. For $q \geq 3$, long-ranged orientational order persists until the melting transition. The width of coexistence region between the liquid and ordered, high-density phase decreases with q up to q = 4.0. The structures of the high-density phases are examined using certain order parameters, distribution functions, and orientational correlation functions. We also find that a fixed simulation cell induces artificial phase transitions that are out of equilibrium. Current fabrication techniques allow for the synthesis of colloidal superballs, and thus the phase behavior of such systems can be investigated experimentally.