C7H12(2+): a prototype hexacoordinate carbonium ion

C(7)H(12)(2+) (1), the prototype hexacoordinate carbonium dication was found to be a viable minimum at the MP2/6-31G** and MP2/cc-pVTZ levels. Structure 1 is a propeller shaped molecule resembling a complex involving a C(2+) with three ethylene molecules resulting in the formation of three two-electron, three-center (2e-3c) bonds. Isomeric structure 2 was found to be 21.8 kcal/mol more stable than structure 1. However, conversion of 1 into 2 through transition structure 3 has a barrier of 5.7 kcal/mol. Related structures 4, 5, and 8 were also located as minima for C(7)H(12)(2+). The isoelectronic boron analogue BC(6)H(12)(+) (10) was also computed to be a minimum at the same level of calculations.     

Photophysical studies of the trans to cis isomerization of the push-pull molecule: 1-(pyridin-4-yl)-2-(N-methylpyrrol-2-yl)ethene (mepepy)

Organic molecules possessing intramolecular charge-transfer properties (D-pi-A type molecules) are of key interest particularly in the development of new optoelectronic materials as well as photoinduced magnetism. One such class of D-pi-A molecules that is of particular interest contains photoswitchable intramolecular charge-transfer states via a photoisomerizable pi-system linking the donor and acceptor groups. Here we report the photophysical and electronic properties of the trans to cis isomerization of 1-(pyridin-4-yl)-2-(N-methylpyrrol-2-yl)ethene ligand (mepepy) in aqueous solution using photoacoustic calorimetry (PAC) and theoretical methods. Density functional theory (DFT) calculations demonstrate a global energy difference between cis and trans isomers of mepepy to be 8 kcal mol(-1), while a slightly lower energy is observed between the local minima for the trans and cis isomers (7 kcal mol(-1)). Interestingly, the trans isomer appears to exhibit two ground-state minima separated by an energy barrier of approximately 9 kcal mol(-1). Results from the PAC studies indicate that the trans to cis isomerization results in a negligible volume change (0.9 +/- 0.4 mL mol(-1)) and an enthalpy change of 18 +/- 3 kcal mol(-1). The fact that the acoustic waves associated with the trans to cis transition of mepepy overlap in frequency with those of a calorimetric reference implies that the conformational transition occurs faster than the approximately 50 ns response time of the acoustic detector. Comparison of the experimental results with theoretical studies provide evidence for a mechanism in which the trans to cis isomerization of mepepy results in the loss of a hydrogen bond between a water molecule and the pyridine ring of mepepy.     

Conformational stability, vibrational assignmenents, barriers to internal rotations and ab initio calculations of 2-aminophenol (d 0 and d3)

The Raman (3700-100 cm(-1)) and infrared (4000-400 cm(-1)) spectra of solid 2-aminophenol (2AP) have been recorded. The internal rotation of both OH and NH2 moieties produce ten conformers with either Cs or C1 symmetry. However, the calculated energies as well as the imaginary vibrational frequencies reduce rotational isomerism to five isomers. The molecular geometry has been optimized without any constraints using RHF, MP2 and B3LYP levels of theory at 6-31G(d), 6-311+G(d) and 6-31++G(d,p) basis sets. All calculations predict 1 (cis; OH is directed towards NH2) to be the most stable conformation except RHF/6-31++G(d,p) basis set. The 1 (cis) isomer is found to be more stable than 8 (trans; OH is away from the NH2 moiety and the NH bonds are out-of-plane) by 1.7 kcal/mol (598 cm(-1)) as obtained from MP2/6-31G(d) calculations. Aided by experimental and theoretical vibrational spectra, cis and trans 2AP are coexist in solution but cis isomer is more likely present in the crystalline state. Aided by MP2 and B3LYP frequency calculations, molecular force fields, simulated vibrational spectra utilizing 6-31G(d) basis set as well as normal coordinate analysis, complete vibrational assignments for HOC6H4NH2 and DOC6H4ND2 have been proposed. Furthermore, we carried out potential surface scan, to determine the barriers to internal rotations of NH2 and OH groups. All results are reported herein and compared with similar molecules when appropriate.     

Trans or (Unusual) Cis Geometry in d(2) Octahedral Dioxo Complexes. A DFT Study

Geometry optimization of the cis and the trans isomers of several octahedral dioxo complexes of d(2) electronic configuration are performed using the gradient-corrected density functional theory (B3LYP and, for some key structures, BP86). With only monodentate sigma donor ligands (ReO(2)(NH(3))(4)(+), 7), the usual energy order is found (i.e., the trans isomer is the most stable). Complexes with a chelating bidentate ligand, OsO(2)(OCH(2)CH(2)O)(NH(3))(2) (10) and ReO(2)(HN=CHCH=NH)(NH(3))(2)(+) (11), are used as models for the experimental complexes 5 and 2 in which the arrangement of the O=M=O unit is trans and cis, respectively. Our calculations actually show an inversion of the relative energy of the two isomers in going from 10 to 11: while the trans isomer is found to be the most stable in 10, the unusual cis diamagnetic isomer is favored by about 29 kcal mol(-)(1) in 11. This result is traced to the geometric and electronic properties of the bidentate ligand, in particular an acute bite angle and good pi acceptor character. In complex 14 with a bipyridine chelating ligand (weaker pi acceptor than diaza-1,4-butadiene in 11), this energy difference is, however, reduced to 7.5 kcal mol(-)(1) (partial geometry optimization).     

Comparative ab initio study of the structures and stabilities of the ethane dication C2H6(2+) and its silicon analogues Si2H6(2+) and CSiH6(2+)

Ab initio MP2/6-311G and QCISD(T)/6-311G levels as well as Gaussian-2 theory were used to perform a comparative study of the structures and stabilities of the ethane dication C(2)H(6)(2+) and its silicon analogues Si(2)H(6)(2+) and CSiH(6)(2+). Similar to previous HF/6-31G results, our present calculations also indicate that the two-electron three-center (2e-3c) bonded carbonium-carbenium structure 1 is more stable than the doubly hydrogen bridged diborane-type structure 2 by about 12 kcal/mol. For the silicon analogue Si(2)H(6)(2+) the calculations, however, indicate that the 2e-3c bonded siliconium-silicenium structure 8 is about 9 kcal/mol less stable than doubly hydrogen bridged structure 9. Similar results were also computed for carbon-silicon mixed CSiH(6)(2+) dication structures. These studies are in agreement with the more electropositive character of silicon compared to carbon. Possible dissociation paths of the minimum structures were also calculated.     

Protonated tert-pentyl dication (C5H122+, isopentane dication)

Structures of the tert-pentyl cation (C(5)H(11)(+)) and its protonated dication (C(5)H(12)(2+), isopentane dication) were studied using ab initio methods at the MP2/cc-pVTZ level. Both C-C and C-H hyperconjugatively stabilized structures 1 and 2 , respectively, were found to be minima on the potential energy surface (PES) of the tert-pentyl cation. Structure 1 was computed to be about as stable as structure 2 (slightly more stable by 0.5 kcal mol(-1)). Inter-conversion between 1 and 2 through transition state 3 has a kinetic barrier of only 1.5 kcal mol(-1). The C-H protonated form (H(3)C)(2)C(+)CH(2)CH(4)(+)4 was found to be the global minimum for the protonated tert-pentyl dication. Charges and (13)C NMR chemical shifts of the dication 4 were calculated and compared to those of monocation 1 to study the effect of the additional charge in the dication.     

On the viability of small endohedral hydrocarbon cage complexes: X@C4H4, X@C8H8, X@C8H14, X@C10H16, X@C12H12, AND X@C16H16

Small hydrocarbon complexes (X@cage) incorporating cage-centered endohedral atoms and ions (X = H(+), H, He, Ne, Ar, Li(0,+), Be(0,+,2+), Na(0,+), Mg(0,+,2+)) have been studied at the B3LYP/6-31G(d) hybrid HF/DFT level of theory. No tetrahedrane (C(4)H(4), T(d)()) endohedral complexes are minima, not even with the very small hydrogen atom or beryllium dication. Cubane (C(8)H(8), O(h)()) and bicyclo[2.2.2]octane (C(8)H(14), D(3)(h)()) minima are limited to encapsulating species smaller than Ne and Na(+). Despite its intermediate size, adamantane (C(10)H(16), T(d)()) can enclose a wide variety of endohedral atoms and ions including H, He, Ne, Li(0,+), Be(0,+,2+), Na(0,+), and Mg(2+). In contrast, the truncated tetrahedrane (C(12)H(12), T(d)()) encapsulates fewer species, while the D(4)(d)() symmetric C(16)H(16) hydrocarbon cage (see Table of Contents graphic) encapsulates all but the larger Be, Mg, and Mg(+) species. The host cages have more compact geometries when metal atoms, rather than cations, are inside. This is due to electron donation from the endohedral metals into C-C bonding and C-H antibonding cage molecular orbitals. The relative stabilities of endohedral minima are evaluated by comparing their energies (E(endo)) to the sum of their isolated components (E(inc) = E(endo) - E(cage) - E(x)) and to their exohedral isomer energies (E(isom) = E(endo) - E(exo)). Although exohedral binding is preferred to endohedral encapsulation without exception (i.e., E(isom) is always exothermic), Be(2+)@C(10)H(16) (T(d)(); -235.5 kcal/mol), Li(+)@C(12)H(12) (T(d)(); 50.2 kcal/mol), Be(2+)@C(12)H(12) (T(d)(); -181.2 kcal/mol), Mg(2+)@C(12)H(12) (T(d)(); -45.0 kcal/mol), Li(+)@C(16)H(16) (D(4)(d)(); 13.3 kcal/mol), Be(+)@C(16)H(16) (C(4)(v)(); 31.8 kcal/mol), Be(2+)@C(16)H(16) (D(4)(d)(); -239.2 kcal/mol), and Mg(2+)@C(16)H(16) (D(4)(d)(); -37.7 kcal/mol) are relatively stable as compared to experimentally known He@C(20)H(20) (I(h)()), which has an E(inc) = 37.9 kcal/mol and E(isom) = -35.4 kcal/mol. Overall, endohedral cage complexes with low parent cage strain energies, large cage internal cavity volumes, and a small, highly charged guest species are the most viable synthetic targets.     

Theoretical study of the C6H3 potential energy surface and rate constants and product branching ratios of the C2H(2Sigma+) + C4H2(1Sigma(g)+) and C4H(2Sigma+) + C2H2(1Sigma(g)+) reactions

Ab initio CCSD(T)cc-pVTZ//B3LYP6-311G(**) and CCSD(T)/complete basis set (CBS) calculations of stationary points on the C(6)H(3) potential energy surface have been performed to investigate the reaction mechanism of C(2)H with diacetylene and C(4)H with acetylene. Totally, 25 different C(6)H(3) isomers and 40 transition states are located and all possible bimolecular decomposition products are also characterized. 1,2,3- and 1,2,4-tridehydrobenzene and H(2)CCCCCCH isomers are found to be the most stable thermodynamically residing 77.2, 75.1, and 75.7 kcal/mol lower in energy than C(2)H + C(4)H(2), respectively, at the CCSD(T)/CBS level of theory. The results show that the most favorable C(2)H + C(4)H(2) entrance channel is C(2)H addition to a terminal carbon of C(4)H(2) producing HCCCHCCCH, 70.2 kcal/mol below the reactants. This adduct loses a hydrogen atom from the nonterminal position to give the HCCCCCCH (triacetylene) product exothermic by 29.7 kcal/mol via an exit barrier of 5.3 kcal/mol. Based on Rice-Ramsperger-Kassel-Marcus calculations under single-collision conditions, triacetylene+H are concluded to be the only reaction products, with more than 98% of them formed directly from HCCCHCCCH. The C(2)H + C(4)H(2) reaction rate constants calculated by employing canonical variational transition state theory are found to be similar to those for the related C(2)H + C(2)H(2) reaction in the order of magnitude of 10(-10) cm(3) molecule(-1) s(-1) for T = 298-63 K, and to show a negative temperature dependence at low T. A general mechanism for the growth of polyyne chains involving C(2)H + H(C[triple bond]C)(n)H --> H(C[triple bond]C)(n+1)H + H reactions has been suggested based on a comparison of the reactions of ethynyl radical with acetylene and diacetylene. The C(4)H + C(2)H(2) reaction is also predicted to readily produce triacetylene + H via barrierless C(4)H addition to acetylene, followed by H elimination.     

Ab initio/GIAO-CCSD(T) study of structures, energies, and 13C NMR chemical shifts of C4H7(+) and C5H9(+) ions: relative stability and dynamic aspects of the cyclopropylcarbinyl vs bicyclobutonium ions

The structures and energies of the carbocations C 4H 7 (+) and C 5H 9 (+) were calculated using the ab initio method. The (13)C NMR chemical shifts of the carbocations were calculated using the GIAO-CCSD(T) method. The pisigma-delocalized bisected cyclopropylcarbinyl cation, 1 and nonclassical bicyclobutonium ion, 2 were found to be the minima for C 4H 7 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level the structure 2 is 0.4 kcal/mol more stable than the structure 1. The (13)C NMR chemical shifts of 1 and 2 were calculated by the GIAO-CCSD(T) method. Based on relative energies and (13)C NMR chemical shift calculations, an equilibrium involving the 1 and 2 in superacid solutions is most likely responsible for the experimentally observed (13)C NMR chemical shifts, with the latter as the predominant equilibrating species. The alpha-methylcyclopropylcarbinyl cation, 4, and nonclassical bicyclobutonium ion, 5, were found to be the minima for C 5H 9 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level ion 5 is 5.9 kcal/mol more stable than the structure 4. The calculated (13)C NMR chemical shifts of 5 agree rather well with the experimental values of C 5H 9 (+).     

High-level ab initio calculations for the four low-lying families of minima of (H2O)20. I. Estimates of MP2/CBS binding energies and comparison with empirical potentials

We report estimates of complete basis set (CBS) limits at the second-order M?ller-Plesset perturbation level of theory (MP2) for the binding energies of the lowest-lying isomers within each of the four major families of minima of (H(2)O)(20). These were obtained by performing MP2 calculations with the family of correlation-consistent basis sets up to quadruple zeta quality, augmented with additional diffuse functions (aug-cc-pVnZ, n=D, T, Q). The MP2/CPS estimates are -200.1 (dodecahedron, 30 hydrogen bonds), -212.6 (fused cubes, 36 hydrogen bonds), -215.0 (face-sharing pentagonal prisms, 35 hydrogen bonds), and -217.9 kcal/mol (edge-sharing pentagonal prisms, 34 hydrogen bonds). The energetic ordering of the various (H(2)O)(20) isomers does not follow monotonically the number of hydrogen bonds as in the case of smaller clusters such as the different isomers of the water hexamer. The dodecahedron lies ca. 18 kcal/mol higher in energy than the most stable edge-sharing pentagonal prism isomer. The TIP4P, ASP-W4, TTM2-R, AMOEBA, and TTM2-F empirical potentials also predict the energetic stabilization of the edge-sharing pentagonal prisms with respect to the dodecahedron, albeit they universally underestimate the cluster binding energies with respect to the MP2/CBS result. Among them, the TTM2-F potential was found to predict the absolute cluster binding energies to within <1% from the corresponding MP2/CBS values, whereas the error for the rest of the potentials considered in this study ranges from 3% to 5%.     

The C4H7+ cation. A theoretical investigation

The potential energy surface of the C4H7+ cation has been investigated with ab initio quantum chemical theory. Extended basis set calculations, including electronic correlation, show that cyclobutyl and cyclopropylcarbinyl cation are equally stable isomers. The saddle point connecting these isomers lies 0.6 kcal/mol above the minima. The global C4H7+ minimum corresponds to the 1-methylallyl cation, which is 9.0 kcal/mol more stable than the cyclobutyl and the cyclopropylcarbinyl cation and 9.5 kcal/mol below the 2-methylallyl cation. These results are in excellent agreement with experimental data.     

Remarkable interplay of redox states and conformational changes in a sterically crowded, cross-conjugated tetrathiafulvalene vinylog

Derivatives of 9-[2-(1,3-dithiol-2-ylidene)ethylidene]thioxanthene have been synthesized using Horner-Wadsworth-Emmons reactions of (1,3-dithiol-2-yl)phosphonate reagents with thioxanthen-9-ylidene-acetaldehyde (5). Further reactions lead to the sterically crowded cross-conjugated "vinylogous tetrathiafulvalene" derivative 9-[2,3-bis-(4,5-dimethyl-1,3-dithiol-2-ylidene)-propylidene]thioxanthene (10). X-ray crystallography, solution electrochemistry, optical spectroscopy, spectroelectrochemistry, and simultaneous electrochemistry and electron paramagnetic resonance spectroscopy, combined with theoretical calculations performed at the B3LYP/6-31G(d) level, elucidate the interplay of the electronic and structural properties in these molecules. For compound 10, multistage redox behavior is observed: the overall electrochemical process can be represented by 10-->10(.+)-->10(2+)-->10(4+) with good reversibility for the 10-->10(.+)-->10(2+) transformations. At the tetracation stage there is the maximum gain in aromaticity at the dithiolium and thioxanthenium rings. Theory predicts that for 10, 10(.+), and 10(2+) the trans isomers are more stable than the cis isomers (by ca. 2-18 kJ mol(-1)), whereas for 10(4+) the cis isomer becomes more stable than the trans isomer (by ca. 25 kJ mol(-1)) [trans and cis refer to the arrangement of the two dithiole moieties with respect to the central ==C(R)--C(H)== fragment]. These data explain the detection in cyclic voltammograms of both trans and cis isomers of 10 and 10(.+) during the reduction of 10(4+) at fast scan rates (>100 mV s(-1)) when the cis-trans isomerization is not completed within the timescale of the experiment. The X-ray structure of the charge-transfer complex (CTC) of 10 with 2,4,5,7-tetranitrofluorene-9-dicyanomethylenefluorene (DTeF) [stoichiometry: 10(.+)(DTeF)(2) (.-)2 PhCl] reveals a twisted conformation of 10(.+) (driven by the bulky thioxanthene moiety) and provides a very rare example of segregated stacking of a fluorene acceptor in a CTC.     

Comparative study of the hypercoordinate ions C7H9+ and C8H9+ by the ab initio/GIAO-CCSD(T) method

A comparative study of the hypercoordinate square-pyramidal carbocations C7H9+ and C8H9+ was performed by the ab initio/GIAO-CCSD(T) method. The structures and 13C NMR chemical shifts of the cations were calculated at the GIAO-CCSD(T)/tzp/dz//MP2/cc-pVTZ level. The bishomo square pyramidal structure 1 was calculated for C7H9+ at the MP2/cc-pVTZ level. The calculated 13C NMR chemical shifts of structure 1 agree extremely well with the experimental values. However, unlike for C7H9+ both the bishomo square pyramidal structure 3 and the trishomocyclopropenium type structure 4 were found to be minima on the potential energy surface of C8H9+. They are very close energetically with cation 3, only 0.7 kcal/mol less stable than cation 4 at the MP2/cc-pVTZ//MP2/cc-pVTZ + ZPE level. Neither structure 3 nor 4 yields NMR spectra that agree with experiment. However, a weighted average of the two reproduces the observed NMR spectrum of C8H9+ (at -80 degrees C) quite well.     

Structures and Stability of HNS_2 Isomers

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｉｔｒｏｇｅｎ ｓｕｌｆｕｒｃｏｎｔａｉｎｉｎｇｃｏｍｐｏｕｎｄｓｈａｖｅａｔ ｔｒａｃｔｅｄｍｕｃｈａｔｔｅｎｔｉｏｎｉｎｔｈｅｆｉｅｌｄｓｏｆｐｏｌｙｍｅｒｉｃｍａｔｅｒｉａｌｓａｎｄａｔｍｏｓｐｈｅｒｅｃｈｅｍｉｓｔｒｙ .1 4 Ａｎｄｓｏｍｅａｖａｉｌａｂｌｅｅｘｐｅｒｉ ｍｅｎｔａｌａｎｄｔｈｅｏｒｅｔｉｃａｌｉｎｆｏｒｍａｔｉｏｎｗａｓｒｅｐｏｒｔｅｄａｂｏｕｔｔｈｅｓｅｃｏｍｐｏｕｎｄｓｗｉｔｈＨ ｅｌｅｍｅｎｔｏｒｏｔｈｅｒｆｕｎｃｔｉｏｎａｌｇｒｏｕｐｓ…     

Ab initio/GIAO-CCSD(T) study of propenoyl (H2C=CH-CO+) and isopentenoyl ((CH3)2C=CH-CO+) cations and their superelectrophilic protonated dications

Structures of superelectrophilic protonated propenoyl (H2C=CH-COH2+) and isopentenoyl ((CH3)2C=CH-COH2+) dications and their parent cations were calculated using ab initio methods at the MP2/6-311+G and MP2/cc-pVTZ levels. Energies were calculated using Gaussian-2 (G2) theory. The alpha-carbon (Calpha) protonated 3 and 7 were found to be the global minima for protonated propenoyl and isopentenoyl dications, respectively. 13C NMR chemical shifts of the cations were also calculated using the GIAO-CCSD(T), GIAO-MP2 and GIAO-SCF methods. 13C NMR chemical shifts of the related tert-butyl cation ((CH3)3C+) and protonated tert-butyl dication ((CH3)2CCH4(2+)) were also computed at the same level to compare and explore the effect of the additional charge in dications.     

Conformational preferences and internal rotation in alkyl- and phenyl-substituted thiourea derivatives

Potential energy surfaces (PES) for rotation about the N-C(sp(3)) or N-C(aryl) bond and energies of stationary points on PES for rotation about the C(sp(2))-N bond are reported for methylthiourea, ethylthiourea, isopropylthiourea, tert-butylthiourea, and phenylurea, using the MP2/aug-cc-pVDZ method. Analysis of alkylthioureas shows that conformations, with alkyl groups cis to the sulfur atom, are more stable (by 0.4-1.5 kcal/mol) than the trans forms. All minima adopt anti configurations with respect to nitrogen pyramidalization, whereas syn configurations are not stationary points on the MP2 potential surface. In contrast, analysis of phenylthiourea reveals that a trans isomer in a syn geometry is the global minimum, whereas a cis isomer in an anti geometry is a local minimum with a relative energy of 2.7 kcal/mol. Rotation about the C(sp(2))-N bond in alkyl and phenyl thioureas is slightly more hindered (9.1-10.2 kcal/mol) than the analogous motion in the unsubstituted molecule (8.6 kcal/mol). The maximum barriers to rotation for the methyl, ethyl, isopropyl, tert-butyl, and phenyl substituents are predicted to be 1.2, 8.9, 8.6, 5.3, and 0.9 kcal/mol, respectively. Corresponding PESs are consistent with the experimental dihedral angle distribution observed in crystal structures. The results of the electronic structure calculations are used to benchmark the performance of the MMFF94 force field. Systematic discrepancies between MMFF94 and MP2 results were improved by modification of selected torsion parameters and one of the van der Waals parameters for sulfur.     

Selection of the cis and trans phosph(III)azane macrocycles [{P(mu-NtBu)}2(1-Y-2-NH-C6H4)]2(Y=O, S)

The 1 : 1 reactions of [ClP(mu-NtBu)]2 with the difunctional aromatic amines 1,2-1-YH-2-NH2-C6H4 in the presence of Et3N give the dimeric phosph(III)azane macrocycles [{P(mu-NtBu)2(1-Y-2-HN-C6H4)]2, predominantly as the cis isomer in the case of Y=O (1.cis) and as the trans isomer for Y=S (2.trans). Model M.O. calculations suggest that the selection of the cis and trans isomers is not thermodynamically controlled. The alternative isomers 1.trans and 2.cis are generated exclusively by the deprotonation of the model intermediates [(1-Y-2-NH2-C6H4)P(mu-NtBu)]2[Y=O (3), S (4)] with nBuLi followed by cyclisation with [ClP(mu-NtBu)]2. The solid-state structures of 1.cis/trans(50 : 50), 2.cis, 3 and 4 are reported.     

Molecular structure of the trans and cis isomers of metal-free phthalocyanine studied by gas-phase electron diffraction and high-level quantum chemical calculations: NH tautomerization and calculated vibrational frequencies

The molecular structure of the trans isomer of metal-free phthalocyanine (H2Pc) is determined using the gas electron diffraction (GED) method and high-level quantum chemical calculations. B3LYP calculations employing the basis sets 6-31G**, 6-311++G**, and cc-pVTZ give two tautomeric isomers for the inner H atoms, a trans isomer having D2h symmetry and a cis isomer having C2v symmetry. The trans isomer is calculated to be 41.6 (B3LYP/6-311++G**, zero-point corrected) and 37.3 kJ/mol (B3LYP/cc-pVTZ, not zero-point corrected) more stable than the cis isomer. However, Hartree-Fock (HF) calculations using different basis sets predict that cis is preferred and that trans does not exist as a stable form of the molecule. The equilibrium composition in the gas phase at 471 degrees C (the temperature of the GED experiment) calculated at the B3LYP/6-311++G** level is 99.8% trans and 0.2% cis. This is in very good agreement with the GED data, which indicate that the mole fraction of the cis isomer is close to zero. The transition states for two mechanisms of the NH tautomerization have been characterized. A concerted mechanism where the two H atoms move simultaneously yields a transition state of D2h symmetry and an energy barrier of 95.8 kJ/mol. A two-step mechanism where a trans isomer is converted to a cis isomer, which is converted into another trans isomer, proceeds via two transition states of C(s) symmetry and an energy barrier of 64.2 kJ/mol according to the B3LYP/6-311++G** calculation. The molecular geometry determined from GED is in very good agreement with the geometry obtained from the quantum chemical calculations. Vibrational frequencies, IR, and Raman intensities have been calculated using B3LYP/6-311++G**. These calculations indicate that the molecule is rather flexible with six vibrational frequencies in the range of 20-84 cm(-1) for the trans isomer. The cis isomer might be detected by infrared matrix spectroscopy since the N-H stretching frequencies are very different for the two isomers.     

Theoretical study of M(PH3)2 complexes of C60, corannulene (C20H10), and sumanene (C21H12) (M = Pd or Pt). Unexpectedly large binding energy of M(PH3)2(C60)

DFT and MP2 to MP4(SDQ) methods were applied to M(PH3)2(C60), Pt(PH3)2(C20H10), and Pt(PH3)2(C21H12) (M = Pd or Pt, C20H10 = corannulene, and C21H12 = sumanene). The binding energy considerably fluctuates around MP2 and MP3 levels but much less upon going from MP3 to MP4(SDQ) in Pt(PH3)2(C2H4), Pt(PH3)2(C20H10), and Pt(PH3)2(C21H12). Also, the MP4(SDQ) method presents a binding energy similar to that of the CCSD(T) method in Pt(PH3)2(C2H4). Thus, it is likely that the MP4(SDQ) method is useful to evaluate binding energies of these complexes. The binding energies of Pt(PH3)2(C20H10) and Pt(PH3)2(C21H12) are evaluated to be 24.9 and 26.1 kcal/mol, respectively, by the MP4(SDQ) method and only +5.8 and -2.6 kcal/mol, respectively, by the DFT(B3LYP) method. These MP4(SDQ)-calculated binding energies of Pt(PH3)2(C20H10) and Pt(PH3)2(C21H12) are similar to that of Pt(PH3)2(C2H4), which strongly suggests that these complexes can be successfully synthesized. The binding energy of Pt(PH3)2(C60) is evaluated to be 44.8 and 45.5 kcal/mol with the ONIOM(MP4(SDQ):UFF) and ONIOM(MP4(SDQ):B3LYP) methods, respectively, and that of the Pd analogue is evaluated to be 39.9 kcal/mol with the ONIOM(MP4(SDQ):UFF) method, whereas the DFT(B3LYP), DFT(BVP86), and DFT(BPW91) methods provide much smaller binding energies. It is noted that these binding energies are much larger than those of the ethylene, corannulene, and sumanene analogues. This difference is reasonably interpreted in terms that the LUMO of C60 is at much lower energy than those of ethylene, corannulene, and sumanene. We investigated also how to separate the high level and the low level regions in the ONIOM calculation of M(PH3)2(C60) and proposed here the reasonable way to evaluate the binding energy of transition-metal complexes of C60.