Joint density-functional theory: ab initio study of Cr2O3 surface chemistry in solution

Advancing electrochemical nitrogen reduction: Efficacy of two-dimensional SiP layered structures with single-atom transition metal catalysts

Researchers are interested in single-atom catalysts with atomically scattered metals relishing the enhanced electrocatalytic activity for nitrogen reduction and 100 % metal atom utilization. In this paper, we investigated 18 transition metals (TM) spanning 3d to 5d series as efficient nitrogen reduction reaction (NRR) catalysts on defective 2D SiPV layered structures through first-principles calculation. A systematic screening identified Mo@SiPV, Nb@SiPV, Ta@SiPV and W@SiPV as superior, demonstrating enhanced ammonia synthesis with significantly lower limiting potentials (-0.25, -0.45, -0.49 and -0.15 V, respectively), compared to the benchmark -0.87 eV for the defective SiP. In addition, the descriptor ΔG*N was introduced to establish the relationship between the different NRR intermediates, and the volcano plot of the limiting potentials were determined for their potential-determining steps (PDS). Remarkably, the limiting voltage of the NRR possesses a good linear relationship with the active center TM atom Ɛd, which is a reliable descriptor for predicting the limiting voltage. Furthermore, we verified the stability (using Ab Initio Molecular Dynamics - AIMD) and high selectivity (UL(NRR)-UL(HER) > -0.5 V) of these four catalysts in vacuum and solvent environments. This study systematically demonstrates the strong catalytic potential of 2D TM@SiPV(TM = Mo, Nb, Ta, W) single-atom catalysts for nitrogen reduction electrocatalysis.

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.

PW-SMD: A Plane-Wave Implicit Solvation Model Based on Electron Density for Surface Chemistry and Crystalline Systems in Aqueous Solution

Electron density-based implicit solvation models are a class of techniques for quantifying solvation effects and calculating free energies of solvation without an explicit representation of solvent molecules. Integral to the accuracy of solvation modeling is the proper definition of the solvation shell separating the solute molecule from the solvent environment, allowing for a physical partitioning of the free energies of solvation. Unlike state-of-the-art implicit solvation models for molecular quantum chemistry calculations, e.g., the solvation model based on solute electron density (SMD), solvation models for systems under periodic boundary conditions with plane-wave (PW) basis sets have been limited in their accuracy. Furthermore, a unified implicit solvation model with both homogeneous solution-phase and heterogeneous interfacial structures treated on equal footing is needed. In order to address this challenge, we developed a high-accuracy solvation model for periodic PW calculations that is applicable to molecular, ionic, interfacial, and bulk-phase chemistry. Our model, PW-SMD, is an extension of the SMD molecular solvation model to periodic systems in water. The free energy of solvation is partitioned into the electrostatic and cavity-dispersion-solvent structure (CDS) contributions. The electrostatic contributions of the solvation shell surrounding solute structures are parametrized based on their geometric and physical properties. In addition, the nonelectrostatic contribution to the solvation energy is accounted for by extending the CDS formalism of SMD to incorporate periodic boundary conditions. We validate the accuracy and robustness of our solvation model by comparing predicted solvation free energies against experimental data for molecular and ionic systems, carved-cluster composite energetic models of solvated reaction energies and barriers on surface systems, and deep-learning-accelerated ab initio molecular dynamics (AIMD). Our developed periodic implicit solvation model shows significantly improved accuracy compared to previous work (namely, solvation models in aqueous solution) and can be applied to simulate solvent effects in a wide range of surface and crystalline materials.

Coverage enhancement accelerates acidic CO2 electrolysis at ampere-level current with high energy and carbon efficiencies

Acidic CO2 electroreduction (CO2R) using renewable electricity holds promise for high-efficiency generation of storable liquid chemicals with up to 100% CO2 utilization. However, the strong parasitic hydrogen evolution reaction (HER) limits its selectivity and energy efficiency (EE), especially at ampere-level current densities. Here we present that enhancing CO2R intermediate coverage on catalysts promotes CO2R and concurrently suppresses HER. We identified and engineered robust Cu6Sn5 catalysts with strong *OCHO affinity and weak *H binding, achieving 91% Faradaic efficiency (FE) for formic acid (FA) production at 1.2 A cm-2 and pH 1. Notably, the single-pass carbon efficiency reaches a new benchmark of 77.4% at 0.5 A cm-2 over 300 hours. In situ electrochemical Fourier-transform infrared spectroscopy revealed Cu6Sn5 enhances *OCHO coverage ~2.8× compared to Sn at pH 1. Using a cation-free, solid-state-electrolyte-based membrane-electrode-assembly, we produce 0.36 M pure FA at 88% FE over 130 hours with a marked full-cell EE of 37%.

Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.

Electrochemical Carbon Dioxide Reduction to Ethylene: From Mechanistic Understanding to Catalyst Surface Engineering

Electrochemical carbon dioxide reduction reaction (CO2RR) provides a promising way to convert CO2 to chemicals. The multicarbon (C2+) products, especially ethylene, are of great interest due to their versatile industrial applications. However, selectively reducing CO2 to ethylene is still challenging as the additional energy required for the C-C coupling step results in large overpotential and many competing products. Nonetheless, mechanistic understanding of the key steps and preferred reaction pathways/conditions, as well as rational design of novel catalysts for ethylene production have been regarded as promising approaches to achieving the highly efficient and selective CO2RR. In this review, we first illustrate the key steps for CO2RR to ethylene (e.g., CO2 adsorption/activation, formation of *CO intermediate, C-C coupling step), offering mechanistic understanding of CO2RR conversion to ethylene. Then the alternative reaction pathways and conditions for the formation of ethylene and competitive products (C1 and other C2+ products) are investigated, guiding the further design and development of preferred conditions for ethylene generation. Engineering strategies of Cu-based catalysts for CO2RR-ethylene are further summarized, and the correlations of reaction mechanism/pathways, engineering strategies and selectivity are elaborated. Finally, major challenges and perspectives in the research area of CO2RR are proposed for future development and practical applications.

Fractional Charge Density Functional Theory and Its Application to the Electro-inductive Effect

We employed the chemical potential equalization principle to demonstrate that fractional electrons are involved in the electro-inductive effect as well as the vibrational Stark effect. By the chemical potential model, we were able to deduce that the frontier molecular orbitals of immobilized molecules can provide valuable insight into these effects. To further understand and quantify these findings, we introduced fractional charge density functional theory (FC-DFT), a canonical ensemble approach for open systems. This method allows for the calculation of electronic energies, nuclear gradients, and the Hessian matrix of fractional electronic systems. To correct the spurious delocalization error commonly found in approximate density functionals for small systems, we imposed the Perdew-Parr-Levy-Balduz (PPLB) condition through linear interpolation of two adjacent integer points (LI-FC-DFT). Although this approach is relatively simple in terms of molecular modeling, the results obtained through LI-FC-DFT calculations predict the same trend seen in experimental reactivity and the frequency change of immobilized molecules.

Adaptive Electric Fields Embedded Electrochemical Barrier Calculations

Electrochemical interfaces are grand canonical ensembles of varying electrons. Simulating them by standard first-principles methods is a challenging task, since the number of electrons (or charge) is fixed in the calculation. Under the constant charge framework, we developed a constant potential simulation method realized by adding an adaptive electric field to a charge neutral system. Electric field is the controlling variable. In addition, we defined an internal reversible hydrogen electrode potential (ϕ_IRHE), which can ensure the model independence of our method. To validate our method, the reaction energies of some electrochemical reactions are calculated, the results are comparable with the computational hydrogen electrode model and experiments. At last, the evolution of transition state structures and charge transfer coefficients of some electrochemical reactions on Ag(111) surface were discussed by our method.

Unravelling the 2e- ORR Activity Induced by Distance Effect on Main-Group Metal InN4 Surface Based on First Principles

The p-electron-dominated main-group metals (Sb, Se, In, etc.) have recently been reported to possess excellent oxygen reduction reaction (ORR) activity by means of heteroatom doping into graphene. However, on these main group metal surfaces, other approaches especially the distance effect to modulate catalytic activity are rarely involved. In this work, the origin of excellent 2e^- ORR catalytic activity of graphene-supported InN₄ moiety by tuning the distance between metallic In atoms is thoroughly investigated by employing the first-principles calculations. Our DFT calculations show that the 2e^- ORR catalytic activity strongly depends on the crystal orbital Hamilton population (COHP) between In and O atoms. This work is useful for the rational design of main group metal single atom electrocatalysts.

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pK_a values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

Advancement of the Homogeneous Background Method for the Computational Simulation of Electrochemical Interfaces

Computational studies of electrochemical interfaces based on density-functional theory (DFT) play an increasingly important role in the present research on electrochemical processes for energy conversion and storage. The homogeneous background method (HBM) offers a straightforward approach to charge the electrochemical system within DFT simulations, but it typically requires the specification of the active fraction of excess electrons based on a certain choice of the electrode-electrolyte boundary location, which can be difficult in the presence of electrode-surface adsorbates or explicit solvent molecules. In this work, we present a methodological advancement of the HBM, both facilitating and extending its applicability. The advanced version requires neither energy corrections nor the specification of the active fraction of excess electrons, providing a versatile and readily available method for the simulation of charged interfaces when adsorbates or explicit solvent molecules are present. Our computational DFT results for Pt(111), Au(111), and Li(100) metal electrodes in high-dielectric-constant solvents demonstrate an excellent agreement in the interfacial charging characteristics obtained from simulations with the advanced HBM in comparison with the (linearized) Poisson-Boltzmann model (PBM).

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.

Ab Initio Simulations of Water/Metal Interfaces

Structures and processes at water/metal interfaces play an important technological role in electrochemical energy conversion and storage, photoconversion, sensors, and corrosion, just to name a few. However, they are also of fundamental significance as a model system for the study of solid-liquid interfaces, which requires combining concepts from the chemistry and physics of crystalline materials and liquids. Particularly interesting is the fact that the water-water and water-metal interactions are of similar strength so that the structures at water/metal interfaces result from a competition between these comparable interactions. Because water is a polar molecule and water and metal surfaces are both polarizable, explicit consideration of the electronic degrees of freedom at water/metal interfaces is mandatory. In principle, ab initio molecular dynamics simulations are thus the method of choice to model water/metal interfaces, but they are computationally still rather demanding. Here, ab initio simulations of water/metal interfaces will be reviewed, starting from static systems such as the adsorption of single water molecules, water clusters, and icelike layers, followed by the properties of liquid water layers at metal surfaces. Technical issues such as the appropriate first-principles description of the water-water and water-metal interactions will be discussed, and electrochemical aspects will be addressed. Finally, more approximate but numerically less demanding approaches to treat water at metal surfaces from first-principles will be briefly discussed.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

Effect of Surface Oxidation on Oxidative Propane Dehydrogenation over Chromia: An Ab Initio Multiscale Kinetic Study

An increasingly utilized way for the production of propene is propane dehydrogenation. The reaction presents an alternative to conventional processes based on petroleum resources. In this work, we investigate theoretically how Cr₂O₃ catalyzes this reaction in oxidative and reducing environments. Although previous studies showed that the reduced catalyst is selective for the non-oxidative dehydrogenation of propane, real operating conditions are oxidative. Herein, we use multiscale modeling to investigate the difference between the oxidized and reduced catalyst and their performance. The complete reaction pathway for propane dehydrogenation, including C-C cracking, formation of side products (propyne, ethane, ethylene, acetylene, and methane), and catalyst coking on oxidized and reduced surfaces of α-Cr₂O₃(0001), is calculated using density functional theory with the Hubbard correction. Parameters describing adsorption, desorption, and surface reactions are used in a kinetic Monte Carlo simulation, which employed industrially relevant conditions (700-900 K, pressures up to 2 bar, and varying oxidants: N₂O, O₂, and none). We observe that over the reduced surface, propene and hydrogen form with high selectivity. When oxidants are used, the surface is oxidized, which changes the reaction mechanism and kinetics. During a much faster reaction, H₂O forms as a coproduct in a Mars-van Krevelen cycle. Additionally, CO₂ is also formed, which represents waste and adversely affects the selectivity. It is shown that the oxidized surface is much more active but prone to the formation of CO₂, while the reduced surface is less active but highly selective toward propene. Moreover, the effect of the oxidant used is investigated, showing that N₂O is preferred to O₂ due to higher selectivity and less catalyst coking. We show that there exists an optimum degree of surface oxidation, where the yield of propene is maximized. The coke, which forms during the reaction, can be burnt away as CO₂ with oxygen.

Ab initiomodelling of interfacial electrochemical properties: beyond implicit solvation limitations

First-principles calculations are an important tool to investigate the complex processes occurring at solid/liquid interfaces which are at the heart of modern technologies. Currently, capturing the whole electrochemical environment at an interface, including the applied potential and solvation, still remains challenging as it necessitates to couple different approaches whose interactions are not fully understood. In this work, a grand canonical density functional theory approach is coupled with solvation models to investigate the electrochemical interfaces under applied potential. We show that a parametrized polarizable continuum model (PCM) which represent solvation in a mean field approach by a continuous polarizable media, possesses catastrophic limitations for the modelling of ionic and charged interfaces. We reveal the origin of PCM instabilities under chemical or electrochemical strong oxidation to be the consequence of a phase transition in the surface Li electronic structure. Thus, PCM undergoes an unphysical response to this phase transition by penetrating within the atomic radius of surface Li atoms. To recover a physical response, an explicit first solvation shell has to be included in addition to the PCM in order to properly describe the electrochemistry of the interface. The Fukui functions show that the first solvation shell becomes involved in the redox process as solvent electron doublet is transferred to the acidic Li⁺. If another explicit solvent layer is added, the interface electrochemical properties become independent of the PCM parameters: in particular capacitance can then be computed from a parameter-free electrochemical approach. This is an important conclusion as the experimental electrochemical capacitance are not easily found and thus the parametrization of the PCM for electrochemical interface can be difficult. This approach can easily be applied to investigate electrochemical properties at the atomic scale and generalized to any electrochemical device for which interfaces play a crucial role.

Electronic structure behavior of PbO2, IrO2, and SnO2 metal oxide surfaces (110) with dissociatively adsorbed water molecules as a function of the chemical potential

A detailed analysis of the electronic structure of three different electrochemical interfaces as a function of the chemical potential (μ) is performed using the grand canonical density functional theory in the joint density functional theory formulation. Changes in the average number of electrons and the density of states are also described. The evaluation of the global softness, which measures the tendency of the system to gain or lose electrons, is straightforward under this formalism. The observed behavior of these quantities depends on the electronic nature of the electrochemical interfaces.

Tackling Solvent Effects by Coupling Electronic and Molecular Density Functional Theory

Solvation effects can have a tremendous influence on chemical reactions. However, precise quantum chemistry calculations are most often done either in vacuum neglecting the role of the solvent or using continuum solvent model ignoring its molecular nature. We propose a new method coupling a quantum description of the solute using electronic density functional theory with a classical grand-canonical treatment of the solvent using molecular density functional theory. Unlike a previous work, both densities are minimized self-consistently, accounting for mutual polarization of the molecular solvent and the solute. The electrostatic interaction is accounted using the full electron density of the solute rather than fitted point charges. The introduced methodology represents a good compromise between the two main strategies to tackle solvation effects in quantum calculation. It is computationally more effective than a direct quantum mechanics/molecular mechanics coupling, requiring the exploration of many solvent configurations. Compared to continuum methods, it retains the full molecular-level description of the solvent. We validate this new framework onto two usual benchmark systems: a water solvated in water and the symmetrical nucleophilic substitution between chloromethane and chloride in water. The prediction for the free energy profiles are not yet fully quantitative compared to experimental data, but the most important features are qualitatively recovered. The method provides a detailed molecular picture of the evolution of the solvent structure along the reaction pathway.

Sulfate, Bisulfate, and Hydrogen Co-adsorption on Pt(111) and Au(111) in an Electrochemical Environment

The co-adsorption of sulfate, bisulfate and hydrogen on Pt(111) and Au(111) electrodes was studied based on periodic density functional calculations with the aqueous electrolyte represented by both explicit and implicit solvent models. The influence of the electrochemical control parameters such as the electrode potential and pH was taken into account in a grand-canonical approach. Thus, phase diagrams of the stable coadsorption phases as a function of the electrochemical potential and Pourbaix diagrams have been derived which well reproduce experimental findings. We demonstrate that it is necessary to include explicit water molecules in order to determine the stable adsorbate phases as the (bi)sulfate adsorbates rows become significantly stabilized by bridging water molecules.

Understanding Reaction Pathways in High Dielectric Electrolytes Using β-Mo2C as a Catalyst for Li-CO2 Batteries

The rechargeable Li-CO₂ battery has attracted considerable attention in recent years because of its carbon dioxide (CO₂) utilization and because it represents a practical Li-air battery. As with other battery systems such as the Li-ion, Li-O₂, and Li-S battery systems, understanding the reaction pathway is the first step to achieving high battery performance because the performance is strongly affected by reaction intermediates. Despite intensive efforts in this area, the effect of material parameters (e.g., the electrolyte, the cathode, and the catalyst) on the reaction pathway in Li-CO₂ batteries is not yet fully understood. Here, we show for the first time that the discharge reaction pathway of a Li-CO₂ battery composed of graphene nanoplatelets/beta phase of molybdenum carbide (GNPs/β-Mo₂C) is strongly influenced by the dielectric constant of its electrolyte. Calculations using the continuum solvents model show that the energy of adsorption of oxalate (C₂O₄^2-) onto Mo₂C under the low-dielectric electrolyte tetraethylene glycol dimethyl ether is lower than that under the high-dielectric electrolyte N,N-dimethylacetamide (DMA), indicating that the electrolyte plays a critical role in determining the reaction pathway. The experimental results show that under the high-dielectric DMA electrolyte, the formation of lithium carbonate (Li₂CO₃) as a discharge product is favorable because of the instability of the oxalate species, confirming that the dielectric properties of the electrolyte play an important role in the formation of the discharge product. The resulting Li-CO₂ battery exhibits improved battery performance, including a reduced overpotential and a remarkable discharge capacity as high as 14,000 mA h g^-1 because of its lower internal resistance. We believe that this work provides insights for the design of Li-CO₂ batteries with enhanced performance for practical Li-air battery applications.

Implicit self-consistent electrolyte model in plane-wave density-functional theory

The ab initio computational treatment of electrochemical systems requires an appropriate treatment of the solid/liquid interfaces. A fully quantum mechanical treatment of the interface is computationally demanding due to the large number of degrees of freedom involved. In this work, we develop a computationally efficient model where the electrode part of the interface is described at the density-functional theory (DFT) level, and the electrolyte part is represented through an implicit solvation model based on the Poisson-Boltzmann equation. We describe the implementation of the linearized Poisson-Boltzmann equation into the Vienna Ab initio Simulation Package, a widely used DFT code, followed by validation and benchmarking of the method. To demonstrate the utility of the implicit electrolyte model, we apply it to study the surface energy of Cu crystal facets in an aqueous electrolyte as a function of applied electric potential. We show that the applied potential enables the control of the shape of nanocrystals from an octahedral to a truncated octahedral morphology with increasing potential.

Interface Properties of the Partially Oxidized Pt(111) Surface Using Hybrid DFT-Solvation Models

This article reports a theoretical-computational effort to model the interface between an oxidized platinum surface and aqueous electrolyte. It strives to account for the impact of the electrode potential, formation of surface-bound oxygen species, orientational ordering of near-surface solvent molecules, and metal surface charging on the potential profile along the normal direction. The computational scheme is based on the DFT/ESM-RISM method to simulate the charged Pt(111) surface with varying number of oxygen adatoms in acidic solution. This hybrid solvation method is known to qualitatively reproduce bulk metal properties like the work function. However, the presented calculations reveal that vital interface properties such as the electrostatic potential at the outer Helmholtz plane are highly sensitive to the position of the metal surface slab relative to the DFT-RISM boundary region. Shifting the relative position of the slab also affects the free energy of the system. It follows that there is an optimal distance for the first solvent layer within the ESM-RISM framework, which could be found by optimizing the position of the frozen Pt(111) slab. As it stands, manual sampling of the position of the slab is impractical and betrays the self-consistency of the method. Based on this understanding, we propose the implementation of a free energy optimization scheme of the relative position of the slab in the DFT-RISM boundary region. This optimization scheme could considerably increase the applicability of the hybrid method.

Elucidating Molecular Interactions in Glycerol Adsorption at the Metal-Water Interface with Density Functional Theory

Glycerol is an extremely versatile platform molecule for chemical and fuel production, as evidenced by successful demonstrations in electrochemical and thermochemical processes, where key catalytic chemistries occur at the solid-liquid interface. Despite the remarkable progress made in enriching the first-principles-based computational tool set to reveal and characterize solvent structures in the past decade, techniques for realistic and efficient molecular-level modeling to study aqueous-phase glycerol chemistry are still far from mature. Many aqueous-phase catalytic systems are deemed too complex for routine modeling because of their highly correlated structures at the heterogeneous solid-liquid interface. This invited feature article merges recent developments in quantum mechanical solvation models and oxygenated hydrocarbon conversion chemistry by revisiting the molecular interactions of adsorbed glycerol and its dehydrogenation intermediates at the water-metal interface. Explicit participation of water through the establishment of water-adsorbate, water-water, and water-metal interactions on Pt(111) was investigated using density functional theory. In periodic models, the adsorption favors networklike structures with adsorbates as nodal points linked by coadsorbed water molecules. We also showed that these adsorption patterns actually preserve the original bond-order-based scaling relationship framework established without the consideration of solvent. This behavior can be exploited to improve computational efficiency for future analysis of catalytic polyol conversions in the aqueous-phase environment.

Grand-canonical approach to density functional theory of electrocatalytic systems: Thermodynamics of solid-liquid interfaces at constant ion and electrode potentials

Properties of solid-liquid interfaces are of immense importance for electrocatalytic and electrochemical systems, but modeling such interfaces at the atomic level presents a serious challenge and approaches beyond standard methodologies are needed. An atomistic computational scheme needs to treat at least part of the system quantum mechanically to describe adsorption and reactions, while the entire system is in thermal equilibrium. The experimentally relevant macroscopic control variables are temperature, electrode potential, and the choice of the solvent and ions, and these need to be explicitly included in the computational model as well; this calls for a thermodynamic ensemble with fixed ion and electrode potentials. In this work, a general framework within density functional theory (DFT) with fixed electron and ion chemical potentials in the grand canonical (GC) ensemble is established for modeling electrocatalytic and electrochemical interfaces. Starting from a fully quantum mechanical description of multi-component GC-DFT for nuclei and electrons, a systematic coarse-graining is employed to establish various computational schemes including (i) the combination of classical and electronic DFTs within the GC ensemble and (ii) on the simplest level a chemically and physically sound way to obtain various (modified) Poisson-Boltzmann (mPB) implicit solvent models. The detailed and rigorous derivation clearly establishes which approximations are needed for coarse-graining as well as highlights which details and interactions are omitted in vein of computational feasibility. The transparent approximations also allow removing some of the constraints and coarse-graining if needed. We implement various mPB models within a linear dielectric continuum in the GPAW code and test their capabilities to model capacitance of electrochemical interfaces as well as study different approaches for modeling partly periodic charged systems. Our rigorous and well-defined DFT coarse-graining scheme to continuum electrolytes highlights the inadequacy of current linear dielectric models for treating properties of the electrochemical interface.

A first-principles roadmap and limits to design efficient supercapacitor electrode materials

A road map to guide researchers to predict the desired properties is presented based on the DFT calculations to allow researchers decide which property of the material they wish to predict or develop and how to choose the proper DFT route to do so.

Theoretical Modeling of Electronic Excitations of Gas-Phase and Solvated TiO2 Nanoclusters and Nanoparticles of Interest in Photocatalysis

The optical absorption spectra of (TiO₂) _n, nanoclusters ( n = 1-20) and nanoparticles ( n = 35, 84) have been calculated from the frequency-dependent dielectric function in the independent particle approximation under the framework of density functional theory. The PBE generalized gradient approach based functional, the so-called PBE+ U method and the PBE0 and PBEx hybrid functionals-containing 25% and 12.5% of nonlocal Fock exchange, respectively-have been used. The simulated spectra have been obtained in the gas phase and in water on previously PBE0 optimized atomic structures. The effect of the solvent has been accounted for by using an implicit water solvation model. For the smallest nanoclusters, the spectra show discrete peaks, whereas for the largest nanoclusters and for the nanoparticles they resemble a continuum absorption band. In the gas phase and for a given density functional, the onset of the absorption (optical gap, O_gap) remains relatively constant for all nanoparticle sizes although it increases with the percentage of nonlocal Fock exchange, as expected. For all tested functionals, the tendency of O_gap in water is very similar to that observed in the gas phase with an almost constant upshift. For comparison, the optical gap has also been calculated at the TD-DFT level with the PBEx functional in the gas phase and in water. Both approaches agree reasonably well although the TD-DFT gap values are lower than those derived from the dielectric-function. Overall, the position of the spectral maxima and the width of the spectra are relatively constant and independent of particle size which may have implications in the understanding of photocatalysis by TiO₂.

High-performance bifunctional porous non-noble metal phosphide catalyst for overall water splitting

Water electrolysis is an advanced energy conversion technology to produce hydrogen as a clean and sustainable chemical fuel, which potentially stores the abundant but intermittent renewable energy sources scalably. Since the overall water splitting is an uphill reaction in low efficiency, innovative breakthroughs are desirable to greatly improve the efficiency by rationally designing non-precious metal-based robust bifunctional catalysts for promoting both the cathodic hydrogen evolution and anodic oxygen evolution reactions. We report a hybrid catalyst constructed by iron and dinickel phosphides on nickel foams that drives both the hydrogen and oxygen evolution reactions well in base, and thus substantially expedites overall water splitting at 10 mA cm-2 with 1.42 V, which outperforms the integrated iridium (IV) oxide and platinum couple (1.57 V), and are among the best activities currently. Especially, it delivers 500 mA cm-2 at 1.72 V without decay even after the durability test for 40 h, providing great potential for large-scale applications.

Force Field for Water over Pt(111): Development, Assessment, and Comparison

Metal/water interfaces are key in many natural and industrial processes, such as corrosion, atmospheric, or environmental chemistry. Even today, the only practical approach to simulate large interfaces between a metal and water is to perform force-field simulations. In this work, we propose a novel force field, GAL17, to describe the interaction of water and a Pt(111) surface. GAL17 builds on three terms: (i) a standard Lennard-Jones potential for the bonding interaction between the surface and water, (ii) a Gaussian term to improve the surface corrugation, and (iii) two terms describing the angular dependence of the interaction energy. The 12 parameters of this force field are fitted against a set of 210 adsorption geometries of water on Pt(111). The performance of GAL17 is compared to several other approaches that have not been validated against extensive first-principles computations yet. Their respective accuracy is evaluated on an extended set of 802 adsorption geometries of H₂O on Pt(111), 52 geometries derived from icelike layers, and an MD simulation of an interface between a c(4 × 6) Pt(111) surface and a water layer of 14 Å thickness. The newly developed GAL17 force field provides a significant improvement over previously existing force fields for Pt(111)/H₂O interactions. Its well-balanced performance suggests that it is an ideal candidate to generate relevant geometries for the metal/water interface, paving the way to a representative sampling of the equilibrium distribution at the interface and to predict solvation free energies at the solid/liquid interface.

Understanding the MXene Pseudocapacitance

MXenes have attracted great attention as next-generation capacitive energy-storage materials, but the mechanisms underlying their pseudocapacitive behavior are not well understood. Here we provide a theoretical description of the surface redox process of Ti₃C₂T _x (T = O, OH), a prototypical MXene, in 1 M H₂SO₄ electrolyte, based on joint density functional theory with an implicit solvation model and the analysis of Gibbs free energy under a constant-electrode potential. From the dependence of the O/OH ratio (or the surface H coverage) and the surface charge on the applied potential, we obtain a clear picture of the capacitive energy-storage mechanism of Ti₃C₂T _x that shows good agreement with previous experimental findings in terms of the integral capacitance and Ti valence change. We find a voltage-dependent redox/double-layer co-charging behavior: the capacitive mechanism is dominated by the redox process, but the electric double-layer charge works against the redox process. This new insight may be useful in improving the capacitance of MXenes.

Polymer-free nanofibers from vanillin/cyclodextrin inclusion complexes: high thermal stability, enhanced solubility and antioxidant property

Vanillin/cyclodextrin inclusion complex nanofibers (vanillin/CD-IC NFs) were successfully obtained from three modified CD types (HPβCD, HPγCD and MβCD) in three different solvent systems (water, DMF and DMAc) via an electrospinning technique without using a carrier polymeric matrix. Vanillin/CD-IC NFs with uniform and bead-free fiber morphology were successfully produced and their free-standing nanofibrous webs were obtained. The polymer-free CD/vanillin-IC-NFs allow us to accomplish a much higher vanillin loading (∼12%, w/w) when compared to electrospun polymeric nanofibers containing CD/vanillin-IC (∼5%, w/w). Vanillin has a volatile nature yet, after electrospinning, a significant amount of vanillin was preserved due to complex formation depending on the CD types. Maximum preservation of vanillin was observed for vanillin/MβCD-IC NFs which is up to ∼85% w/w, besides, a considerable amount of vanillin (∼75% w/w) was also preserved for vanillin/HPβCD-IC NFs and vanillin/HPγCD-IC NFs. Phase solubility studies suggested a 1 : 1 molar complexation tendency between guest vanillin and host CD molecules. Molecular modelling studies and experimental findings revealed that vanillin : CD complexation was strongest for MβCD when compared to HPβCD and HPγCD in vanillin/CD-IC NFs. For vanillin/CD-IC NFs, water solubility and the antioxidant property of vanillin was improved significantly owing to inclusion complexation. In brief, polymer-free vanillin/CD-IC NFs are capable of incorporating a much higher loading of vanillin and effectively preserve volatile vanillin. Hence, encapsulation of volatile active agents such as flavor, fragrance and essential oils in electrospun polymer-free CD-IC NFs may have potential for food related applications by integrating the particularly large surface area of NFs with the non-toxic nature of CD and inclusion complexation benefits, such as high temperature stability, improved water solubility and an enhanced antioxidant property, etc.

Computational Insights into Materials and Interfaces for Capacitive Energy Storage

Supercapacitors such as electric double-layer capacitors (EDLCs) and pseudocapacitors are becoming increasingly important in the field of electrical energy storage. Theoretical study of energy storage in EDLCs focuses on solving for the electric double-layer structure in different electrode geometries and electrolyte components, which can be achieved by molecular simulations such as classical molecular dynamics (MD), classical density functional theory (classical DFT), and Monte-Carlo (MC) methods. In recent years, combining first-principles and classical simulations to investigate the carbon-based EDLCs has shed light on the importance of quantum capacitance in graphene-like 2D systems. More recently, the development of joint density functional theory (JDFT) enables self-consistent electronic-structure calculation for an electrode being solvated by an electrolyte. In contrast with the large amount of theoretical and computational effort on EDLCs, theoretical understanding of pseudocapacitance is very limited. In this review, we first introduce popular modeling methods and then focus on several important aspects of EDLCs including nanoconfinement, quantum capacitance, dielectric screening, and novel 2D electrode design; we also briefly touch upon pseudocapactive mechanism in RuO2. We summarize and conclude with an outlook for the future of materials simulation and design for capacitive energy storage.

Cu metal embedded in oxidized matrix catalyst to promote CO2 activation and CO dimerization for electrochemical reduction of CO2

We propose and validate with quantum mechanics methods a unique catalyst for electrochemical reduction of CO₂ (CO₂RR) in which selectivity and activity of CO and C₂ products are both enhanced at the borders of oxidized and metallic surface regions. This Cu metal embedded in oxidized matrix (MEOM) catalyst is consistent with observations that Cu₂O-based electrodes improve performance. However, we show that a fully oxidized matrix (FOM) model would not explain the experimentally observed performance boost, and we show that the FOM is not stable under CO₂ reduction conditions. This electrostatic tension between the Cu⁺ and Cu⁰ surface sites responsible for the MEOM mechanism suggests a unique strategy for designing more efficient and selective electrocatalysts for CO₂RR to valuable chemicals (HCO_x), a critical need for practical environmental and energy applications.

Antibacterial electrospun zein nanofibrous web encapsulating thymol/cyclodextrin-inclusion complex for food packaging

Thymol (THY)/γ-Cyclodextrin(γ-CD) inclusion complex (IC) encapsulated electrospun zein nanofibrous webs (zein-THY/γ-CD-IC-NF) were fabricated as a food packaging material. The formation of THY/γ-CD-IC (1:1 and 2:1) was proved by experimental (X-ray diffraction (XRD), thermal gravimetric analysis (TGA), ¹H NMR) and computational techniques. THY/γ-CD-IC (2:1) exhibited higher preservation rate and stability than THY/γ-CD-IC (1:1). It is worth mentioning that zein-THY/γ-CD-IC-NF (2:1) preserved much more THY as observed in TGA and stability of THY/γ-CD-IC (2:1) was higher, as shown by a modelling study. Therefore, much more THY was released from zein-THY/γ-CD-IC-NF (2:1) than zein-THY-NF and zein-THY/γ-CD-IC-NF (1:1). Similarly, antibacterial activity of zein-THY/γ-CD-IC-NF (2:1) was higher than zein-THY-NF and zein-THY/γ-CD-IC-NF (1:1). It was demonstrated that zein-THY/γ-CD-IC-NF (2:1) was most effective in inhibiting the growth of bacteria on meat samples. These webs show potential application as an antibacterial food packaging material.

Full atomistic reaction mechanism with kinetics for CO reduction on Cu(100) from ab initio molecular dynamics free-energy calculations at 298 K

A critical step toward the rational design of new catalysts that achieve selective and efficient reduction of CO₂ to specific hydrocarbons and oxygenates is to determine the detailed reaction mechanism including kinetics and product selectivity as a function of pH and applied potential for known systems. To accomplish this, we apply ab initio molecular metadynamics simulations (AIMμD) for the water/Cu(100) system with five layers of the explicit solvent under a potential of -0.59 V [reversible hydrogen electrode (RHE)] at pH 7 and compare with experiment. From these free-energy calculations, we determined the kinetics and pathways for major products (ethylene and methane) and minor products (ethanol, glyoxal, glycolaldehyde, ethylene glycol, acetaldehyde, ethane, and methanol). For an applied potential (U) greater than -0.6 V (RHE) ethylene, the major product, is produced via the Eley-Rideal (ER) mechanism using H₂O + e^- The rate-determining step (RDS) is C-C coupling of two CO, with ΔG^‡ = 0.69 eV. For an applied potential less than -0.60 V (RHE), the rate of ethylene formation decreases, mainly due to the loss of CO surface sites, which are replaced by H*. The reappearance of C₂H₄ along with CH₄ at U less than -0.85 V arises from *CHO formation produced via an ER process of H* with nonadsorbed CO (a unique result). This *CHO is the common intermediate for the formation of both CH₄ and C₂H₄ These results suggest that, to obtain hydrocarbon products selectively and efficiency at pH 7, we need to increase the CO concentration by changing the solvent or alloying the surface.