Compact micro/nano electrohydrodynamic patterning: using a thin conductive film and a patterned template

b The influence of electrostatic heterogeneity on the electric-field-induced destabilization of thin ionic liquid (IL) films is investigated to control spatial ordering and to reduce the lateral dimension of structures forming on the films. Commonly used perfect dielectric (PD) films are replaced with ionic conductive films to reduce the lateral length scales to a sub-micron level in the EHD pattering process. The 3-D spatiotemporal evolution of a thin IL film interface under homogenous and heterogeneous electric fields is numerically simulated. Finite diﬀerences in the spatial directions using an adaptive time step ODE solver are used to solve the 2-D nonlinear thin film equation. The validity of our simulation technique is determined from close agreement between the simulation results of a PD film and the experimental results in the literature. Replacing the flat electrode with the patterned one is found to result in more compact and well-ordered structures particularly when an electrode with square block protrusions is used. This is attributed to better control of the characteristic spatial lengths by applying a heterogeneous electric field by patterned electrodes. The structure size in PD films is reduced by a factor of 4 when they are replaced with IL films, which results in nano-sized features with well-ordered patterns over the domain.


Introduction
Electrohydrodynamic (EHD) destabilization of molten polymer or liquid films has received extensive attention over the past few decades as a unique and interesting approach for creating micron-and submicron-sized features for soft lithography [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] and coatings. 17,18In this process, a thin liquid film is sandwiched between two electrodes and an electric field is applied to the film in the transverse direction (see Fig. 1).The disparity of electrical properties of the film and the bounding fluid (the material which fills the gap between the film and the upper electrode) results in a different electric field in each layer.As a consequence, a net electrostatic force acts at the interface due to the Maxwell stress. 19In the EHD patterning process, thicker regions of the film are subjected to a higher electrostatic force compared to thinner regions, thus the electrostatic force increases the instabilities.In contrast, capillary and viscous forces tend to dampen the instabilities.When the electrostatic force overcomes these damping forces, instabilities grow and pillars (raised columnar structures) form at the interface. 1itial linear stages in the growth of instabilities were investigated through linear stability (LS) analysis. 12][3] The LS analysis shows that the lateral size of pillars decreases by increasing the applied voltage (V), mean initial film thickness (h 0 ), electric permittivity of the film (e) and its conductivity (s).Typically in the EHD patterning process, the thin liquid film is considered as either a perfect dielectric (PD) or a leaky dielectric (LD). 1,6,11,20,21Smaller length scales can be generated by introducing the temperature gradient for combined thermal and electric field induced instability for PD and LD thin films. 22n the PD film there is no free charge and the film is polarized under an applied electric field.In contrast, for the LD film the infinitesimal amount of charge present moves and accumulates at the interface.Both LS and non-linear analysis of EHD patterning demonstrate that the presence of free charges, albeit infinitesimal, reduces the characteristic wavelength l. 4,5,11,20,21 This is attributed to a high electrostatic force arising from the addition of an extra term to the Maxwell stress. 16,235][26][27] Subsequently, a drop is emitted from the cone's tip due to jet breaking.Such a EHD cone-jetting phenomenon, known as electrospraying, has been used extensively as a well-known technique of electrospray mass spectrometry. 28,29n the context of EHD patterning, the formation of the DC Taylor cone has just been reported (to the best of authors' knowledge) for the fabrication of the bifocal microlens using EHD reflowing of pre-patterned polymer films, 30 where the steady state pattern formation is achieved using critical electric field intensity.The Taylor cone shape is observed when the electric field intensity slightly exceeds the critical value of 10 MV m À1 resulting in switching from ordinary single-focus microlenses to the bifocal ones.It should be noted that in the EHD patterning a finite mass is exposed to the electric field, in contrast with electrospraying, and the creeping flow (high viscosity and low velocity) hinders tip streaming and drop emissions.
The electric breakdown of a polymer film and the bounding fluid at a high electric field strength is reported as a limiting factor in EHD patterning especially when the fabrication of nanosized features is needed. 312][33] Employing ionic conductive polymers or simply ionic liquid (IL) films as a bounding fluid enhances the electrostatic pressure acting at the interface while preventing electric breakdown. 31It is shown that in the polymer-ionic liquid bilayer system, the maximum wavelength (l) can be decreased to a submicron level with an increase in the applied voltage and/or a decrease of polymer film thickness.However, it must be noted that thinner films have a higher tendency of dielectric breakdown mainly when the film thickness is less than 50 nm. 315][36][37] To mathematically model the ion distribution within the IL film the Poisson-Nernst-Plank (PNP) model has been used. 23,31,35When the process time is much shorter than the free charge/ion relaxation time the PNP equation can be simplified to the Poisson-Boltzmann equation which relates the electric potential to the charge distribution across the film. 23,31ree ions present in ILs move easily and form a diffuse layer called a double layer (DL) adjacent to the charged surfaces.In the PD and LD models, it is assumed that the thickness of the DL is infinitesimally large and small, respectively.Initially, to bridge the gap between these two limiting cases, an electrokinetic model was developed for drop deformation in an electric field. 38,39ecently an electrokinetic model has also been developed to study EHD patterning and LS analysis is used to investigate the dynamics of polymer films. 23These electrokinetic models apply the linearized Debye-Huckel approximation 40 to study the stability of IL films.In our previous work, a more general model has been developed based on the nonlinear Poisson-Boltzmann equation to investigate the dynamics and stability of IL-PD bilayers. 12,15The EHD patterning process has also been conducted for ion conductive films which resulted in the creation of nano-sized pillars on the film even for high filling ratios (initial film thickness to electrode distance ratio). 16However, the formed pillars were polydispersed in the domain and less likely to have the array of pillars uniformly distributed over the domain.
In order to create well-ordered structures via electrically induced patterning, it is common to use a patterned mask which causes replication of mask on the polymer film. 2,3,7,13,14n this study, a patterned upper electrode is used to electrostatically control the spatial ordering of pillars as well as their sizes in an IL film.First, spatiotemporal evolution of thin PD and IL films in the absence of electrostatic heterogeneity (by using flat electrode) is presented as a baseline.Then the influence of electrostatic heterogeneity on the electric-fieldinduced destabilization of thin IL and PD films is presented.The effect of protrusion shape (strip and square block protrusions) and their size (see Fig. 1(b and c)) on the dynamics, morphological evolution and the pattern formation process of IL films is investigated and the results are compared with the PD base cases.

Mathematical model
A schematic of a thin liquid film resting on a lower electrode and confined with a patterned upper electrode (i.e.mask) is shown in Fig. 1(a).The gap between the film and the mask can be filled with air or a polymer.It is assumed that the material properties of the layers, such as interfacial tension, viscosity and dielectric permittivity, are constant.Gravity effects are neglected and it is assumed that the liquid film is incompressible Newtonian fluid.Constants and parameters used in this study are presented in Table 1.
Continuity and Navier-Stokes (NS) equations are used to describe the thin film dynamics: This journal is © The Royal Society of Chemistry 2016 where f e = Àrf stands for the external body force which consists of the intermolecular and electrostatic forces with the term f being the conjoining/disjoining pressure.The subscript i = 1, 2 represents the ith fluid phase.Boundary conditions are: The kinematic boundary condition for the vertical component of velocity at the interface is: 41 In this study, air is considered as the bounding fluid.Air is considered dynamically inert since it has much lower viscosity than the polymer film.Moreover, the Reynolds number in the EHD patterning process is much less than unity due to the small length scale of the system.The film motion can thus be described by Stokes flow.Combining these assumptions and applying lubrication theory to the thin film dynamics 42 result in the following nonlinear thin film equation: In this equation, subscripts x and y denote the spatial derivatives and subscript t denotes the time derivative.Details of these derivations are available in the literature. 2,15,16,43,44The conjoining/disjoining pressures (f) are the summation of van der Waals, Born repulsion and electrostatic pressures and are defined as The term, f vdW = A/6ph 3 , is the van der Waals pressure and A is the effective Hamaker constant defined as in which 1, 2 and 3 denote the substrate, the liquid film and the bounding fluid.The term f Br = 8B U /(d À h) 9 is the Born repulsion acting at the interface to avoid the contact line singularities when the interface touches the upper electrode.Coefficient B U is calculated by setting the conjoining pressure (f) equal to zero at an equilibrium distance (l 0 ).The electrostatic component of the conjoining pressure (f E ) is found based on the Maxwell stress acting at the fluid interface 19 as: To find the electrostatic pressure in eqn (7), it is required to find the electric potential at each layer.In the IL layer, the ion conservation equation is considered to find the ion/free charge dynamics and distribution as: where r f stands for free charge density, s is conductivity and D is the diffusion coefficient.Using scaling factors of The process time is t p = mh 0 /e 0 e 1 V 0 2 , and the charge relaxation time is t c = e 0 e 1 /s.In this study, t p and t c are in the order of seconds and micro-seconds, respectively.This diminishes the effect of the free ion (charges) dynamics during the pattern formation process and the ion conservation equation (eqn (8)) is simplified to the Poisson-Boltzmann equation: 40 is the Boltzmann distribution of free charges or ions and is zero for air.V ref is the reference potential which is found based on the electroneutrality conditions within the IL film, 40 Ð h 0 r f dz ¼ 0. The term k B T represents the thermal motion where k B is the Boltzmann constant and T (K) is the temperature.The bulk ion number concentration n N = 1000N A M depends on the electrolyte molarity, M (mol L À1 ).
For the IL layer, the Poisson-Boltzmann equation and for air as PD bounding fluid the Laplace equation in the long-wave limit form are given as with the following boundary conditions: The Debye-Huckel approximation was used to simplify and analytically solve the electrostatic governing equations 16,23  (eqn ( 11)).This results in a relation for the electrostatic pressure acting at the interface of the IL film bounded with air as where e r ¼ e 1 e 2 and V s is calculated as: In the presented mathematical model, the following points are worth mentioning: (a) all simulations are terminated at 1 nm (minimum equilibrium distance), which is considered as the limit of the continuum assumption, (b) the analogy of symmetric mono-valent electrolytes is used to analytically relate the electrostatic force and DL thickness in ionic liquids, (c) the height between the two electrodes, and the finite volume of the confined liquid film result in divergence of the hydrodynamic retardation as the electric field intensity diverges which avoids Taylor-cone formation and tip streaming.

Scaling
The thin film equation (eqn (6)) is nondimensionalized using scaling parameters obtained from LS analysis 15,16 as: Then the thin film thickness (i.e.interface height) is scaled with its mean initial thickness, H = h/h 0 .The corresponding wavelength for growth of instabilities is:

Numerical method
To find the spatiotemopral evolution of the thin IL film, the 2D nonlinear thin film equation is solved numerically.Derivatives in spatial coordinates are discretized using finite difference and the resulting set of differential algebraic equations (DAE) are solved in time using an adaptive time step ODE solver.A square 3 Results and discussion

Morphologies in the homogeneous electric field
In this section, an overview of the EHD patterning process under a homogeneous electric field is provided.The effects of ionic conductivity and the film initial filling ratio (mean initial film thickness to electrodes distance ratio) on the shape and size of structures are shown.Snapshots in time of the 3-D patterns for thin PD and IL liquid films (d 1 = d 2 = 100 nm) are presented in Fig. 3. Results show that the addition of ionic conductivity to the film results in an increase in the number of pillars forming at the interface and consequently more compact structures (pillars with smaller size) are created.The number of pillars formed in the IL films is also increased at higher filling ratios.The PD film shows a similar trend below the critical filling ratio of h 0 /d = 0.5.More details regarding the process of pattern formation and the structure size under a homogeneous electric field in ionic conductive films are presented elsewhere.
Base case studies demonstrate that the addition of ionic conductivity to the polymer film leads to more compact but less ordered pillars.Coarsening of the structure occurs randomly over the domain 16 similar to the PD 46,47 films resulting in nonuniform patterns.In PD films, the hexagonally packed pillars typically span about 10 to 15 pillars or they are randomly distributed over the domain. 46To overcome this limitation, a patterned mask is used to electrostatically control the structure size and fluid patterns.The use of a patterned mask, which results in the replication of the mask on the polymer film, has already been investigated for the PD and LD films. 2,3,6,7,9,11,13,14,48abricating controllable submicron sized structures using the EHD patterning process is the main objective of this study and the influence of electrostatic heterogeneity on the EHD pattern formation of thin IL films is investigated.

Conjoining/disjoining pressure
The final shape and size of features in the EHD patterning process are a result of a balance between the electrostatic (as the main driving force), intermolecular, viscous and interfacial force. 49The distribution of conjoining/disjoining pressure components acting at the film interface under a heterogeneous electric field is compared in Fig. 4. The 2-D cross-section height profile of the upper electrode and the patterned film is shown in Fig. 4(a and b), respectively.The Laplace pressure (Fig. 4(c)), the van der Waals pressure (Fig. 4(d)) and electrostatic pressure for the PD (Fig. 4(e)) and IL (Fig. 4(f)) films are all scaled with the electrostatic pressure acting on the PD film at the initial stage (flat interface) when exposed to the homogeneous electric field (d 1 = d 2 = 100 nm and h 0 = 30 nm).
Negative (positive) pressure values in Fig. 4 represent the upward (downward) direction in which the interface is pushed toward the upper (lower) electrode.The electrostatic pressure has a negative value over the domain irrespective to the interface pattern but its magnitude changes significantly with the interface height (Fig. 4(e and f)).The film is exposed to a higher electrostatic force under protrusions (decreased electrode distance) caused by a lower electric potential drop in the bounding layer.Thinner film regions also experience lower electrostatic pressure compared to the thicker regions (located under protrusions), which results in fluid flowing from thinner regions to thicker regions.This confirms the destabilizing role of the electrostatic pressure in the EHD patterning process.
The curvature of the interface determines the Laplace pressure with a convex shape of the interface (in the valleys) resulting in negative values for pressure and vice versa (Fig. 4(c)).Film regions with a higher (lower) interface height experience a pressure toward the lower (upper) electrode.As a consequence, the Laplace pressure tends to flatten the interface, thereby stabilizing the film in contrast to the destabilizing effect of electrostatic pressure.
Short ranged van der Waals forces are considered in the simulation of the EHD pattering process when the film thickness is less than 100 nm.The van der Waals forces depend on the properties and surface energy of materials used for electrodes, films and the bounding fluid.In the present work, the lower electrode surface (apolar substrate) has a higher energy compared to the liquid film leading to a spinodally stable film (negative values for the pressure, shown in Fig. 4(d)) for the initial stage of the EHD patterning process.
Using a patterned electrode applies a heterogeneous pressure at the film interface leading to replication of the upper electrode pattern onto the polymer film (Fig. 4(a and b)).The distribution of scaled electrostatic pressure (F ES ) acting on the IL film at the initial stage (i.e.flat interface) is compared with a PD film as shown in Fig. 4(e) and (f).A smaller electrode distance results in a higher electrostatic pressure in both ILs and PDs; however, the IL films experience higher electrostatic pressure compared to the PD films in the protrusion regions (38 times higher in the case considered here).

Morphologies in the heterogeneous electric field
The effect of a heterogeneous electric field on the EHD patterning of PD and IL films using strip like and block protrusions is  discussed.Effects of the size and shape of the patterns and the initial filling ratio on the EHD pattering process and the morphology of structures are investigated.
3.3.1 Strip-like protrusions.The 3-D snapshots of the thin IL film interface are used to show the change in the structure height over time under a heterogeneous electric field (Fig. 5).The mean initial film thickness, h 0 , is 30 nm and electrode distances d 1 and d 2 are 188 and 100 nm, respectively.Protrusion width W 2 is 0.64 mm and the groove width W 1 is 2.46 mm.The fluid interface height increases toward the upper electrode with the same rate as it decreases downward during the initial linear stage (i) of the process (see Fig. 5(b)).This growth in height results in the formation of undulating 2-D structures mainly under the protrusion region (see Fig. 5c(i)) where the ridges are surrounded by grooves.In the next stage (ii), fragmentation occurs and initially raises cone shaped structures which occurs randomly along the ridges (see Fig. 5c(ii)).Next the raised structures bridge the upper and lower electrodes and their cross-section enlarges to form raised columnar structures (pillars) (see Fig. 5c(iii)).A further increase of the size of pillars due to the collision of neighboring pillars and their coalescence leads to the formation of larger size pillars (see Fig. 5c(iv)).The width of the pillars in this case is less than 0.4 mm which is about 38% less than the width of protrusions on the electrode.
For comparison, a 3-D snapshot of the PD film structure formed under similar conditions is shown in Fig. 5(d) which has much fewer pillars with oval cross-section formed under the protrusions.The cross-section height profiles (along the L 1 and L 2 ) provide more information regarding the size and shape of pillars in the PD and IL films (Fig. 5(e and f)).The center-to-center distance of adjacent pillars in the PD and IL films is around 2.6 and 1.2 mm respectively.A 53% reduction of the pillar width in the IL film, compared to the PD film, is achieved.The cross-section height profile along the X direction (L 2 ) for the PD and IL films shows that the pillars are not necessarily aligned in the x direction.Hence the number of pillars and their shape may vary between the protrusion lines.
When a flat electrode is replaced with the patterned one, two characteristic spatial lengths are identified.The first is the spacing distance that pillars form within each protrusions (l 1 ) and the second characteristic length is the distance of pillars formed in two neighboring protrusions (l 2 ).The first one is controlled by the same parameters that affect the EHD patterning process of pattern-free electrodes such as the applied voltage, the film filling ratio and electric permittivity, etc.But the second one can also be controlled by adjusting the groove width.
To investigate the effect of the patterned electrode on the IL films, the width of protrusions is increased from 0.64 to 2.46 mm.Structures formed at the quasi steady state stage (stage (iv) in Fig. 5(b)) of the EHD process are compared in Fig. 6.The height and applied voltage of protrusions are the same as in Fig. 5. Results for the PD and IL thin films with the mean initial thicknesses of 30 nm in Fig. 6(b and c) and 50 nm in Fig. 6(d and e) are presented.The 30 nm thick IL film forms two rows of pillars under the protrusion at the middle whereas in the PD film just one row of pillars is formed.In addition, the number of pillars formed in one row of the IL film is higher than in the PD film.For the 30 nm film thickness using a patterned electrode reduces the size of pillars in the IL films compared to a PD films.As previously mentioned, when a flat electrode is used, increasing the initial IL film thickness resulted in more compact structures with smaller sizes (see Fig. 3(b and c)).However, increasing the initial film thickness for the patterned electrode case results in the formation of roll-like structures.
The effect of reducing the groove width (w 1 ) from 2.46 to 0.64 mm at a constant width of protrusions (2.46 mm) is investigated and the 3-D snapshots of the structures formed in the PD and IL films are shown in Fig. 7.The total number of pillars for the PD film is 12 whereas it exceeds 46 pillars for the IL film under the same conditions.By decreasing the width of grooves from 2.46 to 0.64 mm the undulating patterns are not formed under the grooves.Therefore the groove width should be smaller than the strip width to avoid secondary pattern formation in the low electric field (high gap) regions.Finally, by using an appropriately patterned electrode, smaller size features are generated and the structures are more compact and well ordered.
3.3.2Square block protrusions.Another common shape of protrusions used in EHD pattering is the square block protrusion. 2,3,48The square shape protrusions are used to create a heterogeneous electrostatic force in both x and y directions to control the two length scales.In order to validate our numerical simulation, the EHD patterning of the PD film is simulated for the same experimental conditions reported in the literature. 48 schematic of the electrode pattern and 3-D snapshots of the pattern formed on the PD and IL films at the quasi-steady stage of the pattering process are shown in Fig. 8(a-c).Both experiments 48 and simulation result in the formation of four pillars with identical height under the square protrusions.Based on our simulation, the pillars formed in the IL film are located exactly at the corners whereas for the PD film they are close to the center.The pillar size is reduced by 42% in the IL film as compared to the PD film.Examining Fig. 8(b) and (c) closely, it can be seen that an equal number of pillars are formed in both PD and IL films for square block protrusions.This differs from the case of strip-like protrusions where more pillars were observed in the IL film.A possible explanation for this result may be the lack of adequate polymer liquid fluid to create more The cross-section height profile of the interface for the PD and IL films is compared in Fig. 8(d).The intra-array spacing of l 1 B 5.9 mm and l 2 B 8.3 mm in the PD film match the experimental observation 48 of l 1 B 6.2 mm and l 2 B 8.8 mm with reasonable accuracy.This provides validation of the fidelity and accuracy of simulation results for the prediction of the pattern formed in the EHD process -at least for this case.
In the IL film, the center to center distance of pillars under the block protrusions is increased to l 1 B 9 mm and the distance between the neighbor blocks is decreased to l 2 B 4.5 mm.To create well ordered structures with a higher number of pillars, the electrode distances are reduced to d 1 = 100 nm, d 2 = 80 nm and an electrode mask with square protrusions (W 1 = W 2 = 0.84 mm) is used.A schematic of the upper electrode and spatiotemporal evolution of PD and IL films are shown in Fig. 9. Tracking the interface height over time (Fig. 9b(i)), it is found that the  pattern formation includes three main stages.First, the interface is pushed toward the upper electrode in regions with a higher electrostatic force (under the protrusions) and castle shaped structures with four corner towers are formed.The cone shaped structures are formed close to the corners of squares as the gradient of electrostatic force is higher at these boundaries (Fig. 9b(ii) and c(ii)).Fluid depletion around boundaries (blue color ring in color print) is the evidence of fluid flow from thinner regions to the thicker ones.The cone shaped structure formed in the IL film is thinner compared to that formed in the PD film and at the corners, whereas in the PD film it is closer to the center.In the second stage, collision and coalescence of neighboring cones in the PD film occurs even before completion of the pillar formation.Hence, in the PD film, only one pillar is formed under the protrusions bridging the two electrodes (Fig. 9b(iii)) while in the IL film, four pillars with separate ring shaped fluid depletion are formed (Fig. 9c(iii)).In the third stage of pattern formation, secondary pillars are formed in the low electric field region (under the grooves) in both PD and IL films (Fig. 9b(iv) and c(iv)).In the IL film the four pillars formed under the protrusions are merged and form a pillar with a larger cross-section (similar to stage ii in the PD film), then the secondary pillars are created.The number of secondary pillars in the IL film is doubled compared to that in the PD film.Reducing the electrode distances d 1 from 188 nm to 100 nm and d 2 from 100 nm to 80 nm results in the formation of secondary pillars under the grooves.This can be used to create a multi-level structure.
For further examining the effect of protrusion size, it is increased to W 2 = 1.26 mm (=1.5 Â l m ) and the groove width is reduced to W 1 = 0.42 mm (=0.5 Â l m ) as illustrated in Fig. 10a(i).As seen in Fig. 9b(iii) and c(iii), using protrusions with the same size as the characteristic wavelength for the growth of instabilities W 2 = l m results in the formation of one (four) pillar(s) under the protrusions in the PD (IL) film.Increasing the protrusion size leads to the formation of four (nine) pillars in the PD (IL) film.The groove width is chosen narrow enough to avoid the formation of secondary pillars (Fig. 10a(ii and iii)).More detailed information regarding the size of pillars and their intra-spacing is shown in Fig. 10a(iv).The center to center distance of pillars formed under the protrusions (l 1 ) is reduced to B0.675 for the PD and 0.43 mm for the IL film.By reducing the groove width, the secondary pattern is no longer formed under the grooves for both PD and IL films.The distances between the neighboring pillars under the grooves (l 2 ) is around 1.0 and 0.78 mm in the PD and IL films, respectively.
The results show that for every n pillars that are formed with a PD film, ffiffi ffi n p þ 1 ð Þ 2 pillars would be created using an IL film under the same conditions.Although the number of pillars formed in the PD is less than the IL film, the size of pillars in the PD film are identical.In the case of the IL film, the pillars located at the center are smaller than those formed at the corners resulting in non-uniformity of pillar sizes.To have a maximum number of pillars with identical shape and size, a new pattern with nine uniform square block protrusions and the same patterning domain of 3.37 Â 3.37 mm 2 as Fig. 9(a) is

Conclusion
EHD patterning using a thin ionic conductive film and patterned electrodes is numerically examined to create compact (submicron level size) and well ordered structures.Dynamics and spatiotemporal evolution of thin PD and IL films, subjected to both homogeneous and heterogeneous electric fields, are investigated.Components of conjoining/disjoining pressure acting at the film interface during the EHD patterning process are compared.Higher electrostatic pressure acts at the interface of the IL films compared to the PD films.However, for both cases the electrostatic pressure distribution pattern over the domain matches the upper electrode pattern.Based on the 2D nonlinear thin film equation, the evolution of a three-dimensional structure is numerically simulated to find the dynamics and the process of pattern formation in the IL films and the results are compared with the PD baseline cases.For the PD base case results, structures formed at the quasi steady stage are in good agreement with the experiments 48 in the literature and found to be within 5%.The effects of protrusion and groove size (width and depth) in strip like and square block patterns and the film initial filling ratio on the size, shape and order of fluid patterns are investigated.When the flat electrode is replaced with the patterned electrode, two characteristic spatial lengths are identified which are controlled by the width of the strip/square block (W 2 ), groove width (W 1 ) and groove depth (d 1 À d 2 ).Results show that depending on the groove width and depth the EHD patterning process has either two or three main stages.In the two-stage process, patterns only form under the protrusions whereas in the three-stage one, patterns (either undulating or secondary pillars) also form under the grooves.The two-stage process can be used for the pattern replication of mask onto the film and the threestage one is utilized when the multilevel structure formation is desired.Furthermore, more pillars are found to be created using the IL film rather than the PD film.By using electrodes with square The key insight our electrostatic model provides is that, the leaky dielectric model may not account for phenomena important for EHD patterning, namely the presence of finite ions in the conducting layers.The Debye length is defined using ionic strengths, which provides a realistic range of free charge concentrations in the conducting fluid layer compared to the previous work in the literature. 23We expect that the results presented here will be useful for the creation of compact submicron sized features for optoelectronic-related applications.

Fig. 1 2
Fig. 1 2-D Schematic representation of (a) thin liquid film sandwiched between two electrodes and a 3-D schematic of the upper electrode with (b) strip-like (c) square block protrusions.l is the wavelength of growing instabilities at the interface.
r kðd À hÞ sinhðkhÞ (14) More detailed information about the derivation of electrostatic pressure is presented in ref. 15, 16, 23.The well known electrostatic pressure acting on the PD film is duplicated for the sake of completeness 2 as domain with a length of 4l m under periodic boundary conditions is used.Term l m is the characteristic wavelength based on the lower electrode distance.Uniform Cartesian grids with four different grid sizes (n x Â n y ) of 91 Â 91, 121 Â 121, 151 Â and 171 Â 171 are examined to find the optimum grid size and the results are shown in Fig.2.The grid size of 121 Â 121 grid points is found to be an acceptable compromise between the accurate interface height growth and the computational cost.Volume conservative random perturbations are used to initiate the film height of the simulations.

Fig. 2
Fig.2Effect of grid size on the nondimensional maximum and minimum interface height variations versus nondimensional time.The PD film with a mean initial thickness of h 0 = 30 nm, V = 20 V.

Fig. 3
Fig. 3 Total number of pillars versus initial filling ratio of the film.3-D snapshots of the patterns formed in PD and IL films with h 0 = 30 and 70 nm.d 1 = d 2 = 100 nm, V = 20 V.

Fig. 6
Fig. 6 (a) 3-D schematic of upper electrode pattern and 3-D snapshots of the interface height of (b and d) PD films, (c and e) IL films.Mean initial film thicknesses, h 0 are (b and c) 30 and (d and e) 50 nm.d 1 = 188 nm, d 2 = 100 nm, W 1 = W 2 = 2.46 mm and V = 20 V.

Fig. 9
Fig. 9 (a) 3-D schematic of the upper electrode pattern and the interface height profile of (b) PD and (c) IL films.(i) Maximum and minimum interface height profile over time and (ii-iv) 3-D snapshots of the formed structure on the thin film from stage (i) to (iii).Nondimensional times are: T = b(i) 5100, b(ii) 20 550, b(iii) 60 021 and c(i) 1000, c(ii) 2000, c(iii) 3000.Mean initial film thickness h 0 is 30 nm. d 1 = 100 nm, d 2 = 80 nm, W 1 = W 2 = 0.84 mm and V = 20 V.

2
pillars are formed in the IL film for every n pillars generated in the PD film.In the PD film, electrode square block mask protrusions with a width of W 2 = 1.5 Â l m and a groove width of W 1 = 0.5 Â l m and in the IL film square block protrusions with a width of W 2 = l m and a groove width of W 1 = 0.25 Â l m result in compact and well ordered fluid pillars.

Table 1
Constants or parameters used in simulations