Wetting characterisation of silicon (1,0,0) surface

Line tension from dual-geometry sessile droplet measurements: Combining contact-angle size-dependence data for axisymmetric and cylindrical droplets to determine the line tension

We perform a thermodynamic analysis of various contributions to the size dependence of the contact-angle cosine for both axisymmetric and cylindrical sessile droplets. This shows that a commonly used method to determine the line tension from the slope of the contact-angle cosine dependence on the three-phase contact-line curvature (for axisymmetric droplets) provides a certain combination of the line tension, the adsorptions at the three interfaces, and the macroscopic contact angle. To extract the contribution related to the contact line in the leading order and determine the line tension, we propose a simple technique using droplet-size dependences of the contact angle for axisymmetric and cylindrical droplets under the same conditions.

Accurate Force Field for Carbon Dioxide-Silica Interactions Based on Density Functional Theory

Fluid-silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional "mixing rules" to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO2 in silica slit pores are also examined.

Effect of molecular branching and surface wettability on solid-liquid surface tension and line-tension of liquid alkane surface nanodroplets

HYPOTHESIS

Surface nanodroplets have important technological applications. Previous experiments and simulations have shown that their contact angle deviates from Young's equation. A modified version of Young's equation considering the three-phase line tension (τ) has been widely used in literature, and a wide range of values for τ are reported. We have recently shown that molecular branching affects the liquid-vapour surface tension γlv of liquid alkanes. Therefore, the wetting behaviour of surface nanodroplets should be affected by molecular branching. This study conducted molecular dynamics (MD) simulations to gain insight into the wetting behaviour of linear and branched alkane nanodroplets on oleophilic and oleophobic surfaces. We aim to examine the Young equation's validity and branching's effect on fundamental properties, including solid-liquid surface tension γsl and line tension τ.

SIMULATIONS

The simulations were performed on a linear alkane, triacontane (C30H62), as well as four of its branched isomers: 2,6,13,17-tetrapropyloctadecane,2,6,9,10,13,17-hexaethyloctadecane, 2,5,7,8,11,12,15-heptaethylhexadecane and 2,3,6,7,10,11-hexapropyldodecane. Nanodroplets with a diameter of approximately 15 nm were released onto the surfaces, and their contact angles were measured. Additionally, using a novel approach, the solid-liquid surface tension (γsl), the validity of Young's equation and line tension for all alkane and surface combinations are determined.

FINDINGS

It was discovered that the calculated γsl, deviated from the theoretical γsl,Young predicted from Young's equation for all alkanes on oleophilic surfaces. However, this deviation was minimal for branched alkanes on the oleophobic surfaces but more significant for the linear alkane. The findings indicated that γsl < 0 for oleophilic surfaces and γsl > 0 for oleophobic surfaces. Moreover, it was observed that |γsl| was lower for branched molecules and decreased as branching increased. Line tension values were then determined through a novel method, showing τ was positive for oleophilic surfaces ranging from 1.30 × 10-10 to 6.27 × 10-11N. On an oleophobic surface, linear alkane shows a negative line tension of -1.15 × 10-10N and branched alkanes up to two orders of magnitude lower values ranging from -2.09 × 10-12 to 2.43 × 10-11N. Line tension values between -1.15 × 10-10 and + 1.1 × 10-10N are calculated for various linear alkane and surface combinations. These findings show the dependence of line tension on the contact angle and branching, demonstrating that for linear alkanes, τ is significant, whereas, for branched alkanes, line tension is smaller or negligible for large contact angles.

Effect of interfacial vibrational coupling on surface wettability and water transport

We report that the atomic-scale vibrational coupling at the solid-fluid interface can substantially alter the interfacial properties such as wettability and fluid slip. The wettability of water droplets on substrates subjected to various vibrational frequencies is studied using molecular dynamics simulation. The contact angle increases (i.e., becomes more hydrophobic) when the oscillation frequency of the substrate matches the intermolecular bending frequency of liquid water. We investigate the underlying mechanism by examining the dynamics of water molecules at the interface and find that the temporal contact between the solid and fluid is shorter when the frequencies match, resulting in weak solid-fluid adsorption. We further report that the vibrational match at the interface reduces wall-fluid friction and enhances water transport through the nanopore. Our findings demonstrate the importance of the atomic-scale vibrational coupling at the solid-fluid interface on the physicochemical behavior of nanodevices and biological nanochannels.

Molecular Dynamics Simulation of Nanodroplets Impacting Stripe-Textured Surfaces

The dynamic behavior of droplets impacting on textured surfaces has an important influence on many engineering applications, such as anti-icing and self-cleaning. However, the mechanism and law of the effect of textured surfaces on the impact behavior of nanodroplets has not been fully revealed yet. In this paper, the molecular dynamics (MD) method is used to model the dynamic behavior of nanodroplets after impacting the solid surface with a striped texture. The influences of texture gap and texture angle on the real contact area, spreading factor, contact time, and bounce velocity of the droplet after impact are also quantitatively analyzed. It is shown that the striped texture produces significant anisotropy in the spreading and contraction behavior of nanodroplets after impact, and the anisotropy is more pronounced on the ridged texture surface than on the grooved texture surface. In addition, we find that the texture gap has little effect on the dynamic behavior of nanodroplets impacting the textured surface. However, as the bottom angle of the texture increases, the real contact area and bounce velocity of the nanodroplet increase significantly, while the contact time and spreading factor decrease. This work further elucidates the characteristics and mechanisms of nanodroplets impacting on stripe-textured surfaces and provides a theoretical basis for the design of nanostructured surfaces in relevant applications.

Energy-Based Interface Detection for Phase Change Processes of Monatomic Fluids in Nanoconfinements

An energy-based liquid-vapor interface detection method is presented using molecular dynamics simulations of liquid menisci confined between two parallel plates under equilibrium and evaporation/condensation conditions. This method defines the liquid-vapor interface at the location where the average kinetic energy of atoms first exceeds the average potential energy imposed by all neighboring molecules. This definition naturally adapts to the location of the menisci relative to the walls and can properly model the behavior of the liquid adsorbed layers. Unlike the density cutoff methods frequently used in the literature that suffer from density layering effects, this new method gives smooth and continuous liquid-vapor interfaces in nanoconfinements. Surface tension values calculated from the equilibrium MD simulations match the Young-Laplace equation better when using the radius of curvatures calculated from this method. Overall, this energy-based liquid-vapor interface detection method can be used in studies of nanoscale phase change processes and other relevant applications.

Deposition and Characterization of RP-ALD SiO2 Thin Films with Different Oxygen Plasma Powers

In this study, silicon oxide (SiO2) films were deposited by remote plasma atomic layer deposition with Bis(diethylamino)silane (BDEAS) and an oxygen/argon mixture as the precursors. Oxygen plasma powers play a key role in the quality of SiO2 films. Post-annealing was performed in the air at different temperatures for 1 h. The effects of oxygen plasma powers from 1000 W to 3000 W on the properties of the SiO2 thin films were investigated. The experimental results demonstrated that the SiO2 thin film growth per cycle was greatly affected by the O2 plasma power. Atomic force microscope (AFM) and conductive AFM tests show that the surface of the SiO2 thin films, with different O2 plasma powers, is relatively smooth and the films all present favorable insulation properties. The water contact angle (WCA) of the SiO2 thin film deposited at the power of 1500 W is higher than that of other WCAs of SiO2 films deposited at other plasma powers, indicating that it is less hydrophilic. This phenomenon is more likely to be associated with a smaller bonding energy, which is consistent with the result obtained by Fourier transformation infrared spectroscopy. In addition, the influence of post-annealing temperature on the quality of the SiO2 thin films was also investigated. As the annealing temperature increases, the SiO2 thin film becomes denser, leading to a higher refractive index and a lower etch rate.

Size dependent influence of contact line pinning on wetting of nano-textured/patterned silica surfaces

Wetting behavior on a heterogeneous surface undergoes contact angle hysteresis as the droplet stabilized at a metastable state with a contact angle significantly different from its equilibrium value due to contact line pinning. However, there is a lack of consensus on how to calculate the influence of pinning forces. In general, the pinning effect can be characterized as (i) microscopic behavior when a droplet is pinned and the contact angle increases/decreases as the droplet volume increases/decreases and (ii) macroscopic behavior as the pinning effects decrease and ultimately, disappear with the increase of the droplet size. The current work studied both behaviors using molecular dynamics (MD) simulation with more than 300 different size water droplets on silica surfaces with three different patterns across two different wetting conditions. Results showed that the contact angle increases linearly with increasing droplet volume through the microscopic behavior, while the droplet is pinned on top of a certain number of patterns. When we normalized the droplet size with the corresponding pattern size, we observed a "wetting similarity" that linear microscopic contact angle variations over different size heterogeneities continuously line up. This shows that the pinning force remains constant and the resulting pinning effects are scalable by the size ratio between the droplet and pattern, independent of the size-scale. The slope of these microscopic linear variations decreases with an increase in the droplet size as observed through the macroscopic behavior. We further found a universal behavior in the variation of the corresponding pinning forces, independent of the wetting condition. In macroscopic behavior, pinning effects become negligible and the contact angle reaches the equilibrium value of the corresponding surface when the diameter of the free-standing droplet is approximately equal to 24 times the size of the surface structure. We found that the pinning effect is scalable with the droplet volume, not the size of the droplet base.

Molecular transportation phenomena of simple liquids through a nanoporous graphene membrane

The mechanism of pressure-driven transport of simple liquid through a nanoporous graphene membrane has been analyzed using nonequilibrium molecular dynamics simulation. In this study, we investigate liquid dynamics properties such as local density, pressure variation, and local viscosity depending on the flow region. With movement of the specular reflection wall at the end of the front and back reservoirs, a pressure difference occurs mainly due to the change in the relative distance between the liquid molecules in the corresponding reservoir. The interfacial pressure difference strongly depends on the intermolecular force of the graphene membrane governed by the layered structure of the simple liquid and the applied flow velocity. The local viscosity was calculated for a nanochannel of simple liquid sheared by graphene walls. The liquid velocity adjacent to the pore edge was considered as the slip velocity, which provides updates in the Sampson flow equation. We observed that the entrance interfacial pressure and higher local viscosity in the vicinity of the graphene membrane, which are associated with the optimized definition of the wall-liquid boundary near the pore edge, play a critical role in the permeation of simple liquids through the nanoporous graphene membrane.

Slip Effects on Ionic Current of Viscoelectric Electroviscous Flows through Different Length Nanofluidic Channels

The pressure driven slip flow of an electrolyte solution is studied through different nanofluidic channel lengths at varying salt concentrations. The viscous-thickening due to the electrostatic interactions within the electric double layer and the reverse ionic transport due to the streaming potential are developed. The influence of the Navier slip boundary condition is described under both electroviscous and viscoelectric effects with a surface charge regulation (CR) model while the observed behavior is compared and validated with molecular dynamic (MD) calculations from multiple studies. Results show that electroviscous and viscoelectric effects decrease transport. Earlier studies at the no slip boundary presented an increase of ionic current by increasing salt concentration and decreasing channel length. In contrast, our study found that the ionic current occurred almost independent of both salt concentration and channel length, except for very short channels and very low salt concentrations, when electroviscous and viscoelectric effects were considered. In the case of the constant slip length condition, ionic conduction was enhanced, but velocity slip developing on surfaces showed significant variation based on the salt concentration and channel length. This is due to the natural CR behavior enhancing the surface charge and consequential near surface electrohydrodynamics as a result of increase in salt concentration and/or decrease of channel length. Considering that the electroviscous effect alone creates up to 70% lower velocity slips than Poiseuille flow predictions, while further including the viscoelectric effect, results in an almost no-slip condition at high salt concentrations and/or short channels. As a result, the ionic current of a viscoelectric electroviscous slip flow is found to be equal to 1/3 of an electroviscous slip flow and to decrease with a decrease in the channel length.

Molecular insight into the formation of adsorption clusters based on the zeta isotherm

This work presents a series of molecular dynamics simulations of argon adsorption on a silicon substrate with different lattice orientations. From the simulation results, the density profiles are discussed and the amount of adsorbed particles is obtained at different pressures. It is found that the solid surface orientation has a great influence on the density distributions and atomic arrangements near the surface. With the collected data, the thermal constants derived from the expression of zeta adsorption isotherms are determined. The calculated isotherms agree well with the simulation results. Also, from a microscopic point of view, the molecular insights show that the structures of the adsorbates are present as clusters with different numbers of particles. The size of the clusters changes with pressure. At a relatively small pressure ratio, most of the clusters consist of a single molecule. As the pressure ratio increases, larger sized clusters appear, forming various cluster-types. The molecular cluster distributions are closely consistent with the basic approximation of the zeta adsorption isotherm. Furthermore, the surface adsorption sites determined from molecular dynamics simulation show good agreement with that predicted by the zeta isotherm model, which reaffirms the effectiveness of the theoretical model. When the isotherm is extended to a pressure ratio greater than unity, a finite amount of adsorption is predicted and the wetting conditions are obtained. Affected by the solid surface orientations, the pressure ratio at wetting for the silicon substrate with the (111) surface plane is larger than those of the (100) and (110) surfaces, indicating that a higher subcooling is required for the wetting transition.

What Is the Value of Water Contact Angle on Silicon?

Silicon is a widely applied material and the wetting of silicon surface is an important phenomenon. However, contradictions in the literature appear considering the value of the water contact angle (WCA). The purpose of this study is to present a holistic experimental and theoretical approach to the WCA determination. To do this, we checked the chemical composition of the silicon (1,0,0) surface by using the X-ray photoelectron spectroscopy (XPS) method, and next this surface was purified using different cleaning methods. As it was proved that airborne hydrocarbons change a solid wetting properties the WCA values were measured in hydrocarbons atmosphere. Next, molecular dynamics (MD) simulations were performed to determine the mechanism of wetting in this atmosphere and to propose the force field parameters for silica wetting simulation. It is concluded that the best method of surface cleaning is the solvent-reinforced de Gennes method, and the WCA value of silicon covered by SiO₂ layer is equal to 20.7° (at room temperature). MD simulation results show that the mechanism of pure silicon wetting is similar to that reported for graphene, and the mechanism of silicon covered by SiO₂ layer wetting is similar to this observed recently for a MOF.

Preparation of Si-SiOx nanoparticles from volatile residue produced by refining of silicon

Residual Si was produced on a furnace wall when upgraded metallurgical grade Si was refined by electron beam melting. It was then recycled to prepare Si-SiO_x nanoparticles with an average size of 100 nm by planetary ball milling. The obtained Si-SiO_x nanoparticles mainly consist of amorphous Si, crystalline Si and amorphous SiO_x, which was confirmed by XRD, FTIR, XPS and TEM. SiO_x is mainly composed of SiO₂ and SiO_1.35. Distilled water used as a grinding aid not only enhances milling efficiency, but also plays a key role in obtaining SiO_x. During refining of upgraded metallurgical grade Si, the deposition pattern of residual Si on furnace wall agrees with model of three-dimension growth. Growth of Si-SiO_x nanoparticles is the mutual effect of distilled water and ball milling. Si-SiO_x nanoparticles were doped into phenolic resin pyrolysis carbon as anode materials for lithium ion batteries, and 10% doping was observed to improve the specific capacity. After 500 cycles, specific capacity of delithiation remained around 550 mA h/g. It suggests the residual Si is a value-added by-product, and it can be recycled as anode materials for lithium ion batteries.

Manipulating thermal resistance at the solid–fluid interface through monolayer deposition

Heat transfer across an interface between a monolayer coated solid substrate and fluid has been extensively analyzed through a series of non-equilibrium molecular dynamics simulations. The effect of the monolayer was studied by varying its atomic mass (m_M) and interaction energy between monolayer particles (ε_MM). Even though the fluid adsorption plays a role in heat transfer at the solid–fluid interface, we found that the interfacial thermal resistance (Kapitza resistance) is highly affected by the insertion of monolayer without any further change in the liquid structure near the solid surface. The Kapitza length monotonically increases with the increase of m_M irrespective of ε_MM. The observations were explained by analysis of the overlap of the phonon spectrum at the interface using vibrational density of states. The effect of the monolayer on the Kapitza length was summarized by a fourth polynomial function that demonstrates the contribution of both m_M and ε_MM with respect to each other on the Kapitza resistance within the parametric range studied.

At the interface between monolayer coated solid substrate and fluid, the effect of interfacial mismatch on Kapitza length due to the monolayer particles has been extensively analyzed through a series of non-equilibrium molecular dynamics simulation.

Multiscale Modeling of the Three-Dimensional Meniscus Shape of a Wetting Liquid Film on Micro-/Nanostructured Surfaces

The design of structured surfaces for increasing the heat flux dissipated during boiling and evaporation processes via enhanced liquid rewetting requires prediction of the liquid meniscus shape on these surfaces. In this study, a general continuum model is developed to predict the three-dimensional meniscus shape of liquid films on micro/nanostructured surfaces based on a minimization of the system free energy that includes solid-liquid van der Waals interaction energy, surface energy, and gravitational potential. The continuum model is validated at the nanoscale against molecular dynamics simulations of water films on gold surfaces with pyramidal indentations, and against experimental measurements of water films on silicon V-groove channels at the microscale. The validated model is used to investigate the effect of film thickness and surface structure depth on the meniscus shape. The meniscus is shown to become more conformal with the surface structure as the film thickness decreases and the structure depth increases. Assuming small interface slope and small variation in film thickness, the continuum model can be linearized to obtain an explicit expression for the meniscus shape. The error of this linearized model is quantitatively assessed and shown to increase with increasing structure depth and decreasing structure pitch. The model developed can be used for accurate prediction of three-dimensional meniscus shape on structured surfaces with micro/nano-scale features, which is necessary for determining the liquid delivery rate and heat flux dissipated during thin-film evaporation. The linearized model is useful for rapid prediction of meniscus shape when the structure depth is smaller than or comparable to the liquid film thickness.

Mimicking wettability alterations using temperature gradients for water nanodroplets

A sessile droplet or a film usually moves from hotter regions to colder regions, due to variations in interfacial tension. This, known as the so-called Marangoni effect, is true for most pure liquids like water for which the surface tension decreases with an increase in temperature. In stark contrast to this existing understanding, we bring forth the coupled effect of wettability and temperature gradients on the dynamics of the three-phase contact line. By simultaneously tracking the dynamic evolution of the three-phase contact line due to the evaporation and diffusion of molecules through molecular dynamics simulations, we explore the coterminous effects of the change of surface tension coefficients and wetting parameters with temperature on sessile droplets residing on surfaces with different wettabilities. We demonstrate, for the very first time, that the inverse Marangoni effect, which is believed to be exclusively observed in mixtures and self-rewetting fluids, is feasible in pure water at scales where inertial effects are negligible. The results of the study find application in electronic chip cooling where by the combined tuning of surface characteristics and Marangoni forces, droplets can be passively transported to warmer regions for efficient thermal management.

Ion-Specific Effects on the Elongation Dynamics of a Nanosized Water Droplet in Applied Electric Fields

We report an all-atom molecular dynamics study of the structures and dynamics of salty water droplets on a silicon surface under the influence of applied electric field. Our simulation results support ion-specific effects on the elongation dynamics of salty nanodroplets, induced by the field. This feature has not been explored up to now in monovalent salts. Nevertheless, the importance of ion-specific effects is widely confirmed in biological and colloidal systems. In particular, the increase of salt concentration enhances the effect of the nature of ions on the wetting properties of droplets. In the presence of electric field (0.05 V Å^-1), a complete spreading is implemented in a short time for different droplets at a concentration of 1 M, and the droplet morphology is stable, observed at long time scales. However, a higher salt concentration of 4 M considerably suppresses the spreading process owing to the increase of surface tension. It was found that the NaCl droplet shows deformation oscillations along the external field, but cannot fully wet the substrate surface. By contrast, the CsCl droplet reaches complete elongation rapidly and adopts a steady strip shape. The KCl droplet undergoes frequent transitions between breakup and connection. Additionally, the droplets can be elongated only when the electric field strength exceeds a threshold value. The dipole orientation of interfacial water and the ionic diffusion exhibit ion-specific dependences, but the hydrogen bond network is scarcely disturbed, excluding a concentration-dependent effect.

Dynamic Chemically Driven Dewetting, Spreading, and Self-Running of Sessile Droplets on Crystalline Silicon

A chemically driven dewetting effect is demonstrated using sessile droplets of dilute hydrofluoric acid on chemically oxidized silicon wafers. The dewetting occurs as the thin oxide is slowly etched by the droplet and replaced by a hydrogen-terminated surface; the result of this is a gradual increase in the contact angle of the droplet with time. The time-varying work of adhesion is calculated from the time-varying contact angle; this corresponds to the changing chemical nature of the surface during dewetting and can be modeled by the well-known logistic (sigmoid) function often used for the modeling of restricted growth, in this case, the transition from an oxidized surface to a hydrogen-terminated silicon surface. The observation of the time-varying contact angle allows one to both measure the etch rate of the silicon oxide and estimate the hydrogenation rate as a function of HF concentration and wafer type. In addition to this, at a certain HF concentration, a self-running droplet effect is observed. In contrast, on hydrogen-terminated silicon wafers, a chemically induced spreading effect is observed using sessile droplets of nitric acid. The droplet spreading can also be modeled using a logistical function, where the restricted growth is the transition from hydrogen-terminated to a chemically induced oxidized silicon surface. The chemically driven dewetting and spreading observed here add to the methods available to study dynamic wetting (e.g., the moving three-phase contact line) of sessile droplets on surfaces. By slowing down chemical kinetics of the wetting, one is able to record the changing profile of the sessile droplet with time and gather information concerning the time-varying surface chemistry. The data also indicates a chemical interface hysteresis (CIH) that is compared to contact angle hysteresis (CAH). The approach can also be used to study the chemical etching and deposition behavior of thin films using liquids by monitoring the macroscopic droplet profile and relating this to the time-varying physical and chemical interface phenomena.

Mass transport and structure of liquid n-alkane mixtures in the vicinity of α-quartz substrates

The asymmetric self-diffusion characteristics of binary n-alkane mixtures in the first adsorption layer near α-quartz substrate were demonstrated and discussed.

Near-surface viscosity effects on capillary rise of water in nanotubes

In this paper, we present an approach for predicting nanoscale capillary imbibitions using the Lucas-Washburn (LW) theory. Molecular dynamics (MD) simulations were employed to investigate the effects of surface forces on the viscosity of liquid water. This provides an update to the modified LW equation that considered only a nanoscale slip length. An initial water nanodroplet study was performed to properly elucidate the wetting behavior of copper and gold surfaces. Intermolecular interaction strengths between water and corresponding solid surfaces were determined by matching the contact angle values obtained by experimental measurements. The migration of liquid water into copper and gold capillaries was measured by MD simulations and was found to differ from the modified LW equation. We found that the liquid layering in the vicinity of the solid surface induces a higher density and viscosity, leading to a slower MD uptake of fluid into the capillaries than was theoretically predicted. The near-surface viscosity for the nanoscale-confined water was defined and calculated for the thin film of water that was sheared between the two solid surfaces, as the ratio of water shear stress to the applied shear rate. Considering the effects of both the interface viscosity and slip length of the fluid, we successfully predicted the MD-measured fluid rise in the nanotubes.