Ab initio lattice dynamics of nonconducting crystals by systematic fragmentation

A systematic method for approximating the ab initio electronic energy of crystal lattices has been improved by the incorporation of long range electrostatic and dispersion interactions. The effect of these long range interactions on the optimization of the crystal structure is reported. The harmonic lattice dynamics have been evaluated to give phonon frequencies and neutron scattering intensities. Exemplary results are reported for diamond, silicon, and α-quartz using Hartree-Fock, M?ller-Plesset perturbation, and coupled-cluster levels of ab initio theory.     

Potential energy surface and MULTIMODE vibrational analysis of C2H3+

A full dimensional, ab initio-based semiglobal potential energy surface for C(2)H(3) (+) is reported. The ab initio electronic energies for this molecule are calculated using the spin-restricted, coupled cluster method restricted to single and double excitations with triples corrections [RCCSD(T)]. The RCCSD(T) method is used with the correlation-consistent polarized valence triple-zeta basis augmented with diffuse functions (aug-cc-pVTZ). The ab initio potential energy surface is represented by a many-body (cluster) expansion, each term of which uses functions that are fully invariant under permutations of like nuclei. The fitted potential energy surface is validated by comparing normal mode frequencies at the global minimum and secondary minimum with previous and new direct ab initio frequencies. The potential surface is used in vibrational analysis using the "single-reference" and "reaction-path" versions of the code MULTIMODE.     

取代对N,N-二甲基苯乙烯基蒽的含时密度泛函研究

运用B3LYP方法和HF单激发态相互作用（CIS）方法分别优化了2-对N,N-二甲基苯乙烯基蒽（2-pNC）和9-对N,N-二甲基苯乙烯基蒽（9-pNC）分子的基态及最低激发单重态几何结构.系统分析了前线分子轨道特征,并探索了电子跃迁机理,应用含时密度泛函理论计算了分子的电子光谱.     

Effects of ethynyl substituents on the electronic structure of cyclobutadiene

The effects of ethynyl substitution on the electronic structure of cyclobutadiene are investigated in this work. Ethynyl substituted cyclobutadienes may be involved in Bergman cyclization reactions and are possible intermediates in the formation of fullerenes and graphitic sheets. Prediction of the electronic structure of cyclobutadiene is challenging for single-reference ab initio methods because of Jahn-Teller distortions and the diradical character of the singlet state. The equation-of-motion spin-flip coupled-cluster with single and double excitations (EOM-SF-CCSD) method accurately describes diradical states and is used to determine vertical and adiabatic singlet-triplet energy splittings in the substituted cyclobutadienes. The adiabatic singlet-triplet gaps decrease upon substituent addition, but the singlet states remain lower in energy. However, the results are affected by spin-contamination of the reference state and deteriorate when an unrestricted HF reference is employed. Additional insights in the electronic structure of cyclobutadienes are obtained by analyzing natural charges and spin densities. The substituents pull the charge out of the cyclobutadiene ring; however, the natural charges and spin densities are found to be nearly independent of the geometry and spin state.     

Ab initio analysis of some fluoride and oxide structures doped with Pr and Yb

Doping effects are analyzed by means of the ab initio perturbed ion (aiPI) method via substitution of the fluorides and oxides of zirconium, hafnium, and thorium. Lattice relaxation is simulated through the calculation of vibrational breathing modes and substitutional effects on the compound are thereby analyzed. In addition, the band gaps (in the pure species) and impurity levels (where substitutional ions Pr⁺⁴ and Yb⁺³ are considered in the doped species) of the fluorides are estimated via the transition energy calculated at the aiPI‐optimized geometry of the pure and doped crystal clusters by means of the configuration interaction with single excitations (CIS) method that accounts for electronic correlation. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Character of the electronic ground state and of charge-transfer excited states in ionic solids: An ab initio cluster model approach

The electronic structure of the ground electronic state and of some special charge-transfer excited states in ionic solids is examined from the ab initio cluster model approach. Different ab initio wave functions, including a frozen orbital approach, the Hartree–Fock self-consistent field, and multireference configuration interaction wave functions, are considered and analyzed using different theoretical techniques. We explicitly consider some alkaline–earth oxides such as CaO, a more difficult case such as A1₂O₃, a transition-metal oxide such as NiO, and a system with a more complicated structure such as KNiF₃. Analysis of ab initio wave functions in terms of valence bond components shows that all these compounds are largely ionic, thus supporting the simple picture arising from the ionic model. However, the nature of the excited states is more complex. Alkaline–earth oxides lowest excited states are essentially described as charge-transfer excitations dominated by a single resonant valence bond structure and the calculated energy difference is comparable to the experimental optical gap. In the case of A1₂O₃, the electronic spectra presents excitonic features and the local charge-transfer excitation excited states provide a reasonable representation of these phenomena. Finally, several different valence bond structures are present in the lowest electronic states of KNiF₃. © 1994 John Wiley & Sons, Inc.     

Interactions between anionic and neutral bromine and rare gas atoms

High-quality, ab initio potential energy functions are obtained for the interaction of bromine atoms and anions with atoms of the six rare gases (Rg) from He to Rn. The potentials of the nonrelativistic (2)Sigma(+) and (2)Pi electronic states arising from the ground-state Br((2)P)-Rg interactions are computed over a wide range of internuclear separations using a spin-restricted version of the coupled cluster method with single and double excitations and noniterative correction to triple excitations [RCCSD(T)] with an extrapolation to the complete basis set limit, from basis sets of d-aug-cc-pVQZ and d-aug-cc-pV5Z quality. These are compared with potentials derived previously from experimental measurements and ab initio calculations. The same approach is used also to refine the potentials of the Br(-)-Rg anions obtained previously [Buchachenko et al., J. Chem. Phys. 125, 064305 (2006)]. Spin-orbit coupling in the neutral species is included both ab initio and via an atomic approximation; deviations between two approaches that are large enough to affect the results significantly are observed only in the Br-Xe and Br-Rn systems. The resulting relativistic potentials are used to compute anion zero electron kinetic energy photoelectron spectra, differential scattering cross sections, and the transport coefficients of trace amounts of both anionic and neutral bromine in the rare gases. Comparison with available experimental data for all systems considered proves a very high precision of the present potentials.     

PdO~(0,±1),PdH~(0,±1)及PdOH分子结构和电子性质的从头计算

利用多种从头计算方法对PdO0,±1,PdH0,±1以及PdOH的结构和电子性质进行了理论研究.优化了PdO分子两个能量相近的态(3Π和3Σ-),结果表明PdO的3Π态较3Σ-态更为稳定,在单双激发耦合簇理论(CCSD)水平下3Σ-态的总能量比3Π态的高出0.286eV.在结构优化的基础上计算了PdO以及PdH的绝热电离能(AIE)和电子亲和能(AEA),计算结果与实验值符合得很好.PdOH的基态为Cs对称性的角型结构Pd—OH(2A′态),另外还优化得到两个2A″态的亚稳结构,分别对应于Pd—OH和O—Pd—H,CCSD水平下两个亚稳态的总能量较基态分别高0.405和2.284eV,优化得到了连接这两个2A″态的过渡态,并计算了相应的反应能垒.     

Rovibrational dynamics of the strontium molecule in the A(1)Σ(u)+, c(3)Π(u), and a(3)Σ(u)+ manifold from state-of-the-art ab initio calculations

State-of-the-art ab initio techniques have been applied to compute the potential energy curves for the electronic states in the A(1)Σ(u)(+), c(3)Π(u), and a(3)Σ(u)(+) manifold of the strontium dimer, the spin-orbit and nonadiabatic coupling matrix elements between the states in the manifold, and the electric transition dipole moment from the ground X(1)Σ(g)(+) to the nonrelativistic and relativistic states in the A+c+a manifold. The potential energy curves and transition moments were obtained with the linear response (equation of motion) coupled cluster method limited to single, double, and linear triple excitations for the potentials and limited to single and double excitations for the transition moments. The spin-orbit and nonadiabatic coupling matrix elements were computed with the multireference configuration interaction method limited to single and double excitations. Our results for the nonrelativistic and relativistic (spin-orbit coupled) potentials deviate substantially from recent ab initio calculations. The potential energy curve for the spectroscopically active (1)0(u)(+) state is in quantitative agreement with the empirical potential fitted to high-resolution Fourier transform spectra [A. Stein, H. Kno?ckel, and E. Tiemann, Eur. Phys. J. D 64, 227 (2011)]. The computed ab initio points were fitted to physically sound analytical expressions, and used in converged coupled channel calculations of the rovibrational energy levels in the A+c+a manifold and line strengths for the A(1)Σ(u)(+)←X(1)Σ(g (+) transitions. Positions and lifetimes of quasi-bound Feshbach resonances lying above the (1)S(0) + (3)P(1) dissociation limit were also obtained. Our results reproduce (semi)quantitatively the experimental data observed thus far. Predictions for on-going and future experiments are also reported.     

Cluster models of ionic crystals: Band gaps

We performed electronic ab initio calculations of ion clusters simulating a piece of ionic crystal. Our main interest here is to estimate the (direct) energy band gap from the cluster energy levels. The model consists of a central cation surrounded by four shells of ions, with additional point charges embedding the cluster. We present results for the series of alkali halides with rock-salt structure. © 1996 John Wiley & Sons, Inc.     

一维酞菁聚合物的从头算晶体轨道研究

用密度泛函自洽场晶体轨道方法计算了几种一维酞菁聚合物及其取代衍生物的能带结构和电子性质,得到这些聚合物均是半导体.研究发现,聚合方式和共轭程度的差异不仅影响能隙的大小,而且影响能带的形状.取代基对聚合物能带形状的影响不大,但取代基的电子特性对能带的位置高低起到一定的作用.     

结构原子组成对Keggin型杂多酸铯盐上丙烷选择氧化性能的影响

系统地研究了具有不同结构原子组成的Keggin型杂多酸铯盐上丙烷选择氧化性能,重点考察了H3+n-xCsxPMo12-nVnO40(n=0-4,x=0-3)系列化合物以及H0.5Cs2.5PW12O40,H1.5Cs2.5PW11VO40和H1.5Cs2.5SiMo12O40样品.应用BET,UV-Vis,FTIR,XRD,TPR和SEM等手段研究了催化剂的物化性质.V取代Mo增加了H3+nPMo12-nVnO40(n=0-4)系列杂多酸的酸量和氧化性能.Cs+对质子的取代调变了样品的表面酸性以及氧化还原性,同时催化剂的热稳定性也随着Cs+取代数的递增而增强,当Cs+的取代数达到2以上时催化剂的热稳定性较好.Cs+取代对H3+n-xCsxPMo12-nVnO40(n=0-4,x=0-3)催化剂上丙烷选择氧化性能有重要影响.随着一级结构中V5+取代数目的增多,可能的Cs+取代数目相应减少;而对丙烯酸的选择性则存在一适宜的Cs+取代数目.Keggin型杂多阴离子的结构原子组成及平衡离子的种类和数量决定了其酸性、氧化还原性和热稳定性,这是影响该类催化剂上丙烷选择氧化性能的首要因素.研究结果表明,以P为中心原子,以Mo和V为配位原子的杂多酸铯盐[x(Cs+)=2.5]具有较好的催化性能.     

四面体型含碳簇合物MC4H12的电子结构计算

用四面体含碳原子簇模型MC₄H₁₂(M=C,B)和量子化学的Gaussian76从头计算方法,对金刚石晶体及其代位硼体系的电子结构进行了无调节参数计算。提出了以激发能概念确定原子簇模型的晶体带隙和以原子内层电子能量校正含杂原子簇杂质能级的观点。所得金刚石晶体带隙值以及硼杂质能级在带隙中的位置与实验值及其他有关计算值都符合得很好。     

An efficient and accurate molecular alignment and docking technique using ab initio quality scoring

An accurate and efficient molecular alignment technique is presented based on first principle electronic structure calculations. This new scheme maximizes quantum similarity matrices in the relative orientation of the molecules and uses Fourier transform techniques for two purposes. First, building up the numerical representation of true ab initio electronic densities and their Coulomb potentials is accelerated by the previously described Fourier transform Coulomb method. Second, the Fourier convolution technique is applied for accelerating optimizations in the translational coordinates. In order to avoid any interpolation error, the necessary analytical formulas are derived for the transformation of the ab initio wavefunctions in rotational coordinates. The results of our first implementation for a small test set are analyzed in detail and compared with published results of the literature. A new way of refinement of existing shape based alignments is also proposed by using Fourier convolutions of ab initio or other approximate electron densities. This new alignment technique is generally applicable for overlap, Coulomb, kinetic energy, etc., quantum similarity measures and can be extended to a genuine docking solution with ab initio scoring.     

Study of the electronic energy band structure of polyethylene using MINDO/3

The electronic structure of polyethylene has been calculated using crystal orbital theory and the latest version (MINDO/3) of the MINDO semi-empirical SCF MO method. The results are in better agreement with the observed photoelectron spectrum than those given by an analogous calculation using MINDO/2, which in turn has been shown superior to treatments based on EHT, CNDO, INDO, or ab initio SCF procedures.     

TD-DFT Study on Pyrazoline Derivatives

1 INTRODUCTION Pyrazoline compounds have been widely used in industry due to their good photoconductivity and high fluorescence quantum productivity[1]. In recent years, pyrazoline derivatives with high vitrification temperature Tg have been found to act as hole trans- fer materials in the electroluminescence device, which has captured the intensive interest of che- mists[2]. But so far, the electron transfer mechanism of these compounds is still unclear[3~7]. Some resear- ches have attr…     

Automatic generation of potential energy and property surfaces of polyatomic molecules in normal coordinates

A procedure for the automatic construction of Born-Oppenheimer (BO) potential energy and molecular property surfaces in rectilinear normal coordinates is presented and its suitability and accuracy when combined with vibrational structure calculations are assessed. The procedure relies on a hierarchical n-mode representation of the BO potential energy or molecular property surface, where the n-mode term of the sequence of potentials/molecular properties includes only the couplings between n or less vibrational degrees of freedom. Each n-mode cut of the energy/molecular property surface is first evaluated in a grid of points with ab initio electronic structure methods. The ab initio data are then spline interpolated and a subsequent polynomial fitting provides an analytical semiglobal representation for use in vibrational structure programs. The implementation of the procedure is outlined and the accuracy of the method is tested on water and difluoromethane. Strategies for improving the proposed algorithm are also discussed.     

Interpolating moving least-squares methods for fitting potential energy surfaces: an application to the H2CN unimolecular reaction

Classical trajectories have been used to compute rates for the unimolecular reaction H2CN-->H+HCN on a fitted ab initio potential energy surface (PES). The ab initio energies were obtained from CCSD(T)/aug-cc-pvtz electronic structure calculations. The ab initio energies were fitted by the interpolating moving least-squares (IMLS) method. This work continues the development of the IMLS method for producing ab initio PESs for use in molecular dynamics simulations of many-atom systems. A dual-level scheme was used in which the preliminary selection of data points was done using a low-level theory and the points used for fitting the final PES were obtained at the desired higher level of theory. Classical trajectories were used on various low-level IMLS fits to tune the fit to the unimolecular reaction under study. Procedures for efficiently picking data points, selecting basis functions, and defining cutoff limits to exclude distant points were investigated. The accuracy of the fitted PES was assessed by comparing interpolated values of quantities to the corresponding ab initio values. With as little as 330 ab initio points classical trajectory rate constants were converged to 5%-10% and the rms error over the six-dimensional region sampled by the trajectories was a few tenths of a kcal/mol.     

8-羟基喹啉铍及其衍生物电子光谱性质的含时密度泛函理论研究

运用密度泛函理论(DFT)B3LYP方法和abinitioHF单激发组态相互作用(CIS)法分别优化了有机金属配合物8-羟基喹啉铍(BeQ₂)及其3种衍生物分子的基态及最低激发单重态几何结构.系统分析了分子结构、前线分子轨道特征和能级分布规律以探索电子跃迁机理.应用含时密度泛函理论(TD-DFT)计算分子的电子光谱,揭示了BeQ₂及其衍生物的发光源于配体中π→π^*电子跃迁,指出通过配体修饰可以有效地影响配合物前线分子轨道分布,调整发光波段,并有效提高电荷转移量.