W4‐17: A diverse and high‐confidence dataset of atomization energies for benchmarking high‐level electronic structure methods

Atomization reactions are among the most challenging tests for electronic structure methods. We use the first‐principles Weizmann‐4 (W4) computational thermochemistry protocol to generate the W4‐17 dataset of 200 total atomization energies (TAEs) with 3σ confidence intervals of 1 kJ mol⁻¹. W4‐17 is an extension of the earlier W4‐11 dataset; it includes first‐ and second‐row molecules and radicals with up to eight non‐hydrogen atoms. These cover a broad spectrum of bonding situations and multireference character, and as such are an excellent benchmark for the parameterization and validation of highly accurate ab initio methods (e.g., CCSD(T) composite procedures) and double‐hybrid density functional theory (DHDFT) methods. The W4‐17 dataset contains two subsets (i) a non‐multireference subset of 183 systems characterized by dynamical or moderate nondynamical correlation effects (denoted W4‐17‐nonMR) and (ii) a highly multireference subset of 17 systems (W4‐17‐MR). We use these databases to evaluate the performance of a wide range of CCSD(T) composite procedures (e.g., G4, G4(MP2), G4(MP2)‐6X, ROG4(MP2)‐6X, CBS‐QB3, ROCBS‐QB3, CBS‐APNO, ccCA‐PS3, W1, W2, W1‐F12, W2‐F12, W1X‐1, and W2X) and DHDFT methods (e.g., B2‐PLYP, B2GP‐PLYP, B2K‐PLYP, DSD‐BLYP, DSD‐PBEP86, PWPB95, ωB97X‐2(LP), and ωB97X‐2(TQZ)). © 2017 Wiley Periodicals, Inc.     

Theoretical investigation of the structures and spectroscopic properties of (H2O4)n (n = 1–4) clusters

Density functional theory and ab initio calculations were performed to elucidate the hydrogen interactions in (H₂O₄)_n (n = 1–4) clusters. The optimized geometries, binding energies, and harmonic vibrational frequencies were predicted at various levels of theory. The trans conformer of the H₂O₄ monomer was predicted to be the most stable structure at the CCSD(T)/aug‐cc‐pVTZ level of theory. The binding energies per H₂O₄ monomer increased in absolute value by 9.0, 10.1, and 11.8 kcal/mol from n = 2 to n = 4 at the MP2/cc‐pVTZ level of theory (after the zero‐point vibrational energy and basis set superposition error corrections). This result implies that the intermolecular hydrogen bonds were stronger in the long‐chain clusters, that is, the formation of the longer chain in the (H₂O₄)_n clusters was more energetically favorable.     

Stability of the chlorinated derivatives of the DNA/RNA nucleobases,purine and pyrimidine toward radical formation via homolytic CCl bond dissociation

The chlorinated derivatives of nucleobases (and nucleosides), as well as those of purine, have well‐established anticancer activity, and in some cases, are also shown to be involved in the link between chronic inflammatory conditions and the development of cancer. In this investigation, the stability of all of the isomeric forms of the chlorinated nucleobases, purine and pyrimidine are investigated from the perspective of their homolytic C?Cl bond dissociation energies (BDEs). The products of these reactions, namely chlorine atom and the corresponding carbon‐centered radicals, may be of importance in terms of potentiating biological damage. Initially, the performance of a wide range of contemporary theoretical procedures were evaluated for their ability to afford accurate C?Cl BDEs, using a recently reported set of 28 highly accurate C?Cl BDEs obtained by means of W1w theory. Subsequent to this analysis, the G3X(MP2)‐RAD procedure (which achieves a mean absolute deviation of merely 1.3 kJ mol^?1, with a maximum deviation of 5.0 kJ mol^?1) was employed to obtain accurate gas‐phase homolytic C?Cl bond dissociation energies for a wide range of chlorinated isomers of the DNA/RNA nucleobases, purine and pyrimidine.     

Assessment of ab initio MP2 and density functionals for characterizing the potential energy profiles of the SN2 reactions at N center

The potential energy profiles of five selected bimolecular nucleophilic substitution (S_N2) reactions at nitrogen (N) center have been reinvestigated with the CCSD(T), G3[MP2,CCSD(T)], MP2, and some density functional methods. The basis sets of 6‐31+G(d,p) and 6‐311+G(3d,2p) are used for the MP2 and density functional calculations. Taking the relative energies at the CCSD(T)/CBS level of theory as benchmarks, we recommend the MP2, B97‐K, B2K‐PLYP, BMK, ωB97X‐D, M06‐2X, M05‐2X, CAM‐B3LYP, M08‐SO, and ωB97X methods to generally characterize the potential energy profiles for the S_N2 reactions at N center. Furthermore, these recommended methods with the relatively small 6‐31+G(d,p) basis set may also be used to perform direct classical trajectory simulations to uncover the dynamic behaviors of the S_N2 reactions at N center. © 2012 Wiley Periodicals, Inc.     

Gibbs energy of activation for thermal isomerization of (1Z)‐acetaldehyde hydrazone and (1Z)‐acetaldehyde N,N‐dimethylhydrazone by Gaussian‐4 theory and CCSD(T)/CBS computations

In this article, we examined the Gibbs energy of activation for the Z/E thermal isomerization reaction of (1Z)‐acetaldehyde hydrazone and (1Z)‐acetaldehyde N,N‐dimethylhydrazone, at 298.15 K in the solvent of cyclohexane. We carried out computations employing both the Gaussian‐4 (G4) theory and the coupled cluster method using both single and double substitutions and triple excitations noniteratively, CCSD(T). The CCSD(T) energy is extrapolated to the complete basis set (CBS). We compared the calculated results to the available experimental observation. It appeared that both G4 and CCSD(T)/CBS computations overestimated the experimental value by as much as about 6 and 12 kcal/mol in the present two cases. We discussed possible sources of error and proposed the experimental kinetic data could be questionable. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Analysis of ToF‐SIMS spectra of poly(2‐vinylpyridine) and poly(4‐vinylpyridine) with density functional theory calculations

The time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS) positive and negative ion spectra of poly(2‐vinylpyridine) (P2VP) and poly(4‐vinylpyridine) (P4VP) were analyzed using density functional theory calculations. Most of the ions from these structural isomers shared the same accurate mass, but had different relative abundance. This could be attributed to the fact that from a thermodynamics perspective, the disparity in the molecular structures can affect the ion stability if we assume that they shared the same mechanistic pathway of formation with similar reaction kinetics. The molecular structures of these ions were assigned, and their stability was evaluated based on calculations using the Kohn‐Sham density functional theory with Becke's 3‐parameter Lee‐Yang‐Parr exchange‐correlation functional and a correlation‐consistent, polarized, valence, double‐zeta basis set for cations and the same basis set with a triple‐zeta for anions. The computational results agreed with the experimental observations that the nitrogen‐containing cations such as C₅H₄N⁺ (m/z = 78), C₈H₇N^+· (m/z = 117), C₈H₈N⁺ (m/z = 118), C₉H₈N⁺ (m/z = 130), C₁₃H₁₁N₂⁺ (m/z = 195), C₁₄H₁₃N₂⁺ (m/z = 209), C₁₅H₁₅N₂⁺ (m/z = 223), and C₂₁H₂₂N₃⁺ (m/z = 316) ions were more favorably formed in P2VP than in P4VP due to higher ion stability because the calculated total energies of these cations were more negative when the nitrogen was situated at the ortho position. Nevertheless, our assumption was invalid in the formation of positive ions such as C₆H₇N⁺˙ (m/z = 93) and C₈H₁₀N⁺ (m/z = 120). Their formation did not necessarily depend on the ion stability. Instead, the transition state chemistry and the matrix effect both played a role. In the negative ion spectra, we found that nitrogen‐containing anions such as C₅H₄N^? (m/z = 78), C₆H₆N^? (m/z = 92), C₇H₆N^? (m/z = 104), C₈H₆N^? (m/z = 116), C₉H₁₀N^? (m/z = 132), C₁₃H₁₁N₂^? (m/z = 195), and C₁₄H₁₃N₂^? (m/z = 209) ions were more favorably formed in P4VP, which is in line with our computational results without exception. We speculate that whether anions would form from P2VP and P4VP is more dependent on the stability of the ions.     

Theoretical investigation on identical anionic halide‐exchange SN2 reaction processes on N‐haloammonium cation NH3X+ (X = F,Cl, Br,and I)

CCSD(T) calculations have been used for identically nucleophilic substitution reactions on N‐haloammonium cation, X^? + NH₃X⁺ (X = F, Cl, Br, and I), with comparison of classic anionic S_N2 reactions, X^? + CH₃X. The described S_N2 reactions are characterized to a double curve potential, and separated charged reactants proceed to form transition state through a stronger complexation and a charge neutralization process. For title reactions X^? + NH₃X⁺, charge distributions, geometries, energy barriers, and their correlations have been investigated. Central barriers ΔE for X^? + NH₃X⁺ are found to be lower and lie within a relatively narrow range, decreasing in the following order: Cl (21.1 kJ/mol) > F (19.7 kJ/mol) > Br (10.9 kJ/mol) > I (9.1 kJ/mol). The overall barriers ΔE relative to the reactants are negative for all halogens: ?626.0 kJ/mol (F), ?494.1 kJ/mol (Cl), ?484.9 kJ/mol (Br), and ?458.5 kJ/mol (I). Stability energies of the ion–ion complexes ΔE_comp decrease in the order F (645.6 kJ/mol) > Cl (515.2 kJ/mol) > Br (495.8 kJ/mol) > I (467.6 kJ/mol), and are found to correlate well with halogen Mulliken electronegativities (R² = 0.972) and proton affinity of halogen anions X^? (R² = 0.996). Based on polarizable continuum model, solvent effects have investigated, which indicates solvents, especially polar and protic solvents lower the complexation energy dramatically, due to dually solvated reactant ions, and even character of double well potential in reactions X^? + CH₃X has disappeared. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Reaction cycles and poisoning in catalysis by gold clusters: A thermodynamics approach

In heterogeneous catalysis, a catalytic process takes place at finite temperature and at finite pressure of the atmosphere of the reactant gases. By applying ab initio atomistic thermodynamics to the model case of free Au₂ and clusters in an atmosphere of O₂ and CO, we derive all the thermodynamically possible reaction paths for the oxidation of CO to CO₂. This analysis lets us explain how gold clusters enable oxidation reactions without breaking the spin‐conservation rule. Furthermore, we identify special cluster + ligands compositions such as reaction intermediates and poisoned species. In particular, a thermodynamically driven poisoning is identified for the catalytic system containing free Au₂, and the experimental (p, T) conditions that avoid its formation are suggested. This implies that for some systems a catalytic cycle can be established, on thermodynamics grounds, only in a defined range of temperatures and pressures. In addition, our predictions for provide the so far most complete interpretation of the available experimental data (Socaciu et al, J. Am. Chem. Soc. 2003). © 2013 Wiley Periodicals, Inc.     

Computational studies of bond dissociation energies, ionization potentials, and heat of formation for NH and NH+. Are hybrid density functional theory methods as accurate as quadratic complete basis set and Gaussian-2 ab initio methods?

Ionization potentials, bond dissociation energies, and heat of formation for NH and NH⁺ molecular species as well as for their elements were computed with highly reliable quadratic complete basis set and Gaussian-2 ab initio methods. The results are compared with experimental results and the assurance of these ab initio approaches is assessed. The same studies were also performed with three hybrid density functional methods (B3LYP, B3P86, and B3PW91) in combination with variously sized basis sets. The computational results are discussed in light of density functional theory reliability for exploring the potential energy of small polar molecular systems. Received: 21 July 1997 / Accepted: 8 December 1997     

Theoretical investigation of LaC3n+ (n = 0, 1, 2) clusters by density functional theory

LaC₃ⁿ⁺ (n=0, 1, 2) clusters have been studied using B3LYP (Becke 3-parameter–Lee-Yang-Parr) density functional method. The basis set is Dunning/Huzinaga valence double zeta for carbon and [2s2p2d] for lanthanum, denoted LANL1DZ. Four isomers are presented for each cluster; two of them are edge binding isomers with C_2v symmetry, the other two are linear chains with C_∞v symmetry. Meanwhile, two spin states for each isomer, that is, singlet and triplet for LaC₃⁺, doublet and quartet for LaC₃ and LaC₃²⁺, respectively, are also considered. Geometries, vibrational frequencies, infrared intensities, and other quantities are reported and discussed. The results indicate that at some spin states; the C_2v symmetry isomers are the dominant structures, while for the other spin states, linear isomers are energetically favored. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 301–307, 1998     

Assessing the role of Hartree‐Fock exchange,correlation energy and long range corrections in evaluating ionization potential,and electron affinity in density functional theory

Accurate determination of ionization potentials (IPs), electron affinities (EAs), fundamental gaps (FGs), and HOMO, LUMO energy levels of organic molecules play an important role in modeling and predicting the efficiencies of organic photovoltaics, OLEDs etc. In this work, we investigate the effects of Hartree Fock (HF) Exchange, correlation energy, and long range corrections in predicting IP and EA in Hybrid Functionals. We observe increase in percentage of HF exchange results in increase of IPs and decrease in EAs. Contrary to the general expectations inclusion of both HF exchange and correlation energy (from the second order perturbation theory MP2) leads to poor prediction. Range separated Hybrid Functionals are found to be more reliable among various DFT Functionals investigated. DFT Functionals predict accurate IPs whereas post HF methods predict accurate EAs. © 2017 Wiley Periodicals, Inc.     

Definitive heat of formation of methylenimine,CH2NH,and of methylenimmonium ion,CH2NH2+, by means of W2 theory

A long‐standing controversy concerning the heat of formation of methylenimine has been addressed by means of the W2 (Weizmann‐2) thermochemical approach. Our best calculated values, ΔH°_f,298(CH₂NH) = 21.1±0.5 kcal/mol and ΔH°_f,298(CH₂NH₂⁺) = 179.4±0.5 kcal/mol, are in good agreement with the most recent measurements but carry a much smaller uncertainty. As a byproduct, we obtain the first‐ever accurate anharmonic force field for methylenimine: upon consideration of the appropriate resonances, the experimental gas‐phase band origins are all reproduced to better than 10 cm^?1. Consideration of the difference between a fully anharmonic zero‐point vibrational energy and B3LYP/cc‐pVTZ harmonic frequencies scaled by 0.985 suggests that the calculation of anharmonic zero‐point vibrational energies can generally be dispensed with, even in benchmark work, for rigid molecules. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1297–1305, 2001     

A comparative computational study of the homolytic and heterolytic bond dissociation energies involved in the activation step of ATRP and SET‐LRP of vinyl monomers

A quantum‐chemical calculation of the homolytic and heterolytic bond dissociation energies of the model compounds of the monomer and dimer is reported. These model compounds include the dormant chloride, bromide, and iodide species for representative activated and nonactivated monomers containing electron‐withdrawing groups as well as for a nonactivated monomer containing an electron‐donor group. Two examples of sulfonyl and N‐halide initiators are also reported. The homolytic inner‐sphere electron‐transfer bond dissociation is known as atom transfer and is responsible for the activation step in ATRP. The heterolytic outer sphere single electron transfer bond dissociation is responsible for the activation step in single electron transfer mediated living radical polymerization (SET‐LRP). The results of this study demonstrated much lower bond dissociation energies for the outer sphere single electron transfer processes. These results explain the higher rate constant of activation, the higher apparent rate constant of propagation, and the lower polymerization temperature for both activated and nonactivated monomers containing electron‐withdrawing groups in SET‐LRP. © 2007 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 1607–1618, 2007     

Theoretical investigation on electronic structure and stability of some inverted cucurbiturils by density functional theory

Structure and stability of diastereoisomers of cucurbit[n]urils (CB[n = 5–10]), the inverted CB[n]s, were investigated by density functional theory (DFT) computations. All the inverted CB[n]s were less stable than their normal CB[n]s and the mono-inverted ones with one inverted glycoluril unit in their structures were more stable than their doubly-inverted isomers. Relative change in dipole moments and molecular electrostatic potentials (MEP) were discussed with the deformation in geometric structure and charge distribution of the normal and inverted CB[n]s.     

Conversion of N2O to N2 on MgO （001） Surface with Vacancy： A DFT Study

The adsorption and decomposition of NzO at regular and defect sites of MgO (001) surface have been studied using cluster models embedded in a large array of point charges (PCs) by DFT/B3LYP method. The results indicate that the MgO (001)surface with oxygen vacancies exhibits high catalytic reactivity toward N2O adsorptive-decomposition. It is different from the regular MgO surface or the surface with magnesium vacancies.Much elongation of O—N bond of N2O after adsorption at oxy-gen vacancy site with O end down shows that O—N bond has been broken with concurrent production of N2, leaving a regu-lar site instead of the original oxygen vacancy site (F center ).The MgO (001) surface with magnesium vacancies hardly ex-hibits catalytic reactivity. It can be concluded that N2O dissoci-ation likely occurs at oxygen vacancy sites of MgO (001) sur-face, which is consistent with the generally accepted viewpoint in the experiments. The potential energy surface (PES) reflects that the dissociation process of N2O does not virtually need to surmount a given energy barrier.     

Density functional theory study of calix[4]arene‐N‐azacrown‐5, calix[4]arene‐N‐phenyl‐azacrown‐5, and their complexes with alkali‐metal cations: Na+, K+, and Rb+

Theoretical studies of 1,3‐alternate‐25,27‐bis(1‐methoxyethyl)calix[4]arene‐azacrown‐5 ( L₁ ), 1,3‐alternate‐25,27‐bis(1‐methoxyethyl)calix[4]arene‐N‐phenyl‐azacrown‐5 ( L₂ ), and the corresponding complexes M⁺/ L of L₁ and L₂ with the alkali‐metal cations: Na⁺, K⁺, and Rb⁺ have been performed using density functional theory (DFT) at B3LYP/6‐31G* level. The optimized geometric structures obtained from DFT calculations are used to perform natural bond orbital (NBO) analysis. The two main types of driving force metal–ligand and cation–π interactions are investigated. The results indicate that intermolecular electrostatic interactions are dominant and the electron‐donating oxygen offer lone pair electrons to the contacting RY* (1‐center Rydberg) or LP* (1‐center valence antibond lone pair) orbitals of M⁺ (Na⁺, K⁺, and Rb⁺). What's more, the cation–π interactions between the metal ion and π‐orbitals of the two rotated benzene rings play a minor role. For all the structures, the most pronounced changes in geometric parameters upon interaction are observed in the calix[4]arene molecule. In addition, an extra pendant phenyl group attached to nitrogen can promote metal complexation by 3D encapsulation greatly. In addition, the enthalpies of complexation reaction and hydrated cation exchange reaction had been studied by the calculated thermodynamic data. The calculated results of hydrated cation exchange reaction are in a good agreement with the experimental data for the complexes. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

气相中2Sr+和2Ba+催化N2O与CO还原反应的理论研究

用密度泛函理论DFT-UB3LYP方法,对2Sr+,2Ba+采用相对论校正赝势基组SDD,对C,N,O采用6-311+G(2d)基组,计算研究了气相中碱土金属离子2Sr+,2Ba+催化N2O(X1∑+)+CI(X1∑+)→N2(X1∑+g)+CO2(X1∑+g)反应的微观机理.优化了二重态势能面上各反应物、中间体和过渡...     

Density functional theory study of tautomerization of 2‐aminothiazole in the gas phase and in solution

The amino/imino tautomeric equilibrium in the isolated, mono‐, di‐, and trihydrate forms and dimer of 2‐aminothiazole, and the effects of hydration or self‐assistance on the transition state structures corresponding to proton transfer from the amino to imino form, have been investigated by the B3LYP method in conjunction with 6‐31+G(d,p) and 6‐311+G(3df,2p) basis sets in the gas phase and in solution. The amino form has been found to be the predominant tautomer. The tautomeric barrier heights for water‐ and self‐assisted tautomerization reactions are significantly lower than that from the amino to imino form by the intramolecular proton transfer, showing the catalytic effect of water molecule(s) and the important role of 2‐aminothiazole itself for intermolecular proton transfer. Comparison between the tautomeric barriers demonstrates that the self‐association tautomerization through the dimerization is the most favorable pathway. Bulk solvent effects have been taken into account using the polarizable continuum model (PCM) of water and CCl₄. The polar medium is favorable for the population of the imino form. The amino/imino equilibrium is also analyzed using the aromaticity index nucleus‐independent chemical shift (NICS); the NICS values for the amino form (about ?10 ppm) are more negative than the imino species (about ?8 ppm), showing that the amino form is more stable. The time‐dependent density functional theory (TDDFT) calculations of electronic absorption spectra suggest that the λ_max of dimer is 255 nm. The oscillator strength of the imino forms is less than the amino form, and increases with the polarity of the solvents. All calculations for the tautomerization of 2‐aminothiazole are in reasonable line with the available experiments. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A theoretical investigation on palladium‐catalyzed one‐pot coupling of aryl iodides,alkynes, and amines through CN bond cleavage for the synthesis of indole derivatives

Detailed mechanism of one‐pot coupling of aryl iodides, alkynes, and amines catalyzed by a palladium complexes is studied using density functional theory. The reaction of three‐component proceeds in four steps after the insertion of alkyne into C? I bond. The four‐component case comprises eight steps incorporating two times of nucleophilic addition. The C? N bond cleavage is rate‐limiting. The reaction sequence is maintained the same. The function of metal Pd is reflected by the coordination ability and electrophility of different oxidation state. The high selectivity of ligand is originated from the style of electron distribution in C? N bond cleavage and amination steps. We also discuss the amination from a bonding orbital matching perspective by means of frontier molecular orbital analysis. © 2014 Wiley Periodicals, Inc.