Structure and Shear Response of Janus Colloid-Polymer Mixtures in Solution

Structure and dynamics of amphiphilic patchy cubes in a nanoslit under shear

Patchy nanocubes are intriguing materials with simple shapes and space-filling and multidirectional bonding properties. Previous studies have revealed various mesoscopic structures such as colloidal crystals in the solid regime and rod-like or fractal-like aggregates in the liquid regime of the phase diagram. Recent studies have also shown that mesoscopic structural properties, such as an average cluster size M and orientational order, in amphiphilic nanocube suspensions are associated with macroscopic viscosity changes, mainly owing to differences in cluster shape among patch arrangements. Although many studies have been conducted on the self-assembled structures of nanocubes in bulk, little is known about their self-assembly in nanoscale spaces or structural changes under shear. In this study, we investigated mixtures of one- and two-patch amphiphilic nanocubes confined in two flat parallel plates at rest and under shear using molecular dynamics simulations coupled with multiparticle collision dynamics. We considered two different patch arrangements for the two-patch particles and two different slit widths H to determine the degree of confinement in constant volume fractions in the liquid regime of the phase diagram. We revealed two unique cluster morphologies that have not been previously observed under bulk conditions. At rest, the size of the rod-like aggregates increased with decreasing H, whereas that of the fractal-like aggregates remained constant. Under weak shear with strong confinement, the rod-like aggregates maintained a larger M than the fractal-like aggregates, which were more rigid and maintained a larger M than the rod-like aggregates under bulk conditions.

Mesoscale simulations of diffusion and sedimentation in shape-anisotropic nanoparticle suspensions

We determine the long-time self-diffusion coefficient and sedimentation coefficient for suspensions of nanoparticles with anisotropic shapes (octahedra, cubes, tetrahedra, and spherocylinders) as a function of nanoparticle concentration using mesoscale simulations. We use a discrete particle model for the nanoparticles, and we account for solvent-mediated hydrodynamic interactions between nanoparticles using the multiparticle collision dynamics method. Our simulations are compared to theoretical predictions and experimental data from existing literature, demonstrating good agreement in the majority of cases. Further, we find that the self-diffusion coefficient of the regular polyhedral shapes can be estimated from that of a sphere whose diameter is the average of their inscribed and circumscribed sphere diameters.

Aggregation of amphiphilic nanocubes in equilibrium and under shear

We investigate the self-assembly of amphiphilic nanocubes into finite-sized aggregates in equilibrium and under shear, using molecular dynamics (MD) simulations and kinetic Monte Carlo (KMC) calculations. These patchy nanoparticles combine both interaction and shape anisotropy, making them valuable models for studying folded proteins and DNA-functionalized nanoparticles. The nanocubes can self-assemble into various finite-sized aggregates ranging from rods to self-avoiding random walks, depending on the number and placement of the hydrophobic faces. Our study focuses on suspensions containing multi- and one-patch cubes, with their ratio systematically varied. When the binding energy is comparable to the thermal energy, the aggregates consist of only few cubes that spontaneously associate/dissociate. However, highly stable aggregates emerge when the binding energy exceeds the thermal energy. Generally, the mean aggregation number of the self-assembled clusters increases with the number of hydrophobic faces and decreases with increasing fraction of one-patch cubes. In sheared suspensions, the more frequent collisions between nanocube clusters lead to faster aggregation dynamics but also to smaller terminal steady-state mean cluster sizes. The results from the MD and KMC simulations are in excellent agreement for all investigated two-patch cases, whereas the three-patch cubes form systematically smaller clusters in the MD simulations compared to the KMC calculations due to finite-size effects and slow aggregation kinetics. By analyzing the rate kernels, we are able to identify the primary mechanisms responsible for (shear-induced) cluster growth and breakup. This understanding allows us to tune nanoparticle and process parameters to achieve desired cluster sizes and shapes.

Self-Assembly of Amphiphilic Cubes in Suspension

We study the self-assembly of amphiphilic cubic colloids using molecular dynamics as well as rejection-free kinetic Monte Carlo simulations. We vary both the number and location of the solvophobic faces (patches) on the cubes at several colloid volume fractions and determine the resulting size and shape distributions of the self-assembled aggregates. When the binding energy is comparable to the thermal energy of the system, aggregates typically consist of only few spontaneously associating/dissociating colloids. Increasing the binding energy (or lowering the temperature) leads to the emergence of highly stable aggregates, e.g., small dimers in pure suspensions of one-patch cubes or large (system-spanning) aggregates in suspensions of multipatch colloids. In mixtures of one- and multipatch cubes, the average aggregation number increases with increasing number of solvophobic faces on the multipatch cubes as well with increasing fraction of multipatch cubes. The resulting aggregate shapes range from elongated rods over fractal objects to compact spheres, depending on the number and arrangement of solvophobic patches on the cubic colloids. Our findings establish the complex self-assembly pathways for a class of building blocks that combine both interaction and shape anisotropy, with the potential of forming hierarchically ordered superstructures.

Diffusion and sedimentation in colloidal suspensions using multiparticle collision dynamics with a discrete particle model

We study self-diffusion and sedimentation in colloidal suspensions of nearly hard spheres using the multiparticle collision dynamics simulation method for the solvent with a discrete mesh model for the colloidal particles (MD+MPCD). We cover colloid volume fractions from 0.01 to 0.40 and compare the MD+MPCD simulations to experimental data and Brownian dynamics simulations with free-draining hydrodynamics (BD) as well as pairwise far-field hydrodynamics described using the Rotne-Prager-Yamakawa mobility tensor (BD+RPY). The dynamics in MD+MPCD suggest that the colloidal particles are only partially coupled to the solvent at short times. However, the long-time self-diffusion coefficient in MD+MPCD is comparable to that in experiments, and the sedimentation coefficient in MD+MPCD is in good agreement with that in experiments and BD+RPY, suggesting that MD+MPCD gives a reasonable description of hydrodynamic interactions in colloidal suspensions. The discrete-particle MD+MPCD approach is convenient and readily extended to more complex shapes, and we determine the long-time self-diffusion coefficient in suspensions of nearly hard cubes to demonstrate its generality.