Open Access Article
Tonghan
Gu
a,
Cao
Zheng
a,
Fan
He
a,
Yunfei
Zhang
a,
Saif A.
Khan
b and
T. Alan
Hatton
*a
aDepartment of Chemical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA. E-mail: tahatton@mit.edu
bNational University of Singapore, Department of Chemical and Bimolecular Engineering, 4 Engineering Drive 4 E5-02-28, 117576 Singapore. E-mail: saifkhan@nus.edu.sg
First published on 28th March 2018
Microfluidic droplets have been applied extensively as reaction vessels in a wide variety of chemical and biological applications. Typically, once the droplets are formed in a flow channel, it is a challenge to add new chemicals to the droplets for subsequent reactions in applications involving multiple processing steps. Here, we present a novel and versatile method that employs a high strength alternating electrical field to tunably transfer chemicals into microfluidic droplets using nanodroplets as chemical carriers. We show that the use of both continuous and cyclic burst square wave signals enables extremely sensitive control over the total amount of chemical added and, equally importantly, the rate of addition of the chemical from the nanodroplet carriers to the microfluidic droplets. An a priori theoretical model was developed to model the mass transport process under the convection-controlled scenario and compared with experimental results. We demonstrate an application of this method in the controlled preparation of gold nanoparticles by reducing chloroauric acid pre-loaded in microfluidic droplets with L-ascorbic acid supplied from miniemulsion nanodroplets. Under different field strengths, L-ascorbic acid is supplied in controllable quantities and addition rates, rendering the particle size and size distribution tunable. Finally, this method also enables multistep synthesis by the stepwise supply of miniemulsions containing different chemical species. We highlight this with a first report of a three-step Au–Pd core–shell nanoparticle synthesis under continuous flow conditions.
Many chemical reaction systems, such as reactions generating solids14–16 or those requiring a rate-controlled addition of one or more chemicals, require the addition of chemicals into generated microdroplets. Multistep preparation in microdroplets has also only been achievable with such techniques.17 The most common addition strategy is to alternately generate microdroplets of similar sizes in pairs carrying two reactive species, and fuse them together by mechanical or electrical stimuli,18,19 which typically requires droplet synchronization.20 Another method is to directly inject chemicals into the droplets through a side channel, where the chemical solution to be added forms a pendant drop at the fluid junction and merges with the microdroplets passing by, building a pathway for mass transfer until the microdroplet leaves and disconnects from this side channel.21,22 Nightingale et al. reported an improvement to this method by introducing gas plugs, which help pinch off the pendent drop to achieve a more precise addition, and applied it to a multistep growth of CdSe quantum dots.23 To enable better control over the chemical infusion rate, Ismagilov et al. fabricated an array of side channels that links to the main reaction channel; the chemical can be added to microdroplets through these side channels at multiple locations along their flow path, and therefore at discrete time points.24 Later studies also showed strategies for improving pressure control as required in this method.25 However, the device still requires fabrication with polydimethylsiloxane (PDMS), which limits its application for many organic solvents.
To precisely control both the amount of chemical added and the rate of addition, we present here a novel method that enables highly controlled material transport into the microfluidic droplets. This work is based on a previous study in which we demonstrated that instead of using an inert oil as the carrier fluid, a W/O miniemulsion can be used as the continuous phase (CP) to generate aqueous monodisperse microdroplets as dispersed phase (DP).26 A W/O miniemulsion (or called nanoemulsion27) is a thermodynamically metastable system with aqueous nanodroplets, typically in the 50–500 nm range, dispersed in an immiscible organic solvent with a stabilizing surfactant. These nanodroplets serve as carriers for chemicals to be supplied to the aqueous microdroplets. In this paper, we enable the controllable transport of nanodroplets into the microdroplets via electrocoalescence, mediated by an applied alternating electrical field. Although electrocolaescence has been applied extensively in droplet microfluidic systems to merge droplets of similar sizes,7 we report here for the first time its application to merge droplets of dramatically different sizes, which requires a much higher electrical field strength (∼10
000 V cm−1) to realize. Since the nanodroplets are three orders of magnitude smaller than the microdroplets, we can easily achieve a nanodroplet-to-microdroplet population ratio of greater than 1 million. Such a large population ratio makes it possible to control the chemical addition rate over a wide range, and the addition is in fact “quasi-continuous”, while a paired fusion of microdroplets of similar sizes or infusion through side channels can only achieve a “one-shot” or “multi-shot” discrete injection. Two different types of signal – continuous and cyclic burst square waves – are used. For the continuous signal, the addition rate is controlled by the field strength while the total quantity of chemical added is controlled by the residence time of the microdroplet under the electrical field. For cyclic burst square wave signals, the control over chemical addition is achieved with the number of cycles and the field-on/field-off periods in each cycle. In this method, the electrical field is applied between two parallel electrodes placed outside the flow channel without any physical contact with the droplet flow. Such non-invasiveness makes it easily implementable to many existing systems, such as in the much-explored area of flow synthesis of nanomaterials.15,18,28–32 Here, we demonstrate an application in the controlled nucleation and growth of gold nanoparticles. The addition of the reducing reagent, L-ascorbic acid, is well-controlled to produce gold nanoparticles of tunable size and size distributions. Finally, this method is extended to achieve a three-step continuous synthesis of Au–Pd core–shell nanoparticles with multiple additions of miniemulsions carrying different chemical species. Our work enhances the operational flexibility of microdroplets as reaction vessels, enabling precisely controlled mass transfer at the sub-millimeter scale.
This paper is organized as follows. First we discuss how chemical addition rate and dosage are controlled by a continuous or cyclic burst electrical signal in section II and III. In section IV, we show how such control capabilities can be leveraged for gold nanoparticle synthesis. Section V demonstrates how the multi-step synthesis of gold–palladium core–shell nanoparticles can be easily achieved with this method. In the last section, an a priori simulation based on the electrostatic theory and Navier–Stokes equation is presented to reveal the physical principles involved in this process.
Continuous square wave signals were used for this experiment. Two different electrical field strengths, 4300 V cm−1 and 11
400 V cm−1, were applied to control the electrocoalescence rate, and their influence on the final gold nanoparticle size and size distribution was investigated. For the low field strength condition, the miniemulsion prepared with 0.17 wt% L-ascorbic acid solution was used. After a residence time of 600 seconds, 75% of the total nanodroplets were consumed (calculated from microdroplet size expansion). This was equivalent to adding 0.0075 mmol of L-ascorbic acid per 1 mL of the HAuCl4 solution (0.01 mmol gold). For the high field condition, we still used a residence time of 600 seconds for consistency, although electrocoalescence was completed in less than 60 seconds. To keep the total amount of added L-ascorbic acid the same, we used the miniemulsion prepared with a 0.13 wt% (75% of 0.17 wt%) L-ascorbic acid solution instead, which was then fully consumed to achieve the same final reaction stoichiometry.
To achieve variable residence times at fixed flow rates of the CP and DP, we fixed the applied electrical field strength at 11400 V cm−1, and varied the tube length under the electric field to achieve residence times of 3, 8, and 60 seconds. A miniemulsion prepared with a higher concentration 0.34 wt% L-ascorbic acid solution was used here, and these three residence times yielded 25%, 50%, and 100% microdroplet size expansion. Equivalently, they resulted in 25%, 50% and 100% nanodroplet consumption respectively, giving final reaction stoichiometries of 0.005, 0.01, and 0.02 mmol L-ascorbic acid per 1 mL of the HAuCl4 solution (containing 0.01 mmol gold).
After electrocoalescence, the reaction mixture was collected in a glass vial. The pink aqueous phase was decanted, washed with ethanol 3 times to remove all of the surfactant, and re-dispersed in deionized water. Due to aggregation, dynamic light scattering (NanoBrook 90Plus PALS) characterization did not give the primary particle information. Thus, UV-vis spectroscopy (Agilent Cary 60) and transmission electron microscopy (TEM, FEI Tecnai F20) were used to characterize the size and size distribution of the primary gold nanoparticles.
The 3-step synthesis of Au–Pd nanoparticles is shown in Fig. 1(b). HAuCl4 aqueous droplets were first generated by the ascorbic acid miniemulsion at the T-junction, with both syringe pumps operated at 30 μL min−1. The microdroplets and the miniemulsion were passed into an 11
400 V cm−1 electrical field for electrocoalescence. Although electrocoalescence was fast at such a high field strength, we allowed an additional two-minute residence time for complete gold nanoparticle growth. Following this crystal growth process, the Pd-precursor emulsion was infused into the system at a rate of 30 μL min−1 through the PVDF Y-junction. This flow rate was low enough to allow the addition of nanodroplets of Na2PdCl4 without breakage of the elongated microdroplet plugs. The nanodroplets and microdroplets were routed under the same electrode couple for the transfer of all Na2PdCl4 into the microdroplet. After that, another stream of ascorbic acid miniemulsion was infused at a rate of 30 μL min−1 through another PVDF Y-junction, followed by another electrocoalescence process to reduce Na2PdCl4. The aqueous phase was then decanted and washed with ethanol three times to remove surfactants and dodecane residue, and re-dispersed in water for UV-Vis and TEM characterization.
000 V cm−1, further increase in field strength does not significantly change the size expansion rate, indicating that we have almost reached a limiting scenario, which will be discussed in detail later in section VI.
Besides addition rate control, dosage control could also be achieved by controlling the residence time. By adjusting the length of the FEP tubing wound under the electrodes, we were able to control the residence time for a droplet to pass through, at constant flow velocity. For example, under 7100 V cm−1, we could allow a 7 second residence time for the passage of the microdroplet for a 20% size expansion, or a 17 second residence time to obtain a 50% size expansion by looking up the required size expansion-time relationship in Fig. 2(b).
In summary, both addition rate and dosage can be tuned with a continuous square wave signal by changing the field strength and residence time. For real-time digital control, however, continuous signals only work for modulation of the addition rate because field strength can be adjusted digitally, while changing the residence time for the modulation of the dosage requires the rearrangement of the tube length or changing electrode size, and is thus not suitable for online control.
000 V cm−1 after amplification.
Still using microdroplet size expansion as a characterization of the rate and dosage of chemical addition, we measured the size expansion of microdroplets as a function of the number of cycles. In these experiments, the tubing length was fixed to give a residence time of 66 seconds. As shown in Fig. 3(b), an increase in the number of cycles or the field-on time in each cycle (correlated to number of periods in each cycle) could increase the dosage. Interestingly, when these two variables were multiplied to form a new variable - the total field-on time, all curves overlapped, and a clear relationship between dosage and total field-on time can be visualized as in Fig. 3(c). Because of the replenishment of nanodroplets near microdroplets during the field-off time, the total field-on time required for microdroplet size expansion was shorter under these conditions than that required for a continuous 10
000 V cm−1 square wave, which has been shown in Fig. 2(b).
Thus, with a fixed residence time, dosage can be tuned digitally by the control of the total field-on time using Fig. 3(c) as a reference without changes to the device geometry, which now provides real-time dosage control capability. In addition, we can tune the field-off time interval between bursts and maintain a fixed total field-on time, which correspondingly results in different total process times, to achieve different addition rates with a fixed dosage.
Under a continuous square wave signal, the field strength was controlled to achieve different rates of L-ascorbic acid addition, and the results follow the hypothesis above very well. By maintaining a reaction stoichiometry of 0.0075 mmol of ascorbic acid to 0.01 mmol of HAuCl4, and using a high field of 11
400 V cm−1 or a low field of 4300 V cm−1, we fused nanodroplets into the microdroplets to achieve fast or slow mass transfer respectively. In the low field case, a long residence time of 600 s was needed, while in the high field case, nanodroplets were completely consumed in less than 60 s (although we let the train of droplets flow in the tube for the entire 600 s for consistency). In both cases, the microdroplets turned pink, showing the formation of gold nanoparticles. A detailed characterization of the UV-vis spectrum of these two samples reveals key differences; as shown in Fig. 4(a), the slow addition sample gave a higher peak absorption wavelength at 560 nm compared to 535 nm for the fast addition sample. Previous research35 shows a positive correlation between the wavelength of the UV-vis absorption peak and the size of the primary particles. This indicates that the slow addition of the reducing agent resulted in a larger mean particle size, as discussed above. We confirmed this inference with TEM characterization, as shown in Fig. 4(c–f). The size of the gold nanoparticles prepared with a fast addition over 60 seconds was 44.6 ± 12.6 nm, while it was 81.2 ± 25.5 nm when L-ascorbic acid was added slowly over 600 seconds. Under an 11
400 V cm−1 continuous square wave signal, dosage of the L-ascorbic acid can be tuned by varying the residence time by changing the length of tubing between the two electrode plates. Under the same high field condition (11
400 V cm−1), we investigated three different residence times 3 s, 8 s, and 60 s (at the same flow speed) to control the added ascorbic acid quantity. We also raised the concentration of L-ascorbic acid to more substantially differentiate these three conditions, resulting in different reaction stoichiometries of 0.005, 0.01, and 0.02 mmol ascorbic acid per 0.01 mmol of HAuCl4. As shown in Fig. 4(b), an increase in the amount of added ascorbic acid results in a higher absorption peak wavelength, due to contributions from different size and shaped nanoparticles in the samples (TEM images and size distribution histograms in the ESI†).
As shown in Fig. 5(a), the prepared Au–Pd nanoparticles formed a black suspension when dispersed in water, which matches previous reports on Au–Pd nanoparticles prepared by batch methods.40Fig. 5(b) shows the UV-vis characterization of the Au nanoparticles collected before Na2PdCl4 addition, and the finally collected Au–Pd nanoparticles after three steps, both prepared with a gold concentration of 10−6 mol mL−1. Sánchez-Ramírez et al.41 and Hu et al.36 studied the difference in optical properties between Au–Pd bimetallic nanoparticles and a mixture of Pd and Au nanoparticles, showing that a complete disappearance of the surface plasmon resonance absorbance is an indication of effective coating. Ideally, the gold surface plasmon resonance peak should be completely eliminated if the coverage of Pd on Au is 100%, and a broad featureless absorption which contributes to the black color of the suspension, should be obtained.36 In our case, the residue peak indicated an insufficient coverage, although the absorbance had been significantly suppressed. The result is also in agreement with the report from Shim et al. that Pd nanoparticles do not have surface plasmon resonance absorption peaks.42 An increase in Na2PdCl4 concentration could result in a further increase in Pd coverage, but also result in more self-nucleated pure Pd nanoparticles.
TEM characterization of the Au–Pd nanoparticles gives a clearer view of the Au core–Pd shell structure. With different electron transmission rates, Pd shell and Au core can be easily distinguished for small particles, which had a shell thickness of about 10 nm, as shown in Fig. 5(d). The Pd-coating was non-uniform, indicating that it formed mostly through secondary nucleation on the gold surface. The morphology is also consistent with previous reports on two-step Au–Pd core–shell nanoparticle synthesis using batch methods.43 For larger particles, however, the relatively thin Pd shell was difficult to observe (Fig. 5(c)) that TEM-EDS (energy dispersive spectrum) was needed to visualize the surface coverage of palladium. Fig. 5(e) and (f) give evidence for effective palladium coating, and a clear Au core–Pd shell structure is seen. Defective edges without thick Pd coating were also observed, and probably contributed to the residual surface plasmon resonance peak in the UV-vis spectrum of the Au–Pd sample (Fig. 5(b)).
| b = 0.643(3Ca)2/3r1 | (1) |
For an elongated microdroplet transported in a circular channel, r1 is the radius of the hemispherical droplet end. Since r1 + b = rtube = 381 μm, we calculate r1 = 368 μm and the lubrication layer had a thickness of only 13 μm. The microdroplet length and drop–drop distance can be calculated from the flow rates and the droplet generation frequency.
With the geometrical parameters determined above, nanodroplet motion can be simulated by numerically solving the Navier–Stokes equation. Under the 50 μL min−1 (CP)/10 μL min−1 (DP) flow condition, for droplets with a diameter smaller than 60 μm, the fluid drag force (calculated based on a characteristic velocity of 1 mm s−1) is stronger than the gravitational force, as calculated from the Stokes law. Thus for nanodroplets, gravitational settling of the nanodroplets can be neglected and the droplets are entrained by dodecane. Treating the miniemulsion as a pseudo-single phase, we determined the nominal viscosity (measured with a Ubbelohde viscometer) and density to be 4 cP and 800 kg m−3, respectively. Thus, the laminar flow field could be calculated numerically, with results shown in Fig. 6(a). Note that we fix the coordinate on the microdroplet, effectively making the wall move in the opposite direction at a speed of 2.2 mm s−1.
![]() | (2) |
![]() | (3) |
From the electrostatic force expressions, it can be inferred that the force is extremely sensitive to drop–drop distance and the nanodroplet size. Thus, a theoretical analysis is complicated by the fact that nanodroplet size and size distribution are not time-invariant since, under a high electrical field, nanodroplet self-coalescence is significant. Self-coalescence of nanodroplets under an alternating electrical field has been reported previously.46 To verify self-coalescence, when we flowed only the miniemulsion continuous phase through this field, the nanodroplets easily self-coalesced into naked eye-visible large droplets in 5–10 seconds. When microdroplets were present, we did not observe the formation of these large droplets due to the continuous consumption of the nanodroplets by the microdroplets, preventing them from growing too large. To in situ characterize the dynamic time-evolution of nanodroplet size and size distribution accurately was beyond our experimental capability, while analysis performed after the collection of self-coalesced droplets was inaccurate due to their easier subsequent coalescence in the collection well. Thus, self-coalescence complicates the theoretical prediction based solely on the electrostatic force expressions, as the actual miniemulsion droplet sizes are unknown and time-dependent.
Details of the assumptions made and the numerical simulation set-up in a commercial finite element software, COMSOL Multiphysics®, are presented in the ESI.† The a priori simulated nanodroplet consumption rate under these physically limiting conditions is shown by the solid line in Fig. 7(a). Note that nanodroplet depletion leads to microdroplet expansion and, at a flow rate ratio of 5 (CP)
:
1 (DP), corresponds to microdroplet size expansion on a one-to-one basis. For example, a 50% nanodroplet depletion corresponds to a 50% microdroplet size expansion. The simulated curve is the theoretical convection-controlled limit of nanoemulsion consumption rate under a continuous signal. Recall that in section 3.2 and Fig. 2(b) we showed that the mass transport rate was barely enhanced when the field strength was sufficiently strong, which corresponds to this theoretical limit. We obtained additional microdroplet size expansion data points for the three highest field strengths, 8600, 10
000 and 11
400 V cm−1, and deduced the rate of nanodroplet depletion based on the one-to-one relationship mentioned before. As shown in Fig. 7(a), once the field was sufficiently high (∼10
000 V cm−1), further increase in the field strength did not contribute to a faster consumption rate, indicating we had approached the convection-controlled limit. The close proximity between the predicted curve and experimental results also strengthens this conclusion and provides a good validation of the numerical model. We further simulated such convection-controlled cases for higher flow rates at a fixed DP-to-CP flow rate ratio of 1
:
5. The droplet flow pattern, including droplet sizes and drop–drop distances, are different under different flow conditions. A higher total flow rate decreases droplet size and drop-drop distance using the same T-junction, which were measured experimentally and used to determine the simulation geometries for each case (see details in ESI†). These additional simulation results are plotted in Fig. 7(b), showing that increasing flow rates greatly enhanced nanodroplet depletion rate by a stronger convection and a shorter drop-drop distance.
000 V cm−1, the electrocoalescence is so fast that the entire process is controlled by the convection of nanodroplets to the microdroplet surfaces as they are dragged by the circulating continuous phase flow, in agreement with theoretical predictions of the process rate under this convection-controlled scenario. This method was applied to a single step synthesis of gold nanoparticles and a multistep flow synthesis of Au–Pd core–shell nanoparticles, with the size and size distribution of the resulting nanoparticles well controlled. Our approach is advantageous for the post-addition of chemicals to microfluidic droplets since there is no requirement for droplet synchronization, as in the droplet pair-fusion method, or surface wettability control, as in the direct injection method, or complicated channel fabrication, as in the arrayed side-channel injection method. The electrical field is applied by external electrodes with no physical contact with the reagents, making this method well suited to many existing aqueous flow synthesis systems.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c8lc00114f |
| This journal is © The Royal Society of Chemistry 2018 |