Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation

First-Principles Studies of the Electronic and Optical Properties of Two-Dimensional Arsenic-Phosphorus (2D As-P) Compounds

In this work, we propose the construction of a two-dimensional system based on the stable phases previously reported for the 2D arsenic and phosphorus compounds, with hexagonal and orthorhombic symmetries. Therefore, we have modeled one hexagonal and three possible orthorhombic structures. To ensure the dynamical stability, we performed phonon spectra calculations for each system. We found that all phases are dynamically stable. To ensure the thermodynamic and mechanical stabilities, we have calculated cohesive energies and elastic constants. Our results show that the criteria for the stabilities are all fulfilled. For these stable structures, we computed the electronic and optical properties from first-principles studies based on density functional theory. The computation of electronic band gaps was performed by using the GW approximation to overcome the underestimation of the results obtained from standard DFT approaches. To study the optical properties, we have computed the dielectric function imaginary part within the BSE approach, which takes into account the excitonic effects and allows us to calculate the exciton binding energies of each system. The study was complemented by the computation of the absorption coefficient. From our calculations, it can be established that the 2D As-P systems are good candidates for several technological applications.

Catalysts with Trimetallic Sites on Graphene-like C2N for Electrocatalytic Nitrogen Reduction Reaction: A Theoretical Investigation

Electrocatalytic nitrogen reduction reaction (NRR) is a green and highly efficient way to replace the industrial Haber-Bosch process. Herein, clusters consisting of three transition metal atoms loaded on C2N as NRR electrocatalysts are investigated using density functional theory (DFT). Meanwhile, Ca was introduced as a promoter and the role of Ca in NRR was investigated. It was found that Ca anchored to the catalyst can act as an electron donor and effectively promote the activation of N2 on M3. In both M3@C2N and M3Ca@C2N (M=Fe, Co, Ni), the limiting potential (UL) is less negative than that of the Ru(0001) surface and has the ability to suppress the competitive hydrogen evolution reaction (HER). Among them, Fe3@C2N is suggested to be the most promising candidate for NRR with high thermal stability, strong N2 adsorption ability, low limiting potential, and good NRR selectivity. The concepts of trimetallic sites and alkaline earth metal promoters in this work provide theoretical guidance for the rational design of atomically active sites in electrocatalytic NRR.

New global minimum conformers for the Pt 19 and Pt 20 clusters: low symmetric species featuring different active sites

CONTEXT

The study of platinum (Pt) clusters and nanoparticles is essential due to their extensive range of potential technological applications, particularly in catalysis. The electronic properties that yield optimal catalytic performance at the nanoscale are significantly influenced by the size and structure of Pt clusters. This research aimed to identify the lowest-energy conformers for Pt18 , Pt19 , and Pt20 species using Density Functional Theory (DFT). We discovered new low-symmetry conformers for Pt19 and Pt20 , which are 3.0 and 1.0 kcal/mol more stable, respectively, than previously reported structures. Our study highlights the importance of using density functional approximations that incorporate moderate levels of exact Hartree-Fock exchange, alongside basis sets of at least quadruple-zeta quality. The resulting structures are asymmetric with varying active sites, as evidenced by sigma hole analysis on the electrostatic potential surface. This suggests a potential correlation between electronic structure and catalytic properties, warranting further investigation.

METHODS

An equivariant graph neural network interatomic potential (NequIP) within the Atomic Simulation Environment suite (ASE) was used to provide initial geometries of the aggregates under study. DFT calculations were performed with the ORCA 5 package, using functional approximations that included Generalized Gradient Approximation (PBE), meta-GGA (TPSS, M06-L), hybrid (PBE0, PBEh), meta-GGA hybrid (TPSSh), and range-separated hybrid ( ω B97x) functionals. Def2-TZVP and Def2-QZVP as well as members of the cc-pwCVXZ-PP family to check basis set convergence were used. QTAIM calculations were performed using the AIMAll suite. Structures were visualized with the AVOGADRO code.

Mechanisms and selectivity of methanol oxidation on PtRuM3/C-MWCNT (M = Fe and Co) electrocatalysts

Methanol oxidation efficiency and resistance to CO poisoning are the most challenging issues associated with direct methanol fuel cells. Much experimental effort has been undertaken, such as generating Pt-based binary and ternary nanoparticles, creating composite substrates, and fabricating nanoparticles with special shapes, to overcome these drawbacks. Our previous experiment showed that ternary PtRuM3/C-MWCNT (M = Fe and Co; C-MWCNT = carbon Vulcan-multiwalled carbon nanotube) electrocatalysts exhibited high methanol oxidation activity and tolerance to CO poisoning. However, reaction mechanisms on ternary PtRuM3/C-MWCNT (M = Fe and Co) electrocatalysts remain unknown. Therefore, this work is devoted to elucidating the problem using density functional theory calculations and thermodynamic models. Our present study showed that methanol oxidation proceeds via four possible reaction pathways on the surface of PtRuM3/C-MWCNTs, where the most favourable one follows a series of steps converting with a thermodynamic barrier of 0.513 eV for applied potentials of U = 0 V and 1.005 V on PtRuFe3/C-MWCNTs and 0.404 eV for U = 0 V and 0.167 eV for U = 1.005 V on PtRuCo3/C-MWCNTs. We also provide physical insights into the interaction between methanol oxidation intermediates and substrates' surface by analysing electronic properties. Our findings support the results of our previous experiment. The results of this study can be useful for rationally designing the anode for fuel cells.

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.

Theoretical calculations and experimental verification of carbon dioxide reduction electrocatalyzed by metalloporphyrin

Metal-functionalized porphyrin-like graphene structures are promising electrocatalysts for carbon dioxide reduction reaction (CO2RR) as their metal centers can modulate activity. Yet, the role of metal center of metalloporphyrins (MTPPs) in CO2 reaction activity is still lacking deep understanding. Here, CO2RR mechanism on MTPPs with five different metal centers (M = Fe, Co, Cu, Zn and Ni) are examined by first-principles calculations. The *COOH formation is the rate determined step on the five MTPP structures, and the CoTPP exhibits the best CO2RR activity while ZnTPP and NiTPP are the worst, which is also verified by our experiment. The CO2RR activity is controlled by adsorption states of intermediates (*CO, *COOH), i.e., chemisorption (e.g., on CoTPP) and physisorption (on ZnTPP and NiTPP) of intermediates will lead to good and poor activity, respectively. The deeper the d-band center of the porphyrin ring complexed metal atom, the weaker bonding of MTPP with CO and COOH. Theoretical calculations and experimental results indicate that MTPPs with Co and Fe centers lead to a reduction in the energy barriers for the two uphill reaction steps in the electrocatalytic CO2 reduction process, thereby enhancing CO2 reduction electrocatalytic activity. Faradaic efficiency of CO is correlated with the reaction energy barrier of the first proton-coupled electron reduction process, displaying a strong linear correlation. This work provides a fundamental understanding of MTPPs used as electrocatalysts for CO2RR.

Identification of K+-determined reaction pathway for facilitated kinetics of CO2 electroreduction

Cations such as K+ play a key part in the CO2 electroreduction reaction, but their role in the reaction mechanism is still in debate. Here, we use a highly symmetric Ni-N4 structure to selectively probe the mechanistic influence of K+ and identify its interaction with chemisorbed CO2-. Our electrochemical kinetics study finds a shift in the rate-determining step in the presence of K+. Spectral evidence of chemisorbed CO2- from in-situ X-ray absorption spectroscopy and in-situ Raman spectroscopy pinpoints the origin of this rate-determining step shift. Grand canonical potential kinetics simulations - consistent with experimental results - further complement these findings. We thereby identify a long proposed non-covalent interaction between K+ and chemisorbed CO2-. This interaction stabilizes chemisorbed CO2- and thus switches the rate-determining step from concerted proton electron transfer to independent proton transfer. Consequently, this rate-determining step shift lowers the reaction barrier by eliminating the contribution of the electron transfer step. This K+-determined reaction pathway enables a lower energy barrier for CO2 electroreduction reaction than the competing hydrogen evolution reaction, leading to an exclusive selectivity for CO2 electroreduction reaction.

Synergistic Effect of Rutile and Brookite TiO2 for Photocatalytic Formic Acid Dehydrogenation

Solar energy is an ideal clean and inexhaustible energy source. Solar-driven formic acid (FA) dehydrogenation is one of the promising strategies to address safety and cost issues related to the storage, transport, and distribution of hydrogen energy. For FA dehydrogenation, the O-H and C-H cleavages are the key steps, and developing a photocatalyst with the ability to break these two bonds is critical. In this work, both density functional theory (DFT) calculation and experimental results confirmed the positive synergistic effect between brookite and rutile TiO2 for O-H and C-H cleavage in HCOOH. Further, brookite TiO2 is beneficial to the generation of the •OH radical and significantly promotes C-H cleavage in formate. Under optimized conditions, the H2 production efficiency of FA dehydrogenation can reach up to 30.4 μmol·mg-1·h-1, which is the highest value compared with similar reported TiO2-based systems and over 1.7 times the reported highest value of Au0.75Pd0.25/TiO2 photocatalysts. More importantly, after more than 42 days (>500 h) of irradiation, the system still demonstrated high H2 production activity, indicating the potential for practical application. This work provides a valuable strategy to improve both the efficiency and stability of photocatalytic FA dehydrogenation under mild conditions.

A CeO2 (100) surface reconstruction unveiled by in situ STEM and particle swarm optimization techniques

The reconstruction of the polar CeO2 (100) surface has been a subject of long-standing debates due to its complexity and the limited availability of experimental data. Herein, we successfully reveal a CeO2 (100)-(4 × 6) surface reconstruction by combining in situ spherical aberration-corrected scanning transmission electron microscopy, density functional theory calculations, and a particle swarm optimization-based algorithm for structure searching. We have further elucidated the stabilizing mechanism of the reconstructed structure, which involves the splitting of the filled Ce(4f) states and the mixing of the lower-lying ones with the O(2p) orbitals, as evidenced by the projected density of states. We also reveal that the surface chemisorption properties toward water molecules, an important step in numerous heterogeneous catalytic reactions, are enhanced. These insights into the distinct properties of ceria surface pave the way for performance improvements of ceria in a wide range of applications.

Synergistic engineering of heterostructure and oxygen vacancy in cobalt hydroxide/aluminum oxyhydroxide as bifunctional electrocatalysts for urea-assisted hydrogen production

Designing inexpensive, high-efficiency and durable bifunctional catalysts for urea oxidation reaction (UOR) and hydrogen evolution reaction (HER) is an encouraging tactic to produce hydrogen with reduced energy expenditure. Herein, oxygen vacancy-rich cobalt hydroxide/aluminum oxyhydroxide heterostructure on nickel foam (denoted as Co(OH)2/AlOOH/NF-100) has been fabricated using one step hydrothermal process. Theoretical calculation and experimental results indicate the electrons transfer from Co(OH)2 to highly active AlOOH results in the interfacial charge redistribution and optimization of electronic structure. Abundant oxygen vacancies in the heterostructure could improve the conductivity and simultaneously serve as the active sites for catalytic reaction. Consequently, the optimal Co(OH)2/AlOOH/NF-100 demonstrates excellent electrocatalytic performance for HER (62.9 mV@10 mA cm-2) and UOR (1.36 V@10 mA cm-2) due to the synergy between heterointerface and oxygen vacancies. Additionally, the in situ electrochemical impedance spectrum (EIS) for UOR suggests that the heterostructured catalyst exhibits rapid reaction kinetics, mass transfer and current response. Importantly, the urea-assisted electrolysis composed of the Co(OH)2/AlOOH/NF-100 manifests a low cell voltage (1.48 V @ 10 mA cm-2) in 1 M KOH containing 0.5 M urea. This work presents a promising avenue to the development of HER/UOR bifunctional electrocatalysts.

Spatial-coupled Ampere-level Electrochemical Propylene Epoxidation over RuO2/Ti Hollow-fiber Penetration Electrodes

The electrochemical propylene epoxidation reaction (PER) provides a promising route for ecofriendly propylene oxide (PO) production, instantly generating active halogen/oxygen species to alleviate chloride contamination inherent in traditional PER. However, the complex processes and unsatisfactory PO yield for current electrochemical PER falls short of meeting industrial application requirements. Herein, a spatial-coupling strategy over RuO2/Ti hollow-fiber penetration electrode (HPE) is adopted to facilitate efficient PO production, significantly improving PER performance to ampere level (achieving over 80 % PO faradaic efficiency and a maximum PO current density of 859 mA cm-2). The synergetic combination of the penetration effect of HPE and the spatial-coupled reaction sequence, enables the realization of ampere-level PO production with high specificity, exhibiting significant potentials for economically viable PER applications.

First-principles study on the heterogeneous formation of environmentally persistent free radicals (EPFRs) over α-Fe2O3(0001) surface: Effect of oxygen vacancy

Metal oxides with oxygen vacancies have a significant impact on catalytic activity for the transformation of organic pollutants in waste-to-energy (WtE) incineration processes. This study aims to investigate the influence of hematite surface oxygen point defects on the formation of environmentally persistent free radicals (EPFRs) from phenolic compounds based on the first-principles calculations. Two oxygen-deficient conditions were considered: oxygen vacancies at the top surface and on the subsurface. Our simulations indicate that the adsorption strength of phenol on the α-Fe2O3(0001) surface is enhanced by the presence of oxygen vacancies. However, the presence of oxygen vacancies has a negative impact on the dissociation of the phenol molecule, particularly for the surface with a defective point at the top layer. Thermo-kinetic parameters were established over a temperature range of 300-1000 K, and lower reaction rate constants were observed for the scission of phenolic O-H bonds over the oxygen-deficient surfaces compared to the pristine surface. The negative effects caused by the oxygen-deficient conditions could be attributed to the local reduction of FeIII to FeII, which lower the oxidizing ability of surface reaction sites. The findings of this study provide us a promising approach to regulate the formation of EPFRs.

Hierarchical NiCo2Se4 Arrays Composed of Atomically Thin Nanosheets: Simultaneous Improvements in Thermodynamics and Kinetics for Electrocatalytic Water Splitting in Neutral Media

The inefficiency of electrocatalysts for water splitting in neutral media stems from a comprehensive impact of poor intrinsic activity, a limited number of active sites, and inadequate mass transport. Herein, hierarchical ultrathin NiCo2Se4 nanosheets are synthesized by the selenization of NiCo2O4 porous nanoneedles. Theoretical and experimental investigations reveal that the intrinsic hydrogen evolution reaction (HER) activity primarily originate from the NiCo2Se4, whereas the high oxygen evolution reaction (OER) performance is related to the NiCoOOH due to the structural reconstruction. The abundant Se and O vacancies introduced by atomically thin nanostructure modulate the electronic structure of NiCo2Se4 and NiCoOOH, thereby improving the intrinsic HER and OER activities, respectively. COMSOL simulation demonstrate the edges of extended nanosheets from the main body significantly promote the charge aggregation, boosting the reduction and oxidation current during HER/OER process. This charge aggregation effect notably exceeds the tip effect for the nanoneedle, highlighting the unique advantage of the hierarchical nanosheet structure. Benefiting from abundant vacancies and unique nanostructure, the hierarchical ultrathin nanosheet simultaneously improve the thermodynamics and kinetics of the electrocatalyst. The optimized samples display an overpotential of 92 mV for HER and 214 mV for OER at 100 mA cm-2, significantly surpassing the performance of currently reported HER/OER catalysts in neutral media.

Introducing a new correlation functional in density functional theory

The correlation functional holds significance in density functional theory as it addresses electron-electron interactions beyond the mean-field approximation, enhancing the accuracy of total energy calculations, electronic excitations, and the prediction of materials properties. There are several expressions to describe this energy, and each of them has a unique set of errors in calculating particular properties of materials. This work offers a new correlation functional by employing the density's dependence on ionization energy. We theoretically derived this functional and combined it with the previously reported ionization energy dependent exchange functional to investigate its effect on the total energy, bond energy, dipole moment, and zero-point energy of 62 molecules. The comparison of this new functional in respect to existing widely used correlation models including QMC, PBE, B3LYP and Chachiyo models shows how well it works in producing accurate results with minimal mean absolute error.

Prediction of potential high-temperature superconductivity in ternary Y-Hf-H compounds under high pressure

Compressed ternary alloy superhydrides are currently considered to be the most promising competitors for high-temperature superconducting materials. Here, the stable stoichiometries in the Y-Hf-H ternary system under pressure are comprehensively explored in theory and four fresh phases are predicted: Pmna-YHfH6 and P4/mmm-YHfH7 at 200 GPa, P4/mmm-YHfH8 at 300 GPa and P-6m2-YHfH18 at 400 GPa. The four Y-Hf-H ternary phases are thermodynamically and dynamically stable at corresponding pressure. In addition, structural features, bonding characteristics, electronic properties, and superconductivity of the four ternary Y-Hf-H phases are systematically calculated and discussed. As the hydrogen content and the density of states of H atoms at the Fermi level increase, the superconducting transition temperatures (Tc) of Y-Hf-H system are significantly enhanced. The P-6m2-YHfH18 with high hydrogen content exhibits a high calculated Tc value of 130 K at 400 GPa.

N2 as an Efficient IR Probe Molecule for the Investigation of Ceria-Containing Materials

Ceria and ceria-based catalysts are very important in redox and acid-base catalysis. Nanoceria have also been found to be important in biomedical applications. To design efficient materials, it is necessary to thoroughly understand the surface chemistry of ceria, and one of the techniques that provides such information about the surface is the vibrational spectroscopy of probe molecules. Although the most commonly used probe is CO, it has some disadvantages when applied to ceria and ceria-based catalysts. CO can easily reduce the material, forming carbonate-like species, and can be disproportionate, thus modifying the surface. Here, we offer a pioneering study of the adsorption of 15N2 at 100 K, demonstrating that dinitrogen can be more advantageous than CO when studying ceria-based materials. As an inert gas, N2 is not able to oxidize or reduce cerium cations and does not form any surface anionic species able to modify the surface. It is infrared and transparent, and thus there is no need to subtract the gas phase spectrum, something that often increases the noise level. Being a weaker base than CO, N2 has a negligible induction effect. By using stoichiometric nano-shaped ceria samples, we concluded that 15N2 can distinguish between surface Ce4+ sites on different, low index planes; with cations on the {110} facets and on some of the edges, Ce4+-15N2 species with IR bands at 2258-2257 cm-1 are formed. Bridging species, where one of the N atoms from the molecule interacts with two Ce4+ cations, are formed on the {100} facets (2253-2252 cm-1), while the interaction with the {111} facets is very weak and does not lead to the formation of measurable amounts of complexes. All species are formed by electrostatic interaction and disappear during evacuation at 100 K. In addition, N2 provides more accurate information than CO on the acidity of the different OH groups because it does not change the binding mode of the hydroxyls.

Substantial Energy Band Modulation of Semiconducting Hexagonal GaTe Quantum Wells by Layer Thickness and Mirror Twin Boundaries

Exploring emerging two-dimensional (2D) van der Waals (vdW) semiconducting materials and precisely tuning their electronic properties at the atomic level have long been recognized as crucial issues for developing their high-end electronic and optoelectronic applications. As a III-VI semiconductor, ultrathin layered hexagonal GaTe (h-GaTe) remains unexplored in terms of its intrinsic electronic properties and band engineering strategies. Herein, we report the successful synthesis of ultrathin h-GaTe layers on a selected graphene/SiC(0001) substrate, via molecular beam epitaxy (MBE). The widely tunable quasiparticle band gaps (∼2.60-1.55 eV), as well as the vdW quantum well states (QWSs) that can be strictly counted by the layer numbers, are well characterized by onsite scanning tunneling microscopy/spectroscopy (STM/STS), and their origins are clearly addressed by density functional theory (DFT) calculations. More intriguingly, distinctive 8|8E and 4|4P (Ga) mirror twin boundaries (MTBs) are identified in the ultrathin h-GaTe flakes, which can induce decreased band gaps and prominent p-doping effects. This work should deepen our understanding on the electronic tunability of 2D III-VI semiconductors by thickness control and line defect engineering, which may hold promise for fabricating atomic-scale vertical and lateral homojunctions toward ultrascaled electronics and optoelectronics.

Borrowed dislocations for ductility in ceramics

The inherent brittleness of ceramics, primarily due to restricted atomic motions from rigid ionic or covalent bonded structures, is a persistent challenge. This characteristic hinders dislocation nucleation in ceramics, thereby impeding the enhancement of plasticity through a dislocation-engineering strategy commonly used in metals. Finding a strategy that continuously generates dislocations within ceramics may enhance plasticity. Here, we propose a "borrowing-dislocations" strategy that uses a tailored interfacial structure with well-ordered bonds. Such an approach enables ceramics to have greatly improved tensile ductility by mobilizing a considerable number of dislocations in ceramic borrowed from metal through the interface, thereby overcoming the challenge associated with direct dislocation nucleation within ceramics. This strategy provides a way to enhance tensile ductility in ceramics.

Creating Asymmetric Fe-N3C-N Sites in Single-Atom Catalysts Boosts Catalytic Performance for Oxygen Reduction Reaction

Fine tuning of the metal site coordination environment of a single-atom catalyst (SAC) to boost its catalytic activity for oxygen reduction reaction (ORR) is of significance but challenging. Herein, we report a new SAC bearing Fe-N3C-N sites with asymmetric in-plane coordinated Fe-N3C and axial coordinated N atom for ORR, which was obtained by pyrolysis of an iron isoporphyrin on polyvinylimidazole (PVI) coated carbon black. The C@PVI-(NCTPP)Fe-800 catalyst exhibited significantly improved ORR activity (E1/2 = 0.89 V vs RHE) than the counterpart SAC with Fe-N4-N sites in 0.1 M KOH. Significantly, the Zn-air batteries equipped with the C@PVI-(NCTPP)Fe-800 catalyst demonstrated an open-circuit voltage (OCV) of 1.45 V and a peak power density (Pmax) of 130 mW/cm2, outperforming the commercial Pt/C catalyst (OCV = 1.42 V; Pmax = 119 mW/cm2). The density functional theory (DFT) calculations revealed that the d-band center of the asymmetric Fe-N3C-N structure shifted upward, which enhances its electron-donating ability, favors O2 adsorption, and supports O-O bond activation, thus leading to significantly promoted catalytic activity. This research presents an intriguing strategy for the designing of the active site architecture in metal SACs with a structure-function controlled approach, significantly enhancing their catalytic efficiency for the ORR and offering promising prospects in energy-conversion technologies.

Iodine Clusters in the Atmosphere I: Computational Benchmark and Dimer Formation of Oxyacids and Oxides

The contribution of iodine-containing compounds to atmospheric new particle formation is still not fully understood, but iodic acid and iodous acid are thought to be significant contributors. While several quantum chemical studies have been carried out on clusters containing iodine, there is no comprehensive benchmark study quantifying the accuracy of the applied methods. Here, we present the first study in a series that investigate the role of iodine species in atmospheric cluster formation. In this work, we have studied the iodic acid, iodous acid, iodine tetroxide, and iodine pentoxide monomers and their dimers formed with common atmospheric precursors. We have tested the accuracy of commonly applied methods for calculating the geometry of the monomers, thermal corrections of monomers and dimers, the contribution of spin-orbit coupling to monomers and dimers, and finally, the accuracy of the electronic energy correction calculated at different levels of theory. We find that optimizing the structures either at the ωB97X-D3BJ/aug-cc-pVTZ-PP or the M06-2X/aug-cc-pVTZ-PP level achieves the best thermal contribution to the binding free energy. The electronic energy correction can then be calculated at the ZORA-DLPNO-CCSD(T0) level with the SARC-ZORA-TZVPP basis for iodine and ma-ZORA-def2-TZVPP for non-iodine atoms. We applied this methodology to calculate the binding free energies of iodine-containing dimer clusters, where we confirm the qualitative trends observed in previous studies. However, we identify that previous studies overestimate the stability of the clusters by several kcal/mol due to the neglect of relativistic effects. This means that their contributions to the currently studied nucleation pathways of new particle formation are likely overestimated.

Shift of the reduction potential of nickel(II) Schiff base complexes in the presence of redox innocent metal ions

With the objective of gaining insight into the modulation of the reduction potential of the Ni(II/I) couple, we have synthesized two mononuclear nickel(II) complexes, NiLen (H2Len = N,N'-bis(3-methoxysalicylidene)-1,2-diamino-2-methylpropane) and NiLpn (H2Lpn = N,N'-bis(3-methoxysalicylidene)-1,3-diamino-2,2-dimethylpropane) of two N2O4 donor ligands and recorded their cyclic voltammograms. Both the nickel complexes show reversible reduction processes for the Ni(II/I) couple in acetonitrile solution but the reduction potential of NiLpn (E1/2 = -1.883 V) is 188 mV more positive than that of NiLen (E1/2 = -2.071 V). In the presence of redox inactive metal ions (Li+, Na+, K+, Mg2+, Ca2+ and Ba2+), the reduction potentials are shifted by 49-331 mV and 99-435 mV towards positive values compared to NiLen and NiLpn, respectively. The shift increases with the decrease of the pKa of the respective aqua-complexes of the metal ion but is poorly co-linear; however, better linearity is found when the shift of the mono- and bi-positive metal ion aqua complexes is plotted separately. Spectrophotometric titrations of these two nickel complexes with the guest metal ions in acetonitrile showed a well-anchored isosbestic point in all cases, confirming the adduct formation of NiLen and NiLpn with the metal ions. Structural analysis of single crystals, [(NiLen)Li(H2O)2]·ClO4 (1), [(NiLpn)Li(H2O)]·ClO4 (2), [(NiLpn)2Na]·BF4 (3) and [(NiLpn)2Ba(H2O)(ClO4)]·ClO4 (4), also corroborates the heterometallic adduct formation. The orbital energies of the optimised heterometallic adducts from which electron transfers originated were calculated in order to explain the observed reduction process. A strong linear connection between the calculated orbital energies and the experimental E1/2 values was observed. According to MEP and 2D vector field plots, the largest shift for divalent metal ions is most likely caused by the local electric field that they impose in addition to Lewis acidity.

A comparative DFT study of HCHO decomposition on different terminations of the Co3O4(110) surface

Density functional theory calculations have been performed to compare the HCHO decomposition on Co3O4(110)-A and (110)-B terminations. The results showed that the energy barriers of the two C-H bond cleavages of HCHO on the (110)-A termination were lower than those on the (110)-B termination, suggesting that the (110)-A termination had stronger HCHO decomposition ability than the (110)-B termination. Electronic structures revealed that the stronger HCHO decomposition ability of the (110)-A termination might be ascribed to the strong covalent bond between HCHO and the (110)-A termination, as well as the higher d-band center of Co3+ ions on the (110)-A termination. Furthermore, we proposed that the preparation of Co3O4 under oxygen-rich growth conditions was beneficial to HCHO decomposition because the (110)-A termination was more stable under oxygen-rich conditions.

First-principles investigations of Fe-based A3BX ceramics with high stiffness and damage tolerance

In the search for high-stiffness and damage-tolerant materials, Fe-based A3BX carbide and nitride anti-perovskites were studied using first-principles calculations. These perovskites were found to be stable in cubic structures, as substantiated by the formation energy, elastic Born stability criterion, and phonon dispersion spectrum analysis. The GGA functional was applied for geometry optimization, and the lattice constants are found to be 3.730 Å, 3.715 Å, 3.832 Å, and 3.828 Å for Fe3AlC, Fe3AlN, Fe3SnC, and Fe3SnN, respectively. Elastic property analysis reveals that all the materials have large elastic moduli, high sound velocities, and high Debye temperatures. Among them, carbides have superior stiffness and quasi-ductile properties, and they can be further improved by applying additional pressure. Preliminary analysis of electronic properties indicates that they are ferromagnetic and metallic compounds. Their high melting temperatures (>2600 K) confirm their potential in high-temperature applications. The lowest thermal conductivity of Fe3SnN suggests its potential in efficient solid-state refrigeration application. Moreover, Fe3SnC is proposed to be a viable damage-tolerant material with good prospects. Under 10 GPa external pressure, it possesses a ductile structure with a Young's modulus of 402.15 GPa and bulk modulus of 280.25 GPa.

Synthesis of a stable crystalline nitrene

Nitrenes are a highly reactive, yet fundamental, compound class. They possess a monovalent nitrogen atom and usually a short life span, typically in the nanosecond range. Here, we report on the synthesis of a stable nitrene by photolysis of the arylazide MSFluindN3 (1), which gave rise to the quantitative formation of the arylnitrene MSFluindN (2) (MSFluind is dispiro[fluorene-9,3'-(1',1',7',7'-tetramethyl-s-hydrindacen-4'-yl)-5',9''-fluorene]) that remains unchanged for at least 3 days when stored under argon atmosphere at room temperature. The extraordinary life span permitted the full characterization of 2 by single-crystal x-ray crystallography, electron paramagnetic resonance spectroscopy, and superconducting quantum interference device magnetometry, which supported a triplet ground state. Theoretical simulations suggest that in addition to the kinetic stabilization conferred by the bulky MSFluind aryl substituent, electron delocalization across the central aromatic ring contributes to the electron stabilization of 2.

Exciton-Driven and Layer-Independent Linear and Nonlinear Optical Properties in NbOCl2

We investigate the electronic structure and linear and nonlinear [second-harmonic generation (SHG)] spectra of the NbOCl2 monolayer, bilayer, and bulk by using a real-time first-principles approach based on many-body theory. First, the interlayer couplings between NbOCl2 layers are very weak, due to the relatively large interlayer distance, saturation of the p orbital of Cl atoms, and high degree of localization of charge density around the Nb atom for both the lowest conduction band and the highest valence band. Second, the quasiparticle gaps and exciton binding energy for the three systems show layer-dependent features and decrease with an increase in layer thickness. Most importantly, the linear and SHG spectra of the NbOCl2 monolayer, bilayer, and bulk are dominated by strong excitonic resonances and exhibit layer-independent features due to the weak interlayer couplings. Our findings demonstrate that excitonic effects should be included in studying the optical properties of not only two-dimensional materials but also layered bulk materials with weak interlayer couplings.

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H2O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data. Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from an experiment. The An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Adsorption in Pyrene-Based Metal-Organic Frameworks: The Role of Pore Structure and Topology

Pore topology and chemistry play crucial roles in the adsorption characteristics of metal-organic frameworks (MOFs). To deepen our understanding of the interactions between MOFs and CO2 during this process, we systematically investigate the adsorption properties of a group of pyrene-based MOFs. These MOFs feature Zn(II) as the metal ion and employ a pyrene-based ligand, specifically 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy). Including different additional ligands leads to frameworks with distinctive structural and chemical features. By comparing these structures, we could isolate the role that pore size, the presence of open-metal sites (OMS), metal-oxygen bridges, and framework charges play in the CO2 adsorption of these MOFs. Frameworks with constricted pore structures display a phenomenon known as the confinement effect, fostering stronger MOF-CO2 interactions and higher uptakes at low pressures. In contrast, entropic effects dominate at elevated pressures, and the MOF's pore volume becomes the driving factor. Through analysis of the CO2 uptakes of the benchmark materials ─some with narrower pores and others with larger pore volumes─it becomes evident that structures with narrower pores and high binding energies excel at low pressures. In contrast, those with larger volumes perform better at elevated pressures. Moreover, this research highlights that open-metal sites and inherent charges within the frameworks of ionic MOFs stand out as CO2-philic characteristics.

A thermodynamic model of the surface hydroxylation of γ-Al2O3

Rational design of γ-alumina-based catalysts relies on an extensive understanding of the distribution of hydroxyl groups on the surface of γ-alumina and their physicochemical properties, which remain unclear and challenging to determine experimentally due to the structural complexity. In this work, by means of DFT and thermodynamic calculations, various hydroxylation modes of γ-alumina (110) and (100) surfaces at different OH coverages were evaluated, based on which a thermodynamic model to reflect the relationship between temperature and the surface structure was established and the stable hydroxylation modes under experimental conditions were predicted. This enables us to identify the experimentally measured IR spectra. The effect of hydroxyl coverages on the surface Lewis acidity was then analyzed, showing that the presence of hydroxyl groups could promote the Lewis acidity of neighboring Al sites. This work provides fundamental insights into the molecular level understanding of the surface properties of γ-alumina and benefits the rational design of alumina-based catalysts.

Insights into the Acidic Site in Manganese Oxide in Terms of the Sulfur and Water Tolerance of Low-Temperature NH3 Selective Catalytic Reduction

A critical constraint impeding the utilization of Mn-based oxide catalysts in NH3 selective catalytic reduction (NH3-SCR) is their inadequate resistance to water and sulfur. This vulnerability primarily arises from the propensity of SO2 to bind to the acidic site in manganese oxide, resulting in the formation of metal sulfate and leading to the irreversible deactivation of the catalyst. Therefore, gaining a comprehensive understanding of the detrimental impact of SO2 on the acidic sites and elucidating the underlying mechanism of this toxicity are of paramount importance for the effective application of Mn-based catalysts in NH3-SCR. Herein, we strategically modulate the acidity of the manganese oxide catalyst surface through the incorporation of Ce and Nb. Comprehensive analyses, including thermogravimetry, NH3 temperature-programmed desorption, in situ diffused reflectance infrared Fourier transform spectroscopy, and density functional theory calculations, reveal that SO2 exhibits a propensity for adsorption at strongly acidic sites. This mechanistic understanding underscores the pivotal role of surface acidity in governing the sulfur resistance of manganese oxide.

In-silico mechanistic analysis of adsorption of Iodinated Contrast Media agents on graphene surface

The study aims at assessing the potential of graphene-based adsorbents to reduce environmental impacts of Iodinated Contrast Media Agents (ICMs). We analyze an extensive collection of ICMs. A modeling approach resting on molecular docking and Density Functional Theory simulations is employed to examine the adsorption process at the molecular level. The study also relies on a Quantitative Structure-Activity Relationship (QSAR) modeling framework to correlate molecular properties with the adsorption energy (Ead) of ICMs, thus enabling identification of the key mechanisms underpinning adsorption and of the key factors contributing to it. A collection of distinct QSAR-based models is developed upon relying on Multiple Linear Regression and a standard genetic algorithm method. Having at our disposal multiple models enables us to take into account the uncertainty associated with model formulation. Maximum Likelihood and formal model identification/discrimination criteria (such as Bayesian and/or information theoretic criteria) are then employed to complement the traditional QSAR modeling phase. This has the advantage of (a) providing a rigorous ranking of the alternative models included in the selected set and (b) quantifying the relative degree of likelihood of each of these models through a weight or posterior probability. The resulting workflow of analysis enables one to seamlessly embed DFT and QSAR studies within a theoretical framework of analysis that explicitly takes into account model and parameter uncertainty. Our results suggest that graphene-based surfaces constitute a promising adsorbent for ICMs removal, π-π stacking being the primary mechanism behind ICM adsorption. Furthermore, our findings offer valuable insights into the potential of graphene-based adsorbent materials for effectively removing ICMs from water systems. They contribute to ascertain the significance of various factors (such as, e.g., the distribution of atomic van der Waals volumes, overall molecular complexity, the presence and arrangement of Iodine atoms, and the presence of polar functional groups) on the adsorption process.

Second Harmonic Generation in β-K2TeW3O12: An Acentric Crystal Designed from Centric Phase via Pressure Modulation

The latent value of nonlinear optical (NLO) crystals applied in solid-state laser equipment necessitates the development of applicable strategies for constructing noncentrosymmetric (NCS) crystals. By modulating the synthetic temperature and pressure to achieve the rearrangement of [TeO3]2- groups, a new NCS tellurium tungstate, β-K2TeW3O12 (β-KTW), with a strong second harmonic generation (SHG) response was synthesized based on its centrosymmetric polymorphic phase α-K2TeW3O12 (α-KTW). Computational calculation reveals that the large SHG response of β-KTW (15 × KH2PO4@1064 and 1.5 × KTiOPO4@1950 nm) could be attributed to the uniform arrangement of the NLO-active [TeO3]2- and [WO6]6- groups. β-KTW also exhibits enlarged birefringence (0.196@1064 nm) and a high laser damage threshold (42.3 MW cm-2), showing great potential as a nonlinear crystalline material. This work also provides a new route for the construction of NLO crystals based on centric structure, i.e., reverse pressure regulation.

Unveiling unconventional CH4-Xe compounds and their thermodynamic properties at extreme conditions

Inert gases (e.g., He and Xe) can exhibit chemical activity at high pressure, reacting with other substances to form compounds of unexpected chemical stoichiometry. This work combines first-principles calculations and crystal structure predictions to propose four unexpected stable compounds of CH4Xe3, (CH4)2Xe, (CH4)3Xe, and (CH4)3Xe2 at pressure ranges from 2 to 100 GPa. All structures are composed of isolated Xe atoms and CH4 molecules except for (CH4)3Xe2, which comprises a polymerization product, C3H8, and hydrogen molecules. Ab initio molecular dynamics simulations indicate that pressure plays a very important role in the different temperature driving state transitions of CH4-Xe compounds. At lower pressures, the compounds follow the state transition of solid-plastic-fluid phases with increasing temperature, while at higher pressures, the stronger Xe-C interaction induces the emergence of a superionic state for CH4Xe3 and (CH4)3Xe2 as temperature increases. These results not only expand the family of CH4-Xe compounds, they also contribute to models of the structures and evolution of planetary interiors.

Sodium catalytic phenylpentazole cracking: a theoretical study

This work presents an in-depth investigation into the cracking reaction mechanism of phenylpentazole (C6H5N5) under the catalytic influence of sodium metal, utilizing density functional theory. The geometries of the reactants, transition states, intermediates, and products are meticulously optimized employing the GGA/PW91/DNP level of theory. Also, a rigorous analysis is undertaken, encompassing various key factors including configuration parameters, Mulliken charges, densities of states, and reaction energies. Three distinct reaction pathways are comprehensively examined, shedding light on the intricate details and intricacies of each pathway. The results show that a remarkable outcome in which the activation energy of the C6H5N5 cracking reaction releases N2, facilitated by catalytic metal Na, reveals a strikingly reduced value of a mere 5.2 kcal mol-1 compared to the previously reported activation energies ranging from 20 to 30 kcal mol-1. Evidently, this significantly lowered barrier can be readily surpassed at typical room temperatures, exhibiting practical applicability. Notably, the alkali metal Na effectively serves as a catalyst, successfully diminishing the activation energy required for N2 production through the pyrolysis of pentazole compounds. This breakthrough discovery provides a theoretical basis for experimental research on the low-temperature cracking of pentazole compounds. It also offers valuable insights for the development and application of new high energy density materials, contributing to the creation of a green and low-carbon circular economic system.

Comparison of lattice thermal conductivity usingab-initioDFT, machine learning interatomic potentials, and temperature dependent effective potential: a case study of hexagonal BN and BP bilayer

Discovering high thermal conductivity materials is essential for various practical applications, particularly in electronic cooling. The significance of two-dimensional (2D) materials lies in their unique properties that emerge due to their reduced dimensionality, making them highly promising for a wide range of applications. Hexagonal boron nitride (BN), both monolayer and bilayer forms, has garnered attention for its fascinating properties. In this work, we focus on bilayer boron phosphide (BP), which is isostructural to its BN analogue. The lattice thermal conductivity of both bilayer BN and BP have been calculated usingab-initiodensity functional theory, machine learning with the moment tensor potential method, and the temperature-dependent effective-potential method (TDEP). The TDEP approach gives more accurate results for both BN and BP materials. The lattice thermal conductivity of bilayer BP is lower than that of bilayer BN at room temperature, attributed to increased phonon anharmonicity. This study highlights the importance of understanding phonon scattering mechanisms in determining the thermal conductivity of 2D materials, contributing to the broader understanding and potential applications of these materials in future technologies.

Phase transition and electrochemical properties of S-functionalized MXene anodes for Li-ion batteries: a first-principles investigation

The advancement of anode materials for achieving high energy storage is a crucial topic for high-performance Li-ion batteries (LIBs). Here, first-principles calculations were used to conduct a thorough and systematic investigation into lithium storage properties of MXenes with new S functional groups as LIB anode materials. Density of states, diffusion energy barriers, open circuit voltages and storage capacities were calculated to comprehensively evaluate the lithium storage properties of S-functionalized MXenes. Based on the computational results, Ti2CS2 and V2CS2 were selected as excellent candidates from ten M2CS2 MXenes. The diffusion energy barriers of M2CS2 within the range of 0.26-0.32 eV are lower than those of M2CO2 and M2CF2, indicating that M2CS2 anodes exhibit faster charge/discharge rates. By examining the stable crystal structures and comparing atomic positions before and after Li adsorptions, structural phase transitions during Li-ion adsorptions could happen for nearly all M2CS2 MXenes. The phase transitions predicted were directly observed using ab initio molecular dynamic simulations. The cycle stability, storage capacity and other lithium storage properties were enhanced by the reversible structural phase transition.

Ni─Co─O─S Derived Catalysts on Hierarchical N-doped Carbon Supports with Strong Interfacial Interactions for Improved Hybrid Water Splitting Performance

Simultaneously improving electrochemical activity and stability is a long-term goal for water splitting. Herein, hierarchical N-doped carbon nanotubes on carbon nanowires derived from PPy are grown on carbon cloth, serving as a support for NiCo oxides/sulfides. The hierarchical electrodes annealed in N2 or H2/N2 display improved intrinsic activity and stability for hydrogen evolution reaction (HER) and glucose oxidation reaction. Compared with Pt/C||Ir/C in alkaline media, the glucose electrolysis assembled with electrodes exhibits a cell voltage of 1.38 V at 10 mA cm-2, durability for >12 h at 50 mA cm-2, and resistance to glucose/gluconic acid poisoning. In addition, electrocatalysts can also be applied in ethanol oxidation reactions. Systematic characterizations reveal the strong interactions between NiCo and N-doped carbon support-induced partial charge transfer at the interface and regulate the local electronic structure of active sites. Density functional theory calculations demonstrate that the synergistic effect between N-doped carbon supports, metallic NiCo, and NiCo oxides/sulfides optimize the adsorption energy of H2O and the H* free energy for HER. The energy barrier of the dehydrogenation of glucose effectively decreased. This work will attract attention to the role of metal-support interactions in enhancing the intrinsic activity and stability of electrocatalysts.

Mott-Schottky Construction Boosted Plasmon Thermal and Electronic Effects on the Ag/CoV-LDH Nanohybrids for Highly-Efficient Water Oxidation

Mott-Schottky construction and plasmon excitation represent two highly-efficient and closely-linked coping strategies to the high energy loss of oxygen evolution reaction (OER), but the combined effect has rarely been investigated. Herein, with Ag nanoparticles as electronic structure regulator and plasmon exciter, Ag/CoV-LDH@G nanohybrids (NHs) with Mott-Schottky heterojunction and notable plasmon effect are well-designed. Combining theoretical calculations with experiments, it is found that the Mott-Schottky construction modulates the Fermi level/energy band structure of CoV-LDH, which in turn leads to lowered d-band center (from -0.89 to -0.93), OER energy barrier (from 6.78 to 1.31 eV), and preeminent plasmon thermal/electronic effects. The thermal effect can offset the endothermic enthalpy change of OER, promote the deprotonation of *OOH, and accelerate electron transfer kinetics. Whereas the electronic effect can increase the density of charge carriers (from 0.70 × 1020 to 1.64 × 1020 cm-3), lower the activation energy of OER (from 30.3 to 17.7 kJ mol-1). Benefiting from these favorable factors, the Ag/CoV-LDH@G NHs show remarkable electrocatalytic performances, with an overpotential of 178 and 263 mV to afford 10 and 100 mA cm-2 for OER, respectively, and a low cell voltage of 1.42 V to drive 10 mA cm-2 for overall water splitting under near-infrared light irradiation.

Electrochemically Induced Ru/CoOOH Synergistic Catalyst as Bifunctional Electrode Materials for Alkaline Overall Water Splitting

Efficient and affordable price bifunctional electrocatalysts based on transition metal oxides for oxygen and hydrogen evolution reactions have a balanced efficiency, but it remains a significant challenge to control their activity and durability. Herein, a trace Ru (0.74 wt.%) decorated ultrathin CoOOH nanosheets (≈4 nm) supported on the surface of nickel foam (Ru/CoOOH@NF) is rationally designed via an electrochemically induced strategy to effectively drive the electrolysis of alkaline overall water splitting. The as-synthesized Ru/CoOOH@NF electrocatalysts integrate the advantages of a large number of different HER (Ru nanoclusters) and OER (CoOOH nanosheets) active sites as well as strong in-suit structure stability, thereby exhibiting exceptional catalytic activity. In particular, the ultra-low overpotential of the HER (36 mV) and the OER (264 mV) are implemented to achieve 10 mA cm-2. Experimental and theoretical calculations also reveal that Ru/CoOOH@NF possesses high intrinsic conductivity, which facilitates electron release from H2O and H-OH bond breakage and accelerates electron/mass transfer by regulating the charge distribution. This work provides a new avenue for the rational design of low-cost and high-activity bifunctional electrocatalysts for large-scale water-splitting technology and expects to help contribute to the creation of various hybrid electrocatalysts.

Mechanistic studies of adsorption and ion exchange of Si(OH)4 molecules on the surface of scorodites

Scorodites are commonly used for arsenic immobilization, and it is also the main component of arsenic bearing tailings. Alkali-activated geopolymers are commonly used to landfill arsenic-bearing minerals. However, there no previous studies have explored the interaction between geopolymer molecules and the surface of scorodite. In this paper, Si(OH)4 as a monomer molecule of geopolymer, the mechanism of adsorption and 'ion exchange' between Si(OH)4 molecule and the surface of scorodite during alkali-activation is studied. Results show that the Fe-terminated scorodite (010) surface has high stability. Si(OH)4 are more easily adsorbed on the hollow site of an Fe-terminated scorodite (010) surface, which is described as chemisorption. Compared with Si(OH)4, NaOH is easier to adsorb on an Fe-terminated scorodite (010) surface. The co-adsorption of NaOH and Si(OH)4 on the Fe-terminated scorodite (010) surface was studied, and also belongs to chemical adsorption. When the hydroxyl binds to the As atom, the adsorbed Si(OH)4 is more likely to undergo an 'ion exchange' reaction with the surface, and the reaction is barrierless. The intermediate As(OH)4 produced by the 'ion exchange' reaction can be deprotonated to form an arsenate molecule, which can occur spontaneously. This work reveals that the interaction mechanism of geopolymer molecules on surface of scorodite.

Facet-selective growth of halide perovskite/2D semiconductor van der Waals heterostructures for improved optical gain and lasing

The tunable properties of halide perovskite/two dimensional (2D) semiconductor mixed-dimensional van der Waals heterostructures offer high flexibility for innovating optoelectronic and photonic devices. However, the general and robust growth of high-quality monocrystalline halide perovskite/2D semiconductor heterostructures with attractive optical properties has remained challenging. Here, we demonstrate a universal van der Waals heteroepitaxy strategy to synthesize a library of facet-specific single-crystalline halide perovskite/2D semiconductor (multi)heterostructures. The obtained heterostructures can be broadly tailored by selecting the coupling layer of interest, and can include perovskites varying from all-inorganic to organic-inorganic hybrid counterparts, individual transition metal dichalcogenides or 2D heterojunctions. The CsPbI2Br/WSe2 heterostructures demonstrate ultrahigh optical gain coefficient, reduced gain threshold and prolonged gain lifetime, which are attributed to the reduced energetic disorder. Accordingly, the self-organized halide perovskite/2D semiconductor heterostructure lasers show highly reproducible single-mode lasing with largely reduced lasing threshold and improved stability. Our findings provide a high-quality and versatile material platform for probing unique optoelectronic and photonic physics and developing further electrically driven on-chip lasers, nanophotonic devices and electronic-photonic integrated systems.

Perdew Festschrift editorial

This Special Issue of the Journal of Chemical Physics is dedicated to the work and life of John P. Perdew. A short bio is available within the issue [J. P. Perdew, J. Chem. Phys. 160, 010402 (2024)]. Here, we briefly summarize key publications in density functional theory by Perdew and his collaborators, followed by a structured guide to the papers contributed to this Special Issue.

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.

Going for Gold(-Standard): Attaining Coupled Cluster Accuracy in Oxide-Supported Nanoclusters

The structure of oxide-supported metal nanoclusters plays an essential role in their sharply enhanced catalytic activity over that of bulk metals. Simulations provide the atomic-scale resolution needed to understand these systems. However, the sensitive mix of metal-metal and metal-support interactions, which govern their structure, puts stringent requirements on the method used, requiring calculations beyond standard density functional theory (DFT). The method of choice is coupled cluster theory [specifically CCSD(T)], but its computational cost has so far prevented its application to these systems. In this work, we showcase two approaches to make CCSD(T) accuracy readily achievable in oxide-supported nanoclusters. First, we leverage the SKZCAM protocol to provide the first benchmarks of oxide-supported nanoclusters, revealing that it is specifically metal-metal interactions that are challenging to capture with DFT. Second, we propose a CCSD(T) correction (ΔCC) to the metal-metal interaction errors in DFT, reaching accuracy comparable to that of the SKZCAM protocol at significantly lower cost. This approach forges a path toward studying larger systems at reliable accuracy, which we highlight by identifying a ground-state structure in agreement with experiments for Au20 on MgO, a challenging system where DFT models have yielded conflicting predictions.

Suzuki-Miyaura Cross-Coupling Reaction Catalyzed by Al12M (M = Be, Al, C, and P) Superatoms with Different Numbers of Valence Electrons

The exploration of low-cost, efficient, environmentally safe, and selective catalysts for the activation of carbon-halogen bonds has become an important and challenging topic in modern chemistry. With the help of density functional theory (DFT), it is found that phenyl bromide (PhBr) can be efficiently chemisorbed by the Al12M (M = Be, Al, C, and P) superatoms via forming highly polarized Al-Br covalent bonds, where the C-Br bonds of PhBr can be effectively activated through the electron transfer from Al12M. The different electronic structures of these four Al12M superatoms pose a substantial effect on their performances on the activation of PhBr and the catalytic mechanisms of the Suzuki-Miyaura (SM) reaction. Among them, the alkali-metal-like superatom Al12P exhibits the best performance for the activation of PhBr. In particular, Al13 and Al12P with open-shell electronic structures exhibit catalytic performances comparable to those of previously reported catalysts for this coupling reaction. Hence, it is highly expected that Al13 and Al12P could be used as novel superatom catalysts for C-C coupling reactions and, therefore, open up new possibilities to use nonprecious superatoms in catalyzing the activation and transformation of carbon-halogen bonds.

Deciphering the Catalytic Mechanism of Peroxidase-like Activity of Iron Sulfide Nanozymes

Iron sulfide nanomaterials represented by FeS2 and Fe3S4 nanozymes have attracted increasing attention due to their biocompatibility and peroxidase-like (POD-like) catalytic activity in disease diagnosis and treatments. However, the mechanism responsible for their POD-like activities remains unclear. Herein, taking the oxidation of 3,3,5,5-tetramethylbenzidine (TMB) by H2O2 on FeS2(100) and Fe3S4(001) surfaces, the catalytic mechanism was investigated in detail using density functional theory (DFT) calculations and experimental characterizations. Our experimental results showed that the catalytic activity of FeS2 nanozymes was significantly higher than that of Fe3S4 nanozymes. Our DFT calculations indicated that the surface iron ions of iron sulfide nanozymes could effectively catalyze the production of HO• radicals via the interactions between Fe 3d electrons and the frontier orbitals of H2O2 in the range of -10 to 5 eV. However, FeS2 nanozymes exhibited higher POD-like activity due to the surface Fe(II) binding to H2O2, forming inner-orbital complexes, which results in a larger binding energy and a smaller energy barrier for the base-like decomposition of H2O2. In contrast, the surface iron ions of Fe3S4 nanozymes bind to H2O2, forming outer-orbital complexes, which results in a smaller binding energy and a larger energy barrier for the base-like decomposition of H2O2. The charge transfer analysis showed that FeS2 nanozymes transferred 0.12 e and Fe3S4 nanozymes transferred 0.05 e from their surface iron ions to H2O2, respectively. The simulations were consistent with the experimental observations that the FeS2 nanozymes had a greater affinity for H2O2 compared to that of Fe3S4 nanozymes. This work provides a theoretical foundation for the rational design and accurate preparation of iron sulfide functional nanozymes.

Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes

Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me2-1.1.1, and 36Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3)2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6-2.0 eV) and hence these complexes belong to the superalkali category. The ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3)2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.

Hydrogen storage in M(BDC)(TED)0.5 metal-organic framework: physical insights and capacities

Finding renewable energy sources to replace fossil energy has been an essential demand in recent years. Hydrogen gas has been becoming a research hotspot for its clean and free-carbon energy. However, hydrogen storage technology is challenging for mobile and automotive applications. Metal-organic frameworks (MOFs) have emerged as one of the most advanced materials for hydrogen storage due to their exceptionally high surface area, ultra-large and tuneable pore size. Recently, computer simulations allowed the designing of new MOF structures with significant hydrogen storage capacity. However, no studies are available to elucidate the hydrogen storage in M(BDC)(TED)0.5, where M = metal, BDC = 1,4-benzene dicarboxylate, and TED = triethylenediamine. In this report, we used van der Waals-dispersion corrected density functional theory and grand canonical Monte Carlo methods to explore the electronic structure properties, adsorption energies, and gravimetric and volumetric hydrogen loadings in M(BDC)(TED)0.5 (M = Mg, V, Co, Ni, and Cu). Our results showed that the most favourable adsorption site of H2 in M(BDC)(TED)0.5 is the metal cluster-TED intersection region, in which Ni offers the strongest binding strength with the adsorption energy of -16.9 kJ mol-1. Besides, the H2@M(BDC)(TED)0.5 interaction is physisorption, which mainly stems from the contribution of the d orbitals of the metal atoms for M = Ni, V, Cu, and Co and the p orbitals of the O, C, N atoms for M = Mg interacting with the σ* state of the adsorbed hydrogen molecule. Noticeably, the alkaline-earth metal Mg strongly enhanced the specific surface area and pore size of the M(BDC)(TED)0.5 MOF, leading to an enormous increase in hydrogen storage with the highest absolute (excess) gravimetric and volumetric uptakes of 1.05 (0.36) wt% and 7.47 (2.59) g L-1 at 298 K and 7.42 (5.80) wt% and 52.77 (41.26) g L-1 at 77 K, respectively. The results are comparable to the other MOFs found in the literature.

Regulation of hydrogen binding energy via oxygen vacancy enables an efficient trifunctional Rh-Rh2O3 electrocatalyst for fuel cells and water splitting

Developing multi-functional electrocatalysts is of great practical significance for fuel cells and water splitting. Herein, Rh-Rh2O3 nanoclusters are prepared and the surface oxygen vacancy content is regulated elaborately by post-treatment. The optimized Rh-Rh2O3/C-400 exhibits superior trifunctional catalytic activity for hydrogen oxidation reaction (HOR), hydrogen evolution reaction (HER) and hydrazine oxidation reaction (HzOR), i.e., the mass activity for HOR is 2.29 mA μgRh-1, and the overpotential for HER and HzOR at 10 mA cm-2 is as low as 12 mV and 31 mV, respectively, superior to the benchmark Pt/C. Rh-Rh2O3/C-400 also displays promising performance in practical devices, with the H2-O2 anion-exchange-membrane fuel cell delivering a peak power density of 0.66 W cm-2, and the hydrazine-assisted water splitting electrolyzer requiring a low electrolysis voltage of 0.161 V at 0.1 A cm-2. The experimental and theoretical investigations discover that the hydrogen binding energy (HBE) is linearly depended on surface oxygen vacancy contents, and the HBE directly determines the catalytic activity for HOR, HER and HzOR. This work not only innovates an efficient Rh-based nanocluster tri-functional electrocatalyst, but also eludicates the intrinsic relationship of surface structure-intermediate adsorption-catalytic activity.

Resolution of the reciprocity between radical species from precursor and closed pore formation in hard carbon for sodium storage

Hard carbon (HC) has emerged as a highly promising anode material for sodium ion batteries, drawing tremendous interest in producing this material with low-cost and easily accessible precursors. The determination of the crucial parameters of precursors influencing the formation of key structures, such as closed pores, in the HC is of paramount importance. Considering the potential role of free radicals in the structural evolution of the precursors, we, for the first time, delve into the impact of radical species on the development of closed pores by electron paramagnetic resonance spectroscopy, with petroleum asphalt as the model system. Our findings reveal that carbon centred radicals, with the g value close to that of the free electron (2.0023), exhibit a propensity to form long-range, well-ordered graphitic structures with lower sodium storage capacity. Conversely, the deliberately incorporated oxygen radicals with the g value over 2.005 require a higher energy for ordering the graphitic structures, leading to the creation of closed pores. As a result, the optimal sample showcases a four-fold increase in plateau capacity for sodium ion storage due to the pore filling process. Our research underscores the pivotal role of employing electron paramagnetic resonance spectroscopy studying the critical structural evolution of functional carbon materials.

Single-Atom palladium engineered cobalt nanocomposite for selective aerobic oxidation of sulfides to sulfoxides

Developing efficient catalysts for the selective oxidation of sulfides to sulfoxides using molecular oxygen as the oxidant is a challenging task. Here, we report a novel catalyst comprising a single atom palladium engineered cobalt nanocomposite (denoted as PdCo@NC-800-0.01) for this reaction. The incorporation of single atom palladium effectively transforms an originally inactive cobalt nanocomposite into a highly efficient and selective catalyst for the oxidation of sulfides. This catalyst PdCo@NC-800-0.01 exhibited outstanding performance in the selective oxidation of sulfides to sulfoxides using O2 as the oxidant in the presence of isobutyraldehyde (IBA) under mild conditions, demonstrating high activity and excellent selectivity for a broad spectrum of sulfides with good tolerance toward various functional groups, including those susceptible to oxidation. Furthermore, the catalyst could be easily recovered and reused up to 10 times without any significant loss in activity and selectivity. Comprehensive characterizations and theoretical calculations revealed that the engineering of cobalt nanocomposite with single atom Pd greatly enhanced the ability to adsorb and activate IBA, leading to the generation of the key acyl radical. This radical then reacted with singlet oxygen 1O2 derived from molecular oxygen, producing reactive oxygen species peroxy radical, which ultimately promoted the catalytic performance.