Controlling the Structure of MoS2 Membranes via Covalent Functionalization with Molecular Spacers.

Restacked two-dimensional (2D) materials represent a new class of membranes for water-ion separations. Understanding the interplay between the 2D membrane's structure and the constituent material's surface chemistry to its ion sieving properties is crucial for further membrane development. Here, we reveal, and tune via covalent functionalization, the structure of MoS2-based membranes. We find features on both the ∼1 nm (interlayer spacing) and ∼100 nm (mesoporous voids between layers) length scales that evolve with the hydration level. The functional groups act as permanent molecular spacers, preventing local impermeability caused by irreversible restacking and promoting the uniform rehydration of the membrane. Molecular dynamics simulations show that the choice of functional group tunes the structure of water within the MoS2 channel and consequently determines the hydrated interlayer spacing. We demonstrate that MoS2 membranes functionalized with acetic acid have consistently ∼92% rejection of Na2SO4 with a flux of ∼1.5 lm-2 hr-1 bar-1.

Unusually strong hydrogen bond cooperativity in particular (H2O)20 clusters.

Drawing upon an intuitive charge-transfer-based picture of hydrogen bonding, we demonstrate that cooperativity effects acting in concert can lead to unusually strong hydrogen bonds in neutral water clusters. The structure, vibrational, and NMR properties of a (H2O)20 pentagonal dodecahedron cluster containing such a strong hydrogen bond were studied using second-order perturbation theory and density functional theory. The hydrogen bond length was found to be shorter than 2.50 Å. A large redshift of over 2000 cm-1 with respect to the isolated water molecule was predicted for the OH stretching frequency of the donor water molecule. A large downfield shift to 13.5 ppm of the isotropic part of the 1H magnetic shielding tensor together with an unusually large shielding anisotropy of 49.9 ppm was obtained. The hydrogen bond energy was calculated using symmetry-adapted perturbation theory and was found to be more than three times stronger than a typical hydrogen bond in liquid water.

Dephasing and Decoherence in Vibrational and Electronic Line Shapes.

Absorption and emission line shapes of vibrational and electronic transitions in liquids are broadened by interactions with the "bath" (in this case, the rotational and translational degrees of freedom of all the molecules in the liquid). If these degrees of freedom are treated classically, the broadening process is often known as dephasing. If, on the other hand, the bath degrees of freedom are instead treated quantum mechanically, there is additional broadening due to what is known in the chemical-physics literature as decoherence. The question addressed in this paper is the relative importance of decoherence (bath quantum effects) and dephasing. We present general developments of this subject for absorption and emission line shapes, discover several new relationships connecting classical and quantum treatments of the bath, and also consider the Stokes shift (difference in peak frequencies in absorption and emission). We next draw some general conclusions by considering a model system whose transition-frequency time-correlation function has only one bath time scale. We then consider a realistic system of the vibrational OH stretch transition of dilute HOD in liquid D2O at room temperature. For this system, we conclude that bath quantum effects are not very important, except for the Stokes shift. More generally, we argue that this is the case for many vibrational and most electronic transitions in room-temperature liquids.

IR Spectroscopy Can Reveal the Mechanism of K+ Transport in Ion Channels.

Ion channels like KcsA enable ions to move across cell membranes at near diffusion-limited rates and with very high selectivity. Various mechanisms have been proposed to explain this phenomenon. Broadly, there is disagreement among the proposed mechanisms about whether ions occupy adjacent sites in the channel during the transport process. Here, using a mixed quantum-classical approach to calculate theoretical infrared spectra, we propose a set of infrared spectroscopy experiments that can discriminate between mechanisms with and without adjacent ions. These experiments differ from previous ones in that they independently probe specific ion binding sites within the selectivity filter. When ions occupy adjacent sites in the selectivity filter, the predicted spectra are significantly redshifted relative to when ions do not occupy adjacent sites. Comparisons between theoretical and experimental peak frequencies will therefore discriminate the mechanisms.

Machine Learning for Vibrational Spectroscopic Maps.

Maps that relate spectroscopic properties of a vibrational mode and collective solvent coordinates have proven useful in theoretical vibrational spectroscopy of condensed-phase systems. It has been realized that the predictive power of such an approach is limited and there is no clear systematic way to improve its accuracy. Here, we propose an adaptation of Δ-machine-learning methodology that goes beyond the spectroscopic maps. The machine-learning part of our approach combines Gaussian process regression used to generate the data set with an artificial neural network used to predict spectroscopic properties of interest. A specific application to the OH-stretch frequencies and transition dipoles of water is presented. Our method approximates these properties about two times more accurately than the spectroscopic-maps-only-based approach. Our results suggest that the new approach may become useful in the study of vibrational spectroscopy of condensed-phase systems.

OH-Stretch Raman Multivariate Curve Resolution Spectroscopy of HOD/H2O Mixtures.

Recently, in an attempt to quantify the role of intermolecular OH stretching vibrational couplings in liquid water, experimental Raman spectra of HOD/H2O mixtures were analyzed using the multivariate curve resolution (Raman-MCR) algorithm. This algorithm allowed for the separation of the HOD solute-correlated spectrum from the spectrum of bulk water. The former spectrum highlights features arising from HOD itself as well as from perturbations it induces on the surrounding H2O molecules. In this work, we apply a mixed quantum-classical methodology developed in our group to simulate the isotropic Raman-MCR spectra of HOD/H2O mixtures. Our results illustrate that intermolecular coupling leads to broadening and a red shift of the OH stretching band, in good agreement with the experiment. Our theoretical analysis provides a molecular-level interpretation of Raman-MCR experiments on HOD/H2O mixtures, suggesting that perturbations affecting the OH stretching vibrational mode of HOD result from intermolecular vibrational coupling to surrounding H2O molecules extending well beyond the first solvation shell.

Mid-IR spectroscopy of supercritical water: From dilute gas to dense fluid.

Mixed quantum-classical methods are commonly used to calculate infrared spectra for condensed-phase systems. These methods have been applied to study water in a range of conditions from liquid to solid to supercooled. Here, we show that these methods also predict infrared line shapes in excellent agreement with experiments in supercritical water. Specifically, we study the OD stretching mode of dilute HOD in H2O. We find no qualitative change in the spectrum upon passing through the near-critical region (Widom line) or the hydrogen-bond percolation line. At very low densities, the spectrum does change qualitatively, becoming rovibrational in character. We describe this rovibrational spectrum from the perspective of classical mechanics and provide a classical interpretation of the rovibrational line shape for both HOD and H2O. This treatment is perhaps more accessible than the conventional quantum-mechanical treatment.

Perspective: Crossing the Widom line in no man's land: Experiments, simulations, and the location of the liquid-liquid critical point in supercooled water.

The origin of liquid water's anomalous behavior continues to be a subject of interest and debate. One possible explanation is the liquid-liquid critical point hypothesis, which proposes that supercooled water separates into two distinct liquids at low temperatures and high pressures. According to this hypothesis, liquid water's anomalies can be traced back to the critical point associated with this phase separation. If such a critical point actually exists, it is located in a region of the phase diagram known as No Man's Land (NML), where it is difficult to characterize the liquid using conventional experimental techniques due to rapid crystallization. Recently, however, experimentalists have managed to explore NML near the proposed location of the Widom line (i.e., the Kanno-Angell line), thereby providing valuable information concerning the liquid-liquid critical point hypothesis. In this perspective, we analyze these experimental results, in conjunction with molecular dynamics simulations based on the E3B3 water model and discuss their implications for the validity of the liquid-liquid critical point hypothesis and the possible location of water's second critical point.

Communication: Diffusion constant in supercooled water as the Widom line is crossed in no man's land.

According to the liquid-liquid critical point (LLCP) hypothesis, there are two distinct phases of supercooled liquid water, namely, high-density liquid and low-density liquid, separated by a coexistence line that terminates in an LLCP. If the LLCP is real, it is located within No Man's Land (NML), the region of the metastable phase diagram that is difficult to access using conventional experimental techniques due to rapid homogeneous nucleation to the crystal. However, a recent ingenious experiment has enabled measurement of the diffusion constant deep inside NML. In the current communication, these recent measurements are compared, with good agreement, to the diffusion constant of E3B3 water, a classical water model that explicitly includes three-body interactions. The behavior of the diffusion constant as the system crosses the Widom line (the extension of the liquid-liquid coexistence line into the one-phase region) is analyzed to derive information about the presence and location of the LLCP. Calculations over a wide range of temperatures and pressures show that the new experimental measurements are consistent with an LLCP having a critical pressure of over 0.6 kbar.

Percolation in supercritical water: Do the Widom and percolation lines coincide?

Hydrogen bonding distinguishes water from simpler fluids. Here, we use classical molecular dynamics to study the percolation transition in the hydrogen bond network of supercritical water. We find that, contrary to some previous work, the percolation line in both the pressure-temperature and density-temperature planes does not coincide with the Widom line. This difference stems from a fundamental distinction between the Widom line, which is thermodynamic in nature, and the percolation transition, which depends only on connectivity. For example, we show that percolation-related quantities collapse onto master curves when plotted with respect to a measure of connectivity rather than thermodynamic variables. We then use the Galam-Mauger formula to understand the properties of the hydrogen bonding network. We find that the percolation transition in supercritical water can shed light on the hydrogen bonding network in room temperature liquid water.

Fermi resonance in OH-stretch vibrational spectroscopy of liquid water and the water hexamer.

Vibrational spectroscopy of water contains a wealth of information about the structure and dynamics of this fascinating substance. Theoretical modeling of fundamental vibrational transitions in condensed water has proven difficult, and in many circumstances, one cannot reach even qualitative agreement with experiment. Due to the ability of water to form hydrogen bonds of various strengths, the OH stretching band spans several hundreds of wave numbers in the spectra, overlapping with the first overtone of the HOH bending band and triggering a resonance between these two vibrations. This effect, known as Fermi resonance, has been traditionally ignored in theoretical condensed-phase simulations due to the additional computational burden and its deemed low importance. Depending on a particular molecular environment, the Fermi resonance manifests itself from small spectral features in the spectra of liquid water to pronounced distinct peaks in the spectra of ice and water clusters. The goal of this work is to illustrate the effects of including the Fermi resonance coupling between the bending overtone and stretching fundamental vibrations in the mixed quantum-classical formalism developed by Skinner and co-workers on the IR and Raman spectra of liquid water and the water hexamer. We show that by adding the Fermi resonance coupling, we are able to reproduce the location of the peak and a shoulder on the red side of the IR spectrum as well as the bimodal structure of the polarized Raman spectrum of liquid water at 300 K. Very good agreement between theory and experiment is achieved for the IR spectra of the water hexamer as well. We suggest that the Fermi resonance should not be ignored if intricate features of spectra are of interest. In spite of these promising results obtained in the region of a spectrum where Fermi resonance is important, further development of spectroscopic maps is needed to improve agreement with the experiment outside of the frequency range affected by the Fermi resonance.

Instantaneous ion configurations in the K+ ion channel selectivity filter revealed by 2D IR spectroscopy.

Potassium channels are responsible for the selective permeation of K+ ions across cell membranes. K+ ions permeate in single file through the selectivity filter, a narrow pore lined by backbone carbonyls that compose four K+ binding sites. Here, we report on the two-dimensional infrared (2D IR) spectra of a semisynthetic KcsA channel with site-specific heavy (13C18O) isotope labels in the selectivity filter. The ultrafast time resolution of 2D IR spectroscopy provides an instantaneous snapshot of the multi-ion configurations and structural distributions that occur spontaneously in the filter. Two elongated features are resolved, revealing the statistical weighting of two structural conformations. The spectra are reproduced by molecular dynamics simulations of structures with water separating two K+ ions in the binding sites, ruling out configurations with ions occupying adjacent sites.

IR spectra of water droplets in no man's land and the location of the liquid-liquid critical point.

No man's land is the region in the metastable phase diagram of water where it is very difficult to do experiments on liquid water because of homogeneous nucleation to the crystal. There are a number of estimates of the location in no man's land of the liquid-liquid critical point, if it exists. We suggest that published IR absorption experiments on water droplets in no man's land can provide information about the correct location. To this end, we calculate theoretical IR spectra for liquid water over a wide range of temperatures and pressures, using our E3B3 model, and use the results to argue that the temperature dependence of the experimental spectra is inconsistent with several of the estimated critical point locations, but consistent with others.

Communication: Vibrational sum-frequency spectrum of the air-water interface, revisited.

Before 2015, heterodyne-detected sum-frequency-generation experiments on the air-water interface showed the presence of a positive feature at low frequency in the imaginary part of the susceptibility. However, three very recent experiments indicate that this positive feature is in fact absent. Armed with a better understanding, developed by others, of how to calculate sum-frequency spectra, we recalculate the spectrum and find good agreement with these new experiments. In addition, we provide a revised interpretation of the spectrum.

Sub- and super-Maxwellian evaporation of simple gases from liquid water.

Non-Maxwellian evaporation of light atoms and molecules (particles) such as He and H2 from liquids has been observed experimentally. In this work, we use simulations to study systematically the evaporation of Lennard-Jones particles from liquid water. We find instances of sub- and super-Maxwellian evaporation, depending on the mass of the particle and the particle-water interaction strength. The observed trends are in qualitative agreement with experiment. We interpret these trends in terms of the potential of mean force and the effectiveness and frequency of collisions during the evaporation process. The angular distribution of evaporating particles is also analyzed, and it is shown that trends in the energy from velocity components tangential and normal to the liquid surface must be understood separately in order to interpret properly the angular distributions.

Evidence for a liquid-liquid critical point in supercooled water within the E3B3 model and a possible interpretation of the kink in the homogeneous nucleation line.

Supercooled water exhibits many thermodynamic anomalies, and several scenarios have been proposed to interpret them, among which the liquid-liquid critical point (LLCP) hypothesis is the most commonly discussed. We investigated Widom lines and the LLCP of deeply supercooled water, by using molecular dynamics simulation with a newly reparameterized water model that explicitly includes three-body interactions. Seven isobars are studied from ambient pressure to 2.5 kbar, and Widom lines are identified by calculating maxima in the coefficient of thermal expansion and the isothermal compressibility (both with respect to temperature). From these data we estimate that the LLCP of the new water model is at 180 K and 2.1 kbar. The oxygen radial distribution function is calculated along the 2 kbar isobar. It shows a steep change in the height of its second peak between 180 and 185 K, which indicates a transition between the high-density liquid and low-density liquid phases and which is consistent with the ascribed location of the critical point. The good agreement of the height of the second peak of the radial distribution function between simulation and experiment at 1 bar, as a function of temperature, supports the validity of the model. The location of the LLCP within the model is close to the kink in the experimental homogeneous nucleation line. We use existing experimental data to argue that the experimental LLCP is at 168 K and 1.95 kbar and speculate how this LLCP and its Widom line might be responsible for the kink in the homogeneous nucleation line.

Low-frequency dynamics of aqueous alkali chloride solutions as probed by terahertz spectroscopy.

Terahertz (far infrared) spectroscopy provides a useful tool for probing both ionic motions in solution and the effect of ionic solutes on the dynamics of the solvent. In this study, we calculate terahertz spectra of aqueous alkali chloride solutions using classical but novel (the water model includes three-body interactions, the ion parameterization is non-standard, and the dipole surface is polarizable) molecular dynamics simulations. The calculated spectra compare reasonably well to experimental spectra. Decomposition of the calculated spectra is used to gain a deeper understanding of the physical phenomena underlying the spectra and the connection to, for instance, the vibrational density of states for the ions. The decomposed results are also used to explain many of the cation-dependent trends observed in the experimental spectra.

Super-Maxwellian helium evaporation from pure and salty water.

Helium atoms evaporate from pure water and salty solutions in super-Maxwellian speed distributions, as observed experimentally and modeled theoretically. The experiments are performed by monitoring the velocities of dissolved He atoms that evaporate from microjets of pure water at 252 K and 4-8.5 molal LiCl and LiBr at 232-252 K. The average He atom energies exceed the flux-weighted Maxwell-Boltzmann average of 2RT by 30% for pure water and 70% for 8.5m LiBr. Classical molecular dynamics simulations closely reproduce the observed speed distributions and provide microscopic insight into the forces that eject the He atoms from solution. Comparisons of the density profile and He kinetic energies across the water-vacuum interface indicate that the He atoms are accelerated by He-water collisions within the top 1-2 layers of the liquid. We also find that the average He atom kinetic energy scales with the free energy of solvation of this sparingly soluble gas. This free-energy difference reflects the steeply decreasing potential of mean force on the He atoms in the interfacial region, whose gradient is the repulsive force that tends to expel the atoms. The accompanying sharp decrease in water density suppresses the He-water collisions that would otherwise maintain a Maxwell-Boltzmann distribution, allowing the He atom to escape at high energies. Helium is especially affected by this reduction in collisions because its weak interactions make energy transfer inefficient.

Two-dimensional infrared spectroscopy of neat ice Ih.

The assignment of the distinct peaks observed in the OH stretch lineshape of ice Ih is controversial. Recent two-dimensional infrared spectroscopic measurements provided new data. The spectra are, however, challenging to interpret and here we provide simulations that help overcome experimental issues as thermal signals and finite pulse duration. We find good agreement with experiment and the difference between H2O and D2O ices is well accounted for. The overall dynamics is demonstrated to be faster than observed for the corresponding liquid water. We find that excitonic cross peaks exist between the dominant exciton peaks. This leads us to conclude that the cross peaks arise due to the formation of delocalized exciton states, which have essentially no directional correlation between their transition dipoles as opposed to what is commonly seen, for example, in isolated water, where the transition dipoles of the eigenstates are perpendicular to each other.