Ab Initio calculation of 19F NMR chemical shielding for alkaline‐earth‐metal fluorides

Gauge‐independent atomic orbital (GIAO) method at Hartree‐Fock (HF) and density functional theory (DFT) levels, respectively, was employed to calculate ¹⁹F NMR chemical shieldings of solid state alkaline‐earth‐metal fluorides MF₂ (M = Mg, Ca, Sr, Ba). The results show that, although the calculated ¹⁹F chemical shieldings tend to be larger than the experimental values, they have a fairly good linear relationship with the observed ones. The calculated results based on different combinations of basis sets show that the B3LYP (a hybrid of DFT with HF) predictions are greatly superior to the HF predictions. When a basis set of metal atom with effective core potential (ECP) has well representation of valence wavefunction, especially wavefunction of d component, and proper definition of core electron number, it can be applied to obtain ¹⁹F chemical shielding which is close to that of all‐electron calculation. The variation of ¹⁹F chemical shielding of alkaline‐earth‐metal fluorides correlates well with the lattice factor A/R².     

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.     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

NLO‐X (X = I–III): New Gaussian basis sets for prediction of linear and nonlinear electric properties

New adjusted Gaussian basis sets are proposed for first and second rows elements (H, B, C, N, O, F, Si, P, S, and Cl) with the purpose of calculating linear and mainly nonlinear optical (L–NLO) properties for molecules. These basis sets are new generation of Thakkar‐DZ basis sets, which were recontracted and augmented with diffuse and polarization extrabasis functions. Atomic energy and polarizability were used as reference data for fitting the basis sets, which were further applied for prediction of L–NLO properties of diatomic, H₂, N₂, F₂, Cl₂, BH, BF, BCl, HF, HCl, CO, CS, SiO, PN, and polyatomic, CH₄, SiH₄, H₂O, H₂S, NH₃, PH₃, OCS, NNO, and HCN molecules. The results are satisfactory for all electric properties tested; dipole moment (µ), polarizability (α), and first hyperpolarizability (β), with an affordable computational cost. Three new basis sets are presented and called as NLO‐I (ADZP), NLO‐II (DZP), and NLO‐III (VDZP). The NLO‐III is the best choice to predict L–NLO properties of large molecular systems, because it presents a balance between computational cost and accuracy. The average errors for β at B3LYP/NLO‐III level were of 8% for diatomic molecules and 14% for polyatomic molecules that are within the experimental uncertainty. © 2014 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     

Electron correlation in the interacting quantum atoms partition via coupled‐cluster lagrangian densities

The electronic energy partition established by the Interacting Quantum Atoms (IQA) approach is an important method of wavefunction analyses which has yielded valuable insights about different phenomena in physical chemistry. Most of the IQA applications have relied upon approximations, which do not include either dynamical correlation (DC) such as Hartree‐Fock (HF) or external DC like CASSCF theory. Recently, DC was included in the IQA method by means of HF/Coupled‐Cluster (CC) transition densities (Chávez‐Calvillo et al., Comput. Theory Chem. 2015, 1053, 90). Despite the potential utility of this approach, it has a few drawbacks, for example, it is not consistent with the calculation of CC properties different from the total electronic energy. To improve this situation, we have implemented the IQA energy partition based on CC Lagrangian one‐ and two‐electron orbital density matrices. The development presented in this article is tested and illustrated with the H₂, LiH, H₂O, H₂S, N₂, and CO molecules for which the IQA results obtained under the consideration of (i) the CC Lagrangian, (ii) HF/CC transition densities, and (iii) HF are critically analyzed and compared. Additionally, the effect of the DC in the different components of the electronic energy in the formation of the T‐shaped (H₂)₂ van der Waals cluster and the bimolecular nucleophilic substitution between F^– and CH₃F is examined. We anticipate that the approach put forward in this article will provide new understandings on subjects in physical chemistry wherein DC plays a crucial role like molecular interactions along with chemical bonding and reactivity. © 2016 Wiley Periodicals, Inc.     

H2O,H2, HF,F2 and F2O nuclear magnetic shielding constants and indirect nuclear spin–spin coupling constants (SSCCs) in the BHandH/pcJ‐n and BHandH/XZP Kohn–Sham limits

Good performance of segmented contracted basis sets XZP, where X = D, T, Q and 5, for obtaining H₂O, H₂, HF, F₂ and F₂O nuclear isotropic shielding constants in the BHandH Kohn–Sham basis set limit was shown. The results of two‐ and three‐parameter complete basis set limit extrapolation schemes were compared with experimental results, earlier literature data and benchmark ab initio results. Similar convergence patterns of shieldings obtained from calculations using general purpose XZP basis sets and from polarization‐consistent basis sets pcS‐n and pcJ‐n, where n = 0, 1, 2, 3 and 4, designed to accurately predict magnetic properties were observed. On the contrary, the SSCCs were more sensitive to the XZP basis set size and generally less accurate than those estimated using pcJ‐n basis set family. The BHandH density functional markedly outperforms B3LYP method in predicting heavy atom shieldings and SSCCs values in the studied systems. Copyright © 2009 John Wiley & Sons, Ltd.     

Enhanced second‐order treatment of electron pair correlation

Ab initio calculations have been performed for F₂, HCCH, H₂O, HF, (HF)₂, and (H₂O)₂, comparing certain electron pair correlation methods, or methods for doubly substituted configurations. In these model systems, the reweighting of substituted configurations that occurs beyond a second‐order perturbative treatment of electron correlation can be partly built into the second‐order analysis in a computationally trivial step. Specific means for doing this are explored, and they offer improvement in certain cases or else very little change. A consistent improvement in the correlation energy when judged against treatment with double substitution coupled cluster theory for the test species is obtained through one of these schemes. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 226–236, 2000     

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.     

Electronic structure and molecular orbital study of the first excited state of the high‐efficiency blue OLED material bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex from ab initio and TD‐B3LYP

Bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex (AlMq₂OH) is used in organic light‐emitting diodes (OLEDs) as an electron transport material and emitting layer. By means of ab initio Hartree–Fock (HF) and density functional theory (DFT) B3LYP methods, the structure of AlMq₂OH was optimized. The frontier molecular orbital characteristics and energy levels of AlMq₂OH have been analyzed systematically to study the electronic transition mechanism in AlMq₂OH. For comparison and calibration, bis(8‐quinolinolato)aluminum(III) hydroxide complex (Alq₂OH) has also been examined with these methods using the same basis sets. The lowest singlet excited state (S₁) of AlMq₂OH has been studied by the singles configuration interaction (CIS) method and time‐dependent DFT (TD‐DFT) using a hybrid functional, B3‐LYP, and the 6‐31G* basis set. The lowest singlet electronic transition (S₀ → S₁) of AlMq₂OH is π → π* electronic transitions and primarily localized on the different quinolate ligands. The emission of AlMq₂OH is due to the electron transitions from a phenoxide donor to a pyridyl acceptor from another quinolate ligand including C → C and O → N transference. Two possible electron transfer pathways are presented, one by carbon, oxygen, and nitrogen atoms and the other via metal cation Al³⁺. The comparison between the CIS‐optimized excited‐state structure with the HF ground‐state structure indicates that the geometric shift is mainly confined to the one quinolate and these changes can be easily understood in terms of the nodal patterns of the highest occupied and lowest unoccupied molecular orbitals. On the basis of the CIS‐optimized structure of the excited state, TD‐B3‐LYP calculations predict an emission wavelength of 499.78 nm. An absorption wavelength at 380.79 nm on the optimized structure of B3LYP/6‐31G* was predicted. They are comparable to AlMq2OH 485 and 390 nm observed experimentally for photoluminescence and UV‐vis absorption spectra of AlMq₂OH solid thin film on quartz, respectively. Lending theoretical corroboration to recent experimental observations and supposition, the reasons for the blue‐shift of AlMq2OH were revealed. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Correlated one‐electron wave functions

Using four basis sets, 6‐311G(d,p), 6‐31+G(d,p), 6‐311++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the acidic H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. By contrast with above the three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen‐bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is a noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD(T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Accurate ab initio potential energy surface and vibration‐rotation energy levels of hydrogen peroxide

The accurate ground‐state potential energy surface of hydrogen peroxide, H₂O₂, has been determined from ab initio calculations using the coupled‐cluster approach in conjunction with the correlation‐consistent basis sets up to septuple‐zeta quality. Results obtained with the conventional and explicitly correlated coupled‐cluster methods were compared. The core–electron correlation, scalar relativistic, and higher‐order valence–electron correlation effects were taken into account. The adiabatic effects were also discussed. The vibration–rotation energy levels of the H₂O₂, D₂O₂, and HOOD isotopologues were predicted, and the experimental vibrational fundamental wavenumbers were reproduced to 1 cm^?1 (“spectroscopic”) accuracy. © 2012 Wiley Periodicals, Inc.     

One‐bond 29Si‐1H spin‐spin coupling constants in the series of halosilanes: benchmark SOPPA and DFT calculations,relativistic effects,and vibrational corrections

A number of most representative second order polarization propagator approach (SOPPA) based wavefunction methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), and density functional theory (DFT) based methods, B3LYP, PBE0, KT2, and KT3, have been benchmarked in the calculation of the one‐bond ²⁹Si‐¹H spin‐spin coupling constants in the series of halosilanes SiH_nX_4?n (X = F, Cl, Br, I), both at the non‐relativistic and full four‐parameter Dirac's relativistic levels taking into account vibrational corrections. At the non‐relativistic level, the wavefunction methods showed much better results as compared with those of DFT. At the DFT level, out of four tested functionals, the Perdew, Burke, and Ernzerhof's PBE0 showed best performance. Taking into account, relativistic effects and vibrational corrections noticeably improves wavefunction methods results, but generally worsens DFT results. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Quantum‐chemical,IR, NMR,and X‐ray diffraction studies on 2‐(4‐chlorophenyl)‐1‐methyl‐ 1H‐benzo[d]imidazole

The title molecule, 2‐(4‐chlorophenyl)‐1‐methyl‐1H‐benzo[d]imidazole (C₁₄H₁₁ClN₂), was prepared and characterized by ¹H NMR, ¹³C NMR, IR, and single‐crystal X‐ray diffraction. The molecular geometry, vibrational frequencies, and gauge including atomic orbital (GIAO) ¹H and ¹³C NMR chemical shift values of the title compound in the ground state have been calculated by using the Hartree‐Fock (HF) and density functional theory (DFT/B3LYP) method with 6‐31G(d) basis sets, and compared with the experimental data. The calculated results show that the optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies and GIAO ¹H and ¹³C NMR chemical shifts show good agreement with experimental values. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6‐31G(d) basis set by applying the Onsager and the polarizable continuum model (PCM). Besides, molecular electrostatic potential (MEP), frontier molecular orbitals (FMO) analysis, and nonlinear optical (NLO) properties of the title compound were investigated by theoretical calculations. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Synthesis of polyarylacetylenes by γ‐ray‐induced polymerization of terminal alkynes. Nanostructures of ortho‐substituted derivatives

Terminal aryl alkynes RC₆H₄C?CH with substituents of different electronic properties and ring position (R = H, 4‐CF₃, 4‐OMe, 2‐CF₃, 2‐OMe, 2‐Me) were exposed to γ‐radiation (50–400 kGy) in organic solvents (hexane, 1,4‐dioxane, ethylacetate, methanol, tetrahydrofuran), at room temperature. The effects arising from substituent, solvent, dilution, and radiation dose allowed to define the conditions suitable for polymerization, which was favored in methanol at increasing dilution of the alkyne. Ortho‐substitution represented the key structural element in the substrate, and the derived polyarylacetylenes were characterized in detail, including gel permeation chromatography, thermal analysis, infrared, NMR, UV–vis, fluorescence, and scanning electron microscope spectroscopy. The results are consistent with the formation of irregular polymers mainly composed of trans‐transoid chains. Controlled aggregation of the polymers by means of an osmosis‐based procedure in solvent/non‐solvent mixtures allowed the formation of nanostructured materials, in particular of hollow nanospheres from THF/water. The methodology sets the basis for the development of γ‐rays‐induced polymerization of alkynes, in a transition metal catalyst‐free environment. © 2012 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem, 2012     

Theoretical study of the electronic spectra of bi‐ and tri‐heteronuclear platinum complexes

The electronic structure and the spectroscopic properties of [Pt(NH₃)₄][Au(CN)₂]₂, [Pt(NH₃)₄][Ag(CN)₂]₂, [Pt(CNCH₃)₄][Pt(CN)₄], and [Pt(CNCH₃)₄][Pd(CN)₄] were studied at the HF, MP2, B3LYP, and PBE levels. In all the complexes, it was found that the nature of the intermetal interactions is consistent with the presence of a high‐ionic contribution (90%) and a dispersion‐type interaction (10%). The absorption spectra of these complexes were calculated by the single‐excitation time‐dependent (TD) method at the HF, B3LYP, and PBE levels. The [Pt(NH₃)₄][M(CN)₂]₂ (M ? Au, Ag) complexes showed a ¹(dσ* → pσ) transition associated with a metal–metal charge transfer. On the other hand, the [Pt(CNCH₃)₄][M(CN)₄] (M ? Pt, Pd) complexes showed a ¹(dσ* → π*) transition associated with a metal‐to‐metal and ligand charge transfer. The values obtained theoretically are in agreement with the experimental range. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Ab initio potential energy surface and vibration‐rotation energy levels of sulfur dioxide

An accurate potential energy surface of sulfur dioxide, SO₂, in its ground electronic state has been determined from ab initio calculations using the coupled‐cluster approach in conjunction with the correlation‐consistent basis sets up to septuple‐zeta quality. The results obtained with the conventional and explicitly correlated coupled‐cluster methods are compared. The role of the core–electron correlation, higher‐order valence–electron correlation, scalar relativistic, and adiabatic effects in determining the structure and dynamics of the SO₂ molecule is discussed. The vibration‐rotation energy levels of the ³²SO₂ and ³⁴SO₂ isotopologues were predicted using a variational approach. It was shown that the inclusion of the aforementioned effects was mandatory to attain the “spectroscopic” accuracy. © 2017 Wiley Periodicals, Inc.     

Calculations of solid‐state 43Ca NMR parameters: A comparison of periodic and cluster approaches and an evaluation of DFT functionals

We present a computational study of magnetic‐shielding and quadrupolar‐coupling tensors of ⁴³Ca sites in crystalline solids. A comparison between periodic and cluster‐based approaches for modeling solid‐state interactions demonstrates that cluster‐based approaches are suitable for predicting ⁴³Ca NMR parameters. Several model chemistries, including Hartree–Fock theory and 17 DFT approximations (SVWN, CA‐PZ, PBE, PBE0, PW91, B3PW91, rPBE, PBEsol, WC, PKZB, BMK, M06‐L, M06, M06‐2X, M06‐HF, TPSS, and TPSSh), are evaluated for the prediction of ⁴³Ca NMR parameters. Convergence of NMR parameters with respect to basis sets of the form cc‐pVXZ (X = D, T, Q) is also evaluated. All DFT methods lead to substantial, and frequently systematic, overestimations of experimental chemical shifts. Hartree–Fock calculations outperform all DFT methods for the prediction of ⁴³Ca chemical‐shift tensors. © 2017 Wiley Periodicals, Inc.