水合铀酰与尿嘧啶及其异构体配位体系的计算研究

运用B3LYP/6-311++G**(RLC ECP)方法研究[UO_2(Uracil)j(H_2O)k]~(2+)(Uijk,i为尿嘧啶6种异构体代号,j+k=5)配位体系的几何结构、振动光谱、结合能等性质,并用极化连续介质模型(PCM)考察了溶剂化效应。结果表明,在U1jk体系中,随着尿嘧啶配体数目的逐渐增加,U-Ouracil配位键和U=O键的键长逐渐伸长,水溶液中U=O键的伸缩振动频率逐渐减小,配离子的总结合能呈增加趋势,且气相中的线性拟合效果较好。在Ui14体系中,U-Ouracil的键长与U-OH_2的键长整体成负相关,与U=O键的伸缩振动频率成正相关,结合能最大的配离子并不是由能量最低的尿嘧啶异构体生成的。电子密度拓扑分析表明U-Ouracil键和U-OH_2键具有离子键性质。原子电荷分析揭示在配位过程中是由配体片段向铀酰发生了电子转移,且尿嘧啶的电荷转移量与该配体数目成负相关,其中Ur6异构体向铀酰离子转移电子数最多。     

水溶液中四种脱氧核苷酸与水合铀酰离子相互作用的理论计算

运用密度泛函理论(DFT)方法 M06-2X,在相对论有效势基组水平下研究水溶液中水合铀酰脱氧核苷酸配合物[UO2(dNMP)(H2O)3]2+(d NMP=单磷酸脱氧核苷酸)的几何结构、能量和电子结构等性质,采用极化连续介质模型考虑水溶剂环境。结果表明,除了脱氧胸苷酸,其它脱氧核苷酸均以磷酸基团P=O键与铀原子形成的配位结构最稳定,4类配离子U=O键长值较接近,但配位键键长有所不同。脱氧腺苷酸配离子的总结合能最大,但脱氧腺苷酸变形能最小。振动光谱分析发现U=O键和P=O键伸缩振动频率均发生红移,脱氧腺苷酸配离子的变化程度最大。电子密度拓扑分析表明U-dNMP键和U-OH2键具有离子键性质。电荷转移分析显示配体片段向铀酰离子发生了电子转移,脱氧腺苷酸配离子中配体转移的电荷数最多。轨道成分分析表明,高占据轨道贡献主要来自脱氧核苷酸配体π轨道,而低阶空轨道基本是U的5f轨道组成,脱氧鸟苷酸配离子的HOMO-LUMO能级差最小。     

金红石TiO2纳米团簇与铀酰相互作用的相对论密度泛函理论计算

采用全电子相对论密度泛函理论方法探索金红石型Ti O_2纳米团簇与铀酰的相互作用。考察金红石团簇模型(包括层数和表面积大小)变化对吸附铀形成复合物结构、吸附作用能等性质的影响,确定2层、表面积为1.1 nm×0.6 nm、包括63个原子的纳米团簇(标记为2L-Ti15)能够合理描述金红石纳米粒子性质的同时,还能节约计算资源。对2L-Ti15-[(UO_2)(H_2O)_3]~(2+)复合物计算表明,纳米团簇和铀酰存在共价键作用;优化得到U-O_(surf)键长0.233~0.238 nm,这一距离在已发现铀酰基配合物U-O距离范围内。在气相条件下,纳米团簇对铀酰吸附反应为放热过程(-3.02 e V);考虑溶剂介质环境的影响,反应则需要吸收少许能量(0.16e V)。U-O_(surf)键的能量分解发现,纳米团簇和铀酰的化学键作用为轨道相互作用主导的(占94%),它的静电吸引略大于Pauli排斥。基于电子密度的QTAIM(quantum theory of atoms in molecule)分析揭示,U-O_(surf)作用是介于离子和共价之间的配位键,其强度高于复合物中的U-OH_2键作用,但比U=O键弱。波函数分析表明,来自纳米团簇的O(2p)贡献HOMO轨道,并混有σ(U=O)成键性质,而LUMO轨道则为Ti(3d)修饰的U(5f)性质,复合物HOMO-LUMO带隙为2.40 e V,相对吸附前的纳米团簇半导体粒子的3.35 e V变窄。从吸收光谱角度而言,复合物体系可能在可见光区域具有更强的捕光性能。     

聚吡咯铀酰氧-酰亚胺配合物结构、振动光谱和基团交换反应的相对论密度泛函理论计算

铀酰氧-酰亚胺混合型配合物的合成为热力学稳定、动力学惰性铀酰（UO₂²⁺）和等电子体系铀酰亚胺（U（NR）₂²⁺）研究建立了纽带. 使用相对论密度泛函理论计算“Pacman”结构配合物[（THF）（OU^VIE）（A₂L）]（E=O、NH、NMe和NPh; A=H和Li; L为八齿氮供体低聚吡咯大环配体）. 优化得到U=O/U=N距离与实验值符合; 键级和电子结构分析显示U=O_exo/U=N有部分三重键特征; 由于受Li离子扰动,U=O_endo强度则介于单键和双键之间. 计算表明O=U=NH和O=U=O成键相近,均具有对称和反对称伸缩振动频率,而取代基Me和Ph耦合作用使得U=N-C吸收峰出现在高频区域. 酰氧和酰亚胺基团交换反应计算发现独特Pacman结构配合物的反应能相对五角双锥型配合物的有所降低,其中=O与=NMe交换反应最易实现; A离子变换不但能调控配合物结构和特征振动谱,还可降低基团交换反应能.     

C6H5—H…X分子间氢键的理论计算

利用密度泛函理论B3LYP方法, 在6-311+G(3df,2p)水平上对C6H5—H…X型分子间氢键进行了几何构型优化、氢键相互作用能、电子密度分布等计算. 其中C6H6为质子供体, HCOH、H2O、NH3、CH2NH和HCN为质子受体. 从电荷布居分析、自然键轨道等角度详细地讨论了C6H5—H…X 体系中, 共轭π键、O和N的不同键型结构对氢键形成的影响以及孤电子对与C—H 反键轨道之间的相互作用(n→σ*)等.     

系列杂氮钛三环中的分子内弱相互作用的紫外光电子能谱和理论研究

首次报道了四种杂氮钛三环的紫外光电子能谱,结合量子化学从头算和和电子密度拓扑分析方法对这些化合物的电子结构和分子内存在的弱相互作用进行了分析研究。通过优化STO-6G(d)基组(C:1s=5.66,ζ2s=1.675,ζ2p=1.66;O:ζ1s=7.67,ζ2s=2.25,ζ2p=2.20;N:ζ1s=6.66,ζ2s=1.91,ζ2p=1.89;H:ζ1s=1.24)对所研究的化合物进行了分子轨道计算,结合计算结果对化合物的紫外光电子能谱进行了解析和指认。实验得到的体现分子内弱相互作用σN-Ti的轨道电子的电离能表明了N,Ti间相互作用的弱成键特性,电子密度拓扑分析显示N,Ti原子间存在体现成键作用的键鞍点,键鞍点处电子密度的拓扑性质也进一步支持了N,Ti间弱成键作用的特性。     

H2CO与双卤分子间卤键的电子密度分析

运用B3LYP和MP2方法在6-311++G(d,p)基组水平上, 对H2CO-XY(XY=F2、Cl2、Br2、ClF、BrF、BrCl)卤键体系进行构型全优化, 得到了O…X—Y型卤键复合物. 结果表明, MP2/6-311++G(d,p)计算结果与实验值较吻合. 并在MP2水平下计算了分子间的相互作用能, 用完全均衡校正CP(counterpoise procedure)方法对基函数重叠误差(BSSE)进行了校正. 利用电子密度拓扑分析方法对卤键复合物的电子密度拓扑性质进行了分析研究.     

UO22+·nH2O和PuO22+·nH2O(n=2,4,5,6)密度泛函理论研究

首先用密度泛函理论(DFT)方法研究了铀酰和钚酰离子的几何与电子结构,计算结果与实验基本符合,表明DFT方法也能用于含铀和钚重原子的化合物计算.然后对铀酰和钚酰水合离子的几何构型、Mulliken集居数分布以及铀酰(钚酰)与配体水分子的结合能进行计算,计算结果表明UO22+·5H2O和PuO22+·5H2O分别为铀酰和钚酰系列水合离子中最稳定的配合物.     

磷酸二甲酯阴离子与二价金属离子相互作用的研究

在B3LYP/6-31G(d, p)水平上,用全电子分别研究了Mg2+、Ca2+、Mn2+、Co2+、Ni2+、Cu2+和Zn2+与磷酸二甲酯(DMP-)“双齿”配位模式下的相互作用,且由4个水分子((H2O)4)“饱和”各金属离子六配位模式下的其它位点。结果表明:这种配位模式下,DMP-对上述金属离子的选择性顺序为:Cu2+>Ni2+>Zn2+>Co2+>Mg2+>Mn2+>Ca2+,和Irving-Williams序列基本一致。且具有未充满d层轨道的金属对配体的变化“响应”明显。通过NBO(自然键轨道)方法对配合物电荷布居分析可得出主族金属离子与配体间以静电相互作用为主,非主族金属离子配合物中存在一定的共价作用的结论,同时关于配合物各部分间键轨道相互作用的分析部分的反映了金属离子与配体的作用实质。     

Au+(1S, 3D)与N2O(1Σ+)反应机理的理论研究

采用密度泛函B3LYP方法, O和N用6-311+G*基组, Au+用赝势基组(8s7p6d)/[6s5p3d], 研究了Au+(1S, 3D)离子和N2O(1Σ+)分子的反应机理. 报道了在基态单重态和激发三重态势能面上各反应物、中间体和过渡态的构型特征及能量. 结果表明, 两个主反应通道Au+(1S)+ N2O(1Σ+)→1NA-Complex-1→1NA-TS1→1NA-Complex-2→1NA-Crossing→[3OAuNN]+和Au+(1S)+ N2O(1Σ+)→1NB-Complex→1NB-Crossing→[AuNN(1Σ+)]++O(3P)都需经过反应交叉势能面, 出现“系间窜越”. 用内禀坐标单点计算垂直激发态的方法确定了势能面交叉点, 并用含时密度泛函TD-B3LYP方法进一步探讨了自旋翻转机理.     

Isomerization barriers in bis(4H-thiopyran) and in bithioxanthenes

The inversion and rotation mechanisms for the isomerization of Feringa’s bithioxanthenes existing in two conformations, up/up and up/down, have been calculated at the B3LYP/6-31G(d) and B3LYP/6-311++G(d,p) levels. The inversion mechanism that maintains the double bond nature of the central bond is a classical one but the rotation mechanisms that break the double bond to form a biradical needs to explore the singlet and triplet states. To do this we have removed the four fused phenyl rings of bithioxanthene and calculated at the CASSCF and CASPT2 levels bis(4H-thiopyran) proving that B3LYP calculations yield reasonable results for the rotation barriers.     

哒嗪裂解反应机理的密度泛函理论研究

利用密度泛函理论(DFT)方法, 以哒嗪为煤的模型化合物, 在(U)B3LYP/6-31G(d,p)水平下计算了标题化合物的键布居数和键裂解能, 并对其热解机理进行了探讨. 在(U)B3LYP/6-311++G(d,p)水平下计算了热解过程中各物种的单点能, 并对总能量进行了校正. 结果表明, 标题物N-N键的布居数和裂解能均最小, 且各键裂解能随温度变化不大, 热解最终产物为HCN和乙炔, NH₃可能是HCN经过二次转换得到的. 采用过渡态理论计算了300~1900 K温度范围内热解的速率常数, 求得lnk与1/T的线性关系.     

Exploring the actinide-actinide bond: theoretical studies of the chemical bond in Ac2, Th2, Pa2, and U2

Multiconfigurational quantum chemical methods (CASSCF/CASPT2) have been used to study the chemical bond in the actinide diatoms Ac2, Th2, Pa2, and U2. Scalar relativistic effects and spin-orbit coupling have been included in the calculations. In the Ac2 and Th2 diatoms the atomic 6d, 7s, and 7p orbitals are the significant contributors to the bond, while for the two heavier diatoms, the 5f orbitals become increasingly important. Ac2 is characterized by a double bond with a 3Sigmag-(0g+) ground state, a bond distance of 3.64. A, and a bond energy of 1.19 eV. Th2 has quadruple bond character with a 3Dg(1g) ground state. The bond distance is 2.76 A and the bond energy (D0) 3.28 eV. Pa2 is characterized by a quintuple bond with a 3Sigmag-(0g+) ground state. The bond distance is 2.37 A and the bond energy 4.00 eV. The uranium diatom has also a quintuple bond with a 7Og (8g) ground state, a bond distance of 2.43 A, and a bond energy of 1.15 eV. It is concluded that the strongest bound actinide diatom is Pa2, characterized by a well-developed quintuple bond.     

Effects of Epoxidation and Nitration on Ballistic Properties of FOX‐7 – A DFT Study

1, 1‐Diamino‐2, 2‐dinitroethylene (FOX‐7) has received increasing attention since it was industrialized in the late 1990s. It has lower sensitivity and comparable performance to RDX. This paper presents ballistic properties of FOX‐7, its mono and dinitro derivatives and their epoxide derivatives computationally. The structures were optimized at the B3LYP/6‐31G(d, p) level and the bond lengths were calculated. The calculated data for FOX‐7 are compatible with the literature one. We have investigated the bond dissociation energies of the molecules. Mulliken electro negativities (χ_M) and chemical hardness (η) were reviewed using Frontier Molecular Orbitals at HF/6‐31G(d, p)//B3LYP/6‐31G(d, p) theoretical level. The detonation performance analyses were done using empirical Kamlet‐Jacobs equations. Additionally, power index values were calculated. All the compounds considered in the present article are powerful candidates for high energy materials.     

Reactions of uranium atoms with ammonia: infrared spectra and quasi-relativistic calculations of the U:NH3, H2N--UH, and HN==UH2 complexes

Ammonia molecules interact with U atoms, and the resulting U:NH3 complex rearranges upon visible irradiation to form the H2N--UH and HN==UH2 molecules in excess argon. These products are identified by functional group frequencies, 15NH3 and ND3 isotopic shifts, and comparison to frequencies calculated by using density functional theory. The N==U pi bond in HN==UH2 is enhanced by partial triple-bond character through N(2p) to U(5f) conjugation, which is comparable to that found in the analogous HN==ThH2 molecule. These products also form complexes with additional ammonia molecules in the matrix. The interesting higher-energy N[triple chemical bond]UH3 complex is not formed.     

Relativistic bond lengthening of UO 2 2+ and UO2

Summary Relativistic calculations on UO₂ [1] have shown that relativity leads to substantial bondlengthening in this compound, in contrast to the bond contraction found almost exclusively for other compounds. The bond lengthening isnot caused by the relativistic expansion of the 5f valence AO of U, which is the primary bond forming orbital on U in UO₂. The origin of the bond lengthening can be traced back to the semi-core resp. subvalence character of the U 6p AO. The valence character of 6p shows up in an increasing depopulation of the 6p upon bond shortening, and hence loss of mass-velocity stabilization. The core character of 6p shows up in large off-diagonal mass-velocity matrix elements 5p|h _MV|6p which are shown to have an overall bond lengthening effect. The larger expansion in UO₂ than in UO ₂ ²⁺ is due to destabilization of U levels in UO₂, caused by repulsion of the two additional 5f electrons.The present analysis corroborates the picture of relativistic bond length effects of Ref. [2].     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

Why are hexavalent uranium cyanides rare while U–F and U–O bonds are common and short?

Relativistic small-core pseudopotential B3LYP and CCSD(T) calculations and frozen-core PW91–PW91 studies are reported for the series UF ₄X ₂ ( X=H, F, Cl, CN, NC, NCO, OCN, NCS and SCN). The bonding in UF ₆ is analyzed and found to have some multiple-bond character, approaching at a theoretical limit a bond order of 1.5. In addition to these s and p orbital interactions, the electrostatic attraction is important. Evidence for p bonding in the other systems studied was also found. The triatomic pseudohalides as well as fluorine and chlorine are in this sense better ligands than cyanide. The –CN group is a s donor and p acceptor, as uranium itself, and hence is unfit to bond to U(VI). The s-bonded UH ₆ is octahedral.     

Theoretical study of the gas-phase thermolysis reaction of alkyl (ethyl,isopropyl, and tert-butyl) N,N-dimethylcarbamates and N,N-diethylcarbamates

Theoretical studies on the thermolysis in the gas phase of alkyl N,N-dialkylcarbamates were carried out using ab initio theoretical methods, at the MP2/6-31G(d), MP2/6-31++G(d,p) and MP2/6-311++G(2d,p)//MP2/6-31G(d) levels. The reactions have two steps: the first one corresponds to the formation of an alkene and a neutral dialkylcarbamic acid intermediate via a six-membered cyclic transition state; the second one is the decarboxylation of this intermediate via a four-membered cyclic transition state, leading to carbon dioxide and the corresponding dialkylamine. The progress of the reactions was followed by means of the Wiberg bond indices. The results indicate that the transition states have character intermediate between reactants and products, and the calculated synchronicities show that the reactions are slightly asynchronous. The bond-breaking processes are more advanced than the bond-forming ones, indicating a bond deficiency in the transition states. The rate constants calculated for all the reactions agree very well with the available experimental data.From the Proceedings of the 28th Congreso de Químicos Teóricos de Expresión Latina (QUITEL 2002)