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.

Modeling the Ultrafast Electron Attachment Dynamics of Solvated Uracil.

Electron attachment to DNA by low energy electrons can lead to DNA damage, so a fundamental understanding of how electrons interact with the components of nucleic acids in solution is an open challenge. In solution, low energy electrons can generate presolvated electrons, epre-, which are efficiently scavanged by pyrimidine nucleobases to form transient negative ions, able to relax to either stable valence bound anions or undergo dissociative electron detachment or transfer to other parts of DNA/RNA leading to strand breakages. In order to understand the initial electron attachment dynamics, this paper presents a joint molecular dynamics and high-level electronic structure study into the behavior of the electronic states of the solvated uracil anion. Both the valence π* and nonvalence epre- states of the solvated uracil system are studied, and the effect of the solvent environment and the geometric structure of the uracil core are uncoupled to gain insight into the physical origin of the stabilization of the solvated uracil anion. Solvent reorganization is found to play a dominant role followed by relaxation of the uracil core.

Probing Heterogeneous Charge Distributions at the α-Al2O3(0001)/H2O Interface.

Unlike metal or semiconductor electrodes, the surface charge resulting from the protonation or deprotonation of insulating mineral oxides is highly localized and heterogeneous in nature. In this work the Stark active C≡N stretch of potassium thiocyanate is used as a molecular probe of the heterogeneity of the interfacial electrostatic potential at the α-Al2O3(0001)/H2O interface. Vibrational sum frequency generation (vSFG) measurements performed in the OH stretching region suggest that thiocyanate species organize interfacial water similarly to halide ions. Changes in the electrostatic potential are then tracked via Stark shifts of the vibrational frequency of the thiocyanate stretch. Our vSFG measurements show that we can simultaneously measure the vSFG response of SCN- ions experiencing charged and neutral surface sites. We assign local potentials of +308 and -154 mV to positively and negatively charged aluminol groups that are present at pH = 4 and pH = 10, respectively. Thiocyanate anions at positively charged surface sites and negatively charged surface sites and those participating in contact ion pairing adopt similar orientations and are oppositely oriented relative to thiocyanate ions near neutral surface sites. All four species followed Langmuir adsorption isotherms. Density functional theory-molecular dynamics (DFT-MD) simulations of SCN- near the neutral α-Al2O3(0001)/H2O interface show that the vSFG response in the C≡N stretch region originates from a SCN-H-O-Al complex, suggesting the surface site specificity of these experiments. To our knowledge this is the first spectroscopic measurement of local potentials associated with a heterogeneously charged surface. The ability to probe the evolution of local charges in situ could provide vital insight into many industrial, electrochemical, and geochemically relevant interfaces.

Aqueous solvation of the chloride ion revisited with density functional theory: impact of correlation and exchange approximations.

The specificity of aqueous halide solvation is fundamental to a wide range of bulk and interfacial phenomena spanning from biology to materials science. Halide polarizability is thought to drive the ion specificity, and if so, it is essential to have an accurate description of the electronic properties of halide ions in water. To this end, the solvation of the chloride anion, Cl- has been reinvestigated with state-of-the-art density functional theory. Specifically, the PBE-D3, PBE0-D3, and SCAN functionals have been employed to probe the impact of correlation and exchange approximations. Anticipating the findings, adding exact exchange improves the electronic structure, but simultaneously significantly reduces the Cl- polarizability, resulting in an over-structured Cl-O radial distribution function (RDF) and longer water H-bond lifetimes to Cl-. SCAN does not yield as much improvement in the energetics of Cl- relative to bulk water, but does result in a smaller reduction of the polarizability and thus a less structured Cl-O RDF, which agrees better with experiment. Special consideration is therefore warranted in assessing the impact of exchange on the energy, charge density, and the charge density response when designing and testing hybrid functionals for aqueous halide solvation.

Bond-Dependent Thole Model for Polarizability and Spectroscopy.

We present a new model for the calculation of molecular polarizabilities from effective atomic polarizabilities. This model is based on the Thole modified dipole interaction model for molecular polarizabilities, where the total polarizability is computed as a sum of effective atomic polarizabilities modified by dipole-dipole interactions. We extend this model by making the atomic polarizabilities explicit functions of the interatomic distances, scaling them by the radius of the volume an atom occupies in a molecule. We use the SCAN functional to show that this model, denoted TholeL, yields accurate molecular polarizabilities with little dependence on the training set. We also demonstrate that the TholeL model yields accurate polarizabilities for configurations far from the ground state structure for a wide range of molecules. Finally, we show that the TholeL model can be used to generate accurate Raman spectra for water, crystalline urea, and urea in water from ab initio molecular dynamics simulations.

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.

Bond Polarizability Model for Sum Frequency Generation at the Al2O3(0001)-H2O Interface.

Sum Frequency Generation (SFG) is a powerful, surface-specific vibrational probe ideally suited to studying buried interfaces; however, insight from theory is often necessary to explain the microscopic origins of the spectral features. To calculate the SFG spectrum at an insulating solid/liquid interface, we develop a flexible polarizability model that takes local dipole interactions into account, rather than assuming additive polarizabilities. We use this model to calculate bond dipoles and polarizabilities that reflect the local geometry of the interface. We apply our method to the Al2O3(0001)-H2O interface, where we reproduce the experimental spectrum and show the two H stretching peaks come from solvent and surface modes separately, not from H2O molecules with different coordination numbers as previously thought. Our work therefore emphasizes the importance of treating both surface and solvent at the same level of theory for accurate spectroscopic calculations.

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.