Detecting Liquid-Liquid Phase Separations Using Molecular Dynamics Simulations and Spectral Clustering.

A stringent test of the accuracy of empirical force fields is reproducing the phase diagram of bulk phases and mixtures. Exploring the phase diagram of mixtures requires the detection of phase boundaries and critical points. In contrast to most solid-liquid transitions, in which a global order parameter (average density) can be used to discriminate between two phases, some demixing transitions entail relatively subtle changes in the local environment of each molecule. In such cases, finite sampling errors and finite-size effects make the identification of trends in local order parameters extremely challenging. Here we analyze one such example, namely a methanol/hexane mixture, and compute several local and global structural properties. We simulate the system at various temperatures and study the structural changes associated with demixing. We show that despite a seemingly continuous transformation between mixed and demixed states, the topological properties of the H-bond network change abruptly as the system crosses the demixing line. In particular, by using spectral clustering, we show that the distribution of cluster sizes develops a fat tail (as expected from percolation theory) in the vicinity of the critical point. We illustrate a simple criterion to identify this behavior, which results from the emergence of large system-spanning clusters from a collection of aggregates. We further tested the spectral clustering analysis on a Lennard-Jones system as a standard example of a system with no H-bonds, and also, in this case, we were able to detect the demixing transition.

Probing the Gold/Water Interface with Surface-Specific Spectroscopy.

Water is an integral component in electrochemistry, in the generation of the electric double layer, and in the propagation of the interfacial electric fields into the solution; however, probing the molecular-level structure of interfacial water near functioning electrode surfaces remains challenging. Due to the surface-specificity, sum-frequency-generation (SFG) spectroscopy offers an opportunity to investigate the structure of water near working electrochemical interfaces but probing the hydrogen-bonded structure of water at this buried electrode-electrolyte interface was thought to be impossible. Propagating the laser beams through the solvent leads to a large attenuation of the infrared light due to the absorption of water, and interrogating the interface by sending the laser beams through the electrode normally obscures the SFG spectra due to the large nonlinear response of conduction band electrons. Here, we show that the latter limitation is removed when the gold layer is thin. To demonstrate this, we prepared Au gradient films on CaF2 with a thickness between 0 and 8 nm. SFG spectra of the Au gradient films in contact with H2O and D2O demonstrate that resonant water SFG spectra can be obtained using Au films with a thickness of ∼2 nm or less. The measured spectra are distinctively different from the frequency-dependent Fresnel factors of the interface, suggesting that the features we observe in the OH stretching region indeed do not arise from the nonresonant response of the Au films. With the newfound ability to probe interfacial solvent structure at electrode/aqueous interfaces, we hope to provide insights into more efficient electrolyte composition and electrode design.

Effect of Functional and Electron Correlation on the Structure and Spectroscopy of the Al2O3(001)-H2O Interface.

Oxide-water interfaces are ubiquitous, with many applications in industry and the environment, yet there is a great deal of controversy over their properties and microscopic structure. This controversy stems, in part, from the unique H-bond networks formed at different surface terminations and mineral compositions. Density functional theory simulations of these interfaces require an accurate description of both the oxide mineral and water in diverse H-bond environments. Thus, herein we simulate the Al2O3(001)-H2O interface using the PBE, PBE-TS, RPBE, SCAN, and HSE06-TS functionals to determine how calculated interfacial properties depend on the choice of functional. We find that the structure of the first few layers of water at the surface is determined by electron correlation in a way that cannot be approximated using semiemipirical van der Waals corrections. Of the functionals investigated, we find that SCAN yields the most accurate interfacial structure, dynamics, and sum frequency generation spectrum. Furthermore, SCAN leads to a reduction in the order of the 2D H-bond network of water at the alumina surface predicted by GGA functionals, leading to a significant decrease in the anisotropy of the diffusion coefficient at the surface. We emphasize the importance of using a functional which accurately describes electron correlation for more complex oxides, such as transition-metal oxides, where electron correlation may play an even greater role in determining the structure and dynamics of the oxide-water interface.

Effect of Ions on H-Bond Structure and Dynamics at the Quartz(101)-Water Interface.

We use ab initio molecular dynamics simulations to study the effect of ions on the structure and dynamics of the quartz(101)-water interface. We study several IA (Na+, Rb+) and IIA (Mg2+, Sr2+) cations, with Cl- as counterion, adsorbed onto acidic, neutral, and basic surface configurations at 300 and 373 K. We find that both cations and anions can bond directly to the surface and perturb the local H-bond network. The adsorbed ions promote the formation of intrasurface H-bonds, as shown by vibrational density of states and orientations of the surface silanols. Both local and global structural correlations of the interfacial H-bond network are studied using a structural definition of the H-bond and a network correlation function. We find the ions' effect on the solvent structure exhibits a complex dependence on specific surface interactions. The structure-making properties of ions are enhanced at the quartz surface, particularly for ions adsorbed without a complete hydration shell, and the structuring effect extends beyond the first solvation shell. The ions have a lesser effect on solvent structure in solution, especially in the presence of counterions. In fact, cations that are the greatest "structure makers" at the surface are the greatest "structure breakers" when in solution with a counterion. Therefore, we find the ions cannot be simply classified as "structure making" or "structure breaking". We discuss the implications of these findings for the effect of ions on the dissolution rate, surface charge, and solvent structure.