Noncovalent interactions that tune the reactivities of the flavins in bifurcating electron transferring flavoprotein.

Bifurcating electron transferring flavoproteins (Bf-ETFs) tune chemically identical flavins to two contrasting roles. To understand how, we used hybrid quantum mechanical molecular mechanical calculations to characterize noncovalent interactions applied to each flavin by the protein. Our computations replicated the differences between the reactivities of the flavins: the electron transferring flavin (ETflavin) was calculated to stabilize anionic semiquinone (ASQ) as needed to execute its single-electron transfers, whereas the Bf flavin (Bfflavin) was found to disfavor the ASQ state more than does free flavin and to be less susceptible to reduction. The stability of ETflavin ASQ was attributed in part to H-bond donation to the flavin O2 from a nearby His side chain, via comparison of models employing different tautomers of His. This H-bond between O2 and the ET site was uniquely strong in the ASQ state, whereas reduction of ETflavin to the anionic hydroquinone (AHQ) was associated with side chain reorientation, backbone displacement, and reorganization of its H-bond network including a Tyr from the other domain and subunit of the ETF. The Bf site was less responsive overall, but formation of the Bfflavin AHQ allowed a nearby Arg side chain to adopt an alternative rotamer that can H-bond to the Bfflavin O4. This would stabilize the anionic Bfflavin and rationalize effects of mutation at this position. Thus, our computations provide insights on states and conformations that have not been possible to characterize experimentally, offering explanations for observed residue conservation and raising possibilities that can now be tested.

All-body concept and quantified limits of cooperativity and related effects in homodromic cyclic water clusters from a molecular-wide and electron density-based approach.

We strongly advocate distinguishing cooperativity from cooperativity-induced effects. From the MOWeD-based approach, the origin of all-body cooperativity is synonymous with physics- and quantum-based processes of electron (e) delocalization throughout water clusters. To this effect, over 10 atom-pairs contribute to the total e-density at a BCP(H,O) between water molecules in a tetramer. Intermolecular all-body e-delocalization, that is, cooperativity, is an energy-minimizing process that fully explains non-additive increase in stability of a water molecule in clusters with an increase in their size. A non-linear change in cooperativity and cooperativity-induced effects, such as (i) structural (e.g., a change in d(O,O)) or topological intra- and intermolecular properties in water clusters (e.g., electron density or potential energy density at bond critical points) is theoretically reproduced by the proposed expression. It predicted the limiting value of delocalized electrons by a H2O molecule in homodromic cyclic clusters to be 1.58e. O-atoms provide the vast majority of electrons that "travel throughout a cluster predominantly on a privileged exchange quantum density highway" (⋅⋅⋅O-H⋅⋅⋅O-H⋅⋅⋅O-H⋅⋅⋅) using Bader's classical bond paths as density bridges linking water molecules. There are, however, additional electron exchange channels that are not seen on molecular graphs as bond paths. A 3D visual representation of the "privileged" and "additional" exchange channels as well as detailed intra- and inter-molecular patterns of e-sharing and (de)localizing is presented. The energy stabilizing contribution made by three O-atoms of neighboring water molecules was found to be large (-597 kcal/mol in cyclic hexamer) and 5 times more significant than that of a classical O-H⋅⋅⋅O intermolecular H-bond.

Mutual influence of non-covalent interactions formed by imidazole: A systematic quantum-chemical study.

Imidazole is a five-membered heterocycle that is part of a number of biologically important molecules such as the amino acid histidine and the hormone histamine. Imidazole has a unique ability to participate in a variety of non-covalent interactions involving the NH group, the pyridine-like nitrogen atom or the π-system. For many biologically active compounds containing the imidazole moiety, its participation in formation of hydrogen bond NH⋯O/N and following proton transfer is the key step of mechanism of their action. In this work a systematic study of the mutual influence of various paired combinations of non-covalent interactions (e.g., hydrogen bonds and π-interactions) involving the imidazole moiety was performed by means of quantum chemistry (PW6B95-GD3/def2-QZVPD) for a series of model systems constructed based on analysis of available x-ray data. It is shown that for considered complexes formation of additional non-covalent interactions can only enhance the proton-donating ability of imidazole. At the same time, its proton-accepting ability can be both enhanced and weakened, depending on what additional interactions are added to a given system. The mutual influence of non-covalent interactions involving imidazole can be classified as weak geometric and strong energetic cooperativity-a small change in the length of non-covalent interaction formed by imidazole can strongly influence its strength. The latter can be used to develop methods for controlling the rate and selectivity of chemical reactions involving the imidazole fragment in larger systems. It is shown that the strong mutual influence of non-covalent interactions involving imidazole is due to the unique ability of the imidazole ring to effectively redistribute electron density in non-covalently bound systems with its participation.

Positively Photo-Responsive Adsorption Over Binary Copper Porphyrin Framework and Graphene Film Sorbents.

Photo-responsive adsorption has emerged as a vibrant area because it provides a promising route to reduce the energy consumption of the traditional adsorption separation. However, the current methodology to fabricate photo-responsive sorbents is still subject to the photo-deforming molecular units. In this study, a new initiative of photo-dissociated electron-hole pairs is proposed to generate amazing adsorption activity, and prove its feasibility. Employing CuPP [PP = 5,10,15,20-tetrakis(4-carboxyphenyl)porphyrin] framework nanosheets compounded with graphene, binary film (BF) sorbents are successfully fabricated. The paradigmatic BF nanostructure brings about efficiently photo-excited electron-hole pairs with durable enough lifetime to meet the needs of microscopic adsorption equilibrium, which ultimately alters the electron density distribution of adsorption surface, and thus markedly modulates the adsorption activity. Therefore, an amazing photo-enhanced adsorption capability for the index gas CO can be gotten. Once exposed to the visible-light at 420 nm, the CO adsorption capacity (0 °C, 1 bar) is risen from 0.23 mmol g-1 in the darkness to 1.66 mmol g-1, changed by + 622%. This is essentially different from majority of current photo-responsive sorbents based on photo-deforming molecular units, of which adsorption capability is only decreased with photo-induction, and the maximum rate of change reported is just -54%.

Regulating Local Electron Density of Cyano Sites in Graphitic Nitride Carbon by Giant Internal Electric Field for Efficient CO2 Photoreduction to Hydrocarbons.

Selective photocatalytic CO2 reduction to high-value hydrocarbons using graphitic carbon nitride (g-C3N4) polymer holds great practical significance. Herein, the cyano-functionalized g-C3N4 (CN-g-C3N4) with a high local electron density site is successfully constructed for selective CO2 photoreduction to CH4 and C2H4. Wherein the potent electron-withdrawing cyano group induces a giant internal electric field in CN-g-C3N4, significantly boosting the directional migration of photogenerated electrons and concentrating them nearby. Thereby, a high local electron density site around its cyano group is created. Moreover, this structure can also effectively promote the adsorption and activation of CO2 while firmly anchoring *CO intermediates, facilitating their subsequent hydrogenation and coupling reactions. Consequently, using H2O as a reducing agent, CN-g-C3N4 achieves efficient and selective photocatalytic CO2 reduction to CH4 and C2H4 activity, with maximum rates of 6.64 and 1.35 µmol g-1 h-1, respectively, 69.3 and 53.8 times higher than bulk g-C3N4 and g-C3N4 nanosheets. In short, this work illustrates the importance of constructing a reduction site with high local electron density for efficient and selective CO2 photoreduction to hydrocarbons.

Enhancing Localized Electron Density over Pd1.4Cu Decorated Oxygen Defective TiO2-x Nanoarray for Electrocatalytic Nitrite Reduction to Ammonia.

Electrocatalytic nitrite (NO2 -) reduction to ammonia (NH3) is a promising method for reducing pollution and aiding industrial production. However, progress is limited by the lack of efficient selective catalysts and ambiguous catalytic mechanisms. This study explores the loading of PdCu alloy onto oxygen defective TiO2-x, resulting in a significant increase in NH3 yield (from 70.6 to 366.4 µmol cm-2 h-1 at -0.6 V vs reversible hydrogen electrode) by modulating localized electron density. In situ and operando studies illustrate that the reduction of NO2 - to NH3 involves gradual deoxygenation and hydrogenation. The process also demonstrated excellent selectivity and stability, with long-term durability in cycling and 50 h stability tests. Density functional theory (DFT) calculations elucidate that the introduction of PdCu alloys further amplified electron density at oxygen vacancies (Ovs). Additionally, the Ti─O bond is strengthened as the d-band center of the Ti 3d rising after PdCu loading, facilitating the adsorption and activation of *NO2. Moreover, the presence of Ovs and PdCu alloy lowers the energy barriers for deoxygenation and hydrogenation, leading to high yield and selectivity of NH3. This insight of controlling localized electron density offers valuable insights for advancing sustainable NH3 synthesis methods.

Defect Engineering Simultaneously Regulating Exciton Dissociation in Carbon Nitride and Local Electron Density in Pt Single Atoms Toward Highly Efficient Photocatalytic Hydrogen Production.

The high exciton binding energy (Eb) and sluggish surface reaction kinetics have severely limited the photocatalytic hydrogen production activity of carbon nitride (CN). Herein, a hybrid system consisting of nitrogen defects and Pt single atoms is constructed through a facile self-assembly and photodeposition strategy. Due to the acceleration of exciton dissociation and regulation of local electron density of Pt single atoms along with the introduction of nitrogen defects, the optimized Pt-MCT-3 exhibits a hydrogen production rate of 172.0 µmol h-1 (λ ≥ 420 nm), ≈41 times higher than pristine CN. The apparent quantum yield for the hydrogen production is determined to be 27.1% at 420 nm. The experimental characterizations and theoretical calculations demonstrate that the nitrogen defects act as the electron traps for the exciton dissociation, resulting in a decrease of Eb from 86.92 to 43.20 meV. Simultaneously, the stronger interaction between neighboring nitrogen defects and Pt single atoms directionally drives free electrons to aggregate around Pt single atoms, and tailors the d-band electrons of Pt, forming a moderate binding strength between Pt atoms and H* intermediates.

Linkage-Mediated Electronic Structure Modulation in Multicomponent Covalent Organic Frameworks for Dramatically Promoted Photocatalytic Hydrogen Evolution.

Linkage chemistry is an essential aspect to covalent organic framework (COF) applications; it is highly desirable to precisely modulate electronic structure mediated directly by linkage for efficient COF-based photocatalytic hydrogen evolution, which however, remains substantially challenging. Herein, as a proof of concept, a collection of robust multicomponent pyrene-based COFs with abundant donor-acceptor (D-A) interactions has been judiciously designed and synthesized through molecularly engineering linkage for photogeneration of hydrogen. Controlled locking and conversion of linkage critically contribute to continuously regulating COFs' electronic structures further to optimize photocatalytic activities. Remarkably, the well-modulated optoelectronic properties turn on the average hydrogen evolution rate from zero to 15.67 mmol g-1 h-1 by the protonated quinoline-linked COF decorated with the trifluoromethyl group (TT-PQCOF-CF3). Using diversified spectroscopy and theoretical calculations, we show that multiple modifications toward linkage synergistically lead to the redistribution of charge on COFs with extended π-conjugation and reinforced D-A effect, making TT-PQCOF-CF3 a promising material with significantly boosted carrier separation and migration. This study provides important guidance for the design of high-performance COF photocatalysts based on the strategy of linkage-mediated electronic structure modulation in COFs.

In Silico Design and Characterization of a New Molecular Electride: Li@Calix[3]Pyrrole.

Electrides, in which anionic electrons are localized independently of the atoms in the compound, have shown promise, especially as catalysts and optoelectronic materials. Here, we present a new computationally designed molecular electride, Li@calix[3]pyrrole (Li@C3P). Electron density and electron localization function analyses unequivocally confirm the existence of localized electride electron density, outside the system, independent of any specific atoms. Non-covalent interaction plots further validate the character of the isolated localized electron, suggesting that the system can be accurately represented by Li+@calix[3]pyrrole ⋅ e-, denoting its distinct charge separation. The remarkable non-linear optical properties of Li@C3P, including average polarizability, α ‾ ${\bar{\alpha }}$ =412.4 au, first hyperpolarizability, ß=4.46×104 au, and second hyperpolarizability, γ ∥ ${{\gamma }_{\parallel }}$ =18.40×106 au, are unparalleled in the previously reported and similar Li@C4P molecular electride. Furthermore, energy decomposition analysis in combination with natural orbital for chemical valence theory sheds light on the mechanism of electron density transfer from Li to the C3P cage, yielding the charge-separated Li@C3P complex. In addition to the electron transfer, a key factor to its electride nature is the electronic structure of the CnP cage, which has its lowest unoccupied molecular orbital located in the void adjacent to the N-H groups at the back of the bowl-shaped CnP cage.

Effects of electron density force to 1.0 and fill to 1.0 on VMAT treatment plans for lung SBRT.

PURPOSE: The aim of this study is to determine the effect of forcing and filling the electron density (ED) to 1.0 of the planning target volume (PTV) overdose distribution in lung SBRT treatment leading to shortening patient treatment time and increasing patient comfort by reducing MU/fraction due to ED manipulation effect. METHODS: In this study, 36 lung SBRT plans of 12 suitable patients who prescribed a total dose of 50 Gy in five fractions were generated with Monaco v.5.10 TPS using the Monte Carlo (MC) algorithm and volumetric modulated arc therapy (VMAT) technique by PTV ED values forcing as well as filling to 1.0 and comparatively assessed. The first group of plans was created by using the patient's original ED, second and third groups of plans were reoptimized by forcing and filling the ED of PTV to 1.0, respectively, therefore acquiring a new dose distribution which lead to comparatively assessment the effects of changes in ED on PTV and OAR doses. RESULTS: Assessment of treatment plans revealed that mean MU/fx numbers were decreased by 76% and 75.25% between Groups 1 and 2, Groups 1 and 3, respectively. The number of segments was also reduced in Group 1 by up to 15% compared with Groups 2 and 3. Maximum HI and CI differences for PTV between Groups 1 and 2 were less than 1% and Groups 1 and 3 were 1.5% which indicates all 3 group plans were comparable in terms of dose distribution within PTV. CONCLUSIONS: Forcing and filling the ED of PTV to 1.0 strategy has provided reduced a number of segments and MU/fx without a significant change in PTV mean and maximum doses, thereby decreasing treatment time and patient discomfort during treatment. This process should be considered in line of a potential number of patients as well as prescribed dose and MU/fx numbers.

Time-Resolved Single-Particle X-ray Scattering Reveals Electron-Density Gradients As Coherent Plasmonic-Nanoparticle-Oscillation Source.

Dynamics of optically excited plasmonic nanoparticles are presently understood as a series of scattering events involving the initiation of nanoparticle breathing oscillations. According to established models, these are caused by statistical heat transfer from thermalized electrons to the lattice. An additional contribution by hot-electron pressure accounts for phase mismatches between theory and experimental observations. However, direct experimental studies resolving the breathing-oscillation excitation are still missing. We used optical transient-absorption spectroscopy and time-resolved single-particle X-ray diffractive imaging to access the electron system and lattice. The time-resolved single-particle imaging data provided structural information directly on the onset of the breathing oscillation and confirmed the need for an additional excitation mechanism for thermal expansion. We developed a new model that reproduces all of our experimental observations. We identified optically induced electron density gradients as the initial driving source.

Revealing the Critical Role of Global Electron Density Transfer in the Reaction Rate of Polar Organic Reactions within Molecular Electron Density Theory.

The critical role of global electron density transfer (GEDT) in increasing the reaction rate of polar organic reactions has been studied within the framework of Molecular Electron Density Theory (MEDT). To this end, the series of the polar Diels-Alder (P-DA) reactions of cyclopentadiene with cyanoethylene derivatives, for which experimental kinetic data are available, have been chosen. A complete linear correlation between the computed activation Gibbs free energies and the GEDT taking place at the polar transition state structures (TSs) is found; the higher the GEDT at the TS, the lower the activation Gibbs free energy. An interacting quantum atoms energy partitioning analysis allows for establishing a complete linear correlation between the electronic stabilization of the electrophilic ethylene frameworks and the GEDT taking place at the polar TSs. This finding supports Parr's proposal for the definition of the electrophilicity ω index. The present MEDT study establishes the critical role of the GEDT in the acceleration of polar reactions, since the electronic stabilization of the electrophilic framework with the electron density gain is greater than the destabilization of the nucleophilic one, making a net favorable electronic contribution to the decrease in the activation energy.

Schiff Bases Derived from Pyridoxal 5'-Phosphate and 2-X-Phenylamine (X = H, OH, SH): Substituent Effects on UV-Vis Spectra and Hydrolysis Kinetics.

Schiff bases are compounds that are widely distributed in nature and have practical value for industry and biomedicine. Another important use of Schiff bases is identifying metal ions and different molecules, including proteins. Their proneness to hydrolysis limits the utilization of Schiff bases to mainly non-aqueous solutions. However, by introducing -OH and -SH substituents to aromatic amine-bearing rings, it is possible to increase the resilience of the Schiff base to destruction in water. The present paper discusses how the hydroxyl or thiol group influences the spectral properties and kinetics of the hydrolysis and formation of Schiff bases derived from pyridoxal 5'-phosphate and aniline, 2-hydroxyaniline, and 2-mercaptoaniline using quantum chemical data. The spectral variation between different imines can be explained by taking into account the geometry and frontier molecular orbital alteration induced by the substituents. The changes in the hydrolysis rate are analyzed using the computed values of local reactivity indices.

Electron Density and Molecular Orbital Analyses of the Nature of Bonding in the η3-CCH Agostic Rhodium Complexes Preceding the C-C and C-H Bond Cleavages.

In our recent work, we revisited C-H and C-C bond activation in rhodium (I) complexes of pincer ligands PCP, PCN, PCO, POCOP, and SCS. Our findings indicated that an η3-Csp2Csp3H agostic intermediate acts as a common precursor to both C-C and C-H bond activation in these systems. We explore the electronic structure and bonding nature of these precleavage complexes using electron density and molecular orbital analyses. Using NBO, IBO, and ESI-3D methods, the bonding in the η3-CCH agostic moiety is depicted by two three-center agostic bonds: Rh-Csp2-Csp3 and Rh-Csp3-H, with all three atoms datively bound to Rh(I). IBO analysis specifically highlights the involvement of three orbitals (CC→Rh and CH→Rh σ donation, plus Rh→CCH π backdonation) in both C-C and C-H bond cleavages. NCIPLOT and QTAIM analyses highlight anagostic (Rh-H) or ß-agostic (Rh-Csp2-H) interactions and the absence of Rh-Csp3 interactions. QTAIM molecular graphs suggest bond path instability under dynamic conditions due to the nearness of line and ring critical points. Several low-frequency and low-force vibrational modes interconvert various bonding patterns, reinforcing the dynamic η3-CCH agostic nature. The kinetic preference for C-H bond breaking is attributed to the smaller reduced mass of C-H vibrations compared to C-C vibrations.

Experimental and Theoretical Characterization of 4π-Electron Möbius Aromatic System of a 1,2-Digermacyclobutadiene†.

We present the visualization of the experimental valence electron-density distribution (EDD) in the isolated 1,2-digermacyclobutadiene ring system, revealing the unique 4π electron-delocalization on the four-membered Ge2C2 ring. A remarkably high Möbius 4π-electron aromatic character in the Ge2C2 ring can be suggested from theoretical calculations, in sharp contrast to the significant antiaromaticity of the all-carbon cyclobutadiene ring.

Precisely Engineering Asymmetric Atomic CoN4 by Electron Donating and Extracting for Oxygen Reduction Reaction.

The development of nonpyrolytic catalysts featuring precisely defined active sites represents an effective strategy for investigating the fundamental relationship between the catalytic activity of oxygen reduction reaction (ORR) catalysts and their local coordination environments. In this study, we have synthesized a series of model electrocatalysts with well-defined CoN4 centers and nonplanar symmetric coordination structures. These catalysts were prepared by a sequential process involving the chelation of cobalt salts and 1,10-phenanthroline-based ligands with various substituent groups (phen(X), where X=OH, CH3, H, Br, Cl) onto covalent triazine frameworks (CTFs). By modulating the electron-donating or electron-withdrawing properties of the substituent groups on the phen-based ligands, the electron density surrounding the CoN4 centers was effectively controlled. Our results demonstrated a direct correlation between the catalytic activity of the CoN4 centers and the electron-donating ability of the substituent group on the phenanthroline ligands. Notably, the catalyst denoted as BCTF-Co-phen(OH), featuring the electron-donating OH group, exhibited the highest ORR catalytic activity. This custom-crafted catalyst achieved a remarkable half-wave potential of up to 0.80 V vs. RHE and an impressive turnover frequency (TOF) value of 47.4×10-3 Hz at 0.80 V vs. RHE in an alkaline environment.

IGMPlot: A program to identify, characterize, and quantify molecular interactions.

We describe the development and features of a program called IGMPlot, which is based on the independent gradient model (IGM) and its local descriptor δ g . The IGM approach analyzes the gradient of the electron density (ED) in a molecular system to identify regions of space where chemical interactions take place. IGMPlot is intended for use by both experimental scientists and theoretical chemists. It is standalone software written in C++, with versions available for multiple platforms. Some key features are: probing and quantifying interactions between two given molecular fragments, determining bond strength (IBSI), estimating the atomic contributions to an intermolecular interaction and preparing data to build 2D and 3D representations of interaction regions. The software has been updated to include new features: critical point analysis of the ED, assessing ED asymmetry of a given bond (PDA) and a new descriptor called q g designed to enhance the IGM- δ g analysis. The program can be found at: http://igmplot.univ-reims.fr.

An information-theoretic approach to basis-set fitting of electron densities and other non-negative functions.

The numerical ill-conditioning associated with approximating an electron density with a convex sum of Gaussian or Slater-type functions is overcome by using the (extended) Kullback-Leibler divergence to measure the deviation between the target and approximate density. The optimized densities are non-negative and normalized, and they are accurate enough to be used in applications related to molecular similarity, the topology of the electron density, and numerical molecular integration. This robust, efficient, and general approach can be used to fit any non-negative normalized functions (e.g., the kinetic energy density and molecular electron density) to a convex sum of non-negative basis functions. We present a fixed-point iteration method for optimizing the Kullback-Leibler divergence and compare it to conventional gradient-based optimization methods. These algorithms are released through the free and open-source BFit package, which also includes a L2-norm squared optimization routine applicable to any square-integrable scalar function.

Bonding in the high spin lithium clusters: Non-nuclear attractors play a crucial role.

The bonding in lithium high-spin clusters contradicts the usual chemical bonding concept since there are no electron pairs between the atoms, and they are bound with parallel spin electrons. Quantum theory of atoms in molecules and interacting quantum atom analysis (IQA) were used to investigate the nature of bonding in the high-spin Li n n + 1 n = 2 - 5 clusters. Our findings demonstrate that the non-nuclear attractors (NNAs) are an essential component of the high-spin lithium clusters and play a key role in keeping them stable. Based on IQA energy terms, an electrostatic destabilizing interaction between the lithium atoms works against the cluster formation. On the other hand, the interactions between lithium atoms and NNA basins are stabilizing and outweigh the lithium-lithium destabilizing effects. In fact, NNAs tend to draw lithium atoms together and stabilize the resulting cluster. The high-spin clusters of lithium can be regarded as electrostatically driven compounds since the electrostatic components are primarily responsible for the stabilizing interactions between NNAs and Li atoms. The only exception is 3 Li2 , which lacks NNA and has a non-repellent lithium-lithium interaction. Indeed, in the 3 Li2 , the interatomic electrostatic component is negligibly small, and the exchange-correlation term leads to a weak bonding interaction.

Axial Chlorine Induced Electron Delocalization in Atomically Dispersed FeN4 Electrocatalyst for Oxygen Reduction Reaction with Improved Hydrogen Peroxide Tolerance.

Atomically dispersed iron sites on nitrogen-doped carbon (Fe-NC) are the most active Pt-group-metal-free catalysts for oxygen reduction reaction (ORR). However, due to oxidative corrosion and the Fenton reaction, Fe-NC catalysts are insufficiently active and stable. Herein, w e demonstrated that the axial Cl-modified Fe-NC (Cl-Fe-NC) electrocatalyst is active and stable for the ORR in acidic conditions with high H2 O2 tolerance. The Cl-Fe-NC exhibits excellent ORR activity, with a high half-wave potential (E1/2 ) of 0.82 V versus a reversible hydrogen electrode (RHE), comparable to Pt/C (E1/2 = 0.85 V versus RHE) and better than Fe-NC (E1/2 = 0.79 V versus RHE). X-ray absorption spectroscopy analysis confirms that chlorine is axially integrated into the FeN4. More interestingly, compared to Fe-NC, the Fenton reaction is markedly suppressed in Cl-Fe-NC. In situ electrochemical impedance spectroscopy reveals that Cl-Fe-NC provides efficient electron transfer and faster reaction kinetics than Fe-NC. Density functional theory calculations reveal that incorporating Cl into FeN4 can drive the electron density delocalization of the FeN4 site, leading to a moderate adsorption free energy of OH* (∆GOH* ), d-band center, and a high onset potential, and promotes the direct four-electron-transfer ORR with weak H2 O2 binding ability compared to Cl-free FeN4, indicating superior intrinsic ORR activity.