Bonding Energy of a Molecular Orbital(Ⅱ)——Ab initio Calculation

Calculation of the bonding energy of a molecular orbital for a series of small molecules has been carried out by using ab initio STO-3G method. The results obtained demonstrate that the concept of the molecular orbital bonding energy is applicable for judging whether a molecular orbital is bonding, nonbonding or antibonding besides Mulliken overlap criterion.

First-principles investigation of cohesive energy and electronic structure in vanadium phosphides

First-principles calculations based on the density-functional theory were employed to study the crystal structure of vanadium phosphide compounds,such as V3P,V2P,VP,VP2 and VP4. Cohesive energy of five types of vanadium phosphide compounds was calculated to assess their structural stability. The charge density distribution and densities of states of vanadium phosphides were discussed to study further their electronic structures. The results show that the structure of metal-rich compounds is considerably more stable than the phosphorus-rich compositions,and covalent bond exists between the V and P atoms of V3P,V2P,VP,VP2 and VP4.

Surface energetic and bonding characteristics of tetrahexahedral platinum nanocrystals enclosed by high-index facets

This paper uses a molecular static approach with a many-body potential to investigate the surface energetic and bonding characteristics of tetrahexahedral platinum nanocrystals enclosed by high-index facets such as {210}, {310}, {410}, {520} and {730}. It mainly focuses on the effect of crystal size and surface Miller index on these characteristics. The results show that the surface energy and dangling bond density increase with decreasing diameter of tetrahexahedral nanocrystals and generally show an order of {210} 〉{730}〉 {520} 〉 {310} 〉 {410}. However, this order is not valid at crystal sizes below 7 nm or so. The results of corresponding surfaces are also presented for comparison.

Bonding Energies and Structural Study of Alkylaluminium

Neutral aluminium alkyls are well known to act as ethylene oligomerization and polymerization catalysts and cocatalysts.On the basis of the full optimization of alkylaluminium compounds with Gaussian 98 program package at the B3LYP/6-31G＊＊ level,the selected structures and bonding energies were investigated extensively.The geometries and bonding energies of AlR3（R = H,CH3,C2H5,C3H7,C4H9） and Al（C2H5）2R＇（R＇ = C2H5,C3H7,C4H9,C5H11,C6H13） were investigated extensively,and we found that,along with the prolongation of carbon chains the terminal C-C bond is shortened gradually until to a constant value of about 0.1532 nm in C4H9;and the bonding energy almost remains constant.The dative bonding of C2H4 to Al（C2H5）3,whose bonding energy is only 15.30 kJ/mol,is very weak.

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 nitroalkane 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.mo1^-1 （for the BDE of tC4Hg-NOz）, 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 mo1^-1.

Thermal Analysis of Thermophysical Data for Equilibrium Pure Fluids

The thermal analysis of precise thermophysical data for pure fluids from electronic databases is developed to investigate the molecular interaction mechanisms and parameters and the structural features of heterogeneities in fluids. The method is based on the series expansion of thermophysical values by powers of the monomer fraction density. Unlike the virial expansion by powers of the total density, the series expansion terms in this method directly reflect properties of the corresponding cluster fractions. The internal energy had been selected among thermophysical properties as the most informative for this method. The thermal analysis of its series expansion coefficients permits to estimate the temperature dependence of the pair bond parameters, the clusters’ bond energies and equilibrium constants, the structural transitions between dominating isomers of clusters. The application of method to different pure fluids, including noble and molecular gases with van der Waals and polar molecular interactions, brings unknown clusters’ characteristics for the fluids under investigation. The thermal analysis of the ordinary and heavy Water vapors points on no trivial isotopic effects. The unpredictable growth of the pair bond energy with temperature in Alkanes points on existence in hydrocarbons of some unknown molecular interaction forces in addition to dispersion forces.

Hydrogen Bonds of C=S, C=Se and C=Te with C-H in Small-Organic Molecule Compounds Derived from the Cambridge Structural Database (CSD)

Considerable interest in hydrogen bonding involving chalcogen has been growing since the IUPAC committee has redefined hydrogen bonding. Not only the focus is on unconventional acceptors, but also on donors not discussed before. It has been mentioned in previous studies that the proton of the H-C group could be involved in hydrogen bonding, but with conventional acceptors. In this study, we explored the ability of hydrogen bond formation of Se, S and Te acceptors with the H-C donor using Cambridge Structural Database in conjunction with Ab Initio calculations. In the CSD, there are respectively 256, 6249 and 11 R1,R2,-C=Se, R1,R2,-C=S and R1,R2,-C=Te structures that form hydrogen bonds, in which the N,N groups are majority. Except for C=S acceptor which can form a hydrogen bond with its C, C group, both C=Se and C=Te acceptors could form a hydrogen bond only with N,C and N,N groups. CSD analysis shows very similar d (norm) around -0.04 Å, while DFT-calculated interaction for N,C and N,N groups are also similar. Both interaction distances derived from CSD analysis and DFT-calculated interaction energies demonstrate that the acceptors form stable complexes with H-CF3. Besides hydrogen bonds, dispersion interactions are forces stabilizing the complexes since their contribution can reach 50%. Analysis of intra-molecular geometries and Ab Initio partial charges show that this bonding stems from resonance induced Cδ+=Xδ- dipoles. In many respects, both C=Se, C=S and C=Te are similar to C=S, with similar d (norm) and calculated interaction strengths.

Theoretical Insights into the Inhibition Performance of Three Neonicotine Derivatives as Novel Type of Inhibitors on Carbon Steel

The adsorption process of new nicotinic derivatives on Fe(110)surface and diffusion of corrosive particles in inhibition film were studied by quantum chemistry and molecular dynamics simulation,and inhibition mechanism of inhibitor was discussed.The results indicated that the main active sites of three inhibitors are located in N atoms on the five membered ring.The inhibitor YM-1 has the strongest activity of electrophilic reaction,and the adsorption process of inhibitor molecules is polycentric chemisorption.The adsorption energy for inhibitors followed the order of YM-1>YM-2>YM-3.The adsorption film YM-1 more effectively impedes the diffusion and migration of the corrosive particles,and also has the best inhibitory effect on Cl-.In addition,the bond energy of chemical bond formed between inhibitor and Fe atom increases as the increase of temperature.The interaction between corrosive particles and inhibitors is dominated by Van der Waals force.

A Relativistic Density Functional Study of the U₂F₆Molecule

All-electronic relativistic density functional theory (DFT) method has been used to study the U2F6 molecule. Results from calculations predict the existence of U2F6 molecule, which has been found to be stable with a multiply bonded U2 unit. The calculations also predict that D3d symmetry of U2F6 is more stable than D3h. The optimized geometries, vibrational frequencies and infrared intensities are reported for D3d symmetry of U2F6 from Becke Three-parameter Lee-Yang-Parr (B3LYP) function with triple zeta valence plus polarization functions basis set (TZP). The bond dissociation energy (BDE) for U-U bond in the U2F6 was obtained using the same method. In addition, the entropies of U2F6 have been investigated at temperature rang from 0 to 3000K in 10 steps using the B3LYP method.

镁合金表面DCPD涂层的制备及其界面结合机制研究

目的研究CaHPO_(4)·2H_(2)O(DCPD)与Mg的界面结合机制,以提高DCPD在镁合金表面的界面结合强度。方法利用电镀法在AZ31镁合金表面制备DCPD涂层,采用SEM、XRD、XPS等对涂层形貌及结构进行表征。同时,运用分子动力学模拟(MD)对DCPD在Mg表面形成机制进行研究,通过统计界面层中不同组分的径向分布函数、密度分布、均力势、总能量等的变化,揭示DCPD/Mg的界面结合能、结合位点及结合方式。结果通过电镀法形成的DCPD涂层形貌为致密的荷花瓣状晶体,主要成分为CaHPO_(4)·2H_(2)O。模拟结果表明,CaHPO_(4)·2H_(2)O的4个晶面(010)、(-120)、(11-1)、(111)、(-120)与Mg的结合能最强(163.63 kJ/mol)。其中起“铆钉”作用的基团是HPO_(4)^(2-)和H_(2)O,结合位点主要为O与Mg,即HPO_(4)^(2-)和H_(2)O通过静电作用及范德华力与Mg形成Mg-HPO_(4)^(2-)和Mg-H_(2)O偶极对。研究发现,形成的偶极对中HPO_(4)^(2-)及H_(2)O的配位数分别为0.75和1.16,Mg-H_(2)O的解离能更大,结构更稳定。结论提出改善DCPD/Mg结合强度的方法,电镀前可将镁合金置于NH4H_(2)PO_(4)溶液中浸泡片刻,促进CaHPO_(4)·2H_(2)O(-120)晶面的形成。

Molecular dynamics study on the relationships of modeling,structural and energy properties with sensitivity for RDX-based PBXs

In this paper,a primary model is established for MD(molecular dynamics) simulation for the PBXs(polymer-bonded explosives) with RDX(cyclotrimethylene trinitramine) as base explosive and PS as polymer binder.A series of results from the MD simulation are compared between two PBX models,which are represented by PBX1 and PBX2,respectively,including one PS molecular chain having 46 repeating units and two PS molecular chains with each having 23 repeating units.It has been found that their structural,interaction energy and mechanical properties are basically consistent between the two models.A systematic MD study for the PBX2 is performed under NPT conditions at five different temperatures,i.e.,195 K,245 K,295 K,345 K,and 395 K.We have found that with the temperature increase,the maximum bond length(L max) of RDX N N trigger bond increases,and the interaction energy(E N-N) between two N atoms of the N-N trigger bond and the cohesive energy density(CED) decrease.These phenomena agree with the experimental fact that the PBX becomes more sensitive as the temperature increases.Therefore,we propose to use the maximum bond length L max of the trigger bond of the easily decomposed and exploded component and the interaction energy E N-N of the two relevant atoms as theoretical criteria to judge or predict the relative degree of heat and impact sensitivity for the energetic composites such as PBXs and solid propellants.

氢键调控提升聚硫脲介电储能特性

为理解氢键对聚硫脲(polythiourea,PTU)介电储能特性的影响,该文采用实验和分子模拟方法,研究脂环族PTU(AcPTU)和芳香族PTU(ArPTU)中氢键作用模式和强度对介电常数、击穿场强和储能特性的影响。利用红外光谱和核磁共振氢谱表征化学组成和分子构型,利用宽频介电谱、击穿和电滞回线测试获得介电储能特性,基于H-N公式提取偶极极化特征参数,利用热刺激去极化电流法提取了陷阱深度,采用Multiwfn软件识别氢键作用模式和强度。结果显示,氢键作用模式和强度可以通过改变偶极极化和陷阱特性影响介电性能。脂环结构可以增多由反式/反式构型硫脲组装成的双氢键阵列,增强其介电松弛强度,进而提高介电常数;苯环结构可以增强π氢键数目和作用强度,加深由苯环碳原子和硫原子等氢键受体贡献的深陷阱能级,进而提高击穿场强。AcPTU具有高介电常数,归因于脂环结构增强了双氢键阵列的介电松弛强度;ArPTU的击穿场强高达855 MV/m,是因为苯环结构强化了氢键的电荷捕获效应。此外,当电场强度为300 MV/m时,AcPTU的储能密度和效率分别达到2.62 J/cm3和93.31%。该研究为提升本征型聚合物介电储能特性提供了理论依据。

Simulative Calculation of Mechanical Property, Binding Energy and Detonation Property of TATB/Fluorine-polymer PBX

Molecular dynamics （MD） method was used to simulate 1,3,5-triamino-2,4,6-trinitrobenzene （TATB） coated with fluorine containing polymers. The mechanical properties and binding energies of PBXs were obtained. It was found that when the number of chain monomers of fluorine containing polymers was the same, the elasticity of TATB/F2314 was increased more greatly than others and the binding energy of TATB/F2311 was the largest among four PBXs. Detonation heat and velocity of such four PBXs were calculated according to theoretical and empirical formulas. The results show that the order of detonation heat is TATB〉TATB/PVDF〉TATB/F2311〉TATB/ F2314 〉 TATB/PCTFE while the order of detonation velocity is TATB/PVDF 〈 TATB/F2311 〈 TATB/F2314 〈 TATB/PCTFE 〈TATB.

Simulation investigations in the binding energy and mechanical properties of HMX-based polymer-bonded explosives

The molecular simulations of the well-known high explosive β-HMX (cyclotetramethylene tetranitramine) anditsfluorine containing polymer-bonded explosives(PBXs) were carried out with the combination method of quantum mechanics, molecular mechanics and molecular dynamics. The atomic cluster model, containing the β-HMX molecule and the polymer molecule whose chain dimension was about the same as β-HMX’s, was fully optimized by AM1 and PM3 semi-empirical molecular orbital and molecular mechanical methods using COMPASS and PCFF force field. Then the calculated binding energy is found to be linearly correlated to each other. Molecular dynamics simulations using COM- PASS force field were performed for β-HMX crystal and the PBXs involving β-HMX and a series of fluorine containing polymers. Their elastic coefficients, moduli and Poisson’s ratios were calculated. It is found that the mechanical prop- erties of β-HMX can be effectively improved by blending with fluorine containing polymers in small amounts.

Relativistic density functional investigation of UX_6(X=F,Cl,Br and I)

The molecular structures and the vibrational frequencies of uranium hexahalides UX6 （X=F, Cl, Br and I） molecules are investigated by using local density approximation （LDA） and generalised gradient approximation （GGA） functions （BP, BLYP and RPBE） in combination with two different relativistic methods （scalar and scalar＋spin-orbit relativistic effects）. The calculated results show that the differences are trivial between scalar and scalar＋spin-orbit relativistic methods. The vibrational frequencies are also compared with existing experimental values, and overall, the RPBE approach gives the smallest error. The bond dissociation energies （BDEs） of UX6 are computed by using the RPBE function, thereby obtaining exact vibrational frequencies. In addition, the calculated magnitudes of the spin orbit effect on the BDE of UX6 （X=F, Cl, Br, and I） are found to be approximately -0.3198, 0.3218, -0.3609 and -0.4415 eV, respectively.

XeO_(3)和XeOF_(2)与三种不同杂化方式N原子的Xe…N非共价作用的理论研究

本文采用密度泛函理论和二级微扰理论方法研究了XeO_(3)和XeOF_(2)与一系列不同杂化含氮给体配合物之间的相互作用.含氮给体与XeO_(3)和XeOF_(2)之间的相互作用能在6.5~19.9 kcal/mol之间.相互作用强弱顺序为sp^(3)>sp^(2)>sp,且XeO_(3)比XeOF_(2)高.约化密度梯度图证明,对于sp^(2)和sp^(3)杂化的给体,空间位阻效应在相互作用中起着一定的影响.能量成分分析发现,Xe……N作用主要的吸引部分是静电相互作用.在XeO_(3)配合物中,极化能的占比大于色散,而在XeOF_(2)配合物中则相反.除表面静电势(V_(s))外,其他5个相互作用参数之间的线性相关系数均较大.

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5-Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1 ： 1 5-hydroxytryptamine-water （5-HT-H20） complexes are studied at the MP2/6-311 ＋ ＋ G（d,p） level. Natural bond orbital （NBO）, quantum theory of atoms in molecules （QTAIM） analyses and the localized molecular orbital energy decomposition analysis （LMO-EDA） were performed to explore the nature of the hydrogen-bonding interactions in these complexes. Various types of hydro- gen bonds （H-bonds） are formed in these 5-HT-H20 complexes. The intermolecular C4H55HT＇＂Ow H-bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5-HT-H20 complexes. H-bond in which nitrogen atom of amino in 5-HT acted as proton donors was stronger than other H-bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5-HT-H20 complexes.

Deconstructing Vibrational Motions on the Potential Energy Surfaces of Hydrogen-Bonded Complexes

Internal vibrations underlie transient structure formation,spectroscopy,and dynamics.However,at least two challenges exist when aiming to elucidate the contributions of vibrational motions on the potential energy surfaces.One is the acquisition of well-resolved experimental infrared spectra,and the other is the development of efficient theoretical methodologies that reliably predict band positions,relative intensities,and substructures.Here,we report size-specific infrared spectra of ammonia clusters to address these two challenges.Unprecedented agreement between experiment and state-of-the-art quantum simulations reveals that the vibrational spectra are mainly contributed by proton-donor ammonia.A striking Fermi resonance observed at approximately 3210 and 3250 cm^(−1)originates from the coupling of NH symmetric stretch fundamentals with overtones of free and hydrogen-bonded NH bending,respectively.These novel,intriguing findings contribute to a better understanding of vibrational motions in a large variety of hydrogen-bonded complexes with orders of magnitude improvements in spectral resolution,efficiency,and specificity.

Chalcogen Bonds in Small-Organic Molecule Compounds Derived from the Cambridge Structural Database (CSD)

Growing interest in non-covalent interactions involving chalcogen atoms has been ascribed to their importance in crystal engineering, molecular recognition and macromolecular edifices. The present study is dealing with chalcogen bonds involving divalent Sulphur, Selenium and Tellurium atoms, acting as sigma-hole donors, in small-molecule compounds using the Cambridge Structural Database (CSD) in conjunction with ab initio calculations. Results derived from CSD surveys and computational study revealed that nucleophiles formed complexes with the chalcogen-bond donors R1-X-R2 (X = S, Se or Te). The main forces stabilizing the complexes were chalcogen bonds, enhanced by dispersion interactions. Complexation pattern and energetics show that nucleophile bonding at divalent S, Se and Te atoms is a relatively strong and directed interaction. The bond consists of a charge transfer from a nucleophile atom lone pair to an X-R1 or X-R2 antibonding orbital.

双原子分子键长对振动热容量的影响

在低温下双原子分子的振动自由度对气体的热容量没有贡献,当温度高于双原子分子振动特征温度的时候,振动自由度才被激发,这个双原子分子的振动是对原子间相互作用的简谐近似。不过,文献中给出的简谐近似判据,都是力学结果,而且这些结果都是定性的。本研究给出统计物理判据:简谐近似的有效性是温度必须远低于分子键能特征温度。