Molecular Dynamics Studies on the Condensation Coefficient of Water

Effects of Organic Surface Contamination on the Mass Accommodation Coefficient of Water: A Molecular Dynamics Study

The mass accommodation coefficient (MAC), a parameter that quantifies the possibility of a phase change to occur at a liquid-vapor interface, can strongly affect the evaporation and condensation rates at a liquid surface. Due to the various challenges in experimental determination of the MAC, molecular dynamics (MD) simulations have been widely used to study the MAC on liquid surfaces with no impurities or contaminations. However, experimental studies show that airborne hydrocarbons from various sources can adsorb on liquid surfaces and alter the liquid surface properties. In this work, therefore, we study the effects of organic surface contamination, which is immiscible with water, on the MAC of water by equilibrium and nonequilibrium MD simulations. The equilibrium MD simulation results show that the MAC decreases almost linearly with increasing surface coverage of the organic contaminants. With the MAC determined from EMD simulations, the nonequilibrium MD simulation results show that the Schrage equation, which has been proven to be accurate in predicting the evaporation/condensation rates on clean liquid surfaces, is also accurate in predicting the condensation rate at contaminated water surfaces. The key assumption about the molecular velocity distribution in the Schrage analysis is still valid for condensing vapor molecules near contaminated water surfaces. We also find that under nonequilibrium conditions the adsorption of the water vapor molecules on the organic surface results in an adsorption vapor flux near the contaminated water surface. When the water surface is almost fully covered by the model organic contaminants, the adsorption flux dominates over the water condensation flux and leads to a false prediction of the MAC from the Schrage equation.

Mechanism, Kinetics, and Potential of Mean Force of Evaporation of Water from Aqueous Sodium Chloride Solutions of Varying Concentrations

The mechanism, kinetics, and potential of mean force of evaporation of water from aqueous NaCl solutions are investigated through both unbiased molecular dynamics simulations and also biased simulations using the umbrella sampling method. The results are obtained for aqueous solutions of three different NaCl concentrations ranging from 0.6 to 6.0 m and also for pure water. The rate of evaporation is found to decrease in the presence of ions. It is found that the process of evaporation of a surface water molecule from ionic solutions can be triggered through its collision with another water or chloride ion. Such collisions provide the additional kinetic energy that is required for evaporation. However, when the collision takes place with a Cl^- ion, the evaporation of the escaping water also involves a collision with water in the vicinity of the ion at the same time along with the ion-water collision. These two collisions together provide the required kinetic energy for escape of the evaporating water molecule. Thus, the mechanism of evaporation process of ionic solutions can be more complex than that of pure water. The potential of mean force (PMF) of evaporation is found to be positive and it increases with increasing ion concentration. Also, no barrier in the PMF is found to be present for the condensation of water from vapor phase to the surfaces of the solutions. A detailed analysis of the unsuccessful evaporation attempts by surface water molecules is also made in the current study.

Evaporation in nano/molecular materials

Evaporation is a physical phenomenon with fundamental significance to both nature and technology ranging from plant transpiration to DNA engineering. Various analytical and empirical relationships have been proposed to characterize evaporation kinetics at macroscopic scales. However, theoretical models to describe the kinetics of evaporation from nano and sub-nanometer (molecular) confinements are absent. On the other hand, the fast advancements in technology concentrated on development of nano/molecular-scale devices demand appropriate models that can accurately predict physics of phase-change in these systems. A thorough understanding of the physics of evaporation in nano/molecular materials is, thus, of critical importance to develop the required models. This understanding is also crucial to explain the intriguing evaporation-related phenomena that only take place when the characteristic length of the system drops to several nanometers. Here, we comprehensively review the underlying physics of evaporation phenomenon and discuss the effects of nano/molecular confinement on evaporation. The role of liquid-wall interface-related phenomena including the effects of disjoining pressure and flow slippage on evaporation from nano/molecular confinements are discussed. Different driving forces that can induce evaporation in small confinements, such as heat transfer, pressure drop, cavitation and density fluctuations are elaborated. Hydrophobic confinement induced evaporation and its potential application for synthetic ion channels are discussed in detail. Evaporation of water as molecular clusters rather than isolated molecules is discussed. Despite the lack of experimental investigations on evaporation at nanoscale, there exist an extensive body of literature that have applied different simulation techniques to predict the phase change behavior of liquids in nanoconfinements. We infer that exploring the effect of electrostatic interactions and flow slippage to enhance evaporation from nanoconduits is an interesting topic for future endeavors. Further future studies could be devoted to developing nano/molecular channels with evaporation-based gating mechanism and utilization of 2D materials to tune energy barrier for evaporation leading to enhanced evaporation.

Molecular dynamics study on characteristics of reflection and condensation molecules at vapor-liquid equilibrium state

The kinetic boundary condition (KBC) represents the evaporation or condensation of molecules at the vapor–liquid interface for molecular gas dynamics (MGD). When constructing the KBC, it is necessary to classify molecular motions into evaporation, condensation, and reflection in molecular-scale simulation methods. Recently, a method that involves setting the vapor boundary and liquid boundary has been used for classifying molecules. The position of the vapor boundary is related to the position where the KBC is applied in MGD analyses, whereas that of the liquid boundary has not been uniquely determined. Therefore, in this study, we conducted molecular dynamics simulations to discuss the position of the liquid boundary for the construction of KBCs. We obtained some variables that characterize molecular motions such as the positions that the molecules reached and the time they stayed in the vicinity of the interface. Based on the characteristics of the molecules found from these variables, we investigated the valid position of the liquid boundary. We also conducted an investigation on the relationship between the condensation coefficient and the molecular incident velocity from the vapor phase to the liquid phase. The dependence of the condensation coefficient on the incident velocity of molecules was confirmed, and the value of the condensation coefficient becomes small in the low-incident-velocity range. Furthermore, we found that the condensation coefficient in the non-equilibrium state shows almost the same value as that in the equilibrium state, although the corresponding velocity distribution functions of the incident velocity significantly differ from each other.

Thermal transport across the interface between liquid n-dodecane and its own vapor: A molecular dynamics study

There are two possible thermal transport mechanisms at liquid-gas interfaces, namely, evaporation/condensation (i.e., heat transfer by liquid-vapor phase change at liquid surfaces) and heat conduction (i.e., heat exchange by collisions between gas molecules and liquid surfaces). Using molecular dynamics (MD) simulations, we study thermal transport across the liquid-vapor interface of a model n-dodecane (C₁₂H₂₆) under various driving force conditions. In each MD simulation, we restrict the thermal energy to be transferred across the liquid-vapor interface by only one mechanism. In spite of the complex intramolecular interactions in n-dodecane molecules, our modeling results indicate that the Schrage relationships, which were shown to give accurate predictions of evaporation and condensation rates of monatomic fluids, are also valid in the prediction of evaporation and condensation rates of n-dodecane. In the case of heat conduction at the liquid-vapor interface of n-dodecane, the interfacial thermal conductance obtained from MD simulations is consistent with the prediction from the kinetic theory of gases. The fundamental understanding of thermal transport mechanisms at liquid-gas interfaces will allow us to formulate appropriate boundary conditions for continuum modeling of heating and evaporation of small fuel droplets.

A unified relationship for evaporation kinetics at low Mach numbers

We experimentally realized and elucidated kinetically limited evaporation where the molecular gas dynamics close to the liquid–vapour interface dominates the overall transport. This process fundamentally dictates the performance of various evaporative systems and has received significant theoretical interest. However, experimental studies have been limited due to the difficulty of isolating the interfacial thermal resistance. Here, we overcome this challenge using an ultrathin nanoporous membrane in a pure vapour ambient. We demonstrate a fundamental relationship between the evaporation flux and driving potential in a dimensionless form, which unifies kinetically limited evaporation under different working conditions. We model the nonequilibrium gas kinetics and show good agreement between experiments and theory. Our work provides a general figure of merit for evaporative heat transfer as well as design guidelines for achieving efficient evaporation in applications such as water purification, steam generation, and thermal management.

Evaporation plays a key role in applications such as cooling and desalination. Here, the authors experimentally demonstrated a unifying relationship between dimensionless flux and driving potential for evaporation kinetics under different working conditions.

Evaporation Mass Flux: A Predictive Model and Experiments

Evaporation is a fundamental and core phenomenon in a broad range of disciplines including power generation and refrigeration systems, desalination, electronic/photonic cooling, aviation systems, and even biosciences. Despite its importance, the current theories on evaporation suffer from fitting coefficients with reported values varying in a few orders of magnitude. Lack of a sound model impedes simulation and prediction of characteristics of many systems in these disciplines. Here, we studied evaporation at a planar liquid-vapor interface through a custom-designed, controlled, and automated experimental setup. This experimental setup provides the ability to accurately probe thermodynamic properties in vapor, liquid, and close to the liquid-vapor interface. Through analysis of these thermodynamic properties in a wide range of evaporation mass fluxes, we cast a predictive model of evaporation based on nonequilibrium thermodynamics with no fitting parameters. In this model, only the interfacial temperatures of liquid and vapor phases along with the vapor pressure are needed to predict evaporation mass flux. The model was validated by the reported study of an independent research group. The developed model provides a foundation for all liquid-vapor phase change studies including energy, water, and biological systems.

Evaporation of nanoscale water on a uniformly complete wetting surface at different temperatures

The evaporation of nanoscale water films on surfaces affects many processes in nature and industry. Using molecular dynamics (MD) simulations, we show the evaporation of a nanoscale water film on a uniformly complete wetting surface at different temperatures. With the increase in temperature, the growth of the water evaporation rate becomes slow. Analyses show that the hydrogen bond (H-bond) lifetimes and orientational autocorrelation times of the outermost water film decrease slowly with the increase in temperature. Compared to a thicker water film, the H-bond lifetimes and orientational autocorrelation times of a monolayer water film are much slower. This suggests that the lower evaporation rate of the monolayer water film on a uniformly complete wetting surface may be caused by the constriction of the water rotation due to the substrate. This finding may be helpful for controlling nanoscale water evaporation within a certain range of temperatures.

Modeling of Evaporation from Nanopores with Nonequilibrium and Nonlocal Effects

Evaporation from nanopores is of fundamental interest in nature and various industrial applications. We present a theoretical framework to elucidate evaporation and transport within nanopores by incorporating nonequilibrium effects due to the deviation from classical kinetic theory. Additionally, we include the nonlocal effects arising from phase change in nanoporous geometries and the self-regulation of the shape and position of the liquid-vapor interface in response to different operating conditions. We then study the effects of different working parameters to determine conditions suitable for maximizing evaporation from nanopores.

Nonequilibrium kinetic boundary condition at the vapor-liquid interface of argon

A boundary condition for the Boltzmann equation (kinetic boundary condition, KBC) at the vapor-liquid interface of argon is constructed with the help of molecular dynamics (MD) simulations. The KBC is examined at a constant liquid temperature of 85 K in a wide range of nonequilibrium states of vapor. The present investigation is an extension of a previous one by Ishiyama, Yano, and Fujikawa [Phys. Rev. Lett. 95, 084504 (2005)] and provides a more complete form of the KBC. The present KBC includes a thermal accommodation coefficient in addition to evaporation and condensation coefficients, and these coefficients are determined in MD simulations uniquely. The thermal accommodation coefficient shows an anisotropic behavior at the interface for molecular velocities normal versus tangential to the interface. It is also found that the evaporation and condensation coefficients are almost constant in a fairly wide range of nonequilibrium states. The thermal accommodation coefficient of the normal velocity component is almost unity, while that of the tangential component shows a decreasing function of the density of vapor incident on the interface, indicating that the tangential velocity distribution of molecules leaving the interface into the vapor phase may deviate from the tangential parts of the Maxwell velocity distribution at the liquid temperature. A mechanism for the deviation of the KBC from the isotropic Maxwell KBC at the liquid temperature is discussed in terms of anisotropic energy relaxation at the interface. The liquid-temperature dependence of the present KBC is also discussed.

Insight into the molecular mechanism of water evaporation via the finite temperature string method

The process of water's evaporation at its liquid/air interface has proven challenging to study experimentally and, because it constitutes a rare event on molecular time scales, presents a challenge for computer simulations as well. In this work, we simulated water's evaporation using the classical extended simple point charge model water model, and identified a minimum free energy path for this process in terms of 10 descriptive order parameters. The measured free energy change was 7.4 kcal/mol at 298 K, in reasonable agreement with the experimental value of 6.3 kcal/mol, and the mean first-passage time was 1375 ns for a single molecule, corresponding to an evaporation coefficient of 0.25. In the observed minimum free energy process, the water molecule diffuses to the surface, and tends to rotate so that its dipole and one O-H bond are oriented outward as it crosses the Gibbs dividing surface. As the water molecule moves further outward through the interfacial region, its local density is higher than the time-averaged density, indicating a local solvation shell that protrudes from the interface. The water molecule loses donor and acceptor hydrogen bonds, and then, with its dipole nearly normal to the interface, stops donating its remaining hydrogen bond. At that point, when the final, accepted hydrogen bond is broken, the water molecule is free. We also analyzed which order parameters are most important in the process and in reactive trajectories, and found that the relative orientation of water molecules near the evaporating molecule, and the number of accepted hydrogen bonds, were important variables in reactive trajectories and in kinetic descriptions of the process.

Raman Thermometry Measurements of Free Evaporation from Liquid Water Droplets

Recent theoretical and experimental studies of evaporation have suggested that on average, molecules in the higher-energy tail of the Boltzmann distribution are more readily transferred into the vapor during evaporation. To test these conclusions, the evaporative cooling rates of a droplet train of liquid water injected into vacuum have been studied via Raman thermometry. The resulting cooling rates are fit to an evaporative cooling model based on Knudsen's maximum rate of evaporation, in which we explicitly account for surface cooling. We have determined that the value of the evaporation coefficient (gamma(e)) of liquid water is 0.62 +/- 0.09, confirming that a rate-limiting barrier impedes the evaporation rate. Such insight will facilitate the formulation of a microscopic mechanism for the evaporation of liquid water.

Isotope Fractionation of Water during Evaporation without Condensation

The microscopic events engendering liquid water evaporation have received much attention over the last century, but remain incompletely understood. We present measurements of isotope fractionation occurring during free molecular evaporation from liquid microjets and show that the isotope ratios of evaporating molecules exhibit dramatic differences from equilibrium vapor values, strong variations with the solution deuterium mole fraction, and a clear temperature dependence. These results indicate the existence of an energetic barrier to evaporation and that the evaporation coefficient of water is less than unity. These new insights into water evaporation promise to advance our understanding of the processes that control the formation and lifetime of clouds in the atmosphere.