The Surface of Ionic Liquids in Water: From an Ionic Tug of War to a Quasi-Ordered Two-Dimensional Layer.

The sustainable development encompasses the search for new materials for energy storage, gas capture, separation, and solvents in industrial processes that can substitute conventional ones in an efficient and clean manner. Ionic liquids (ILs) emerged and have been advanced as alternative materials for such applications, but an obstacle is their hygroscopicity and the effects on their physical properties in the presence of humidity. Several industrial processes depend on the aqueous interfacial properties, and the main focus of this work is the water/IL interface. The behavior of the aqueous ionic liquids at the water-vacuum interface is representative for their water interfacial properties. Using X-ray photoelectron spectroscopy in combination with molecular dynamics simulations we investigate four aqueous IL systems, and provide molecular level insight on the interfacial behaviour of the ionic liquids, such as ion-pair formation, orientation and surface concentration. We find that ionic liquids containing a chloride anion have a lowered surface enrichment due to the low surface propensity of chloride. In contrast, the ionic liquids containing a bistriflimide anion are extremely surface-enriched due to cooperative surface propensity between the cations and anions, forming a two-dimensional ionic liquid on the water surface at low concentrations.

An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation.

Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.

Microscopic origins of conductivity in molten salts unraveled by computer simulations.

Molten salts are crucial materials in energy applications, such as batteries, thermal energy storage systems or concentrated solar power plants. Still, the determination and interpretation of basic physico-chemical properties like ionic conductivity, mobilities and transference numbers cause debate. Here, we explore a method for determination of ionic electrical mobilities based on non-equilibrium computer simulations. Partial conductivities are then determined as a function of system composition and temperature from simulations of molten LiFαClßIγ (with α + ß + γ = 1). High conductivity does not necessarily coincide with high Li+ mobility for molten LiFαClßIγ systems at a given temperature. In salt mixtures, the lighter anions on average drift along with Li+ towards the negative electrode when applying an electric field and only the heavier anions move towards the positive electrode. In conclusion, the microscopic origin of conductivity in molten salts is unraveled here based on accurate ionic electrical mobilities and an analysis of the local structure and kinetics of the materials.

Systematically improved melting point prediction: a detailed physical simulation model is required.

Accurate prediction of fundamental properties such as melting points using direct physical simulation is challenging. Here, we investigate the melting point (Tm) of alkali halides that are often considered to be the simplest category of salts. Popular force fields that have been examined for this task leave considerable room for improvement. Recently we introduced a new force field for alkali halides (WBK) as part of the Alexandria project, featuring explicit polarisation and distributed charges. This new force field significantly improves the prediction of a large set of physicochemical properties and in this contribution we show that the same is valid for the prediction of Tm. For reference, we calculated Tm using a non-polarisable force field by Joung and Cheatham (JC), and compare our results to existing literature data on the widely used Tosi-Fumi (TF) parameters. In contrast to the predictions of the WBK model, the JC force field consistently overestimates the experimental Tm, while the accuracy of the TF model strongly depends on the investigated salt. Our results show that the inclusion of more realistic physics into a force field opens up the possibility to accurately describe many physicochemical properties over a large range of temperatures, even including phase transitions.

Molten alkali halides - temperature dependence of structure, dynamics and thermodynamics.

The renewed interest in molten salts in the energy industry fuels the need of a thorough understanding of their physicochemical properties. Alkali halide melts are perhaps the simplest ionic liquids, but they are used as electrolytes in batteries or for thermal energy storage. Although their structure is considered to be well documented and understood, a systematic evaluation of experimental structural data reveals significant discrepancies, while there is only limited experimental information on dynamic properties. Here, we investigate structure, dynamics and thermodynamic properties of pure alkali halide melts using state-of-the-art simulation models at different temperatures. The simulations provide a consistent picture of the structure of alkali halide melts with coordination numbers that lie in between experimental numbers. The simulations reveal a strengthening of the cation-anion bonds with increasing temperature that, somewhat counter-intuitively, coincides with faster dynamics in the melts. The thermodynamic analysis unveils that structure breaking proceeds on the picosecond timescale through an associative substitution mechanism as signified by a negative entropy of activation. The results on ion pair lifetimes contribute to an improved understanding of the microscopic origin of dynamical properties, such as e.g. conductivity of salt melts. The structural analysis provided here contributes to a more coherent picture of the coordination numbers in alkali halides than what is currently available from experimental data.

Strong enrichment of atmospherically relevant organic ions at the aqueous interface: the role of ion pairing and cooperative effects.

Surface affinity, orientation and ion pairing are investigated in mixed and single solute systems of aqueous sodium hexanoate and hexylammonium chloride. The surface sensitive X-ray photoelectron spectroscopy technique has been used to acquire the experimental results, while the computational data have been calculated using molecular dynamics simulations. By comparing the single solute solutions with the mixed one, we observe a non-linear surface enrichment and reorientation of the organic ions with their alkyl chains pointing out of the aqueous surface. We ascribe this effect to ion paring between the charged functional groups on the respective organic ion and hydrophobic expulsion of the alkyl chains from the surface in combination with van der Waals interactions between the alkyl chains. These cooperative effects lead to a substantial surface enrichment of organic ions, with consequences for aerosol surface properties.

Phase-Transferable Force Field for Alkali Halides.

A longstanding goal of computational chemistry is to predict the state of materials in all phases with a single model. This is particularly relevant for materials that are difficult or dangerous to handle or compounds that have not yet been created. Progress toward this goal has been limited, as most work has concentrated on just one phase, often determined by particular applications. In the framework of the development of the Alexandria force field, we present here new polarizable force fields for alkali halides with Gaussian charge distributions for molecular dynamics simulations. We explore different descriptions of the van der Waals interaction, like the commonly applied 12-6 Lennard-Jones (LJ), and compare it to "softer" ones, such as the 8-6 LJ, Buckingham, and a modified Buckingham potential. Our results for physicochemical properties of the gas, liquid, and solid phases of alkali halides are compared to experimental data and calculations with reference polarizable and nonpolarizable force fields. The new polarizable force field that employs a modified Buckingham potential predicts the tested properties for gas, liquid, and solid phases with a very good accuracy. In contrast to reference force fields, this model reproduces the correct crystal structures for all alkali halides at low and high temperature. Seeing that experiments with molten salts may be tedious due to high temperatures and their corrosive nature, the models presented here can contribute significantly to our understanding of alkali halides in general and melts in particular.

Shifted equilibria of organic acids and bases in the aqueous surface region.

Acid-base equilibria of carboxylic acids and alkyl amines in the aqueous surface region were studied using surface-sensitive X-ray photoelectron spectroscopy and molecular dynamics simulations. Solutions of these organic compounds were examined as a function of pH, concentration and chain length to investigate the distribution of acid and base form in the surface region as compared to the aqueous bulk. Results from these experiments show that the neutral forms of the studied acid-base pairs are strongly enriched in the aqueous surface region. Moreover, we show that for species with at least four carbon atoms in their alkyl-chain, their charged forms are also found to be abundant in the surface region. Using a combination of XPS and MD results, a model is proposed that effectively describes the surface composition. Resulting absolute surface concentration estimations show clearly that the total organic mole fractions in the surface region change drastically as a function of solution pH. The origin of the observed surface phenomena, hydronium/hydroxide concentrations in the aqueous surface region and why standard chemical equations, used to describe equilibria in dilute bulk solution are not valid in the aqueous surface region, are discussed in detail. The reported results are of considerable importance especially for the detailed understanding of properties of small aqueous droplets that can be found in the atmosphere.

Anomalous surface behavior of hydrated guanidinium ions due to ion pairing.

Surface affinity of aqueous guanidinium chloride (GdmCl) is compared to that of aqueous tetrapropylammonium chloride (TPACl) upon addition of sodium chloride (NaCl) or disodium sulfate (Na2SO4). The experimental results have been acquired using the surface sensitive technique X-ray photoelectron spectroscopy on a liquid jet. Molecular dynamics simulations have been used to produce radial distribution functions and surface density plots. The surface affinities of both TPA+ and Gdm+ increase upon adding NaCl to the solution. With the addition of Na2SO4, the surface affinity of TPA+ increases, while that of Gdm+ decreases. From the results of MD simulations it is seen that Gdm+ and SO42- ions form pairs. This finding can be used to explain the decreased surface affinity of Gdm+ when co-dissolved with SO42- ions. Since SO42- ions avoid the surface due to the double charge and strong water interaction, the Gdm+-SO42- ion pair resides deeper in the solutions' bulk than the Gdm+ ions. Since TPA+ does not form ion pairs with SO42-, the TPA+ ions are instead enriched at the surface.

Ethanol Solvation in Water Studied on a Molecular Scale by Photoelectron Spectroscopy.

Because of the amphiphilic properties of alcohols, hydrophobic hydration is important in the alcohol-water system. In the present paper we employ X-ray photoelectron spectroscopy (XPS) to investigate the bulk and surface molecular structure of ethanol-water mixtures from 0.2 to 95 mol %. The observed XPS binding energy splitting between the methyl C 1s and hydroxymethyl C 1s groups (BES_[CH3-CH2OH]) as a function of the ethanol molar percentage can be divided into different regions: one below 35 mol % with higher values (about 1.53 eV) and one starting at 60 mol % up to 95 mol % with 1.49 eV as an average value. The chemical shifts agree with previous quantum mechanics/molecular mechanics (QM/MM) calculations [ Löytynoja , T. ; J. Phys. Chem. B 2014 , 118 , 13217 ]. According to these calculations, the BES_[CH3-CH2OH] is related to the number of hydrogen bonds between the ethanol and the surrounding molecules. As the ethanol concentration increases, the average number of hydrogen bonds decreases from 2.5 for water-rich mixtures to 2 for pure ethanol. We give an interpretation for this behavior based on how the hydrogen bonds are distributed according to the mixing ratio. Since our experimental data are surface sensitive, we propose that this effect may also be manifested at the interface. From the ratio between the XPS C 1s core lines intensities we infer that below 20 mol % the ethanol molecules have their hydroxyl groups more hydrated and possibly facing the solution's bulk. Between 0.1 and 14 mol %, we show the formation of an ethanol monolayer at approximately 2 mol %. Several parameters are derived for the surface region at monolayer coverage.

Surface Partitioning in Organic-Inorganic Mixtures Contributes to the Size-Dependence of the Phase-State of Atmospheric Nanoparticles.

Atmospheric particulate matter is one of the main factors governing the Earth's radiative budget, but its exact effects on the global climate are still uncertain. Knowledge on the molecular-scale surface phenomena as well as interactions between atmospheric organic and inorganic compounds is necessary for understanding the role of airborne nanoparticles in the Earth system. In this work, surface composition of aqueous model systems containing succinic acid and sodium chloride or ammonium sulfate is determined using a novel approach combining X-ray photoelectron spectroscopy, surface tension measurements and thermodynamic modeling. It is shown that succinic acid molecules are accumulated in the surface, yielding a 10-fold surface concentration as compared with the bulk for saturated succinic acid solutions. Inorganic salts further enhance this enrichment due to competition for hydration in the bulk. The surface compositions for various mixtures are parametrized to yield generalizable results and used to explain changes in surface tension. The enhanced surface partitioning implies an increased maximum solubility of organic compounds in atmospheric nanoparticles. The results can explain observations of size-dependent phase-state of atmospheric nanoparticles, suggesting that these particles can display drastically different behavior than predicted by bulk properties only.

Water at Interfaces.

The interfaces of neat water and aqueous solutions play a prominent role in many technological processes and in the environment. Examples of aqueous interfaces are ultrathin water films that cover most hydrophilic surfaces under ambient relative humidities, the liquid/solid interface which drives many electrochemical reactions, and the liquid/vapor interface, which governs the uptake and release of trace gases by the oceans and cloud droplets. In this article we review some of the recent experimental and theoretical advances in our knowledge of the properties of aqueous interfaces and discuss open questions and gaps in our understanding.

Acid-base speciation of carboxylate ions in the surface region of aqueous solutions in the presence of ammonium and aminium ions.

The acid-base speciation of surface-active carboxylate ions in the surface region of aqueous solutions was studied with synchrotron-radiation-based photoelectron spectroscopy. The protonated form was found at an extraordinarily large fraction compared to that expected from the bulk pH. When adding salts containing the weak acid NH4(+) to the solution, the fraction of the acidic form at the surface increases, and to a much greater extent than expected from the bulk pH of the solution. We show that ammonium ions also are overrepresented in the surface region, and propose that the interaction between the surface-active anionic carboxylates and cationic ammonium ions creates a carboxylate-ammonium bilayer close to the surface, which increases the probability of the protonation of the carboxylate ions. By comparing the situation when a salt of the less volatile amine diethanolamine is used, we also show that the observed evaporation of ammonia that occurs after such an event only affects the equilibrium marginally.

Electron Beam-Induced Writing of Nanoscale Iron Wires on a Functional Metal Oxide.

Electron beam-induced surface activation (EBISA) has been used to grow wires of iron on rutile TiO2(110)-(1 × 1) in ultrahigh vacuum. The wires have a width down to ∼20 nm and hence have potential utility as interconnects on this dielectric substrate. Wire formation was achieved using an electron beam from a scanning electron microscope to activate the surface, which was subsequently exposed to Fe(CO)5. On the basis of scanning tunneling microscopy and Auger electron spectroscopy measurements, the activation mechanism involves electron beam-induced surface reduction and restructuring.

Electron beam induced surface activation of ultrathin porphyrin layers on Ag(111).

We demonstrate how a focused electron beam can be used to chemically activate porphyrin layers on Ag(111) such that they become locally reactive toward the decomposition of iron pentacarbonyl, Fe(CO)5. This finding considerably expands the scope of electron beam induced surface activation (EBISA) and also has implications for electron beam induced deposition (EBID). The influence of the porphyrin layer thickness on both processes is studied in detail using scanning tunneling microscopy (STM) and scanning electron microscopy (SEM) as well as Auger electron spectroscopy (AES) and scanning Auger microscopy (SAM). While a closed monolayer of porphyrin molecules does exhibit some activity toward Fe(CO)5 decomposition after electron irradiation, a growth enhancement is found for bi- and multilayer films. This is attributed to a partial quenching of activated centers in the first layer due to the close proximity of the silver substrate. In addition, we demonstrate that the catalytic decomposition of gaseous Fe(CO)5 on Ag(111) can be effectively inhibited by introducing a densely packed monolayer of 2H-tetraphenylporphyrin (2HTPP) molecules.

Defects in oxygen-depleted titanate nanostructures.

The identification of defects and their controlled generation in titanate nanostructures is a key to their successful application in photoelectronic devices. We comprehensively explored the effect of vacuum annealing on morphology and composition of Na(2)Ti(3)O(7) nanowires and protonated H(2)Ti(3)O(7) nanoscrolls using a combination of scanning electron microscopy, Auger and Fourier-transform infrared (FT-IR) spectroscopy, as well as ab initio density functional theory (DFT) calculations. The observation that H(2)Ti(3)O(7) nanoscrolls are more susceptible to electronic reduction and annealing-induced n-type doping than Na(2)Ti(3)O(7) nanowires is attributed to the position of the conduction band minimum. It is close to the vacuum level and, thus, favors the Fermi level-induced compensation of donor states by cation vacancies. In agreement with theoretical predictions that suggest similar formation energies for oxygen and sodium vacancies, we experimentally observed the annealing induced depletion of sodium from the surface of the nanowires.

Generation of clean iron nanocrystals on an ultra-thin SiO(x) film on Si(001).

Upon exposure to Fe(CO)(5), the formation of pure cubic Fe nanocrystals with dimensions up to ~75 nm is reported on ultra-thin SiO(x) films (thickness ≈ 0.5 nm) on Si(001), which have been prepared in situ under UHV conditions. The active centers for initial decomposition of Fe(CO)(5) resulting in the growth of the Fe clusters are proposed to be SiO sites. After nucleation at these sites, further crystal growth is observed due to autocatalytic dissociation of Fe(CO)(5) at room temperature. The density of the Fe clusters can be increased by irradiating the surface with a focused electron beam (15 keV) prior to gas exposure. The formation of the active SiO sites upon electron irradiation is attributed to oxygen desorption via the Knotek-Feibelman mechanism.

Tetraphenylporphyrin picks up zinc atoms from a silver surface.

We demonstrate that adsorbed meso-tetraphenylporphyrin molecules can coordinate Zn atoms that are pre-deposited on an Ag(111) surface, forming a complex that is identical to directly deposited tetraphenylporphyrinato-zinc(II); this reaction, which we studied with XPS, is the first example of an oxidative dissolution of a metal by a large organic ligand under ultrahigh vacuum conditions.