Dielectric response of thin water films: a thermodynamic perspective

The surface of a polar liquid presents a special environment for the solvation and organization of charged solutes, which differ from bulk behaviors in important ways. These differences have motivated many attempts to understand electrostatic response at aqueous interfaces in terms of a spatially varying dielectric permittivity, typically concluding that the dielectric constant of interfacial water is significantly lower than in the bulk liquid. Such analyses, however, are complicated by the potentially nonlocal nature of dielectric response over the short length scales of interfacial heterogeneity. Here we circumvent this problem for thin water films by adopting a thermodynamic approach. Using molecular simulations, we calculate the solvent''s contribution to the reversible work of charging a parallel plate capacitor. We find good agreement with a simple dielectric continuum model that assumes bulk dielectric permittivity all the way up to the liquid''s boundary, even for very thin (∼1 nm) films. This comparison requires careful attention to the placement of dielectric boundaries between liquid and vapor, which also resolves apparent discrepancies with dielectric imaging experiments.

Free energy calculations from molecular simulations reveal that water''s interfacial dielectric response is well-described by bulk properties.     

Hydrogen bonds in liquid water studied by photoelectron spectroscopy

The authors report on photoelectron emission spectroscopy measurements of the oxygen 1s orbital of liquid water, using a liquid microjet in ultrahigh vacuum. By suitably changing the soft x-ray photon energy, within 600-1200 eV, the electron probing depth can be considerably altered as to either predominantly access the surface or predominantly bulk water molecules. The absolute probing depth in liquid water was inferred from the evolution of the O1s signal and from comparison with aqueous salt solution. The presence of two distinctive components in the core-level photoelectron spectrum, with significantly different binding energies, is revealed. The dominant contribution, at a vertical binding energy of 538.1 eV, was found in bulk and surface sensitive spectra. A weaker component at 536.6 eV binding energy appears to be present only in bulk water. Hartree-Fock calculations of O1s binding energies in different geometric arrangements of the water network are presented to rationalize the experimental distribution of O1s electron binding energies.     

Excess-Electron Attachment and Ionization of Aqueous Uridine Monophosphate Anion

We applied quantum mechanics/classical mechanics simulations to study excess-electron attachment and ionization of uridine monophosphate anion (dUMP\begin{document}$^-$\end{document}) in explicit aqueous solutions. We calculated vertical electron affinities (VEAs), adiabatic electron affinities (AEAs), vertical detachment energies (VDEs), vertical ionization energies (VIEs), and adiabatic ionization energies (AIEs) of the 40 structures obtained from molecular dynamic trajectory. The excess-electron and hole distributions were analyzed in electron attachment and ionization of aqueous dUMP\begin{document}$^-$\end{document}. The converged mean VEA (-0.31 eV) and AEA (2.13 eV) suggest that excess-electron can easily attach to dUMP\begin{document}$^-$\end{document}. The mean vertical (-0.50 e) and adiabatic (-0.62 e) excess-electron on uracil reveal that main excess-electrons are localized on nucleobases at the most snapshots. The distributions at several special snapshots demonstrate the excess-electron delocalization over nucleobases/ribose or ribose/phosphate group after the structural relaxations of dUMP\begin{document}$^{2-}$\end{document} dianion. The VDE value (2.78 eV) indicates that dUMP\begin{document}$^{2-}$\end{document} dianion could be very stable. Moreover, the mean VIE is 8.13 eV which is in agreement with the previous calculation using solvation model. The hole distributions on uracil suggest that the nucleobases are easily ionized after the irradiation of high-energy rays. In vertical ionizations, the holes would be delocalized over uracil and ribose at several snapshots. Observing the adiabatic hole distributions, it can be found that electrons on phosphate group and holes on nucleobases can be transferred to ribose at the special snapshots in the structural relaxation of neutral species.     

Effect of solvation on the ionization of guanine nucleotide: A hybrid QM/EFP study

Ionization of nucleobases is affected by their biological environment, which includes both the effect of adjacent nucleotides as well as the presence of water around it. Guanine and its nucleotide have the lowest ionization potentials among the various DNA bases. Therefore, the threshold of ionization is dependent on that of guanine and its characterization is crucial to the prediction of interaction of light with DNA. We investigate the effect of solvation on the vertical ionization energies (VIEs) of guanine and its nucleotide. In this work, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) approach with effective fragment potential as the MM method of choice and equation‐of‐motion coupled‐cluster for ionization potential with singles and doubles (EOM‐IP‐CCSD) as the QM method. The performance of the hybrid scheme with respect to the full QM method shows an accuracy of 0.02–0.04 eV. The lowest few ionizations of the nucleotide are found to be from different parts of the moiety, that is, the nucleic acid base, phosphate, or sugar, and these ionization energies are very closely spaced giving rise to a very complicated spectrum. Furthermore, microsolvation has large effects on these ionizations and can lead to red or blue shift depending on the position of the water molecule. Even a single water molecule can change the order of ionized states in the nucleotide. The VIEs of the bulk solvated chromophores are predicted and compared to existing experimental spectra. The predominant role of polarization in the solvatochromic shift is noticed. © 2017 Wiley Periodicals, Inc.     

Assessing long-range contributions to the charge asymmetry of ion adsorption at the air–water interface

Anions generally associate more favorably with the air–water interface than cations. In addition to solute size and polarizability, the intrinsic structure of the unperturbed interface has been discussed as an important contributor to this bias. Here we assess quantitatively the role that intrinsic charge asymmetry of water''s surface plays in ion adsorption, using computer simulations to compare model solutes of various size and charge. In doing so, we also evaluate the degree to which linear response theory for solvent polarization is a reasonable approach for comparing the thermodynamics of bulk and interfacial ion solvation. Consistent with previous works on bulk ion solvation, we find that the average electrostatic potential at the center of a neutral, sub-nanometer solute at the air–water interface depends sensitively on its radius, and that this potential changes quite nonlinearly as the solute''s charge is introduced. The nonlinear response closely resembles that of the bulk. As a result, the net nonlinearity of ion adsorption is weaker than in bulk, but still substantial, comparable to the apparent magnitude of macroscopically nonlocal contributions from the undisturbed interface. For the simple-point-charge model of water we study, these results argue distinctly against rationalizing ion adsorption in terms of surface potentials inherent to molecular structure of the liquid''s boundary.

Cations and anions have different affinities for the air-water interface. The intrinsic orientation of surface molecules suggests such an asymmetry, but the bias is dominated by solvent response that is spatially local and significantly nonlinear.

Counter to expectations from conventional theories of solvation, there is a large body of both computational and experimental evidence indicating that small ions can adsorb to the air–water interface.^1–9 Implications across the biological, atmospheric and physical sciences have inspired efforts to understand the microscopic driving forces for ions associating with hydrophobic interfaces in general.^10–21 A particular emphasis has been placed on understanding ion specificity, i.e., why some ions exhibit strong interfacial affinity while others do not. Empirical trends indicate that ion size and polarizability are important factors, as could be anticipated from conventional theory. More surprisingly, the sign of a solute''s charge can effect a significant bias, with anions tending to adsorb more favorably than cations.Here we examine the microscopic origin of this charge asymmetry in interfacial ion adsorption. We specifically assess whether the thermodynamic preference can be simply and generally understood in terms of long-range biases that are intrinsic to an aqueous system surrounded by vapor. By “long-range” and “nonlocal” we refer to macroscopically large scales, i.e., collective forces that are felt at arbitrarily long distance. Such a macroscopically long-range bias is expected from the air–water interface due to its average polarization, and by some measures the bias is quite strong. By contrast, “local” contributions comprise the entire influence of a solute''s microscopic environment, including electrostatic forces from molecules that are many solvation shells away – any influence that decays over a sub-macroscopic length scale.The importance of macroscopically nonlocal contributions has been discussed extensively in the context of ion solvation in bulk liquid water, which we review in Section 1 as a backdrop for interfacial solvation. The notion that such contributions strongly influence charge asymmetry of solvation at the air–water interface has informed theoretical approaches and inspired criticism of widely used force fields for molecular simulation.^22,23 A full understanding of their role in interfacial adsorption, however, is lacking.In the course of this study, we will also evaluate the suitability of dielectric continuum theory (DCT) to describe the adsorption process. DCT has provided an essential conceptual framework for rationalizing water''s response to electrostatic perturbations. But a more precise understanding of its applicability is needed, particularly for the construction of more elaborate models (e.g., with heterogeneous polarizability near interfaces^24–26) and for the application of DCT to evermore complex (e.g., nanoconfined^27,28) environments.     

Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk

The effect of hydration on the vertical ionization energy (VIE) of thymine was characterized using equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods. We considered several microsolvated clusters as well as thymine solvated in bulk water. The VIE in bulk water was computed by averaging over solvent-solute configurations obtained from equilibrium molecular dynamics trajectories at 300 K. The effect of microsolvation was analyzed and contrasted against the combined effect of the first solvation shell in bulk water. Microsolvation reduces the ionization energy (IE) by about 0.1 eV per water molecule, while the first solvation shell increases the IE by 0.1 eV. The subsequent solvation lowers the IE, and the bulk value of the solvent-induced shift of thymine's VIE is approximately -0.9 eV. The combined effect of the first solvation shell was explained in terms of specific solute-solvent interactions, which were investigated using model structures. The convergence of IE to the bulk value requires the hydration sphere of approximately 13.5 ? radius. The performance of the EOM-IP-CCSD/EFP scheme was benchmarked against full EOM-IP-CCSD using microhydrated structures. The errors were found to be less than 0.01-0.02 eV. The relative importance of the polarization and higher multipole moments in EFP model was also investigated.     

Low Energy Electron Attachment by Some Chlorosilanes

In this paper, the rate coefficients (k) and activation energies (E_a) for SiCl₄, SiHCl₃, and Si(CH₃)₂(CH₂Cl)Cl molecules in the gas phase were measured using the pulsed Townsend technique. The experiment was performed in the temperature range of 298–378 K, and carbon dioxide was used as a buffer gas. The obtained k depended on temperature in accordance with the Arrhenius equation. From the fit to the experimental data points with function described by the Arrhenius equation, the activation energies (E_a) were determined. The obtained k values at 298 K are equal to (5.18 ± 0.22) × 10⁻¹⁰ cm³·s⁻¹, (3.98 ± 1.8) × 10⁻⁹ cm³·s⁻¹ and (8.46 ± 0.23) × 10⁻¹¹ cm³·s⁻¹ and E_a values were equal to 0.25 ± 0.01 eV, 0.20 ± 0.01 eV, and 0.27 ± 0.01 eV for SiHCl₃, SiCl₄, and Si(CH₃)₂(CH₂Cl)Cl, respectively. The linear relation between rate coefficients and activation energies for chlorosilanes was demonstrated. The DFT/B3LYP level coupled with the 6-31G(d) basis sets method was used for calculations of the geometry change associated with negative ion formation for simple chlorosilanes. The relationship between these changes and the polarizability of the attaching center (α_centre) was found. Additionally, the calculated adiabatic electron affinities (AEA) are related to the α_centre.     

Solvated electrons at the water-air interface: surface versus bulk signal in low kinetic energy photoelectron spectroscopy

Time-resolved photoelectron spectroscopy at low kinetic energies (?5 eV) is applied to dilute iodide solutions with different surface and bulk contributions. The results indicate a pronounced surface sensitivity. Signals assigned to solvated electrons near the liquid surface decay rapidly on a sub-ps timescale. In contrast to the literature, a long-lived surface solvated electron at 1.6 eV binding energy is not observed.     

Stereoselective β-mannosylations and β-rhamnosylations from glycosyl hemiacetals mediated by lithium iodide

Stereoselective β-mannosylation is one of the most challenging problems in the synthesis of oligosaccharides. Herein, a highly selective synthesis of β-mannosides and β-rhamnosides from glycosyl hemi-acetals is reported, following a one-pot chlorination, iodination, glycosylation sequence employing cheap oxalyl chloride, phosphine oxide and LiI. The present protocol works excellently with a wide range of glycosyl acceptors and with armed glycosyl donors. The method doesn''t require conformationally restricted donors or directing groups; it is proposed that the high β-selectivities observed are achieved via an S_N2-type reaction of α-glycosyl iodide promoted by lithium iodide.

Breaking from the current paradigm, a highly selective synthesis of hard-to-make β-mannosides and β-rhamnosides from simple glycosyl hemi-acetals has been achieved without using conformational restrictions.     

Time-resolved photoelectron spectroscopy of solvated electrons in aqueous NaI solution

Time-resolved photoelectron spectroscopy was used to study the energetics and dynamics of solvated electrons in aqueous solution. Solvated electrons are generated by ultrafast photodetachment in a 100 mM aqueous NaI solution. Initially, an ensemble of strongly bound ("cold") solvated electrons and an ensemble of weakly bound ("hot") electrons in an unequilibrated solvent environment are observed. We report an ultrafast recombination channel for the "hot" electrons with a rate of (800 fs)(-1) which is in competition with thermalization occurring with a rate of (1.1 ps)(-1). The thermalized electrons recombine with the iodide radical with a rate of (22 ps)(-1). About 35% of the thermalized electrons escape geminate recombination and form free, solvated electrons. The vertical detachment energy for the solvated electron is determined to be 3.40 eV. No indication for a surface-bound electron at lower binding energies was observed.     

Electron impact ionisation and UV absorption study of α- and β-pinene

We have performed a coordinated set of experiments to measure the electron impact ionisation and UV photoabsorption cross sections of α- and β-pinene. The adiabatic ionisation energies of α- and β-pinene were derived from experiment and found to be 8.3 and 8.6 eV which compared well with high-level quantum chemical calculations (G3MP2) yielding values of 8.29 and 8.41 eV. Additionally, vertical ionisation energies of 8.62 and 8.96 eV were calculated using an OVGF method. UV photoabsorption cross sections were measured using a high-resolution synchrotron radiation source and electronic states interpreted on the basis of the TD quantum chemical methods.     

Study by Optical Spectroscopy of Bismuth Emission in a Nanosecond-Pulsed Discharge Created in Liquid Nitrogen

Time-resolved optical emission spectroscopy of nanosecond-pulsed discharges ignited in liquid nitrogen between two bismuth electrodes is used to determine the main discharge parameters (electron temperature, electron density and optical thickness). Nineteen lines belonging to the Bi I system and seven to the Bi II system could be recorded by directly plunging the optical fibre into the liquid in close vicinity to the discharge. The lack of data for the Stark parameters to evaluate the broadening of the Bi I lines was solved by taking advantage of the time-resolved information supported by each line to determine them. The electron density was found to decrease exponentially from 6.5 ± 1.5 × 10¹⁶ cm⁻³ 200 ns after ignition to 1.0 ± 0.5 × 10¹⁶ cm⁻³ after 1050 ns. The electron temperature was found to be 0.35 eV, close to the value given by Saha’s equation.     

Ultraschnelle Dynamik von Strukturen und Schwingungen

Aqueous systems display ultrafast structural fluctuations in the time domain between 10^–14 and 10^–11 s, due to the limited strength of hydrogen bonds. Multidimensional ultrafast infrared spectroscopy allows for a time‐resolved observation of such elementary dynamics in bulk water, DNA oligomers and phospholipids. Both the couplings of different vibrational degrees of the hydrated system and bulk water dynamics can be mapped directly. Structural fluctuations of bulk water occur in the time range between 5×10^–14 and 10^–12 s, the breaking and reformation of hydrogen bonds on a time scale of 10^–12 s. At the interface between the water hydration shell and DNA or phospholipids, specific hydration geometries exist for periods of some 10^–11 s with minor fluctuations of their structure. However, a long‐lived 'memory of water' is absent.     

Charge referencing issues in XPS of insulators as evidenced in the case of Al‐Si‐N thin films

Accurate charge referencing in XPS of insulating specimens is a delicate issue. This difficulty is illustrated in the case of Al‐Si‐N composite thin films deposited by reactive magnetron sputtering with variable composition from pure aluminum nitride to pure silicon nitride. The samples were mounted with Au‐coated metallic clamps. Argon sputter cleaning was required to remove a surface native oxide before analysis. For charge referencing implanted argon atoms from the sputter gas and a small amount of gold re‐deposited from the metallic clamps onto the specimen surface during sputter cleaning were evaluated. For the argon atoms, a surprisingly large chemical shift (～1 eV) and a significant peak broadening (0.6 eV) of the Ar 2p_3/2 photoelectron line were found with varying the Si content of the films. This could be related to chemical and structural changes of the Al‐Si‐N films. Hence implanted argon could not be used for charge referencing of Al‐Si‐N samples. In contrast to the implanted argon, the Au 4f_7/2 line width of the gold re‐deposited onto the sample surface did not depend on the Si content of Al‐Si‐N films. A constant energy shift (～1.2 eV) of the Au 4f_7/2 line as compared with bulk gold was, however, found, which was related to the size of gold particles formed on the insulating films. Therefore gold could be reliably used to study chemical shifts of sample‐relevant species in Al‐Si‐N films, but the absolute binding energies of Al 2p, Si 2p and N 1s photoelectrons could not be determined. Copyright © 2011 John Wiley & Sons, Ltd.     

Electron binding energies of aqueous alkali and halide ions: EUV photoelectron spectroscopy of liquid solutions and combined ab initio and molecular dynamics calculations

Photoelectron spectroscopy combined with the liquid microjet technique enables the direct probing of the electronic structure of aqueous solutions. We report measured and calculated lowest vertical electron binding energies of aqueous alkali cations and halide anions. In some cases, ejection from deeper electronic levels of the solute could be observed. Electron binding energies of a given aqueous ion are found to be independent of the counterion and the salt concentration. The experimental results are complemented by ab initio calculations, at the MP2 and CCSD(T) level, of the ionization energies of these prototype ions in the aqueous phase. The solvent effect was accounted for in the electronic structure calculations in two ways. An explicit inclusion of discrete water molecules using a set of snapshots from an equilibrium classical molecular dynamics simulations and a fractional charge representation of solvent molecules give good results for halide ions. The electron binding energies of alkali cations computed with this approach tend to be overestimated. On the other hand, the polarizable continuum model, which strictly provides adiabatic binding energies, performs well for the alkali cations but fails for the halides. Photon energies in the experiment were in the EUV region (typically 100 eV) for which the technique is probing the top layers of the liquid sample. Hence, the reported energies of aqueous ions are closely connected with both structures and chemical reactivity at the liquid interface, for example, in atmospheric aerosol particles, as well as fundamental bulk solvation properties.     

Low-energy electron distributions from the photoionization of liquid water: a sensitive test of electron mean free paths

The availability of accurate mean free paths for slow electrons (<50 eV) in water is central to the understanding of many electron-driven processes in aqueous solutions, but their determination poses major challenges to experiment and theory alike. Here, we describe a joint experimental and theoretical study demonstrating a novel approach for testing, and, in the future, refining such mean free paths. We report the development of Monte-Carlo electron-trajectory simulations including elastic and inelastic electron scattering, as well as energy loss and secondary-electron production to predict complete photoelectron spectra of liquid water. These simulations are compared to a new set of photoelectron spectra of a liquid-water microjet recorded over a broad range of photon energies in the extreme ultraviolet (20–57 eV). Several previously published sets of scattering parameters are investigated, providing direct and intuitive insights on how they influence the shape of the low-energy electron spectra. A pronounced sensitivity to the escape barrier is also demonstrated. These simulations considerably advance our understanding of the origin of the prominent low-energy electron distributions in photoelectron spectra of liquid water and clarify the influence of scattering parameters and the escape barrier on their shape. They moreover describe the reshaping and displacement of low-energy photoelectron bands caused by vibrationally inelastic scattering. Our work provides a quantitative basis for the interpretation of the complete photoelectron spectra of liquids and opens the path to fully predictive simulations of low-energy scattering in liquid water.

Our study reveals the detailed influence of elastic and inelastic mean-free paths on the complete photoelectron spectra of liquid water, including the low-energy electron distributions and the reshaping of the primary photoelectron bands.     

Phase Equilibria for Quaternary Systems of Water + Methanol or Ethanol + Ethyl Methyl Carbonate + Octane at <Emphasis Type="Italic">T</Emphasis> = 298.15 K

Liquid–liquid equilibrium (LLE) data were measured at 298.15 K and atmospheric pressure for the quaternary [water + methanol + ethyl methyl carbonate (EMC) + octane] and (water + ethanol + EMC + octane) systems, and relevant ternary (water + methanol + EMC), (water + ethanol + EMC), (water + ethanol + octane) and (water + EMC + octane) systems. The modified and extended UNIQUAC activity coefficient models were employed to correlate the experimental LLE results, which were also correlated by the Othmer–Tobias equation. The separation factors were calculated from the LLE data. The extracting selectivity of methanol and ethanol from aqueous solutions using pure EMC and mixed solvent (EMC + octane) was studied.     

Halogen Interactions in Halogenated Oxindoles: Crystallographic and Computational Investigations of Intermolecular Interactions

X-ray structural determinations and computational studies were used to investigate halogen interactions in two halogenated oxindoles. Comparative analyses of the interaction energy and the interaction properties were carried out for Br···Br, C-H···Br, C-H···O and N-H···O interactions. Employing Møller–Plesset second-order perturbation theory (MP2) and density functional theory (DFT), the basis set superposition error (BSSE) corrected interaction energy (E_int(BSSE)) was determined using a supramolecular approach. The E_int(BSSE) results were compared with interaction energies obtained by Quantum Theory of Atoms in Molecules (QTAIM)-based methods. Reduced Density Gradient (RDG), QTAIM and Natural bond orbital (NBO) calculations provided insight into possible pathways for the intermolecular interactions examined. Comparative analysis employing the electron density at the bond critical points (BCP) and molecular electrostatic potential (MEP) showed that the interaction energies and the relative orientations of the monomers in the dimers may in part be understood in light of charge redistribution in these two compounds.     

Simple and Equipment-Free Paper-Based Device for Determination of Mercury in Contaminated Soil

This work presents a simple and innovative protocol employing a microfluidic paper-based analytical device (µPAD) for equipment-free determination of mercury. In this method, mercury (II) forms an ionic-association complex of tetraiodomercurate (II) ion (HgI₄²⁻_(aq)) using a known excess amount of iodide. The residual iodide flows by capillary action into a second region of the paper where it is converted to iodine by pre-deposited iodate to liberate I_2(g) under acidic condition. Iodine vapor diffuses across the spacer region of the µPAD to form a purple colored of tri-iodide starch complex in a detection zone located in a separate layer of the µPAD. The digital image of the complex is analyzed using ImageJ software. The method has a linear calibration range of 50–350 mg L⁻¹ Hg with the detection limit of 20 mg L⁻¹. The method was successfully applied to the determination of mercury in contaminated soil and water samples which the results agreed well with the ICP-MS method. Three soil samples were highly contaminated with mercury above the acceptable WHO limits (0.05 mg kg⁻¹). To the best of our knowledge, this is the first colorimetric µPAD method that is applicable for soil samples including mercury contaminated soils from gold mining areas.     

Computational Analysis and Biological Activities of Oxyresveratrol Analogues,the Putative Cyclooxygenase-2 Inhibitors

Polyphenols are a large family of naturally occurring phytochemicals. Herein, oxyresveratrol was isolated from ethanolic crude extracts of Artocarpus lacucha Buch.-Ham., and chemically modified to derive its lipophilic analogues. Biological screening assays showed their inhibitory potency against cyclooxygenase-2 (COX-2) with very low cytotoxicity to the MRC-5 normal cell lines. At the catalytic site of COX-2, docking protocols with ChemPLP, GoldScore and AutoDock scoring functions were carried out to reveal hydrogen bonding interactions with key polar contacts and hydrophobic pi-interactions. For more accurate binding energetics, COX-2/ligand complexes at the binding region were computed in vacuo and implicit aqueous solvation using M06-2X density functional with 6-31G+(d,p) basis set. Our computational results confirmed that dihydrooxyresveratrol (4) is the putative inhibitor of human COX-2 with the highest inhibitory activity (IC₅₀ of 11.50 ± 1.54 µM) among studied non-fluorinated analogues for further lead optimization. Selective substitution of fluorine provides a stronger binding affinity; however, lowering the cytotoxicity of a fluorinated analogue to a normal cell is challenging. The consensus among biological activities, ChemPLP docking score and the binding energies computed at the quantum mechanical level is obviously helpful for identification of oxyresveratrol analogues as a putative anti-inflammatory agent.