A new all-round density functional based on spin states and S[sub N]2 barriers

PtIV-Containing Hexaplatinate(II) [PtIVPtII6O6(AsO2(CH3)2)6]2- and Hexapalladate(II) [PtIVPdII6O6(AsO2(CH3)2)6]2

The first PtIV-containing discrete polyoxoplatinate(II) [PtIVPtII6O6(AsO2(CH3)2)6]2- (Pt7) and polyoxopalladate(II) [PtIVPdII6O6(AsO2(CH3)2)6]2- (PtPd6) have been prepared and characterized in the solid state, in solution, and in the gas phase. The molecular structures of the noble metal-oxo clusters Pt7 and PtPd6 comprise a central, octahedral PtIVO6 hetero group surrounded by six square-planar MO4 (M = PtII, PdII) units, which are capped by six dimethylarsinate ligands. The polyanions were prepared under simple one-pot aqueous solution conditions by reacting H2Pt(OH)6 with either K2PtCl4 or Pd(NO3)2 in sodium dimethylarsinate buffer (pH 7) at 80 °C. Catalytic studies were performed on Pt7 supported on SBA15-apts for o-xylene hydrogenation at 300 °C and 90 bar H2 pressure and indicated excellent activity and recyclability with low activation temperature.

Structure and stability of the sH binary hydrate cavity and host-guest versus guest-guest interactions therein: A DFT approach

The intrinsic ability of clathrate hydrates to encage gaseous molecules is explored. Encapsulation ability depends on the cavity size and the type of guest gaseous species in addition to the physical parameters, temperature and pressure. Here we have reported the structure, stability and nature of interaction in dissimilar guest occupied sH hydrate cavity. Diatomic gas molecules and small polyatomic hydrocarbons are considered as guests. The irregular icosahedron (5¹² 6⁸ ) cavity of sH hydrate is the host. Different thermodynamic parameters for the guest molecules encapsulation were calculated using three different hybrid DFT functionals, B3LYP, M05-2X, M06, and moreover using dispersion correction (PBE0-D3). With the consideration of large H-bonded systems the 6-31G* and cc-pVTZ basis sets were used for two set of computations. To disclose the nature of interaction between the host-guest systems as well as the interaction between the guest molecules inside the host the non-covalent interaction (NCI) indices and energy decomposition analysis (EDA) were done. Impact of host-guest and guest-guest interactions are discussed.

One-electron self-interaction error and its relationship to geometry and higher orbital occupation

Density Functional Theory (DFT) sees prominent use in computational chemistry and physics; however, problems due to the self-interaction error (SIE) pose additional challenges to obtaining qualitatively correct results. As an unphysical energy an electron exerts on itself, the SIE impacts most practical DFT calculations. We conduct an in-depth analysis of the one-electron SIE in which we replicate delocalization effects for simple geometries. We present a simple visualization of such effects, which may help in future qualitative analysis of the one-electron SIE. By increasing the number of nuclei in a linear arrangement, the SIE increases dramatically. We also show how molecular shape impacts the SIE. Two- and three-dimensional shapes show an even greater SIE stemming mainly from the exchange functional with some error compensation from the one-electron error, which we previously defined [D. R. Lonsdale and L. Goerigk, Phys. Chem. Chem. Phys. 22, 15805 (2020)]. Most tested geometries are affected by the functional error, while some suffer from the density error. For the latter, we establish a potential connection with electrons being unequally delocalized by the DFT methods. We also show how the SIE increases if electrons occupy higher-lying atomic orbitals; seemingly one-electron SIE free methods in a ground are no longer SIE free in excited states, which is an important insight for some popular, non-empirical density functional approximations (DFAs). We conclude that the erratic behavior of the SIE in even the simplest geometries shows that robust DFAs are needed. Our test systems can be used as a future benchmark or contribute toward DFT development.

Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Benchmarking Approximate Methods for Adiabatic Spin-State Energy Differences in Fe(II) Complexes

The CASPT2+δMRCI composite approach reported in a companion paper has been extended and used to provide high-quality reference data for a series of adiabatic spin gaps (defined as ΔE = E_quintet - E_singlet) of [Fe^IIL₆]²⁺ complexes (L = CNH, CO, NCH, NH₃, H₂O), either at nonrelativistic level or including scalar relativistic effects. These highly accurate data have been used to evaluate the performance of various more approximate methods. Coupled-cluster theory with singles, doubles, and perturbative triples, CCSD(T), is found to agree well with the new reference data for Werner-type complexes but exhibits larger underestimates by up to 70 kJ/mol for the π-acceptor ligands, due to appreciable static correlation in the low-spin states of these systems. Widely used domain-based local CCSD(T) calculations, DLPNO-CCSD(T), are shown to depend very sensitively on the cutoff values used to construct the localized domains, and standard values are not sufficient. A large number of density functional approximations have been evaluated against the new reference data. The B2PLYP double hybrid gives the smallest deviations, but several functionals from different rungs of the usual ladder hierarchy give mean absolute deviations below 20 kJ/mol. This includes the B97-D semilocal functional, the PBE0* global hybrid with 15% exact-exchange admixture, as well as the local hybrids LH07s-SVWN and LH07t-SVWN. Several further functionals achieve mean absolute errors below 30 kJ/mol (M06L-D4, SSB-D, B97-1-D4, LC-ωPBE-D4, LH12ct-SsirPW92-D4, LH12ct-SsifPW92-D4, LH14t-calPBE-D4, LHJ-HFcal-D4, and several further double hybrids) and thereby also still overall outperform CCSD(T) or uncorrected CASPT2. While exact-exchange admixture is a crucial factor in favoring high-spin states, the present evaluations confirm that other aspects can be important as well. A number of the better-performing functionals underestimate the spin gaps for the π-acceptor ligands but overestimate them for L = NH₃, H₂O. In contrast to a previous suggestion, non-self-consistent density functional theory (DFT) computations on top of Hartree-Fock orbitals are not a promising path to produce accurate spin gaps in such complexes.

Computational Prediction of Tc-99 NMR Chemical Shifts in Technetium Complexes with Radiopharmaceutical Applications

The Tc-99m nucleus is the most used nuclide in radiopharmaceuticals designed for imaging diagnosis. The metal can exist in nine distinct oxidation states and forms distinct coordination complexes with a variety of chelating agents and geometries. These complexes are usually characterized through Tc-99 NMR that is very sensitive to the Tc coordination sphere. Therefore, predicting Tc-99 NMR might be useful to assist experimentalists in structural characterization. In the present study, we propose three computational protocols for predicting Tc-99 NMR chemical shifts based on density functional theory calculations using relativistic and nonrelativistic Hamiltonians: the relativistic Model 1, the nonrelativistic Model 2, and the empirical nonrelativistic Model 3. In Models 2 and 3, the NMR-DKH basis set was used for all atoms, including the Tc, for which it was developed here. All models were applied for a set of 41 Tc-complexes with metal oxidation states 0, I, and V, for which the Tc-99 chemical shift was available experimentally. The mean absolute deviation and the mean relative deviation were 67 ppm and 4.8% (Model 1), 92 ppm and 6.2% (Model 2), and 65 ppm and 4.9% (Model 3), respectively. Last, the effect of the explicit solvent was evaluated for the [TcO₂(en)₂]⁺─Tc(V) complex. The calculated results for the Tc-99 NMR chemical shift at SO-ZORA-SSB-D/TZ2P-ZORA/COSMO//TPSS/def2-SVP/IEF-PCM(UFF) show that the inclusion of 14 water molecules (first solvation shell) together with the implicit solvation model leads to an absolute deviation of only 7 ppm (0.3%) from the experimental value, indicating that the solvent effects play a key role in predicting Tc-99 NMR.

Discovery of Polythioplatinate(II) [Pt3S2(SO3)6]10- and Study of Its Solution and Catalytic Properties

We have discovered the first polythioplatinate(II), [Pt^II₃S₂(SO₃)₆]^10- (1), which was synthesized in aqueous basic medium (pH 11) by hydrothermal heating at 150 °C. Polyanion 1 comprises a discrete, triangular assembly of three Pt²⁺ ions linked by two μ₃-sulfido ligands, and their square-planar coordination geometry is completed by two terminal S-bound sulfito ligands. Polyanion 1 was isolated as a hydrated sodium salt, Na₁₀[Pt^II₃(μ₃-S)₂(SO₃)₆]·22H₂O (Na-1), which was characterized in the solid state by single-crystal X-ray diffraction, Fourier-transform infrared spectroscopy, thermogravimetric analysis, X-ray photoelectron spectra, and elemental analysis, in solution by ¹⁹⁵Pt NMR and atomic absorption spectroscopy, and in the gas phase by electrospray ionization mass spectrometry. Density functional theory calculations were performed, and the ¹⁹⁵Pt NMR chemical shifts of 1 were computed theoretically and shown to match well with the experimental data. Furthermore, the discrete title polyanion 1 was immobilized on mesoporous SBA-15 support and used as a precatalyst for the hydrogenation of o-xylene.

Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics

Spin-state energetics of transition metal complexes remain one of the most challenging targets for electronic structure methods. Among single-reference wave function approaches, local correlation approximations to coupled cluster theory, most notably the domain-based local pair natural orbital (DLPNO) approach, hold the promise of bringing the accuracy of coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T), to molecular systems of realistic size with acceptable computational cost. However, recent studies on spin-state energetics of iron-containing systems raised doubts about the ability of the DLPNO approach to adequately and systematically approximate energetics obtained by the reference-quality complete active space second-order perturbation theory with coupled-cluster semicore correlation, CASPT2/CC. Here, we revisit this problem using a diverse set of iron complexes and examine several aspects of the application of the DLPNO approach. We show that DLPNO-CCSD(T) can accurately reproduce both CASPT2/CC and canonical CCSD(T) results if two basic principles are followed. These include the consistent use of the improved iterative (T₁) versus the semicanonical perturbative triple corrections and, most importantly, a simple two-point extrapolation to the PNO space limit. The latter practically eliminates errors arising from the default truncation of electron-pair correlation spaces and should be viewed as standard practice in applications of the method to transition metal spin-state energetics. Our results show that reference-quality results can be readily achieved with DLPNO-CCSD(T) if these principles are followed. This is important also in view of the applicability of the method to larger single-reference systems and multinuclear clusters, whose treatment of dynamic correlation would be challenging for multireference-based approaches.

Chalcogen bonds: Hierarchical ab initio benchmark and density functional theory performance study

We have performed a hierarchical ab initio benchmark and DFT performance study of D2 Ch•••A- chalcogen bonds (Ch = S, Se; D, A = F, Cl). The ab initio benchmark study is based on a series of ZORA-relativistic quantum chemical methods [HF, MP2, CCSD, CCSD(T)], and all-electron relativistically contracted variants of Karlsruhe basis sets (ZORA-def2-SVP, ZORA-def2-TZVPP, ZORA-def2-QZVPP) with and without diffuse functions. The highest-level ZORA-CCSD(T)/ma-ZORA-def2-QZVPP counterpoise-corrected complexation energies (ΔECPC ) are converged within 1.1-3.4 kcal mol-1 and 1.5-3.1 kcal mol-1 with respect to the method and basis set, respectively. Next, we used the ZORA-CCSD(T)/ma-ZORA-def2-QZVPP (ΔECPC ) as reference data for analyzing the performance of 13 different ZORA-relativistic DFT approaches in combination with the Slater-type QZ4P basis set. We find that the three-best performing functionals are M06-2X, B3LYP, and M06, with mean absolute errors (MAE) of 4.1, 4.2, and 4.3 kcal mol-1 , respectively. The MAE for BLYP-D3(BJ) and PBE amount to 8.5 and 9.3 kcal mol-1 , respectively.

Targeting structural features of viral genomes with a nano-sized supramolecular drug

RNA targeting is an exciting frontier for drug design. Intriguing targets include functional RNA structures in structurally-conserved untranslated regions (UTRs) of many lethal viruses. However, computational docking screens, valuable in protein structure targeting, fail for inherently flexible RNA. Herein we harness MD simulations with Markov state modeling to enable nanosize metallo-supramolecular cylinders to explore the dynamic RNA conformational landscape of HIV-1 TAR untranslated region RNA (representative for many viruses) replicating experimental observations. These cylinders are exciting as they have unprecedented nucleic acid binding and are the first supramolecular helicates shown to have anti-viral activity in cellulo: the approach developed in this study provides additional new insight about how such viral UTR structures might be targeted with the cylinder binding into the heart of an RNA-bulge cavity, how that reduces the conformational flexibility of the RNA and molecular details of the insertion mechanism. The approach and understanding developed represents a new roadmap for design of supramolecular drugs to target RNA structural motifs across biology and nucleic acid nanoscience.

Spin-state dependence of exchange-correlation holes

Applications of density-functional theory (DFT) in computational chemistry rely on an approximate exchange-correlation (xc) functional. However, existing approximations can fail dramatically for open-shell molecules, in particular for transition-metal complexes or radicals. Most importantly, predicting energy differences between different spin-states with approximate exchange-correlation functionals remains extremely challenging. Formally, it is known that the exact xc functional should be spin-state dependent, but none of the available approximations feature such an explicit spin-state dependence [C. R. Jacob and M. Reiher, Int. J. Quantum Chem., 2012, 112, 3661-3684]. Thus, to find novel approximations for the xc functional for open-shell systems, the development of spin-state dependent xc functionals appears to be a promising avenue. Here, we set out to shed light on the spin-state dependence of the xc functional by investigating the underlying xc holes, which we extract from configuration interaction calculations for model systems. We analyze the similarities and differences between the xc holes of the lowest-energy singlet and triplet states of the dihydrogen molecule, the helium atom, and the lithium dimer. To shed further light on the spin-state dependence of these xc holes we also discuss exact conditions that can be derived from the spin structure of the reduced two-electron density matrix. Altogether, our results suggest several possible routes towards the construction of explicitly spin-state dependent approximations for the xc functional.

Metal-Organic Frameworks of Cu(II) Constructed from Functionalized Ligands for High Capacity H2 and CO2 Gas Adsorption and Catalytic Studies

Two Cu(II)-based metal-organic frameworks (MOFs) having paddle-wheel secondary building units (SBUs), namely, 1_Me and 1_ipr, were synthesized solvothermally using two new bent di-isophthalate ligands incorporating different substituents. The MOFs showed high porosity (BET surface area, 2191 m²/g for 1_Me and 1402 m²/g for 1_ipr). For 1_Me, very high CO₂ adsorption (98.5 wt % at 195 K, 42.9 wt % at 273 K, 23.3 wt % at 298 K) at 1 bar was found, while for 1_ipr, it was significantly less (14.3 wt % at 298 K and 1 bar, 54.4 wt % at 298 K at 50 bar). 1_Me exhibited H₂ uptake of 3.2 wt % at 77 K and 1 bar of pressure, which compares well with other benchmark MOFs. For 1_ipr, the H₂ uptake was found to be 2.54 wt % under similar experimental conditions. The significant adsorption of H₂ and CO₂ for 1_Me could be due to the presence of micropores as well as unsaturated metal sites in these MOFs besides the presence of substituents that interact with the gas molecules. The experimental adsorption behavior of the MOFs could be justified by theoretical calculations. Additionally, catalytic conversions of CO₂ and CS₂ into useful chemicals like cyclic carbonates, cyclic trithiocarbonates, and cyclic dithiocarbonates could be achieved.

Iron(II) and Iron(III) Spin Crossover: Toward an Optimal Density Functional

Spin crossover (SCO) plays a major role in biochemistry, catalysis, materials, and emerging technologies such as molecular electronics and sensors, and thus accurate prediction and design of SCO systems is of high priority. However, the main tool for this purpose, density functional theory (DFT), is very sensitive to applied methodology. The most abundant SCO systems are Fe(II) and Fe(III) systems. Even with average good agreement, a functional may be significantly more accurate for Fe(II) or Fe(III) systems, preventing balanced study of SCO candidates of both types. The present work investigates DFT's performance for well-known Fe(II) and Fe(III) SCO complexes, using various design types and customized versions of GGA, hybrid, meta-GGA, meta-hybrid, double-hybrid, and long-range-corrected hybrid functionals. We explore the limits of DFT performance and identify proficient Fe(II)-Fe(III)-balanced functionals. We identify and quantify remarkable differences in the DFT description of Fe(II) and Fe(III) systems. Most functionals become more accurate once Hartree-Fock exchange is adjusted to 10-17%, regardless of the type of functionals involved. However, this typically introduces a clear Fe(II)-Fe(III) bias. The most accurate functionals measured by mean absolute errors <10 kJ/mol are CAMB3LYP-17, B3LYP*, and B97-15 with 15-17% Hartree-Fock exchange, closely followed by CAMB3LYP and CAMB3LYP-15, OPBE, rPBE-10, and B3P86-15. While GGA functionals display a small Fe(II)-Fe(III) bias, they are generally inaccurate, except the O exchange functional. Hybrid functionals (including B2PLYP double hybrids and meta hybrids) tend to favor HS too much in Fe(II) vs Fe(III), which is important in many studies where the oxidation state of iron can vary, e.g. rational SCO design and studies of catalytic processes involving iron. The only functional with a combined bias <5 kJ/mol and a decent MAE (15 kJ/mol) is our customized PBE0-12 functional. Alternatively one has to sacrifice Fe(II)-Fe(III) balance to use the best functionals for each group separately. We also investigated the precision (measured as the standard deviation of errors) and show that the target accuracy for iron SCO is 10 kJ/mol for accuracy and 5 kJ/mol for precision, and DFT is probably not going to break this limit in the near future. Importantly, all four types of functional behavior (accurate/precise, accurate/imprecise, inaccurate/precise, inaccurate/imprecise) are observed. More generally, our work illustrates the importance not only of overall accuracy but also of balanced accuracy for systems likely to occur in context.

The role of spin states in the catalytic mechanism of the intra- and extradiol cleavage of catechols by O2

Iron-dependent enzymes and biomimetic iron complexes can catalyze the ring cleavage of very inert, aromatic compounds. The mechanisms of these transformations and the factors that lead either to extradiol cleavage or intradiol cleavage have not been fully understood. By using density functional theory we have elucidated the mechanism of the catalytic cycle for two biomimetic complexes, and explained the difference in the experimentally obtained products.

Evaluation of the Factors Impacting the Accuracy of 13C NMR Chemical Shift Predictions using Density Functional Theory-The Advantage of Long-Range Corrected Functionals

The various factors influencing the accuracy of ¹³C NMR calculations using density functional theory (DFT), including the basis set, exchange-correlation (XC) functional, and isotropic shielding calculation method, are evaluated. A wide selection of XC functionals (over 70) were considered, and it was found that long-range corrected functionals offer a significant improvement over the other classes of functionals. Based on a thorough study, it is recommended that for calculating NMR chemical shifts (δ) one should use the CSGT method, the COSMO solvation model, and the LC-TPSSTPSS exchange-correlation functional in conjunction with the cc-pVTZ basis set. A selection of problems in natural product identification are considered in light of the newly recommended level of theory.

A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions

We present the updated and extended GMTKN55 benchmark database for more accurate and extensive energetic evaluation of density functionals and other electronic structure methods with detailed guidelines for method users.

Spinning around in Transition-Metal Chemistry

The great diversity and richness of transition metal chemistry, such as the features of an open d-shell, opened a way to numerous areas of scientific research and technological applications. Depending on the nature of the metal and its environment, there are often several energetically accessible spin states, and the progress in accurate theoretical treatment of this complicated phenomenon is presented in this Account. The spin state energetics of a transition metal complex can be predicted theoretically on the basis of density functional theory (DFT) or wave function based methodology, where DFT has advantages since it can be applied routinely to medium-to-large-sized molecules and spin-state consistent density functionals are now available. Additional factors such as the effect of the basis set, thermochemical contributions, solvation, relativity, and dispersion, have been investigated by many researchers, but challenges in unambiguous assignment of spin states still remain. The first DFT studies showed intrinsic spin-state preferences of hybrid functionals for high spin and early generalized gradient approximation functionals for low spin. Progress in the development of density functional approximations (DFAs) then led to a class of specially designed DFAs, such as OPBE, SSB-D, and S12g, and brought a very intriguing and fascinating observation that the spin states of transition metals and the S_N2 barriers of organic molecules are somehow intimately linked. Among the many noteworthy results that emerged from the search for the appropriate description of the complicated spin state preferences in transition metals, we mainly focused on the examination of the connection between the spin state and the structures or coordination modes of the transition metal complexes. Changes in spin states normally lead only to changes in the metal-ligand bond lengths, but to the best of our knowledge, the dapsox ligand showed the first example of a transition-metal complex where a change in spin state leads also to changes in the coordination, switching between pentagonal-bipyramidal and capped-octahedron. Moreover, we have summarized the results of the thorough study that corrected the experimental assignment of the nature of the recently synthesized Sc³⁺ adduct of [Fe^IV(O)(TMC)]²⁺ (TMC = 1,4,8,11-tetramethylcyclam) and firmly established that the Sc³⁺-capped iron-oxygen complex corresponds to high-spin Fe^III. Last, but not least, we have provided deeper insight and rationalization of the observation that unlike in metalloenzymes, where the Fe^IV-oxo is usually observed with high spin, biomimetic Fe^IV-oxo complexes typically have a intermediate spin state. Energy decomposition analyses on the trigonal-bypiramidal (TBP) and octahedral model systems with ammonia ligands have revealed that the interaction energy of the prepared metal ion in the intermediate spin state is much smaller for the TBP structure. This sheds light on the origin of the intermediate spin state of the biomimetic TBP Fe^IV-oxo complexes.

Reactivity of an FeIV-Oxo Complex with Protons and Oxidants

High-valent Fe-OH species are often invoked as key intermediates but have only been observed in Compound II of cytochrome P450s. To further address the properties of non-heme FeIV-OH complexes, we demonstrate the reversible protonation of a synthetic FeIV-oxo species containing a tris-urea tripodal ligand. The same protonated FeIV-oxo species can be prepared via oxidation, suggesting that a putative FeV-oxo species was initially generated. Computational, Mössbauer, XAS, and NRVS studies indicate that protonation of the FeIV-oxo complex most likely occurs on the tripodal ligand, which undergoes a structural change that results in the formation of a new intramolecular H-bond with the oxido ligand that aids in stabilizing the protonated adduct. We suggest that similar protonated high-valent Fe-oxo species may occur in the active sites of proteins. This finding further argues for caution when assigning unverified high-valent Fe-OH species to mechanisms.

Electron Transport in a Dioxygenase-Ferredoxin Complex: Long Range Charge Coupling between the Rieske and Non-Heme Iron Center

Dioxygenase (dOx) utilizes stereospecific oxidation on aromatic molecules; consequently, dOx has potential applications in bioremediation and stereospecific oxidation synthesis. The reactive components of dOx comprise a Rieske structure Cys₂[2Fe-2S]His₂ and a non-heme reactive oxygen center (ROC). Between the Rieske structure and the ROC, a universally conserved Asp residue appears to bridge the two structures forming a Rieske-Asp-ROC triad, where the Asp is known to be essential for electron transfer processes. The Rieske and ROC share hydrogen bonds with Asp through their His ligands; suggesting an ideal network for electron transfer via the carboxyl side chain of Asp. Associated with the dOx is an itinerant charge carrying protein Ferredoxin (Fdx). Depending on the specific cognate, Fdx may also possess either the Rieske structure or a related structure known as 4-Cys-[2Fe-2S] (4-Cys). In this study, we extensively explore, at different levels of theory, the behavior of the individual components (Rieske and ROC) and their interaction together via the Asp using a variety of density function methods, basis sets, and a method known as Generalized Ionic Fragment Approach (GIFA) that permits setting up spin configurations manually. We also report results on the 4-Cys structure for comparison. The individual optimized structures are compared with observed spectroscopic data from the Rieske, 4-Cys and ROC structures (where information is available). The separate pieces are then combined together into a large Rieske-Asp-ROC (donor/bridge/acceptor) complex to estimate the overall coupling between individual components, based on changes to the partial charges. The results suggest that the partial charges are significantly altered when Asp bridges the Rieske and the ROC; hence, long range coupling through hydrogen bonding effects via the intercalated Asp bridge can drastically affect the partial charge distributions compared to the individual isolated structures. The results are consistent with a proton coupled electron transfer mechanism.

Complexes of Ni(ii) and Cu(ii) with small peptides: deciding whether to deprotonate

Infrared multiple photon dissociation (IRMPD) spectroscopy differentiates two binding modes (iminol versus charge solvated) for Ni(ii) bound to model peptides.

Density Functional Calculations of E2 and SN2 Reactions: Effects of the Choice of Method, Algorithm, and Numerical Accuracy

Herein we provide a detailed account on how the potential energy surfaces of the E2 and SN2 reactions of X(-) + CH3CH2X (X = F, Cl) depend on various methodological and technical choices in density functional calculations. We cover a choice of density functionals (OLYP, various M06-types, and the new SSB-D), basis sets (up to quintuple- and quadruple-ζ for Gaussian- and Slater-type orbitals, respectively, plus polarization and diffuse functions), and other aspects of the computations (among others: nonrelativistic versus zeroth-order regular approximation relativistic; numerical integration accuracy; all-electron versus frozen core; self-consistent field (SCF) versus post-SCF). The program codes ADF and NWChem are used for calculations with Slater- and Gaussian-type basis sets, respectively. The fluoride systems (X = F) appear to not only depend extremely sensitively on the basis set size (especially the presence of diffuse functions) but also on other technical settings, especially in the case of hybrid meta-generalized gradient approximation functionals. This work complements a recent contribution (Y. Zhao, D. G. Truhlar, J. Chem. Theory Comput. 2010, 6, 1104) and provides recommendations for density functionals, basis sets, and technical settings.

A faux hawk fullerene with PCBM-like properties

Reaction of C60, C6F5CF2I, and SnH(n-Bu)3 produced, among other unidentified fullerene derivatives, the two new compounds 1,9-C60(CF2C6F5)H (1) and 1,9-C60(cyclo-CF2(2-C6F4)) (2). The highest isolated yield of 1 was 35% based on C60. Depending on the reaction conditions, the relative amounts of 1 and 2 generated in situ were as high as 85% and 71%, respectively, based on HPLC peak integration and summing over all fullerene species present other than unreacted C60. Compound 1 is thermally stable in 1,2-dichlorobenzene (oDCB) at 160 °C but was rapidly converted to 2 upon addition of Sn2(n-Bu)6 at this temperature. In contrast, complete conversion of 1 to 2 occurred within minutes, or hours, at 25 °C in 90/10 (v/v) PhCN/C6D6 by addition of stoichiometric, or sub-stoichiometric, amounts of proton sponge (PS) or cobaltocene (CoCp2). DFT calculations indicate that when 1 is deprotonated, the anion C60(CF2C6F5)- can undergo facile intramolecular SNAr annulation to form 2 with concomitant loss of F-. To our knowledge this is the first observation of a fullerene-cage carbanion acting as an SNAr nucleophile towards an aromatic C-F bond. The gas-phase electron affinity (EA) of 2 was determined to be 2.805(10) eV by low-temperature PES, higher by 0.12(1) eV than the EA of C60 and higher by 0.18(1) eV than the EA of phenyl-C61-butyric acid methyl ester (PCBM). In contrast, the relative E1/2(0/-) values of 2 and C60, -0.01(1) and 0.00(1) V, respectively, are virtually the same (on this scale, and under the same conditions, the E1/2(0/-) of PCBM is -0.09 V). Time-resolved microwave conductivity charge-carrier yield × mobility values for organic photovoltaic active-layer-type blends of 2 and poly-3-hexylthiophene (P3HT) were comparable to those for equimolar blends of PCBM and P3HT. The structure of solvent-free crystals of 2 was determined by single-crystal X-ray diffraction. The number of nearest-neighbor fullerene-fullerene interactions with centroid···centroid (⊙···⊙) distances of ≤10.34 Å is significantly greater, and the average ⊙···⊙ distance is shorter, for 2 (10 nearest neighbors; ave. ⊙···⊙ distance = 10.09 Å) than for solvent-free crystals of PCBM (7 nearest neighbors; ave. ⊙···⊙ distance = 10.17 Å). Finally, the thermal stability of 2 was found to be far greater than that of PCBM.

Accurate calculation of31P NMR chemical shifts in polyoxometalates

We search for an optimal DFT/basis set combination for computationally reproducing³¹P NMR chemical shifts in polyoxometalates.

Pair-potential approach to accurate dispersion energies between group 12 (Zn, Cd, Hg) clusters

Dispersion interactions between group 12 (Zn, Cd, Hg) metal clusters are studied at the CCSD(T) level with triple-ζ basis sets. We use atomic orientation dependent C6 dispersion coefficients derived from simple model systems to calculate the intermolecular dispersion energy between larger metal clusters. By using an atomic pair-potential model, we are able to obtain highly accurate results comparable to the CCSD(T) level of theory. For all the systems studied, our method results in a robust description of the dispersion energy at long distances, and it is able to outperform the examined dispersion corrected DFT calculations.

Structural and electronic properties of the Pt(n)-PAH complex (n = 1, 2) from density functional calculations

A detailed density functional study of the Pt atom and the Pt dimer adsorption on a polyaromatic hydrocarbon (PAH) is presented. The preferred adsorption site for a Pt atom is confirmed to be the bridge site. Upon adsorption of a single Pt atom, however, it is found here that the electronic ground state changes from the triplet state (5d(9)6s(1) configuration) to the closed-shell singlet state (5d(10)6s(0) configuration), which consequently will affect the catalytic activity of Pt when single Pt atoms bind to a carbon surface. The preferred adsorption site for the Pt dimer in the upright configuration is the hollow site. In contrast to the adsorption of a single Pt atom, the formation of a Pt-C bond in the adsorption of a Pt dimer is not accompanied by a change in the spin state, so the most stable electronic state is still the triplet state. While the atomic charge on the Pt atoms and dimers (in parallel configuration) in the Ptn-PAH complex is positive, a negative charge is found on the upper Pt atom for the upright configuration, indicating that single layers of Pt atoms will have a different catalytic activity as compared to Pt clusters on a carbon surface. Comparing the Pt-C bond length and the charge transfer on different sites, the magnitude of the charge transfer decreases with bond elongation, indicating that the catalytic activity of the Pt atom and dimer can be changed by modifying its chemical surroundings. The adsorption energy for the Pt dimer on a PAH surface is larger than that for two individual Pt atoms on the surface indicating that aggregation of Pt atoms on the PAH surface is favorable.

Benchmark study on the smallest bimolecular nucleophilic substitution reaction: H⁻+CH₄-->CH₄+H⁻

We report here a benchmark study on the bimolecular nucleophilic substitution (S(N)2) reaction between hydride and methane, for which we have obtained reference energies at the coupled cluster toward full configuration-interaction limit (CC-cf/CBS). Several wavefunction (HF, MP2, coupled cluster) and density functional methods are compared for their reliability regarding these reference data.