pubs.acs.org/Langmuir © 2010 American Chemical Society
High-Resolution Direct Patterning of Gold Nanoparticles by the Microfluidic Molding Process Michael T. Demko,*,†,‡ Jim C. Cheng,†,§ and Albert P. Pisano†,‡,§ †
Berkeley Sensor & Actuator Center (BSAC), ‡Department of Mechanical Engineering, and Department of Electrical Engineering & Computer Sciences, University of California at Berkeley, Berkeley, California 94720, United States
§
Received June 3, 2010. Revised Manuscript Received August 6, 2010 A novel microfluidic molding process was used to form microscale features of gold nanoparticles on polyimide, glass, and silicon substrates. This technique uses permeation pumping to pattern and concentrate a nanoparticle ink inside microfluidic channels created in a porous polymer template in contact with a substrate. The nanoparticle ink is self-concentrated in the microchannels, resulting in dense, close-packed nanoparticle features. The method allows for better control over the structure of printed features at a resolution that is comparable to inkjet printing, and is purely additive with no residual layers or etching required. The process uses low temperatures and pressures and takes place in an ambient environment. After patterning, the gold nanoparticles were sintered into continuous and conductive gold traces.
Introduction Precise patterning of nanoparticles is critical for a number of different applications, including low temperature electrode deposition,1 optical coatings and photonics,2,3 biosensors,4 catalysts,5 and MEMS applications.6 Several methods are available for patterning nanoparticles, the most popular of which is inkjet printing.7,8 Inkjet printing is attractive due to its simplicity, high throughput, and low material loss. However, patterning with inkjet printing is limited to a resolution of around 20-50 μm with current printers,8 with higher resolution possible by adding complexity to the substrate prior to printing.9 Electrohydrodynamic printing has been proposed to increase the resolution beyond the limits of inkjet printing, achieving a line resolution as small as 700 nm.10 Both inkjet and electrohydrodynamic printing, however, do not allow precise control over the structure of the printed lines, often resulting in lines with scalloped edges or nonuniform width, and offer only limited control over the height of the printed features.1,8,11,12 Recently, nanoimprint lithography has been proposed as a means of decreasing the feature size of patterned nanoparticles while allowing more precise control over the *To whom correspondence should be addressed. E-mail:
[email protected]. (1) Huang, D.; Liao, F.; Molesa, S.; Redinger, D.; Subramanian, V. J. Electrochem. Soc. 2003, 150, G412–G417. (2) Nakanishi, H.; Bishop, K. J. M.; Kowalczyk, B.; Nitzan, A.; Weiss, E. A.; Tretiakov, K. V.; Apodaca, M. M.; Klajn, R.; Stoddart, J. F.; Grzybowski, B. A. Nature 2009, 460, 371–375. (3) Maier, S. A.; Atwater, H. A. J. Appl. Phys. 2005, 98, 011101. (4) Alivisatos, P. Nat. Biotechnol. 2004, 22, 47–52. (5) Haruta, M. Cattech 2002, 6, 102–115. (6) Fuller, S. B.; Wilhelm, E. J.; Jacobson, J. M. J. Microelectromech. Syst. 2002, 11, 54–60. (7) Ridley, B. A.; Nivi, B.; Jacobson, J. M. Science 1999, 286, 746–749. (8) Sirringhaus, H.; Kawase, T.; Friend, R. H.; Shimoda, T.; Inbasekaran, M.; Wu, W.; Woo, E. P. Science 2000, 290, 2123–2126. (9) Wang, J. Z.; Zheng, Z. H.; Li, H. W.; Huck, W. T. S.; Sirringhaus, H. Nat. Mater. 2004, 3, 171–176. (10) Park, J.-U.; Hardy, M.; Kang, S. J.; Barton, K.; Adair, K.; Mukhopadhyay, D. K.; Lee, C. Y.; Strano, M. S.; Alleyne, A. G.; Georgiadis, J. G.; Ferreira, P. M.; Rogers, J. A. Nat. Mater. 2007, 6, 782–789. (11) Ahmad, Z.; Rasekh, M.; Edirisinghe, M. Macromol. Mater. Eng. 2010, 295, 315–319. (12) Samarasinghe, S. R.; Pastoriza-Santos, I.; Edirisinghe, M. J.; Reece, M. J.; Liz-Marzan, L. M. Gold Bull. 2006, 39, 48–53. (13) Ko, S. H.; Park, I.; Pan, H.; Grigoropoulos, C. P.; Pisano, A. P.; Luscombe, C. K.; Frechet, J. M. J. Nano Lett. 2007, 7, 1869–1877. (14) Ko, S. H.; Park, I.; Pan, H.; Misra, N.; Rogers, M. S.; Grigoropoulos, C. P.; Pisano, A. P. Appl. Phys. Lett. 2008, 92, 154102.
16710 DOI: 10.1021/la1022533
structure of the printed lines.13-15 In this fabrication method, the nanoparticle inks are patterned by pressing with an elastomer mold and the particles dried into their final shape. While the resolution of this method is improved over inkjet printing, there exists a residual layer on the substrate that must be etched away after patterning, and control over the height of features can be frustrated by capillary interactions between the mold and the drying ink, especially along the length of longer features. As an alternative to nanoimprint lithography, nanoparticle self-assembly methods based on capillary filling of photoresist templates have been proposed.16 While these can produce high aspect ratio features with smooth edges, the photoresist must be etched away in subsequent processing steps without removing the particles themselves, which can be technically challenging or nonfeasible in some instances. Here, we demonstrate a nanoparticle patterning method based on permeation pumping17-19 that concentrates nanoparticles in selective regions inside a vapor-permeable polymer mold. This method is completely additive (no etching required) and allows for control over the structure of the patterned lines, including smooth edges and control over the height of the patterned features. The resolution obtained is comparable to that obtained with inkjet printing. Long, continuous lines of gold nanoparticles were patterned over large areas. After patterning, the gold nanoparticles were sintered into conductive traces. This method is compatible with the idea of ambient environment roll-to-roll processing and works on a wide variety of substrates with a wide variety of nanoparticle inks.
Experimental Method Gold nanoparticles encapsulated in a hexanethiol monolayer were synthesized using a two-phase reduction method following the method of Brust et al. and subsequently encapsulated in a hexanethiol self-assembled monolayer and dispersed in an (15) Park, I.; Ko, S. H.; Pan, H.; Grigoropoulos, C. P.; Pisano, A. P.; Frechet, J. M. J.; Lee, E. S.; Jeong, J. H. Adv. Mater. 2008, 20, 489. (16) Choi, S.; Park, I.; Hao, Z.; Holman, H. Y. N.; Pisano, A. P.; Zohdi, T. I. Langmuir 2010, 26, 4661–4667. (17) Randall, G. C.; Doyle, P. S. Proc. Natl. Acad. Sci. U.S.A. 2005, 102, 10813– 10818. (18) Leng, J.; Lonetti, B.; Tabeling, P.; Joanicot, M.; Ajdari, A. Phys. Rev. Lett. 2006, 96, 084503. (19) Moreau, P.; Dehmoune, J.; Salmon, J. B.; Leng, J. Appl. Phys. Lett. 2009, 95, 033108.
Published on Web 10/01/2010
Langmuir 2010, 26(22), 16710–16714
Demko et al.
Article nanoparticle ink began to evaporate, self-concentrating the nanoparticles in the fluid. Eventually, the nanoparticle ink became so concentrated that the ink could no longer flow in the microchannel. At this point, the remaining solvent evaporated and the nanoparticles in the slurry began to close-pack. Some volume reduction of the nanoparticle features occurred at this point, as will be discussed later, but cohesion of the liquid in the solvent ensured that the nanoparticles remained in continuous traces, and the volume reduction occurred uniformly only in directions perpendicular to free surfaces. When the ink was completely dried in the microchannel, the system was cooled and the polymer template removed. Following the nanoparticle patterning, the nanoparticles were sintered by heating in an oven. The nanoparticles were placed in an oven and the temperature increased to 220 C for 8 h. The elevated temperature allows for oxidation of the hexanethiol monolayer and sintering of the gold nanoparticles. The long bake time was necessary to ensure the maximum conductivity of the gold, allowing the diffusion of the oxidized thiols from the interior of the gold traces.
Theoretical Calculations
Figure 1. Microfluidic molding process. alpha-terpineol solvent.20 The specific nanoparticle synthesis method used here has been documented in previous publications.14 The ink was diluted to 3.5 wt % in alpha-terpineol. A master template was made using standard photolithography using SU8-2007 resist on a silicon wafer. The master template consisted of lines that are 15 μm in width, 9.5 μm in height, and 1 cm in length. PDMS (Sylgard 184 - Dow) prepolymer and cross-linker were mixed in a 10:1 ratio, poured over the silicon wafer and cured. After curing, the PDMS elastomer template was removed and cleaned by sonication in ethanol for 10 min and then dried. The resulting PDMS template was porous to solvent vapors. The nanoparticles were patterned on a substrate using the microfluidic molding self-assembly method, which is illustrated schematically in Figure 1. A polyimide, glass, or silicon substrate was placed in a custom-built press, coated with fresh solvent free of nanoparticles and positioned below the patterned PDMS elastomer template. The polymer template was placed on the top side of the press and the fresh solvent patterned by pressing. The pressure used was just enough to completely conform the elastomer mold to the substrate, and did not exceed 14 kPa. This filled the template channels completely with fluid, and excluded the fluid from all other areas. Any residual fluid trapped between the mold and substrate in areas without patterned features was simply evaporated through the porous polymer template in subsequent steps and, since this fluid contained no nanoparticles, left no residual layer. Next, the nanoparticle ink was added to designated filling ports at the ends of the channels. The temperature of the system was raised to 75 C to evaporate the solvent through the porous polymer mold. The evaporating fluid drew the nanoparticle ink into the microfluidic channels, which replaced the fluid lost through evaporation. As the ink was brought into the channel, the solvent of the (20) Brust, M.; Walker, M.; Bethell, D.; Schiffrin, D. J.; Whyman, R. J. Chem. Soc.-Chem. Commun. 1994, 801–802.
Langmuir 2010, 26(22), 16710–16714
The patterning of the nanoparticles in the microfluidic channels can be separated into two distinct stages: filling and drying. In the filling stage, the nanoparticle ink is drawn into the channel and self-concentrated by the evaporating solvent. This continues until the nanoparticle ink saturates at some point in the channel. After this point, part of the channel is saturated with nanoparticles and sustains no additional fluid flow, and part continues to allow fluid flow and self-concentration of the nanoparticle ink. The channel continues to fill with nanoparticles until the entire length of the channel is saturated with particles, at which time the nanoparticle patterning is complete. Each of the two stages was modeled by considering a onedimensional channel. This model is appropriate as long as the length of the channel is much longer than the width and height of the channel to ensure that lateral diffusion across the channel is much stronger than diffusion along the length of the channel, causing a nearly uniform concentration at any given cross section. The velocity in the channel is approximated as a slug flow, and an effective diffusivity used to account for the Taylor dispersion due to the shear-enhanced diffusion of nanoparticles across streamlines in the parabolic velocity distribution. This one-dimensional channel is open to a reservoir of nanoparticle ink at a constant concentration at one end, and closed at the other. For the filling stage, the velocity of the fluid in the channel is found by assuming a constant rate of evaporation of the solvent through the porous walls of the polymer template. This mass flux due to evaporation is represented by a volume flow rate per unit of surface area of the porous polymer template, q00 . To a first approximation, q00 is assumed constant. This approximation ignores the dependence of the evaporation rate on the nanoparticle concentration, which is valid for a sufficiently low concentration nanoparticle ink and not valid for areas with highly concentrated ink toward the end of the channel. Given that the velocity of the fluid at the closed end of the channel must be zero, the velocity as a function of length along the channel is vðxÞ ¼
q00 Pp ðL0 - xÞ A
ð1Þ
where Pp is the part of the perimeter of the microchannel consisting of the porous polymer, A is the cross sectional area of the one-dimensional channel, and L0 is the total length of the channel. The equation for the velocity shows a characteristic length scale equal to the length of the channel, and a characteristic time scale given by ts = A/(q00 Pp). The dimensional variables are DOI: 10.1021/la1022533
16711
Article
Demko et al.
nondimensionalied based on these characteristic scales, resulting in an equation for the velocity given by vðxÞ ¼
L0 ð1 - xÞ ts
ð2Þ
where the quantities with bars above are nondimensional equivalents of the corresponding dimensional variable without the bar above. An equation for the nanoparticle concentration φ(x,t) can be found by considering the flux of nanoparticles by both advection and diffusion. In one dimension, the advective diffusion equation is δφ δ δ δφ þ ðvφÞ ¼ D δt δx δx δx
ð3Þ
where D(x) is the effective diffusivity of the nanoparticles in the one-dimensional system, and will be discussed in detail later. Equation 3 is nondimensionalized based on the constant concentration of the ink in the reservoir φ0 and the length and time scales given above, and eq 2 is inserted. This results in δφ δφ Dts δ2 φ ts dD δφ ¼ φ þ ðx - 1Þ þ þ δt δx L20 δx2 L20 dx δx
ð4Þ
At x = 0, the nanoparticle ink concentration is equal to the concentration of the ink in the reservoir (φh = 1), which is taken to be a constant. The nanoparticle flux is zero at the end of the microcapillary, at x = 1. At t = 0, there are no nanoparticles in the channel, so φh = 0. The magnitude of the diffusivity of nanoparticles in solvent can be estimated using the Einstein-Stokes equation for diffusion of spherical particles through a liquid with low Reynolds number. This equation is given by D0 ¼
kB T 6πηr
ð5Þ
where η is the fluid viscosity, r is the particle radius, kB is the Boltzmann constant, and T is the temperature. However, the value of the diffusivity must be modified to include shearenhanced diffusion caused by velocity gradients perpendicular to the direction of the flow. These velocity gradients enhance the diffusion of nanoparticles as compared to the slug flow assumed here. A correction factor is therefore used to account for this enhanced diffusion21 DðxÞ ¼
1 ½vðxÞrh 2 48 D0
ð6Þ
where rh is the hydraulic radius of the microchannel. The partial differential equation for the nanoparticle concentration was solved numerically using an explicit Euler method. Additionally, in the case where diffusion is negligible throughout the channel, the last two terms in eq 4 are negligible and an analytical solution to the above partial differential equation exists. This solution is given by φðx, tÞ ¼
1 uð1 - expð - tÞ - xÞ 1-x
ð7Þ
(21) Probstein, R. F. Physiochemical Hydrodynamics; Wiley: New York, 1994; p 90.
16712 DOI: 10.1021/la1022533
Figure 2. Nanoparticle ink concentrations versus length and time in the filling stage for ts = 31 s. The dotted black line shows the analytical solution for the case in which D = 0 in the limit as t f ¥.
where u(ζ) is the unit step function. The numerical solution for concentration as a function of position for various times is plotted together with the analytical solution for concentration versus position in the limit of very long times (see Figure 2). The effect of diffusion is seen to cause a smearing of the concentration profile at the end of the channel, and is observed to be stronger near the entrance of the channel due to the stronger velocity gradients and corresponding shear-enhanced diffusion, allowing the concentration gradient to be much steeper at the end of the microchannel. The filling stage continues until the maximum concentration of nanoparticles is reached at some point in the channel, stopping the flow of fluid and nanoparticles in that region. At this point, the drying stage begins and the effective length of the channel with respect to nanoparticle self-concentration and fluid flow, L(t), can no longer be considered a constant, but is rather a function of time. The dimensionless velocity can be modeled as L0 LðtÞ -x vðx, tÞ ¼ ts L0
! ð8Þ
The effective diffusivity is the same as in eq 6, except it is now a function of time as well as position. The nondimensionalized equation for the concentration of nanoparticles can be modified to ! δφ LðtÞ δφ Dts δ2 φ ts δD δφ þ ¼ φ þ xþ δt L0 δx L20 δx2 L20 δx δx
ð9Þ
Again, the nanoparticle ink concentration is taken as constant at the ink reservoir, and the nanoparticle flux at the effective length of the channel is zero. The initial condition is chosen to match the concentration profile from the filling stage when the channel first clogs with nanoparticles. This equation was solved numerically using an explicit Euler method, and the length of the wet channel was recorded as a function of time. Additionally, for the case in which diffusion is negligible everywhere, an analytical solution exists for the length as a function of time. In dimensionless form, this is given by LðtÞ φ0 ¼ exp - t L0 φf
ð10Þ
where φf is the saturation concentration of nanoparticles. The numerical solution to eq 9 for the length of the channel as a Langmuir 2010, 26(22), 16710–16714
Demko et al.
Article
Figure 3. Position of the maximum concentration in the channel as a function of time.
function of time was found to closely follow the analytical solution without diffusion given in eq 10. Note that, after a long time, the effective length of the channel approaches zero, and the entire channel fills uniformly with nanoparticles.
Results To characterize the filling process, optical images were taken using a Canon 5D Mark II camera with a Canon EF 100 mm f/2.8 macro lens with a working distance of 8 cm. Continuous video was taken of the filling process, and the video was subsequently split up into individual frames for analysis. Images were analyzed using Matlab (MathWorks) to convert the images to grayscale and subtract the background from the images. ImageJ was used to analyze the average intensity of the light in the images along the microchannels. The light intensity in the images is proportional to the nanoparticle concentration, with darker areas containing more nanoparticles. Since the location of the maximum light intensity is easily tracked and corresponds identically to the location of the maximum nanoparticle concentration, the location of the maximum light intensity was tracked and plotted as a function of time and compared to the analytical model (see Figure 3). The time scale of the filling process was found by minimizing the least squared error between the theory and data. For the experimental data here, the time scale was found to be approximately ts = 31 s. Once the unsintered gold nanoparticles were successfully patterned, images of the resulting patterns were obtained. An optical microscope (Olympus BH-2) was used to confirm that the pattern was transferred over the entire area of the die. The optical image of the result for the glass and polyimide substrates are shown in Figure 4. A scanning electron microscope (LEO 1550, Zeiss) was used to show the fine details of the resulting pattern and measure the final dimensions of the resulting printed features (see Figure 5). The width of the printed lines, measured using an SEM, were 11 μm. The height of the printed lines were measured by imaging with a 90 tilt, which showed a uniform height of approximately 6.5 μm over the entire area (see Figure 6a). Additionally, the features did show some sag in the center, known in the literature as the rabbitear effect (see Figure 6b). The volume reduction and resulting rabbit-ear height profile was to be expected. At the end of the filling stage, the nanoparticle ink becomes so concentrated that it can no longer flow in the channel, yet additional volume reduction results from further solvent evaporation and possibly some oxidation and removal of thiol groups on the gold nanoparticles. The nanoparticles pack densely together, leaving a gap between the particles and the top of the mold. The presence of this gap causes capillary adhesion between the remaining liquid solvent Langmuir 2010, 26(22), 16710–16714
Figure 4. Optical micrograph of a nanoparticle pattern on (a) glass and (b) polyimide.
Figure 5. SEM image of the top of a line of patterned gold nanoparticles on a polyimide substrate.
Figure 6. SEM image of (a) the side of a patterned line of nanoparticles and (b) a cut cross-section of the patterned particle line showing a slight rabbit-ear effect.
and the solid mold to form the rabbit-ear profile in the final dry structure. However, this effect is small, as the height of the sag in the center of the features was limited to approximately 15% of the total height of the feature. DOI: 10.1021/la1022533
16713
Article
Figure 7. SEM image of the side of a sintered gold line, showing a different height and profile from the unsintered nanoparticles.
Figure 8. Average line resistance as a function of distance. Each point is the average of three measurements on different lines, and the error bars represent the standard deviation of the three measurements.
The gold nanoparticles were then sintered, and the resulting gold traces imaged (see Figure 7). The average height of the sintered gold particles was approximately 2.5 μm. The feature heights exhibited a small amount of waviness over short distances, but the average height of the traces was constant over the entire area. The profile of the sintered gold traces was significantly different from that of the unsintered particles. The reflow of the gold during sintering changed the shape profile from having rabbit-ears to having a dome shape. After the reflow, the average line width was around 16 μm. The sintered gold lines were electrically characterized. The resistance of the lines as a function of distance was measured (see Figure 8). The resulting measurements were linear over large
16714 DOI: 10.1021/la1022533
Demko et al.
distances, showing that the conductivity and the cross-sectional area of the sintered gold are constant over the entire printed area. This result further proves that the initial nanoparticle patterns were uniform over large areas. The cross-sectional area of the gold wires was approximated using a circular segment with the average height and width of the patterned lines. Using the measured resistance per unit length and the calculated cross-sectional area of the gold traces, the average resistivity of the lines was measured to be 6.7 10-7 Ω m. The gold lines are approximately 4% of the conductivity of bulk gold. This value is slightly lower than other values reported in the literature for thermally annealed gold nanoparticles, which may be caused by the larger height of the printed features and the corresponding difficulty for the oxidized hexanethiol monolayer to migrate out of the sintered gold.1,13,15,22-24
Conclusions The proposed method for patterning nanoparticles allows for excellent control over the structure of the printed features, including both line widths and heights. In contrast to other nanoimprint-based patterning methods, this process is designed to form uniform nanoparticle patterns over large areas and leave no residual layer that must be etched away later. Additionally, the nanoparticles can be patterned at low temperatures and pressures in an ambient environment. The process can be done very quickly and the method is compatible with a roll-to-roll processing methodology. Finally, the proposed technology has the potential to be scaled down to achieve higher resolution, while retaining control over the structure of the printed features. Acknowledgment. The authors thank Heng Pan and Costas P. Grigoropoulos for providing the gold nanoparticles, and Timothy P. Brackbill for his assistance with the imaging of the nanoparticle filling process. This work was supported by NSF Grant No. CMMI-0825189. (22) Bieri, N. R.; Chung, J.; Haferl, S. E.; Poulikakos, D.; Grigoropoulos, C. P. Appl. Phys. Lett. 2003, 82, 3529–3531. (23) Wu, Y.; Li, Y.; Liu, P.; Gardner, S.; Ong, B. S. Chem. Mater. 2006, 18, 4627–4632. (24) Cui, W.; Lu, W.; Zhang, Y.; Lin, G.; Wei, T.; Jiang, L. Colloids Surf. A: Physicochem. Eng. Aspects 2010, 358, 35–41.
Langmuir 2010, 26(22), 16710–16714