No evidence of gas-liquid coexistence in dipolar hard spheres

Percolation transition and bimodal density distribution in hydrogen fluoride

Hydrogen-bond networks in associating fluids can be extremely robust and characterize the topological properties of the liquid phase, as in the case of water, over its whole domain of stability and beyond. Here, we report on molecular dynamics simulations of hydrogen fluoride (HF), one of the strongest hydrogen-bonding molecules. HF has more limited connectivity than water but can still create long, dynamic chains, setting it apart from most other small molecular liquids. Our simulation results provide robust evidence of a second-order percolation transition of HF's hydrogen bond network occurring below the critical point. This behavior is remarkable as it underlines the presence of two different cohesive mechanisms in liquid HF, one at low temperatures characterized by a spanning network of long, entangled hydrogen-bonded polymers, as opposed to short oligomers bound by the dispersion interaction above the percolation threshold. This second-order phase transition underlines the presence of marked structural heterogeneity in the fluid, which we found in the form of two liquid populations with distinct local densities.

Multicore-based ferrofluids in zero field: initial magnetic susceptibility and self-assembly mechanisms

The necessity to improve magnetic building blocks in magnetic nano-structured soft materials stems from a fascinating potential these materials have in bio-medical applications and nanofluidics. Along with practical reasons, the interplay of magnetic and steric interactions on one hand, and entropy, on the other, makes magnetic soft matter fundamentally challenging. Recently, in order to tailor magnetic response of the magnetic particle suspensions, the idea arose to replace standard single-core nanoparticles with nano-sized clusters of single-domain nanoparticles (grains) rigidly bound together by solid polymer matrix - multicore magnetic nanoparticles (MMNPs). To pursue this idea, a profound understanding of the MMNP interactions and self-assembly is required. In this work we present a computational study of the MMNP suspensions and elucidate their self-assembly and magnetic susceptibility. We show that depending on the magnetic moment of individual grains the suspensions exhibit qualitatively distinct regimes. Firstly, if the grains are moderately interacting, they contribute to a significant decrease of the remanent magnetisation of MMNPs and as such to a decrease of the magnetic susceptibility, this way confirming previous findings. If the grains are strongly interacting, instead, they serve as anchor points and support formation of grain clusters that span through several MMNPs, leading to MMNP cluster formation and a drastic increase of the initial magnetic response. Both the topology of the clusters and their size distribution in MMNP suspensions is found to be notably different from those formed in conventional magnetic fluids or magnetorheological suspensions.

Inhomogeneous steady shear dynamics of a three-body colloidal gel former

We investigate the stationary flow of a colloidal gel under an inhomogeneous external shear force using adaptive Brownian dynamics simulations. The interparticle forces are derived from the Stillinger-Weber potential, where the three-body term is tuned to enable network formation and gelation in equilibrium. When subjected to the shear force field, the system develops remarkable modulations in the one-body density profile. Depending on the shear magnitude, particles accumulate either in quiescent regions or in the vicinity of maximum net flow, and we deduce this strong non-equilibrium response to be characteristic of the gel state. Studying the components of the internal force parallel and perpendicular to the flow direction reveals that the emerging flow and structure of the stationary state are driven by significant viscous and structural superadiabatic forces. Thereby, the magnitude and nature of the observed non-equilibrium phenomena differ from the corresponding behavior of simple fluids. We demonstrate that a simple power functional theory reproduces accurately the viscous force profile, giving a rationale of the complex dynamical behavior of the system.

The influence of anisotropy on the microstructure and magnetic properties of dipolar nanoplatelet suspensions

In recent years, there has been an increasing interest in magnetic nanoparticles with non-spherical shapes. This is largely due to their broad span of tuneable properties, which allow for tailoring of the colloidal properties by altering the magnetic anisometry or shape anisotropy of the nanoparticle. Although extensive research has gone into novel synthesis methods, the theoretical and analytical treatment of magnetic colloidal suspensions still predominantly focuses on spherical particles. This paper explores the microstructure and initial static magnetic susceptibility of systems of anisometric dipolar magnetic nanoplatelets in order to understand the applicability of dipolar sphere-based theories and models for such systems. We find that the microstructure as characterized by the particle distribution and magnetic clustering of platelets diverges significantly from that of spheres both quantitatively and qualitatively. We find lower initial static magnetic susceptibilities in nanoplatelet systems than in comparable suspensions of dipolar spheres. At lower values of the magnetic coupling constant, this can be accounted for by applying corrections to the volume fraction. However, this approach is less accurate for systems with stronger magnetic interactions. By providing predictions of and explanations for the observed effects, we aim to facilitate the use of magnetic nanoplatelet suspensions in the broad range of applications.

Emergent vortices and phase separation in systems of chiral active particles with dipolar interactions

Using Brownian dynamics (BD) simulations we investigate the self-organization of a monolayer of chiral active particles with dipolar interactions. Each particle is driven by both, translational and rotational self-propulsion, and carries a permanent point dipole moment at its center. The direction of the translational propulsion for each particle is chosen to be parallel to its dipole moment. Simulations are performed at high dipolar coupling strength and a density below that related to motility-induced phase separation in simple active Brownian particles. Despite this restriction, we observe a wealth of phenomena including formation of two types of vortices, phase separation, and flocking transitions. To understand the appearance and disappearance of vortices in the many-particle system, we further investigate the dynamics of simple ring structures under the impact of self-propulsion.

How chains and rings affect the dynamic magnetic susceptibility of a highly clustered ferrofluid

The dynamic magnetic susceptibility, χ(ω), of a model ferrofluid at a very low concentration (volume fraction, approximately 0.05%), and with a range of dipolar coupling constants (1≤λ≤8), is examined using Brownian dynamics simulations. With increasing λ, the structural motifs in the system change from unclustered particles, through chains, to rings. This gives rise to a nonmonotonic dependence of the static susceptibility χ(0) on λ and qualitative changes to the frequency spectrum. The behavior of χ(0) is already understood, and the simulation results are compared to an existing theory. The single-particle rotational dynamics are characterized by the Brownian time, τ_{B}, which depends on the particle size, carrier-liquid viscosity, and temperature. With λ≤5.5, the imaginary part of the spectrum, χ^{''}(ω), shows a single peak near ω∼τ_{B}^{-1}, characteristic of single particles. With λ≥5.75, the spectrum is dominated by the low-frequency response of chains. With λ≥7, new features appear at high frequency, which correspond to intracluster motions of dipoles within chains and rings. The peak frequency corresponding to these intracluster motions can be computed accurately using a simple theory.

Monodisperse patchy particle glass former

Glass formers are characterized by their ability to avoid crystallization. As monodisperse systems tend to rapidly crystallize, the most common glass formers in simulations are systems composed of mixtures of particles with different sizes. Here, we make use of the ability of patchy particles to change their local structure to propose them as monodisperse glass formers. We explore monodisperse systems with two patch geometries: a 12-patch geometry that enhances the formation of icosahedral clusters and an 8-patch geometry that does not appear to strongly favor any particular local structure. We show that both geometries avoid crystallization and present glassy features at low temperatures. However, the 8-patch geometry better preserves the structure of a simple liquid at a wide range of temperatures and packing fractions, making it a good candidate for a monodisperse glass former.

Self-diffusion in bidisperse systems of magnetic nanoparticles

In the present paper, we study the self-diffusion of aggregating magnetic particles in bidisperse ferrofluids. We employ density functional theory (DFT) and coarse-grained molecular dynamics (MD) simulations to investigate the impact of granulometric composition of the system on the cluster self-diffusion. We find that the presence of small particles leads to the overall increase of the self-diffusion rate of clusters due the change in cluster size and composition.

Field-induced anti-nematic and biaxial ordering in binary mixtures of discotic mesogens and spherical magnetic nanoparticles

Using computer simulations we explore the equilibrium structure and response to external stimuli of complex magnetic hybrids consisting of magnetic particles in discotic liquid crystalline matrices. We show that the anisotropy of the liquid crystalline matrix (either in the nematic or in the columnar phase) promotes the collective orientational ordering of self-assembled magnetic particles. Upon applying an external homogeneous magnetic field in an otherwise isotropic state, the magnetic particles self-assemble into linear-rodlike-chains. At the same time structural changes occur in the matrix. The matrix transforms from an isotropic to a non-conventional anti-nematic state in which the symmetry axis of the discs is, on average, perpendicular to the magnetic field. In addition, a stable biaxial nematic state is found upon applying an external field to an otherwise uniaxial discotic nematic state. These observed morphologies constitute an appealing alternative to binary mixtures of rigid rod-disc system and indicate that non-trivial biaxial ordering can be obtained in the presence of a uniaxial external stimulus.

Chain Formation and Phase Separation in Ferrofluids: The Influence on Viscous Properties

Ferrofluids have attracted considerable interest from researchers and engineers due to their rich set of unique physical properties that are valuable for many industrial and biomedical applications. Many phenomena and features of ferrofluids' behavior are determined by internal structural transformations in the ensembles of particles, which occur due to the magnetic interaction between the particles. An applied magnetic field induces formations, such as linear chains and bulk columns, that become elongated along the field. In turn, these structures dramatically change the rheological and other physical properties of these fluids. A deep and clear understanding of the main features and laws of the transformations is necessary for the understanding and explanation of the macroscopic properties and behavior of ferrofluids. In this paper, we present an overview of experimental and theoretical works on the internal transformations in these systems, as well as on the effect of the internal structures on the rheological effects in the fluids.

Perturbation approaches for describing dipolar fluids and electrolyte solutions

This work proposes perturbation approaches for describing dipolar fluids as well as model and aqueous electrolyte solutions. The electrostatic pair potentials are split into short- and long-ranged contributions, whereas a third order perturbation expansion is applied for the short-ranged potentials. This circumvents the problem of divergent correlation integrals. The dipolar perturbation terms are represented through a [2,1]-Padé approximation to resum the poorly convergent series. For the remaining charge-charge and charge-dipole contributions, we present a new approximant, which provides a (quasi)linear dependence of the Helmholtz energy. The underlying correlation integrals are adjusted to results from molecular simulations. The long-ranged contribution to the electrostatic interactions is treated through an analytic expression developed by Rodgers and Weeks [J. Chem. Phys. 131, 244108 (2010)]. Theoretical predictions of our perturbation theory are compared to results from a widely used integral equation theory, namely, the mean spherical approximation, and we find that our perturbation theory provides much more accurate results. Furthermore, the theory shows some quantities in rather good agreement with reference data, namely, Helmholtz energies, internal energies, and densities at higher densities of solutions. Limitations of the approach, however, are observed for several other partial molar quantities, such as the mean activity coefficient.

Dynamical self-assembly of dipolar active Brownian particles in two dimensions

Based on Brownian Dynamics (BD) simulations, we study the dynamical self-assembly of active Brownian particles with dipole-dipole interactions, stemming from a permanent point dipole at the particle center. The propulsion direction of each particle is chosen to be parallel to its dipole moment. We explore a wide range of motilities and dipolar coupling strengths and characterize the corresponding behavior based on several order parameters. At low densities and low motilities, the most important structural phenomenon is the aggregation of the dipolar particles into chains. Upon increasing the particle motility, these chain-like structures break, and the system transforms into a weakly correlated isotropic fluid. At high densities, we observe that the motility-induced phase separation is strongly suppressed by the dipolar coupling. Once the dipolar coupling dominates the thermal energy, the phase separation disappears, and the system rather displays a flocking state, where particles form giant clusters and move collective along one direction. We provide arguments for the emergence of the flocking behavior, which is absent in the passive dipolar system.

Self-assembly of Pseudo-Dipolar Nanoparticles at Low Densities and Strong Coupling

Nanocolloids having directional interactions are highly relevant for designing new self-assembled materials easy to control. In this article we report stochastic dynamics simulations of finite-size pseudo-dipolar colloids immersed in an implicit dielectric solvent using a realistic continuous description of the quasi-hard Coulombic interaction. We investigate structural and dynamical properties near the low-temperature and highly-diluted limits. This system self-assembles in a rich variety of string-like configurations, depicting three clearly distinguishable regimes with decreasing temperature: fluid, composed by isolated colloids; string-fluid, a gas of short string-like clusters; and string-gel, a percolated network. By structural characterization using radial distribution functions and cluster properties, we calculate the state diagram, verifying the presence of string-fluid regime. Regarding the string-gel regime, we show that the antiparallel alignment of the network chains arises as a novel self-assembly mechanism when the characteristic interaction energy exceeds the thermal energy in two orders of magnitude, u_d/k_BT ≈ 100. This is associated to relevant structural modifications in the network connectivity and porosity. Furthermore, our results give insights about the dynamically-arrested nature of the string-gel regime, where we show that the slow relaxation takes place in minuscule energy steps that reflect local rearrangements of the network.

Arrested dynamics of the dipolar hard sphere model

We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (soft core) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down of the dynamics of the translational and orientational degrees of freedom: at low (η = 0.2) and intermediate (η = 0.4) volume fractions, both dynamics are strongly coupled and become simultaneously arrested upon cooling. At high concentrations (η≥ 0.6), the translational dynamics shows the features of an ordinary glass transition, either by compressing or cooling down the system, but with the orientations remaining ergodic, thus indicating the existence of partially arrested states. In this density regime, but at lower temperatures, the relaxation of the orientational dynamics also freezes. The physical scenario provided by the simulations is discussed and compared against results obtained with the self-consistent generalized Langevin equation theory, and both provide a consistent description of the dynamical arrest transitions in the system. Our results are summarized in an arrested states diagram which qualitatively organizes the simulation data and provides a generic picture of the glass transitions of a dipolar fluid.

The structure of clusters formed by Stockmayer supracolloidal magnetic polymers

Unlike Stockmayer fluids, that prove to undergo gas-liquid transition on cooling, the system of dipolar hard or soft spheres without any additional central attraction so far has not been shown to have a critical point. Instead, in the latter, one observes diverse self-assembly scenarios. Crosslinking dipolar soft spheres into supracolloidal magnetic polymer-like structures (SMPs) changes the self-assembly behaviour. Moreover, aggregation in systems of SMPs strongly depends on the constituent topology. For Y- and X-shaped SMPs, under the same conditions in which dipolar hard spheres would form chains, the formation of very large loose gel-like clusters was observed (E. Novak et al., J. Mol. Liq. 271, 631 (2018)). In this work, using molecular dynamics simulations, we investigate the self-assembly in suspensions of four topologically different SMPs --chains, rings, X and Y-- whose monomers interact via Stockmayer potential. As expected, compact drop-like clusters are formed by SMPs in all cases if the central isotropic attraction is introduced, however, their shape and internal structure turn out to depend on the SMPs topology.

Field-responsive colloidal assemblies defined by magnetic anisotropy

Particle dispersions provide a promising tool for the engineering of functional materials that exploit self-assembly of complex structures. Dispersion made from magnetic colloidal particles is a great choice; they are biocompatible and remotely controllable among many other advantages. However, their dominating dipolar interaction typically limits structural complexity to linear arrangements. This paper shows how a magnetostatic equilibrium state with noncollinear arrangement of the magnetic moments, as reported for ferromagnetic Janus particles, enables the controlled self-organization of diverse structures in two dimensions via constant and low-frequency external magnetic fields. Branched clusters of staggered chains, compact clusters, linear chains, and dispersed single particles can be formed and interconverted reversibly in a controlled way. The structural diversity is a consequence of both the inhomogeneity and the spatial extension of the magnetization distribution inside the particles. We draw this conclusion from calculations based on a model of spheres with multiple shifted dipoles. The results demonstrate that fundamentally new possibilities for responsive magnetic materials can arise from interactions between particles with a spatially extended, anisotropic magnetization distribution.

Self-assembly of magnetic colloids with shifted dipoles

The self-assembly of colloidal magnetic Janus particles with a laterally displaced (or shifted), permanent dipole in a quasi-two-dimensional system is studied using Brownian dynamics simulations. The rate of formation of clusters and their structures are quantified for several values of dipolar shift from the particle center, which is nondimensionalized using the particle's radius so that it takes values ranging from 0 to 1, and examined under different magnetic interaction strengths relative to Brownian motion. For dipolar shifts close to 0, chain-like structures are formed, which grow at long times following a power law, while particles of shift higher than 0.2 generally aggregate in ring-like clusters that experience limited growth. In the case of shifts between 0.4 and 0.5, the particles tend to aggregate in clusters of 3 to 6, while for all shifts higher than 0.6 clusters rarely contain more than 3 particles due to the antiparallel dipole orientations that are most stable at those shifts. The strength of the magnetic interactions hastens the rate at which clusters are formed; however, the effect it has on cluster size is lessened by increases in the shift of the dipoles. These results contribute to better understand the dynamics of magnetic Janus particles and can help the synthesis of functionalized materials for specific applications such as drug delivery.

Speeding up the study of diluted dipolar systems

We study the regimes of a diluted dipolar system through Monte Carlo numerical simulations in the NVT ensemble. To accelerate the dynamics, several approximations and speed-up algorithms are proposed and tested. In particular, it turns out that "cluster move Monte Carlo" algorithm speeds-up to two decades faster than the traditional Monte Carlo, depending on temperature and density. We find simple-fluid, chain-fluid, ring-fluid, gel, and antiparallel columnar regimes, which are studied and characterized through positional, orientational, and thermodynamical observables.

Phase diagrams of mixtures of dipolar rods and discs

Self-assembly of binary mixtures that contain anisotropic, interacting colloidal particles have been proposed as a way to create new, multi-functional materials. We simulate binary mixtures of dipolar rods and dipolar discs in two-dimensions using discontinuous molecular dynamics to determine how the assembled structures of these mixtures differ from those seen in single component systems. Two different binary mixtures are investigated: a mixture of an equal number of dipolar rods and dipolar discs ("equal number"), and a mixture where the area fraction of dipolar rods is equal to the area fraction of dipolar discs ("equal area"). Phase boundaries between fluid, string-fluid, and "gel" phases are calculated and compared to the phase boundaries of the pure components. Looking deeper at the underlying structure of the mixture reveals a complex interplay between the rods and discs and the formation of states where the two components are in different phases. The mixtures exhibit phases where both rods and discs are in the fluid phase, where rods form a string-fluid while discs remain in the fluid phase, a rod string-fluid coexisting with a disc string-fluid, a "gel" that consists primarily of rods while the discs form either a fluid or string-fluid phase, and a "gel" that contains both rods and discs. Our results give insight into the general assembly pathway of binary mixtures, and how complex aggregates can be created by varying the mixture composition, strength of interaction between the two components, and the temperature. By manipulating the properties of one of the components it should be possible to fabricate bifunctional, thermally responsive self-assembled materials.

How to simulate patchy particles⋆

Patchy particles is the name given to a large class of systems of mesoscopic particles characterized by a repulsive core and a discrete number of short-range and highly directional interaction sites. Numerical simulations have contributed significantly to our understanding of the behaviour of patchy particles, but, although simple in principle, advanced simulation techniques are often required to sample the low temperatures and long time-scales associated with their self-assembly behaviour. In this work we review the most popular simulation techniques that have been used to study patchy particles, with a special focus on Monte Carlo methods. We cover many of the tools required to simulate patchy systems, from interaction potentials to biased moves, cluster moves, and free-energy methods. The review is complemented by an educationally oriented Monte Carlo computer code that implements all the techniques described in the text to simulate a well-known tetrahedral patchy particle model.

Temperature- and field-induced structural transitions in magnetic colloidal clusters

Magnetic colloidal clusters can form chain, ring, and more compact structures depending on their size. In the present investigation we examine the combined effects of temperature and external magnetic field on these configurations by means of extensive Monte Carlo simulations and a dedicated analysis based on inherent structures. Various thermodynamical, geometric, and magnetic properties are calculated and altogether provide evidence for possibly multiple structural transitions at low external magnetic field. Temperature effects are found to overcome the ordering effect of the external field, the melted stated being associated with low magnetization and a greater compactness. Tentative phase diagrams are proposed for selected sizes.

Free energy calculations for rings and chains formed by dipolar hard spheres

We employ a method based on Monte Carlo grand-canonical simulations to precisely calculate partition functions of non-interacting chains and rings formed by dipolar hard spheres (DHS) at low temperature. The extended low temperature region offered by such cluster calculations, compared to what had been previously achieved with standard simulations, opens up the possibility of exploring a part of the DHS phase diagram which was inaccessible before. The reported results offer the unique opportunity of verifying well-established theoretical models based on the ideal gas of cluster approximation in order to clarify their range of validity. They also provide the basis for future studies in which cluster-cluster interactions will be included.

Modified mean-field theory of the magnetic properties of concentrated, high-susceptibility, polydisperse ferrofluids

The effects of particle-size polydispersity on the magnetostatic properties of concentrated ferrofluids are studied using theory and computer simulation. The second-order modified mean-field (MMF2) theory of Ivanov and Kuznetsova [Phys. Rev. E 64, 041405 (2001)1063-651X10.1103/PhysRevE.64.041405] has been extended by calculating additional terms of higher order in the dipolar coupling constant in the expansions of the initial magnetic susceptibility and the magnetization curve. The theoretical predictions have been tested rigorously against results from Monte Carlo simulations of model monodisperse, bidisperse, and highly polydisperse ferrofluids. Comparisons have been made between systems with the same Langevin susceptibility and the same saturation magnetization. In all cases, the new theoretical magnetization curve shows better agreement with simulation data than does the MMF2 theory. As for the initial susceptibility, MMF2 theory is most accurate for the monodisperse model, while the new theory works best for polydisperse systems with a significant proportion of large particles. These results are important for the analysis and characterization of recently synthesized polydisperse ferrofluids with record-breaking values of the initial magnetic susceptibility.

Reversible and irreversible aggregation of magnetic liposomes

Understanding stabilization and aggregation in magnetic nanoparticle systems is crucial to optimizing the functionality of these systems in real physiological applications. Here we address this problem for a specific, yet representative, system. We present an experimental and analytical study on the aggregation of superparamagnetic liposomes in suspension in the presence of a controllable external magnetic field. We study the aggregation kinetics and report an intermediate time power law evolution and a long time stationary value for the average aggregate diffusion coefficient, both depending on the magnetic field intensity. We then show that the long time aggregate structure is fractal with a fractal dimension that decreases upon increasing the magnetic field intensity. By scaling arguments we also establish an analytical relation between the aggregate fractal dimension and the power law exponent controlling the aggregation kinetics. This relation is indeed independent on the magnetic field intensity. Despite the superparamagnetic character of our particles, we further prove the existence of a population of surviving aggregates able to maintain their integrity after switching off the external magnetic field. Finally, we suggest a schematic interaction scenario to rationalize the observed coexistence between reversible and irreversible aggregation.

Self-assembly behaviour of hetero-nuclear Janus dumbbells

We investigate the fluid structure and self-assembly of a system of Janus dumbbells by means of aggregation-volume-bias Monte Carlo simulations and Simulated Annealing techniques. In our approach, Janus dumbbells model asymmetric colloidal particles constituted by two tangent (touching) spheres (labelled as h and s) of different sizes and interaction properties: specifically, the h spheres interact with all other spheres belonging to different dumbbells via hard-sphere potentials, whereas two s spheres interact via a square-well potential. By introducing a parameter α ∈ [0,2] that controls the size ratio between the h and s spheres, we are able to investigate the overall phase behaviour of Janus dumbbells as a function of α. In a previous paper (O'Toole et al., Soft Matter, 2017, 13, 803) we focused on the region where the s sphere is larger than the h sphere (α > 1), documenting the presence of a variety of phase behaviours. Here we investigate a different regime of size ratios, predominantly where the hard sphere is larger than (or comparable to) the attractive one. Under these conditions, we observe the onset of many different self-assembled super-structures. Depending on the specific value of α we document the presence of spherical clusters (micelles) progressively evolving into more exotic structures including platelets, filaments, networks and percolating fluids, sponge structures and lamellar phases. We find no evidence of a gas-liquid phase separation for α ≤ 1.1, since under these conditions it is pre-empted by the development of self-assembled phases.

How cube-like must magnetic nanoparticles be to modify their self-assembly?

Systems whose magnetic response can be finely tuned using control parameters, such as temperature and external magnetic field strength, are extremely desirable, functional materials. Magnetic nanoparticles, in particular suspensions thereof, offer opportunities for this controllability to be realised. Cube-like particles are particularly mono-disperse examples that, together with their favourable packing behaviour, make them of significant interest for study. Using a combination of analytical calculations and molecular dynamics we have studied the self-assembly of permanently magnetised dipolar superballs. The superball shape parameter was varied in order to interpolate the region between the already well-studied sphere system and that of the recently investigated cube. Our findings show that as a superball particle becomes more cubic the chain to ring transition, observed in the ground state of spherical particles, occurs at an increasingly larger cluster size. This effect is mitigated, however, by the appearance of a competing configuration; asymmetric rings, a conformation that we show superballs can readily adopt.

Simulation study on the structural properties of colloidal particles with offset dipoles

A major research theme in materials science is determining how the self-assembly of new generations of colloidal particles of complex shape and surface charge is guided by their interparticle interactions. In this paper, we describe results from quasi-2D Monte Carlo simulations of systems of colloidal particles with offset transversely-oriented extended dipole-like charge distributions interacting via an intermediate-ranged Yukawa potential. The systems are cooled slowly through an annealing procedure during which the temperature is lowered in discrete steps, allowing the system to equilibrate. We perform ground state calculations for two, three and four particles at several shifts of the dipole vector from the particle center. We create state diagrams in the plane spanned by the temperature and the area fraction outlining the boundaries between fluid, string-fluid and percolated states at various values of the shift. Remarkably we find that the effective cooling rate in our simulations has an impact on the structures formed, with chains being more prevalent if the system is cooled quickly and cyclic structures more prevalent if the system is cooled slowly. As the dipole is further shifted from the center, there is an increased tendency to assemble into small cyclic structures at intermediate temperatures. These systems further self-assemble into open lattice-like arrangements at very low temperatures. The novel structures identified might be useful for photonic applications, new types of porous media for filtration and catalysis, and gel matrices with unusual properties.

Sedimentation equilibrium of magnetic nanoparticles with strong dipole-dipole interactions

Langevin dynamics simulation is used to study the suspension of interacting magnetic nanoparticles (dipolar spheres) in a zero applied magnetic field and in the presence of a gravitational (centrifugal) field. A particular emphasis is placed on the equilibrium vertical distribution of particles in the infinite horizontal slab. An increase in the dipolar coupling constant λ (the ratio of dipole-dipole interaction energy to thermal energy) from zero to seven units causes an increase in the particle segregation coefficient by several orders of magnitude. The effect of anisotropic dipole-dipole interactions on the concentration profile of particles is the same as that of the isotropic van der Waals attraction modeled by the Lennard-Jones potential. In both cases, the area with a high-density gradient separating the area with high and low particle concentration is formed on the profiles. Qualitative difference between two potentials manifests itself only in the fact that in the absence of a gravitational field the dipole-dipole interactions do not lead to the "gas-liquid" phase transition: no separation of the system into weakly and highly concentrated phases is observed. At high particle concentration and at large values of λ, the orientational ordering of magnetic dipoles takes place in the system. Magnetic structure of the system strongly depends on the imposed boundary conditions. Spontaneous magnetization occurs in the infinite horizontal slab (i.e., in the rectangular cell with two-dimensional periodic boundary conditions). Replacement of the infinite slab by the finite-size hard-wall vertical cylinder leads to the formation of azimuthal (vortex-like) order. The critical values of the coupling constant corresponding to the transition into an ordered state are very close for two geometries.

Reconfiguring active particles by electrostatic imbalance

Active materials represent a new class of condensed matter in which motile elements may collectively form dynamic, global structures out of equilibrium. Here, we present a general strategy to reconfigure active particles into various collective states by introducing imbalanced interactions. We demonstrate the concept with computer simulations of self-propelled colloidal spheres, and experimentally validate it in a two-dimensional (2D) system of metal-dielectric Janus colloids subjected to perpendicular a.c. electric fields. The mismatched, frequency-dependent dielectric responses of the two hemispheres of the colloids allow simultaneous control of particle motility and colloidal interactions. We realized swarms, chains, clusters and isotropic gases from the same precursor particle by changing the electric-field frequency. Large-scale polar waves, vortices and jammed domains are also observed, with the persistent time-dependent evolution of their collective structure evoking that of classical materials. This strategy of asymmetry-driven active self-organization should generalize rationally to other active 2D and three-dimensional (3D) materials.

Orientational order and translational dynamics of magnetic particle assemblies in liquid crystals

Implementing extensive molecular dynamics simulations we explore the organization of magnetic particle assemblies (clusters) in a uniaxial liquid crystalline matrix comprised of rodlike particles. The magnetic particles are modelled as soft dipolar spheres with diameter significantly smaller than the width of the rods. Depending on the dipolar strength coupling the magnetic particles arrange into head-to-tail configurations forming various types of clusters including rings (closed loops) and chains. In turn, the liquid crystalline matrix induces long range orientational ordering to these structures and promotes their diffusion along the director of the phase. Different translational dynamics are exhibited as the liquid crystalline matrix transforms either from isotropic to nematic or from nematic to smectic state. This is caused due to different collective motion of the magnetic particles into various clusters in the anisotropic environments. Our results offer a physical insight for understanding both the structure and dynamics of magnetic particle assemblies in liquid crystalline matrices.

Self-assembly of three-dimensional ensembles of magnetic particles with laterally shifted dipoles

We consider a model of colloidal spherical particles carrying a permanent dipole moment which is laterally shifted out of the particles' geometrical centres, i.e. the dipole vector is oriented perpendicular to the radius of the particles. Varying the shift δ from the centre, we analyse ground state structures for two, three and four hard spheres, using a simulated annealing procedure. We also compare earlier ground state results. We then consider a bulk system at finite temperatures and different densities. Using molecular dynamics simulations, we examine the equilibrium self-assembly properties for several shifts. Our results show that the shift of the dipole moment has a crucial impact on both the ground state configurations as well as the self-assembled structures at finite temperatures.

Self-assembly of colloidal magnetic particles: energy landscapes and structural transitions

The self-assembly of colloidal magnetic particles is of particular interest for the rich variety of structures it produces and the potential for these systems to be reconfigurable.

Temperature-dependent dynamic correlations in suspensions of magnetic nanoparticles in a broad range of concentrations: a combined experimental and theoretical study

We study the effects of temperature and concentration on the dynamic spectra of polydisperse magnetic nanoparticle suspensions.

Phase separation and self-assembly in a fluid of Mickey Mouse particles

Recent developments in the synthesis of colloidal particles allow for control over shape and inter-particle interaction. One example, among others, is the so-called "Mickey Mouse" (MM) particle for which the self-assembly properties have been previously studied yielding a stable cluster phase together with elongated, tube-like structures. Here, we investigate under which conditions a fluid of Mickey Mouse particles can yield phase separation and how the self-assembly behaviour affects the gas-liquid coexistence. We vary the distance between the repulsive and the attractive lobes (bond length), and the interaction range, and follow the evolution of the gas-liquid (GL) coexistence curve. We find that upon increasing the bond length distance the binodal line shifts to lower temperatures, and that the interaction range controls the transition between phase separation and self-assembly of clusters. Upon further reduction of the interaction range and temperature, the clusters assume an increasingly ordered tube-like shape, ultimately matching the one previously reported in literature. These results are of interest when designing particle shape and particle-particle interaction for self-assembly processes.

Self-organization of magnetic moments in dipolar chains with restricted degrees of freedom

Equilibrium behavior of a single chain of dipolar spheres is investigated by the method of molecular dynamics in a wide range of the dipolar coupling constant λ. Two cases are considered: rodlike and flexible chains. In the first case, particle centers are immovably fixed on one axis, but their magnetic moments retain absolute orientational freedom. It has been found that at λ≳1.5 particle moments are chiefly aligned parallel to the chain axis, but the total moment of the chain continuously changes its sign with some mean frequency, which exponentially decreases with the growth of λ. Such behavior of the rodlike chain is analogous to the Néel relaxation of a superparamagnetic particle with a finite energy of magnetic anisotropy. In the flexible chain particles are able to move in the three-dimensional space, but the distance between centers of the first-nearest neighbors never exceeds a given limiting value r(max). If r(max)≃d (d is the particle diameter) then the most probable shape of the chain of five or more particles at λ≳6 is that of a ring. The behavior of chains with r(max)≥2d is qualitatively different: At λ≃4 long chains collapse into dense quasispherical globules and at λ≳8 these globules take toroidal configuration with a spontaneous azimuthal ordering of magnetic dipoles. With the increase of r(max) to larger values (r(max)>10d) globules expand and break down to form separate rings.

Tunable structures of mixtures of magnetic particles in liquid-crystalline matrices

We investigate the self-organization of a binary mixture of similar sized rods and dipolar soft spheres by means of Monte-Carlo simulations. We model interparticle interactions by employing anisotropic Gay-Berne, dipolar and soft-sphere interactions. In the limit of vanishing magnetic moments we obtain a variety of fully miscible liquid crystalline phases including nematic, smectic and lamellar phases. For the magnetic mixture, we find that the liquid crystalline matrix supports the formation of orientationally ordered ferromagnetic chains. Depending on the relative size of the species the chains align parallel or perpendicular to the director of the rods forming uniaxial or biaxial nematic, smectic and lamellar phases. As an exemplary external perturbation we apply a homogeneous magnetic field causing uniaxial or biaxial ordering to an otherwise isotropic state.

Behaviour of magnetic Janus-like colloids

We present a theoretical study of Janus-like magnetic particles at low temperature. To describe the basic features of the Janus-type magnetic colloids, we put forward a simple model of a spherical particle with a dipole moment shifted outwards from the centre and oriented perpendicular to the particle radius. Using direct calculations and molecular dynamics computer simulations, we investigate the ground states of small clusters and the behaviour of bigger systems at low temperature. In both cases the important parameter is the dipolar shift, which leads to different ground states and, as a consequence, to a different microscopic behaviour in the situation when the thermal fluctuations are finite. We show that the head-to-tail orientation of dipoles provides a two-particle energy minima only if the dipoles are not shifted from the particle centres. This is one of the key differences from the system of shifted dipolar particles (sd-particles), in which the dipole was shifted outwards radially, studied earlier (Kantorovich et al 2011 Soft Matter 7 5217-27). For sd-particles the dipole could be shifted out of the centre for almost 40% before the head-to-tail orientation was losing its energetic advantage. This peculiarity manifests itself in the topology of the small clusters in the ground state and in the response of the Janus-like particle systems to an external magnetic field at finite temperatures.

Hierarchical self-assembly of colloidal magnetic particles into reconfigurable spherical structures

Colloidal self-assembly has enormous potential as a bottom-up means of structure fabrication. Here we demonstrate hierarchical self-assembly of rationally designed charge-stabilised colloidal magnetic particles into ground state structures that are topologically equivalent to a snub cube and a snub dodecahedron, the only two chiral Archimedean solids, for size-selected clusters. These spherical structures open up in response to an external magnetic field and demonstrate controllable porosity. Such features are critical to their applications as functional materials.

Temperature-induced structural transitions in self-assembling magnetic nanocolloids

With the help of a unique combination of density functional theory and computer simulations, we discover two possible scenarios, depending on concentration, for the hierarchical self-assembly of magnetic nanoparticles on cooling.

Structure factor of model bidisperse ferrofluids with relatively weak interparticle interactions

In the present manuscript we develop a theoretical approach to describe the pair correlation function of bidisperse magnetic dipolar hard- and soft-spheres. We choose bidisperse system as the first step to allow for polydispersity when studying thermodynamics of magnetic fluids. Using diagram technique we calculate the virial expansion of the pair correlation function up to the first order in density and fourth order in the dipolar strength. Even though, the radial distribution functions are extremely sensitive to the steric potential, we show that the behaviour of the isotropic centre-centre structure factor is almost indifferent to the type of the short-range repulsion. We extensively compare our theoretical results to the data of molecular dynamics simulations, which helps us to understand the range of validity of the virial expansion both on density and magnetic dipolar strength. We also investigate the influence of the granulometric composition on the height, width, and position of the structure factor first peak in order to clarify whether it is possible to extract structural information from experimentally measured small angle neutron scattering intensities.