几个小分子静电势的密度泛函理论计算

运用多种密度泛函理论(DFT)方法和从头算(ab initio)方法,研究了具有代表性的一系列分子的静电势,选择QCISD计算出的分子静电势为标准静电势,通过比较多种DFT法和QCISD从头算,以及它们之间的均方根偏差和静电势差值曲线图,结果表明B3LYP-DFT法与QCISD-ab initio法计算结果较吻合,暗示了B3LYP-DFT法在得到分子静电势时是一个有用的工具,尤其对于从头计算难以进行研究的大分子体系.     

Crystal chemistry of molybdenum phosphides from density functional theory calculations

The structures, bonding and relative thermodynamic stabilities of the crystalline molybdenum phosphides Mo₃P, Mo₈P₅, Mo₄P₃, MoP, MoP₂, MoP₄, have been investigated by density functional theory calculations. The elastic stiffness coefficients for MoP and MoP₂ have been predicted.     

磺胺甲恶唑的振动光谱和密度泛函理论研究

磺胺甲恶唑(Sulfamethoxazole, SMX)是一种磺胺类广谱抗生素.基于密度泛函理论(DFT)的B3LYP方法,在B3LYP/6-311++G (d, p)基组水平上对SMX初始构型进行结构优化,得出该分子最稳定结构构型.同时,给出了SMX的键长、键角和二面角等空间结构参数值,确定该分子空间构型为非平面结构.相同基组水平下进行频率计算,通过频率校正因子0.9630修正并绘制了SMX分子的红外光谱和拉曼光谱图.发现光谱谱段主要分布在3600-2800 cm^-1和1700-300 cm^-1频率段,并且理论计算值和实验测定值有较好吻合,表明DFT计算方法是可取和可信的.借助GaussView 6.0软件和计算得到的势能分布结果对频率范围内每一个振动模式进行了指认归属.此外,采用Multiwfn和VMD程序计算并分析SMX分子的静电势和前线分子轨道,预测了该分子的反应活性位点.该结果为SMX分子的振动光谱检测提供基础数据,为其结构鉴定提供技术参考.     

Partition density functional theory and its extension to the spin-polarized case

Partition density functional theory (PDFT) [P. Elliott, K. Burke, M.H. Cohen, and A. Wasserman, Phys. Rev. A 82 (2), 024501 (2010)] is a formally exact method for obtaining molecular properties from Kohn–Sham calculations on isolated fragments. Here, we express the partition energy of PDFT as an implicit functional of the molecular spin-densities for a given choice of fragmentation, and use the principle of von Barth and Hedin to formulate the spin-decomposed version of PDFT. We introduce a partition energy functional of the spin-up and spin-down electronic densities and derive the associated polarized partition potentials, which are found to be global quantities that influence every fragment in the molecule. Along with the formal theory, we propose a simplified approach to computing the spin-partition potentials, and illustrate its utility and accuracy with two simple examples. Finally, we propose a viable approach to including external electric and magnetic fields in the framework of spin-PDFT.     

Classical density functional theory meets Monte-Carlo simulations

ABSTRACT

Typically the quality of an approximate density functional is evaluated by a direct comparison of its predictions in a given test case to exact data obtained by computer simulations. An important example for such an approach is the test of equilibrium structure of a simple fluid as measured by the pair distribution function g(r) or the cavity correlation function y(r). However, the combination of exact density profiles and the analytical structure of density functional theory allows one to determine and potentially improve the quality of a functional in a more sophisticatedway.     

有限长扶手椅形单壁碳纳米管分子静电势的理论研究

在混合密度泛函B3LYP理论下,用3-21G基函数对有限长扶手椅形单壁碳纳米管(4,4)、(5,5)和(6,6)的构型进行优化和分子静电势计算.结果表明:除近核区域为正常的正电势外,碳纳米管结构模型的管内和管外为负电势区域;在碳纳米管结构模型的管内,管心处均出现负电势的最小值,且负电势的绝对值随着碳纳米管的曲率降低而增大,管心轴线上静电势的变化随碳纳米管的曲率降低而减少,带电粒子流比较容易通过纳米管.     

On the relationship between the density functional formalism and the potential distribution theory for nonuniform fluids

It is shown that the variational principle for the grand potential of a nonuniform fluid as a functional of the singlet density yields the potential distribution theory for the equilibrium density. We derive the explicit form that the functional takes for a system of hard rods, and propose an approximate one for hard spheres. Attractive interactions are also considered in mean-field approximation. In all cases the pair direct correlation function of the nonuniform system is obtained and the density gradient expansion of the free energy is investigated.     

A 3D non-local density functional theory for any pore geometry

A general framework of classical non-local density functional theory (NLDFT) is presented, in order to consider the adsorption of spherical molecules in porous materials of any geometry. Fluid-fluid interactions and fluid-solid interactions can be repulsive or attractive. Some techniques that have been developed for the computation of weighted densities of hard-spheres are extended to attractive ones, in order to deal with an arbitrary pore geometry. This way, the computation method introduced in this work is validated by a comparison with analytical results for simple cases, and is directly applied to more complex systems. Density distributions depending on multi-dimensional effects are presented, and some radial distribution functions are recovered from NLDFT computations. Finally, the case of attractive continuous curved walls is detailed, which represents a large variety of real systems (e.g. micro and mesoporous silica, zeolites, carbonaceous nanoporous materials, etc.). With the new way of computation proposed, a general solution is presented, valid for any shape of continuous pore surface, by considering mathematical properties of discrete geometry due to the discretisation of the computational space with FFT computations.     

Toward ab initio density functional theory for nuclei

We survey approaches to non-relativistic density functional theory (DFT) for nuclei using progress toward ab initio DFT for Coulomb systems as a guide. Ab initio DFT starts with a microscopic Hamiltonian and is naturally formulated using orbital-based functionals, which generalize the conventional ‘local density plus gradients’ form. The orbitals satisfy single-particle equations with multiplicative (local) potentials. The DFT functionals can be developed starting from internucleon forces using wavefunction-based methods or by Legendre transform via effective actions. We describe known and unresolved issues for applying these formulations to the nuclear many-body problem and discuss how ab initio approaches can help improve empirical energy density functionals.     

基于密度泛函理论的叶绿素A性质的研究

在液相环境中,采用密度泛函理论(DFT)、含时密度泛函理论(TD-DFT)、Multiwfn波函数分析软件,在pbepbe/6-311g(d)基组水平上,计算并分析了叶绿素A的结构、紫外光谱和电子—空穴分布,结果表明:pbepbe/6-311g(d)方法是计算叶绿素A紫外吸收光谱更精确的方法;叶绿素A分子的吡咯环与取代基相互作用的过程中,吡咯环Ⅳ受侧链"尾巴"的影响最大;理论计算的紫外光谱与实验数据吻合较好,其中635.71 nm和446.87 nm处的两个吸收峰可认为是叶绿素A的特征吸收峰;侧链或取代基团在叶绿素A激发过程中是给电子体,卟啉"头"既是电子供体,也是电子受体.     

Vibrational assignment of isoxazole aided by rovibrational data and density functional theory

IR and Raman spectra of isoxazole and D₃-isoxazole are reassigned with the aid of density functional theory calculations at the B3LYP/6-311+G(d,p) level. Harmonic wavenumber values of all non-CH stretch modes give an r.m.s. deviation from experimental values of 3 cm⁻¹ when they are uniformly scaled by 0.9818. Anharmonic values obtained using a second-order perturbative approach show an r.m.s. deviation of 7.6 cm⁻¹. Results of equal quality are obtained for D₃-isoxazole. Predicted Raman and IR band strengths match experimental spectra very closely. Detailed re-examination of rovibrational constants and inertial defect data from high resolution IR spectra, in comparison with values from the anharmonic frequency calculation, provides confirmation of the vibrational assignments. ‘Dark state’ fundamentals ν₁₂ (900.2 cm⁻¹) and ν₁₅ (866 cm⁻¹) are detected through the perturbations they cause in nearby bands.     

基于密度泛函理论的CO2吸附微观机理比较研究

为探究吸附法捕获CO₂过程中的微观机理和吸附剂材料间的作用关系,基于密度泛函理论方法,综合比较了典型吸附剂包括煤基官能团、Fe、限域离子液体、Na₂CO₃、SrTiO₃与CO₂的吸附过程和差异性.根据不同计算策略,着重分析比较了吸附能、结构优化参数、吸附构型以及原子分布等参数.结果表明,化学吸附中CO₂分子与吸附面呈平行关系时通常吸附能最大;在一种材料的同类型官能团中,吸附能大小与氧原子的数量呈正相关关系;吸附过程中C-O键的伸长活化会生成一种重要的中间产物CO₂^-.提出在探寻CO₂吸附材料时可以在含氧原子较多的官能团、活性金属表面等方面进一步探究.最后对基于密度泛函理论的CO₂的吸附机理的进一步研究方向进行了展望.     

基于密度泛函理论的三聚氰胺结构及振动光谱研究

利用密度泛函(DFT)算法,对三聚氰胺分子的空间结构进行了优化;然后分别用MP2/6-31G和DFT/DGTIVP两种算法计算了它的拉曼光谱与红外光谱,给出光谱强度图;对比了两种算法的拉曼光谱图和其实验光谱图,结果显示了很好的一致性;给出了三聚氰胺分子中的各原子间键长,键角等空间结构参数;并对三聚氰胺分子在550～3 800cm-1区间的振动谱做了指认。上述工作将有助于食品中三聚氰胺含量测量技术的研究。     

密度泛函理论下氯普鲁卡因分子的能级结构和光谱计算

本文基于密度泛函理论,采用B3LYP方法,在6-31G(d,p)基组上对麻醉剂氯普鲁卡因的分子结构进行几何优化,在此基础上以乙醇为溶剂计算分子的前20个激发态,所有计算在Gaussian 09W-D01中进行。利用Multiwfn3.7软件绘制红外光谱图,并对其分子振动进行分析;利用Origin 2018 64Bit软件和Multiwfn3.7软件相结合绘制紫外光谱图,并计算空穴-电子来分析分子的激发态性质;通过计算前线轨道来预测氯普鲁卡因分子的活性位点。结果表明,在所计算得到的激发态中,由基态到第2、3、6、10激发态为局域激发,由基态到第19激发态为电荷转移激发。氯普鲁卡因乙氨基上的N22为亲电反应位点,苯环上的碳原子和脂基上的氧原子为亲核反应位点。本研究对更好的了解氯普鲁卡因分子的反应机理和在医学上的麻醉活性提供理论参考。     

A perturbation density functional theory for hydrogen bonding cyclic molecules

Using the framework of Wertheim's thermodynamic perturbation theory, a new polyatomic density functional theory is developed to account for the intermolecular association of cyclic molecules in interfacial systems. To test the theory, Monte Carlo simulations in the canonical ensemble were performed for the specific case of an associating triatomic ring with one association site next to a hard wall. The theory and simulation results were found to be in good agreement.     

Spontaneous spin polarization in rubrene studied by density functional theory calculations

We theoretically studied the spontaneous spin polarization properties of organic molecule rubrene by using density functional theory calculations. Our investigations show that normally nonmagnetic molecule rubrene could be spin polarized by spinless-hole injection. Magnetic moment of the molecule increases linearly with the extra hole charge amount only when the injected hole charges reach a certain value. The spin density resides predominantly on the carbon atoms in the tetracene backbone of rubrene molecule and also the bond lengths change differently due to the injected charge. Spontaneous spin polarization can be explained as the preferably filling of the spin-splitted carbon p_z orbitals near the Fermi energy for the injected charge.     

Computational studies of CO and CO: Density functional theory and time-dependent density functional theory

Potential energy curves, equilibrium interatomic distances, term energies and harmonic vibration frequencies for the 16 lowest states of neutral carbon monoxide and the six lowest states of singly ionized carbon monoxide are calculated by density functional theory (DFT) and linear-response time-dependent density functional (LR-TDDFT) theory. The results are compared with experimental data. The two theories, DFT and LR-TDDFT, are described briefly.     

基于密度泛函理论的单层氢化石墨烯特性分析

本文采用密度泛函理论的第一性原理方法,研究了不同尺寸H-graphene的稳定性、HOMO-LU-MO能隙以及电子激发态.研究结果表明,对于C_(16)H_(10)、C_(30)H_(14)、C_(48)H_(18)、C_(70)H_(22)、C_(96)H_(26)、C_(126)H_(30)计算的比结合能,C_(126)H_(30)相比C_(16)H_(10)的比结合能增长23.9%,且比结合能随着H-graphene尺寸扩大而增加,意味着稳定性不断提高.通过对HOMO-LUMO能隙分析发现,在较小尺寸的H-graphene中,由于量子效应起主要作用,因此出现了较大的HOMO-LUMO能隙,且随着H-graphene团簇尺寸的增加,能隙逐渐缩小可以看出,对于无限大的H-graphene团簇中,HOMO-LUMO能隙无限趋近于零(相当于零带隙),其电子性质与纯石墨烯相似.通过分析C_(16)H_(10)、C_(30)H_(14)、C_(48)H_(18)、C_(70)H_(22)激发态以及了吸收光谱,发现随着尺寸的扩大,吸收光谱发生红移,为石墨烯在电子器件领域的应用提供理论基础.