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     

Reassessment of methyl rotation barriers and conformations by correlated quantum chemistry methods

Internal rotations of the methyl group in ortho‐substituted and 2,6‐disubstituted toluenes in their ground state have been investigated by means of various ab initio quantum chemistry methods. Computed barriers at the Hartree‐Fock (HF) level using medium sized basis sets agreed reasonably with experimental results in the case of the studied ortho‐substituted toluenes. However, this agreement worsens when using very large basis sets. Furthermore, the determination of the conformation and barriers of more weakly hindered methyl groups, that is, for 2,6‐dihalogenotoluenes or toluene itself, necessitates high level correlated computations, because of a possible failure of HF calculations in this case. Density functional theory (DFT) techniques required, in several cases, much more extended basis sets than the post‐HF Møller‐Plesset perturbation (MP2, MP4) ones, to insure the convergence of the computed barriers. Non‐negligible variations of the computed barriers when using different DFT functionals are observed for some systems. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 2093–2100, 2003     

New basis sets for the evaluation of interaction‐induced electric properties in hydrogen‐bonded complexes

Interaction‐induced static electric properties, that is, dipole moment, polarizability, and first hyperpolarizability, of the CO? (HF)_n and N₂? (HF)_n, n = 1–9 hydrogen‐bonded complexes are evaluated within the finite field approach using the Hartree–Fock, density functional theory, Møller–Plesset second‐order perturbation theory, and coupled cluster methods, and the LPol‐n (n = ds, dl, fs, fl) basis sets. To compare the performance of the different methods with respect to the increase of the complex size, we consider as model systems linear chains of the complexes. We analyze the results in terms of the many‐body and cooperative effects. © 2012 Wiley Periodicals, Inc.     

BSSE‐free description of the formamide dimers

The different configurations (linear, zig‐zag, and cyclic) of formamide dimers have been studied at the level of both Hartree–Fock (HF) and second order Møller–Plesset perturbation theory (MP2). The widely used a posteriori Boys–Bernardi “counterpoise” (CP) correction scheme has been compared with our a priori methods utilizing the “chemical Hamiltonian approach” (CHA). The appropriate interaction energies have been calculated in six different basis sets (6‐31G, 6‐31G**, DZV, DZP, TZV, and cc‐pVDZ). © 2001 John Wiley & Sons, Inc. Int J Quant Chem, 2001     

3He NMR studies on helium–pyrrole,helium–indole,and helium–carbazole systems: a new tool for following chemistry of heterocyclic compounds

The ³He nuclear magnetic shieldings were calculated for free helium atom and He–pyrrole, He–indole, and He–carbazole complexes. Several levels of theory, including Hartree–Fock (HF), Second‐order Møller‐Plesset Perturbation Theory (MP2), and Density Functional Theory (DFT) (VSXC, M062X, APFD, BHandHLYP, and mPW1PW91), combined with polarization‐consistent pcS‐2 and aug‐pcS‐2 basis sets were employed. Gauge‐including atomic orbital (GIAO) calculated ³He nuclear magnetic shieldings reproduced accurately previously reported theoretical values for helium gas. ³He nuclear magnetic shieldings and energy changes as result of single helium atom approaching to the five‐membered ring of pyrrole, indole, and carbazole were tested. It was observed that ³He NMR parameters of single helium atom, calculated at various levels of theory (HF, MP2, and DFT) are sensitive to the presence of heteroatomic rings. The helium atom was insensitive to the studied molecules at distances above 5 Å. Our results, obtained with BHandHLYP method, predicted fairly accurately the He–pyrrole plane separation of 3.15 Å (close to 3.24 Å, calculated by MP2) and yielded a sizable ³He NMR chemical shift (about ?1.5 ppm). The changes of calculated nucleus‐independent chemical shifts (NICS) with the distance above the rings showed a very similar pattern to helium‐3 NMR chemical shift. The ring currents above the five‐membered rings were seen by helium magnetic probe to about 5 Å above the ring planes verified by the calculated NICS index. Copyright © 2014 John Wiley & Sons, Ltd.     

Theoretical study on the molecular structures of toluene,para‐fluorotoluene,para‐chlorotoluene,and 4‐methylpyridine and their sixfold internal rotational barriers

Two kinds of sixfold internal rotational configurations of toluene, para‐fluorotoluene, para‐chlorotoluene, and 4‐methylpyridine were calculated using Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and Beck's three parameter hybrid functional using the LYP correlation functional (B3LYP) theory methods with various high‐level basis sets. Structures and energies were compared for different configurations. Calculations indicate that the orthogonal configuration has a local minimum while the planar configuration is a transition structure. Furthermore, geometries of the orthogonal and the planar configurations are quite similar, except for a methyl CH bond. Sixfold internal rotational barriers were calculated from the energy difference of two different configurations. For the calculated results, HF methods underestimated the rotational barriers, but MP2 calculations overestimated them. However, the density functional theory (DFT) method is a reliable method since the calculated internal rotational barriers are similar to the experimental ones. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 772–778, 2000     

Encapsulation of the interstellar abundant H3+ in a C60 fullerene

The possible encapsulation of the interstellar abundant H₃⁺ ion inside a C₆₀ fullerene cage has been examined by using the Hartree‐Fock (HF) and the second order Møller‐Plesset perturbation (MP2) methods both with the 6‐31G** basis set. It was found that H₃⁺ forms various stable endohedral complexes inside the cage. Six configurations have been examined among which four were stable compared with the separated initial species, the dissociated H₂ + H⁺ inside the cage being the most stable. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical study of silicon–sulfur clusters (SiS2)n− (n=1–6)

The possible geometrical structures and relative stability of silicon–sulfur clusters (SiS₂) (n=1–6) are explored by means of density functional theory (DFT) quantum chemical calculations. We also compare DFT with second‐order Møller–Plesset (MP2) and Hartree–Fock (HF) methods. The effects of polarization functions, diffuse functions, and electron correlation are included in MP2 and B3LYP quantum chemical calculations, and B3LYP is effective in larger cluster structure optimization, so we can conclude that the DFT approach is useful in establishing trends. The electronic structures and vibrational spectra of the most stable geometrical structures of (SiS₂)_n⁻ are analyzed by B3LYP. As a result, the regularity of the (SiS₂)_n⁻ cluster growing is obtained, and the calculation may predict the formation mechanism of the (SiS₂)_n⁻ cluster. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 280–290, 2001     

Analytical gradients of the second‐order Møller–Plesset energy using Cholesky decompositions

An algorithm for computing analytical gradients of the second‐order Møller–Plesset (MP2) energy using density fitting (DF) is presented. The algorithm assumes that the underlying canonical Hartree–Fock reference is obtained with the same auxiliary basis set, which we obtain by Cholesky decomposition (CD) of atomic electron repulsion integrals. CD is also used for the negative semidefinite MP2 amplitude matrix. Test calculations on the weakly interacting dimers of the S22 test set (Jure?ka et al., Phys. Chem. Chem. Phys. 2006, 8, 1985) show that the geometry errors due to the auxiliary basis set are negligible. With double‐zeta basis sets, the error due to the DF approximation in intermolecular bond lengths is better than 0.1 pm. The computational time is typically reduced by a factor of 6–7. © 2013 Wiley Periodicals, Inc.     

Evaluation of DFT methods for computing the interaction energies of homomolecular and heteromolecular dimers of monosubstituted benzene

We present density functional theory (DFT) interaction energies for the sandwich and T‐shaped conformers of substituted benzene dimers. The DFT functionals studied include TPSS, HCTH407, B3LYP, and X3LYP. We also include Hartree–Fock (HF) and second‐order Møller–Plesset perturbation theory calculations (MP2), as well as calculations using a new functional, P3LYP, which includes PBE and HF exchange and LYP correlation. Although DFT methods do not explicitly account for the dispersion interactions important in the benzene–dimer interactions, we find that our new method, P3LYP, as well as HCTH407 and TPSS, match MP2 and CCSD(T) calculations much better than the hybrid methods B3LYP and X3LYP methods do. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Møller–Plesset theory for atomic ground state energies

Mø–Plesset theory, in which electron correlation energy is calculated by perturbation techniques, is used in second order to calculate energies of the ground states of atoms up to neon. The unrestricted Hartree–Fock (UHF) Hamiltonian is used as the unperturbed system and the technique is then described as unrestricted Mø–Plesset to second order (UMP2). Use of large Gaussian basis sets suggests that the limiting UMP2 energies with a complete basis of s, p, and d functions account for 75–84% of the correlation energy. Preliminary estimates of the contributions of basis functions with higher angular quantum numbers indicate that full UMP2 limits give even more accurate total energies.     

Ab initio pseudopotential study of M2As− and M2Br+ (M = Cu,Ag, Au)

The equilibrium geometries and the vibration frequencies of M₂As⁻ and M₂Br⁺ (M = Cu, Ag, Au) are calculated at the Hartree–Fock (HF) and the second‐order Møller–Plesset (MP2) levels with pseudopotentials. The calculated results indicate that the species have a bent structure (C_2v). The electron correlation corrections on the geometrical structure are investigated at the MP2 level, the bond angles are reduced by 10°–20° for considered species. The electron correlation effects on the geometry of the Au₂As⁻ are studied particularly at MP2, MP3, MP4, CCSD and CCSD(T) levels. Comparing the species containing Ag and Au, the relativistic effects slightly short the bond lengths of the species. The bonding possibility of the Au₂As⁻ is predicted. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 80: 38–43, 2000     

Theoretical study on intermolecular interaction of epoxyethane dimer

Ab initio calculations at Hartree–Fock and fourth‐order Mø ller–Plesset (MP4) correlation correction levels with 6‐31G* basis set have been performed on the epoxyethane dimer. Dimer binding energies have been corrected for the basis set superposition error (BSSE) and the zero‐point energy. The greatest corrected dimer binding energy is −8.36 kJ/mol at the MP4/6‐31G*//HF/6‐31G* level. The natural bond orbital analysis has been performed to trace the origin of the weak interactions that stabilize dimer. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 94–98, 2000     

Ab initio study of Rg2I− (Rg = Ar,Kr, Xe)

The equilibrium geometries, vibrational frequencies, and dissociation energies of rare gas iodine clusters Rg₂I^?(Rg = Ar, Kr, Xe) were calculated at the Hartree–Fock (HF), second‐order Møller–Plesset (MP2), the coupled cluster method with single and double excitation and a noniterative correction for triple excitations method [CCSD(T)] levels. The title species have bent C_2v structure of about 60° angle. The electron correlation effects and relativistic effects on the geometry and stability were investigated at CCSD(T) level. Both effects stabilize title species. The calculated electron affinities are in good agreement with the experimental results available. The effect of high angular momentum functions (g and h) was studied. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical Studies on the Mechanistic Insertions of Singlet Methylene and Halomethylene into Polar S‒H Bonds of Methanethiol,Ethanethiol, 1-Propanethiol and 2-Propanethiol

Abstract

The insertion of ¹CH₂, ¹CHF, ¹CF₂,¹CHCl, and ¹CCl₂ into primary and secondary polar S?H bonds of methanethiol, ethanethiol, 1-propanethiol, and 2-propanethiol have been investigated at HF (Hartree–Fock), MP2 (Møller–Plesset Perturbation Theory), and DFT (density functional theory) levels using 6-31G (d, p) basis set. The insertions follow a two-step mechanism. The potential energy surface exploration identifies ylide-like structures undergoing the 1,2-sigmatropic hydrogen shift concertedly giving the thioether product. But ¹CF₂ forms weak complexes involving 1,2-hydrogen shift. The barrier height in the concerted mechanism varies with the type of carbene moiety and S?H bond. The initial interaction seems to be a function of HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energies of thiols and carbenes, respectively. The natural bond orbital charge analysis showing a net charge flow from the substrate thiol to the reagent carbene at the transition states corroborates this fact.     

Triazene proton affinities: A comparison between density functional,Hartree–Fock,and post-Hartree–Fock methods

The consistency of three density functional computational implementations (DMol, DGauss, and deMon) are compared with high-level Hartree–Fock and Møller–Plesset (MP) calculations for triazene (HN?NNH₂) and formyl triazene (HN?NNHCOH). Proton affinities on all electronegative sites are investigated as well as the geometries of the neutral and protonated species. Density functional calculations employing the nonlocal gradient corrections show agreement with MP calculations for both proton affinities and geometries of neutral and protonated triazenes. Local spin density approximation DMol calculations using numerical basis sets must employ an extended basis to agree with other density functional codes using analytic Gaussian basis sets. The lowest energy conformation of triazene was found to be nonplanar; however, the degree of nonplanarity, as well as some bond lengths, is dependent on the basis set, electron correlation treatment, and methods used for the calculation. © 1994 by John Wiley & Sons, Inc.

1 This article is a U.S. Government work and, as such, is in the public domain in the United States of America.

     

Structural Characteristics of Graphane‐Type C and BN Nanostructures by Periodic Local MP2 Approach

The structural characteristics of fully‐hydrogenated carbon and boron nitride mono‐ and multilayer slabs, together with nanotubes derived from the slabs, are investigated mainly by means of periodic local second‐order Møller–Plesset perturbation (LMP2) calculations and the results are compared with Hartree–Fock (HF), density functional theory (DFT), and dispersion function‐augmented DFT (DFT‐D) obtained ones. The investigated systems are structurally analogous to (111) and (110) slabs of diamond, where the hydrogenated (111) slab of diamond corresponds to the experimentally known graphane. Multilayering of monolayers and nanotubes is energetically favorable at the LMP2 level for both C and BN, while HF and DFT are not able to reproduce this behavior for CH systems. The work highlights the importance of utilizing methods capable of properly describing weak interactions in the investigation of dispersively‐bound systems such as the multilayered graphanes and the corresponding nanotubes.     

Multiply Borylated Arenes: X‐ray Crystal Structure Analyses and Quantum Chemical Calculations

Optimised synthesis procedures and results of X‐ray single crystal structure analyses for 4‐(dibromoboryl)toluene, 1, 3‐bis(dibromoboryl)benzene, 1, 4‐bis(dibromoboryl)benzene, and 1, 3, 5‐tris(dibromoboryl)benzene are reported. These compounds have also been studied by Hartree‐Fock (HF), density functional theory (DFT), and Mßller‐Plesset second‐order perturbation (MP2) methods in combination with the polarized double‐ζ valence (SVP) and polarized triple‐ζ valence (TZVP) basis sets of Ahlrichs and coworkers. A comparison of the quantum chemical results for optimised geometries and computed NMR chemical shifts with experiment is presented to test the quality of the various methods for this class of compounds. All DFT methods tested yield optimised geometries within the experimental error bars of 3σ for bond lengths, whereas larger deviations among the methods are observed for computed NMR chemical shifts. This calibration recommends the B3LYP/SVP combination as a reliable and computationally efficient level of theory to assess the structures and absolute and relative ¹H‐, ¹³C‐ and ¹¹B NMR shift values of borylated aromatic compounds in future investigations.     

An approximate potential energy surface for HeI2 collisions

We have calculated the interaction potential for HeI₂ in T-shaped geometries using Hartree—Fock and Møller—Plesset third-order perturbation t