Dynamic Monte Carlo simulations of inhomogeneous colloidal suspensions

Structural and dynamical equilibrium properties of hard board-like particles in parallel confinement

We performed Monte Carlo and dynamic Monte Carlo simulations to model the diffusion of monodispersed suspensions composed of impenetrable cuboidal particles, specifically hard board-like particles (HBPs), in the presence of parallel hard walls. The impact of the walls was investigated by adjusting the size of the simulation box while maintaining constant packing fractions, fixed at η = 0.150, for systems consisting of HBPs with prolate, dual-shaped, and oblate geometries. We observed that increasing the distance between the walls led to the recovery of an isotropic bulk phase, while local particle organization near the walls remained stable. Due to their shape, oblate HBPs exhibit more efficient anchoring at wall surfaces compared to prolate shapes. The formation of nematic-like particle assemblies near the walls, confirmed by theoretical calculations based on density functional theory, significantly influenced local particle dynamics. This effect was particularly pronounced to the extent that a modest portion of cuboids near the walls tended to diffuse exclusively in planes parallel to the confinement, even more efficiently than observed in the bulk regions.

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.

Brownian dynamics simulations of hard rods in external fields and with contact interactions

We propose a simulation method for Brownian dynamics of hard rods in one dimension for arbitrary continuous external force fields. It is an event-driven procedure based on the fragmentation and mergers of clusters formed by particles in contact. It allows one to treat particle interactions in addition to the hard-sphere exclusion as long as the corresponding interaction forces are continuous functions of the particle coordinates. We furthermore develop a treatment of sticky hard spheres as described by Baxter's contact interaction potential.

Active microrheology of colloidal suspensions of hard cuboids

By performing dynamic Monte Carlo simulations, we investigate the microrheology of isotropic suspensions of hard-core colloidal cuboids. In particular, we infer the local viscoelastic behavior of these fluids by studying the dynamics of a probe spherical particle that is incorporated in the host phase and is dragged by an external force. This technique, known as active microrheology, allows one to characterize the microscopic response of soft materials upon application of a constant force, whose intensity spans here three orders of magnitude. By tuning the geometry of cuboids from oblate to prolate as well as the system density, we observe different responses that are quantified by measuring the effective friction perceived by the probe particle. The resulting friction coefficient exhibits a linear regime at forces that are much weaker and larger than the thermal forces, whereas a nonlinear, force-thinning regime is observed at intermediate force intensities.

Long-time relaxation dynamics in nematic and smectic liquid crystals of soft repulsive colloidal rods

Understanding the relaxation dynamics of colloidal suspensions is crucial to identifying the elements that influence the mobility of their constituents, assessing their macroscopic response across the relevant time and length scales, and thus disclosing the fundamentals underpinning their exploitation in formulation engineering. In this work, we specifically assess the impact of long-ranged ordering on the relaxation dynamics of suspensions of soft repulsive rodlike particles, which are able to self-organize into nematic and smectic liquid-crystalline phases. Rods are modeled as soft repulsive spherocylinders with a length-to-diameter ratio L^{★}=5, interacting via the truncated and shifted Kihara potential. By performing dynamic Monte Carlo simulations, we analyze the effect of translational and orientational order on the diffusion of the rods along the relevant directions imposed by the morphology of the background phases. To provide a clear picture of the resulting dynamics, we assess its dependence on temperature, which can dramatically determine the response time of the system relaxation and the self-diffusion coefficients of the rods. The computation of the van Hove correlation functions allows us to identify the existence of rods that diffuse significantly faster than the average and whose concentration can be accurately adjusted by a suitable choice of temperature.

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.

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.