Cohesive energy of 3d transition metals: Density functional theory atomic and bulk calculations

DFT investigation of geometrical, vibrational, elastic, electronic, optical, and thermoelectric properties of aluminum pnictogens compounds

The aim of this study is to investigate the geometrical, vibrational, elastic, electronic, optical, and thermoelectric characteristics of aluminum pnictides in monolayer square lattice and bilayer hexagonal phases (s- and h-AlX; X = N, P, As) using first principles. The s- and h-AlX materials are mechanically, energetically, and dynamically stable, through phonon dispersion and elastic properties investigations. It was observed that s-AlX materials exhibited both direct and indirect bandgaps, whereas h-AlX materials exhibited indirect bandgap behavior. The energy bandgap values for s- and h-AlX materials measured between 0.79 eV and 3.49 eV for the PBE functional, and between 1.49 eV and 4.74 eV for the HSE06 functional. The effective mass, mobility and relaxation time of electron carriers as well as hole carriers from the band structure of s- and h-AlX are examined to gain a better perception into these materials. The AlP monolayer square lattice phase has the highest mobility and relaxation time of 266129.60 cm2V-1s-1 and 740369.83 fs among entire s- and h-AlX materials. The optical characteristics of s- and h-AlX materials are examined in the existence of field polarizations. The thermoelectric properties of the AlX materials are assessed for temperature dependent. Our investigated results expose that AlP/AlP and AlAs/AlAs are the proficient thermoelectric materials at room temperature in the considered sequence. The present investigation shows that the s- and h-AlX materials are mostly active in the UV region of electromagnetic spectrum, and may find applications in UV-photodetectors and UV-protectant materials.

Design a Functional Graphene with Decoration of Dual Transition Metal Dopants for Hydrogen Evolution Electrocatalysis

Since hydrogen is a promising alternative to fossil fuels due to its high energy density and environmental friendliness, water electrolysis for hydrogen production has received widespread attentions wherein the development of active and stable catalytic materials is a key research direction. This article designs a dual transition metal doped functional graphene for hydrogen evolution reaction via density functional theory calculations. Among varied combinations, 16 candidates are screened out that are expected to be stable as reflected by the criterion of formation energy Ef<0 and active due to its free energy of hydrogen adsorption ▵GH within the window of ±0.3 eV. Considering its feasibility in structural modification and electronic adjustment due to the strong dd orbital couplings, the homogeneous dual-atom moiety delivers improved performance toward hydrogen evolution in comparison with the single-atom counterpart. Owing to the good resistance of electrochemical dissolution, the work figures out the potential combinations of Cu2C3N3, Rh2C6, Rh2C3N3 and Rh2N6 endowed with the ▵GH values of -0.03, 0.12, -0.21, and 0.06 eV, respectively, being comparable to the benchmark Pt materials. Therefore, this study provides a new direction for the experimental synthesis of highly active carbon-based electrocatalysts and highlights the well-tuning ability posed by the dual-atom interaction.

Systematic investigation of structural, magneto-electronic, mechanical, thermophysical, optical and thermoelectric properties of Hf2VZ (Z = Ga, In, Tl) inverse Heusler alloy for spintronics applications

The structural stability, magneto-electronic, mechanical, thermodynamic, thermoelectric and optical, characteristics of the Hf2VZ (Z = Ga, In, Tl) Heusler alloy are revealed and understood by a comprehensive investigation employing density functional theory simulations. The stability of these alloys in F-43m phase is confirmed by structural optimizations and cohesive energies, which also provide the equilibrium lattice parameters. Compared to generalized gradient approximations, modified Becke-Johnson methods were more effective in determining the electrical structure and ground state attributes. Hf2VZ (Z = Ga, In, Tl) is predicted to have half-metallic ferromagnetic characteristics with indirect spin-up gaps based on the band structure analysis and density of state calculations. Stability of these compounds is determined by calculating the elastic constants indicating the ductile nature of these alloys. The quasi-harmonic Debye model is used to predict the effects of temperature and pressure on thermodynamic characteristics, conveying the alloys' thermodynamic stability. To estimate the thermoelectric performance of these materials, we compute electrical conductivities and Seebeck coefficients. The optical parameters like absorption coefficient, optical conductivity, dielectric constants etc., were determined to show the photo-voltaic applications of these alloys. Hence, the finding will lead to future research on developing new types of Hafnium based Heusler alloys for spintronics, thermoelectric and optoelectronics 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.

Stable hydrogen evolution reaction at high current densities via designing the Ni single atoms and Ru nanoparticles linked by carbon bridges

Continuous and effective hydrogen evolution under high current densities remains a challenge for water electrolysis owing to the rapid performance degradation under continuous large-current operation. In this study, theoretical calculations, operando Raman spectroscopy, and CO stripping experiments confirm that Ru nanocrystals have a high resistance against deactivation because of the synergistic adsorption of OH intermediates (OHad) on the Ru and single atoms. Based on this conceptual model, we design the Ni single atoms modifying ultra-small Ru nanoparticle with defect carbon bridging structure (UP-RuNiSAs/C) via a unique unipolar pulse electrodeposition (UPED) strategy. As a result, the UP-RuNiSAs/C is found capable of running steadily for 100 h at 3 A cm-2, and shows a low overpotential of 9 mV at a current density of 10 mA cm-2 under alkaline conditions. Moreover, the UP-RuNiSAs/C allows an anion exchange membrane (AEM) electrolyzer to operate stably at 1.95 Vcell for 250 h at 1 A cm-2.

Assessing r2SCAN meta-GGA functional for structural parameters, cohesive energy, mechanical modulus, and thermophysical properties of 3d, 4d, and 5d transition metals

The recent development of accurate and efficient semilocal density functionals on the third rung of Jacob's ladder of density functional theory, such as the revised regularized strongly constrained and appropriately normed (r2SCAN) density functional, could enable rapid and highly reliable prediction of the elasticity and temperature dependence of thermophysical parameters of refractory elements and their intermetallic compounds using the quasi-harmonic approximation (QHA). Here, we present a comparative evaluation of equilibrium cell volumes, cohesive energy, mechanical moduli, and thermophysical properties (Debye temperature and thermal expansion coefficient) for 22 transition metals using semilocal density functionals, including the local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE) and PBEsol generalized gradient approximations (GGAs), and the r2SCAN meta-GGA. PBEsol and r2SCAN deliver the same level of accuracies for structural, mechanical, and thermophysical properties. PBE and r2SCAN perform better than LDA and PBEsol for calculating cohesive energies of transition metals. Among the tested density functionals, r2SCAN provides an overall well-balanced performance for reliably computing cell volumes, cohesive energies, mechanical properties, and thermophysical properties of various 3d, 4d, and 5d transition metals using QHA. Therefore, we recommend that r2SCAN could be employed as a workhorse method to evaluate thermophysical properties of transition metal compounds and alloys in high throughput workflows.

Theoretical Insight into the Special Synergy of Bimetallic Site in Co/MoC Catalyst to Promote N2 -to-NH3 Conversion

The catalytic mechanisms of nitrogen reduction reaction (NRR) on the pristine and Co/α-MoC(001) surfaces were explored by density functional theory calculations. The results show that the preferred pathway is that a direct N≡N cleavage occurs first, followed by continuous hydrogenations. The production of second NH3 molecule is identified as the rate-limiting step on both systems with kinetic barriers of 1.5 and 2.0 eV, respectively, indicating that N2 -to-NH3 transformation on bimetallic surface is more likely to occur. The two components of the bimetallic center play different roles during NRR process, in which Co atom does not directly participate in the binding of intermediates, but primarily serves as a reservoir of H atoms. This special synergy makes Co/α-MoC(001) have superior activity for ammonia synthesis. The introduction of Co not only facilitates N2 dissociation, but also accelerates the migration of H atom due to the antibonding characteristic of Co-H bond. This study offers a facile strategy for the rational design and development of efficient catalysts for ammonia synthesis and other reactions involving the hydrogenation processes.

Spin-Resolved Band Structure of Hoffman Clathrate [Fe(pz)2Pt(CN)4] as an Essential Tool to Predict Optical Spectra of Metal-Organic Frameworks

Paramount spin-crossover properties of the 3D-Hoffman metalorganic framework (MOF) [Fe(pz)₂Pt(CN)₄] are generally described on the basis of the ligand field theory, which provides adequate insight into theoretical and simulation analysis of spintronic complexes. However, the ligand field approximation does not take into account the 3D periodicity of the actual complex lattice and surface effects and therefore cannot predict a full-scale periodic structure without utilizing more advanced methods. Therefore, in this paper, the electronic properties of the exemplar MOF were analyzed from the band structure perspective in low-spin (LS) and high-spin (HS) states. The density-of-states spectra determined for both spin-up and spin-down electrons of Fe d⁶ orbitals indicate spin-orbital splitting and delocalization for HS due to spin polarization in the iron atom ligand field. Presence of the surface states in the real crystal causes a red shift of the metal-metal charge transfer (MMCT) and metal-ligand charge transfer (MLCT) peaks for both HS and LS states. The addition of residual water molecules and disorder among the pyrazine rings reveal additional influences on the positions of the pyrazine band and, therefore, on the absorption spectra of the crystal. The results show a magnification of the peak correlated with the MLCT in the HS state and a significant red shift of the LS characteristic absorption band. The presented approach involving band structure analysis delivers a more complete image of the electronic properties of the [Fe(pz)₂Pt(CN)₄] crystalline network and can be a landmark for insightful studies of other MOFs.

Effective Approaches for Designing Stable M-Nx /C Oxygen-Reduction Catalysts for Proton-Exchange-Membrane Fuel Cells

The large-scale commercialization of proton-exchange-membrane fuel cells (PEMFCs) is extremely limited by their costly platinum-group metals (PGMs) catalysts, which are used for catalyzing the sluggish oxygen reduction reaction (ORR) kinetics at the cathode. Among the reported PGM-free catalysts so far, metal-nitrogen-carbon (M-N_x /C) catalysts hold a great potential to replace PGMs catalysts for the ORR due to their excellent initial activity and low cost. However, despite tremendous progress in this field in the past decade, their further applications are restricted by fast degradation under practical conditions. Herein, the theoretical fundamentals of the stability of the M-N_x /C catalysts are first introduced in terms of thermodynamics and kinetics. The primary degradation mechanisms of M-N_x /C catalysts and the corresponding mitigating strategies are discussed in detail. Finally, the current challenges and the prospects for designing highly stable M-N_x /C catalysts are outlined.

Review on Magnetism in Catalysis: From Theory to PEMFC Applications of 3d Metal Pt-Based Alloys

The relationship between magnetism and catalysis has been an important topic since the mid-20th century. At present time, the scientific community is well aware that a full comprehension of this relationship is required to face modern challenges, such as the need for clean energy technology. The successful use of (para-)magnetic materials has already been corroborated in catalytic processes, such as hydrogenation, Fenton reaction and ammonia synthesis. These catalysts typically contain transition metals from the first to the third row and are affected by the presence of an external magnetic field. Nowadays, it appears that the most promising approach to reach the goal of a more sustainable future is via ferromagnetic conducting catalysts containing open-shell metals (i.e., Fe, Co and Ni) with extra stabilization coming from the presence of an external magnetic field. However, understanding how intrinsic and extrinsic magnetic features are related to catalysis is still a complex task, especially when catalytic performances are improved by these magnetic phenomena. In the present review, we introduce the relationship between magnetism and catalysis and outline its importance in the production of clean energy, by describing the representative case of 3d metal Pt-based alloys, which are extensively investigated and exploited in PEM fuel cells.

The evolution mechanism of an FeMo alloy catalyst for growth of single-walled carbon nanotubes

Adding small fractions of Mo to Fe nanoparticles (NPs) can reduce the melting point of FeMo NPs to lower than that of Fe NPs to prolong the lifetime of the alloy catalyst which in turn promotes the quality of catalytically synthesized single-walled carbon nanotubes (SWCNTs). In this study, we reveal the mechanism of the above-mentioned abnormal melting behavior by employing molecular dynamics simulations. Our results indicate that the bond length between the Fe atoms and the number of bonds between the Mo atoms play an important role in reducing the melting point of the FeMo NPs. This study provides useful insight into the evolution mechanism of the alloy catalyst for the growth of SWCNTs.

First-principle investigation of CO, CH4, and CO2 adsorption on Cr-doped graphene-like hexagonal borophene

It is important for life safety and scientific research to design new sensing materials for detecting CO, CH₄, and CO₂ from the environment. We theoretically designed a new Cr-doped graphene-like hexagonal borophene (CrB₆) as potential sensor material for these gases. Carrying out first-principle density-functional calculations, we calculated the adsorption energy, band structure, adsorption distance, charge transfer, charge density difference, density of states, and partial density of states of CO, CH₄, and CO₂ gas molecules absorbed on CrB₆ monolayer. The calculated results show that the adsorption behavior of CO is different from those of CH₄ and CO₂. CO adsorbed on CrB₆ monolayer prefers chemisorption with the adsorption energy of - 2.59 eV while CH₄ and CO₂ adsorbed on CrB₆ monolayer prefer physisorption with the adsorption energy of - 0.72 and - 0.69 eV. As a result, the different adsorption behaviors have significant influence on the band structures and density of states of CrB₆ monolayer. We hope that our results can help experimentalists synthesize better sensor materials based on hexagonal borophene.

Highly ordered nanoarrays catalysts embedded in carbon nanotubes as highly efficient and robust air electrode for flexible solid-state rechargeable zinc-air batteries

The development of multicomponent materials is the most efficient and successful way for creating advanced multifunctional catalysts. Herein, the bimetal FeCo nanoarrays enclosed N-CNTs have a high surface on carbon cloth support, which promotes efficient electron transport and prevents nanoparticle aggregation. Taking advantage of the high-level use of active material and fast charge transfer, the developed electrocatalyst exhibits excellent multifunctional electrocatalyst such as oxygen reduction reaction (ORR), oxygen evolution reaction (OER), and hydrogen evolution reaction (HER). The N-CNTs@MOF FeCo nanoarrays @CC exhibit higher activity than reference catalysts including MOF FeCo nanoarrays@CC, FeCo nanoarrays@CC, and CC. Interestingly, the synthesized multifunctional catalyst, which serves as the air electrode in zinc-air batteries with liquid electrolytes as well as solid-state gel electrolytes possesses outstanding charging-discharge performance and long service life. This study provides enormous potential for the real implementation of portable, even wearable, and efficient rechargeable batteries in the future.

An orbital principle to design P2-NaxMO2 cathode materials for sodium-ion batteries

Layered oxide materials are regarded to be the most promising high-performance cathode materials for sodium-ion batteries owing to their high working voltage and facile synthesis. Here, we study the influences of 3d transition metals on the cohesive energies, structural changes and operating voltages of P2-Na_xMO₂ during discharge based on first-principles calculations. Our results confirm that the performances of P2-Na_xMO₂ are associated with the chemical properties of the transition metals. In addition to this, we disclose that the involved orbitals of the 3d transition metal also greatly impact the electrochemical performance of the P2-Na_xMO₂ material during discharge according to the analysis of electronic structures. The jumps in the working voltage and volume during discharge are closely related to the occupation of the e_g and t_2g orbitals. Therefore, it is necessary to ensure that the discharge or charge process is carried out in one degenerate orbital to avoid jumps in the voltage and volume of the material. Our results could shed a light on the subsequent design of layered oxide cathodes with high cycle stability and a smooth voltage curve.

The influence of the oxygen vacancies on the Pt/TiO2 single-atom catalyst-a DFT study

The titanium dioxide (TiO₂) surface is suitable as a substrate for single-atom catalysts (SACs) for oxygen reduction reaction (ORR). As a common defect on TiO₂, oxygen vacancies may have a significant impact on the adsorption and activity of the adatoms. This work aims to investigate whether titanium dioxide containing surface oxygen vacancies is more suitable as a base material for SACs. This paper calculates the changes in the adsorption energy of the Pt atom and the energy of the d-band center on the perfect surface and the surface containing oxygen vacancies. Concerning the perfect surface, the surface containing oxygen vacancies fixes the Pt atom more firmly and increases the center energy of the d-band of Pt, thereby improving the performance of the Pt atom as SACs. Consequently, the (110) surface of rutile TiO₂ with oxygen vacancies may be the best substrate for SACs.

Ab initio DFT simulation of electronic and magnetic properties of Tin+1 and FeTin clusters

We report a computational investigation of the electronic and magnetic properties of neutral Ti_n+1 and FeTi_n (n = 1-10) clusters using ab initio calculations based on density functional theory (DFT) within the generalized gradient approximation (GGA). The best structures for Ti_n+1 and FeTi_n clusters are planar for size n < 5, while from n = 5, they showed a compact three-dimensional cage structure. For the best structures of the FeTi_n clusters, the Fe atoms favor the peripheral position with the highest coordination with the neighboring Ti atoms. The evolution as a function of the size of the average binding energies (Eb/atom) and HOMO-LUMO gaps of Ti_n+1 and FeTi_n (n = 1-10) clusters are studied. The stability results show that the Ti_n+1 clusters have relatively higher stability than the FeTi_n cluster with the same size. In addition, the vertical ionization potentials and electron affinities, chemical hardness, and atomic magnetic moment of Ti_n+1 and FeTi_n (n = 1-10) clusters are also investigated.

Transition-metal decorated graphdiyne monolayer as an efficient sensor toward phosphide (PH3) and arsine (AsH3)

Graphdiyne (GDY), a two-dimensional (2D) carbon, uniquely possesses mixed sp–sp2 hybridization, uniform nano-sized porous structure, semiconducting character, and excellent electrical conductivity.

Stability and Catalytic Performance of Single-Atom Supported on Ti2 CO2 for Low-Temperature CO Oxidation: A First-Principles Study

Based on first-principles calculations, the potential of Ti₂ CO₂ monolayer (MXene) as a single-atom catalyst (SAC) support for 3d transition metal (TM) atoms (Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn) is studied for CO oxidation. We first screen the support effect according to the stability of a single metal atom and find that Sc and Ti supported on Ti₂ CO₂ have stronger adsorption energy than the cohesive energy of their bulk counterparts and therefore, we selected Sc and Ti supported on Ti₂ CO₂ for further catalytic reactions. The stability and the potential catalytic reactivity are verified by electronic structure and charge transfer analysis. Both Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms are considered in this study, and lower energy barriers of 0.002 and 0.37 eV were found in the ER mechanism compared to the LH mechanism, which are 0.25 and 0.34 eV for Sc and Ti catalysts, respectively. Moreover, kinetic ER and LH mechanisms are favorable for both Sc- and Ti/Ti₂ CO₂ because of the comparable energy barrier to other metals and SAC supported on 2D materials. However, Ti/Ti₂ CO₂ catalyst is thermodynamically unfavorable. Based on these calculations, we propose that Sc supported on Ti₂ CO₂ is the best catalyst for CO-oxidation. The current study not only broadens the scope of the single-atom Sc catalyst but also extends the consideration of MXene support for catalyst optimization.

On the Origins of Some Spectroscopic Properties of "Purple Iron" (the Tetraoxoferrate(VI) Ion) and Its Pourbaix Safe-Space

In this work, the authors attempt to interpret the visible, infrared and Raman spectra of ferrate(VI) by means of theoretical physical-inorganic chemistry and historical highlights in this field of interest. In addition, the sacrificial decomposition of ferrate(VI) during water treatment will also be discussed together with a brief mention of how Rayleigh scattering caused by the decomposition of Fe^VIO₄²⁻ may render absorbance readings erroneous. This work is not a compendium of all the instrumental methods of analysis which have been deployed to identify ferrate(VI) or to study its plethora of reactions, but mention will be made of the relevant techniques (e.g., Mössbauer Spectroscopy amongst others) which support and advance this overall discourse at appropriate junctures, without undue elaboration on the foundational physics of these techniques.

Linking Bi-Metal Distribution Patterns in Porous Carbon Nitride Fullerene to Its Catalytic Activity toward Gas Adsorption

Immobilization of two single transition metal (TM) atoms on a substrate host opens numerous possibilities for catalyst design. If the substrate contains more than one vacancy site, the combination of TMs along with their distribution patterns becomes a design parameter potentially complementary to the substrate itself and the bi-metal composition. By means of DFT calculations, we modeled three dissimilar bi-metal atoms (Ti, Mn, and Cu) doped into the six porphyrin-like cavities of porous C₂₄N₂₄ fullerene, considering different bi-metal distribution patterns for each binary complex, viz. Ti_xCu_z@C₂₄N₂₄, Ti_xMn_y@C₂₄N₂₄, and Mn_yCu_z@C₂₄N₂₄ (with x, y, z = 0-6). We elucidate whether controlling the distribution of bi-metal atoms into the C₂₄N₂₄ cavities can alter their catalytic activity toward CO₂, NO₂, H₂, and N₂ gas capture. Interestingly, Ti₂Mn₄@C₂₄N₂₄ and Ti₂Cu₄@C₂₄N₂₄ complexes showed the highest activity and selectively toward gas capture. Our findings provide useful information for further design of novel few-atom carbon-nitride-based catalysts.

First-principles calculations of hybrid inorganic-organic interfaces: from state-of-the-art to best practice

The computational characterization of inorganic-organic hybrid interfaces is arguably one of the technically most challenging applications of density functional theory. Due to the fundamentally different electronic properties of the inorganic and the organic components of a hybrid interface, the proper choice of the electronic structure method, of the algorithms to solve these methods, and of the parameters that enter these algorithms is highly non-trivial. In fact, computational choices that work well for one of the components often perform poorly for the other. As a consequence, default settings for one materials class are typically inadequate for the hybrid system, which makes calculations employing such settings inefficient and sometimes even prone to erroneous results. To address this issue, we discuss how to choose appropriate atomistic representations for the system under investigation, we highlight the role of the exchange-correlation functional and the van der Waals correction employed in the calculation and we provide tips and tricks how to efficiently converge the self-consistent field cycle and to obtain accurate geometries. We particularly focus on potentially unexpected pitfalls and the errors they incur. As a summary, we provide a list of best practice rules for interface simulations that should especially serve as a useful starting point for less experienced users and newcomers to the field.

Role of initial stage nitridation on the mechanical properties of an α-Fe(100) nanofilm in NH3

Nitrogen is one of the most significant non-native interstitial elements that is present in the structure of Fe. Initial stage nitridation dramatically influences the mechanical properties of steel, especially for micro to nanoscale applications, but is not yet fully understood. By means of reactive force field molecular dynamics (ReaxFF MD) simulations, the initial stage of the nitridation process of nanofilm Fe, as well as its role on the mechanical properties of the material, were investigated. To clarify the temperature effect, nitridation was simulated in the range of 500-900 K, demonstrating that the adsorption of both N and H atoms into Fe was enhanced by thermal actuation. Corresponding tension test simulations were performed, manifesting that the Fe nanofilm nitrided at 600 K presents the highest yield stress. Further analysis shows that there is a competitive mechanism between the inward diffusion of N atoms that enhances the strength and simultaneous adsorption of H atoms, which leads to brittleness of the material as the temperature increases. Hence, an intermediate temperature could lead to optimal mechanical properties due to the balance of improving the strength while controlling the brittleness of the material. To probe the deformation mechanism, evolutions of partial dislocation and twin boundary at plasticity beginning for pure Fe and the nitrided Fe nanofilm are discussed. The present results show the nitridation strengthening technology of Fe in NH3 and its related microscale mechanism, which may theoretically support the technical design and improvement in the properties of steel.

London Dispersion Corrections to Density Functional Theory for Transition Metals Based on Fitting to Experimental Temperature-Programmed Desorption of Benzene Monolayers

Standard implementations of generalized gradient approximation (GGA)-based density functional theory (DFT) describe well strongly bound molecules and solids but fail to describe long-range London dispersion or van der Waals (vdW) attraction interactions that are important in molecular crystals and two-dimensional solids. To provide accurate values for the vdW distance and energies for the metals Cu, Ag, Au, Ni, Pd, and Pt, we determined empirical vdW corrections to Perdew, Burke, and Ernzerhof (PBE) DFT by fitting the experimental adsorption enthalpies measured by temperature-programmed desorption (TPD) from benzene monolayers by Campbell and co-workers ( J. Phys. Chem. C 2016, 120, 25161-25172). Benzene physisorbed to these metals without chemical reaction; therefore, we consider the bonding to be vdW. We use the low gradient form for the vdW corrections, E_vdW-LG = -C_6LG/[R⁶ + R_vdwLG⁶] with just two parameters per atom (C_6LG and R_vdwLG). This LG form leads to negligible changes in bond distances and angles, so adjusting the parameters should not sacrifice accuracy for the bonding interactions. We demonstrate that the parameters fitted to benzene also describe well the physisorption enthalpies for other hydrocarbons (naphthalene, cyclohexane, methane, ethane, and propane) on Pt. We also report low gradient vdW correction parameters for the noble gases that fit the equilibrium lattice parameter and heat of vaporization of the crystals.

Theoretical understanding of the properties of stepped iron surfaces with van der Waals interaction corrections

The stepped surfaces in nanoscale zero-valent iron (nZVI) play an essential role for environmental application. However, there is still currently a deficiency in the atomic understanding of stepped surface properties due to the limitation of the computational methodology. In this study, stepped Fe(210) and (211) surfaces were theoretically investigated using density functional theory (DFT) computations in terms of the flat Fe(110) surface. Our results suggest that the consideration of van der Waals (vdW) interaction correction is beneficial for the DFT study on Fe-based systems. The DF-cx method is found to be the most promising vdW correction method. The DF-cx results reveal that the stepped Fe(210) and Fe(211) surfaces experience significant surface relaxation and abnormal trends in their work function. Their electronic properties and reactivities of the surface atoms are strongly affected by the Fe coordination numbers and the strong adsorption strengths of oxygen on the surfaces are dependent on both the coordination number of the adsorbed atoms and the geometry of the adsorption sites.

The role of an unintentional carbon dopant in resolving the controversial conductivity aspects in BiFeO3

Light elements like carbon may enter unintentionally into a material during material processing owing to their ubiquitous nature, and may significantly influence its observed electronic and magnetic properties. In the present work, the energetics and kinetics of carbon impurity related defects in BiFeO₃ (BFO) are studied using first principles calculations in order to gain insight into the ongoing controversial aspects of conductivity of BFO. The results suggest that oxygen deficient conditions provide a favorable chemical environment for incorporation of carbon in BFO. Calculations based on the formation energy predict that carbon can spontaneously occupy interstitials, O, and Fe sites in BFO (where it is found to introduce impurity induced shallow acceptor type states at an energy of 0.05 eV above the valence band maximum). Carbon occupying cationic sites (C_Bi and C_Fe) tends to ionize their vacancies (V_Bi and V_Fe), resulting in the formation of a CO₃ cluster, whereas it induces localized electron traps with energy levels composed of impurity states near the center of the band gap (0.9 eV above the valence band maximum) when occupying interstitial sites in BFO. An understanding of the migration of C impurity in BFO is developed, which suggests the favorable incorporation of carbon impurity via a vacancy mechanism. In order to confirm the theoretical results, experimental studies are carried out where BFO and carbon doped BFO (BCFO) thin films are grown by the pulsed laser deposition technique. Polycrystalline pure phase (R3c) thin films of BFO and BCFO are obtained. The presence of defect states in the deposited thin films is optically analyzed by the photoluminescence (PL) technique. In order to highlight the critical role of carbon in modifying the electrical conductivity of BFO, a BCFO/BFO/ITO based p-i-n heterojunction is prepared. The electrical characteristics depict remarkable rectifying characteristics, thus suggesting the p-type nature of carbon dopant in otherwise intrinsic BFO.

2D PC3 as a promising thermoelectric material

In the present article, we report the thermoelectric properties of monolayer PC3 for the first time. The structural, vibrational, electronic and thermoelectric properties of PC3 are investigated in detail using a combination of density functional and Boltzman transport theory, and are compared to the carbon (graphene) and phosphorous (phosphorene) analogues. The results show that the PC3 monolayer is dynamically stable and robust upon oxygen contact as well. Also, PC3 is found to be an indirect band gap semiconductor in comparison to the zero gap carbon (graphene) and direct gap phosphorous (phosphorene) analogues. The effect of axial strains is also investigated on the electronic and thermoelectric properties of PC3. The present work reveals monolayer PC3 to be an excellent thermoelectric material with significant thermoelectric performance (ZT ∼ 1) for a large scale operating temperature range of 200-1200 K.

First principles analysis of surface dependent segregation in bimetallic alloys

Stability is an important aspect of alloys, and proposed alloys may be unstable due to unfavorable atomic interactions. Segregation of an alloy may occur preferentially at specific exposed surfaces, which could affect the alloy's structure since certain surfaces may become enriched in certain elements. Using density functional theory (DFT), we modeled surface segregation in bimetallic alloys involving all transition metals doped in Pt, Pd, Ir, and Rh. We not only modeled common (111) surfaces of such alloys, but we also modeled (100), (110), and (210) facets of such alloys. Segregation is more preferred for early and late transition metals, with middle transition metals being most stable within the parent metal. We find these general trends in segregation energies for the parent metals: Pt > Rh > Pd > Ir. A comparison of different surfaces suggests no consistent trends across the different parent hosts, but segregation energies can vary up to 2 eV depending on the exposed surface. We also developed a statistical model to predict surface-dependent segregation energies. Our model is able to distinguish segregation at different surfaces (as opposed to generic segregation common in previous models), and agrees well with the DFT data. The present study provides valuable information about surface-dependent segregation and helps explain why certain alloy structures occur (e.g. core-shell).

Tailoring the capability of carbon nitride (C3N) nanosheets toward hydrogen storage upon light transition metal decoration

To nurture the full potential of hydrogen (H₂) as a clean energy carrier, its efficient storage under ambient conditions is of great importance. Owing to the potential of material-based H₂ storage as a promising option, we have employed here first principles density functional theory calculations to study the H₂ storage properties of recently synthesized C₃N monolayers. Despite possessing fascinating structural and mechanical properties C₃N monolayers weakly bind H₂ molecules. However, our van der Waals corrected simulations revealed that the binding properties of H₂ on C₃N could be enhanced considerably by suitable Sc and Ti doping. The stabilities of Sc and Ti dopants on a C₃N surface has been verified by means of reaction barrier calculations and ab initio molecular dynamics simulations. Upon doping with C₃N, the existence of partial positive charges on both Sc and Ti causes multiple H₂ molecules to bind to the dopants through electrostatic interactions with adsorption energies that are within an ideal range. A drastically high H₂ storage capacity of 9.0 wt% could be achieved with two-sided Sc/Ti doping that ensures the promise of C₃N as a high-capacity H₂ storage material.

Perspective: Chemical Information Encoded in Electron Density

In this perspective, we review the chemical information encoded in electron density and other ingredients used in semilocal functionals. This information is usually looked at from the functional point of view: the exchange density or the enhancement factor are discussed in terms of the reduced density gradient. However, what parts of a molecule do these 3D functions represent? We look at these quantities in real space, aiming to understand the electronic structure information they encode and provide an insight from the quantum chemical topology (QCT). Generalized gradient approximations (GGAs) provide information about the presence of chemical interactions, whereas meta-GGAs can differentiate between the different bonding types. By merging these two techniques, we show new insight into the failures of semilocal functionals owing to three main errors: fractional charges, fractional spins, and non-covalent interactions. We build on simple models. We also analyze the delocalization error in hydrogen chains, showing the ability of QCT to reveal the delocalization error introduced by semilocal functionals. Then, we show how the analysis of localization can help understand the fractional spin error in alkali atoms, and how it can be used to correct it. Finally, we show that the poor description of GGAs of isodesmic reactions in alkanes is due to 1,3-interactions.

Evaluation and comparison of classical interatomic potentials through a user-friendly interactive web-interface

Classical empirical potentials/force-fields (FF) provide atomistic insights into material phenomena through molecular dynamics and Monte Carlo simulations. Despite their wide applicability, a systematic evaluation of materials properties using such potentials and, especially, an easy-to-use user-interface for their comparison is still lacking. To address this deficiency, we computed energetics and elastic properties of variety of materials such as metals and ceramics using a wide range of empirical potentials and compared them to density functional theory (DFT) as well as to experimental data, where available. The database currently consists of 3248 entries including energetics and elastic property calculations, and it is still increasing. We also include computational tools for convex-hull plots for DFT and FF calculations. The data covers 1471 materials and 116 force-fields. In addition, both the complete database and the software coding used in the process have been released for public use online (presently at http://www.ctcms.nist.gov/∼knc6/periodic.html) in a user-friendly way designed to enable further material design and discovery.

Chemical states of 3d transition metal impurities in a liquid lead–bismuth eutectic analyzed using first principles calculations

Interactions of 4s and 3d orbitals with the 6p band largely determine the chemical states of 3d transition-metal impurities in a liquid LBE.

Fe– and Co–P4-embedded graphenes as electrocatalysts for the oxygen reduction reaction: theoretical insights

Fe–P₄-embedded graphene exhibits high catalytic activity for the ORR in alkaline media.

The design of d-character Dirac cones based on graphene

We introduce a new framework for designing a transition metal (TM) d-electrons dominant Dirac cone spectrum based on the hybridization between graphene and a modulated TM d impurity band. The obtained Dirac cone behaves like a 'copy' from graphene, insensitive to the TM coverage and order. First-principles calculations reveal such a system of Mn intercalated epitaxial graphene on SiC(0 0 0 1), dubbed manganosine. The robustness of the Dirac cone is discussed in terms of the possible imperfection of Mn atoms. The mechanism at work is expected to be rather general and may open the door to designing new d- or f-character Dirac systems.

First principles computational study on the electrochemical stability of Pt-Co nanocatalysts

Using density functional theory (DFT) calculations, we identify the thermodynamically stable configurations of Pt-Co alloy nanoparticles of varying Co compositions and particle sizes. Our results indicate that the most thermodynamically stable structure is a shell-by-shell configuration where the Pt atom only shell and the Co only shell alternately stack and the outermost shell consists of a Pt skin layer. DFT calculations show that the structure has substantially higher dissolution potential of the outermost Pt shell compared with pure Pt nanoparticles of approximately the same size. Furthermore, our DFT calculations also propose that the shell-by-shell structure shows much better oxygen reduction reaction (ORR) activity than conventional bulk or nanoparticles of pure Pt. These novel catalyst properties can be changed when the surfaces are adsorbed with oxygen atoms via selective segregation followed by the electrochemical dissolution of the alloyed Co atoms. However, these phenomena are thermodynamically not plausible if the chemical potentials of oxygen are controlled below a certain level. Therefore, we propose that the shell-by-shell structures are promising candidates for highly functional catalysts in fuel cell applications.