Correlated polarization propagator calculations of static polarizabilities

We present a systematic comparison of the correlation contribution at the level of the second-order polarization propagator approximation (SOPPA ) and MP 2 to the static dipole polarizability of (1) Be, BeH^?, BH, CH⁺, MgH^?, AIH, SiH⁺, and GeH⁺; (2) BH₃, CH₄, NH₃, H₂O, HF, BF, and F₂; and (3) N₂, CO, CN^?, HCN, C₂H₂, and HCHO . Fairly extended basis sets were used in the calculations. We find that the agreement with experimental values is improved in SOPPA and MP .2 over the results at the SCF level. The signs and magnitudes of the correlation contribution in SOPPA are similar to those obtained in analytical derivative MP 2 calculations. However, it is not possible to say, in general, which method gives the largest correlation contribution or the best agreement with experiment, nor is it possible to make a priori prediction of the sign of the correlation contribution. For the first group of molecules, which have a quasi-degenerate ground state, additional CCDPPA and CCSDPPA calculations were performed and compared with polarizabilities obtained as analytical/numerical derivatives of the CCD and CCSD energies. The CCSDPPA results were found to be in better agreement with other calculations than were the SOPPA results, demonstrating the necessity of using methods based on infinite-order perturbation theory for these systems. © 1994 John Wiley & Sons, Inc.     

Estimation of isotropic nuclear magnetic shieldings in the CCSD(T) and MP2 complete basis set limit using affordable correlation calculations

A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH₂O, H₂O, HF, F₂, HCN, SiH₄ and H₂S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T)), with affordable pcS‐2 basis set and corresponding complete basis set results, estimated from calculations with the family of polarization‐consistent pcS‐n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and shieldings obtained with the significantly smaller basis sets pcS‐2 and aug‐cc‐pVTZ‐J for the selected set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS‐2, MP2/pcS‐2 and DFT/CBS calculations with pcS‐n basis sets. The proposed method leads to a fairly accurate estimation of nuclear magnetic shieldings and considerable saving of computational efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

Second-order polarization propagator calculations of dynamic dipole polarizabilities and C6 coefficients

The frequency-dependent dipole polarizability, α(E), is calculated using the second-order polarization propagator approximation (SOPPA ). We have shown how to express α(E) as a function of E² and thus obtained a form of α(E) that can be used to compute C₆-coefficients without invoking complex arithmetic. For He we find that SOPPA recovers a large fraction of the correlation contribution for all frequencies, whereas for H₂ where the correlation contributions are much smaller and also basis set-dependent, we find a less definite trend of SOPPA relative to RPA .     

Correlated calculations of indirect nuclear spin-spin coupling constants using second-order polarization propagator approximations: SOPPA and SOPPA(CCSD)

We present correlated calculations of the indirect nuclear spin-spin coupling constants of HD, HF, H₂O, CH₄, C₂H₂, BH, AlH, CO and N₂ at the level of the second-order polarization propagator approximation (SOPPA) and the second-order polarization propagator approximation with coupled-cluster singles and doubles amplitudes – SOPPA(CCSD). Attention is given to the effect of the so-called W ₄ term, which has not been included in previous SOPPA spin-spin coupling constant studies of these molecules. Large sets of Gaussian basis functions, optimized for the calculation of indirect nuclear spin-spin coupling constants, were used instead of the in general rather small basis sets used in previous studies. We find that for nearly all couplings the SOPPA(CCSD) method performs better than SOPPA. Received: 6 July 1998 / Accepted: 8 September 1998 / Published online: 23 November 1998     

Molecular orbital studies of the NMR hydrogen bond shift

Magnetic shielding constants are calculated for the protons in XOH and XOH…OH₂ (XH, CH₃, NH₂, OH and F) molecules using a slightly extended set of atomic functions modified by gauge factors. These results are used to determine theoretical values for the NMR hydrogen bond shifts in the XOH…OH₂ systems. Such theoretical data are consistent with the few available experimental data. An analysis of the theoretical results reveals that there are three major types of shielding contribution to the NMR hydrogen bond shift; (a) a deshielding change due to the variation of the local currents on the hydrogen bonded proton; (b) a reduction in shielding from currents localized on the oxygen atom of the proton donor; (c) a deshielding contribution from currents induced on the oxygen atom of the proton acceptor. Except for the water dimer, contributions (a), (b) and (c) are of comparable importance for changes in isotropic shielding. For (H₂O)₂ contributions (a) and (c) are somewhat more important than contribution (b). Contribution (c) is almost totally responsible for the changes in the anistropies of the shielding tensors associated with the hydrogen bonded protons. The proton shielding anisotropy changes which occur on hydrogen bond formation are generally much larger than the corresponding variations in the isotropic values of the shielding tensors. This suggests that proton magnetic shielding anisotropies may be more sensitive measures of features of hydrogen bonding than are isotropic proton shielding constants.     

Calculation of magnetic properties of HF,H2O,NH3, and CH4 molecules using a longitudinal gauge for the vector potential

A transformation of the transverse Coulomb vector potential was implemented to calculate molecular magnetic properties via the random-phase approximation (RPA) within the framework of a “longitudinal gauge.” In this gauge, the diamagnetic contribution to magnetic susceptibility is a tensor with equal diagonal components as in atoms, irrespective of molecular symmetry, whereas diagonal and average diamagnetic contributions to the nuclear magnetic shielding are the same as in the Coulomb gauge. Near-Hartree–Fock magnetic susceptibility and nuclear magnetic shielding tensors were evaluated for a set of small molecules, HF, H₂O, NH₃, and CH₄, employing extended Gaussian basis sets. The peculiar features of the longitudinal gauge, and the fulfillment of a series of sum rules involving the virial operator, which must be satisfied to guarantee gauge invariance of total magnetic tensors, were exploited to check the degree of convergence of theoretical values and to estimate the corresponding Hartree–Fock limit for the properties. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 31–45, 1998     

Study on Pair Correlation Energies of Isoelectronic Systems F^－, HF and H2F^＋

The intrapair and interpair correlation energies of F-, HF and H2F^ systems are calculated and analyzed using MP2-OPT2 method of MELD program with cc-PVSZ^* basis set. From the analysis of pair correlation energies of these isoelectronlc sysoterns, it is found that the 1sF^2 pair correlation energy is trans-ferable in these three isociectronic systems. According to the definition of pair correlation contribution of one electron pair to a system, the pair correlation contribution values of these three systems are calculated. The correlation contribution values of inner electron pairs and H—F bonding electron pair in HF molecule with those in H2F^ system are compared. The results indicate that the bonding effect of a molecule is one of the im-portant factors to influence electron correlation energy of the system. The comparison of correlation energy contributions in-cluding triple and quadruple excitations with those only includ-ing singles and doubles calculated with 6-311 G(d) basis set shows that the higher.excitation correlation energy contribution gives more than 2 % of the total correlation energy for these sys-tems.     

Unexpected geometrical effects on paramagnetic spin-orbit and spin-dipolar 2J(FF) couplings

The second-rank tensor character of the paramagnetic spin-orbit and spin-dipolar contributions to nuclear spin-spin coupling constants is usually ignored when NMR measurements are carried out in the isotropic phase. However, in this study it is shown that isotropic (2)J(FF) couplings strongly depend on the relative orientation of the C-F bonds containing the coupling nuclei and the eigenvectors of such tensors. Predictions about such effect are obtained using a qualitative approach based on the polarization propagator formalism at the RPA, and results are corroborated performing high-level ab initio spin-spin coupling calculations at the SOPPA(CCSD)/EPR-III//MP2/EPR-III level in a model system. It is highlighted that no calculations at the RPA level were carried out in this work. The quite promising results reported in this paper suggest that similar properties are expected to hold for the second-rank nuclear magnetic shielding tensor.     

Calculation of optimum geometries and force fields by the CNDO/force method

CNDO/Force calculations have been performed on a series of molecules, H₂CO, F₂CO, CF₄, CHF₃, CH₂F₂ and CH₃F. The optimum geometries and force fields are reported. It is found that the method can successfully predict the geometries of polyatomic molecules. The bending force constants and interaction force constants are, in general, comparable with experimental values both with respect to sign and magnitude. The stretching force constants have higher values than the experimental force constants. However, the trend in stretching force constants of a series of molecules is comparable with that of the corresponding experimental values.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants: 1. Three‐membered rings

The carbon–carbon indirect nuclear spin–spin coupling constants in cyclopropane, aziridine and oxirane were investigated by means of ab initio calculations at the RPA, SOPPA and DFT/B3LYP levels. We found that the carbon–carbon couplings are by far dominated by the Fermi contact term. Our best SOPPA and DFT results are in a very good agreement with each other and with the experimental values, whereas calculations at the RPA level of theory strongly overestimate the carbon–carbon couplings. Significant differences in the basis set dependence of the calculated carbon–carbon coupling constants obtained with either wavefunction method, RPA or SOPPA, or the density functional method, DFT/B3LYP, are observed. The SOPPA results depend much more strongly on the quality of the basis set than the results of DFT/B3LYP calculations. The medium‐sized core‐valence basis sets cc‐pCVTZ and even cc‐pCVDZ were found to perform fairly well at the SOPPA level for the one‐bond carbon–carbon couplings investigated here. Copyright © 2002 John Wiley & Sons, Ltd.     

The analysis of NMR J-couplings of saturated and unsaturated compounds by the localized second order polarization propagator approach method

Calculations of NMR J-coupling with polarization propagators are not invariant under unitary transformations at second order level of approach, second order polarization propagator approach (SOPPA). They are only invariant at first order or random phase level of approach (RPA). We performed "localized" SOPPA (Loc-SOPPA), calculations of J-couplings applying two different schemes for the localization of molecular orbitals(LMO): Foster-Boys and Pipek-Mezey. We show here that results of such Loc-SOPPA calculations are different though not much: they are less than 6% different in the worst case. Therefore it is possible to apply them with confidence in the analysis of the transmission of different coupling mechanisms within the molecule. We are able now to get reliable information on what LMOs are the most important (and so which are not important) for a given J-coupling in a molecule. This information can then be used for selecting which are the paths that should be described with the highest possible accuracy for that J-coupling calculation. A few unsaturated compounds are analyzed: ethene, trans-difluoroethene or DiF-ethene, and imine. It is shown that different lone pairs (of p(z) or p(x/y) type) are responsible for the vicinal F-F J-coupling in DiF-ethene; and also the fact that the main LP contributor is not the same for the fermi contact and the spin-dipolar mechanisms. We also studied phosphorous containing compounds such as phosphine and cis-propylene phosphine. In both cases the analysis of the main LMO contributing to one-bond P-H coupling and through-space P-C coupling were performed. The above mentioned unsaturated molecular systems have quasiinstability problems that arise at RPA level of approach. We show here that they are mostly originated in the antibonding π? LMO, corresponding to the C=C or C=N double bonds. We performed the analysis of the origin of quasiinstabilities for the SD mechanism. The contribution of each kind of excitation terms to SOPPA calculations were considered, meaning the main contributions by single and double excitations. It is shown that one can get more than 97% of the total electron correlation contribution when including terms that mainly contain single excitations (though double-excitation matrix elements should still be calculated).     

Molecular Vibrations and Mean Square Amplitudes. Supplement III Mean Amplitudes of Vibration for Inorganic Six-Atomic Molecules

Review of mean amplitudes of vibration for inorganic six-atomic molecules from spectroscopic calculations and electrondiffraction-data. The work contains new results of spectroscopic mean amplitudes and force constants for a number of molecules, viz.: S₆, N₂O₄, B₂F₄, N₂H₄, N₂F₄, P₂I₄ and SOF₄.     

From CCSD(T)/aug-cc-pVTZ-J to CCSD(T) complete basis set limit isotropic nuclear magnetic shieldings via affordable DFT/CBS calculations

It is shown that a linear correlation exists between nuclear shielding constants for nine small inorganic and organic molecules (N(2), CO, CO(2), NH(3), CH(4), C(2)H(2), C(2)H(4), C(2)H(6) and C(6)H(6)) calculated with 47 methods (42 DFT methods, RHF, MP2, SOPPA, SOPPA(CCSD), CCSD(T)) and the aug-cc-pVTZ-J basis set and corresponding complete basis set results, estimated from calculations with the family of polarization-consistent pcS-n basis sets. This implies that the remaining basis set error of the aug-cc-pVTZ-J basis set is very similar in DFT and CCSD(T) calculations. As the aug-cc-pVTZ-J basis set is significantly smaller, CCSD(T)/aug-cc-pVTZ-J calculations allow in combination with affordable DFT/pcS-n complete basis set calculations the prediction of nuclear shieldings at the CCSD(T) level of nearly similar accuracy as those, obtained by fitting results obtained from computationally demanding pcS-n calculations at the CCSD(T) limit. A significant saving of computational efforts can thus be achieved by scaling inexpensive CCSD(T)/aug-cc-pVTZ-J calculations of nuclear isotropic shieldings with affordable DFT complete basis set limit corrections.     

Gauge invariant calculations of nuclear magnetic shielding constants using the continuous transformation of the origin of the current density approach. II. Density functional and coupled cluster theory

The quantum mechanical current density induced in a molecule by an external magnetic field is invariant to translations of the coordinate system. This fundamental symmetry is exploited to formally annihilate the diamagnetic contribution to the current density via the approach of "continuous transformation of the origin of the current density-diamagnetic zero" (CTOCD-DZ). The relationships obtained by this method for the magnetic shielding at the nuclei are intrinsically independent of the origin of the coordinate system for any approximate computational scheme relying on the algebraic approximation. The authors report for the first time an extended series of origin-independent estimates of nuclear magnetic shielding constants using the CTOCD-DZ approach at the level of density functional theory (DFT) with four different types of functionals and unrelaxed coupled cluster singles and doubles linear response (CCSD-LR) theory. The results obtained indicate that in the case of DFT the procedure employed is competitive with currently adopted computational methods allowing for basis sets of gauge-including atomic orbitals, whereas larger differences between CTOCD-DZ and common origin CCSD-LR results are observed due to the incomplete fulfillment of hypervirial relations in standard CCSD-LR theory. It was found furthermore that the unrelaxed CCSD-LR calculations predict larger correlation corrections for the shielding constants of almost all nonhydrogen atoms in their set of molecules than the usual relaxed energy derivative CCSD calculations. Finally the results confirm the excellent performance of Keal and Tozer's third functional, in particular, for the multiply bonded systems with a lot of electron correlation, but find also that the simple local density functional gives even better results for the few singly bonded molecules in their study where correlation effects are small.     

Nonempirical calculations of the one‐bond 29Si–13C spin–spin coupling constants taking into account relativistic and solvent corrections

The computational study of the one‐bond ²⁹Si–¹³C spin–spin coupling constants has been performed at the second‐order polarization propagator approximation (SOPPA) level in the series of 60 diverse silanes with a special focus on the main factors affecting the accuracy of the calculation including the level of theory, the quality of the basis set, and the contribution of solvent and relativistic effects. Among three SOPPA‐based methods, SOPPA(MP2), SOPPA(CC2), and SOPPA(CCSD), the best result was achieved with SOPPA(CCSD) when used in combination with Sauer's basis set aug‐cc‐pVTZ‐J characterized by the mean absolute error of calculated coupling constants against the experiment of ca 2 Hz in the range of ca 200 Hz. The SOPPA(CCSD)/aug‐cc‐pVTZ‐J method is recommended as the most accurate and effective computational scheme for the calculation of ¹J(Si,C). The slightly less accurate but essentially more economical SOPPA(MP2)/aug‐cc‐pVTZ‐J and/or SOPPA(CC2)/aug‐cc‐pVTZ‐J methods are recommended for larger molecular systems. It was shown that solvent and relativistic corrections do not play a major role in the computation of the total values of ¹J(Si,C); however, taking them into account noticeably improves agreement with the experiment. The rovibrational corrections are estimated to be of about 1 Hz or 1–1.5% of the total value of ¹J(Si,C). Copyright © 2014 John Wiley & Sons, Ltd.     

Verification of DFT-predicted hydrogen storage capacity of VC3H3 complex using molecular dynamics simulations

Density functional theory (DFT) and Fourth‐order Möller–Plesset (MP4) perturbation theory calculations are performed to examine the possibility of hydrogen storage in V‐capped VC₃H₃ complex. Stability of bare and H₂ molecules adsorbed V‐capped VC₃H₃ complex is verified using DFT and MP4 method. Thermo‐chemistry calculations are carried out to estimate the Gibbs free corrected averaged H₂ adsorption energy which reveals whether H₂ adsorption on V‐capped VC₃H₃ complex is energetically favorable, at different temperatures. We use different exchange and correlation functionals employed in DFT to see their effect on H₂ adsorption energy. Molecular dynamic (MD) simulations are performed to confirm whether this complex adsorbs H₂ molecules at a finite temperature. We elucidate the correlation between H₂ adsorption energy obtained from density functional calculations and retaining number of H₂ molecules on VC₃H₃ complex during MDs simulations at various temperatures. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Core-valence correlation effects for molecules containing first-row atoms. Accurate results using effective core polarization potentials

The accuracy of employing effective core polarization potentials (CPPs) to account for the effects of core-valence correlation on the spectroscopic constants and dissociation energies of the molecules B₂, C₂, N₂, O₂, F₂, CO, CN, CH, HF, and C₂H₂ has been investigated by comparison to accurate all-electron benchmark calculations. The results obtained from the calculations employing CPPs were surprisingly accurate in every case studied, reducing the errors in the calculated valence D _e values from a maximum of nearly 2.5 kcal/mol to just 0.3 kcal/mol. The effects of enlarging the basis set and using higher-order valence electron correlation treatments were found to have only a small influence on the core-valence correlation effect predicted by the CPPs. Thus, to accurately recover the effects of intershell correlation, effective core polarization potentials such as the ones used in the present work provide an attractive alternative to carrying out computationally demanding calculations where the core electrons are explicitly included in the correlation treatment. Received: 11 May 1998 / Accepted: 27 July 1998 / Published online: 28 October 1998     

Binding energy of gas molecule with two pyrazine molecules as organic linker in metal�Corganic framework: its theoretical evaluation and understanding of determining factors

We explored the interactions of gas molecules such as H₂, CH₄, C₂H₄, C₂H₆, CO₂, and CS₂ sandwiched by two pyrazine (Pz) molecules, which were employed as a model of organic linker in the Hofmann-type metal?Corganic framework (MOF). The MP2.5/aug-cc-pVTZ method was employed here, because this method presents almost the same binding energy as that calculated by the CCSD(T)/aug-cc-pVDZ with MP2.5-evaluated basis set extension effects to aug-cc-pVTZ basis set. The binding energy of the gas molecule increases in the order H₂?<?CH₄?<?CO₂?<?C₂H₄????C₂H₆?<?CS₂. The energy decomposition analysis of the interaction energy indicates that the electrostatic term presents the largest contribution to the interaction energy at the Hartree?CFock level. However, the dispersion interaction provides dominant contribution to the total binding energy at correlated level. We newly found a linear correlation between the z-component of polarizability of gas molecules and dispersion energy, where the z-axis was taken to be perpendicular to two Pz rings. These results are useful for understanding and predicting the binding energy of the gas molecule with the organic linkers of MOF.     

Eclipsed and staggered conformations of (SIH3)2F+: An ab initio study

Ab initio calculations at 6–31G**, 6–31++G**, and MP2/6–31G** levels were performed on disilyl–fluoronium, (SiH₃)₂F⁺, with the SiH₃ group eclipsed or staggered. Optimized geometries, total energies, dipole moments, atomic charges, electronic density, and vibrational frequencies were computed. The results were compared with calculated structural parameters and vibrational frequencies of H₃SiF, H₂SiF⁺, H₂SiF^?, and H₄SiF⁺ ions. The basis-set effects were studied. Several thermochemistry parameters—ZPE, thermal energy, rotational constants, and entropies—were also calculated. © 1994 John Wiley & Sons, Inc.     

Ab initio study of the molecular structures, force fields, and vibrational spectra of alkaline metal fluoride dimers MM′F2 (M, M′=Li, Na, K)

Ab initio calculations of the equilibrium geometrical parameters, force constants, and IR vibration frequencies and intensities of Li₂F₂, Na₂F₂, K₂F₂, LiNaF₂, LiKF₂, and NaKF₂ are reported. The calculations use the Hartree-Fock-Roothaan method and second-order Möller-Plesset perturbation theory along with configuration interaction theory including singly and doubly excited configurations and corrections for quartic excitations with basis sets of grouped Gaussian functions: Li — (9s3p1d/4s3p1d), Na — (12s8p1d/6s4p1d), K — (14s11p3d/9s8p3d), F — (9s5p1d/4s2p1d). According to the results of calculations, the equilibrium structures of the molecules are planar cyclic structures of D_2h (for M₂F₂) and C_2v (for MM′F₂) symmetries. The linear configurations M-F-M′-F (of C_∞v symmetry) are 70–190 kJ/mole less stable than the cyclic ones; for all molecules except M-F-K-F, these configurations correspond to local minima on the potential energy surface. The role of correlation effects in ab initio calculations of the geometry, force fields, and IR characteristics of molecules with highly polar chemical bonds is discussed. The theoretical force fields of the molecules are represented in canonical form in a system of redundant natural vibrational coordinates. The force fields of MM′F₂ molecules are studied. The results of the ab initio calculations are compared with the experimental structural data and vibrational spectra available in the literature.