Gaussian basis sets of quadruple zeta quality for potassium through xenon: application in CCSD(T) atomic and molecular property calculations

Correlation consistent basis sets designed for density functional theory: Third-row atoms (Ga-Br)

The correlation consistent basis sets (cc-pVnZ with n = D, T, Q, 5) for the Ga-Br elements have been redesigned, tuning the sets for use for density functional approximations. Steps to redesign these basis sets for an improved correlation energy recovery and efficiency include truncation of higher angular momentum functions, recontraction of basis set coefficients, and reoptimization of basis set exponents. These redesigned basis sets are compared with conventional cc-pVnZ basis sets and other basis sets, which are, in principle, designed to achieve systematic improvement with respect to increasing basis set size. The convergence of atomic energies, bond lengths, bond dissociation energies, and enthalpies of formation to the Kohn-Sham limit is improved relative to other basis sets where convergence to the Kohn-Sham limit is typically not observed.

Electron group localization in atoms and molecules

Partitioning atomic and molecular charge densities in non overlapping chemically significant regions is a challenging problem for quantum chemists. The present method aims to build a tool enabling the determination of ``good boundaries' with the help of elementary statistical methods or information theory. This is done by minimizing an objective function with respect to the boundaries of the localization regions, the choice of this function being guided by a clarity requirement. With the sum of the indexes ofdispersion ($\Sigma D$) or the mutual information ($MI$) as objective function,the method yield partitions in good agreement with the Aufbau rules for Li-Rn atoms and with Lewis's pairing model for molecules.

Ag2O versus Cu2O in the Catalytic Isomerization of Coordinated Diaminocarbenes to Formamidines: A Theoretical Study

DFT theoretical calculations for the Ag₂O-induced isomerization process of diaminocarbenes to formamidines, coordinated to Mn(I), have been carried out. The reaction mechanism found involves metalation of an N-H residue of the carbene ligand by the catalyst Ag₂O and the formation of a key transition state showing a μ-η²:η² coordination of the formamidinyl ligand between manganese and silver, which allows a translocation process of Mn(I) and silver(I) ions between the carbene carbon atom and the nitrogen atom, before the formation of the formamidine ligand is completed. Calculations carried out using Cu₂O as a catalyst instead of Ag₂O show a similar reaction mechanism that is thermodynamically possible, but highly unfavorable kinetically and very unlikely to be observed, which fully agrees with experimental results.

A detailed analysis of the spin-crossover reaction of H2S binding to heme and the six-coordinated FeP(Im)-HS- porphyrin complex

The potential energy surfaces of the H₂S binding to iron-porphyrin (FeP) with the imidazole (Im) ligand via intersystem crossings are investigated by using density functional theory. The minimum energy intersystem crossing point (MEISCP) between the quintet and triplet states (MEISCP_TQ) for the Fe(II)P(Im)-H₂S complex is located at a Fe-S distance of 3.39 Å with only 1.1 kcal/mol above the quintet state minimum. The second spin-crossover point, where a change from the triplet to the singlet state occurs, comes at a much shorter Fe-S distance of 2.79 Å, and the MEISCP_ST is located at 3.7 kcal/mol above the triplet state minimum. The nature of the chemical bonding along the Fe-S reaction coordinate from the ground state singlet to the quintet state along the path to the separated species is analyzed. An inspection of the vibrational modes reveals that the largest contribution to the triplet-quintet transition around the quintet and triplet state minimum comes from the symmetric shrinking of the pyrrole units of the porphyrin ring, indicating that the related reaction coordinate plays a main role in the intersystem crossing. The fully optimized structures of the Fe(II)P(Im)-HS^- complex corresponding to three different spin multiplicities (M = 1, 3, 5) are characterized by a bent Fe-H-S conformation. The binding of the hydrosulfide anion to Fe(II)P(Im) in the quintet state induces a 0.2 Å displacement of the Fe atom out of the nitrogen porphyrin (N_pyr) plane. The fully optimized structure of the ground state of Fe(II)P(Im)-HS^- agrees well with experimental data for the corresponding heme models.

In search of stable singlet metalla-N-heterocyclic carbenes (MNHCs): a contribution from theory

Theoretical studies predict that the stability of the singlet state of metalla-N-heterocyclic carbenes (MNHCs) is strongly influenced not only by the substituents at the α-nitrogen atoms and the nature of ligands at the transition metal center but also by the substituents at the carbenic backbone. All the MNHCs were found to have a stable singlet ground state and the computed ΔE_S-T values for some of the MNHCs were found to be significantly large and lie within the range of experimentally known carbenes (31.0-84.0 kcal mol^-1). Furthermore, they were found to have superior σ-donation ability to conventional NHCs. The calculated proton affinities, gallium pyramidalization, pK_a and nucleophilicity index values for the MNHCs were found to be in good agreement with their calculated electron donation ability.

Performance of density functional theory for describing hetero-metallic active-site motifs for methane-to-methanol conversion in metal-exchanged zeolites

Methane-to-methanol conversion (MMC) can be facilitated with high methanol selectivities by copper-exchanged zeolites. There are however two open questions regarding the use of these zeolites to facilitate the MMC process. The first concerns the possibility of operating the three cycles in the stepwise MMC process by these zeolites in an isothermal fashion. The second concerns the possibility of improving the methanol yields by systematic substitution of some copper centers in these active sites with other earth-abundant transition metals. Quantum-mechanical computations can be used to compare methane activation by copper oxide species and analogous mixed-metal systems. To carry out such screening, it is important that we use theoretical methods that are accurate and computationally affordable for describing the properties of the hetero-metallic catalytic species. We have examined the performance of 47 exchange-correlation density functionals for predicting the relative spin-state energies and chemical reactivities of six hetero-metallic [M-O-Cu]²⁺ and [M-O₂ -Cu]²⁺ , (where MCo, Fe, and Ni), species by comparison with coupled cluster theory including iterative single, double excitations as well as perturbative treatment of triple excitations, CCSD(T). We also performed multireference calculations on some of these systems. We considered two types of reactions (hydrogen addition and oxygen addition) that are relevant to MMC. We recommend the use of τ-HCTH and OLYP to determine the spin-state energy splittings in the hetero-metallic motifs. ωB97, ωB97X, ωB97X-D3, and MN15 performed best for predicting the energies of the hydrogen and oxygen addition reactions. In contrast, local, and semilocal functionals do poorly for chemical reactivity. Using [Fe-O-Cu]²⁺ as a test, we see that the nonlocal functionals perform well for the methane CH activation barrier. In contrast, the semilocal functionals perform rather poorly. © 2018 Wiley Periodicals, Inc.

Distinguishing Bonds

The energy change per electron in a chemical or physical transformation, ΔE/n, may be expressed as Δχ̅ + Δ(VNN + ω)/n, where Δχ̅ is the average electron binding energy, a generalized electronegativity, ΔVNN is the change in nuclear repulsions, and Δω is the change in multielectron interactions in the process considered. The last term can be obtained by the difference from experimental or theoretical estimates of the first terms. Previously obtained consequences of this energy partitioning are extended here to a different analysis of bonding in a great variety of diatomics, including more or less polar ones. Arguments are presented for associating the average change in electron binding energy with covalence, and the change in multielectron interactions with electron transfer, either to, out, or within a molecule. A new descriptor Q, essentially the scaled difference between the Δχ̅ and Δ(VNN + ω)/n terms, when plotted versus the bond energy, separates nicely a wide variety of bonding types, covalent, covalent but more correlated, polar and increasingly ionic, metallogenic, electrostatic, charge-shift bonds, and dispersion interactions. Also, Q itself shows a set of interesting relations with the correlation energy of a bond.

Unifying General and Segmented Contracted Basis Sets. Segmented Polarization Consistent Basis Sets

We propose a method, denoted P-orthogonalization, for converting a general contracted basis set to a computationally more efficient segmented contracted basis set, while inheriting the full accuracy of the general contracted basis set. The procedure can be used for any general contracted basis set to remove the redundancies between general contracted functions in terms of primitive functions. The P-orthogonalization procedure is used to construct a segmented contracted version of the polarization consistent basis sets, which are optimized for density functional theory calculations. Benchmark calculations show that the new pcs-n basis sets provide uniform error control of the basis set incompleteness for molecular systems composed of atoms from the first three rows in the periodic table (H-Kr) and for different exchange-correlation functionals. The basis set errors at a given zeta quality level are lower than other existing basis sets, and the pcs-n basis sets are furthermore shown to be among the computationally most efficient. The pcs-n basis sets are available in qualities ranging from (unpolarized) double-zeta to pentuple zeta quality and should therefore be well suited for both routine and benchmark calculations using density functional theory methods in general.