Throwing Tetrahedral Dice

Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies

Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.

Probing n-Octyl α-d-Glycosides Using Deuterated Water in the Lyotropic Phase by Deuterium NMR

The lyotropic phase behavior of four common and easily accessible glycosides, n-octyl α-d-glycosides, namely, α-Glc-OC₈, α-Man-OC₈, α-Gal-OC₈, and α-Xyl-OC₈, was investigated. The presence of normal hexagonal (H_I), bicontinuous cubic (V_I), and lamellar (L_α) phases in α-Glc-OC₈ and α-Man-OC₈ including their phase diagrams in water reported previously was verified by deuterium nuclear magnetic resonance (²H NMR), via monitoring the D₂O spectra. Additionally, the partial binary phase diagrams and the liquid crystal structures formed by α-Gal-OC₈ and α-Xyl-OC₈ in D₂O were constructed and confirmed using small- and wide-angle X-ray scattering and ²H NMR. The average number of bound water molecules (n_b) per headgroup in the L_α phase was determined by the systematic measurement of the quadrupolar splitting of D₂O over a wide range of molar ratio values (glycoside/D₂O), especially at high glucoside composition. The number of bound water molecules bound to the headgroup was found to be around 1.5-2.0 for glucoside, mannoside, and galactoside, all of which possesses four OH groups. In the case of xyloside, which has only three OH groups, the bound water content is ∼2.0. Our findings confirmed that the bound water content of all n-octyl α-d-glycosides studied is lower compared to the number of possible hydrogen bonding sites possibly due to the fact that most of the OH groups are involved in intralayer interaction that holds the lipid assembly together.

How many water molecules are needed to solvate one?

Many efforts undertaken to study the solvation process have led to general theories that may describe mean properties, but are unable to provide a detailed understanding at the molecular level. Remarkably, the basic question of how many solvent molecules are necessary to solvate one solute molecule is still open. By exploring several water aggregates of increasing complexity, in this contribution we employ semiclassical spectroscopy to determine on quantum dynamical grounds the minimal network of surrounding water molecules to make the central one display the same vibrational features of liquid water. We find out that double-acceptor double-donor tetrahedral coordination constituting the standard picture is necessary but not sufficient, and that particular care must be reserved for the quantum description of the combination band due to the coupling of the central monomer bending mode with network librations. It is actually our ability to investigate the combination band with a quantum-derived approach that allows us to answer the titular question. The minimal structure eventually responsible for proper solvation is made of a total of 21 water molecules and includes two complete solvation shells, of which the whole first one is tetrahedrally coordinated to the central molecule.

How quantum spectroscopic simulations can explain water solvation by comparison with experimental spectra.

Molecular Dynamics Simulation of Calcium-Silicate-Hydrate for Nano-Engineered Cement Composites-A Review

With the continuous research efforts, sophisticated predictive molecular dynamics (MD) models for C-S-H have been developed, and the application of MD simulation has been expanded from fundamental understanding of C-S-H to nano-engineered cement composites. This paper comprehensively reviewed the current state of MD simulation on calcium-silicate-hydrate (C-S-H) and its diverse applications to nano-engineered cement composites, including carbon-based nanomaterials (i.e., carbon nanotube, graphene, graphene oxide), reinforced cement, cement–polymer nanocomposites (with an application on 3D printing concrete), and chemical additives for improving environmental resistance. In conclusion, the MD method could not only compute but also visualize the nanoscale behaviors of cement hydrates and other ingredients in the cement matrix; thus, fundamental properties of C-S-H structure and its interaction with nanoparticles can be well understood. As a result, the MD enabled us to identify and evaluate the performance of new advanced nano-engineered cement composites.

Contribution of Discharge Excited Atomic N, N2 *, and N2 + to a Plasma/Liquid Interfacial Reaction as Suggested by Quantitative Analysis

Electric-discharge nitrogen comprises three main types of excited nitrogen species-atomic nitrogen (N_atom ), excited nitrogen molecules (N₂ *), and nitrogen ions (N₂ ⁺ ) - which have different lifetimes and reactivities. In particular, the interfacial reaction locus between the discharged nitrogen and the water phase produces nitrogen compounds such as ammonia and nitrate ions (denoted as N-compounds generically); this is referred to as the plasma/liquid interfacial (P/L) reaction. The N_atom amount was analyzed quantitatively to clarify the contribution of N_atom to the P/L reaction. We focused on the quantitative relationship between N_atom and the produced N-compounds, and found that both N₂ * and N₂ ⁺ , which are active species other than N_atom , contributed to P/L reaction. The production of N-compounds from N₂ * and N₂ ⁺ was enhanced upon UV irradiation of the water phase, but the production of N-compounds from N_atom did not increase by UV irradiation. These results revealed that the P/L reactions starting from N_atom and those starting from N₂ ^* and N₂ ⁺ follow different mechanisms.

Pressure response of the THz spectrum of bulk liquid water revealed by intermolecular instantaneous normal mode analysis

The radial distribution functions of liquid water are known to change significantly their shape upon hydrostatic compression from ambient conditions deep into the kbar pressure regime. It has been shown that despite their eye-catching changes, the fundamental locally tetrahedral fourfold H-bonding pattern that characterizes ambient water is preserved up to about 10 kbar (1 GPa), which is the stability limit of liquid water at 300 K. The observed increase in coordination number comes from pushing water molecules into the first coordination sphere without establishing an H-bond, resulting in roughly two such additional interstitial molecules at 10 kbar. THz spectroscopy has been firmly established as a powerful experimental technique to analyze H-bonding in aqueous solutions given that it directly probes the far-infrared lineshape and thus the prominent H-bond network mode around 180 cm^-1. We, therefore, set out to assess pressure effects on the THz response of liquid water at 10 kbar in comparison to the 1 bar (0.1 MPa) reference, both at 300 K, with the aim to trace back the related lineshape changes to the structural level. To this end, we employ the instantaneous normal mode approximation to rigorously separate the H-bonding peak from the large background arising from the pronounced librational tail. By exactly decomposing the total molecular dynamics into hindered translations, hindered rotations, and intramolecular vibrations, we find that the H-bonding peak arises from translation-translation and translation-rotation correlations, which are successively decomposed down to the level of distinct local H-bond environments. Our utmost detailed analysis based on molecular pair classifications unveils that H-bonded double-donor water pairs contribute most to the THz response around 180 cm^-1, whereas interstitial waters are negligible. Moreover, short double-donor H-bonds have their peak maximum significantly shifted toward higher frequencies with respect to such long H-bonds. In conjunction with an increasing relative population of these short H-bonds versus the long ones (while the population of other water pair classes is essentially pressure insensitive), this explains not only the blue-shift of the H-bonding peak by about 20-30 cm^-1 in total from 1 bar to 10 kbar but also the filling of the shallow local minimum of the THz lineshape located in between the network peak and the red-wing of the librational band at 1 bar. Based on the changing populations as a function of pressure, we are also able to roughly estimate the pressure-dependence of the H-bond network mode and find that its pressure response and thus the blue-shifting are most pronounced at low kbar pressures.

On the complex structural diffusion of proton holes in nanoconfined alkaline solutions within slit pores

The hydroxide anion OH⁻(aq) in homogeneous bulk water, that is, the solvated proton hole, is known to feature peculiar properties compared with excess protons solvated therein. In this work, it is disclosed that nanoconfinement of such alkaline aqueous solutions strongly affects the key structural and dynamical properties of OH⁻(aq) compared with the bulk limit. The combined effect of the preferred hypercoordinated solvation pattern of OH⁻(aq), its preferred perpendicular orientation relative to the confining surfaces, the pronounced layering of nanoconfined water and the topology of the hydrogen bond network required for proton hole transfer lead to major changes of the charge transport mechanism, in such a way that the proton hole migration mechanism depends exquisitely on the width of the confined space that hosts the water film. Moreover, the anionic Zundel complex, which is of transient nature in homogeneous bulk solutions, can be dynamically trapped as a shallow intermediate species by suitable nanoconfinement conditions.

Confined liquids can exhibit very different properties compared with the bulk. Here, the authors show that the migration mechanism of the hydroxide anion in water is greatly affected by bi-dimensional nano-confinement and strongly depends on the width of the confined space.

Water bilayers on ZnO(101̄0) surfaces: data-driven structural search

We demonstrate an approach for the use of data science methods for structural search for high-stability atomic structures in ab initio simulation, via the analysis of a large set of candidate structures.

Neural network molecular dynamics simulations of solid–liquid interfaces: water at low-index copper surfaces

Molecular dynamics simulation of the water–copper interface have been carried out using high-dimensional neural network potential based on density functional theory.

Structure and dynamics of H2O vis-á-vis phenylalanine recognition at a DPPC lipid membrane via interfacial H-bond types: insights from polarized FT-IRRAS and ADMP simulations

Preferential and enantioselective interactions of L-/D-Phenylalanine (L-Phe and D-Phe) and butoxycarbonyl-protected L-/D-Phenylalanine (LPA and DPA) as guest with 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (L-DPPC) as host were tapped by using real time Fourier transform infrared reflection absorption spectroscopy (FT-IRRAS). Polarization-modulated FT-IRRAS of DPPC monolayers above the phenylalanine modified subphases depicted fine structure/conformation differences under considerations of controlled 2D surface pressure. Selective molecular recognition of D-enantiomer over L-enantiomer driven by the DPPC head group via H-bonding and electrostatic interactions was evident spectroscopically. Accordingly, binding constants (K) of 145, 346, 28, and 56 M(-1) for LPA, DPA, L-Phe, and D-Phe, respectively, were estimated. The real time FT-IRRAS water bands were strictly conformation sensitive. The effect of micro-solvation on the structure and stability of the 1:1 diastereomeric L-lipid⋯, LPA/DPA and L-lipid⋯, (L/D)-Phe adducts was investigated with the aid of Atom-centered Density Matrix Propagation (ADMP), a first principle quantum mechanical molecular dynamics approach. The phosphodiester fragment was the primary site of hydration where specific solvent interactions were simulated through single- and triple- "water-phosphate" interactions, as water cluster's "tetrahedral dice" to a "trimeric motif" transformation as a partial de-clusterization was evident. Under all the hydration patterns considered in both static and dynamic descriptions of density functional theory, L-lipid/D-amino acid enantiomer adducts continued to be stable structures while in dynamic systems, water rearranged without getting "squeezed-out" in the process of recognition. In spite of the challenging computational realm of this multiscale problem, the ADMP simulated molecular interactions complying with polarized vibrational spectroscopy unraveled a novel route to chiral recognition and interfacial water structure.

Changes in the hydrogen-bonding strength of internal water molecules and cysteine residues in the conductive state of channelrhodopsin-1

Water plays an essential role in the structure and function of proteins, particularly in the less understood class of membrane proteins. As the first of its kind, channelrhodopsin is a light-gated cation channel and paved the way for the new and vibrant field of optogenetics, where nerve cells are activated by light. Still, the molecular mechanism of channelrhodopsin is not understood. Here, we applied time-resolved FT-IR difference spectroscopy to channelrhodopsin-1 from Chlamydomonas augustae. It is shown that the (conductive) P2(380) intermediate decays with τ ≈ 40 ms and 200 ms after pulsed excitation. The vibrational changes between the closed and the conductive states were analyzed in the X-H stretching region (X = O, S, N), comprising vibrational changes of water molecules, sulfhydryl groups of cysteine side chains and changes of the amide A of the protein backbone. The O-H stretching vibrations of "dangling" water molecules were detected in two different states of the protein using H2 (18)O exchange. Uncoupling experiments with a 1:1 mixture of H2O:D2O provided the natural uncoupled frequencies of the four O-H (and O-D) stretches of these water molecules, each with a very weakly hydrogen-bonded O-H group (3639 and 3628 cm(-1)) and with the other O-H group medium (3440 cm(-1)) to moderately strongly (3300 cm(-1)) hydrogen-bonded. Changes in amide A and thiol vibrations report on global and local changes, respectively, associated with the formation of the conductive state. Future studies will aim at assigning the respective cysteine group(s) and at localizing the "dangling" water molecules within the protein, providing a better understanding of their functional relevance in CaChR1.

Regioselective long-range proton transfer in new rifamycin antibiotics: a process in which crown ethers act as stronger Brønsted bases than amines

Water-mediated proton transfer in six new derivatives of 3-formylrifamycin SV that contain crown, aza-crown, and benzo-crown ether rings were investigated by FTIR and NMR spectroscopy. (1)H-(1)H COSY couplings provide evidence for the formation of zwitterionic structures of the aza-crown and crown ether derivatives of rifamycin, in which a proton from one of the phenolic groups is transferred to tertiary and secondary nitrogen atoms. The increased intensity of the continuous absorption in the mid-infrared region together with the NMR data indicate proton transfer from the phenol group of the rifamycin core to the cavity of the benzo-crown ether ring. This proton transfer is achieved by formation of hydronium (H3O(+)) or Zundel ions (H5O2(+)), which form intermolecular hydrogen bonds with the oxygen atoms of the crown ether. DFT calculations are in agreement with the spectroscopic data and allow visualization of the structures of all new rifamycin derivatives, characterized by different intramolecular protonation sites.

Water structure and solvation of osmolytes at high hydrostatic pressure: pure water and TMAO solutions at 10 kbar versus 1 bar

Solvation structures of trimethylamine N-oxide change drastically due to the increase in the hydrostatic pressure.

Hydrogen bonding and molecular orientation at the liquid–vapour interface of water

We discover two hydrogen bonding motifs in liquid water and a related asymmetry linked with molecular orientation at its liquid–vapour interface.

A statistical mechanical theory of proton transport kinetics in hydrogen-bonded networks based on population correlation functions with applications to acids and bases

Extraction of relaxation times, lifetimes, and rates associated with the transport of topological charge defects in hydrogen-bonded networks from molecular dynamics simulations is a challenge because proton transfer reactions continually change the identity of the defect core. In this paper, we present a statistical mechanical theory that allows these quantities to be computed in an unbiased manner. The theory employs a set of suitably defined indicator or population functions for locating a defect structure and their associated correlation functions. These functions are then used to develop a chemical master equation framework from which the rates and lifetimes can be determined. Furthermore, we develop an integral equation formalism for connecting various types of population correlation functions and derive an iterative solution to the equation, which is given a graphical interpretation. The chemical master equation framework is applied to the problems of both hydronium and hydroxide transport in bulk water. For each case it is shown that the theory establishes direct links between the defect's dominant solvation structures, the kinetics of charge transfer, and the mechanism of structural diffusion. A detailed analysis is presented for aqueous hydroxide, examining both reorientational time scales and relaxation of the rotational anisotropy, which is correlated with recent experimental results for these quantities. Finally, for OH(-)(aq) it is demonstrated that the "dynamical hypercoordination mechanism" is consistent with available experimental data while other mechanistic proposals are shown to fail. As a means of going beyond the linear rate theory valid from short up to intermediate time scales, a fractional kinetic model is introduced in the Appendix in order to describe the nonexponential long-time behavior of time-correlation functions. Within the mathematical framework of fractional calculus the power law decay ∼t(-σ), where σ is a parameter of the model and depends on the dimensionality of the system, is obtained from Mittag-Leffler functions due to their long-time asymptotics, whereas (stretched) exponential behavior is found for short times.

Dissecting the THz spectrum of liquid water from first principles via correlations in time and space

Solvation of molecules in water is at the heart of a myriad of molecular phenomena and of crucial importance to understanding such diverse issues as chemical reactivity or biomolecular function. Complementing well-established approaches, it has been shown that laser spectroscopy in the THz frequency domain offers new insights into hydration from small solutes to proteins. Upon introducing spatially-resolved analyses of the absorption cross section by simulations, the sensitivity of THz spectroscopy is traced back to characteristic distance-dependent modulations of absorption intensities for bulk water. The prominent peak at approximately 200 cm(-1) is dominated by first-shell dynamics, whereas a concerted motion involving the second solvation shell contributes most significantly to the absorption at about 80 cm(-1) approximately 2.4 THz. The latter can be understood in terms of an umbrella-like motion of two hydrogen-bonded tetrahedra along the connecting hydrogen bond axis. Thus, a modification of the hydrogen bond network, e.g., due to the presence of a solute, is expected to affect vibrational motion and THz absorption intensity at least on a length scale that corresponds to two layers of solvating water molecules. This result provides a molecular mechanism explaining the experimentally determined sensitivity of absorption changes in the THz domain in terms of distinct, solute-induced dynamical properties in solvation shells of (bio)molecules--even in the absence of well-defined resonances.

Proton transfer 200 years after von Grotthuss: insights from ab initio simulations

In the last decade, ab initio simulations and especially Car-Parrinello molecular dynamics have significantly contributed to the improvement of our understanding of both the physical and chemical properties of water, ice, and hydrogen-bonded systems in general. At the heart of this family of in silico techniques lies the crucial idea of computing the many-body interactions by solving the electronic structure problem "on the fly" as the simulation proceeds, which circumvents the need for pre-parameterized potential models. In particular, the field of proton transfer in hydrogen-bonded networks greatly benefits from these technical advances. Here, several systems of seemingly quite different nature and of increasing complexity, such as Grotthuss diffusion in water, excited-state proton-transfer in solution, phase transitions in ice, and protonated water networks in the membrane protein bacteriorhodopsin, are discussed in the realms of a unifying viewpoint.

Coherence loss in light backscattering by random media with nanoscale nonuniformities

An experimental technique for measuring time-resolved coherence loss and destruction of backscattered wave packets in random media is described. The results of such measurements, performed with a modified Michelson interferometer, contain rich information about the characteristics of media nonuniformities. Experimental data for model nanosuspensions are compared with theoretical expressions developed in the paper which include the effects of Mie-type resonant scattering. We attribute one such observed effect to enhanced ineleastic optical transitions near the surface of nonmetallic nanoparticles. The inverse problem of characterization of multiscattering random media by backscattering is also considered.

The Liquid Water−Benzene System

The 500 MHz NMR spectra of water-benzene solution near saturation at 303.15, 323.15, and 343.15 K indicate that there is a proton-proton exchange between the water and benzene molecules. In the solution water appears to be present as a dimer attached to the benzene pi cloud on one side of each of the two (initially degenerate) fundamental energy levels, as predicted by the Jahn-Teller effect. This view is reinforced by the fact that one of its hydrogen atoms hovers above one of the carbon atoms and the other three are spread upward around the C6 axis of the benzene molecule. It is also supported by the calculated NMR spectra. Both effects are responsible for the change in the NMR spectra of the water molecules from a single line into four AB signals.

Water adsorption on ZnO(101̄0): from single molecules to partially dissociated monolayers

Static and dynamic density functional calculations have been used to study the structure and energetics of water adsorbed on the main cleavage plane of ZnO. In the single molecule limit we find that molecular adsorption is strongly preferred. The water binding energy increases for higher coverages due to an almost isotropic attractive water-water interaction which leads to clustering and formation of monolayer islands in the low water coverage regime. A thermodynamic analysis further shows that the full water monolayer is clearly the most stable phase until water starts to desorb. The water monolayer is even more stabilized by a partial dissociation of the water molecules, yielding as most stable configuration a (2x1) superstructure where every second water molecule is cleaved. The dissociation barrier for this process is very small which allows for an auto-dissociation of the water molecules even at low temperatures as observed experimentally. Finally we find that the energy cost involved to form [1210]-oriented domain boundaries between (2x1) patches with different orientation is almost negligible which explains the abundance of such domain boundaries in STM images.

Functional waters in intraprotein proton transfer monitored by FTIR difference spectroscopy

Much progress has been made in our understanding of water molecule reactions on surfaces, proton solvation in gas-phase water clusters and proton transfer through liquids. Compared with our advanced understanding of these physico-chemical systems, much less is known about individual water molecules and their cooperative behaviour in heterogeneous proteins during enzymatic reactions. Here we use time-resolved Fourier transform infrared spectroscopy (trFTIR) and in situ H2(18)O/H2(16)O exchange FTIR to determine how the membrane protein bacteriorhodopsin uses the interplay among strongly hydrogen-bonded water molecules, a water molecule with a dangling hydroxyl group and a protonated water cluster to transfer protons. The precise arrangement of water molecules in the protein matrix results in a controlled Grotthuss proton transfer, in contrast to the random proton migration that occurs in liquid water. Our findings support the emerging paradigm that intraprotein water molecules are as essential for biological functions as amino acids.

Dual-Scale Modeling of Benzene Adsorption onto Ni(111) and Au(111) Surfaces in Explicit Water

We present a multiscale modeling approach for studying interactions of organic molecules with metal surfaces in explicit water. The approach is based on combining adsorption energies of isolated molecules on transition metal surfaces calculated by ab initio density functional methods and classical molecular dynamics simulations using atomistically detailed force fields. The interaction of benzene with Ni(111) and Au(111) surfaces was studied. It is shown that a strong affinity of water for the hydrophilic surfaces makes benzene adsorption on Au thermodynamically unfavorable, while on Ni there is no preference. The work presented here serves as a first step in modeling the interactions of larger organic molecules with metal surfaces.

New information on water interfacial structure revealed by phase-sensitive surface spectroscopy

A phase-sensitive sum-frequency vibrational spectroscopic technique is developed to study interfacial water structure of water/quartz interfaces. Measurements allow deduction of both real and imaginary parts of the surface nonlinear spectral response, revealing an unprecedentedly detailed picture of the net polar orientations of the water species at the interface. The orientations of the icelike and liquidlike species appear to respond very differently to the bulk pH change indicating the existence of different surface sites on quartz with different deprotonation pK values.

How does water boil?

Insight into the boiling of water is obtained from molecular dynamics simulations. The process is initiated by the spontaneous formation of small vacuum cavities in liquid water. By themselves, these defects are very short lived. If, however, several cavities occur at close distances, they are likely to merge into larger vacuum holes. At the liquid-vapor interfaces, single or small groups of water molecules tend to leave the liquid surface. Once the system is propagated beyond the transition state, these evaporation events outnumber the competing reintegration into the hydrogen-bonded network.