Wavefunction theory and density functional theory analysis of ground and excited electronic states of TaB and WB

Several low-lying electronic states of TaB and WB molecules were studied using ab initio multireference configuration interaction (MRCI), Davidson corrected MRCI (MRCI+Q), and coupled cluster singles doubles and perturbative triples [CCSD(T)] methods. Their full potential energy curves (PECs), equilibrium electron configurations, equilibrium bond distances (res), dissociation energies (Des), excitation energies (Tes), harmonic vibrational frequencies (ωes), and anharmonicities (ωexes) are reported. The MRCI dipole moment curves (DMCs) of the first 5 electronic states of both TaB and WB are also reported and the equilibrium dipole moment (μ) values are compared with the CCSD(T) μ values. The most stable 13Π (1σ22σ23σ11π3) and 15Δ (1σ22σ23σ11π21δ1) electronic states of TaB lie close in energy with ∼62 kcal mol-1De with respect to the Ta(4F) + B(2P) asymptote. However, spin-orbit coupling effects make the 15Δ0+ state the true ground state of TaB. The ground electronic state of WB (16Π) has the 1σ22σ13σ11π31δ2 electron configuration and is followed by the excited 16Σ+ and 14Δ states. Finally, the MRCI De, re, ωe, and ωexe values of the 13Π state of TaB and 16Π and 14Δ states of WB are used to assess the density functional theory (DFT) errors on a series of exchange-correlation functionals that span multiple-rungs of the Jacob's ladder of density functional approximations (DFA).

Ab initio electronic structure analysis of ground and excited states of HfN0,

High-level ab initio electronic structure analysis of third-row transition metal (TM)-based diatomic species is challenging and has been perpetually lagging. In this work, fourteen and eighteen electronic states of HfN and HfN+ respectively are studied, employing multireference configuration interaction (MRCI) and coupled cluster singles doubles and perturbative triples [CCSD(T)] theories under larger correlation-consistent basis sets. Their potential energy curves (PECs), energetics, and spectroscopic parameters are reported. Core electron correlation effects on their properties are also investigated. Chemical bonding patterns of several low-lying electronic states are introduced based on the equilibrium electron configurations. The ground state of HfN (X2Σ+) has the 1σ22σ23σ11π4 electronic configuration, and the ionization of the 3σ1 electron produces the ground state of HfN+ (X1Σ+). Ground states of both HfN and HfN+ are triple bonded in nature and bear 124.86 and 109.10 kcal mol-1 binding energies with respect to their ground state fragments. The findings of this work agree well with the limited experimental literature available and provide useful reference values for future experimental analysis of HfN and HfN+.

Gas-phase and solid-state electronic structure analysis and DFT benchmarking of HfCO

Ab initio multi-reference configuration interaction (MRCI) and coupled cluster singles doubles and perturbative triples [CCSD(T)] levels of theory were used to study ground and excited electronic states of HfCO. We report potential energy curves, dissociation energies (De), excitation energies, harmonic vibrational frequencies, and chemical bonding patterns of HfCO. The 3Σ- ground state of HfCO has an 1σ22σ21π2 electron configuration and a ∼30 kcal mol-1 dissociation energy with respect to its lowest-energy fragments Hf(3F) + CO(X1Σ+). We further evaluated the De of its isovalent HfCX (X = S, Se, Te, Po) series and observed that they increase linearly from the lighter HfCO to the heavier HfCPo with the dipole moment of the CX ligand. The same linear relationship was observed for TiCX and ZrCX. We utilized the CCSD(T) benchmark values of De, excitation energy, and ionization energy (IE) values to evaluate density functional theory (DFT) errors with 23 exchange-correlation functionals spanning GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, and double-hybrid functional families. The global GGA hybrid B3LYP and range-separated hybrid ωB97X performed well at representing the ground state properties of HfCO (i.e., De and IE). Finally, we extended our DFT analysis to the interaction of a CO molecule with a Hf surface and observed that the surface chemisorption energy and the gas-phase molecular dissociation energy are very similar for some DFAs but not others, suggesting moderate transferability of the benchmarks on these molecules to the solid state.

Transition metal oxide complexes as molecular catalysts for selective methane to methanol transformation: any prospects or time to retire?

Transition metal oxides have been extensively used in the literature for the conversion of methane to methanol. Despite the progress made over the past decades, no method with satisfactory performance or economic viability has been detected. The main bottleneck is that the produced methanol oxidizes further due to its weaker C-H bond than that of methane. Every improvement in the efficiency of a catalyst to activate methane leads to reduction of the selectivity towards methanol. Is it therefore prudent to keep studying (both theoretically and experimentally) metal oxides as catalysts for the quantitative conversion of methane to methanol? This perspective focuses on molecular metal oxide complexes and suggests strategies to bypass the current bottlenecks with higher weight on the computational chemistry side. We first discuss the electronic structure of metal oxides, followed by assessing the role of the ligands in the reactivity of the catalysts. For better selectivity, we propose that metal oxide anionic complexes should be explored further, while hydrophylic cavities in the vicinity of the metal oxide can perturb the transition-state structure for methanol increasing appreciably the activation barrier for methanol. We also emphasize that computational studies should target the activation reaction of methanol (and not only methane), the study of complete catalytic cycles (including the recombination and oxidation steps), and the use of molecular oxygen as an oxidant. The titled chemical conversion is an excellent challenge for theory and we believe that computational studies should lead the field in the future. It is finally shown that bottom-up approaches offer a systematic way for exploration of the chemical space and should still be applied in parallel with the recently popular machine learning techniques. To answer the question of the title, we believe that metal oxides should still be considered provided that we change our focus and perform more systematic investigations on the activation of methanol.

Electronic structure of the ground and excited states of neutral and charged silicon hydrides, SiHx0/+/-, x = 1-4

The ground and excited electronic states of the titled species are investigated with multi-reference configuration interaction and diffuse basis sets. We found that in addition to the valence orbitals, the inclusion of the 4s, 4p, and especially 3d orbitals (although with minimal population) of silicon in the active space of the reference complete active space self-consistent field wavefunction are necessary for the proper convergence of the calculations. We also demonstrate that the aug-cc-pVTZ basis set provides quite accurate results compared to both larger basis sets and basis set limit results at a lower computational cost. The excited states involve excitations within the 3s and 3p orbitals of silicon (especially for the mono- and di-hydrides), followed by excitations from the Si-H bonding orbitals to either silicon valence or Rydberg (4s, 4p) orbitals. The number of electronic states per energy unit decrease as we add hydrogen atoms, and the first excited state of SiH₄ is at 9.0 eV and leads to SiH₃ + H. All species have stable ground state structures with all hydrogen atoms bound to silicon, except for SiH₄⁺ and SiH₄^-. The former dissociates to SiH₂⁺ + H₂, while the latter loses an electron or can dissociate forming H₂ as well.

Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf-O or Hf-B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X⁴Σ^-) and HfO(X¹Σ⁺) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D_e) from HfO to HfPo and the singlet-triplet energy gap (ΔE_S-T) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange-correlation functionals spanning multiple rungs of "Jacob's ladder." When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.

Electronic Structure of RhO2+, Its Ammoniated Complexes (NH3)1-5RhO2+, and Mechanistic Exploration of CH4 Activation by Them

High-level electronic structure calculations are initially performed to investigate the electronic structure of RhO²⁺. The construction of potential energy curves for the ground and low-lying excited states allowed the calculation of spectroscopic constants, including harmonic and anharmonic vibrational frequencies, bond lengths, spin-orbit constants, and excitation energies. The equilibrium electronic configurations were used for the interpretation of the chemical bonding. We further monitored how the Rh-O bonding scheme changes with the gradual addition of ammonia ligands. The nature of this bond remains unaffected up to four ammonia ligands but adopts a different electronic configuration in the pseudo-octahedral geometry of (NH₃)₅RhO²⁺. This has consequences in the activation mechanism of the C-H bond of methane by these complexes, especially (NH₃)₄RhO²⁺. We show that the [2 + 2] mechanism in the (NH₃)₄RhO²⁺ case has a very low energy barrier comparable to that of a radical mechanism. We also demonstrate that methane can coordinate to the metal in a similar fashion to ammonia and that knowledge of the electronic structure of the pure ammonia complexes provides qualitative insights into the optimal reaction mechanism.

Radical abstraction vs. oxidative addition mechanisms for the activation of the SH, OH, and CH bonds using early transition metal oxides

Quantum chemical calculations are performed to study the S-H, O-H, and C-H bond activation of H2S, H2O, and CH4 by bare and ligated ZrO+ and NbO+ units. These representative oxides bear low energy oxo and higher energy oxyl units. S-H and C-H bonds are readily activated by metal oxyl states (radical mechanism), but the O-H bond is harder to activate with either the oxyl or oxo states. Our results suggest that known practices for the C-H bond activation can be applied to S-H, but not to O-H bonds. The identified trends are rationalized in terms of the HS-H, HO-H, and H3C-H dissociation energies to the homolytic or heterolytic fragments. We also found that these dissociation energies drop to about half after coordination of H2S or H2O to the metal oxide unit. In addition, chlorine ligands are shown to stabilize the higher energy oxyl states of the transition metal oxygen unit enhancing the reactivity of the formed complexes.

Ab initio investigation of the ground and excited states of RuO+,0,- and their reaction with water

High-level quantum chemical calculations on RuO0,± elucidate the electronic structure of their low-lying electronic states. For thirty-two states, we report the electronic configurations, bond lengths, vibrational frequencies, spin-orbit splittings, and excitation energies. The electronic states of RuO can be generated from those of RuO+ by adding one electron to the σ non-bonding orbital closely resembling the 5s atomic orbital of Ru. The ground states for RuO and RuO- are clearly identified as 5Δ and 4Δ, but the two states (4Δ and 2Π) compete for RuO+. The difficulty of calculations is revealed by our small binding energies compared to the experimental values. In addition, we studied the reaction of the three species with water in their ground and selected low-lying electronic states. We found a consistent decrease of the activation energy barriers and higher exothermicity as we add electrons to the system. RuO- is found to facilitate the reaction for both kinetic and thermodynamic reasons.

Weak-field ligands enable inert early transition metal oxides to convert methane to methanol: the case of ZrO

Zirconium monoxide, ZrO, was studied by multi-reference configuration interaction (MRCI) and coupled cluster methods using large basis sets in conjunction with effective core potentials. Complete potential energy curves were constructed and bonding patterns are proposed for several electronic states. Numerical results include accurate equilibrium bond lengths, harmonic vibrational frequencies, anharmonicities, excitation energies, dipole moments, and binding energies for both ground and excited states. The application of a ZrO unit as the catalytic center for methane activation is explored through the reaction ZrO + CH₄→ Zr + CH₃OH. Optimal density functional structures combined with single-point MRCI energy calculations are obtained for the complete reaction pathway. It is found that the lower energy singlet and triplet multiplicities (oxo states) favor the [2+2] mechanism and the higher energy quintets (oxyl states) favor the radical mechanism, which is overall more efficient in producing methanol. We finally suggest proper ligands that stabilize the oxyl states. These include halogens or other weak-field ligands, which finally convert the inert early transition metal oxide units to efficient methane-to-methanol catalysts.

O-H and C-H Bond Activations of Water and Methane by RuO2+ and (NH3)RuO2+: Ground and Excited States

Investigation of the ground and excited states of RuO²⁺ is carried out using multireference quantum chemical methodologies. The electronic structure is explored in detail, and accurate spectroscopic constants for 12 states are reported. Although ruthenium belongs to the same group as iron, the ground state of RuO²⁺ is ¹Σ⁺ with a strong oxo character as opposed to the ³Δ of FeO²⁺ with primarily oxyl character. To see the effect of the different electronic structure of RuO²⁺ on the O-H and C-H bond activation processes, we studied its reaction with one water or methane molecule. Reaction energies and activation barriers are given for six low-lying electronic states of singlet, triplet, and quintet spin multiplicities. It is found that the higher-energy quintet state (⁵Σ⁺) provides the lowest activation energies and is the same state responsible for the C-H activation for FeO²⁺ complexes. The reason is attributed to its weaker metal-oxygen bond (longer bond length), which is "prepared" to be activated at the same time with the O-H and C-H bonds. The effect of an ammonia ligand in the chemical activity is also discussed.

Electronic and geometric structure analysis of neutral and anionic metal nitric chalcogens: The case of MNX series (M=Li, Na, Be and X=O, S, Se, Te)

Coupled cluster and multireference configuration approaches are employed to study the electronic and geometric structures of mono-coordinated complexes of lithium, sodium, and beryllium with nitric oxide and its isovalent NS, NSe, and NTe species. Ground and low-lying excited states were examined for both linear-bonded and side-bonded isomers. We show that the ionic M⁺ NX^- (M=Li, Na, Be and X=O, S, Se, Te) picture is a more natural representation and can account for the symmetry of the low-lying electronic states as Σ^- , Δ, and Σ⁺ , the smaller excitation energies and the larger binding energies for heavier X. An additional electron binds to the positively charged Li and Na terminal creating stable anions. The electron affinity (EA) of LiNX and NaNX species is in the 0.5-0.8 eV range. Despite the negative EA of beryllium and the very small EA of NO, the BeNO molecule has an EA of ~1.0 eV, which is increased to ~1.5 eV for the heavier BeNX species. This is attributed to the fact that the additional electron goes to the beryllium end for BeNO but to a π(MN)π*(NX) orbital of the rest species. Our accurate results contradict previous findings and serve as a guide for future experimental studies. © 2019 Wiley Periodicals, Inc.

Quantum chemical calculations on NbO and its reaction with methane: ground and excited electronic states

Our high-level calculations show that high-spin NbO electronic states facilitate the methane to methanol transformation via a very efficient radical mechanism, as opposed to the [2+2] mechanism observed for the rest of the low-lying states.

Ab initio calculations on the ground and excited electronic states of neutral and charged palladium monoxide, PdO0,+,−

Multi-reference configuration interaction and coupled cluster calculations were carried out for the ground and several low-lying excited electronic states for PdO, PdO⁺, and PdO⁻. The photoelectron spectrum peaks of PdO were assigned.