DFT investigation of the mechanism of CH2CO + O(3P) reaction

A detailed theoretical survey of the potential energy surface (PES) for the CH₂CO + O(³P) reaction is carried out at the QCISD(T)/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) level. The geometries, vibrational frequencies, and energies of all stationary points involved in the reaction are calculated at the B3LYP/6‐311+G(d,p) level. More accurate energy information is provided by single‐point calculations at the QCISD(T)/6‐311+G(3df,2p) level. Relationships of the reactants, transition states, intermediates, and products are confirmed by the intrinsic reaction coordinate (IRC) calculations. The results suggest that P1(CH₂+CO₂) is the most important product. This study presents highlights of the mechanism of the title reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio quantum chemical studies of the reactions of CF3CFHO2 with HO2

The potential energy surface (PES) for the CF₃CFHO₂+HO₂ reaction has been theoretically investigated using the DFT [B3LYP/6‐311G(d,p)] and B3LYP/6‐311++G(3df,3pd)//B3LYP/6‐311G(d,p) levels of theory. Both singlet and triplet PESs are investigated. The reaction mechanism on the triplet surface is simple. It is revealed that the formation of CF₃CFHOOH+³O₂ is the dominant channel on the triplet surface. On the basis of the ab initio data, the total rate constants for the reaction CF₃CFHO₂+HO₂ in the T = 210–500 K range have been computed using conventional transition state theory with Wigner's tunneling correction and have been fitted by a rate constant expression as k = 1.04 ×10^?12(cm³ molecule^?1 s^?1) exp (700.33/T). Calculated transition state rate constants with Wigner's tunneling correction for the reaction CF₃CFHO₂+HO₂ are in good agreement with the available experimental values. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Ab initio study of the OH + CH2O reaction: The effect of the OH··OCH2 complex on the H‐abstraction kinetics

Kinetics for the reaction of OH radical with CH₂O has been studied by single‐point calculations at the CCSD(T)/6‐311+G(3df, 2p) level based on the geometries optimized at the B3LYP/6‐311+G(3df, 2p) and CCSD/6‐311++G(d,p) levels. The rate constant for the reaction has been computed in the temperature range 200–3000 K by variational transition state theory including the significant effect of the multiple reflections above the OH··OCH₂ complex. The predicted results can be represented by the expressions k₁ = 2.45 × 10^‐21 T^2.98 exp (1750/T) cm³ mol^?1 s^?1 (200–400 K) and 3.22 × 10^‐18 T^2.11 exp(849/T) cm³ mol^?1 s^?1 (400–3000 K) for the H‐abstraction process and k₂ = 1.05 × 10^‐17 T^1.63 exp(?2156/T) cm³ mol^?1 s^?1 in the temperature range of 200–3000 K for the HO‐addition process producing the OCH₂OH radical. The predicted total rate constants (k₁ + k₂) can reproduce closely the recommended kinetic data for OH + CH₂O over the entire range of temperature studied. © 2006 Wiley Periodicals, Inc. Int J Chem Kinet 38: 322–326, 2006     

Theoretical predictions of the structure,gas‐phase acidity,and aromaticity of tetrathiosquaric acid

Results of ab initio self‐consistent‐field and density functional theory calculations of the gas‐phase structure, acidity (free energy of deprotonation, ΔG⁰), and aromaticity of tetrathiosquaric acid (3,4‐dithiohydroxy‐3‐cyclobutene‐1,2‐dithione, H₂C₄S₄) are reported. The global minimum found on the potential energy surface of tetrathiosquaric acid presents a planar conformation. The ZZ isomer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very close in energy. The optimized geometric parameters exhibit a bond length equalization relative to reference compounds, cyclobutanedithione, and cyclobutenedithiol. The computed aromatic stabilization energy by homodesmotic reaction is −18.4 (MP2(fu)/6‐311+G**//RHF/6‐311+G**) and −15.1 kcal/mol (B3LYP//6‐311+G**// B3LYP/6‐311+G**). The aromaticity of tetrathiosquaric acid is indicated by the calculated diamagnetic susceptibility exaltation (Λ) −11.77 (CSGT(IGAIM)‐RHF/6‐311+G**// RHF/6‐311+G**) and −18.08 (CSGT(IGAIM)‐B3LYP/6‐311+G**// B3LYP/6‐311+G**). Thus, tetrathiosquaric acid fulfils the geometric, energetic and magnetic criteria of aromaticity. The most reliable theoretical gas‐phase acidities are $\Delta G^{0}_{1(298\mathrm{K})}=303.7$ and $\Delta G^{0}_{2(298\mathrm{K})}=394.1$ kcal/mol. Hence, tetrathiosquaric acid is a stronger acid than squaric acid (3,4‐dihydroxy‐3‐ cyclobutene‐1,2‐dione, H₂C₄O₄). Comparisons of the computed results of tetrathiosquaric acid with squaric acid have also been made. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 443–449, 2000     

Complexation of Phenol and Thiophenol by Amine N‐Oxides: Isothermal Titration Calorimetry and ab Initio Calculations

To develop a new solvent‐impregnated resin (SIR) system for removal of phenols from water, the complex formation of dimethyldodecylamine N‐oxide (DMDAO), trioctylamine N‐oxide (TOAO), and tris(2‐ethylhexyl)amine N‐oxide (TEHAO) with phenol (PhOH) and thiophenol (PhSH) is studied. To this end we use isothermal titration calorimetry (ITC) and quantum chemical modeling (on B3LYP/6‐311G(d,p)‐optimized geometries: B3LYP/6‐311+G(d,p), B3LYP/6‐311++G(2d,2p), MP2/6‐311+G(d,p), and spin component scaled (SCS) MP2/6‐311+G(d,p); M06‐2X/6‐311+G(d,p)//M06‐2X/6‐311G(d,p), MP2 with an extrapolation to the complete basis set limit (MP2/CBS), as well as CBS‐Q). The complexes are analyzed in terms of structural (e.g., bond lengths) and electronic elements (e.g., charges). Furthermore, complexation and solvent effects (in benzene, toluene, and mesitylene) are investigated by ITC measurements, yielding binding constants K, enthalpies ΔH⁰, Gibbs fre energies ΔG⁰, and entropies ΔS⁰ of complex formation, and stoichiometry N. The ITC measurements revealed strong 1:1 complex formation between both DMDAO–PhOH and TOAO–PhOH. The binding constant (K=1.7–5.7×10⁴ M ^?1) drops markedly when water‐saturated toluene was used (K=5.8×10³ M ^?1), and π–π interaction with the solvent is shown to be relevant. Quantum mechanical modeling confirms formation of stable 1:1 complexes with linear hydrogen bonds that weaken on attachment of electron‐withdrawing groups to the amine N‐oxide moiety. Modeling also showed that complexes with PhSH are much weaker than those with PhOH, and in fact too weak for ITC determination. CBS‐Q incorrectly predicts equal or even higher binding enthalpies for PhSH than for PhOH, which invalidates it as a benchmark for other calculations. Data from the straightforward SCS‐MP2 method without counterpoise correction show very good agreement with the MP2/CBS values.     

Does CH prefer a C2v rather than a Cs structure?

The closely related C_s ( 1 ) and C_2v ( 3 ) structures of CH have been reinvestigated at many ab initio levels using MP2/6-31G** and MP2/6-311 + + G(2df, 2pd) geometries. The largest basis sets employed were 6-311G(3df, 2p), 6-311 + + G(3df, 3pd), and the Dunning “correlation consistent” polarized triple-split valence basis set (cc-pVTZ). Electron correlation was probed at the MP4 level, but the QCISD method was also used with the largest basis sets. While electron correlation favors 3 over 1 by about 2 kcal/mol, the correlated relative energies with all basis sets employed range from 0.36–1.03 kcal/mol in favor of 1 . The best estimate of this difference, 0.86 kcal/mol, is essentially identical with the (scaled) zero-point energy difference, 0.84 kcal/mol, favoring 3 over 1 . These results indicate that 1 and 3 have almost exactly the same energy at 0 K. Our best value for the dissociation energy of CH is 42.0 kcal/mol [QCISD(T)/6-311 + + G(3df, 3pd)//MP2(fu)/6-311 + + G(2df, 2pd), corrected to 298 K], which agrees very well with the experimental value. © 1992 by John Wiley & Sons, Inc.     

Theoretical study of the gas‐phase SN2 reactions of X− with CH3OY (X,Y = Cl,Br, I)

The gas‐phase nucleophilic substitution reactions at saturated oxygen X^? + CH₃OY (X, Y = Cl, Br, I) have been investigated at the level of CCSD(T)/6‐311+G(2df,p)//B3LYP/6‐311+G(2df,p). The calculated results indicate that X^? preferably attacks oxygen atom of CH₃OY via a S_N2 pathway. The central barriers and overall barriers are respectively in good agreement with both the predictions of Marcus equation and its modification, respectively. Central barrier heights (ΔH and ΔH) correlate well with the charges (Q) of the leaving groups (Y), Wiberg bond orders (BO) and the elongation of the bonds (O? Y and O? X) in the transition structures. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical investigation on multinuclear NMR chemical shifts of some tris(trifluoromethyl)boron complexes

Tris(trifluoromethyl)boron complexes have unusual properties and may find applications in many fields of chemistry, biology, and physics. To gain insight into their NMR properties, the isotropic ¹¹B, ¹³C, and ¹⁹F NMR chemical shifts of a series of tris(trifluoromethyl)boron complexes were systematically studied using the gauge‐included atomic orbitals (GIAO) method at the levels of B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G* and B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p). Solvent effects were taken into account by polarizable continuum models (PCM). The calculated results were compared with the experimental values. The reason that the structurally inequivalent fluorine atoms in a specific species give a same chemical shift in experimental measurements is attributed to the fast rotation of CF₃ group around the B? C(F₃) bond because of the low energy barrier. The calculated ¹¹B, ¹³C(F₃), and ¹⁹F chemical shifts are in good agreement with the experimental measurements, while the deviations of calculated ¹³C(X, X = O, N) chemical shifts are slightly large. For the latter, the average absolute deviations of the results from B3LYP/6‐311 + G(d,p)//B3LYP/6‐311 + G(d,p) are smaller than those from B3LYP/6‐31 + G(d,p)//B3LYP/6‐31G*, and the inclusion of PCM reduces the deviation values. The calculated ¹⁹F and ¹¹B chemical shieldings of (CF₃)₃BCO are greatly dependent on the optimized structures, while the influence of structural parameters on the calculated ¹³C chemical shieldings is minor. Copyright © 2009 John Wiley & Sons, Ltd.     

A DFT study of hydride transfers to the carbonyl oxygen of DDQ

Hydride‐transfer reactions between benzylic substrates and 2,3‐dichloro‐5,6‐dicyano‐1,4‐benzoquinone (DDQ) were investigated by DFT (density functional theory) calculations. The lowest unoccupied molecular orbital of DDQ has the largest extension on two carbonyl oxygens, which comes from two‐step mixing of antisymmetric orbitals of fragment π MOs. Transition‐state (TS) geometries and activation energies of reactions of four benzylic substrates R²? CH₂‐para‐C₆H₄? R¹ (R¹, R² = H and/or OCH₃) with DDQ were calculated. M06‐2X/6‐311(+*)G* was found to be a practical computational method, giving energies and geometries similar to those of M06‐2X/6‐311++G(3df,2pd) and wB97xD/6‐311++G(3df,2pd). For toluene (R¹ = R² = H), an initiation‐propagation model was suggested, and the calculated kinetic isotope effect k(H)/k(D) = 5.0 with the tunnel correction at the propagating step is in good agreement with the experimental value 5.2. A reaction of para‐MeO? C₆H₄? CH₂(OMe) + DDQ + (H₂O)₁₄ → para‐MeO? C₆H₄? C(?O)H + HOMe + DDQH₂ + (H₂O)₁₃ was investigated by M06‐2X/6‐311(+*)G*. Four elementary processes were found and the hydride transfer (TS1) is the rate‐determining step. The hydride transfer was promoted by association with the water cluster. The size of the water cluster, (H₂O)_n, at TS1 was examined. Three models of n = 14, 20, and 26 were found to give similar activation energies. Metal‐free neutral hydride transfers from activated benzylic substrates to DDQ were proposed to be ready processes both kinetically and thermodynamically. © 2015 Wiley Periodicals, Inc.     

Ab initio study on the reaction mechanism of ozone with bromine atom

Ab initio calculations of the potential energy surface (PES) for the Br+O₃ reaction have been performed using the MP2, CCSD(T), and QCISD(T) methods with 6‐31G(d), 6‐311G(d), and 6‐311+G(3df). The reaction begins with a transition state (TS) when the Br atom attacks a terminal oxygen of ozone, producing an intermediate, the bromine trioxide (M), which immediately dissociates to BrO+O₂. The geometry optimizations of the reactants, products, and intermediate and transition states are carried out at the MP2/6‐31G(d) level. The reaction potential barrier is 3.09 kcal/mol at the CCSD(T)/6‐311+G(3df)//MP2 level, which shows that the bromine atom trends intensively to react with the ozone. The comparison of the Br+O₃ reaction with the F+O₃ and Cl+O₃ reactions indicates that the reactions of ozone with the halogen atoms have the similar reaction mechanism. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Conformers and intramolecular hydrogen bonding of the oxalic acid monomer and its anions

Density functional theory, B3LYP/6‐31G** and B3LYP/6‐311+G(2d,p), and ab initio MP2/6‐31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISD/6‐31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJ/mol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar results, and that the conformer L1 (C_2h) with the double‐interfunctional‐groups hydrogen bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJ/mol for interfunctional groups, 17 kJ/mol for weak interfunctional groups, 24 kJ/mol for intra‐COOH in (COOH)₂, and 60 kJ/mol for interfunctional groups in (COOH)COO⁻¹ ion if calculated using the B3LYP/6‐311+G(2d,p) method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 541–551, 2000     

Theoretical studies on a possible synthesis reaction pathway on N8 (CS) clusters

The potential energy surfaces of N₈ clusters were investigated by density functional theory (DFT) and a possible synthesis reaction pathway for N₈ (C_S) was suggested. The species involved were fully optimized up to the B3LYP/6‐311+G* level of theory. Relative energies were further calculated at the QCISD/6‐311+G*//B3LYP/6‐311+G* level. The reaction rate constants of these steps from the 1 (N₅⁺?N₃^?, complex, C_S) to 2 (N₈, C_S), 2 (N₈, C_S) to 3 (N₈, C_S), 3 (N₈, C_S) to 4 (N₈, D_2d), and 4 (N₈, D_2d) to 5 (N₈, C_S) reactions were predicted by the VTST theory. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1334–1339, 2001     

[C5Li5]Mgn[C5Li5] (n?=?2–8): novel sandwich complexes containing –Mg–Mg– chain

A density functional theory (DFT) investigation on novel sandwich-type D ₅ [C₅Li₅]Mg _n [C₅Li₅] (n = 2–8) complexes containing –Mg–Mg– chain has been performed in this work. The equilibrium geometries, electronic structures, vibrational frequencies, and stabilities of these complexes are researched by B3LYP and BP86 methods at 6-311+G(d) levels of theory. The Mg _n ²⁺ sandwich complexes with D ₅ symmetry are all true minima on the potential energy surface. NBO analyses for the series of complexes reveal that the Mg–Mg bond is a weak σ covalent bond. There are mainly electrostatic interactions between C₅Li₅ ⁻ ligands and Mg _n ²⁺(n = 2–8) nuclear in these complexes. The NICS and NICS_zz computed with GIAO-B3LYP/6-311+G(d) indicates that the C₅Li₅ ⁻ rings in the series of complexes are aromatic. These novel complexes turn out to be strongly thermodynamically favored in the gas phases and may be targeted in future experiments to expand the structural domain of sandwich-type complexes.     

Exploration on the structure,stability, and isomerization of planar CnB5 (n = 1−7) clusters

The structures, stabilities, nature of bonding, and potential energy surfaces of low‐energy isomers of planar C_nB₅ (n = 1?7) have been systematically explored at the CCSD(T)/6‐311+G(d)//B3LYP/6‐311+G(d) level. Incremental binding energy (IBE) and second order energy difference (Δ²E) analyses demonstrate that C_nB₅ clusters with even n have relatively higher stability. The nature of bonding in these clusters is discussed based on valence molecular orbital (VMO), and Mayer bond order (MBO). Hückel (4n + 2) rule and nucleus‐independent chemical shift (NICS) values suggest that the ground states of C₃B₅, C₄B₅, and C₇B₅ have π aromaticity. VMO, electron localization function (ELF), adaptive natural density partitioning (AdNDP), and NICS analyses reveal the double aromaticity of C₃B₅ cation. CB₅ and C₃B₅ are stable both thermodynamically and kinetically based on isomerization analysis. In addition, the simulated IR spectra are expected to be helpful for future experimental studies of these clusters.     

Polynuclear sandwich derivatives of [10]annulene: A quantum chemical study

Density-functional theory (DFT) calculations at B3LYP/6-311+G(df,p) and M06/6-311+G(df,p) levels of theory predict the stability of new polynuclear sandwich complexes [Cu₅(C₁₀H₁₀)₂]⁺ and Ni₅(C₁₀H₁₀)₂ based on flat [10]annulene cycles.     

Conformational study of the structure of free 12-thiacrown-4 and some of its cation metal complexes

Conformational search of 12-thiacrown-4, 12t4, was performed using the CONFLEX method and the MMFF94S force field whereby 156 conformations were predicted. Optimized geometries of the 156 predicted conformations were calculated at the HF, B3LYP, CAM-B3LYP, M06, M06L, M062x and M06HF levels using the 6-311G** basis set. The correlation energy was recovered at the MP2 level using the same 6-311G** basis set. Optimized geometries at the MP2/6-311G** level and G3MP2 energies were calculated for some of the low energy conformations. The D ₄ conformation was predicted to be the ground state conformation at all levels of theory considered in this work. Comparison between the dihedral angles of the predicted conformations indicated that for the stability of 12t4, a SCCS dihedral angle of 180° requirement is more important than a gauche CSCC dihedral angle requirement. Conformational search was performed also for the 12t4?CAg⁺, Bi³⁺, Cd²⁺, Cu⁺ and Sb³⁺ cation metal complexes using the CONFLEX method and the CAChe-augmented MM3 and MMFF94S force fields. Conformations with relative energies less than 10?kcal/mol at the MP2/6-31+G*//HF/6-31+G* level, with double zeta quality basis set on the metal cations, were considered for computations at the same levels as those used for free 12t4, using also the 6-311G** basis set. The cc-pVTZ-pp basis set was used for the metal cations. The predicted ground state conformations of the 12t4?CAg⁺, Bi³⁺, Cd²⁺, Cu⁺ and Sb³⁺ cation metal complexes are the C ₄, C ₄, C ₄, C _2v and C ₄ conformations, respectively. This is in agreement with the experimental X-ray data for the 12t4?CAg⁺ and Cd²⁺ cation metal complexes, but experimentally by X-ray, the 12t4?CBi³⁺ and Cu⁺ cation metal complexes have C _s and C ₄ structures, respectively.     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.