系列金化合物[X-{Au(PMe_3)}_2]的结构、电子光谱和非线性光学性质

采用密度泛函理论(DFT)B3LYP和BhandHLYP方法,对6个不同金化合物[X-{Au(PMe3)}2]的几何结构、电子光谱以及极化率和三阶非线性光学(NLO)效应进行了计算分析.本文采用的计算方法和基组适合所研究的对象,并发现用PMe3代替PPh3对分子整体结构影响不大.另外,分子桥连部分的空间效应对极化率有明显影响,对三阶极化率影响却不大.关于6个分子的三阶NLO性质,由于X部分的电子性质及共轭程度不同,分子1a的三阶极化率γ值最小,分子2a的γ值最大.通过分析电子光谱和对应的分子轨道组成可知,Au在分子1a中做电子给体,而在分子2a～6a中做电子受体,表明Au在此类化合物中对NLO性质的贡献不同.     

N-亚苄基苯胺衍生物的非线性光学有限场/AM1计算

利用有限场/AM1方法对具有D-π-A结构的N-亚苄基苯胺衍生物的一阶极化率、非线性光学二阶和三阶极化率进行了计算,并考察和讨论了不同取代基对分子的非线性光学性质的影响。     

N-亚苄基苯胺衍生物的非线性光学有限场/AM1计算

利用有限场／ＡＭ１方法对具有Ｄ－π－Ａ结构的Ｎ－亚苄基苯胺衍生物的一阶极化率、非线性光学二阶和三阶极化率进行了计算,并考察和讨论了不同取代基对分子的非线性光学性质的影响。     

X…H2O（X=Li，Na，K）非线性光学性质的从头算理论

在MP2水平上采用6-311G基组计算了van der Waals复合物X…H2O(X=Li, Na, K)的非线性光学性质(μ, α, β), 讨论了基组效应和电子相关效应对计算结果的贡献. 在MP2/6-311++G(2df, 2pd)水平上计算得到的三个复合物分子X(X=Li, Na, K)•••H2O的非线性光学性质. 结果表明, 三种复合物分子均具有巨大的一阶超极化率, 其中最外层电子的弥散特性对一阶超极化率有很大的影响.     

取代苯体系的二阶非线性光学性质:动力学李代数方法

提出了一种动力学李代数方法来研究取代苯体系的非线性光学性质. 对于给定的PPP模型(Pariser-Parr-Pople)哈密顿量, 生成了一个动力学李代数. 依据这些代数元构造出演化算子作为群参数的函数, 通过求解一组非线性微分方程能够得到这些群参数. 再按照统计力学中的密度算子公式给出取代苯分子体系偶极矩的统计平均值. 于是导出二阶极化率的表达式. 与其他量子力学计算结果比较, 表明这种动力学李代数方法在预言有机共轭分子的非线性光学性质上同样有用.     

C50富勒烯及其二聚物C100、C101的光学性质

运用B3LYP方法在6-31G*基组水平上对C50富勒烯以及它的两个不同二聚物C100、C101的几何构型进行了全优化. 在优化所得构型的基础上, 采用TDB3LYP方法在3-21G*基组水平上对其激发态性质、电子吸收光谱进行了研究, 根据计算得到的态态间跃迁偶极矩和跃迁能等数据, 结合使用态求和公式进一步计算得到了它们不同光学过程中的三阶非线性极化率. 结果表明, 当C50富勒烯二聚以后, 其电子吸收光谱的最大波长吸收峰发生了明显的红移, 三阶非线性极化率有了较大的提高. 其中, [5,5]-[5,5]哑铃型二聚物C101有着比[2+2]闭环型二聚物C100更大的三阶非线性极化率.     

阴离子自由基TCNQ-及其铜盐CuTCNQ非线性光学性质的理论研究

采用量子化学方法计算分析了阴离子自由基TCNQ-及其铜盐CuTCNQ单体、二聚体的非线性光学(NLO)性质.结果表明,将Cu+引入TCNQ-分子,极化率值减小,Cu+的位置不同对分子NLO系数的贡献不同;二聚体分子的极化率与其结合能相关,同时分子的几何构型影响二阶NLO系数.由前线分子轨道组成得到,电子在前线分子轨道间...     

溶致变色法测定Schiff碱类化合物非线性二阶极化率张量分量β~xxx和估算二阶极化率向量分量β~x

提出了一个以二能级模型为基础,结合溶致变色法计算非线性光学分子二阶极化率向量分量β_x的新方法.使用溶致变色法测定了新的Schiff碱化合物的二阶极化率张量分量β_(xxx),并对其非线性光学特性进行了评估.     

一水甲酸锂晶体二阶非线性光学系数的理论计算

采用有限场方法和含时耦合微扰方法,在6-31+G*基组水平上并考虑了电子相关效应和色散效应的影响,计算了一水甲酸锂晶体的宏观二阶非线性光学极化率和二阶非线性光学系数。其中非线性光学系数d31和d32与Roberts报道的相符,而另一个系数d33则比文献值大。计算结果还表明,在非共振情况下,电子相关效应对非线性光学极化率的影响远远超过了色散效应,同时表明Roberts报道的非线性光学系数比Singh等人报道的更为合理。     

尿素分子团簇线性和一阶非线性光学极化率的理论计算

检验分子晶体非线性光学极化率的可加性,对于确认有向气体模型(Oriented-Gas Model) 具有重要的科学意义,对探索新型实用的分子非线性光学材料也有实际的意义。本文对尿素晶体中分子簇非线性光学极化率的可加性做了详尽的理论研究。首次做到了尿素分子最近邻排布团簇的非线性光学极化率的从头算,从而扩展了前人的工作。研究结果有效地确认了有向气体模型。同时,我们还提出了壳层结构模型 (Shell-Structure Model),用以解决有向气体模型中的缺陷,可实现更为有效可靠地计算晶体宏观非线性光学系数。     

一些密度泛函方法预测电子亲和势

选取了杂化泛函B3LYP, B3PW91, O3LYP, PBE0, 以及与之相对应的GGA泛函BLYP, BPW91, OLYP和PBE, 还选取了能更好地兼顾强相互作用和弱相互作用的X3LYP泛函和在预测NMR的化学位移有较好表现的OPBE泛函, 以及两种meta-GGA泛函VSXC和TPSS, 共12种泛函, 详细地考察了这些泛函在预测EA方面的准确性.     

[Li…X]e-[1](X=FH,OH2, NH3)的光电性质从头算

在MP2理论水平上采用6-311G基组系列计算了一价阴离子van der Waals复合物[Li…X]e-[1](X=FH, OH2, NH3)的偶极矩(μ)、平均极化率(α)以及平均一阶超极化率(β), 讨论了基组效应和电子相关效应对计算结果的影响, 比较了价电子对复合物一阶超极化率的贡献. 在MP4(SDQ)/6-311++G(2df, 2pd)水平上计算得到[Li…FH]e-[1]的μ=2.5633 a.u., α=1.0476×10³ a.u., β=1.0948×10⁵ a.u.；[Li…OH2]e-[1] 的μ=2.3204 a.u., α=1.2201×10³ a.u., β=2.1410×10⁵ a.u.；[Li…NH3]e-[1]的μ=2.4687 a.u., α=1.4817×10³ a.u., β=3.4040×10⁵ a.u.. 计算结果表明, 三种一价阴离子复合物分子均具有非常大的一阶超极化率, 而一个价电子对复合物的一阶超极化率的贡献超过1.0×10⁵ a.u..     

Time-dependent density functional theory for nonlinear properties of open-shell systems

This paper presents response theory based on a spin-restricted Kohn-Sham formalism for computation of time-dependent and time-independent nonlinear properties of molecules with a high spin ground state. The developed approach is capable to handle arbitrary perturbations and constitutes an efficient procedure for evaluation of electric, magnetic, and mixed properties. Apart from presenting the derivation of the proposed approach, we show results from illustrating calculations of static and dynamic hyperpolarizabilities of small Si(3n+1)H(6n+3) (n=0,1,2) clusters which mimic Si(111) surfaces with dangling bond defects. The results indicate that the first hyperpolarizability tensor components of Si(3n+1)H(6n+3) have an ordering compatible with the measurements of second harmonic generation in SiO2/Si(111) interfaces and, therefore, support the hypothesis that silicon surface defects with dangling bonds are responsible for this phenomenon. The results exhibit a strong dependence on the quality of basis set and exchange-correlation functional, showing that an appropriate set of diffuse functions is required for reliable predictions of the first hyperpolarizability of open-shell compounds.     

Computational study of static first hyperpolarizability of donor-acceptor substituted (E)-benzaldehyde phenylhydrazone

Various hybrid functionals (B3LYP, B97-2, PBE0, BMK, BH&HLYP, CAM-B3LYP, and LC-ωPBE) implemented in density functional theory were applied to give estimate of static first hyperpolarizabilty (β(0)) of (E)-benzaldehyde phenylhydrazone designated as (E)-BPH. Against those of MP2 computations as a function of the underlying density functional, good agreement was obtained with the BH&HLYP and CAM-B3LYP functionals. The LC-ωPBE functional and the B3LYP, PBE0, B97-2, and BMK functionals underestimated and overestimated β(0), respectively. The basis set effect on the calculated β(0) was also investigated. It turned out that the 6-311+G(2d,p) basis set provided excellent converged value of β(0). On the basis of the calculated results, we investigated the substituent effect on β(0) of donor-acceptor (D-A) substituted (E)-BPH systematically by using the BH&HLYP and CAM-B3LYP computations with the 6-311+G(2d,p) basis set. We proposed a Zwitterion structure to explain the calculated trend in the substituent effect and the enhanced hyperpolarizability of type II compounds (A-(E)-BPH-D) than type I compounds (D-(E)-BPH-A). Natural bonding orbital analysis carried out at BH&HLYP/6-311+G(2d,p)//B3LYP/6-31G(2df,p) level of theory substantiated the claim.     

The role of the basis set and the level of quantum mechanical theory in the prediction of the structure and reactivity of cisplatin

In this article, we conducted an extensive ab initio study on the importance of the level of theory and the basis set for theoretical predictions of the structure and reactivity of cisplatin [cis‐diamminedichloroplatinum(II) (cDDP)]. Initially, the role of the basis set for the Pt atom was assessed using 24 different basis sets, including three all‐electron basis sets (ABS). In addition, a modified all‐electron double zeta polarized basis set (mDZP) was proposed by adding a set of diffuse d functions onto the existing DZP basis set. The energy barrier and the rate constant for the first chloride/water exchange ligand process, namely, the aquation reaction, were taken as benchmarks for which reliable experimental data are available. At the B3LYP/mDZP/6‐31+G(d) level (the first basis set is for Pt and the last set is for all of the light atoms), the energy barrier was 22.8 kcal mol^?1, which is in agreement with the average experimental value, 22.9 ± 0.4 kcal mol^?1. For the other accessible ABS (DZP and ADZP), the corresponding values were 15.4 and 24.5 kcal mol^?1, respectively. The ADZP and mDZP are notably similar, raising the importance of diffuse d functions for the prediction of the kinetic properties of cDDP. In this article, we also analyze the ligand basis set and the level of theory effects by considering 36 basis sets at distinct levels of theory, namely, Hartree‐Fock, MP2, and several DFT functionals. From a survey of the data, we recommend the mPW1PW91/mDZP/6‐31+G(d) or B3PW91/mDZP/6‐31+G(d) levels to describe the structure and reactivity of cDDP and its small derivatives. Conversely, for large molecules containing a cisplatin motif (for example, the cDDP‐DNA complex), the lower levels B3LYP/LANL2DZ/6‐31+G(d) and B3LYP/SBKJC‐VDZ/6‐31+G(d) are suggested. At these levels of theory, the predicted energy barrier was 26.0 and 25.9 kcal mol^?1, respectively, which is only 13% higher than the actual value. © 2012 Wiley Periodicals, Inc.     

Irreversibility,decay, and asymptotic dynamics

Ab initio molecular-orbital theory has been used to study the 1,3-sigmatropic hydrogen rearrangements: propene → propene, formic acic → formic acid, and vinyl alcohol → acetaldehyde. Fully optimized structures of stable molecules and transition states have been determined using gradient procedures and the 4-31G basis set. Improved energies have been obtained using a variety of techniques with basis sets up to the size of double-ζ plus polarization (DZP ) and electron correlation up to the CEPA /DZP level. Although both polarization functions and electron correlation lead to a lowering of the calculated barriers, the values remain substantial for all three rearrangements.     

Molecular energies and properties from density functional theory: Exploring basis set dependence of Kohn—Sham equation using several density functionals

The performance of four commonly used density functionals (VWN, BLYP, BP91, and Becke's original three-parameter approximation to the adiabatic connection formula, referred to herein as the adiabatic connection method or ACM) was studied with a series of six Gaussian-type atomic basis sets [DZP, 6–31G**, DZVP, TZVP, TZ2P, and uncontracted aug-cc-pVTZ (UCC)]. The geometries and dipole moments of over 100 first-row and second-row molecules and reaction energies of over 300 chemical reactions involving such molecules were computed using each of the four density functionals in combination with each of the six basis sets. The results were compared to experimentally determined values. Based on overall mean absolute theory versus experiment errors, it was found that ACM is the best choice for predictions of both energies of reaction [overall mean absolute theory versus experiment error (MATvEE) of 4.7 kcal/mol with our most complete (UCC) basis set] and molecular geometries (overall MATvEE of 0.92 pm for bond distances and 0.88° for bond angles with the UCC basis set). For routine calculations with moderate basis sets (those of double-ζ type: DZP, 6–31G**, and DZVP) the DZVP basis set was, on average, the best choice. There were, however, examples of reactions where significantly larger basis sets were required to achieve reasonable accuracy (errors ≤ 5 kcal/mol). For dipole moments, ACM, BP91, and BLYP performed comparably (overall MATvEE of 0.071, 0.067, and 0.059 debye, respectively, with the UCC basis set). Basis sets that include additional polarization functions and diffuse functions were found to be important for accurate density functional theory predictions of dipole moments. © 1997 by John Wiley & Sons, Inc.     

Density functional theory studies on the dissociation energies of metallic salts: relationship between lattice and dissociation energies

The formation and physicochemical properties of polymer electrolytes strongly depend on the lattice energy of metal salts. An indirect but efficient way to estimate the lattice energy through the relationship between the heterolytic bond dissociation and lattice energies is proposed in this work. The heterolytic bond dissociation energies for alkali metal compounds were calculated theoretically using the Density Functional Theory (DFT) of B3LYP level with 6‐311+G(d,p) and 6‐311+G(2df,p) basis sets. For transition metal compounds, the same method was employed except for using the effective core potential (ECP) of LANL2DZ and SDD on transition metals for 6‐311+G(d,p) and 6‐311+G(2df,p) calculations, respectively. The dissociation energies calculated by 6‐311+G(2df,p) basis set combined with SDD basis set were better correlated with the experimental values with average error of ca. ±1.0% than those by 6‐311+G* combined with the LANL2DZ basis set. The relationship between dissociation and lattice energies was found to be fairly linear (r>0.98). Thus, this method can be used to estimate the lattice energy of an unknown ionic compound with reasonably high accuracy. We also found that the dissociation energies of transition metal salts were relatively larger than those of alkaline metal salts for comparable ionic radii. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 827–834, 2001     

Basis set dependence of the vibrational wavenumbers of the out-of-plane modes of conjugated π-electron ring systems

Basis set dependence of the vibrational wavenumbers of out-of-plane modes calculated at the MP2 level of ab initio molecular orbital theory is examined for benzene, p-benzoquinodimethane, p-benzoquinone, furan, and thiophene. Various polarization functions up to (3df,p) are used in combination with the 6-31G and 6-311G basis sets. It is shown that, especially in the case of normal modes with alternate out-of-plane motions of the carbon atoms (such as the ν₄ (b_2g) mode of benzene), the calculated wavenumbers depend strongly on the exponents (_d) of the d functions on the carbon atoms. It is therefore necessary to include d functions with an optimum exponent (_d0.4) on the carbon atoms to obtain reasonable out-of-plane vibrational force fields. In a few cases (such as the ν₁₆ (a₂) mode of furan), inclusion of a set of f functions on the carbon atoms has some effects on the calculated wavenumbers of out-of-plane modes. However, unless the basis set contains an optimum set of d functions, inclusion of a set of f functions does not improve the agreement between the observed and calculated vibrational wavenumbers. As a case with an exaggerated effect of basis set, it is shown that the wavenumber of the ν₄₁ (b_2g) mode of the planar optimized structure of p-benzoquinodimethane is calculated to be imaginary by using the 6-311G(d,p) or 6-311G(df,p) basis set at the MP2 level. For all the molecules treated in the present study, reasonable out-of-plane vibrational force fields are obtained by using the 6-31G(2df,p) and 6-311G(2df,p) basis sets.