Search

Among lattice configurations of densely packed hard ellipses, Monte Carlo simulations are used to identify the so-called parallel and diagonal lattices as the two favorable states. The free energies of these two states are computed for several system sizes employing the Einstein crystal method. An accurate calculation of the free energy difference between the two states reveals the parallel lattice as the state with the lowest free energy. The origin of the entropic difference between the two states is further elucidated by assessing the roles of the translational and rotational degrees of freedom. [ABSTRACT FROM AUTHOR]

We investigate the two-dimensional behavior of colloidal patchy ellipsoids specifically designed to follow a two-step assembly process from the monomer state to mesoscopic liquid-crystal phases via the formation of the so-called bent-core units at the intermediate stage. Our model comprises a binary mixture of ellipses interacting via the Gay-Berne potential and decorated by surface patches, with the binary components being mirror-image variants of each other-referred to as left-handed and right-handed ellipses according to the position of their patches. The surface patches are designed so as in the first stage of the assembly the monomers form bent-cores units, i.e., V-shaped dimers with a specific bent angle. The Gay-Berne interactions, which act between the ellipses, drive the dimers to subsequently form the characteristic phase observed in bent-core liquid crystals. We numerically investigate-by means of both Molecular Dynamics and Monte Carlo simulations-the described two-step process: we first optimize a target bent-core unit and then fully characterize its state diagram in temperature and density, defining the regions where the different liquid crystalline phases dominate., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We employ classical density-functional theory to investigate the phase diagram of an assembly of mutually penetrable, parallel ellipsoids interacting via the generalized exponential model of index four (GEM-4) pair potential. We show that the crystal phases of the system are obtained from those of the spherically symmetric GEM-4 model by rescaling the lattice vectors. Performing this rescaling in combination with an arbitrary rotation of the lattice leads to infinitely many different structures with the same free energy, thereby implying their infinite degeneracy. These findings generalize to non-zero temperature the results formerly obtained by us [Pini et al., J. Chem. Phys. 153, 164901 (2020)] for the ground state of a similar system of ellipsoids interacting via a Gaussian potential. According to the mean-field free-energy functional used here, our conclusions apply to soft-core potentials both when they form cluster crystals as the GEM-4 and when they form single-occupancy crystals as the Gaussian itself., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

We investigate the phase ordering (pattern formation) of systems of two-dimensional core-shell particles using Monte Carlo (MC) computer simulations and classical density functional theory (DFT). The particles interact via a pair potential having a hard core and a repulsive square shoulder. Our simulations show that on cooling, the liquid state structure becomes increasingly characterized by long wavelength density modulations and on further cooling forms a variety of other phases, including clustered, striped, and other patterned phases. In DFT, the hard core part of the potential is treated using either fundamental measure theory or a simple local density approximation, whereas the soft shoulder is treated using the random phase approximation. The different DFTs are benchmarked using large-scale grand-canonical-MC and Gibbs-ensemble-MC simulations, demonstrating their predictive capabilities and shortcomings. We find that having the liquid state static structure factor S(k) for wavenumber k is sufficient to identify the Fourier modes governing both the liquid and solid phases. This allows us to identify from easier-to-obtain liquid state data the wavenumbers relevant to the periodic phases and to predict roughly where in the phase diagram these patterned phases arise., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Phase separation plays a key role in determining the self-assembly of biological and soft-matter systems. In biological systems, liquid-liquid phase separation inside a cell leads to the formation of various macromolecular aggregates. The interaction among these aggregates is soft, i.e., they can significantly overlap at a small energy cost. From a computer simulation point of view, these complex macromolecular aggregates are generally modeled by soft particles. The effective interaction between two particles is defined via the generalized exponential model of index n, with n = 4. Here, using molecular dynamics simulations, we study the phase separation dynamics of a size-symmetric binary mixture of ultrasoft particles. We find that when the mixture is quenched to a temperature below the critical temperature, the two components spontaneously start to separate. Domains of the two components form, and the equal-time order parameter reveals that the domain sizes grow with time in a power-law manner with an exponent of 1/3, which is consistent with the Lifshitz-Slyozov law for conserved systems. Furthermore, the static structure factor shows a power-law decay with an exponent of 4, consistent with the Porod law., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)