配合物[M(CO)3(PPh2py)2](M=Fe,Ru)异构体的理论研究

对PPh2py配合物[M(CO)3(PPh2py)2](M=Fe, Ru)的三种构型的异构体1-6进行了研究. 其中PPh2py以两个P原子与M配位形成HH构型1(Fe)和4(Ru), 以一个P和一个N原子与M配位形成HT构型2(Fe)和5(Ru), 以两个N原子与M配位形成HH’构型3(Fe)和6(Ru). 结果表明, (1) PPh2py中P原子对HOMO轨道的贡献最大, PPh2py作为电子给体时易以P原子与金属原子结合. (2)从分子能量和相互作用能数据表明, 配合物中HH构型最稳定, HH'构型最不稳定, 这与合成产物为HH构型的结果一致. (3) 键长和Wiberg键级均表明P—M键比N—M键结合力强. P、M原子间存在σ键, 而N、Fe原子间仅存在nN→n*M或nN→σ*M-P的电荷转移作用. (4) HH构型中M对HOMO的贡献最大, PPh2py向M的电荷转移最强, 使M的负电荷最大, 故HH构型最易作为电子给体以M原子与第二个金属配位形成双核配合物.     

(Mg3N2)n(n=1～4)团簇结构与性质的密度泛函研究

用密度泛函理论(DFT)的杂化密度泛函B3LYP方法在6-31G*基组水平上对(Mg3N2)n(n=1～4)团簇各种可能的构型进行几何结构优化,预测了各团簇的最稳定结构.并对最稳定结构的振动特性、成键特性、电荷特性和稳定性等进行了理论分析.结果表明:(Mg3N2)n=1～4团簇易形成笼状结构,其最稳定构型中N原子配位数以3、4较多见;团簇主要由Mg-N键组成,Mg-N键长为0.194～0.218nm,Mg-Mg 键长为0.262～0.298 nm;N原子的平均自然电荷为-2.06 e,Mg原子的平均自然电荷为 1.37 e;(Mg3N2)2团簇有相对较高的动力学稳定性.     

硅基稀土掺杂团簇MSi_7~q(M=Eu,Sm,Yb;q=0,-1)结构预测与光电子能谱研究

采用随机踢球模型结合密度泛函理论,在PBEPBE/RE/SDD/Si/6-311+G(d)水平下研究了中性和阴性的硅基稀土掺杂团簇MSi_7~q(M=Eu,Sm,Yb;q=0,-1)的几何结构、稳定性及电子和磁学性质.计算结果表明,阴性团簇的基态结构是在五角双锥的双锥侧面外法向方向加入一个Si原子而形成的3D结构,并且稀土原子M处于五角双锥的顶点;中性团簇的最低能结构是一个畸变的双帽八面体,并且M原子处于八面体的赤道面上.SmSi_7~-团簇在这3种稀土掺杂的团簇中具有最高的平均结合能和掺杂能,是这3种稀土掺杂团簇中最稳定的一种.Si_7团簇是非磁性团簇,但是当M原子(M=Eu,Sm,Yb)掺入其中时,由于镧系元素独特的原子磁性,使其变成了磁性团簇.此外,还模拟了各团簇前几种低能异构体的光电子能谱.     

碱土金属叠氮化合物(MgN_6)_n(n=1～5)团簇结构与性质的密度泛函理论研究

利用密度泛函理论在B3LYP/6-311G*水平上对碱土金属叠氮化合物(MgN6)n(n=1～5)团簇各种可能构型进行了几何优化,预测了各团簇的最稳定结构。并对最稳定结构的成键特性、电荷分布、振动特性及稳定性进行理论研究。结果表明,叠氮化合物中叠氮基以直线型存在,MgN6团簇最稳定结构为直线型;(MgN6)2团簇最稳定结构为Mg2N2四元环平面结构;(MgN6)n(n=3～5)团簇最稳定结构是由2个叠氮基与2个Mg原子首先构成近似菱形,再由近似菱形延伸形成的链状结构。叠氮基中间的N原子显示正电性,两端的N原子显示负电性,且与Mg直接作用的N原子负电性更强,金属Mg原子和N原子之间形成很强的离子键。(MgN6)n(n=1～5)团簇最稳定结构的IR光谱分为4个部分,其最强振动峰均位于2209～2313cm-1,振动模式为叠氮基中N-N键的反对称伸缩振动。稳定性分析显示,(MgN6)3和(MgN6)5团簇相对于其他团簇较为稳定。     

C59XH(X=N,B)自由基加成区域选择性的理论研究

用半经验的AMl方法,对C59XHCl2n(X=N,B;n=1～2)和C60H2Cl2n(n=1～2)的异构体进行几何构型全优化和振动频率计算,结合密度泛函B3LYP／6—31G^*单点能计算确定各异构体的相对稳定性．对比C59XH(X=N,B)和C60H2的H2加成方式,计算结果表明H2或Cl2加在碳笼官能化部分的邻近位置在能量上都是有利的：C59NH和C59BH自由基多加成物区域选择性的差别可归因于N原子和B原子电子性质的不同;立体效应是导致H2和Cl2加成方式不同的主要原因．     

AlmN2－ (m=1～8)团簇的结构与稳定性

用密度泛函理论(DFT)的B3LYP方法, 在6-311G*水平上对AlmN2－ (m =1～8)团簇的几何构型、电子结构、振动频率、电荷分布与成键方式进行了理论研究. 结果表明, AlmN2－ 团簇的基态结构有两种基本构型, 一种是以N—N键为核心, 周围与Al原子配位形成的； 另一种是由两个AlnN(n＜m)分子碎片通过共用Al原子或Al—Al键相互结合形成的. 对AlnN分子碎片相互结合形成的基态结构亲和能讨论得到, m为奇数的AlmN2－团簇比m为偶数的稳定.     

羰基硼化合物(BCO)n(n==1～12)的理论研究

运用HF/3-21G方法和密度泛函理论(DFT)的B3LYP/6-31G*方法, 对羰基硼化合物(BCO)n(n=1~12)的各种可能结构进行了优化, 对在B3LYP/6-31G*水平上得到的几何构型、电子态、结合能、振动频率、核独立化学位移(NICS)、能量二次差分和热力学性质进行了理论研究, 得到了(BCO)n(n=1~12)结构的稳定性信息. 十二种基态结构都是端配位(μ1-CO), (BCO)n(n=1~3, 5, 6) 的基态是线型或平面结构, (BCO)n(n=4, 7~12)的基态是笼状结构; B—C平均键能呈现奇偶交替现象, 偶数的结构比奇数稳定; 能量二次差分得到同样的结论；羰基的振动频率与实验值非常吻合; 热力学性质的研究对实验具有重要的指导意义.     

钛,锆和铪羰基化合物M(CO)_n(M=Ti,Zr,Hf;n=4～7)的理论研究

利用密度泛函方法对标题化合物的平衡几何、热化学及振动频率进行了理论预测,发现这3种金属原子都有相似的M(CO)n(n=4～7)结构.全局最低构型对M(CO)7都是单态C3v戴帽八面体7S-1,对M(CO)6都是三重态D3d畸变八面体6T(而对应的单重态M(CO)5仅比它低不到21 kJ·mol-1).对M(CO)n(n=5,4)都是三重态6S-1,其构型分别为从6T中移去1个或2个CO基的衍生物5T和4T.此外,五重态的D3h的三角双锥M(CO)5和单态的Td四面体M(CO)4以及能量更高的含有C和O同时与金属成键的独特配位CO基的M(CO)6和M(CO)3也被发现.最后,给出M(CO)7→M(CO)6+CO反应的离解能.并讨论了金属18价电子的Ti(CO)7存在的可能性.     

AlmN2 (m=1～8)团簇结构与稳定性的DFT研究

用密度泛函理论(DFT)的B3LYP方法，在6-31G*水平上对AlmN2(m=1～8) 团簇的几何构型、电子结构、振动频率和热力学性质进行了理论研究. 结果表明， AlmN2团簇的基态结构有两种基本构型， m≤2的结构是以N－N键为核心周围与Al原子相配位形成的; m >2的结构是由两个AlnN（n< m）分子碎片通过共用Al原子或Al－Al键相互结合形成的，这为较快找到AlmN2团簇基态结构提供了一条有效途径. 通过对基态结构能量二次差分的讨论，得到了m为偶数的AlmN2团簇比m为奇数的稳定.     

AlCn及AlC+n(n=1-4)原子簇的理论研究

用ab initio能量解析梯度法，在UHF（RHF）／3－21G水平上优化得到AlCn，AlCn+（n＝1～3）的四十三个构型，在RHF（UHF）／3－21G水平上优化得到AlC4，AlC4+的二十六个构型，CISD能量．从能量角度所研究的AlCn，AlCn+（n=1～4）构型中，最稳定的构型均是线性构型，且Al全都在端点上，同文献报导的理论和实验结果相一致．此外，还研究了原子簇的离子化能，原子平均结合能以及原子簇的碎片化通道和碎片化能，并计算了上述最稳定构型的谐振动光谱常数．     

Theoretical Study on the Electronic Structures of Small TinNm （n ＋ m = 5, 6） Clusters

45 isomers of TinNm （n ＋ m = 5, 6） clusters, including linear, some planar and some stero configurations, have been predicted by density functional theory method. For five-atom clusters Ti3N2 and Ti2N3, the most stable structures are trigonal bipyramid in D3h symmetry, and for TiaN cluster, the isomer with one nitrogen atom occupying the center of quasi-tetrahedron is the most stable. In the isomers of Ti4N2 and Ti3N3, the planar networks are more stable, but for Ti2N4, the six-membered ring configuration is the most favorable. Most linear structures can form weak-strong bonds alternately with higher energy. As regards to planar structures, the more Ti-N bonds are formed, the more stable they will be; for stero closed polyhedral isomers, their energies are lower.     

Density functional theoretical study of a series of binary azides M(N3)n (n = 3, 4)

Geometrical structures of a series of binary azides M(N3)n (M = elements in groups 3 and 13 (n = 3) and in groups 4 and 14 (n = 4)) were investigated at the B3LYP/6-311+G level of theory. Our calculations found that binary group 3 triazides M(N3)3 (M = Sc, Y, La) and binary group 4 tetraazides M(N3)4 (M = Ti, Zr, Hf) turn out to be stable with all frequencies real having a similar linear M-N-NN structural feature, as previously reported for M(N3)4 (M = Ti, Zr, Hf). However, binary azides of group 13 M(N3)3 (M = B, Al, Ga, In, Tl) and group 14 elements M(N3)4 (C, Si, Ge, Sn, Pb) with bent M-N-NN bond angles differ obviously from binary group 3 and 4 azides in geometrical structure. These facts are mainly explained by the difference in electronic density overlap between the central atom and the alpha-N atoms of the azido groups. Two lone-pair electrons on the sp hybridization alpha-N atoms in the binary group 3 and 4 azides donate electron density into two empty d orbitals of the central transition metal atom and a pair of valence bonding electrons, resulting in the alpha-N atoms acting as a tridentate ligand. The sp2 hybridization alpha-N atoms of the binary group 13 and 14 azides only give one valence electron to form one valence bonding electron pair acting virtually as monodentate donors.     

第ⅢA族金属叠氮多聚体结构和性质的理论研究

在DFT-B3LYP/SDD水平上计算研究了第A族金属叠氮多聚体(Me₂MN₃)_n(n=1～3,M=Ga,Al)的结构和性质.多聚体(Me₂MN₃)₂,3的各优化构型均为环状,通过一子体系叠氮基的α-N和另一子体系的金属Ga或Al相连.二聚体(Me₂MN₃)₂中含M2N2平面四元环结构,三聚体(Me₂MN₃)₃具有结合能相近的扭船式和椅式两种构象,均含M₃N₃六元环结构.与单体相比,多聚体的几何参数变化较大.报道了它们在不同温度下的热力学性质,发现叠氮二甲基镓和铝体系以二聚体形式存在.     

Geometry and stability of BenCm (n=1-10; m=1, 2, ..., to 11-n) clusters

Density functional theory B3PW91/6-31+G* calculations on BenCm (n=1-10; m=1, 2, ..., to 11-n) clusters have been carried out to examine the effect of cluster size, relative composition, binding energy per atom, HOMO-LUMO gap, vertical ionization potential, and electron affinity on their relative stabilities. The most stable planar cyclic conformations of these clusters always show at least a set of two carbon atoms between two beryllium atoms, while structures where beryllium atoms cluster together, or allow the intercalation of one carbon atom between two of them, generally seem to be the least stable ones. Clusters containing 1, 2, and 3 beryllium atoms (Be2C8, Be3C6, Be2C6, BeC6, Be2C4, BeC4, Be2C2, and BeC2) are identified as clusters of "magic numbers" in terms of their high binding energy per atom, high HOMO-LUMO gap, vertical ionization potential, and second difference in energy per beryllium atom.     

Geometric, electronic, and bonding properties of AuNM (N = 1-7, M = Ni, Pd, Pt) clusters

Employing first-principles methods, based on density functional theory, we report the ground state geometric and electronic structures of gold clusters doped with platinum group atoms, Au(N)M (N = 1-7, M = Ni, Pd, Pt). The stability and electronic properties of Ni-doped gold clusters are similar to that of pure gold clusters with an enhancement of bond strength. Due to the strong d-d or s-d interplay between impurities and gold atoms originating in the relativistic effects and unique properties of dopant delocalized s-electrons in Pd- and Pt-doped gold clusters, the dopant atoms markedly change the geometric and electronic properties of gold clusters, and stronger bond energies are found in Pt-doped clusters. The Mulliken populations analysis of impurities and detailed decompositions of bond energies as well as a variety of density of states of the most stable dopant gold clusters are given to understand the different effects of individual dopant atom on bonding and electronic properties of dopant gold clusters. From the electronic properties of dopant gold clusters, the different chemical reactivity toward O(2), CO, or NO molecule is predicted in transition metal-doped gold clusters compared to pure gold clusters.     

原子簇NixFe和NiFex(x=1-5)电子性质的DFT研究

More than one hundred models were designed to reflect the local structure and electronic property of Ni-Fe amorphous alloys. After calculating by DFF method, a series of configurations of clusters NixFe and NiFex （x = 1 - 5） were gained. The configurations, which possessed the lowest energies and non-imaginary frequencies, were considered the most stable optimized structures. The catalytic activity, charge and magnetic properties were analyzed and discussed. The different Fe content changed the catalytic properties of clusters through altering the value of Fermi level of every cluster. However the density of state （DOS） nearby Fermi level and average 3d orbital population of atom Ni, which were also important properties related to the catalytic activation, were little changed. Based on the Fermi level, the activity of catalyst toward hydrogenation reaction would be considered best when the ratio of Ni to Fe was close to 1. The Fermi level of clusters was far distant to the level of nitrogen in singlet state. It would be the reason why the reaction condition in ammonia synthesis and nitrogen fixation process was rigorous. When Fe atom contents were higher than 75% （NiFe3）, the electrons transferred from atom Fe to Ni, but when the ratio was decreased, the transfer was reversed. The ratio of atoms of local structure also played an important role in the aspect of electron transition. On the average 3d orbital population of atom Fe, the average magnetic moments of Fe atoms in clusters were calculated. When Fe atom contents were 50% nearly, the average magnetic moment achieved the highest point.     

Theoretical study of the adsorption of H on Si n clusters, (n=3-10)

A recently proposed local Fukui function is used to predict the binding site of atomic hydrogen on silicon clusters. To validate the predictions, an extensive search for the more stable SinH (n=3-10) clusters has been done using a modified genetic algorithm. In all cases, the isomer predicted by the Fukui function is found by the search, but it is not always the most stable one. It is discussed that in the cases where the geometrical structure of the bare silicon cluster suffers a considerable change due to the addition of one hydrogen atom, the situation is more complicated and the relaxation effects should be considered.     

Cyano-bridged CuII-MV (M = Mo, W) bimetallic layered complexes exhibiting novel honeycomblike structures

Self-assembly of a new precursor [Cu(L)](ClO4)2 (1) (L = macrocyclic ligand) with octacyanometalates [M(CN)8]3- (M = Mo, W) produced two-dimensional cyano-bridged Cu(II)-M(V) bimetallic assemblies [Cu(L)]3n[M(CN)8]2n.6nH2O [M = Mo (2), W (3)] with novel honeycomblike structures, characterized by spectroscopic data, single-crystal X-ray diffraction studies, and magnetic measurements. The crystallographic determination reveals that compounds 2 and 3 are isostructural and crystallize in the triclinic system (P). The Cu atom in a distorted octahedral environment experiences a tetragonal elongation of apical nitrogen atoms exhibiting average Cu-Nax lengths of 2.566 Angstroms for 2 and 2.593 Angstroms for 3, which accounts for the Jahn-Teller effect of a Cu(II) ion. The Cu-NC angles are magnetically important, ranging from 135.7 to 159.2 degrees. Three types of L in the crystal lattice are observed, which are dependent on the relative positions of the pendant hydroxyl groups with respect to the CuN4 basal plane. The positions are correlated with hydrogen bonding of OH groups to neighboring atoms. The magnetic data indicate that ferromagnetic and antiferromagnetic interactions between Cu(II) and M(V) through the CN linkage coexist.     

Ab initio finite field (hyper)polarizability computations on stoichiometric gallium arsenide clusters GanAsn (n=2-9)

We report reliable ab initio finite field (hyper)polarizability values at Hartree-Fock and second order Moller-Plesset perturbation theory (MP2) levels of theory for different geometrical configurations of small gallium arsenide clusters Ga(n)As(n) with n=2-5. We relied on all-electron basis sets and pseudopotentials suitable for (hyper)polarizability calculations. In each case, we used structures that have been established in the literature after we optimized their geometries at B3LYP/cc-pVTZ-PP level of theory. Our results suggest that the first order hyperpolarizability (beta) is much more sensitive to the special geometric features than the second order hyperpolarizability (gamma). For the most stable configurations up to ten atoms the second order hyperpolarizability at MP2 level of theory varies between 15 x 10(4) and 32 x 10(4) e(4)a0 (4)Eh(-3). In addition, we examined the polarizability per atom evolution versus the cluster size for Ga(n)As(n) with n=2-9. Our work extends earlier theoretical studies which were limited to eight atoms and exposes that the polarizability/atom of the most stable stoichiometric configurations up to Ga(9)As(9) continues the monotonic downward trend with increasing size. Lastly, from the methodological point of view, our analysis shows that apart from polarizabilities, augmented pseudopotentials yield reliable first and second hyperpolarizability values as well.     

A density functional study on nitrogen-doped carbon clusters CnN3- (n=1-8)

Using molecular graphics software, we designed numerous models of CnN3- (n=1-8). Geometry optimization and calculation on vibration frequency were carried out by the B3LYP density functional method. After comparison of structure stability, we found that the structures of ground-state CN3- and C2N3- are bent chains with a nitrogen atom at either end, whereas when n=3-8, the ground-state clusters show three branches, each with a nitrogen atom located at the end. When n=5-8, the longest branch of CnN3- is polyacetylenelike. When n=5 or 7, the longest branch is connected to the central sp2 carbon in a nonlinear manner. The CnN3- (n=1-8) with an even number of carbon atoms are more stable than those with odd numbers, matching the peak pattern observed in laser-induced mass spectra of CnN3-. The trend of such odd/even alternation is explained based on concepts of bonding characteristics, electron affinities, and incremental binding energies.