Benchmark calculations of 29Si–1H spin–spin coupling constants across double bond

Benchmark calculations of geminal and vicinal ²⁹Si–¹H spin–spin coupling constants across double bond in three reference alkenylsilanes have been carried out at both DFT and SOPPA levels in comparison with experiment. At the former, four density functionals, B3LYP, B3PW91, PBE0 and KT3, were tested in combination with five representative basis sets. At the latter, three main SOPPA‐based methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), were examined in combination with the same series of basis sets. On the whole, the wavefunction methods showed much better results as compared to DFT, with the most efficient combination of SOPPA/cc‐pVTZ‐su2 characterized by a mean absolute error of only 0.4 Hz calculated for a set of nine coupling constants in three compounds with a sample span of around 40 Hz. Copyright © 2012 John Wiley & Sons, Ltd.     

Prediction of water's isotropic nuclear shieldings and indirect nuclear spin–spin coupling constants (SSCCs) using correlation‐consistent and polarization‐consistent basis sets in the Kohn–Sham basis set limit

Density functional theory (DFT) was used to estimate water's isotropic nuclear shieldings and indirect nuclear spin–spin coupling constants (SSCCs) in the Kohn–Sham (KS) complete basis set (CBS) limit. Correlation‐consistent cc‐pVxZ and cc‐pCVxZ (x = D, T, Q, 5, and 6), and their modified versions (ccJ‐pVxZ, unc‐ccJ‐pVxZ, and aug‐cc‐pVTZ‐J) and polarization‐consistent pc‐n and pcJ‐n (n = 0, 1, 2, 3, and 4) basis sets were used, and the results fitted with a simple mathematical formula. The performance of over 20 studied density functionals was assessed from comparison with the experiment. The agreement between the CBS DFT‐predicted isotropic shieldings, spin–spin values, and the experimental values was good and similar for the modified correlation‐consistent and polarization‐consistent basis sets. The BHandH method predicted the most accurate ¹H, ¹⁷O isotropic shieldings and ¹J(OH) coupling constant (deviations from experiment of about ? 0.2 and ? 1 ppm and 0.6 Hz, respectively). The performance of BHandH for predicting water isotropic shieldings and ¹J(OH) is similar to the more advanced methods, second‐order polarization propagator approximation (SOPPA) and SOPPA(CCSD), in the basis set limit. Copyright © 2008 John Wiley & Sons, Ltd.     

Non‐empirical calculations of NMR indirect carbon‐carbon coupling constants. Part 5: Bridged bicycloalkanes

All possible J(C,C) of the bicarbocyclic frameworks together with J(C,H) and J(H,H) at bridgeheads in the series of six bridged bicycloalkanes, bicyclo[1.1.0]butane, bicyclo[2.1.0]pentane, bicyclo[3.1.0]hexane, bicyclo[2.2.0]hexane, bicyclo[3.2.0]heptane and bicyclo[3.3.0]octane, were calculated at the SOPPA level with correlation consistent Dunning sets cc‐pVTZ‐Cs augmented with inner core s‐functions and locally dense Sauer sets aug‐cc‐pVTZ‐J augmented with tight s‐functions and rationalized in terms of the multipath coupling mechanism and hybridization effects explaining many interesting structural trends. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Relativistic effects in the one‐bond spin–spin coupling constants involving selenium

One‐bond spin–spin coupling constants involving selenium of seven different types, ¹ J(Se,X), X = ¹H, ¹³C, ¹⁵ N, ¹⁹ F, ²⁹Si, ³¹P, and ⁷⁷Se, were calculated in the series of 14 representative compounds at the SOPPA(CCSD) level taking into account relativistic corrections evaluated both at the RPA and DFT levels of theory in comparison with experiment. Relativistic corrections were found to play a major role in the calculation of ¹ J(Se,X) reaching as much as almost 170% of the total value of ¹ J(Se,Se) and up to 60–70% for the rest of ¹ J(Se,X). Scalar relativistic effects (Darwin and mass‐velocity corrections) by far dominate over spin–orbit coupling in the total relativistic effects for all ¹ J(Se,X). Taking into account relativistic corrections at both random phase approximation and density functional theory levels essentially improves the agreement of theoretical results with experiment. The most ‘relativistic’ ¹ J(Se,Se) demonstrates a marked Karplus‐type dihedral angle dependence with respect to the mutual orientation of the selenium lone pairs providing a powerful tool for stereochemical analysis of selenoorganic compounds. Copyright © 2014 John Wiley & Sons, Ltd.     

Improving on equilibrium isotropic nuclear shielding constants

Isotropic nuclear shielding constants at the equilibrium molecular structure σ_eq and zero‐point vibrational corrections (ZPVCs) to σ_eq are evaluated using the B3LYP/aug‐cc‐pVTZ level of theory, as well as the KT2/aug‐cc‐pVTZ level of theory. Various scaling factors and systematic corrections are obtained by linear regression to experimental shielding constants. Comparisons of the scaled and systematically corrected equilibrium and vibrationally averaged shielding constants reveal that, at the 99% confidence level, the ZPVCs via second‐order perturbation theory do not improve the agreement of B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ calculated shielding constants with experiment. This holds true when the same analysis is applied to CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq of Teale et al. [Journal of Chemical Physics 2013, 138, 024111]. In addition, at the 99% confidence level, B3LYP/aug‐cc‐pVTZ and KT2/aug‐cc‐pVTZ scaled and systematically corrected shielding constants are found to be statistically no different from CCSD(T)/aug‐cc‐pCV[TQ]Z calculated σ_eq. The use of scaling factors and systematic corrections could thus provide a cheaper but yet reasonably accurate alternative for the study of nuclear shielding constants of larger systems.     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Ab initio calculations of the Ar–ethane intermolecular potential energy surface using bond function basis sets

The intermolecular potential energy surface (PES) of argon with ethane has been studied by ab initio calculations at the levels of second‐order Møller–Plesset perturbation (MP2) theory and coupled‐cluster theory with single, double, and noniterative triple configurations (CCSD(T)) using a series of augmented correlation‐consistent basis sets. Two sets of bond functions, bf1 (3s3p2d) and bf2 (6s6p4d2f), have been added to the basis sets to show a dramatic and systematic improvement in the convergence of the entire PES. The PES of Ar–ethane is characterized by a global minimum at a near T‐shaped configuration with a well depth of 0.611 kcal mol^?1, a second minimum at a collinear configuration with a well depth of 0.456 kcal mol^?1, and a saddle point connecting the two minima. It is shown that an augmented correlation‐consistent basis set with a set of bond functions, either bf1 or bf2, can effectively produce results equivalent to the next larger augmented correlation‐consistent basis set, that is, aug‐cc‐pVDZ‐bf1 ≈ aug‐cc‐pVTZ, aug‐cc‐pVTZ‐bf1 ≈ aug‐cc‐pVQZ. Very importantly, the use of bond functions improves the PES globally, resulting accurate potential anisotropy. Finally, MP2 method is inadequate for accurate calculations, because it gives a potentially overestimated well depth and, more seriously, a poor potential anisotropy. © 2012 Wiley Periodicals, Inc.     

Structural and vibrational analysis of the OH torsional motion in difluorohydroxyborane

In this work, we present a complete structural and vibrational analysis of the OH torsional motion in difluorohydroxyborane (BF₂OH) at the HF/aug‐cc‐pVTZ, MP2(full)/aug‐cc‐pVTZ, and CCSD/aug‐cc‐pVTZ theory levels. After full relaxation of the geometry, the equilibrium structure is found in a planar conformation of C_s symmetry. The difference in the two BF distances suggests the existence of a nonbonded interaction between the fluorine and oxygen atoms. The structural and energetic variation of BF₂OH as a function of the OH torsional angle is considered. The torsional barrier, at the CCSD/aug‐cc‐pVTZ level, and including the effect of the zero‐point energy of the remaining vibrations, is found 2,728 cm^?1. In addition, an anharmonic Hamiltonian for the OH torsional mode is presented and variationally solved. To simplify the treatment and to classify the energy levels, BF₂OH is classified under a G₄ nonrigid group accounting for the inversion symmetry of the molecule and the interchange of the fluorine atoms. The computed torsional energy levels exhibit a very small inversion splitting. The torsional spectrum is simulated considering the dipole moment components along the principal axes of inertia as an explicit function of the torsional coordinate. We observe three dominant bands in the spectrum formed by doublets corresponding to ν₉ = 0 → 1, 2 transitions. The fundamental is an a‐type, Franck–Condon, transition. This is the strongest and appears at 466.80 cm^?1 with relative intensity 0.4312. The ν₉ = 0 → 2 bands correspond to doublets of b‐ and c‐type, i.e., Herzberg–Teller transitions. These are two overlapping bands found at 890.92 and 890.94 cm^?1 with intensity 0.2207 for the b‐type band and 0.2193 for the c‐type band. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Computational thermochemistry of glycolaldehyde

We use a variant of the focal point analysis to refine estimates of the relative energies of the four low‐energy torsional conformers of glycolaldehyde. The most stable form is the cis‐cis structure which enjoys a degree of H‐bonding from hydroxyl H to carbonyl O; here dihedral angles τ₁ (O?C? C? O) and τ₂ (C? C? O? H) both are zero. We optimized structures in both CCSD(T)/aug‐cc‐pVDZ and aug‐cc‐pVTZ; the structures agree within 0.01 Å for bond lengths and 1.0 degrees for valence angles, but the larger basis brings the rotational constants closer to experimental values. According to our extrapolation of CCSD(T) energies evaluated in basis sets ranging to aug‐cc‐pVQZ the trans‐trans form (180°, 180°) has a relative energy of 12.6 kJ/mol. The trans‐gauche conformer (160°, ±75°) is situated at 13.9 kJ/mol and the cis‐trans form (0°, 180°) at 18.9 kJ/mol. Values are corrected for zero point vibrational energy by MP2/aug‐cc‐pVTZ frequencies. Modeling the vibrational spectra is best accomplished by MP2/aug‐cc‐pVTZ with anharmonic corrections. We compute the Watsonian parameters that define the theoretical vibrational‐rotational spectra for the four stable conformers, to assist the search for these species in the interstellar medium. Six transition states are located by G4 and CBS‐QB3 methods as well as extrapolation using energies for structures optimized in CCSD(T)/aug‐cc‐pVDZ structures. We use two isodesmic reactions with two well‐established thermochemical computational schemes G4 and CBS‐QB3 to estimate energy enthalpy and Gibbs energy of formation as well as the entropy of the gas phase system. Our extrapolated electronic energies of species appearing in the isodesmic reactions produce independent values of thermodynamic quantities consistent with G4 and CBS‐QB3. © 2013 Wiley Periodicals, Inc.     

Structural trends of 77Se?1H spin–spin coupling constants and conformational behavior of 2‐substituted selenophenes

Experimental measurements and second‐order polarization propagator approach (SOPPA) calculations of ⁷⁷Se? ¹H spin–spin coupling constants together with theoretical energy‐based conformational analysis in the series of 2‐substituted selenophenes have been carried out. A new basis set optimized for the calculation of ⁷⁷Se? ¹H spin–spin coupling constants has been introduced by extending the aug‐cc‐pVTZ‐J basis for selenium. Most of the spin–spin coupling constants under study, especially vicinal ⁷⁷Se? ¹H couplings, demonstrated a remarkable stereochemical behavior with respect to the internal rotation of the substituent in the 2‐position of the selenophene ring, which is of major importance in the stereochemical studies of the related organoselenium compounds. Copyright © 2009 John Wiley & Sons, Ltd.     

Theoretical investigation of the structures and spectroscopic properties of (H2O4)n (n = 1–4) clusters

Density functional theory and ab initio calculations were performed to elucidate the hydrogen interactions in (H₂O₄)_n (n = 1–4) clusters. The optimized geometries, binding energies, and harmonic vibrational frequencies were predicted at various levels of theory. The trans conformer of the H₂O₄ monomer was predicted to be the most stable structure at the CCSD(T)/aug‐cc‐pVTZ level of theory. The binding energies per H₂O₄ monomer increased in absolute value by 9.0, 10.1, and 11.8 kcal/mol from n = 2 to n = 4 at the MP2/cc‐pVTZ level of theory (after the zero‐point vibrational energy and basis set superposition error corrections). This result implies that the intermolecular hydrogen bonds were stronger in the long‐chain clusters, that is, the formation of the longer chain in the (H₂O₄)_n clusters was more energetically favorable.     

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.     

Difluorobenzenes revisited: an experimental and theoretical study of spin–spin coupling constants for 1,2‐, 1,3‐, and 1,4‐difluorobenzene

The experimental spin–spin coupling constants (SSCCs) for 1,3‐ and 1,4‐difluorobenzene have been determined anew, and found to be consistent with previously determined values. SSCCs for 1,2‐, 1,3‐, and 1,4‐difluorobenzene have been analyzed by comparing them with the coupling constants computed using the second‐order polarization propagator approximation (SOPPA) and the equation‐of‐motion coupled cluster singles and doubles method (EOM‐CCSD). Eighty experimental values have been analyzed using SOPPA calculations, and a subset of 40 values using both SOPPA and EOM‐CCSD approaches. One‐bond coupling constants ¹J(C? C) and ¹J(C? F) are better described by EOM‐CCSD, whereas one‐bond ¹J(C? H) values are better described by SOPPA. An empirical equation is presented which allows for the prediction of unknown coupling constants from computed SOPPA values. A similar approach may prove useful for predicting coupling constants in larger systems. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of the nitrile and methyl torsional vibrations of n‐propyl cyanide in its S0 electronic state

This work presents a structural and vibrational theoretical study of n‐propyl cyanide as a function of the nitrile and methyl torsional modes. A potential energy hypersurface is built at the MP4(SDQ)/aug‐cc‐pVTZ//MP2/aug‐cc‐pVTZ theory level. The equilibrium structure is found in a gauche conformation. Another minimum is found for the trans form. The maximum appears in a cis conformation. For the first time, the interconversion barriers between the different forms are calculated. A two‐dimensional anharmonic vibrational Hamiltonian is built for the nitrile and methyl torsional modes. We find the vibrational energy levels to organize in two stacks associated to the gauche and trans forms. Fundamental frequencies of 113.12 and 220.54 cm^?1 are predicted for the nitrile and methyl torsions in the equilibrium, gauche, conformer. In addition, we find symmetry allowed transitions between the gauche and trans energy levels stacks. The lowest transition is predicted to appear at 24.49 cm^?1. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

First example of the correlated calculation of the one‐bond tellurium–carbon spin–spin coupling constants: Relativistic effects,vibrational corrections,and solvent effects

This work reports on the comprehensive calculation of the NMR one‐bond spin–spin coupling constants (SSCCs) involving carbon and tellurium, ¹J(¹²⁵Te,¹³C), in four representative compounds: Te(CH₃)₂, Te(CF₃)₂, Te(C?CH)₂, and tellurophene. A high‐level computational treatment of ¹J(¹²⁵Te,¹³C) included calculations at the SOPPA level taking into account relativistic effects evaluated at the 4‐component RPA and DFT levels of theory, vibrational corrections, and solvent effects. The consistency of different computational approaches including the level of theory of the geometry optimization of tellurium‐containing compounds, basis sets, and methods used for obtainig spin–spin coupling values have also been discussed in view of reproducing the experimental values of the tellurium–carbon SSCCs. Relativistic corrections were found to play a major role in the calculation of ¹J(¹²⁵Te,¹³C) reaching as much as almost 50% of the total value of ¹J(¹²⁵Te,¹³C) while relativistic geometrical effects are of minor importance. The vibrational and solvent corrections account for accordingly about 3–6% and 0–4% of the total value. It is shown that taking into account relativistic corrections, vibrational corrections and solvent effects at the DFT level essentially improves the agreement of the non‐relativistic theoretical SOPPA results with experiment. © 2016 Wiley Periodicals, Inc.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants. Part 4: Bicycloalkanes

A systematic study of the one‐bond and long‐range J(C,C), J(C,H) and J(H,H) in the series of nine bicycloalkanes was performed at the SOPPA level with special emphasis on the coupling transmission mechanisms at bridgeheads. Many unknown couplings were predicted with high reliability. Further refinement of SOPPA computational scheme adjusted for better performance was carried out using bicyclo[1.1.1]pentane as a benchmark to investigate the influence of geometry, basis set and electronic correlation. The calculations performed demonstrated that classical ab initio SOPPA applied with the locally dense Dunning's sets augmented with inner core s‐functions used for coupled carbons and Sauer's sets augmented with tight s‐functions used for coupled hydrogens performs perfectly well in reproducing experimental values of different types of coupling constants (the estimated reliability is ca 1–2 Hz) in relatively large organic molecules of up to 11 carbon atoms. Additive coupling increments were derived for J(C,C), J(C,H) and J(H,H) based on the calculated values of coupling constants within SOPPA in the model bicycloalkanes, in reasonably good agreement with the known values obtained earlier on pure empirical grounds. Most of the bridgehead couplings in all but one bicycloalkane appeared to be essentially additive within ca 2–3 Hz while bicyclo[1.1.1]pentane demonstrated dramatic non‐additivity of ?14.5 Hz for J(C,C), +16.6 Hz for J(H,H) and ?5.5 Hz for J(C,H), in line with previous findings. Non‐additivity effects in the latter compound established at the SOPPA level should be attributed to the through‐space non‐bonded interactions at bridgeheads due to the essential overlapping of the bridgehead rear lobes which provides an additional and effective non‐bonding coupling path for the bridgehead carbons and their protons in the bicyclopentane framework. Copyright © 2003 John Wiley & Sons, Ltd.     

Highly accurate relativistic universal Gaussian basis set for Dirac–Fock–Breit calculations

Ion mobilities of H₂O⁺ drifting in helium are calculated and compared with experiment. These calculations employ global potential energy surfaces of the H₂O⁺–He complex, which in the present case were calculated ab initio at the unrestricted MP2 level of theory using a basis set of aug‐cc‐pVTZ quality, and treating the ion as a rigid body. Details are presented of the general characteristics of both the ground and first‐excited electronic states of the complex. Although only the ground‐state surface was used for the mobility calculations, the ab initio determination of the ground state necessitated the inclusion of the first‐excited state owing to the presence of a crossing between the two. This crossing is also described. Mobilities calculated from the global surfaces are in good agreement with experiment. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

IR Low‐Temperature Matrix and ab Initio Study on β‐Alanine Conformers

For the first time the argon‐matrix low‐temperature FTIR spectra of β‐alanine are recorded. They reveal a quite complicated spectral pattern which suggests the presence of several β‐alanine conformers in the matrix. To interpret the spectra, the eighteen β‐alanine conformers, stable in the gas phase, are estimated at the B3LYP and MP2 levels combined with the aug‐cc‐pVDZ. Ten low‐energy structures are reoptimized at the QCISD/aug‐cc‐pVDZ and B3LYP and MP2 levels by using the aug‐cc‐pVTZ basis sets. Assignment of the experimental spectra is undertaken on the basis of the calculated B3LYP/aug‐cc‐pVDZ anharmonic IR frequencies as well as careful estimation of the conformer population. The presence of at least three β‐alanine conformers is demonstrated. The detailed analysis of IR spectra points to the possible presence of five additional β‐alanine conformers.