Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Modeling the surface of fast-cured polymer droplet lenses for precision fabrication

Open Access Open Access

Abstract

Optical lenses with diameter in the millimeter range have found important commercial use in smartphone cameras. Although these lenses are typically made by molding, recent demonstration of fast-cured polymer droplets by inkjet printing has gained interest for cost-effective smartphone microscopy. In this technique, the surface of a fast-cured polydimethylsiloxane droplet obtains dynamic equilibrium via the interplay of surface tension, gravity, thermalization, and a steep viscosity hike. The nature of surface formation involves multiple physical and chemical domains, which represent significant challenges in modeling with the Young–Laplace theory, assuming constant surface tension and viscosity. To overcome these challenges, we introduce the concept of effective surface tension, which allows fast-cured polymer droplets to be modeled as normal liquid droplets with constant viscosity.

© 2018 Optical Society of America

1. INTRODUCTION

Optical lenses with up to 5 mm diameter and focal length are important in smartphone camera modules, which generally have a minimum working distance of 10–20 cm and a corresponding sampling resolution of 100 μm per pixel. A single additional external millimeter lens attached to the camera window can achieve 1 μm resolution, which effectively turns a smartphone camera into a microscope [1]. One simple idea among many others is to fabricate low-cost lenses by free-formed liquid surfaces. In fact, the development of free-formed curable-liquid lenses started in 1995 with moisture-cured polyurethane for surface doming [2]. In 2000, Komatsu et al. deposited high-surface-tension UV-cured resins on a surface to form micro droplets, and subsequently curing to form a microlens array. In 2014, Lee et al. developed a method of creating lenses by inverting the substrate with a deposited resin droplet and relying on gravity to produce significant curvature [3]. In these techniques, the liquid surface has already reached the endpoint of shape evolution and thus represents a steady-state curvature due to the long duration between dispensing and curing. Thus, the final surface curvature can be well predicted by the Young–Laplace (Y-L) model. In 2015, we demonstrated that heat-accelerated curing of polydimethylsiloxane (PDMS) can effectively stop droplet surface evolution, while it still possesses significant surface curvature. This single-step method has enabled high-throughput, high-quality, low-cost lens fabrication [4,5]. Since then, a variety of methods sharing the same or similar principles have been developed [69]. Most recently, DotLens smartphone nanocolorimetry has been demonstrated for lead sensing in drinking water [10]. In these techniques, it is crucial to understand the relationship between droplet surface geometry and fabrication parameters. To date, only empirical methods have been employed [11]. As mentioned earlier, the challenge of modeling surface curvature arises from the rapid viscosity hike during the heat-accelerated curing process, which violates the major assumption in the Y-L model. On the other hand, time-domain numerical modeling is difficult due to the multiple physical and chemical domain interactions of heat transfer, surface tension, viscosity and phase change, and of course, gravity [12]. Nevertheless, our observations suggest that polymer droplet lenses produced by fast curing does appear like normal liquid droplets as shown in Fig. 1. This insight has motivated us to introduce the concept of effective surface tension (EST) for a fast-cured droplet [13]. In other words, EST is the surface tension of the cured droplet as if it were still in the liquid form, and can be employed in the Y-L model.

 figure: Fig. 1.

Fig. 1. 50 μL droplets on surfaces. (a) Water on glass slide (left) and cover slip glass (right) exhibits different contact angles due to different surface energies; (b) fast-cured PDMS on the same glass slide at 100°C (left) and 200°C (right) exhibits different contact angles due to different curing speeds.

Download Full Size | PDF

2. MODELING LIQUID DROPLETS

To quantify the lensing effect of a droplet situated between an air–solid interface, we first consider the interfacial energy determined by the Young’s equation γsaγslγlacosθc=0, where γsa is the solid surface free energy, γsl is the solid–liquid interfacial tension, γla is the liquid–air surface tension, and θc is the contact angle. Second, its curvature at equilibrium can be determined by solving the Y-L equation Δp=γla¯·n^=γla(Rx1+Ry1) via surface energy minimization [14], where Δp is Laplace pressure, n^ is unit normal pointing out of surface, and Rx and Ry are radii of curvature in the axes parallel to surface. However, the Y-L equation has an analytical solution only in the simplest cases and does not consider other influential forces such as gravity; thus, the interface shape is usually calculated numerically. A widely used numerical modeling software is Surface Evolver [15], which simulates steady-state droplet shape by surface energy minimization subject to different forces and constraints in addition to surface tension, e.g., gravity [1618]. Once the droplet surface is determined, lens parameters such as focal length (f) and numerical aperture (NA) can be calculated. Here we consider a droplet having the shape of a spherical cap resting on a flat substrate. The droplet shape is controlled by the Bond number Bo=(ρg/γla)×rh, where ρ is liquid density (ρPDMS=960kg/m3), g is gravitational constant, r is radius, h is height, and γ is liquid surface tension. For small droplets (less than 5 μL), the Bond number is less than 1, and surface tension forces dominate, so droplets can be accurately regarded as a spherical cap described by V=16πh(3r2+h2), where V is volume, h is height, and r is radius.

Here, water droplets of different volumes were deposited on glass surface with θc=90°±2.6°, as shown in Fig. 2(a) (inset, lower row). The generated model is shown for comparison (inset, upper row). The measured and simulated droplet height versus radius are shown in Fig. 2(a) (bottom line), and the corresponding volumes of the droplet show a good fit to the spherical cap equation for up to 10 μL before the gravitational effect starts to flatten the droplet [19]. For a different glass surface (Corning Gorilla Glass 3) that exhibits θc=110°±3°, the measured and simulated droplet variables are shown in Fig. 2(a) (top line). For fixed surface tension, curvature is highly dependent on volume.

 figure: Fig. 2.

Fig. 2. (a) Deposited and modeled measurements of pure water droplets on glass. (i) 0.05 μL, (ii) 0.1 μL, (iii) 0.2 μL, (iv) 0.5 μL, (v) 1 μL, (vi) 2 μL, (vii) 5 μL, (viii) 10 μL, (ix) 50 μL, (x) 100 μL, (xi) 200 μL, and (xii) 500 μL. (b) Droplet curvature versus volume.

Download Full Size | PDF

The curvature of the droplet was obtained by imaging and parabolic curve fitting from the droplet vertex. The numerical curvature was obtained as 1/R for a parabola in the form of 2y=1Rx2, and the results are shown in Fig. 2(b). The focal length can be obtained with the lensmaker’s equation, simplified as f=R/(n1) for a plano-convex lens. Smaller droplets possess larger curvatures that result in larger spherical aberrations; therefore, we defined the aperture as 1/3 the lens diameter. The curvature obtained this way typically deviates less than 5% from the curvature at the vertex.

In Fig. 2(a), as the droplet volume increases, its radius and height both increase as expected; however, the height eventually plateaus due to gravity, and additional volume creates only lateral expansion. The maximum puddle height that a liquid can attain on a surface is related to the energy minimization of surface tension and gravity, and its value can be calculated by Eq. (1):

hmax=2γlaρg(1cosθc),
where hmax is maximum puddle height that a liquid puddle can attain for an infinitely large droplet, γla is liquid surface tension, θc is contact angle, ρ is liquid density, and g is gravitational constant [20]. For water droplets on the selected substrates, hmax is, respectively, 3.83 mm (θc=90°) and 4.43 mm (θc=110°) [21]. For a lower surface tension liquid such as PDMS (γla=20mN/m) on glass, the contact angle rarely exceeds 30°, which yields negligible curvature and thus optical power, as shown in Fig. 3. For this reason, accelerated curing of curable liquid polymers was developed as a method to stop the surface evolution before the equilibrium is reached [11].

 figure: Fig. 3.

Fig. 3. (a) Maximum puddle height for water and liquid PDMS at different contact angles.

Download Full Size | PDF

3. MODELING CURABLE LIQUID POLYMER DROPLETS

Traditional methods of droplet modeling are applicable to slowly cured PDMS if the droplet is allowed to reach static equilibrium. For example, at 80°C, PDMS requires 1 min to cure, which ensures that static equilibrium is attained or nearly attained, as shown in Fig. 4(a). However, this model is oversimplified when a PDMS droplet undergoes fast curing on a heated substrate, while the surface evolution stops at a point of dynamic equilibrium, as shown in Fig. 4(b).

 figure: Fig. 4.

Fig. 4. (a) 80°C: droplet attains static equilibrium before curing; normal low-temperature curing of a PDMS droplet can be modeled by its static equilibrium state. (b) 200°C: accelerated curing prevents droplet from reaching equilibrium; fast-cured PDMS causes material to undergo chemical phase change in a state of dynamic equilibrium. (c) Trend in geometry versus volume and temperature; varying substrate temperature and volume creates predictable changes in droplet radius, height, and curvature. Scale bars are 5 mm.

Download Full Size | PDF

Despite this, observations indicate that fast-cured polymer droplets retain the appearance of a liquid droplet, and an empirical trend is observed between droplet geometries (radius, height, contact angle, and curvature) and its volume and substrate temperature, as shown in Fig. 4(c). For example, a 100 μL PDMS droplet shows decreasing radius and increasing height with increasing temperatures due to increased contact angle, similar to water droplets deposited on increasingly hydrophobic substrates. Droplets cured at 200°C show increasing radius and height (until plateau) with increasing volumes, also similar to water droplets with increasing volumes.

A series of cured PDMS lenses has been fabricated by a range of droplet volumes (12 values) and temperatures (13 values). In Fig. 5(a), the droplet radius was measured: for a fixed temperature, the radius increases with volume as expected due to the larger droplet. For a fixed volume, the radius decreases with temperature as the droplet cures before spreading out; this also indicates a height increase to offset for constant volume. In Fig. 5(b), the droplet height was measured: for a fixed temperature, the height increases with volume, but effectively plateaus in a fashion similar to the maximum puddle height equation. For example, at 300°C, the maximum height that can be achieved with very large volumes was 3.61 mm, while at 140°C the maximum height was 2.35 mm. Beyond these heights, the volume increase causes only lateral expansion. For a selected volume, the temperature increase causes height increase before it plateaus at high temperature.

 figure: Fig. 5.

Fig. 5. Fast-cured PDMS droplets exhibit controllable geometry with changing volume and temperature. (a) Droplet radius, (b) height, and (c) contact angle, respectively, plotted versus substrate temperatures and volumes.

Download Full Size | PDF

In Fig. 5(c), the droplet contact angle was measured: for a fixed temperature, the contact angle is independent of the droplet volume, as surface energy is not related to volume in the Young’s equation. For a fixed volume, the contact angle displays itself as a function of substrate temperature, which is another mode of surface energy as defined by the Young’s equation. This is a unique phenomenon observed in fast-cured polymer droplets. The relationship between contact angle and temperature can be fitted by a sigmoidal function, shown in Eq. (2). The ability to tune the droplet contact angle by changing temperature makes it possible to obtain a wide range of usable PDMS droplet curvatures for lensing, overcoming the inherent problem of small contact angles in slowly cured PDMS droplets:

θc=f(T)=a1+ebT+c=1801+e0.0185T+2.614.
Two observations point to the similarities of fast-cured liquid PDMS droplets to normal liquid droplets. First, the droplets attain a maximum puddle height by which any volumetric increase results in strictly lateral expansion [Fig. 5(b)]. Second, the contact angle depends solely on temperature, i.e., surface energy [Fig. 5(c)], similar to the Young’s equation. Thus, we propose a droplet model that views fast-cured polymer droplets as normal liquid droplets with an EST, and circumvent the complex curing process. The EST is calculated by measuring the maximum puddle height attained by droplets at every temperature, and its contact angle at that temperature. Upon reorganization of Eq. (1), the EST is calculated by Eq. (3):
γeff=gρ2hmax2(1cosθc)=gρ2hmax2(1cos(1801+e0.0185T+2.614)).
In general, surface tension in a normal liquid decreases with increasing temperatures due to reduction in cohesive forces caused by increase in molecular thermal activity. In contrast, EST of a fast-cured polymer droplet increases with temperature due to the polymer cross-linkage reaction rate, as shown in Fig. 6(a). By using these EST values and modeling droplets in Surface Evolver, the results for different volumes and temperatures are shown in Fig. 6(b). The radius, height, and curvature of empirical and modeled droplets exhibit a good fit for all volumes and temperature combinations, as shown in Figs. 7(a) and 7(b). Particularly, it can be seen that the same contact angle can be obtained at different volume–temperature combinations, as shown in Fig. 7(c). A comprehensive graph for different variable combination modeling is shown in Fig. 8(a). The error between the modeled and experimental results are shown in Fig. 8(b) and was <10% in most cases, which was largely attributed to images where edges were not clearly defined. Some data points with volume <1μL at low temperatures yielded results that were hard to image and analyze due to their very small flat geometry, and were omitted.

 figure: Fig. 6.

Fig. 6. (a) EST of fast-cured PDMS versus temperature; lower values have higher error due to shorter height. (b) Geometry of fast-cured PDMS droplets modeled as normal liquid droplets with EST.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. Droplet curvatures for different printing variables. Droplets marked with blue circles are the same. The 50 μL 120°C droplet has similar curvature to 75 μL 200°C droplet (blue diamond). (a) 50 μL volume, various temperatures. (b) 200°C temperature, various volumes.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. Comprehensive results for (a) curvature versus printing variables, and (b) error between modeled and actual droplets. (*) No data as droplets too small and flat.

Download Full Size | PDF

4. CONCLUSION

Curable polymers such as PDMS possess low surface tension and exhibit small contact angles on the air–solid interface, leading to negligible surface curvature. However, heat-induced fast-curing causes accelerated solidification in a state of dynamic equilibrium results in significant surface curvature and lensing effect. Due to the dynamic nature of the fast-curing process, the Y-L model, which assumes constant surface tension, is not applicable. Here, we model fast-cured polymer droplets as normal liquid droplets that possess an EST, a static surface tension value that the cured droplets appear to exhibit. Using experimentally extracted EST, the final cured droplet surface can be modeled accurately using a standard surface evolution program.

Funding

National Science Foundation (NSF) (1558240, 1643391).

REFERENCES

1. Y. Sung, F. Campa, and W. C. Shih, “Open-source do-it-yourself multi-color fluorescence smartphone microscopy,” Biomed. Opt. Express. 8, 5075–5086 (2017). [CrossRef]  

2. W. Emmerling, T. Podola, L. Unger, and M. Majolo, “Moisture-curing alkoxysilane-terminated polyurethanes,” U.S. patent 5,525,654 (10 September 1996).

3. W. Lee, A. Upadhya, P. Reece, and T. G. Phan, “Fabricating low cost and high performance elastomer lenses using hanging droplets,” Biomed. Opt. Express 5, 1626–1635 (2014). [CrossRef]  

4. Y. Sung, J. Jeang, C. H. Lee, and W. C. Shih, “Fabricating optical lenses by inkjet printing and heat-assisted in situ curing of polydimethylsiloxane for smartphone microscopy,” J. Biomed. Opt. 20, 047005 (2015). [CrossRef]  

5. W. C. Shih and Y. Sung, “Fabrication of lenses by droplet formation on a pre-heated surface,” U.S. patent 9,995,851 (12 June 2018).

6. X. Meng, H. Huang, K. Yan, X. Tian, W. Yu, H. Cui, Y. Kong, L. Xue, C. Liu, and S. Wang, “Smartphone based hand-held quantitative phase microscope using the transport of intensity equation method,” Lab Chip 17, 104–109 (2017). [CrossRef]  

7. S. Ekgasit, N. Kaewmanee, P. Jangtawee, C. Thammacharoen, and M. Donphoongpri, “Elastomeric PDMS planoconvex lenses fabricated by a confined sessile drop technique,” ACS Appl. Mater. Interfaces 8, 20474–20482 (2016). [CrossRef]  

8. T. Kamal, R. Watkins, Z. Cen, J. Rubinstein, G. Kong, and W. M. Lee, “Design and fabrication of a passive droplet dispenser for portable high resolution imaging system,” Sci. Rep. 7, 41482 (2017). [CrossRef]  

9. R. Amarit, A. Kopwitthaya, P. Pongsoon, U. Jarujareet, K. Chaitavon, S. Porntheeraphat, S. Sumriddetchkajorn, and T. Koanantakool, “High-quality large-magnification polymer lens from needle moving technique and thermal assisted moldless fabrication process,” PloS One 11, e0146414 (2016). [CrossRef]  

10. H. Nguyen, Y. Sung, K. O’Shaughnessy, X. Shan, and W.-C. Shih, “Smartphone nanocolorimetry for on-demand lead detection and quantitation in drinking water,” Anal. Chem. 90, 11517–11522 (2018). [CrossRef]  

11. Y. Sung, J. Garan, H. Nguyen, Z. Hu, and W. C. Shih, “Automated batch characterization of inkjet-printed elastomer lenses using a LEGO platform,” Appl. Opt. 56, 7346–7350 (2017). [CrossRef]  

12. J. White, M. Santos, M. Rodríguez-Valverde, and S. Velasco, “Numerical study of the most stable contact angle of drops on tilted surfaces,” Langmuir 31, 5326–5332 (2015). [CrossRef]  

13. D. Zang, Z. Chen, Y. Zhang, K. Lin, X. Geng, and B. P. Binks, “Effect of particle hydrophobicity on the properties of liquid water marbles,” Soft Matter 9, 5067–5073 (2013). [CrossRef]  

14. H. Liu and G. Cao, “Effectiveness of the Young-Laplace equation at nanoscale,” Sci. Rep. 6, 23936 (2016). [CrossRef]  

15. K. A. Brakke, “The surface evolver,” Exp. Math. 1, 141–165 (1992). [CrossRef]  

16. M. Santos and J. White, “Theory and simulation of angular hysteresis on planar surfaces,” Langmuir 27, 14868–14875 (2011). [CrossRef]  

17. M. Santos, S. Velasco, and J. White, “Simulation analysis of contact angles and retention forces of liquid drops on inclined surfaces,” Langmuir 28, 11819–11826 (2012). [CrossRef]  

18. G. Soligno, M. Dijkstra, and R. van Roij, “The equilibrium shape of fluid-fluid interfaces: derivation and a new numerical method for Young’s and Young-Laplace equations,” J. Chem. Phys. 141, 244702 (2014). [CrossRef]  

19. Y. Yonemoto and T. Kunugi, “Theoretical model of droplet wettability on a low-surface-energy solid under the influence of gravity,” Sci. World J. 2014, 647694 (2014). [CrossRef]  

20. F. Brochard-Wyart and D. Quéré, Capillary and Wetting Phenomena-Drops, Bubbles, Pearls, Waves (Translated by A. Reisinger, Springer, 2002).

21. L. R. White, “On deviations from Young’s equation,” J. Chem. Soc., Faraday Trans. 1 73, 390–398 (1977). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1.
Fig. 1. 50 μL droplets on surfaces. (a) Water on glass slide (left) and cover slip glass (right) exhibits different contact angles due to different surface energies; (b) fast-cured PDMS on the same glass slide at 100°C (left) and 200°C (right) exhibits different contact angles due to different curing speeds.
Fig. 2.
Fig. 2. (a) Deposited and modeled measurements of pure water droplets on glass. (i) 0.05 μL, (ii) 0.1 μL, (iii) 0.2 μL, (iv) 0.5 μL, (v) 1 μL, (vi) 2 μL, (vii) 5 μL, (viii) 10 μL, (ix) 50 μL, (x) 100 μL, (xi) 200 μL, and (xii) 500 μL. (b) Droplet curvature versus volume.
Fig. 3.
Fig. 3. (a) Maximum puddle height for water and liquid PDMS at different contact angles.
Fig. 4.
Fig. 4. (a) 80°C: droplet attains static equilibrium before curing; normal low-temperature curing of a PDMS droplet can be modeled by its static equilibrium state. (b) 200°C: accelerated curing prevents droplet from reaching equilibrium; fast-cured PDMS causes material to undergo chemical phase change in a state of dynamic equilibrium. (c) Trend in geometry versus volume and temperature; varying substrate temperature and volume creates predictable changes in droplet radius, height, and curvature. Scale bars are 5 mm.
Fig. 5.
Fig. 5. Fast-cured PDMS droplets exhibit controllable geometry with changing volume and temperature. (a) Droplet radius, (b) height, and (c) contact angle, respectively, plotted versus substrate temperatures and volumes.
Fig. 6.
Fig. 6. (a) EST of fast-cured PDMS versus temperature; lower values have higher error due to shorter height. (b) Geometry of fast-cured PDMS droplets modeled as normal liquid droplets with EST.
Fig. 7.
Fig. 7. Droplet curvatures for different printing variables. Droplets marked with blue circles are the same. The 50 μL 120°C droplet has similar curvature to 75 μL 200°C droplet (blue diamond). (a) 50 μL volume, various temperatures. (b) 200°C temperature, various volumes.
Fig. 8.
Fig. 8. Comprehensive results for (a) curvature versus printing variables, and (b) error between modeled and actual droplets. (*) No data as droplets too small and flat.

Equations (3)

Equations on this page are rendered with MathJax. Learn more.

hmax=2γlaρg(1cosθc),
θc=f(T)=a1+ebT+c=1801+e0.0185T+2.614.
γeff=gρ2hmax2(1cosθc)=gρ2hmax2(1cos(1801+e0.0185T+2.614)).
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.