Two-Steps Versus One-Step Solidification Pathways of Binary Metallic Nanodroplets

The solidification of AgCo, AgNi, and AgCu nanodroplets is studied by molecular dynamics simulations in the size range of 2-8 nm. All these systems tend to phase separate in the bulk solid with surface segregation of Ag. Despite these similarities, the simulations reveal clear differences in the solidification pathways. AgCo and AgNi already separate in the liquid phase, and they solidify in configurations close to equilibrium. They can show a two-step solidification process in which Co-/Ni-rich parts solidify at higher temperatures than the Ag-rich part. AgCu does not separate in the liquid and solidifies in one step, thereby remaining in a kinetically trapped state down to room temperature. The solidification mechanisms and the size dependence of the solidification temperatures are analyzed, finding qualitatively different behaviors in AgCo/AgNi compared to AgCu. These differences are rationalized by an analytical model.

Self-Organization of Large Alumina Platelets and Silica Nanoparticles by Heteroaggregation and Sedimentation: Toward an Alternative Shaping of Nacre-Like Ceramic Composites

Nacre-like ceramic composites are of importance in a wide range of applications, because of their mechanical properties, combining high mechanical strength and high fracture toughness. Those mechanical properties are the result of strongly aligned platelets glued in a matrix. Different methods exist to shape such a "brick-and-mortar" hierarchical structure. In this paper, we propose to use the phenomenon of heteroaggregation between silica nanoparticles and large alumina platelets. Experimental and numerical results show that silica nanoparticles can adsorb on alumina platelets with good distribution. This adsorption promotes the deagglomeration of alumina that can self-organize in layers by sedimentation. This phenomenon can be exploited to shape alumina-silica nacre-like composites.

Brownian dynamics simulations of one-patch inverse patchy particles

Inverse patchy particles are promising colloids to develop new architectures in ceramic materials based on their self-assembly. Nonetheless, a good understanding of their aggregation is required. Several previous studies have shown that the behavior of ceramic colloids can be well described by the DLVO interaction potential. In the present paper, we develop new coarse-grained Brownian dynamics simulations, where particles are represented by an assembly of beads interacting via DLVO interactions, whose parameters can be directly linked to experimental characterization. First, the validity of the simulations is proved by studying the heteroaggregation of homogeneously charged particles. Then, simulations are applied to one-patch inverse patchy particles to study the effect of the patch size. They show that the smaller the patch, the more elongated the aggregates. Simulations are also performed to understand the role of the Debye screening length in the particular case of large patches and they show that aggregation leads always to compact aggregates.

Role of Electrostatic Interactions in Oil-in-Water Emulsions Stabilized by Heteroaggregation: An Experimental and Simulation Study

Oil-in-water emulsion stabilization by heteroaggregation of hydrophilic particles without a surfactant is of importance in a wide range of applications; however, the stabilization mechanism is little described. To shed light on the early stage of the stabilization mechanism, a model system composed of an oil wax phase dispersed in water with oppositely charged colloidal particles is studied experimentally and numerically. Experiments show that the colloids do not penetrate deeply in the oil phase, suggesting that adsorption of the colloidal particles on the wax droplets is mainly due to electrostatic interactions. Experiments and Brownian dynamics simulations show also that when oppositely charged colloidal particles are present in the emulsion, a multilayer coating of heteroaggregated colloidal particles is formed around the wax droplets. This protective coating is expected to prevent from the oil droplet coalescence and therefore to stabilize the emulsion.

Computer simulations of heteroaggregation with large size asymmetric colloids

HYPOTHESIS

Hetero-aggregation of inorganic colloids is influenced by numerous parameters, which dictate the suspension properties. When particles are different in size, the suspension can be either stable or unstable according to concentration of components, ionic strength, and pH. Experimentally, understanding the role of each parameter is sometimes difficult because parameters cannot easily be modified independently. Numerical simulations are thus very useful to discriminate between different effects.

SIMULATIONS

Brownian dynamics simulations are used here to study the heteroaggregation of dilute suspensions composed of two populations of colloids with large size asymmetry. Special attention is paid to the effect of small-particle concentration, surface potentials, and ionic strength.

FINDINGS

The simulation results show that hetero-aggregation can be tuned by modifying these different parameters, and that the resulting aggregate structures depend more on the surface properties of small particles than on those of large particles. The simulations shed light on a further parameter crucially influencing hetero-aggregation, i.e. the mobility of small particles when adsorbed on large ones. The present results rationalize numerous experimental observations reported in the literature and can be used as reference to explain future experimental observations.

Early Dynamics and Stabilization Mechanisms of Oil-in-Water Emulsions Containing Colloidal Particles Modified with Short Amphiphiles: A Numerical Study

Emulsions stabilized by mixtures of particles and amphiphilic molecules are relevant for a wide range of applications, but their dynamics and stabilization mechanisms on the colloidal level are poorly understood. Given the challenges to experimentally probe the early dynamics and mechanisms of droplet stabilization, Brownian dynamics simulations are developed here to study the behavior of oil-in-water emulsions stabilized by colloidal particles modified with short amphiphiles. Simulation parameters are based on an experimental system that consists of emulsions obtained with octane as the oil phase and a suspension of alumina colloidal particles modified with short carboxylic acids as the continuous aqueous medium. The numerical results show that attractive forces between the colloidal particles favor the formation of closely packed clusters on the droplet surface or of a percolating network of particles throughout the continuous phase, depending on the amphiphile concentration. Simulations also reveal the importance of a strong adsorption of particles at the liquid interface to prevent their depletion from the droplet surface when another droplet approaches. Strongly adsorbed particles remain immobile on the droplet surface, generating an effective steric barrier against droplet coalescence. These findings provide new insights into the early dynamics and mechanisms of stabilization of emulsions using particles and amphiphilic molecules.

Shear viscosity in hard-sphere and adhesive colloidal suspensions with reverse non-equilibrium molecular dynamics

We employ the reverse non-equilibrium molecular dynamics method (RNEMD) of Müller-Plathe [Phys. Rev. E, 1999, 59, 4894] to calculate the shear viscosity of colloidal suspensions within the stochastic rotation dynamics-molecular dynamics (SRD-MD) simulation method. We examine the influence of different coupling schemes in SRD-MD on the colloidal volume fraction ϕ_c dependent viscosity from the dilute limit up to ϕ_c = 0.3. Our results demonstrate that the RNEMD method is a robust and reliable method for calculating rheological properties of colloidal suspensions. To obtain quantitatively accurate results beyond the dilute regime, the hydrodynamic interactions between the effective fluid particles in the SRD and the MD colloidal particles must be carefully considered in the coupling scheme. We benchmark the method by comparing with the hard sphere suspension case, and then calculate relative viscosities for colloids with mutually attractive interactions. We show that the viscosity displays a sharp increase at the onset of aggregation of the colloidal particles with increasing volume fraction and attraction.

Interdiffusion and crystallization of oppositely charged colloids

Innovative way of controlling colloidal heteroaggregation.

Aggregation of binary colloidal suspensions on attractive walls

The adsorption of colloidal particles from a suspension on a solid surface is of fundamental importance to many physical and biological systems. In this work, Brownian Dynamics simulations are performed to study the aggregation in a suspension of oppositely charged colloidal particles in the presence of an attractive wall. For sufficiently strong attractions, the wall alters the microstructure of the aggregates so that B2 (CsCl-type) structures are more likely obtained instead of B1 (NaCl-type) structures. The probability of forming either B1 or B2 crystallites depends also on the inverse interaction range κa. Suspensions with small κa are more likely to form B2 crystals than suspensions with larger κa, even if the energetic stability of the B2 phase decreases with decreasing κa. The mechanisms underlying this aggregation and crystallization behaviour are analyzed in detail.

Computation of shear viscosity of colloidal suspensions by SRD-MD

The behaviour of sheared colloidal suspensions with full hydrodynamic interactions (HIs) is numerically studied. To this end, we use the hybrid stochastic rotation dynamics-molecular dynamics (SRD-MD) method. The shear viscosity of colloidal suspensions is computed for different volume fractions, both for dilute and concentrated cases. We verify that HIs help in the collisions and the streaming of colloidal particles, thereby increasing the overall shear viscosity of the suspension. Our results show a good agreement with known experimental, theoretical, and numerical studies. This work demonstrates the ability of SRD-MD to successfully simulate transport coefficients that require correct modelling of HIs.

Compact and ordered colloidal clusters from assembly-disassembly cycles: a numerical study

Brownian dynamics simulations are used to investigate the assembly of attractive colloids whose interaction potential well is periodically changed over time. Our system is composed of spherical, mono-disperse, highly charged, alumina particles whose interactions are modeled by the DLVO theory. The depth of the potential well is periodically changed by varying the ionic strength of the liquid medium. The simulations show that, with a right choice of some key parameters, a potential well depth alternating between low and higher values allows a faster aggregation into more compact and also more ordered colloidal clusters. This result is quantified by the computation of two relevant coordination parameters during the aggregation. This finding may help elucidate the assembly of colloidal particles in complex biological processes (e.g. biomineralization) and could be useful for the development of photonic crystals from attractive colloidal particles.

Numerical and experimental study of suspensions containing carbon blacks used as conductive additives in composite electrodes for lithium batteries

Suspensions of carbon blacks and spherical carbon particles are studied experimentally and numerically to understand the role of the particle shape on the tendency to percolation. Two commercial carbon blacks and one lab-synthesized spherical carbon are used. The percolation thresholds in suspensions are experimentally determined by two complementary methods: impedance spectroscopy and rheology. Brownian dynamics simulations are performed to explain the experimental results taking into account the fractal shape of the aggregates in the carbon blacks. The results of Brownian dynamics simulations are in good agreement with the experimental results and allow one to explain the experimental behavior of suspensions.

Study of the B1-B2 transition in colloidal clusters

The possible mechanisms for the B1 (NaCl-type) to B2 (CsCl-type) transition in crystalline colloidal clusters of equally sized particles are studied by means of two computational techniques: metadynamics and nudged elastic band calculations. The system is modelled by a screened Coulomb potential. Different interaction ranges are considered. The transition from a perfect NaCl cubic cluster to a full CsCl cluster is forced by metadynamics, revealing a transition path with intermediate metastable configurations in which planes are shifted one by one. The presence of metastable configurations in the transition path, corresponding to a certain number of NaCl planes turned into CsCl, has clear analogies with the known Hyde and O'Keeffe mechanism for ionic crystals, with some important differences due to finite-size effects. These comprise the fact that the transition starts by shifting a surface plane by means of a row-by-row mechanism that has no analog in bulk crystals. The energy barriers between the local minima in the transition path are calculated, showing that the barriers strongly depend on the screening length, in such a way that the B1 metastable phase can have very long lifetimes when the interaction is sufficiently long-ranged.

Kinetically driven ordered phase formation in binary colloidal crystals

The aggregation of binary colloids of the same size and balanced charges is studied by Brownian dynamics simulations for dilute suspensions. It is shown that, under appropriate conditions, the formation of colloidal crystals is dominated by kinetic effects leading to the growth of well-ordered crystallites of the sodium-chloride (NaCl) bulk phase. These crystallites form with very high probability even when the cesium-chloride (CsCl) phase is more stable thermodynamically. Global optimization searches show that this result is not related to the most favorable structures of small clusters, which are either amorphous or of the CsCl structure. The formation of the NaCl phase is related to the specific kinetics of the crystallization process, which takes place by a two-step mechanism. In this mechanism, dense fluid aggregates form at first and then crystallization follows. It is shown that the type of short-range order in these dense fluid aggregates determines which phase is finally formed in the crystallites. The role of hydrodynamic effects in the aggregation process is analyzed by stochastic rotation dynamics - molecular dynamics simulations, and we find that these effects do not play a major role in the formation of the crystallites.

Brownian dynamics simulations of colloidal suspensions containing polymers as precursors of composite electrodes for lithium batteries

Dilute aqueous suspensions of silicon nanoparticles and sodium carboxymethylcellulose salt (CMC) are studied experimentally and numerically by brownian dynamics simulations. The study focuses on the adsorption of CMC on silicon and on the aggregation state as a function of the suspension composition. To perform simulations, a coarse-grained model has first been developed for the CMC molecules. Then, this model has been applied to study numerically the behavior of suspensions of silicon and CMC. Simulation parameters have been fixed on the basis of experimental characterizations. Results of brownian dynamics simulations performed with our model are found in qualitative good agreement with experiments and allow a good description of the main features of the experimental behavior.

Oppositely charged model ceramic colloids: numerical predictions and experimental observations by confocal laser scanning microscopy

Fluorescent silica and alumina-like spherical particles with almost equal sizes are synthesized. Dilute aqueous suspensions are prepared with various ratios of those colloidal particles that exhibit opposite surface charges. These suspensions undergo heteroaggregation for a wide range of compositions. The structure of the formed aggregates is analyzed by means of confocal microscopy. The experimental results are compared to those of Brownian dynamics simulations in which the interactions between colloids are modeled by the DLVO potential. Good agreement between experiments and simulations is obtained.

Granulating titania powder by colloidal route using polyelectrolytes

A new, convenient, and inexpensive approach to process and granulate titania powders by a chemical route is proposed. It is based on the use of a formulation that includes a polyanion such as poly(sodium 4-styrenesulfonate) (PSS). Such a polyelectrolyte is most often considered to achieve dispersion of oxide powders in water. Basically, it adsorbs onto the surface of particles and induces electrical and/or steric interactions between particles in the suspension, which prevents agglomeration and rapid sedimentation. The advantages of polyelectrolytes in ceramic processing is well documented in the literature to produce low viscosity suspensions that are further used to form ceramic parts. In the case of TiO2 powders, such aqueous dispersions were obtained by adding small quantities of PSS. However, when exploring the behavior of mixtures containing lower contents of dispersant, we have discovered that, well below the optimum concentration required to get stable dispersions, the polyelectrolyte can act as a binder for titania particles. This can confer cohesion to the agglomerates, which can be processed to form large size (e.g., millimeter size) spheres. This phenomenon takes place when the oxide surface carries both positive and negative electrical charges and can be explained on a simple basis involving surface chemistry. For the optimum concentration of PSS that disperses titania, a polycation such as chitosan should be added to get spheres. This simple technique is expected to receive increasing attention due its potentialities and strong advantages with respect to other granulation techniques, such as spray-drying, which are energy consuming.

Heteroaggregation between Al2O3 submicrometer particles and SiO2 nanoparticles: experiment and simulation

The aggregation process of a two-component dilute system (3 vol %), made of alumina submicrometer particles and silica nanoparticles, is studied by Brownian dynamics simulations. Alumina and silica particles have very different sizes (diameters of 400 and 25 nm, respectively). The particle-particle interaction potential is of the DLVO form. The parameters of the potential are extracted from the experiments. The simulations show that the experimentally observed aggregation phenomena between alumina particles are due to the silica-alumina attraction that induces an effective driving force for alumina-alumina aggregation. The experimental data for silica adsorption on alumina are very well reproduced.

Agglomeration of alumina submicronparticles by silica nanoparticles: Application to processing spheres by colloidal route

In aqueous media, heterocoagulation between submicronic alumina (400 nm) and nanometric silica (25 nm) leads to the adsorption of silica on the alumina surface. By controlling the coverage rate of alumina particles, this adsorption destabilizes the suspension that leads to a very porous network of agglomerated particles. This work shows that the structure is all the more open as the density of charge carried by the two oxides is high and the ionic strength in the suspension low. From such a flocculated suspension, a new colloidal process to fabricate ceramic spheres is proposed which is based on a size increase of agglomerates. Under a controlled rotation of the vessel, electrostatic attraction between the surface charges of opposite polarity induces a size increase of agglomerates until the formation of spheres occurs. It has been shown that the mechanism of growth is poisoned by species adsorbed such as ions. Nevertheless, this new process proves very promising because it leads to a narrow size distribution of spheres by colloidal way, which can be subsequently consolidated by sintering, with a smooth surface.

Scaling and universality of self-organized patterns on unstable vicinal surfaces

We propose a unified treatment of the step bunching instability during epitaxial growth. The scaling properties of the self-organized surface pattern are shown to depend on a single parameter, the leading power in the expansion of the biased diffusion current in powers of the local surface slope. We demonstrate the existence of universality classes for the self-organized patterning appearing in models and experiments.