Differential Capacitance of Ionic Liquid Confined between Metallic Interfaces.

We present here a detailed analysis of the electric double layer at the gold electrode/[BMIM][BF4] interface using a polarizable model for the electrode, based on our recent approach to include image charges [Geada et al. Nat. Commun. 2018, 9, 716]. A double bell (camel) shape is obtained for the differential capacitance, where the inclusion of metal polarization allows for a higher density of ions in the double layer, particularly around the maxima, thereby increasing the capacitance. The charging mechanism differs for the positive and negative electrodes, with counterion adsorption prevailing at the anode and co-ion desorption prevailing at the cathode. The charging mechanism is predominantly governed by the BF4 anions, serving as counterions and co-ions at the anode and cathode, respectively. Within the considered range of potentials, only minor changes are observed in the dynamical properties, specifically in the diffusion coefficients. Notably, it is interesting to observe that bulk properties are restored at a shorter distance from the gold surface in the case of the anode compared to the cathode.

Molecular dynamics simulations of electrified interfaces including the metal polarisation.

Understanding electrified interfaces requires an accurate description of the electric double layer which also takes into account the metal polarisation. Here we present a simple approach to the molecular dynamics simulation of electrified interfaces which combines fixed charges and a core-shell model for the description of the polarisable electron density on the metal electrode. The approach has been applied to the Au(111) surface in contact with a NaCl aqueous electrolyte solution in order to calculate the differential capacitance and to gain a detailed picture of the charging mechanism. Metal polarisation enhances the interfacial capacitance with a difference between the cathode and anode. In particular, we find that the influence of the metal polarisation on the electric double layer depends on the ion's solvation shell structure and, for the investigated interface, is more important at the cathode, where it modifies the sodium ion distribution.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment.

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Modeling of minimal systems based on ATP-Zn coordination for chemically fueled self-assembly.

Following nature's example, there is currently strong interest in using adenosine 5'-triphosphate (ATP) as a fuel for the self-assembly of functional materials with transient/non-equilibrium behaviours. These hold great promise for applications, e.g. in catalysis and drug delivery. In a recent seminal work [Maiti et al., Nat. Chem., 2016, 8, 725], binding of ATP to the metallosurfactant zinc hexadecyl-1,4,7-triazacyclononane ([ZnC16 TACN]2+) was exploited to produce ATP-fueled transient vesicles. Crucial to the complex formation is the ability of ATP to bind to the metal ion. As a first step to unveil the key elements underlying this process, we investigate the interaction of ATP with Zn2+ and with methyl-1,4,7-triazacyclononane ([ZnCH3 TACN]2+), using all-atom molecular dynamics simulations. The free energy landscape of the complex formation is sampled using well-tempered metadynamics with three collective variables, corresponding to the coordination numbers of Zn2+ with the oxygen atoms of the three phosphate groups. We find that the structure of the ternary complex is controlled by direct triphosphate coordination to zinc, with a minor role played by the interactions between ATP and CH3 TACN which, however, may be important for the build-up of supramolecular assemblies.

Charging of Dielectric Surfaces in Contact with Aqueous Electrolytes─the Influence of CO2.

The charge state of dielectric surfaces in aqueous environments is of fundamental and technological importance. Here, we study the influence of dissolved molecular CO2 on the charging of three chemically different surfaces (SiO2, Polystyrene, Perfluorooctadecyltrichlorosilane). We determine their charge state from electrokinetic experiments. We compare an ideal, CO2-free reference system to a system equilibrated against ambient CO2 conditions. In the reference system, the salt-dependent decrease of the magnitudes of ζ-potentials follows the expectations for a constant charge scenario. In the presence of CO2, the starting potential is lower by some 50%. The following salt-dependent decrease is weakened for SiO2 and inverted for the organic surfaces. We show that screening and pH-driven charge regulation alone cannot explain the observed effects. As an additional cause, we tentatively suggest dielectric regulation of surface charges due to a diffusively adsorbed thin layer of molecular CO2. The formation of such a dynamic layer, even at the hydrophilic and partially ionized silica surfaces, is supported by a minimal theoretical model and results from molecular simulations.

Lower degree of dissociation of pyruvic acid at water surfaces than in bulk.

Understanding the acid/base behavior of environmentally relevant organic acids is of key relevance for accurate climate modelling. Here we investigate the effect of pH on the (de)protonation state of pyruvic acid at the air-water interface and in bulk by using the analytical techniques surface-specific vibrational sum frequency generation and attenuated total reflection spectroscopy. To provide a molecular interpretation of the observed behavior, simulations are carried out using a free energy perturbation approach in combination with electronic structure-based molecular dynamics. In both the experimental and theoretical results we observe that the protonated form of pyruvic acid is preferred at the air-water interface. The increased proton affinity is the result of the specific microsolvation at the interface.

Effects of shear flow on the structure and dynamics of ionic liquids in a metallic nanoconfinement.

It has been shown that a weak shear can induce crystallisation in a disordered, glassy state. In this study, we use molecular dynamics simulations in order to investigate the out-of-equilibrium properties of [BMIM][BF4] confined between metal slabs. In particular, we want to understand the extent to which the shear flow modifies the interfacial properties. In particular, the questions we address here are (i) is the shear able to promote the crystalline phase in [BMIM][BF4]? (ii) Can, as a consequence of shear flow, a solid-like layer develop at the interface with a metallic surface? (iii) What are the tribological properties of nanoconfined [BMIM][BF4]? We find that the system behaves quite differently from the ideal linear Couette flow. Indeed, the portion of fluid closer to the shearing slabs behaves as a disordered, solid-like layer, which, under the investigated conditions extends to a few nanometres. The linear velocity regime is only recovered in the central region of the ionic liquid slab. The formation of such a solid-like glassy rather than crystalline layer is in agreement with recent mechanical impedance measurements performed on nano-confined ionic liquids.

On the origin of controlled anisotropic growth of monodisperse gold nanobipyramids.

We elucidate the crucial role of the cetyl trimethylammonium bromide (CTAB) surfactant in the anisotropic growth mechanism of gold nano-bipyramids, nano-objects with remarkable optical properties and high tunability. Atomistic molecular dynamics simulations predict different surface coverages of the CTAB (positively charged) heads and their (bromide) counterions as function of the gold exposed surfaces. High concentration of CTAB surfactant promotes formation of gold nanograins in solution that work as precursors for the smooth anisotropic growth of more elongated nano-bipyramidal objects. Nanobipyramids feature higher index facets with respect to nanorods, allowing higher CTAB coverages that stabilize their formation and leading to narrower inter-micelles channels that smooth down their anisotropic growth. Absorption spectroscopy and scanning electron microscopy confirmed the formation of nanograins and demonstrated the importance of surfactant concentration on driving the growth towards nano-bipyramids rather than nanorods. The outcome explains the formation of the monodisperse bipyramidal nano-objects, the origin of their controlled shapes and sizes along with their remarkable stability.

Role of pH in the synthesis and growth of gold nanoparticles using L-asparagine: a combined experimental and simulation study.

The use of biomolecules as capping and reducing agents in the synthesis of metallic nanoparticles constitutes a promising framework to achieve desired functional properties with minimal toxicity. The system's complexity and the large number of variables involved represent a challenge for theoretical and experimental investigations aiming at devising precise synthesis protocols. In this work, we use L-asparagine (Asn), an amino acid building block of large biomolecular systems, to synthesise gold nanoparticles (AuNPs) in aqueous solution at controlled pH. The use of Asn offers a primary system that allows us to understand the role of biomolecules in synthesising metallic nanoparticles. Our results indicate that AuNPs synthesised in acidic (pH 6) and basic (pH 9) environments exhibit somewhat different morphologies. We investigate these AuNPs via Raman scattering experiments and classical molecular dynamics simulations of zwitterionic and anionic Asn states adsorbing on (111)-, (100)-, (110)-, and (311)-oriented gold surfaces. A combined analysis suggests that the underlying mechanism controlling AuNPs geometry correlates with amine's preferential adsorption over ammonium groups, enhanced upon increasing pH. Our simulations reveal that Asn (both zwitterionic and anionic) adsorption on gold (111) is essentially different from adsorption on more open surfaces. Water molecules strongly interact with the gold face-centred-cubic lattice and create traps, on the more open surfaces, that prevent the Asn from diffusing. These results indicate that pH is a relevant parameter in green-synthesis protocols with the capability to control the nanoparticle's geometry, and pave the way to computational studies exploring the effect of water monolayers on the adsorption of small molecules on wet gold surfaces.

Role of image charges in ionic liquid confined between metallic interfaces.

The peculiar properties of ionic liquids in confinement have not only become essential for energy storage, catalysis and tribology, but still pose fundamental questions. Recently, an anomalous liquid-solid phase transition has been observed in atomic force microscopy experiments for 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]), the transition being more pronounced for metallic surfaces. Image charges have been suggested as the key element driving the anomalous freezing. Using atomistic molecular dynamics simulations, we investigate the impact of image charges on structure, dynamics and thermodynamics of [BMIM][BF4] confined between gold electrodes. Our results not only unveil a minor role played by the metal polarisation, but also provide a novel description of the interfacial layer. Although no diffuse layer can be defined in terms of the electrostatic potential, long range effects are clearly visible in the dynamical properties up to 10 nanometers away from the surface, and are expected to influence viscous forces in the experiments.

Surface Charges at the CaF2 /Water Interface Allow Very Fast Intermolecular Vibrational-Energy Transfer.

We investigate the dynamics of water in contact with solid calcium fluoride, where at low pH, localized charges can develop upon fluorite dissolution. We use 2D surface-specific vibrational spectroscopy to quantify the heterogeneity of the interfacial water (D2 O) molecules and provide information about the sub-picosecond vibrational-energy-relaxation dynamics at the buried solid/liquid interface. We find that strongly H-bonded OD groups, with a vibrational frequency below 2500 cm-1 , display very rapid spectral diffusion and vibrational relaxation; for weakly H-bonded OD groups, above 2500 cm-1 , the dynamics slows down substantially. Atomistic simulations based on electronic-structure theory reveal the molecular origin of energy transport through the local H-bond network. We conclude that strongly oriented H-bonded water molecules in the adsorbed layer, whose orientation is pinned by the localized charge defects, can exchange vibrational energy very rapidly due to the strong collective dipole, compensating for a partially missing solvation shell.

The puzzling issue of silica toxicity: are silanols bridging the gaps between surface states and pathogenicity?

BACKGROUND: Silica continues to represent an intriguing topic of fundamental and applied research across various scientific fields, from geology to physics, chemistry, cell biology, and particle toxicology. The pathogenic activity of silica is variable, depending on the physico-chemical features of the particles. In the last 50 years, crystallinity and capacity to generate free radicals have been recognized as relevant features for silica toxicity. The 'surface' also plays an important role in silica toxicity, but this term has often been used in a very general way, without defining which properties of the surface are actually driving toxicity. How the chemical features (e.g., silanols and siloxanes) and configuration of the silica surface can trigger toxic responses remains incompletely understood. MAIN BODY: Recent developments in surface chemistry, cell biology and toxicology provide new avenues to improve our understanding of the molecular mechanisms of the adverse responses to silica particles. New physico-chemical methods can finely characterize and quantify silanols at the surface of silica particles. Advanced computational modelling and atomic force microscopy offer unique opportunities to explore the intimate interactions between silica surface and membrane models or cells. In recent years, interdisciplinary research, using these tools, has built increasing evidence that surface silanols are critical determinants of the interaction between silica particles and biomolecules, membranes, cell systems, or animal models. It also has become clear that silanol configuration, and eventually biological responses, can be affected by impurities within the crystal structure, or coatings covering the particle surface. The discovery of new molecular targets of crystalline as well as amorphous silica particles in the immune system and in epithelial lung cells represents new possible toxicity pathways. Cellular recognition systems that detect specific features of the surface of silica particles have been identified. CONCLUSIONS: Interdisciplinary research bridging surface chemistry to toxicology is progressively solving the puzzling issue of the variable toxicity of silica. Further interdisciplinary research is ongoing to elucidate the intimate mechanisms of silica pathogenicity, to possibly mitigate or reduce surface reactivity.

Heterogeneous Interactions between Gas-Phase Pyruvic Acid and Hydroxylated Silica Surfaces: A Combined Experimental and Theoretical Study.

The adsorption of gas-phase pyruvic acid (CH3COCOOH) on hydroxylated silica particles has been investigated at 296 K using transmission Fourier transform infrared (FTIR) spectroscopy and theoretical simulations. Under dry conditions (<1% relative humidity, RH), both the trans-cis (Tc) and trans-trans (Tt) pyruvic acid conformers are observed on the surface as well as the (hydrogen bonded) pyruvic acid dimer. The detailed surface interactions were further understood through ab initio molecular dynamics simulations. Under higher relative humidity conditions (above 10% RH), adsorbed water competes for surface adsorption sites. Adsorbed water is also observed to change the relative populations of the different adsorbed pyruvic acid configurations. Overall, this study provides valuable insights into the interaction of pyruvic acid with hydroxylated silica surfaces on the molecular level from both experimental and theoretical analyses. Furthermore, these results highlight the importance of the environment (relative humidity and coadsorbed water) in the adsorption, partitioning, and configurations of pyruvic acid at the surface.

Dynamical heterogeneities of rotational motion in room temperature ionic liquids evidenced by molecular dynamics simulations.

Room temperature ionic liquids (RTILs) have been shown to exhibit spatial heterogeneity or structural heterogeneity in the sense that they form hydrophobic and ionic domains. Yet studies of the relationship between this structural heterogeneity and the ∼picosecond motion of the molecular constituents remain limited. In order to obtain insight into the time scales relevant to this structural heterogeneity, we perform molecular dynamics simulations of a series of RTILs. To investigate the relationship between the structures, i.e., the presence of hydrophobic and ionic domains, and the dynamics, we gradually increase the size of the hydrophobic part of the cation from ethylammonium nitrate (EAN), via propylammonium nitrate (PAN), to butylammonium nitrate (BAN). The two ends of the organic cation, namely, the charged Nhead-H group and the hydrophobic Ctail-H group, exhibit rotational dynamics on different time scales, evidencing dynamical heterogeneity. The dynamics of the Nhead-H group is slower because of the strong coulombic interaction with the nitrate counter-ionic anions, while the dynamics of the Ctail-H group is faster because of the weaker van der Waals interaction with the surrounding atoms. In particular, the rotation of the Nhead-H group slows down with increasing cationic chain length, while the rotation of the Ctail-H group shows little dependence on the cationic chain length, manifesting that the dynamical heterogeneity is enhanced with a longer cationic chain. The slowdown of the Nhead-H group with increasing cationic chain length is associated with a lower number of nitrate anions near the Nhead-H group, which presumably results in the increase of the energy barrier for the rotation. The sensitivity of the Nhead-H rotation to the number of surrounding nitrate anions, in conjunction with the varying number of nitrate anions, gives rise to a broad distribution of Nhead-H reorientation times. Our results suggest that the asymmetry of the cations and the larger excluded volume for longer cationic chain are important for both the structural heterogeneity and the dynamical heterogeneities. The observed dynamical heterogeneities may affect the rates of chemical reactions depending on where the reactants are solvated in ionic liquids and provide an additional guideline for the design of RTILs as solvents.

A set-up for simultaneous measurement of second harmonic generation and streaming potential and some test applications.

We present a measurement cell that allows simultaneous measurement of second harmonic generation (SHG) and streaming potential (SP) at mineral-water interfaces with flat specimen that are suitable for non-linear optical (NLO) studies. The set-up directly yields SHG data for the interface of interest and can also be used to obtain information concerning the influence of flow on NLO signals from that interface. The streaming potential is at present measured against a reference substrate (PTFE). The properties of this inert reference can be independently determined for the same conditions. With the new cell, for the first time the SHG signal and the SP for flat surfaces have been simultaneously measured on the same surface. This can in turn be used to unambiguously relate the two observations for identical solution composition. The SHG test of the cell with a fluorite sample confirmed previously observed differences in NLO signal under flow vs. no flow conditions in sum frequency generation (SFG) investigations. As a second test surface, an inert ("hydrophobic") OTS covered sapphire-c electrolyte interface was studied to verify the zeta-potential measurements with the new cell. For this system we obtained combined zeta-potential/SHG data in the vicinity of the point of zero charge, which were found to be proportional to each other as expected. Furthermore, on the accessible time scales of the SHG measurements no effects of flow, flow velocity and stopped flow occurred on the interfacial water structure. This insensitivity to flow for the inert surface was corroborated by concomitant molecular dynamics simulations. Finally, the set-up was used for simultaneous measurements of the two properties as a function of pH in automated titrations with an oxidic surface. Different polarization combinations obtained in two separate titrations, yielded clearly different SHG data, while under identical conditions zeta-potentials were exactly reproduced. The polarization combination that is characteristic for dipoles perpendicular to the surface scaled with the zeta-potentials over the pH-range studied, while the other did not. The work provides an advanced approach for investigating liquid/surface interactions which play a major role in our environment. The set-up can be upgraded for SFG studies, which will allow more detailed studies on the chemistry and the water structure at a given interface, but also the combined study of specific adsorption including kinetics in combination with electrokinetics. Such investigations are crucial for the basic understanding of many environmental processes from aquatic to atmospheric systems.

Atypical titration curves for GaAl12 Keggin-ions explained by a joint experimental and simulation approach.

Although they have been widely used as models for oxide surfaces, the deprotonation behaviors of the Keggin-ions (MeAl127+) and typical oxide surfaces are very different. On Keggin-ions, the deprotonation occurs over a very narrow pH range at odds with the broad charging curve of larger oxide surfaces. Depending on the Me concentration, the deprotonation curve levels off sooner (high Me concentration) or later (for low Me concentration). The leveling off shows the onset of aggregation before which the Keggin-ions are present as individual units. We show that the atypical titration data previously observed for some GaAl12 solutions in comparison to the originally reported data can be explained by the presence of Ga2Al11 ions. The pKa value of aquo-groups bound to octahedral Ga was determined from ab initio molecular dynamics simulations relative to the pure GaAl12 ions. Using these results within a surface complexation model, the onset of deprotonation of the crude solution is surprisingly well predicted and the ratio between the different species is estimated to be in the proportion 20 (Ga2Al11) : 20 (Al13) : 60 (GaAl12).

A Microscopic Interpretation of Pump-Probe Vibrational Spectroscopy Using Ab Initio Molecular Dynamics.

What happens when extra vibrational energy is added to water? Using nonequilibrium molecular dynamics simulations, also including the full electronic structure, and novel descriptors, based on projected vibrational density of states, we are able to follow the flow of excess vibrational energy from the excited stretching and bending modes. We find that the energy relaxation, mostly mediated by a stretching-stretching coupling in the first solvation shell, is highly heterogeneous and strongly depends on the local environment, where a strong hydrogen bond network can transport energy with a time scale of 200 fs, whereas a weaker network can slow down the transport by a factor 2-3.

Increased Acid Dissociation at the Quartz/Water Interface.

As shown by a quite significant amount of literature, acids at the water surface tend to be "less" acid, meaning that their associated form is favored over the conjugated base. What happens at the solid/liquid interface? In the case of the silica/water interface, we show how the acidity of adsorbed molecules can instead increase. Using a free energy perturbation approach in combination with electronic structure-based molecular dynamics simulations, we show how the acidity of pyruvic acid at the quartz/water interface is increased by almost two units. Such increased acidity is the result of the specific microsolvation at the interface and, in particular, of the stabilization of the deprotonated form by the silanols on the quartz surface and the special interfacial water layer.

Insight into induced charges at metal surfaces and biointerfaces using a polarizable Lennard-Jones potential.

Metallic nanostructures have become popular for applications in therapeutics, catalysts, imaging, and gene delivery. Molecular dynamics simulations are gaining influence to predict nanostructure assembly and performance; however, instantaneous polarization effects due to induced charges in the free electron gas are not routinely included. Here we present a simple, compatible, and accurate polarizable potential for gold that consists of a Lennard-Jones potential and a harmonically coupled core-shell charge pair for every metal atom. The model reproduces the classical image potential of adsorbed ions as well as surface, bulk, and aqueous interfacial properties in excellent agreement with experiment. Induced charges affect the adsorption of ions onto gold surfaces in the gas phase at a strength similar to chemical bonds while ions and charged peptides in solution are influenced at a strength similar to intermolecular bonds. The proposed model can be applied to complex gold interfaces, electrode processes, and extended to other metals.