Search

We theoretically study the conformational and dynamical properties of semiflexible active polar ring polymers under linear shear flow. A ring is described as a continuous semiflexible Gaussian polymer with a tangential active force of a constant density along its contour. The linear but non-Hermitian equation of motion is solved using an eigenfunction expansion, which yields activity-independent, but shear-rate-dependent, relaxation times and activity-dependent frequencies. As a consequence, the ring's stationary-state properties are independent of activity, and its conformations and rheological properties are equal to those of a passive ring under shear. The presence of characteristic time scales by relaxation and the activity-dependent frequencies give rise to a particular dynamical behavior. A tank-treading-like motion emerges for long relaxation times and high activities, specifically for stiff rings. In the case of very flexible polymers, the relaxation behavior dominates over activity contributions suppressing tank-treading. Shear strongly affects the crossover from a tank-treading to a relaxation-dominated dynamics, and the ring polymer exhibits tumbling motion at high shear rates. This is reflected in the tumbling frequency, which displays two shear-rate dependent regimes, with an activity-dependent plateau at low shear rates followed by a power-law regime with increasing tumbling frequency for high shear rates., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

We present Brownian dynamics simulation results of a flexible linear polymer with excluded-volume interactions under shear flow in the presence of active noise. The active noise strongly affects the polymer's conformational and dynamical properties, such as the stretching in the flow direction and compression in the gradient direction, shear-induced alignment, and shear viscosity. In the asymptotic limit of large activities and shear rates, the power-law scaling exponents of these quantities differ significantly from those of passive polymers. The chain's shear-induced stretching at a given shear rate is reduced by active noise, and it displays a non-monotonic behavior, where an initial polymer compression is followed by its stretching with increasing active force. The compression of the polymer in the gradient direction follows the relation ∼Wi Pe -3/4 as a function of the activity-dependent Weissenberg number Wi Pe , which differs from the scaling observed in passive systems ∼Wi Pe -1/2 . The flow-induced alignment at large Péclet numbers Pe ≫ 1, where Pe is the Péclet number, and large shear rates Wi Pe ≫ 1 displays the scaling behavior Wi Pe -1/2 , with an exponent differing from the passive value -1/3. Furthermore, the polymer's zero-shear viscosity displays a non-monotonic behavior, decreasing in an intermediate activity regime due to excluded-volume interactions and increasing again for large Pe. Shear thinning appears with increasing Weissenberg number with the power-laws Wi Pe -1/2 and Wi Pe -3/4 for passive and active polymers, respectively. In addition, our simulation results are compared with the results of an analytical approach, which predicts quantitatively similar behaviors for the various aforementioned physical quantities.

The conformational and dynamical properties of isolated flexible active polar linear polymers (APLPs) are studied analytically. The APLPs are modeled as Gaussian bead-spring linear chains augmented by tangential active forces, both in a discrete and continuous representation. The polar forces lead to linear non-Hermitian equations of motion, which are solved by an eigenfunction expansion in terms of a biorthogonal basis set. Our calculations show that the polymer conformations are independent of activity. However, tangential propulsion strongly impacts the polymer dynamics and yields an active ballistic regime as well as an activity-enhanced long-time diffusive regime, which are both absent in passive systems. The polar forces imply a coupling of modes in the eigenfunction representation, in particular with the translational mode, with a respective strong influence on the polymer dynamics. The total polymer mean-square displacement on scales smaller than the radius of gyration is determined by the active internal dynamics rather than the collective center-of-mass motion, in contrast to active Brownian polymers at large Péclet numbers, reflecting the distinct difference in the propulsion mechanism.

Systems comprised of self-steering active Brownian particles are studied via simulations for a minimal cognitive flocking model. The dynamics of the active Brownian particles is extended by an orientational response with limited maneuverability to an instantaneous visual input of the positions of neighbors within a vision cone and a cut-off radius. The system exhibits large-scale self-organized structures, which depend on selected parameter values, and, in particular, the presence of excluded-volume interactions. The emergent structures in two dimensions, such as worms, worm-aggregate coexistence, and hexagonally close-packed structures, are analysed and phase diagrams are constructed. The analysis of the particle's mean-square displacement shows ABP-like dynamics for dilute systems and the worm phase. In the limit of densely packed structures, the active diffusion coefficient is significantly smaller and depends on the number of particles in the cluster. Our analysis of the cluster-growth dynamics shows distinct differences to processes in systems of short-range attractive colloids in equilibrium. Specifically, the characteristic time for the growth and decay of clusters of a particular size is longer than that of isotropically attractive colloids, which we attribute to the non-reciprocal nature of the directed visual perception. Our simulations reveal a strong interplay between ABP-characteristic interactions, such as volume exclusion and rotational diffusion, and cognitive-based interactions and navigation.

The conformational and dynamical properties of isolated semiflexible active polar ring polymers are investigated analytically. A ring is modeled as a continuous Gaussian polymer exposed to tangential active forces. The analytical solution of the linear non-Hermitian equation of motion in terms of an eigenfunction expansion shows that ring conformations are independent of activity. In contrast, activity strongly affects the internal ring dynamics and yields characteristic time regimes, which are absent in passive rings. On intermediate timescales, flexible rings show an activity-enhanced diffusive regime, while semiflexible rings exhibit ballistic motion. Moreover, a second active time regime emerges on longer timescales, where rings display a snake-like motion, which is reminiscent to a tank-treading rotational dynamics in shear flow, dominated by the mode with the longest relaxation time.

The properties of microswimmer dumbbells composed of pusher-puller pairs are investigated by mesoscale hydrodynamic simulations employing the multiparticle collision dynamics approach for the fluid. An individual microswimmer is represented by a squirmer, and various active-stress combinations in a dumbbell are considered. The squirmers are connected by a bond, which does not impose any geometrical restriction on the individual rotational motion. Our simulations reveal a strong influence of the squirmers' flow fields on the orientation of their propulsion directions, their fluctuations, and the swimming behavior of a dumbbell. The properties of pusher-puller pairs with an equal magnitude of the active stresses depend only weakly on the stress magnitude. This is similar to dumbbells of microswimmers without hydrodynamic interactions. However, for non-equal stress magnitudes, the active stress implies strong orientational correlations of the swimmers' propulsion directions with respect to each other, as well as the bond vector. The orientational coupling is most pronounced for pairs with large differences in the active-stress magnitude. The alignment of the squirmers' propulsion directions with respect to each other is preferentially orthogonal in dumbbells with a strong pusher and weak puller, and antiparallel in the opposite case when the puller dominates. These strong correlations affect the active motion of dumbbells, which is faster for strong pushers and slower for strong pullers.

Semiflexible polymers comprised of active Brownian particles (ABPOs) exhibit intriguing activity-driven conformational and dynamical features. Analytically, the generic properties of ABPOs can be obtained in a mean-field description applying the Gaussian semiflexible polymer model. In this article, we derive a path integral representation of the stationary-state distribution function of such ABPOs, based on the stationary-state distribution function of the normal mode amplitudes following from the Langevin equation of motion. The path integral includes characteristic semiflexible polymer contributions from entropy and bending energy, with activity dependent coefficients, and, in addition, activity-induced torsional and higher order correlations along the polymer contour. Focusing on a semiflexible polymer approximation, we determine various properties such as the tangent-vector correlation function, effective persistence length, and the mean-square end-to-end distance. The latter reflects the characteristic features of ABPOs, and good quantitative agreement is obtained with the full solution for larger activities, specifically for flexible polymers. Moreover, the approximation indicates the relevance of torsional and higher order contour correlations for the ABPO conformations. In general, the ABPO path integral illustrates how colored noise (active fluctuations) affects semiflexible polymer conformations in comparison to white noise thermal fluctuations.

The conformational and dynamical properties of active Brownian polymers embedded in a fluid depend on the nature of the driving mechanism, e.g., self-propulsion or external actuation of the monomers. Implementations of self-propelled and actuated active Brownian polymers in a multiparticle collision (MPC) dynamics fluid are presented, which capture the distinct differences between the two driving mechanisms. The active force-free nature of self-propelled monomers requires adaptations of the MPC simulation scheme, with its streaming and collision steps, where the monomer self-propulsion velocity has to be omitted in the collision step. Comparison of MPC simulation results for active polymers in dilute solution with results of Brownian dynamics simulations accounting for hydrodynamics via the Rotne-Prager-Yamakawa tensor confirm the suitability of the implementation. The polymer conformational and dynamical properties are analyzed by the static and dynamic structure factor, and the scaling behavior of the latter with respect to the wave number and time dependence are discussed. The dynamic structure factor displays various activity-induced temporal regimes, depending on the considered wave number, which reflect the persistent diffusive motion of the whole polymer at small wave numbers, and the activity-enhanced internal dynamics at large wave numbers. The obtained simulation results are compared with theoretical predictions.

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

The conformational and dynamical properties of a passive polymer embedded in a bath of active Brownian particles (ABPs) are studied by Langevin dynamics simulations. Various activities and ABP concentrations below and above the critical values for motility-induced phase separation (MIPS) are considered. In a homogeneous ABP fluid, the embedded polymer swells with increasing bath activity, with stronger swelling for larger densities. The polymer dynamics is enhanced, with the diffusion coefficient increasing by a power-law with increasing activity, where the exponent depends on the ABP concentration. For ABP concentrations in the MIPS regime, we observe a localization of the polymer in the low-density ABP phase associated with polymer collapse for moderate activities and a reswelling for high activities accompanied by a preferred localization in the high-density ABP phase. Localization and reswelling are independent of the polymer stiffness, with stiff polymers behaving similarly to flexible polymers. The polymer collapse is associated with a slowdown of its dynamics and a significantly smaller center-of-mass diffusion coefficient. In general, the polymer dynamics can only partially be described by an effective (bath) temperature. Moreover, the properties of a polymer embedded in a homogeneous active bath deviate quantitatively from those of a polymer composed of active monomers, i.e., linear chains of ABPs; however, such a polymer exhibits qualitatively similar activity-dependent features.

The properties of semiflexible polymers tethered by one end to an impenetrable wall and exposed to oscillatory shear flow are investigated by mesoscale simulations. A polymer, confined in two dimensions, is described by a linear bead-spring chain, and fluid interactions are incorporated by the Brownian multiparticle collision dynamics approach. At small strain, the polymers follow the applied flow field. However, at high strain, we find a strongly nonlinear response with major conformational changes. Polymers are stretched along the flow direction and exhibit U-shaped conformations while following the flow. As a consequence of confinement in the half-space, frequency doubling in the time-dependent polymer properties appears along the direction normal to the wall.

The surface distribution of flagella in peritrichous bacterial cells has been traditionally assumed to be random. Recently, the presence of a regular grid-like pattern of basal bodies has been suggested. Experimentally, the manipulation of the anchoring points of flagella in the cell membrane is difficult, and thus, elucidation of the consequences of a particular pattern on bacterial locomotion is challenging. We analyze the bundle formation process and swimming properties of Bacillus subtilis-like cells considering random, helical, and ring-like arrangements of flagella by means of mesoscale hydrodynamics simulations. Helical and ring patterns preferentially yield configurations with a single bundle, whereas configurations with no clear bundles are most likely for random anchoring. For any type of pattern, there is an almost equally low probability to form V-shaped bundle configurations with at least two bundles. Variation of the flagellum length yields a clear preference for a single major bundle in helical and ring patterns as soon as the flagellum length exceeds the body length. The average swimming speed of cells with a single or two bundles is rather similar, and approximately [Formula: see text] larger than that of cells of other types of flagellar organization. Considering the various anchoring patterns, rings yield the smallest average swimming speed independent of the type of bundle, followed by helical arrangements, and largest speeds are observed for random anchoring. Hence, a regular pattern provides no advantage in terms of swimming speed compared to random anchoring of flagella, but yields more likely single-bundle configurations.

Correction for 'Hydrodynamic interactions in squirmer dumbbells: active stress-induced alignment and locomotion' by Judit Clopés et al., Soft Matter, 2020, 16, 10676-10687, DOI: .

Hydrodynamic interactions are fundamental for the dynamics of swimming self-propelled particles. Specifically, bonds between microswimmers enforce permanent spatial proximity and, thus, enhance emergent correlations by microswimmer-specific flow fields. We employ the squirmer model to study the swimming behavior of microswimmer dumbbells by mesoscale hydrodynamic simulations, where the squirmers' rotational motion is geometrically unrestricted. An important aspect of the applied particle-based simulation approach-the multiparticle collision dynamics method-is the intrinsic account for thermal fluctuations. We find a strong effect of active stress on the motility of dumbbells. In particular, pairs of strong pullers exhibit orders of magnitude smaller swimming efficiency than pairs of pushers. This is a consequence of the inherent thermal fluctuations in combination with the strong coupling of the squirmers' rotational motion, which implies non-exponentially decaying auto- and cross-correlation functions of the propulsion directions, and active stress-dependent characteristic decay times. As a consequence, specific stationary-state relative alignments of the squirmer propulsion directions emerge, where pullers are preferentially aligned in an antiparallel manner along the bond vector, whereas pushers are preferentially aligned normal to the bond vector with a relative angle of approximately 60° at weak active stress, and one of the propulsion directions is aligned with the bond at strong active stress. The distinct differences between dumbbells comprised of pusher or pullers suggest means to control microswimmer assemblies for future microbot applications.

Adding shape and interaction anisotropy to a colloidal particle offers exquisitely tunable routes to engineer a rich assortment of complex-architected structures. Inspired by the hierarchical self-assembly concept with block copolymers and DNA liquid crystals and exploiting the unique assembly properties of DNA, we report here the construction and self-assembly of DNA-based soft-patchy anisotropic particles with a high degree of modularity in the system's design. By programmable positioning of thermoresponsive polymeric patches on the backbone of a stiff DNA duplex with linear and star-shaped architecture, we reversibly drive the DNA from a disordered ensemble to a diverse array of long-range ordered multidimensional nanostructures with tunable lattice spacing, ranging from lamellar to bicontinuous double-gyroid and double-diamond cubic morphologies, through the alteration of temperature. Our results demonstrate that the proposed hierarchical self-assembly strategy can be applied to any kind of DNA nanoarchitecture, highlighting the design principles for integration of self-assembly concepts from the physics of liquid crystals, block copolymers, and patchy colloids into the continuously growing interdisciplinary research field of structural DNA nanotechnology.

Active matter agents consume internal energy or extract energy from the environment for locomotion and force generation. Already, rather generic models, such as ensembles of active Brownian particles, exhibit phenomena, which are absent at equilibrium, particularly motility-induced phase separation and collective motion. Further intriguing nonequilibrium effects emerge in assemblies of bound active agents as in linear polymers or filaments. The interplay of activity and conformational degrees of freedom gives rise to novel structural and dynamical features of individual polymers, as well as in interacting ensembles. Such out-of-equilibrium polymers are an integral part of living matter, ranging from biological cells with filaments propelled by motor proteins in the cytoskeleton and RNA/DNA in the transcription process to long swarming bacteria and worms such as Proteus mirabilis and Caenorhabditis elegans, respectively. Even artificial active polymers have been synthesized. The emergent properties of active polymers or filaments depend on the coupling of the active process to their conformational degrees of freedom, aspects that are addressed in this article. The theoretical models for tangentially and isotropically self-propelled or active-bath-driven polymers are presented, both in the presence and absence of hydrodynamic interactions. The consequences for their conformational and dynamical properties are examined, with emphasis on the strong influence of the coupling between activity and hydrodynamic interactions. Particular features of emerging phenomena in semi-dilute systems, induced by steric and hydrodynamic interactions, are highlighted. Various important, yet theoretically unexplored, aspects are featured, and future challenges are discussed.

Microswimmers such as E. coli bacteria accumulate and exhibit an intriguing dynamics near walls, governed by hydrodynamic and steric interactions. Insight into the underlying mechanisms and predominant interactions demand a detailed characterization of the entrapment process. We employ a mesoscale hydrodynamics simulation approach to study entrapment of an E. coli-type cell at a no-slip wall. The cell is modeled by a spherocylindrical body with several explicit helical flagella. Three stages of the entrapment process can be distinguished: the approaching regime, where a cell swims toward a wall on a nearly straight trajectory; a scattering regime, where the cell touches the wall and reorients; and a surface-swimming regime. Our simulations show that steric interactions may dominate the entrapment process, yet, hydrodynamic interactions slow down the adsorption dynamics close to the boundary and imply a circular motion on the wall. The locomotion of the cell is characterized by a strong wobbling dynamics, with cells preferentially pointing toward the wall during surface swimming.

Activity and autonomous motion are fundamental in living and engineering systems. This has stimulated the new field of 'active matter' in recent years, which focuses on the physical aspects of propulsion mechanisms, and on motility-induced emergent collective behavior of a larger number of identical agents. The scale of agents ranges from nanomotors and microswimmers, to cells, fish, birds, and people. Inspired by biological microswimmers, various designs of autonomous synthetic nano- and micromachines have been proposed. Such machines provide the basis for multifunctional, highly responsive, intelligent (artificial) active materials, which exhibit emergent behavior and the ability to perform tasks in response to external stimuli. A major challenge for understanding and designing active matter is their inherent nonequilibrium nature due to persistent energy consumption, which invalidates equilibrium concepts such as free energy, detailed balance, and time-reversal symmetry. Unraveling, predicting, and controlling the behavior of active matter is a truly interdisciplinary endeavor at the interface of biology, chemistry, ecology, engineering, mathematics, and physics. The vast complexity of phenomena and mechanisms involved in the self-organization and dynamics of motile active matter comprises a major challenge. Hence, to advance, and eventually reach a comprehensive understanding, this important research area requires a concerted, synergetic approach of the various disciplines. The 2020 motile active matter roadmap of Journal of Physics: Condensed Matter addresses the current state of the art of the field and provides guidance for both students as well as established scientists in their efforts to advance this fascinating area.

The conformational and dynamical properties of active polymers in solution are determined by the nature of the activity. Here, the behavior of polymers with self-propelled, active Brownian particle-type monomers differs qualitatively from that of polymers with monomers driven externally by colored-noise forces. We present simulation and theoretical results for polymers in solution in the presence of external active noise. In simulations, a semiflexible bead-spring chain is considered, in analytical calculations, a continuous linear wormlike chain. Activity is taken into account by independent monomer or site velocities, with orientations changing in a diffusive manner. In simulations, hydrodynamic interactions (HIs) are taken into account by the Rotne-Prager-Yamakawa tensor or by an implementation of the active polymer in the multiparticle-collision-dynamics approach for fluids. To arrive at an analytical solution, the preaveraged Oseen tensor is employed. The active process implies a dependence of the stationary-state properties on HIs via the polymer relaxation times. With increasing activity, HIs lead to an enhanced swelling of flexible polymers, and the conformational properties differ substantially from those of polymers with self-propelled monomers in the presence of HIs, or free-draining polymers. The polymer mean-square displacement is enhanced by HIs. Over a wide range of timescales, hydrodynamics leads to a subdiffusive regime of the site mean-square displacement for flexible active polymers, with an exponent of 5/7, larger than that of the Rouse (1/2) and Zimm (2/3) models of passive polymers.

Ensembles of actuated colloids are excellent model systems to explore emergent out-of-equilibrium structures, complex collective dynamics, and design rules for the next generation materials. Here, we demonstrate that ferromagnetic microparticles suspended at an air-water interface and energized by an external rotating magnetic field spontaneously form dynamic ensembles of synchronized spinners in a certain range of the excitation field parameters. Each spinner generates strong hydrodynamic flows, and collective interactions of the multiple spinners promote a formation of dynamic lattices. On the basis of experiments and simulations, we reveal structural transitions from liquid to nearly crystalline states in this novel active spinner material and demonstrate that dynamic spinner lattices are reconfigurable, capable of self-healing behavior and that the transport of embedded inert cargo particles can be remotely tuned by the parameters of the external excitation field. Our findings provide insights into the behavior of active spinner materials with reconfigurable structural order and tunable functionalities., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

The rotational diffusive motion of a self-propelled, attractive spherical colloid immersed in a solution of self-avoiding polymers is studied by mesoscale hydrodynamic simulations. A drastic enhancement of the rotational diffusion by more than an order of magnitude in the presence of activity is obtained. The amplification is a consequence of two effects, a decrease of the amount of adsorbed polymers by active motion and an asymmetric encounter with polymers on the squirmer surface, which yields an additional torque and random noise for the rotational motion. Our simulations suggest a way to control the rotational dynamics of squirmer-type microswimmers by the degree of polymer adsorption and system heterogeneity.

The emergent fluctuating hydrodynamics of a viscoelastic fluid modeled by the multiparticle collision dynamics (MPC) approach is studied. The fluid is composed of flexible, Gaussian phantom polymers that interact by local momentum-conserving stochastic MPCs. For comparison, the analytical solution of the linearized Navier-Stokes equation is calculated, where viscoelasticity is taken into account by a time-dependent shear relaxation modulus. The fluid properties are characterized by the transverse velocity autocorrelation function in Fourier space as well as in real space. Various polymer lengths are considered-from dumbbells to (near-)continuous polymers. Viscoelasticity affects the fluid properties and leads to strong correlations, which overall decay exponentially in Fourier space. In real space, the center-of-mass velocity autocorrelation function of individual polymers exhibits a long-time tail, independent of the polymer length, which decays as t -3/2 , similar to a Newtonian fluid, in the asymptotic limit t → ∞. Moreover, for long polymers, an additional power-law decay appears at time scales shorter than the longest polymer relaxation time with the same time dependence, but negative correlations, and the polymer length dependence L -1/2 . Good agreement is found between the analytical and simulation results.

The conformational and dynamical properties of active self-propelled filaments/polymers are investigated in the presence of hydrodynamic interactions by both, Brownian dynamics simulations and analytical theory. Numerically, a discrete linear chain composed of active Brownian particles is considered, analytically, a continuous linear semiflexible polymer with active velocities changing diffusively. The force-free nature of active monomers is accounted for-no Stokeslet fluid flow induced by active forces-and higher order hydrodynamic multipole moments are neglected. Hence, fluid-mediated interactions are assumed to arise solely due to intramolecular forces. The hydrodynamic interactions (HI) are taken into account analytically by the preaveraged Oseen tensor, and numerically by the Rotne-Prager-Yamakawa tensor. The nonequilibrium character of the active process implies a dependence of the stationary-state properties on HI via the polymer relaxation times. In particular, at moderate activities, HI lead to a substantial shrinkage of flexible and semiflexible polymers to an extent far beyond shrinkage of comparable free-draining polymers; even flexible HI-polymers shrink, while active free-draining polymers swell monotonically. Large activities imply a reswelling, however, to a less extent than for non-HI polymers, caused by the shorter polymer relaxation times due to hydrodynamic interactions. The polymer mean square displacement is enhanced, and an activity-determined ballistic regime appears. Over a wide range of time scales, flexible active polymers exhibit a hydrodynamically governed subdiffusive regime, with an exponent significantly smaller than that of the Rouse and Zimm models of passive polymers. Compared to simulations, the analytical approach predicts a weaker hydrodynamic effect. Overall, hydrodynamic interactions modify the conformational and dynamical properties of active polymers substantially.

The stress of a fluid on a confining wall is given by the mechanical wall forces, independent of the nature of the fluid being passive or active. At thermal equilibrium, an equation of state exists and stress is likewise obtained from intrinsic bulk properties; even more, stress can be calculated locally. Comparable local descriptions for active systems require a particular consideration of active forces. Here, we derive expressions for the stress exerted on a local volume of a systems of spherical active Brownian particles (ABPs). Using the virial theorem, we obtain two identical stress expressions, a stress due to momentum flux across a hypothetical plane, and a bulk stress inside of the local volume. In the first case, we obtain an active contribution to momentum transport in analogy to momentum transport in an underdamped passive system, and we introduce an active momentum. In the second case, a generally valid expression for the swim stress is derived. By simulations, we demonstrate that the local bulk stress is identical to the wall stress of a confined system for both, non-interacting ABPs as well as ABPs with excluded-volume interactions. This underlines the existence of an equation of state for a system of spherical ABPs. Most importantly, our calculations demonstrated that active stress is not a wall (boundary) effect, but is caused by momentum transport. We demonstrate that the derived stress expression permits the calculation of the local stress in inhomogeneous systems of ABPs.

The properties of a semiflexible polymer with fixed ends exposed to oscillatory shear flow are investigated by simulations. The two-dimensionally confined polymer is modeled as a linear bead-spring chain, and the interaction with the fluid is described by the Brownian multiparticle collision dynamics approach. For small shear rates, the tethering of the ends leads to a more-or-less linear oscillatory response. However, at high shear rates, we found a strongly nonlinear reaction, with a polymer (partially) wrapped around the fixation points. This leads to an overall shrinkage of the polymer. Dynamically, the location probability of the polymer center-of-mass position is largest on a spatial curve resembling a limaçon, although with an inhomogeneous distribution. We found shear-induced modifications of the normal-mode correlation functions, with a frequency doubling at high shear rates. Interestingly, an even-odd asymmetry for the Cartesian components of the correlation functions appears, with rather similar spectra for odd x - and even y -modes and vice versa. Overall, our simulations yielded an intriguing nonlinear behavior of tethered semiflexible polymers under oscillatory shear flow.

The conformational and dynamical properties of semiflexible active Brownian ring polymers are investigated analytically. A ring is described by the Gaussian semiflexible polymer model accounting for the finite contour length. Activity is implemented by a Gaussian, non-Markovian stochastic process resembling either an external nonthermal force or a local self-propulsion velocity as for an active Ornstein-Uhlenbeck particle. Specifically, the fluctuation spectrum of normal-mode amplitudes is analyzed. At elevated activities, flexible (tension) modes dominate over bending modes even for semiflexible rings, corresponding to enhanced conformational fluctuations. The fluctuation spectrum exhibits a crossover from a quadratic to a quartic dependence on the mode number with increasing mode number, originating from intramolecular tension, but the relaxation behavior is either dominated by intra-polymer processes or the active stochastic process. A further increase in activity enhances fluctuations at large length scales at the expense of reduced fluctuations at small scales. Conformationally, the mean square ring diameter exhibits swelling qualitatively comparable to liner polymers. The ring's diffusive dynamics is enhanced, and the mean square displacement shows distinct activity-determined regimes, consecutively, a ballistic, a subdiffusive, and a diffusive regime. The subdiffusive regime disappears gradually with increasing activity.

The effect of a nonspherical particle shape on the dynamics in crowded solutions presents a significant challenge for a comprehensive understanding of interaction and structural relaxation in biological and soft matter. We report that small deviations from a spherical shape induce a nonmonotonic contribution to the crowding effect on the short-time cage diffusion compared with spherical systems, using molecular dynamics simulations with mesoscale hydrodynamics of a multiparticle collision dynamics fluid in semidilute systems with volume fractions smaller than 0.35. We show that the nonmonotonic effect due to anisotropy is caused by the combination of a reduced relative mobility over the entire concentration range and a looser and less homogeneous cage packing of nonspherical particles. Our finding stresses that nonsphericity induces new complexity, which cannot be accounted for in effective sphere models, and is of great interest in applications such as formulations as well as for the fundamental understanding of soft matter in general and crowding effects in living cells in particular.

The spatiotemporal dynamics in systems of active self-propelled particles is controlled by the propulsion mechanism in combination with various direct interactions, such as steric repulsion and hydrodynamics. These direct interactions are typically anisotropic, and come in different "flavors", such as spherical and elongated particle shapes, pusher and puller flow fields, etc. The combination of the various aspects is expected to lead to new emergent behavior. However, it is a priori not evident whether shape and hydrodynamics act synergistically or antagonistically to generate motility-induced clustering (MIC) and phase separation (MIPS). We employ a model of prolate spheroidal microswimmers-called squirmers-in quasi-two-dimensional confinement to address this issue by mesoscale hydrodynamic simulations. For comparison, non-hydrodynamic active Brownian particles (ABPs) are considered to elucidate the contribution of hydrodynamic interactions. For spherical particles, the comparison between ABPs and hydrodynamic-squirmer ensembles reveals a suppression of MIPS due to hydrodynamic interactions. Yet, our analysis shows that dynamic clusters exist, with a broad size distribution. The fundamental difference between ABPs and squirmers is attributed to an increased reorientation of squirmers by hydrodynamic torques during their collisions. In contrast, for elongated squirmers, hydrodynamics interactions enhance MIPS. The transition to a phase-separated state strongly depends on the nature of the swimmer's flow field-with an increased tendency toward MIPS for pullers, and a reduced tendency for pushers. Thus, hydrodynamic interactions show opposing effects on MIPS for spherical and elongated microswimmers, and details of the propulsion mechanism of biological microswimmers may be very important to determine their collective behavior.

The conformational and rheological properties of active filaments/polymers exposed to shear flow are studied analytically. Using the continuous Gaussian semiflexible polymer model extended by the activity, we derive analytical expressions for the dependence of the deformation, orientation, relaxation times, and viscosity on the persistence length, shear rate, and activity. The model yields a Weissenberg-number dependent shear-induced deformation, alignment, and shear thinning behavior, similarly to the passive counterpart. Thereby, the model shows an intimate coupling between activity and shear flow. As a consequence, activity enhances the shear-induced polymer deformation for flexible polymers. For semiflexible polymers/filaments, a nonmonotonic deformation is obtained because of the activity-induced shrinkage at moderate and swelling at large activities. Independent of stiffness, activity-induced swelling facilitates and enhances alignment and shear thinning compared to a passive polymer. In the asymptotic limit of large activities, a polymer length- and stiffness-independent behavior is obtained, with universal shear-rate dependencies for the conformations, dynamics, and rheology.

Adaptive hydrogels, often termed smart materials, are macromolecules whose structure adjusts to external stimuli. Responsive micro- and nanogels are particularly interesting because the small length scale enables very fast response times. Chemical cross-links provide topological constraints and define the three-dimensional structure of the microgels, whereas their porous structure permits fast mass transfer, enabling very rapid structural adaption of the microgel to the environment. The change of microgel structure involves a unique transition from a flexible, swollen finite-size macromolecular network, characterized by a fuzzy surface, to a colloidal particle with homogeneous density and a sharp surface. In this contribution, we determine, for the first time, the structural evolution during the microgel-to-particle transition. Time-resolved small-angle x-ray scattering experiments and computer simulations unambiguously reveal a two-stage process: In a first, very fast process, collapsed clusters form at the periphery, leading to an intermediate, hollowish core-shell structure that slowly transforms to a globule. This structural evolution is independent of the type of stimulus and thus applies to instantaneous transitions as in a temperature jump or to slower stimuli that rely on the uptake of active molecules from and/or exchange with the environment. The fast transitions of size and shape provide unique opportunities for various applications as, for example, in uptake and release, catalysis, or sensing.

The structural and dynamical properties of ultra-soft colloids-star polymers-exposed to a uniform external force field are analyzed by applying the multiparticle collision dynamics technique, a hybrid coarse-grain mesoscale simulation approach, which captures thermal fluctuations and long-range hydrodynamic interactions. In the weak-field limit, the structure of the star polymer is nearly unchanged; however, in an intermediate regime, the radius of gyration decreases, in particular transverse to the sedimentation direction. In the limit of a strong field, the radius of gyration increases with field strength. Correspondingly, the sedimentation coefficient increases with increasing field strength, passes through a maximum, and decreases again at high field strengths. The maximum value depends on the functionality of the star polymer. High field strengths lead to symmetry breaking with trailing, strongly stretched polymer arms and a compact star-polymer body. In the weak-field-linear response regime, the sedimentation coefficient follows the scaling relation of a star polymer in terms of functionality and arm length.

Dispersions of particles with short-range attractive and long-range repulsive interactions exhibit rich equilibrium microstructures and a complex phase behavior. We present theoretical and simulation results for structural and, in particular, short-time diffusion properties of a colloidal model system with such interactions, both in the dispersed-fluid and equilibrium-cluster phase regions. The particle interactions are described by a generalized Lennard-Jones-Yukawa pair potential. For the theoretical-analytical description, we apply the hybrid Beenakker-Mazur pairwise additivity (BM-PA) scheme. The static structure factor input to this scheme is calculated self-consistently using the Zerah-Hansen integral equation theory approach. In the simulations, a hybrid simulation method is adopted, combing molecular dynamics simulations of colloids with the multiparticle collision dynamics approach for the fluid, which fully captures hydrodynamic interactions. The comparison of our theoretical and simulation results confirms the high accuracy of the BM-PA scheme for dispersed-fluid phase systems. For particle attraction strengths exceeding a critical value, our simulations yield an equilibrium cluster phase. Calculations of the mean lifetime of the appearing clusters and the comparison with the analytical prediction of the dissociation time of an isolated particle pair reveal quantitative differences pointing to the importance of many-particle hydrodynamic interactions for the cluster dynamics. The cluster lifetime in the equilibrium-cluster phase increases far stronger with increasing attraction strength than that in the dispersed-fluid phase. Moreover, significant changes in the cluster shapes are observed in the course of time. Hence, an equilibrium-cluster dispersion cannot be treated dynamically as a system of permanent rigid bodies.

Colloidal particles subject to an external periodic forcing exhibit complex collective behavior and self-assembled patterns. A dispersion of magnetic microparticles confined at the air-liquid interface and energized by a uniform uniaxial alternating magnetic field exhibits dynamic arrays of self-assembled spinners rotating in either direction. Here, we report on experimental and simulation studies of active turbulence and transport in a gas of self-assembled spinners. We show that the spinners, emerging as a result of spontaneous symmetry breaking of clock/counterclockwise rotation of self-assembled particle chains, generate vigorous vortical flows at the interface. An ensemble of spinners exhibits chaotic dynamics due to self-generated advection flows. The same-chirality spinners (clockwise or counterclockwise) show a tendency to aggregate and form dynamic clusters. Emergent self-induced interface currents promote active diffusion that could be tuned by the parameters of the external excitation field. Furthermore, the erratic motion of spinners at the interface generates chaotic fluid flow reminiscent of 2D turbulence. Our work provides insight into fundamental aspects of collective transport in active spinner materials and yields rules for particle manipulation at the microscale., Competing Interests: The authors declare no conflict of interest., (Copyright © 2017 the Author(s). Published by PNAS.)

Hydrodynamic fluctuations in simple fluids under shear flow are demonstrated to be spatially correlated, in contrast to the fluctuations at equilibrium, using mesoscopic hydrodynamic simulations. The simulation results for the equal-time hydrodynamic correlations in a multiparticle collision dynamics (MPC) fluid in shear flow are compared with the explicit expressions obtained from fluctuating hydrodynamics calculations. For large wave vectors k, the nonequilibrium contributions to transverse and longitudinal velocity correlations decay as k^{-4} for wave vectors along the flow direction and as k^{-2} for the off-flow directions. For small wave vectors, a crossover to a slower decay occurs, indicating long-range correlations in real space. The coupling between the transverse velocity components, which vanishes at equilibrium, also exhibits a k^{-2} dependence on the wave vector. In addition, we observe a quadratic dependency on the shear rate of the nonequilibrium contribution to pressure.

The morphology of core-shell microgels under different swelling conditions and as a function of the core-shell thickness ratio is systematically characterized by mesoscale hydrodynamic simulations. With increasing hydrophobic interaction of the shell polymers, we observe drastic morphological changes from a core-shell structure to an inverted microgel, where the core is turned to the outside, or a microgel with a patchy surface of core polymers directly exposed to the environment. We establish a phase diagram of the various morphologies. Moreover, we characterize the polymer and microgel conformations. For sufficiently thick shells, the changes of the shell size upon increasing hydrophobic interactions are well described by the Flory-Rehner theory. Additionally, this theory provides a critical line in the phase diagram separating core-shell structures from the distinct two other phases. The appearing new phases provide a novel route to nano- and microscale functionalized materials.

The intramolecular dynamics of flexible and semiflexible polymers in response to active noise is studied theoretically. The active noise may either originate from interactions of a passive polymer with a bath of active Brownian particles or the polymer itself is comprised of active Brownian particles. We describe the polymer by the continuous Gaussian semiflexible-polymer model, taking into account the finite polymer extensibility. Our analytical calculations predict a strong dependence of the polymer dynamics on the activity. In particular, active semiflexible polymers exhibit a crossover from a bending elasticity-dominated dynamics at weak activity to that of flexible polymers at strong activity. The end-to-end vector correlation function decays exponentially for times longer than the longest polymer relaxation time. Thereby, the polymer relaxation determines the decay of the correlation function for long and flexible polymers. For shorter and stiffer polymers, the relaxation behavior of individual active Brownian particles dominates the decay above a certain activity. The diffusive dynamics of a polymer is substantially enhanced by the activity. Three regimes can be identified in the mean square displacement for sufficiently strong activities: an activity-induced ballistic regime at short times, followed by a Rouse-type polymer-specific regime for any polymer stiffness, and free diffusion at long times, again determined by the activity.

We investigate the conformational properties of polymers in ionic microgels in the presence of salt ions by molecular dynamics simulations and analytical theory. A microgel particle consists of coarse-grained linear polymers, which are tetra-functionally crosslinked. Counterions and salt ions are taken into account explicitly, and charge-charge interactions are described by the Coulomb potential. By varying the charge interaction strength and salt concentration, we characterize the swelling of the polyelectrolytes and the charge distribution. In particular, we determine the amount of trapped mobile charges inside the microgel and the Debye screening length. Moreover, we analyze the polymer extension theoretically in terms of the tension blob model taking into account counterions and salt ions implicitly by the Debye⁻Hückel model. Our studies reveal a strong dependence of the amount of ions absorbed in the interior of the microgel on the electrostatic interaction strength, which is related to the degree of the gel swelling. This implies a dependence of the inverse Debye screening length κ on the ion concentration; we find a power-law increase of κ with the Coulomb interaction strength with the exponent 3 / 5 for a salt-free microgel and an exponent 1 / 2 for moderate salt concentrations. Additionally, the radial dependence of polymer conformations and ion distributions is addressed.

We analyze the dynamics of polymers in a microgel system under different swelling conditions. A microgel particle consists of coarse-grained linear polymers which are tetra-functionally crosslinked and undergoes conformational changes in response to the external stimuli. Here, a broad range of microgel sizes, extending from tightly collapsed to strongly swollen particles, is considered. In order to account for hydrodynamic interactions, the microgel is embedded in a multiparticle collision dynamics fluid while hydrophobic attraction is modelled by an attractive Lennard-Jones potential and swelling of ionic microgels is described through the Debye-Hückel potential. The polymer dynamics is analyzed in terms of the monomer mean square displacement and the intermediate scattering function S(q, t). The scattering function decays in a stretched-exponential manner, with a decay rate exhibiting a crossover from a collective diffusive dynamics at low magnitudes of the wavevector q to a hydrodynamic-dominated dynamics at larger q. There is little difference between the intermediate scattering functions of microgels under good solvent conditions and strongly swollen gels, but strongly collapsed gels exhibit a faster decay at short times and hydrodynamic interactions become screened. In addition, we present results for the dynamics of the crosslinks, which exhibit an unexpected, semiflexible polymer-like dynamics.

In the dense and crowded environment of the cell cytoplasm, an individual protein feels the presence of and interacts with all surrounding proteins. While we expect this to strongly influence the short-time diffusion coefficient D s of proteins on length scales comparable to the nearest-neighbor distance, this quantity is difficult to assess experimentally. We demonstrate that quantitative information about D s can be obtained from quasi-elastic neutron scattering experiments using the neutron spin echo technique. We choose two well-characterized and highly stable eye lens proteins, bovine α-crystallin and γ B -crystallin, and measure their diffusion at concentrations comparable to those present in the eye lens. While diffusion slows down with increasing concentration for both proteins, we find marked variations that are directly linked to subtle differences in their interaction potentials. A comparison with computer simulations shows that anisotropic and patchy interactions play an essential role in determining the local short-time dynamics. Hence, our study clearly demonstrates the enormous effect that weak attractions can have on the short-time diffusion of proteins at concentrations comparable to those in the cellular cytosol.

A wide spectrum of Peritrichous bacteria undergo considerable physiological changes when they are inoculated onto nutrition-rich surfaces and exhibit a rapid and collective migration denoted as swarming. Thereby, the length of such swarmer cells and their number of flagella increases substantially. In this article, we investigated the properties of individual E. coli-type swarmer cells confined between two parallel walls via mesoscale hydrodynamic simulations, combining molecular dynamics simulations of the swarmer cell with the multiparticle particle collision dynamics approach for the embedding fluid. E. coli-type swarmer cells are three-times longer than their planktonic counter parts, but their flagella density is comparable. By varying the wall separation, we analyze the confinement effect on the flagella arrangement, on the distribution of cells in the gap between the walls, and on the cell dynamics. We find only a weak dependence of confinement on the bundle structure and dynamics. The distribution of cells in the gap changes from a geometry-dominated behavior for very narrow to fluid-dominated behavior for wider gaps, where cells are preferentially located in the gap center for narrower gaps and stay preferentially next to one of the walls for wider gaps. Dynamically, the cells exhibit a wide spectrum of migration behaviors, depending on their flagella bundle arrangement, and ranges from straight swimming to wall rolling.

We propose a hydrodynamic model for a spheroidal microswimmer with two tangential surface velocity modes. This model is analytically solvable and reduces to Lighthill's and Blake's spherical squirmer model in the limit of equal major and minor semi-axes. Furthermore, we present an implementation of such a spheroidal squirmer by means of particle-based mesoscale hydrodynamics simulations using the multiparticle collision dynamics approach. We investigate its properties as well as the scattering of two spheroidal squirmers in a slit geometry. Thereby we find a stable fixed point, where two pullers swim cooperatively forming a wedge-like conformation with a small constant angle.

The conformational properties of flexible and semiflexible polymers exposed to active noise are studied theoretically. The noise may originate from the interaction of the polymer with surrounding active (Brownian) particles or from the inherent motion of the polymer itself, which may be composed of active Brownian particles. In the latter case, the respective monomers are independently propelled in directions changing diffusively. For the description of the polymer, we adopt the continuous Gaussian semiflexible polymer model. Specifically, the finite polymer extensibility is taken into account, which turns out to be essential for the polymer conformations. Our analytical calculations predict a strong dependence of the relaxation times on the activity. In particular, semiflexible polymers exhibit a crossover from a bending elasticity-dominated dynamics to the flexible polymer dynamics with increasing activity. This leads to a significant activity-induced polymer shrinkage over a large range of self-propulsion velocities. For large activities, the polymers swell and their extension becomes comparable to the contour length. The scaling properties of the mean square end-to-end distance with respect to the polymer length and monomer activity are discussed.

The dynamical properties of a flexible dumbbell composed of active Brownian particles are analytically analyzed. The dumbbell is considered as a simplified description of a linear active polymer. The two beads are independently propelled in directions which change in a diffusive manner. The relaxation behavior of the internal degree of freedom is tightly coupled to the dumbbell activity. The latter dominates the dynamics for strong propulsion. As is shown, limitations in bond stretching strongly influence the relaxation behavior. Similarly, under shear flow, activity determines the relaxation and tumbling behavior at strong propulsion. Moreover, shear leads to a preferred alignment and consequently to shear thinning. Thereby, a different power-law dependence on the shear rate compared to passive dumbbells under flow is found.

The friction and diffusion coefficients of rigid spherical colloidal particles dissolved in a fluid are determined from velocity and force autocorrelation functions by mesoscale hydrodynamic simulations. Colloids with both slip and no-slip boundary conditions are considered, which are embedded in fluids modeled by multiparticle collision dynamics with and without angular momentum conservation. For no-slip boundary conditions, hydrodynamics yields the well-known Stokes law, while for slip boundary conditions the lack of angular momentum conservation leads to a reduction of the hydrodynamic friction coefficient compared to the classical result. The colloid diffusion coefficient is determined by integration of the velocity autocorrelation function, where the numerical result at shorter times is combined with the theoretical hydrodynamic expression for longer times. The suitability of this approach is confirmed by simulations of sedimenting colloids. In general, we find only minor deviations from the Stokes-Einstein relation, which even disappear for larger colloids. Importantly, for colloids with slip boundary conditions, our simulation results contradict the frequently assumed additivity of local and hydrodynamic diffusion coefficients.

Polyelectrolyte gels are known to undergo significant conformational changes in response to external stimuli such as pH, temperature, or the dielectric constant. Specifically, an increase of the degree of ionization associated with an increasing number of counterions leads to swelling of the network. For a macroscopically large gel, which is electrostatically neutral in its interior, swelling is no longer governed by electrostatic interactions, but rather by the osmotic pressure of counterions. However, this electrostatic neutrality is typically violated for nanogels, because counterions are free to leave a gel particle. Although nanogel-swelling exhibits similar features as swelling of micro- and macrogels, another mechanism has to be relevant. Here, we use molecular dynamics simulations and scaling theory to unravel the structural properties of nanogels upon changing the electrostatic interactions. We demonstrate that the swelling of nanogels is governed by screened electrostatic interactions without a relevant contribution by the counterion osmotic pressure.

The nonequilibrium dynamical behavior and structure formation of end-functionalized semiflexible polymer suspensions under flow are investigated by mesoscale hydrodynamic simulations. The hybrid simulation approach combines the multiparticle collision dynamics method for the fluid, which accounts for hydrodynamic interactions, with molecular dynamics simulations for the semiflexible polymers. In equilibrium, various kinds of scaffold-like network structures are observed, depending on polymer flexibility and end-attraction strength. We investigate the flow behavior of the polymer networks under shear and analyze their nonequilibrium structural and rheological properties. The scaffold structure breaks up and densified aggregates are formed at low shear rates, while the structural integrity is completely lost at high shear rates. We provide a detailed analysis of the shear- rate-dependent flow-induced structures. The studies provide a deeper understanding of the formation and deformation of network structures in complex materials.