Calculation of Electrochemical Energy Levels in Water Using the Random Phase Approximation and a Double Hybrid Functional

Emerging Atomistic Modeling Methods for Heterogeneous Electrocatalysis

Heterogeneous electrocatalysis lies at the center of various technologies that could help enable a sustainable future. However, its complexity makes it challenging to accurately and efficiently model at an atomic level. Here, we review emerging atomistic methods to simulate the electrocatalytic interface with special attention devoted to the components/effects that have been challenging to model, such as solvation, electrolyte ions, electrode potential, reaction kinetics, and pH. Additionally, we review relevant computational spectroscopy methods. Then, we showcase several examples of applying these methods to understand and design catalysts relevant to green hydrogen. We also offer experimental views on how to bridge the gap between theory and experiments. Finally, we provide some perspectives on opportunities to advance the field.

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.

Unraveling the origin of reductive stability of super-concentrated electrolytes from first principles and unsupervised machine learning

Developing electrolytes with excellent electrochemical stability is critical for next-generation rechargeable batteries. Super-concentrated electrolytes (SCEs) have attracted great interest due to their high electrochemical performances and stability. Previous studies have revealed changes in solvation structures and shifts in lowest unoccupied molecular orbitals from solvents to anions, promoting the formation of an anion-derived solid-electrolyte-interphase (SEI) in SCE. However, a direct connection at the atomic level to electrochemical properties is still missing, hindering the rational optimization of electrolytes. Herein, we combine ab initio molecular dynamics with the free energy calculation method to compute redox potentials of propylene carbonate electrolytes at a range of LiTFSI concentrations, and moreover employ an unsupervised machine learning model with a local structure descriptor to establish the structure-property relations. Our calculation indicates that the network of TFSI- in SCE not only helps stabilize the added electron and renders the anion more prone to reductive decomposition, but also impedes the solvation of F- and favors LiF precipitation, together leading to effective formation of protective SEI layers. Our work provides new insights into the solvation structures and electrochemistry of concentrated electrolytes which are essential to electrolyte design in batteries.

Automated Workflow for Computation of Redox Potentials, Acidity Constants and Solvation Free Energies Accelerated by Machine Learning

Fast evolution of modern society stimulates intense development of new materials with novel functionalities in energy and environmental applications. Due to rapid progress of computer science, computational design of materials with target properties has recently attracted lots of interests. Accurate and efficient calculation of fundamental thermodynamic properties, including redox potentials, acidity constants, and solvation free energies, is of great importance for selection and design of desirable materials. Free energy calculation based on ab initio molecular dynamics (AIMD) can predict these properties with high accuracy at complex environments, however being impeded by high computational costs. To address this issue, this work develops an automated scheme that combines iterative training of machine learning potentials (MLPs) and free energy calculation, and deomonstrates that these thermodynamic properties can be computed by ML accelerated MD with ab initio accuracy and much longer time scale at cheaper costs, improving poor statistics and convergence of numerical integration by AIMD. Our automated scheme lays the foundation for computational chemistry-assisted materials design.

The energetics of electron and proton transfer to CO2 in aqueous solution

The electrocatalytic reduction of CO₂ is considered an effective method to reduce CO₂ emissions and achieve electrical/chemical energy conversion. It is crucial to determine the reaction mechanism so that the key reaction intermediates can be targeted and the overpotential lowered. The process involves the interaction with the electrode surface and with species, including the solvent, at the electrode-electrolyte interface, and it is therefore not easy to separate catalytic contributions of the electrode from those of the electrolyte. We have used density functional theory-based molecular dynamics to calculate the Gibbs free energy of the proton and electron transfer reactions corresponding to each step in the electroreduction of CO₂ to HCOOH in aqueous media. The results show thermodynamic pathways consistent with the mechanism proposed by Hori. Since electrodes are not included in this work, differences between the calculated results and the experimental observations can help determine the catalytic contribution of the electrode surface.

Doubly Hybrid Functionals Close to Chemical Accuracy for Both Finite and Extended Systems: Implementation and Test of XYG3 and XYGJ-OS

While being widely used to understand the chemical reactions in heterogeneous catalysis or other multidisciplinary systems, a great challenge that semilocal and hybrid density functional approximations (DFAs) are facing is to deliver a uniformly accurate description for both finite and extended systems. Herein, we perform reliable and well-converged periodic calculations of two doubly hybrid approximations (DHAs), XYG3 and XYGJ-OS, and demonstrate that the good accuracy of DHAs achieved for molecules is transferable to the semiconductors and insulators. Such an accuracy is not only for energetic properties but also for the first- and second-order response properties, which is general for different kinds of chemical environments, including simple cubic bulks, perovskite-type transition metal oxides like TiO2, and heterogeneous systems like CO adsorption on the NaCl(100) surface. The present finding has strengthened the predictive power of DFT, which not only will inspire the future development of the top-rung DFAs but also will boost their applications in multidisciplinary studies with high accuracy and efficiency.

Tuning the Performance of Aqueous Organic Redox Flow Batteries via First-Principles Calculations

Aqueous organic redox flow batteries have many appealing properties in the application of large-scale energy storage. The large chemical tunability of organic electrolytes shows great potential to improve the performance of flow batteries. Computational studies at the quantum-mechanics level are very useful for guiding experiments, but in previous studies, explicit water interactions and thermodynamic effects were ignored. Here, we applied the computational electrochemistry method based on ab initio molecular dynamics and thermodynamic integration to calculate redox potentials of quinones and their derivatives. The calculated results are in excellent agreement with experimental data. We mixed side chains to tune their reduction potentials and found that solvation interactions and entropy effects play a significant role in side-chain engineering. On the basis of our calculations, we proposed several high-performance negative and positive electrolytes. Our first-principles study paves the way toward the development of large-scale and sustainable electrical energy storage.

Computing Surface Acidity Constants of Proton Hopping Groups from Density Functional Theory-Based Molecular Dynamics: Application to the SnO2(110)/H2O Interface

Proton transfer at metal oxide/water interfaces plays an important role in electrochemistry, geochemistry, and environmental science. The key thermodynamic quantity to characterize this process is the surface acidity constant. An ab initio method that combines density functional theory-based molecular dynamics (DFTMD) and free energy perturbation theory has been established for computing surface acidity constants. However, it involves a reversible proton insertion procedure in which frequent proton hopping, e.g., for strong bases and some oxide surfaces (e.g., SnO₂), can cause instability issues in electronic structure calculation. In the original implementation, harmonic restraining potentials are imposed on all O-H bonds (denoted by the V_rH scheme) to prevent proton hopping and thus may not be applicable for systems involving spontaneous proton hopping. In this work, we introduce an improved restraining scheme with a repulsive potential V_rep to compute the surface acidities of systems in which proton hopping is spontaneous and fast. In this V_rep scheme, a Buckingham-type repulsive potential V_rep is applied between the deprotonation site and all other protons in DFTMD simulations. We first verify the V_rep scheme by calculating the pK_a values of H₂O and aqueous HS^- solution (i.e., strong conjugate bases) and then apply it to the SnO₂(110)/H₂O interface. It is found that the V_rep scheme leads to a prediction of the point of zero charge (PZC) of 4.6, which agrees well with experiment. The intrinsic individual pK_a values of the terminal five-coordinated Sn site (Sn_5cOH₂) and bridge oxygen site (Sn₂O_brH⁺) are 4.4 and 4.7, respectively, both being almost the same as the PZC. The similarity of the two pK_a values indicates that dissociation of terminal water has almost zero free energy at this proton hopping interface (i.e., partial water dissociation), as expected from the acid-base equilibrium on SnO₂.

Computational Ag/AgCl Reference Electrode from Density Functional Theory-Based Molecular Dynamics

We have developed a scheme to compute the standard potential of the Ag/AgCl reference electrode using density functional theory-based molecular dynamics, similar to the computational standard hydrogen electrode (SHE) developed by Cheng, Sulpizi, and Sprik [J. Chem. Phys. 2009, 131, 154504], with which our new computational reference electrode was compared. We have obtained a similar value of the potential of the Ag/AgCl electrode versus SHE to the experiment. The newly developed computational reference electrode will be extended to nonaqueous solvents in the future, where it will be used to predict standard equilibrium potentials to be compared with experimental data.

Assessment of Density Functional Theory in Predicting Interaction Energies between Water and Polycyclic Aromatic Hydrocarbons: from Water on Benzene to Water on Graphene

The interactions of water with polycyclic aromatic hydrocarbons, from benzene to graphene, are investigated using various exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. The accuracy of the different functionals is assessed through comparisons with random phase approximation (RPA) and coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] calculations. Diffusion Monte Carlo (DMC) data reported in the literature are also used for comparison. Relatively large variations are found in interaction energies predicted by different DFT models, with GGA functionals underestimating the interaction strength for configurations with the water oxygen pointing toward the aromatic molecules. The meta-GGA B97M-rV and range-separated hybrid, meta-GGA ωB97M-V functionals provide nearly quantitative agreement with CCSD(T) values for the water-benzene, water-coronene, and water-circumcoronene dimers, while RPA and DMC predict interaction energies that differ by up to ∼1 kcal/mol and ∼0.4 kcal/mol from the corresponding CCSD(T) values, respectively. Similar trends among GGA, meta-GGA, and hybrid functionals are observed for larger polycyclic aromatic hydrocarbons. By performing absolutely localized molecular orbital energy decomposition analyses (ALMO-EDA), it is found that, independently of the number of carbon atoms and exchange-correlation functional, the dominant contributions to the interaction energies between water and polycyclic aromatic hydrocarbon molecules are the electrostatic and dispersion terms while polarization and charge transfer effects are negligibly small. Calculations carried out with GGA and meta-GGA functionals indicate that, as the number of carbon atoms increases, the interaction energies slowly converge to the corresponding values obtained for an infinite graphene sheet.

A very sensitive and highly selective organic selector in CNTs composite chemiresistive for efficient differentiation of organic amine vapours

With the call for the IoE (Internet of Everything), stable and efficient electric noses/tongues have become the most critical part of the sensor network. Identifying target gases efficiently and rapidly at ambient air becomes a focus on sensor research. We designed a chemiresistive sensor based on a composite of a specific selector and single-walled carbon nanotubes (SWNTs) for the detection and differentiation of organic amine vapours in air (25 ℃, 55% RH). The synergetic combination of F4-TCNQ (2,3,5,6-Tetrafluoro-7,7',8,8'-tetracyanoquinodimethane) and SWCNTs could modulate the electrical properties of sensor leading to the enhancement of response up to ppb-level for primary amine vapor detection. Different from traditional chemiresistive sensor, this sensing materials exhibit unique differences in response to different types of amines thought different mechanisms. We have proven the practical possibilities through the detection of the simulated complexed environmental atmosphere in industrial production. Furthermore, we explored the working mechanism of high-performance sensors, which could provide theoretical guidance for sensor design for more commercial applications. This study provided a simple, convenient, and highly efficient practical method for organic amine detection at ambient air for real-life applications.

Dynamics of the Bulk Hydrated Electron from Many-Body Wave-Function Theory

The structure of the hydrated electron is a matter of debate as it evades direct experimental observation owing to the short life time and low concentrations of the species. Herein, the first molecular dynamics simulation of the bulk hydrated electron based on correlated wave‐function theory provides conclusive evidence in favor of a persistent tetrahedral cavity made up by four water molecules, and against the existence of stable non‐cavity structures. Such a cavity is formed within less than a picosecond after the addition of an excess electron to neat liquid water, with less regular cavities appearing as intermediates. The cavities are bound together by weak H−H bonds, the number of which correlates well with the number of coordinated water molecules, each type of cavity leaving a distinct spectroscopic signature. Simulations predict regions of negative spin density and a gyration radius that are both in agreement with experimental data.

Determining Potentials of Zero Charge of Metal Electrodes versus the Standard Hydrogen Electrode from Density-Functional-Theory-Based Molecular Dynamics

We develop a computationally efficient scheme to determine the potentials of zero charge (PZC) of metal-water interfaces with respect to the standard hydrogen electrode. We calculate the PZC of Pt(111), Au(111), Pd(111) and Ag(111) at a good accuracy using this scheme. Moreover, we find that the interface dipole potentials are almost entirely caused by charge transfer from water to the surfaces, the magnitude of which depends on the bonding strength between water and the metals, while water orientation hardly contributes at the PZC conditions.

Electronic structure of aqueous solutions: Bridging the gap between theory and experiments

Predicting the electronic properties of aqueous liquids has been a long-standing challenge for quantum mechanical methods. However, it is a crucial step in understanding and predicting the key role played by aqueous solutions and electrolytes in a wide variety of emerging energy and environmental technologies, including battery and photoelectrochemical cell design. We propose an efficient and accurate approach to predict the electronic properties of aqueous solutions, on the basis of the combination of first-principles methods and experimental validation using state-of-the-art spectroscopic measurements. We present results of the photoelectron spectra of a broad range of solvated ions, showing that first-principles molecular dynamics simulations and electronic structure calculations using dielectric hybrid functionals provide a quantitative description of the electronic properties of the solvent and solutes, including excitation energies. The proposed computational framework is general and applicable to other liquids, thereby offering great promise in understanding and engineering solutions and liquid electrolytes for a variety of important energy technologies.

Nuclear Quantum Effects on Aqueous Electron Attachment and Redox Properties

Nuclear quantum effects (NQEs) on the reduction and oxidation properties of small aqueous species (CO₂, HO₂, and O₂) are quantified and rationalized by first-principles molecular dynamics and thermodynamic integration. Vertical electron attachment, or electron affinity, and detachment energies (VEA and VDE) are strongly affected by NQEs, decreasing in absolute value by 0.3 eV going from a classical to a quantum description of the nuclei. The effect is attributed to NQEs that lessen the solvent response upon oxidation/reduction. The reduction of solvent reorganization energy is expected to be general for small solutes in water. In the thermodynamic integral that yields the free energy of oxidation/reduction, these large changes enter with opposite sign, and only a small net effect (0.1 eV) remains. This is not obvious for CO₂, where the integrand is strongly influenced by NQEs due to the onset of interaction of the reduced orbital with the conduction band of the liquid during thermodynamic integration. We conclude that NQEs might not have to be included in the computation of redox potentials, unless high accuracy is needed, but are important for VEA and VDE calculations.

Free Energies of Redox Half-Reactions from First-Principles Calculations

Quantitative prediction of the energetics of redox half-reactions is still a challenge for modern computational chemistry. Here, we propose a simple scheme for reliable calculations of vertical ionization and attachment energies, as well as of redox potentials of solvated molecules. The approach exploits linear response approximation in the context of explicit solvent simulations with spherical boundary conditions. It is shown that both vertical ionization energies and vertical electron affinities, and, consequently redox potentials, exhibit linear dependence on the inverse radius of the solvation sphere. The explanation of the linear dependence is provided, and an extrapolation scheme is suggested. The proposed approach accounts for the specific short-range interactions within hybrid DFT and effective fragment potential approach as well as for the asymptotic system-size effects. The computed vertical ionization energies and redox potentials are in excellent agreement with the experimental values.

Redox potentials of aryl derivatives from hybrid functional based first principles molecular dynamics

We report the redox potentials of a set of organic aryl molecules, including quinones, juglone, tyrosine and tryptophan, calculated using a first principles molecular dynamics (FPMD) based method.