Hybrid Density Functional Theory Methods Research Articles

Predicting decomposition temperatures and carbonization feasibility of heterocycles with their decomposition mechanism

Pyrolysis of heterocyclic compounds such as Meldrum’s acid derivatives was studied with two key mechanisms analysing their energy barrier diagrams by the hybrid density functional theory method M06-2X. The feasibility of these decomposition reactions was calculated specifically for Meldrum’s acid, monomethyl Meldrum’s acid, dimethyl Meldrum’s acid (DMMA), methylene Meldrum’s acid (MeMA) and (dimethylamino)methylene Meldrum’s acid (DMAMeMA), based on previous experimental analysis. Their initial decomposition products (ketene or substituted ketenes with CO2 and acetone) were found endothermic requiring activation energies (∼180–––230 kJ/mol) at the decomposition temperatures of 280, 299, 301, 503 and 414 K respectively. The transition state theory calculations show MeMA and DMAMeMA are thermally more stable than other derivatives with low values for the rate constants and high energy barriers to produce initial decomposition products, suggesting they can produce fine carbon materials during carbonization. DMMA is found the best option for obtaining heteroatoms-doped carbon material formation.

Prediction of Feasibility of Polaronic OER on (110) Surface of Rutile TiO2.

The polaronic effects at the atomic level hold paramount significance for advancing the efficacy of transition metal oxides in applications pertinent to renewable energy. The lattice-distortion mediated localization of photoexcited carriers in the form of polarons plays a pivotal role in the photocatalysis. This investigation focuses on rutile TiO2, an important material extensively explored for solar energy conversion in artificial photosynthesis, specifically targeting the generation of green H2 through photoelectrochemical (PEC) H2O splitting. By employing Hubbard-U corrected and hybrid density functional theory (DFT) methods, we systematically probe the polaronic effects in the catalysis of oxygen evolution reaction (OER) on the (110) surface of rutile TiO2. Theoretical understanding of polarons within the surface, coupled with simulations of OER at distinct titanium (Ti) and oxygen (O) active sites, reveals diverse polaron formation energies within the lattice sites with strong preference for bulk and surface bridge (Ob) oxygen sites. Moreover, we provide the evidence for the facilitative role of polarons in OER. We find that hole polarons situated at the equatorial oxygen sites near the Ti-active site, along with bridge site hole polarons distal from the Ob active site yield a small reduction in OER overpotential by ~0.06 eV and ~0.12 eV, respectively. However, subsurface, equatorial, and bridge site hole polarons significantly reduce the Ti-active site OER overpotential by ~0.4 eV through the peroxo-type oxygen pathway. We also observe that the presence of hole polarons stabilizes the *OH, *O, and *OOH intermediate species compared to the scenario without hole polarons. Overall, this study provides a detailed mechanistic insight into polaron-mediated OER, offering a promising avenue for improving the catalytic activity of transition metal oxide-based photocatalysts catering to renewable energy requisites.

DLPNO-CCSD(T) and DFT study of the acetate-assisted C-H activation of benzaldimine at [RuCl2(p-cymene)]2: the relevance of ligand exchange processes at ruthenium(II) complexes in polar protic media.

To gain mechanistic insights into the acetate-assisted cyclometallations of arylimines promoted by [RuCl2(p-cymene)]2 in polar protic media, DFT geometry optimizations (with M06 and ωB97X-D3 functionals and the cc-pVDZ-PP[Ru] basis set) followed by DLPNO-CCSD(T)/CBS energy evaluations were performed using benzaldimine as a model substrate and methanol as the solvent (with CPCM or SMD models). The calculation results show that coordination of the imine to an acetate ruthenium precursor is followed by anion (chloride or acetate) dissociation as the rate-determining step of the process. H-Bonding of two explicit MeOH to the anion reduces the calculated activation energy to ca. 23 kcal mol-1, in good agreement with the experimental half-life at room temperature. Subsequent AMLA/CMD C-H activation of the intermediate cationic complexes is a faster, reversible process. Alternative reaction pathways involving neutral diacetate ruthenium complexes offer AMLA/CMD transition state structures of lower energy but are precluded due to higher energy barriers for the initial ligand exchange processes at ruthenium. Solvent assistance accelerates the final chloride/acetate exchange processes on the cycloruthenate intermediates, particularly when compression in the condensed phase is taken into consideration. The performance of six DFT functionals (with the aug-pVTZ-PP[Ru] basis set) was assessed using the DLPNO-CCSD(T)/CBS reference energies. Neutral diacetate ruthenium complexes were incorrectly predicted as being kinetically relevant when using hybrid DFT methods (PBE0-D3(BJ), M06-2X or ωB97M-V). Good agreement between the calculated barrier heights and our benchmark energy results was obtained by using double-hybrid DFT methods. PWPB95 with D3(BJ) or D4 dispersion energy corrections was found to be the most accurate (ΔG≠ MUE of ca. 1 kcal mol-1). This study may aid our understanding of and help with further experimental investigations of synthetically useful carboxylate-assisted C-H bond functionalizations involving (N,C)-cyclometallated (p-cymene)Ru(II) intermediate complexes in sustainable polar protic solvents.

In silico modelling of chelate stabilized tetrylene derivatives.

The steric and electronic effects of specific ligands can play crucial roles in stabilizing unsaturated tetrylene species. In this work, hybrid density functional theory (DFT) methods, quantum theory of atoms in molecules (QTAIM) investigations and natural bond orbital (NBO) calculations are employed to evaluate the stabilization of low-valent E(ii) centers (E = Si, Ge, Sn, Pb) through the chelating effect generated by an electron-rich ligand containing the P[double bond, length as m-dash]C-P[double bond, length as m-dash]X moiety (X = O or S). Based on several types of analyses, such as the bond dissociation energy (BDE) or the interplay between attractive (i.e., charge-transfer) and repulsive (i.e., Pauli-exchange) effects, we highlight that the stabilization energy induced by chelation is up to ca. 70 kcal mol-1 for silylenes, yet slightly decreases within the heavier analogues. Moreover, it is emphasized that chelate-stabilized silylenes can form highly stable hybrid metal-metalloid complexes with transition metals (e.g., gold). Due to push-pull effects occurring in the X→Si(ii)→Au fragment, the Si(ii)→Au bonding is significantly stronger than the X→Au, P(sp2)→Au or π(C[double bond, length as m-dash]P)→Au donor-acceptor bonds, which are potentially formed by the electron-rich P[double bond, length as m-dash]C-P[double bond, length as m-dash]X unit with the AuCl fragment. These findings are supported by energy decomposition analysis (EDA) calculations.

DAPD Set of Pd-Containing Diatomic Molecules: Accurate Molecular Properties and the Great Lengths to Obtain Them.

In the present study, we obtained reliable bond energy, bond length, and zero-point vibrational frequency for a set of diatomic Pd species (the DAPD set). It includes PdH, Pd2, and PdX (X = B, C, N, O, F, Al, Si, P, S, and Cl). Our highest-level protocol (W4X-L) represents scalar and spin-orbit relativistic, valence- and inner-valence correlated, extrapolated CCSDTQ(5) energy. The DAPD set of molecules is challenging for computational chemistry methods in different manners; for Pd2, the spin-orbit contribution to the bond energy is fairly large, whereas for PdC and PdSi, the post-CCSD(T) correlation components are considerable. The diverse range of requirements represents a significant challenge for lower-level methods. While density functional theory (DFT) methods generally yield good agreements for bond lengths and vibrational frequencies, large deviations are found for bond energies. In general, hybrid DFT methods are more accurate than nonhybrid functionals, but the agreement in individual cases varies. This illustrates the critical role that new high-quality reference data would play in the continual development of lower-cost methods.

The Nature of the Chemical Bonds of High-Valent Transition-Metal Oxo (M=O) and Peroxo (MOO) Compounds: A Historical Perspective of the Metal Oxyl-Radical Character by the Classical to Quantum Computations.

This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).

The stability of oxygen-centered radicals and its response to hydrogen bonding interactions.

The stability of various alkoxy/aryloxy/peroxy radicals, as well as TEMPO and triplet dioxygen (3 O2 ) has been explored at a variety of theoretical levels. Good correlations between RSEtheor and RSEexp are found for hybrid DFT methods, for compound schemes such as G3B3-D3, and also for DLPNO-CCSD(T) calculations. The effects of hydrogen bonding interactions on the stability of oxygen-centered radicals have been probed by addition of a single solvating water molecule. While this water molecule always acts as a H-bond donor to the oxygen-centered radical itself, it can act as a H-bond donor or acceptor to the respective closed-shell parent.

Streamlining and Optimizing Strategies of Electrostatic Parameterization.

Accurate characterization of electrostatic interactions is crucial in molecular simulation. Various methods and programs have been developed to obtain electrostatic parameters for additive or polarizable models to replicate electrostatic properties obtained from experimental measurements or theoretical calculations. Electrostatic potentials (ESPs), a set of physically well-defined observables from quantum mechanical (QM) calculations, are well suited for optimization efforts due to the ease of collecting a large amount of conformation-dependent data. However, a reliable set of QM ESP computed at an appropriate level of theory and atomic basis set is necessary. In addition, despite the recent development of the PyRESP program for electrostatic parameterizations of induced dipole-polarizable models, the time-consuming and error-prone input file preparation process has limited the widespread use of these protocols. This work aims to comprehensively evaluate the quality of QM ESPs derived by eight methods, including wave function methods such as Hartree-Fock (HF), second-order Møller-Plesset (MP2), and coupled cluster-singles and doubles (CCSD), as well as five hybrid density functional theory (DFT) methods, used in conjunction with 13 different basis sets. The highest theory levels CCSD/aug-cc-pV5Z (a5z) and MP2/aug-cc-pV5Z (a5z) were selected as benchmark data over two homemade data sets. The results show that the hybrid DFT method, ωB97X-D, combined with the aug-cc-pVTZ (a3z) basis set, performs well in reproducing ESPs while taking both accuracy and efficiency into consideration. Moreover, a flexible and user-friendly program called PyRESP_GEN was developed to streamline input file preparation. The restraining strengths, along with strategies for polarizable Gaussian multipole (pGM) model parameterizations, were also optimized. These findings and the program presented in this work facilitate the development and application of induced dipole-polarizable models, such as pGM models, for molecular simulations of both chemical and biological significance.

QM evaluation of the intramolecular aromatic π-π interactions of Ir(I) complex transition states

Major π-π interactions within Ir(I) complex Cp*Ir(I)PPh3 were quantum mechanically characterized using the major hybrid DFT, MP2 and CCSD(T) methods. Almost all the neighborhood π-π interactions are found to be slightly repulsive, and MP2 interaction energies range from -0.49 to 19.16 kcal/mol. The representative A-B π-π pair (See Figure 2) is further addressed using DFT methods. Interaction energies of the A-B π-π pair calculated from 6 hybrid DFT methods ranges from 0.52 kcal/mol by the SVWN method to 8.23 kcal/mol by the O3LYP method. The interaction energies calculated with unrestricted open shell MP2 theory can be up to 25% higher than those computed using restricted open shell MP2 method. Comparison between the DFT results and the benchmark CCSD(T) results shows that the DFT calculations of π-π interaction are reasonable but still need to be interpreted with caution. This study helps to pinpoint intramolecular aromatic π-π interactions of transition state for alkyl CH bond activation of saturated hydrocarbons and rational molecular design.

ωB97X-3c: A composite range-separated hybrid DFT method with a molecule-optimized polarized valence double-ζ basis set.

A new composite density functional theory (DFT) method is presented. It is based on ωB97X-V as one of the best-performing density functionals for the GMTKN55 thermochemistry database and completes the family of "3c" methods toward range-separated hybrid DFT. This method is consistently available for all elements up to Rn (Z = 1-86). Its further key ingredients are a polarized valence double-ζ (vDZP) Gaussian basis set, which was fully optimized in molecular DFT calculations, in combination with large-core effective core potentials and a specially adapted D4 dispersion correction. Unlike most existing double-ζ atomic orbital sets, vDZP shows only small basis set superposition errors (BSSEs) and can compete with standard sets of triple-ζ quality. Small residual BSSE effects are efficiently absorbed by the D4 damping scheme, which overall eliminates the need for an explicit treatment or empirical corrections for BSSE. Thorough tests on a variety of thermochemistry benchmark sets show that the new composite method, dubbed ωB97X-3c, is on par with or even outperforms standard hybrid DFT methods in a quadruple-zeta basis set at a small fraction of the computational cost. Particular strengths of this method are the description of non-covalent interactions and barrier heights, for which it is among the best-performing density functionals overall.

Gas phase H+, H3O+ and NH4+ affinities of oxygen-bearing volatile organic compounds; DFT calculations for soft chemical ionisation mass spectrometry.

Quantum chemistry calculations were performed using the density functional theory, DFT, to understand the structures and energetics of organic ions relevant to gas phase ion chemistry in soft chemical ionisation mass spectrometry analytical methods. Geometries of a range of neutral volatile organic compound molecules and ions resulting from protonation, the addition of H3O+ and the addition of NH4+ were optimised using the B3LYP hybrid DFT method. Then, the total energies and the normal mode vibrational frequencies were determined, and the total enthalpies of the neutral molecules and ions were calculated for the standard temperature and pressure. The calculations were performed for several feasible structures of each of the ions. The proton affinities of several benchmark molecules agree with the accepted values within ±4 kJ mol-1, indicating that B3LYP/6-311++G(d,p) provides chemical accuracy for oxygen-containing volatile organic compounds. It was also found that the binding energies of H3O+ and NH4+ to molecules correlate with their proton affinities. The results contribute to the understanding of ligand switching ion-molecule reactions important for secondary electrospray ionisation, SESI, and selected ion flow tube, SIFT, mass spectrometries.

Molecular Recognition of FDA-Approved Small Molecule Protein Kinase Drugs in Protein Kinases.

Protein kinases are key enzymes that catalyze the covalent phosphorylation of substrates via the transfer of the γ-phosphate of ATP, playing a crucial role in cellular proliferation, differentiation, and various cell regulatory processes. Due to their pivotal cellular role, the aberrant function of kinases has been associated with cancers and many other diseases. Consequently, competitive inhibition of the ATP binding site of protein kinases has emerged as an effective means of curing these diseases. Decades of intense development of protein kinase inhibitors (PKIs) resulted in 71 FDA-approved PKI drugs that target dozens of protein kinases for the treatment of various diseases. How do FDA-approved protein kinase inhibitor PKI drugs compete with ATP in their own binding pocket? This is the central question we attempt to address in this work. Based on modes of non-bonded interactions and their calculated interaction strengths by means of the advanced double hybrid DFT method B2PLYP, the molecular recognition of PKI drugs in the ATP-binding pockets was systematically analyzed. It was found that (1) all the FDA-approved PKI drugs studied here form one or more hydrogen bond(s) with the backbone amide N, O atoms in the hinge region of the ATP binding site, mimicking the adenine base; (2) all the FDA-approved PKI drugs feature two or more aromatic rings. The latter reach far and deep into the hydrophobic regions I and II, forming multiple CH-π interactions with aliphatic residues L(3), V(11), A(15), V(36), G(51), L(77) and π-π stacking interactions with aromatic residues F(47) and F(82), but ATP itself does not utilize these regions extensively; (3) all FDA-approved PKI drugs studied here have one thing in common, i.e., they frequently formed non-bonded interactions with a total of 12 residues L(3),V(11), A(15), K(17), E(24),V(36),T(45), F(47), G(51), L(77), D(81) and F(82) in the ATP binding. Many of those 12 commonly involved residues are highly conserved residues with important structural and catalytic functional roles. K(17) and E(24) are the two highly conserved residues crucial for the catalytic function of kinases. D(81) and F(82) belong to the DFG motif; T(45) was dubbed the gate keeper residue. F(47) is located on the hinge region and G(51) sits on the linker that connects the hinge to the αD-helix. It is this targeting of highly conserved residues in protein kinases that led to promiscuous PKI drugs that lack selectivity. Although the formation of hydrogen bond(s) with the backbone of the hinge gives PKI drugs the added binding affinity and the much-needed directionality, selectivity is sacrificed. That is why so many FDA-approved PKI drugs are known to have multiple targets. Moreover, off-target-mediated toxicity caused by a lack of selectivity was one of the major challenges facing the PKI drug discovery community. This work suggests a road map for future PKI drug design, i.e., targeting non-conserved residues in the ATP binding pocket to gain better selectivity so as to avoid off-target-mediated toxicity.

Cyclo[18]carbon‐A new class of electron acceptor for organic solar cells applications

The all-carboatomic ring, cyclo[18]carbon (C18), has the potential to act as an electron acceptor due to promising electronic and optical properties. In this study, we first illustrated the geometrical, electronic, and excited-state properties of C18 using various hybrid and long-range corrected density functional theory (DFT) methods. Further, we studied the nature of intermolecular interactions between dimers of C18 to gain insights into packing configurations of cyclo[18] in dimer and trimer configurations. Also, using the state-of-the-art DFT methods, we have reported the thorough characterization of the lowest excited-state (i.e. charge-transfer state) in various donor-acceptor model complexes based on pentacene and C18. We established an interplay between the molecular packing of C18 and pentacene molecules on the energy of charge transfer state. All these results could help in the designing of more efficient organic solar cells.

Fullerenes Pose a Strain on Hybrid Density Functional Theory.

The computational modeling of fullerenes plays a fundamental role in designing low-dimension carbon nanostructures. Nevertheless, the relative energies of fullerenes larger than C20 and C24 have not been comprehensively examined by means of highly accurate ab initio methods, for example, the CCSD(T) method. Here we report such an investigation for a diverse set of 29 C40 isomers. We calculate the energies of the C40 fullerenes using the G4(MP2) composite ab initio method, which approximates the CCSD(T) energy in conjunction with a triple-ζ-quality basis set (CCSD(T)/TZ). The CCSD(T)/TZ isomerization energies span 43.1-763.3 kJ mol-1. We find a linear correlation (R2 = 0.96) between the CCSD(T)/TZ isomerization energies and the fullerene pentagon signatures (P1 index), which reflect the strain associated with fused pentagon-pentagon rings. Using the reference CCSD(T)/TZ isomerization energies, we examine the relationship between the percentage of exact Hartree-Fock (HF) exchange in hybrid density functional theory (DFT) methods and the pentagon-pentagon strain energies. We find that the performance of hybrid DFT methods deteriorates with the pentagon-pentagon strain energy. This deterioration in performance becomes more pronounced with the inclusion of high amounts of HF exchange. For example, for B3LYP (20% HF exchange), the root-mean-square deviation (RMSD) relative to G4(MP2) increases from 8.9 kJ mol-1 for the low-strain isomers (P1 = 11) to 18.0 kJ mol-1 for the high-strain isomers (P1 > 13). However, for BH&HLYP (50% HF exchange) the RMSD increases from 23.0 (P1 = 11) to 113.2 (P1 > 13) kJ mol-1. A similar trend is observed for the M06/M06-2X pair of functionals. Namely, for M06 (27% HF exchange) the RMSD increases from 0.8 (P1 = 11) to 21.0 (P1 > 13) kJ mol-1, whereas for M06-2X (54% HF exchange) the RMSD increases from 16.7 (P1 = 11) to 77.7 (P1 > 13) kJ mol-1. Overall, we find that the strain associated with pentagon adjacency is an inherently challenging problem for hybrid DFT methods involving high amounts of HF exchange and that there is an inverse relationship between the optimal percentage of HF exchange and the pentagon-pentagon strain energy. For example, for BLYP the optimal percentages of HF exchange are 13% (P1 = 11), 10% (P1 = 12), 7.5% (P1 = 13), and 6% (P1 > 13).

Quantum chemical study of the hydrolysis of oxidized endogenous psychedelic N,N-dimethyltryptamine

The metabolism of biogenic N,N-dimethyltryptamine (DMT) starts with its enzymatic oxidation to indole-3-N,N-dimethylethaniminium cation (imDMT+) by monoamine oxidase (MAO). Following this reaction, a non-enzymatic hydrolysis occurs, yielding indole-3-acetaldehyde (IAL) and dimethylamine (DMA). In our in silico study, using a hybrid DFT method at the M06-2X/aug-cc-pVDZ level, we examined the necessity of explicit water molecules to accurately model the hydrolysis. Our calculations revealed that despite an implicit continuum solvent model, one or two explicit solvent molecules should be added to obtain meaningful activation barriers. Thus, we demonstrate that the hydrolysis can occur in a human cell without enzymatic assistance, paving the way for deeper understanding of DMT metabolism.

Removal of Heavy Metal Ions from Water Using Functionalized Carbon Nanosheet: A Density-Functional Theory Study

Background: Currently, a clean environment and human health have been regarded as the most important challenges of humankind; therefore, the ability of carbon-based nanomaterials as a new type of adsorbents in the removal of various pollutants from aqueous solutions has caused the widespread attention of research groups. Objectives: Carbon nanosheet (CNS) and its amine-functionalized derivatives were considered to remove heavy metal ions (HMIs) (zinc(II) ion [Zn2+]/cadmium(II) ion [Cd2+]/mercury(II) ion [Hg2+]) from water. Methods: This theoretical study in the framework of density-functional theory (DFT) was performed to gain insight into the HMIs (Zn2+/Cd2+/Hg2+) adsorption and removal by CNS and {(NH2)n-CNS (n = 1 - 2)}. The calculation level was hybrid DFT methods, such as wB97XD/lanl2dz. Results: Based on DFT calculations, more negative chemical potential value (-1.827 eV), high global softness (0.964 eV), low highest occupied molecular orbital (NH2)2-CNS lowest unoccupied molecular orbital-HMI gap, and more negative adsorption energy (range: -75 to -93 kJ mol-1) in (NH2)2-CNS demonstrated that this compound as a suitable adsorbent removed HMIs (Zn2+/Cd2+/Hg2+) from water than other nanosheets. The HMIs adsorption was confirmed by the natural bond orbital and quantum theory of atoms in molecules. Conclusions: The CNS and amine-functionalized structures have a good ability for HMIs (Zn2+/Cd2+/Hg2+) removal due to high chemical potential, electrophilicity, and adsorption energy.

The structure of anionic NbH n − (n = 2–15) clusters and their maximum hydrogen capacity

The geometrical structures of the lowest energy states of anionic niobium-doped hydrogen clusters are determined based on Crystal structure AnaLYsis by the Particle Swarm Optimization (CALYPSO) isomer search method. The geometries and relative stabilities are determined by the B3LYP hybrid density functional theory (DFT) method. The calculated results show that a niobium atom can absorb up to twelve hydrogen atoms in the anionic clusters. In exploring the stability, electronic properties, and photoelectron spectra of NbH n − clusters, the NbH12 − cluster is found to possesses a relatively high binding energy, high second-order differential energy, and a large HOMO–LUMO gap. The hydrogen storage density of the NbH12 − cluster is calculated as 11.5 wt%. The present findings enrich the database of transition metal doped hydrogen clusters and the insights are important in the pursuit of efficient hydrogen energy storage materials.

Calculation of vertical and adiabatic ionization potentials for some benzaldehydes using hybrid DFT, multilevel G3B3 and MP2 methods

The adiabatic ionization potential (AIP) and vertical ionization potential (VIP) for eight benzaldehydes were calculated using DFT (B3LYP, CAM-B3LYP and ωB97XD), MP2, G3B3 and CBS-QB3 methods. Basis set dependence was checked using four different basis sets 6-311G(d,p), 6-31+G(d,p), 6-311++G(2df,2p) and cc-pVTZ. Computed results were compared with the experimental data. Error metrics and statistical tools were used to validate the methods tested. Results showed an excellent agreement between the computed results and the experimental ones. Nature of the substituted groups and solvent effect on AIP values were investigated. Changes in the geometrical structures of cationic species were compared with neutral ones.

Determining the Band Alignment at the BiVO4|NiOOH Interface Using the Hybrid DFT Technique

It is important to understand the band offsets between semiconductors, which are crucial to determine the direction of electron transfer at the interfaces. Two methods are normally used to determine the direction from the first principles: alternating slabs put in contact (without empty spaces between them) and separate calculations for each material surface in the presence of vacuum spaces. The first method may introduce distortions due to insufficient epitaxial match, which may lead to bandgap changes, and the second may neglect electron transfer at the interface, which may be important in systems exhibiting very different average electronegativities. This can also imply a spill of electronic density into the vacuum spaces, which will not be present at real interfaces. Herein, both approaches were used to study the BiVO4/NiOOH interface, and the results were compared; the results are here relevant for photoelectrochemistry. The method is based on hybrid Density Functional Theory methods which give for the bulk phases Bandgap values that agree with the experimental ones (in one case, a value reflecting the theoretical value). The distances between the (hybrid DFT-derived) band positions and the corresponding profiles of the Hartree electrostatic potential were transferred to the interfaces. This helps determining the appropriate positions of the valence and conduction bands (as has been suggested by C. G. Van de Walle & R. M. Martin, Phys. Rev. B 1987, 35, 8154). It is ensured that the interfaces are nonpolar (Tasker’s criterion: P.W. Tasker, J. Phys. C: Solid State Phys. 1979, 12, 4977).

Symmetry Dilemma of Doubly Hybrid Density Functionals for Equilibrium Molecular Property Calculations.

In electronic structure theory, when an approximate wavefunction tends to artifactually break the symmetry of the exact Hamiltonian, the corresponding method is referred to as having a "symmetry dilemma" problem. Such types of artifacts were often reported when Hartree-Fock (HF) and the low-level post-HF methods were used, while the traditional Kohn-Sham density functional theory (KS-DFT) methods were usually found to be more resistant to this breakdown. In this work, we present a systematic study on the reliability of the doubly hybrid (DH) DFT methods for several violable cases. Almost all the commonly used B2PLYP-type (bDH) functionals are shown to have a severe "symmetry dilemma" problem and yield dramatically unreliable molecular properties, such as dipole moment, vibrational frequency, and static polarizability at the equilibrium geometry. A one-parameter bDH functional model study demonstrates that such a problem is a combined effect of the inappropriate portion of the HF exchange (over 50%) for the self-consistent field (SCF) calculation and the augmentation of the second-order perturbative contribution. It is remarkable that the XYG3-type (xDH) functionals show a good capability to resist the artifactual symmetry breaking and yield reliable molecular properties when the same critical cases are calculated. In the xDH, there are two functionals of different purposes, namely, the SCF functional and the energy functional, which have different amounts of the HF exchange and different portions of the correlation contributions. The success of the xDHs can be attributed to this flexibility in xDH construction to avoid using an improperly large portion of the HF exchange in the SCF functional. The insights gained in this work are of significance for the development of an improved DH functional.