Structure and vibrational spectrum of the hydrogen-bonded system between 4-tert-butylphenol and N-bases: Ab initio and DFT studies

The structural and vibrational characteristics of the hydrogen-bonded system between 1,5,7-triazabicyclo [4.4.0] dec-5-ene (TBD) and 4-tert-butylphenol have been investigated employing ab initio and DFT calculations at different basis sets. The calculations show that the optimized structure of the studied system is cyclic. The corrected values of the dissociation energy for the hydrogen-bonded complex have been calculated in order to estimate its stability. The influence of the hydrogen bonding on the properties of the monomers (TBD and 4-tert-butylphenol) has been investigated. The hydrogen bonding between TBD and 4-tert-butylphenol leads to changes in the structural (bond lengths and angles) and vibrational (vibrational frequencies and infrared intensities) characteristics of the monomers. It was established that the TBD molecule is considerably deformed upon hydrogen bonding, while the deformation of the 4-t-BuPhOH is smaller. In agreement with the experiment, the calculations show that the stretching O-H vibration from 4-tert-butylphenol is shifted to lower frequency upon hydrogen bonding. The predicted frequency shift Deltanu(O-H) (-338cm(-1)) is in very good agreement with the experimentally observed (-351cm(-1)). In the same time the IR intensity of the nu(O-H) increases dramatically in the hydrogen-bonded system.     

Cluster and periodic DFT calculations of MgO/Pd(CO) and MgO/Pd(CO)(2) surface complexes

The bonding and vibrational properties of Pd(CO) and Pd(CO)(2) complexes formed at the (100) surface of MgO have been investigated using the gradient-corrected DFT approach and have been compared to the results of infrared and thermal desorption experiments performed on ultrathin MgO films. Two complementary approaches have been used for the calculation of the electronic properties: the embedded cluster method using localized atomic orbital basis sets and supercell periodic calculations using plane waves. The results show that the two methods provide very similar answers, provided that sufficiently large supercells are used. Various regular and defect adsorption sites for the Pd(CO) and Pd(CO)(2) have been considered: terraces, steps, neutral and charged oxygen vacancies (F and F(+) centers), and divacancies. From the comparison of the computed and experimental results, it is concluded that the most likely site where the Pd atoms are stabilized and where carbonyl complexes are formed are the F(+) centers, paramagnetic defects consisting of a single electron trapped in an anion vacancy.     

Vibrational spectroscopic (FT-IR and FT-Raman), first-order hyperpolarizablity, HOMO, LUMO, NBO, Mulliken charges and structure determination of 2-bromo-4-chlorotoluene

The FT-IR and FT-Raman spectra of 2-bromo-4-chlorotoluene (2B4CT) molecule have been recorded in the region 4000-400 cm(-1) and 3500-50 cm(-1), respectively. Optimized geometrical structures, harmonic vibrational frequencies, intensities, reduced mass, force constants and depolarization ratio have been computed by the B3 based (B3LYP) density functional methods using 6-31+G(d,p) and 6-311++G(d,p) basis sets. The observed FT-IR and FT-Raman vibrational frequencies are analysed and compared with theoretically predicted vibrational frequencies. The geometries and normal modes of vibration obtained from DFT method are in good agreement with the experimental data. The Mulliken charges, the natural bonding orbital (NBO) analysis, the values of electric dipole moment (μ) and the first-order hyperpolarizability (β) of the investigated molecule were computed using DFT calculations. The calculated HOMO and LUMO energies show that charge transfer occurs within molecule. The influences of bromine atom, chlorine atom and methyl group on the geometry of benzene and its normal modes of vibrations have also been discussed.     

Ab initio and DFT studies of the vibrational spectra of hydrogen-bonded PhOH...(H2O)4 complexes

The vibrational characteristics (vibrational frequencies and infrared intensities) for the hydrogen-bonded complex of phenol with four water molecules PhOH...(H2O)4 (structure 4A) have been predicted using ab initio and DFT (B3LYP) calculations with 6-31G(d,p) basis set. The changes in the vibrational characteristics from free monomers to a complex have been calculated. The ab initio and B3LYP calculations show that the observed four intense bands at 3299, 3341, 3386 and 3430 cm(-1) can be assigned to the hydrogen-bonded OH stretching vibrations in the complex PhOH...(H2O)4 (4A). The complexation leads to very large red shifts of these vibrations and very strong increase in their IR intensity. The predicted red shifts for these vibrations with B3LYP/6-31G(d,p) calculations are in very good agreement with the experimentally observed. It was established that the phenolic OH stretching vibration is the most sensitive to the hydrogen bonding. The predicted red-shift with the B3LYP/6-31G(d,p) calculations for the most stable ring structure 4A (-590 cm(-1)) is in better agreement with the experimentally observed than the red-shift, predicted with SCF/6-31G(d,p) calculations. The magnitude of the wavenumber shift is indicative of relatively strong OH...H hydrogen-bonded interaction. The complexation between phenol and four water molecules leads to strong increase of the IR intensity of the phenolic OH stretching vibration (up to 38 times).     

Molecular structures,bond energies,and bonding analysis of group 11 cyanides TM(CN) and isocyanides TM(NC) (TM = Cu,Ag, Au)

We report on quantum chemical calculations at the DFT (BP86/TZP) and ab initio (CCSD(T)/III+) levels of the title compounds. The geometries, vibrational spectra, heats of formation, and homolytic and heterolytic bond dissociation energies are given. The calculated bond length of Cu-CN is in reasonable agreement with experiment. The theoretical geometries for CuNC and the other group 11 cyanides and isocyanides which have not been measured as isolated species provide a good estimate for the exact values. The theoretical bond dissociation energies and heats of formation should be accurate with an error limit of +/-5 kcal/mol. The calculation of the vibrational spectra shows that the C-N stretching mode of the cyanides, which lies between 2170 and 2180 cm(-)(1), is IR inactive. The omega(1)(C-N) vibrations of the isocyanides are shifted by approximately 100 cm(-)(1) to lower wavenumbers. They are predicted to have a very large IR intensity. The nature of the metal-ligand interactions was investigated with the help of an energy partitioning analysis in two different ways using the charged fragments TM(+) + CN(-) (TM = transition metal) and the neutral fragments TM(*) + CN(*) as bonding partners. The calculations suggest that covalent interactions are the driving force for the formation of the TM-CN and TM-NC bonds, but the finally formed bonds are better described in terms of interactions between TM(+) and CN(-), which have between 73% and 80% electrostatic character. The contribution of the pi bonding is rather small. The lower energy of the metal cyanides than that of the isocyanides comes from the stronger electrostatic interaction between the more diffuse electron density at the carbon atom of the cyano ligand and the positively charged nucleus of the metal.     

FTIR and FT-Raman spectra, assignments, ab initio HF and DFT analysis of xanthine

The molecular vibrations of xanthine were investigated in polycrystalline sample, at room temperature by Fourier transform infrared (FTIR) and FT-Raman spectroscopies. The spectra of the molecule have been recorded in the regions 4000-50 cm(-1) and 3500-100 cm(-1), respectively. Theoretical information on the optimized geometry, harmonic vibrational frequencies, infrared and Raman intensities were obtained by means of ab initio Hartree-Fock (HF) and density functional theory (DFT) gradient calculations with complete relaxation in the potential energy surface using 6-311++G(d,p) basis set. The vibrational frequencies which were determined experimentally from the spectral data are compared with those obtained theoretically from ab initio and DFT calculations. A close agreement was achieved between the observed and calculated frequencies by refinement of the scale factors. The infrared and Raman spectra were also predicted from the calculated intensities. Thermodynamic properties like entropy, heat capacity, zero point energy have been calculated for the molecule. Unambiguous vibrational assignment of all the fundamentals was made using the potential energy distribution (PED).     

Molecular orbital cluster model study of bonding and vibrations of CO adsorbed on MgO surface

The interaction of CO with the MgO(100) surface has been investigated by means of all electron cluster model calculations. The CO molecule is bound on the Mg²⁺ site of MgO with a chemisorption energy of about 0.2 eV. The binding mechanism is electrostatic in nature and arises almost entirely from the interaction of the weak electric field generated by the ionic surface and the CO charge distribution, with negligible contributions from chemical effects as the CO σ donation. When CO is bound through carbon, its vibrational frequency increases with respect to the gas-phase value. This shift, Δ, has been analyzed and decomposed into the sum of different contributions. It is found that the positive Δω does not arise entirely from the field–dipole interaction but is due, in part, to the increase in Pauli repulsion occurring when the C? O molecule vibrates in the presence of the surface “wall.” A stronger electrostatic interaction, bringing the CO adsorbate closer to the surface, increases this wall effect and results in a more pronounced positive ω shift. It is also found that the two CO orientations exhibit opposite shifts in ω_e, thus, the two orientations can be distinguished, in principle, by IR spectroscopy. The analysis of our ab initio cluster wave functions gives a very different picture than the standard view of the metal–CO bond as arising from σ donation and π back donation.     

Theoretical study of the changes in the vibrational characteristics arising from the hydrogen bonding between Vitamin C (L-ascorbic acid) and H2O

The vibrational characteristics (vibrational frequencies, infrared intensities and Raman activities) for the hydrogen-bonded system of Vitamin C (L-ascorbic acid) with five water molecules have been predicted using ab initio SCF/6-31G(d,p) calculations and DFT (BLYP) calculations with 6-31G(d,p) and 6-31++G(d,p) basis sets. The changes in the vibrational characteristics from free monomers to a complex have been calculated. The ab initio and BLYP calculations show that the complexation between Vitamin C and five water molecules leads to large red shifts of the stretching vibrations for the monomer bonds involved in the hydrogen bonding and very strong increase in their IR intensity. The predicted frequency shifts for the stretching vibrations from Vitamin C taking part in the hydrogen bonding are up to -508 cm(-1). The magnitude of the wavenumber shifts is indicative of relatively strong OH...H hydrogen-bonded interactions. In the same time the IR intensity and Raman activity of these vibrations increase upon complexation. The IR intensity increases dramatically (up to 12 times) and Raman activity increases up to three times. The ab initio and BLYP calculations show, that the symmetric OH vibrations of water molecules are more sensitive to the complexation. The hydrogen bonding leads to very large red shifts of these vibrations and very strong increase in their IR intensity. The asymmetric OH stretching vibrations of water, free from hydrogen bonding are less sensitive to the complexation than the hydrogen-bonded symmetric OH stretching vibrations. The increases of the IR intensities for these vibrations are lower and red shifts are negligible.     

Density functional theory study on molecular structure and vibrational IR spectra of isocytosine

Density functional theory with the combined Becke3-LYP exchange-correlation energy functional [DFT(B3-LYP) method] using the 6-31G(d, p) basis set is applied to predict molecular parameters (geometries, rotational constants, dipole moments) and vibrational IR spectra (harmonic wavenumbers, absolute intensities) of six tautomers of the isocytosine molecule. The results are compared with the corresponding data calculated at the conventional ab initio Hartree-Fock (HF) level using the same basis set and with available experimental data. Calculations show that (a) three amino tautomers are slightly nonplanar species with, evidently, a distorted amino group, (b) the DFT (B3-LYP)/6-31G(d, p) method predicts better molecular parameters, than do the HF calculations, and (c) the DFT(B3-LYP)-calculated vibrational IR spectra of isocytosine agree well with the available recorded IR spectra, and they show marked improvement over the IR spectra predicted at the HF/6-31G(d, p) level. Tautomeric stabilities of isocytosine are discussed on the basis of computed electronic energies by the DFT(B3-LYP) and ab initio approaches [including the MP2 and MP4(SDQ) calculations of electronic energies] and predicted zero-point vibrational energies by DFT(B3-LYP) and HF methods. This relative energies at 0 K of the tautomeric forms of isocytosine predicted by both conventional ab initio and DFT(B3-LYP) methods correlate well with the experimental data, showing the predominance of the aminohydroxy tautomer of isocytosine for an isolated molecule. © 1997 John Wiley & Sons, Inc.     

Charging of Au atoms on TiO2 thin films from CO vibrational spectroscopy and DFT calculations

Au atoms have been deposited on oxidized and reduced TiO2 thin films grown on Mo(110). The gold binding sites and the occurrence of Au-TiO2 charge transfer were identified by measuring infrared spectra as a function of temperature and substrate preparation. The results have been interpreted by slab model DFT calculations. Au binds weakly to regular TiO2 sites (De < 0.5 eV) where it remains neutral, and diffuses easily even at low temperature until it gets trapped at strong binding sites such as oxygen vacancies (De = 1.7 eV). Here, a charge transfer from TiO2 to Au occurs. Au(delta-)CO complexes formed on oxygen vacancies easily lose CO (De = 0.4 eV), and the CO stretching frequency is red-shifted. On nondefective surfaces, CO adsorption induces a charge transfer from Au to TiO2 with formation of strongly bound Audelta+CO complexes (De = 2.4 eV); the corresponding CO frequency is blue-shifted with respect to free CO. We propose possible mechanisms to reconcile the observed CO desorption around 380 K with the unusually high stability of Au-CO complexes formed on regular sites predicted by the calculations. This implies: (a) diffusion of AuCO complexes above 150 K; (b) formation of gold dimers when the diffusing AuCO complex encounters a Au atom bound to an oxygen vacancy (reduced TiO2) or a second AuCO unit (oxidized TiO2); and (c) CO desorption from the resulting dimer, occurring around 350-400 K.     

Vibrational analyses for CHFClBr and CDFClBr based on high level ab initio calculations

Anharmonicity corrections to the harmonic vibrational spectra of CHFClBr and its deuterated isotopomer were computed by means of variational and perturbational approaches. A comparison of both methods is provided. Based on CCSD(T)/aug-cc-pVTZ electronic structure calculations excellent agreement with experimental data was obtained. Absolute mean deviations are in the range of about 4 cm(-1) for the fundamental modes, while slightly larger values of about 7 cm(-1) were found for the first vibrational overtones. In addition, vibrationally averaged structural parameters are provided for both molecules. The calculations will serve as a future starting point for parity-violation effects in vibrational transitions in these chiral molecules.     

Quantum cascade laser spectroscopy and photoinduced chemistry of Al-(CO)n clusters in helium nanodroplets

Helium nanodroplet isolation and a tunable quantum cascade laser are used to probe the fundamental CO stretch bands of aluminum carbonyl complexes, Al-(CO)(n) (n ≤ 5). The droplets are doped with single aluminum atoms via the resistive heating of an aluminum wetted tantalum wire. The downstream sequential pick-up of CO molecules leads to the rapid formation and cooling of Al-(CO)(n) clusters within the droplets. Near 1900 cm(-1), rotational fine structure is resolved in bands that are assigned to the CO stretch of a linear (2)Π(1/2) Al-CO species and the asymmetric and symmetric CO stretch vibrations of a planar C(2v) Al-(CO)(2) complex in a (2)B(1) electronic state. Bands corresponding to clusters with n ≥ 3 lack resolved rotational fine structure; nevertheless, the small frequency shifts from the n = 2 bands indicate that these clusters consist of an Al-(CO)(2) core with additional CO molecules attached via van der Waals interactions. A second n = 2 band is observed near the CO stretch of Al-CO, indicating a local minimum on the n = 2 potential consisting of an "unreacted" (Al-CO)-CO cluster. The line width of this band is ～0.3 cm(-1), which is about 30 times broader than the transitions within the Al-CO band. The additional broadening is consistent with a homogeneous mechanism corresponding to a rapid vibrational excitation induced reaction within the (Al-CO)-CO cluster to form the covalently bonded Al-(CO)(2) complex. Ab initio CCSD(T) calculations and natural bond orbital (NBO) analyses are carried out to investigate the nature of the bonding in the n = 1, 2 complexes. The NBO calculations show that both π-donation (from the occupied aluminum p orbital into a π* antibonding CO orbital) and σ-donation (from CO into the empty aluminum p orbitals) play a significant role in the bonding, analogous to transition-metal carbonyl complexes. The large red shift observed for the CO stretch vibrations is consistent with this bonding analysis.     

Experimental and theoretical studies of the vibrational spectra of cis-1-bromo-2-fluoroethene

The gas-phase infrared spectrum of cis-1-bromo-2-fluoroethene has been studied at low resolution in the range 200-6500 cm(-1), leading to a complete assignment of the fundamentals, except the lowest vibrational mode nu9 predicted at 167 cm(-1). The remaining vibrational structure has been mainly interpreted in terms of first overtone or two quanta combination bands. Isotopic (79/91)Br shift has been observed only in the nu8 fundamental. The equilibrium structure and the quadratic force field have been investigated theoretically at CCSD(T) level of theory employing Dunning's correlation consistent triple-zeta basis set. Cubic and semidiagonal quartic force field have been calculated using second-order M?ller-Plesset perturbation theory and Ahlrich' split valence (SV) contracted basis set. After a minor scaled quantum mechanical (SQM) adjustment of the quadratic force constants, the vibrational analysis, based on the second-order perturbation theory, has been carried out with the calculated force constants.     

Gas-phase IR spectroscopy of anionic iron carbonyl clusters

The first gas-phase vibrational spectra are presented for several anionic iron carbonyl clusters, ranging in size from Fe(CO)4- to Fe5(CO)14- in the CO-stretching region (1600-2100 cm-1). The experimental spectra provide some immediate structural information about the clusters in the form of low-wavenumber (1750-1850 cm-1) bands marking the presence of bridging carbonyl ligands (mu2-COs). Supporting DFT calculations are presented for the smaller clusters (<3 Fe atoms) and give good agreement with the experimental data, allowing structural assignments for these cases. The Fe2(CO)7- spectrum suggests a structure lacking bridging carbonyl ligands, in agreement with the DFT results. For the case of Fe2(CO)8-, there are two possible structures based on the calculations, both with and without bridging carbonyls. The presence of a low-frequency band ( approximately 1770 cm-1) in the experimental spectrum conclusively demonstrates the existence of the bridged form. The ramifications of these data for metal-metal bonding in the clusters are also considered.     

A combined theoretical and experimental study of the reaction products of laser-ablated thorium atoms with CO: first identification of the CThO, CThO(-), OthCCO, OTh(eta(3)-CCO), and Th(CO)(n) (n = 1-6) molecules

Laser-ablated thorium atoms have been reacted with CO molecules during condensation in excess neon. Absorptions at 617.7 and 812.2 cm(-1) are assigned to Th-C and Th-O stretching vibrations of the CThO molecule. Absorptions at 2048.6, 1353.6, and 822.5 cm(-1) are assigned to the OThCCO molecule, which is formed by CO addition to CThO and photochemical rearrangement of Th(CO)(2). The OThCCO molecule undergoes further photoinduced rearrangement to OTh(eta(3)-CCO), which is characterized by C-C, C-O, and Th-O stretching vibrations at 1810.8, 1139.2, and 831.6 cm(-1). The Th(CO)(n) (n = 1-6) complexes are formed on deposition or on annealing. Evidence is also presented for the CThO(-) and Th(CO)(2)(-) anions, which are formed by electron capture of neutral molecules. Relativistic density functional theory (DFT) calculations of the geometry structures, vibrational frequencies, and infrared intensities strongly support the experimental assignments. It is found that CThO is an unprecedented actinide-containing carbene molecule with a triplet ground state and an unusual bent structure ( angleCThO = 109 degrees ). The OThCCO molecule has a bent structure while its rearranged product OTh(eta(3)-CCO) is found to have a unique exocyclic structure with side-bonded CCO group. We also find that both Th(CO)(2) and Th(CO)(2)(-) are, surprisingly, highly bent, with the angleC-Th-C bond angle being close to 50 degrees; the unusual geometries are the result of extremely strong Th-to-CO back-bonding, which causes significant three-centered bonding among the Th atom and the two C atoms.     

Gas-phase infrared and ab initio study of the unstable CF3CNO molecule and its stable furoxan ring dimer

The unstable trifluoroacetonitrile N-oxide molecule, CF3CNO, has been generated in high yield in the gas phase from CF3BrC=NOH and studied for the first time by gas-phase mid-infrared spectroscopy. Cold trapping of this molecule followed by slow warming forms the stable ring dimer, bis(trifluoromethyl)furoxan, also investigated by gas-phase infrared spectroscopy. The spectroscopy provides an investigation into the vibrational character of the two molecules, the assignments supported by calculations of the harmonic vibrational frequencies using in the case of CF3CNO both ab initio (CCSD(T)) and density functional theory (B3LYP) and B3LYP for the ring dimer. The ground-state structures of both molecules were investigated at the B3LYP level of theory, with CF3CNO further investigated using coupled-cluster. The CCSD(T) method suggests a slightly bent (C(s)) structure for CF3CNO, while the B3LYP method (with basis sets ranging from 6-311G(d) to cc-pVTZ) suggests a close-to-linear or linear CCNO chain. The CCN bending potential in CF3CNO was explored at the CCSD(T)(fc)/cc-pVTZ level, with the results suggesting that CF3CNO exhibits strong quasi-symmetric top behavior with a barrier to linearity of 174 cm(-1). Since both isomerization and dimerization are feasible loss processes for this unstable molecule, the relative stability of CF3CNO with respect to the known cyanate (CF3OCN), isocyanate (CF3NCO), and fulminate (CF3ONC) isomers and the mechanism of the dimerization process to the ring furoxan and other isomers were studied with density functional theory.     

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.     

Theoretical studies of UO2(H2O)n(2+), NpO2(H2O)n(+), and PuO2(H2O)n(2+) complexes (n=4-6) in aqueous solution and gas phase

Extensive ab initio calculations both in gas phase and solution have been carried out to study the equilibrium structure, vibrational frequencies, and bonding characteristics of various actinyl (UO2(2+), NpO2(+), and PuO2(2+)) and their hydrated forms, AnO2(H2O)n(z+) (n=4, 5, and 6). Bulk solvent effects were studied using a continuum method. The geometries were fully optimized at the coupled-cluster singles + doubles (CCSD), density-functional theory (DFT), and M?ller-Plesset (MP2) level of theories. In addition vibrational frequencies have been obtained at the CCSD as well as MP2/DFT levels. The results show that both the short-range and long-range solvent effects are important. The combined discrete-continuum model, in which the ionic solute and the solvent molecules in the first and second solvation shells are treated quantum mechanically while the solvent is simulated by a continuum model, can predict accurately the bonding characteristics. Moreover, our values of solvation free energies suggest that five- and six-coordinations are equally preferred for UO2(2+), and five-coordinated species are preferred for NpO2(+) and PuO2(2+). On the basis of combined quantum-chemical and continuum treatments of the hydrated complexes, we are able to determine the optimal cavity radii for the solvation models. The coupled-cluster computations with large basis sets were employed for the vibrational spectra and equilibrium geometries both of which compare quite favorably with experiment. Our most accurate computations reveal that both five- and six-coordination complexes are important for these species.     

Has AuF7 been made?

Quantum chemical calculations at DFT (BP86, B3LYP, BHLYP), MP2, CCSD, and CCSD(T) levels have been carried out on various fluoro complexes of gold in oxidation states +V through +VII to evaluate the previously claimed existence of AuF7. The calculations indicate clearly that elimination of F2 from AuF7 is a strongly exothermic reaction with a low activation barrier. This is inconsistent with the reported stability of AuF7 up to room temperature. A reported experimental vibrational frequency at 734 cm(-1) for AuF7 could not be verified computationally. It is concluded that the reported observation of AuF7 was probably erroneous. As the calculations indicate also an extremely large electron affinity and little stability for AuF6, Au(V) remains the highest well-established gold oxidation state.     

Theoretical predictions of nuclear magnetic resonance parameters in a novel organo-xenon species: chemical shifts and nuclear quadrupole couplings in HXeCCH

We calibrate the methodology for the calculation of nuclear magnetic resonance (NMR) properties in novel organo-xenon compounds. The available state-of-the-art quantum-chemical approaches are combined and applied to the HXeCCH molecule as the model system. The studied properties are (129)Xe, (1)H, and (13)C chemical shifts and shielding anisotropies, as well as (131)Xe and (2)H nuclear quadrupole coupling constants. The aim is to obtain, as accurately as currently possible, converged results with respect to the basis set, electron correlation, and relativistic effects, including the coupling of relativity and correlation. This is done, on one hand, by nonrelativistic correlated ab initio calculations up to the CCSD(T) level and, on the other hand, for chemical shifts and shielding anisotropies by the leading-order relativistic Breit-Pauli perturbation theory (BPPT) with correlated ab initio and density-functional theory (DFT) reference states. BPPT at the uncorrelated Hartree-Fock level as well as the corresponding fully relativistic Dirac-Hartree-Fock method are found to be inapplicable due to a dramatic overestimation of relativistic effects, implying the influence of triplet instability in this multiply bonded system. In contrast, the fully relativistic second-order Moller-Plesset perturbation theory method can be applied for the quadrupole coupling, which is a ground-state electric property. The performance of DFT with various exchange-correlation functionals is found to be inadequate for the nonrelativistic shifts and shielding anisotropies as compared to the CCSD(T) results. The relativistic BPPT corrections to these quantities can, however, be reasonably predicted by DFT, due to the improved triplet excitation spectrum as compared to the Hartree-Fock method, as well as error cancellation within the five main BPPT contributions. We establish three computationally feasible models with characteristic error margins for future calculations of larger organo-xenon compounds to guide forthcoming experimental NMR efforts. The predicted (129)Xe chemical shift in HXeCCH is in a novel range for this nucleus, between weakly bonded or solvated atomic xenon and xenon in the hitherto characterized molecules.