We report computer simulations of two-dimensional convex hard superellipse particle phases vs. particle shape parameters including aspect ratio, corner curvature, and sidewall curvature. Shapes investigated include disks, ellipses, squares, rectangles, and rhombuses, as well as shapes with non-uniform curvature including rounded squares, rounded rectangles, and rounded rhombuses. Using measures of orientational order, order parameters, and a novel stretched bond orientational order parameter, we systematically identify particle shape properties that determine liquid crystal and crystalline phases including their coarse boundaries and symmetry. We observe phases including isotropic, nematic, tetratic, plastic crystals, square crystals, and hexagonal crystals (including stretched variants). Our results catalog known benchmark shapes, but include new shapes that also interpolate between known shapes. Our results indicate design rules for particle shapes that determine two-dimensional liquid, liquid crystalline, and crystalline microstructures that can be realized via particle assembly.
more »
« less
Anisotropic particle multiphase equilibria in nonuniform fields
We report a method to predict equilibrium concentration profiles of hard ellipses in nonuniform fields, including multiphase equilibria of fluid, nematic, and crystal phases. Our model is based on a balance of osmotic pressure and field mediated forces by employing the local density approximation. Implementation of this model requires development of accurate equations of state for each phase as a function of hard ellipse aspect ratio in the range k = 1–9. The predicted density profiles display overall good agreement with Monte Carlo simulations for hard ellipse aspect ratios k = 2, 4, and 6 in gravitational and electric fields with fluid–nematic, fluid–crystal, and fluid–nematic–crystal multiphase equilibria. The profiles of local order parameters for positional and orientational order display good agreement with values expected for bulk homogeneous hard ellipses in the same density ranges. Small discrepancies between predictions and simulations are observed at crystal–nematic and crystal–fluid interfaces due to limitations of the local density approximation, finite system sizes, and uniform periodic boundary conditions. The ability of the model to capture multiphase equilibria of hard ellipses in nonuniform fields as a function of particle aspect ratio provides a basis to control anisotropic particle microstructure on interfacial energy landscapes in diverse materials and applications.
more »
« less
- Award ID(s):
- 2113594
- PAR ID:
- 10504143
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 159
- Issue:
- 12
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We generate and study dense positionally and/or orientationally disordered, including jammed, monodisperse packings of hard convex lens-shaped particles (lenses). Relatively dense isotropic fluid configurations of lenses of various aspect ratios are slowly compressed via a Monte Carlo method based procedure. Under this compression protocol, while ‘flat’ lenses form a nematic fluid phase (where particles are positionally disordered but orientationally ordered) and ‘globular’ lenses form a plastic solid phase (where particles are positionally ordered but orientationally disordered), ‘intermediate’, neither ‘flat’ nor ‘globular’, lenses do not form either mesophase. In general, a crystal solid phase (where particles are both positionally and orientationally ordered) does not spontaneously form during lengthy numerical simulation runs. In correspondence to those volume fractions at which a transition to the crystal solid phase would occur in equilibrium, a ‘downturn’ is observed in the inverse compressibility factor versus volume fraction curve beyond which this curve behaves essentially linearly. This allows us to estimate the volume fraction at jamming of the dense non-crystalline packings so generated. These packings are nematic for ‘flat’ lenses and plastic for ‘globular’ lenses, while they are robustly isotropic for ‘intermediate’ lenses, as confirmed by the calculation of the τ order metric, among other quantities. The structure factors S ( k ) of the corresponding jammed states tend to zero as the wavenumber k goes to zero, indicating they are effectively hyperuniform ( i.e. , their infinite-wavelength density fluctuations are anomalously suppressed). Among all possible lens shapes, ‘intermediate’ lenses with aspect ratio around 2/3 are special because they are those that reach the highest volume fractions at jamming while being positionally and orientationally disordered and these volume fractions are as high as those reached by nematic jammed states of ‘flat’ lenses and plastic jammed states of ‘globular’ lenses. All of their attributes, taken together, make such ‘intermediate’ lens packings particularly good glass-forming materials.more » « less
-
Assembling different shaped particles into ordered microstructures is an open challenge in creating multifunctional particle-based materials and devices. Here, we report the two-dimensional (2D) AC electric field mediated assembly of different shaped colloidal particles into amorphous, liquid crystalline, and crystalline microstructures. Particle shapes investigated include disks, ellipses, squares, and rectangles, which show how systematic variations in anisotropy and corner curvature determine the number and type of resulting microstructures. AC electric fields induce dipolar interactions to control particle positional and orientational order. Microstructural states are determined via particle tracking to compute order parameters, which agree with computer simulations and show how particle packing and dipolar interactions together produce each structure. Results demonstrate how choice of particle shape and field conditions enable kinetically viable routes to assemble nematic, tetratic, and smectic liquid crystal structures as well as crystals with stretched 4- and 6-fold symmetry. Results show it is possible to assemble all corresponding hard particle phases, but also show how dipolar interactions influence and produce additional microstructures. Our findings provide design rules for the assembly of diverse microstructures of different shaped particles in AC electric fields, which could enable scalable and reconfigurable particle-based materials, displays, and printing technologies.more » « less
-
Monte Carlo simulations were used to study the influence of particle aspect ratio on the kinetics and phase behavior of hard gyrobifastigia (GBF). First, the formation of a highly anisotropic nucleus shape in the isotropic-to-crystal transition in regular GBF is explained by the differences in interfacial free energies of various crystal planes and the nucleus geometry predicted by the Wulff construction. GBF-related shapes with various aspect ratios were then studied, mapping their equations of state, determining phase coexistence conditions via interfacial pinning, and computing nucleation free-energy barriers via umbrella sampling using suitable order parameters. Our simulations reveal a reduction of the kinetic barrier for isotropic–crystal transition upon an increase in aspect ratio, and that for highly oblate and prolate aspect ratios, an intermediate nematic phase is stabilized. Our results and observations also support two conjectures for the formation of the crystalline state from the isotropic phase: that low phase free energies at the ordering phase transition correlate with low transition barriers and that the emergence of a mesophase provides a steppingstone that expedites crystallization.more » « less
-
Lenarčič, Jadran; Husty, Manfred (Ed.)The multidirectional transmission characteristics of a five-bar linkage can be visualized by plotting Jacobian-defined velocity ellipses inside its workspace. The orientation, size, and aspect ratio of these ellipses indicate directional force and velocity multiplication from the actuators to the end-effector. Our broader goal is approximate dimensional synthesis to achieve desired ellipses. On a workspace bound, the minor axis of a velocity ellipse collapses while the major axis aligns tangential to the bound. Interior to the workspace, ellipses vary with continuity. Therefore, the shape of a workspace bound influences the interior ellipses. The workspace bounds of a five-bar linkage are formed from segment of four-bar coupler curves (the locus of endpoint positions while the five-bar is held in output singularity conditions) and circular segments. Therefore, interior ellipses can be influenced by the path synthesis of four-bar linkages that represent the five-bar situated with certain links held colinear (the output singularity conditions). This paper details the synthesis of these four-bar coupler curves for forming the workspace bounds of a five-bar in order to influence its interior ellipses. Our approach employs saddle graphs that detail the connectivity of critical points over an optimization function.more » « less