Crossover from normal to anomalous diffusion in field-aligned dipolar systems
Using molecular dynamics simulations we investigate the translational dynamics of particles with dipolar interactions in homogenous external fields. For a broad range of concentrations, we find that the anisotropic, yet normal diffusive behavior characterizing weakly coupled systems becomes anomalous both parallel and perpendicular to the field at sufficiently high dipolar coupling and field strength. After the ballistic regime, chain formation first yields cage-like motion in all directions, followed by transient, mixed diffusive-superdiffusive behavior resulting from cooperative motion of the chains. The enhanced dynamics disappears only at higher densities close to crystallization.
There is a strong current interest in colloidal soft-matter systems displaying anomalous dynamic behavior beyond simple diffusion. A paradigm example are colloidal hard spheres which exhibit a (repulsive) glass transition in near-perfect agreement with Mode Coupling theory Pusey87 (). Recent attention focuses, on one hand, on the anomalous slow dynamics in a broad variety of systems such as low-density colloidal gels Lu08 (), cells Lebold09 (), driven granular fluids Fiege09 (), and disordered solids Kurzidim09 (). In many cases the reduced mobility is accompanied by dynamic heterogeneities and cooperative motion as recently shown, e.g., in colloids Kob97 (); Chaud07 () and in liquid crystals Patti09 (). On the other hand, there is increasing interest in systems exhibiting enhanced mobility such as rods between obstacles Hoefling08 () and activated systems such as catalytic nanorods Dhar06 () or self-propelled bacteria Golestian09 ().
In this letter we discuss the translational dynamics of dipolar systems which, as we will demonstrate, can exhibit both, slow and enhanced mobility. Prominent examples of dipole-coupled systems are suspensions of magnetic nanoparticles (ferrocolloids), polarizable colloids Leu09 () or colloidal ”molecules” (e.g. dumbells) with partial charges Mano03 (), but also electric field-driven living cells Gupta09 (). The strong current attention received by dipolar systems (or, more generally, systems with directional interactions) is born by their importance for self-assembly processes towards novel classes of equilibrium and non-equilibrium materials Gangwal10 (); Zhang07 (). In zero field, the anisotropic interactions between permanent dipoles can induce network formation Tlusty00 () including the possibility of an equilibrium gel Blaak07 (). On the contrary, the dynamics in an external field appears far less understood, although this situation is clearly most relevant for many practical applications of dipolar systems Odenbach (). Previous experimental and theoretical studies indicate that the single-particle motion under field can be strongly anisotropic Ilg05 () and even sub-diffusive (in lateral directions) due to ”caging” of the particles in the chains formed at sufficiently high coupling strength Furst00 (); Toussaint04 (); Belovs06 (). However, only recently dynamic light-scattering experiments Mertelj09 () reveal, in addition, strong evidence for superdiffusive motion and dynamic heterogeneity. Motivated by these exciting and seemingly somewhat contradictory findings, we have performed extensive molecular dynamics (MD) simulations of a simple dipolar model system. Analyzing separately the microscopic motion parallel and perpendicular to the field for a broad range of parameters we find, for the first time, a coupling strength-induced crossover from normal to anomalous dynamics. The latter is characterized by a step-wise behavior involving sub-diffusion, directed motion and dynamic heterogeneities.
Our model fluid consists of dipolar soft spheres (DSS) of diameter with embedded, permanent point dipole moments . The total pair interaction at distance consists of a (truncated and shifted) soft-sphere repulsion, , and the long-range, anisotropic dipole-dipole potential, . We employ constant-temperature MD simulations with particles. Thus, we focus on the direct (conservative) effects on the dynamics rather than on those stemming from a solvent, in accordance with other recent studies (see, e.g., Blaak07 ()). The dipolar interactions are treated using the Ewald method with conducting boundaries (see Jordanov09 () for further details). Simulations are carried out at a reduced density (which is typical for the experiments in Mertelj09 ()) and a reduced temperature (with and being Boltzmann’s constant and true temperature, respectively). The dipole moment is chosen such that the resulting coupling parameter has values , consistent with the range considered in recent ferrofluid experiments Mertelj09 (); klokkenburg06 (). The fluids are subject to a homogeneous magnetic field with fields strengths up to the (still realistic klokkenburg06 ()) value .
One key quantity of our analysis is the translational mean-squared displacement (MSD), . Representative MSDs for a moderate coupling parameter () are shown in Fig. 1a).
For we differ between the MSDs parallel and perpendicular to the field direction. This is sensible due to the structural anisotropy resulting from the formation of elongated clusters along the field, see Fig. 1b). Despite these (small) clusters, all of systems have ”normal” translational dynamics in the sense that there is an initial ballistic regime, , after which the MSD directly goes over into a diffusive regime . The latter, Einstein-like behavior allows us to define diffusion constants (for ) and , () in the standard way. Of particular interest is the anisotropy ratio , results for which are plotted in Fig. 1c). At the coupling strengths considered, motion perpendicular to is ”faster” than along the field, i.e., (consistent with earlier simulation results Ilg05 ()), despite the fact that the clusters form in -direction. Thus the lateral translational fluctuations of the particles are still so large that their diffusive behavior resembles disks than rods (where one expects ). However, as Fig. 1c) clearly indicates, there is a crossover towards different dynamic behavior as increases beyond .
Typical MSD’s for this strongly coupled regime are plotted in Fig. 2 along with corresponding simulation snapshots ().
At zero field, the particles arrange into a globally isotropic network of worm-like chains [Fig. 2c)], as it is typical for dilute, strongly interacting dipolar systems Tlusty00 (); Blaak07 (). The corresponding MSD still behaves ”normal” in the sense described above, consistent with MD results for dipolar dumbbells Blaak07 () at comparable values of . With field, long chains form which become more and more straight upon increase of . We first consider the case [see Figs. 2b) and e)]. After the initial ballistic regime, the motion of a particle significantly slows down due to the confinement by its nearest neighbors, yielding a sub-diffusive regime in the MSD characterized by a time dependence with . This holds not only in -direction (where sub-diffusion occurs for with ), but also perpendicular to (, ), indicating pronounced spatial-temporal correlations in all spatial directions. We characterize these correlations by the distinct parts of the van Hove correlation functions and with and being the longitudinal and lateral distance, respectively. Results are plotted in Figs. 3a) and b). For the functions coincide with the corresponding static pair distribution functions.
The data for at clearly reflect long-range correlations, and thus, a ”trapping” of particles along the chain, which starts to die out only at times beyond the sub-diffusive regime of the MSD (i.e., ). The exponent within the sub-diffusive regime is close to the value of found in the MSD of interior beads in linear polymer chains Baschnagel00 () prior to crossover to free diffusion. In perpendicular directions [see Fig. 3b)], the correlation ”hole” at reflects that positions directly besides a particle are disfavored due to the repulsive dipolar interactions associated to this configuration. At later times, however, we observe development of a peak at , indicating that the original position of a particle is replaced by another particle. These strong and long-living correlations explain the sub-diffusive behavior of the perpendicular MSD, characterized by an exponent . A similar exponent () has been found in experiments probing lateral diffusion in magnetorheological systems Furst00 (); Toussaint04 () involving chains of induced (and thus, perfectly aligned) dipoles. This similarity already indicates that the strong suppression of orientational fluctuations, i.e., a strong field, is an important ingredient for the observed sub-diffusive behavior in our present system consisting of permanent dipoles. Indeed, the MSD’s at the somewhat weaker field [see Fig. 2a),d)] do not reveal a sub-diffusive regime. Instead, one merely observes a slowing down (, ) as compared to the initial ballistic regime.
However, both anisotropic systems ( and ) display, at , a further intriguing feature not present in weakly coupled systems considered in Fig. 1. Specifically, we observe in Figs. 2a)-b) a ”kink” in the MSDs, which separate the intermediate time range discussed above and a later time regime where the MSDs grow faster than linearly in time. This behavior strongly contrasts that in other disordered systems (see, e.g., Baschnagel00 ()), where the sub-diffusive regime, if present, changes directly into true diffusion (). In our present system, the initial exponent within the ”fast” regime (after the kink) lies between and ; however, these values decrease in time and also strongly depend on , as does the time related to the kink. We conclude that there is a (transient) superposition of diffusive and ballistic () behavior, with a tendency towards normal diffusion (visible particularly in the lateral part) at large times. In fact, at the low density considered () we did not reach the true diffusive regime particularly for the longitudinal part in the total simulation time (contrary to what is found at larger , see Fig. 5). To test possible finite-size effects we have repeated the simulations with particles, the changes [see Fig. 2b)], however, being only marginal.
A direct illustration of the superposition of diffusive and ballistic motion is provided by typical trajectories plotted in the bottom of Fig. 3. In -direction [3c)], one clearly recognizes small random-walk-like movements of the particles around a mean path consisting of essentially straight stretches, the latter giving rise to the ”fast” regime of the parallel MSD [see Fig. 2]. Similar random movement around a mean path occurs perpendicular to [Fig. 3d)], where, however, the stretches are shorter and, consequently, the growth of the MSD (after the kink) is somewhat less pronounced. Moreover, in both directions, the motion of two particles which are nearest neighbors at is correlated over an extremely long time, consistent with the behavior of (and ) shown in the top of Fig. 3. These observations suggest that, within the fast regime, the particles move collectively as chains (indeed, at the large field considered, there are essentially no unbounded particles). Moreover, the stretches in the trajectories indicate that the entire chains can move quasi-independently (i.e., ballistically) over a finite time interval before they are slowed down due to interactions with particles in other chains.
The occurrence of cooperative motion of the particles as chains is further confirmed by the behavior of the self parts of the van Hove functions, and measuring the probability that a particle moves over a certain distance in a time . Exemplary results for time (i.e., within the ”fast” regime) are shown in Fig. 4.
Particularly the parallel part strongly deviates from the usual Gaussian shape characterizing diffusive systems. Similar deviations occur at other times (and in the perpendicular part), yielding large values of the ”non-Gaussian”-parameter (see insets in Fig. 4). At , specifically, the buckled structure of [Fig. 4a)] suggests that the system consists of several portions (”populations”) of particles. To better identify these populations we have resolved the instantaneous data for with respect to the actual particles in the system as illustrated in Fig. 4c),d). One clearly recognizes that particles related to one ”population” in belong to the same chain in the actual configuration. Also, comparison of Figs. 4c) and d) shows that chains which are fast in -direction move more slowly in perpendicular directions and vice versa.
Taken altogether, we interpret the ”kink” in the MSDs (and thus, the superdiffusive behavior), as a consequence of the motion of the entire chains, which over certain times, interact only weakly with particles of other chains. This interpretation prompts the question what would happen at higher densities, where chain-chain interactions should be more prominent. To this end, we plot in Fig. 5 results for MSD’s in the range . The simulations have been performed at the same strong coupling conditions as before, yielding essentially percolated chains in -direction.
An inspection of the data in Fig. 5 reveals that the ”kink” indicating onset of (transient) superdiffusive motion along survives up to in field direction; at higher one observes normal dynamics where the initial ballistic regime directly changes into diffusion. Perpendicular to , the disappearance of the kink occurs already ; at plateau-like regions in the MSD appear. Indeed, at high densities one expects Tao91 () formation of body-centered-tetragonal structures with neighboring chains being out of registry, a tendency which is visible in the snapshot in Fig. 5c). Clearly, this lateral ordering induces additional dynamic effects not present at small .
To summarize, our MD simulation results demonstrate that chain formation of strongly coupled, dilute or moderately dense dipolar fluids in external fields gives rise to anomalous translational dynamics. The latter involves both, sub-diffusive behavior at intermediate times and mixed diffusive-superdiffusive behavior at later times, consistent with recent ferrofluid experiments Mertelj09 (). According to our analysis in the superdiffusive regime, the particles move collectively as chains, a behavior which appears as dynamic heterogeneity on the single-particle level. From a fundamental point of view, a particularly intriguing result of our study is that the systems may display enhanced dynamics already in equilibrium, contrary to active (bio-) materials such as self-propelled bacteria Golestian09 () and catalytic nanorods Dhar06 (), which are intrinsically out of equilibrium. From an application perspective, we expect our results to be relevant not only for ferrofluids, but for the much broader class of colloidal systems where external fields induce directed self-assembly, examples being patchy colloids such as metallodielectric (Janus) particles Gangwal10 () and living cells suspensions Gupta09 (). Clearly more detailed experiments are needed to confirm our findings, such as the ”kink” separating the dynamic regimes in the MSD. On the theoretical side, points deserving further attention are the role of orientational fluctuations at moderate fields, as well as solvent (hydrodynamic) effects which are so far neglected. Work in these directions is in progress.
Acknowledgements.Financial support from the German Science Foundation via grant KL1215/6 is gratefully acknowledged.
- (1) P. N. Pusey and W. van Megen, Phys. Rev. Lett. 59, 2083 (1987).
- (2) P. J. Lu, E. Zaccarelli, F. Ciulla, A. B. Schofield, F. Sciortino, D. A. Weitz, Nature 453, 499 (2008).
- (3) T. Lebold, C. Jung, J. Michaelis, and C. Bräuchle, Nanolett. 9, 2877 (2009).
- (4) A. Fiege, T. Aspelmeier, and A. Zippelius, Phys. Rev. Lett. 102, 098001 (2009).
- (5) J. Kurzidim, D. Coslovich, and G. Kahl, Phys. Rev. Lett. 103, 138303 (2009).
- (6) W. Kob, C. Donati, S. J. Plimpton, P. H. Poole, and S. C. Glotzer, Phys. Rev. Lett. 79, 2827 (1997); C. Donati, J. F. Douglas, W. Kob, S. J. Plimpton, P. H. Poole, and S. C. Glotzer, Phys. Rev. Lett. 80, 2338 (1997)
- (7) P. Chaudhuri, L. Berthier, and W. Kob, Phys. Rev. Lett. 99, 060604 (2007).
- (8) A. Patti, D. El Masri, R. van Roij, and M. Dijkstra, Phys. Rev. Lett. 103, 248304 (2009); M. P. Lettinga and E. Grelet, Phys. Rev. Lett. 99, 197802 (2007).
- (9) F. Höfling, E. Frey, and T. Franosch, Phys. Rev. Lett. 101, 120605 (2008).
- (10) P. Dhar, Th. M. Fischer, Y. Wang, T. E. Mallouk, W. F. Paxton, and A. Sen, Nanolett. 6, 66 (2006).
- (11) R. Golestanian, Phys. Rev. Lett. 102, 188305 (2009).
- (12) M. E. Leunissen, H. R. Vutukuri, and A. van Blaaderen, Adv. Mat. 21, 3116 (2009).
- (13) V. N. Manohoran, M. T. Elsesser, and D. J. Pine, Science 301, 483 (2003).
- (14) S. Gupta et al., Langmuir 26, 3441 (2009).
- (15) S. Gangwal, A. Pawar, I. Kretzschmar, and O. D. Velev, Soft Matter 6, 1413 (2010).
- (16) Z. Zhang, Z. Tang, N. A. Kotov, and S. C. Glotzer, Nanolett. 7, 1670 (2007).
- (17) T. Tlusty and S. A. Safran, Science 290, 1328 (2000).
- (18) R. Blaak, M. A. Miller, and J.-P. Hansen, Europhys. Lett. 78, 26002 (2007).
- (19) S. Odenbach (Ed.), Colloidal Magnetic Fluids: Basics, Development and Application of Ferrofluids, Lect. Notes Phys. 763 (Springer, Berlin Heidelberg 2009).
- (20) P. Ilg and M. Kröger, Phys. Rev. E 72, 031504 (2005).
- (21) E. M. Furst and A. P. Gast, Phys. Rev. E 62, 6916 (2000).
- (22) R. Toussaint, G. Helgesen, and E. G. Flekkoy, Phys. Rev. Lett. 93, 108304 (2004).
- (23) M. Belovs and A. Cebers, Phys. Rev. E 73, 021507 (2006).
- (24) A. Mertelj, L. Cmok, and M. Copic, Phys. Rev. E 79, 041402 (2009).
- (25) J. Jordanovic and S. H. L. Klapp, Phys. Rev. E 79, 021405 (2009).
- (26) M. Klokkenburg et al., Phys. Rev. Lett. 97, 185702 (2006).
- (27) J. Baschnagel, C. Bennemann, W. Paul, and K. Binder, J. Phys.: Condensed Matter 12, 6365 (2000).
- (28) R. Tao and J. M. Sun, Phys. Rev. Lett. 67, 398 (1991).