Novel method of self‐interaction corrections in density functional calculations

It is demonstrated that the commonly applied self‐interaction correction (SIC) used in density functional theory does not remove all self‐interaction. We present as an alternative a novel method that, by construction, is totally free from self‐interaction. The method has the correct asymptotic 1/r dependence. We apply the new theory to localized f electrons in praseodymium and compare with the old version of SIC, the local density approximation (LDA) and with an atomic Hartree–Fock calculation. The results show a lowering of the f level, a contraction of the f electron cloud and a lowering of the total energy by 13 eV per 4 f electron compared to LDA. The equilibrium volume of the new SIC method is close to the ones given by LDA and the older SIC method and is in good agreement with experiment. The experimental cohesive energy is in better agreement using the new SIC method, both compared to LDA and another SIC method. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 247–252, 2001     

Crystal structures and thermodynamic investigations of NaSc(BH4)4 from first‐principles calculations

The double cation borohydride NaSc(BH₄)₄ has a total H₂ content of 12.67 wt.% and has been suggested as a potential candidate for hydrogen storage applications. This study reports first‐principles calculations of the structure and reaction thermodynamics of NaSc(BH₄)₄. The calculations indicate that NaSc(BH₄)₄ is decomposed into a mixture of ScB₂, NaBH₄, and Na₂(B₁₀H₁₀) with H₂ release of 9.3 wt.% at 118 K at a partial pressure of H₂ of 100 bar. Reactant compositions that can destabilize NaSc(BH₄)₄ were evaluated. This effort identified four destabilization reactions that are predicted to have reaction thermodynamics for hydrogen release within the temperature range of 78–109 K. Even though the reactions conclusively produce undesired compounds, such as refractory materials or kinetically stable B₁₂H₁₂‐containing species, the thermodynamic study suggests a direction for improving the thermodynamics of double cation borohydride‐based systems being actively considered for hydrogen storage applications. © 2012 Wiley Periodicals, Inc.     

First‐principles calculation of electronic structure of V‐doped anatase TiO2

The energetic and electronic structures of V‐doped anatase TiO₂ have been investigated systematically by the GGA+U approach, including replacement of Ti by V in the absence and presence of oxygen vacancies and the presence of an interstitial site. It was found that V should exist as a V⁴⁺ ion in the replacement of Ti in the anatase lattice, the electron transitions of which to the conduction band from V 3d states are responsible for the experimentally observed visible light absorption. The influence of V dopant concentration on the electronic and magnetic properties is also discussed, such as the influence of the U value in systems containing oxygen vacancies and spin flip phenomena for interstitial V‐doping.     

A first‐principles study on the existence and structures of the lighter alkaline‐earth pernitrides

The lighter alkaline‐earth pernitrides BeN₂, MgN and CaN₂ have been structurally predicted by a series of density‐functional (GGA/PBE/PAW) electronic‐structure calculations. Despite their crystal chemistry clearly pointing towards the formulation M²⁺N₂^2? with an N? N distance of 1.26 Å, all phases turn out as metallic compounds which are exothermic with respect to the elements. The M²⁺ coordination numbers are a simple function of the cationic radius. The bulk moduli are about three times smaller than those of the noble‐metal pernitrides, a consequence of the smaller anionic charge in the former phases. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

First‐principles calculations on the four phases of BaTiO3

The calculations based on linear combination of atomic orbitals basis functions as implemented in CRYSTAL09 computer code have been performed for cubic, tetragonal, orthorhombic, and rhombohedral modifications of BaTiO₃ crystal. Structural and electronic properties as well as phonon frequencies were obtained using local density approximation, generalized gradient approximation, and hybrid exchange‐correlation density functional theory (DFT) functionals for four stable phases of BaTiO₃. A comparison was made between the results of different DFT techniques. It is concluded that the hybrid PBE0 [J. P. Perdew, K. Burke, M. Ernzerhof, J. Chem. Phys. 1996, 105, 9982.] functional is able to predict correctly the structural stability and phonon properties both for cubic and ferroelectric phases of BaTiO₃. The comparative phonon symmetry analysis in BaTiO₃ four phases has been made basing on the site symmetry and irreducible representation indexes for the first time. © 2012 Wiley Periodicals, Inc.     

First‐Principles Calculations on Electronic Structures of N/V‐Doped and N‐V‐Dodoped Anatase TiO2 (101) Surfaces

The energetic and electronic properties of N/V‐doped and N‐V‐codoped anatase TiO₂ (101) surfaces are investigated by first‐principles calculations, with the aim to elucidate the relationship between the electronic structure and the photocatalytic performance of N‐V‐codoped TiO₂. Several substitutional and interstitial configurations for the N and/or V impurities in the bulk phase and on the surface are studied, and the relative stability of different doping configurations is compared by the impurity formation energy. Systematic calculations reveal that N and V impurities can be encapsulated by TiO₂ to form stable structures as a result of strong N‐V interactions both in the bulk and the surface model. Through analyzing and comparing the electronic structures of different doping systems, the synergistic doping effects are discussed in detail. Based on these discussions, we suggest that N_OV_Ti codoping cannot only narrow the band gap of anatase TiO₂, but also forms impurity states, which are propitious for the separation of photoexcited electron–hole pairs. In the case of N_OV_Ti‐codoped TiO₂ (101) surfaces, this phenomenon is especially prominent. Finally, a feasible synthesis route for N_OV_Ti codoping into anatase TiO₂ is proposed.     

A combinatorial study of inverse Heusler alloys by first‐principles computational methods

In continuation of our recent combinatorial work on 810 X₂YZ full Heusler alloys, a computational study of the same class of materials but with the inverse (XY)XZ crystal structure has been performed on the basis of first‐principles (GGA) total‐energy calculations using pseudopotentials and plane waves. The predicted enthalpies of formation evidence 27 phases to be thermochemically stable against the elements and the regular X₂YZ type. A chemical‐bonding study yields an inherent tendency for structural distortion in a majority of these alloys, and we predict the existence of the new tetragonal phase Fe₂CuGa (P4₂/ncm; a = 5.072 Å, c = 7.634 Å; c/a ≈ 1.51) with a saturation moment of μ = 4.69 μ_B per formula unit. Thirteen more likewise new, isotypical phases are predicted to show essentially the same behavior. Six phases turn out to be the most stable in the inverse tetragonal arrangement. The course of the magnetic properties as a function of the valence‐electron concentration is analyzed using a Slater‐Pauling approach. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Interaction potential models for bulk Zns,Zns nanoparticle,and Zns nanoparticle‐PMMA from first‐principles

An ab initio derived transferable polarizable force‐field has been developed for Zinc sulphide (ZnS) nanoparticle (NP) and ZnS NP‐PMMA nanocomposite. The structure and elastic constants of bulk ZnS using the new force‐field are within a few percent of experimental observables. The new force‐field show remarkable ability to reproduce structures and nucleation energies of nanoclusters (Zn₁S₁‐Zn₁₂S₁₂) as validated with that of the density functional theory calculations. A qualitative agreement of the radial distribution functions of Zn? O, in a ZnS nanocluster‐PMMA system, obtained using molecular mechanics molecular dynamics (MD) and ab initio MD (AIMD) simulations indicates that the ZnS–PMMA interaction through Zn? O bonding is explained satisfactorily by our force‐field. © 2015 Wiley Periodicals, Inc.     

Design of low band gap polymers employing density functional theory—hybrid functionals ameliorate band gap problem

Band gaps in solids and excitation energies in finite systems are underestimated significantly if estimated from differences between eigenvalues obtained within the local spin density approximation (LSDA). In this article we present results on 20 small- and medium-sized π-systems which show that HOMO–LUMO energy differences obtained with the B3LYP, B3P86, and B3PW91 functionals are in good agreement with vertical excitation energies from UV-absorption spectra. The improvement is a result of the use of the exact Hartree–Fock exchange with hybrid methods. Negative HOMO energies and negative LUMO energies do not provide good estimates for IPs and EAs. In contrast to Hartree–Fock theory, where IPs are approximated well and EAs are given poorly, DFT hybrid methods underestimate IPs and EAs by about the same amount. LSDA yields reasonable EAs but poor IPs. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1943–1953, 1997     

A DFT study on 1,4‐dihydro‐1,4‐azaborinine annulated linear polyacenes: Absorption spectra,singlet‐triplet energy gap,aromaticity, and HOMO–LUMO energy modulation

Linear polyacene (LPA) mimics containing multiple heterocycles have been computationally designed by annulating 1,4‐dihydro‐1,4‐azaborinine moieties to benzene (aB₁–aB₅), naphthalene (aN₁–aN₅), anthracene (aA₁–aA₅), and tetracene (aT₁–aT₅) cores. DFT studies conducted on them using M06L/6‐311++G(d,p) method reveal a perfect planar structure for all and suggest the utilization of nitrogen lone pairs for aromatic π‐electron delocalization. The computed values of aromaticity indices such as HOMA, NICS, and dehydrogenation energy (E _dh) of heterocycles support strong aromatic character for each six‐membered ring in the LPA mimics. On the basis of the minimum value of the molecular electrostatic potential (V _min) observed on each LPA unit in the LPA mimics, the extended delocalization of π‐electrons is verified. The energetic parameter E _dh showed strong linear correlation with HOMA, NICS and V _min parameters, which strongly supports the multidimensional character of aromaticity in LPA mimics. The electronic property modification is shown by the theoretical absorption spectra data and singlet‐triplet energy gap (ΔE _ST). The bandgap and ΔE _ST tunings are achieved for LPA mimics by selecting appropriate number of azaborinine type units and the size of LPA core used for annulation. © 2017 Wiley Periodicals, Inc.     

First‐principles calculations on the energetics,electronic structures and magnetism of SrFeO2

The electronic and magnetic properties of SrFeO₂ with different magnetic configurations have been calculated via the plane‐wave pseudopotential density functional theory method, using the experimental lattice parameters. The results give an antiferromagnetic ground state for SrFeO₂ with an absolute magnetic moment agreeing very well with the experimental report. In comparison with the counterparts whose magnetic moments are parallel to the c axis, the structures with spin moments parallel to the a (or b) axis exhibit no observable preference in total energy, but show different density distributions of the Fe 3d and Fe 3d states. The square‐planar crystal field splits the Fe 3d orbitals into a high‐level d, a low d, and intermediate d_xy and d_xz or d_yz components. The exchange splitting is larger than the crystal‐field splitting, resulting in the high‐spin Fe 3d states. Referred to the triplet O₂, the O‐vacancy formation energy from SrFeO₃ to SrFeO₂ has been deduced as well, along with its dependence on the temperature and O₂ partial pressure. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

Transition metal atoms on oxide supports density functional calculations

The adsorption properties of Cu, Ag, Ni, and Pd atoms on O^2?, F, and F⁺ sites of MgO, CaO, SrO, and BaO (001) surfaces have been studied by means of density functional calculations. The examined clusters were embedded in the simulated Coulomb fields that closely approximate the Madelung fields of the host surfaces. The adsorption properties have been analyzed with reference to the basicity and energy gap of the oxide support in addition to orbital interactions. While the free Ni d⁹s¹ triplet ground state is preserved on adsorption on the O^2? sites of MgO, CaO, and SrO surfaces, it is no longer preserved on the O^2? site of BaO. For all adsorbates considered, adsorption is found to be stronger on F⁺ sites compared with regular O^2? sites. While on the O^2? site, Pd and Ni form the most stable complexes, on the F site, Pd and Cu form the most stable counterparts. On the F⁺ site, the single valence electron of Cu and Ag atoms couples with the unpaired electron of the vacancy forming a covalent bond. As a result, the adsorption energies of these atoms on the F⁺ site are stronger than those on the F and O^2? sites. The adsorption properties correlate linearly with the basicity and energy gap of the oxide support in addition to orbital interactions. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Electronic structure,magnetic properties,and mixed valence character of Ce2Ni3Si5 from first principles calculations

Cerium intermetallic compounds exhibit anomalous physical properties such as heavy fermion and Kondo behaviors. Here, an ab initio study of the electronic structure, magnetic properties, and mixed valence character of Ce₂Ni₃Si₅ using density functional theory (DFT) is presented. Two theoretical methods, including pure Perdew–Burke–Ernzerhof (PBE) and PBE + U , are used. In this study, Ce³⁺ and Ce⁴⁺ are considered as two different constituents in the unit cell. The formation energy calculations on the DFT level propose that Ce is in a stable mixed valence of 3.379 at 0 K. The calculated electronic structure shows that Ce₂Ni₃Si₅ is a metallic compound with a contribution at the Fermi level from Ce 4f and Ni 3d states. With the inclusion of the effective Hubbard parameter (U _eff), the five valence electrons of 5 Ce³⁺ ions are distributed only on Ce³⁺ 4f orbitals. Therefore, the occupied Ce³⁺ 4f band is located in the valence band (VB) while Ce⁴⁺ 4f orbitals are empty and Located at the Fermi level. The calculated magnetic moment in Ce₂Ni₃Si₅ is only due to cerium (Ce³⁺) in good agreement with the experimental results. The U _eff value of 5.4 eV provides a reasonable magnetic moment of 0.981 for the unpaired electron per Ce³⁺ ion. These results may serve as a guide for studying present mixed valence cerium‐based compounds. © 2017 Wiley Periodicals, Inc.     

Temperature effects,energy shifts,and band entropies of Si semiconductors

Understanding the electronic properties of silicon semiconductors is important for the preparation of high-performance semiconductor materials. We calculated the band entropies, electronic structures, and bonding properties of a silicon semiconductor using density functional theory and the binding-energy and bond-charge model. The relationship between Si energy and temperature was studied using the tight binding (TB) approximation and bond-order-length-strength (BOLS) theory (BOLS-TB), with the Si (111) surface as an example. The specific binding energies and bonding properties of Si atoms in different surface atomic layers are discussed by analyzing the X-ray photoelectron spectra of the Si (111) surface at 953 and 1493 K. This study improves our understanding of how surface properties reflect local bonding states and deepens our understanding of how atomic-relaxation-derived Hamiltonian perturbations and temperature influence the binding energy of the surface region. It also contributes to the development of Si-based semiconductor materials by providing new ideas and methods.     

A first‐principles study of diatomic NiAl: Ground state,structure, and spectroscopic constants

A computational study of diatomic NiAl is reported. Molecular properties evaluated include the equilibrium bond length (r_e), equilibrium stretching frequency (ω_e), doublet‐quartet energy splitting, and nickel‐aluminum bond strength. Several interesting conclusions have resulted from this research. First, convergence in calculated properties is smoother with recently reported correlation consistent basis sets than earlier basis sets for Ni and Al. Second, with the exception of bond strength, basis set limit properties extrapolated using correlation basis sets are in agreement with reported data. Third, this research suggests that caution may be needed with regard to the use of DFT for developing interatomic potentials for larger scale simulations. For example, B97‐1 showed better agreement with reported r_e for ²NiAl than B3LYP. However, the situation was reversed for the calculation of ω_e. With respect to bond strength, the situation is unclear due to the scatter among experiment and calculations. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

The Influence of Defects on Mo‐Doped TiO2 by First‐Principles Studies

TiO₂ doped with transition metals shows improved photocatalytic efficiency. Herein the electronic and optical properties of Mo‐doped TiO₂ with defects are investigated by DFT calculations. For both rutile and anatase phases of TiO₂, the bandgap decreases continuously with increasing Mo doping level. The 4d electrons of Mo introduce localized states into the forbidden band of TiO₂, and this shifts the absorption edge into the visible‐light region and enhances the photocatalytic activity. Since defects are universally distributed in TiO₂ or doped TiO₂, the effect of oxygen deficiency due to oxygen vacancies or interstitial Mo atoms is systemically studied. Oxygen vacancies associated with the Mo dopant atoms or interstitial Mo will reduce the spin polarization and magnetic moment of Mo‐doped TiO₂. Moreover, oxygen deficiency has a negative impact on the improved photocatalytic activity of Mo‐doped TiO₂. The current results indicate that substitutional Mo, interstitial Mo, and oxygen vacancy have different impacts on the electronic/optical properties of TiO₂ and are suited to different applications.     

Kinetic energy density for orbital‐free density functional calculations by axiomatic approach

An axiomatic approach is herein used to determine the physically acceptable forms for general D‐dimensional kinetic energy density functionals (KEDF). The resulted expansion captures most of the known forms of one‐point KEDFs. By statistically training the KEDF forms on a model problem of noninteracting kinetic energy in 1D (six terms only), the mean relative accuracy for 1000 randomly generated potentials is found to be better than the standard KEDF by several orders of magnitudes. The accuracy improves with the number of occupied states and was found to be better than for a system with four occupied states. Furthermore, we show that free fitting of the coefficients associated with known KEDFs approaches the exactly analytic values. The presented approach can open a new route to search for physically acceptable kinetic energy density functionals and provide an essential step toward more accurate large‐scale orbital free density functional theory calculations.     

Sterically Less‐Hindered Half‐Titanocene(IV) Phenoxides: Ancillary‐Ligand Effect on Mono‐, Bis‐, and Tris(2‐Alkyl‐/arylphenoxy) Titanium(IV) Chloride Complexes

A series of mono‐, bis‐, and tris(phenoxy)–titanium(IV) chlorides of the type [Cp*Ti(2‐R? PhO)_nCl_3?n] (n=1–3; Cp*=pentamethylcyclopentadienyl) was prepared, in which R=Me, iPr, tBu, and Ph. The formation of each mono‐, bis‐, and tris(2‐alkyl‐/arylphenoxy) series was authenticated by structural studies on representative examples of the phenyl series including [Cp*Ti(2‐Ph? PhO)Cl₂] ( 1 PhCl2 ), [Cp*Ti(2‐Ph? PhO)₂Cl] ( 2 PhCl ), and [Cp*Ti(2‐Ph? PhO)₃] ( 3 Ph ). The metal‐coordination geometry of each compound is best described as pseudotetrahedral with the Cp* ring and the 2‐Ph? PhO and chloride ligands occupying three leg positions in a piano‐stool geometry. The mean Ti? O distances, observed with an increasing number of 2‐Ph? PhO groups, are 1.784(3), 1.802(4), and 1.799(3) Å for 1 PhCl2 , 2 PhCl , and 3 Ph , respectively. All four alkyl/aryl series with Me, iPr, tBu, and Ph substituents were tested for ethylene homopolymerization after activation with Ph₃C⁺[B(C₆F₅)₄]^? and modified methyaluminoxane (7% aluminum in isopar E; mMAO‐7) at 140 °C. The phenyl series showed much higher catalytic activity, which ranged from 43.2 and 65.4 kg (mmol of Ti?h)^?1, than the Me, iPr, and tBu series (19.2 and 36.6 kg (mmol of Ti?h)^?1). Among the phenyl series, the bis(phenoxide) complex of 2 PhCl showed the highest activity of 65.4 kg (mmol of Ti?h)^?1. Therefore, the catalyst precursors of the phenyl series were examined by treating them with a variety of alkylating reagents, such as trimethylaluminum (TMA), triisobutylaluminum (TIBA), and methylaluminoxane (MAO). In all cases, 2 PhCl produced the most catalytically active alkylated species, [Cp*Ti(2‐Ph? PhO)MeCl]. This enhancement was further supported by DFT calculations based on the simplified model with TMA.     

Density Functional Investigations of Electronic Structure and Dehydrogenation Reactions of Al‐ and Si‐Substituted Magnesium Hydride

The effect on the hydrogen storage attributes of magnesium hydride (MgH₂) of the substitution of Mg by varying fractions of Al and Si is investigated by an ab initio plane‐wave pseuodopotential method based on density functional theory. Three supercells, namely, 2×2×2, 3×1×1 and 5×1×1 are used for generating configurations with varying amounts (fractions x=0.0625, 0.1, and 0.167) of impurities. The analyses of band structure and density of states (DOS) show that, when a Mg atom is replaced by Al, the band gap vanishes as the extra electron occupies the conduction band minimum. In the case of Si‐substitution, additional states are generated within the band gap of pure MgH₂—significantly reducing the gap in the process. The reduced band gaps cause the Mg? H bond to become more susceptible to dissociation. For all the fractions, the calculated reaction energies for the stepwise removal of H₂ molecules from Al‐ and Si‐substituted MgH₂ are much lower than for H₂ removal from pure MgH₂. The reduced stability is also reflected in the comparatively smaller heats of formation (ΔH_f) of the substituted MgH₂ systems. Si causes greater destabilization of MgH₂ than Al for each x. For fractions x=0.167 of Al, x=0.1, 0.167 of Si (FCC) and x=0.0625, 0.1 of Si (diamond), ΔH_f is much less than that of MgH₂ substituted by a fraction x=0.2 of Ti (Y. Song, Z. X. Guo, R. Yang, Mat. Sc. & Eng. A 2004 , 365, 73). Hence, we suggest the use of Al or Si instead of Ti as an agent for decreasing the dehydrogenation reaction and energy, consequently, the dehydrogenation temperature of MgH₂, thereby improving its potential as a hydrogen storage material.     

Density Functional Characterization of the Electronic Structure and Visible‐Light Absorption of Cr‐Doped Anatase TiO2

To evaluate the electronic and optical properties of Cr‐doped anatase TiO₂, three possible Cr‐doped TiO₂ models, including Cr at a Ti site (model I), Cr at a Ti site with an oxygen vacancy compensation (model II), and an interstitial Cr site (model III), are studied by means of first principles density functional theory calculations. In model I, the splitting behavior of the Cr 3d states and the insulating properties are successfully depicted by the GGA+U method, from which it is proposed that Cr at a Ti site should exist as Cr⁴⁺ instead of the generally believed Cr³⁺. As a result, the electron transitions between these impurity states, the conduction band (CB), and the valence band (VB), as well as the d–d transitions between occupied and unoccupied Cr 3d states, provide a reasonable explanation for the experimentally observed major and minor absorption bands. In models II and III, the impurity states and associated optical transition processes—as well as the corresponding electron configurations—are examined.