Dispersion Interactions with Density-Functional Theory: Benchmarking Semiempirical and Interatomic Pairwise Corrected Density Functionals

Ultralow lattice thermal conductivity and high thermoelectric performance of monolayer KCuTe: a first principles study

Monolayer KCuTe is a new-type of two-dimensional (2D) semiconductor material with high carrier mobility and large power energy conversion efficiencies, suggesting its potential application in thermoelectric (TE) and photoelectric fields. Based on the density functional theory (DFT) and semiclassical Boltzmann transport equation, the electronic and phonon transport properties of monolayer KCuTe are systematically studied. Our results show that it possesses an ultralow lattice thermal conductivity value of nearly ∼0.13 W m⁻¹ K⁻¹ at 300 K, mainly attributed to its small phonon group velocity, large Grüneisen parameters, and strong phonon–phonon scattering. Furthermore, the intralayer opposite phonon vibrations greatly restrict the heat transport. Monolayer KCuTe shows an ideal direct band gap of ∼1.21 eV, and a high twofold degeneracy appearing at the Γ point gives a high Seebeck coefficient of ∼2070 μV K⁻¹, leading to high TE performance. Using the transport coefficients together with constant electron relaxation time, the figure of merit (ZT) can reach 2.71 at 700 K for the p-type doping, which is comparable to the well-known TE material SnSe (2.6 ± 0.3 at 935 K). Our theoretical studies may provide perspectives to TE applications of monolayer KCuTe and stimulate further experimental synthesis.

The excellent thermoelectric performance of monolayer KCuTe is discovered by first-principles study for the first time.

Two-Dimensional Li-Based Ternary Chalcogenides for Photocatalysis

Motivated by fundamental interest and practical applications, the investigations of two-dimensional photocatalysts are fascinating subjects in clean energy. Herein, we propose that two-dimensional Li-based ternary chalcogenides LiXY₂ (X = Al, Ga, In; Y = S, Se, Te) have intrinsic polarization and direct band gaps. Our results show that LiXY₂ materials possess optical absorption spectra covering the visible and ultraviolet range. We show that these materials possess extremely high electron mobility (∼10³ cm² V^-1 s^-1), providing great potential in overall water splitting. Furthermore, LiAlS₂ and LiGaS₂ can facilitate overall water splitting regardless of their energy gaps because of the large differences of surface electronic potentials of LiXY₂. Importantly, it is feasible to exfoliate the layered LiAlTe₂ from its bulk counterpart in experiments. Our findings open an exotic pathway to realizing promising photocatalytic applications in two-dimensional ternary materials.

Theoretical Study of Methane Storage in Cu24(m-BDC)24

Calculations on the Cu₂₄(m-BDC)₂₄ (m-BDC = 1,3-benzenedicarboxylate) polyoxometalate (POM) cage with 0, 12, 24, and 40 methane molecules inside were made using the M06 exchange/correlation functional. During filling of the cage with 40 CH₄ molecules, the 12 strongest binding CH₄ molecules are those to the coordination unsaturated sites (CUS) to the inwardly directed Cu(+2) centers via agostic interactions. The next 12 CH₄ molecules are less tightly bound followed by the next 16 CH₄ molecules with average binding energies of 8.27, 7.88, and 7.36 kcal/mol per CH₄, respectively. A section of the Cu₂₄(m-BDC)₂₄ cage was taken with the formula Cu₄(m-BDC)(BC)₆ (BC = benezenecarboxylate) in order to estimate zero-point, thermal, and entropy corrections of the larger cage. Estimating free energies at 1 bar, the Cu₂₄(m-BDC)₂₄ POM is predicted to lose 16, 12, and 12 CH₄ molecules at 67, 123, and 171 °C, respectively. The 40CH₄@Cu₂₄(m-BDC)₂₄ cage, which is isostructural to the main cavity of HKUST-1 with 40 CH₄ molecules inside, is predicted to have a loading of 224 cm³(STP) cm^-3 at 1 bar.

Comparative study of density functionals for the description of lithium-graphite intercalation compounds

The lack of description of van der Waals interactions in layered materials such as graphite and binary graphite intercalation compounds remains a main drawback of conventional density functional theory. Two fundamentally different approaches to overcome this issue are the employment of semiempirical dispersion correction scheme such as Grimme dispersion correction or nonlocal density functionals. We carefully compare these two approaches for the description of the geometric structure and the thermodynamic stability of pure graphite and Li-GICs at different lithium concentrations and stages. Based on the computed formation energies, we also evaluate the lithium-graphite intercalation potential. We find that PBE-D3(BJ) accurately reproduces the lattice parameters and the interlayer binding energy of graphite, although it underestimates the thermodynamic stability of stage-II Li-GICs mainly due to overbinding of carbon atoms in pure graphite. The nonlocal van der Waals functionals optB88-vdW, optB86b-vdW, and revPBE-vdW show a good agreement with experiments concerning stability of Li-GICs of different stages, although they overestimate the van der Waals interactions in graphite. The experimentally determined decreasing step-function behavior of Li-graphite intercalation potential can be qualitatively reproduced only by employing the revPBE van der Waals functional, whereas the other density functionals fail in the description. © 2019 Wiley Periodicals, Inc.

Encapsulation of anticancer drug doxorubicin inside dendritic macromolecular cavities: First-principles benchmarks

By using first-principles approaches based on Density Functional Theory, we explore the possibility of using dendritic macromolecular structures as carriers of the doxorubicin anticancer drug. In particular, we consider macromolecular cavities of different sizes composed of phenylene-, thiophene-, phenyl-cored thiophen- and thioazole-based dendrimers. The comparison between the optimized molecular geometries of the monomers and of the host-guest complexes reveals that only slight structural changes are observed in doxorubicin upon complexation. Also, the encapsulation energies for the host-guest complexes suggest that these systems are of potential use for pharmacology applications in vivo. The interaction of the guest doxorubicin with the macromolecular cavities exploits different types of weak intermolecular forces including σ, π and hydrogen bond interactions. The electronic structure of these complexes is discussed, with particular emphasis placed on the role of the charge distribution and the nature of the frontier molecular orbitals in the encapsulation process. Spectroscopic properties of these complexes are derived to facilitate their detection in laboratory and in vivo. These include IR vibrational frequencies, absorption wavelengths and relative oscillator strengths for the main transitions in the UV-Vis spectrum.

New Family of Two-Dimensional Ternary Photoelectric Materials

Screening unique two-dimensional (2D) materials with high mobility and applicable band gaps is motivated by not only the interest in basic science but also the practical applications for photoelectric materials. In this work, we have systematically studied a new family of 2D ternary quintuple layers (QLs), named ABC (A = Na, K, and Rb; B = Cu, Ag, and Au; C = S, Se, and Te). Our results indicate that the QLs of KCuTe, KAgS, KAgSe, KAuTe, RbCuTe, RbAgSe, and RbAgTe host direct band gaps. Moreover, KCuTe, RbCuTe, and RbAgTe QLs show extremely high mobilities of ∼10⁴ cm² V^-1 s^-1. Interestingly, the linear scaling between exciton binding energy and quasiparticle band gap for ABC QLs exhibits an unexpected deviation with the 1/4 law. In addition, KAgSe, KAgS, RbAgSe, and RbAgTe show outstanding power energy conversion efficiencies of up to 21.5%, suggesting that they are good potential donor materials. Our results provide many potential candidates for applications in photoelectric materials, which may be realized in experiments due to the possible exfoliation from their parent compounds.

Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study

We present the first, to the best of our knowledge, ab initio molecular dynamics (AIMD) investigation on three aqueous solutions where an abasic cyclic nucleotide model is solvated in the presence of distinct cations (i.e., Na⁺, K⁺ and Mg²⁺). We elucidate the typical modalities of interaction between those ionic species and the nucleotide moiety by first-principles numerical simulations, starting from an inner-shell binding configuration on a time scale of 100 ps (total simulation time of ∼600 ps). Whereas the strong "structure-maker" Mg²⁺ is permanently bound to one of the two oxygen atoms of the phosphate group of the nucleotide model, Na⁺ and K⁺ show binding times τ_b of 65 ps and 10-15 ps, respectively, thus reflecting their chemical nature in aqueous solutions. Furthermore, we qualitatively relate these findings to approximate free-energy barriers of the cations' unbinding obtained by means of exploratory well-tempered metadynamics. With the aim of shedding light on the features of commonly employed force-fields (FFs), classical MD simulations (almost 200 trajectories with a total simulation time of ∼18 μs) using the biomolecular AMBER FF are also reported. By choosing several combinations of the parametrization for the water environment (i.e., TIP3P, SPC/E and OPC) and cations (i.e., Joung-Cheatham, Li-Merz 12-6 and Li-Merz 12-6-4), we found significant differences in the radial distribution functions and residence times compared to the ab initio results. The Na⁺ and K⁺ ions wrongly show quasi-identical radial distribution functions and the Li & Merz 12-6-4 Lennard-Jones parameters for Mg²⁺ were found to be essential in quickly reaching the binding state consistent with AIMD.

Theoretical insight into the regioselective ring-expansions of bicyclic aziridinium ions

Transient bicyclic aziridinium ions are known to undergo ring-expansion reactions, paving the way to functionalized nitrogen-containing heterocycles. In this study, the regioselectivity observed in the ring-expansion reactions of 1-azoniabicyclo[n.1.0]alkanes was investigated from a computational viewpoint to study the ring-expansion pathways of two bicyclic systems with different ring sizes. Moreover, several nucleophiles leading to different experimental results were investigated. The effect of solvation was taken into account using both explicit and implicit solvent models. This theoretical rationalization provides valuable insight into the observed regioselectivity and may be used as a predictive tool in future studies.

Assessment of DFT for endohedral complexes' dipole moment: PNO-LCCSD-F12 as a reference method

We present a systematic evaluation of the performance of a wide range of exchange-correlation functionals and related dispersion correction schemes for the computation of dipole moments of endohedral complexes, formed through the encapsulation of an AB molecule (AB = LiF, HCl) inside carbon nanotubes (CNTs) of different diameter. The consistency and accuracy of (i) generalized gradient approximation, (ii) meta GGA, (iii) global hybrid, and (iv) range-separated hybrid density functionals are assessed. In total, 37 density functionals are tested. The results obtained using the highly accurate pair natural orbitals based explicitly correlated local coupled cluster singles doubles (PNO-LCCSD-F12) method of Werner and co-workers [Schwilk et al., J. Chem. Theory Comput., 2017, 13, 3650; Ma et al., J. Chem. Theory Comput., 2017, 13, 4871] with the aug-cc-pVTZ basis set serve as a reference. The static electric dipole moment is computed via the finite field response or, when possible, as the expectation value of the dipole operator. Among others, it is shown that functionals belonging to the class of range-separated hybrids, provide results closest to the coupled cluster reference data. In particular, the ωB97X as well as the M11 functional may be considered as a promising choice for computing electric properties of noncovalent endohedral complexes. On the other hand, the worst performance was found for the functionals which do not include the Hartree-Fock exchange. The analysis of both the coupled cluster and the DFT results indicates a strong coupling of dispersion and polarization that may also explain why lower level DFT methods, as well as Hartree-Fock and MP2, cannot yield dipole moments beyond a qualitative agreement with the higher order reference data. Interestingly, the much smaller and less systematically constructed basis sets of Pople of moderate size provide results of accuracy at least comparable with the extended Dunning's aug-cc-pVTZ basis set.

Tailoring van der Waals dispersion interactions with external electric charges

van der Waals (vdW) dispersion interactions strongly impact the properties of molecules and materials. Often, the description of vdW interactions should account for the coupling with pervasive electric fields, stemming from membranes, ionic channels, liquids, or nearby charged functional groups. However, this quantum-mechanical effect has been omitted in atomistic simulations, even in widely employed electronic-structure methods. Here, we develop a model and study the effects of an external charge on long-range vdW correlations. We show that a positive external charge stabilizes dispersion interactions, whereas a negative charge has an opposite effect. Our analytical results are benchmarked on a series of (bio)molecular dimers and supported by calculations with high-level correlated quantum-chemical methods, which estimate the induced dispersion to reach up to 35% of intermolecular binding energy (4 kT for amino-acid dimers at room temperature). Our analysis bridges electrostatic and electrodynamic descriptions of intermolecular interactions and may have implications for non-covalent reactions, exfoliation, dissolution, and permeation through biological membranes.

Stability of 2',3' and 3',5' cyclic nucleotides in formamide and in water: a theoretical insight into the factors controlling the accumulation of nucleic acid building blocks in a prebiotic pool

Synthesis of the first RNAs represents one of the cornerstones of the emergence of life. Recent studies demonstrated powerful scenarios of prebiotic synthesis of cyclic nucleotides in aqueous and formamide environments. This raised a question about their thermodynamic stability, a decisive factor determining their accumulation in a prebiotic pool. Here we performed ab initio molecular dynamics simulations at various temperatures in formamide and water to study the relative stabilities of the 2',3' and 3',5' isomers of cyclic nucleotides. The computations show that in an aqueous environment 2',3' cyclic nucleotides are more stable than their 3',5' counterparts at all temperatures up to the boiling point. In contrast, in formamide higher temperatures favor the accumulation of the 3',5' cyclic form, whereas below about 400 K the 2',3' cyclic form becomes more stable. The latter observation is consistent with a formamide-based origin scenario, suggesting that 3',5' cyclic nucleotides accumulated at higher temperatures subsequently allowed oligomerization reactions after fast cooling to lower temperatures. A statistical analysis of the geometrical parameters of the solutes indicates that thermodynamics of cyclic nucleotides in aqueous and formamide environments are dictated by the floppiness of the molecules rather than by the ring strain of the cyclic phosphodiester linkages.

Theoretical study of gas and solvent phase stability and molecular adsorption of noncanonical guanine bases on graphene

The gas and solvent phase stability of noncanonical (Gua)_n nucleobases is investigated in the framework of dispersion-corrected density functional theory (DFT). The calculated results strongly support the high tendency for the dimerization of (Gua)_n bases in both gas and solvent phases. An interplay between intermolecular and bifurcated H-bonds is suggested to govern the stability of (Gua)_n bases which bears a correlation with the description of dispersion correction terms employed in the DFT calculations. For example, a higher polarity is predicted for (Gua)_n bases by the dispersion-corrected DFT in contrast to the non-polar nature of (Gua)₃ and (Gua)₄ predicted by the hybrid meta-GGA calculations. This distinct variation becomes significant under physiological conditions as polar (Gua)_n is likely to exhibit greater stabilization in the gas phase compared to solvated (Gua)_n. Graphene acting as a substrate induces modification in base configurations via maximization of π-orbital overlap between the base and substrate. In solvent, the substrate-induced effects are further heightened with lowering of the dipole moments of (Gua)_n as also displayed by the corresponding isosurface of the electrostatic potential. The graphene-induced stability in both gas and solvent phases appears to fulfill one of the prerequisite criteria for molecular self-assembly. The DFT results therefore provide atomistic insights into the stability and molecular assembly of free-standing noncanonical (Gua)_n nucleobases which can be extended to understanding the self-assembly process of functional biomolecules on 2D materials for potential biosensing applications.

Modeling the archetype cysteine protease reaction using dispersion corrected density functional methods in ONIOM-type hybrid QM/MM calculations; the proteolytic reaction of papain

A proteolytic reaction of papain with a simple peptide model substrate N-methylacetamide has been studied. Our aim was twofold: (i) we proposed a plausible reaction mechanism with the aid of potential energy surface scans and second geometrical derivatives calculated at the stationary points, and (ii) we investigated the applicability of the dispersion corrected density functional methods in comparison with the popular hybrid generalized gradient approximations (GGA) method (B3LYP) without such a correction in the QM/MM calculations for this particular problem. In the resting state of papain the ion pair and neutral forms of the Cys-His catalytic dyad have approximately the same energy and they are separated by only a small barrier. Zero point vibrational energy correction shifted this equilibrium slightly to the neutral form. On the other hand, the electrostatic solvation free energy corrections, calculated using the Poisson-Boltzmann method for the structures sampled from molecular dynamics simulation trajectories, resulted in a more stable ion-pair form. All methods we applied predicted at least a two elementary step acylation process via a zwitterionic tetrahedral intermediate. Using dispersion corrected DFT methods the thioester S-C bond formation and the proton transfer from histidine occur in the same elementary step, although not synchronously. The proton transfer lags behind (or at least does not precede) the S-C bond formation. The predicted transition state corresponds mainly to the S-C bond formation while the proton is still on the histidine Nδ atom. In contrast, the B3LYP method using larger basis sets predicts a transition state in which the S-C bond is almost fully formed and the transition state can be mainly featured by the Nδ(histidine) to N(amid) proton transfer. Considerably lower activation energy was predicted (especially by the B3LYP method) for the next amide bond breaking elementary step of acyl-enzyme formation. Deacylation appeared to be a single elementary step process in all the methods we applied.

Role of surface adsorption in tuning the properties of black phosphorus

The synergetic effect of O₂and H₂O during the oxidation of black phosphorus (BP) at the atomic level is revealed, and the effects of H₂O and/or O₂on the properties of BP are also investigated.

A variable-temperature X-ray diffraction and theoretical study of conformational polymorphism in a complex organic molecule (DTC)

Two conformational crystal polymorphs of 3-diethylamino-4-(4-methoxyphenyl)-1,1-dioxo-4H-1λ⁶,2-thiazete-4-carbonitrile (DTC) have been analyzed in the 100 K-room temperature range by single crystal X-ray diffraction and high quality DFT calculations. DTC has strongly polar as well as aliphatic substituents but no hydrogen bonding groups, and thus qualifies as a test molecule for the relative importance of electrostatic vs. dispersion–repulsion terms. The two polymorphs have the same P2₁/n space group and differ by a flipping of the –OCH₃ group, the two conformations being almost equi-energetic and separated by a low barrier. The system is monotropic in the observed temperature range with nearly identical thermal expansion coefficients and energy–temperature slopes, one phase consistently predicted to be more stable in agreement with the relative ease of appearance. Energy decompositions show that the electrostatic term is dominant and stabilizes with decreasing temperature. Dispersion and repulsion show the expected behavior, the former becoming more stabilizing at lower temperature in contrast with increasing repulsion at higher density. Absolute values and trends are very similar in the two phases, explaining the small total energy difference. Geometrical analyses of intermolecular contacts using fingerprint plots, as well as the study of molecular dipole moments as a function of T in the framework of the Quantum Theory of Atoms in Molecules, reveal more details of phase stability.

A small conformational change in the asymmetric unit has a significant effect on how non-covalent interactions determine (i) the crystal packing and (ii) the effect of T on the relative balance of electrostatics and dispersion–repulsions.

Valence electronic structure of [EMIM][BF4] ionic liquid: photoemission and DFT+D study

The ultraviolet photoelectron spectrum (UPS) of the [EMIM][BF₄] ionic liquid was recorded and analyzed. Together with the gas-phase UPS spectrum of the [EMIM][BF₄] vapor and ab initio calculation methods, detailed insight into the electronic structure of this simple ionic liquid is possible. The low binding energy tail in the UPS spectrum is about 7.4 eV, in agreement with previous estimations of the HOMO–LUMO gap of the [EMIM][BF₄] ion-pair. The bulk ab initio calculations are able to explain most of the features in the spectrum. However, DFT consistently lacks accuracy in the description of the top of the valence band. The dispersion corrected PBE calculation (PBE-D3) did offer very good agreement with the experimental structure, but the recently-developed vdW-DF functionals C09, optPBE, optB88 and CX were found to offer the best agreement in terms of the electronic structure.

The ultraviolet photoelectron spectrum (UPS) of the [EMIM][BF₄] ionic liquid was recorded and compared to previously measured vapor phase UPS spectrum.

Conformational and structural stability of the single molecule and hydrogen bonded clusters of para aminobenzoic acid in the gas and solution phases

The structures of α- and β-para aminobenzoic acid are deconstructed into their hydrogen bonding molecular structural building blocks, where they are analysed usingab initioquantum mechanical calculations of their conformation and cluster stability in solution.

Crystalline packing in pentacene-like organic semiconductors

Combining Hirshfeld surface analysis with single molecule electrostatic property calculations allows rationalizing crystalline packing motifs of organic semiconductors.

Non-covalent interactions for carbonaceous materials: impacts of doping, curving and their combination

Non-covalent interactions and their influencing factors are addressed after studying the binding of ion pairs with graphene (doped, curved and combination).

Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals

We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm-3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.

First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications

Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.

Validation of density functionals for pancake-bonded π-dimers; dispersion is not enough

π-Stacking pancake bonding between radicals poses special challenges to density functional theories (DFTs) due to their shorter than van der Waals contact distances, their multireference singlet ground states and the concurrently important dispersion interactions.

Ionic diffusion and proton transfer in aqueous solutions of alkali metal salts

We report on a series of ab initio molecular dynamics investigations on LiCl, NaCl, and KCl aqueous solutions under the effect of static electric fields.

Enhanced hydrogen adsorption on Li-coated B12C6N6

The hydrogen storage property of Li-coated B₁₂C₆N₆ is investigated by density functional theory calculations. B₁₂C₆N₆ is an electron deficient fullerene. Li atoms can be strongly bound to this cage by donating their valance electrons to the virtual 2p orbitals of carbon in the cluster. The binding energy (-2.90 eV) is much larger than the cohesive energy (1.63 eV) of bulk Li, and it prevents the Li atoms from aggregation. The coated Li atoms have large positive charges and the adsorbed hydrogen molecules can be moderately polarized by the Li⁺ ions. The computation shows that each Li atom coated on B₁₂C₆N₆ can hold 2-3 H₂ molecules with adsorption energies in the range of 0.21-0.24 eV/H₂. The B₁₂C₆N₆Li₈ can adsorb 16 H₂ and achieve a gravimetric hydrogen density of 8.63 wt. %. The present results indicate that alkali-metal atoms coated on electron deficient fullerenes can serve as hydrogen storage materials that can operate at ambient temperatures with high recycling storage capacity.

Failure of Density Functional Dispersion Correction in Metallic Systems and Its Possible Solution Using a Modified Many-Body Dispersion Correction

Previous density functional dispersion corrections to density functional theory lead to an unphysical description of metallic systems, as exemplified by alkali and alkaline earth compounds. We demonstrate that it is possible to remedy this limitation by including screening effects into the form of interacting smeared-out dipoles in the many-body expansion of the interaction. Our new approach, called the coupled fluctuating smeared dipole model, describes equally well noncovalent systems, such as molecular pairs and crystals, and metallic systems.

Benchmarking Density Functional Theory Based Methods To Model NiOOH Material Properties: Hubbard and van der Waals Corrections vs Hybrid Functionals

NiOOH has recently been used to catalyze water oxidation by way of electrochemical water splitting. Few experimental data are available to rationalize the successful catalytic capability of NiOOH. Thus, theory has a distinctive role for studying its properties. However, the unique layered structure of NiOOH is associated with the presence of essential dispersion forces within the lattice. Hence, the choice of an appropriate exchange-correlation functional within Density Functional Theory (DFT) is not straightforward. In this work, we will show that standard DFT is sufficient to evaluate the geometry, but DFT+U and hybrid functionals are required to calculate the oxidation states. Notably, the benefit of DFT with van der Waals correction is marginal. Furthermore, only hybrid functionals succeed in opening a bandgap, and such methods are necessary to study NiOOH electronic structure. In this work, we expect to give guidelines to theoreticians dealing with this material and to present a rational approach in the choice of the DFT method of calculation.