Speeding up Dynamics by Tuning the Noncommensurate Size of Rodlike Particles in a Smectic Phase

Curvature-directed anchoring and defect structure of colloidal smectic liquid crystals in confinement

Rod-like objects at high packing fractions can exhibit liquid crystalline ordering. By controlling how the rods align near a boundary, i.e. the anchoring, the defects of a liquid crystal can be selected and tuned. For smectic phases, the rods break rotational and translational symmetry by forming lamellae. Smectic defects thereby include both discontinuities in the rod orientational order (disclinations), as well as in the positional order (dislocations). In this work, we use experiments and simulations to uncover the geometrical conditions necessary for a boundary to set the anchoring of a confined, particle-resolved, smectic liquid crystal. We confine a colloidal smectic within elliptical wells of varying size and shape for a smooth variation of the boundary curvature. We find that the anchoring depends upon the local boundary curvature, with an anchoring transition observed at a critical radius of curvature approximately twice the rod length. Surprisingly, the critical radius of curvature for an anchoring transition holds across a wide range of rod lengths and packing fractions. The anchoring controls the defect structure. By analyzing topological charges and networks composed of maximum density (rod centers) and minimum density (rod ends), we quantify disclinations and dislocations formed with varying confinement geometry. Circular confinements, characterized by planar anchoring, promote disclinations, whereas elliptical confinements, featuring antipodal regions of homeotropic anchoring, promote long-range smectic order and dislocations. Our findings demonstrate how geometrical constraints can control the anchoring and defect structures of liquid crystals-a principle that is applicable from molecular to colloidal length scales.

Preferential ordering of incommensurate-length guest particles in a smectic host

Using density functional theory, we study the preferential ordering of rod-like guest particles immersed in a smectic host fluid. Within a model of perfectly aligned rods and assuming that the guest particles do not perturb the smectic host fluid, simple excluded-volume arguments explain that guest particles that are comparable in length to the host particles order in phase with the smectic host density layering, whereas guest particles that are considerably shorter or longer order in antiphase. The corresponding free-energy minima are separated by energetic barriers on the order of the thermal energy kBT, suggesting that guest particles undergo hopping-type diffusion between adjacent smectic layers. Upon introducing a slight orientational mismatch between the guest particles and the perfectly aligned smectic host, an additional, smaller free-energy barrier emerges for a range of intermediate guest-to-host length ratios, which splits the free-energy minimum into two. Guest particles in this range occupy positions intermediate between in-phase ordering and in-antiphase ordering. Finally, we use Kramers' theory to identify slow, fast, and intermediate diffusive regimes for the guest particles as a function of their length. Our model is in qualitative agreement with experiment and simulation and provides an alternative, complementary explanation in terms of a free-energy landscape for the intermediate diffusive regime, which was previously hypothesized to result from temporary caging effects [M. Chiappini, E. Grelet, and M. Dijkstra, Phys. Rev. Lett. 124, 087801 (2020)]. We argue that our simple model of aligned rods captures the salient features of incommensurate-length guest particles in a smectic host if a slight orientational mismatch is introduced.

Insight into the Viscoelasticity of Self-Assembling Smectic Liquid Crystals of Colloidal Rods from Active Microrheology Simulations

The rheology of colloidal suspensions is of utmost importance in a wide variety of interdisciplinary applications in formulation technology, determining equally interesting questions in fundamental science. This is especially intriguing when colloids exhibit a degree of long-range positional or orientational ordering, as in liquid crystals (LCs) of elongated particles. Along with standard methods, microrheology (MR) has emerged in recent years as a tool to assess the mechanical properties of materials at the microscopic level. In particular, by active MR one can infer the viscoelastic response of a soft material from the dynamics of a tracer particle being dragged through it by external forces. Although considerable efforts have been made to study the diffusion of guest particles in LCs, little is known about the combined effect of tracer size and directionality of the dragging force on the system's viscoelastic response. By dynamic Monte Carlo simulations, we apply active MR to investigate the viscoelasticity of self-assembling smectic (Sm) LCs consisting of rodlike particles. In particular, we track the motion of a spherical tracer whose size is varied within a range of values matching the system's characteristic length scales and being dragged by constant forces that are parallel, perpendicular, or at 45° to the nematic director. Our results reveal a uniform value of the effective friction coefficient as probed by the tracer at small and large forces, whereas a nonlinear, force-thinning regime is observed at intermediate forces. However, at relatively weak forces the effective friction is strongly determined by correlations between the tracer size and the structure of the host fluid. Moreover, we also show that external forces forming an angle with the nematic director provide additional details that cannot be simply inferred from the mere analysis of parallel and perpendicular forces. Our results highlight the fundamental interplay between tracer size and force direction in assessing the MR of Sm LC fluids.

Unconventionally fast transport through sliding dynamics of rodlike particles in macromolecular networks

Transport of rodlike particles in confinement environments of macromolecular networks plays crucial roles in many important biological processes and technological applications. The relevant understanding has been limited to thin rods with diameter much smaller than network mesh size, although the opposite case, of which the dynamical behaviors and underlying physical mechanisms remain unclear, is ubiquitous. Here, we solve this issue by combining experiments, simulations and theory. We find a nonmonotonic dependence of translational diffusion on rod length, characterized by length commensuration-governed unconventionally fast dynamics which is in striking contrast to the monotonic dependence for thin rods. Our results clarify that such a fast diffusion of thick rods with length of integral multiple of mesh size follows sliding dynamics and demonstrate it to be anomalous yet Brownian. Moreover, good agreement between theoretical analysis and simulations corroborates that the sliding dynamics is an intermediate regime between hopping and Brownian dynamics, and provides a mechanistic interpretation based on the rod-length dependent entropic free energy barrier. The findings yield a principle, that is, length commensuration, for optimal design of rodlike particles with highly efficient transport in confined environments of macromolecular networks, and might enrich the physics of the diffusion dynamics in heterogeneous media.

Elastic response of colloidal smectic liquid crystals: Insights from microscopic theory

Elongated colloidal rods at sufficient packing conditions are known to form stable lamellar or smectic phases. Using a simplified volume-exclusion model, we propose a generic equation of state for hard-rod smectics that is robust against simulation results and is independent of the rod aspect ratio. We then extend our theory by exploring the elastic properties of a hard-rod smectic, including the layer compressibility (B) and bending modulus (K_{1}). By introducing weak backbone flexibility we are able to compare our predictions with experimental results on smectics of filamentous virus rods (fd) and find quantitative agreement between the smectic layer spacing, the out-of-plane fluctuation strength, as well as the smectic penetration length λ=sqrt[K_{1}/B]. We demonstrate that the layer bending modulus is dominated by director splay and depends sensitively on lamellar out-of-plane fluctuations that we account for on the single-rod level. We find that the ratio between the smectic penetration length and the lamellar spacing is about two orders of magnitude smaller than typical values reported for thermotropic smectics. We attribute this to the fact that colloidal smectics are considerably softer in terms of layer compression than their thermotropic counterparts while the cost of layer bending is of comparable magnitude.

Topological fine structure of smectic grain boundaries and tetratic disclination lines within three-dimensional smectic liquid crystals

Observing and characterizing the complex ordering phenomena of liquid crystals subjected to external constraints constitutes an ongoing challenge for chemists and physicists alike. To elucidate the delicate balance appearing when the intrinsic positional order of smectic liquid crystals comes into play, we perform Monte-Carlo simulations of rod-like particles in a range of cavities with a cylindrical symmetry. Based on recent insights into the topology of smectic orientational grain boundaries in two dimensions, we analyze the emerging three-dimensional defect structures from the perspective of tetratic symmetry. Using an appropriate three-dimensional tetratic order parameter constructed from the Steinhardt order parameters, we show that those grain boundaries can be interpreted as a pair of tetratic disclination lines that are located on the edges of the nematic domain boundary. Thereby, we shed light on the fine structure of grain boundaries in three-dimensional confined smectics.

Topology of Orientational Defects in Confined Smectic Liquid Crystals

We propose a general formalism to characterize orientational frustration of smectic liquid crystals in confinement by interpreting the emerging networks of grain boundaries as objects with a topological charge. In a formal idealization, this charge is distributed in pointlike units of quarter-integer magnitude, which we identify with tetratic disclinations located at the end points and nodes. This coexisting nematic and tetratic order is analyzed with the help of extensive Monte Carlo simulations for a broad range of two-dimensional confining geometries as well as colloidal experiments, showing how the observed defect networks can be universally reconstructed from simple building blocks. We further find that the curvature of the confining wall determines the anchoring behavior of grain boundaries, such that the number of nodes in the emerging networks and the location of their end points can be tuned by changing the number and smoothness of corners, respectively.

Microrheology of colloidal suspensions via dynamic Monte Carlo simulations

Understanding the rheology of colloidal suspensions is crucial in the formulation of a wide selection of industry-relevant products, such as paints, foods and inks. To characterise the viscoelastic behaviour of these soft materials, one can analyse the microscopic dynamics of colloidal tracers diffusing through the host fluid and generating local deformations and stresses. This technique, referred to as microrheology, links the bulk rheology of fluids to the microscopic dynamics at the particle scale. If tracers are subjected to external forces, rather than freely diffusing, it is called active microrheology. Motivated by the impact of microrheology in providing information on local structure in complex systems such as colloidal glasses, active matter or biological systems, we have extended the dynamic Monte Carlo (DMC) technique to investigate active microrheology in colloidal suspensions. The original DMC theoretical framework, able to accurately describe the Brownian dynamics of colloids at equilibrium, is here reconsidered and expanded to describe the effects of an external force pulling a tracer embedded in isotropic colloidal suspensions at different densities. To this end, we studied the dynamics of a spherical tracer dragged by a constant external force through a bath of spherical and rod-like particles of comparable size. We could extract valuable details on its effective friction coefficient, being constant at small and large values of the external force, but otherwise displaying a nonlinear behaviour that indicates the occurrence of a force-thinning regime. Our DMC simulation results are in excellent quantitative agreement with past Langevin dynamics simulations and theoretical works for the bath of spherical colloids. The bath of rod-like particles is studied in the isotropic phase, and displays an example where DMC is more convenient than Brownian or Langevin dynamics, in this case, in dealing with particle rotation.

When bigger is faster: A self-Van Hove analysis of the enhanced self-diffusion of non-commensurate guest particles in smectics

We investigate the anomalous dynamics in smectic phases of short host rods where, counter-intuitively, long guest rod-shaped particles diffuse faster than the short host ones due to their precise size mismatch. In addition to the previously reported mean-square displacement, we analyze the time evolution of the self-Van Hove functions G(r, t), as this probability density function uncovers intrinsic heterogeneous dynamics. Through this analysis, we show that the dynamics of the host particles parallel to the director becomes non-Gaussian and therefore heterogeneous after the nematic-to-smectic-A phase transition, even though it exhibits a nearly diffusive behavior according to its mean-squared displacement. In contrast, the non-commensurate guest particles display Gaussian dynamics of the parallel motion, up to the transition to the smectic-B phase. Thus, we show that the self-Van Hove function is a very sensitive probe to account for the instantaneous and heterogeneous dynamics of our system and should be more widely considered as a quantitative and complementary approach of the classical mean-squared displacement characterization in diffusion processes.

Particle-resolved topological defects of smectic colloidal liquid crystals in extreme confinement

Confined samples of liquid crystals are characterized by a variety of topological defects and can be exposed to external constraints such as extreme confinements with nontrivial topology. Here we explore the intrinsic structure of smectic colloidal layers dictated by the interplay between entropy and an imposed external topology. Considering an annular confinement as a basic example, a plethora of competing states is found with nontrivial defect structures ranging from laminar states to multiple smectic domains and arrays of edge dislocations, which we refer to as Shubnikov states in formal analogy to the characteristic of type-II superconductors. Our particle-resolved results, gained by a combination of real-space microscopy of thermal colloidal rods and fundamental-measure-based density functional theory of hard anisotropic bodies, agree on a quantitative level.

Helicoidal dynamics of biaxial curved rods in twist-bend nematic phases unveiled by unsupervised machine learning techniques

Uniaxial rods in a nematic phase diffuse preferentially in the direction parallel to the nematic director n[over ̂]. The nematic director field n[over ̂](r) of a chiral twist-bend nematic (N_{TB}) phase of achiral banana-shaped particles, recently discovered experimentally, displays a heliconical twist of given handedness and periodicity. Using simulations, we investigate the long-time macroscopic diffusion in N_{TB} phases, and find that the predilection of curved rods to diffuse in the direction of the twisting n[over ̂](r) yields a fascinating chiral dynamics along helices, even though achiral curved rods display Brownian motion with a nontrivial rototranslational coupling. We devise a machine learning protocol to characterize the helicoidal particle trajectories, finding that their pitch and radius are determined by the pitch and conical angle of the N_{TB} phase thereby connecting its structural and dynamical properties.

Dynamic Monte Carlo simulations of inhomogeneous colloidal suspensions

The dynamic Monte Carlo (DMC) method is an established molecular simulation technique for the analysis of the dynamics in colloidal suspensions. An excellent alternative to Brownian dynamics or molecular dynamics simulation, DMC is applicable to systems of spherical and/or anisotropic particles and to equilibrium or out-of-equilibrium processes. In this work, we present a theoretical and methodological framework to extend DMC to the study of heterogeneous systems, where the presence of an interface between coexisting phases introduces an additional element of complexity in determining the dynamic properties. In particular, we simulate a Lennard-Jones fluid at the liquid-vapor equilibrium and determine the diffusion coefficients in the bulk of each phase and across the interface. To test the validity of our DMC results, we also perform Brownian Dynamics simulations and unveil an excellent quantitative agreement between the two simulation techniques.