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
Effect of particle anisotropy on the thermodynamics and kinetics of ordering transitions in hard faceted particles
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
- Award ID(s):
- 1907369
- PAR ID:
- 10426412
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 4
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 044502
- 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
-
We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems.more » « less
-
Effect of Filler Aspect Ratio on Stiffness and Conductivity in Phase‐Changing Particulate CompositesAbstract Variable stiffness in elastomers can be achieved through the introduction of low melting point alloy particles, such as Field's metal (FM), enabling on‐demand switchable elasticity and anisotropy in response to thermal stimulus. Because the FM particles are thermally transitioned between solid and liquid phases, it is beneficial for the composite to be electrically conductive so the stiffness may be controlled via direct Joule heating. While FM is highly conductive, spherical particles contribute to a high percolation threshold. In this paper, it is shown that the percolation threshold of FM particulate composites can be reduced with increasing particles aspect ratio. Increasing the aspect ratio of phase‐changing fillers also increases the rigid‐to‐soft modulus ratio of the composite by raising the elastic modulus in the rigid state while preserving the low modulus in the soft state. The results indicate that lower quantities of high aspect ratio FM particles can be used to achieve both electrical conductivity and stiffness‐switching via a single solution and without introducing additional conductive fillers. This technique is applied to enable a highly stretchable, variable stiffness, and electrically conductive composite, which, when patterned around an inflatable actuator, allows for adaptable trajectories via selective softening of the surface materials.more » « less
-
Liquid crystal elastomers (LCEs) are made of liquid crystal molecules integrated with rubber-like polymer networks. An LCE exhibits both the thermotropic property of liquid crystals and the large deformation of elastomers. It can be monodomain or polydomain in the nematic phase and transforms to an isotropic phase at elevated temperature. These features have enabled various new applications of LCEs in robotics and other fields. However, despite substantial research and development in recent years, thermomechanical coupling in polydomain LCEs remains poorly studied, such as their temperature-dependent mechanical response and stretch-influenced isotropic-nematic phase transition. This knowledge gap severely limits the fundamental understanding of the structure-property relationship, as well as future developments of LCEs with precisely controlled material behaviors. Here, we construct a theoretical model to investigate the thermomechanical coupling in polydomain LCEs. The model includes a quasi-convex elastic energy of the polymer network and a free energy of mesogens. We study the working conditions where a polydomain LCE is subjected to various prescribed planar stretches and temperatures. The quasi-convex elastic energy enables a “mechanical phase diagram” that describes the macroscopic effective mechanical response of the material, and the free energy of mesogens governs their first-order nematic-isotropic phase transition. The evolution of the mechanical phase diagram and the order parameter with temperature is predicted and discussed. Unique temperature-dependent mechanical behaviors of the polydomain LCE that have never been reported before are shown in their stress-stretch curves. These results are hoped to motivate future fundamental studies and new applications of thermomechanical LCEs.more » « less