Search

A comprehensive set of single-component and binary isotherms were collected for ethanol/water adsorption into the siliceous forms of 185 known zeolites using grand-canonical Monte Carlo simulations. Using these data, a systematic analysis of ideal/real adsorbed-solution theory (IAST/RAST) was conducted and activity coefficients were derived for ethanol/water mixtures adsorbed in different zeolites based on RAST. It was found that activity coefficients of ethanol are close to unity while activity coefficients of water are larger in most zeolites, indicating a positive excess free energy of the mixture. This observation can be attributed to water/ethanol interactions being less favorable than water/water interactions in the single-component adsorption of water at comparable loadings. The deviation from ideal behavior can be highly structure-dependent but no clear correlation with pore diameters was identified. Our analysis also demonstrates the following: (1) accurate unary isotherms in the low-loading regime are critical for obtaining physically sensible activity coefficients; (2) the global regression scheme to solve for activity model parameters performs better than fitting activity models to activity coefficients calculated locally at each binary state point; and (3) including the dependence on adsorption potential offers only a minor benefit for describing binary adsorption at the lowest fugacities. Finally, the Margules activity model was found incapable of capturing the non-ideal adsorption behavior over the entire range of fugacities and compositions in all zeolites, but for conditions typical of solution-phase adsorption, RAST predictions using zeolite-specific or even bulk Margules parameters provide an improved description compared to IAST. [ABSTRACT FROM AUTHOR]

A linker extension strategy for generating metal-organic frameworks (MOFs) with superior moisture-capturing properties is presented. Applying this design approach involving experiment and computation results in MOF-LA2-1 {[Al(OH)(PZVDC)], where PZVDC2- is (E)-5-(2-carboxylatovinyl)-1H-pyrazole-3-carboxylate}, which exhibits an approximately 50% water capacity increase compared to the state-of-the-art water-harvesting material MOF-303. The power of this approach is the increase in pore volume while retaining the ability of the MOF to harvest water in arid environments under long-term uptake and release cycling, as well as affording a reduction in regeneration heat and temperature. Density functional theory calculations and Monte Carlo simulations give detailed insight pertaining to framework structure, water interactions within its pores, and the resulting water sorption isotherm.

The speed of sound in bubbly water is an important parameter in the wave equations governing pressure–density relations for turbulent multi-phase flow simulations. Recent molecular simulation results indicate that, for bubbles that are thermodynamically stable at finite volume conditions, the derivative of total pressure P with density ρ has a negative sign, complicating the interpretation of the speed of sound. We show that such a negative compressibility is thermodynamically consistent in a single-component two-phase model at finite volume, and identify an empirically derived equation of state to illustrate that this observation is not an artifact of small simulation length scales. To reconcile this thermodynamic relation with measurements of sound propagation, we decompose the derivative ∂P/∂ρ for bubbly water into its constituent phases to identify absorptive and transmissive contributors, both with an equation of state and using molecular simulations. We find that the speed of sound in the liquid phase remains real-valued while the bubble attenuates sound, giving a negative system compressibility. The inclusion of N2 molecules in molecular simulations illustrates that these observations are robust and hold also for mixtures. From these simulations, we also compute scattering functions for bubbly systems to identify oscillations associated with the speed of sound. Finally, the spherical harmonic modes of bubble oscillations are analyzed in the context of resonance with propagating waves. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 106 and 4.5 × 109 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L−1. Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models' differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models' thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. We also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state. [ABSTRACT FROM AUTHOR]

Mechanisms of CaCO3 nucleation from solutions that depend on multistage pathways and the existence of species far more complex than simple ions or ion pairs have recently been proposed. Herein, we provide a tightly coupled theoretical and experimental study on the pathways that precede the initial stages of CaCO3 nucleation. Starting from molecular simulations, we succeed in correctly predicting bulk thermodynamic quantities and experimental data, including equilibrium constants, titration curves, and detailed x-ray absorption spectra taken from the supersaturated CaCO3 solutions. The picture that emerges is in complete agreement with classical views of cluster populations in which ions and ion pairs dominate, with the concomitant free energy landscapes following classical nucleation theory.

Computational high-throughput screening was carried out to assess a large number of experimentally reported metal-organic frameworks (MOFs) and zeolites for their utility in hexane isomer separation. Through the work, we identified many MOFs and zeolites with high selectivity (SL+M > 10) for the group of n-hexane, 2-methylpentane, and 3-methylpentane (linear and monobranched isomers) versus 2,2-dimethylbutane and 2,3-dimethylbutane (dibranched isomers). This group of selective sorbents includes VICDOC (Fe2(BDP)3), a MOF with triangular pores that is known to exhibit high isomer selectivity and capacity. For three of these structures, the adsorption isotherms for a 10-component mixture of hexane and heptane isomers were calculated. Subsequent simulations of column breakthrough curves showed that the DEYVUA MOF exhibits a longer process cycle time than VICDOC MOF or MRE zeolite, which are previously reported, high-performing materials, illustrating the importance of capacity in designing MOFs for practical applications. Among the identified candidates, we synthesized and characterized a MOF in a new copper form with high predicted adsorbent capacity (qL+M > 1.2 mol/L) and moderately high selectivity (SL+M ≈ 10). Finally, we examined the role of pore shape in hexane isomer separations, especially of triangular-shaped pores. We show through the potential energy surface and three-dimensional siting analyses that linear alkanes do not populate the corners of narrow triangular channels and that structures with nontriangular pores can efficiently separate hexane isomers. Detailed thermodynamic analysis illustrates how differences in the free energy of adsorption contribute to shape-selective separation in nanoporous materials.