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.     

Accurate prediction of vertical electronic transitions of Ni(II) coordination compounds via time dependent density functional theory

Time dependent density functional theory calculations are completed for five Ni(II) complexes formed by polydentate peptides to predict the electronic absorption spectrum. The ligands examined were glycyl‐glycyl‐glycine (GGG), glycyl‐glycyl‐glycyl‐glycine (GGGG), glycyl‐glycyl‐histidine (GGH), glycyl‐glycyl‐cysteine (GGC), and triethylenetetramine (trien). Fifteen functionals and two basis sets were tested. On the basis of the mean absolute percent deviation (MAPD), the ranking among the functionals is: HSE06 ∼ MPW1PW91 ∼ PBE0 > ω‐B97x‐D ∼ B3P86 ∼ B3LYP ∼ CAM‐B3LYP > PBE ∼ BLYP ∼ BP86 > TPSS > TPSSh > BHandHLYP > M06 ≫ M06‐2X. Concerning the basis sets, the triple‐ζ def2‐TZVP performs better than the double‐ζ LANL2DZ. With the functional HSE06 and basis set def2‐TZVP the MAPD with respect to the experimental λ_max is 1.65% with a standard deviation of 1.26%. The absorption electronic spectra were interpreted in terms of vertical excitations between occupied and virtual MOs based on Ni‐d atomic orbitals. The electronic structure of the Ni(II) species is also discussed.     

DSD‐PBEP86‐NL and DOD‐PBEP86‐NL functionals for noncovalent interactions: Basis set effects and tentative applications to large noncovalent systems

Basis set effects on the DSD‐PBEP86‐NL and DOD‐PBEP86‐NL functionals for noncovalent interactions have been extensively studied in this work. The cc‐pVXZ (X = D, T, Q, 5, 6) and augmented aug‐cc‐pVXZ (X = D, T, Q) basis sets are systematically tested without counterpoise (CP) corrections against the well‐known S66 database. Additionally, the basis sets of def2‐TZVPP and def2‐TZVPPD are also examined. Based on our computations, the performances of the aug‐cc‐pVQZ, cc‐pV5Z, and cc‐pV6Z basis sets are very approximate to those obtained with the def2‐QZVP basis set for both the DSD‐PBEP86‐NL and DOD‐PBEP86‐NL functionals. Note that the short‐range attenuation parameters for these two functionals were directly optimized using the def2‐QZVP basis set without CP corrections against the S66 database. Generally speaking, the cc‐pVXZ (X = D, T, Q), aug‐cc‐pVXZ (X = D, T, Q), def2‐TZVPP, and def2‐TZVPPD basis sets favor half CP correction for these two functionals. Nevertheless, the aug‐cc‐pVQZ basis set already performs well without any CP correction, especially for the DOD‐PBEP86‐NL functional. With respect to accuracy and computational cost, the cc‐pVTZ and def2‐TZVPP basis sets with half CP corrections are recommended for these two functionals to evaluate interaction energies of large noncovalent complexes.     

Accurate Coulomb-fitting basis sets for H to Rn

A series of auxiliary basis sets to fit Coulomb potentials for the elements H to Rn (except lanthanides) is presented. For each element only one auxiliary basis set is needed to approximate Coulomb energies in conjunction with orbital basis sets of split valence, triple zeta valence and quadruple zeta valence quality with errors of typically below ca. 0.15 kJ mol(-1) per atom; this was demonstrated in conjunction with the recently developed orbital basis sets of types def2-SV(P), def2-TZVP and def2-QZVPP for a large set of small molecules representing (nearly) each element in all of its common oxidation states. These auxiliary bases are slightly more than three times larger than orbital bases of split valence quality. Compared to non-approximated treatments, computation times for the Coulomb part are reduced by a factor of ca. 8 for def2-SV(P) orbital bases, ca. 25 for def2-TZVP and ca. 100 for def2-QZVPP orbital bases.     

The confirmation of accurate combination of functional and basis set for transition‐metal dimers: Fe2, Co2, Ni2, Ru2, Rh2, Pd2, Os2, Ir2, and Pt2

Eight kinds of density functionals named B3LYP, PBE1PBE, B1B95, BLYP, BP86, G96PW91, mPWPW91, and SVWN along with two different valence basis sets (LANL2DZ and CEP‐121g) are employed to study the transition‐metal dimers for the elements of group VIII. By comparing the equilibrium bond distances, vibrational frequencies, and dissociation energies of the ground state of these dimers with the available experimental values and theoretical data, we show that the “pure” DFT methods (G96PW91, BLYP, and BP86) with great‐gradient approximation always give better results relative to the hybrid HF/DFT schemes (B3LYP, PBE1PBE, and B1B95). The striking case found by us is that the G96PW91 functional, which is not tested in previous systemic studies, always predicts the dissociation energy to be well. The Ru₂ and Os₂ dimers are sensitive to not only the functionals employed but also the valence basis sets adopted. The natural bond orbital population is analyzed, and the molecular orbitals of the unpaired electrons are determined. Furthermore, our results indicate that the s and d orbitals of these dimers always hybridize with each other except for Rh₂ and Pt₂ molecules. And by analyzing the electron configuration of the bonding atom, the dissociation limit of the ground state is obtained. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

First‐principle modelling of forsterite surface properties: Accuracy of methods and basis sets

The seven main crystal surfaces of forsterite (Mg₂SiO₄) were modeled using various Gaussian‐type basis sets, and several formulations for the exchange‐correlation functional within the density functional theory (DFT). The recently developed pob‐TZVP basis set provides the best results for all properties that are strongly dependent on the accuracy of the wavefunction. Convergence on the structure and on the basis set superposition error‐corrected surface energy can be reached also with poorer basis sets. The effect of adopting different DFT functionals was assessed. All functionals give the same stability order for the various surfaces. Surfaces do not exhibit any major structural differences when optimized with different functionals, except for higher energy orientations where major rearrangements occur around the Mg sites at the surface or subsurface. When dispersions are not accounted for, all functionals provide similar surface energies. The inclusion of empirical dispersions raises the energy of all surfaces by a nearly systematic value proportional to the scaling factor s of the dispersion formulation. An estimation for the surface energy is provided through adopting C₆ coefficients that are more suitable than the standard ones to describe O? O interactions in minerals. A 2 × 2 supercell of the most stable surface (010) was optimized. No surface reconstruction was observed. The resulting structure and surface energy show no difference with respect to those obtained when using the primitive cell. This result validates the (010) surface model here adopted, that will serve as a reference for future studies on adsorption and reactivity of water and carbon dioxide at this interface. © 2015 Wiley Periodicals, Inc.     

Evaluating Shielding‐Based Ring‐Current Models by Using the Gauge‐Including Magnetically Induced Current Method

Magnetically induced current densities and integrated ring‐current strength susceptibilities have been calculated at the density functional theory (DFT) level for a test set consisting of 17 ring‐shaped molecules using the gauge‐including magnetically induced current (GIMIC) method. Reliable values for the ring‐current strengths have been obtained by performing numerical integration of the current‐density susceptibility passing a cut plane perpendicularly to the molecular ring. The current densities and ring current strengths were calculated at the DFT level using the B3LYP functional and def2‐TZVP basis sets. Current densities and ring‐current strengths have also been calculated at the Hartree‐Fock self‐consistent field (HF‐SCF) level using Dunning’s aug‐cc‐pVTZ basis sets, which allow a direct comparison with ring‐current strengths that have previously been estimated using ring‐current models based on magnetic shielding calculations. Current density calculations at both levels of theory show that the magnetic shielding based ring‐current models are not a very accurate means to estimate the magnetically induced ring current strengths, whereas they provide qualitatively the correct aromaticity trends for the studied molecules.     

Influence of the exchange‐correlation functional in all‐electron calculations of the vibrational frequencies of corundum (α‐Al2O3)

The equilibrium structural parameters, high‐ and low‐frequency dielectric tensors, Born effective charges, and Γ‐point vibrational frequencies of bulk Al₂O₃ corundum are calculated by using the periodic, ab initio program CRYSTAL, which adopts an all‐electron Gaussian‐type basis set. The effect of basis set and the performance of three different functionals, i.e., LDA, PW91, and B3LYP, are discussed. The mean absolute deviation from the measured frequencies is as small as 7 cm^?1 for both the LDA and B3LYP functionals, indicating that these functionals perform extremely well in this case. The mean absolute deviation increases to 18 cm^?1 when the PW91 functional is used. All three functionals reproduce the equilibrium geometry of corundum to a high level of accuracy, with LDA and B3LYP outperforming PW91 slightly. The comparison of the current all‐electron calculations with previous plane‐wave, pseudo‐potential calculations shows an overall similar performance. The results of isotopic substitution for both Al and O are also presented. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

A comparison between plane wave and Gaussian-type orbital basis sets for hydrogen bonded systems: formic acid as a test case

The formic acid molecule, its dimers, and its molecular crystal are adopted as test systems to compare results obtained with plane wave (PW) basis sets and norm-conserving pseudopotentials to all-electron Gaussian-type orbital (GTO) calculations. The CPMD and CRYSTAL06 codes, respectively, are applied with the PBE, PW91, and BLYP density functionals. Hydrogen bonding is the leading interaction in the dimers and the crystal. In the latter, dispersive and weak C-H...O interactions are also relevant. Irrespective of the adopted functional, for all considered structures PW and GTO results converge smoothly as a function of the quality of the adopted basis sets to the same values for structures, energies of interaction, and harmonic vibrational features. To achieve a high level of mutual agreement the use of GTO basis sets of at least of triple-zeta quality including one set of polarization functions and PW basis sets with a kinetic energy cutoff higher than 110 Ry is recommended. Pros and cons of both approaches for studying molecular crystals are also discussed.     

A DFT study on Cu(I) coordination in Cu‐ZSM‐5: Effects of the functional choice and tuning of the ONIOM approach

The coordination of Cu⁺ at the T1 and T7 positions of the M7 ring of Cu‐ZSM‐5, and the interaction of NO with coordinated Cu⁺ were investigated by means of DFT/ONIOM calculations. The B3LYP, BLYP, PBE1PBE, PBE, M06, and M062X functionals with the def2‐TZVP (def2‐QZVP for Cu) basis set were used in the high‐level part of ONIOM calculations, with the HF/3‐21G, B3LYP/LANL2DZ, M06/LANL2DZ, and M062X/LANL2DZ methods in the low‐level part. The ability of suitable combinations of the above methods to reproduce (i) the crystallographic structure of purely siliceous ZSM‐5, (ii) the tendency of Cu⁺ to be twofold or fourfold coordinated by framework oxygen atoms of Cu‐ZSM‐5, and (iii) the interaction energy and the N? O stretching frequency of adsorbed nitrogen oxide are discussed, showing that different results are obtained depending on the adopted computational approach. With reference to above properties, some considerations about the employment of the ONIOM approximations are also included. © 2015 Wiley Periodicals, Inc.     

Stacking of Metal Chelates with Benzene: Can Dispersion‐Corrected DFT Be Used to Calculate Organic–Inorganic Stacking?

CCSD(T)/CBS energies for stacking of nickel and copper chelates are calculated and used as benchmark data for evaluating the performance of dispersion‐corrected density functionals for calculating the interaction energies. The best functionals for modeling the stacking of benzene with the nickel chelate are M06HF‐D3 with the def2‐TZVP basis set, and B3LYP‐D3 with either def2‐TZVP or aug‐cc‐pVDZ basis set, whereas for copper chelate the PBE0‐D3 with def2‐TZVP basis set yielded the best results. M06L‐D3 with aug‐cc‐pVDZ gives satisfying results for both chelates. Most of the tested dispersion‐corrected density functionals do not reproduce the benchmark data for stacking of benzene with both nickel (no unpaired electrons) and copper chelate (one unpaired electron), whereas a number of these functionals perform well for interactions of organic molecules.     

Analytic projection from plane‐wave and PAW wavefunctions and application to chemical‐bonding analysis in solids

Quantum‐chemical computations of solids benefit enormously from numerically efficient plane‐wave (PW) basis sets, and together with the projector augmented‐wave (PAW) method, the latter have risen to one of the predominant standards in computational solid‐state sciences. Despite their advantages, plane waves lack local information, which makes the interpretation of local densities‐of‐states (DOS) difficult and precludes the direct use of atom‐resolved chemical bonding indicators such as the crystal orbital overlap population (COOP) and the crystal orbital Hamilton population (COHP) techniques. Recently, a number of methods have been proposed to overcome this fundamental issue, built around the concept of basis‐set projection onto a local auxiliary basis. In this work, we propose a novel computational technique toward this goal by transferring the PW/PAW wavefunctions to a properly chosen local basis using analytically derived expressions. In particular, we describe a general approach to project both PW and PAW eigenstates onto given custom orbitals, which we then exemplify at the hand of contracted multiple‐ζ Slater‐type orbitals. The validity of the method presented here is illustrated by applications to chemical textbook examples—diamond, gallium arsenide, the transition‐metal titanium—as well as nanoscale allotropes of carbon: a nanotube and the fullerene. Remarkably, the analytical approach not only recovers the total and projected electronic DOS with a high degree of confidence, but it also yields a realistic chemical‐bonding picture in the framework of the projected COHP method. © 2013 Wiley Periodicals, Inc.     

Auxiliary basis sets for density‐fitting second‐order Møller–Plesset perturbation theory: Weighted core‐valence correlation consistent basis sets for the 4d elements Y–Pd

Auxiliary basis sets (ABS) specifically matched to the cc‐pwCVnZ‐PP and aug‐cc‐pwCVnZ‐PP orbital basis sets (OBS) have been developed and optimized for the 4d elements Y‐Pd at the second‐order Møller‐Plesset perturbation theory level. Calculation of the core‐valence electron correlation energies for small to medium sized transition metal complexes demonstrates that the error due to the use of these new sets in density fitting is three to four orders of magnitude smaller than that due to the OBS incompleteness, and hence is considered negligible. Utilizing the ABSs in the resolution‐of‐the‐identity component of explicitly correlated calculations is also investigated, where it is shown that i‐type functions are important to produce well‐controlled errors in both integrals and correlation energy. Benchmarking at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations level indicates impressive convergence with respect to basis set size for the spectroscopic constants of 4d monofluorides; explicitly correlated double‐ζ calculations produce results close to conventional quadruple‐ζ, and triple‐ζ is within chemical accuracy of the complete basis set limit. © 2013 Wiley Periodicals, Inc.     

HF and MP2 calculations on CN−, N2, AlF,SiO, PN,SC, ClB,and P2 using correlated molecular wave functions

Contracted basis sets of double zeta valence quality plus polarization functions (DZP) and augmented DZP basis sets, which were recently constructed for the first‐ and second‐row atoms, are applied to study the electronic ground states of the diatomic molecules CN^?, N₂, AlF, SiO, PN, SC, ClB, and P₂. At the Hartree–Fock (HF) and/or Møller–Plesset second‐order (MP2) levels, total and molecular orbital energies, dissociation energies, bond lengths, harmonic vibrational frequencies, and dipole moments are calculated and compared with available experimental data and with the results obtained from correlation consistent polarized valence basis sets of Dunning's group. For N₂, calculations of polarizabilities at the HF and MP2 levels with the sets presented above are also done and compared with results reported in the literature. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Ab initio investigation of ground‐state properties of group‐12 fluorides

The performance of wavefunction‐based correlation methods in theoretical solid‐state chemistry depends on reliable Hartree–Fock (HF) results for infinitly extended systems. Therefore, we optimized basis sets of valence‐triple‐ζ quality based on HF calculations for the periodic system of group‐12‐metal difluorides. Scalar‐relativistic effects were included in the case of the metal‐ions by applying small‐core pseudopotentials. To assess the quality of the proposed basis sets, the structural parameters, bulk moduli as well as cohesive and lattice energies of the systems were evaluated at the HF and the density functional theory levels. In addition to these two mean‐field approaches and to assess further employment of our basis sets to wavefunction‐based correlation methods we performed periodic local MP2 computations. Finally, the possibilities of pressure induced structural phase transitions occurring in the ZnF₂, CdF₂, and HgF₂ were investigated. © 2014 Wiley Periodicals, Inc.     

MPI/OpenMP hybrid parallel algorithm for resolution of identity second‐order Møller–Plesset perturbation calculation of analytical energy gradient for massively parallel multicore supercomputers

A massively parallel algorithm of the analytical energy gradient calculations based the resolution of identity Møller–Plesset perturbation (RI‐MP2) method from the restricted Hartree–Fock reference is presented for geometry optimization calculations and one‐electron property calculations of large molecules. This algorithm is designed for massively parallel computation on multicore supercomputers applying the Message Passing Interface (MPI) and Open Multi‐Processing (OpenMP) hybrid parallel programming model. In this algorithm, the two‐dimensional hierarchical MP2 parallelization scheme is applied using a huge number of MPI processes (more than 1000 MPI processes) for acceleration of the computationally demanding O (N ⁵) step such as calculations of occupied–occupied and virtual–virtual blocks of MP2 one‐particle density matrix and MP2 two‐particle density matrices. The new parallel algorithm performance is assessed using test calculations of several large molecules such as buckycatcher C₆₀@C₆₀H₂₈ (144 atoms, 1820 atomic orbitals (AOs) for def2‐SVP basis set, and 3888 AOs for def2‐TZVP), nanographene dimer (C₉₆H₂₄)₂ (240 atoms, 2928 AOs for def2‐SVP, and 6432 AOs for cc‐pVTZ), and trp‐cage protein 1L2Y (304 atoms and 2906 AOs for def2‐SVP) using up to 32,768 nodes and 262,144 central processing unit (CPU) cores of the K computer. The results of geometry optimization calculations of trp‐cage protein 1L2Y at the RI‐MP2/def2‐SVP level using the 3072 nodes and 24,576 cores of the K computer are presented and discussed to assess the efficiency of the proposed algorithm. © 2017 Wiley Periodicals, Inc.     

Orbital signatures of methyl in L-alanine

Molecular orbital signatures of the methyl substituent in L-alanine have been identified with respect to those of glycine from information obtained in coordinate and momentum space, using dual space analysis. Electronic structural information in coordinate space is obtained using ab initio (MP2/TZVP) and density functional theory (B3LYP/TZVP) methods, from which the Dyson orbitals are simulated based on the plane wave impulse approximation into momentum space. In comparison to glycine, relaxation in geometry and valence orbitals in L-alanine is found as a result of the attachment of the methyl group. Five orbitals rather than four orbitals are identified as methyl signatures. That is, orbital 6a in the core shell, orbitals 11a and 12a in the inner valence shell, and orbitals 19a and 20a in the outer valence shell. In the inner valence shell, the attachment of methyl to glycine causes a splitting of its orbital 10a' into orbitals 11a and 12a of L-alanine, whereas in the outer valence shell the methyl group results in an insertion of an additional orbital pair of 19a and 20a. The frontier molecular orbitals, 24a and 23a, are found without any significant role in the methylation of glycine.     

All‐electron LCAO calculations of the LiF crystal phonon spectrum: Influence of the basis set,the exchange‐correlation functional,and the supercell size

For the first time the convergence of the phonon frequencies and dispersion curves in terms of the supercell size is studied in ab initio frozen phonon calculations on LiF crystal. Helmann–Feynman forces over atomic displacements are found in all‐electron calculations with the localized atomic functions (LCAO) basis using CRYSTAL06 program. The Parlinski–Li–Kawazoe method and FROPHO program are used to calculate the dynamical matrix and phonon frequencies of the supercells. For fcc lattice, it is demonstrated that use of the full supercell space group (including the supercell inner translations) enables to reduce essentially the number of the displacements under consideration. For Hartree–Fock (HF), PBE and hybrid PBE0, B3LYP, and B3PW exchange‐correlation functionals the atomic basis set optimization is performed. The supercells up to 216 atoms (3 × 3 × 3 conventional unit cells) are considered. The phonon frequencies using the supercells of different size and shape are compared. For the commensurate with supercell k ‐points the best agreement of the theoretical results with the experimental data is found for B3PW exchange‐correlation functional calculations with the optimized basis set. The phonon frequencies at the most non‐commensurate k ‐points converged for the supercell consisting of 4 × 4 × 4 primitive cells and ensures the accuracy 1–2% in the thermodynamic properties calculated (the Helmholtz free energy, entropy, and heat capacity at the room temperature). © 2009 Wiley Periodicals, Inc. J Comput Chem 2009