Ahmed
Hafez
,
Qi
Liu
* and
J. Carlos
Santamarina
Earth Science and Engineering, KAUST, Thuwal 23955-6900, Saudi Arabia. E-mail: qi.liu@kaust.edu.sa
First published on 4th June 2021
Self-assembly is ubiquitous at all scales in nature. Most studies have focused on the self-assembly of micron-scale and nano-scale components. In this study, we explore the self-assembly of millimeter-scale magnetic particles in a bubble-column reactor to form 9 different structures. Two component systems (N–N and S–S particles) assemble faster than one-component systems (all particles have N–S poles) because they have more numerous bonding pathways. In addition, two-components add control to process initiation and evolution, and enable the formation of complex structures such as squares, tetrahedra and cubes. Self-assembly is collision-limited, thus, the formation time increases with the total number of bonds required to form the structure and the injected power. The dimensionless Mason number captures the interplay between hydrodynamic forces and magnetic interactions: self-assembly is most efficient at intermediate Mason numbers (the system is quasi-static at low Mason numbers with limited chances for particle interaction; on the other hand, hydrodynamic forces prevail over dipole–dipole interactions and hinder bonding at high Mason numbers). Two strategies to improve yield involve (1) the inclusion of pre-assembled nucleation templates to prevent the formation of incorrect initial structures that lead to kinetic traps, and (2) the presence of boundaries to geometrically filter unwanted configurations and to overcome kinetic traps through particle–wall collisions. Yield maximization involves system operation at an optimal Mason number, the inclusion of nucleation templates and the use of engineered boundaries (size and shape).
All self-assembly systems involve: (1) elemental particles (e.g. molecules, DNA tiles, PDMS plates, 3D-printed tiles), (2) their interactions (e.g. electrical bonds, capillarity, magnetic), (3) a gas or liquid medium, (4) boundaries including interfaces, and (5) an energy source to sustain relative displacements that can lead to self-assembly e.g. mechanical vibration, heat, external field.12,16–27
Kinetic traps are intermediate states with local energy minima that hinder the evolution towards the most stable target configuration.28,29 Reversible bonds may help to overcome kinetic traps and correct errors during assembly.2,10,12 Nucleation templates – as in heterogenous nucleation – and geometrical constraints can guide self-assembly and improve yield;11,30,31 in fact, boundaries aid the self-assembly of lipid vesicles at the molecular scale32 and of 3D printed particles at the macroscale.27,29
Multi-component building blocks can form more complex structures than single-component systems and exhibit a wider range of properties;33,34 furthermore, multi-component systems allow for precise process initiation control. Salient examples include the self-assembly of DNA from thousands of distinct elements23,24 and the self-assembly of proteins and peptides to form supramolecular functional structures.33
Most of the work on self-assembly has focused on micron-scale and nano-scale particles controlled by electrical bonds and Brownian motion. Very few studies investigate the self-assembly of millimeter-size particles, where gravity, inertia and hydrodynamic interactions prevail. Yet, self-assembly can become a manufacturing technique for millimeter-scale objects that are inconvenient for robotic manipulation,18,29,35,36 especially when formation times are not a limiting factor.1 Furthermore, millimeter-scale self-assembly can provide valuable insight relevant to micron and nano-scale processes.
This study focuses on the self-assembly of millimeter-scale magnetic particles and seeks to understand concepts that apply to self-assembly at all scales, including assembly kinetics and the role of boundaries.
![]() | ||
Fig. 2 Experimental study. (A) Rod and tripod particles. Small disc-shaped magnets are fixed at the ends of each particle. The magnetic pole facing the particle is the hidden pole; brown and green magnets have exposed poles with opposite polarities. (B) Large bubble column reactor. (C) Possible initial nucleation structures for tripods and rods with probability of formation P (refer to Table S1, ESI†). |
These particle shapes and end-polarity allow us to explore 9 different structures (Fig. 1). One-component circular quadrants involve particles with N–S poles (Fig. 1A). On the other hand, two-component systems involve particles with the same end polarity, either N–N or S–S (Fig. 1B–G). The two-component three-dimensional tetrahedra and cubes combine tripods with three exposed S-poles and rods with two exposed N-poles (Fig. 1H and I).
Particles are suspended in a glycerol–water mixture (26% by volume glycerol, density ρf = 1.04 g ml−1, viscosity = 2 cP) inside a bubble-column reactor (length = width = 190 mm, height = 395 mm, Fig. 2B). Compressed air enters the bottom of the column through nine evenly-spaced needles (needle internal diameter = 0.66 mm, length = 14 mm). A high-resolution flow-meter (Brooks Instruments 1250-55) continuously monitors the air flow rate while a video camera (Sony) records the assembly process.
The test protocol used for all two-component tests conducted in this study consists of four steps: (1) fill the column with the viscous fluid, (2) start air injection, (3) release the first group of particles into the fluid and (4) add the second group. Particle image velocimetry measurements show that the average particle velocity increases with increasing air injection rate (see Fig. S1 for details, ESI†).37 We use a similar but smaller bubble column reactor (cylinder diameter = 94 mm) to study self-assembly in confined environments. The self-assembly of two-component systems can only start after the two components are mixed in the suspension. There is a weak interaction between the exposed S-pole in one rod and the hidden N-pole in the other, however, these incidental bonds readily break.
![]() | ||
Fig. 3 Assembly of one-component (blue) and two-component (red) circles. Cumulative distribution function for the time required to form the circles. The one-component circle is made of 4 identical N–S particles; the two-component circle is made of two N–N particles and two S–S particles (refer to Fig. 1A and B). Points: experimental data. Dashed lines: log-normal distribution. Mean ± standard deviation – one component: 16.8 ± 12.1 min. Two components: 10.2 ± 12.6 min. Air injection rate = 3.7 l min−1. |
Let's consider the assembly pathways for the four particles involved in one- and two-component circles. In the one-component system, the N-pole in one particle can bind to the S-poles in any of the three other particles. However, the N-pole in two-component N–N particles can bind to four S-poles in S–S particles. The complete “bonding graphs” for both systems show that the two-component system has a total of 128 possible assembly pathways, while the one-component circle has only 48 possible pathways (Fig. S2, ESI†). Hence, the two-component circle experiences more efficient self-assembly (Fig. 3).
Formation times ranged from 40 seconds to more than 10 hours. Fig. 5A shows that the formation time increases with the total number of bonds Nt; in fact, the higher the number of bonds Nt, the higher the number of intermediate structures and kinetic traps the system encounters during self-assembly. We repeat the experiments for each target structure ten times and define the product yield as the ratio of the number of trials that result in successful error-free structures within 10 hours to the total number of 10 trials attempted. The yield is 100% for the 2D shapes (circles Nt = 4, triangles Nt = 6 and squares Nt = 8) but significantly lower for the 3D shapes (tetrahedra Nt = 12 and cubes Nt = 24, Fig. 5B).
The initial “nucleation structure” plays a critical role in the formation time and yield of complex structures. In fact, incorrect initial nucleation structures are responsible for the low yield in tetrahedra and cubes. Fig. 2C shows four equally stable nucleation structures formed by tripod and rod particles during cube assembly. We can conclude that (1) the nucleation of the correct initial structure has a low probability P = 12.5% – Table S1 (ESI†), and (2) alternative initial structures are equally stable and particles remain bonded in a kinetic trap that would require a high input energy transient to overcome.
We perform additional experiments where we introduce a pre-assembled nucleation template made of four tripods and four rods into the reactor (structure 4 in Fig. 2C), in addition to the remaining four tripods and eight rods. In this case, the yield improves significantly from 0.2 to 0.8 (10 trials each; air injection rate = 2.0 l min−1 for trials with and without template).
Some experiments involve the particles needed to generate one structure only; this initial condition minimizes kinetic traps. In fact, the outcome is not unique when experiments involve the particles needed to form two structures (Fig. 1E–G). For example, the 6 vertices and 6 rods can either assemble into two separate triangles or conjoined triangles.
![]() | (1) |
We explore the assembly kinetics of 8 tripods and 12 rods in the large bubble column reactor (note: the asymptotic structure is the cube – Fig. 1I). We subject the system to five different air injection rates and repeat each experiment three times (Movies S5 and S6, ESI†). Results in Fig. 6 show that the number of bonds increases asymptotically with time, in agreement with previous studies.40 A hyperbolic saturation curve adequately describes the data:
![]() | (2) |
![]() | ||
Fig. 6 The effect of kinetic energy on bond creation. Target structure: cube made of 8 tripods and 12 rods (large bubble-column reactor). The air injection rate varies from 2.6 l min−1 to 9.7 l min−1. A hyperbolic saturation curve fits the data, where t* is the characteristic time to form nbt/2 of the bonds. We repeat experiments three times at each air injection rate. Inset: The inverse of the characteristic time (1/t*) as a function of the air injection rate. Points: experimental data. Dashed lines: upper and lower bound fits using eqn (4). |
These results resemble the effect of heat on crystal nucleation and growth rate.41,42 The classical nucleation theory estimates the nucleation rate J [# events m−3 s−1] as a function of the number of molecules near the nucleation site and the probability that molecules have enough energy to overcome the formation barrier:43,44
![]() | (3) |
An analogous analysis applies to millimeter-scale magnetic particles. There are two competing trends. At low flow rates, the assembly rate Ja [# events m−3 s−1] depends on the particle collision rate jc [# events m−3 s−1] which is proportional to the air injection rate q [m3 s−1] and particle concentration n [# particles m−3]; let's adopt an asymptotically linear relation jc = α·n·q. At high flow rates, the jc-collisions will not render stable bonds if the kinetic energy resulting from the hydrodynamic force K [J] exceeds the binding energy arising from magnetic dipole interactions U [J]. If we assume that particle velocities follow a Maxwell–Boltzmann distribution,45 then the probability that K > U is βexp(−U/K) where the β-factor from the Maxwell–Boltzmann distribution has a weak dependence on K and U, and the kinetic energy K = 1/2·m·vp2 depends on the air injection rate q. The two trends combine and anticipate an optimal air injection rate (Fig. 6, inset):
![]() | (4) |
We test the effect of confinement by varying the confinement factor CF defined as the ratio of the target structure's characteristic size to the diameter of the reactor (experiments: we varied the triangle side length and used the small reactor in all tests). Fig. 8 shows that the formation time and yield for the self-assembly of triangles have an optimum value at an intermediate confinement factor CF. An under-confined system (low CF) results in unrestricted structural geometries, as well as low particle volumetric concentration and low particle collision rate (eqn (4)). Assembly time and yield improve as the confinement factor increases. Particle–wall collision frequency increases with increasing confinement factor, and causes more frequent breakage of partial structures. While disassembly allows recovery from kinetic traps, frequent structure-wall collisions in an over-confined reactor cause continuous fragmentation and delay structure formation. Clearly, an optimum confinement factor minimizes formation time and maximizes yield. (Note: air injection favors particle movement in the vertical direction and adds an additional constrain or bias – see Fig. S1 and Movies S1 through S6, ESI†).
We can anticipate that the reactor shape also affects the evolving morphology of the assembled structure. In fact, cubic nanoparticles self-assemble into cubic superlattices in under-confined systems46 and into spherical structures in confining emulsion droplets.47 Similarly, copolymers assemble into thin films when confined between parallel plates, into cylinders when confined in pores with square and triangular cross-sections, and bend into helices and stacked toroids in cylindrical pores with circular cross-sections.48 Furthermore, different structures emerge under moving rather than static boundary conditions.49 This suggests the possibility of feedback-controlled boundary conditions to optimize self-assembly.
The self-assembly of millimeter-scale particles is collision-limited. The optimal kinetic energy for self-assembly is a trade-off between collision rate, bonding probability and debonding out of energy traps. Self-assembly is most efficient at intermediate Mason numbers.
Formation time increases and yield decreases with structural complexity. We explored two possible strategies to improve yield. Pre-assembled initial nucleation templates prevent the formation of incorrect initial structures that lead to kinetic traps; therefore, nucleation templates result in significant yield improvements. Equally important is the role of boundaries, including interfaces: boundaries act as geometric filters that exclude unwanted configurations, constrain assembly pathways and help overcome kinetic traps through particle–wall collisions. Geometric boundaries can be engineered (size and shape) to optimize self-assembly.
J [# events m−3 s−1] | Nucleation rate |
j [# events per s] | Attachment rate |
J a [# events m−3 s−1] | Assembly rate |
j c [# events per s] | Collision rate |
K [J] | Kinetic energy of a particle |
k B [J K−1] | Boltzmann constant |
Mn [ ] | Mason number = hydrodynamic drag/magnetic dipole–dipole interaction |
n [1 m−3] | Particle concentration |
n b [ ] | Number of bonds formed at time t during self-assembly |
n bt [ ] | Asymptotic number of bonds the particles form during self-assembly |
q [m3 s−1] | Air injection flow rate |
T [K] | Temperature |
t* [s] | Characteristic time to form nbt/2 bonds |
〈t〉 [min] | Average formation time |
U [J] | Magnetic potential energy for dipole–dipole interaction |
v p [m s−1] | Average particle velocity |
α [1 m−3] | Proportionality factor between collision rate and injection flow rate |
β [ ] | Factor in the Maxwell–Boltzmann distribution |
ΔG* [J] | Formation Gibbs free energy for nucleation |
ρ f [kg m−3] | Fluid density |
ρ p [kg m−3] | Particle bulk density |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sm00588j |
This journal is © The Royal Society of Chemistry 2021 |