Evaporation-driven colloidal cluster assembly using droplets on superhydrophobic fractal-like structures

Colloidal particles are assembled into structures of various shapes – from flattened arrangements to three-dimensional spheroids – by confining them in droplets and manipulating their evaporation over a microfabricated surface with fractal-like pillars that adopts superhydrophobic character when properly coated. Covering a wide range of initial particle concentrations, we demonstrate its dominant role over other parameters such as the ambient humidity. As featured in:


Introduction
Agglomerating small solid particles in clusters of a particular shape and size has been an intriguing challenge of fundamental soft matter research for decades. Precise knowledge of the underlying mechanisms and the governing factors of, for example, grain size and morphology, is a necessity for many fields of modern microtechnology, such as spray drying of colloidals, 1,2 (3D) printing, 3 coating processes, 4 and drug delivery. 5 Handling micro-or nanoparticles in a dry state is challenging for many applications, due to adhesion forces and electrostatic charges. Often, it is easier to first suspend the particles in a liquid and then gain control over the evaporation of the suspending phase. At controlled evaporation, the particles can agglomerate into clusters of bigger size, which are easier to handle. Evaporating suspensions is a powerful tool for particle deposition, as was observed in works as early as Brown's, deepening his understanding of the motion now named after him. 6 In those early works, Brown already observed ''. . .those currents from centre to circumference [of the droplet], at first hardly perceptible, then more obvious, and at last very rapid, which constantly exist in drops exposed to the air. . .''. More than a century later, Deegan et al. 7 explained the phenomenon, which is now popularly known as the coffee-stain effect: particles are being transported to the rim of a sessile evaporating droplet by ''currents from centre to circumference'', and form a stain-like deposit, once the evaporation is complete. A crucial condition for this phenomenon to occur, is a stationary droplet perimeter at a contact angle below 901 with respect to the substrate. This leads to a replenishment of liquid at the outer rim in order to maintain a minimal liquid-air interface. The droplet perimeter, also described as the (three-phase) contact line, remains stationary or pinned as long as surface defects or surface-adhered particles prevent the contact line motion. 8 However, if the contact line of the droplet is unpinned and moves freely across a pristine substrate, the coffee-stain effect will be mostly suppressed. 9,10 In order to agglomerate particle clusters from suspension droplets via evaporation, all particles should ideally end up agglomerated around one point, e.g. at the centre of the droplet. To avoid the irreversible transport of particles away from the droplet centre, an unpinned contact line and high contact angles are therefore essential requirements for a particle cluster formation. Both requirements can be met by evaporating the droplets on superhydrophobic substrates. [11][12][13][14][15][16][17] On a superhydrophobic substrate, water-based droplets minimize their contact area with the solid due to the considerably higher solid-liquid surface energy than the droplet's liquid-gas energy. The increased energy difference is often induced by a chemical treatment of the solid substrate (see Section 3.1 for more details). To reduce the contact between liquid and solid as far as possible, the substrate is often decorated with micro-or nano-scopic asperities, which lead to a composite solid-liquidgas interface underneath the droplet. The droplet then lays on top of a forest of microstructures, like a fakir, laying on a bed of needles. [18][19][20] Such a droplet features high macroscopic contact angles with the substrate and little or no contact line pinning, as long as the pattern of microstructures is homogeneous and does not suffer from impurities or irregularities.
This superhydrophobic state of almost floating droplets is well described by the Cassie-Baxter model. 12 Consequently, a droplet following this model is said to be in a Cassie-Baxter (CB) state. In contrast, when the droplet is fully engulfing the microstructures, the model developed by Wenzel 11 applies, and the state in which the droplet is impaled on the asperities is called the Wenzel (W) state.
CB states typically involve larger surface energies than W states. This translates to a certain degree of metastability and larger contact angles in CB states than observed in W states. Droplets transition between the higher energy CB to the lower energy W by impaling on the asperities. 21,22 However, this transition is rarely spontaneous and requires a certain input of energy due to the complex energy landscape. 23 Impaling can happen due to, for example, evaporation and a consequent pressure increase inside the droplet, 21 or due to mechanical force, pushing the droplet down into the substrate. 24 Reentrant features on the asperities 14,25,26 can help to prevent the droplet from impaling.
Water droplets on superhydrophobic substrates are therefore ideal vessels for (non-volatile) colloidals, such as e.g. particles, polymers, and biological matter, due to the reduced contact with the solid substrate and the reduced interaction with the gas phase. But when the droplet evaporates, such ideal containers can become rather complex systems. Water evaporation leads to a continuous increase in volume concentration of non-volatile contents. As the droplet's volume reduces, nonvolatile content can eventually aggregate in different structures depending on their nature. 27,28 As a result, generating particle clusters from sessile suspension droplets is a simple way (due to the low energy input) to achieve agglomeration of particles on a very small length scale.
Despite the abundant literature on such systems, most of these studies are typically very narrow in the parameter space explored, and therefore there is a lack of understanding on which are the most important parameters to control the shape of the resulting colloidal structures via droplet evaporation.
In our study, we demonstrate that, using a robust superhydrophobic substrate and a stable colloidal solution, one can obtain a vast variety of particle assemblies by simply controlling the initial particle concentration (and to a lesser extent, humidity). Our results show two extreme kinds of cluster shapes: spherically-shaped (3D) clusters and flattened (2D) clusters, which we classify quantitatively based on their shape. On the one hand, 3D clusters offer high particle packing fraction and feature low surface area to volume ratios. On the other hand, 2D clusters offer direct access to the majority of particles in the deposit, and are of special interest for 2D printing applications. 29

Methods
The performed experiments involve the evaporation of sessile droplets on a superhydrophobic substrate. The experiments take place inside a closed (not pressurized) chamber to minimize environmental influence. See Fig. 1 for a schematic of the set-up. The substrate is placed in line with the CCD camera (Ximea MQ013MG-ON) to capture the side view of the experiment. To deposit a droplet on top of the superhydrophobic substrate, we use a threaded-plunger syringe (glass syringe from Hamilton, Model 1750 LT Threaded Plunger SYR) with a blunt cut needle with a nominal outer diameter of 0.21 mm (also supplied from Hamilton, 33 gauge, Metal Hub Needle).
In order to gently deposit the droplet, the needle is brought close to the substrate. We then form a droplet at the needle tip by slowly turning the threaded plunger, until the droplet touches the substrate and the needle at the same time. Once the droplet volume is big enough, we force the detachment of the droplet by retracting the needle upwards. The droplet consequently comes to rest on top of the substrate. With this method, used for all experiments, we reach droplet volumes between 3 ml and 4 ml. Our analysis show that such variation in initial volume does not affect the results.
The deposited droplets are either purified water droplets (Milli-Q s IQ 7000 Water Purification System), or water droplets with suspended fluorescent polystyrene particles of 0.98 AE 0.04 mm in diameter. We use suspensions of different particle concentrations, ranging from 0.0009 vol% to 1.66 vol% (the error margin for all concentrations stated in this work are AE10%). The fluorescent polystyrene particles (distributed by Microparticles GmbH) are functionalized with sulphate groups that convey electrostatic stability to the suspension with no need of additional surfactants. Prior to the experiments, the particle solutions are washed twice by diluting, centrifuging and extracting the supernatant. After the washing process, the particles stay homogeneously suspended for longer than 24 hours. The suspensions are stored in the fridge and are immersed in an ultrasound bath for 10 minutes prior to each experiment. In every experiment, the evaporation process is captured from the side view with the CCD camera, with a frame rate of 0.5 frames per second (see Fig. 2 for sample images). From these images, we extract contact angles and droplet volumes as functions of time, as well as aspect ratios of particle clusters, after the liquid has completely evaporated. We categorize the clusters based on their morphology (see Section 4.2). For this classification, we also rely on scanning electron microscopy (SEM), see for example Fig. 5.
The role of the chamber is to shield the droplets against spurious air currents or sudden changes in temperature and humidity. For a number of experiments, we also use a humidity controller (HGC 30 from DataPhysics) to keep the humidity in the chamber at a constant value. The values for humidity and temperature within the chamber were monitored throughout all experiments.

Fractal-like substrates
Typical superhydrophobic substrates have two components which render them superhydrophobic: 26,30,31 a mechanical component, consisting of microscopic structures, that enlarge the surface area and enhance air entrapment, and a chemical component, consisting of a low-surface-energy coating. The substrates in our research were fabricated with oxide-only corner lithography and an anisotropic etching process of silicon, 32,33 which delivers SiO 2 substrates with microscopic fractal-like structures (see Fig. 3). The fractal-like structures consist of stacked octahedra, resulting from the crystalline structure of silicon. The intricate shape of the fractal-like asperities features several re-entrant geometries (see Fig. 3c). These ''overhanging'' steps are known to considerably improve the stability of the floating CB state of droplets. 14,25,34,35 The base width of one fractal-like structure is 5 mm, the structures are ordered in a hexagonal pattern on the substrate and the pitch between structures is 12 mm. The actual surface area of this substrate is 13 times bigger than the projected one.
Following Young's law of surface tensions, 36,37 the high surface energy of SiO 2 renders substrates of this material generally hydrophilic. Decorating the surface with microscopic structures, and therefore adding a roughness, will increase the wettability further, 38 much like wetting a porous medium. However, applying a low surface energy coating to the surface renders the substrate superhydrophobic. To achieve this, we applied a vapour deposition protocol to coat our fractal-like

Experimental results
Since the evaporation process is the driving force for particle agglomeration, we first analyse the droplet evaporation process. Evaporating sessile droplets either shrink with a constant contact radius of the wetted area on the solid substrate (CCR mode), or with a constant contact angle, while the wetted area is decreasing (CCA mode). 41,42 A droplet evaporating in a CCA mode maintains an unpinned contact line, with the radius of the wetted area decreasing due to decreasing droplet volume. If the contact line pins (and the droplet enters a CCR mode), the contact angle will need to decrease due to the continued volume loss. If the contact line de-pins again, the contact angle will jump back up to its original receding value. We do not observe significant pinning events on our substrates throughout the evaporation of the droplets, as is shown in Fig. 2 and its inset, where we can see that the change in contact angle as a function of time is negligible until the last instants. Our experiments therefore follow, for the majority of their duration, the CCA mode with contact angles of E1551 due to the superhydrophobicity of the substrate, which prevented the droplet perimeter from pinning.
Both modes of evaporation have been extensively studied. In this analysis, we will make use of the celebrated work by Popov, 40 who investigated spatial deposition patterns of particles during the evaporation of sessile suspension droplets in CCR mode. However, since the derivation of evaporation rates is valid for arbitrary contact angles, this model has also been adapted for cases with droplets in CCA mode, e.g. by Nguyen et al. 43 and Dash and Garimella 15 and we implement their analytic approach in the following.
The evaporation rate can be described as a function of contact radius and contact angle: where V is the droplet volume, t is time, D is the diffusion coefficient of vapour in air, c s the saturation concentration of water in air, H r the relative humidity, y is the contact angle of the droplet with the substrate, r the liquid density, R c the base radius of the droplet and f (y) a contact angle dependent function for arbitrary contact angles, as defined by Popov: 40 where t is the integration variable. We evaluated the integral in eqn (2) numerically from t = 0 to t = t max , choosing a large enough t max value, such that the integration becomes independent of the choice. Assuming the droplet's surface adopts a spherical cap shape, we can calculate the volume as a function of y and of R c of the droplet. From this, we extract the base radius as a function of volume: We want to emphasize that for all our experiments, the contact angle stays approximately constant and well above 1501 throughout 90% of the droplet's lifetime. We therefore consider the functions f (y) and g(y), (eqn (2) and (4)) also constant for the majority of the duration of the experiments.
By numerically integrating eqn (2), and combining the result with eqn (3) into eqn (1), we obtain the following relation for the volume of the droplet in time: with In order to make a prediction of the droplet volume evolution, we require empirical input from the initial droplet volume, the initial contact angle and the relative humidity in the chamber during the experiment. However, we noticed that the measured relative humidity does not seem to be representative of the far field vapour concentration felt by the droplet, and therefore we choose to use H r as a fitting parameter. We find that our measured values are consistently off for all experiments by a factor of roughly B1.2, which might be caused by the fixed position of the humidity sensor.
In Fig. 4 we present the results of the measurements, expressed as volume loss over time, based on eqn (5). Fig. 4a shows data from droplets with and without added particles, conducted at similar relative humidities. As the plots overlap within the experimental error, we can safely conclude that the presence of particles at these initial concentrations does not influence the evaporation process. This holds even for the last stages, when the particle concentration reaches its highest values.
In Fig. 4b, we show experiments with pure water droplets, conducted at different relative humidities. An increase in relative humidity reflects in a decrease in evaporation rate, manifested in smaller slopes of the curves. In order to test the validity of the modelled evaporation rate for all experiments, we write eqn (5) in dimensionless form as where Ṽ = V/V 0 andt R eff,0 is the effective radius of the droplet at time t = 0. Note that since g(y) and f (y) are only functions of y, they can be considered constant for the majority of the droplet's lifetime.
In Fig. 4c, we show the same experimental data from panels (a) and (b) in dimensionless form. Eqn (7) is plotted making use of the mean y found in the experiments. The experimental data from the different experimental settings clearly collapses onto a single curve, matching also the theoretical predictions from eqn (7), and therefore we can conclude that our experiments are all governed by a diffusion-limited regime.

Morphological categories
By varying the initial particle concentration of the evaporating droplets within more than three decades, we find a variety of shapes in the resulting particle clusters. We categorize them in three groups, depending on shape and formation, which are shown in Fig. 5. The two extreme shapes are on the one hand thin, flattened clusters (''pancakes'') and on the other hand, spheroid-shaped clusters with a smoothly curved surface, which we will describe as ''spheroids''. In between these two extreme cases, we find clusters with irregular shapes which we will refer to as ''buckled clusters''. This latter category of agglomerates experiences a sudden buckling instability in a late stage of evaporation.
In order to better quantify the classification, we introduce the aspect ratio AR of the final cluster as AR = d v /d h , calculated from the side view images, where d v is the maximum vertical extend of the cluster (normal to the substrate) and d h the maximum horizontal one (parallel to the substrate).  (7), which is plotted making use of the mean y found in the experiments. All plots collapse, indicating that our experiments are governed by a diffusion-limited evaporation regime which can be described by the analytical solution by Popov. 40 Fig. 5 Overview of the three different categories of cluster shapes. All panels present a side view schematic, a top view SEM image of a typical cluster and a zoom-in view of the cluster surface. The particles are 0.98 mm in diameter. (a) Flattened clusters, called pancakes, with low aspect ratios and no considerable height extent. The microscopic fractal-like structures are in some areas visible through the cluster. Note the perfectly circular shape. Experimental details: F 0 = 0.0009%, H r E 78%. (b) Buckled clusters with folded rims and a hollowed core (not necessarily visible from the side view), which we show in the schematic with a dashed line. Experimental details: F 0 = 0.009%, H r E 60%. (c) Clusters with an intact core, called spheroids, which present a more regular shape and a smooth cluster surface. We find extensive patches of hexagonally ordered particles at the cluster surface. Experimental details: F 0 = 0.045%, H r E 78%.
The main control parameter is the initial particle concentration or initial particle packing fraction F 0 , defined as the fraction of volume taken up in the droplet by the solid particles. As the evaporation process takes the place, the particle packing fraction F(t) will increase until its maximum, which will change depending on the final particle arrangement. A perfectly ordered arrangement of the particles in close packing would yield a F| max of approximately 0.74. Unfortunately, due to the irregular shapes of the agglomerates and inaccuracies in the total number of particles per droplet, we are unable to compute their final packing fraction with sufficient precision. Fig. 6 a shows the resulting aspect ratio AR as a function of the initial particle concentration F 0 . Starting with the lowest AR (F 0 = 0.0009%, about 8 Â 10 4 particles in 4 ml), at the left side of the figure, we find the pancake category, with aspect ratios well below 0.4.
Increasing the initial particle concentration to the range from 0.0045% to 0.009%, we start to find buckled clusters, which results in AR values of up to 0.9 in this range of packing fractions. As we continue further along the x-axis of the plot, towards the largest initial particle concentrations (F 0 Z 1.66%), buckled clusters and spheroids, are both found in a similar range of AR values. Buckled clusters present a large variety of shapes, which explains their larger dispersion in AR. Note however, that the AR measurements are based only on a single side view image of each object, and therefore such measurements should be interpreted with caution. While spheroids and buckled clusters are found with similar probabilities for initial concentrations in the range from 0.009% up to 0.5%, regardless the humidity values, we do see a sharp dependence on the evaporation rate for the largest initial concentrations. This is shown in Fig. 6b, where spheroids are obtained with the highest probability at high humidity (slower evaporation rates), while buckled clusters are found with their larger dispersion in AR for lower humidity values (faster evaporation rate).
The adherence of the clusters to the substrates is an interesting property and worth mentioning: whereas the flattened pancakes stick to the substrates, clusters that consolidate a more 3D shape detach easily from the substrate after evaporation. Simple tilting of the substrate out of the horizontal plane will be enough to displace the clusters from their original position. Such lack of adherence was also observed by Wooh et al. 17 for titanium oxide nanoparticles over soot-template superhydrophobic substrates.
In summary, our analysis indicates that the initial particle concentration in the droplet is a (if not the most) crucial factor in determining the final shape of the agglomerate. There are two interesting transitions that we will discuss in the following section: First, the fact that only flattened clusters are observed below packing fractions of 0.004%. Second, the puzzling transition from buckled clusters to stable spheroids, which does not occur for a certain critical initial particle concentration, but rather for a combination of F 0 and humidity.

Discussion
Our results show that capillary water droplets (i.e. where the dominant length scale is smaller or comparable to the capillary length) evaporate following a CCA mode on surfaces decorated Fig. 6 (a) Characterization of the clusters according to aspect ratio AR = d v /d h , and the initial particle concentration F 0 . Note the logarithmic scaling of the axes. The inset schematics and the denominations correspond to Fig. 5. With increasing F 0 the final clusters present with higher aspect ratios. The distinction between buckles clusters and spheroids is not always sharp, as can be seen in the overlap of particle concentrations, where both categories co-exist. The data for high F 0 is plotted again in panel (b), as AR = d v /d h over relative humidity H r . Note here the linear scaling of the axes. We can see a clear trend for buckled clusters to present for lower H r than spheroids. This trend is not as clear in lower F 0 . Therefore, we assume that a sufficiently high initial particle concentration is crucial for the influence of humidity to become quantifiable.
with fractal-like micropillars. The evaporation process follows a purely diffusive model, based on a modified version of the Popov model. 40 The presence of particles does neither affect the droplet evaporation, nor the contact line motion. Even the droplets containing the highest initial particle concentration still evaporate following a diffusion-limited process (Fig. 4a), and the contact line still follows a smooth CCA motion. As an indication of smooth contact line motion, we want to highlight the fact that no particles can be found outside the bounds of the final cluster. The suspended solid phase in our droplets is therefore completely recovered in a single agglomerate (supraparticle) via the evaporation of the suspending liquid phase. This is an important difference to previous studies on supraparticles on microstructured substrates 14,20,44 in which -after reaching a certain particle concentration in the droplet -particles were deposited on top of pillars in the wake of the receding droplet perimeter, as a consequence of a regular stick-and-slip contact line motion.
We do not observe any stick-and-slip motion of the contact line (also see Fig. 2, inset) and instead, the droplet's movement across the substrate is smooth, until a final pinning event occurs at the end of the evaporation process. This final pinning event also limits the 3D cluster agglomeration. Droplets with the lowest initial particle concentration (F 0 = 0.0009%) do not consolidate a robust particle cluster before the final pinning event occurs. The particle deposit then adopts a flattened circular shape (as seen in Fig. 5a), where the perimeter is determined by the final pinning event, however the height keeps decreasing further after the pinning happened.
The evaporation in CCA mode, with the droplet's base radius steadily decreasing until the final pinning event, gives an indication, that the droplet remains in a Cassie-Baxter state during the whole evaporation process. SEM close-ups of pancake clusters (see Fig. 7) show that the particles are partially engulfing the fractal-like structures. Since there is no evidence for particles at the bottom of the substrate between the fractallike structures, the liquid seems to never fully penetrate the microstructured substrate. The pancake deposits indicate that the droplet gets partially impaled on the microstructures in a late stage of the evaporation. The particles, wrapping the fractal-like structures only on the top half, then indicate the liquid penetration level. As the final pinning event determines the extent of the pancake perimeter, final pinning and partial impalement possibly happen at the same time and are actually part of the same phenomenon: droplets in this small size range must be partially impaled by the structure, 22 which eventually hinders the contact line motion. Wrapping the top half of the microstructures locks the pancakes in their final drying position, which also leads to these clusters to not easily detach, but to stick to the solid substrate.
As the initial concentration of particles increases from F 0 = 0.0009% to 0.0045%, we observe a transition in cluster shapes: particles begin to consolidate three-dimensional clusters, without influence of the final pinning event.
Buckled clusters appear with variable probabilities for all initial concentrations above F 0 = 0.0045%. In comparison to the pancakes, these clusters maintain a 3D shape, even though their height can decrease considerably during the buckling event. This leads to folded rims and hollow cores, as can be seen in Fig. 5b. Buckling instabilities in evaporating droplets have often been described as mechanical shell instabilities. [45][46][47][48] Predictions for critical buckling conditions can be found for spray drying 49 and Leidenfrost drying; 45,46 as well as for the dissolution of water-based suspension droplets in oil. 50,51 SEM images of buckled clusters (e.g. Fig. 5b) indicate that the buckling in the experiments of this study also has its origin in a shell formation at the interface of the droplet. 52,53 However, a rigorous quantification of the mechanical shell (in-)stability during the evaporation process is outside of the scope of this work.
All 3D clusters (buckled or unbuckled) reveal highly-ordered particle patches at their surface (see the zoom-in views in Fig. 5), similar to those previously found for clusters of similar size. 20 We interpret this as evidence of a predominant particle accumulation at the droplet surface (or in its vicinity). To confirm this we compare the particle diffusive time scale t B = r p 2 /D p -where r p is the particle radius and D p its Brownian diffusivity -with the time scale at which the interface approaches the particles. To construct the time scale of the interface we use the inwards radial velocity of the interface U surf (t) -which increases non-linearly with time -as well as the particle radius for a length scale, which yields t surf = r p /U surf (t). Comparing the two time scales yields a maximum value of t B /t surf B 0.5, from which we interpret that the interfacial motion is always slower than (but close to) the particles' Brownian motion. Particle-particle ordering, which leads to the hexagonally-ordered facets we see at the cluster surfaces, is therefore not directly influenced by the motion of the interface. However, the comparison of time scales reaches values close to O(1). Therefore, the inwards motion of the interface is fast enough to have a sweeping effect on the bulk particles in its way. With time, the particle concentration increases in the region ahead of the moving droplet interface. This picture is consistent with the presence of a loose shell being formed during the evaporation close to the droplet's interface. The hexagonally-ordered facets on the clusters surface originate in this loose shell, which bends and eventually buckles. 52 For initial particle concentrations above F 0 = 0.045%, some of the clusters assemble into spheroids, for which no buckling event can be detected. Considering that buckling is based on a mechanical shell instability, the particle clustering of spheroids has to be mechanically stable enough to withstand rising tensions. In this study, where spheroids only occur for F 0 Z 0.045%, experiments suggest that a higher initial particle concentration favours the agglomeration into spheroids. With a higher initial particle concentration, the sweeping effect of the interface on the particles in the bulk will lead to a shell, which can extend further into the bulk and therefore convey mechanical stability.
In principle, buckling events, i.e. unstable shells, are likely to occur for all the initial particle concentrations explored. However, for experiments with high initial particle concentrations F 0 Z 0.5%, we observe a strong influence of the relative humidity, see panel (b) in Fig. 6, which is absent for lower concentrations. Spheroids are then found with the highest probability when the humidity is increased from approximately 60% to about 80%. This increase corresponds to a proportional decrease in the evaporation rate, and consequently in the time scale ratio t B /t surf to values closer to 0.3. We conclude from this result that the sweeping effect, too, substantially decreases, which only has a significant impact in those cases, where the inter-particle distance is shortest (i.e. for highest particle concentrations).
Other studies on evaporation-driven colloidal agglomerates have obtained other cluster shapes, like the ''doughnut'' shapes in Rastogi et al., 16 which we have not observed in our experiments. It is difficult to point to one specific reason to explain such a difference, but we can definitively point out that their doughnut shapes were found for heavier particles (silica is almost twice as dense as our polystyrene particles), much larger initial concentrations and, as explicitly mentioned by the authors, an early pinning of the contact line. Nonetheless, the colloidal agglomerates found by Rastogi et al. 16 at their lowest particle concentrations seem very similar to the spheroidal clusters in our work. Wooh et al. 17 also occasionally found doughnut shapes in their titanium oxide particle clusters, for small droplets in which the contact line motion cannot keep up with the fast evaporation. Our clusters however are very different from those by Wooh et al. 17 In their case, the final particle packing is rather low due to the complex inter-particle structures resulting from the attractive interaction among the particles.
In summary, the data in Fig. 6 shows that, for high initial particle concentrations, a slower evaporation leads to a higher packing of particles inside the shell, conveying stability to the cluster. Interestingly, such a dependence on the evaporation rate has also been observed in spray drying systems for silica nano-particle solutions with particle sizes in the range of 10 nm, with initial volume fractions in the range of 2%, 49 but with completely different evaporation rates. Furthermore, in a similar system as ours, Sperling et al. 53 directly observed with confocal microscopy such shell formation, using fused silica nano-particle solutions at slightly higher initial particle concentrations, as the ones used here (around 1.7% volume fraction).

Conclusions
We have shown that, using a stable droplet suspension and a robust superhydrophobic substrate, we can obtain a large variety of particle agglomerates depending only on the initial particle fraction within the droplet. The fractal-like microstructures employed have proved to be a powerful tool yielding high contact angles and minimized pinning of evaporating water droplets. Although wetting transitions from CB to W states were not observed, final pinning events for minimal droplet sizes seem to be unavoidable, 22,35 and prevent the formation of spheroid-shaped particle clusters with low number of particles. We distinguish three different categories of cluster shapes based on their aspect ratio (i.e. their three-dimensional character), for a wide range of initial particle packing fractions in the droplet, spanning over 3 decades. From low to high initial particle concentration we found: pancakes, flattened clusters which engulf the upper part of the micropillars, buckled clusters, which adopt their shape due to a sudden collapse of the cluster's shell structure, and spheroids, unbuckled 3D clusters with smooth surfaces.
The controlled assembly of shell structures by exploiting a simple process, such as the evaporation of particle suspensions, is a very exciting goal. Applications such as spray drying, 1 drug delivery, 5 or advanced detection techniques 54,55 utilize shell geometries down to the nano-scale. The ability to control the shell thickness and stability is crucial, as well as insights about stratification 56 and/or cavity formation on the inside of the cluster.

Conflicts of interest
There are no conflicts to declare.