# Continuous Flow Controlled Synthesis of Gold Nanoparticles Using Pulsed Mixing Microfluidic System

Authors: Guojun Liu; Xuhao Yang; Yan Li; Zhigang Yang; Wen Hong; JianFang Liu
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

To prepare the gold nanoparticles (AuNPs) with uniform sizes, fine morphology, and good monodispersity, a pulsed mixing microfluidic system based on PZT actuation was presented. The system includes PZT micropump and Y type micromixer. By adjusting voltage (entrance flow rate), pulsed frequency, phase, and other parameters, a variety of mixing modes can be achieved, so as to realize the controllable synthesis of nanoparticles in a certain range. By numerical simulation and analysis, the channel section size, entrance angle, and pulse frequency were optimized. Based on the optimized structure and working parameters, the test prototype has been manufactured in lab, and the related synthesis tests of AuNPs were carried out. The test results indicate that AuNPs with uniform morphology and good monodispersity can be synthesized using the system with the section size (0.4 mm × 0.4 mm), the entrance channel angle (60°) under condition of the pulsed frequency (300 Hz), and the entrance flow rate (4 mL/min). The average diameter and its standard deviation of AuNPs synthesized were 21.6 nm, 4.83 nm, respectively. The research work above can be applied to the fields such as the controlled synthesis of noble metal nanoparticles, biomedicine, and microchemical system.


---

## Body

## 1. Introduction

Gold nanoparticles (AuNPs) have received tremendous growing attention during the recent few years because of their unique characteristics [1–4] and have already found many kinds of applications in catalysis [5–7], sensing [8–10], diagnosis [11–13], biomedicine [1, 14, 15], drug delivery [16–18], and so forth. In addition, it also has been applied to a chemical analysis with Localized Surface Plasmon (LSP) [19] and Surface Enhanced Raman Spectroscopy (SERS) [20, 21]; both the average size of the gold nanoparticles and the size distribution are critically important factors in such applications [22].AuNPs which have been synthesized by various physical and chemical methods mainly include UV irradiation, laser ablation, plasma synthesis, citrate synthesis, turkevich method, wet chemical synthesis, and chemical reduction method [1]. The shapes of obtained AuNPs are multifarious, such as triangles [23, 24], nanorods [25], and stars [26]. The most of AuNPs synthesis process was achieved by reducing the gold salt in the water phase [27]. On the basis of the method, microfluidic mixer has been used for the synthesis of gold nanoparticles by scholars. The continuous flow methods using micromixers have been proven to be successful alternative methods with a possibility to synthesize large quantities with better control on the distribution, shape, and size of the particles [28, 29]. However, some of reports still cannot realize a narrow distribution of particle size. The reason for that is the lack of rapid mixability at a laminar flow condition with a microchannel.Mixing is a basic unit operation for chemical reactions and micromixing can extremely influence the selectivity, yield, and quality of final products in various chemical processes [30–32]. In this paper, a pulsed mixing microfluidic system for the synthesis of AuNPs is developed. Firstly, the mechanism of gold nanoparticles formation and pulsed micromixing is introduced. Secondly, the structure and working parameters of Y type micromixer are optimized after simulation analysis. Finally, the gold nanoparticles are prepared and the test results are discussed.


## 2. Mechanism of Controlled Synthesis of AuNPs

### 2.1. Mechanism of AuNPs Formation

Based on chemical reduction mechanism, the liquid-phase synthesis method of AuNPs includes three processes, namely, precipitation of gold atoms, nucleation, and growth of crystal nuclei. The synthesis mechanism can be described by Figure1 as follows. In the initial stage of mixing of the reducing agent (such as sodium citrate) and precursor (such as chloroauric acid), as the precursor is reduced, the concentration of gold atoms continues increasing (stage I: precipitation of gold atoms). When the concentration of gold atoms exceeds the critical supersaturation (stage II: nucleation), the gold atoms begin to gather and form crystal nuclei. Also, with the formation and growth of crystal nuclei, the concentration of gold atoms begins to decrease. After a large amount of gold atoms are consumed, the concentration of gold atoms gradually declines down below the critical supersaturation (stage III: growth of crystal nuclei), the number of crystal nuclei no longer increases, and the growth of nuclei will dominate the whole process. When the concentration of gold atoms decreases to the saturation level, crystal nuclei stop growing and the AuNPs are finally prepared. As can be seen from the above reaction processes, stages I and II are particularly important, whose mixing effect largely affects the amount and size uniformity of AuNPs prepared by the liquid-phase synthesis method. To obtain the AuNPs with uniform morphology and small size deviation, besides the appropriate adjustment of ingredients proportion, concentration, and so forth of two-phase agents, it is crucial to regulate the mixed reaction efficiency. Only when the mixing effect can be controlled freely, will it become possible to realize the controlled synthesis of AuNPs. So the controllability of microfluidic reaction process is an important condition to achieve the AuNPs with uniform particle size, regular morphology.Figure 1
Schematic illustration of synthesis mechanism of gold nanoparticles.


### 2.2. Mechanism of Pulsed Micromixing

Pulsed micromixing is an active mixing method with high mixing efficiency [22, 33] and its working mechanism is shown in Figure 2. Pulsed micromixing system includes two functional parts, namely, driving source (PZT micropump) and Y type micromixer. Two PZT micropumps with same structure and working performance provide the mixing power for Y type micromixer. And two different kinds of fluids are pumped into the mixer with the changing flow rate, which fluctuates from high to low periodically by different frequency and phase. Therefore, the pulsing fluid fields will cause rapid stretching and folding of the fluid interfaces and then realize a variety of mixed modes, such as the laminar mixing, chaotic advection, and turbulent mixing. As illustrated in Figure 2, when two sine wave input signals have a phase difference of 180 degrees, solution A (such as chloroauric acid) and solution B (such as sodium citrate) achieve pulsating crossover type input in the micromixer. Thus, it increases the contact area and the contact opportunity of two kinds of liquid, shortens the mixing time, and improves the efficiency of mixing.Figure 2
Schematic diagram of pulsed micromixing.


## 2.1. Mechanism of AuNPs Formation

Based on chemical reduction mechanism, the liquid-phase synthesis method of AuNPs includes three processes, namely, precipitation of gold atoms, nucleation, and growth of crystal nuclei. The synthesis mechanism can be described by Figure1 as follows. In the initial stage of mixing of the reducing agent (such as sodium citrate) and precursor (such as chloroauric acid), as the precursor is reduced, the concentration of gold atoms continues increasing (stage I: precipitation of gold atoms). When the concentration of gold atoms exceeds the critical supersaturation (stage II: nucleation), the gold atoms begin to gather and form crystal nuclei. Also, with the formation and growth of crystal nuclei, the concentration of gold atoms begins to decrease. After a large amount of gold atoms are consumed, the concentration of gold atoms gradually declines down below the critical supersaturation (stage III: growth of crystal nuclei), the number of crystal nuclei no longer increases, and the growth of nuclei will dominate the whole process. When the concentration of gold atoms decreases to the saturation level, crystal nuclei stop growing and the AuNPs are finally prepared. As can be seen from the above reaction processes, stages I and II are particularly important, whose mixing effect largely affects the amount and size uniformity of AuNPs prepared by the liquid-phase synthesis method. To obtain the AuNPs with uniform morphology and small size deviation, besides the appropriate adjustment of ingredients proportion, concentration, and so forth of two-phase agents, it is crucial to regulate the mixed reaction efficiency. Only when the mixing effect can be controlled freely, will it become possible to realize the controlled synthesis of AuNPs. So the controllability of microfluidic reaction process is an important condition to achieve the AuNPs with uniform particle size, regular morphology.Figure 1
Schematic illustration of synthesis mechanism of gold nanoparticles.


## 2.2. Mechanism of Pulsed Micromixing

Pulsed micromixing is an active mixing method with high mixing efficiency [22, 33] and its working mechanism is shown in Figure 2. Pulsed micromixing system includes two functional parts, namely, driving source (PZT micropump) and Y type micromixer. Two PZT micropumps with same structure and working performance provide the mixing power for Y type micromixer. And two different kinds of fluids are pumped into the mixer with the changing flow rate, which fluctuates from high to low periodically by different frequency and phase. Therefore, the pulsing fluid fields will cause rapid stretching and folding of the fluid interfaces and then realize a variety of mixed modes, such as the laminar mixing, chaotic advection, and turbulent mixing. As illustrated in Figure 2, when two sine wave input signals have a phase difference of 180 degrees, solution A (such as chloroauric acid) and solution B (such as sodium citrate) achieve pulsating crossover type input in the micromixer. Thus, it increases the contact area and the contact opportunity of two kinds of liquid, shortens the mixing time, and improves the efficiency of mixing.Figure 2
Schematic diagram of pulsed micromixing.


## 3. The Structure and Working Parameters Optimization of Y Type Micromixer

### 3.1. Modeling of Y Type Micromixer

To optimize the structure of mixer, Y type micromixer model is established as shown in Figure3. Through the channel entrances 1 and 2, the two kinds of fluid are, respectively, pumped into the micromixer, complete mixing, or reaction in the mainstream channel of micromixer and then are output from the outlet. The main structure dimensions of mixer include entrance channel length S, angle between two flow entrances θ, length of main mixing channel L, width of the channel W, and depth of flow channel H. The origin of coordinate system is set in the center of the junction of two entrance channels, with the X direction along the mainstream channel.Figure 3
Simulation model of Y type micromixer.


### 3.2. Evaluation of Mixing Efficiency

To evaluate the mixing efficiency quantitatively, an important index parameter, namely, the degree of mixing (σ), is defined [34]:(1)σ=1.0-1c∞∑i=1Nci-c∞2N,where N is the number of sampling points inside the cross section, ci is the mass fraction at sampling point i, and c∞ is the optimal mixing mass fraction. The value of σ should be initially zero and tend monotonically towards a large value during the mixing progresses. Larger σ indicates better mixing efficiency. When the mixing degree σ=0, there is no mixture between two phases of fluid. When the mixing degree σ=1, two phases of fluid are in a real state of complete mixing.


### 3.3. Optimization of the Structure Parameters

#### 3.3.1. Analysis of Effects of Section Size on the Mixing Effect

In experimental test, the flow rate of PZT micropump is fixed on a specified value (in this work, the flow rate of PZT micropump is 4 mL/min). At the entrance of mixer, the velocity amplitudes of two phase flow vary with size of section. In simulation analysis, four sets of structure parameters are selected as shown in Table1, and the angle between the entrance channels is initially selected 90°. According to the different sizes of cross section, the four volume concentration distributions in the flow channel are analyzed by computational fluid dynamics (CFD) simulations, as shown in Figure 4. Analyzing the concentration distribution diagram of Figure 4(a), it can be clearly seen that, after one and a half cycles of pulsed mixing, the mixing effect of mixed solution is not obvious, and the boundary of two solution is still clear. The mixing effect has been improved a little in Figure 4(b), after one and a half pulsed mixing cycles. As for the Figure 4(c), the mixing is almost completed. In Figure 4(d), the field near the two sides of channel wall still has clear color boundary. Observing the concentration distribution diagram, it can be concluded that the mixing effect of Figure 4(c) is optimal. The simulation results indicate that when the entrance flow rate and working frequency are ascertained, there exist corresponding optimum values of the section size in order to obtain the best mixing effect. When the output volume of PZT micropump per pulse can be filled with the whole cross section, the microsolution presents a bullet type distribution in the flow channel. And the shorter the length of bullet is, the larger the contact area of two solutions is. Of course, the mixing effect will be better. While the output volume per pulse just can be filled with microchannel exactly, the mixing effect is the best. If it fails to fill the entire cross section, as shown in Figure 4(d), the mixing state belongs to the chaotic convection mode, and the mixing effect decreases compared with Figure 4(c). So when the flow rate and working frequency of micropumps are fixed, the size of cross section has an optimum value. Accordingly, the mixing degrees versus different section size and different section location are shown in Table 2. Based on the simulation results, comprehensively considering the fabrication process, structure, and saving the dosage of chemical reagent and other factors, the final section size after optimization is determined as 0.4 mm × 0.4 mm.Table 1
Section size of microchannel and the corresponding velocity amplitude of entrance.
W
×
H (mm)
0.2 × 0.2
0.3 × 0.3
0.4 × 0.4
0.5 × 0.5
Velocity amplitudeV (m/s)
1.667
0.741
0.417
0.267Table 2
Mixing degree of the fluid for different section size and location.
W
×
H (mm)
Mixing degreeσ
X
=
9
X
=
10
X
=
11
0.3 × 0.3
0.6341
0.6753
0.7892
0.4 × 0.4
0.8446
0.8761
0.9189
0.5 × 0.5
0.8537
0.8928
0.9187Figure 4
Concentration distributions of flow for different cross sections.
(a)
(b)
(c)
(d)


#### 3.3.2. Analysis of Effects of Different Entrance Channel Angle on the Mixing Effect

The angle of entrance channel is another important structure parameter influencing the mixing effect. In this work, the effects of different entrance channel angle on the mixing effect are stimulated by CFD software, when section size is determined as 0.4 mm × 0.4 mm with other conditions unchanged. The simulation results of volume concentration distribution in the mixing channel are shown in Figure5. The volume concentration distributes in bullet shape, and the mixing effect by layer morphology is better than other angles due to more contact and diffusion areas in boundary. Table 3 lists the mixing degree value of mainstream channel with different entrance channel angles and different section locations. The simulation results in Table 3 show that the mixing degree increases first and then decreases with the increase of angle between entrance channels. When the angle equals 60°, the maximum value of mixing degree is achieved, and three mixing degree values are all greater than 0.9. Therefore, 60° is finally selected as the optimal angle of entrance channels.Table 3
Mixing degree of the fluid phase 1 for different flow angles and locations.
X (mm)
Mixing degreeσ
θ
=
3
0
°
θ
=
6
0
°
θ
=
9
0
°
θ
=
12
0
°
θ
=
15
0
°
θ
=
18
0
°
X
=
9
0.8367
0.9104
0.8446
0.8118
0.7048
0.6610
X
=
10
0.8557
0.9397
0.8761
0.8371
0.7449
0.6922
X
=
11
0.8429
0.9084
0.9189
0.7898
0.6761
0.6261Figure 5
Concentration distributions of flow for different flow angle.


### 3.4. Optimization of the Working Parameter

For the pulsed micromixing system, the working frequency and entrance flow rate are the most important working parameters which directly affect the mixing effect. Since the testing system adopts the mature driving technology of PZT micropump, the optimization of working parameters focuses on the research of working frequency. The specific frequency in the range of 50~350 Hz is studied and analyzed for its influence on the mixing effect, assuming that the entrance flow rate is 4 mL/min. The simulation results of CFD and the relationship curves between the mixing degree and working frequency are shown in Figure6. Analyzing the distribution graphs and curves in Figure 6, it can be concluded that when the flow rate remains constant (4 mL/min), the mixing degree increases first and then decreases with the increase of working frequency. When the frequency is 300 Hz, the mixing degree reaches the maximum value at 0.97, which realizes a nearly complete mixture. Under this condition, the mixing effect is optimal.Figure 6
Concentration distributions of flow for different working frequencies and relationship curves between working frequency and mixing degree in different sections of microchannel.
(a)
(b)


## 3.1. Modeling of Y Type Micromixer

To optimize the structure of mixer, Y type micromixer model is established as shown in Figure3. Through the channel entrances 1 and 2, the two kinds of fluid are, respectively, pumped into the micromixer, complete mixing, or reaction in the mainstream channel of micromixer and then are output from the outlet. The main structure dimensions of mixer include entrance channel length S, angle between two flow entrances θ, length of main mixing channel L, width of the channel W, and depth of flow channel H. The origin of coordinate system is set in the center of the junction of two entrance channels, with the X direction along the mainstream channel.Figure 3
Simulation model of Y type micromixer.


## 3.2. Evaluation of Mixing Efficiency

To evaluate the mixing efficiency quantitatively, an important index parameter, namely, the degree of mixing (σ), is defined [34]:(1)σ=1.0-1c∞∑i=1Nci-c∞2N,where N is the number of sampling points inside the cross section, ci is the mass fraction at sampling point i, and c∞ is the optimal mixing mass fraction. The value of σ should be initially zero and tend monotonically towards a large value during the mixing progresses. Larger σ indicates better mixing efficiency. When the mixing degree σ=0, there is no mixture between two phases of fluid. When the mixing degree σ=1, two phases of fluid are in a real state of complete mixing.


## 3.3. Optimization of the Structure Parameters

### 3.3.1. Analysis of Effects of Section Size on the Mixing Effect

In experimental test, the flow rate of PZT micropump is fixed on a specified value (in this work, the flow rate of PZT micropump is 4 mL/min). At the entrance of mixer, the velocity amplitudes of two phase flow vary with size of section. In simulation analysis, four sets of structure parameters are selected as shown in Table1, and the angle between the entrance channels is initially selected 90°. According to the different sizes of cross section, the four volume concentration distributions in the flow channel are analyzed by computational fluid dynamics (CFD) simulations, as shown in Figure 4. Analyzing the concentration distribution diagram of Figure 4(a), it can be clearly seen that, after one and a half cycles of pulsed mixing, the mixing effect of mixed solution is not obvious, and the boundary of two solution is still clear. The mixing effect has been improved a little in Figure 4(b), after one and a half pulsed mixing cycles. As for the Figure 4(c), the mixing is almost completed. In Figure 4(d), the field near the two sides of channel wall still has clear color boundary. Observing the concentration distribution diagram, it can be concluded that the mixing effect of Figure 4(c) is optimal. The simulation results indicate that when the entrance flow rate and working frequency are ascertained, there exist corresponding optimum values of the section size in order to obtain the best mixing effect. When the output volume of PZT micropump per pulse can be filled with the whole cross section, the microsolution presents a bullet type distribution in the flow channel. And the shorter the length of bullet is, the larger the contact area of two solutions is. Of course, the mixing effect will be better. While the output volume per pulse just can be filled with microchannel exactly, the mixing effect is the best. If it fails to fill the entire cross section, as shown in Figure 4(d), the mixing state belongs to the chaotic convection mode, and the mixing effect decreases compared with Figure 4(c). So when the flow rate and working frequency of micropumps are fixed, the size of cross section has an optimum value. Accordingly, the mixing degrees versus different section size and different section location are shown in Table 2. Based on the simulation results, comprehensively considering the fabrication process, structure, and saving the dosage of chemical reagent and other factors, the final section size after optimization is determined as 0.4 mm × 0.4 mm.Table 1
Section size of microchannel and the corresponding velocity amplitude of entrance.
W
×
H (mm)
0.2 × 0.2
0.3 × 0.3
0.4 × 0.4
0.5 × 0.5
Velocity amplitudeV (m/s)
1.667
0.741
0.417
0.267Table 2
Mixing degree of the fluid for different section size and location.
W
×
H (mm)
Mixing degreeσ
X
=
9
X
=
10
X
=
11
0.3 × 0.3
0.6341
0.6753
0.7892
0.4 × 0.4
0.8446
0.8761
0.9189
0.5 × 0.5
0.8537
0.8928
0.9187Figure 4
Concentration distributions of flow for different cross sections.
(a)
(b)
(c)
(d)


### 3.3.2. Analysis of Effects of Different Entrance Channel Angle on the Mixing Effect

The angle of entrance channel is another important structure parameter influencing the mixing effect. In this work, the effects of different entrance channel angle on the mixing effect are stimulated by CFD software, when section size is determined as 0.4 mm × 0.4 mm with other conditions unchanged. The simulation results of volume concentration distribution in the mixing channel are shown in Figure5. The volume concentration distributes in bullet shape, and the mixing effect by layer morphology is better than other angles due to more contact and diffusion areas in boundary. Table 3 lists the mixing degree value of mainstream channel with different entrance channel angles and different section locations. The simulation results in Table 3 show that the mixing degree increases first and then decreases with the increase of angle between entrance channels. When the angle equals 60°, the maximum value of mixing degree is achieved, and three mixing degree values are all greater than 0.9. Therefore, 60° is finally selected as the optimal angle of entrance channels.Table 3
Mixing degree of the fluid phase 1 for different flow angles and locations.
X (mm)
Mixing degreeσ
θ
=
3
0
°
θ
=
6
0
°
θ
=
9
0
°
θ
=
12
0
°
θ
=
15
0
°
θ
=
18
0
°
X
=
9
0.8367
0.9104
0.8446
0.8118
0.7048
0.6610
X
=
10
0.8557
0.9397
0.8761
0.8371
0.7449
0.6922
X
=
11
0.8429
0.9084
0.9189
0.7898
0.6761
0.6261Figure 5
Concentration distributions of flow for different flow angle.


## 3.3.1. Analysis of Effects of Section Size on the Mixing Effect

In experimental test, the flow rate of PZT micropump is fixed on a specified value (in this work, the flow rate of PZT micropump is 4 mL/min). At the entrance of mixer, the velocity amplitudes of two phase flow vary with size of section. In simulation analysis, four sets of structure parameters are selected as shown in Table1, and the angle between the entrance channels is initially selected 90°. According to the different sizes of cross section, the four volume concentration distributions in the flow channel are analyzed by computational fluid dynamics (CFD) simulations, as shown in Figure 4. Analyzing the concentration distribution diagram of Figure 4(a), it can be clearly seen that, after one and a half cycles of pulsed mixing, the mixing effect of mixed solution is not obvious, and the boundary of two solution is still clear. The mixing effect has been improved a little in Figure 4(b), after one and a half pulsed mixing cycles. As for the Figure 4(c), the mixing is almost completed. In Figure 4(d), the field near the two sides of channel wall still has clear color boundary. Observing the concentration distribution diagram, it can be concluded that the mixing effect of Figure 4(c) is optimal. The simulation results indicate that when the entrance flow rate and working frequency are ascertained, there exist corresponding optimum values of the section size in order to obtain the best mixing effect. When the output volume of PZT micropump per pulse can be filled with the whole cross section, the microsolution presents a bullet type distribution in the flow channel. And the shorter the length of bullet is, the larger the contact area of two solutions is. Of course, the mixing effect will be better. While the output volume per pulse just can be filled with microchannel exactly, the mixing effect is the best. If it fails to fill the entire cross section, as shown in Figure 4(d), the mixing state belongs to the chaotic convection mode, and the mixing effect decreases compared with Figure 4(c). So when the flow rate and working frequency of micropumps are fixed, the size of cross section has an optimum value. Accordingly, the mixing degrees versus different section size and different section location are shown in Table 2. Based on the simulation results, comprehensively considering the fabrication process, structure, and saving the dosage of chemical reagent and other factors, the final section size after optimization is determined as 0.4 mm × 0.4 mm.Table 1
Section size of microchannel and the corresponding velocity amplitude of entrance.
W
×
H (mm)
0.2 × 0.2
0.3 × 0.3
0.4 × 0.4
0.5 × 0.5
Velocity amplitudeV (m/s)
1.667
0.741
0.417
0.267Table 2
Mixing degree of the fluid for different section size and location.
W
×
H (mm)
Mixing degreeσ
X
=
9
X
=
10
X
=
11
0.3 × 0.3
0.6341
0.6753
0.7892
0.4 × 0.4
0.8446
0.8761
0.9189
0.5 × 0.5
0.8537
0.8928
0.9187Figure 4
Concentration distributions of flow for different cross sections.
(a)
(b)
(c)
(d)


## 3.3.2. Analysis of Effects of Different Entrance Channel Angle on the Mixing Effect

The angle of entrance channel is another important structure parameter influencing the mixing effect. In this work, the effects of different entrance channel angle on the mixing effect are stimulated by CFD software, when section size is determined as 0.4 mm × 0.4 mm with other conditions unchanged. The simulation results of volume concentration distribution in the mixing channel are shown in Figure5. The volume concentration distributes in bullet shape, and the mixing effect by layer morphology is better than other angles due to more contact and diffusion areas in boundary. Table 3 lists the mixing degree value of mainstream channel with different entrance channel angles and different section locations. The simulation results in Table 3 show that the mixing degree increases first and then decreases with the increase of angle between entrance channels. When the angle equals 60°, the maximum value of mixing degree is achieved, and three mixing degree values are all greater than 0.9. Therefore, 60° is finally selected as the optimal angle of entrance channels.Table 3
Mixing degree of the fluid phase 1 for different flow angles and locations.
X (mm)
Mixing degreeσ
θ
=
3
0
°
θ
=
6
0
°
θ
=
9
0
°
θ
=
12
0
°
θ
=
15
0
°
θ
=
18
0
°
X
=
9
0.8367
0.9104
0.8446
0.8118
0.7048
0.6610
X
=
10
0.8557
0.9397
0.8761
0.8371
0.7449
0.6922
X
=
11
0.8429
0.9084
0.9189
0.7898
0.6761
0.6261Figure 5
Concentration distributions of flow for different flow angle.


## 3.4. Optimization of the Working Parameter

For the pulsed micromixing system, the working frequency and entrance flow rate are the most important working parameters which directly affect the mixing effect. Since the testing system adopts the mature driving technology of PZT micropump, the optimization of working parameters focuses on the research of working frequency. The specific frequency in the range of 50~350 Hz is studied and analyzed for its influence on the mixing effect, assuming that the entrance flow rate is 4 mL/min. The simulation results of CFD and the relationship curves between the mixing degree and working frequency are shown in Figure6. Analyzing the distribution graphs and curves in Figure 6, it can be concluded that when the flow rate remains constant (4 mL/min), the mixing degree increases first and then decreases with the increase of working frequency. When the frequency is 300 Hz, the mixing degree reaches the maximum value at 0.97, which realizes a nearly complete mixture. Under this condition, the mixing effect is optimal.Figure 6
Concentration distributions of flow for different working frequencies and relationship curves between working frequency and mixing degree in different sections of microchannel.
(a)
(b)


## 4. Experiment

### 4.1. Reagents

HAuCl4·4H2O (precursor), Na3C6H5O7 (reducing agent), and CTAB (cationic surface active agent) were purchased from Sinopharm Chemical Reagent Corporation (Shanghai, China). The other reagents were of the highest grades available commercially. The water used in the experiments was of Milli-Q grade.


### 4.2. Apparatus

UV-Vis absorption spectrophotometer (UV-2501 PC, SHIMADZU, Japan) was used to record the characteristic absorption spectrum of gold colloid. The size and morphology of synthesized AuNPs were characterized by high resolution transmission electron microscopy (TEM TECNAI-12, Philips, Holland). The PZT micropump was fabricated using mature MEMS process in lab [35]. The output flow rate of micropump can be adjusted in the range of 0~6 mL/min, and the working frequency of micropump can be adjusted in the range of 10~400 Hz. The fabrication of Y type micromixer used a computer-numerically controlled (CNC) machining process and a polydimethylsiloxane (PDMS) replication process [4]. The mixing channel and reaction channel were packaged using PDMS and glass. Also, the design and manufacture of micromixer all adopted the optimized structural parameters above. The constant temperature heating platform (BHW-05A shanghai, China) was used to provide the stable reaction temperature. In addition, precision electronic balance and power supply controller were used in experiment test.Figure7 shows the scheme of the experimental setup, the prototype photo of PZT micropumps, and testing photo of controlled synthesis of gold nanoparticles using pulsed mixing microfluidic system. The detailed experimental procedure is as follows. First, all glassware used was clean in freshly prepared aqua regia solution (HCI : HNO3, 3 : 1) and rinsed with deionized water three times before used. Before the synthesis of AuNPs, the chambers of PZT micropump and the channels of micromixer were cleaned by deionized water. The gold precursor, reducing agent, and surfactant were, respectively, dissolved in deionized water to prepare the reaction reagents of the corresponding concentration [36], which include 100 mL HAuCl4 of 1 mM, 100 mL Na3C6H5O7 of 1.5 mM, and 100 mL CTAB of 200 mM. The temperature control platform was adjusted to keep the constant temperature of 100°C. With adjustments of the working voltage, the flow rate of PZT micropumps was fixed at 4 mL/min. To achieve the crossing pulsing micromixing of two reaction reagents, the phase of input signal of PZT micropumps was regulated to maintain 180°. After turning on the power of PZT micropumps, the gold precursor (HAuCl4) and the reducing agent (Na3C6H5O7) were pumped into the micromixer and growth channel by the crossing pulsing state. At the outlet of microchannel, the solution which experienced the “pulsed mixing” and passed through the “growth flow channel” was collected. At the same time, in order to avoid the agglomeration of AuNPs, the corresponding amount of CTAB was also added to the collected samples one time in proportion (nCTAB/nHAuCl4 = 200 : 1). After about ten minutes, the synthesized AuNPs reached stable state, and then the size and morphology of AuNPs can be characterized using the relevant testing apparatus. Note that the experiments were conducted in the frequency of 100, 200, and 300 Hz, respectively.Figure 7
(a) Experimental scheme; (b) the prototype photo of PZT micropumps; (c) testing photo of controlled synthesis of gold nanoparticles using pulsed mixing microfluidic system.
(a)
(b)
(c)


### 4.3. Measurement

To monitor the formation of gold nanostructures, UV-Vis spectroscopy was used to characterize the resulting samples in real-time by UV-Vis absorption spectrophotometer. The size and morphology of AuNPs synthesized based on different frequencies were characterized by transmission electron microscope (TEM). The size distribution of synthesized AuNPs was calculated and analyzed by measuring 300 nanoparticles (six groups) in each sample based on TEM micrographs with ImageJ software.


## 4.1. Reagents

HAuCl4·4H2O (precursor), Na3C6H5O7 (reducing agent), and CTAB (cationic surface active agent) were purchased from Sinopharm Chemical Reagent Corporation (Shanghai, China). The other reagents were of the highest grades available commercially. The water used in the experiments was of Milli-Q grade.


## 4.2. Apparatus

UV-Vis absorption spectrophotometer (UV-2501 PC, SHIMADZU, Japan) was used to record the characteristic absorption spectrum of gold colloid. The size and morphology of synthesized AuNPs were characterized by high resolution transmission electron microscopy (TEM TECNAI-12, Philips, Holland). The PZT micropump was fabricated using mature MEMS process in lab [35]. The output flow rate of micropump can be adjusted in the range of 0~6 mL/min, and the working frequency of micropump can be adjusted in the range of 10~400 Hz. The fabrication of Y type micromixer used a computer-numerically controlled (CNC) machining process and a polydimethylsiloxane (PDMS) replication process [4]. The mixing channel and reaction channel were packaged using PDMS and glass. Also, the design and manufacture of micromixer all adopted the optimized structural parameters above. The constant temperature heating platform (BHW-05A shanghai, China) was used to provide the stable reaction temperature. In addition, precision electronic balance and power supply controller were used in experiment test.Figure7 shows the scheme of the experimental setup, the prototype photo of PZT micropumps, and testing photo of controlled synthesis of gold nanoparticles using pulsed mixing microfluidic system. The detailed experimental procedure is as follows. First, all glassware used was clean in freshly prepared aqua regia solution (HCI : HNO3, 3 : 1) and rinsed with deionized water three times before used. Before the synthesis of AuNPs, the chambers of PZT micropump and the channels of micromixer were cleaned by deionized water. The gold precursor, reducing agent, and surfactant were, respectively, dissolved in deionized water to prepare the reaction reagents of the corresponding concentration [36], which include 100 mL HAuCl4 of 1 mM, 100 mL Na3C6H5O7 of 1.5 mM, and 100 mL CTAB of 200 mM. The temperature control platform was adjusted to keep the constant temperature of 100°C. With adjustments of the working voltage, the flow rate of PZT micropumps was fixed at 4 mL/min. To achieve the crossing pulsing micromixing of two reaction reagents, the phase of input signal of PZT micropumps was regulated to maintain 180°. After turning on the power of PZT micropumps, the gold precursor (HAuCl4) and the reducing agent (Na3C6H5O7) were pumped into the micromixer and growth channel by the crossing pulsing state. At the outlet of microchannel, the solution which experienced the “pulsed mixing” and passed through the “growth flow channel” was collected. At the same time, in order to avoid the agglomeration of AuNPs, the corresponding amount of CTAB was also added to the collected samples one time in proportion (nCTAB/nHAuCl4 = 200 : 1). After about ten minutes, the synthesized AuNPs reached stable state, and then the size and morphology of AuNPs can be characterized using the relevant testing apparatus. Note that the experiments were conducted in the frequency of 100, 200, and 300 Hz, respectively.Figure 7
(a) Experimental scheme; (b) the prototype photo of PZT micropumps; (c) testing photo of controlled synthesis of gold nanoparticles using pulsed mixing microfluidic system.
(a)
(b)
(c)


## 4.3. Measurement

To monitor the formation of gold nanostructures, UV-Vis spectroscopy was used to characterize the resulting samples in real-time by UV-Vis absorption spectrophotometer. The size and morphology of AuNPs synthesized based on different frequencies were characterized by transmission electron microscope (TEM). The size distribution of synthesized AuNPs was calculated and analyzed by measuring 300 nanoparticles (six groups) in each sample based on TEM micrographs with ImageJ software.


## 5. Results and Discussion

### 5.1. Test Results

Table4 summarizes the values characterized for AuNPs synthesized at different frequency. (where λm is the wavelength of absorption maximum, αm is the absorption maximum, PWHM is peak width at half maximum, nm¯ is an average diameter of nanoparticles, and σ is the standard deviation of particles size).Table 4
Three groups of experimental results.
Number
f/Hz
λ
m/nm
α
m
PWHM/nm
n
m
¯/nm
σ/nm
1
100
534
0.852
96
25.8
6.08
2
200
532
1.0191
90
28.6
5.47
3
300
524
0.6804
84
21.6
4.83


### 5.2. Analysis of Effect of Pulsing Frequency

Figure8 shows the UV-Vis spectra of three kinds of AuNPs synthesized under different mixing frequencies, when the entrance flow rate is fixed at 4 mL/min. And the corresponding frequencies of the samples (numbers 1, 2, and 3) are 100, 200, and 300 Hz, respectively. As shown in Figure 8, with the increase of pulsing mixing frequency, the absorbance of sample first increases and then decreases. Accordingly, the absorbance maximum of samples rises from 0.852 to 1.091 and then drops to 0.6804. While the absorbance maximum reflects the concentration of particles size, which indicates that the concentration of AuNPs of sample number 2 is the largest, the concentration of AuNPs of sample number 3 is the minimum. As shown in Table 4 and Figure 8, it is also observed that, with the increase of mixing frequency, a continuous blue-shift of the wavelength of absorption maximum occurs. To be more specific, the corresponding wavelength shifts from 534 nm to 532 nm and then blue-shifted to 524 nm. Note that samples (numbers 1, 2, and 3) of synthesized AuNPs showed the typical and unique optical characteristics of gold colloid, and surface plasmon resonance (SPR) by AuNPs was observed at 500~550 nm [37, 38], which further confirms the synthesis of AuNPs. Analyzing the UV-Vis spectra of prepared AuNPs in detail, we can know that the PWHM values of three samples are 96 nm (number 1), 90 nm (number 2), and 84 nm (number 3), respectively. Actually, PWHM value is the peak width at half of absorption maximum, which can characteristic the monodispersity of AuNPs. The lower the PWHM value is, the narrower the size distribution of nanoparticles is [39]. Therefore, it can be inferred that particle diameter deviation of sample number 3 is the lowest, and number 1 is the largest. The inferences are confirmed further through the TEM observations.Figure 8
UV-Vis absorption spectra of gold colloids at different frequencies.Figure9 shows the TEM images and the size distribution histogram of synthesized AuNPs, photos of gold of their corresponding gold colloids. As the photos of Figure 9 show, it can be clearly observed that the color of sample number 1 is pink, number is wine red, and number 3 is light pink, which is in good agreements with the analysis results of UV-Vis absorption spectra presented above. The higher the intensity of absorbance is, the larger the concentration of nanoparticles is. Analyzing the TEM images, we can find that most of AuNPs synthesized by the mixing frequency of 100 Hz are of irregular shape, and the particle monodispersity is poorer. The mean particle size and standard deviation of sample number 1 are 25.8±6.8 nm. When the mixing frequency was set at 200 Hz, the obtained nanoparticles showed good monodispersity, but the morphology of AuNPs still does not reach the expected level. The spherical or quasispherical particles account for a higher ratio, and the mean particle size and standard deviation are 28.6±5.47 nm. With the frequency increasing to 300 Hz, the resulting nanoparticles show better sphericity, most of which are spherical and, meanwhile, have a good monodispersion. The corresponding mean particle size and standard deviation are 21.6±4.83 nm, which indicates that the AuNPs with smaller diameter and narrower size distribution are prepared successfully. The experimental results obtained above indicate that AuNPs synthesized under the mixing frequency of 300 Hz have the optimal size uniformity, morphology, and monodispersity, which are in good accordance with the results of previous simulation optimization. Also, the simulation analysis above showed that mixing in case of 300 Hz had the highest mixing degree. That is to say, the mixing effect under this mixing condition was the best. Exploring the deeper causes of test results, we consider that when the entrance flow rate keeps a certain value, different structure parameters (such as the section size and angle) and different working parameters will lead to kinds of different mixing effects. In the case all parameters are best matched, micropulsing mixing system will achieve the best mixing effects. In this condition, the two solutions layer at the Y type mixing channel became thinner, and the contact areas and opportunities of reaction reagents increased. Then, more fully mixing and reaction will be achieved. As Figure 1 shows, perfect mixing and reaction will directly shorten the nuclear reaction time (stage II: nucleation). At the same time, the total quantity of AuNPs will be ascertained quickly and then reaction entered the next stage (stage III: growth of crystal nucleus). Since two solutions have realize full mixing before entering the growth channel of nanoparticles, the growth environment of nanoparticles is equal, which ensures the final synthesis of nanoparticles with more uniform size and morphology.Figure 9
TEM micrographs and histograms of gold nanoparticles at different frequencies ((a)–(c) inf=100 Hz, f=200 Hz, and f=300 Hz conditions).
(a)
(b)
(c)


## 5.1. Test Results

Table4 summarizes the values characterized for AuNPs synthesized at different frequency. (where λm is the wavelength of absorption maximum, αm is the absorption maximum, PWHM is peak width at half maximum, nm¯ is an average diameter of nanoparticles, and σ is the standard deviation of particles size).Table 4
Three groups of experimental results.
Number
f/Hz
λ
m/nm
α
m
PWHM/nm
n
m
¯/nm
σ/nm
1
100
534
0.852
96
25.8
6.08
2
200
532
1.0191
90
28.6
5.47
3
300
524
0.6804
84
21.6
4.83


## 5.2. Analysis of Effect of Pulsing Frequency

Figure8 shows the UV-Vis spectra of three kinds of AuNPs synthesized under different mixing frequencies, when the entrance flow rate is fixed at 4 mL/min. And the corresponding frequencies of the samples (numbers 1, 2, and 3) are 100, 200, and 300 Hz, respectively. As shown in Figure 8, with the increase of pulsing mixing frequency, the absorbance of sample first increases and then decreases. Accordingly, the absorbance maximum of samples rises from 0.852 to 1.091 and then drops to 0.6804. While the absorbance maximum reflects the concentration of particles size, which indicates that the concentration of AuNPs of sample number 2 is the largest, the concentration of AuNPs of sample number 3 is the minimum. As shown in Table 4 and Figure 8, it is also observed that, with the increase of mixing frequency, a continuous blue-shift of the wavelength of absorption maximum occurs. To be more specific, the corresponding wavelength shifts from 534 nm to 532 nm and then blue-shifted to 524 nm. Note that samples (numbers 1, 2, and 3) of synthesized AuNPs showed the typical and unique optical characteristics of gold colloid, and surface plasmon resonance (SPR) by AuNPs was observed at 500~550 nm [37, 38], which further confirms the synthesis of AuNPs. Analyzing the UV-Vis spectra of prepared AuNPs in detail, we can know that the PWHM values of three samples are 96 nm (number 1), 90 nm (number 2), and 84 nm (number 3), respectively. Actually, PWHM value is the peak width at half of absorption maximum, which can characteristic the monodispersity of AuNPs. The lower the PWHM value is, the narrower the size distribution of nanoparticles is [39]. Therefore, it can be inferred that particle diameter deviation of sample number 3 is the lowest, and number 1 is the largest. The inferences are confirmed further through the TEM observations.Figure 8
UV-Vis absorption spectra of gold colloids at different frequencies.Figure9 shows the TEM images and the size distribution histogram of synthesized AuNPs, photos of gold of their corresponding gold colloids. As the photos of Figure 9 show, it can be clearly observed that the color of sample number 1 is pink, number is wine red, and number 3 is light pink, which is in good agreements with the analysis results of UV-Vis absorption spectra presented above. The higher the intensity of absorbance is, the larger the concentration of nanoparticles is. Analyzing the TEM images, we can find that most of AuNPs synthesized by the mixing frequency of 100 Hz are of irregular shape, and the particle monodispersity is poorer. The mean particle size and standard deviation of sample number 1 are 25.8±6.8 nm. When the mixing frequency was set at 200 Hz, the obtained nanoparticles showed good monodispersity, but the morphology of AuNPs still does not reach the expected level. The spherical or quasispherical particles account for a higher ratio, and the mean particle size and standard deviation are 28.6±5.47 nm. With the frequency increasing to 300 Hz, the resulting nanoparticles show better sphericity, most of which are spherical and, meanwhile, have a good monodispersion. The corresponding mean particle size and standard deviation are 21.6±4.83 nm, which indicates that the AuNPs with smaller diameter and narrower size distribution are prepared successfully. The experimental results obtained above indicate that AuNPs synthesized under the mixing frequency of 300 Hz have the optimal size uniformity, morphology, and monodispersity, which are in good accordance with the results of previous simulation optimization. Also, the simulation analysis above showed that mixing in case of 300 Hz had the highest mixing degree. That is to say, the mixing effect under this mixing condition was the best. Exploring the deeper causes of test results, we consider that when the entrance flow rate keeps a certain value, different structure parameters (such as the section size and angle) and different working parameters will lead to kinds of different mixing effects. In the case all parameters are best matched, micropulsing mixing system will achieve the best mixing effects. In this condition, the two solutions layer at the Y type mixing channel became thinner, and the contact areas and opportunities of reaction reagents increased. Then, more fully mixing and reaction will be achieved. As Figure 1 shows, perfect mixing and reaction will directly shorten the nuclear reaction time (stage II: nucleation). At the same time, the total quantity of AuNPs will be ascertained quickly and then reaction entered the next stage (stage III: growth of crystal nucleus). Since two solutions have realize full mixing before entering the growth channel of nanoparticles, the growth environment of nanoparticles is equal, which ensures the final synthesis of nanoparticles with more uniform size and morphology.Figure 9
TEM micrographs and histograms of gold nanoparticles at different frequencies ((a)–(c) inf=100 Hz, f=200 Hz, and f=300 Hz conditions).
(a)
(b)
(c)


## 6. Conclusions

Continuous synthesis of AuNPs using pulsed mixing microfluidic system has been carried out. The structure and working parameters for Y type micromixer have also been optimized, which includes the optimal analysis of section size, entrance channel angle, and pulse frequency on the mixing effect. According to the optimized structural parameters and working parameters, a corresponding pulsed mixing system prototype has been manufactured, and AuNPs with relatively uniform size and morphology have been successfully synthesized. Through the research above, the following conclusions are drawn:(1)
The pulsed micromixing system actuated PZT micropumps are competent to implement continuous, high-efficiency, controlled synthesis of AuNPs, and the prepared nanoparticles have relatively uniform size, morphology, and good monodispersity.(2)
The structure parameters (e.g., the section size and entrance angle) and working parameters (e.g., frequency, voltage, and phase of electrical signal) are significant important factors affecting the particles size and morphology by using pulsed micromixing system for synthesizing AuNPs. Optimizing the structure and operating parameters are an important technological link in the design and fabrication of pulsing micromixing system. The higher the mixing degree of mixing system is, the better the mixing and reaction effect is, and the particle size of nanoparticles prepared is smaller; size distribution of nanoparticles is also narrower.The preparation of AuNPs in liquid phase is a complex chemical reaction process. Many effect facts such as the concentration of the reactants, pH value, temperature, UV light, and ultrasonic vibration will affect the particle size and morphology of nanoparticles. In the future, the pulsed micromixing system will be further optimized, and the more comprehensive, systematic controlled preparation of nanoparticles will be carried out, on the basis of considering many influence factors.



---

*Source: 160819-2015-11-16.xml*