Surface rumpling of cubic CaTiO3 from density functional theory

In this paper, the structure of cubic CaTiO3 （001） surfaces with CaO and TiO2 terminations has been studied from density functional calculations. It has been found that the Ca atom has the largest relaxation for both kinds of terminations, and the rumpling of the CaO-terminated surface is much larger than that of TiO2-terminated surface. Also we have found that the metal atom relaxes much more prominently than the O atom does in each layer. The CaO-terminated surface is slightly more energetically favourahle than the TiO2-terminated surface from the analysis of the calculated surface energy.     

稳态Cu-Zr二十面体团簇电子结构的密度泛函研究

采用第一原理对以Cu为心的低能稳态Cu_nZ_(r13-n)(n=6,7,8,9)二十面体团簇的电子结构进行计算,结果表明:同一化学组分下,以Cu为心的Cu-Zr二十面体团簇中出现的同类原子聚集现象可以增强团簇的稳定性,降低费米能级(EF)上的电子数N(EF),这为低能稳态团簇拥有较小的N(EF)提供了深层次的理论解释.进一步的差分电子密度与Mulliken布居分析得知,Cu-Zr二十面体中共价键与离子键共存,成键态与反键态共存,且团簇在形成时壳层Zr与中心Cu原子是电子的提供者,壳层Cu是电子的获得者.该电荷转移方向是金属玻璃中以Cu为心的Cu-Zr二十面体团簇普遍遵循的规律,不随团簇的化学序参数及化学组分的变化而变化.计算的红外振动谱为实验上准确表征不同二十面体原子团提供了一种新的思路.     

Microscopic study of neutron-induced fission process of 239Pu via zero- and finite-temperature density functional theory

To study the neutron-induced fission of \begin{document}$ ^{239} $\end{document}Pu, potential energy surface (PES) calculations were performed using zero and finite-temperature density functional theory (FT-DFT) with the Skyrme force. The energy of the incident neutron was simulated by the temperature of the FT-DFT. The variations of the least-energy fission path, fission barrier, total kinetic energy, scission line, and mass distribution of fission fragments with the incident neutron energy were analyzed. It was learned that an increase in the temperature lowers the barrier height, the isomeric-state energy, and the ridge between symmetric and asymmetric fission valleys. Additionally, the gaps of the single particle levels become smaller with an increase in the temperature. As the temperature increases, the pre-fission region shrinks, and the scission occurs at smaller deformation around the symmetric fission channel. At low temperatures, the pairing correlations in the collective space are similar to those in zero-temperature DFT, and when the temperature is \begin{document}$ T> $\end{document} 0.3 MeV, the pairing gaps decrease rapidly. Two different methods were used to calculate the fission yields of the neutron-induced fission \begin{document}$ ^{239} $\end{document}Pu (n, f) with different incident neutron energies, in the framework of time-dependent generator coordinate method (TDGCM). One way to calculate the fission yield of \begin{document}$ ^{239} $\end{document}Pu (n, f) is to solve the collective equation of the TDGCM by using the PES from the FT-DFT with the corresponding temperature. The other involves using the PES from the zero-temperature DFT and adjusting the initial collective energy of the wave packet in the TDGCM according to the incident neutron energy. For the cases of the lower incident neutron energies, these two methods gave similar results and reproduced the experimental peak and width of fission fragment distribution. However, for the highest incident neutron energy considered in this study, the results from the TDGCM using the PES from zero-temperature DFT deviated explicitly from the experimental data, whereas those obtained by using the PES from FT-DFT remained close to the experimental data. This indicated that, with the increase in the incident neutron energy, the shell structure of the compound nuclei changed explicitly; thus, it may not be effective to use the PES from zero-temperature to perform the fission dynamic calculation.     

Li掺杂8-羟基喹啉铝的密度泛函理论研究

采用密度泛函理论研究了Li原子掺杂8-羟基喹啉铝(Alq₃)分子的几何构型、 前线分子轨道及电子转移特性. 研究结果表明, Li原子掺杂Alq₃后, Li原子与Alq₃的O, N原子键合, 形成电子转移复合物. Li原子将部分电子转移到Alq₃的吡啶环上, 在Alq₃的带隙内形成施主能级, 这种n型掺杂结构有效地提高了电子的传输效率; 但过多的Li原子的掺杂会使Alq₃分解, 从而减弱其电子传输能力. 为使Alq₃的电子传输能力达到最高, Li原子的掺杂应保持在2:1左右的比例.     

氧化镁纳米管团簇电子结构的密度泛函研究

用密度泛函理论（DFT）的杂化密度泛函B3LYP方法在6-31G（d）基组水平上对MgO纳米管团簇的二元环双管、三元环、三元环双管三种构型共21个团簇进行优化,对各构型的平均结合能、能隙、平均原子电荷以及总电荷密度进行了理论研究. 结果表明,平均结合能和配位数呈线性关系;随着纳米管的生长,团簇的稳定性增加,其中以三元环纳米管最为稳定;生长过程中发生原子间的电荷转移现象,预测出至无限长时的平均原子电荷分别为1298,1270,1306;混合离子共价键始终存在于MgO纳米管团簇之中. 关键词： 氧化镁 纳米管团簇 密度泛函理论 电子结构     

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.     

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.     

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.     

Energy density functional analysis of the fission properties of 240Pu: The effect of pairing correlations

We have calculated the potential energy surfaces for ²⁴⁰Pu up to the scission point using the density functional theory with different pairing strengths to investigate the effect of pairing correlations on its fission properties. An enhancement in the pairing correlations lowers the barrier heights, isomeric state, and ridge between the symmetric and asymmetric fission valleys significantly. Moreover, it weakens the microscopic shell structure around the Fermi surface, shrinks the scission frontiers, especially for the symmetric and very asymmetric fission regions, and lifts the total kinetic energies (TKEs) for the symmetric fission region. It is also emphasized that the microscopic calculation qualitatively reproduces the trend of the distribution of the measured TKEs, especially for the positions of the peaks at begin{document}$A_{rm{frag}}simeq132$end{document} and begin{document}$A_{rm{frag}}simeq108$end{document}.     

Kinetic component of the correlation energy density functional: a quantitative description from information theory

In the framework of information theory, a new method to determine T _c , the kinetic energy component of the correlation energy density functional for atoms, is presented. This approach is based on Shannon entropy and information energy that are obtained by fractional occupation probabilities of natural atomic orbitals. It is indicated that the calculated Shannon entropy using discrete probabilities is an increasing function while information energy is a decreasing function of the number of electrons. An expression is proposed with explicit dependence on the Shannon entropy or information energy and atomic number for the purpose. Applications of formulas for estimation of T _c values for neutral atoms up to Xe and their first positive and negative ions are then examined and validity of the proposed approach is numerically verified.     

基于密度泛函理论的C_(24)H_(38)O_4分子外场效应研究

各种环境毒物危害着人类的健康,塑化剂更是破坏了食品安全.为研究外电场对塑化剂主要成分之一C_(24)H_(38)O_4(邻苯二甲酸二辛酯,dioctyl phthalate,DOP)的影响,采用密度泛函理论B3LYP方法在6-311G(d,p)基组水平上优化了不同静电场0—0.0125 a.u.(0—6.4278×10~9 V/m)作用下DOP分子的基态几何结构,在此基础上利用同样的方法计算了DOP分子的电偶极矩和分子总能量,最后利用含时密度泛函理论在同一基组下研究了不同外电场对DOP分子紫外-可见(UV-vis)吸收光谱产生的影响,并与实验测得的光谱图进行了比较.结果表明:分子几何构型与电场大小呈现强烈的依赖关系,分子偶极矩随着外电场的增强先减小后增加,而分子总能量随着外电场的增强先增加而后急剧减小;DOP分子激发态的振子强度随着外电场的增强而减小,UV-vis吸收峰显著红移.     

Electronic structures and magnetoelectric properties of tetragonal BaFeO3: an ab initio density functional theory study

First-principles calculations have been performed to investigate the ground state electronic properties of BaFeO₃ (BFO). Local spin density approximation (LSDA) plus U (LSDA+U) treatment modified the metallic behaviour to insulated one with a band gap of 4.12eV. The spontaneous polarization was found to be 89.3\muC/cm² with Berry phase scheme in terms of the modern theory of polarization. Fe-3d e_g were split into two singlet states (d_z² and d _x²-y²}), and Fe-3d t_2g were split into one doublet states(d_xz and d_yz) and one singlet states(d_xy) after Fe and O displaced along the c axis. Meanwhile the occupation numbers of d_z², d_yz, d_xz and O_T p_z (on the top of Fe) were increased at the expense of those in xy plane. Our results showed that it was the sensitivity of hybridization to ferroelectric distortions, not just the total change of hybridization, that produced the possibility of ferroelectricity. Moreover, the increasing occupation numbers of O_T p_z and Fe d_z² favoured the 180^o coupling between Fe-3d e_g and Fe-3d t_2g, leading to ferromagnetic ordering, which has been confirmed by the increase of magnetic moment by 0.13μ_B per formula unit in the polarized direction. Hence, the magnetization can be altered by the reversal of external electric field.     

碳分子线C5在激光场中的含时密度泛函理论研究

运用含时密度泛函理论和分子动力学相结合的方法, 研究了C₅分子线在强激光场中的电离激发.研究发现, 当考虑激光强度对C₅分子线激发的影响时, 激光强度越强, 分子吸收的能量越多, 电离也越早, 最终电离的电子也越多, 而且沿激光极化方向的偶极矩的变化及峰值也越大. 关于激光极化方向对C₅分子线激发的影响的研究表明, 当激光极化方向沿着C₅分子线轴向时, 分子的电离大大增强, x方向的激光脉冲仅能激发起x方向的偶极振荡, 而y方向的激光脉冲仅能激发起y方向的偶极振荡, 而且x方向的激光脉冲激发的偶极振荡强. 研究还表明, 当激光极化方向沿着C₅分子线轴向时, 尽管由于电离增强而导致C₅分子线C–C键振动的同步性变差, 但在两种激光极化方向情况下, C₅分子线的振动模式与中性C₅分子线的振动模式相同. 关键词： 含时密度泛函理论 分子动力学 分子电离 碳分子线     

用密度泛函理论研究氮化硅新相的电子结构、光学性质和相变

运用第一性原理赝势方法,对氮化硅新相(六方Pˉ6和Pˉ6′相)的电子结构、光学性质和相变过程进行分析,研究能带结构、介电函数谱、反射谱和能量损失函数的变化机理.研究发现,β→Pˉ6→δ相变是可行的,在室温下β→Pˉ6和Pˉ6→δ相变的临界压强分别为42.9和47.7 GPa;相界的斜率为正值表明Pˉ6→δ相变过程伴随着晶胞体积的塌缩;Pˉ6和Pˉ6′相分别属于直接带隙和间接带隙半导体,能隙宽度分别为4.98和4.01 eV;得到了两相的零频介电常数;反射谱表明,两相的强反射峰均位于真空紫外线区域,因此可以用作紫外光屏蔽或紫外探测材料;在可见光区域,两相表现为近似透明.     

Ab-initio density functional theory study of a WO3 NH3-sensing mechanism

WO 3 bulk and various surfaces are studied by an ab-initio density functional theory technique.The band structures and electronic density states of WO 3 bulk are investigated.The surface energies of different WO 3 surfaces are compared and then the (002) surface with minimum energy is computed for its NH 3 sensing mechanism which explains the results in the experiments.Three adsorption sites are considered.According to the comparisons of the energy and the charge change between before and after adsorption in the optimal adsorption site O 1c,the NH 3 sensing mechanism is obtained.     

质子与羟基碰撞的含时密度泛函理论研究

采用将含时密度泛函理论和分子动力学非绝热耦合的方法,研究了不同入射速度下质子与羟基碰撞的反应动力学.计算了碰撞前后质子动能和羟基动能的变化及羟基电子和质子的运动.计算结果表明,质子沿垂直羟基分子轴方向入射时,质子与羟基碰撞后,质子被反弹且动能损失并俘获了羟基中氧的一部分电子,而丢失部分电子的羟基则获得动能以伸缩振动的形式向计算边界平动.随着入射质子的初动能增加,质子从羟基中俘获的电子增多,碰撞后羟基的键长变长,羟基振动变强而伸缩振动频率降低.此外,还发现质子的入射方向对碰撞过程的激发动力学有很大的影响.质子从不同的方向入射时,质子的入射初动能越大,其损失的动能越多且损失的动能与入射初动能呈线性关系,而入射方向对质子动能损失的影响很小.在质子入射初动能较低(小于25 eV)的情况下,羟基获得的动能与质子入射初动能呈线性关系且与入射方向无关;在质子入射初动能较高(大于25 eV)时,当质子沿羟基分子轴方向入射时,羟基动能的增量远大于质子沿垂直于羟基分子轴方向入射时羟基动能的增量.     

Structure and electronic structure of S-doped graphitic C3N4 investigated by density functional theory

The structures of the heptazine-based graphitic C3N4 and the S-doped graphitic C3N4 are investigated by using the density functional theory with a semi-empirical dispersion correction for the weak long-range interaction between layers.The corrugated structure is found to be energetically favorable for both the pure and the S-doped graphitic C3N4.The S doptant is prone to substitute the N atom bonded with only two nearest C atoms.The band structure calculation reveals that this kind of S doping causes a favorable red shift of the light absorption threshold and can improve the electroconductibility and the photocatalytic activity of the graphitic C3N4.     

Microscopic study of higher-order deformation effects on the ground states of superheavy nuclei around 270Hs

We investigate the effects of higher-order deformations \begin{document}$\beta_\lambda$\end{document} (\begin{document}$\lambda=4,6,8,$\end{document} and 10) on the ground state properties of superheavy nuclei (SHN) near the doubly magic deformed nucleus \begin{document}$^{270}{\rm{Hs}}$\end{document} using the multidimensionally-constrained relativistic mean-field (MDC-RMF) model with five effective interactions: PC-PK1, PK1, NL3^*, DD-ME2, and PKDD. The doubly magic properties of \begin{document}$^{270}{\rm{Hs}}$\end{document} include large energy gaps at \begin{document}$N=162$\end{document} and \begin{document}$Z=108$\end{document} in the single-particle spectra. By investigating the binding energies and single-particle levels of \begin{document}$^{270}{\rm{Hs}}$\end{document} in the multidimensional deformation space, we find that, among these higher-order deformations, the deformation \begin{document}$\beta_6$\end{document} has the greatest impact on the binding energy and influences the shell gaps considerably. Similar conclusions hold for other SHN near \begin{document}$^{270}{\rm{Hs}}$\end{document}. Our calculations demonstrate that the deformation \begin{document}$\beta_6$\end{document} must be considered when studying SHN using MDC-RMF.     

Ultraviolet laser ionization studies of 1-fluoronaphthalene clusters and density functional theory calculations

This paper studies supersonic jet-cooled 1-fluoronaphthalene(1FN) clusters by ultraviolet(UV) laser ionization at 281 nm in a time-of-flight mass spectrometer.The(1FN) + n(n=1-3) series cluster ions are observed where the signal intensity decreases with increasing cluster size.The effects of sample inlet pressures and ionization laser fluxes to mass spectral distribution are measured.Using density functional theory calculations,it obtains a planar geometric structure of 1FN dimer which is combined through two hydrogen bonds.The mass spectra indicate that the intensity of 1FN trimer is much weaker than that of 1FN dimer and this feature is attributed to the fact that the dimer may form the first "shell" in geometric structure while the larger clusters are generated based on this fundamental unit.