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

Modification of metals and ligands in two-dimensional conjugated metal-organic frameworks for CO2 electroreduction: A combined density functional theory and machine learning study

Electrochemical carbon dioxide reduction reaction (CO2RR) is a promising technology to establish an artificial carbon cycle. Two-dimensional conjugated metal-organic frameworks (2D c-MOFs) with high electrical conductivity have great potential as catalysts. Herein, we designed a range of 2D c-MOFs with different transition metal atoms and organic ligands, TMNxO4-x-HDQ (TM = Cr∼Cu, Mo, Ru∼Ag, W∼Au; x  = 0, 2, 4; HDQ = hexadipyrazinoquinoxaline), and systematically studied their catalytic performance using density functional theory (DFT). Calculation results indicated that all of TMNxO4-x-HDQ structures possess good thermodynamic and electrochemical stability. Notably, among the examined 37 MOFs, 6 catalysts outperformed the Cu(211) surface in terms of catalytic activity and product selectivity. Specifically, NiN4-HDQ emerged as an exceptional electrocatalyst for CO production in CO2RR, yielding a remarkable low limiting potential (UL) of -0.04 V. CuN4-HDQ, NiN2O2-HDQ, and PtN2O2-HDQ also exhibited high activity for HCOOH production, with UL values of -0.27, -0.29, and -0.27 V, respectively, while MnN4-HDQ, and NiO4-HDQ mainly produced CH4 with UL values of -0.58 and -0.24 V, respectively. Furthermore, these 6 catalysts efficiently suppressed the competitive hydrogen evolution reaction. Machine learning (ML) analysis revealed that the key intrinsic factors influencing CO2RR performance of these 2D c-MOFs include electron affinity (EA), electronegativity (χ), the first ionization energy (Ie), p-band center of the coordinated N/O atom (εp), the radius of metal atom (r), and d-band center (εd). Our findings may provide valuable insights for the exploration of highly active and selective CO2RR electrocatalysts.

Machine learning-driven shortening the screening process towards high-performance nitrogen reduction reaction electrocatalysts with four-step screening strategy

The urgent need to prepare clean energy by environmentally friendly and efficient methods, which has led to widespread attention on electrocatalytic nitrogen reduction reaction (NRR) for ammonia production. At present, single atom catalytic nitrogen reduction has become the earliest promising method for industrial production due to its high atomic utilization rate, high selectivity, high controllability, and high stability. However, how to quickly screen catalysts with high catalytic efficiency and selectivity in single-atom catalysts (SACs) remains a challenge. Herein, the 29 SACs are constructed from C6N2 nanosheets doped with transition metals (TM@C6N2), which are analyzed for stability, adsorption performance, NRR catalytic activity, electronic properties, and competitiveness using first-principles calculations. The results show that Mo@C6N2 and Re@C6N2 exhibit the most outstanding catalytic performances, with limiting potentials (UL) of -0.29 and -0.31 V, respectively, in the solvent model. Machine learning is used to derive descriptors from the intrinsic features to predict the free energy changes for the potential-determining step. The importance of features is calculated, with the first ionisation energy (IE1) being the most significant influencing factor. Based on the guidance of machine learning and considering that IE1 is related to the ability of metal atoms to donate electrons, a four-step screening strategy using the Integrated Crystal Orbital Hamilton Populations (ICOHP) to screen catalysts instead of the traditional five-step screening not only improves the screening efficiency but also obtains completely consistent screening results. This work presents a new approach to predicting the catalytic performance of SACs and provides new insights into the influence of intrinsic properties on catalytic activity.

Computational screening of single-atom transition metals on boron-rich boron nitride nanosheets for efficient hydrogen evolution catalysis in all pH range

Low-cost and high-efficiency catalysts are of crucial importance for the electrocatalytic hydrogen evolution reaction (HER). Two-dimensional (2D) boron nitride (B-N) compounds formed by the combination of boron and nitrogen atoms of group III and V elements are promising candidates for electrocatalytic HER due to their significant electronic properties. Hence, an electrocatalyst is computer-aided designed with isolated single atoms of 3d, 4d, and 5d transition metals supported on 2D B-N (B2N, B5N3, and B7N5) monolayers to fabricate single-atom catalysts (SACs) with an excellent HER performance. Moreover, pH modulations are considered to improve the HER activity theoretically based on first-principles calculation. Our results indicate that B-N compounds surface doping with transition metal atoms can effectively enhance the HER catalytic performance over a wide range of pH. Among all SACs studied, Co-, Ti-, V-, Nb-, Ru-, Tc-, Zr-, and Os-embedded B2N, Sc-, Cr-, Mn-, Ti-, and Y-embedded B5N3, and Sc- and Mn-embedded B7N5 have excellent catalytic activity under acidic conditions, while Mo-, Ir-, Re-, Ta-, and W-embedded B2N and Ti- and Fe-embedded B7N5 show high catalytic activity under alkaline conditions. Interestingly, Hf@B2N and V@B5N3 systems exhibit promising catalytic activity under acidic, neutral, and alkaline conditions. Our work offers cost-effective candidates with a wide pH range HER performance for exploring ideal electrocatalysts.

Electrical monitoring of single-event protonation dynamics at the solid-liquid interface and its regulation by external mechanical forces

Detecting chemical reaction dynamics at solid-liquid interfaces is important for understanding heterogeneous reactions. However, there is a lack of exploration of interface reaction dynamics from the single-molecule perspective, which can reveal the intrinsic reaction mechanism underlying ensemble experiments. Here, single-event protonation reaction dynamics at a solid-liquid interface are studied in-situ using single-molecule junctions. Molecules with amino terminal groups are used to construct single-molecule junctions. An interfacial cationic state present after protonation is discovered. Real-time electrical measurements are used to monitor the reversible reaction between protonated and deprotonated states, thereby revealing the interfacial reaction mechanism through dynamic analysis. The protonation reaction rate constant has a linear positive correlation with proton concentration, whereas the deprotonation reaction rate constant has a linear negative correlation. In addition, external mechanical forces can effectively regulate the protonation reaction process. This work provides a single-molecule perspective for exploring interface science, which will contribute to the development of heterogeneous catalysis and electrochemistry.

Structure-Stability Relation of Single-Atom Catalysts under Operating Conditions of CO2 Reduction

Single-atom catalysts (SACs) have exhibited exceptional atomic efficiency and catalytic performance in various reactions but suffer poor stability. Understanding the structure-stability relation is the prerequisite for stability optimization but has been rarely explored due to complexity of the degradation process and reaction environments. Herein, we successfully established the structure-stability relation of N-doped carbon-supports SACs (MN4 SACs) under working conditions of CO2 reduction, by using advanced constant-potential density functional theory calculations. Systematic mechanism investigation that considered different factors identifies the key role of initial hydrogen adsorption on the coordination N atom in catalytic stability, where the feasibility of the adsorption eventually determines the leaching of the metal atom. On this basis, a simple descriptor consisting of electron number and electronegativity is constructed, realizing accurate and rapid prediction of the stability of SACs. Furthermore, strategies via modifying the local geometric structure to improve the stability without changing the active centers are proposed accordingly, which are supported by related experiments. These findings fill the current void in understanding SAC stability under practical working conditions, potentially advancing the widespread application of SACs in sustainable energy conversion systems.

Optimizing QLED Performance and Stability via the Surface Modification of PEDOT:PSS Experimental Insights and DFT Calculations

The presence of the acidic and weak ionic conductor polystyrenesulfonate (PSS) in poly(3,4-ethylenedioxythiophene:PSS (PEDOT:PSS) leads to degradation and limits the charge transfer within quantum dot light-emitting diodes (QLEDs). Two-step solvent treatment resulted in a 40% reduction of PSS, which could be attributed to ethylene glycol (EG) attenuating the ionic interactions between PSS and PEDOT via interacting with PSS through hydrogen bonding. Methanol dissolved the predominant PSS and EG from the surface. The redshift of the peak representing the symmetrical vibration of Cα═Cβ in the Raman spectrum confirmed the conformation of benzoid structure to quinoid structure after the surface treatment. This conformation was attributed to the extension of the conjugation length and the reduction of the energy barrier within the PEDOT chain. This resulted in the improved conductivity and charge hopping of the PEDOT:PSS, which was also proven using density functional theory (DFT) calculations. Reducing the insulating and acidic PSS improved the electroluminescence performance and extended the operational lifetime of the QLEDs. The tris(dimethylamino)phosphine-based InP QLEDs exhibited an external quantum efficiency (EQE) of 6.4%, that value is comparable to those of tris(trimethylsilyl)phosphine-based QLEDs, and operational lifetime (T50) of 125.6 h.

Reaction modelling of hydrogen evolution on nickel phosphide catalysts: density functional investigation

Nickel phosphides (NixPy), particularly Ni2P, are promising catalysts for the acidic hydrogen evolution reaction (HER). Using density functional theory (DFT), we model HER at the potential of zero charge (PZC), incorporating solvation effects via an explicit water cluster and implicit surrounding solvent. Comparing the Volmer, Tafel, and Heyrovsky steps under saturated hydrogen coverage on Ni2P(0001) terminations, we find that the Ni3P2 (pristine) surface termination prefers the Volmer-Volmer-Tafel (VVT) pathway with activation energy (Ea) of 0.57 eV. Conversely, the Ni3P2 + 4P (reconstructed) surface favors the Volmer-Heyrovsky (VH) pathway with Ea = 0.60 eV. For the pristine surface termination, the differential gas-phase hydrogen adsorption free energies (ΔGdiff) correlate with the Volmer and Tafel step reaction energies, and a linear Bell-Evans-Polanyi relationship for the calculated activation and reaction energies validates the usefulness of the ΔGdiff descriptor for the Volmer step under PZC conditions. Nickel atoms play a crucial role in H2 production on both pristine and reconstructed surfaces, suggesting that modifications of the Ni sites can be used for catalyst design. Our findings highlight the importance of considering surface reconstruction and solvation effects on the HER catalytic performance.

Rigidly Axial O Coordination-Induced Spin Polarization on Single Ni-N4-C Site by MXene Coupling for Boosting Electrochemical CO2 Reduction to CO

Regulating the spin states in transition-metal (TM)-based single-atom catalysts (SACs), such as the TM-Nx-C configurations, is crucial for improving the catalytic activity. However, the role of spin in single Ni atoms facilitating the electrochemical CO2 reduction reaction (CO2RR) has been largely overlooked. Using first-principles simulations, we investigated the electrocatalytic performance of Ni-N4-C SACs vertically stacked on the O-terminated MXene nanosheets for the CO2RR. The terminated O atoms on MXene axially interact with the Ni atom due to significant charge transfer between them. Unlike the pure Ni-N4 site, which lacks spin polarization, the newly formed Ni-N4O configuration breaks the spin degeneracy of Ni d orbitals, dramatically lifting the energy level of spin-down d orbitals relative to that of spin-up d orbitals. As a result, the d electrons of Ni in the two spin channels are rearranged, leading to large net spin moments of 1.4 μB. Compared to the Ni-N4 site, the partially filled minority-spin dz2 orbitals of Ni on Ni-N4O weaken the occupied d-π* orbitals between Ni and *COOH, significantly stabilizing the key intermediate. The detailed reaction mechanisms and energetics show that four MXenes, namely, Hf3C2, Zr3C2, Hf2C, and Zr2C, can induce a large spin on the Ni site, thereby improving catalytic activity for CO2 reduction to CO, with a lower onset potential of about -0.75 V vs SHE compared to pure Ni SACs (-1.17 V) according to the potential-constant model with an explicit solvent environment.

Integrating Active Learning and DFT for Fast-Tracking Single-Atom Alloy Catalysts in CO2-to-Fuel Conversion

Electrocatalytic carbon dioxide reduction (CO2RR) technology enables the conversion of excessive CO2 into high-value fuels and chemicals, thereby mitigating atmospheric CO2 concentrations and addressing energy scarcity. Single-atom alloys (SAAs) possess the potential to enhance the CO2RR performance by full utilization of atoms and breaking linear scaling relationships. However, quickly screening high-performance metal portfolios of SAAs remains a formidable challenge. In this study, we proposed an active learning (AL) framework to screen high-performance catalysts for CO2RR to yield fuels such as CH4 and CH3OH. After four rounds of AL iterations, the ML model attained optimal prediction performance with the test set R2 of approximately 0.94 and successful prediction was achieved for the binding free energy of *CHO, *COH, *CO, and *H on 380 catalyst surfaces with an accuracy within 0.20 eV. Subsequent analysis of the SAA catalysts' activity, selectivity, and stability led to the discovery of eight previously unexplored SAA catalysts for CO2RR. Notably, the surface activity of Ti@Cu(100), Au@Pt(100), and Ag@Pt(100) shines prominently. Utilizing DFT calculations, we elucidated the complete reaction pathway of the CO2RR on the surfaces of these catalysts, confirming their high catalytic activity with limiting potentials of -0.11, -0.34, and -0.46 eV, respectively, which are significantly lower than those of pure Cu catalysts. The results showcase the exceptional predictive prowess of AL, providing a valuable reference for the design of CO2RR catalysts.

Nonbonding Metal-Metal Interaction in Metal-Nitrogen-Carbon Single-Atom Catalysts Boosts CO Electroreduction

Metal-nitrogen-carbon single-atom catalysts (SACs) have recently emerged as selective electrocatalysts for the reduction of CO2 to CO, but their ability to further electroreduce CO is poor. Here, based on constant-potential density functional theory simulations, we predict that Co-N-M (M = Fe, Co) SACs with nonbonding metallic centers bridged by a common nitrogen atom can catalyze four-electron reduction of CO to methanol at an ultralow overpotential of 220-310 mV. We show that the metal atoms in the SACs are terminated by H species which prevent the formation of σ bonding between CO and the metal atoms. Thanks to the nonbonding electrostatic repulsion between Co and its adjacent M atom, the Co dxz band is broadened and shifted toward the Fermi level, leading to enhanced dxz - 2π* interaction that gives rise to stable CO adsorption and promotes its active and selective reduction. This work offers an alternative strategy to tackle the challenge of CO electroreduction on SACs and highlights the role of nonbonding metal-metal interactions in modulating adsorption properties of SACs.

Unveiling the role of dopants in boosting CuS supercapacitor performance: insights from first-principles calculations

Transition metal sulfides have become famous in high energy density supercapacitor materials owing to their rich redox and high conductivity. While their development has achieved a breakthrough in terms of capacitance, there is little knowledge from the theoretical perspective on how dopants play a role in enhancing their capacitances. In this work, pseudocapacitance and quantum capacitance were evaluated through first-principles calculation to describe their role in transition metal sulfide, which here is represented by copper sulfide (CuS). The resulting quantum capacitance (CQ) was calculated in both the bulk and surface of CuS to determine which structure has a greater effect on the capacitance of the system. It was observed that the dopant increased CQ in the bulk system, which is different from the CQ of surface structures. Meanwhile, K+ ions were introduced on the surface structure to calculate transfer charge and work function shift, thus determining pseudocapacitance. All dopant types were able to increase the pseudocapacitance value, with Fe doping showing the highest capacitance of 111 F g-1, which is higher than that of the pristine structure (47 F g-1). The role of the dopant is discussed in detail in this work. Our results suggest that the increased capacitance of doped TMS materials was originated not only from the geometrical perspective but also from the higher pseudocapacitance value. Quantum capacitance, alternatively, could also contribute to the system when the dopant occurs in the bulk rather than only in the surface structure. This work may open a new perspective on how dopants play a role in increasing supercapacitor performance.

Effect of Strain Engineering on the Spin State of the Ni-N4/C Single-Atom Catalyst and Its Consequence in Electrocatalysis

Strain engineering is an effective strategy to improve the activity of catalysts, especially for flexible carbon-based materials. Nitrogen-coordinated single atomic metals on a carbon skeleton (M-Nx/C) are of interest in catalytic electroreduction reactions due to their high activity and atomic utilization. However, the effect of strain on the structure-activity relationship between the electrochemical activity and the electronic and geometric structures of Ni-Nx/C remains unclear. Here, we found that by applying tensile strain on the Ni-N4/C, the spin state of the single atom can be changed from a low-spin to a high-spin state. Moreover, the energy gap between the highest occupied d orbital of Ni and the lowest unoccupied molecular orbital of the adsorbed species narrowed. With an increasing strain rate, the catalytic activity of O2 and CO2 electroreduction can be improved. Especially for the 2e- O2 reduction, the implicit solvent model, constant-potential method, and microkinetic model were used to verify the positive effect of suitable stretching on the catalytic activity from thermodynamic and kinetic viewpoints. This work can reveal the relationship between strain, spin state, and the catalytic activity of Ni-Nx/C.

Machine learning-assisted dual-atom sites design with interpretable descriptors unifying electrocatalytic reactions

Low-cost, efficient catalyst high-throughput screening is crucial for future renewable energy technology. Interpretable machine learning is a powerful method for accelerating catalyst design by extracting physical meaning but faces huge challenges. This paper describes an interpretable descriptor model to unify activity and selectivity prediction for multiple electrocatalytic reactions (i.e., O2/CO2/N2 reduction and O2 evolution reactions), utilizing only easily accessible intrinsic properties. This descriptor, named ARSC, successfully decouples the atomic property (A), reactant (R), synergistic (S), and coordination effects (C) on the d-band shape of dual-atom sites, which is built upon our developed physically meaningful feature engineering and feature selection/sparsification (PFESS) method. Driven by this descriptor, we can rapidly locate optimal catalysts for various products instead of over 50,000 density functional theory calculations. The model's universality has been validated by abundant reported works and subsequent experiments, where Co-Co/Ir-Qv3 are identified as optimal bifunctional oxygen reduction and evolution electrocatalysts. This work opens the avenue for intelligent catalyst design in high-dimensional systems linked with physical insights.

Electrochemical Stability of MXenes in Water Based on Constant Potential AIMD Simulations

MXene has been recently explored as promising electrocatalytic materials to accelerate the electrocatalytic process for hydrogen evolution, but their dynamic stability under electrochemical conditions remains elusive. Here we performed first-principle ab initio molecular dynamics calculations to reveal the electrochemical stability of Ti2CTx MXene in different aqueous environments. The results revealed the high vulnerability of the pure and vacancy-defected Ti2CO2 MXene towards water attack, leading to surface oxidation of MXene under neutral electrochemical condition that formed adsorbed oxygen species to Ti and dissociated proton in solution. The surface oxidation of Ti2CO2 could be prevented in the acid condition or in the neutral condition under the negative potential. Differently, the fully F- or OH-functionalized Ti2CF2 and Ti2C(OH)2 as well as the mixed functionalized Ti2C(O0.5OH0.5)2 and Ti2CO1.12F0.88 are highly stable under various electrochemical conditions, which can effectively prevent close contact between water and surface Ti atoms via electronic repulsion or steric hindrance. These findings provide atomic level understanding of the aqueous stability of MXene and provide useful strategies to prevent degradation and achieve highly stable MXenes.

Surface Termination (-O, -F or -OH) and Metal Doping on the HER Activity of Mo2CTx MXene

Two-dimensional MXenes have recently garnered significant attention as electrocatalytic materials for hydrogen evolution reaction (HER). However, previous theoretical studies mainly focused on the effect of pure functional groups while neglecting hybrid functional groups that are commonly observed in experiments. Herein, we investigated the hybrid functionalized Mo2CTx MXene (T=-O, -F or -OH) to probe the HER properties. In binary O/F co-functionalization, the presence of F groups would attenuate the H adsorption and lead to the enhanced HER activity than the fully O-terminated Mo2CO2. However, the surface HER activity of ternary O/F/OH functionalized Mo2CTx is not satisfactory owing to the relatively weak H adsorption capacity. To further enhance the catalytic activity, modification was performed by introducing another metal element into its lattice structure. The doped metal (Fe, Co, Ni, Cu) exhibits reduced charge transfer to O compared to Mo atoms, leading to enhanced H adsorption and improved overall activity. The synergistic effect of hybrid functionalization and TM modification provides useful guidance for achieving feasible Mo2CTx candidates with high HER performance, which can be applied to the electrocatalytic applications of other MXenes.

H-assisted CO2 dissociation on PdnPt(4-n)/In2O3 catalysts: a density functional theory study

CO2 hydrogenation into valuable chemical compounds can effectively address the issues of greenhouse gas emissions and energy scarcity. The activation and dissociation processes of CO2 are crucial for its reduction reactions, but the effects of *H adatoms on the C-O cleavage are still confusing. This study investigates the H-assisted CO2 dissociation pathways on the PdnPt(4-n)/In2O3 (n = 0-4) catalysts via DFT calculation. Initially, the adsorption properties of *H2, *COOH, and *HCOO species are calculated. Then, two H-assisted CO2 dissociation channels, i.e., *CO2 + *H → *COOH → *CO + *OH and *CO2 + *H → *HCOO → *CHO + *O, are studied. Results show that Pt and Pd promote the CO2 hydrogenation and C-O bond cleavage reactions, respectively. In comparison to CO2 direct dissociation, the COOH-mediated and HCOO-mediated channels facilitate and impede the C-O bond cleavage, respectively. Overall, the Pd3Pt/In2O3 constituent is suggested for the H-assisted CO2 dissociation reaction. The electronic effects of the PdnPt(4-n) bimetals adjust the stabilities of the intermediates and barriers of the elementary steps, and the interactions between PdnPt(4-n) and In2O3 provide extra sites for the adsorbates and reaction steps. This study reveals the effects of *H on the C-O bond dissociation processes and provides useful insight into designing PdPt/In2O3 catalysts for CO2 hydrogenation reactions.

Stabilized Cu0 -Cu1+ dual sites in a cyanamide framework for selective CO2 electroreduction to ethylene

Electrochemical reduction of carbon dioxide to produce high-value ethylene is often limited by poor selectivity and yield of multi-carbon products. To address this, we propose a cyanamide-coordinated isolated copper framework with both metallic copper (Cu0) and charged copper (Cu1+) sites as an efficient electrocatalyst for the reduction of carbon dioxide to ethylene. Our operando electrochemical characterizations and theoretical calculations reveal that copper atoms in the Cuδ+NCN complex enhance carbon dioxide activation by improving surface carbon monoxide adsorption, while delocalized electrons around copper sites facilitate carbon-carbon coupling by reducing the Gibbs free energy for *CHC formation. This leads to high selectivity for ethylene production. The Cuδ+NCN catalyst achieves 77.7% selectivity for carbon dioxide to ethylene conversion at a partial current density of 400 milliamperes per square centimeter and demonstrates long-term stability over 80 hours in membrane electrode assembly-based electrolysers. This study provides a strategic approach for designing catalysts for the electrosynthesis of value-added chemicals from carbon dioxide.

Grand Canonical Ensemble Approaches in CP2K for Modeling Electrochemistry at Constant Electrode Potentials

In electrochemical experiments, the number of electrons of the electrode immersed in the electrolyte is usually variable. Additionally, the numbers of adsorbed substances on the surface of the electrode, the solvent molecules, and counter charge ions in the near-surface region can also vary. Treating electrochemical solid-liquid interfaces with the typical fixed electron number density functional theory (DFT) approach tends to be a challenge. This can be addressed by using grand canonical ensemble approaches. We present the implementation of two grand canonical ensemble approaches in the open-source computational chemistry software CP2K that go beyond the existing canonical ensemble paradigm. The first approach is based on implicit solvent models and explicit atomistic solute (electrode with/without adsorbed species) models, and includes two recent developments: (a) grand canonical self-consistent field (GC-SCF) method (J. Chem. Phys. 2017, 146, 114104) allowing the electron number of the system to fluctuate naturally and accordingly with the experimental electrode potential, (b) planar counter charge (J. Chem. Phys. 2019, 150, 041722, Phys. Rev. B 2003, 68, 245416) salt model completely screening the net charge of the electrode model. In contrast with previous studies, in our implementation, the work function (WF) (absolute electrode potential if the potential drop at the electrolyte-vacuum interface is omitted) is the constrained quantity during an SCF optimization instead of the Fermi energy. The chemical potential of electrons (negative WF) is a natural variable of the grand potential in the GC ensemble of electronic states, and this method can easily achieve stable SCF convergence and obtain an electronic structure that precisely corresponds to a user-specified WF. The second approach referred to as the GC DFT molecular dynamics (DFT-MD) simulation scheme (Phys. Rev. Lett. 2002, 88, 213002, J. Chem. Phys. 2005, 122, 234505, J. Am. Chem. Soc. 2004, 126 (12), 3928-3938) is based on fully explicit modeling the solvent molecules and the ions and is used to calculate the electron chemical potential corresponding to an equilibrium electrochemical half-reaction (M(n+m)+ + ne- ⇌ Mm+) which involves DFT-MD, by allowing the number of electrons to vary during the DFT-MD simulation process. This opens the way for forefront electrochemical calculations in CP2K for a broad range of systems.

Plasmonic Pd-Sb nanosheets for photothermal CH4 conversion to HCHO and therapy

Photothermal catalysis effectively increases catalytic activity by using the photothermal effect of metal nanomaterials; however, the combination of strong light absorption and high catalytic performance remains a challenge. Here, we demonstrate hexagonal ~5-nanometer-thick palladium antimony (chemical formula as Pd8Sb3) nanosheets (NSs) that exhibit strong light absorption within full spectral and localized surface plasmon resonance (LSPR) effects in the visible region. Such LSPR features lead to strong photothermal effects, and Pd8Sb3 NSs aqueous dispersion enables enhanced photothermal methane (CH4) conversion to formaldehyde (HCHO) under full-spectrum light irradiation at 1.7 watts per square centimeter, leading to selectivity of ~98.7%, productivity of ~665 millimoles per gram of catalyst, ~700 times higher than that of Pd NSs. Mechanism investigations suggest that different radicals were generated on Pd8Sb3 (·OH) and Pd NSs (·O2-), where Pd8Sb3 NSs displays stronger adsorption strength to CH4 and facilitates CH4 oxidation to HCHO. Besides, the strong light absorption ability of Pd8Sb3 NSs enables photothermal therapy for breast cancer.

Revisiting the activity origin of the PtAu24(SR)18 nanocluster for enhanced electrocatalytic hydrogen evolution by combining first-principles simulations with the experimental in situ FTIR technique

Thiolate-protected metal nanoclusters (NCs) have been widely used in various electrocatalytic reactions, yet the dynamic evolution of metal NCs during electrocatalysis has been rarely explored and the activity origin remains largely ambiguous. Herein, using a PtAu24(SCH3)18 NC as a prototype model, we combined advanced first-principles calculations and attenuated total reflection surface-enhanced infrared spectroscopy (ATR-SEIRAS) to re-examine its active site and reaction dynamics in the hydrogen evolution reaction (HER). It has been previously assumed that the central Pt is the only catalytic center. However, differently, we observed the spontaneous desorption of thiolate ligands under moderate potential, and the dethiolated PtAu24 exhibits excellent HER activity, which is contributed not only by the central Pt atom but also by the exposed bridged Au sites. Particularly, the exposed Au exhibits high activity even comparable to Pt, and the synergistic effect between them makes dethiolated PtAu24 an extraordinary HER electrocatalyst, even surpassing the commercial Pt/C catalyst. Our predictions are further verified by electrochemical activation experiments and in situ FTIR (ATR-SEIRAS) characterization, where evident adsorption of Au-H* and Pt-H* bonds is monitored. This work detected, for the first time, the Au-S interfacial dynamics of the PtAu24 nanocluster in electrocatalytic processes, and quantitatively evaluated the essential catalytic role of the exposed Au sites that has been largely overlooked in previous studies.

Long-range interactions driving neighboring Fe-N4 sites in Fenton-like reactions for sustainable water decontamination

Actualizing efficient and sustainable environmental catalysis is essential in global water pollution control. The single-atom Fenton-like process, as a promising technique, suffers from reducing potential environmental impacts of single-atom catalysts (SACs) synthesis and modulating functionalized species beyond the first coordination shell. Herein, we devised a high-performance SAC possessing impressive Fenton-like reactivity and extended stability by constructing abundant intrinsic topological defects within carbon planes anchored with Fe-N4 sites. Coupling atomic Fe-N4 moieties and adjacent intrinsic defects provides potent synergistic interaction. Density functional theory calculations reveal that the intrinsic defects optimize the d-band electronic structure of neighboring Fe centers through long-range interactions, consequently boosting the intrinsic activity of Fe-N4 sites. Life cycle assessment and long-term steady operation at the device level indicate promising industrial-scale treatment capability for actual wastewater. This work emphasizes the feasibility of synergistic defect engineering for refining single-atom Fenton-like chemistry and inspires rational materials design toward sustainable environmental remediation.

Solar energy powered electrochemical reduction of CO2 on In2O3 nanosheets with high energy conversion efficiency at a large current density

Electrochemical CO2 reduction (ECO2R) to value-added chemicals offers a promising approach to both mitigate CO2 emission and facilitate renewable energy conversion. We demonstrate a solar energy powered ECO2R system operating at a relatively large current density (57 mA cm-2) using In2O3 nanosheets (NSs) as the cathode and a commercial perovskite solar cell as the electricity generator, which achieves the high solar to formate energy conversion efficiency of 6.6 %. The significantly enhanced operative current density with a fair solar energy conversion efficiency on In2O3 NSs can be ascribed to their high activity and selectivity for formate production, as well as the fast kinetics for ECO2R. The Faradic efficiencies (FEs) of formate In2O3 NSs are all above 93 %, with the partial current density of formate ranging from 2.3 to 342 mA cm-2 in a gas diffusion flow cell, which is among the widest for formate production on In-based catalysts. In-situ Raman spectroscopy and density functional theory simulations reveal that the exceptional performances of formate production on In2O3 NSs originates from the presence of abundant low coordinated edge sites, which effectively promote the selective adsorption of *OCHO while inhibiting *H adsorption.

Disentangling heterogeneous thermocatalytic formic acid dehydrogenation from an electrochemical perspective

Heterogeneous thermocatalysis of formic acid dehydrogenation by metals in solution is of great importance for chemical storage and production of hydrogen. Insightful understanding of the complicated formic acid dehydrogenation kinetics at the metal-solution interface is challenging and yet essential for the design of efficient heterogeneous formic acid dehydrogenation systems. In this work, formic acid dehydrogenation kinetics is initially studied from a perspective of electrochemistry by decoupling this reaction on Pd catalyst into two short-circuit half reactions, formic acid oxidation reaction and hydrogen evolution reaction and manipulating the electrical double layer impact from the solution side. The pH-dependences of formic acid dehydrogenation kinetics and the associated cation effect are attributed to the induced change of electric double layer structure and potential by means of electrochemical measurements involving kinetic isotope effect, in situ infrared spectroscopy as well as grand canonical quantum mechanics calculations. This work showcases how kinetic puzzles on some important heterogeneous catalytic reactions can be tackled by electrochemical theories and methodologies.

Best practices of modeling complex materials in electrocatalysis, exemplified by oxygen evolution reaction on pentlandites

Pentlandites are natural ores with structural properties comparable to that of [FeNi] hydrogenases. While this class of transition-metal sulfide materials - (Fe,Ni)9S8 - with a variable Fe : Ni ratio has been proven to be an active electrode material for the hydrogen evolution reaction, it is also discussed as electrocatalyst for the alkaline oxygen evolution reaction (OER), corresponding to the bottleneck of anion exchange membrane electrolyzers for green hydrogen production. Despite the experimental evidence for the use of (Fe,Ni)9S8 as an OER catalyst, a detailed investigation of the elementary reaction steps, including consideration of adsorbate coverages and limiting steps under anodic polarizing conditions, is still missing. We address this gap in the present manuscript by gaining atomistic insights into the OER on an Fe4.5Ni4.5S8(111) surface through density functional theory calculations combined with a descriptor-based analysis. We use this system to introduce best practices for modeling this rather complex material by pointing out hidden pitfalls that can arise when using the popular computational hydrogen electrode approach to describe electrocatalytic processes at the electrified solid/liquid interface for energy conversion and storage.

Diffusion of ferrocene through vanadyl phosphate by density functional theory

Here, we employed the nudged elastic band (NEB) method to simulate the diffusion of ferrocene through vanadyl phosphate (VOPO4), with a focus on understanding the diffusion pathways arising from the complex structure of ferrocene. We systematically evaluated a total of 36 potential diffusion paths, categorizing them into three groups based on their directional orientation: 15 paths between V sites along the [110] direction, 15 paths from V to P sites along the [100] direction, and 6 paths between P sites also along the [110] direction. Our analysis revealed that the energy barriers for diffusion along the [110] direction typically ranged between 0.25 and 0.35 eV, which are notably higher than those observed for pathways along the [100] direction, where the energy barriers ranged from 0.11 to 0.20 eV. To further elucidate the complex deformation of ferrocene during diffusion, we established four key measures to characterize the structural conformation: the angle of the axis of the ferrocene molecule relative to the [010] direction within the (001) plane, the dihedral angle between the two cyclopentadienyl rings, the orientation angle of the -CH bonds with respect to the [001] direction, and the angle between two -CH bonds from the two cyclopentadienyl rings.

Revealing Atomic Structure of Hybrid Octacalcium-Phosphate Derivative

Octacalcium phosphate (OCP), which is one of the bioactive calcium phosphates, can incorporate various organic molecules in its crystal lattice, forming the organic-inorganic hybrid derivatives. However, detailed atomic arrangements of OCP hybridized with organic molecules such as dicarboxylate are still unknown, although many years have passed since the first discovery of the materials systems. In the present study, some black-box optimization methods combined with first-principles calculations were used to theoretically identify the most stable atomic structure of the OCP with the incorporation of malonate ions as a typical case study. The results showed that the calculated interplanar spacing on the (100) plane of the most stable structure agrees well with experimental data, by taking account of implicit solvent of aqueous solution. An underlying mechanism that realizes the bridging feature of the incorporated malonate ions between the apatitic layers is also discussed. The present methodology can pave the way to accurately explore reliable atomic structures of such complicated organic-inorganic hybrid biomaterials with high structural degrees of freedom.

Atomic-Dispersed Cu Catalysts for Electrochemical Nitrate Reduction: Coordination Engineering and Fundamental Insights

The development of Cu-based atomic dispersed catalysts with tailored coordination environments represents a significant step forward in enhancing the electrocatalytic reduction of nitrate to ammonia. By precisely modulating the electronic structures of Cu active centers, the binding strength of the *NO3 intermediates is successfully tuned, thereby substantially improving the catalytic activity toward electrochemical nitrate reduction reaction (eNO3RR). This study reveals that the N4-coordinated Cu single-atom catalyst (Cu-SAC) exhibits superior performance due to its robust interaction with coordinating atoms. Notably, this optimized catalyst achieves a low limiting potential of -0.38 V, while the dual-atom system further reduces this value to -0.32 V, demonstrating exceptional activity. Detailed electronic structure analysis, including the examination of d-band centers, Bader charges, and projected density of states (PDOS), provides a comprehensive understanding of the origin of this high activity. Specifically, the high and concentrated density of states near the Fermi level plays a crucial role in facilitating the electrocatalytic nitrate reduction process. This work not only offers crucial insights into the underlying mechanisms of eNO3RR but also provides valuable guidelines for the rational design of highly efficient electrocatalysts for this important reaction.

Doping-Induced Enhancement of Hydrogen Evolution at MoS2 Electrodes

Rate theory and DFT calculations of hydrogen evolution reaction (HER) on MoS2 with Co, Ni and Pt impurities show the significance of dihydrogen (H2*) complex where both hydrogen atoms are interacting with the surface. Stabilization of such a complex affects the competing Volmer-Heyrovsky (direct H2 release) and Volmer-Tafel (H2* intermediate) pathways. The resulting evolution proceeds with a very small overpotential for all dopants ( η ${\eta }$ =0.1 to 0.2 V) at 25 % edge substitution, significantly reduced from the already low η ${\eta }$ =0.27 V for the undoped edge. At full edge substitution, Co-MoS2 remains highly active ( η ${\eta }$ =0.18 V) while Ni- and Pt-MoS2 are deactivated ( η ${\eta }$ =0.4 to 0.5 V) due to unfavorable interaction with H2*. Instead of the single S-vacancy, the site of intrinsic activity in the basal plane was found to be the undercoordinated central Mo-atom in threefold S-vacancy configurations, enabling hydrogen evolution with η ${\eta }$ =0.52 V via a H2* intermediate. The impurity atoms interact favorably with the intrinsic sulfur vacancies on the basal plane, stabilizing but simultaneously deactivating the triple vacancy configuration. The calculated shifts in overpotential are consistent with reported measurements, and the dependence on doping level may explain variations in experimental observations.

Theoretical and Experimental Studies of Molecular Interactions between Engineered Graphene and Phosphate Ions for Graphene-Based Phosphate Sensing

Fundamental understanding of the interactions of nanoscale materials with molecules of interest is essential for the development of electronic devices, such as sensors. In particular, structures and molecular interaction properties of engineered graphenes are still largely unexplored, despite these materials' great potential to be used as molecular sensors. As an example of end user application, the detection of phosphorus in the form of phosphate in a soil environment is important for soil fertility and plant growth. However, due to the lack of an affordable technology, it is currently hard to measure the amount of phosphate directly in the soil; therefore, suitable sensor technologies need to be developed for phosphate sensors. In this work, pristine graphene and several modified graphene materials (oxygenated graphene, graphene with vacancies, and curved graphene) were studied as candidates for phosphate sensor materials using density functional theory (DFT) calculations. Our calculations showed that both pristine graphene and functionalized graphene were able to adsorb phosphate species strongly. In addition, these graphene nanomaterials showed selectivity of adsorption of phosphate with respect to nitrate, with stronger adsorption energies for phosphate. Furthermore, our calculations showed significant changes in electrical conductivities of pristine graphene and functionalized graphenes after phosphate species adsorption, in particular, on graphene with oxygen (hydroxyl and epoxide) functional groups. Experimental measurements of electrical resistivity of graphene before and after adsorption of dihydrogen phosphate showed an increase in resistivity upon adsorption of phosphate, consistent with the theoretical predictions. Our results recommend graphene and functionalized graphene-based nanomaterials as good candidates for the development of phosphate sensors.

Integrated electrocatalytic synthesis of ammonium nitrate from dilute NO gas on metal organic frameworks-modified gas diffusion electrodes

The electrocatalytic conversion of NO offers a promising technology for not only removing the air pollutant but also synthesizing valuable chemicals. We design an integrated-electrocatalysis cell featuring metal organic framework (MOF)-modified gas diffusion electrodes for simultaneous capture of NO and generation of NH4NO3 under low-concentration NO flow conditions. Using 2% NO gas, the modified cathode exhibits a higher NH4+ yield and Faradaic efficiency than an unmodified cathode. Notably, the modified cathode shows a twofold increase in NH4+ production with 20 ppm NO gas supply. Theoretical calculations predict favorable transfer of adsorbed NO from the adsorption layer to the catalyst layer, which is experimentally confirmed by enhanced NO mass transfer from gas to electrolyte across the modified electrode. The adsorption layer-modified anode also exhibits a higher NO3- yield for NO electro-oxidation compared to the unmodified electrode under low NO concentration flow. Among various integrated-cell configurations, a single-chamber setup produces a higher NH4NO3 yield than a double-chamber setup. Furthermore, a higher NO utilization efficiency is obtained with a single-gasline operation mode, where the NO-containing gas flows sequentially from the cathode to the anode.

Origin of copper dissolution under electrocatalytic reduction conditions involving amines

Cu dissolution has been identified as the dominant process that causes cathode degradation and losses even under cathodic conditions involving methylamine. Despite extensive experimental research, our fundamental and theoretical understanding of the atomic-scale mechanism for Cu dissolution under electrochemical conditions, eventually coupled with surface restructuring processes, is limited. Here, driven by the observation that the working Cu electrode is corroded using mixtures of acetone and methylamine even under reductive potential conditions (-0.75 V vs. RHE), we employed Grand Canonical density functional theory to understand this dynamic process under potential from a microscopic perspective. We show that amine ligands in solution directly chemisorb on the electrode, coordinate with the metal center, and drive the rearrangement of the copper surface by extracting Cu as adatoms in low coordination positions, where other amine ligands can coordinate and stabilize a surface copper-ligand complex, finally forming a detached Cu-amine cationic complex in solution, even under negative potential conditions. Calculations predict that dissolution would occur for a potential of -1.1 V vs. RHE or above. Our work provides a fundamental understanding of Cu dissolution facilitated by surface restructuring in amine solutions under electroreduction conditions, which is required for the rational design of durable Cu-based cathodes for electrochemical amination or other amine involving reduction processes.

High vacancy formation energy boosts the stability of structurally ordered PtMg in hydrogen fuel cells

Alloys of platinum with alkaline earth metals promise to be active and highly stable for fuel cell applications, yet their synthesis in nanoparticles remains a challenge due to their high negative reduction potentials. Herein, we report a strategy that overcomes this challenge by preparing platinum-magnesium (PtMg) alloy nanoparticles in the solution phase. The PtMg nanoparticles exhibit a distinctive structure with a structurally ordered intermetallic core and a Pt-rich shell. The PtMg/C as a cathode catalyst in a hydrogen-oxygen fuel cell exhibits a mass activity of 0.50 A mgPt-1 at 0.9 V with a marginal decrease to 0.48 A mgPt-1 after 30,000 cycles, exceeding the US Department of Energy 2025 beginning-of-life and end-of-life mass activity targets, respectively. Theoretical studies show that the activity stems from a combination of ligand and strain effects between the intermetallic core and the Pt-rich shell, while the stability originates from the high vacancy formation energy of Mg in the alloy.

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.

Water splitting at imine-linked covalent organic frameworks

Covalent organic frameworks (COFs) are a promising class of metal-free catalysts, offering a high structural and functional variety. Here, we systematically study imine-linked COFs with donor (D) and acceptor (A) groups using density functional theory (DFT). Using water splitting as a model reaction, we analyze the effects of protonation of the catalyst, the orientation of the imine linkage leading to different constitutional isomers, and solvation. In agreement with experimental results, we show that protonation decreases the band gap. In addition, COFs in which the donor is closer to the nitrogen atom of the imine group (DNCA) have lower band gaps than those in which the donor is closer to the carbon atom (DCNA). Three different D/A COFs are compared in this work, for which energies for the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER) and corresponding electrochemical overpotentials are computed. We show that reaction energies are very similar for DCNA and DNCA COFs. The differences in hydrogen evolution rates between the constitutional isomers observed experimentally in (photocatalytic) HER (Yang et al., Nat. Commun., 2022, 13, 6317), are proposed to be at least in part a consequence of differences in charge distribution.

Achieving Exceptional Volumetric Desalination Capacity Using Compact MoS2 Nanolaminates

Capacitive deionization (CDI) has emerged as a promising technology for freshwater recovery from low-salinity brackish water. It is still inapplicable in specific scenarios (e.g., households, islands, or offshore platforms) due to too low volumetric adsorption capacities. In this study, a high-density semi-metallic molybdenum disulfide (1T'-MoS2) electrode with compact architecture obtained by restacking of exfoliated nanosheets, which achieve high capacitance up to ≈277.5 F cm-3 under an ultrahigh scan rate of 1000 mV s-1 with a lower charge-transfer resistance and nearly tenfold higher electrochemical active surface area than the 2H-MoS2 electrode, is reported. Furthermore, 1T'-MoS2 electrode demonstrates exceptional volumetric desalination capacity of 65.1 mgNaCl cm-3 in CDI experiments. Ex situ X-ray diffraction (XRD) reveal that the cation storage mechanism with the dynamic expansion of 1T'-MoS2 interlayer to accommodate cations such as Na+, K+, Ca2+, and Mg2+, which in turn enhances the capacity. Theoretical analysis unveils that 1T' phase is thermodynamically preferable over 2H phase, the ion hydration and channel confinement also play critical role in enhancing ion adsorption. Overall, this work provides a new method to design compact 2D-layered nanolaminates with high-volumetric performance for CDI desalination.

Correcting implicit solvation at metal/water interfaces through the incorporation of competitive water adsorption

Conventional continuum solvation models are ubiquitous in computational catalysis, including for describing metal/water interfaces, which are relevant to both solution-phase heterogeneous catalysis and electrocatalysis. Nonetheless, we find that such continuum models qualitatively fail to describe both the adsorption free energy and conformational preference for many organic molecules at such interfaces, largely due to the failure of continuum models to incorporate the role of competitive water adsorption. We develop a simple phenomenological model that accounts for competitive water adsorption and show that the model, when used in conjunction with continuum solvation, provides a dramatic improvement in the description of both adsorption and conformational preference. The model is also extended to additionally incorporate the influence of applied potential at the electrode surface, thus facilitating computationally efficient applications to scenarios including 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.

FeN4-Embedded Graphene as a Highly Sensitive and Selective Single-Atom Sensor for Reaction Intermediates of Electrochemical CO2 Reduction

Exploring effective ways to detect intermediates during the electrochemical CO2 reduction reaction (CO2RR) process is pivotal for understanding reaction pathways and underlying mechanisms. Recently, two-dimensional FeN4-embedded graphene has received increasing attention as a promising catalyst for CO2RR. Here, by means of density functional theory computations combined with the non-equilibrium Green's function (NEGF) method, we proposed a detection device to evaluate the performance of FeN4-embedded graphene in intermediates detection during the CO2RR process. Our results reveal that the four key intermediates, including *COOH, *OCHO, *CHO, and *COH, can be chemisorbed on FeN4-embedded graphene with high adsorption energies and appropriate charge transfer. The computed current-voltage (I-V) characteristics and transmission spectra suggest that the adsorption of these intermediates induces significant type-dependent changes in currents and transmission coefficients of FeN4-embedded graphene. Remarkably, the FeN4-embedded graphene is more sensitive to *COOH and *COH than to *OCHO and *CHO within the entire bias window. Consequently, our theoretical study indicates that the FeN4-embedded graphene can effectively detect the key intermediates during the CO2RR process, providing a practical scheme for identifying catalytic reaction pathways and elucidating underlying reaction mechanisms.

Parameter-Fitting-Free Continuum Modeling of Electric Double Layer in Aqueous Electrolyte

Electric double layers (EDLs) play fundamental roles in various electrochemical processes. Despite the extensive history of EDL modeling, there remain challenges in the accurate prediction of its structure without expensive computation. Herein, we propose a predictive multiscale continuum model of EDL that eliminates the need for parameter fitting. This model computes the distribution of the electrostatic potential, electron density, and species' concentrations by taking the extremum of the total grand potential of the system. The grand potential includes the microscopic interactions that are newly introduced in this work: polarization of solvation shells, electrostatic interaction in parallel plane toward the electrode, and ion-size-dependent entropy. The parameters that identify the electrode and electrolyte materials are obtained from independent experiments in the literature. The model reproduces the trends in the experimental differential capacitance with multiple electrode and nonadsorbing electrolyte materials (Ag(110) in NaF, Ag(110) in NaClO4, and Hg in NaF), which verifies the accuracy and predictiveness of the model and rationalizes the observed values to be due to changes in electron stability. However, our calculation on Pt(111) in KClO4 suggests the need for the incorporation of electrode/ion-specific interactions. Sensitivity analyses confirmed that effective ion radius, ion valence, the electrode's Wigner-Seitz radius, and the bulk modulus of the electrode are significant material properties that control the EDL structure. Overall, the model framework and findings provide insights into EDL structures and predictive capability at low computational cost.

Electroreduction of Captured CO2 on Silver Catalysts: Influence of the Capture Agent and Proton Source

In the context of carbon reutilization, the direct electroreduction of captured CO2 (c-CO2RR) appears as an appealing approach since it avoids the energetically costly separation of CO2 from the capture agent. In this process, CO2 is directly reduced from its captured form. Here, we investigate the influence of the capture agent and proton source on that reaction from a combination of theory and experiment. Specifically, we consider methoxide-captured CO2, NH3-captured CO2, and bicarbonate on silver electrocatalysts. We show that the proton source plays a key role in the interplay of the chemistries for the electroreduction of protons, free CO2, and captured CO2. Our density functional theory calculations, including the influence of the potential, demonstrate that a proton source with smaller pKa improves the reactivity for c-CO2RR, but also increases the selectivity toward the hydrogen evolution reaction (HER) on silver surfaces. Since c-CO2RR requires an additional chemical protonation step, the influence of the proton source is stronger than that of the HER. However, c-CO2RR cannot compete with the HER on Ag, Experimentally, the dominant product observed is H2 with low amounts of CO being produced. Through a rotating cylinder electrode cell of well-defined mass-transport properties, we conclude that although methanol solvent exhibits a lower HER activity, HER remains dominant over c-CO2RR. Our work suggests that methoxide is a potential alternative capture agent to NH3 for direct reduction of captured CO2, though challenges in catalyst design, particularly in reducing the onset potential of c-CO2RR to surpass the HER, remain to be addressed.

Molecular insights into the interactions between PEG carriers and drug molecules from Celastrus hindsii: a multi-scale simulation study

Efficient drug delivery is crucial for the creation of effective pharmaceutical treatments, and polyethylene glycol (PEG) carriers have been emerged as promising candidates for this purpose due to their bio-compatibility, enhancement of drug solubility, and stability. In this study, we utilized molecular simulations to examine the interactions between PEG carriers and selected drug molecules extracted from Celastrus hindsii: Hindsiilactone A, Hindsiiquinoflavan B, Maytenfolone A, and Celasdin B. The simulations provided detailed insights into the binding affinity, stability, and structural properties of these drug molecules when complexed with PEG carriers. A multi-scale approach combining density functional theory (DFT), extended tight-binding (xTB), and molecular dynamics (MD) simulations was conducted to investigate both unbound and bound states of PEG/drug systems. The results from DFT and xTB calculations revealed that the unbound complex has an unfavorable binding free energy, primarily due to negative contributions of delta solvation free energy and entropy. The MD simulations provided more detailed insights into the interactions between PEG and drug molecules in water solutions. By integrating the findings from the multi-scale simulations, a comprehensive picture of the unbound and bound states of PEG and drug systems were obtained. This information is valuable for understanding the molecular mechanisms governing the binding of drugs in PEG-based delivery platforms, and it contributes to the rational design and optimization of these systems.

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.

Computational Design of an Electro-Organocatalyst for Conversion of CO2 into Long Chain Aldehydes

Density functional theory calculations employing a hybrid implicit/explicit solvation method were used to demonstrate that an electro-organocatalyst designed in our previous work for reducing CO2 to formaldehyde could also be capable of coupling formaldehyde to form long chain aldehydes. The catalytic activity is enabled by an electron-rich vicinal enediamine (>N-C═C-N<) backbone that activates formaldehyde by reversing the polarity on the carbon atom, enabling it to act as a nucleophile in the subsequent aldol addition step. The catalyst then enables reductive dehydroxylation of the aldol addition product by facilitating outer-sphere electron transfer. The optimal pH as well as the limiting potential and formaldehyde concentration are identified and related to the kinetic balance between several rate limiting steps. Finally, the optimal conditions for coupling with the CO2 reduction cycle are discussed, demonstrating that the electro-organocatalyst is capable of efficiently converting CO2 into aldehyde products with a turnover frequency (per carbon atom) on the order of 0.1-1 s-1.

Potential-dependent activities in interpreting the reaction mechanism of dual-metal atom catalysts for Li-CO2 batteries

CO2 electrochemistry has been considered as a promising cathode reaction for energy storage due to its high theoretical energy density, high electrochemical potential, and ability to fix CO2. However, the low efficiency and poor reversibility of Li-CO2 evolution significantly impede the applications of Li-CO2 batteries. Herein, first-principles calculations were employed to investigate the 21 M1M2N4C dual-atom catalysts and explore the catalytic mechanism for the Li-CO2 evolution reaction. Among these dual-atom catalysts, the MoMoN4C shows the highest adsorption interaction with CO2 due to its high d-center and d-p orbital coupling. The effects of dual-atom sites on the catalytic activities and selectivities were investigated by searching the possible reaction pathways toward the battery-discharging processes in the ether electrolyte with the help of implicit constant electrode potential simulations. The compared results show that the Li-CO2 discharging process was limited by the rate-determining reactions involving *Li + CO2 → *LiCO2 and *LiC2O4@ + Li+ + e- → *CO + Li2CO3, and these processes on graphene are relatively sluggish due to the low onset potential range of -2 to -2.36 V vs. SHE. By contrast, The optimized onset potentials of -1.15 to -1.31 V vs. SHE were obtained at the MoMoN4C active site. Furthermore, the MoMoN4C active site shows a lower energy barrier for the decomposition of *Li2CO3 than the pure graphene, which reveals the MoMoN4C active site with excellent CO2 activation ability can reduce the polarization of the discharging reactions and energy barrier for the CO bond cleavage. This work provides deep insight into the Li-CO2 evolution mechanisms and guides the design of advanced dual-atom catalysts for highly reversible Li-CO2 batteries.

Cu Based Dilute Alloys for Tuning the C2+ Selectivity of Electrochemical CO2 Reduction

Electrochemical CO2 reduction is a promising technology for replacing fossil fuel feedstocks in the chemical industry but further improvements in catalyst selectivity need to be made. So far, only copper-based catalysts have shown efficient conversion of CO2 into the desired multi-carbon (C2+) products. This work explores Cu-based dilute alloys to systematically tune the energy landscape of CO2 electrolysis toward C2+ products. Selection of the dilute alloy components is guided by grand canonical density functional theory simulations using the calculated binding energies of the reaction intermediates CO*, CHO*, and OCCO* dimer as descriptors for the selectivity toward C2+ products. A physical vapor deposition catalyst testing platform is employed to isolate the effect of alloy composition on the C2+/C1 product branching ratio without interference from catalyst morphology or catalyst integration. Six dilute alloy catalysts are prepared and tested with respect to their C2+/C1 product ratio using different electrolyzer environments including selected tests in a 100-cm2 electrolyzer. Consistent with theory, CuAl, CuB, CuGa and especially CuSc show increased selectivity toward C2+ products by making CO dimerization energetically more favorable on the dominant Cu facets, demonstrating the power of using the dilute alloy approach to tune the selectivity of CO2 electrolysis.

Ab Initio Studies of Electrocatalytic CO2 Reduction for Small Cu Cluster Supported on Polar Substrates

Small Cu clusters are excellent candidates for the electrocatalytic reduction of carbon dioxide (CO2RR), and their catalytic performance is expected to be significantly influenced by the interaction between the substrate and cluster. In this study, we systematically investigate the CO2RR for a Cu3 cluster anchored on Janus MoSX (X = Se, Te) substrates using density functional theory calculations. These substrates feature a broken vertical mirror symmetry, which generates spontaneous out-of-plane polarization and offers two distinct polar surfaces to support the Cu3 cluster. Our findings reveal that the CO2RR performance on the Cu3 cluster is strongly influenced by the polarization direction and strength of the MoSX (X = Se, Te) substrates. Notably, the Cu3 cluster supported on the S-terminated MoSTe surface (Cu3(S)@MoSTe) demonstrates the highest CO2RR activity, producing methane. These results underscore the pivotal role of substrate polarization in modulating the binding strength of reactants and reaction intermediates, thereby enhancing the CO2RR efficiency.

Unlocking the Potential of Bi2S3-Derived Bi Nanoplates: Enhanced Catalytic Activity and Selectivity in Electrochemical and Photoelectrochemical CO2 Reduction to Formate

Various electrocatalysts are extensively examined for their ability to selectively produce desired products by electrochemical CO2 reduction reaction (CO2RR). However, an efficient CO2RR electrocatalyst doesn't ensure an effective co-catalyst on the semiconductor surface for photoelectrochemical CO2RR. Herein, Bi2S3 nanorods are synthesized and electrochemically reduced to Bi nanoplates that adhere to the substrates for application in the electrochemical and photoelectrochemical CO2RR. Compared with commercial-Bi, the Bi2S3-derived Bi (S-Bi) nanoplates on carbon paper exhibit superior electrocatalytic activity and selectivity for formate (HCOO-) in the electrochemical CO2RR, achieving a Faradaic efficiency exceeding 93%, with minimal H2 production over a wide potential range. This highly selective S-Bi catalyst is being employed on the Si photocathode to investigate the behavior of electrocatalysts during photoelectrochemical CO2RR. The strong adhesion of the S-Bi nanoplates to the Si nanowire substrate and their unique catalytic properties afford exceptional activity and selectivity for HCOO- under simulated solar irradiation. The selectivity observed in electrochemical CO2RR using the S-Bi catalyst correlates with that seen in the photoelectrochemical CO2RR system. Combined pulsed potential methods and theoretical analyses reveal stabilization of the OCHO* intermediate on the S-Bi catalyst under specific conditions, which is critical for developing efficient catalysts for CO2-to-HCOO- conversion.

The Selectivity Origins in Ag-Catalyzed CO2 Electroreduction

Ag exhibits high selectivity of electrochemical CO2 reduction (CO2R) toward C1 products, while the hydrogenation involving the concerted proton-electron transfer (CPET) or sequential electron-proton transfer (SEPT) mechanism is still in debate. Toward a better understanding of the Ag-catalyzed electrochemical CO2R, we employed a microkinetic model based on the Marcus electron transfer theory to thoroughly investigate the selectivity of C1 products of electrochemical CO2R over the Ag(111) surface. We found that at an acidic condition of pH = 1.94, formate is the main product when U < -0.94 V via the CPET mechanism, whereas CO becomes the primary product when U > -0.94 V via the SEPT mechanism. Conversely, at an alkaline condition of pH = 13.95, formate is the main product following the SEPT mechanism. Our findings provide novel insights into the influence of external factors (applied potential and pH) on the product selectivity and hydrogenation mechanism of electrochemical CO2R.

In-situ dynamic catalysis electrode electrolyte interphase enabling Mg2+ insertion in 2D metal-organic polymer for high-capacity and long-lifespan magnesium batteries

Rechargeable magnesium battery is regarded as the promising candidate for the next generation of high-specific-energy storage systems. Nevertheless, issues related to severe Mg-Cl dissociation at the electrolyte-electrode interface impede the insertion of Mg2+ into most materials, leading to severe polarization and low utilization of Mg-storage electrodes. In this study, a metal-organic polymer (MOP) Ni-TABQ (Ni-coordinated tetramino-benzoquinone) with superior surface catalytic activity is proposed to achieve the high-capacity Mg-MOP battery. The layered Ni-TABQ cathode, featuring a unique 2D π-d linear conjugated structure, effectively reduces the dissociation energy of MgxCly clusters at the Janus interface, thereby facilitating Mg2+ insertion. Due to the high utilization of active sites, Ni-TABQ achieves high capacities of 410 mAh/g at 200 mA g-1, attributable to a four-electron redox process involving two redox centers, benzoid carbonyls, and imines. This research highlights the importance of surface electrochemical processes in rechargeable magnesium batteries and paves the way for future development in multivalent metal-ion batteries.

Alkali cation-induced cathodic corrosion in Cu electrocatalysts

The reconstruction of Cu catalysts during electrochemical reduction of CO2 is a widely known but poorly understood phenomenon. Herein, we examine the structural evolution of Cu nanocubes under CO2 reduction reaction and its relevant reaction conditions using identical location transmission electron microscopy, cyclic voltammetry, in situ X-ray absorption fine structure spectroscopy and ab initio molecular dynamics simulation. Our results suggest that Cu catalysts reconstruct via a hitherto unexplored yet critical pathway - alkali cation-induced cathodic corrosion, when the electrode potential is more negative than an onset value (e.g., -0.4 VRHE when using 0.1 M KHCO3). Having alkali cations in the electrolyte is critical for such a process. Consequently, Cu catalysts will inevitably undergo surface reconstructions during a typical process of CO2 reduction reaction, resulting in dynamic catalyst morphologies. While having these reconstructions does not necessarily preclude stable electrocatalytic reactions, they will indeed prohibit long-term selectivity and activity enhancement by controlling the morphology of Cu pre-catalysts. Alternatively, by operating Cu catalysts at less negative potentials in the CO electrochemical reduction, we show that Cu nanocubes can provide a much more stable selectivity advantage over spherical Cu nanoparticles.