skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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. null (Ed.)
    The method of surface grating decay has been used to measure surface diffusion in the glasses of two rod-like molecules posaconazole (POS) and itraconazole (ITZ). Although structurally similar antifungal medicines, ITZ forms liquid-crystalline phases while POS does not. Surface diffusion in these systems is significantly slower than in the glasses of quasi-spherical molecules of similar volume when compared at the glass transition temperature T g . Between the two systems, ITZ has slower surface diffusion. These results are explained on the basis of the near-vertical orientation of the rod-like molecules at the surface and their deep penetration into the bulk where mobility is low. For molecular glasses without extensive hydrogen bonds, we find that the surface diffusion coefficient at T g decreases smoothly with the penetration depth of surface molecules and the trend has the double-exponential form for the surface mobility gradient observed in simulations. This supports the view that these molecular glasses have a similar mobility vs. depth profile and their different surface diffusion rates arise simply from the different depths at which molecules are anchored. Our results also provide support for a previously observed correlation between the rate of surface diffusion and the fragility of the bulk liquid. 
    more » « less
  2. 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
  3. The surface diffusion kinetics on a Pd77.5Cu6Si16.5 metallic glass is measured using a scratch smoothing method in the range of 107–57 K below the glass transition temperature. Within this temperature range, the surface diffusion coefficients are determined to vary between (8.66 ± 0.80) × 10−19 and (5.90 ± 0.60) × 10−18 m2 s−1. The corresponding activation energy is 0.93 ± 0.18 eV, which is about half the value for bulk diffusion. These measurements also corroborate the correlation between enhanced surface diffusion and liquid fragility in glasses. 
    more » « less
  4. 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
  5. Surface diffusion has been measured in the glass of an organic semiconductor, MTDATA, using the method of surface grating decay. The decay rate was measured as a function of temperature and grating wavelength, and the results indicate that the decay mechanism is viscous flow at high temperatures and surface diffusion at low temperatures. Surface diffusion in MTDATA is enhanced by 4 orders of magnitude relative to bulk diffusion when compared at the glass transition temperature T g . The result on MTDATA has been analyzed along with the results on other molecular glasses without extensive hydrogen bonds. In total, these systems cover a wide range of molecular geometries from rod-like to quasi-spherical to discotic and their surface diffusion coefficients vary by 9 orders of magnitude. We find that the variation is well explained by the existence of a steep surface mobility gradient and the anchoring of surface molecules at different depths. Quantitative analysis of these results supports a recently proposed double-exponential form for the mobility gradient: log  D( T, z) = log  D v ( T) + [log  D 0 − log  D v ( T)]exp(− z/ξ), where D( T, z) is the depth-dependent diffusion coefficient, D v ( T) is the bulk diffusion coefficient, D 0 ≈ 10 −8  m 2 /s, and ξ ≈ 1.5 nm. Assuming representative bulk diffusion coefficients for these fragile glass formers, the model reproduces the presently known surface diffusion rates within 0.6 decade. Our result provides a general way to predict the surface diffusion rates in molecular glasses. 
    more » « less