[Si,N,N,P]体系异构体结构与稳定性

在B3LYP/6-311G(d)和QCISD(t)/6-311+G(2df)(单点)水平下计算得到了[Si,N,N,P]体系势能面上的一些驻点, 研究了异构体之间的异构化过程, 并讨论了该体系异构体的结构与稳定性. 结果表明在[Si,N,N,P]体系中的11个异构体中只有含有Si-P交叉成键的类蝶型四元环异构体SiNPN(E1, ²A″)是稳定的, 其他异构体由于易于重排或是解离, 稳定性较低. 研究同时给出了E1的电子与几何结构、振动频率及指认、偶极矩和转动常数等相关数据. 分析结果表明通过中间体E3(SiNPN)的反应SiN(²Π) + PN(¹Σ)→E1是由SiN(²Π)和PN(¹Σ)碎片产生E1的最有利途径, 因此E1可成为星际观测的目标. G2级别下此反应的反应焓变为215.25 kJ/mol (298.15 K), 而异构体E1在298.15 K下的标准生成焓为457.99 kJ/mol.     

[C, N, N, P]体系异构体结构与稳定性

在B3LYP/6-311G(d)和QCISD(t)/6-311+G(2df)水平下计算得到了[C, N, N, P]体系的9个异构体、20个过渡态及一些相关的解离碎片, 并研究了异构体之间的异构化过程, 讨论了该体系异构体的结构与稳定性. 结果表明在得到的9个异构体中, 具有2A′电子态的、非线型链状NCNP, NCPN, CNPN和CNNP 4个异构体是较稳定的. 其余异构体由于易于重排或是解离, 稳定性较低. 研究同时给出了这4个异构体的振动频率、转动常数、偶极矩、第一绝热电离能和绝热电子亲合能等相关数据, 为可能的实验和星际光谱观测提供相关的理论数据. 得到的计算结果与等价电子体系SiNNP进行了比较.     

HBO2异构体的结构和相对稳定性

用ab initio方法在MP2/6 311＋＋G(d,p)水平下优化得到了HBO2体系的若干异构体和过渡态，并在QCISD(t)/6 311＋＋G(3df,2p)//MP2/6 311＋＋G(d,p)水平下进行了单点能量校正.对计算结果的分析表明，无论是在热力学还是在动力学上，具有链状结构的HOBO异构体(E1)是势能面上最稳定的结构，并对E1的电子结构进行了分析；另一具有C2v对称性的HBO(O)结构的异构体(E2)的能量比E1高381.72 kJ•mol－1，由于E2处于一个较深的势垒中，因此是比较稳定的，可以推断，在适合的实验中应该可以观察到异构体E2.     

N8H8链状异构体结构与稳定性的理论研究

采用密度泛函理论(DFT), 在B3LYP/6-311＋＋G(d,p)基组上计算得到了21种N₈H₈链状异构体, 并研究了这些异构体间可能的互变异构情况. 为了得到更为精确的能量信息, 计算了QCISD(T)/6-311G(d,p)基组水平上各物质的能量. 所得的21种异构体分为4类(4种类型链状化合物): A为直链, B有一个支链, C有2个支链, D有3个支链; D类只有一种, A类稳定构型2种, B类稳定构型12种, C类稳定构型6种; 相对稳定的分别为: B2-1构型, B2-3构型和C23-2构型. 我们研究发现N₈H₈链状异构体中含有明显N＝N双键特征有利于化合物稳定性的提高.     

HPOS体系异构体结构与稳定性

在MP2／6－311＋＋G(d,p)和QCISD（T）／6－311＋＋G（3df,2p)（单点）水 平下计算得到了HPOS体系势能面上18个异构体和25个过渡态及解离碎片等驻点，并 分析了这些异构体的结构及异构化过程，讨论了可能的解离方式。在得到的异构体 中，有8个异构体是动力学较稳定的，它们是dis-HOPS,trans-HOPS,trans-HSPO, cis-HSPO,HP(O)S(Cs),trans-HPSO,cis-HPSO和HP（O）S（C1）。这些异构体在实 验中应该可以观测到。理论研究表明，P与S原子较强的超价能力在降低异构体能量 ，提高异构体动力学稳定性方面起到了关键的作用。得到的计算结果与HPO2， HPS2，HNOS等价电子相同的体系进行了比较。     

HAsO~2异构体结构、相对稳定性与体系势能面

在MP2/6-311++G(d,p)和QCISD(T)/6-311++G(3df,2p)(单点)水平下计算得到了包括9个异构体和10个过滤态的HAsO~2体系势能面。在势能面上，异构体cis-HOAsO(E1)的能量是最低的，其次是trans-HOAsO(E2)和HAsO(O)(C~2~V,E3)，能量分别比cis-HOAsO高13.15和192.74kJ/mol。根据体系的势能面，异构体E1，E3及cis-HOOAs(E6),trans-HOOAs(E5)具有一定的动力学稳定性，在实验中应该可以观测到。AsH和O~2反应的第一步产物将会异构化为具有较高动力学稳定性的异构体E3；而OH和AsO反应可直接生成E1。计算结果与HPO~2,HPS~2,HNO~2,HNS~2等价电子相同的分子的势能面进行了比较。     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

质子化甜菜碱阳离子与双三氟甲磺酰亚胺阴离子相互作用的理论研究

利用密度泛函理论(DFT)方法, 在B3LYP/6-311＋G(d,p)水平下, 对质子化甜菜碱阳离子与双三氟甲磺酰亚胺阴离子形成的气态阴阳离子对([Hbet][Tf₂N])进行理论研究, 通过几何结构优化和频率分析得到势能面上的六个稳定构型. 计算结果表明[Hbet]^＋和[Tf₂N]^－之间能够形成较强的氢键相互作用, 氢键相互作用的稳定化能主要来源于[Tf₂N]^－中O, N原子的孤对电子和[Hbet]^＋中参与形成氢键的O-H, C-H反键轨道之间的相互作用. 自然布局分析(NPA)给出气态[Hbet][Tf₂N]离子对中阴阳离子间的电荷转移比传统离子液体中电荷转移的数值小. AIM (atoms in molecules)分析得到[Hbet]^＋和[Tf₂N]^－之间的氢键相互作用以静电作用为主. 通过实验和理论结果相比较, 初步探讨影响离子液体熔点及其对金属离子选择性溶解的结构因素.     

新型取代基炔类分子的合成及其与B(C6F5)3的反应

合成了5种不同取代基的炔类化合物Mes2HSiC≡CPh(1,Mes=2,4,6-Me3C6H2)、[tBuC(NAr)2]GeC≡CPh(2,Ar=2,6-iPr2C6H3)、[PhC(NtBu)2]SnC≡CPPh2(3)、[HC(CMe)2(NAr)2]Sn C≡CPPh2(4)和[HC(CMe)2(NAr)2]ZnC≡CPPh2(5),研究了这些化合物与B(C6F5)3的反应.在与B(C6F5)3的反应中,1和2均发生1,1-碳硼化反应生成烯烃化合物(Ph)(Mes2HSi)C=C(C6F5)B(C6F5)2 (6)和{[tBuC(NAr)2]Ge}(Ph)C=C(C6F5)B(C6F5)2 (7), 7是一种GeⅡ/B松散Lewis酸碱对化合物;3~5则都发生B(C6F5)3与配体金属基的位置交换、进而配体金属基转换键合PPh2的反应,分别生成新颖的分子内双性离子炔烃化合物[PhC(NtBu)2]SnP(Ph2)C≡CB(C6F5)3 (8)、[HC(CMe)2(NAr)2]SnP(Ph2)C≡CB(C6F5)3(9)、[HC(CMe)2(NAr)2]ZnP(Ph2)C≡CB(C6F5)3 (10).文中还讨论了反应机理.     

HPS2异构体结构与相对稳定性

在MP2/6-311++G(d,p)和QCISD(t)/6-311++G(3df,2p)(单点)水平下计算得到了包括9个异构体和15个过渡态的HPS2势能面.其中,异构体trans-HSPS(E1)的能量最低, 其次是cis-HSPS(E2)和HPS(S)(C2_V,E3),能量分别比trans-HSPS(E1)高3.43和14.17 kJ/mol.计算结果表明,异构体E1,E2,E3和立体的三元环结构的异构体HP(S)S(C_S,E4), 及与E4共存的trans-HPSS(E5)和cis-HPSS(E6)具有一定的稳定性.在QCISD(t)/6-311++G(3df,2p)//MP2/6-311++G(d,p)并包含零点能水平下,PH和S₂反应生成的E6和E5分别越过65.75和71.73 kJ/mol的势垒就可以异构化为具有较高动力学稳定性的产物E4, 计算结果对实验认定的产物是cis-HPSS(E6)的结论进行了修正.     

Homoaromaticity in carbene intermediates

Analyses of isodesmic reactions for singlet 7-carbenanorbornene (4S), 8-carbena-endo-tricyclo[3.2.1.0(2,4)]octane (5S), 3-carbenabicyclo[3.1.0]hexane (3S), 2-carbenanorbornene 6S, and 2-carbenabicyclooctadiene 7S at the B3LYP/6-311+G(3df, 2p)//B3LYP/6-31G level provide stabilization energies of 13.83, 13. 50, 3.00, -2.22, and -3.01 kcal/mol, respectively. The C7 carbene in 4S and the C8 carbene center in 5S are strongly bent toward the double bond and cyclopropane ring, respectively, in contrast to their related triplets, 4T and 5T and parent hydrocarbons. The geometric change for 3S compared to 3T or to parent bicyclo[3.1. 0]hexane is minimal. Comparison of the stability of 6S with 2-carbenanorbornane and the geometry of 6S with that of 6T and also with the singlet and triplet 2-carbenanorbornane suggests very modest bridging. The stabilization energy and geometry of 2-carbenabicyclooctadiene 7 resemble antihomoaromatic bicyclooctadienyl cation 9 rather than the related homoaromatic bicyclooctadienyl anion 8. The diamagnetic susceptibility exhaltations (Deltachi) of 3, 4, 5, and 7, calculated at the B3LYP/6-311+G(2d,p)//B3LYP/6-31G level, are -0.7, 22.7, 26.0, and -10.3 cgs-ppm, respectively. The singlet-triplet energy differences, DeltaE(TS), for carbenacyclohexane, carbenacyclopentane, 3, 4, 5, 6, 7, and 2-carbenabicyclo[3.2.1]oct-3-ene at the B3LYP/6-311+G(3df, 2p)//B3LYP/6-31G level are 3.4, 10.3, 8.8, 27.1, 25.9, 12.7, -4.0, and -0.9 kcal/mol, respectively.     

Theoretical study on structures and stabilities of [H,Ge,C,N

Theoretical investigations are performed for the first time on the simplest hydrogenated germanium cyanide [H,Ge,C,N], whose analogs [H,C(2),N] and [H,Si,C,N] have been detected in space and laboratory, respectively. The detailed potential energy surfaces in both singlet and triplet states are constructed at the CCSD(T)/6-311+G(3df,2p)//B3LYP/6-31G(d)+ZPVE level, including 18 minimum isomers and 26 interconversion transition states. The former three low-lying and kinetically stabilized isomers are HGeCN (1)1 (0.0 kcal/mol), HGeNC (1)2 (5.1 kcal/mol), and cyclic cCHNGe(1)7 (11.1 kcal/mol). In addition, five isomers HCNGe (1)3 (33.8), HNCGe (1)5 (29.8), cNHCGe (1)8 (37.9), HGeCN (3)1 (30.1), and HNCGe (3)5 (26.5) each have considerable barriers, despite their high energies. Future laboratory characterization and astrophysical detection of the eight [H,Ge,C,N] isomers, especially the former three low-lying species (1)1, (1)2, and (1)7, are highly recommended. The accurate spectroscopic data at the QCISD/6-311G(d,p) level are provided. For some species, the CBS-QB3 calculations are also performed. Wherever possible, comparisons with the analogous [H,C(2),N] and [H,Si,C,N] are made on the structural, energetic, and bonding properties.     

Theoretical study on germanium cyanide radical GeCN and its ions

A detailed theoretical study is performed on the hitherto unknown germanium cyanide radical and its ions. The (2)Pi state GeCN lies 5.0 kcal/mol lower than the (2)Pi state GeNC at the coupled-cluster theory including single and double excitations and perturbative inclusion of triple excitations [CCSD(T)]/6-311++G(3df)//quadratic configuration interaction with single and double excitations (QCISD)/6-311G(d)+zero-point vibrational energy (ZPVE) level. For interconversion between them, two electronic state pathways (2)A(') and (2)A(") are located, with the latter being 0.7 kcal/mol more favorable than the former. On the (2)A(") path, the GeCN-->GeNC and GeNC-->GeCN conversion barriers are 14.5 and 9.5 kcal/mol, respectively. The detailed singlet and triplet potential-energy surfaces of both the cationic and anionic GeCN species are also investigated. On the ground-state electronic hypersurface, singlet GeNC(+) is 4.6 kcal/mol more stable than singlet GeCN(+), whereas triplet GeNC(-) is 10.0 kcal/mol less stable than triplet GeCN(-). The relative energy difference between the GeCN(0,+/-) and GeNC(0,+/-) can be well correlated with the number of vacant orbitals on the Ge atom. The stability of the neutral and ionic CGeN and cyclic cGeCN is also discussed. The predicted structures, spectroscopies, ionization, and affinity energies as well as the Renner-Teller properties are expected to provide reliable estimates for future characterization of the potential GeCN and GeNC radicals as well as their ionic counterparts both in the laboratory and in the interstellar space.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Role of singlet diradicals in reactions of 2-carbenabicyclo[3.2.1]octa-3,6-diene

The generation of 2-carbenabicyclo[3.2.1]octa-3,6-diene (1) results in the formation of C(8)H(8) hydrocarbons endo-6-ethynylbicyclo[3.1.0]hex-2-ene (4), semibullvalene (5), and 5-ethynyl-1,3-cyclohexadiene (6), and C(8)H(10) hydrocarbons bicyclo[3.2.1]octa-2,6-diene (7), tricyclo[3.2.1.0(4,6)]oct-2-ene (8), and tetracyclo[3.3.0.0(2,8)0(4,6)]octane (9). Focus is placed on three mechanistic pathways for the formation of the C(8)H(10) hydrocarbon fraction: (a) abstraction of hydrogen by triplet carbene 1T to produce an equilibrating set of monoradicals, (b) interconversion of triplet carbene 1T into tricyclic triplet diradical 19T and tetracyclic triplet diradical 20T, and (c) interconversion of singlet 1S with analogous singlet diradical 19S and 20S. Ab initio calculations at the (U)B3LYP/6-311+G(3df,2p)//(U)B3LYP/6-31G(d,p) and broken spin symmetry UBS B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d,p) levels rule out choices (a) and (b) and are consistent with the singlet diradical process.     

Molecular structure of the octatetranyl anion, C8H-: a computational study

The equilibrium molecular structure of the octatetranyl anion, C8H(-), which has been recently detected in two astronomical environments, is investigated with the aid of both ab initio post-Hartree-Fock and density functional theory (DFT) calculations. The model chemistry adopted in this study was selected after a series of benchmark calculations performed on molecular acetylene for which accurate gas-phase structural data are available. Geometry optimizations performed at the CCSD/6-311+G(2d,p), QCISD/6-311+G(2d,p), and MP4(SDQ)/6-311+G(2d,p) levels of theory yield for C8H(-) an interesting polyyne-type structure that defies the chemical formula displaying a simple alternation of triple and single carbon-carbon bonds, [:C[triple bond]C-C[triple bond]C-C[triple bond]C-C[triple bond]CH](1-). In the optimized geometry of C8H(-), as one proceeds from the naked carbon atom on one side of the chain to the CH unit on the opposite side of the chain, the short (formally triple) carbon-carbon bonds decrease in length from 1.255 to 1.213 A whereas the long (formally single) carbon-carbon bonds increase (albeit only slightly) in length from 1.362 to 1.378 A (CCSD results). In striking contrast, both MP2 and DFT (B3LYP and PBE0) calculations fail in reproducing the pattern of the carbon-carbon bond lengths obtained with the CCSD, QCISD, and MP4 methods. The structures of three shorter n-even chains, C(n)H(-) (n = 2, 4, and 6), along with those of four n-odd compounds (n = 3, 5, 7, and 9) are also investigated at the CCSD/6-311+G(2d,p) level of theory.     

Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

Neutral homoaromaticity in some heterocyclic systems

[reaction: see text] Neutral homoaromaticity has been evaluated in heterocyclic systems related to the bicyclo[3.2.1]octane skeleton with replacement of CH(2) at C-2 in bicyclo[3.2.1]octa-3,6-diene with X = BH, AlH, Be, Mg, O, S, PH, NH (12); replacement of CH at C-3 in bicyclo[3.2.1]octa-3,6-dien-2-yl anion with PH, S, NH, O (13); and replacement at C-2 and C-3 with N and O (14). Stabilization energies (SE) are evaluated using density functional theory and homodesmotic equations at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level for series 12-14. Stabilization energies are compared with diamagnetic susceptibility exaltations, Lambda, CSGT-B3LYP/6-31G(d)//B3LYP/6-31G(d), and nucleus-independent chemical shifts (NICS), GIAO-B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d). Analysis of frontier orbitals and geometries, B3LYP/6-31G(d)//B3LYP/6-31G(d), and proton affinities of 2-azabicyclo[3.2.1]octa-3,6-diene, pyrrole, and related model systems, B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d), provide complementary evidence supporting the division of the substrates evaluated into antihomoaromatic (12, X = BH, AlH, and Be), nonhomoaromatic (12, X = O, S, NH, PH), and homoaromatic (13, X = S, PH, NH, O and 14 X = ON), with 12 (X = Mg) appearing as transitional between anti- and nonhomoaromatic.     

Quadricyclane radical cation rearrangements: a computational study of the transformations to 1,3,5-cycloheptatriene and norbornadiene

An alternative skeletal rearrangement of the quadricyclane radical cation (Q*+) explains the side products formed in the one-electron oxidation to norbornadiene. First, the bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cation, with an activation energy of 14.9 kcal mol(-1), is formed. Second, this species can further rearrange to 1,3,5-cycloheptatriene through two plausible paths, that is, a multistep mechanism with two shallow intermediates and a stepwise path in which the bicyclo[3.2.0]hepta-2,6-diene radical cation is an intermediate. The multistep rearrangement has a rate-limiting step with an estimated activation energy of 16.5 kcal mol(-1), which is 2.8 kcal mol(-1) lower in energy than the stepwise mechanism. However, the lowest activation energy is found for the Q*+ cycloreversion to norbornadiene that has a transition structure, in close correspondence with earlier studies, and an activation energy of 10.1 kcal mol(-1), which agrees well with the experimental estimate of 9.3 kcal mol(-1). The computational estimates of activation energies were done using the CCSD(T)/6-311+G(d,p) method with geometries optimized on the B3LYP/6-311+G(d,p) level, combined with B3LYP/6-311+G(d,p) frequencies.     

Absolute thermodynamic measurements of alkali metal cation interactions with a simple dipeptide and tripeptide

Absolute bond dissociation energies (BDEs) of glycylglycine (GG) and glycylglycylglycine (GGG) to sodium and potassium cations and sequential bond energies of glycine (G) in Na+G2 were determined experimentally by threshold collision-induced dissociation (TCID) in a guided ion beam tandem mass spectrometer. Experimental results showed that the binding energies follow the order of Na+ > K+ and M+GGG > M+GG > M+G. Theoretical calculations at the B3LYP/6-311+G(d) level show that all complexes had charge-solvated structures (nonzwitterionic) with either [CO,CO] bidentate or [N,CO,CO] tridentate coordination for M+GG complexes, [CO,CO,CO] tridentate or [N,CO,CO,CO] tetradentate coordination for M+GGG complexes, and [N,CO,N,CO] tetradentate coordination for Na+G2. Ab initio calculations at three different levels of theory (B3LYP, B3P86, and MP2(full) using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d) level) show good agreement with the experimental bond energies. This study demonstrates for the first time that TCID measurements of absolute BDEs can be successfully extended to biological molecules as complex as a tripeptide.