A Density Functional Investigation of Fluorinated Silicon Clusters

Density functional calculations have been carried out on a series of fluorinated empty cages X_nF_n(n=2–20) with X?Si, Ge, and Sn. It indicates that the fullerene‐like cage structure with pentagons turns out to be the most stable with n increasing, and the stability of the X_nF_n isomers increases with the number of five‐membered rings. The HOMO‐LUMO gap for Ge (n=6, 10) cages is found to be even larger than the values for Si cages, though in bulk Ge has a smaller band gap than Si. Moreover, calculation of the Gibbs free energy of oligomerization reaction of SiF→1/n (SiF)_n showed that this reaction is exothermic even at 900 K, indicating the favorability of their formation from the SiF monomer.     

Density Functional Calculations of C–NO_2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nit- roalkane energetic materials have been calculated by using the three hybrid density functional theory (DFT) methods B3LYP, B3PW91 and B3P86 with 6-31g** and 6-311g** basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g** and 6-311g** basis sets can obtain satisfactory bond dissociation energies (BDEs), which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g** method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal·mol-1 (for the BDE of tC4H9–NO2), which further proves the reliability of B3P86/6-311g** method. In addition, it is noted that the BDEs of C–NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mol-1.     

Density Functional Theoretical Study of Polynitrogen Compounds N5+Y− (Y=B(CF3)4, BF4, PF6 and B(N3)4)

The structures and the stabilities of polynitrogen compounds N₅⁺Y^? [Y=B(CF₃)₄, BF₄, PF₆, and B(N₃)₄], as the potential high energy density compounds, have been investigated at the B3LYP/6‐31G(d,p) and B3LYP/6‐311+G(d,p) levels. On the basis of our geometry optimization calculations, the structural properties of the N₅⁺Y^? compounds are discussed, and it is found that the combination of the N₅⁺ cation and the Y^? anions leads to distortion of the structures of the Y^? anions. Based on the TS calculations for the N₂‐loss dissociations of the N₅⁺Y^? compounds, the stabilities of these compounds are discussed, and the following conclusion can be drawn that among the four compounds, N₅⁺B(CF₃)₄^? is the most stable one and N₅⁺B(N₃)₄^? is the most unstable, and the relative stability of these compounds is always consistent using different basis sets. From these discussions, it is revealed that there are close correlations between the stuctrual distortions of the Y^? anions and the stabilities of the N₅⁺Y^? compounds, and between the nitrogen content in the compounds and the stabilities of the N₅⁺Y^? compounds.     

过渡金属混合簇Nb2Rh2的密度泛函研究

采用密度泛函方法(DFT)研究了过渡金属混合簇Nb_mRh_n (m, n≤2) 的结构、稳定性规律及它们的成键情况. 结果表明, Nb—Nb键较强, Rh—Rh键较弱, 而Nb—Rh键的强度则介于两者之间. 在Nb₂Rh₂直线和折线构型中, 金属键有强弱交替的现象. Nb₂Rh₂的各种构型在自旋多重度较小时稳定.     

原子团簇Ge11同分异体的密度泛函研究

用分子图形软件设计出多种Ge_11原子团簇模型,使用B3LYP密度泛函方法进行几何构型优化和振动频率计算,比较了14种同分异构体的总能量,得到了新的基态构型。锗原子团簇的大部分原子以三、四、五和六配位成键。在Ge11原子团簇中,双角反四棱柱衍生出来的构型能量较低,它是设计大分子锗原子团簇初始模型的重要结构单元。由三棱柱演变的构型能量居中。带心结构和共边构型带有高配位原子,其总能量较高,是不稳定的结构。     

Density Functional Theory Study of HCN Adsorption on Small Gold Clusters

A theoretical study was carried out on the adsorption of hydrocyanic acid on small Aun （n ≤ 7） clusters using density functional methods. For HCN adsorption on gold clusters, no dependence was found with respect to the even-odd alternation in relation to the number of gold atoms in the cluster. The HCN molecule is adsorbed at simple adsorption sites （1-fold coordination）, perpendicular to the adsorption site. The largest adsorption energy is only about 74.61 kJ·mol^-1, which indicates that the HCN molecule does not decompose and the C-N bond retains triple bond, and that the C-H and C-N stretching frequencies are only weakly perturbed. The adsorbed C-N and C-H stretching frequencies are blue- and red-shifted compared with the values of free HCN, respectively.     

Theoretical Study of the N-NO2 Bond Dissociation Energies for Energetic Materials with Density Functional Theory

The N-NO2 bond dissociation energies （BDEs） for 7 energetic materials were computed by means of accurate density functional theory （B3LYP, B3PW91 and B3P86） with 6-31G＊＊ and 6-311G＊＊ basis sets. By comparing the computed energies and experimental results, we find that the B3P86/6-311G＊＊ method can give good results of BDE, which has the mean absolute deviation of 1.30kcal/mol. In addition, substituent effects were also taken into account. It is noted that the Hammett constants of substituent groups are related to the BDEs of the N-NO2 bond and the bond dissociation energies of the energetic materials studied decrease when increasing the number of NO2 group.     

Density Functional Calculations on a Double Hydrogen-bonded Dimer

Density functional theory (DFT) calculations on a double hydrogen-bonded dimer of o-hydroxybenzoic acid were carried out at the B3LYP/6-31G^* level. The optimized geometry of the dimer closely resembles that of the crystal. The calculated results show that the total energy of the dimer is much lower than the sum energies of the two monomers, and the average strength of the double hydrogen bonds is about 38.37 kJ/mol. In order to probe the origin of the interactions in the dimer, natural bond orbital analyses were performed. The thermodynamic properties of the title compound at different temperatures have also been calculated on the basis of vibrational analyses and ΔGT, the change of Gibbs free energy for the aggregation from monomer to the dimmer, is 26.47 kJ/mol at 298.15 K and 0.1MPa, implying the spontaneous process of forming the dimer. The correlation graphics of Sm^0 Hm^0 and temperatures is depicted.     

烃类化合物碳氢键解离焓的密度泛函理论研究

通过比较10种密度泛函方法对烃类化合物碳氢键解离焓的计算精度, 发现新型密度泛函BMK方法具有最高的计算精度. 利用该方法计算了包含饱和链烃,、不饱和链烃、脂环烃和芳香烃在内的172个烃类化合物的碳氢键解离焓,计算均方根误差仅为7.95 kJ•mol^－1, 线性拟合常数为0.985. 通过自然键轨道法分析发现, 烃类物质的碳氢键解离焓与母体的碳氢键杂化轨道成分p%, 自由基奇电子轨道杂化成分p%及自由基的自旋密度三个参数之间存在较好的定量关系. 此外, 饱和链烷烃及不饱和链烃的碳氢键解离焓与碳氢键键长之间也存在较好的线性关系.     

1, 4-Pyrone Effects on O-H Bond Dissociation Energies of Catechols in Flavonoids: A Density Functional Theory Study

In recent years, there has been growing interest in selecting efficient antioxidants with low toxicity to reduce the damage of free radicals. Among these antioxidants, flavonoids have been paid much attention, owing to their excellent antioxidative and pharmacological activities1. Up to now, many efforts have been given to summarize the structure-activity relationships (SAR) for flavonoids. It has been widely accepted that two structural factors are critical for flavonoids to enhance the…     

The Structural and Electronic Properties of Tin-Based Heteroatom Clusters Studied by the Density Functional Theory

The focus of this study is on compound clusters and, due to the existence of many phases with different structural properties, tin-based materials have been chosen as the reference case. The clusters considered below are of two types: in the first case the clusters have the skeleton of the pure tin clusters and are doped with oxygen and aluminum atoms with composition Sn _x Y _y with Y = Al, O, x = 1, 10 and y = 1, 2. In the second case the clusters have a rutile lattice with a columnar or a spherical shape and a size up to 80 atoms and are doped with a number of aluminum atoms up to 20. The calculations are based on the Density Functional Theory (DFT) and the results describe the cluster structure, its binding energy and the density of states (DOS). The general indication of the calculations is that the additive coordinates outside, rather than inside, the pristine skeleton with the formation of hybrid bonds with properties similar to the ones of the host atoms. Also conspicuous effects of hybridization are observed in the electronic structure and, due to these effects, the binding energy may decrease with respect to the one in the undoped clusters.     

新型异三聚体超分子的密度泛函计算及其晶体结构、热性能研究

运用密度泛函B3LYP方法,在6-31G^*水平上设计优化了丁二酮肟与苯甲酸通过四重氢键构筑的异三聚体超分子、此异三聚体的总能量比三个单体总能量之和低144．0816kJ／mol,四重氢键的平均键能约为36．0204kJ／mol,属中高强度;298．15K时吉布斯自由能变化△GT=-41．70kJ／mol,显示形成三聚体的反应可自发进行．在此设想指导下,实验合成出相关异三聚体,并得到了晶体结构,与预测的结构极为相似．热性能分析显示,此三聚体不是三个单体的简单叠加,而是一个不同于其单体的全新化合物．理论计算数据与实验结果的一致性,证明了理论计算的预见性．     

X-H (X = C, N, O, Si, P, S) 键离解能计算的密度泛函方法研究

使用了不同密度泛函方法计算X-H (X = C, N, O, Si, P, S) 键离解能，并分析不同密度泛函方法的计算精度。研究发现大多数密度泛函方法包括B3LYP, B3P86, B3PW91, G96LYP, PBE1PBE，和BH&HLYP都明显低估键离解能13-25 kJ/mol。该现象与是否使用无限基组无关，因为即使使用无限基组键离解能仍然被低估。因此密度泛函方法不适合用于键离解能的估算。其中B3P86方法的偏差最小。进一步分析表明，使用限制性开壳层计算并无任何优势，在大多数情况下非限制性开壳层计算实际上比限制性开壳层计算要好。最后，我们发现了密度泛函方法对键离解能的低估是系统的，因此建议利用校准后的UDFT/6-311++G(d, p)方法计算化学键离解能。     

碳氟键均裂解离能的密度泛函理论研究

运用6种密度泛函方法(B3LYP, B3P86, B3PW91, PBE1PBE, MPW1B95, MPW1K)对15个含氟有机化合物的碳氟键均裂解离能进行理论计算, 得到的理论值与实验值比较, 发现B3P86方法用于碳氟键均裂解离能的计算相对可靠. 使用验证后的理论方法对含氟杂环有机化合物和卤氟烃中的碳氟键均裂解离能进行了预测和分析, 并进一步讨论了α-取代基效应以及Hammett型取代基效应对碳氟键均裂解离能的影响.     

密度泛函方法研究噻唑橙类菁染料的光谱性质

4-(二乙基氨基)丁基取代的三甲川噻唑橙(DEAB-TO3)具有极低的本底荧光可用于核酸检测。本文运用密度泛函理论研究了4-(二乙基氨基)丁基取代的一甲川噻唑橙(DEAB-TO1)和DEAB-TO3光谱性质。基态和激发态几何优化显示激发态构型高度扭曲;光谱计算和轨道分析得出第一激发态是暗态,出现了扭曲的分子内电荷转移;势能曲线计算可知,DEAB-TO1和DEAB-TO3有着极低的能隙和旋转能垒。以上结果解释了本底荧光极弱的实验现象。     

大体系分区密度泛函计算方法

提出一种新的对大体系进行分区密度泛函计算的方法。将大体系划分为若干较小的区,每个小区是一个相对独立的量子力学子体系,计及其它区势场的影响和电子的Pauli排斥作用,可以进行相对独立的密度泛函计算。对各子体系求解单电子方程：(F^k F^kp)C^k=S^kC^kε^k K=A,B,C,…式中F^k,C^k,S^k,ε^k分别为子体系K的Fbck矩阵、轨道系数矩阵、基组重叠矩阵和本征值矩阵,F^kb起强制属于不同子体系的占据轨道之间保持正交的作用。得到的轨道是分区定域化的,汇总各区的计算结果得出整个体系的电子结构信息。通过对一些较大分子的计算,考察了几种因素对分区计算精度的影响。结果表明,提出的方法是可行的,通过控制各区基组的大小,可以基本消除基组截断误差,得到精确的计算结果。对于足够大的体系,本方法是一种线性标度算法;和文献报道的相关方法比较,更容易用于对体系的某些区域进行特别的计算研究。     

C24和B12N12团簇的结构与稳定性

本文采用量子化学密度泛函理论的B3LYP/6-31G*方法,对C24和B12N12团簇的12种异构体进行了优化,并对它们的几何构型、振动频率、核独立化学位移(NICS)和结合能进行了理论探讨, 比较了C24和B12N12团簇结构的稳定性。研究表明:C24团簇的最稳定几何构型为类石墨结构d,B12N12团簇的最稳定结构为4/6笼状结构g。C24异构体的稳定性大小顺序为d > b > f > c > a > e。B12N12团簇异构体稳定性大小顺序为a > f> c> d > e >b。     

Density functional study on highly energetic organic azides with empirical formula Cn(N3)m

The structural, electronic, and thermodynamical properties of C_n(N₃)_m (n = 1–7) (m = 4, 6) organic azides have been investigated using Density Functional Calculations. The ground state structures of organic azides were compared with C_nH_m (n = 1–7) (m = 4, 6) cumulenes which shows their higher relative stability. The stability and reactivity of organic azides were analyzed by calculating the HOMO-LUMO gap, binding energies, and harmonic frequencies of the azides. The binding energy and formation energy of C_n(N₃)_m (n = 1–7) (m = 4, 6) organic azides suggest their energetic stability. The structural analysis of the azide group in C_n(N₃)_m (n = 1–7) (m = 4, 6) organic azide shows a tendency to stabilize at a maximum separation between functional azide groups. Ionization potential, electron affinities, and global hardness have been computed for C_n(N₃)_m (n = 1–7) (m = 4, 6) organic azides and the odd–even alternation rule was observed. The molecular dynamic simulation performed at 300 K for 1 fs confirms organic azide's structural stability at room temperature, except for C₄(N₃)₄, and the members of their family can be synthesized.     

Reactivity Analysis in Diamond Surfaces with a Density Functional Calculation

The surface states of different diamond surfaces are studied using total and partial density of states (DOS) curves, and are related qualitatively to the reactivity of these surfaces, which are important in the process of diamond growth. The calculations combined atomic and functional density approaches with MARVIN and LMTO–ASA codes, respectively. In the atomic calculation, the interatomic potentials are as follows: the b parameter in the Morse potential is 0.5523 Å, the A and B parameters in the nonbonding Lennards–Jonnes potential are 639.6258 eV Ź² and 3.632 eV Å⁶, and the three-body bending potential K₃, is 0.7797 eV rad². To validate these results, the elastic constants were evaluated, finding a good agreement with the experiment. With these potentials, a slab, for each of the diamond surfaces, of 40 carbon atoms with periodic conditions in two dimensions was optimized. The output coordinates were used for DOS calculations. These results were later verified with a surface-cluster calculation of HOMO and LUMO frontier orbitals. These were calculated using a 9-carbon cluster with the DGauss code. In a nonrelaxed surface, two surface states are identified: the first is an occupied state placed at the center of the gap, and the other, adjacent to a valence band maximum, is an empty state of p character and, therefore, potentially able to participate in chemical interactions. In the relaxation process of the surface, the surface states become narrower as the valences of the surface carbons are saturated; in this case the isolated p state participates in dangling bonds. With the monohydrogenation, the surface state placed at the center of the gap of the relaxed surface, becomes a subsurface state, that is, the highest density is not at the surface layer, but in inner layers. As a consequence, the reactivity diminishes. Therefore, it is possible to conclude that the study of surface states could give predictive information about the reactivity of surfaces.     

Performance of Density Functional Theory for Transition Metal Oxygen Bonds

The accuracy of density functional theory (DFT) limits predictions in theoretical catalysis, and strong chemical bonds between transition metals and oxygen pose a particular challenge. We benchmarked 30 diverse density functionals against the bond dissociation enthalpies (BDE) of the 30 MO and 30 MO⁺ diatomic systems of all the 3d, 4d, and 5d metals, to test universality across the d-block as required in comparative studies. Seven functionals, B98, B97-1, B3P86, B2PLYP, TPSSh, B3LYP, and B97-2, display mean absolute errors (MAE) <30 kJ/mol. In contrast, many commonly used functionals such as PBE and RPBE overestimate M−O bonding by +30 kJ/mol and display MAEs from 48–76 kJ/mol. RPBE and OPBE reduce the over-binding of PBE but remain very inaccurate. We identify a linear relationship (p-value 7.6 ⋅ 10⁻⁵) between the precision and accuracy of DFT, i. e. inaccurate functionals tend to produce larger, unpredictable random errors. Some functionals commonly deviate from this relationship: Thus, M06-2X is very precise but not very accurate, whereas B3LYP* and MN15-L are more accurate but less precise than M06-2X. The best-performing hybrids have 10–30 % HF exchange, but this can be relieved by double hybrids (B2PLYP). Most functionals describe trends well, but errors comparing 5d to 4d/3d are ∼10 kJ/mol larger than group-wise errors, due to uncertainties in the spin-orbit coupling corrections for effective core potentials, affecting e. g. Pt/Pd or Au/Ag comparisons.