Research on Hydrophobic Properties of Grating Structure on Monocrystalline Silicon Fabricated Using Micromachining

In this paper, micromilling was used to process the surface of single crystal silicon, and six different structural parameters of the grating structure were tested to get the contact angle in a different wettability. )e contact angle obtained by the experiment was compared with the theoretical values of theWenzel and Cassie model optimized based on the characteristics of micromilling. )e relationship between structural parameters and wettability was verified. First, the contact angle of the parallel grating structure was greater than the vertical direction, which was influenced by the solid-liquid interface tension. Second, the hydrophobicity of the specimen was in good agreement with the predicted trend of the optimized prediction model C when the width of the convex is reduced. )e support of the theoretical model to the experimental results is instructive to the construction of the structure. In addition, the molecular dynamics were used to verify the hydrophobicity of grating structures from a molecular point of view.


Introduction
e wetting property of any solid material's surface refers to the spreading ability of a given liquid on the solid surface, which is characterized by the contact angle.e contact angle at the solid-liquid interface is a result of the balance between the cohesion and adhesion forces acting at the solid-liquid-gas three-phase contact line [1][2][3].It is also an important technical index of hydrophobic performance.In general, the surface is called hydrophobic when it possesses water contact angle (WCA) equal or higher than 90 °, while it is considered hydrophilic if the WCA is below 90 °. ere are two cases of extreme wetting and nonwetting, where the WCA can approach 0 °or 180 °, respectively, called superhydrophilicity and superhydrophobicity [4][5][6][7].
e hydrophobic surface of monocrystalline silicon has great application value in solar energy for its heat transfer performance, which brings great convenience to people's production and life.In addition, hydrophobic surfaces have important application value in many industrial fields, including self-cleaning surfaces, anti-icing/ frosting systems, and anti-biofouling coatings, which are of great value to the development of society [8][9][10].e multilevel structure and low surface energy materials, composed of micro-and nanoscale structures, can make the surface of any material hydrophobic or superhydrophobic [11][12][13].
Some materials, such as quartz glass, monocrystalline silicon, and diamond, are widely used in high-tech fields such as civil and national defense for their high hardness, high strength, and stable chemical properties in the context of normal temperature.Among them, the monocrystalline silicon has the unique advantages of high thermal conductivity, low density, low expansion coefficient, and excellent optical performance.It also plays an important role in solar energy technology as the monocrystalline silicon is the main material for the production of solar cells [14].
At present, chemical and physical treatment methods account for the majority of the various methods for preparing hydrophobic surfaces.Razavi et al. [15] used oxidation to make microstructures on copper.Wang et al. [16] used ion etching to make microstructures for Si.However, neither of them can control the structural morphology.Chen et al. [17] developed a rapid one-step electrodepositing process to fabricate the superhydrophobic cathodic surface on the copper plate.Liu et al. [18] used the electrodeposition process to fabricate the hydrophobic surface on Cu.Yang et al. [19] and Gurav et al. [20] developed some new coatings using a sol-gel method.Zorba et al. [21,22] used femtosecond laser to tailor the wetting response of silicon surfaces.However, the parameters of those physically and chemically patterned structures are not controllable, and the ambient environments of the experiment are very specific [23].As a result, it is difficult to study the rules of variation between surface structure and wettability accurately.e processing is very expensive, which means the large-scale production is difficult [24].At the same time, the size of the milling tool is reduced gradually with the study of the micromachining mechanism.It is possible to study the wettability of micron-sized microstructures by micromilling [25].e controllability and repeatability of the machining process is good.Besides, the high-speed micromilling technology is widely used in machinery manufacturing industry.is technology is known for high processing precision, high processing efficiency, simple preparation process, and relatively low preparation cost and arbitrary shape of three-dimensional structure processing [26].In this context, the structure parameters and morphology of any surface can be determined accurately with micromilling technology. is processing method ensures to study the rules of variation between surface structure and wettability accurately.
e processing is also simplified.Moreover, the theoretical models cannot be applied to all structures, which resulted in a large deviation between theoretical and experimental values of contact angle [27].
In this paper, the structure of the monocrystalline silicon surface was changed and the appropriate grating structure was selected to obtain hydrophobicity, which improves the thermal conductivity and solar energy efficiency and enhances the anticorrosion ability [28].First, a theoretical model suitable for micromilling characteristics was established to study the effect of micron particle size on the contact angle of droplets resting on the silicon surfaces.Second, the micrograting structure was fabricated by high-speed micromilling [29,30].ird, the contact angle measured after machining was compared with the analysis under molecular dynamics simulation, and the influence of microstructure size parameters on the contact angle was studied.

Basic Principle of Hydrophobicity.
e wetting phenomenon reflects the equilibrium relation of three-phase interfacial tension between solid, liquid, and gas.In 1805, the functional correspondence between the static contact angle of a smooth surface and the interfacial tension is derived for the first time by Young [31], which is based on the thermodynamic equilibrium theory.And this functional relationship is named "Young's equation."Young's equation is proposed under the condition of neglecting the line tension and the gravity of water droplets.It is used to measure the ideal contact angle of the absolutely smooth surface.Contact angle θ with a smooth surface is calculated as follows [32]: where σ sv , σ sl , and σ lv are the interfacial tensions of the solidgas phase, solid-liquid phase, and gas-liquid phase, respectively, which are shown in Figure 1.
Actually, the surface of natural solids has a certain degree of roughness.e Wenzel [33] model and the Cassie-Baxter [34] model are both proposed to estimate the effect of surface roughness.
e Wenzel model shows the case that the droplet maintains contact with the entirety of the rough surface, which results in the increase of the interfacial contact area.In contrast, the Cassie model indicates the case that the droplet skips between the peaks of the roughness and leaves the interspace of solid-liquid and liquid-gas interfaces below it.e structures are shown in Figure 2.
e contact angle θ w of the Wenzel model is calculated as follows: e contact angle θ c of the Cassie-Baxter model is calculated as follows: In formulas ( 2) and ( 3), r and f are the ratio of the actual solid-liquid contact area to the horizontal projection area of the droplet, which was named after roughness factors [35].
e Cassie model plays an important role in enhancing the hydrophobic properties of the microstructure surface.As a result, the Cassie model is an ideal model to describe the superhydrophobicity of the solid surfaces towards the liquid (water) [36].

Construction of Prediction Model Based on Parameters.
e schematic diagram of the grating structure is shown in Figure 3, where a is the convex platform width of the grating structure, b represents the spacing width of the grating structure, h is the height of convex, and l represents the length of the contact surface of the droplet on the surface.
e precise prediction model suitable for micromilling characteristics was established on the basis of Gibbs free energy.
e surface area of the convex platform under droplets is Gibbs free energy can be associated with the contact angle due to the interfacial free energy, and surface tension is equal in value.In this paper, the prediction models are named prediction model W and prediction model C, respectively.
e contact areas of solid-liquid, gas-liquid, and solidgas interfaces of prediction model W are as follows: 2 Advances in Materials Science and Engineering According to equation (1), the relationship between surface tension and contact areas can be expressed as follows: erefore, the relationship between contact angle and geometric parameters in prediction model W can be shown as follows: Wenzel model Cassie model

Advances in Materials Science and Engineering
In contrast, the contact areas of prediction model C are different from prediction model W: us, the relationship between surface tension and the contact areas of prediction model C can be represented as follows: According to equation ( 8), the relationship between the contact angle and geometric parameters can be obtained as follows: Figure 4 shows the relationship between structural parameters and contact angle, which were derived from equations ( 6) and (9).
It can be seen that when the height and spacing width of the convex are constant, the contact angle of prediction model W increases with the increase of convex platform width, while the contact angle of prediction model C decreases.Moreover, when the height and platform width of convex are fixed, the contact angle of prediction model W goes up and so does of prediction model C.

Experimental Method
e experimental material was selected to be polished monocrystalline silicon with 100 crystal surface.
e diameter of it was 25 mm, and the thickness was 400 μm [37].However, the silicon wafer was too thin to fix it on the worktable by the traditional clamping method.As a result, a copper with the size of 25 mm * 25 mm * 25 mm was fixed on the vise.
en, the upper surface of the copper was smoothed with a 3 mm diameter milling cutter.At last, the silicon wafer was pasted on the copper smoothly and tightly.
e workpieces were machined on high-speed milling center, as shown in Figure 5, with a double-edged micromilling cutter.e milling cutter was a diameter of 0.2 mm with diamond coating.e experiments were carried out under the condition that the spindle speed was 48000 r/min and constant feed speed was 6 mm/min [38].15 grating grooves of different sizes were processed on the intermediate surface of the workpieces to form the grating structure.
After processing, these samples were cleaned in an ultrasonic cleaner containing ethanol for 5 min to remove the burr and impurities, followed by a 10% HF aqueous treatment to remove the oxide grown on the surface of these samples.en, the samples were air-dried naturally to reduce the external errors for the measurement of contact angle.All the apparent morphology of samples was observed with VHX-5000 Superhigh magnification lens zoom 3D microscope (VHX-5000, Keyence, Japan).e contact angle of the grating structure was measured by an optical contact angle measuring instrument (OCA, DataPhysics, Germany).
e middle of the grating structure was selected as the measuring point of the contact angle to ensure that the experimental data have the same measuring point.e contact angle of each specimen was measured five times to ensure the reliability of these data [39].

Result and Analysis
During processing, both the spacing width b and the height of convex h were fixed.e geometrical parameters of the design of the grating structure are shown in Table 1. Figure 6 shows the obtained grating structure surface for Si.
Figure 7 shows the obtained contact angles for the sample.
e contact angle of the unprocessed monocrystalline silicon is 56.3 °.Besides, the theoretical contact angle is calculated based on equations ( 6) and ( 9).All of the data are shown in Table 2.
According to Table 2, the variation curve of each contact angle is shown in Figure 8. Comparing Table 2 and Figure 8, it can be seen that the contact angle of the six samples is greatly improved than the original specimen.
From Figure 8, it can be seen that a huge mismatch in the value of the contact angle between the vertical direction and the parallel direction.In these six specimens, the contact angle in the parallel direction is always greater than the vertical direction, and the average difference is 51.85 °.It is clear that the wetting phenomenon is related to the three-phase equilibrium of the solid, liquid, and gas.In this paper, there is no convex platform constraining the droplets on both sides [40].
e solid-liquid interfacial tension is lower in the direction vertical to the groove than the parallel direction.e droplet will tend to spread completely on the solid in the vertical direction as (σ sv − σ sl ) is larger than σ lv .As a result, the hydrophobic anisotropy will occur in the grating structure.e sphere is easier to form in the parallel direction, which will cause an increase in the contact angle.
e trend of the contact angle in the parallel direction was in good agreement with the value of prediction model C. Figure 7(b) shows that there is a gap between the droplet and the structure of surface, which means this structure matches prediction model C [41].
e established model predicts the contact angle accurately.In prediction model C, the contact angle decreased with the increase of the convex platform width and so does the contact angle of the parallel direction.
It can be seen from Figure 8 that the contact angle in two directions increases as the width of the convex changes from 20 μm to 60 μm.Among these specimens, the hydrophobic property of sample 1 was the best.erefore, the smaller the convex platform width, the larger the contact angle, and the hydrophobic performance of the structure was stronger.

4
Advances in Materials Science and Engineering Based on the analysis above, the contact angle of sample 1 was the largest and the hydrophobicity was the best.e further analysis on sample 1 was carried out.
Figure 9 shows the surface microstructure and size of sample 1, which was observed under 200 times of the VHX-5000 ultradeep 3D microscope.Comparing Figures 9 and 10, it can see that there are some burrs on the edge of sample 1.As the tensile stress of monocrystalline silicon is much greater than its shear stress, curling is more easily to occur in the milling process of monocrystalline silicon.
en, the cracks and breakage will result from it.As a result, the experimental value of the contact angle is larger than the value of the theoretical model [42].However, the overall gap is not so big, and the average D value is 4.3 °. e overall contact angles achieve the desired results.

Model Establishment.
Monocrystalline silicon is an octagonal diamond crystal structure with anisotropy.Combining with experimental, the 100 plane of single crystal silicon was selected for MD simulation.e grating micro-nanostructure model was constructed by removing the crystal cell of the model from the built model.Figure 11 is the simulation model in this experiment.e upper water molecules were arranged evenly, and the lower silicon took four atomic layer thicknesses.
In this simulation, the initial temperature of the system was 298 K, and the integration step was 1.0 fs. e periodic boundary conditions were applied in the X and Y directions of the simulation space.Besides, the direction of Z used the specular reflection boundary condition, and the simulation time was 800 ps [43].e SPC model was used to simulate the droplets.Furthermore, the Lennard-Jones potential was used between the oxygen atoms, and the Tersoff potential was used to describe Si and Si in single crystal silicon.Lennard-Jones potential was used to describe the interaction between droplet and silicon substrate.Lennard-Jones potential energy was calculated as follows [44,45]: where ε and σ are energy parameters and length parameters, respectively.Advances in Materials Science and Engineering e calculation formula of the Tersoff potential function for monocrystalline silicon was as follows: where u ij is the potential function between atomic i and atomic j, f C (r ij ) is the truncated function, f R (r ij ) is the repulsive function, b ij is the low order function, f A (r ij ) is the attraction function, and r ij is the distance between atomic i and atomic j.       e hydrophobicity of the grating structure was verified by the molecular structure.

Conclusion
(1) A theoretical model of contact angle suitable for milling characteristics was established in this paper.e characteristics of the grating structure with micromilling were considered, and the accuracy of the theoretical model was guaranteed, which makes the experimental data more reliable.
(2) e contact angle in the parallel direction was always greater than that in the vertical direction.e droplets were mutually constrained in two   Advances in Materials Science and Engineering directions, making the hydrophobic state more stable.e hydrophobic anisotropy of the grating structure can be beneficial to the special wetting conditions.
(3) e relationship between the contact angle and the size parameter of the grating structure was determined based on experimental data.For the grating structure with constant spacing width and convex height, the smaller the convex platform width, the larger the contact angle.e grating structure designed in this paper is simple, which can also ensure the improvement of the hydrophobicity of the material.ere is good guidance for the construction of the hydrophobic surface on the grating structure.(4) e difference between the experimental value and the theoretical value of the contact angle was due to the existence of some collapse edges on the edge of the convex.Overall, the contact angles of the grating structure achieved the desired results.

Data Availability
e data used to support the findings of this study were supplied by Ziyang Cao under license and so cannot be made freely available.Requests for access to these data should be made to the corresponding author upon request via mail.

Figure
Figure Variation curve of contact angle between measurement and theory.

6
Advances in Materials Science and Engineering e potential energy parameters of the SPC model are shown in Table3[46].
Figure 12(a)  shows the droplet state under the surface of normal monocrystalline silicon.According to the previous simulation, it can be seen that the contact angle is 87 °.
Figure 12(b) shows the water droplet state under the grating structure.It is obvious that the contact angle of the grating structure is 124.4 °.

Table 1 :
Processing parameters of the grating structure.

Table 2 :
eoretical and practical values of contact angle.

Table 3 :
Parameter values of the SPC model.