Improved description of the structure of molecular and layered crystals: ab initio DFT calculations with van der Waals corrections

Theoretical study of the catalytic hydrodeoxygenation of furan, methylfuran and benzofurane on MoS2

Herein, we have studied the direct deoxygenation (DDO) (without prior hydrogenation) of furan, 2-methylfuran and benzofuran on the metal edge of MoS2 with a vacancy created under pressure of dihydrogen. For the three molecules, we found that the desorption of the water molecule for the regeneration of the vacancy is the most endothermic. Based on the thermodynamic and kinetic aspects, the reactivity order of the oxygenated compounds is furan ≈ 2-methylfuran > benzofuran, which is in agreement with literature. We present the key stages of the mechanisms and highlight the effects of substituents.

S-Scheme Boron Phosphide/MoS2 Heterostructure with Excellent Light Conversion Ability for Solar Cells and Water Splitting Photocatalysts

Monolayer molybdenum disulfide (MoS2) with a suitable direct band gap and strong optical absorption is very attractive for utilization in solar cells and photocatalytic water splitting. Nevertheless, the broader utilization of MoS2 is impeded by its low carrier mobility and limited responsiveness to infrared light. To overcome these challenges, we constructed a variety of stackings for the boron phosphide (BP)/MoS2 van der Waals heterostructure (vdWH), all of which display S-scheme band alignments except for the AC' stacking. The constituent BP monolayer has superior carrier mobility and strong infrared and visible light response, which makes up for the shortcomings of MoS2. The study revealed that the AB stacking exhibits a remarkable power conversion efficiency of 22.27%, indicating its significant application prospect in solar cells. Additionally, the AB stacking also exhibits a promising application prospect in photocatalytic water splitting due to its suitable band structure, S-scheme band alignment, strong optical adsorption characteristic, high solar-to-hydrogen efficiency, and robust built-in electric field. Meanwhile, applying uniaxial tensile strains along the x-axis direction is more beneficial for photocatalytic water splitting. Hence, the AB-stacked BP/MoS2 vdWH shows significant potential for use in both solar cells and photocatalytic water splitting. This work paves the way for exploring the application of S-scheme heterostructures in solar energy conversion systems.

Voltage-Driven Fluorine Motion for Novel Organic Spintronic Memristor

Integrating tunneling magnetoresistance (TMR) effect in memristors is a long-term aspiration because it allows to realize multifunctional devices, such as multi-state memory and tunable plasticity for synaptic function. However, the reported TMR in different multiferroic tunnel junctions is limited to 100%. This work demonstrates a giant TMR of -266% in La0.6Sr0.4MnO3(LSMO)/poly(vinylidene fluoride)(PVDF)/Co memristor with thin organic barrier. Different from the ferroelectricity-based memristors, this work discovers that the voltage-driven florine (F) motion in the junction generates a huge reversible resistivity change up to 106% with nanosecond (ns) timescale. Removing F from PVDF layer suppresses the dipole field in the tunneling barrier, thereby significantly enhances the TMR. Furthermore, the TMR can be tuned by different polarizing voltage due to the strong modification of spin-polarization at the LSMO/PVDF interface upon F doping. Combining of high TMR in the organic memristor paves the way to develop high-performance multifunctional devices for storage and neuromorphic applications.

Van der Waals polarity-engineered 3D integration of 2D complementary logic

Vertical three-dimensional integration of two-dimensional (2D) semiconductors holds great promise, as it offers the possibility to scale up logic layers in the z axis1-3. Indeed, vertical complementary field-effect transistors (CFETs) built with such mixed-dimensional heterostructures4,5, as well as hetero-2D layers with different carrier types6-8, have been demonstrated recently. However, so far, the lack of a controllable doping scheme (especially p-doped WSe2 (refs. 9-17) and MoS2 (refs. 11,18-28)) in 2D semiconductors, preferably in a stable and non-destructive manner, has greatly impeded the bottom-up scaling of complementary logic circuitries. Here we show that, by bringing transition metal dichalcogenides, such as MoS2, atop a van der Waals (vdW) antiferromagnetic insulator chromium oxychloride (CrOCl), the carrier polarity in MoS2 can be readily reconfigured from n- to p-type via strong vdW interfacial coupling. The consequential band alignment yields transistors with room-temperature hole mobilities up to approximately 425 cm2 V-1 s-1, on/off ratios reaching 106 and air-stable performance for over one year. Based on this approach, vertically constructed complementary logic, including inverters with 6 vdW layers, NANDs with 14 vdW layers and SRAMs with 14 vdW layers, are further demonstrated. Our findings of polarity-engineered p- and n-type 2D semiconductor channels with and without vdW intercalation are robust and universal to various materials and thus may throw light on future three-dimensional vertically integrated circuits based on 2D logic gates.

Comparative analysis of surface phase diagrams in aqueous environment: Implicit vs explicit solvation models

Identifying the stable surface phases under a given electrochemical conditions serves as the basis for studying the atomistic mechanism of reactions at solid/water interfaces. In this work, we systematically compare the performance of the two main approaches that are used to capture the impact of an aqueous environment, implicit and explicit solvent, on surface energies and phase diagrams. As a model system, we consider the magnesium/water interface with (i) Ca substitution and (ii) proton and hydroxyl adsorption. We show that while the implicit solvent model is computationally very efficient, it suffers from two shortcomings. First, the choice of the implicit solvent parameters significantly influences the energy landscape in the vicinity of the surface. The default parameters benchmarked on solvation in water underestimate the energy of the dissolved Mg ion and lead to spontaneous dissolution of the surface atom, resulting in large differences in the surface energetics. Second, in systems containing a charged surface and a solvated ion, the implicit solvent model may not converge to the energetically stable ionic charge state but remain in a high-energy metastable configuration, representing the neutral charge state of the ion. When these two issues are addressed, surface phase diagrams that closely match the explicit water results can be obtained. This makes the implicit solvent model highly attractive as a computationally-efficient surrogate model to compute surface energies and phase diagrams.

Atomic-Layer Controlled Transition from Inverse Rashba-Edelstein Effect to Inverse Spin Hall Effect in 2D PtSe2 Probed by THz Spintronic Emission

2D materials, such as transition metal dichalcogenides, are ideal platforms for spin-to-charge conversion (SCC) as they possess strong spin-orbit coupling (SOC), reduced dimensionality and crystal symmetries as well as tuneable band structure, compared to metallic structures. Moreover, SCC can be tuned with the number of layers, electric field, or strain. Here, SCC in epitaxially grown 2D PtSe2 by THz spintronic emission is studied since its 1T crystal symmetry and strong SOC favor SCC. High quality of as-grown PtSe2 layers is demonstrated, followed by in situ ferromagnet deposition by sputtering that leaves the PtSe2 unaffected, resulting in well-defined clean interfaces as evidenced with extensive characterization. Through this atomic growth control and using THz spintronic emission, the unique thickness-dependent electronic structure of PtSe2 allows the control of SCC. Indeed, the transition from the inverse Rashba-Edelstein effect (IREE) in 1-3 monolayers (ML) to the inverse spin Hall effect (ISHE) in multilayers (>3 ML) of PtSe2 enabling the extraction of the perpendicular spin diffusion length and relative strength of IREE and ISHE is demonstrated. This band structure flexibility makes PtSe2 an ideal candidate to explore the underlying mechanisms and engineering of the SCC as well as for the development of tuneable THz spintronic emitters.

Repulsive Lateral Interaction of Water Molecules at the Initial Stages of Adsorption in Microporous AlPO4-11 According to 27Al NMR and DFT

Lateral (adsorbate-adsorbate) interactions between adsorbed molecules affect various physical and chemical properties of microporous adsorbents and catalysts, influencing their functional properties. In this work, we studied the hydration of microporous AlPO4-11 aluminophosphate, which has an unusually ordered structure upon adsorption of water vapor, and according to 27Al NMR data, only tetrahedrally or octahedrally coordinated Al sites are present in the AlPO4-11. These 27Al NMR data are consistent with the results of density functional theory (DFT) calculations of hydrated AlPO4-11, which revealed the presence of a strong repulsive lateral interaction at the initial stage of adsorption, suppressing the adsorption of water on neighboring (separated by one -O-P-O- bridge) Al crystallographic sites. As a result, of all the different aluminum sites, only half of the Al1 sites adsorb two water molecules and acquire octahedral coordination.

Phase transitions of choline dihydrogen phosphate: A vibrational spectroscopy and periodic DFT study

Choline dihydrogen phosphate, [Chol][H2PO4], is a proton-conducting ionic plastic crystal exhibiting a complicated sequence of phase transitions. Here, we address the argument in the literature around the thermal properties of [Chol][H2PO4] using Raman and infrared microspectroscopy. The known structure of the low-temperature crystal, which contains the anti-conformer of [Chol]+ and hydrogen-bonded dimers of anions, was used to do periodic density functional theory calculations of the vibrational frequencies. Raman spectra indicate that the solid-solid transition at 20 °C is linked to a conformational change to the gauche [Chol] conformer with a concurrent local rearrangement of the anions. The distinct bands of lattice modes in the low-frequency range of the Raman spectra vanish at the 20 °C transition. Given the ease with which metastable crystals can be produced, Raman mappings demonstrate that a sample of [Chol][H2PO4] at ambient temperature can contain a combination of anti- and gauche conformers. Heating to 120 °C causes continuous changes in the local environment of anions rather than melting as suggested by a recent calorimetric investigation of [Chol][H2PO4]. The monotonic change in vibrational spectra is consistent with earlier observations of a very small entropy of fusion and no abrupt jump in the temperature dependence of ionic conductivity along the phase transitions of [Chol][H2PO4].

A Study on the Role of Electric Field in Low-Temperature Plasma Catalytic Ammonia Synthesis via Integrated Density Functional Theory and Microkinetic Modeling

Low-temperature plasma catalysis has shown promise for various chemical processes such as light hydrocarbon conversion, volatile organic compounds removal, and ammonia synthesis. Plasma-catalytic ammonia synthesis has the potential advantages of leveraging renewable energy and distributed manufacturing principles to mitigate the pressing environmental challenges of the energy-intensive Haber-Bosh process, towards sustainable ammonia production. However, lack of foundational understanding of plasma-catalyst interactions poses a key challenge to optimizing plasma-catalytic processes. Recent studies suggest electro- and photoeffects, such as electric field and charge, can play an important role in enhancing surface reactions. These studies mostly rely on using density functional theory (DFT) to investigate surface reactions under these effects. However, integration of DFT with microkinetic modeling in plasma catalysis, which is crucial for establishing a comprehensive understanding of the interplay between the gas-phase chemistry and surface reactions, remains largely unexplored. This paper presents a first-principles framework coupling DFT calculations and microkinetic modeling to investigate the role of electric field on plasma-catalytic ammonia synthesis. The DFT-microkinetic model shows more consistent predictions with experimental observations, as compared to the case wherein the variable effects of plasma process parameters on surface reactions are neglected. In particular, predictions of the DFT-microkinetic model indicate electric field can have a notable effect on surface reactions relative to other process parameters. A global sensitivity analysis is performed to investigate how ammonia synthesis pathways will change in relation to different plasma process parameters. The DFT-microkinetic model is then used in conjunction with active learning to systematically explore the complex parameter space of the plasma-catalytic ammonia synthesis to maximize the amount of produced ammonia while inhibiting reactions dissipating energy, such as the recombination of H2 through gas-phase H radicals and surface-adsorbed H. This paper demonstrates the importance of accounting for the effects of electric field on surface reactions when investigating and optimizing the performance of plasma-catalytic processes.

Decoding Van der Waals Impact on Chirality Transfer in Perovskite Structures: Density Functional Theory Insights

Chiral organic-inorganic perovskites exhibit unique physicochemical properties driven by the symmetry of monovalent organic cations. However, an atomistic understanding of how chiral cations transfer their chirality to the inorganic framework and the role played by van der Waals (vdW) interactions in this process is still incomplete. In this work, we report a theoretical investigation, based on density functional theory calculations within the Perdew-Burke-Ernzerhof (PBE) formulation for the exchange-correlation functional, into the role of the vdW interactions in the chirality transfer process. For that, we selected several vdW corrections, namely, Grimme (D2, D3, D3(BJ)), Tkatchenko-Scheffler (TS, TS+SCS, TS+HSI), density-dependent energy correction (dDsC), and many-body scattering (MBD) energy method correction. For the chiral perovskite systems, we selected a set of chiral organic-inorganic perovskites with several dimensions, namely, from zero-dimensional to three-dimensional, each having enantiomers with R and S configurations. Based on a statistical treatment of the relative errors of all lattice parameters with respect to experimental data, we found that D3, D3(BJ), TS, TS+SCS, TS+HSI, and MBD vdW are the most accurate corrections to describe the equilibrium structural properties of chiral perovskites using the PBE method. We identify chirality-induced sequential asymmetries of distorted octahedrons and propose angular descriptors to quantify them, where the orientations of these distortions depend on the R or S nature of the chiral cations. Furthermore, we demonstrate the importance of accurate vdW interactions in precisely describing these asymmetric distortions. By means of binding energies and charge-transfer analysis, we show that the impact of vdW corrections on the charge distribution leads to a subtle strengthening of hydrogen bonds between chiral cations and inorganic octahedra, resulting in an increase in the binding energy. Finally, we identified that the Rashba-Dresselhaus effect in two-dimensionality is refined by vdW interactions.

Computational understanding of Na-LTA for ethanol-water separation

There is a growing demand for high purity ethanol as an electronic chemical. The conventional distillation process is effective for separating ethanol from water but consumes a significant amount of energy. Selective membrane separation using the LTA-type molecular sieve has been introduced as an alternative. The density functional theory simulation indicates that aluminum (Al) sites are evenly distributed throughout the framework, while sodium (Na+) ions are preferentially located in the six-membered ring. The movement of ethanol molecules can cause Na+ ions to be transported towards the eight-membered ring, hindering the passage of ethanol through the channel. In contrast, the energy barrier for water molecules passing through the channel occupied by Na+ ions is significantly lower, leading to a high level of selectivity for ethanol-water separation.

Investigating the role of dispersion corrections and anharmonic effects on the phase transition in SrZrS3: A systematic analysis from AIMD free energy calculations

A thermally driven needle-like (NL) to distorted perovskite (DP) phase transition in SrZrS3 was investigated by means of ab initio free energy calculations accelerated by machine learning. As a first step, a systematic screening of the methods to include long-range interactions in semilocal density functional theory Perdew-Burke-Ernzerhof calculations was performed. Out of the ten correction schemes tested, the Tkatchenko-Scheffler method with iterative Hirshfeld partitioning method was found to yield the best match between calculated and experimental lattice geometries, while predicting the correct order of stability of NL and DP phases at zero temperature. This method was then used in free energy calculations, performed using several approaches, so as to determine the effect of various anharmonicity contributions, such as the anisotropic thermal lattice expansion or the thermally induced internal structure changes, on the phase transition temperature (TNP→DP). Accounting for the full anharmonicity by combining the NPT molecular dynamics data with thermodynamic integration with harmonic reference provided our best estimate of TNL→DP = 867 K. Although this result is ∼150 K lower than the experimental value, it still provides an improvement by nearly 300 K compared to the previous theoretical report by Koocher et al. [Inorg. Chem. 62, 11134-11141 (2023)].

Nonredox trivalent nickel catalyzing nucleophilic electrooxidation of organics

A thorough comprehension of the mechanism behind organic electrooxidation is crucial for the development of efficient energy conversion technology. Here, we find that trivalent nickel is capable of oxidizing organics through a nucleophilic attack and electron transfer via a nonredox process. This nonredox trivalent nickel exhibits exceptional kinetic efficiency in oxidizing organics that possess the highest occupied molecular orbital energy levels ranging from -7.4 to -6 eV (vs. Vacuum level) and the dual local softness values of nucleophilic atoms in nucleophilic functional groups, such as hydroxyls (methanol, ethanol, benzyl alcohol), carbonyls (formamide, urea, formaldehyde, glucose, and N-acetyl glucosamine), and aminos (benzylamine), ranging from -0.65 to -0.15. The rapid electrooxidation kinetics can be attributed to the isoenergetic channels created by the nucleophilic attack and the nonredox electron transfer via the unoccupied eg orbitals of trivalent nickel (t2g6eg1). Our findings are valuable in identifying kinetically fast organic electrooxidation on nonredox catalysts for efficient energy conversions.

Chemical-Dealloying-Derived PtPdPb-Based Multimetallic Nanoparticles: Dimethyl Ether Electrocatalysis and Fuel Cell Application

In this work, we report a novel multimetallic nanoparticle catalyst composed of Pt, Pd, and Pb and its electrochemical activity toward dimethyl ether (DME) oxidation in liquid electrolyte and polymer electrolyte fuel cells. Chemical dealloying of the catalyst with the lowest platinum-group metal (PGM) content, Pt2PdPb2/C, was conducted using HNO3 to tune the catalyst activity. Comprehensive characterization of the chemical-dealloying-derived catalyst nanoparticles unambiguously showed that the acid treatment removed 50% Pb from the nanoparticles with an insignificant effect on the PGM metals and led to the formation of smaller-sized nanoparticles. Electrochemical studies showed that Pb dissolution led to structural changes in the original catalysts. Chemical-dealloying-derived catalyst nanoparticles made of multiple phases (Pt, Pt3Pb, PtPb) provided one of the highest PGM-normalized power densities of 118 mW mgPGM-1 in a single direct DME fuel cell operated at low anode catalyst loading (1 mgPGM cm-2) at 70 °C. A possible DME oxidation pathway for these multimetallic catalysts was proposed based on an online mass spectrometry study and the analysis of the reaction products.

Tuning the electronic and optical properties of small organic acenedithiophene molecular crystals for photovoltaic applications: First principles calculations

Periodic density functional theory was employed to investigate the impact of chemical modifications on the properties of π-conjugated acenedithiophene molecular crystals. Here, we highlight the importance of the β-methylthionation effect, the position of the sulfur atoms of the thiacycle group and their size, and the number of central benzene rings in the chemical modification strategy. Our results show that the introduction of the methylthio groups at the β-positions of the thiophene and the additional benzene ring at the center of the BDT crystal structure are a promising strategy to improve the performance of organic semiconductors, as observed experimentally. We found that β-MT-ADT exhibits large charge carrier mobility, which is in good agreement with the experimental results and comparable to that of rubrene. In addition, the electronic and optical properties of these ambipolar materials suggest promising performances with β-MT-ADT > ADT >β-MT-NDT > NDT > BEDT-BDT >β-MT-BDT > BDT. Moreover, functionalization with thiacycle-fused sulfur atoms of different sizes and numbers improve the properties of BDT but is still less efficient than the methylthionation effect. Overall, our findings suggest a promising molecular modification strategy for possibly high performance ambipolar organic semiconducting materials.

Atomistic Wear Mechanisms in Diamond: Effects of Surface Orientation, Stress, and Interaction with Adsorbed Molecules

Despite its unrivaled hardness, diamond can be severely worn during the interaction with others, even softer materials. In this work, we calculate from first-principles the energy and forces necessary to induce the atomistic wear of diamond and compare them for different surface orientations and passivation by oxygen, hydrogen, and water fragments. The primary mechanism of wear is identified as the detachment of the carbon chains. This is particularly true for oxidized diamond and diamonds interacting with silica. A very interesting result concerns the role of stress, which reveals that compressive stresses can highly favor wear, making it even energetically favorable.

Nanotribological Properties of Oxidized Diamond/Silica Interfaces: Insights into the Atomistic Mechanisms of Wear and Friction by Ab Initio Molecular Dynamics Simulations

Controlling friction and wear at silica-diamond interfaces is crucial for their relevant applications in tribology such as micro-electromechanical systems and atomic force microscopes. However, the tribological performance on diamond surfaces is highly affected by the working environment where atmospheric gases are present. In this work, we investigate the effects of adsorbed oxygen on the friction and wear of diamond surfaces sliding against silica by massive ab initio molecular dynamics simulations. Different surface orientations, O-coverages, and tribological conditions are considered. The results suggest that diamond surfaces with full oxygen passivation are very effective in preventing surface adhesion, and as a result present extremely low friction and wear. At low oxygen coverage, Si-O-C bond formation was observed as well as atomistic wear initiated from C-C bond breaking at extreme pressure. The analysis of electronic structures of the configurations resulting from key tribochemical reactions clarifies the mechanisms of friction reduction and atomistic wear. Overall, our accurate in silico experiments shed light on the influence of adsorbed oxygen on the tribological properties and wear mechanisms of diamond against silica.

Stacking-Order-Dependent Excitonic Properties Reveal Interlayer Interactions in Bulk ReS2

Rhenium disulfide, a member of the transition metal dichalcogenide family of semiconducting materials, is unique among 2D van der Waals materials due to its anisotropy and, albeit weak, interlayer interactions, confining excitons within single atomic layers and leading to monolayer-like excitonic properties even in bulk crystals. While recent work has established the existence of two stacking modes in bulk, AA and AB, the influence of the different interlayer coupling on the excitonic properties has been poorly explored. Here, we use polarization-dependent optical measurements to elucidate the nature of excitons in AA and AB-stacked rhenium disulfide to obtain insight into the effect of interlayer interactions. We combine polarization-dependent Raman with low-temperature photoluminescence and reflection spectroscopy to show that, while the similar polarization dependence of both stacking orders indicates similar excitonic alignments within the crystal planes, differences in peak width, position, and degree of anisotropy reveal a different degree of interlayer coupling. DFT calculations confirm the very similar band structure of the two stacking orders while revealing a change of the spin-split states at the top of the valence band to possibly underlie their different exciton binding energies. These results suggest that the excitonic properties are largely determined by in-plane interactions, however, strongly modified by the interlayer coupling. These modifications are stronger than those in other 2D semiconductors, making ReS2 an excellent platform for investigating stacking as a tuning parameter for 2D materials. Furthermore, the optical anisotropy makes this material an interesting candidate for polarization-sensitive applications such as photodetectors and polarimetry.

Stability and Semiconducting Behavior of Magnesium Polysulfides by Nonequivalent sp3 Hybridizations

The Mg/S battery has attracted enormous interest in recent years due to its high theoretical capacity, low cost, and high security. However, the understanding of many intermediate magnesium polysulfides in the Mg/S battery remains elusive. Combining extensive structural search and first-principles calculations, we investigate the phase stability, structural character, and electronic structure of magnesium polysulfides in a wide range from MgS to MgS8. The pyrite-type MgS2 (space group: Pa3̅) is predicted to be stable. Five magnesium polysulfides, MgSx (x = 3, 4, 5, 6, and 8), are found to be metastable, with formation enthalpies slightly above the convex hull. S2 dimer, "V"-like S3, and highly distorted Sx chains are found for the polysulfides with bond lengths close to or slightly longer than S8 and bond angles similar to S8. A wide range of band gaps (0.77-2.82 eV) are revealed for the polysulfides due to the contribution of the nonequivalent sp3 hybridization of the S atoms in Sx2-. Our results can help to further understand the electrochemical process in the Mg/S battery.

Influence of Cations on HCOOH and CO Formation during CO2 Reduction on a PdMLPt(111) Electrode

Understanding the role of cations in the electrochemical CO2 reduction (CO2RR) process is of fundamental importance for practical application. In this work, we investigate how cations influence HCOOH and CO formation on PdMLPt(111) in pH 3 electrolytes. While only (a small amount of adsorbed) CO forms on PdMLPt(111) in the absence of metal cations, the onset potential of HCOOH and CO decreases with increasing cation concentrations. The cation effect is stronger on HCOOH formation than that on CO formation on PdMLPt(111). Density functional theory simulations indicate that cations facilitate both hydride formation and CO2 activation by polarizing the electronic density at the surface and stabilizing *CO2-. Although the upshift of the metal work function caused by high coverage of adsorbates limits hydride formation, the cation-induced electric field counterbalances this effect in the case of *H species, sustaining HCOOH production at mild negative potentials. Instead, at the high *CO coverages observed at very negative potentials, surface hydrides do not form, preventing the HCOOH route both in the absence and presence of cations. Our results open the way for a consistent evaluation of cationic electrolyte effects on both activity and selectivity in CO2RR on Pd-Pt catalysts.

The Role of Self-Catalysis Induced by Co Doping in Nonaqueous Li-O2 Batteries

This work systematically studies the product self-catalysis of in situ electrochemical cobalt doping of Li2O2 and reveals its potential mechanism for improving the performance of lithium-oxygen (Li-O2) batteries. Theoretical calculations demonstrate that the discharge products contain substituted and interstitial Co impurities, which serve as active sites to promote the formation of Li3O4 crystallization, thus switching the nucleation mechanism from the main discharge product Li2O2 to Li3O4. This Co-doping behavior leads to the thermodynamically favorable and dynamically stable formation of Li3O4 crystals during the discharge process. Through systematic investigation of the structural, energetic, electronic, diffusive, and catalytic properties of the Co-doped Li2O2 and Li3O4 compounds, we found that Li3O4 has better charge/mass transport and a lower overpotential for the Li3O4 formation/decomposition reaction. Consequently, this work elucidates that Co doping provides a simple and effective approach for increasing the proportion of Li3O4, which can significantly improve the Li-O2 battery performance.

Toward functionalization of ZnO nanotubes and monolayers with 5-aminolevulinic acid drugs as possible nanocarriers for drug delivery: a DFT based molecular dynamic simulation

We have investigated the interactions between a 5-aminolevulinic acid (ALA) drug and ZnO nanostructures including ZnO monolayers and ZnO nanotubes (ZnONTs) using density functional theory (DFT) calculations. In the context of the dispersion corrected Perdew-Burke-Ernzerhof (PBE) approach, the energetics, charge transfer, electronic structure and equilibrium geometries have been estimated. As ALA is adsorbed onto/into the ZnONTs and on the ZnO monolayer with interaction energies (Eint) of -2.55/-2.75 eV and -2.51 eV, respectively, the calculated Eint values and bonding distances (∼2 Å) reveal that the interaction type is chemisorption. The ZnO nanostructures showed promising performance in the ALA drug functionalization, taking into account the interaction energy values. The band gap almost remains unchanged for both of the substrates under consideration after ALA adsorption, and the semiconductor properties of the substrates are preserved, according to the analyzed density of states (DOSs) spectra. The interaction nature of the ALA-ZnO nanostructures according to the atom in molecule (AIM) analysis was found to be polar attraction with partial covalent bonding between O and Zn. Our DFT based molecular dynamic (MD) simulation results demonstrate that, in the aqueous solution, ALA moves toward the interior sidewall of the ZnONTs and ZnO nanosheet surface and binds to the Zn atom through its O (carbonyl/hydroxyl groups) and N atoms and the hydroxyl H atom was dissociated and binds to the O atom of the ZnO surface. However, in the case of ALA adsorption onto the outer surface of ZnONTs, only the O atoms of carbonyl groups bind to the Zn atom and the structure of the drug remains undestroyed during the adsorption. The current findings shed light on the polar drug adsorption/encapsulation behavior on/into ZnO nanostructures, which may encourage further use of ZnO-based nanomaterials in the field of drug delivery and bio-functionalized nanomaterials.

How the Temperature and Composition Govern the Structure and Band Gap of Zr-Based Chalcogenide Perovskites: Insights from ML Accelerated AIMD

The effects of temperature and composition on the structural and electronic properties of chalcogenide perovskite (CP) materials AZrX3 (A = Ba, Sr, Ca; X = S, Se) in the distorted perovskite (DP) phase are investigated using ab initio molecular dynamics (AIMD) accelerated by machine-learned force fields. Long-range van der Waals (vdW) interactions, incorporated into the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional using the DFT-D3 scheme, are found to be crucial for achieving correct predictions of structural parameters. Our calculations show that the distortion of the DP structure with respect to the parent cubic (C) phase, realized in the form of interoctahedral tilting, decreases with the increasing size of the A cations. The tendency for a gradual transformation of the DP-to-C phase with increasing temperature is shown to be strongly composition-dependent. The transformation temperature decreases with the size of cation A and increases with the size of anion X. Thus, within the range of the temperatures considered here (300-1200 K), a complete transformation is observed only for BaZrS3 (∼600 K) and BaZrSe3 (∼900 K). The computed band gap of CPs is shown to monotonically decrease with increasing temperature, and the magnitude of this decrease is found to be proportional to the extent of the thermally induced changes in the internal structure. Diverse factors affecting the magnitude of band gaps of CP materials are analyzed.

Electrosorption, Desorption, and Oxidation of Perfluoroalkyl Carboxylic Acids (PFCAs) via MXene-Based Electrocatalytic Membranes

MXenes exhibit excellent conductivity, tunable surface chemistry, and high surface area. Particularly, the surface reactivity of MXenes strongly depends on surface exposed atoms or terminated groups. This study examines three types of MXenes with oxygen, fluorine, and chlorine as respective terminal atoms and evaluates their electrosorption, desorption, and oxidative properties. Two perfluorocarboxylic acids (PFCAs), perfluorobutanoic acid (PFBA) and perfluorooctanoic acid (PFOA) are used as model persistent micropollutants for the tests. The experimental results reveal that O-terminated MXene achieves a significantly higher adsorption capacity of 215.9 mg·g^-1 and an oxidation rate constant of 3.9 × 10^-2 min^-1 for PFOA compared to those with F and Cl terminations. Electrochemical oxidation of the two PFCAs (1 ppm) with an applied potential of +6 V in a 0.1 M Na₂SO₄ solution yields >99% removal in 3 h. Moreover, PFOA degrades about 20% faster than PFBA on O-terminated MXene. The density functional theory (DFT) calculations reveal that the O-terminated MXene surface yielded the highest PFOA and PFBA adsorption energy and the most favorable degradation pathway, suggesting the high potential of MXenes as highly reactive and adsorptive electrocatalysts for environmental remediation.

Structural, elastic, electronic, optical and vibrational properties of single-layered, bilayered and bulk molybdenite MoS2-2H

In recent years, transition metal dichalcogenides have received great attention since they can be prepared as two-dimensional semiconductors, presenting heterodesmic structures incorporating strong in-plane covalent bonds and weak out-of-plane interactions, with an easy cleavage/exfoliation in single or multiple layers. In this context, molybdenite, the mineralogical name of molybdenum disulfide, MoS₂, has drawn much attention because of its very promising physical properties for optoelectronic applications, in particular a band gap that can be tailored with the material's thickness, optical absorption in the visible region and strong light-matter interactions due to the planar exciton confinement effect. Despite this wide interest and the numerous experimental and theoretical articles in the literature, these report on just one or two specific features of bulk and layered MoS₂ and sometimes provide conflicting results. For these reasons, presented here is a thorough theoretical analysis of the different aspects of bulk, monolayer and bilayer MoS₂ within the density functional theory (DFT) framework and with the DFT-D3 correction to account for long-range interactions. The crystal chemistry, stiffness, and electronic, dielectric/optical and phonon properties of single-layered, bilayered and bulk molybdenite have been investigated, to obtain a consistent and detailed set of data and to assess the variations and cross correlation from the bulk to single- and double-layer units. The simulations show the indirect-direct transition of the band gap (K-K' in the first Brillouin zone) from the bulk to the single-layer structure, which however reverts to an indirect transition when a bilayer is considered. In general, the optical properties are in good agreement with previous experimental measurements using spectroscopic ellipsometry and reflectivity, and with preliminary theoretical simulations.

Structural, Electronic and Optical Properties of Some New Trilayer Van de Waals Heterostructures

Constructing two-dimensional (2D) van der Waals (vdW) heterostructures is an effective strategy for tuning and improving the characters of 2D-material-based devices. Four trilayer vdW heterostructures, BP/BP/MoS₂, BlueP/BlueP/MoS₂, BP/graphene/MoS₂ and BlueP/graphene/MoS₂, were designed and simulated using the first-principles calculation. Structural stabilities were confirmed for all these heterostructures, indicating their feasibility in fabrication. BP/BP/MoS₂ and BlueP/BlueP/MoS₂ lowered the bandgaps further, making them suitable for a greater range of applications, with respect to the bilayers BP/MoS₂ and BlueP/MoS₂, respectively. Their absorption coefficients were remarkably improved in a wide spectrum, suggesting the better performance of photodetectors working in a wide spectrum from mid-wave (short-wave) infrared to violet. In contrast, the bandgaps in BP/graphene/MoS₂ and BlueP/graphene/MoS₂ were mostly enlarged, with a specific opening of the graphene bandgap in BP/graphene/MoS₂, 0.051 eV, which is much larger than usual and beneficial for optoelectronic applications. Accompanying these bandgap increases, BP/graphene/MoS₂ and BlueP/graphene/MoS₂ exhibit absorption enhancement in the whole infrared, visible to deep ultraviolet or solar blind ultraviolet ranges, implying that these asymmetrically graphene-sandwiched heterostructures are more suitable as graphene-based 2D optoelectronic devices. The proposed 2D trilayer vdW heterostructures are prospective new optoelectronic devices, possessing higher performance than currently available devices.

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. Our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

Binuclear spin-crossover [Fe(bt)(NCS)2]2(bpm) complex: A study using first principles calculations

The spin-crossover [Fe(bt)(NCS)₂]₂(bpm) complex is studied using spin-polarized density functional theory within the generalized gradient approximation, the Hubbard U and the weak van der Waals interactions in conjunction with the projector augmented wave method in its molecular and periodic arrangements. It is shown that the considered complex has three magnetic configurations [high spin state (HS)-HS, HS-low spin state (LS), and LS-LS] corresponding to those observed experimentally after two transition temperatures T_c ⁽¹⁾ of 163 K and T_c ⁽²⁾ of 197 K. For the HS-HS magnetic state, we found that the two Fe centers are antiferromagnetically coupled for both molecular and periodic structures in good agreement with the experimental observations. Our results show that the computed total energy difference between the magnetic state configurations of the considered Fe₂ complex is significantly smaller compared to those reported in the literature for other mono- or binuclear compounds.

Anharmonic Correction to Free Energy Barriers from DFT-Based Molecular Dynamics Using Constrained Thermodynamic Integration

For the calculation of anharmonic contributions to free energy barriers, constrained thermodynamic λ-path integration (λ-TI) from a harmonic reference force field to density functional theory is presented as an alternative to the established Blue Moon ensemble method (ξ-TI), in which free energy gradients along the reaction coordinate ξ are integrated. With good agreement in all cases, the λ-TI method is benchmarked against the ξ-TI method for several reactions, including the internal CH₃ group rotation in ethane, a nucleophilic substitution of CH₃Cl, a retro-Diels-Alder reaction, and a proton transfer in zeolite H-SSZ-13. An advantage of λ-TI is that one can use virtually any reference state to compute anharmonic contributions to reaction free energies or free energy barriers. This is particularly relevant for catalysis, where it is now possible to compute anharmonic corrections to the free energy of a transition state relative to any reference, for example, the most stable state of the active site and the reactants in the gas phase. This is in contrast to ξ-TI, where free energy barriers can only be computed relative to an initial state with all reactants coadsorbed. Finally, the Bennett acceptance ratio method combined with λ-TI is demonstrated to reduce the number of required integration grid points with tolerable accuracy, favoring thus λ-TI over ξ-TI in terms of computational efficiency.

Highly Stable Garnet Fe2 Mo3 O12 Cathode Boosts the Lithium-Air Battery Performance Featuring a Polyhedral Framework and Cationic Vacancy Concentrated Surface

Lithium-air batteries (LABs), owing to their ultrahigh theoretical energy density, are recognized as one of the next-generation energy storage techniques. However, it remains a tricky problem to find highly active cathode catalyst operating within ambient air. In this contribution, a highly active Fe₂ Mo₃ O₁₂ (FeMoO) garnet cathode catalyst for LABs is reported. The experimental and theoretical analysis demonstrate that the highly stable polyhedral framework, composed of FeO octahedrons and MO tetrahedrons, provides a highly effective air catalytic activity and long-term stability, and meanwhile keeps good structural stability. The FeMoO electrode delivers a cycle life of over 1800 h by applying a simple half-sealed condition in ambient air. It is found that surface-rich Fe vacancy can act as an O₂ pump to accelerate the catalytic reaction. Furthermore, the FeMoO catalyst exhibits a superior catalytic capability for the decomposition of Li₂ CO₃ . H₂ O in the air can be regarded as the main contribution to the anode corrosion and the deterioration of LAB cells could be attributed to the formation of LiOH·H₂ O at the end of cycling. The present work provides in-depth insights to understand the catalytic mechanism in air and constitutes a conceptual breakthrough in catalyst design for efficient cell structure in practical LABs.

Synergistic promotions between CO2 capture and in-situ conversion on Ni-CaO composite catalyst

The integrated CO₂ capture and conversion (iCCC) technology has been booming as a promising cost-effective approach for Carbon Neutrality. However, the lack of the long-sought molecular consensus about the synergistic effect between the adsorption and in-situ catalytic reaction hinders its development. Herein, we illustrate the synergistic promotions between CO₂ capture and in-situ conversion through constructing the consecutive high-temperature Calcium-looping and dry reforming of methane processes. With systematic experimental measurements and density functional theory calculations, we reveal that the pathways of the reduction of carbonate and the dehydrogenation of CH₄ can be interactively facilitated by the participation of the intermediates produced in each process on the supported Ni-CaO composite catalyst. Specifically, the adsorptive/catalytic interface, which is controlled by balancing the loading density and size of Ni nanoparticles on porous CaO, plays an essential role in the ultra-high CO₂ and CH₄ conversions of 96.5% and 96.0% at 650 °C, respectively.

Bi-C monolayer as a promising 2D anode material for Li, Na, and K-ion batteries

Electrode materials with high electrochemical efficiency are required for battery technology that can be used to store renewable energy. Bismuth (Bi) has shown great potential as an electrode material for metal ion batteries due to its large volumetric capacity and reasonable operating potential. However, the cycling performance deteriorates due to the drastic volume changes that occur during alloying and dealloying. Herein, we design a 2D Bi-C metal sheet using density functional theory and investigate the feasibility of this nanosheet for alkali metal ion batteries. The predicted metallic Bi-C monolayer (ML) are highly stable and show sound electrode performance. Moreover, alkali metal atoms exhibit high diffusivities on both sides (Bi and C sides) with low energy barriers of 0.252/0.201, 0.217/0.169, and 0.179/0.136 eV for Li, Na, and K ions, respectively. Furthermore, the Bi-C ML shows high theoretical storage capacities of (485 mA h g^-1) for Li and Na and (364 mA h g^-1) for K and low open-circuit voltage of 0.12, 0.24, and 0.32 V for Li, Na, and K ions, respectively. These exciting findings show that the predicted Bi-C ML can be used as an anode material for Li-, Na- and K-ion batteries.

Theoretical exploration of the structural, electronic and optical properties of g-C3N4/C3N heterostructures

Integration of graphene-like carbon nitride materials is essential for nanoelectronic applications. Using density-functional theory (DFT), we systematically investigate the structural, electronic and optical properties of a s-triazine-based g-C₃N₄/C₃N heterostructure under different modified conditions. The g-C₃N₄/C₃N van der Waals heterostructure (vdWH) formed has an indirect bandgap with type-II band alignment and the band structures can be tuned from type-II band alignment to type-I band alignment by applying biaxial strains and external electric fields (E_field). Compared to single transition metal (TM) atoms at g-C₃N₄/C₃N surfaces, the TM atoms anchored in the interlayer region exhibit more stability, and the corresponding bandgaps are changed from 0.19 eV to 0.61 eV. In addition, the g-C₃N₄/C₃N heterostructure has a strong absorption coefficient in the ultraviolet-visible light region along the x direction. It is found that compressive strain has a large influence on the absorption coefficient of the g-C₃N₄/C₃N system. With the increased compressive strain, the absorption spectra in the visible light region disappeared. Tensile strain has a slight effect on the absorption range, but causes a red shift of the absorption spectrum. In comparison, the light absorption coefficient of the g-C₃N₄/C₃N system remains almost unchanged under the E_field conditions. In summary, the formation of a s-triazine-based g-C₃N₄/C₃N heterostructure has shown potential for applications in nanoelectronic and optoelectronic devices.

Effect of 3d Transition Metal Atom Intercalation Concentration on the Electronic and Magnetic Properties of Graphene/MoS2 Heterostructure: A First-Principles Study

The electronic and magnetic properties of graphene/MoS₂ heterostructures intercalated with 3d transition metal (TM) atoms at different concentrations have been systematically investigated by first principles calculations. The results showed that all the studied systems are thermodynamically stable with large binding energies of about 3.72 eV-6.86 eV. Interestingly, all the TM-intercalated graphene/MoS₂ heterostructures are ferromagnetic and their total magnetic moments increase with TM concentration. Furthermore, TM concentration-dependent spin polarization is obtained for the graphene layer and MoS₂ layer due to the charge transfer between TM atoms and the layers. A significant band gap is opened for graphene in these TM-intercalated graphene/MoS₂ heterostructures (around 0.094 eV-0.37 eV). With the TM concentration increasing, the band gap of graphene is reduced due to the enhanced spin polarization of graphene. Our study suggests a research direction for the manipulation of the properties of 2D materials through control of the intercalation concentration of TM atoms.

Influence of strain on an ultrafast phase transition

The flexibility of 2D materials combined with properties highly sensitive to strain makes strain engineering a promising avenue for manipulation of both structure and function. Here we investigate the influence of strain, associated with microstructural defects, on a photo-induced structural phase transition in Td-WTe₂. Above threshold photoexcitation of uniform, non-strained, samples result in an orthorhombic Td to a metastable orthorhombic 1T* phase transition facilitated by shear displacements of the WTe₂ layers along the b axis of the material. In samples prepared with wrinkle defects WTe₂ continue its trajectory through a secondary transition that shears the unit cell along the c axis towards a metastable monoclinic 1T' phase. The time scales and microstructural evolution associated with the transition and its subsequent recovery to the 1T* phase is followed in detail by a combination of ultrafast electron diffraction and microscopy. Our findings show how local strain fields can be employed for tailoring phase change dynamics in ultrafast optically driven processes with potential applications in phase change devices.

High-Throughput Computational Screening of Two-Dimensional Semiconductors

Two-dimensional (2D) materials have attracted great attention mainly due to their unique physical properties and ability to fulfill the demands of future nanoscale devices. By performing high-throughput first-principles calculations combined with a semiempirical van der Waals dispersion correction, we have screened 73 direct- and 183 indirect-gap 2D nonmagnetic semiconductors from nearly 1000 monolayers according to the criteria for thermodynamic, mechanical, dynamic, and thermal stabilities and conductivity type. We present the calculated lattice constants, formation energy, Young's modulus, Poisson's ratio, shear modulus, anisotropic effective mass, band structure, band gap, ionization energy, electron affinity, and simulated scanning tunnel microscopy for each candidate meeting our criteria. The resulting 2D semiconductor database (2DSdb) can be accessed via the Web site https://materialsdb.cn/2dsdb/index.html. The 2DSdb provides an ideal platform for computational modeling and design of new 2D semiconductors and heterostructures in photocatalysis, nanoscale devices, and other applications. Further, a linear fitting model was proposed to evaluate band gap, ionization energy, and electron affinity of 2D semiconductors from the density functional theory (DFT) calculated data as initial input. This model can be as precise as hybrid DFT but with much lower computational cost.

Selective sensing properties and enhanced ferromagnetism in CrI3monolayer via gas adsorption

Recent fabrication of chromium triiodide (CrI₃) monolayers has raised potential prospects of developing two-dimensional (2D) ferromagnetic materials for spintronic device applications. The low Curie temperature has stimulated further interest for improving the ferromagnetic stability of CrI₃monolayer. Here, based on density functional theory calculations, we investigated the adsorption energy, charge transfer, electronic and magnetic properties of gases (CO, CO₂, N₂, NH₃, NO, NO₂, O₂, and SO₂) adsorption on the CrI₃monolayer. It is found that CrI₃is sensitive to the NH₃, NO, and NO₂adsorption due to the high adsorption energy and large charge transfer. The electrical transport results show that the conductivity of CrI₃monolayer is significantly reduced with the adsorption of N-based gases, suggesting that CrI₃exhibits superior sensitivity and selectivity toward N-based gases. In addition, the ferromagnetic stability and Curie temperature (T_C) of CrI₃monolayer can be effectively enhanced by the adsorption of magnetic gases (NO, NO₂, O₂). This work not only demonstrates that CrI₃monolayer can be used as a promising candidate for gas sensing, but also brings further interest to tune the electronic and magnetic properties of 2D ferromagnetic materials via gas adsorption.

Effect of temperature on CO oxidation over Pt(111) in two-dimensional confinement

Confined catalysis between a two-dimensional (2D) cover and metal surfaces has provided a unique environment with enhanced activity compared to uncovered metal surfaces. Within this 2D confinement, weakened adsorption and lowered activation energies were observed using surface science experiments and density functional theory (DFT) calculations. Computationally, the role of electronic and mechanical factors responsible for the improved activity was deduced only from static DFT calculations. This demands a detailed investigation on the dynamics of reactions under 2D confinement, including temperature effects. In this work, we study CO oxidation on a 2D graphene covered Pt(111) surface at 90 and 593 K using DFT-based ab initio molecular dynamics simulations starting from the transition state configuration. We show that CO oxidation in the presence of a graphene cover is substantially enhanced (2.3 times) at 90 K. Our findings suggest that 2D confined spaces can be used to enhance the activity of chemical reactions, especially at low temperatures.

Investigating the role of structural water on the electrochemical properties of α-V2O5 through density functional theory

The α polymorph of V₂O₅ is one of the few known cathodes capable of reversibly intercalating multivalent ions such as Mg, Ca, Zn and Al, but suffers from sluggish diffusion kinetics. The role of H₂O within the electrolyte and between the layers of the structure in the form of a xerogel/aerogel structure, though, has been shown to lower diffusion barriers and lead to other improved electrochemical properties. This density functional theory study systematically investigates how and why the presence of structural H₂O within α-V₂O₅ changes the resulting structure, voltage, and diffusion kinetics for the intercalation of Li, Na, Mg, Ca, Zn, and Al. We found that the coordination of H₂O molecules with the ion leads to an improvement in voltage and energy density for all ions. This voltage increase was attributed to the extra host sites for electrons present with H₂O, thus leading to a stronger ionization of the ion and a higher voltage. We also found that the increase in interlayer distance and a potential "charge shielding" effect drastically changes the electrostatic environment and the resulting diffusion kinetics. For Mg and Ca, this resulted in a decrease in diffusion barrier from 1.3 eV and 2.0 eV to 0.89 eV and 0.4 eV, respectively. We hope that our study motivates similar research regarding the role of water in both V₂O₅ xerogels/aerogels and other layered transition metal oxides.

Lithium stabilizes square-two-dimensional metal sheets: a computational exploration

Based on the M₄-square-containing M₄Li₂ (M = Al, Ga, In, Tl, Ge, Sn, Pb, Sb, Bi, Cu, Ag, Au, and Hg) clusters, we computationally designed two-dimensional (2D) M₂Li sheets consisting of M₄-square motifs. The four M₂Li-I (M = Sb, Bi, Ag, and Au) monolayers with Li square sublayer sandwiched between two M square sublayers (P4/mmm space group) were confirmed to be stable (high cohesive energies, positive vibrational frequencies, moderate Young's moduli, and structural integrity during first-principles molecular dynamics simulations at 500 K), and the particle swarm optimization (PSO) method identified these constructed monolayers as the global minima in the 2D space. The three M₂Li-I (M = Sb, Bi, and Ag) monolayers demonstrated a half-auxetic behavior. Ag₂Li-I could well activate CO₂ and convert it into HCOOH by following the path * → *CO₂ → *OCHO → *HCOOH → *+HCOOH. Particularly, Ag₂Li-I shows great promise as an electrocatalyst for CO₂ reduction as its limiting potential is as low as 0.40 (0.27) V without (with) considering the solvent effect. Our theoretical explorations reveal that lithium can stabilize the square metal monolayers, and the stable square binary metal sheets exhibit diverse mechanical and electrochemical properties, which can be used in the fields of mechanics and electrochemical catalysis.

Vanadium Carbide (V4C3) MXene as an Efficient Anode for Li-Ion and Na-Ion Batteries

Li-ion batteries (LIBs) and Na-ion batteries (SIBs) are deemed green and efficient electrochemical energy storage and generation devices; meanwhile, acquiring a competent anode remains a serious challenge. Herein, the density-functional theory (DFT) was employed to investigate the performance of V₄C₃ MXene as an anode for LIBs and SIBs. The results predict the outstanding electrical conductivity when Li/Na is loaded on V₄C₃. Both Li_2xV₄C₃ and Na_2xV₄C₃ (x = 0.125, 0.5, 1, 1.5, and 2) showed expected low-average open-circuit voltages of 0.38 V and 0.14 V, respectively, along with a good Li/Na storage capacity of (223 mAhg^-1) and a good cycling performance. Furthermore, there was a low diffusion barrier of 0.048 eV for Li_0.0625V₄C₃ and 0.023 eV for Na_0.0625V₄C₃, implying the prompt intercalation/extraction of Li/Na. Based on the findings of the current study, V₄C₃-based materials may be utilized as an anode for Li/Na-ion batteries in future applications.

Electrocatalytic activity of a β-Sb two-dimensional surface for the hydrogen evolution reaction

Hydrogen energy is considered to be one of the most promising clean energy sources. The development of highly active, low-cost catalysts, and good stability is essential for hydrogen production. Herein, the catalytic activity of a two-dimensional β-Sb surface doped with main-group elements (N, P, As, O, S, Se, and Te) for the hydrogen evolution reaction (HER) was investigated by density functional theory, and the catalytic activity of the β-Sb monolayer can be improved by doping group VIA atoms. The catalytic activity of Se@Sb and O@Sb structures at the doping concentration of 2.78% and the S@Sb structure at the doping concentration of 5.56% may be as good as the Pt(111) surface, while keeping energetically stable. In addition, the catalytic performance could be optimized under biaxial strain. Further analysis suggests that the activity is caused by hole states in the lone pair electrons, which are created by the group VIA atom dopants. And our work also reveals that the density of states at the Fermi level could be an appropriate descriptor of the hydrogenation Gibbs free energy. This work not only proposes a novel non-platinum HER catalyst but also provides physical foundations for further application on antimonene-based catalysts.

Performance of Dispersion-Inclusive Density Functional Theory Methods for Energetic Materials

Molecular crystals of energetic materials (EMs) are denser than typical molecular crystals and are characterized by distinct intermolecular interactions between nitrogen-containing moieties. To assess the performance of dispersion-inclusive density functional theory (DFT) methods, we have compiled a data set of experimental sublimation enthalpies of 31 energetic materials. We evaluate the performance of three methods: the semilocal Perdew-Burke-Ernzerhof (PBE) functional coupled with the pairwise Tkatchenko-Scheffler (TS) dispersion correction, PBE with the many-body dispersion (MBD) method, and the PBE-based hybrid functional (PBE0) with MBD. Zero-point energy contributions and thermal effects are described using the quasi-harmonic approximation (QHA), including explicit treatment of thermal expansion, which we find to be non-negligible for EMs. The lattice energies obtained with PBE0+MBD are the closest to experimental sublimation enthalpies with a mean absolute error of 9.89 kJ/mol. However, the state-of-the-art treatment of vibrational and thermal contributions makes the agreement with experiment worse. Pressure-volume curves are also examined for six representative materials. For pressure-volume curves, all three methods provide reasonable agreement with experimental data with mean absolute relative errors of 3% or less. Most of the intermolecular interactions typical of EMs, namely nitro-amine, nitro-nitro, and nitro-hydrogen interactions, are more sensitive to the choice of the dispersion method than to the choice of the exchange-correlation functional. The exception is π-π stacking interactions, which are also very sensitive to the choice of the functional. Overall, we find that PBE+TS, PBE+MBD, and PBE0+MBD do not perform as well for energetic materials as previously reported for other classes of molecular crystals. This highlights the importance of testing dispersion-inclusive DFT methods for diverse classes of materials and the need for further method development.

Defect Healing of MAPbI3 Perovskite Single Crystal Surface by Benzylamine

Controlling the surface traps in metal halide perovskites (MHPs) is essential for device performance, stability, and commercialization. Here, a facile approach is introduced to passivate the methylammonium lead iodide (MAPbI3) perovskite single crystal (PSC) surface defects by benzylamine (BA) ligand treatment, and the natural crystallographic (100) facets surface of PSC is chosen as the research platform to provide a deeper understanding of the passivation process. The confocal photoluminescence (PL) results show that the pristine three-dimensional (3D) MAPbI3 PSC surface with a symmetric emission spectrum is normally converted to a pure two-dimensional (2D) BA2PbI4, and also forms a quasi-2D Ruddlesden–Popper perovskite (RPP) BA2MAn−1PbnI3n+1 (n = 2, 3, 4, … ∞) after BA exchange with cation defects. The blue shift in the PL peak, as well as the extended exciton lifetimes of time-resolved photoluminescence (TRPL), indicate the realization of surface defect passivation. Additionally, changes in surface morphology are also investigated. The reaction starts with the formation of small, layered crystallites over the surface; as time elapses, the layered crystallites spread and merge in contact with each other, eventually resulting in smooth features. Our findings present a simple approach for MAPbI3 PSC surface defect passivation, which aims to advance MHP optimization processes toward practical perovskite device applications.

Metal dimers embedded vertically in defect-graphene as gas sensors: a first-principles study

A highly symmetric structure of metal dimers embedded in defect-graphene (M₂⊥gra) in a perpendicular manner was designed. Five M₂⊥gra (M = Co, Ni, Rh, Ir and Pt) monolayers were identified to be stable by density functional theory (DFT) calculations. To investigate the capability of those new structures as gas sensors, the adsorption behavior of ten gas molecules (O₂, N₂, CO, CO₂, NO, NO₂, NH₃, H₂O, H₂S and SO₂) on M₂⊥gra was explored, and the charge transfer, magnetism changes, etc. of these adsorption systems were analyzed. The Ni₂⊥gra can be used as a gas sensor for O₂ at 500 K by the analysis of electronic and magnetic properties. At room temperature, the Pt₂⊥gra is expected to be an excellent CO₂ gas selector due to its high selectivity, sensitivity and short recovery time (1.04 × 10^-12 s). The electronic and magnetic coupling between the metal atoms in the vertical metal dimers plays an important role in sensing gas molecules. Our work paves a new way to design metal-dimer-based nanomaterials.

Understanding the CH4 Conversion over Metal Dimers from First Principles

Inspired by the advantages of bi-atom catalysts and recent exciting progresses of nanozymes, by means of density functional theory (DFT) computations, we explored the potential of metal dimers embedded in phthalocyanine monolayers (M₂-Pc), which mimics the binuclear centers of methane monooxygenase, as catalysts for methane conversion using H₂O₂ as an oxidant. In total, 26 transition metal (from group IB to VIIIB) and four main group metal (M = Al, Ga, Sn and Bi) dimers were considered, and two methane conversion routes, namely *O-assisted and *OH-assisted mechanisms were systematically studied. The results show that methane conversion proceeds via an *OH-assisted mechanism on the Ti₂-Pc, Zr₂-Pc and Ta₂-Pc, a combination of *O- and *OH-assisted mechanism on the surface of Sc₂-Pc, respectively. Our theoretical work may provide impetus to developing new catalysts for methane conversion and help stimulate further studies on metal dimer catalysts for other catalytic reactions.

Synergy of cobalt vacancies and iron doping in cobalt selenide to promote oxygen electrode reactions in lithium-oxygen batteries

Electronic structural engineering plays a key role in the design of high-efficiency catalysts. Here, to achieve optimal electronic states, introduction of exotic Fe dopant and Co vacancy into CoSe₂ nanosheet (denoted as Fe-CoSe₂-V_Co) is presented. The obtained Fe-CoSe₂-V_Co demonstrates excellent catalytic activity as compared to CoSe₂. Experimental results and density functional theory (DFT) calculations confirm that Fe dopant and Co defects cause significant electron delocalization, which reduces the adsorption energy of LiO₂ intermediate on the catalyst surface, thereby obviously improving the electrocatalytic activity of Fe-CoSe₂-V_Co towards oxygen redox reactions. Moreover, the synergistic effect between Co vacancy and Fe dopant is able to optimize the microscopic electronic structure of Co ion, further reducing the energy barrier of oxygen electrode reactions on Fe-CoSe₂-V_Co. And the lithium-oxygen batteries (LOBs) based on Fe-CoSe₂-V_Co electrodes demonstrate a high Coulombic efficiency (CE) of about 72.66%, a large discharge capacity of about 13723 mA h g^-1, and an excellent cycling life of about 1338 h. In general, the electronic structure modulation strategy with the reasonable introduction of vacancy and dopant is expected to inspire the design of highly efficient catalysts for various electrochemical systems.