While experiments and continuum models have provided a relatively good understanding of the evaporation of macroscopic water droplets, elucidating how sessile nanodroplets evaporate is an open question critical for advancing nanotechnological applications where nanodroplets can play an essential role. Here, using molecular dynamics simulations, we find that evaporating nanodroplets, in contrast to their macroscopic counterparts, are not always in thermal equilibrium with the substrate and that the vapor concentration on the nanodroplet surface does not reach a steady state. As a result, the evaporative behavior of nanodroplets is significantly different. Regardless of hydrophobicity, nanodroplets do not follow conventional evaporation modes but instead exhibit dynamic wetting behavior characterized by huge, non-equilibrium, isovolumetric fluctuations in the contact angle and contact radius. For hydrophilic nanodroplets, the evaporation rate, controlled by the vapor concentration, decays exponentially over time. Hydrophobic nanodroplets follow stretched exponential kinetics arising from the slower thermalization with the substrate. The evaporative half-lifetime of the nanodroplets is directly related to the thermalization time scale and therefore increases monotonically with the hydrophobicity of the substrate. Finally, the evaporative flux profile along the nanodroplet surface is highly nonuniform but does not diverge at the contact line as the macroscopic continuum models predict.

Flexible nanoscale confinement is critical to understanding the role that bending fluctuations play on biological processes where soft interfaces are ubiquitous or to exploit confinement effects in engineered systems where inherently flexible 2D materials are pervasively employed. Here, using molecular dynamics simulations, we compare the phase behavior of water confined between flexible and rigid graphene sheets as a function of the in-plane density, ρ_2D. We find that both cases show commensurate mono-, bi-, and trilayered states; however, the water phase in those states and the transitions between them are qualitatively different for the rigid and flexible cases. The rigid systems exhibit discontinuous transitions between an (n)-layer and an (n+1)-layer state at particular values of ρ_2D, whereas under flexible confinement, the graphene sheets bend to accommodate an (n)-layer and an (n+1)-layer state coexisting in equilibrium at the same density. We show that the flexible walls introduce a very different sequence of ice phases and their phase coexistence with vapor and liquid phases than that observed with rigid walls. We discuss the applicability of these results to real experimental systems to shed light on the role of flexible confinement and its interplay with commensurability effects.

We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm^-3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.

The mobility of radiocesium in the environment is largely mediated by cation exchange in micaceous clays, in particular Illite-a non-swelling clay mineral that naturally contains interlayer K⁺ and has high affinity for Cs⁺. Although exchange of interlayer K⁺ for Cs⁺ is nearly thermodynamically nonselective, recent experiments show that direct, anhydrous Cs⁺-K⁺ exchange is kinetically viable and leads to the formation of phase-separated interlayers through a mechanism that remains unclear. Here, using classical atomistic simulations and density functional theory calculations, we identify a molecular-scale positive feedback mechanism in which exchange of the larger Cs⁺ for the smaller K⁺ significantly lowers the migration barrier of neighboring K⁺, allowing exchange to propagate rapidly once initiated at the clay edge. Barrier lowering upon slight increase in layer spacing (∼0.7 Å) during Cs⁺ exchange is an example of "chemical-mechanical coupling" that likely explains the observed sharp exchange fronts leading to interstratification. Interestingly, we find that these features are thermodynamically favored even in the absence of a heterogeneous layer charge distribution.

Developing accurate ab initio molecular dynamics (AIMD) models that capture both electronic reorganization and nuclear quantum effects associated with hydrogen bonding is key to quantitative understanding of bulk water and its anomalies as well as its role as a universal solvent. For condensed phase simulations, AIMD has typically relied on the generalized gradient approximation (GGA) of density functional theory (DFT) as the underlying model chemistry for the potential energy surface, with nuclear quantum effects (NQEs) sometimes modeled by performing classical molecular dynamics simulations at elevated temperatures. Here we show that the properties of liquid water obtained from the meta-GGA B97M-rV functional, when evaluated using accelerated path integral molecular dynamics simulations, display accuracy comparable to a computationally expensive dispersion-corrected hybrid functional, revPBE0-D3. We show that the meta-GGA DFT functional reproduces bulk water properties including radial distribution functions, self-diffusion coefficients, and infrared spectra with comparable accuracy of a much more expensive hybrid functional. This work demonstrates that the underlying quality of a good DFT functional requires evaluation with quantum nuclei and that high-temperature simulations are a poor proxy for properly treating NQEs.

When driven out of equilibrium, many diverse systems can form complex spatial and dynamical patterns, even in the absence of attractive interactions. Using kinetic Monte Carlo simulations, we investigate the phase behavior of a binary system of particles of dissimilar size confined between semiflexible planar surfaces, in which the nanoconfinement introduces a non-local coupling between particles, which we model as an activation energy barrier to diffusion that decreases with the local fraction of the larger particle. The system autonomously reaches a cyclical non-equilibrium state characterized by the formation and dissolution of metastable micelle-like clusters with the small particles in the core and the large ones in the surrounding corona. The power spectrum of the fluctuations in the aggregation number exhibits 1/f noise reminiscent of self-organized critical systems. We suggest that the dynamical metastability of the micellar structures arises from an inversion of the energy landscape, in which the relaxation dynamics of one of the species induces a metastable phase for the other species.

The VV10 and rVV10 nonlocal correlation functionals are consistently implemented and assessed, with the goal of determining if the rVV10 nonlocal correlation functional can replace the VV10 nonlocal correlation functional in the recently developed B97M-V density functional, to give the B97M-rV density functional. Along the way, four density functionals are simultaneously tested: VV10, rVV10, B97M-V, and B97M-rV. An initial assessment is carried out across the S22 data set, and the short-range damping variable, b, is varied for all four density functionals in order to determine the sensitivity of the functionals to the empirical parameter. The results of this test indicate that a value of b = 6 (fortuitously the same as that in B97M-V) is suitable for B97M-rV. The functionals are then compared across an extensive database of interaction energies, and it is demonstrated that B97M-rV either matches or outperforms B97M-V for all of the tests considered. Finally, the optimization of b across the S22 data set is extended to two range-separated hybrid density functionals, ωB97X-V and ωB97M-V, and a value of b = 6.2 is recommended for both ωB97X-rV and ωB97M-rV.