Mesophases of soft-sphere aggregates (original) (raw)
Related papers
Soft spheres make more mesophases
Europhysics Letters (EPL), 2007
We use both mean-field methods and numerical simulation to study the phase diagram of classical particles interacting with a hard-core and repulsive, soft shoulder. Despite the purely repulsive interaction, this system displays a remarkable array of aggregate phases arising from the competition between the hard-core and shoulder length scales. In the limit of large shoulder width to core size, we argue that this phase diagram has a number of universal features, and classify the set of repulsive shoulders that lead to aggregation at high density. Surprisingly, the phase sequence and aggregate size adjusts so as to keep almost constant inter-aggregate separation.
From compact to fractal crystalline clusters in concentrated systems of monodisperse hard spheres
Soft Matter, 2012
We address the crystallization of monodisperse hard spheres in terms of the properties of finitesize crystalline clusters. By means of large scale event-driven Molecular Dynamics simulations, we study systems at different packing fractions φ ranging from weakly supersaturated state points to glassy ones, covering different nucleation regimes. We find that such regimes also result in different properties of the crystalline clusters: compact clusters are formed in the classical-nucleation-theory regime (φ ≤ 0.54), while a crossover to fractal, ramified clusters is encountered upon increasing packing fraction (φ ≥ 0.56), where nucleation is more spinodal-like. We draw an analogy between macroscopic crystallization of our clusters and percolation of attractive systems to provide ideas on how the packing fraction influences the final structure of the macroscopic crystals. In our previous work (Phys. Rev. Lett., 106, 215701, 2011), we have demonstrated how crystallization from a glass (at φ > 0.58) happens via a gradual (many-step) mechanism: in this paper we show how the mechanism of gradual growth seems to hold also in super-saturated systems just above freezing showing that static properties of clusters are not much affected by dynamics.
Mesoscopic theory for soft-matter systems
2009
A mesoscopic description of various systems containing spherical charged particles in solvent inducing effective attraction is developed by a systematic coarse-graining procedure. For weak ordering the theory can be reduced to the Landau-Ginzburg or the Landau-Brazovskii field theory, depending on the form of the effective interactions between particles. Within the framework of this theory we obtain and discuss the A-line and the universal sequence of phases: disordered, bcc, hexagonal, lamellar, inverted hexagonal, inverted bcc, disordered, for increasing density well below the close-packing density. The sequence of phases agrees with experimental observations and with simulations of many self-assembling systems. In addition to the above phases, more complex phases may appear depending on the interaction potentials. For a particular form of the short-range attraction long-range repulsion potential we find the bicontinuous gyroid phase (Ia3d symmetry) that may be related to a network forming cluster of colloids in a mixture of colloids and nonadsorbing polymers.
Free cooling phase-diagram of hard-spheres with short- and long-range interactions
The European Physical Journal Special Topics, 2014
We study the stability, the clustering and the phase-diagram of free cooling granular gases. The systems consist of mono-disperse particles with additional non-contact (long-range) interactions, and are simulated here by the event-driven molecular dynamics algorithm with discrete (short-range shoulders or wells) potentials (in both 2D and 3D). Astonishingly good agreement is found with a mean field theory, where only the energy dissipation term is modified to account for both repulsive or attractive non-contact interactions. Attractive potentials enhance cooling and structure formation (clustering), whereas repulsive potentials reduce it, as intuition suggests. The system evolution is controlled by a single parameter: the non-contact potential strength scaled by the fluctuation kinetic energy (granular temperature). When this is small, as expected, the classical homogeneous cooling state is found. However, if the effective dissipation is strong enough, structure formation proceeds, before (in the repulsive case) non-contact forces get strong enough to undo the clustering (due to the ongoing dissipation of granular temperature). For both repulsive and attractive potentials, in the homogeneous regime, the cooling shows a universal behaviour when the (inverse) control parameter is used as evolution variable instead of time. The transition to a non-homogeneous regime, as predicted by stability analysis, is affected by both dissipation and potential strength. This can be cast into a phase diagram where the system changes with time, which leaves open many challenges for future research.
Crystallization and dynamical arrest of attractive hard spheres
The Journal of Chemical Physics, 2009
Crystallization of hard spheres interacting with a square well potential was investigated by numerical simulations using so-called Brownian cluster dynamics. The phase diagram was determined over a broad range of volume fractions. The crystallization rate was studied as a function of the interaction strength expressed in terms of the second virial coefficient. For volume fractions below about 0.3 the rate was found to increase abruptly with increasing attraction at the binodal of the metastable liquid-liquid phase separation. The rate increased until a maximum was reached after which it decreased with a power law dependence on the second virial coefficient. Above a critical percolation concentration, a transient system spanning network of connected particles was formed. Crystals were formed initially as part of the network, but eventually crystallization led to the breakup of the network. The lifetime of the transient gels increased very rapidly over a small range of interaction energies. Weak attraction destabilized the so-called repulsive crystals formed in pure hard sphere systems and shifted the coexistence line to higher volume fractions. Stronger attraction led to the formation of a denser, so-called attractive, crystalline phase. Nucleation of attractive crystals in the repulsive crystalline phase was observed close to the transition.
From sticky-hard-sphere to Lennard-Jones-type clusters
Physical review. E, 2018
A relation M_{SHS→LJ} between the set of nonisomorphic sticky-hard-sphere clusters M_{SHS} and the sets of local energy minima M_{LJ} of the (m,n)-Lennard-Jones potential V_{mn}^{LJ}(r)=ɛ/n-m[mr^{-n}-nr^{-m}] is established. The number of nonisomorphic stable clusters depends strongly and nontrivially on both m and n and increases exponentially with increasing cluster size N for N≳10. While the map from M_{SHS}→M_{SHS→LJ} is noninjective and nonsurjective, the number of Lennard-Jones structures missing from the map is relatively small for cluster sizes up to N=13, and most of the missing structures correspond to energetically unfavorable minima even for fairly low (m,n). Furthermore, even the softest Lennard-Jones potential predicts that the coordination of 13 spheres around a central sphere is problematic (the Gregory-Newton problem). A more realistic extended Lennard-Jones potential chosen from coupled-cluster calculations for a rare gas dimer leads to a substantial increase in th...
Phase behaviour and gravity-directed self assembly of hard convex spherical caps
Soft matter, 2017
We investigate the phase behaviour and self-assembly of convex spherical caps using Monte Carlo simulations. This model is used to represent the main features observed in experimental colloidal particles with mushroom-cap shape [Riley et al., Langmuir, 2010, 26, 1648]. The geometry of this non-centrosymmetric convex model is fully characterized by the aspect ratio χ* defined as the spherical cap height to diameter ratio. We use NPT Monte Carlo simulations combined with free energy calculations to determine the most stable crystal structures and the phase behaviour of convex spherical caps with different aspect ratios. We find a variety of crystal structures at each aspect ratio, including plastic and dimer-based crystals; small differences in chemical potential between the structures with similar morphology suggest that convex spherical caps have the tendency to form polycrystalline phases rather than crystallising into a single uniform structure. With the exception of plastic cryst...
The Journal of Chemical Physics, 2006
Reversible aggregation of spheres is simulated using a novel method in which clusters of bound spheres diffuse collectively with a diffusion coefficient proportional to their radius. It is shown that the equilibrium state is the same as with other simulation techniques, but with the present method more realistic kinetics are obtained. The behavior as a function of volume fraction and interaction strength was tested for two different attraction ranges. The binodal and the percolation threshold were determined. The cluster structure and size distribution close to the percolation threshold were found to be consistent with the percolation model. Close to the binodal phase separation occurred through the growth of spherical dense domains, while for deep quenches a system spanning network is formed that coarsens with a rate that decreases with increasing attraction. We found no indication for arrest of the coarsening.