Physical foundation of the fluid particle dynamics method for colloid dynamics simulation

Heterogeneous solvent dissipation coupled with particle rearrangement in shear-thinning non-Brownian suspensions

Dense non-Brownian suspensions exhibit significant shear thinning, although a comprehensive understanding of the full scope of this phenomenon remains elusive. This study numerically reveals intimate heterogenous coupled dynamics between many-body particle motions and solvent hydrodynamics in shear-thinning non-Brownian suspensions. In our simulation systems, we do not account for frictional contact forces, reflecting experimental conditions under low shear rates where shear thinning occurs, while hydrodynamic interactions are directly incorporated using the Smoothed Profile Method. We demonstrate the spatially correlated viscous dissipation and particle motions; they share the same characteristic length, which decreases with increasing shear rate. We further show that, at lower shear rates, significant particle density changes are induced against the incompressibility of the solvent, suggesting the cooperative creation and annihilation of gaps and flow channels. We discuss that hydrodynamic interactions may substantially restrict particle rearrangements even in highly dense suspensions, influencing the quantitative aspects of macroscopic rheology.

Microrheology of active suspensions

We study the microrheology of active suspensions through direct hydrodynamic simulations using model pusher-like microswimmers. We demonstrate that the friction coefficient of a probe particle is notably reduced by hydrodynamic interactions (HIs) among a moving probe and the swimmers. When a swimmer approaches a probe from the rear (front) side, the repulsive HIs between them are weakened (intensified), which results in a slight front-rear asymmetry in swimmer orientation distribution around the probe, creating a significant additional net driving force acting on the probe from the rear side. The present drag-reduction mechanism qualitatively differs from that of the viscosity-reduction observed in sheared bulk systems and depends on probing details. This study provides insights into our fundamental knowledge of hydrodynamic effects in active suspensions and serves as a practical example illuminating distinctions between micro- and macrorheology measurements.

Impact of Hydrodynamic Interactions on the Kinetic Pathway of Protein Folding

Protein folding is a fundamental process critical to cellular function and human health, but it remains a grand challenge in biophysics. Hydrodynamic interaction (HI) plays a vital role in the self-organization of soft and biological materials, yet its role in protein folding is not fully understood despite folding occurring in a fluid environment. Here, we use the fluid particle dynamics method to investigate many-body hydrodynamic couplings between amino acid residues and fluid motion in the folding kinetics of a coarse-grained four-α-helices bundle protein. Our results reveal that HI helps select fast folding pathways to the native state without being kinetically trapped, significantly speeding up the folding kinetics compared to its absence. First, the directional flow along the protein backbone expedites protein collapse. Then, the incompressibility-induced squeezing flow effects retard the accumulation of non-native hydrophobic contacts, thus preventing the protein from being trapped in local energy minima during the conformational search of the native structure. We also find that the significance of HI in folding kinetics depends on temperature, with a pronounced effect under biologically relevant conditions. Our findings suggest that HI, particularly the short-range squeezing effect, may be crucial in avoiding protein misfolding.

Hydrodynamic Effects on the Collapse Kinetics of Flexible Polyelectrolytes

Understanding the collapse kinetics of polyelectrolytes (PEs) is crucial for comprehending various biological and industrial phenomena. Despite occurring in an aqueous environment, previous computational studies have overlooked the influence of hydrodynamic interactions (HIs) facilitated by fluid motion. Here, we directly compute the Navier-Stokes equation to investigate the collapse kinetics of a highly charged flexible PE. Our findings reveal that HI accelerates PE collapse induced by hydrophobicity and multivalent salt. In the case of hydrophobicity, HI induces long-range collective motion of monomers, accelerating the coarsening of local clusters through either Brownian-coagulation-like or evaporation-condensation-like processes, depending on the strength of hydrophobicity with respect to electrostatic interaction. Regarding multivalent salt, HI does not affect the condensation dynamics of multivalent ions but facilitates quicker movement of local dipolar clusters along the PE, thereby expediting the collapse process. These results provide valuable insights into the underlying mechanisms of HI in PE collapse kinetics.

Mechanical Slowing Down of Network-Forming Phase Separation of Polymer Solutions

Phase separation is a fundamental phenomenon leading to spatially heterogeneous material distribution, which is critical in nature, biology, material science, and industry. In ordinary phase separation, the minority phase always forms droplets. Contrary to this common belief, even the minority phase can form a network structure in viscoelastic phase separation (VPS). VPS can occur in any mixture with significant mobility differences between their components and is highly relevant to soft matter and biomatter. In contrast to classical phase separation, experiments have shown that VPS in polymer solutions lacks self-similar coarsening, resulting in the absence of a domain-coarsening scaling law. However, the underlying microscopic mechanism of this behavior remains unknown. To this end, we perform fluid particle dynamics simulations of bead-spring polymers, incorporating many-body hydrodynamic interactions between polymers through a solvent. We discover that polymers in the dense-network-forming phase are stretched and store elastic energy when the deformation speed exceeds the polymer dynamics. This self-generated viscoelastic stress mechanically interferes with phase separation and slows its dynamics, disrupting self-similar growth. We also highlight the essential role of many-body hydrodynamic interactions in VPS. The implications of our findings may hold importance in areas such as biological phase separation, porous material formation, and other fields where network structures play a pivotal role.

Impact of Inverse Squeezing Flow on the Self-Assembly of Oppositely Charged Colloidal Particles under Electric Field

Hydrodynamic interactions (HIs) play a critical role in the self-organization of colloidal suspensions and biological solutions. However, their roles have remained elusive particularly for charged soft matter systems. Here we consider the role of HIs in the self-assembly of oppositely charged colloidal particles, which is a promising candidate for electrical tunable soft materials. We employ the fluid particle dynamics method to consider many-body HIs and the coupling between the colloid, ion, and fluid motions. We find that, under a constant electric field, oppositely charged colloidal particles form clusters and percolate into a gel network, unlike bundlelike aggregates aligned in the field direction observed by Brownian dynamics simulations neglecting HIs. We reveal that the cluster-forming tendency originates from the incompressibility-induced "inverse squeezing flow" effect that dramatically slows down the disaggregation of attached colloids. Our findings indicate that the HI selects a unique kinetic pathway to the nonequilibrium colloidal self-assembly.

A dissipative particle dynamics model for studying dynamic phenomena in colloidal rod suspensions

A dissipative particle dynamics (DPD) model is developed and demonstrated for studying dynamics in colloidal rod suspensions. The solvent is modeled as conventional DPD particles, while individual rods are represented by a rigid linear chain consisting of overlapping solid spheres, which interact with solvent particles through a hard repulsive potential. The boundary condition on the rod surface is controlled using a surface friction between the solid spheres and the solvent particles. In this work, this model is employed to study the diffusion of a single colloid in the DPD solvent and compared with theoretical predictions. Both the translational and rotational diffusion coefficients obtained at a proper surface friction show good agreement with calculations based on the rod size defined by the hard repulsive potential. In addition, the system-size dependence of the diffusion coefficients shows that the Navier-Stokes hydrodynamic interactions are correctly included in this DPD model. Comparing our results with experimental measurements of the diffusion coefficients of gold nanorods, we discuss the ability of the model to correctly describe dynamics in real nanorod suspensions. Our results provide a clear reference point from which the model could be extended to enable the study of colloid dynamics in more complex situations or for other types of particles.

Power-law coarsening in network-forming phase separation governed by mechanical relaxation

A space-spanning network structure is a basic morphology in phase separation of soft and biomatter, alongside a droplet one. Despite its fundamental and industrial importance, the physical principle underlying such network-forming phase separation remains elusive. Here, we study the network coarsening during gas-liquid-type phase separation of colloidal suspensions and pure fluids, by hydrodynamic and molecular dynamics simulations, respectively. For both, the detailed analyses of the pore sizes and strain field reveal the self-similar network coarsening and the unconventional power-law growth more than a decade according to ℓ ∝ t1/2, where ℓ is the characteristic pore size and t is the elapsed time. We find that phase-separation dynamics is controlled by mechanical relaxation of the network-forming dense phase, whose limiting process is permeation flow of the solvent for colloidal suspensions and heat transport for pure fluids. This universal coarsening law would contribute to the fundamental physical understanding of network-forming phase separation.

Transient coarsening and the motility of optically heated Janus colloids in a binary liquid mixture

A gold-capped Janus particle suspended in a near-critical binary liquid mixture can self-propel under illumination. We have immobilized such a particle in a narrow channel and carried out a combined experimental and theoretical study of the non-equilibrium dynamics of a binary solvent around it - lasting from the very moment of switching illumination on until the steady state is reached. In the theoretical study we use both a purely diffusive and a hydrodynamic model, which we solve numerically. Our results demonstrate a remarkable complexity of the time evolution of the concentration field around the colloid. This evolution is governed by the combined effects of the temperature gradient and the wettability, and crucially depends on whether the colloid is free to move or is trapped. For the trapped colloid, all approaches indicate that the early time dynamics is purely diffusive and characterized by composition layers travelling with constant speed from the surface of the colloid into the bulk of the solvent. Subsequently, hydrodynamic effects set in. Anomalously large nonequilibrium fluctuations, which result from the temperature gradient and the vicinity of the critical point of the binary liquid mixture, give rise to strong concentration fluctuations in the solvent and to permanently changing coarsening patterns not observed for a mobile particle. The early time dynamics around initially still Janus colloids produces a force which is able to set the Janus colloid into motion. The propulsion due to this transient dynamics is in the direction opposite to that observed after the steady state is attained.

Influence of Hydrodynamic Interactions on Colloidal Crystallization

One of the biggest unresolved problems in crystallization phenomena is the significant discrepancy in the nucleation rate between experiments and simulations even for the simplest liquid, i.e., the hard-sphere system. A popular explanation for this discrepancy is the neglect of hydrodynamic interactions (HI) in simulation studies. By comparing simulations with and without HI, we show that the long-time diffusive dynamics of the colloids is slowed down more rapidly by hydrodynamic lubrication effects with increasing volume fraction. We find that the kinetics of both nucleation and growth are controlled by this long-time diffusion and that it is possible to account for most of the effects of HI by rescaling with this timescale. Therefore, we conclude that HI is not the primary cause of the accelerated nucleation rates observed in experiments.

Complex dynamical interplay between solid particles and flow in driven granular suspensions

Granular materials immersed in a fluid are ubiquitous in daily life, industry, and nature. They include food processing, pastes, cosmetics, paints, concretes, cements, muds, wet sands, snows, landslides, and lava flow. They are known to exhibit a rich variety of complex rheological behavior, but the role of a fluid component in such behavior has remained poorly understood due to the nonlocal and many-body nature of hydrodynamic interactions between solid particles mediated by the fluid. We address this fundamental problem by comparing the microrheological response of athermal granular suspensions with and without hydrodynamic interactions to an externally driven probe particle by numerical simulations. We find that the presence of the fluid drastically increases the drag coefficient of the probe particle by more than one order of magnitude near the jamming transition. We reveal that this is a consequence of the nontrivial long-range nature of hydrodynamic interactions, which originates from unlimited cumulative transmission of near-field hydrodynamic interactions due to the incompressibility of both fluid and solid particles. Force chain formation of solid particles is dynamically coupled with hydrodynamic flow, leading to strong spatiotemporal fluctuations of flow pattern and nonlinear rheological response. Our study reveals essential roles of hydrodynamic interactions in complex rheological behavior of dense granular suspensions under an external drive.

Bridging the gap between molecular dynamics and hydrodynamics in nanoscale Brownian motions

Through molecular dynamics simulations, we examined the hydrodynamic behavior of the Brownian motion of fullerene particles based on molecular interactions. The solvation free energy and velocity autocorrelation function (VACF) were calculated by using the Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) potentials for the solute-solvent and solvent-solvent interactions and by changing the size of the fullerene particles. We also measured the diffusion constant of the fullerene particles and the shear viscosity of the host fluid, and then the hydrodynamic radius aHD was quantified from the Stokes-Einstein relation. The aHD value exceeds that of the gyration radius of the fullerene when the solvation free energy exhibits largely negative values using the LJ potential. In contrast, aHD is similar to the size of bare fullerene when the solvation free energy is positive using the WCA potential. Furthermore, the VACF of the fullerene particles is directly comparable with the analytical expressions utilizing the Navier-Stokes equations both in incompressible and compressible forms. A hydrodynamic long-time tail t-3/2 is demonstrated for timescales longer than the kinematic time of the momentum diffusion over the particle size. However, the VACF at shorter timescales deviates from the hydrodynamic description, particularly for smaller fullerene particles and for the LJ potential. This occurs even though the compressible effect is considered when characterizing the decay of the VACF around the sound-propagation timescale over the particle size. These results indicate that the nanoscale Brownian motion is influenced by the solvation structure around the solute particles originating from the molecular interaction.

Hydrodynamic simulations of charge-regulation effects in colloidal suspensions

Self-organization of charged soft matter is of crucial importance in biology. However, it is an extremely complex phenomenon due to dynamical couplings between the hydrodynamic flow, ions, and charges of soft matter. For colloidal suspensions, the coupling between the former two has already been studied by numerical simulations, with the colloid surface charge being fixed. However, the self-organization of colloids and/or the application of an external electric field make the electrostatic environment of each colloid inhomogeneous in both space and time. Thus, this leads to inhomogenisation of the surface charge of each colloid under ionisation equilibrium conditions. This effect is known as "charge regulation" and is of great importance in various electrostatic and electrokinetic phenomena not only in colloid suspensions but also in solutions of biomolecules. However, there has so far been no success in taking the charge regulation effect into account in numerical simulations of colloidal electrokinetics. Here, we extend the fluid particle dynamics (FPD) method to incorporate the charge regulation effect. We present a theoretical formulation of the method and its application to two types of problems, where charge regulation plays an important role: (i) cluster formation of colloid particles and (ii) a single colloid particle under an external field. By these examples, we show not only the importance of considering charge-regulation effects in the self-organization of charged systems but also the applicability of our simulation method to more complex problems of charged soft matter systems.