Dilute gel networks vs. clumpy gels in colloidal systems with a competition between repulsive and attractive interactions.

Using Brownian dynamics simulations we study gel-forming colloidal systems. The focus of this article lies on the differences of dense and dilute gel networks in terms of structure formation both on a local and a global level. We apply reduction algorithms and observe that dilute networks and dense gels differ in the way structural properties like the thickness of strands emerge. We also analyze the percolation behavior and find that two different regimes of percolation exist which might be responsible for structural differences. In dilute networks we confirm that solidity is mainly a consequence of pentagonal bipyramids forming in the network. In dense gels, tetrahedral structures also influence solidity.

Structural and dynamical properties of dilute gel networks in colloid-polymer mixtures.

The competition of short-ranged depletion attraction and long-ranged repulsion between colloidal particles in colloid-polymer mixtures leads to the formation of heterogeneous gel-like structures. Our special focus will be on the states where the colloids arrange in thin strands that span the whole system and that we will refer to as dilute gel networks. These states occur at low packing fractions for attractions that are stronger than those at both the binodal line of the equilibrium gas-liquid phase separation and the directed percolation transition line. By using Brownian dynamics simulations, we explore the formation, structure, and aging dynamics of dilute gel networks. The essential connections in a dilute gel network are determined by constructing reduced networks. We compare the observed properties to those of clumpy gels or cluster fluids. Our results demonstrate that both the structure and the (often slow) dynamics of the stable or meta-stable heterogeneous states in colloid-polymer mixtures possess distinct features on various length and time scales and thus are richly diverse.

Event-chain Monte Carlo simulations of the liquid to solid transition of two-dimensional decagonal colloidal quasicrystals.

In event-chain Monte Carlo simulations, we model colloidal particles in two dimensions that interact according to an isotropic short-ranged pair potential which supports the two typical length scales present in decagonal quasicrystals. We investigate the assembled structures as we vary the density and temperature. Our special interest is related to the transition from quasicrystal to liquid. In contrast to the KTHNY melting theory for quasicrystals which predicts an intermediate pentahedratic phase, we find a one-step first-order melting transition. However, we discover that the slow relaxation of phasonic flips, i.e. rearrangements of the particles due to additional degrees of freedom in quasicrystals, changes the positional correlation functions, to the extent that structures with long-range orientational correlations, but exponentially decaying positional correlations, are observed.

Growth of two-dimensional decagonal colloidal quasicrystals.

The growth of quasicrystals, i.e. structures with long-range positional order but no periodic translational symmetry, is more complex than the growth of periodic crystals. By employing Brownian dynamics simulations in two dimensions for colloidal particles that interact according to an isotropic pair potential with two incommensurate lengths, we study the growth of quasicrystalline structures by sequentially depositing particles at their surface. We quantify the occurrence of quasicrystalline order as a function of the temperature and the rate of added particles. In addition, we explore defects like local triangular order or gaps within the quasicrystalline structure. Furthermore, we analyze the shapes of the surfaces in grown structures which tend to build straight lines along the symmetry axes of the quasicrystal. Finally, we identify phasonic flips which are rearrangements of the particles due to additional degrees of freedom. The number of phasonic flips decreases with the distance to the surface.

Dislocation-free growth of quasicrystals from two seeds due to additional phasonic degrees of freedom.

We explore the growth of two-dimensional quasicrystals, i.e., aperiodic structures that possess long-range order, from two seeds at various distances and with different orientations by using dynamical phase-field crystal calculations. We compare the results to the growth of periodic crystals from two seeds. There, a domain border consisting of dislocations is observed in case of large distances between the seed and large angles between their orientation. Furthermore, a domain border is found if the seeds are placed at a distance that does not fit to the periodic lattice. In the case of the growth of quasicrystals, we only observe domain borders for large distances and different orientations. Note that all distances do inherently not match to a perfect domain wall-free quasicrystalline structure. Nevertheless, we find dislocation-free growth for all seeds at a small enough distance and for all seeds that approximately have the same orientation. In periodic structures, the stress that occurs due to incommensurate distances between the seeds results in phononic strain fields or, in the case of too large stresses, in dislocations. In contrast, in quasicrystals an additional phasonic strain field can occur and suppress dislocations. Phasons are additional degrees of freedom that are unique to quasicrystals. As a consequence, the additional phasonic strain field helps to distribute the stress and facilitates the growth of dislocation-free quasicrystals from multiple seeds. In contrast, in the periodic case the growth from multiple seeds most likely leads to a structure with multiple domains. Our work lays the theoretical foundations for growing perfect quasicrystals from different seeds and is therefore relevant for many applications.

Detection of phonon and phason modes in intrinsic colloidal quasicrystals by reconstructing their structure in hyperspace.

Phasons are additional degrees of freedom which occur in quasicrystals alongside the phonons known from conventional periodic crystals. The rearrangements of particles that are associated with a phason mode are hard to interpret in physical space. We reconstruct the quasicrystal structure by an embedding into extended higher-dimensional space, where phasons correspond to displacements perpendicular to the physical space. In dislocation-free decagonal colloidal quasicrystals annealed with Brownian dynamics simulations, we identify thermal phonon and phason modes. Finite phononic strain is pinned by phasonic excitations even after cooling down to zero temperature. For the phasonic displacements underlying the flip pattern, the reconstruction method gives an approximation within the limits of a multi-mode harmonic ansatz, and points to fundamental limitations of a harmonic picture for phasonic excitations in intrinsic colloidal quasicrystals.

Directed percolation identified as equilibrium pre-transition towards non-equilibrium arrested gel states.

The macroscopic properties of gels arise from their slow dynamics and load-bearing network structure, which are exploited by nature and in numerous industrial products. However, a link between these structural and dynamical properties has remained elusive. Here we present confocal microscopy experiments and simulations of gel-forming colloid-polymer mixtures. They reveal that gel formation is preceded by continuous and directed percolation. Both transitions lead to system-spanning networks, but only directed percolation results in extremely slow dynamics, ageing and a shrinking of the gel that resembles synaeresis. Therefore, dynamical arrest in gels is found to be linked to a structural transition, namely directed percolation, which is quantitatively associated with the mean number of bonded neighbours. Directed percolation denotes a universality class of transitions. Our study hence connects gel formation to a well-developed theoretical framework, which now can be exploited to achieve a detailed understanding of arrested gels.

Asymptotic densities of ballistic Lévy walks.

We propose an analytical method to determine the shape of density profiles in the asymptotic long-time limit for a broad class of coupled continuous-time random walks which operate in the ballistic regime. In particular, we show that different scenarios of performing a random-walk step, via making an instantaneous jump penalized by a proper waiting time or via moving with a constant speed, dramatically effect the corresponding propagators, despite the fact that the end points of the steps are identical. Furthermore, if the speed during each step of the random walk is itself a random variable, its distribution gets clearly reflected in the asymptotic density of random walkers. These features are in contrast with more standard nonballistic random walks.

Growth modes of quasicrystals.

The growth of quasicrystals, i.e., aperiodic structures with long-range order, seeded from the melt is investigated using a dynamical phase field crystal model. Depending on the thermodynamic conditions, two different growth modes are detected, namely defect-free growth of the stable quasicrystal and a mode dominated by phasonic flips which are incorporated as local defects into the grown structure such that random tilinglike ordering emerges. The latter growth mode is unique to quasicrystals and can be verified in experiments on one-component mesoscopic systems.

Crystallization induced by multiple seeds: dynamical density functional approach.

Using microscopic dynamical density functional theory, we calculate the dynamical formation of polycrystals by following the crystal growth around multiple crystalline seeds imposed to an undercooled fluid. Depending on the undercooling and the size ratio as well as the relative crystal orientation of two neighboring seeds, three possibilities of the final state emerge, namely no crystallization at all, formation of a monocrystal, or two crystallites separated by a curved grain boundary. Our results, which are obtained for two-dimensional hard disk systems using a fundamental-measure density functional, shed new light on the particle-resolved structure and growth of polycrystalline material in general.

Rhombic preordering on a square substrate.

A competition of incommensurate symmetries occurs whenever a system is forced to conform to an ordering that is different from the intrinsically preferred structure of the system itself. As a model system of such a competition, we study the rivalry between the triangular ordering of hard disks and the square symmetry induced by a periodic square substrate. By using density functional theory as well as Monte Carlo computer simulations, we determine the full phase behavior for the case of one particle per minimum. We observe a rhombic preordering structure preceding the hexagonal solid as a direct consequence of the competing symmetries. Furthermore, the square-rhombic transition is reentrant with increasing substrate interaction. Our predictions can be verified in experiments of colloids in laser fields.

A colloidal model system with tunable disorder: solid-fluid transition and discontinuities in the limit of zero disorder.

We study a colloidal model system where disorder can be continuously tuned from no disorder--corresponding to a system that can crystallize--to large disorder where geometrical frustration occurs. The model system consists of colloidal particles with screened electrostatic repulsion. They can only move on single lines which are parallel and equidistant to each other. We introduce disorder by modulating the particle line density. The system exhibits a solid-to-fluid transition which we study by the structure factor and the temporal evolution of the mean-square distance of nearest neighbors on neighboring lines. A determining feature is the occurrence of discontinuities when disorder is tuned to zero. We observe that the peak height of the pair correlation function in the solid phase does not extrapolate to the value of the perfect crystal. Similarly, the mean interaction energy and the screening length at which the solid-fluid transition occurs seem to be discontinuous when the limit of zero disorder is approached.

Archimedean-like colloidal tilings on substrates with decagonal and tetradecagonal symmetry.

Two-dimensional colloidal suspensions subjected to laser interference patterns with decagonal symmetry can form an Archimedean-like tiling phase where rows of squares and triangles order aperiodically along one direction (J. Mikhael et al., Nature 454, 501 (2008)). In experiments as well as in Monte Carlo and Brownian dynamics simulations, we identify a similar phase when the laser field possesses tetradecagonal symmetry. We characterize the structure of both Archimedean-like tilings in detail and point out how the tilings differ from each other. Furthermore, we also estimate specific particle densities where the Archimedean-like tiling phases occur. Finally, using Brownian dynamics simulations we demonstrate how phasonic distortions of the decagonal laser field influence the Archimedean-like tiling. In particular, the domain size of the tiling can be enlarged by phasonic drifts and constant gradients in the phasonic displacement. We demonstrate that the latter occurs when the interfering laser beams are not ideally adjusted.

Brownian particles in random and quasicrystalline potentials: how they approach the equilibrium.

We study the Brownian motion of an ensemble of single colloidal particles in a random square and a quasicrystalline potential when they start from non-equlibrium. For both potentials, Brownian dynamics simulations reveal a widespread subdiffusive regime before the diffusive long-time limit is reached in thermal equilibrium. We develop a random trap model based on a distribution for the depths of trapping sites that reproduces the results of the simulations in detail. Especially, it gives analytic formulas for the long-time diffusion constant and the relaxation time into the diffusive regime. Aside from detailed differences, our work demonstrates that quasicrystalline potentials can be used to mimic aspects of random potentials.

Photon channelling in foams.

Experiments by Gittings, Bandyopadhyay and Durian (Europhys. Lett. 65, 414 (2004)) demonstrate that light possesses a higher probability to propagate in the liquid phase of a foam due to total reflection. The authors term this observation photon channelling which we investigate in this article theoretically. We first derive a central relation in the work of Gitting et al. without any free parameters. It links the photon's path-length fraction f in the liquid phase to the liquid fraction epsilon. We then construct two-dimensional Voronoi foams, replace the cell edges by channels to represent the liquid films and simulate photon paths according to the laws of ray optics using transmission and reflection coefficients from Fresnel's formulas. In an exact honeycomb foam, the photons show superdiffusive behavior. It becomes diffusive as soon as disorder is introduced into the foams. The dependence of the diffusion constant on channel width and refractive index is explained by a one-dimensional random-walk model. It contains a photon channelling state that is crucial for the understanding of the numerical results. At the end, we shortly comment on the observation that photon channelling only occurs in a finite range of epsilon.