Divide-and-conquer calculations for clean surfaces and surface adsorption

The divide-and-conquer density-functional approach is implemented for electronic structure calculations of clean surfaces and surface adsorption. The method divides a semi-infinite solid-state system into subsystems and determines each subsystem density separately. We have applied the method to study two low index Li surfaces and the chemisorption of a monolayer of atomic hydrogen on the Li(100) surface. Surface and adsorption energies converge well with respect to the number of buffer atoms and surface layers. Our converged energies are in very good agreement with accurate theoretical and available experimental results. Received: 8 January 1997 / Accepted: 14 January 1997     

Grid-free DFT implementation of local and gradient-corrected XC functionals

Following an approach to density functional theory calculations based on the matrix representation of operators, we implemented a scheme as an alternative to traditional grid-based methods. These techniques allow integrals over exchange-correlation operators to be evaluated through matrix manipulations. Both local and gradient-corrected functionals can be treated in a similar way. After deriving all the required expressions, selected examples with various functionals are given. Received: 7 March 1998 / Accepted: 21 May 1998 / Published on line: 6 August 1998     

Geometry and tautomerism of sapphyrin—DFT studies

The structure and total energy have been investigated, applying the density functional theory (DFT), for idealized sapphyrin derivatives, created by replacements of meso-phenyl groups of 5,10,15,20-tetraphenylsapphyrin or β-alkyl groups of decaalkylsapphyrin with hydrogens or two methyl groups, respectively. The following species have been considered: the dication of plain sapphyrin (SapH₅²⁺) 12,13-dimethylsapphyrin (12,13-DMSapH₃), 12,13-dimethylsapphyrin dication (12,13-DMSapH₅²⁺), and the 10,15-dimethylsapphyrin dication (10,15-DMSapH₅²⁺). To address the issue of tautomerization the calculations have been carried out for the six feasible distribution patterns of NH protons among the five nitrogen centers of the planar sapphyrin P0_i (i=1–6). The B3LYP/6-31G optimized bond lengths and angles of sapphyrin or (dimethylated sapphyrin) skeletons are in satisfactory agreement with X-ray crystallographic values determined for decaalkylsapphyrin dications. The calculated total energies, using the B3LYP/6-31G**//3-21G approach demonstrate that the plain sapphyrin and 12,13-dimethylsapphyrin dications favor the planar geometry of the macrocycle. Localization of the methyl groups at the 10,15-meso positions of sapphyrin reversed the order of stability leading to the energy preference for the inverted structure. Thus, this steric factor related to the meso substitution seems to be instrumental in the C pyrrolic ring inversion. The DFT calculation also provided the relative stability frame for the complex tautomeric equilibria of SapH₃ which involve six different forms. Their relative stability decreases in the range: {25-NH, 26-N, 27-NH, 28-NH, 29-N} P0₃>({25-NH, 26-N, 27-NH, 28-N, 29-NH} P0₁>{25-NH, 26-NH, 27-N, 28-NH, 29-N} P0₅>{25-N, 26-NH, 27-NH, 28-NH, 29-N} P0₂>{25-NH, 26-N, 27-N, 28-NH, 29-NH} P0₆>{25-N, 26-N, 27-NH, 28-NH, 29-NH} P0₄. The stability order can be qualitatively related to the number and the nature of the cis NH interactions present in the sapphyrin core.     

A DFT analysis of the vibrational spectra of nitrobenzene

Raman and FTIR, spectra of nitrobenzene, nb, and its isotopomers, nb-¹⁵N, nb-¹³C₆ and nb-d₅, were obtained and the fundamental vibrational modes assigned with the aid of a B3LYP/6-311+G** calculation, without the need for scaling of the force constants. The changes in vibrational coupling between the nitro and benzene groups upon certain isotopic substitutions are well modelled by the calculation, which is able to reproduce the isotopic shifts in frequencies for the nitro vibrations, as well as changes in IR intensities.     

DFT Study of Benzofuroxan Synthesis Mechanism from 2-Nitroaniline via Sodium Hypochlorite

The oxidative cyclization reaction of 2-nitroaniline via sodium hypochlorite to yield benzo-furoxan is investigated by the hybrid density functional theory B3LYP/6-31G(d,p) method. Solvent effects are estimated with the polarizable continuum model to optimize structures. The title reaction is predicted to undergo two pathways, each of which is a stepwise process.Path A includes four steps, namely oxidization, H-attack, hydrolysis, and cyclization. Path B involves the nucleophilic attack of OH^- to the H atom of the N-H bond and the proton transfer to the N atom of amino group leading to the cleavage of the N-H single bond in the amino group. The calculated results indicate that path A is favored mechanism for the title reaction. Furthermore, it is rational for one water molecule serving as a bridge to assist in the hydrolysis step of Path A and our calculations exhibit that this process is the rate-determining step.     

A DFT variational approach to Hooke''s atom based on local-scaling transformations

The local-scaling transformation version of density functional theory, LS-DFT, is employed in order to construct energy functionals for Hooke's atom. The components of the energy are analyzed and the resulting exchange and correlation potentials are compared with the exact ones. In addition, the representation of the exact one-particle density in terms of the various components of the total energy density is discussed.     

基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

Capture and separation of CO2 on BC3 nanosheets: A DFT study

In order to reduce the greenhouse effect caused by the rapid increase of CO₂ concentration in the atmosphere, it is necessary to develop more efficient, controllable, and highly sensitive adsorbing materials. In this study, the adsorption behavior of CO₂ on BC₃ nanosheets under an external electric field was explored based on density functional theory (DFT). It was found that CO₂ experienced a transition from physisorption to chemisorption in the electric field range of 0.0060-0.0065 a.u.. In addition, the adsorption/desorption of CO₂ is reversible and can be precisely controlled by switching on/off at the electric field of 0.0065 a.u.. The selective adsorption of CO₂/H₂/CH₄ by BC₃ can also be used to realize gas separation and purification under different electric fields. This study highlighted the potential application of BC₃ nanosheets as a high-performance, controllable material for CO₂ capture, regeneration, and separation in an electric field.     

Approximation of the molecular electrostatic potential in a gaussian density functional method

The recently developed Asymptotic Density Model (ADM) [6, 9] is here implemented in the density functional framework using the program deMon-KS [13]. While the original implementation divided the atoms into a core shell and a valence shell, the present version allows for an arbitrary number of shells making it therefore more flexible and, as shown with benzene, potentially more accurate. Moreover, since this method is derived through Poisson's equation, an expression for the electronic charge density is also obtained. However, the present discussion will restrict itself to the electrostatic potential. Finally, even though this method requires parametrization, it is shown that the parameters obtained for homonuclear diatomic species, and used as is in molecular calculations, yield satisfactory results. Indeed, the ADM reproduces almost all basic features of the MEP for all molecules presented here, (water, ammonia, ethylene, acetylene, hydrogen cyanide, carbon monoxide, benzene, nitrous acid). Received: 5 July 1996 / Accepted: 12 November, 1996     

Reactivity of the carbon-carbon double bond towards nucleophilic additions. A DFT analysis

The global and local electrophilicity indexes have been used to characterize the reactivity pattern of the CC double bond towards nucleophilic addition reactions. A wide family of molecules including ketones, esters, anhydrides, nitriles and nitrocompounds containing appropriate substitution on the CC double bond have been classified within an unique scale of reactivity. The predictive capability of the theoretical model is tested against a series of benzylidenemalononitriles and substituted α-nitrostilbenes.     

Vibrational analysis of small species in the liquid phase by a combined density functional theory and polarizable continuum method

The vibrational frequencies of small organic molecules and inorganic ions were predicted both in vacuum and in the liquid phase at the B3LYP/6-311G++(d,p) level of theory, within the harmonic approximation. The solvent effect was introduced as an electrostatic influence of the bulk by means of the integral equation formalism of the polarizable continuum model. The results show that the application of this continuum solvation model can reduce the overall error in the theoretical predictions. In order to improve the quality of the results, two different scaling procedures were applied. The frequencies obtained for the continuum, when compared with vacuum calculations, show a better linear correlation with experimental data, which increases the efficiency of the scaling procedures. A further reduction of the errors is obtained by partitioning the overall set of frequencies according to the nature of the normal modes and the charge of the solutes. The quality of the theoretical predictions, especially for the stretching modes of vibration of nonionic chemical species, is here noteworthy.     

DFT Study of Alkali Metal Atom Adsorption on Defect-Free MgO(001) Surface

The adsorption of isolated alkali metal atoms (Li, Na, K, Rb, and Cs) on defect-free sur-face of MgO(001) has been systemically investigated with density functional theory using a pseudopotential plane-wave approach. The adsorption energy calculated is about -0.72 eV for the lithium on top of the surface O site and about one third of this value for the other alkali metals. The relatively strong interaction of Li with the surface O can be explained by a more covalent bonding involved, evidenced by results of both the projected density of states and the charge density difference. The bonding mechanism is discussed in detail for all alkali metals.     

低标度XYG3双杂化密度泛函的开发与测评

以XYG3为代表的双杂化泛函代表了目前密度泛函方法的最高精度,可以胜任复杂化学环境中不同类型的相互作用的精准描述.然而由于XYG3采用二级微扰相关能(PT2)的形式引入未占轨道信息,其计算标度随体系尺寸呈5次方增长,同时内存需求呈4次方增长,这极大削弱了XYG3方法在复杂大体系研究中的可用性.本文介绍了一种结合局域密度拟合方法(Local Resolution-of-Identity, RI-LVL)与拉普拉斯变换(Laplace Transformation, LT)的低标度算法.采用MPI与OpenMP的混合并行设计,实现了可兼顾计算精度、计算效率和数值稳定性的低标度XYG3(LT-XYG3)计算功能.在一系列水簇及ISOL22测试集的测评中, LT-XYG3有效地降低了双杂化泛函的计算成本与内存消耗,展示了良好的可靠性与应用前景.     

Generalized Christoffel–Darboux formulae and the frontier Kohn–Sham molecular orbitals

The Christoffel–Darboux formula for classical orthogonal polynomials is generalized to arbitrary sets of orthogonal functions in three dimensions, yielding an explicit link between frontier Kohn–Sham molecular orbitals and the Kohn–Sham density matrix. Methods using this result could significantly accelerate Kohn–Sham density functional theory calculations, as only a subset of the Kohn–Sham equations would need to be addressed. The result can also be seen as an explicit justification for the utility of frontier molecular orbital theory.     

DFT Study on Homolytic Dissociation Enthalpies of C-I Bonds

The C-I bond dissociation enthalpies （BDE） of various organic iodides were calculated using high-level theoretical methods including MP2 and CCSD（T） with extrapolated basis set as well as a number of density functional theory methods. After systematic evaluation of the theoretical results against available experimental C-I BDEs, it was found that the MPW LYPIM method gave the lowest root mean square error. We, therefore, used this method to examine the substituent effects on different categories of C（sp3）-I and C（sp2）-I bonds. Fur thermore, the remote substituent effects on the C-I BDEs of substituted iodobenzenes and substituted （iodomethyl）benzenes were also investigated at the same level. The C-I BDEs of typical heteroaromatic iodides including five-membered and six-membered heterocyclic iodides were also examined.     

DFT characterization of the mechanism for Staudinger/aza-Wittig tandem organocatalysis

A computational simulation with DFT calculations and microkinetic modeling is carried out on a complete catalytic cycle, involving Staudinger ligation, aza-Wittig condensation and phosphine oxide recycle, for non-truncated substrates and catalyst. The Staudinger reaction produces phosphazenes (R₃P?=?NR), also known as iminophosphoranes, from phosphines and organic azides. Electrophilic carbonyl groups react with phospazenes to produce imines and phosphine oxides. Recently the Staudinger reaction and the aza-Wittig condensation have been combined to spawn intramolecular tandems producing cyclic molecules of great pharmaceutical interest (e.g. benzoxazoles). The release of diatomic nitrogen combined with the formation of phosphine oxide represents the driving force of the reaction. The implementation of in situ recycling of the exhausted phosphine oxide into the Staudinger/aza-Wittig tandem improves the scopes and the applicability of the reaction, transforming it into a powerful and versatile synthetic tool.     

甲烷与含氮有机杂环化合物吸附作用的DFT研究

根据煤中氮的存在形式,提炼出11种具有不同杂化方式及含氮量的有机环状化合物.应用密度泛函理论(DFT)模拟方法对甲烷在这11种含氮化合物上的吸附模型进行结构优化,并结合吸附能、电荷分布及静电势能面的分析,从微观角度考察了甲烷与这些化合物之间的相互作用.计算结果表明:甲烷与含氮化合物的相互作用能在3.81-6.82kJ·mol-1范围内,且通过氢键和静电力相互作用;当化合物中氮的杂化方式为sp2时,其与甲烷的作用能大于sp3杂化方式的;当化合物中的氮含量增加时,可以提供更多的甲烷吸附位点.     

The C−NO2 bond dissociation energies of some nitroaromatic compounds: DFT study

The C−NO₂ bond dissociation energies in nitrobenzene; 3-amino-nitrobenze; 4-amino-nitrobenze; 1,3-dinitrobenzene; 1,4-dinitrobenzene; 2-methyl-nitrobenzene; 4-methyl-nitrobenzene; and 1,3,5-trinitrobenzene nitroaromatic molecules, are computed using B3LYP, B3PW91, B3P86 three-parameter hybrid Density Functional Theory (DFT) methods in conjunction with 6-31G^** basis set. By comparing the computed energies and experimental ones, it is found that B3P86/6-31G^** is not capable of predicting the satisfactory bond dissociation energy (BDE). The BDEs computed with both B3LYP/6-31G^** and B3PW91/6-31G^** for the nitroaromatic molecules are closer to the experimental ones than those obtained with B3P86/6-31G^**. But, when compared with the experimental one, the BDE from the B3LYP/6-31G^** has the maximum deviation, which is completely outside our desired target accuracy for chemical predictions (less than 2.00 kcal mol⁻¹). Therefore, we suggest B3PW91/6-31G^** method as a reliable method of computing the BDE for removal of the nitrogen dioxide group in the nitroaromatic compounds. In addition, the C−NO₂ BDEs for 2,4,6-trinitrotoluene (TNT), triaminotrinitrobenzene (TATB), diaminotrinitrobenzene (DATB), and picramide are studied with B3PW91/6-31G^** method.