Quantum chemical calculations of bond dissociation energies for COOH scission and electronic structure in some acids

Quantum chemical calculations are performed to investigate the equilibrium C-COOH bond distances and the bond dissociation energies(BDEs) for 15 acids.These compounds are studied by utilizing the hybrid density functional theory(DFT)(B3LYP,B3PW91,B3P86,PBE1PBE) and the complete basis set(CBS-Q) method in conjunction with the 6311G** basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work.Comparisons between the computational results and the experimental values reveal that CBS-Q method,which can produce reasonable BDEs for some systems in our previous work,seems unable to predict accurate BDEs here.However,the B3P86 calculated results accord very well with the experimental values,within an average absolute error of 2.3 kcal/mol.Thus,B3P86 method is suitable for computing the reliable BDEs of C-COOH bond for carboxylic acid compounds.In addition,the energy gaps between the highest occupied molecular orbital(HOMO) and the lowest unoccupied molecular orbital(LUMO) of studied compounds are estimated,based on which the relative thermal stabilities of the studied acids are also discussed.     

Structures and electronic properties of SimN8-m（0 〈 m 〈 8） clusters： a density functional theory study

The geometries, electronic structures and related properties of SimN8-m（0 〈 m 〈 8） clusters are studied using density functional theory （DFT） with hybrid functional B3LYP. The calculated results reveal several trends. For any stoichiometric clusters, the lowest energy isomers with an alteration of N and Si atoms are favourable in energy if the numbers of Si and N atoms are large enough to form ... Si N-Si-N... alternative chains. The bond lengths of single Si-N bonds are very close to the corresponding values of the bulk and other SiN clusters. The geometries for N-rich and Si4N4 clusters are planar structures, but three-dimensional structures are favourable in energy for Si-rich clusters. With the increase of m, the isotropic polarizability and average polarizability increase, the total binding energies generally decrease, the HOMO-LUMO gap and vertical ionization potential oscillate with increasing number of valence electrons, and their values with even valence electrons are larger than those with odd valence electrons. The atomic charges, IR and Raman properties are also reported.     

Molecular dynamics simulations of displacement cascades in Fe-10%Cr systems

Molecular dynamics simulations of the displacement cascades in Fe-10%Cr systems are used to simulate the primary knocked-on atom events of the irradiation damage at temperatures 300,600,and 750 K with primary knockedon atom energies between 1 and 15 keV.The results indicate that the vacancies produced by the cascade are all in the central region of the displacement cascade.During the cascade,all recoil Fe and Cr atoms combine with each other to form Fe-Cr or Fe-Fe interstitial dumbbells as well as interstitial clusters.The number and the size of interstitial clusters increase with the energy of the primary knocked-on atom and the temperature.A few large clusters consist of a large number of Fe interstitials with a few Cr atoms,the rest are Fe-Cr clusters with small and medium sizes.The interstitial dumbbells of Fe-Fe and Fe-Cr are in the 111 and 110 series directions,respectively.     

Theoretical study of molecular hydrogen and spiltover hydrogen storage on two-dimensional covalent-organic frameworks

Molecular hydrogen and spiltover hydrogen storages on five two-dimensional （2D） covalent-organic frameworks （COFs） （PPy-COF, TP-COF, BTP-COF, COF-18 A, and HHTP-DPB COF） are investigated using the grand canonical Monte Carlo （GCMC） simulations and the density functional theory （DFT）, respectively. The GCMC simulated results show that HHTP-DPB COF has the best performance for hydrogen storage, followed by BTP-COF, TP-COF, COF-18 A, and PPy-COE However, their adsorption amounts at room temperature are all too low to meet the uptake target set by US Department of Energy （US-DOE） and enable practical applications. The effects of pore size, surface area, and isosteric heat of hydrogen on adsorption amount are considered, which indicate that these three factors are all the important factors for determining the H2 adsorption amount. The chemisorptions of spiltover hydrogen atoms on these five COFs represented by the cluster models are investigated using the DFT method. The saturation cluster models are constructed by considering all possible adsorption sites for these cluster models. The average binding energy of a hydrogen atom and the saturation hydrogen storage density are calculated. The large average binding energy indicates that the spillover process may pro- ceed smoothly and reversibly. The saturation hydrogen storage density is much larger than the physisorption uptake of H2 molecules at 298 K and 100 bar （1 bar = 105 Pa）, and is close to or exceeds the 2010 US-DOE target of 6 wt% for hydrogen storage. This suggests that the hydrogen storage capacities of these COFs by spillover may be significantly enhanced. Thus 2D COFs studied in this paper are suitable hydrogen storage media by spillover.     

Assignment of Photoelectron Spectra of MC2 （M= V, Cr, Fe, and Co）

Hybrid density functional theory （DFT） calculations are performed to study MC2 （M= V, Cr, Fe and Co） clusters in the neutral and anionic charge states. We find that the equilibrium geometries of MC2 and their anions are all cyclic structures with C2v symmetry, which agrees well with the previous theoretical studies. The Mulliken charge and spin populations of MC2 clusters and their anions are also calculated, and it is found that the electron charge transformations from anions to neutral molecules mainly take place on the M atoms. Time-dependent DFT is used to calculate the excited states, and a theoretical assignment for the features in the experimental photoelectron spectrum is given, which are in good agreement with the available experimental data.     

Lattice Dynamical Simulation of Guest-Host Interaction in N2 Clathrate Hydrate

We perform the lattice dynamical simulation studies of hydrate host lattice interacting with Xe, Ar, and N2 atoms/molecules. The calculated results show that the well-defined peaks （2.0meV and 3.8meV） and another peak （6.2 me V） are assigned to the vibrations of N2 molecules in large and small cages, respectively. It is confirmed that the double N2 molecule occupancies of large cage lead to filling of the mode gap between the small cage and the large cage.     

Influence of Ni on Cu precipitation in Fe Cu Ni ternary alloy by an atomic study

The early aging Cu precipitations in Fe-3%Cu and Fe-3%Cu-4%Ni ternary alloys are investigated by molecular dynamics （MD） simulations. The results show that the average size of Cu clusters in Fe-3%Cu-4%Ni alloy is larger than that in Fe-3%Cu alloy. The diffusion of Cu is accelerated by Ni according to the mean square displacement （MSD）. Furthermore, the whole formation process of Cu-rich clusters is analyzed in detail, and it is found that the presence of Ni promotes small Cu-rich clusters to be combined into big ones. Ni atoms prefer to stay at the combination positions of small clusters energetically due to a large number of the first nearest neighbor Cu-Ni interactions, which is verified by first-principles calculations based on density functional theory （DFT）.     

Molecular Dynamics Simulations of Helium Behaviour in Titanium Crystals

Molecular dynamics simulations are performed to investigate the behaviour of helium atoms in titanium at a temperature of 30OK. The nucleation and growth of helium bubble has been simulated up to 50 helium atoms. The approach to simulate the bubble growth is to add helium atoms one by one to the bubble and let the system evolve. The titanium cohesion is based on the tight binding scheme derived from the embedded atom method, and the helium-titanium interaction is characterized by fitted potential in the form of a Lennard-Jones function. The pressure in small helium bubbles is approximately calculated. The simulation results show that the pressure will decrease with the increasing bubble size, while increase with the increasing helium atoms. An analytic function about the quantitative relationship of the pressure with the bubble size and number of helium atoms is also fitted.     

Computation of Bond Dissociation Energies for Removal of Nitrogen Dioxide Groups in Certain Aliphatic Nitro Compounds

Bond dissociation energies for removal of nitrogen dioxide groups in 10 aliphatic nitro compounds, including nitromethane, nitroethylene, nitroethane, dinitromethane, 1-nitropropane, 2-nitropropane, 1-nitrobutane, 2-methyl-2-nitropropane, nitropentane, and nitrohexane, are calculated using the highly accurate complete basis set （CBS-Q） and the three hybrid density functional theory （DFT） methods B3LYP, B3PW91 and B3P86 with 6-31G^＊＊ basis set. By comparing the computed bond dissociation energies and experimental results, we find that the B3LYP/6-31G^＊＊ and B3PW91/6-31G^＊＊ methods are incapable of predicting the satisfactory bond dissociation energy （BDE）. However, B3P86/6-31G^＊＊ and CBS-Q computations are capable of giving the calculated BDEs, which are in extraordinary agreement with the experimental data. Nevertheless, since CBS-Q computational demands increase rapidly with the number of containing atoms in molecules, larger molecules soon become prohibitively expensive. Therefore, we suggest to take the B3P86/6-31G^＊＊ method as a reliable method of computing the BDEs for removal of the NO2 groups in the aliphatic nitro compounds.     

Structural, electronic and vibrational properties of indium oxide clusters

Geometric, electronic and vibrational properties of the most stable and energetically favourable configurations of indium oxide clusters In_mO_n (1≤m, n≤4) are investigated using density functional theory. The lowest energy geometries prefer the planar arrangement of the constituent atoms with a trend to maximize the number of ionic In-O bonds. Due to the charge transfer from In to O atoms, the electrostatic repulsion occurs between the atoms with the same kind of charge. The minimization of electrostatic repulsion and the maximization of In-O bond number compete between each other and determine the location of the isometric total energy. The most stable linear In-O-In-O structure of In₂O₂ cluster is attributed to the reduced electrostatic repulsive energy at the expense of In-O bond number, while the lowest energy rhombus-like structure of In₂O₃ cluster reflects the maximized number of In-O bonds. Furthermore, the vibrational frequencies of the lowest energy clusters are calculated and compared with the available experimental results. The energy gap and the charge density distribution for clusters with varying oxygen/indium ratio are also discussed.     

四碘甲状腺素团簇结构与光谱性质的密度泛函理论研究

在B3LYP/Lanl2mb基组水平上,利用密度泛函理论(DFT)优化了四碘甲状腺素团簇的几何结构.基于该团簇的几何结构下,其吸收和发射光谱的研究使用相同的基组水平并采用极化连续介质模型(PCM)下用含时密度泛函理论(TDDFT).研究结果表明,优化所得甲状腺素团簇的几何结构对称性为C_1;在基态稳定结构基础上,得出其输运性质,即甲状腺素团簇为p型输运材料;通过含时密度泛函理论,在优化好的基态结构基础上,又计算了它的溶剂效应,进一步得出该分子在水溶剂中的吸收光谱和发射光谱特性.     

Density Functional Approach Based on Numerically Obtained Bridge Functional

The Ornstein Zernike equation is solved with the Rogers Young approximation for bulk hard sphere fluidand Lennard-Jones fluid for several state points. Then the resulted bulk fluid radial distribution function combinedwith the test particle method is employed to determine numerically the function relationship of bridge functional as afunction of indirect correlation function. It is found that all of the calculated points from different phase space statepoints for a same type of fluid collapse onto a same smooth curve. Then the numerically obtained curve is used tosubstitute the analytic expression of the bridge functional as a function of indirect correlation function required in themethodology [J. Chem. Phys. 112 (2000) 8079] to deterrnine the density distribution of non-uniform hard spherefluid and Lennard Jones fluid. The good agreement of theoretical predictions with the computer simulation data isobtained. The present numerical procedure incorporates the knowledge of bulk fluid radial distribution function intothe constructing of the density functional approximation and makes the original methodology more accurate and moreflexible for various interaction potential fluid.     

三碘甲状腺素团簇结构与光谱性质的密度泛函理论研究

我们利用密度泛函理论(DFT)，在B3LYP/Lan12mb基组水平上，得到了三碘甲状腺素团簇的几何和电子结构.在此基础上，利用含时密度泛函理论(TDDFT)，使用相同的基组和采用极化连续介质模型(PCM)，对其溶剂效应下的吸收光谱进行研究．研究结果表明，优化所得三碘甲状腺素团簇的对称性为C₁;在基态稳定结构基础上，研究了该分子的红外和拉曼分子振动谱特性，同时研究了其输运性质，即三碘甲状腺素团簇为p型输运材料;通过含时密度泛函理论，在优化好的基态结构基础上，又计算了它的溶剂效应，进一步得出该分子在水溶剂中的吸收光谱特性.     

Quantum Computing without Magic,by Zdzislaw Meglicki

We will discuss the key concepts in density functional theory (DFT), how it can be used to model experimental data, and consider how the synergy between DFT and experiment can give significant insights. The discussion will centre on the scanning tunnelling microscope (STM) and surface problems, tracking the author's personal interest, though the general principles are widely applicable.     

Augmented Kierlik--Rosinberg Fundamental Measure Functional and Extension of Fundamental Measure Functional to Inhomogeneous Non-hard Sphere Fluids

From point of view of weighted density procedure, it is guessed that a Percus-Yevick （PY） compressibility excess free energy density, appearing in the Kierlik Rosinberg type fundamental measure functional （KR-FMF） and expressed in terms of scaled particle variables, can be substituted by a corresponding expression dictated by a more accurate Mansoori Carnahan-Starling Leland （MCSL） equation of state, while retaining the original weighting functions; it is numerically indicated that the resultant undesirable non-self-consistency between the PY type weighting function and MCSL type excess free energy density had no bad effect on the performance of the resultant augmented KR-ffMF which, on the one hand, preserves the exact low-density limit of the original KR-FMF and holds a high degree of pressure self-consistency, on the other hand, improves significantly, as expected, the predictions of density profile of hard sphere fluid at single hard wall contact location and its vicinity, and of the bulk hard sphere second order direct correlation function （DCF）, obtained from functional differentiation. The FMF is made applicable to inhomogeneous non-hard sphere fluids by supplementing a functional perturbation expansion approximation truncated at the lowest order with summation of higher order terms beyond the lowest term calculated by the FMF for an effective hard sphere fluid; the resultant extended FMF only needs second order DCF and pressure of the fluid considered at coexistence state as inputs, consequently is applicable whether the considered temperature is above critical point or below critical point. The extended MCSL-augmented KR-FMF is found to be endowed with an excellent performance for predictions of density profile and surface tension by comparing the present predictions of these two quantities with available computer simulation data for inhomogeneous hard core attractive Yukawa fluid and Lennard-3ones fluid.     

密度泛函理论研究吲哚并咔唑同分异构体结构,芳香性和光谱性质

在B3LPY/6-31G(d, p)基组水平上,利用密度泛函理论(DFT)优化了吲哚并咔唑五种同分异构体的几何和电子结构.基于这五种同分异构体的几何结构下,其吸收和发射光谱的研究使用相同的基组水平并采用极化连续介质模型(PCM)下用含时密度泛函理论(TD-DFT)计算.由于三种近似线性分子(吲哚并[2,3-a]咔唑、吲哚并[2,3-b]咔唑和吲哚并[3,2-b]咔唑)的电荷转移跃迁的振荡强度较大,这些异构体的发射光谱存在明显差异;基于吲哚并[2,3-c]咔唑和吲哚并[3,2-a]咔唑的构型特征,这两种同分异构体的发射谱具有高能量.比较计算结果表明,吲哚并[2,3-b]咔唑在这些分子中的振荡强度最大.这是因为当吲哚并咔唑的五种同分异构体的结构从基态变为激发态时,这些分子的分子轨道(MO)能级不同.由计算结果还得出,这五个分子适用于P型传输材料,并且每个分子的三个苯环均具有共轭效应.     

Augmented Kierlik-Rosinberg Fundamental Measure Functional andExtension of Fundamental Measure Functional to Inhomogeneous Non-hard Sphere Fluids

From point of view of weighted density procedure, it isguessed that a Percus-Yevick (PY) compressibility excess free energydensity, appearing in the Kierlik--Rosinberg type fundamentalmeasure functional (KR-FMF) and expressed in terms of scaledparticle variables, can be substituted by a corresponding expressiondictated by a more accurate Mansoori-Carnahan-Starling-Leland(MCSL) equation of state, while retaining the original weightingfunctions; it is numerically indicated that the resultantundesirable non-self-consistency between the PY type weightingfunction and MCSL type excess free energy density had no bad effecton the performance of the resultant augmented KR-FMF which, on theone hand, preserves the exact low-density limit of the originalKR-FMF and holds a high degree of pressure self-consistency, on theother hand, improves significantly, as expected, the predictions ofdensity profile of hard sphere fluid at single hard wall contactlocation and its vicinity, and of the bulk hard sphere second orderdirect correlation function (DCF), obtained from functionaldifferentiation. The FMF is made applicable to inhomogeneousnon-hard sphere fluids by supplementing a functional perturbationexpansion approximation truncated at the lowest order with summationof higher order terms beyond the lowest term calculated by the FMFfor an effective hard sphere fluid; the resultant extended FMF onlyneeds second order DCF and pressure of the fluid considered atcoexistence state as inputs, consequently is applicable whether theconsidered temperature is above critical point or below criticalpoint. The extended MCSL-augmented KR-FMF is found to be endowedwith an excellent performance for predictions of density profile andsurface tension by comparing the present predictions of these twoquantities with available computer simulation data for inhomogeneoushard core attractive Yukawa fluid and Lennard-Jones fluid.     

Electronic Excitations in Pyrrole: A Test Case for Determination of Chromophores in the Chromogenic Effects of Neurotoxic Hydrocarbons by Time-Dependent Density Functional Theory and Single-Excitation Configuration Interaction Methods

Time-dependent density functional theory (TD-DFT) and single-excitation configuration interaction (CIS) calculations on the electronic excitations in pyrrole have been performed to examine the reliability of these first-principles electronic structure methods in predicting electronic excitation spectraof pyrrole-containing compounds. Both the TD-DFT and CIS calculations led to satisfactory results when compared to available experimental data, particularly for low-lying excited states. The TD-DFT and CIS calculations provide lower and upper limits of the excitation energies, respectively, for low-lying singlet excited states. These results suggest that these methods can be used for the prediction of the excitation spectra, particularly the excitation energies for low-lying excited states, of chromophores responsible for the chromogenic effects of neurotoxic hydrocarbons, which are believed to be substituted pyrroles and their adducts with proteins. As an example of a practical application, the spectrum of the widely used 2,5-dimethylpyrrole has been calculated. It is shown that the 2,5-dimethylpyrrole molecule does not have an absorption in the region of the visible spectrum (400-700 nm), suggesting that the absorption observed at 530 nm and the color of 2,5-dimethylpyrrole is due to another species, probably a product of possible 2,5-dimethylpyrrole autoxidation. This suggests that the conclusions from previously reported experimental studies of biochemical reactions of neurotoxic γ-diketones need to be reexamined in terms of the relationship of chromogenicity to neurotoxicity.     

Universality principle and the development of classical density functional theory

The universality principle of the free energy density functional and the ‘test particle' trick by Percus are combined to construct the approximate free energy density functional or its functional derivative. Information about the bulk fluid radial distribution function is integrated into the density functional approximation directly for the first time in the present methodology. The physical foundation of the present methodology also applies to the quantum density functional theory.     

团簇LaO的理论研究

采用密度泛函理论研究LaO团簇体系。中性分子LaO的基态是两重态(2Σ),阴离子LaO-和阳离子LaO 的基态都是单重态(1Σ)。使用不同的方法计算团簇LaO的电子亲和能和电离能。计算结果表明用BLYP方法和弥散极化基组计算结果和实验数据吻合较好。用含时密度泛函理论计算团簇LaO的低能激发态,从理论上归属LaO-的光电子能谱的谱峰和LaO的吸收光谱的谱峰。计算得到与实验一致的结果。