Murat Barisik

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

Water droplets on bare silicon surfaces are studied to examine the wetting behaviour as a function of the surface energy and to parameterise water–silicon interactions in order to recover the hydrophobic behaviour measured by experiments. Two different wetting regimes characterised by a critical interaction strength value are observed. At a threshold value of the water–silicon interaction parameter, water molecules start penetrating into the first layer of silicon surface under thermally vibrating walls, resulting in two distinct wetting behaviours. Fixed (cold) silicon walls do not exhibit the two different wetting characteristics. Size effects are studied for nano-scale droplets, and line tension influence is observed depending on the surface wettability. Decrease in the droplet size increases the contact angle values for the low wetting cases, while contact angles decrease for smaller droplets on the high wetting surfaces. Considering the line tension effects and droplet size, ϵSi–O for water–silicon interactions to recover the hydrophobic behaviour of silicon surfaces is estimated to be 12.5% of the value predicted using the Lorentz–Berthelot mixing rule.

Boundary treatment effects on molecular dynamics simulations of interface thermal resistance

Molecular Dynamics simulations of heat conduction in liquid Argon confined in Silver nano-channels are performed subject to three different thermal conditions. Particularly, different surface temperatures are imposed on Silver domains using a thermostat in all and limited number of solid layers, resulting in heat flux in the liquid domain. Alternatively, energy is injected and extracted from solid layers to create a NVE liquid Argon system, which corresponds to heat flux specification. Imposition of a constant temperature region in the solid domain results in an unphysical temperature jump, indicating the presence of an artificial thermal resistance induced by the thermostat. Thermal resistance analyses for the components of each case are performed to distinguish the artificial and interface thermal resistance effects. Constant wall temperature simulations are shown to exhibit superposition of the artificial and interface thermal resistance values at the liquid/solid interface, while applying thermostat on wall layers sufficiently away from the liquid/solid interface results in consistent predictions of the interface thermal resistance. Injecting and extracting energy from each solid layer eliminates the artificial resistance. However, the method cannot directly specify a desired temperature difference between the two solid domains.

Pressure dependence of Kapitza resistance at gold/water and silicon/water interfaces

We conducted non-equilibrium molecular dynamics simulations to investigate Kapitza length at solid/liquid interfaces under the effects of bulk liquid pressures. Gold and silicon were utilized as hydrophilic and hydrophobic solid walls with different wetting surface behaviors, while the number of confined liquid water molecules was adjusted to obtain different pressures inside the channels. The interactions of solid/liquid couples were reparameterized accurately by measuring the water contact angle of solid substrates. In this paper, we present a thorough analysis of the structure, normal stress, and temperature distribution of liquid water to elucidate thermal energy transport across interfaces. Our results demonstrate excellent agreement between the pressures of liquid water in nano-channels and published thermodynamics data. The pressures measured as normal stress components were characterized using a long cut-off distance reinforced by a long-range van der Waals tail correction term. To clarify the effects of bulk liquid pressures on water structure at hydrophilic and hydrophobic solid surfaces, we defined solid/liquid interface spacing as the distance between the surface and the peak value of the first water density layer. Near the gold surface, we found that interface spacing and peak value of first water density layer were constant and did not depend on bulk liquid pressure; near the silicon surface, those values depended directly upon bulk liquid. Our results reveal that the pressure dependence of Kapitza length strongly depends on the wettability of the solid surface. In the case of the hydrophilic gold surface, Kapitza length was stable despite increasing bulk liquid pressure, while it varied significantly at the hydrophobic silicon surface.

Scale effects in gas nano flows

Most previous studies on gas transport in nano-scale confinements assume dynamic similarity with rarefied gas flows, and employ kinetic theory based models. This approach is incomplete, since it neglects the van der Waals forces imposed on gas molecules by the surfaces. Using three-dimensional molecular dynamics (MD) simulations of force driven gas flows, we show the significance of wall force field in nano-scale confinements by defining a new dimensionless parameter (B) as the ratio of the wall force-penetration length to the channel height. Investigation of gas transport in different nano-channels at various Knudsen numbers show the importance of wall force field for finite B values, where the dynamic similarity between the rarefied and nano-scale gas flows break down. Comparison of MD results employing molecularly structured three-dimensional walls versus reflection of gas molecules from a two-dimensional planar surface with Maxwell distribution show that the nano-confinement effects cannot be resolved...

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.

Electric field controlled transport of water in graphene nano-channels

Motivated by electrowetting-based flow control in nano-systems, water transport in graphene nano-channels is investigated as a function of the applied electric field. Molecular dynamics simulations are performed for deionized water confined in graphene nano-channels subjected to opposing surface charges, creating an electric field across the channel. Water molecules respond to the electric field by reorientation of their dipoles. Oxygen and hydrogen atoms in water face the anode and cathode, respectively, and hydrogen atoms get closer to the cathode compared to the oxygen atoms near the anode. These effects create asymmetric density distributions that increase with the applied electric field. Force-driven water flows under electric fields exhibit asymmetric velocity profiles and unequal slip lengths. Apparent viscosity of water increases and the slip length decreases with increased electric field, reducing the flow rate. Increasing the electric field above a threshold value freezes water at room temperature.

Atomic density effects on temperature characteristics and thermal transport at grain boundaries through a proper bin size selection

This study focuses on the proper characterization of temperature profiles across grain boundaries (GBs) in order to calculate the correct interfacial thermal resistance (ITR) and reveal the influence of GB geometries onto thermal transport. The solid-solid interfaces resulting from the orientation difference between the (001), (011), and (111) copper surfaces were investigated. Temperature discontinuities were observed at the boundary of grains due to the phonon mismatch, phonon backscattering, and atomic forces between dissimilar structures at the GBs. We observed that the temperature decreases gradually in the GB area rather than a sharp drop at the interface. As a result, three distinct temperature gradients developed at the GB which were different than the one observed in the bulk solid. This behavior extends a couple molecular diameters into both sides of the interface where we defined a thickness at GB based on the measured temperature profiles for characterization. Results showed dependence on the selection of the bin size used to average the temperature data from the molecular dynamics system. The bin size on the order of the crystal layer spacing was found to present an accurate temperature profile through the GB. We further calculated the GB thickness of various cases by using potential energy (PE) distributions which showed agreement with direct measurements from the temperature profile and validated the proper binning. The variation of grain crystal orientation developed different molecular densities which were characterized by the average atomic surface density (ASD) definition. Our results revealed that the ASD is the primary factor affecting the structural disorders and heat transfer at the solid-solid interfaces. Using a system in which the planes are highly close-packed can enhance the probability of interactions and the degree of overlap between vibrational density of states (VDOS) of atoms forming at interfaces, leading to a reduced ITR. Thus, an accurate understanding of thermal characteristics at the GB can be formulated by selecting a proper bin size.

Wetting of chemically heterogeneous striped surfaces: Molecular dynamics simulations

Using molecular dynamics simulations, we thoroughly investigated the wetting behaviors of a chemically heterogeneous striped substrate patterned with two different wetting materials, face-centered cubic gold and face-centered cubic silver. We analyzed the density distributions, normal stress distributions, surface tensions, and contact angles of a water droplet placed on the substrates at different heterogeneities. We found that the density and stress profile of the water droplet near the substrate-water interface were noticeably affected by altering the gold and silver contents in the substrate. Specifically, a greater portion of gold (more wetting) or smaller portion of silver (less wetting) in the substrate composition induced higher densities and higher normal stresses in the vicinity of the substrate surface. Also, it was observed that the surface tensions at liquid-vapor interface and solid-vapor interface were not largely impacted by the change of the substrate composition while the solid-liquid surface tension decreased exponentially with increasing fraction of gold. Most importantly, we found that contact angle of a nanometer-sized water droplet resting on the chemically heterogeneous striped substrate does not show linear dependence on corresponding surface fractions like that predicted by Cassie-Baxter model at the macro-scale. Consequently, we proposed a method for successfully predicting the contact angle by including the critical effects of the substrate heterogeneity on both surface tensions and line tension at the three-phase contact line of the water droplet and the chemically striped substrate.

Heat Conduction and Interface Thermal Resistance in Liquid Argon Filled Silver and Graphite Nanochannels

Heat conduction between two parallel solid walls separated by liquid argon is investigated using three-dimensional molecular dynamics (MD) simulations. Liquid argon molecules confined in silver and graphite nano-channels are examined separately. Heat flux and temperature distribution within the nano-channels are calculated by maintaining a fixed temperature difference between the two solid surfaces. Temperature profiles are linear sufficiently away from the walls, and heat transfer in liquid argon obeys the Fourier law. Temperature jump due to the interface thermal resistance (i.e., Kapitza length) is characterized as a function of the wall temperature. MD results enabled development of a phenomenological model for the Kapitza length, which is utilized as the coefficient of a Navier-type temperature jump boundary condition using continuum heat conduction equation. Analytical solution of this model results in successful predictions of temperature distribution in liquid-argon confined in silver and graphite nano-channels as thin as 7 nm and 3.57 nm, respectively.Copyright