Consistent gaussian basis sets of double‐ and triple‐zeta valence with polarization quality of the fifth period for solid‐state calculations

Consistent basis sets of double‐ and triple‐zeta valence with polarization quality for the fifth period have been derived for periodic quantum‐chemical solid‐state calculations with the crystalline‐orbital program CRYSTAL. They are an extension of the pob‐TZVP basis sets, and are based on the full‐relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and on the def2‐SVP and def2‐TZVP valence basis of the Ahlrichs group. We optimized orbital exponents and contraction coefficients to supply robust and stable self‐consistent field (SCF) convergence for a wide range of different compounds. The computed crystal structures are compared to those obtained with standard basis sets available from the CRYSTAL basis set database. For the applied hybrid density functional PW1PW, the average deviations of calculated lattice constants from experimental references are smaller with pob‐DZVP and pob‐TZVP than with standard basis sets. © 2018 Wiley Periodicals, Inc.     

Consistent Gaussian basis sets of triple‐zeta valence with polarization quality for solid‐state calculations

Consistent basis sets of triple‐zeta valence with polarization quality for main group elements and transition metals from row one to three have been derived for periodic quantum‐chemical solid‐state calculations with the crystalline‐orbital program CRYSTAL. They are based on the def2‐TZVP basis sets developed for molecules by the Ahlrichs group. Orbital exponents and contraction coefficients have been modified and reoptimized, to provide robust and stable self‐consistant field (SCF) convergence for a wide range of different compounds. We compare results on crystal structures, cohesive energies, and solid‐state reaction enthalpies with the modified basis sets, denoted as pob‐TZVP, with selected standard basis sets available from the CRYSTAL basis set database. The average deviation of calculated lattice parameters obtained with a selected density functional, the hybrid method PW1PW, from experimental reference is smaller with pob‐TZVP than with standard basis sets, in particular for metallic systems. The effects of basis set expansion by diffuse and polarization functions were investigated for selected systems. © 2012 Wiley Periodicals, Inc.     

Minimally augmented Karlsruhe basis sets

We propose an extension of the basis sets proposed by Ahlrichs and coworkers at Karlsruhe (these basis sets are designated as the second-generation default or “def2” basis sets in the Turbomole program). The Karlsruhe basis sets are very appealing because they constitute balanced and economical basis sets of graded quality from partially polarized double zeta to heavily polarized quadruple zeta for all elements up to radon (Z = 86). The extension consists of adding a minimal set of diffuse functions to a subset of the elements. This yields basis sets labeled minimally augmented or with “ma” as a prefix. We find that diffuse functions are not quite as important for the def2 basis sets as they are for Pople basis sets, but they are still necessary for good results on barrier heights and electron affinities. We provide assessments and validations of this extension for a variety of data sets and representative cases. We recommend the new ma-TZVP basis set for general-purpose applications of density functional theory.     

dSH扭曲程序的建立以及在两种基组(NBO, LFMO)下对扭曲分子的能量分解

介绍了d_SH 扭曲法产生的背景, 为d_SH扭曲法研究Mills-Nixon效应提供了一个自动化程序和方法. 通过在程序中引入随机函数模拟并代替了人工调节构型参数的过程, 提高了效率并拓宽了程序的应用范围. 利用此程序算得三元苯并分子C₆B₃H₃ 和C₆O₃的一系列的扭曲结构, 并利用自然键轨道(NBO)和定域片断分子轨道(LFMO)两种基组分别进行了NBO能量分析和Morokuma SCF能量分解, 并以动态的方式分析了扭曲过程中垂直共振能及其各分量的变化趋势, 比较了不同基组和不同分析方法下的扭曲的驱动力和阻力. 由分析知, 由于NBO基组的非绝对定域性和NBO能量分析方法的一次性对角化直接导致NBO基组及其能量分析方法在d_SH扭曲能量分析中的不合理.     

Core correlating basis functions for elements 31–118

Gaussian functions for correlation of all core shells of elements from Z = 31 to Z = 118 have been optimized in relativistic singles and doubles CI calculations, performed on the shell of highest angular momentum for each principal quantum number. The SCF functions were derived from the double-zeta, triple-zeta, and quadruple-zeta basis sets previously optimized by the author. Only those Gaussian functions that are not represented in the SCF basis sets were optimized. The functions are available from the Dirac program web site, .     

吡啶甲酸分子平衡构型的ab initio理论研究

本文用能量梯度法TEXAS从头计算程序,选取4-21G(,N原子含极化d轨道)基组,对α-,β-吡啶甲酸分子的平衡几何构型进行了全优化计算,二套基组的理论优化表明,三个吡啶甲酸分子都具有平面的平衡构型,其计算键长经4-21G 基组的标准经验校正值校正后,所得的理论预测平衡几何构型与X射线单晶衍射实验结果有很好的一致性.对部分键长以及羧基与吡啶华平面的二面角中出现的差异进行了讨论.     

Analytical GIAO and hybrid-basis integral derivatives: application to geometry optimization of molecules in strong magnetic fields

Analytical integral evaluation is a central task of modern quantum chemistry. Here we present a general method for evaluating differentiated integrals over standard Gaussian and mixed Gaussian/plane-wave hybrid orbitals. The main idea is to have a representation of basis sets that is flexible enough to enable differentiated integrals to be reinterpreted as standard integrals over modified basis functions. As an illustration of the method, we report a very simple implementation of Hartree-Fock level geometrical derivatives in finite magnetic fields for gauge-origin independent atomic orbitals, within the London program. As a quantum-chemical application, we optimize the structure of helium clusters and some well-known covalently bound molecules (water, ammonia and benzene) subject to strong magnetic fields.     

Some aspects of parallelization of two-electron integrals in molecular orbital programs

Evaluation of two-electron integrals forms a substantial part of the CPU time for any ab initio molecular orbital program. This part of the package, “MICROMOL”, is parallelized. However, this parallelization leads to only sublinear speedups (typically 3 on a 4-node machine). In view of these results, the task of development of an efficient program for two-electron integrals suitable for the parallel environment has been taken up. The program is written in FORTRAN considering specific symmetry features and application of rigorous bounds. This program is further parallelized with a good load balancing strategy. The molecules used as the test cases are: trans-butadiene, benzene, nitrobenzene, naphtalene and cytosine, with 3G and 4–31G basis sets. The results indicate that the parallel version of this program gives a typical speedup of 3.6 for a 3G basis set and approximately 3.4 for a 4–31G basis set for all the molecules tested. The sequential version of this program is ～1.2 times faster than the sequential version of MICROMOL, whereas the parallel version is ～1.4 times faster than the parallelized MICROMOL.     

Mechanism for large first hyperpolarizabilities of phosphonic acid stilbene derivatives

This paper presents calculations of dipole moments (mu), static polarizabilities (alpha), and first hyperpolarizabilities (beta) of phosphonic acid stilbene derivatives calculated in the framework of density functional theory. These calculations were performed using a finite field approach implemented in the density functional program ALLCHEM and were of an all-electron type using local exchange-correlation functional and specially designed basis sets. The molecular structures have been fully optimized using the semiempirical program MSINDO. Some of the investigated stilbenes have been synthesized very recently while others are described for the first time. Donor and acceptor groups of these analogues have been modified and the influence of these changes on the first hyperpolarizabilities has been investigated. This work demonstrates that the nonlinear optical response beta of these compounds increases dramatically when the acceptor moiety is displaced by analogues containing alkali metal groups. A general mechanism for the design of novel nonlinear optical materials with large first hyperpolarizabilities is described.     

Second-order Møller-Plesset theory with linear R12 terms (MP2-R12) revisited: auxiliary basis set method and massively parallel implementation

Ab initio electronic structure approaches in which electron correlation explicitly appears have been the subject of much recent interest. Because these methods accelerate the rate of convergence of the energy and properties with respect to the size of the one-particle basis set, they promise to make accuracies of better than 1 kcal/mol computationally feasible for larger chemical systems than can be treated at present with such accuracy. The linear R12 methods of Kutzelnigg and co-workers are currently the most practical means to include explicit electron correlation. However, the application of such methods to systems of chemical interest faces severe challenges, most importantly, the still steep computational cost of such methods. Here we describe an implementation of the second-order M?ller-Plesset method with terms linear in the interelectronic distances (MP2-R12) which has a reduced computational cost due to the use of two basis sets. The use of two basis sets in MP2-R12 theory was first investigated recently by Klopper and Samson and is known as the auxiliary basis set (ABS) approach. One of the basis sets is used to describe the orbitals and another, the auxiliary basis set, is used for approximating matrix elements occurring in the exact MP2-R12 theory. We further extend the applicability of the approach by parallelizing all steps of the integral-direct MP2-R12 energy algorithm. We discuss several variants of the MP2-R12 method in the context of parallel execution and demonstrate that our implementation runs efficiently on a variety of distributed memory machines. Results of preliminary applications indicate that the two-basis (ABS) MP2-R12 approach cannot be used safely when small basis sets (such as augmented double- and triple-zeta correlation consistent basis sets) are utilized in the orbital expansion. Our results suggest that basis set reoptimization or further modifications of the explicitly correlated ansatz and/or standard approximations for matrix elements are necessary in order to make the MP2-R12 method sufficiently accurate when small orbital basis sets are used. The computer code is a part of the latest public release of Sandia's Massively Parallel Quantum Chemistry program available under GNU General Public License.     

An augmented effective core potential basis set for the calculation of molecular polarizabilities

Calculations of molecular polarizabilities require basis sets capable of accurately describing the responses of the electrons to an external perturbation. Unfortunately, basis sets that yield suitable quantitative results have traditionally been all-electron sets with large numbers of primitives, making their use computationally intractable even for moderately sized systems. We present a systematic augmentation of the effective core potential basis set of Stevens et al. [J Chem Phys 81, 12 (1984), Can J Chem 70, 612 (1992)] for 39 main group elements based on the procedure used to construct diffuse and polarization functions in the well-known Sadlej basis sets [Collec Czech Chem Comm 53, 1995 (1988)]. Representative calculations have been performed and we have shown that results to within 1% of all-electron calculations using the Sadlej basis set can be obtained for <1-35% of the computational cost using this new basis set.     

ERKALE—A flexible program package for X‐ray properties of atoms and molecules

ERKALE is a novel software program for computing X‐ray properties, such as ground‐state electron momentum densities, Compton profiles, and core and valence electron excitation spectra of atoms and molecules. The program operates at Hartree–Fock or density‐functional level of theory and supports Gaussian basis sets of arbitrary angular momentum and a wide variety of exchange‐correlation functionals. ERKALE includes modern convergence accelerators such as Broyden and ADIIS and it is suitable for general use, as calculations with thousands of basis functions can routinely be performed on desktop computers. Furthermore, ERKALE is written in an object oriented manner, making the code easy to understand and to extend to new properties while being ideal also for teaching purposes. © 2012 Wiley Periodicals, Inc.