High acidity of polycyano azatriquinanes—theoretical prediction by the DFT calculations

It is shown by the B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d) calculations that the hexacyano derivative of aza-acepentalene is an extremely powerful superacid both in the gas phase and in DMSO as evidenced by the ΔH_acid = 255.1  kcal mol⁻¹ and pK_a (DMSO) = −26.5. Its synthesis is strongly recommended, in particular, since the related conjugate base hexachloro aza-acepentalenide anion was prepared recently.     

Theoretical study of intermolecular proton transfer reaction in isolated 5-hydroxyisoxazole–water complexes

A systematic investigation in isolated 5-hydroxyisoxazole–water complexes (5-HIO · (H₂O)_nn = 1–3) is performed at the DFT level, employing B3LYP/6-31G(d, p) basis set. Single-point energy calculations are also performed at the MP2 level using B3LYP/6-31G(d, p) optimized geometries and the 6-311++G(d, p) basis set. The computational results show that the keto tautomer K₂ is the most stable isomer in the gas phase, and the tautomer K₁ to be the next most stable tautomer. Hydrogen bonding between HIO and the water molecule(s) will dramatically lower the barrier by a concerted multiple proton transfer mechanism. The proton transfer process of 3WE_cis ↔ 3WK₁ and 2WE_trans ↔ 2WK₂ is found to be more efficient in two tautomerization, and the barrier heights are 7.03 and 14.15 kcal/mol at B3LYP/6-31G(d, p) level, respectively. However, the proton transfer reaction between E_cis and K₁ cannot happen without solvent-assisted.     

Dimeric water embedded within a hydrophobic cavity of tetra-(p-tert-butyl)thiacalix[4]arene

X-ray crystallographic analysis and density functional B3LYP/6-31G(d) calculation confirm that dimeric water is embedded within a hydrophobic cavity of tetra-(p-tert-butyl)thiacalix[4]arene and stabilized by hydrogen bondings of aromatic π?H₂O(1) and methyl?H₂O(2) in the dimeric water inclusion complex with binding energies of 1.4 and 0.9 kcal mol⁻¹ respectively, and by hydrogen bonding formed between H₂O(2) and four phenolic OH groups from an adjacent tetra-(p-tert-butyl)thiacalix[4]arene, with binding energy of 3.8-4.2 kcal mol⁻¹.     

Preparation and evaluation of 2-azinyl-2H-benzotriazoles as bidentate ligands: Synthesis and characterization of [2-(2-pyridynyl)-2H-benzotriazole](bpy)2Ru

Five 2-azinyl-2H-benzotriazoles (azinyl = 2-pyridinyl, 2-pyrazinyl, 2-pyrimidinyl, 6-methoxy-3-pyridazinyl, 5-methyl-2-pyridinyl were prepared and characterized as bidentate ligands. The electronic structure of these and related heterocycles was investigated spectroscopically and computationally (TD-DFT). They were tested at the B3LYP/6-31++G(d, p)//B3LYP/6-31G(d, p) level of theory as ligands for MgH₂, which permitted the elucidation of trends in complex formation, its geometry as a function of the ring structure, and the number and position of the nitrogen atoms in the azine ring. A Ru²⁺ complex 7a-Ru with 2-pyridinyl-2H-benzotriazole (7a) and two bpy ligands was prepared and characterized structurally, spectroscopically and electrochemically. The results were compared to those for similar complexes and discussed in the context of computational results for MgH₂ complexes.     

Are closed clusters expected from the (n + 1) skeletal electron pairs rule in alanes and gallanes? A DFT structural study of AnHn + 2 (A = Al, Ga, and n = 4-6)

It has been suggested recently that the alanes Al_nH_n + 2 can be treated by the polyhedral skeletal electron pair theory (PSEPT) of Wade and Mingos (W-M) as it was successful for their borane congeners such as B_nH_n + 2, well known as the deprotonated B_nH_n²⁻. To do so, the neutral Al_nH_n + 2 have been considered as Al_nH_n²⁻ + 2H⁺. The additional hydrogens donate their electrons to the Al_nH_n polyhedral framework and according to the n + 1 electron pairs rule; these clusters should have closo-polyhedral structures. In this work the homologous gallanes, the structures and stabilities of Ga_nH_n + 2 are studied at high levels of calculational theory and we investigated the applicability of the W-M rule to the alanes and gallanes A_nH_n + 2 (n = 4-6; A = Al, Ga). It will be shown that the presence of bridging hydrogen atoms reduces the compactness of the corresponding polyhedron and so these species do not have the closed structures. The computations were performed at B3LYP/6-311+G(d,p), BPW91/6-311G(d,p) and B3LYP/6-311+G(3df,2p) levels of theory. Our interest in these compounds includes their potential use as hydrogen storage species and future clean sources of energy.     

The isomers of phenalene and their singlet and triplet states: A Hartree–Fock and density functional computational investigation

Each of the isomers of phenalene, 1H-, 2H-, 3aH-, and 9bH-phenalene, as well as the cation, neutral radical, and anion in the phenalenyl system, have been examined at the Hartree–Fock 6-31G(d) and the density functional B3LYP/6-31G(d) levels of theory. The structures and properties of the phenalenes were determined as both singlets and triplets. While the data indicate that both the 9bH- and the 1H-isomer will exist as ground state singlets, the 2H- and 3aH-phenalenes are predicted to exist as ground state triplets; only the synthesis of the 1H-isomer has been reported in the literature. Structurally, 1H- and 2H-phenalene are planar systems, 3aH- and 9bH-phenalene are non-planar systems puckered at the saturated carbon, and the cation, neutral radical, and anion of phenalenyl are planar D_3h systems.     

The phonon-assisted proton transfer between hydrazinium cations in (N2H5)2HMF6·2H2O (M∈{Ga, Al, Fe}) type compounds. FT IR and quantum theoretical study

FT IR spectra of a series of compounds with a general formula (N₂H₅)₂HMF₆·2H₂O (where M∈{Ga, Al, Fe}) were recorded at variable temperatures (from ∼100 to 300 K, at 10 K intervals). The appearance of the spectral region of ν(N-N) modes due to hydrazinium cations further supports the conclusions regarding the N₂H₅⁺?H⁺?N₂H₅⁺ hydrogen bond potential well based on Raman spectroscopic data [J. Raman Spectrosc. 28 (1997) 315]. The appearance of two bands corresponding to the ν(N-N) modes in the low temperature FT IR spectra that merge into one upon heating is a clear evidence of a symmetric potential well through which a phonon-assisted proton transfer (PAPT) occurs at higher temperatures. Ab initio MP2/6-311++G(2d,p) quantum chemical study of the proton transfer potential within the N₂H₅⁺?H⁺?N₂H₅⁺ cluster confirmed its double-minimum character. The first-order saddle point found on the MP2/6-311++G(2d,p) potential energy hypersurface corresponds to a centrosymmetric structure (C_2h symmetry), with the proton placed at the inversion center. The potential energy curve along the tunnelling coordinate was calculated by the intrinsic reaction coordinate (IRC) methodology, leading to an adiabatic PT barrier height of 3.94 kcal mol⁻¹ and a tunneling rate of 1.98 s⁻¹. The corresponding MP4(SDTQ)/6-311++G(2d,p)//MP2/6-311++G(2d,p) value of the adiabatic PT barrier height is 4.26 kcal mol⁻¹.     

Infrared and Raman spectra, conformational stability and vibrational assignment of 1-chloro-1-silacyclopentane

Infrared and Raman spectra (3500-60 cm⁻¹) of gas and/or liquid and solid 1-chloro-1-silacyclopentane (c-C₄H₈SiClH) have been recorded and the vibrational data indicate the presence of a single conformer with no symmetry which is consistent with the twisted form. Ab initio calculations with a variety of basis sets up to MP2(full)/aug-cc-pVTZ predict the envelope-axial and envelope-equatorial conformers are saddle points with nearly the same energies but much lower in energy than the planar conformer. Density functional theory calculations by the B3LYP method predicts slightly lower energies for the two envelope forms and considerably lower for the planar form. By utilizing the MP2(full)/6-31G(d) calculations the force constants, frequencies, infrared intensities, band contours, Raman activities, and depolarization values have been obtained to support the vibrational assignment. Estimated r₀ structural parameters have been obtained from adjusted MP2(full)/6-311 + G(d, p) calculations. These experimental and theoretical results are compared to the corresponding quantities of some other five-membered rings.     

The theoretical research of the medium effect on the vibrational spectrum and the energy of intramolecular coordination O → Si in (benzoyloxymethyl)trifluorosilane molecule

The characteristics of (benzoyloxymethyl)trifluorosilane C₆H₅C(O)OCH₂SiF₃ containing a five-membered heterocycle closed by intramolecular coordination O → Si bond (Ia) and its most stable acyclic isomer (Ib) have been calculated by HF, MP2(Full) non-empirical methods, and DFT(B3LYP) using 6-311G(d) and 6-311 + G(2d,p) basis sets. The (C₈H₁₈, C₆H₆, (C₄H₉)₂O, CHCl₃, (CH₂)₄O, CH₂Cl₂, CH₃CN) medium effect on the energy and structural characteristics, dipole moments, and vibrational spectra of Ia and Ib isomers was calculated by the DFT(B3LYP)/6-311 + G(2d,p) method in the Onsager SCRF model approximation. The DFT(B3LYP)/6-311 + G(2d,p) calculation reasonably reproduces the medium effect on coordination energy, geometry, dipole moments, and band frequencies in the vibrational spectrum of Ia.     

Improper hydrogen-bonded cyclohexane C-Hax?Yax contacts: experimental evidence from H NMR spectroscopy of suitable axial cyclohexane models

B3LYP/6-31+G(d,p) calculations predicted the presence of improper hydrogen-bonded C-H_ax?Y_ax contacts of different strength in cyclohexane derivatives;¹ it was predicted that the addition of an appropriate bridging fragment X_ax between an axial substituent Y₁ and a cyclohexane carbon would strengthen the improper hydrogen-bonded contact C-H_ax?Y₁ when the X_ax-Y₁ bond vector bisects the cyclohexane ring. To support the theoretical predictions with experimental evidence for this effect, several 2-substituted adamantane analogues with suitable improper H-bonded C-H_ax?O contacts of different strength were synthesized, as models of the corresponding cyclohexane derivatives, and their ¹H NMR spectra were recorded at 298 K. The ¹H NMR signal separation within the cyclohexane ring γ-CH₂s is increased when the B3LYP/6-31+G(d,p)-calculated strength of the H-bonded C-H_ax?O=C_ax contact interaction is increased.     

Theoretical modelling of the epoxidation of vinylallenes to give cyclopentenones

Epoxidation of vinylallenes (1,2,4-pentatrienes) can lead to cyclopent-2-enones. Various experimental results suggest that these reactions are concerted and that vinylallene oxides undergo a concerted and thermal ring-closing reaction to give cyclopentenones. In 1977, this view has been supported by the epoxidation of a non-racemic 4-methyl-2,3,5-hexatriene to give a non-racemic 2,5-dimethyl-2-cyclopenten-1-one. Indeed, for the cyclization of (3E)-4-methyl-2,3-epoxy-2,3,5-hexatriene or (4E)-4-methyl-3,4-epoxy-2,3,5-hexatriene into 2,5-dimethyl-2-cyclopenten-1-one, a transition structure for the concerted rearrangement was located, and IRC calculations showed it linked together. The activation barrier predicted at the B3LYP/6-311++G(3df,3pd)//B3LYP/6-31G(d) level of theory is 80.9 or 124.2 kJ mol⁻¹, respectively and the reaction is exothermic by 127.2 or 89.2 kJ mol⁻¹, respectively.     

Experimental and quantum chemical study of pyrrole self-association through N-H?π hydrogen bonding

An FT-IR study of pyrrole self-association in CCl₄ solutions was carried out. According to the IR measurements, pyrrole forms self-associated dimeric species via N-H?π hydrogen bonding. This was also confirmed by quantum chemical calculations for pyrrole monomer and dimer at B3LYP/6-31++G(d,p) level of theory. A T-shaped minimum was located on B3LYP/6-31++G(d,p) PES of pyrrole dimer characterized with a hydrogen bond of an N-H?π type, with centers-of-mass separation of monomeric units of 4.520 Å, H?π distance of 2.475 Å, the interplanar angle between the two monomeric units being 72.9°. The anharmonic vibrational frequency shift upon dimer formation calculated on the basis of 1D DFT vibrational potentials is in excellent agreement with the experimental data (84 vs. 87 cm⁻¹). Harmonic vibrational analysis predicts somewhat smaller shift (68 cm⁻¹). On the basis of NIR spectroscopic data, anharmonicity constants for the 2ν(N-H) and 2ν(N-H?π) vibrational transitions were calculated. The orientational dynamics of monomeric and self-associated pyrrole species was studied within the framework of the transition dipole moment time correlation function formalism. The period of essentially free rotation in the condensed phase reduces from 0.05 ps for the monomeric pyrrole to 0.02 ps for the proton-donor molecule within the dimer.     

The substituent effects of the leaving groups on the aminolysis of phenyl acetates: DFT studies

The density functional theory B3LYP/6-31G(d, p) method is employed to study the mechanism of aminolysis reaction of p-substituted phenyl acetates (CH₃C(O)OC₆H₄X, X = H, NH₂, and NO₂) with ammonia in the gas phase. Two reaction pathways are considered: the concerted process and the stepwise pathway through neutral intermediates. The substituent effects of the leaving groups on the reactivity of phenyl acetates are discussed. The solvent effect of acetonitrile on the title reaction is also assessed by the polarizable continuum model (CPCM model) at B3LYP/6-31++G(d, p) level of theory. The calculated results show that the activation barriers of the concerted pathways are lower than those of the rate-controlling steps of the stepwise processes for all the three aminolysis reactions. This aminolysis of phenyl acetates is more favorable for X = NO₂ than for X = H and NH₂ in the gas phase and in acetonitrile.     

Proton affinities of ketones, vicinal diketones and α-keto esters: a computational study

The proton affinities of seven different ketones, vicinal diketones, and α-keto esters (acetophenone, 2,2,2-trifluoroacetophenone, 2,3-butanedione, 1-phenyl-1,2-propanedione, methyl pyruvate, ethyl benzoylformate, and ketopantolactone) have been evaluated theoretically using the conventional ab initio HF and several post-HF methods (MP2, MP4, CCSD), density functional methods with the B3LYP hybrid functional, as well as some ab initio model chemistries [CBS-4M, G2(MP2), and G3(MP2)//B3LYP]. The chemical compounds studied are frequently used substrates in the asymmetric hydrogenation over chirally modified platinum catalysts where the protonation properties of the chiral modifier and the substrates are of great interest. In most cases, the proton affinities (PAs) evaluated with the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) methods are in good agreement with the existing experimental ones. However, the previously reported PA of 2,3-butanedione seems to be too high by 10-15 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP model chemistries predict proton affinities that are systematically higher and lower than the experimental PAs, respectively. If proton affinities are evaluated as the average of the PAs calculated with these two theoretical methods a very good agreement with the experimental results is obtained. The mean absolute deviation (MAD) from experiment of this combination method for the PAs of 13 test molecules is 4.0 kJ mol⁻¹. For 9 molecules composed only of first-row atoms the MAD is 2.5 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP methods provide significant savings in computational time and disk space compared to the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) models. Therefore, it is suggested that if no experimental or highly accurate theoretical data is available (due to computational cost), the proton affinities of similar compounds as investigated in this paper, can be evaluated with the combination method. For the studied molecules, this method gives the following PAs (in kJ mol⁻¹): 788 (2,3-butanedione, exptl 802); 798 (2,2,2-trifluoroacetophenone, exptl 799); 811 (ketopantolactone); 813 (methyl pyruvate); 825 (1-phenyl-1,2-propanedione); 862 (acetophenone, exptl 861); 865 (ethyl benzoylformate).     

DFT and ab initio theoretical study for the CF3S + CO reaction

The potential energy surface for the reaction of CF₃S with CO is calculated at the G4//B3LYP/6-311++G(d,p) level of theory. The results show that F-abstraction and addition-elimination mechanisms are involved, and the latter one is dominant thermodynamically and kinetically. The dominant channel is the reactant addition to form CF₃SCO, and then decomposes to CF₃ + OCS. While the direct F-abstraction channel and CF₃SCO isomerization channel are not significant for the title reaction due to higher barriers involved. The comparisons among four reactions of CX₃Y + CO (X = H, F; and Y = O, S) are made to imply the similar and different properties and reactivities of the same family elements and the F- and S-substituted derivatives.     

Ab initio/DFT theory and multichannel RRKM study on the mechanisms and kinetics for the CH3S + CO reaction

The potential energy surface for the reaction of CH₃S with CO was calculated at the G3MP2//B3LYP/6-311++G(d,p) level. The rate constants for feasible channels leading to several products were calculated by TST and multichannel-RRKM theory. The results show that addition–elimination mechanism is dominant, while hydrogen abstraction mechanism is uncompetitive. The major channel is the addition of CO to CH₃S leading to an intermediate CH₃SCO which then decomposes to CH₃ + OCS. In the temperature range of 200–3000 K, the overall rate constants are positive temperature dependence and pressure independence, and it can be described by the expression as k = 1.10 × 10⁻¹⁶T^1.57exp(−3359/T) cm³ molecule⁻¹ s⁻¹. At temperature between 208 and 295 K, the calculated rate constants are in good agreement with the experimental upper limit data. At T = 1000 and 2000 K, the major product is CH₃ + OCS at lower pressure; while at higher pressure, the stabilization of IM1 is dominant channel.     

Synthesis, crystal structure and DFT calculations on 2,6-diisopropylphenylcopper; its use in the preparation of dichloro-2,6-diisopropylphenylphosphine

The homoleptic aryl copper reagent [Cu₄Dipp₄] (Dipp = 2,6-diisopropylphenyl) has been prepared and structurally characterized by a single-crystal X-ray diffraction study. Its tetrameric structure differs in significant details from that of the previously reported [Cu₄Tripp₄] (Tripp = 2,4,6-triisopropylphenyl). The electronic structure of the cluster has been probed through B3LYP/6311G(2d,p)//B3LYP/6-31G calculations on [Cu₄Ph₄] constrained to D_2d symmetry. The utility of the new copper reagent is demonstrated by the preparation of pure DippPCl₂, for which the crystal structure is also reported.     

Origin of the stereo- and regioselectivities in the Diels-Alder reactions of azaphospholes: a DFT investigation

Computations of the Diels-Alder (DA) reactions of azaphosphole representative namely, thiazolo[3,2-d][1,4,2]diazaphosphole with 1,3-butadiene and isoprene at the density functional theory level reveal concerted mechanisms via asynchronous transition states. The activation energies (B3LYP/6-311++G**// B3LYP/6-311G**), 16-19 kcal mol⁻¹, are much smaller than the value (32.57 kcal mol⁻¹) calculated for the DA reaction of the non-phosphorus analogue, imidazo[2,1-b]thiazole with 1,3-butadiene. An electron-withdrawing group at the 3-position of the dienophile enhances both stereo- and regioselectivities, which agree nicely with the experimental values. Inclusion of solvent effect (PCM model) reveals that the stereo- and regioselectivities are not affected appreciably. The relative stabilities of the transition structures corresponding to the endo/exo stereoisomers and meta (P/Me, 1:3)/para (P/Me, 1:4) regioisomers have been rationalized on the basis of the secondary molecular orbital interactions.     

Computational study of keto–enol equilibria of tropolone in gas and aqueous solution phase

Density functional theory calculations employing the B3LYP exchange-correlation functional, as well as Hartree–Fock computations, were performed on 2-hydroxy-2,4,6-cycloheptatrien-1-one (tropolone) and 3,5- and 3,6-cycloheptadiene-1,2-dione in gas and aqueous solution phases in order to determine the equilibrium constant for keto to enol interconversion of the isomers of C₇H₆O₂. Two standard basis sets were used throughout, namely 6-311++G^∗∗ and aug-cc-pVDZ. Solvent effects were modelled using two different self-consistent reaction field approaches – the Onsager dipole and the polarizable continuum models (PCM). In addition, the G3 method was used for calculations on species in the gas phase. Molecular geometries were fully optimized at each model chemistry, and it was found that the two keto isomers are always higher in energy than the enol form. From the results of B3LYP/6-311++G^∗∗ calculations of the difference in Gibbs free energy in the gas phase and using PCM, the relative pK values for the 2-hydroxy-2,4,6-cycloheptatrien-1-one ? 3,5- and 3,6-cycloheptadiene-1,2-dione system are 13.75 (g), 15.78 (g) and 13.05 (aq) and 13.45 (aq), respectively. That equilibrium is tilted almost exclusively in the direction of tropolone is due to resonance stabilization of the enol as a result of aromaticity, and is most easily understood on the basis of elementary Hückel theory.     

The infrared and Raman spectra, ab initio calculations and spectral assignments of cyclopropylmethyl dichlorosilane (c-C3H5)SiCl2CH3

Raman spectra of cyclopropylmethyl dichlorosilane (c-C₃H₅)SiCl₂CH₃ as a liquid were recorded at 293 K and polarization data were obtained. Additional Raman spectra were recorded at various temperatures between 293 and 163 K, and intensity changes of certain bands with temperature were detected. No crystallization was ever obtained in the Raman cryostat in spite of extensive annealing. The infrared spectra have been studied as a vapour, as an amorphous solid at 78 K and as a liquid in the range 600-100 cm⁻¹. No infrared bands present in the vapour or liquid seemed to vanish upon cooling, and the sample never formed crystals on the CsI window of an infrared cryostat.The compound exists a priori in two conformers, syn and gauche, and the experimental results suggest an equilibrium in which the gauche conformer has 1.64 kJ mol⁻¹ lower enthalpy than syn in the liquid, leading to 20% syn at ambient temperature. Most of the syn bands were situated close to the corresponding gauche bands and it was difficult to obtain reliable ΔH values.B3LYP calculations with various basis sets and the CBS-QB3 and G2 and G3 models were employed, yielding the conformational enthalpy difference ΔH (syn-gauche) between 2.6 and 3.4 kJ mol⁻¹. Infrared and Raman intensities, polarization ratios and vibrational frequencies for the syn and gauche conformers were calculated. Instead of scaling the calculated wavenumbers in the harmonic approximation, calculations from B3LYP/cc-pVTZ were derived in the anharmonic approximation. In most cases these values were in good agreement with the experimental results for 38 observed modes of the gauche and 8 modes of the syn conformer with a deviation of ca. 1%.