Ab initio molecular dynamics simulation of the solvation and transport of hydronium and hydroxyl ions in water

Tailoring the Superatomic Characteristics and Optical Behavior of Metal-Free Boron Clusters via Ligand Engineering

It is of great importance to understand how the number and type of ligands influence the properties of clusters through ligand engineering, as this knowledge is crucial for the rational design and optimization of functional materials. Herein, the geometrical structures, binding energies, and electronic properties of nonmetallic Bn (n = 20 and 40) clusters with CO, PEt3, F, NO2, and CN ligands are systematically explored based on density functional theory (DFT) calculations. Our findings demonstrate that the CO ligand acts as an electron donor when attached to these two boron clusters, in contrast to their role as electron acceptors in interactions with metal oxide and metal chalcogenide clusters. This emphasizes the necessity of considering the intrinsic properties of the host cluster when modifying with ligands. Moreover, it was observed that substituting PEt3 with F, NO2, or CN converted the B20 cluster from an electron acceptor to an electron donor, thereby demonstrating the versatility in tuning the redox characteristics of boron clusters by selecting appropriate ligands. Intriguingly, the attachment of the PEt3, F, NO2, and CN ligands to B20 can significantly modulate the electronic properties of B20 to realize the formation of metal-free superalkali (B20(PEt3)n, n = 3-5) and superhalogen (B20F, B20NO2, and B20CN) clusters. Furthermore, the structure, stability, and optical absorption of the charge transfer complex B20(PEt3)3+B20F were analyzed. This complex has been identified as an efficient material for harvesting visible light. Our findings provide insights into the effects of ligand variations on boron cluster functionalities, offering a new perspective for the design of advanced materials with tailored cluster properties through ligand engineering.

Anisotropic Effects in Local Anodic Oxidation Nanolithography on Silicon Surfaces: Insights from ReaxFF Molecular Dynamics

Fully understanding the anisotropic effect of silicon surface orientations in local anodic oxidation (LAO) nanolithography processes is critical to the precise control of oxide quality and rate. This study used ReaxFF MD simulations to reveal the surface anisotropic effects in the LAO through the analysis of adsorbed species, atomic charge, and oxide growth. Our results show that the LAO behaves differently on silicon (100), (110), and (111) surfaces. Specifically, the application of an electric field significantly increases the quantity of surface-adsorbed -OH2 while reducing -OH on the (111) surface, and results in a higher charge on a greater number of Si atoms on the (100) surface. Moreover, the quantity of surface-adsorbed -OH plays a pivotal role in influencing the oxidation rate, as it directly correlates with an increased formation rate of Si-O-Si bonds. During bias-induced oxidation, the (111) surface appears with a high initial oxidation rate among three surfaces, while the (110) surface underwent increased oxidation at higher electric field strengths. This conclusion is based on the analysis of the evolution of Si-O-Si bond number, surface elevation, and oxide thickness. Our findings align well with prior theoretical and experimental studies, providing deeper insights and clear guidance for the fabrication of high-performance nanoinsulator gates using LAO nanolithography.

Characterization of the Coordination and Solvation Dynamics of Solvated Systems─Implications for the Analysis of Molecular Interactions in Solutions and Pure H2O

The characterization of solvation shells of atoms, ions, and molecules in solution is essential to relate solvation properties to chemical phenomena such as complex formation and reactivity. Different definitions of the first-shell coordination sphere from simulation data can lead to potentially conflicting data on the structural properties and associated ligand exchange dynamics. The definition of a solvation shell is typically based on a given threshold distance determined from the respective solute-solvent pair distribution function g(r) (i.e., GC). Alternatively, a nearest neighbor (NN) assignment based on geometric properties of the coordination complex without the need for a predetermined cutoff criterion, such as the relative angular distance (RAD) or the modified Voronoi (MV) tessellation, can be applied. In this study, the effect of different NN algorithms on the coordination number and ligand exchange dynamics evaluated for a series of monatomic ions in aqueous solution, carbon dioxide in aqueous and dichloromethane solutions, and pure liquid water has been investigated. In the case of the monatomic ions, the RAD approach is superior in achieving a well separated definition of the first solvation layer. In contrast, the MV algorithm provides a better separation of the NNs from a molecular point of view, leading to better results in the case of solvated CO2. When analyzing the coordination environment in pure water, the cutoff-based GC framework was found to be the most reliable approach. By comparison of the number of ligand exchange reactions and the associated mean ligand residence times (MRTs) with the properties of the coordination number autocorrelation functions, it is shown that although the average coordination numbers are sensitive to the different definitions of the first solvation shell, highly consistent estimates for the associated MRT of the solvated system are obtained in the majority of cases.

Transition metal single-atoms supported on hexagonal ZnIn2S4 monolayers for the hydrogen evolution reaction

Herein, we report a series of 5d transition metal (TM) single atoms supported on ZIS as promising catalysts for the hydrogen evolution reaction using first-principles calculations. The binding behaviors of TMs with the ZIS surface in single-atom catalyst formation are analysed using the adsorption energy (Eads), partial density of states (PDOS), charge density difference (CDD), and crystal orbital Hamilton population (COHP). The TM@ZIS (TM = Ta, W, Re, Os, Ir, and Pt) shows excellent hydrogen evolution performance with the Gibbs free energy (ΔGH*) values from -0.120 to 0.128 eV. The Tafel and Heyrovsky reaction mechanisms to drive H2 formation are also identified.

Nuclear Quantum Effects in Hydroxide Hydrate Along the H-Bond Bifurcation Pathway

Path integral (PI) simulations are used to explore nuclear quantum effects (NQEs) in hydroxide hydrate and its perdeuterated isotopomer along the H-bond bifurcation pathway. Toward this, a new potential energy surface using the symmetric gradient domain machine learning method with ab initio data at the CCSD(T)/aug-cc-pVTZ level is built. From PI umbrella sampling (US) simulations, free energy profiles along the bifurcation coordinate are explored as a function of temperature. At ambient temperature, the bifurcation barrier is increased upon inclusion of NQEs. At low temperatures in the deep tunneling regime, the barrier is strongly decreased and flattened. These trends are examined, and the role of the O-O distance is also investigated through two-dimensional US simulations.

Transition from vehicle to Grotthuss proton transfer in a nanosized flask: cryogenic ion spectroscopy of protonated p-aminobenzoic acid solvated with D2O

Proton transfer (PT) is one of the most ubiquitous reactions in chemistry and life science. The unique nature of PT has been rationalized not by the transport of a solvated proton (vehicle mechanism) but by the Grotthuss mechanism in which a proton is transported to the nearest proton acceptor along a hydrogen-bonded network. However, clear experimental evidence of the Grotthuss mechanism has not been reported yet. Herein we show by infrared spectroscopy that a vehicle-type PT occurs in the penta- and hexahydrated clusters of protonated p-aminobenzoic acid, while Grotthuss-type PT is observed in heptahydrated clusters, indicating a change in the PT mechanism depending on the degree of hydration. These findings emphasize the importance of the usually ignored vehicle mechanism as well as the degree of hydration. It highlights the possibility of controlling the PT mechanism by the number of water molecules in chemical and biological environments.

Effective boosting of halogenated α, β-unsaturated C4-dicarbonyl electrocatalytic hydrodehalogenation by 1 T'-MoS2/Ti3C2T2 (T = O, OH, F) heterojunctions: A theoretical study

Halogenated α, β-unsaturated C4-dicarbonyl (X-BDA), a novel family of high-toxicity ring cleavage products, is produced during the disinfection of phenolic compounds. The technique of electrocatalytic hydrodehalogenation (ECH) is efficient in rupturing carbon-halogen bonds and generating useful chemicals. This study used first principles to examine the ECH reaction mechanism of X-BDA and the subsequent hydrogenation reaction of the toxic derivative BDA over the 1 T'-MoS2/Ti3C2T2 (T = O, OH, F) catalysts. The catalytic activity of Ti3C2T2 (T = O, OH, F) catalysts decreases gradually with -OH, -F, -O functional group. The loading of 1 T'-MoS2 onto the Ti3C2T2 surface improves the stability and selectivity of Ti3C2T2. In particular, 1 T'-MoS2/Ti3C2(OH)2 is most conducive to the ECH reaction of X-BDA via a direct-indirect continuous reduction process. It exhibits excellent removal capability towards Cl-BDA, with decreasing reactivity in the order of the Cl-, Br-, and I-BDA. The material offers a solution to the challenging dechlorination issue. The dehalogenated product BDA can be hydrogenated to produce 1,4-butanedial, 1,4-butanediol, and 1,4-butenediol. Three valuable chemicals can be obtained by exerting an applied potential of - 0.65 V. This work suggests that the formation of heterojunction catalyst may lead to new strategies to improve ECH for environmental remediation applications.

Comparison of classical and ab initio simulations of hydronium and aqueous proton transfer

Proton transport in aqueous systems occurs by making and breaking covalent bonds, a process that classical force fields cannot reproduce. Various attempts have been made to remedy this deficiency, by valence bond theory or instantaneous proton transfers, but the ability of such methods to provide a realistic picture of this fundamental process has not been fully evaluated. Here we compare an ab initio molecular dynamics (AIMD) simulation of an excess proton in water to a simulation of a classical H3O+ in TIP3P water. The energy gap upon instantaneous proton transfer from H3O+ to an acceptor water molecule is much higher in the classical simulation than in the AIMD configurations evaluated with the same classical potential. The origins of this discrepancy are identified by comparing the solvent structures around the excess proton in the two systems. One major structural difference is in the tilt angle of the water molecules that accept an hydrogen bond from H3O+. The lack of lone pairs in TIP3P produces a tilt angle that is too large and generates an unfavorable geometry after instantaneous proton transfer. This problem can be alleviated by the use of TIP5P, which gives a tilt angle much closer to the AIMD result. Another important factor that raises the energy gap is the different optimal distance in water-water vs H3O+-water H-bonds. In AIMD the acceptor is gradually polarized and takes a hydronium-like configuration even before proton transfer actually happens. Ways to remedy some of these problems in classical simulations are discussed.

Understanding molecular and electrochemical charge transfer: theory and computations

Electron, proton, and proton-coupled electron transfer (PCET) are crucial elementary processes in chemistry, electrochemistry, and biology. We provide here a gentle overview of retrospective and currently developing theoretical formalisms of chemical, electrochemical and biological molecular charge transfer processes, with examples of how to bridge electron, proton, and PCET theory with experimental data. We offer first a theoretical minimum of molecular electron, proton, and PCET processes in homogeneous solution and at electrochemical interfaces. We illustrate next the use of the theory both for simple electron transfer processes, and for processes that involve molecular reorganization beyond the simplest harmonic approximation, with dissociative electron transfer and inclusion of all charge transfer parameters. A core example is the electrochemical reduction of the S2O82- anion. This is followed by discussion of core elements of proton and PCET processes and the electrochemical dihydrogen evolution reaction on different metal, semiconductor, and semimetal (say graphene) electrode surfaces. Other further focus is on stochastic chemical rate theory, and how this concept can rationalize highly non-traditional behaviour of charge transfer processes in mixed solvents. As a second major area we address ("long-range") chemical and electrochemical electron transfer through molecular frameworks using notions of superexchange and hopping. Single-molecule and single-entity electrochemistry are based on electrochemical scanning probe microscopies. (In operando) scanning tunnelling microscopy (STM) and atomic force microscopy (AFM) are particularly emphasized, with theoretical notions and new molecular electrochemical phenomena in the confined tunnelling gap. Single-molecule surface structure and electron transfer dynamics are illustrated by self-assembled thiol molecular monolayers and by more complex redox target molecules. This discussion also extends single-molecule electrochemistry to bioelectrochemistry of complex redox metalloproteins and metalloenzymes. Our third major area involves computational overviews of molecular and electronic structure of the electrochemical interface, with new computational challenges. These relate to solvent dynamics in bulk and confined space (say carbon nanostructures), electrocatalysis, metallic and semiconductor nanoparticles, d-band metals, carbon nanostructures, spin catalysis and "spintronics", and "hot" electrons. Further perspectives relate to metal-organic frameworks, chiral surfaces, and spintronics.

Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations

The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions.

Microscopic Mechanism of Proton Transfer in Pure Water under Ambient Conditions

Water molecules and the associated proton transfer (PT) are prevalent in chemical and biological systems and have been a hot research topic. Spectroscopic characterization and ab initio molecular dynamics (AIMD) simulations have previously revealed insights into acidic and basic liquids. Presumably, the situation in the acidic/basic solution is not necessarily the same as in pure water; in addition, the autoionization constant for water is only 10^-14 under ambient conditions, making the study of PT in pure water challenging. To overcome this issue, we modeled periodic water box systems containing 1000 molecules for tens of nanoseconds based on a neural network potential (NNP) with quantum mechanical accuracy. The NNP was generated by training a dataset containing the energies and atomic forces of 17 075 configurations of periodic water box systems, and these data points were calculated at the MP2 level that considers electron correlation effects. We found that the size of the system and the duration of the simulation have a significant impact on the convergence of the results. With these factors considered, our simulations showed that hydronium (H₃O⁺) and hydroxide (OH^-) ions in water have distinct hydration structures, thermodynamic and kinetic properties, e.g., the longer-lasting and more stable hydrated structure of OH^- ions than that of H₃O⁺, as well as a significantly higher free energy barrier for the OH^--associated PT than that of H₃O⁺, leading the two to exhibit completely different PT behaviors. Given these characteristics, we further found that PT via OH^- ions tends not to occur multiple times or between many molecules. In contrast, PT via H₃O⁺ can synergistically occur among multiple molecules and prefers to adopt a cyclic pattern among three water molecules, while it occurs mostly in a chain pattern when more water molecules are involved. Therefore, our studies provide a detailed and solid microscopic explanation for the PT process in pure water.

Ultrafast Proton Transfer Pathways Mediated by Amphoteric Imidazole

Imidazole, being an amphoteric molecule, can act both as an acid and as a base. This property enables imidazole, as an essential building block, to effectively facilitate proton transport in high-temperature proton exchange membrane fuel cells and in proton channel transmembrane proteins, enabling those systems to exhibit high energy conversion yields and optimal biological function. We explore the amphoteric properties of imidazole by following the proton transfer exchange reaction dynamics with the bifunctional photoacid 7-hydroxyquinoline (7HQ). We show with ultrafast ultraviolet-mid-infrared pump-probe spectroscopy how for imidazole, in contrast to expectations based on textbook knowledge of acid-base reactivity, the preferential reaction pathway is that of an initial proton transfer from 7HQ to imidazole, and only at a later stage a transfer from imidazole to 7HQ, completing the 7HQ tautomerization reaction. An assessment of the molecular distribution functions and first-principles calculations of proton transfer reaction barriers reveal the underlying reasons for our observations.

The effect of solvent on determining highest occupied molecular orbital energies of semiconducting organic molecules: Insight from a combined computational approach

Although cyclic voltammetry (CV) measurements in solution have been widely used to determine the highest occupied molecular orbital energy (E_HOMO ) of semiconducting organic molecules, an understanding of the experimentally observed discrepancies due to the solvent used is lacking. To explain these differences, we investigate the solvent effects on E_HOMO by combining density functional theory and molecular dynamics calculations for four donor molecules with a common backbone moiety. We compare the experimental E_HOMO values to the calculated values obtained from either implicit or first solvation shell theories. We find that the first solvation shell method can capture the E_HOMO variation arising from the functional groups in solution, unlike the implicit method. We further applied the first solvation shell method to other semiconducting organic molecules measured in solutions for different solvents. We find that the E_HOMO obtained using an implicit method is insensitive to solvent choice. The first solvation shell, however, produces E_HOMO values that are sensitive to solvent choices and agrees with published experimental results. The solvent sensitivity arises from a hierarchy of three effects: (1) the solute electronic state within a surrounding dielectric continuum, (2) ambient temperature or solvent atoms changing the solute geometry, and (3) electronic interactions between the solute and solvents. The implicit method, on the other hand, only captures the effect of a dielectric environment. Our findings suggest that E_HOMO obtained by CV measurements should account for the influence of solvent when the results are reported, interpreted, or compared to other molecules.

In situ electrochemical Raman spectroscopy and ab initio molecular dynamics study of interfacial water on a single-crystal surface

The dynamics and chemistry of interfacial water are essential components of electrocatalysis because the decomposition and formation of water molecules could dictate the protonation and deprotonation processes on the catalyst surface. However, it is notoriously difficult to probe interfacial water owing to its location between two condensed phases, as well as the presence of external bias potentials and electrochemically induced reaction intermediates. An atomically flat single-crystal surface could offer an attractive platform to resolve the internal structure of interfacial water if advanced characterization tools are developed. To this end, here we report a protocol based on the combination of in situ Raman spectroscopy and ab initio molecular dynamics (AIMD) simulations to unravel the directional molecular features of interfacial water. We present the procedures to prepare single-crystal electrodes, construct a Raman enhancement mode with shell-isolated nanoparticle, remove impurities, eliminate the perturbation from bulk water and dislodge the hydrogen bubbles during in situ electrochemical Raman experiments. The combination of the spectroscopic measurements with AIMD simulation results provides a roadmap to decipher the potential-dependent molecular orientation of water at the interface. We have prepared a detailed guideline for the application of combined in situ Raman and AIMD techniques; this procedure may take a few minutes to several days to generate results and is applicable to a variety of disciplines ranging from surface science to energy storage to biology.

Tracking the Delocalized Proton in Concerted Proton Transfer in Bulk Water

A solvated proton in water is often characterized as a charge or structural defect, and it is important to track its evolution on-the-fly in certain dynamics simulations. Previously, we introduced the proton indicator, a pseudo-atom, whose position approximates the location of the excess proton modeled as a structural defect. The proton indicator generally yields a smooth trajectory of a hydrated proton diffusing in aqueous solutions, including in the events of stepwise proton transfer via the Grotthuss mechanism; however, the proton indicator did not perform well in the events of concerted proton transfer, for which it occasionally yielded large position displacements between two successive time steps. To overcome this hurdle, we develop a new algorithm of a proton indicator with greatly enhanced performance for concerted proton transfer in bulk water. A protocol is proposed to exhaustively explore the hydrogen-bonding network of the water wires over which the excess proton is delocalized and to properly account for the contributions of the water molecules in this network as the geometry evolves. The new proton indicator (called Indicator 2.0) is assessed in dynamics simulations of an excess proton in bulk water and in specially constructed model systems of more complex architectures. The results demonstrate that the new indicator yields a smooth trajectory in both stepwise and concerted proton transfers.

Platinum-Water Interaction Induced Interfacial Water Orientation That Governs the pH-Dependent Hydrogen Oxidation Reaction

Understanding the electrode-water interface structure in acid and alkali is crucial to unveiling the underlying mechanism of pH-dependent hydrogen oxidation reaction (HOR) kinetics. In this work, we construct the explicit Pt(111)-H₂O interface models in both acid and alkali to investigate the relationship between the HOR mechanism and electrode-electrolyte interface structure using ab initio molecular dynamics and density functional theory. We find that the interfacial water orientation in the outer Helmholtz layer (OHP) induced by the Pt-water interaction governs the pH-dependent HOR kinetics on Pt(111). In alkali, the strong Pt-interfacial water electrostatic interaction behaves as a narrow OHP, which increases the proportion of "H-down" interfacial water and leads to less adsorbed water entering the inner Helmholtz plane (IHP), decreasing the work function of Pt(111). Furthermore, the more "H-down" interfacial water stabilizes the H_ad adsorption, prevents H_ad desorption, and suppresses the Volmer step of HOR by forming the solvated [H_ad···H₂O···H₂O] complex. Our work provided a visualized molecular-level mechanism to understand the nature of pH-dependent HOR kinetics.

Computational Insight into Defective Boron Nitride Supported Double-Atom Catalysts for Electrochemical Nitrogen Reduction

Designing highly selective and efficient double-atom electrocatalysts (DACs) is essential for achieving a superior nitrogen-reduction reaction (NRR) performance. Herein, we explored the defective boron nitride–supported cage-like double-atom catalysts to rummage the qualified NRR catalysts. Based on a systematic evaluation of the stability, N2 adsorption, NRR selectivity and activity of 10 DACs of TM1-TM2@VB-BN, we predicted Ru-Ti@VB-BN to be the NRR candidate with a limiting potential of −0.40 V. Compared to the corresponding single-atom catalysts, the introduction of Ti/Mo modulates the d-band center of the active metal atom, which improves the NRR performance. Moreover, the magnetic Ru-Ti dimer can facilitate the transfer of charge to molecular N2, ensuring a significant activation of the inert N≡N bond. This research not only opens up new avenues for designing boron nitride–supported DACs for NRR, but also deepens the understanding of DACs in N2 activation.

Proton-transfer spectroscopy beyond the normal-mode scenario

A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single [Formula: see text] cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.

Hydrogen Bond Thermodynamics in Aqueous Acid Solutions: A Combined DFT and Classical Force-Field Approach

The thermodynamics of hydrogen bonds in aqueous and acidic solutions significantly impacts the kinetics and thermodynamics of acid reaction chemistry. We utilize in this work a multiscale approach, combining density functional theory (DFT) with classical molecular dynamics (MD) to model hydrogen bond thermodynamics in an acidic solution. Using thermodynamic cycles, we split the solution phase free energy into its gas phase counterpart plus solvation free energies. We validate this DFT/MD approach by calculating the aqueous phase hydrogen bond free energy between two water molecules (H₂O-···-H₂O), the free energy to transform an H₃O⁺ cation into an H₅O₂⁺ cation, and the hydrogen bond free energy of protonated water clusters (H₃O⁺-···-H₂O and H₅O₂⁺-···-H₂O). The computed equilibrium hydrogen bond free energy of H₂O-···-H₂O is remarkably accurate, especially considering the large individual contributions to the thermodynamic cycle. Turning to cations, we find the ion to be more stable than H₃O⁺ by roughly 1-2 k_BT. This small free energy difference allows for thermal fluctuation between the two idealized motifs, consistent with spectroscopic and simulation studies. Lastly, hydrogen bonding free energies between either H⁺ cation and H₂O in solution were found to be stronger than between two H₂O, though much less so than in vacuum because of dielectric screening in solution. Altogether, our results suggest the DFT/MD approach is promising for application in modeling hydrogen bonding and proton transfer thermodynamics in condensed phases.

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200-300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm-1, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm-1 and agrees in location and width with previous experimental mid-infrared spectra.

Inherent Promotion of Ionic Conductivity via Collective Vibrational Strong Coupling of Water with the Vacuum Electromagnetic Field

Hydrogen bonding interactions among water molecules play a critical role in chemical reactivity, dynamic proton mobility, static dielectric behavior, and the thermodynamic properties of water. In this study, we demonstrate the modification of ionic conductivity of water through hybridization with a vacuum electromagnetic field by strongly coupling the O─H stretching mode of H₂O to a Fabry-Perot cavity mode. The hybridization generates collective vibro-polaritonic states, thereby enhancing the proton conductivity by an order of magnitude at resonance. In addition, the dielectric constants increase at resonance in the coupled state. The findings presented herein reveal how a vacuum electromagnetic environment can be engineered to control the ground-state properties of water.

Molecular Mechanism of Autodissociation in Liquid Water: Ab Initio Molecular Dynamics Simulations

Autodissociation in liquid water is one of the most important processes in various topics of physical chemistry, such as acid-base chemistry. Molecular simulations have elucidated most of the molecular mechanisms at the atomic level, yet quantitative analysis to compare with experiments using the potential of mean force (PMF) remains a hurdle, including the definition of reaction coordinates and the accuracy of liquid structures by ab initio molecular dynamics (AIMD) simulations with density functional theory (DFT) methods. Here, we perform AIMD simulations with the revPBE-D3 exchange-correlation functional to compute the PMF profiles of autoionization, or proton transfer (PT), in liquid water. For the quantitative analysis with physically meaningful reaction coordinates, we employ a PT coordinate, donor-acceptor (OH^--H₃O⁺) distance, and hydrogen (H)-bond number. The one-dimensional (1D) PMF profile along the PT coordinate shows no local minimum in the product state of PT (OH^- and H₃O⁺), which is necessary to accurately compute the acid dissociation constant (or pK_a). On the other hand, the 2D PMF profiles along the PT coordinate and donor-acceptor distance show local minima in the product state and reaction barriers, and the computed pK_w is comparable to the experiment. In addition, the 2D PMF profiles along the PT coordinate and the H-bond number reveal the molecular mechanism of the H-bond rearrangement concomitant with PT, in which the H-bond breaking before PT is slightly preferable. These findings indicate that an accurate evaluation of pK_a by MD simulations requires the donor-acceptor distance in addition to the conventional PT coordinate.

Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling

Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water has been a challenging task from ab initio molecular dynamics (AIMD). The difficulty mainly comes from a lack of accurate and efficient exchange-correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly-constrained and appropriately normed (SCAN) meta-GGA functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange-correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice in describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. Our AIMD simulations indicate hydroxide diffuses slower than hydronium in water, which is consistent with the experiments.

Grotthuss Molecular Dynamics Simulations for Modeling Proton Hopping in Electrosprayed Water Droplets

Excess protons in water exhibit unique transport properties because they can rapidly hop along H-bonded water wires. Considerable progress has been made in unraveling this Grotthuss diffusion mechanism using quantum mechanical-based computational techniques. Unfortunately, high computational cost tends to restrict those techniques to small systems and short times. Molecular dynamics (MD) simulations can be applied to much larger systems and longer time windows. However, standard MD methods do not permit the dissociation/formation of covalent bonds, such that Grotthuss diffusion cannot be captured. Here, we bridge this gap by combining atomistic MD simulations (using Gromacs and TIP4P/2005 water) with proton hopping. Excess protons are modeled as hydronium ions that undergo H₃O⁺ + H₂O → H₂O + H₃O⁺ transitions. In accordance with ab initio MD data, these Grotthuss hopping events are executed in "bursts" with quasi-instantaneous hopping across one or more waters. The bursts are separated by regular MD periods during which H₃O⁺ ions undergo Brownian diffusion. The resulting proton diffusion coefficient agrees with the literature value. We apply this Grotthuss MD technique to highly charged water droplets that are in a size regime encountered during electrospray ionization (5 nm radius, ∼17,000 H₂O). The droplets undergo rapid solvent evaporation and occasional H₃O⁺ ejection, keeping them at ca. 81% of the Rayleigh limit. The simulated behavior is consistent with phase Doppler anemometry data. The Grotthuss MD technique developed here should be useful for modeling the behavior of various proton-containing systems that are too large for high-level computational approaches. In particular, we envision future applications related to electrospray processes, where earlier simulations used metal cations while in reality excess protons dominate.

Revealing the Anticorrelation Behavior Mechanism between the Grotthuss and Vehicular Diffusions for Proton Transport in Concentrated Acid Solutions

In this study, we performed reactive molecular dynamics simulations to characterize proton solvation and transport in concentrated hydrochloric acid solutions. The correlation contribution to the total proton mean square displacement is found to be negative for all acid concentrations, indicating the anticorrelation between the Grotthuss and vehicular diffusions. For the vehicular diffusion, the hydronium ions tend to move freely toward the lone pair side independent of acid concentrations, whereas for the Grotthuss diffusion, the proton hopping direction is limited to one of the hydrogen-bonded water molecules on the opposite side of the lone pair region, which are specifically oriented with respect to the neighboring hydronium ion at higher acid concentrations. This result is justified by our findings of the higher fraction of proton rattling with the single hopping event and longer hydrogen bond lifetimes at higher acid concentrations. However, the angular distribution for both the vehicular and Grotthuss diffusions is found to be rather broad and comparable for all acid concentrations, and thus, the anticorrelation shows a minimal dependence on the acid concentration. Our results reveal that the anticorrelation behavior between the vehicle and Grotthuss diffusions is attributed to the amphiphilic nature of hydronium ions and thus is independent of the acid concentrations in solutions.

Modulating the growth of chemically deposited ZnO nanowires and the formation of nitrogen- and hydrogen-related defects using pH adjustment

ZnO nanowires (NWs) grown by chemical bath deposition (CBD) have received great interest for nanoscale engineering devices, but their formation in aqueous solution containing many impurities needs to be carefully addressed. In particular, the pH of the CBD solution and its effect on the formation mechanisms of ZnO NWs and of nitrogen- and hydrogen-related defects in their center are still unexplored. By adjusting its value in a low- and high-pH region, we show the latent evolution of the morphological and optical properties of ZnO NWs, as well as the modulated incorporation of nitrogen- and hydrogen-related defects in their center using Raman and cathodoluminescence spectroscopy. The increase in pH is related to the increase in the oxygen chemical potential (μ _O), for which the formation energy of hydrogen in bond-centered sites (H_BC) and V_Zn-N_O-H defect complexes is found to be unchanged, whereas the formation energy of zinc vacancy (V_Zn) and zinc vacancy-hydrogen (V_Zn-nH) complexes steadily decreases as shown from density-functional theory calculations. Revealing that these V_Zn-related defects are energetically favorable to form as μ _O is increased, ZnO NWs grown in the high-pH region are found to exhibit a higher density of V_Zn-nH defect complexes than ZnO NWs grown in the low-pH region. Annealing at 450 °C under an oxygen atmosphere helps tuning the optical properties of ZnO NWs by reducing the density of H_BC and V_Zn-related defects, while activating the formation of V_Zn-N_O-H defect complexes. These findings show the influence of pH on the nature of Zn(ii) species, the electrostatic interactions between these species and ZnO NW surfaces, and the formation energy of the involved defects. They emphasize the crucial role of the pH of the CBD solution and open new possibilities for simultaneously engineering the morphology of ZnO NWs and the formation of nitrogen- and hydrogen-related defects.

Controlling Hydronium Diffusivity in Model Proton Exchange Membranes

Fuel-cell-based proton exchange membranes (PEMs) show great potential as cost-effective and clean energy conversion devices. In our recent work, we found that for the low-hydrated model PEMs with a inhomogeneous water distribution and a sulfonate anionic functional end group (SO₃^-), the H₃O⁺ reacts with SO₃^- according to SO₃^- + H₃O⁺ ↔ SO₃H + H₂O, indicating that the anions in PEMs become active participants in the hydronium diffusion. In this work, we use fully atomistic ab initio molecular dynamics simulations to elucidate the optimal conditions that would promote the participation of SO₃^- in the hydronium diffusion mechanism by increasing the H₃O⁺/SO₃^- reactivity, thus increasing the hydronium diffusivity along the cell. The results presented in this work allow us to suggest a set of design rules for creating novel, highly conductive PEMs operating at high temperatures under a nonuniform water distribution using a linker/anion with a relatively high pK_a such as (CH₂)₂SO₃. We expect that the discovery of these key design principles will play an important role in the synthesis of high-performing materials for emerging PEM-based fuel cell technologies.

The Ionic Product of Water in the Eye of the Quantum Cluster Equilibrium

The theoretical description of water properties continues to be a challenge. Using quantum cluster equilibrium (QCE) theory, we combine state-of-the-art quantum chemistry and statistical thermodynamic methods with the almost historical Clausius-Clapeyron relation to study water self-dissociation and the thermodynamics of vaporization. We pay particular attention to the treatment of internal rotations and their impact on the investigated properties by employing the modified rigid-rotor-harmonic-oscillator (mRRHO) approach. We also study a novel QCE parameter-optimization procedure. Both the ionic product and the vaporization enthalpy yield an astonishing agreement with experimental reference data. A significant influence of the mRRHO approach is observed for cluster populations and, consequently, for the ionic product. Thermodynamic properties are less affected by the treatment of these low-frequency modes.

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.

Photoelectron spectra of water and simple aqueous solutions at extreme conditions

Determining the electronic structure of aqueous solutions at extreme conditions is an important step towards understanding chemical bonding and reactions in water under pressure (P) and at high temperature (T)....

Molecular dynamics simulations of proton conducting mediums containing phosphoric acid

An anhydrous proton conducting electrolyte is a key material to realise medium-temperature fuel cells that can drastically simplify heat radiation systems in transportation applications. However, practical applications are limited by...

Resolving the Structural Debate for the Hydrated Excess Proton in Water

It has long been proposed that the hydrated excess proton in water (aka the solvated "hydronium" cation) likely has two limiting forms, that of the Eigen cation (H₉O₄⁺) and that of the Zundel cation (H₅O₂⁺). There has been debate over which of these two is the more dominant species and/or whether intermediate (or "distorted") structures between these two limits are the more realistic representation. Spectroscopy experiments have recently provided further results regarding the excess proton. These experiments show that the hydrated proton has an anisotropy reorientation time scale on the order of 1-2 ps. This time scale has been suggested to possibly contradict the picture of the more rapid "special pair dance" phenomenon for the hydrated excess proton, which is a signature of a distorted Eigen cation. The special pair dance was predicted from prior computational studies in which the hydrated central core hydronium structure continually switches (O-H···O)* special pair hydrogen-bond partners with the closest three water molecules, yielding on average a distorted Eigen cation with three equivalent and dynamically exchanging distortions. Through state-of-art simulations it is shown here that anisotropy reorientation time scales of the same magnitude are obtained that also include structural reorientations associated with the special pair dance, leading to a reinterpretation of the experimental results. These results and additional analyses point to a distorted and dynamic Eigen cation as the most prevalent hydrated proton species in aqueous acid solutions of dilute to moderate concentration, as opposed to a stabilized or a distorted (but not "dancing") Zundel cation.

Characterizing and Contrasting Structural Proton Transport Mechanisms in Azole Hydrogen Bond Networks Using Ab Initio Molecular Dynamics

Imidazole and 1,2,3-triazole are promising hydrogen-bonded heterocycles that conduct protons via a structural mechanism and whose derivatives are present in systems ranging from biological proton channels to proton exchange membrane fuel cells. Here, we leverage multiple time-stepping to perform ab initio molecular dynamics of imidazole and 1,2,3-triazole at the nanosecond time scale. We show that despite the close structural similarities of these compounds, their proton diffusion constants vary by over an order of magnitude. Our simulations reveal the reasons for these differences in diffusion constants, which range from the degree of hydrogen-bonded chain linearity to the effect of the central nitrogen atom in 1,2,3-triazole on proton transport. In particular, we uncover evidence of two "blocking" mechanisms in 1,2,3-triazole, where covalent and hydrogen bonds formed by the central nitrogen atom limit the mobility of protons. Our simulations thus provide insights into the origins of the experimentally observed 10-fold difference in proton conductivity.

p-Aminobenzoic acid protonation dynamics in an evaporating droplet by ab initio molecular dynamics

Protonation equilibria are known to vary from the bulk to microdroplet conditions, which could induce many chemical and physical phenomena. Protonated p-aminobenzoic acid (PABA + H⁺) can be considered a model for probing the protonation dynamics in an evaporating droplet, as its protonation equilibrium is highly dependent on the formation conditions from solution via atmospheric pressure ionization sources. Experiments using diverse experimental techniques have shown that protic solvents allow formation of the O-protomer (PABA protonated in the carboxylic acid group) stable in the gas phase, while aprotic solvents yield the N-protomer (protonated in the amino group) that is the most stable protomer in solution. In this work, we explore the protonation equilibrium of PABA solvated by different numbers of water molecules (n = 0 to 32) using ab initio molecular dynamics. For n = 8-32, the protonation is either at the NH₂ group or in the solvent network. The solvent network interacts with the carboxylic acid group, but there is no complete proton transfer to form the O-protomer. For smaller clusters, however, solvent-mediated proton transfers to the carboxylic acid were observed, both via the Grotthuss mechanism and the vehicle or shuttle mechanism (for n = 1 and 2). Thermodynamic considerations allowed a description of the origins of the kinetic trapping effect, which explains the observation of the solution structure in the gas phase. This effect likely occurs in the final evaporation steps, which are outside the droplet size range covered by previous classical molecular dynamics simulations of charged droplets. These results may be considered relevant in determining the nature of the species observed in the ubiquitous ESI based mass spectrometry analysis, and in general for droplet chemistry, explaining how protonation equilibria are drastically changed from bulk to microdroplet conditions.

High-throughput screening of single metal atom anchored on N-doped boron phosphide for N2 reduction

Developing eco-friendly and highly-efficient catalysts for the electrochemical nitrogen reduction reaction (NRR) under ambient conditions to replace the energy-intensive and environment-polluting Haber-Bosch process is of great significance, while remaining a long-standing challenge in the field of energy conversion today. Herein, through the first principles high-throughput screening, we systematically investigated the catalytic activity of a series of single metal atom immobilized on N-doped boron phosphide (N₃-BP) for N₂ reduction, denoted as MN₃-BP. In particular, a "four-step" screening strategy, involving the structural stability, N₂ chemisorption, low energy cost, as well as good selectivity, was adopted for the stringent screening of the promising MN₃-BP candidates for NRR. Our results unveil that among these candidates, MoN₃-BP eventually stands out, benefiting from its high selectivity and activity, as well as accompanying a considerably favorable limiting potential of -0.25 V for NRR. More impressively, the NRR activity origin of various candidates was revealed by the descriptor φ and ICOHP. Overall, our work not only accelerates the discovery of SACs for converting N₂ into sustainable NH₃ but also provides an exciting impetus for the rational design of NRR catalysts with high stability, high activity, and high selectivity.

Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism

The properties of the water network in concentrated HCl acid pools in nanometer-sized reverse non-ionic micelles were probed with TeraHertz absorption, dielectric relaxation spectroscopy, and reactive force field simulations capable of describing proton hopping mechanisms. We identify that only at a critical micelle size of W0=9 do solvated proton complexes form in the water pool, accompanied by a change in mechanism from Grotthuss forward shuttling to one that favors local oscillatory hopping. This is due to a preference for H+ and Cl- ions to adsorb to the micelle interface, together with an acid concentration effect that causes a "traffic jam" in which the short-circuiting of the hydrogen-bonding motif of the hydronium ion decreases the forward hopping rate throughout the water interior even as the micelle size increases. These findings have implications for atmospheric chemistry, biochemical and biophysical environments, and energy materials, as transport of protons vital to these processes can be suppressed due to confinement, aggregation, and/or concentration.

Electron binding energies and Dyson orbitals of OnH2n+1 +,0,- clusters: Double Rydberg anions, Rydberg radicals, and micro-solvated hydronium cations

Ab initio electron propagator methods are employed to predict the vertical electron attachment energies (VEAEs) of OH₃ ⁺(H₂O)_n clusters. The VEAEs decrease with increasing n, and the corresponding Dyson orbitals are diffused over exterior, non-hydrogen bonded protons. Clusters formed from OH₃ ^- double Rydberg anions (DRAs) and stabilized by hydrogen bonding or electrostatic interactions between ions and polar molecules are studied through calculations on OH₃ ^-(H₂O)_n complexes and are compared with more stable H^-(H₂O)_n+1 isomers. Remarkable changes in the geometry of the anionic hydronium-water clusters with respect to their cationic counterparts occur. Rydberg electrons in the uncharged and anionic clusters are held near the exterior protons of the water network. For all values of n, the anion-water complex H^-(H₂O)_n+1 is always the most stable, with large vertical electron detachment energies (VEDEs). OH₃ ^-(H₂O)_n DRA isomers have well separated VEDEs and may be visible in anion photoelectron spectra. Corresponding Dyson orbitals occupy regions beyond the peripheral O-H bonds and differ significantly from those obtained for the VEAEs of the cations.

The hopping mechanism of the hydrated excess proton and its contribution to proton diffusion in water

In this work, a series of analyses are performed on ab initio molecular dynamics simulations of a hydrated excess proton in water to quantify the relative occurrence of concerted hopping events and "rattling" events and thus to further elucidate the hopping mechanism of proton transport in water. Contrary to results reported in certain earlier papers, the new analysis finds that concerted hopping events do occur in all simulations but that the majority of events are the product of proton rattling, where the excess proton will rattle between two or more waters. The results are consistent with the proposed "special-pair dance" model of the hydrated excess proton wherein the acceptor water molecule for the proton transfer will quickly change (resonate between three equivalent special pairs) until a decisive proton hop occurs. To remove the misleading effect of simple rattling, a filter was applied to the trajectory such that hopping events that were followed by back hops to the original water are not counted. A steep reduction in the number of multiple hopping events is found when the filter is applied, suggesting that many multiple hopping events that occur in the unfiltered trajectory are largely the product of rattling, contrary to prior suggestions. Comparing the continuous correlation function of the filtered and unfiltered trajectories, we find agreement with experimental values for the proton hopping time and Eigen-Zundel interconversion time, respectively.

OH- and H3O+ Diffusion in Model AEMs and PEMs at Low Hydration: Insights from Ab Initio Molecular Dynamics

Fuel cell-based anion-exchange membranes (AEMs) and proton exchange membranes (PEMs) are considered to have great potential as cost-effective, clean energy conversion devices. However, a fundamental atomistic understanding of the hydroxide and hydronium diffusion mechanisms in the AEM and PEM environment is an ongoing challenge. In this work, we aim to identify the fundamental atomistic steps governing hydroxide and hydronium transport phenomena. The motivation of this work lies in the fact that elucidating the key design differences between the hydroxide and hydronium diffusion mechanisms will play an important role in the discovery and determination of key design principles for the synthesis of new membrane materials with high ion conductivity for use in emerging fuel cell technologies. To this end, ab initio molecular dynamics simulations are presented to explore hydroxide and hydronium ion solvation complexes and diffusion mechanisms in the model AEM and PEM systems at low hydration in confined environments. We find that hydroxide diffusion in AEMs is mostly vehicular, while hydronium diffusion in model PEMs is structural. Furthermore, we find that the region between each pair of cations in AEMs creates a bottleneck for hydroxide diffusion, leading to a suppression of diffusivity, while the anions in PEMs become active participants in the hydronium diffusion, suggesting that the presence of the anions in model PEMs could potentially promote hydronium diffusion.

Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane

To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. The results reported here are in reasonable agreement with experiment.

Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy

The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H⁺(aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H⁺⋅(H₂O)_n cluster ions isolated in the gas phase shed light on the chemical nature of H⁺(aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared-infrared (IR-IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H⁺(aq).

Proton transfer in bulk water using the full adaptive QM/MM method: integration of solute- and solvent-adaptive approaches

The quantum mechanical/molecular mechanical (QM/MM) method is a hybrid molecular simulation technique that increases the accessibility of local electronic structures of large systems. The technique combines the benefit of accuracy found in the QM method and that of cost efficiency found in the MM method. However, it is difficult to directly apply the QM/MM method to the dynamics of solution systems, particularly for proton transfer. As explained in the Grotthuss mechanism, proton transfer is a structural interconversion between hydronium ions and solvent water molecules. Hence, when the QM/MM method is applied, an adaptive treatment, namely on-the-fly revisions on molecular definitions, is required for both the solute and solvent. Although several solvent-adaptive methods have been proposed, a full adaptive framework, which is an approach that also considers adaptation for solutes, remains untapped. In this paper, we propose a new numerical expression for the coordinates of the excess proton and its control algorithm. Furthermore, we confirm that this method can stably and accurately simulate proton transfer dynamics in bulk water.

Water orientation and dynamics in the closed and open influenza B virus M2 proton channels

The influenza B M2 protein forms a water-filled tetrameric channel to conduct protons across the lipid membrane. To understand how channel water mediates proton transport, we have investigated the water orientation and dynamics using solid-state NMR spectroscopy and molecular dynamics (MD) simulations. 13C-detected water 1H NMR relaxation times indicate that water has faster rotational motion in the low-pH open channel than in the high-pH closed channel. Despite this faster dynamics, the open-channel water shows higher orientational order, as manifested by larger motionally-averaged 1H chemical shift anisotropies. MD simulations indicate that this order is induced by the cationic proton-selective histidine at low pH. Furthermore, the water network has fewer hydrogen-bonding bottlenecks in the open state than in the closed state. Thus, faster dynamics and higher orientational order of water molecules in the open channel establish the water network structure that is necessary for proton hopping.

Probing the Local Reaction Environment During High Turnover Carbon Dioxide Reduction with Ag-Based Gas Diffusion Electrodes

Discerning the influence of electrochemical reactions on the electrode microenvironment is an unavoidable topic for electrochemical reactions that involve the production of OH⁻ and the consumption of water. That is particularly true for the carbon dioxide reduction reaction (CO₂RR), which together with the competing hydrogen evolution reaction (HER) exert changes in the local OH⁻ and H₂O activity that in turn can possibly affect activity, stability, and selectivity of the CO₂RR. We determine the local OH⁻ and H₂O activity in close proximity to a CO₂‐converting Ag‐based gas diffusion electrode (GDE) with product analysis using gas chromatography. A Pt nanosensor is positioned in the vicinity of the working GDE using shear‐force‐based scanning electrochemical microscopy (SECM) approach curves, which allows monitoring changes invoked by reactions proceeding within an otherwise inaccessible porous GDE by potentiodynamic measurements at the Pt‐tip nanosensor. We show that high turnover HER/CO₂RR at a GDE lead to modulations of the alkalinity of the local electrolyte, that resemble a 16 m KOH solution, variations that are in turn linked to the reaction selectivity.