3,4-二硒方酸芳香性和气相酸性的理论研究

在ＲＨＦ／６－３１１Ｇ水平优化得到３,４－二硒方酸（３,４－二氢硒基－３－环丁烯－１,２－二酮）３ 种平面构象异构体的平衡几何构型。进一步用ＭＰ２（ｆｕｌｌ）／６－３１１Ｇ／／ＲＨＦ／／６－３１１Ｇ方法计算单点能量,发现ＺＺ型异构体是能量最低构象,且ＺＺ和ＺＥ型能量非常接近。用优化的最稳定构象ＺＺ型异构体在ＲＨＦ／６－３１１Ｇ／／ＲＨＦ／６－３１１Ｇ、ＲＨＦ／６－３１１＋ Ｇ／／ＲＨＦ／６－３１１＋ Ｇ、ＭＰ２（ｆｕｌｌ）／６－３１１＋ Ｇ／／ＲＨＦ／６－３１１＋ Ｇ 和Ｂ３ＬＹＰ／６－３１１＋ Ｇ／／Ｂ３ＬＹＰ／６－３１１＋ Ｇ水平计算其气相酸性（ΔＧ°）和同键反应芳香性稳定化能（ＨＡＳＥ）。用基团加和法（ｇｒｏｕｐ ｉｎｃｒｅｍ ｅｎｔ ａｐ－ｐｒｏａｃｈ ） 在 ＲＨＦ／６－３１１ ＋ Ｇ／／ＲＨＦ／６－３１１ ＋ Ｇ 和 Ｂ３ＬＹＰ／６－３１１ ＋ Ｇ／／Ｂ３ＬＹＰ／６－３１１＋ Ｇ水平计算其磁化率增量（Λ）。计算结果指出标题化合物的同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

3,4-二硒方酸芳香性和气相酸性的理论研究

在ＲＨＦ／６－３１１Ｇ水平优化得到３,４－二硒方酸（３,４－二氢硒基－３－环丁烯－１,２－二酮）３ 种平面构象异构体的平衡几何构型。进一步用ＭＰ２（ｆｕｌｌ）／６－３１１Ｇ／／ＲＨＦ／／６－３１１Ｇ方法计算单点能量,发现ＺＺ型异构体是能量最低构象,且ＺＺ和ＺＥ型能量非常接近。用优化的最稳定构象ＺＺ型异构体在ＲＨＦ／６－３１１Ｇ／／ＲＨＦ／６－３１１Ｇ、ＲＨＦ／６－３１１＋ Ｇ／／ＲＨＦ／６－３１１＋ Ｇ、ＭＰ２（ｆｕｌｌ）／６－３１１＋ Ｇ／／ＲＨＦ／６－３１１＋ Ｇ 和Ｂ３ＬＹＰ／６－３１１＋ Ｇ／／Ｂ３ＬＹＰ／６－３１１＋ Ｇ水平计算其气相酸性（ΔＧ°）和同键反应芳香性稳定化能（ＨＡＳＥ）。用基团加和法（ｇｒｏｕｐ ｉｎｃｒｅｍ ｅｎｔ ａｐ－ｐｒｏａｃｈ ） 在 ＲＨＦ／６－３１１ ＋ Ｇ／／ＲＨＦ／６－３１１ ＋ Ｇ 和 Ｂ３ＬＹＰ／６－３１１ ＋ Ｇ／／Ｂ３ＬＹＰ／６－３１１＋ Ｇ水平计算其磁化率增量（Λ）。计算结果指出标题化合物的同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

1,2-二硫方酸的气相酸性和芳香性

在ＲＨＦ／６－３１１Ｇ、ＲＨＦ／６－３１１＋Ｇ　和Ｂ３ＬＹＰ／６－３１１＋Ｇ　水平,优化得到１,２－二硫方酸（３,４．二羟基－３－环丁烯－１,２－二硫酮）三种平面构象异构体的平衡几何构型．进一步用ＭＰＺ（ｆｕｌｌ）／６－３１１＋Ｇ　、／／ＲＨＦ／６－３１１＋Ｇ　方法计算三种异构体的单点能量,发现ＺＺ型异构体是能量最低构象,且ＺＺ和ＺＥ型能量非常接近．用优化的最稳定构象ＺＺ型异构体,在ＲＨＦ／６－３１１Ｇ　／／ＲＨＦ／６－３１１Ｇ　、ＲＨＦ／６．３１１＋Ｇ　／／ＲＨＦ／６－３１１＋Ｇ　、ＭＰＺ（ｆｕｌｌ）／６．３１１＋Ｇ　／／ＲＨＦ／６－３１１＋Ｇ　和Ｂ３ＬＹＰ／６－３１１＋Ｇ　斤Ｂ３ＬＹＰ／６－３１１＋Ｇ　水平,计算其气相酸性（△Ｇ＿２９８）,并用同键反应方法在同样水平计算其芳香性稳定化能．用基团加和法（ｇｒｏｕｐｉｎｃｒｅｍｅｎｔａｐｐｒｏａｃｈ）在ＲＨＦ／６．３１１＋Ｇ　／／ＲＨＦ／６－３１１＋Ｇ　和Ｂ３ＬＹＰ／６－３１１＋Ｇ　／／Ｂ３ＬＴｈ／６－３１１＋Ｇ　水平计算其磁化率增量（　）．计算结果指出,标题化合物的键长发生了平均化,芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香     

1, 2-二硒方酸气相酸性和芳香性的理论研究

在RHF/6-311G^*^*水平优化得到1,2-二硒方酸(3,4-二羟基-3-环丁烯-1,2-二硒酮)三种平面构象异构体的平衡几何构型。进一步用MP2(full)/6-311G^*//RHF/6-311G^*^*方法计算三种异构体的单点能量,发现ZZ型异构体是能量最低构象,且ZZ和ZE型能量非常接近。用优化的最稳定构象ZZ型异构体在RHF/6-311G^*^*//RHF/6-311G^*^*,RHF/6-311+G^*^*//RHF/6-311+G^*^*,MP2(full)/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其气相酸性[ΔGⅲ~(~2~9~8~K~)]和同键反应芳香性稳定化能(HASE)。用基团加和法(groupincrementapproach)在RHF/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其磁化率增量(Λ)。计算结果指出标题化合物的键长发生了平均化,同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

3,4—二硫方酸的从头算研究

在６－３１Ｇ°水平上对３，４－二硫方酸（３，４－二巯基－３－环丁烯－１，２－二酮）的３种平面构象异构体进行ＳＣＦ计算，结果表明，ＺＺ型异构体最稳定，ＸＥ型次之，从等键反应能量分析，３，４－二太酸的稳定性，与苯作探讨其芳香性，关在６－３１Ｇ水平上计算了３种构象的振动频率。     

1,2-二硒-3，4-二硫方酸的从头计算

在RHF／STO－3G、STO－3G和3－21G（仅对ZE型）水平上，用abinitioSCF方法优化得到1，2－二硒－3，4－二硫方酸（3，4－二巯基－3－环丁烯－1，2－二硒酮）三种平面异构体的平衡构型，进一步用3－21G／／STO－3G（对ZZ和EE型）方法计算总能量，发现三种平面导构体中ZZ型是最稳定构象，ZE型次之，并与1，2－二硒方酸和3，4－二硫方酸从头算结果作了比较，用abinitio数值方法在STO－3G水平上计算了三种异构体的谐振动频率．     

3，4－二硫方酸的从头算研究

在６－３１Ｇ水平上对３，４－二硫方酸（３，４－二巯基－３－环丁烯－１，２－二酮）的３种平面构象异构体进行ＳＣＦ计算．结果表明，ＺＺ型异构体最稳定，ＺＥ型次之，从等键反应能量分析３，４－二硫方酸的稳定性，与苯作比较探讨其芳香性．并在６－３１Ｇ水平上计算了３种构象的振动频率．     

1,2-二硒方酸的从头计算

在ＲＨＦ／ＳＴＯ－３Ｇ和ＳＴＯ－３Ｇ＊水平上用ａｂｉｎｉｔｉｏＳＣＦ方法优化得到１，２－二硒方酸（３，４－二羟基－３－环丁烯－１，２－二硒酮）三种平面异构体的平衡几何构型，发现三种平面异构体中ＺＺ型是能量最低构象。用ａｂｉｎｉｔｉｏ数值方法在ＲＨＦ／ＳＴＯ－３Ｇ＊水平上计算了三种平面异构体的谐振动频率。     

二苯二硫和二苄二硫润滑性能的量子化学研究

The molecular geometries optimization and electronic structures of diphenyl disulfide (DPDS) and dibenzyl disulfide (DBDS) compounds were investigated by density functional theory (DFT) and ab initio method at the 6-31G basis set level. The active atoms and bonds of reaction were provided by frontier molecular orbital theory. The molecular orbital parameters of DPDS and DBDS compounds and iron atom cluster were calculated by using density functional theory. The interaction pattern between the organic disulfide compounds and iron atom cluster was discussed based on the approximate rule of orbital energy. Some parameters characterizing the action strength between the organic disulfide compounds and iron atom cluster, including the bonding strength, reactive strength and static action strength, were analyzed by using frontier electron density, super de-localizability, net atomic charge and the interaction energy of chemical adsorption as criteria. The results indicate that S-S chemical bond and C-S chemical bond of the compounds are inclined to be broken when DPDS and DBDS interact with the metal. The anti wear ability order of DPDS and DBDS compounds is DPDS>DBDS, and the extreme pressure ability order of DPDS and DBDS compounds is DBDS>DPDS, and the prediction results based on quantum chemistry calculations are in good accordance with the friction and wear test results.     

2-溴丙酸气相热消除反应的机理

用密度泛函（DFT）方法,在B3LYP/6-31G**水平上对2-溴丙酸气相消除反应机理进行了研究.计算表明,反应主要是通过半极化五元环结构过渡态进行的,羧基上的氢原子协助溴原子离去,羧基氧原子帮助稳定过渡态.在B3LYP/6-311＋＋G(3df,3pd)水平上对B3LYP/6-31G**优化的几何构型进行了单点能计算,计算所得反应的速度控制步骤的活化能为189.461 kJ•mol－1,偏离实验值((180.3±3.4) kJ•mol－1)5.08％.     

2-溴乙酸气相热消除反应的机理探讨

分别利用密度泛函、从头算分子轨道方法及半经验分子轨道AM 1法对 2 溴乙酸气相热消除反应的机理进行了理论研究 .结果表明 ,虽然用三种不同的方法得到的结果有一定差异 (如反应势能面上各驻点的几何构型、反应的活化能 ) ,但研究得到的反应机理是一致的 ,即反应分两步进行 .其中 ,具有五元环过渡态的第一步反应为速控步骤 ,与实验工作者所建议的机理相符 .     

NiTi过渡金属团簇的结构和芳香性的理论研究

在B3PW91/6-311+G(d)计算水平上, 计算并讨论了Ni₄Ti₂, [Ni₄Ti₂]²⁺, [Ni₄Ti₂]^2-与Ni₄Ti₄, [Ni₄Ti₄]²⁺, [Ni₄Ti₄]^2-团簇的几何结构和芳香性. 在构型优化过程中得到了Ni₄Ti₂(D_4h), [Ni₄Ti₂]²⁺(D_4h), [Ni₄Ti₂]^2-(D_4h)和Ni₄Ti₄(D_2h)4个稳定构型, 发现当引入上下2个Ti原子后, Ni₄环成为了平面正方形构型. 核无关化学位移(NICS)计算结果表明, Ni₄Ti₂(D_4h)与Ni₄Ti₄(D_2h)的NICS值为正, 而[Ni₄Ti₂]²⁺(D_4h)和[Ni₄Ti₂]^2-(D_4h)的NICS值为负, 且[Ni₄Ti₂]^2-(D_4h)的NICS值更负. 同时还发现, 由s与d轨道参与形成的反磁性环流是引起[Ni₄Ti₂]²⁺(D_4h)和[Ni₄Ti₂]^2-(D_4h)具有较大芳香性的主要原因; 其中Ti原子主要提供d_z²与s轨道, 而Ni原子主要利用其d_z²与d_x²-y²轨道形成正方形环, 它们之间构成了球状的d轨道环流, 且[Ni₄Ti₂]²⁺(D_4h)和[Ni₄Ti₂]^2-(D_4h)中还有非常明显的π轨道环流.     

Theoretical predictions of the structure, gas-phase acidity and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure, acidity (free energy of deprotonation, ΔG^o), and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid (3,4-dithiohydroxy-3-cyclobutene-1,2-diselenone, H₂C₄Se₂S₂) are reported. The global minimum found on the potential energy surface of 1,2-diseleno-3,4-dithiosquaric 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, cyclobutanediselenone, and cyclobutenedithiol. The computed aromatic stabilization energy (ASE) by homodesmotic reaction (Eq 1) is −20.1 kcal/mol (MP2(fu)/6-311+G** //RHF/6-311+G**) and −14.9 kcal/mol (B3LYP//6-311+G**//B3LYP/6-311+G**). The aromaticity of 1,2-diseleno-3,4-dithiosquaric acid is indicated by the calculated diamagnetic susceptibility exaltation (Λ) −17.91 (CSGT(IGAIM)-RHF/6-311+G**//RHF/6-311+G**) and −31.01 (CSGT(IGAIM)-B3LYP/6-311+G**//B3LYP/6-311+G**). Thus, 1,2-diseleno-3,4-dithiosquaric acid fulfils the geometric, energetic and magnetic criteria of aromaticity. The calculated theoretical gas-phase acidity is ΔG^o _1(298K)=302.7 kcal/mol and ΔG^o _2(298K)=388.4 kcal/mol. Hence, 1,2-diseleno-3,4-dithiosquaric acid is a stronger acid than squaric acid(3,4-dihydroxy-3-cyclobutene-1,2-dione, H₂C₄O₄). Received: 11 April 2000 / Accepted: 7 July 2000 / Published online: 27 September 2000     

杂硼原子簇XB6+ (X=C,Si, Ge,Sn, Pb)的稳定性和芳香性

利用从头算MP2方法和密度泛函理论B3LYP和B3PW91方法, 研究了杂硼原子簇XB₆⁺ (X=C, Si, Ge, Sn, Pb)的结构、稳定性及化学键合情况. 对C, Si, Ge, B使用6-311+G(d)基组, 对Sn和Pb使用LANL2DZ赝势基组. 研究结果表明, 具有C_s对称性的假平面XB₆⁺ (X=C, Si, Ge, Sn, Pb)结构是势能面上的全域极小点, 其稳定性要高于C_6v对称性的锥形结构和C₂对称性的假锥形结构. 在B3LYP水平上, 对这些异构体的势能面的极小点进行了自然键轨道(NBO)的分析; 对最稳定构型的最高占据分子轨道(HOMO)和最低空轨道(LUMO)能级差、分子轨道(MO)和核独立化学位移(NICS)进行了计算和讨论. 分析了杂原子和硼原子间、相邻硼原子间的键合情况, 讨论了最稳定构型的芳香性质.     

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