skip to main content


Title: Molecular Dynamics Investigation of Liquid and Vapor Interactions Near an Evaporating Interface: A Theoretical Genetics Perspective
Abstract

Understanding phase transition between the liquid and gaseous states has gained significant interest, and has been ubiquitously observed in many places ranging from natural systems to water–energy nexus and thermal management applications. Phase transition phenomena at liquid–vapor interfaces are greatly governed by intermolecular‐level kinetics, which requires the use of empirical parameters in continuum‐level relations to explain the discrete nature of molecular particles. Despite its significance, it has been a great challenge to find detailed expressions of empirical parameters such as accommodation coefficients, which represent the probabilities for phase transition of liquid or vapor molecules at the interface. Here, direct statistical measurements of accommodation coefficients are reported by tracking the trajectories of liquid and vapor molecules in molecular simulations. The measurements reveal that evaporation and condensation coefficients are different by ≈50%, whereas they have been assumed to be equal in most previous studies. Then, the indirect measurement method is studied from a perspective of theoretical genetics based on the diffusion approximation. A good agreement between two approaches suggests that diffusion approximation can contribute to provide empirical parameters with a cost‐effective method.

 
more » « less
Award ID(s):
1752147
NSF-PAR ID:
10456591
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Advanced Theory and Simulations
Volume:
3
Issue:
7
ISSN:
2513-0390
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Rotational-translational decoupling, in which translational motion is apparently enhanced over rotational motion in violation of Stokes-Einstein (SE) and Debye-Stokes-Einstein (DSE) predictions, has been observed in materials near their glass transition temperatures (Tg). This has been posited to result from ensemble averaging in the context of dynamic heterogeneity. In this work, ensemble and single molecule experiments are performed in parallel on a fluorescent probe in high molecular weight polystyrene near its Tg. Ensemble results show decoupling onset at approximately 1.15Tg, increasing to over three orders of magnitude at Tg. Single molecule measurements also show a high degree of decoupling, with typical molecules at Tgshowing translational diffusion coefficients nearly 400 times higher than expected from SE/DSE predictions. At the single molecule level, higher degree of breakdown is associated with particularly mobile molecules and anisotropic trajectories, providing support for anomalous diffusion as a critical driver of rotational-translational decoupling and SE/DSE breakdown.

     
    more » « less
  2. null (Ed.)
    Abstract. Mass accommodation is an essential process for gas–particle partitioning oforganic compounds in secondary organic aerosols (SOA). The massaccommodation coefficient is commonly described as the probability of a gasmolecule colliding with the surface to enter the particle phase. It is oftenapplied, however, without specifying if and how deep a molecule has topenetrate beneath the surface to be regarded as being incorporated into thecondensed phase (adsorption vs. absorption). While this aspect is usuallynot critical for liquid particles with rapid surface–bulk exchange, it canbe important for viscous semi-solid or glassy solid particles to distinguishand resolve the kinetics of accommodation at the surface, transfer acrossthe gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective massaccommodation coefficient αeff that depends on penetrationdepth and is a function of surface accommodation coefficient, volatility,bulk diffusivity, and particle-phase reaction rate coefficient. Applicationof αeff in the traditional Fuchs–Sutugin approximation ofmass-transport kinetics at the gas–particle interface yields SOApartitioning results that are consistent with a detailed kinetic multilayermodel (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Modelfor Simulating Aerosol Interactions and Chemistry, MOSAIC;Zaveri et al., 2014) but deviate substantially from earlier modelingapproaches not considering the influence of penetration depth and relatedparameters. For highly viscous or semi-solid particles, we show that the effective massaccommodation coefficient remains similar to the surface accommodationcoefficient in the case of low-volatility compounds, whereas it can decrease byseveral orders of magnitude in the case of semi-volatile compounds. Such effectscan explain apparent inconsistencies between earlier studies deriving massaccommodation coefficients from experimental data or from molecular dynamicssimulations. Our findings challenge the approach of traditional SOA models using theFuchs–Sutugin approximation of mass transfer kinetics with a fixed massaccommodation coefficient, regardless of particle phase state and penetrationdepth. The effective mass accommodation coefficient introduced in this studyprovides an efficient new way of accounting for the influence of volatility,diffusivity, and particle-phase reactions on SOA partitioning in processmodels as well as in regional and global air quality models. While kineticlimitations may not be critical for partitioning into liquid SOA particlesin the planetary boundary layer (PBL), the effects are likely important foramorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT)as well as in the PBL at low relative humidity and low temperature. 
    more » « less
  3. Abstract The reactive uptake of N 2 O 5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N 2 O 5 in water. The free energy profile highlights that N 2 O 5 is selectively adsorbed to the liquid–vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction–diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction–diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N 2 O 5 in more complex solutions. 
    more » « less
  4. The phase transition from subcritical to supercritical conditions, referred to as transcritical behavior, significantly impacts the evaporation and fuel–air mixing in high-pressure liquid-fuel propulsion systems. Transcritical behavior is characterized as a transition from classical two-phase evaporation to a single-phase gas-like diffusion regime as surface tension and latent heat of vaporization reduce. However, the interfacial behavior represented by the surface tension coefficient and evaporation rate during this transition which are crucial inputs for Computational Fluid Dynamics (CFD) simulations of practical transcritical fuel spray is still missing. This study aims at developing new evaporation rate and surface tension models for transcritical n-dodecane droplets using molecular dynamics (MD) simulations irrespective of the droplet size. As MD simulations are primarily limited to the nanoscale, the new models can bridge the gap between MD and continuum simulations and enable the direct application of these findings to microscopic droplets. A new characteristic timescale, i.e., “undroplet time,” is defined which marks the transition from classical two-phase evaporation to single-phase gas-like diffusion behavior. The undroplet time indicates the onset of droplet core disintegration and penetration of nitrogen molecules into the droplet, which occurs after the vanishment of the surface tension. By normalizing the time with respect to the undroplet time, the rate of surface tension decay, evaporation rate, and the rate of droplet mass depletion become independent of the droplet size. Calculation of pairwise correlation coefficients for the entire MD results shows that both surface tension coefficient and evaporation rate are strongly correlated with the background temperature, while pressure and droplet size play a less significant role past the critical point. Therefore, new models for surface tension coefficient and evaporation rate spanning from sub- to supercritical conditions are developed as a function of background pressure and temperature, which can be used in continuum simulations. The identified phase change behavior based on the undroplet time shows a good agreement with the phase change regime maps obtained using microscale experiments and nanoscale MD predictions. 
    more » « less
  5. Abstract

    Atmospheric water harvesting (AWH)—producing fresh water via collecting moisture from air—enables sustainable water delivery without geographical and hydrologic limitations. However, the fundamental design principle to prepare materials that can convert the water vapor in the air to collectible liquid water is still mostly unknown. Here, a super moisture‐absorbent gel, which is composed of hygroscopic polypyrrole chloride penetrating in hydrophilicity‐switchable polymeric network of polyN‐isopropylacrylamide, is shown. Based on such design, a high‐efficiency water production by AWH has been achieved in a broad range of relative humidity. The synergistic effect enabled by the molecular level integration of hygroscopic and hydrophilicity‐switchable polymers in a network architecture presents controllable interaction between the gel and water molecules, simultaneously realizing efficient vapor capturing, in situ water liquefaction, high‐density water storage and fast water releasing under different weather conditions. Being an effective method to regulate migration of water molecules, such design represents a novel strategy to improve the AWH, and it is also fundamental to other water management systems for environmental cooling, surficial moisturizing and beyond.

     
    more » « less