for atomization energies and geometries , density functional theory ( dft ) constitutes a cost - effective alternative to wavefunction - based _ ab initio _ methods , being capable of accuracies on the order of a few kcal / mol and a few picometer , respectively , if basis sets of polarized triple zeta quality are employed@xcite .    as a result of this astonishing success , especially considering the modest computational cost of dft methods , a variety of new applications are now being explored@xcite , including , very recently , the calculation of molecular anharmonic potential energy surfaces @xcite .
anharmonic potential energy surfaces ( see refs.@xcite for initial pioneering studies at the scf and cisd level ) allow for the prediction of molecular vibration - rotation spectra that can be compared directly to experiment , obviating the need to make estimations or approximations for the effect of anharmonicity .
ccsd(t)/@xmath0 anharmonic force fields can achieve accuracies on the order of 10 @xmath1 or better for fundamentals ( e.g. @xcite and references therein ) , and even greater accuracy is achievable if still larger basis sets and corrections for inner - shell correlation are considered@xcite .
+ dft anharmonic force fields for small molecules have recently been the subject of two validation studies@xcite .
our own validation study@xcite suggested that , for fundamental frequencies , an rms accuracy of about 18 @xmath1 can be attained when using the b97 - 1 functional with sufficiently large basis sets .
the question arises as to how capable dft is to treat medium - sized organic molecules in this fashion . particularly relevant here is the study the study of miani et al.@xcite on the fundamental frequencies of benzene , which employed the b3lyp@xcite functional with the tz2p ( triple - zeta plus double polarization ) basis set .
this paper , as well as a more recent study of furan , pyrrole and thiophene@xcite , in fact yielded even more accurate results on these medium - sized organic systems than suggested by the small - molecule validation studies cited above .
+ the azabenzenes are obtained by systematically replacing ch moieties in benzene by nitrogen atoms .
numerous experimental spectra of these compounds are available in the literature , ranging from mono - azabenzene ( pyridine ) to 1,2,4,5-tetraazabenzene ( s - tetrazine ) .
pentazine is quite elusive@xcite ; the equally elusive n@xmath2 ( hexazine or hexaazabenzene ) molecule has been the subject of extensive theoretical studies@xcite : at the ccsd(t)/cc - pvtz level , the hexazine ( hexaazabenzene ) structure undergoes distortion from the idealized @xmath3 ring to a @xmath4 geometry , which is a local minimum situated some 23 kcal / mol above the @xmath5 diazide global minimum@xcite .
straka recently made the interesting suggestion that cyclic n@xmath2 may form very stable m(@xmath6-n@xmath2 ) complexes with m = ti , zr , hf , th@xcite .    in practice ,
nine molecules are relevant to the present study : the parent benzene molecule , pyridine , pyridazine , pyrimidine , pyrazine , 1,2,3-triazine , 1,2,4-triazine , 1,2,5-triazine ( _ sym_-triazine ) and 1,2,4,5-tetrazine ( s - tetrazine ) .
all molecules are displayed in figure 1 .
+ azabenzene skeletons serve as building blocks in nature , e.g. , pyridine in pyridixol ( vitamin b@xmath2 ) , pyridazine and pyrimidine in pteridine , itself found in folic acid ( vitamin vitamin b@xmath7 ) and riboflavin ( vitamin b@xmath8 ) .
moreover the four dna bases are derivates of pyrimidine ( c , t ) and of purine ( a , g )  itself a pyrimidine - imidazole fused ring system . for these reasons ,
azabenzene - like compounds figure in drug design , e.g. the reverse transcriptase inhibitors azt ( contains pyrimidine ) or nevirapine ( contains pyridine ) .
derivatives of pyridazine have been found to have potential therapeutic or plant growth inhibitory effects .
+ quite different applications include melamine - based plastics ( melamine is a derivative of s - triazine ) and the potential use of s - tetrazine in molecular data storage applications ( one of the relaxation mechanisms of the @xmath9 state leads to 2 hcn+n@xmath8.@xcite ) .
+ innes , ross and moomaw ( irm throughout this paper ) published a compilation and critical review of experimental vibration spectroscopic data current to 1988@xcite .
one of us ( with c. van alsenoy ) carried out a harmonic - only b3lyp / cc - pvtz study of the vibrational frequencies and geometries.@xcite the limitations of this approach are obvious .
+ in the present contribution , we will consider dft _
anharmonic _ force fields and anharmonic vibrational spectra  which are directly comparable to experiment  for the azabenzenes , and demonstrate their power as a spectroscopic assignment tool for medium - sized organic molecules .
in our validation study@xcite , we considered a wide variety of exchange - correlation functionals , as well as convergence in terms of the basis set and the integration grids ( both kohn - sham , ks , and coupled perturbed kohn - sham ) .    especially for organic molecules , we found the best performance to be delivered by the b97 - 1 functional , which is handy s reparametrization@xcite of becke s 1997 hybrid functional@xcite .
satisfactory basis set convergence was generally found to be achieved with the tz2p basis set@xcite .
this is therefore the functional / basis set combination employed for the present study .
we found@xcite results for anharmonic force field calculations to be exceedingly dependent on the ks integration grid , much less so on the cpks grid .
for the present study , we ended up using a ( 200,974 ) grid , that is , the direct product of a 200-point euler - maclaurin radial grid with a 974-point lebedev angular grid . for the cpks steps
, we employed a ( 75,194 ) grid , which considerably reduced the overall computational cost .
neither grid was pruned .
the quartic force fields were calculated by numerical central differences ( in rectilinear normal coordinates ) of analytical second derivatives , using a stand - alone driver program adapted from the cadpac electronic structure program system@xcite .
the actual dft second derivative calculations were done using the gaussian 98 rev . a11 program package@xcite .
the step size in a numerical derivatives calculation always represents a compromise between discretization error and roundoff error .
based on experience , we determined the step size as a functions of both the absolute value of the harmonic frequency associated with the normal coordinate involved , and the associated reduced mass ( @xmath10 ) : @xmath11 in order to reduce roundoff error as much as possible , the ks and cpks equations were basically converged to machine precision .    to calculate this factor , both the reduced mass of the mode @xmath12 and the eigenvalues of the normalisation
matrix @xmath13 are needed .
the last part of equation 1 , depending on the harmonic frequencies @xmath14 , additionally ensures that each displaced geometry has approximately the same energy difference to the minimum geometry .
+ in this manner , we obtain a complete cubic force field , as well as all the diagonal and semidiagonal quartic force constants .
these are sufficient for second - order rovibrational perturbation theory analyses@xcite , which were carried out using the spectro@xcite and polyad@xcite programs developed in the cambridge group and at weizmann , respectively .
the parent molecule of our series , benzene ( figure 1a ) , has been the subject of a few previous anharmonic force field studies , two at the scf / dzp level@xcite , and a very recent one at the b3lyp / tz2p level@xcite .
an extensive experimental literature exists on the subject ( see refs.@xcite for reviews ) : high resolution data are available for many of the fundamentals , and all the assignments can be regarded as conclusive . ( an analysis of the b3lyp / cc - pvtz normal modes of benzene and the azabenzenes in terms of pulay s redundant internal coordinates can be found in table 17 of ref.@xcite . )    experimentally obtaining harmonic frequencies for a molecule this large , even with such high symmetry , is a nearly impossible task . in table
[ tab1 ] we compare computed harmonic frequencies with selected prior calculations ( in particular ccsd(t)/@xmath0 data@xcite ) , as well as several sets of experimentally derived data .
the first such set are the ` @xmath15 ' estimates of handy , murray and amos ( hma)@xcite , themselves obtained as averages of an empirical estimate by goodman , ozkabak , and thakur ( got)@xcite and their own combination of experimental fundamentals with the scf / dzp computed anharmonicities of maslen et al.@xcite .
the second set are derived from the experimental fundamentals and a new scf / dzp anharmonic analysis by handy and willetts ( hw)@xcite . the third set is derived from the experimental fundamentals and the b3lyp / tz2p anharmonicities@xcite , and is expected to be the most reliable . + both sets of dft numbers compare about equally well with the ccsd(t ) data .
the c - h stretching frequencies are underestimated by about 20 @xmath1 , while all other frequencies are being reproduced very accurately .
( as an aside , we note that the c  h underestimation problem also occurs@xcite at the mller - plesset perturbation theory ( mp2 ) level . )
although the c  h stretches are reproduced marginally worse by the b97 - 1 functional ( 3 @xmath1 ) compared to the reference ccsd(t)/@xmath0 values , the overall performance is improved : the mean absolute error decreases from 9.3 ( b3lyp ) to 6.5 ( b97 - 1 ) @xmath1 , and the rms deviation from 10.6 ( b3lyp ) to 9.6 ( b97 - 1 ) .
the smaller decrease for the rms error ( compared to the mean absolute error ) can be attributed to the greater proportional weight given to larger errors .
when excluding the c - h stretching frequencies , the rms errors decrease to 8.4 ( b3lyp ) and 4.3 @xmath1 ( b97 - 1 ) . considering the computational cost of density functional theory ,
the b3lyp results are already of spectacular compared to the reference _ ab initio _ data .
the b97 - 1 functional , however , reduces the rms error of the non c - h stretching frequencies by an additional 50% .
+ both functionals yield excellent agreement with experiment for the fundamental frequencies ( table [ tab2 ] ) . the experimental data were taken from the earlier compilation of got@xcite and from more recent gas - phase measurements@xcite .
the two sets of experimental data agree very well with each other , except for @xmath16 and @xmath17 which are involved in strong fermi resonances . to obtain the b97 - 1/tz2p fundamentals , we had to take the following fermi resonances into consideration : @xmath18 and @xmath19 ( with the deperturbed frequency @xmath19 at 3018 @xmath1 and the combination at 2947 @xmath20 ) , @xmath21 and @xmath17 ( unperturbed @xmath22 at 3040 and the combination band at 3082 @xmath20 ) , @xmath23 and @xmath24 ( with the unperturbed frequencies at 3064 and 3020 @xmath20 ) , and @xmath25 and @xmath16 ( with @xmath16 ( unperturbed ) at 1600 , @xmath25 at 1603 @xmath20 ) .
+ interestingly , because of the small differences in the harmonic frequencies between the b3lyp / tz2p and the b97 - 1/tz2p calculations and the different force field , now new fermi resonances are predicted .
furthermore , the fundamentals affected by these resonances are the ones which differ the most when using the b3lyp or b97 - 1 functional .
because of the fermi resonances , the difference between different functionals and methods becomes more important than on the harmonic level of approximation .
the results of both functionals compared to experiment are , however , very similar , with mean absolute errors ( compared to the latest experimental results of miani et al.@xcite ) of 9.4 ( b3lyp ) and 8.3 ( b97 - 1 ) @xmath1 and rms errors of 17.0 and 13.7 @xmath1 , respectively .
again , all frequencies are very well described with the exception of the c - h stretches , which are underestimated by 1630 @xmath1 ( discounting @xmath24 which is severely perturbed by resonances ) . without the c
- h stretches , the mean absolute errors are reduced to 3.7 ( b3lyp ) and 3.9 ( b97 - 1 ) @xmath1 and the rms errors to 4.8 ( b3lyp ) and 4.5 @xmath1 ( b97 - 1 ) .
hence , while clearly better for harmonic frequencies , the b97 - 1 functional is only slightly better for predicting the fundamental frequencies of benzene .
nevertheless , the accuracy obtained by density functional theory , which yields rms errors of less than 5 @xmath1 for non c - h stretching frequencies , is very good .
+ in the third column , we added the b97 - 1/tz2p anharmonic corrections to the ccsd(t)/ano4321 harmonic frequencies of ref.@xcite , in order to see how such a ` hybrid ' ab initio - dft approach would perform .
the principal improvement seen is for the c  h stretches , which are thus all bought into the range of the experimental values , except for @xmath26 , which is in resonance with @xmath23 as discussed above .
when discounting c  h stretching frequencies , the improvement compared to a pure b97 - 1/tz2p calculation is quite modest , to 3.1 @xmath1 for the mean absolute error and to 4.1 @xmath1 for the rms error .
the corresponding error statistics including c
h stretching frequencies change rather more significantly , to a mean absolute error of 4.0 @xmath1 and an rms error of 7.4 @xmath1 .
+ summarizing , b97 - 1/tz2p is likely to be a useful tool for analysis of vibrational spectra ( and verification of their assignments ) of aromatic organic molecules in general and of the azabenzene series in particular .
pyridine being the simplest azabenzene ( figure 1b ) and the closest to the parent molecule , we may anticipate similar accuracy as for benzene .
as expected for this chemically important species , many experimental spectra are available , in both liquid and gas phases .
most of the experimental results and the latest experimental assignments are summarized by klots@xcite .
the computed and observed vibrational frequencies are presented in table [ tab3 ] .
as expected , the c - h stretching frequencies are not well described . here , they are also heavily perturbed , and , considering performance for the corresponding bands in benzene , our method may not be sufficiently accurate to assist in the experimental assignment . thus , although we will attempt to elaborate on the c - h stretching frequencies for most azabenzenes , the results are more tentative .
additionally , a small change in the original unperturbed fundamental frequencies will cause large changes in the perturbed frequencies .
nevertheless , the two asymmetric stretches @xmath27 and @xmath19 do not follow the general trend of being underestimated by 2040 @xmath1 . while the perturbed @xmath27 is close to the experimental value , @xmath19 differs by 70 @xmath1 .
especially the latter appears too large to be accounted for merely by deficiencies in our dft calculation . here , the severe fermi resonance seems to be responsible for this assignment , as its deperturbed value is at 3021 @xmath1 .
this mode is resonating with ( deperturbed frequencies given ) @xmath27 at 3053 @xmath1 , @xmath28 at 3015 @xmath1 and @xmath29 at 3054 @xmath1 , with the perturbed @xmath28 now appearing at 3023 @xmath1 .
it appears that the latter value is more in line with experiment , and that perhaps this assignment might be more reasonable . nevertheless , some caution is appropriate as the c  h harmonic stretching frequencies are not as well described by dft as the remaining modes .
large basis set ccsd(t ) harmonic frequencies would be helpful here , but would require inordinate amounts of cpu time .
+ another picture emerges for the non c - h frequency modes . here
, we primarily compare to the medium resolution raman vapor data of klots ( klots2)@xcite : their assignments seem to be in line with most other experiments , although they differ from the most recent data obtained by partal et al.@xcite . comparing to these latter inelastic neutron scattering ( ins ) data ,
however , amounts to comparing apples and oranges , as the ins data are based on further refinement of a dft computed _ harmonic _ force field by maximizing agreement between simulated and observed ins spectral .
the pseudo - harmonic frequencies from the refined force field correspond neither to true harmonic nor to true fundamental frequencies  in effect , they are neither fish nor fowl .
+ we can confirm the assignment of klots@xcite of @xmath30 to the lower value of 1053 @xmath1 .
difficult assignments seem to result from the fermi resonances for both @xmath31 and @xmath32 , since these are the only fundamentals which give errors of 10 @xmath1 or larger .
the deperturbed mode of @xmath31 at 1578 @xmath1 resonates with @xmath33 at 1593 @xmath1 resulting in perturbed modes at 1575 and 1596 @xmath1 .
klots may have misassigned the latter band to @xmath31 rather than @xmath33 .
@xmath32 , on the other hand , seems to be very difficult to assign from experimental data because of its depolarized character .
the problem is further exacerbated as @xmath32 is located at the shoulder of @xmath34 .
although @xmath32 is involved in a fermi resonance , it is probably not the source of the disagreement with , since its unperturbed frequency at 1241 @xmath1 resonates with @xmath35 at 1282 @xmath1 , resulting in bands at 1237 and 1289 @xmath1 . with the exception of these two strongly perturbed bands ( @xmath31 and @xmath32 ) ,
the agreement between experiment and theory can only be described as stunning for bands other than c - h stretches . for the latter , probably only a full ccsd(t ) force field
can resolve the assignment .
even when including @xmath31 and @xmath32 ( but excluding all c - h stretching frequencies ) , the mean absolute error for the computed values compared to the experimental data of klots@xcite is 3.6 @xmath1 and the rms error is 5.3 @xmath1 ; without @xmath31 and @xmath32 they are reduced to 2.7 and 3.8 @xmath1 , respectively .
+ the double - harmonic infrared intensities are in good agreement with the b3lyp / cc - pvtz values@xcite and in reasonable agreement with the experimental numbers . even here ,
the c - h stretches differ , probably not only because of anharmonic contributions and fermi resonances , but also because of the poor description of these modes by dft .
still , the excellent performance seen for benzene results is repeated for pyridine .
pyridazine ( figure 1c ) has been of particular interest in the last couple of years @xcite .
a large body of experimental data is available , generally measured in the gas phase but when certain modes were unavailable , the authors of ref.@xcite substituted their own liquid or solid phase measurements in order to get a full complement of frequencies .
all of the spectra have been assigned with the aid of scaled hartree - fock , mp2 , blyp or b3lyp harmonic force fields using small basis sets .
some of them have their force fields fitted to match the experimental data . despite the large amount of data ,
huge discrepancies exist , as shown in table [ tab4 ] .
this seems to be especially true for the out - of - plane modes where experimental assignments and values seem almost arbitrary . for @xmath36 ,
for example , the experimental numbers range from 765 to 949 @xmath1 , and even the two most recent numbers appear at opposite ends of this range .
these assignments have been done using different methods and scaling techniques , since no reliable theoretical data were available at the time .
+ the same problem with the c - h stretches that was encountered for benzene and pyridine occurs for pyridazine as well ; the situation is further complicated since the experimental data differs and thus we can not assign the frequencies of the experiment . despite the discrepancies between various experimental numbers , all the c - h stretches appear plausible within the error range of our method . + as the ins data set dramatically varies from the remainder of the experimental data , and
because it did not appear to be reliable for pyridine , it was excluded from our analysis , although it is reported in table [ tab4 ] for completeness .
+ since such a large amount of experimental data is available , it is best to discuss each mode separately , starting with the in - plane modes .    * for @xmath37 ,
all experiments except the first of ref.@xcite agree on an assignment to the 1570 @xmath1 band .
our results , however , would seem to indicate that the assignment around 1555 @xmath1 might be more plausible .
* the same occurs for @xmath31 .
experimental values of refs.@xcite cluster around 1440 @xmath1 , with refs .
@xcite around 1415 @xmath1 .
our calculation clearly favors the former assignment . * for @xmath24
, significant differences between experiments of refs.@xcite and the computed value of 1145 @xmath1 can be found .
the experimental results are thus untenable in light of this huge discrepancy of more than a hundred @xmath1 , and only ref.@xcite proposes an assignment of 1160 @xmath1 , which still appears to be on the high side . * for @xmath34 ,
most experiments have assignments around the computed 1157 @xmath1 ; however , ref.@xcite assigned this mode to 1120 @xmath1 because of their previous assignment .
our calculated @xmath24 and @xmath34 nearly coincide , * for @xmath38 , all experimental datasets are in agreement with each other and our calculations . considering the generally small anharmonicities for this type of vibrations , the discrepancy of 40 @xmath1 between the computed @xmath39 and the observed @xmath40 seems to be a bit on the large side .
the experimental assignment may have been complicated by the band s position in an ir band envelope going from 960 to 980 @xmath1@xcite . * as with @xmath24 , the experimental assignment for @xmath41 of vazquez et al.@xcite differs from the other obtained fundamentals , and is the only one that can be confirmed by our calculations . * for @xmath16 , only refs.@xcite differ from the computed values , while for @xmath42 , the experimental assignment of stidham and tucci@xcite both in the ir and the raman phase appear doubtful . @xmath43 agrees nicely with all experiments . * for @xmath32 and @xmath22 ,
the same assignment problems that we experienced for @xmath24 and @xmath34 seem to have occurred ; while most experiments assign values above 1100 @xmath1 to @xmath32 , only vazquez et al.@xcite seem to give the correct assignment .
this results , however , in the value of @xmath22 thus ending up a bit on the low side , and only the value of ref.@xcite ends up close to our calculated number . the same seems to have happened to @xmath30 , where only ref.@xcite is in line with our calculations .
hence , for the in - plane modes , the data set of vazquez et al.@xcite appears to be most reliable , excluding @xmath37 and @xmath34 ( which appear to be dubious assignments ) and perhaps @xmath40 and @xmath22 .
* for the out - of - plane modes , the experimental assignment appears to be even more problematic . for @xmath44 ,
all experiments differ by 20 @xmath1 , and are either found at 998 or 842 @xmath1 .
both assignments appear to be unlikely , since the harmonic b97 - 1 frequency is at 981 @xmath1 , which would put the fundamental too low for 998 and definitely way too high for assignment to the 842 @xmath1 band .
* for @xmath45 , all experimental assignments are around the computed 751 @xmath1 , while for @xmath46 , all the experiments appear to yield assignments higher than our computed fundamental of 362 @xmath1 . * based on our data , all the @xmath47 symmetry modes ( @xmath48 through @xmath49 ) should probably be re - examined , since several experimental results differ drastically from each other and from our calculations .
a discussion of these modes is almost impossible , and with the exception of refs .
@xcite for @xmath49 , all experimental assignments lie beyond the likely error bars of the theoretical values .    on the whole , the pyridazine vibrational spectra definitely merit further investigation , as the many experimental and theoretical datasets available are considerably at variance with each other .
thus , it is impossible to give an error estimate for the method used based on this molecule .
however if our computational results prove to be as accurate as for benzene and pyridine , the experimental spectra have to be reexamined . based on this study
, one might wonder whether the addition of more nitrogens or the n - n bond is causing the deterioration of the calculated values .
pyrimidine ( figure 1d ) , has been studied less in recent years .
notable are the ir experiments compiled by irm@xcite , the ir experiments of billes et al.@xcite , and the inelastic neutron scattering work of navarro et al.@xcite . here , we compare mainly to the ir experiments , since the ins data for neither pyridine nor pyridazine agree well with our computed values and the other experiments , for reasons outlined above
. nevertheless , for pyrimidine the ins data seem a lot closer to the ir data than for the previous two azabenzenes ( table [ tab5 ] ) .
the values for the @xmath50 c - h stretches are in agreement with experiment , although this seems fortuitous .
otherwise , all in - plane modes are within 8 @xmath1 of experiment , with the sole exception of @xmath42 .
interestingly , both @xmath38 and @xmath41 undergo strong fermi type 1 and 2 resonances ; their resonance matrices are displayed in table[tab6 ] with their respective eigensolutions . for @xmath38 , very small changes in the unperturbed fundamentals would result in the value ending up not at 1071 @xmath1 , but at the experimental assignment of 1065 @xmath1 .
the assignment of the two resonating modes appears to be purely arbitrary . as for @xmath41 ,
the resonating overtone of @xmath46 matches the experimental number exactly .
since both modes are in resonance , the overtone might have `` borrowed intensity '' ( as described in ref.@xcite ) from @xmath41 , and thus the assignment .
+ when looking at the out - of - plane modes in table [ tab5 ] , only @xmath48 and @xmath22 yield different results compared to the first set of experimental results .
both discrepancies have been mentioned before by the harmonic - only study@xcite , and @xmath22 has a very small intensity .
moreover , by using scaled hf frequencies , wiberg@xcite has suggested that the assignment for @xmath48 might also be incorrect .
the irm@xcite assignments for modes 5 , 6 , 13 and 21 appear more plausible than from ref.@xcite . neglecting the problematic modes 10 and 20 and reassigning mode 9 to 677 @xmath1 , the errors , including the c - h stretches , are 6.3 ( mean absolute ) and 11.4 @xmath1 ( rms ) , and without them are 4.1 and 5.3 @xmath1 , respectively .
thus this method again appears to be very reliable for predicting fundamental frequencies .      as with pyrimidine
, only few experiments have been done for pyrazine ( figure 1e ) .
the vibrational assignments are displayed in table [ tab7 ] .
in ref.@xcite , assignments were done in the gas phase for the ir active modes and in a melt for the raman spectrum .
overall , the agreement is again very good for modes other than c - h stretches .
if the previously mentioned tendency to underestimate the c - h stretches is taken into account , all assignments appear plausible , with the possible exception of the assignment of @xmath36 from ref.@xcite .
it should be noted that all of the c - h stretches are perturbed by fermi resonances . for the in - plane modes
, especially @xmath30 appears to lie beyond the error bars of our method .
it has been mentioned by the harmonic - only study in ref.@xcite that this ` kekul mode ' is problematic and it was concluded that this mode has significant multireference character , since it corresponds to the dissociation into 2 hcn + c@xmath8h@xmath8 .
however , hybrid density functionals are usually capable to describe nondynamical correlation to a certain extent , since for example the atomization energy of ozone is reasonably well reproduced . for this mode ,
the assignment of ref.@xcite appears to be out of question .
+ for all the other modes , most the experiments scatter only by about 5 @xmath1 , with this agreement suggesting that their assignments and values are pretty accurate .
the out - of - plane modes are generally quite well described by theory .
note that the two @xmath51 modes since they are neither ir nor raman active ; for this reason they are not assigned in ref.@xcite .
the ins values are somewhat disappointing for both frequencies , since the chief advantage of this method is the ability to predict such bands .
another discrepancy is the first @xmath52 mode ( @xmath41 ) , although the second experiment of ref.@xcite is closer to our values . +
this leads to mean absolute and rms errors compared to ref.@xcite ( ref.@xcite ) of 7.8 ( 7.9 ) and 10.4 ( 11.1 ) @xmath1 including the c - h stretches , and of 5.4 ( 5.1 ) and 6.9 ( 6.5 ) @xmath1 when excluding them .
if the problematic @xmath53 and @xmath41 are excluded , the mean absolute and rms errors for the non c - h stretches reduces to 4.1 ( 4.1 ) and 5.4 ( 5.5 ) @xmath1 , respectively .
thus , the general rms discrepancy of 5 to 7 @xmath1 is obtained for pyrazine as well .
+      1,3,5-triazine ( figure 1f ) is the only azabenzene with a degenerate point group symmetry ( @xmath54 ) . because of this ,
the @xmath55 modes ( @xmath31 and @xmath24 ) are both ir and raman inactive . the ins @xcite value for @xmath56 coincides exactly with our calculation ( table [ tab8 ] ) ; the discrepancy of 124 @xmath1 for @xmath24 is outlandish even for the ins data and this band is probably misassigned .
the inelastic neutron scattering data@xcite are clearly inferior to the other experimental datasets .
this is hardly surprising as this method relies on fitting and scaling a harmonic force field which is obtained from dft in small split - valence basis sets .
already for the harmonics , this introduces errors on the order of about 50 @xmath1 . the error for the c - h stretches
is expected to be even larger .
it would be interesting to reanalyze their data based on the present anharmonic force fields , although this still would not resolve the problems associated with the c - h stretches , for which large basis set coupled cluster frequencies would be desirable .
the experimental ir / raman data of daunt et al.@xcite and lancaster et al.@xcite are very similar to those of irm@xcite and are not listed in table [ tab8 ] . with the exception of the first @xmath57 @xmath17 and the c - h stretches , which are again underestimated , all computed values are within 11 @xmath1 of experiment again , the errors are similar to the other azabenzenes ; including the c - h stretches , the mean absolute error is 7.6 @xmath1 and the rms error 11 @xmath1 . excluding the stretches ( and @xmath17 ) ,
these numbers reduce to 5.1 ( 3.4 ) and 7.5 ( 5.0 ) @xmath1 , respectively .      for 1,2,3-triazine ( figure 1 g )
, only one ir and raman spectrum in kbr has been experimentally determined@xcite : the assignment was assisted by scaled mp2/6 - 31 g * harmonic frequencies .
the spectrum looks very similar to that of pyrimidine ; even the same fermi resonances are present .
+ at first glance at the spectra in table [ tab9 ] , the agreement between theory and experiment is disappointing .
for very few of the modes is there agreement between the spectra .
however , the authors did not consider _ any _ fermi resonances in the lower spectrum ( although they noted that one might exist , involving @xmath42 ) , and simply assumed the ordering of the fundamentals to be identical to that of the mp2/6 - 31 g * harmonics .
this however leads to obvious misassignments .
for example , @xmath34 and @xmath43 lie almost on top of each other , while @xmath58 borrows intensity from @xmath34 because of the fermi resonance , as can be seen in table [ tab10 ] . if the resonance were not as strong as computed
specifically , if the unperturbed fundamental were shifted only a few @xmath1 towards higher values , the perturbed @xmath58 , presently calculated at 1056 @xmath1 , might easily correspond to the band at 1064 @xmath1 , assigned to @xmath34 in ref.@xcite .
the same actually applies to @xmath22 , which has a strong fermi resonance with @xmath59 + @xmath49 , that has , however , been calculated at a lower value of 644 @xmath1 .
we have also reassigned the first @xmath50 c - h stretch ( mode 1 ) , since the discrepancy between experiment and theory appears too high even for a c - h stretch .
needless to say , this mode also has a strong fermi resonance with the @xmath60 overtone , which has its largest component for the eigenvector with associated eigenvalue at at 3100 @xmath1 , again explaining the experimental value perfectly .
the only other value which is still in question would be a small peak at 1597 @xmath1 , which might correspond to an impurity ; a band was found that was assigned to acetylene , a fragmentation product of 1,2,3-triazine .
the value of 318 @xmath1 , which has been reported for @xmath49 is actually just a scaled hartree - fock harmonic frequency .
+ therefore , a dramatic change in the errors is obtained by reassigning the experimental spectrum . including the c - h stretches in the error analysis , the mean absolute error decreases from 17.7 to 4.7 @xmath1 and the rms error from 27.9 to 6.1 @xmath1 . by excluding the c - h stretches , the mean absolute and rms errors improve from 14.6 and 22.1 @xmath1 , respectively , to 4.1 and 5.2 @xmath1 .
thus , accuracy consistent with the preceding azabenzenes can be obtained following some reassignments in the experimental spectrum .
this also corroborates our conclusion concerning pyridazine that neither substitution of c - h by n nor n - n bonds seem to be responsible for the discrepancies between experimental and computed vibrational spectra .      because 1,2,4-triazine ( figure 1h ) is highly unstable ,
only one `` full '' experimental spectrum has been measured as a liquid film@xcite , although more recently two high - resolution ir spectra have been recorded for four frequencies@xcite . for the c - h stretches ( table [ tab11 ] ) ,
possibly @xmath61 and @xmath53 may have been swapped , since they both take part in a large resonance matrix block .
overall , ten frequencies and combination bands interact forming this large resonance block that include modes 1 to 3 .
hence , it is not surprising that a correct assignment based merely on experimental data would be problematic . with this swap ,
all three c - h stretches are underestimated by the usual 10 - 30 @xmath1 .
only @xmath36 appears questionable in view of our results .
however , @xmath41 has a fermi resonance with the combination of @xmath42 and @xmath30 , and the resonance , eigenvector and eigenvalue matrices are displayed in table [ tab12 ] .
the fermi resonances displayed in these tables are similar to pyrimidine , 1,2,3-triazine and 1,2,4-triazine .
assuming that both of the fermi resonances are slightly stronger than predicted , since @xmath30 is underestimated by 3 @xmath1 compared to the high - resolution ir data , the new @xmath42 + @xmath30 could well end up lower than the predicted 1143 @xmath1 , explaining the value of 1136 @xmath1 that would then not be assigned to any mode .
the assignment of @xmath36 in the experiment would then correspond to @xmath48 , in agreement with the calculated spectrum .
the assignment of @xmath49 also appears to be incorrect .
+ of the four modes obtained by high - resolution ir , the newer data for @xmath17 validates our computational value of 1044 @xmath1 rather than the 1050 @xmath1 predicted by the older experiment . on the other hand , the low - resolution ir value of 769 @xmath1 for mode 19 is closer to the hi - res data compared to our predicted 779 @xmath1 .
these new experiments show the usefulness of the theoretical value , since its resolution is apparently even comparable to that obtained by low - resolution ir data . + without making any reassignments to the experimental spectrum , the mean absolute error from experiment is 13.6 @xmath1 , and the rms error is 19.2 @xmath1 including the c - h stretches .
this appears to be quite large in comparison to the other azabenzenes . with
the abovementioned suggested reassignments , these errors change to 7.5 @xmath1 and 10.0 @xmath1 , respectively . excluding the c - h stretches
, we once again end up with a mean absolute error of 5.1 @xmath1 and an rms error of 5.8 @xmath1 .      for 1,2,4,5-tetrazine ( figure 1i ) , ccsd(t ) harmonic frequencies are available in the literature@xcite .
unfortunately , they were calculated with all electrons correlated in a basis set ( 6 - 311 g * * ) that is only minimal in the inner - shell orbitals .
the resulting errors in harmonic frequencies can easily exceed 20 @xmath1 for small molecules@xcite , several times larger than the actual effect of neglecting core correlation ( typically less than 10 @xmath1 in hcnof systems@xcite ) .
we have recomputed the ccsd(t)/6 - 311 g * * frequencies correlating only valence electrons ( table [ tab13 ] ) . for the in - plane modes , differences with the all - electron calculations
are quite minor ( 6 @xmath1 or less ) .
the out - of - plane modes , however , are drastically affected , up to 60 @xmath1 ( table [ tab13 ] , especially the @xmath62 and the second b@xmath63 frequencies ) .
this adds yet another stanza to the long litany of rejoinders in the literature ( e.g.@xcite ) against the practice of correlating inner - shell electrons in basis sets unsuitable for the purpose .
+ we have recalculated the harmonic frequencies at the ccsd(t ) level ( frozen core ) with basis sets of @xmath0 quality , specifically dunning s popular cc - pvtz @xcite and the almlf - taylor atomic natural orbital@xcite basis set used in ref.@xcite for benzene .
the ccsd(t)/cc - pvtz and ccsd(t)/ano4321 frequencies ( table [ tab13 ] ) are in excellent agreement with each other , the largest difference being 6 @xmath1 for the lowest frequency ( @xmath64 ) .
( this is very unlike the case of benzene@xcite , where the two @xmath52 modes are hypersensitive to the basis set and differ significantly even between cc - pvtz and ano4321 . ) .
the two corresponding modes in s - tetrazine display appreciable differences between 6 - 311 g * and cc - pvtz but barely between cc - pvtz and ano4321 . ) in contrast , the ccsd(t)/6 - 311 g * * frequencies differs significantly from both data sets , particularly for the two @xmath52 modes ( 39 and 33 @xmath1 w.r.t .
ano4321 ) , followed by the lowest @xmath65 mode ( a.k.a . ` kekul mode ' , 27 @xmath1 ) and the @xmath51 mode ( 25 @xmath1 , but a relative error of 7% ) .
so even for a valence - only ccsd(t ) calculation , the 6 - 311 g * * basis set is insufficient , and the rms deviation from the ccsd(t)/ano4321 results , 16.4 @xmath1 ( 17.3 when excluding the c - h stretching frequencies ) , is comparable to that of the much more cost - effective b97 - 1 calculations , 18.1 ( 15.1 ) @xmath1 .
turning now to the b97 - 1/tz2p harmonic frequencies , the most notable differences with the ccsd(t)/ano4321 data are for the ch stretches ( underestimates of 34 @xmath66 for both ) and the kekul mode ( overestimated by 37 @xmath66 ) .
possibly , self - interaction error would adversely affect@xcite the presently used functional s description of the potential surface along the kekul mode , and it would be interesting how a self - interaction corrected functional would perform .
unfortunately no implementation of such a functional is available to the authors , let alone an implementation including analytical derivatives .
exclusive of this mode and the ch stretches , the rms deviation between the ccsd(t)/ano4321 and b97 - 1/tz2p frequencies is about 11 @xmath1 , compared to just 4 @xmath1 for benzene .    as expect ,
agreement between b97 - 1/tz2p fundamentals and experiment is likewise compromised ( table [ tab14 ] ) .
we note that s - tetrazine is both the least stable and the most ` inorganic ' molecule of the series , and that it has a low - lying excited state at less than 2 ev@xcite . if we combine the ccsd(t)/ano4321 harmonic frequencies with the b97 - 1/tz2p anharmonic corrections , the picture brightens somewhat .
compared to the assignment of irm@xcite , the rms deviation amounts to 9.8 @xmath1 if the symmetric c  h stretch is excluded .
the peculiarly large splitting of 76 @xmath1 between irm s symmetric and asymmetric ch stretches is impossible to reconcile with any of our calculations , all of which suggest a very small splitting on the order of just 3 @xmath1 , in agreement with the earlier film ir and raman study of sigworth and pace ( sp)@xcite .
( their band origin of 3089 @xmath1 for the symmetric ch stretch is quite compatible with our calculations : irm s value of 3010 @xmath1 was taken from single - vibronic - level fluorescence spectra@xcite .
an older gas - phase study by franks , merer , and innes ( fmi)@xcite found both ch stretches at 3090 @xmath1 ; the symmetric one in the ir , the antisymmetric one in the raman spectrum . ) our calculations do not reveal a resonance interaction affecting @xmath61 that is severe enough to make a large ch symmetric / antisymmetric spllitting plausible .
we hence conclude that @xmath61 of irm is erroneous .
one other irm assignment which is difficult to reconcile with our best calculations concerns the highest @xmath67 mode ( calculated : 902 , irm : 929 @xmath1 ) for which the franks , merer , and innes@xcite value of 904 @xmath1 nearly perfectly matches our calculation .
it would appear that even the somewhat degraded performance , exhibited by our model for this molecule , still ought to be quite useful for resolving assignment issues or guiding high - resolution spectroscopic measurements .
we have assessed the performance of modern density functional theory ( particularly , b97 - 1/tz2p ) for the anharmonic vibrational spectra of a series of related medium - sized molecules .
the c - h stretching frequencies are consistently underestimated by about 20 - 30 @xmath1 ; disregarding them , the level of agreement with experiment that can be achieved is quite astonishing ( on the order of 6 @xmath1 rms deviation for fundamentals ) .
somewhat poorer accuracy is achieved for pyridazine and particularly for s - tetrazine . for several systems ( particularly 1,2,3-triazine ) , our calculations strongly suggest revised assignments of the observed frequencies . for pyridazine , unfortunately ,
the experimental results are mutually contradictory and our analysis inconclusive . for s - tetrazine
, basically all predicted fundamentals agree less well with experiment ; the more ` inorganic ' character of the molecule may have adversely affected the performance of our calculations .
nevertheless , even here a rms deviation of about 15 @xmath1 for dft and 10 @xmath1 for the combined ccsd(t)/dft results is obtained .
+ overall , b97 - 1/tz2p quartic force fields combined with second - order rovibrational perturbation theory shows great promise for the assignment of vibrational spectra of medium - sized organic molecules .
a. d. boese is grateful for financial support from the feinberg graduate school .
this research was supported by the lise meitner - minerva center for computational quantum chemistry , of which j. m. l. martin is a member , and the helen and martin a. kimmel center for molecular design .
b97 - 1/tz2p quartic force fields in normal coordinates of the molecules studied are available for download ( in spectro format ) at the url http://theochem.weizmann.ac.il/web/papers/azabenzenes.html .    99 boese , a. d. ; handy , n. c. ; _ j. chem .
phys . _ * 114 * 5497 ( 2001 ) .
boese , a. d. ; martin , j. m. l. ; handy , n. c. ; _ j. chem .
phys . _ * 2003 * , _ 119 _ , 3005 .
colwell , s. m. ; handy , n. c. ; lee , a. m. ; _ phys .
* 1996 * , _ 53 _ , 1316 .
jamorski , c. ; casida , m. e. ; salahub , d. r. ; _ j. chem .
phys . _ * 1996 * , _ 104 _ , 5134 .
bauernschmitt , r. ; ahlrichs , r. ; _ chem .
lett . _ * 1996 * , _ 256 _ , 454 .
van caille , c. ; amos , r. d. ; _ chem .
_ * 2000 * , _ 317 _ , 159 .
amos , r. d. ; _ chem .
* 2002 * , _ 364 _ , 612 .
dressler , s. ; thiel , w.;_chem .
lett . _ * 1997 * , _ 273 _ , 71 .
clabo jr . , d. a. ; allen , w. d. ; remington , r. b. ; yamaguchi , y. ; schaefer iii , h. f. ; _ chem .
_ * 1988 * , _ 123 _ , 117
. allen , w. d. ; yamaguchi , y. ; a. g. csszr , clabo , d. a. ; remington , r. b. ; schaefer iii , h. f. ; _ chem .
phys . _ * 1990 * , _ 145 _ , 427 .
martin , j. m. l. ; lee , t. j. ; taylor , p. r. ; franois , j. p. ; _ j. chem .
phys . _ * 1995 * , _ 103 _ , 2589 .
martin , j. m. l. ; lee , t. j. ; and taylor , p. r. ; _ j. chem .
phys . _ * 1998 * , _ 108 _ , 676 .
martin , j. m. l. ; _ chem .
lett . _ * 1998 * , _ 292 _ , 411 .
neugebauer , j. ; hess , b. a. ; _ j. chem .
phys . _ * 2003 * , _ 118 _ , 7215 .
boese , a. d. ; martin , j. m. l. ; submitted .
miani , a. ; cane , e. ; palmieri , p. ; trombetti , a. ; handy , n. c. ; _ j. chem .
_ * 2000 * , _ 112 _ , 248 .
becke , a. d. ; _ j. chem .
phys . _ * 1993 * , _ 98 _ , 5648 .
lee , c. ; yang , w. ; parr , r. g. ; _ phys .
b _ * 37 * , 785 ( 1988 ) burcl , r. ; handy , n. c. ; carter , s.;_spectrochim .
acta a _ * 2003 * , _ 59 _ , 1881 .
politzer , p. ; murray , j. s. ; seminario , j. m. ; miller , r. s. ; _ j. mol .
( theochem ) _ * 1992 * , _ 94 _ , 155 .
saxe , p. ; schaefer iii , h. f. ; _ j. am .
soc . _ * 1983 * , _ 105 _ , 1760 .
tobita , m. ; bartlett , r. j. ; _ j. phys .
* 2001 * , _ 105 _ , 4107 and references therein ; see also http://www.qtp.ufl.edu/~bartlett/downloads/polynitrogen.pdf .
straka , m. ; _ chem .
lett . _ * 2002 * , _ 358 _ , 531 .
scheiner , a. c. ; schaefer iii , h. f. ; _ j. chem . phys . _ * 1987 * , _ 87 _ , 3539 m. schtz , hutter , j. ; lthi , h. p. ; _ j. chem . phys . _ * 1995 * , _ 103 _ , 7048 .
stanton , j. f. ; and gauss , j. ; _ j. chem .
phys . _ * 1996 * , _ 104 _ , 9859 .
rubio , m. ; roos , b. o. ; _ mol .
phys . _ * 1999 * , _ 96 _ , 603 .
nooijen , m. ; _ j. phys .
a _ * 2000 * , _ 104 _ , 4553 .
innes , k. k. ; ross , i. g. ; moomaw , w. r. ; _ j. mol .
spectrosc . _ * 1988 * , _ 132 _ , 492 .
martin , j. m. l. ; van alsenoy , c. ; _ j. phys .
chem . _ * 1996 * , _ 100 _ , 6973 .
hamprecht , f. a. ; cohen , a. j. ; tozer , d. j. ; and handy , n. c. ; _ j. chem .
phys . _ * 1998 * , _ 109 _ , 6264 .
becke , a. d. ; _ j. chem .
phys . _ * 1997 * , _ 107 _ , 8554 dunning , t. h. ; _ j. chem .
phys . _ * 1971 * , _ 55 _ , 716 . the cambridge analytic derivatives package ( cadpac ) , issue 6.5 , cambridge , 1998 ; developed by amos , r. d. , with contributions from alberts , i. l. ; andrews , j. s. ; colwell , s. m. ; handy , n. c. ; jayatikala , d. ; knowles , p.j . ; kobayashi , r.;laidig , k. e. ; laming , g. ; lee , a. m. ; maslen , p. e. ; murray , c. w. ; palmieri , p. ; rice , j. e. ; simandiras , e. d. ; stone , a. j. ; su , m .- d . ; tozer , d. j. frisch , m.j . ; trucks , g.w . ; schlegel , h.b .
; scuseria , g.e . ; robb , m.a . ; cheeseman , j.r . ; zakrzewski , v.g . ; j.a .
montgomery jr . ,
stratmann , r. e. ; burant , j. c. ; dapprich , s. ; millam , j. m. ; daniels , a.d . ;
kudin , k.n . ; strain , m.c . ; farkas , o. ; tomasi , j. ; barone , v. ; cossi , m. ; cammi , r. ; mennucci , b. ; pomelli , c. ; adamo , c. ; clifford , s. ; ochterski , j. ; petersson , g.a . ; ayala , p.y . ; cui , q. ; morokuma , k. ; malick , d.k . ; rabuckm , a.d . ; raghavachari , k. ; foresman , j.b . ; cioslowski , j. ; ortiz , j.v . ; baboul , a.g . ; stefanov , b.b . ; liu , g. ; liashenko , a. ; piskorz , p. ; komaromi , i. ; gomperts , r. ; martin , r.l . ; fox , d.j . ; keith , t. ; m.a .
al - laham , peng , c.y . ; nanayakkara , a. ; gonzalez , c. ; challacombe , m. ; gill , p.m.w . ; johnson , b. ; chen , w. ; wong , m.w . ; andres , j.l . ; gonzalez , c. ; m. head - gordon , replogle , e.s . ; and pople , j.a . ; _
gaussian 98 , revision a. 11 _ ( gaussian , inc . , pittsburgh , pa , 1998 ) .
mills , i. m. ; in : molecular spectroscopy : modern research , eds .
k. narahari rao , c. w. mathews , vol .
1 , academic press , new york 1972 , pp . 115 - 140 .
d. a. clabo jr . ,
allen , w. d. ; remington , r. b. ; yamaguchi , y. ; schaefer iii , h. f. ; _ chem .
phys . _ * 1988 * , _ 123 _ , 187 .
allen , w. d. ; yamaguchi , y. ; csaszar , a. g. ; clabo , d. a. ; remington , r. b. ; schaefer iii , h. f. ; _ chem .
phys . _ * 1990 * , _ 145 _ , 427 .
schneider , w. ; and thiel , w. ; _ chem .
lett . _ * 1989 * , _ 157 _ , 367 .
hoy , a. r. ; mills , i. m. ; and strey , g. ; _ mol .
_ * 1972 * , _ 24 _ , 1265 .
gaw , j. f. ; willets , a. ; green , w. h. ; handy , n. c. ; in _ advances in molecular vibration and collision dynamics _ ; bowman , j. m. , ed . ;
jai : grennwich , ct , 1990 .
martin , j. m. l. ; `` polyad : a vibrational perturbation theory program including arbitrary resonance matrices '' ( weizmann institute of science , reovot , 1997 ) .
maslen , p. e. ; handy , n. c. ; amos , r. d. ; jayatilaka , d. ; _ j. chem .
phys . _ * 1992 * , _ 97 _ , 4233 .
handy , n. c. ; willets , a. ; _ spectrochim .
acta a _ * 1997 * , _ 53 _ , 1169 .
goodman , l. ; ozkabak , a. g. ; thakur , s. n. ; _ j. phys .
chem . _ * 1991 * , _ 95 _ , 9044 .
martin , j. m. l. ; taylor , p. r. ; and lee , t. j. ; _ chem .
_ * 1997 * , _ 275 _ , 414 .
handy , n. c. ; murray , c. w. ; and amos , r. d. ; _ j. phys .
chem . _ * 1993 * , _ 97 _ , 4392 .
handy , n. c. ; maslen , p. e. ; amos , r. d. ; andrews , j. s. ; murray , c. w. ; and laming , g. j. ; _ chem .
lett . _ * 1992 * , _ 197 _ , 506 .
klots , t. d. ; _ spectrochim .
acta a _ * 1998 * , _ 54 _ , 1481 .
partal , f. ; fernandez gomez , m. ; lopez gonzalez , j. j. ; navarro , a. ; kearley , g. h. ; _ chem .
_ * 2000 * , _ 261 _ , 239 .
navarro , a. ; vazquez , j. ; montejo , m. ; lopez gonzalez , j. j. ; kearley , g. j. ; _ chem .
lett . _ * 2002 * , _ 361 _ , 483 .
billes , f. ; mikosch , h. ; and holly , s. ; _ j. mol .
( theochem ) _ * 1998 * , _ 423 _ , 225 .
vazquez , j. ; lopez gozalez , j. j. ; marquez , f. ; boggs , j. e. ; _ j. raman spectrosc . _
* 1998 * , _ 29 _ , 547 .
ozono , y. ; nibu , y. ; shimada , h. ; and shimada , r. ; _ bull .
jpn . _ * 1986 * , _ 59 _ , 2997 .
stidham , h. d. ; and tucci , j. v. ; _ spectrochim .
acta a _ * 1967 * , _ 23 _ , 2233 .
navarro , a. ; fernandez gomez , m. ; lopez gonzalez , j. j. ; fernandez liencres , m. p. ; martinez torres , e. m. ; tomkinson , j. ; kearley , g. j. ; _ j. phys .
* 1999 * , _ 103 _ , 5833 .
burcl , r. ; carter , s. ; and handy , n. c. ; _ chem .
lett . _ * 2003 * , _ 380 _ , 237 .
wiberg , k. b. ; _ j. mol .
struct . _ * 1990 * , _ 244 _ , 61 .
arenas , j. f. ; lopez navarrete , j. t. ; otero , j. c. ; marcos , j. i. ; _ j. chem .
faraday trans .
2 _ * 1985 * , _ 81 _ , 405 .
navarro , a. ; lopez gonzalez , j. j. ; kearley , g. j. ; tomkinson , j. ; parker , s. f. ; and sivia , d. s. ; _ chem .
phys . _ * 1995 * , _ 200 _ , 395 .
kearley , g. j. ; tomkinson , j. ; navarro , a. ; lopez gonzalez , j. j. ; fernandez gomez , m. ; _ chem .
* 1997 * , _ 216 _ , 323 .
daunt , s. j. ; shurvell , h. f. ; _ j. mol .
spectrosc . _ * 1976 * , _ 62 _ , 373 .
lancaster , j. e. ; stamm , r. f. ; colthup , n. b. ; _ spectrochim .
acta _ * 1961 * , _ 17 _ , 155 .
fischer , g. ; nwankwoala , a. u. ; olliff , m. p. ; scott , a. p. ; _ j. mol .
spectrosc . _ * 1993 * , _ 161 _ , 388 .
neunhoffer , h. ; hennig , h. ; _ chemische berichte _ * 1968 * , _ 101 _ , 3952 .
palmer , m. h. ; maier , r. r. j. ; hegelund , f. ; newnham , d. a. ; _ j. mol
. spectrosc . _ * 1998 * , _ 192 _ , 331 .
hegelund , f. ; larsen , r. w. ; nelander , b. ; christen , d. ; palmer , m. h. ; _ j. mol .
spectrosc . _ * 2003 * , _ 217 _ , 9 .
martin , j. m. l. ; schwenke , d. w. ; lee , t. j. ; taylor , p. r. ; _ j. chem .
phys . _ * 1996 * , _ 104 _ , 4657 .
martin , j. m. l. ; _ chem .
lett . _ * 1995 * , _ 242 _ , 343 .
little can be added to this statement : `` using valence - optimized basis sets and including core - correlation is not only a waste of computer time , but a potential source of problems , as it can substantially increase bsse .
this point is not well appreciated : the prevailing view appears to be that no harm can come of correlating the core when the basis set is inadequate .
this is not so . ''
( taylor , p. r. ; in _ lecture notes in quantum chemistry _ ; roos , b. o. , ed . ; springer : berling , 1992 ; p. 406 ) . dunning jr . , t. h. ; _ j. chem .
_ * 90 * , 1007 ( 1989 ) almlf , j. ; taylor , p. r. ; _ j. chem .
phys . _ * 86 * , 4070 ( 1987 ) patchkovskii , s. ; ziegler , t. ; _ j. chem .
phys . _ * 116 * , 7806 ( 2002 ) sigworth , w. d. ; pace , e. l. ; _ spectrochim .
acta a _ * 1971 * , _ 27 _ , 747 brumbaugh , d. v. ; innes , k. k. ; _ chem .
phys . _ * 1981 * , _ 59 _ , 413 .
franks , l. a. ; merer , a. j. ; innes , k. k. ; _ j. mol .
spectrosc . _
* 1968 * , _ 26 _ , 458 .
dilella , d. p. ; and stidham , h. d. ; _ j. raman spectrosc . _ * 1980 * , _ 9 _ , 90 .
wiberg , k. b. ; walters , v. a. ; colson , s. d. ; _ j. phys .
chem . _ * 1984 * , _ 88 _ , 6067 ; wong , k. n. ; colson , s. d. ; _ j. mol .
spectrosc . _ * 1984 * , _ 104 _ , 129 ; walters , v. a. ; snavely , d. a. ; colson , d. a. ; wiberg , k. b. ; wong , k. n. ; _ j. phys .
chem . _ * 1986 * , _ 90 _ , 592 .
kakiuti , y. ; akiyama , m. ; saito , n. ; saito , h. ; _ j. mol .
spectrosc . _ * 1976 * , _ 61 _ , 164 .
califano , s. ; crawford jr .
, b. _ spectrochim . acta _ * 1960 * , _ 16 _ , 900 .
sverdlov , l. m. ; kovner , m. a. ; krainov , e. p. , eds . ; j. wiley and sons : new york , toronto , 1974 ; p. 533
ff .         & b97 - 1 & b3lyp@xcite & ccsd(t)@xcite & hma@xmath68@xcite & hw@xmath69@xcite & mcpth@xmath70@xcite +   + a@xmath71 & 3189 & 3192 & 3210 & 3198 & 3191 & 3218 + & 1008 & 1012 & 1003 & 1001 & 1008 & 1003 + a@xmath72 & 1380 & 1392 & 1380 & 1378 & 1367 & 1392 + b@xmath72 & 1012 & 1017 & 1009 & 1000 & 990 & 1012 + & 718 & 723 & 709 & 712 & 718 & 717 + e@xmath72 & 3165 & 3168 & 3183 & 3182 & 3174 & 3210 + & 1628 & 1635 & 1637 & 1623 & 1607 & 1645 + & 1195 & 1201 & 1194 & 1185 & 1178 & 1197 + & 617 & 624 & 611 & 610 & 613 & 617 + e@xmath71 & 863 & 864 & 865 & 856 & 847 & 861 + a@xmath73 & 686 & 686 & 687 & 680 & 686 & 683 + b@xmath74 & 3155 & 3159 & 3173 & 3173 & 3174 & + & 1022 & 1031 & 1020 & 1016 & 1024 & 1030 +   + b@xmath73 & 1330 & 1333 & 1326 & 1313 & 1318 & 1338 + & 1170 & 1178 & 1163 & 1158 & 1167 & 1163 + e@xmath73 & 987 & 992 & 985 & 978 & 967 & 987 + & 408 & 412 & 406 & 402 & 1058 & 1056 + e@xmath74 & 3180 & 3183 & 3200 & 3186 & 3181 & 3212 + & 1509 & 1519 & 1509 & 1503 & 1494 & 1522 + & 1055 & 1060 & 1056 & 1048 & 1058 & 1057 +   +   +   +     & mode & b97 - 1/tz2p & b3lyp / tz2p & ccsd(t)/b97 - 1 & got@xcite & mcp@xcite +   + a@xmath71 & 1 & 3048 & 3051 & 3069 & 3074 & 3073 + & 2 & 992 & 995 & 987 & 993 & 993 + a@xmath72 & 3 & 1348 & 1351@xmath75 & 1348 & 1350 & 1350 + b@xmath72 & 7 & 987 & 997 & 984 & 990 & 993 + & 8 & 707 & 708 & 698 & 707 & 702 + e@xmath72 & 15 & 3032@xmath75 & 3028 & 3050 & 3057 & 3057 + & 16 & 1611@xmath75 & 1613@xmath75 & 1620 & 1601 & 1610 + & 17 & 1178 & 1181 & 1177 & 1178 & 1178 + & 18 & 613 & 615@xmath75 & 607 & 608 & 608 + e@xmath71 & 11 & 841 & 846 & 843 & 847 & 847 + a@xmath73 & 4 & 672 & 677 & 673 & 674 & 674 + b@xmath74 & 5 & 3004@xmath75 & 2988@xmath75 & 3022 & 3057 & 3057 + & 6 & 1006 & 1015 & 1004 & 1010 & 1014 +   + b@xmath73 & 9 & 1309 & 1305 & 1305 & 1309 & 1309 + & 10 & 1156 & 1163 & 1149 & 1150 & 1148 + e@xmath73 & 19 & 964 & 972 & 962 & 967 & 968 + & 20 & 400 & 403 & 398 & 398 & 398 + e@xmath74 & 12 & 3031@xmath75 & 3023@xmath75 & 3051 & 3064 & 3048 + & 13 & 1486 & 1484 & 1486 & 1484 & 1484 + & 14 & 1045 & 1038 & 1046 & 1038 & 1038 + & & 15.2 ( 4.2 ) & 17.0 ( 4.8 ) + & & 12.2 ( 2.8 ) & 13.7 ( 4.5 ) + & & 7.8 ( 4.5 ) & 7.4 ( 4.1 ) +   +     & & & & & ins + & mode & harmonic&fundamental & ref.@xcite & ref.@xcite & ref.@xcite & ref.@xcite & ref.@xcite & ref.@xcite & ref.@xcite +   + a@xmath76 & 1 & 3193 & 3061@xmath75 & 3068 & & 3082 & 3086 & 3064 & 3070 & 3086 + & 2 & 3166 & 3034@xmath75 & 3056 & & 3052 & 3071 & 3052 & 3053 & 3071 + & 3 & 1604 & 1558 & 1570 & 1555 & 1568 & 1570 & 1572 & 1570 & 1547 + & 4 & 1477 & 1446 & 1415 & 1440 & 1418 & 1444 & 1417 & 1441 & 1465 + & 5 & 1180 & 1145 & & 1340 & 1337 & 1160 & 1347 & 1352 & 1188 + & 6 & 1171 & 1157 & 1159 & 1153 & 1154 & 1119 & 1160 & 1150 & 1177 + & 7 & 1088 & 1064 & 1061 & 1055 & 1061 & 1061 & 1063 & 1063 & 1075 + & 8 & 1006 & 987 & 963 & 960 & 969 & 968 & 964 & 963 & 998 + & 9 & 676 & 669 & 629 & 622 & 622 & 668 & 630 & 632 & 670 + b@xmath8 & 14 & 3180 & 3056@xmath75 & 3085 & & 3082 & 3079 & 3083 & 3080 & 3079 + & 15 & 3162 & 3052@xmath75 & 3056 & & 3080 & 3057 & 3052 & 3041 & 3057 + & 16 & 1601 & 1556 & 1563 & 1540 & 1572 & 1563 & 1566 & 1564 & 1571 + & 17 & 1438 & 1406 & 1446 & 1408 & 1412 & 1413 & 1450 & 1401 & 1436 + & 18 & 1311 & 1284 & 1283 & 1283 & 1281 & 1281 & 1283 & 1287 & 1306 + & 19 & 1081 & 1058 & 1131 & 1112 & 1131 & 1049 & 1129 & 1113 & 1075 + & 20 & 1056 & 1039 & & 1058 & 1058 & 1027 & 1032 & 1052 & 1049 + & 21 & 629 & 622 & 664 & 663 & 673 & 622 & 667 & 660 & 628 +   + b@xmath76 & 22 & 981 & 962 & & & 985 & 987 & 842 & 986 & 998 + & 23 & 762 & 751 & 760 & 760 & 745 & 745 & 759 & 755 & 782 + & 24 & 370 & 362 & 369 & 372 & 369 & 376 & 370 & 370 & 375 + a@xmath8 & 10 & 1019 & 1002 & & & 963 & 1025 & 938 & 970 & 1039 + & 11 & 943 & 926 & & & 765 & 945 & 861 & 786 & 949 + & 12 & 775 & 763 & & & 785 & 729 & 753 & 775 & 793 + & 13 & 372 & 365 & & & 377 & 367 & 410 & 363 & 372 +   +     & & & & ins + & mode & harmonic&fundamental & irm@xcite & ref.@xcite & ref.@xcite +   + a@xmath76 & 1 & 3194 & 3074@xmath75 & 3074 & 3082 & 3074 + & 2 & 3161 & 3045@xmath75 & 3052 & 3053 & 3052 + & 3 & 3147 & 3025@xmath75 & 3038 & 3053 & 3038 + & 4 & 1607 & 1565 & 1570 & 1572 & 1580 + & 5 & 1433 & 1404 & 1398 & 1465 & 1390 + & 6 & 1160 & 1139 & 1147 & 1155 & 1139 + & 7 & 1077 & 1071@xmath75 & 1065 & 1065 & 1080 + & 8 & 1004 & 989 & 991 & 969 & 972 + & 9 & 691 & 685@xmath75 & 678 & 679 & 665 + b@xmath8 & 12 & 3151 & 3040@xmath75 & 3086 & 3047 & 3086 + & 13 & 1605 & 1562 & 1568 & 1569 & 1565 + & 14 & 1494 & 1464 & 1466 & 1411 & 1470 + & 15 & 1393 & 1365 & 1370 & 1356 & 1376 + & 16 & 1252 & 1225 & 1225 & 1224 & 1225 + & 17 & 1211 & 1173 & 1159 & 1158 & 1136 + & 18 & 1091 & 1072 & 1071 & 1071 & 1021 + & 19 & 630 & 623 & 623 & 621 & 628 +   + b@xmath76 & 20 & 1026 & 1006 & 980 & 1033 & 1070 + & 21 & 981 & 962 & 955 & 980 & 980 + & 22 & 822 & 809 & 811 & 804 & 826 + & 23 & 736 & 723 & 721 & 719 & 722 + & 24 & 347 & 340 & 344 & 347 & 344 + a@xmath8 & 10 & 1001 & 982 & [ 927 ] & 960 & 1020 + & 11 & 404 & 396 & 399 & 398 & 344 + & & 5.3 & 21.2 ( 9.9@xmath68 ) & +   +   +     deperturbed & @xmath77 & @xmath45 + @xmath46 + & 1060.0 & + & -10.9 & 1060.1 +   + deperturbed & @xmath59 & 2 @xmath46 + & 682.7 & + & 3.7 & 680.0 +   +   +   +   + perturbed & @xmath77 & @xmath45 + @xmath46 + coefficient of & 1049.2 & 1070.9 + @xmath77 & -0.7085 & -0.7057 + @xmath45 + @xmath46 & -0.7057 & -0.7085 +   + perturbed & @xmath59 & 2 @xmath46 + coefficient of & 677.4 & 685.3 + @xmath59 & -0.5760 & 0.8175 + 2 @xmath46 & 0.8175 & -0.5760 +     & & & & ins + & mode & harmonic&fundamental & irm@xcite & ref.@xcite & ref.@xcite & ref.@xcite +   + a@xmath78 & 1 & 3168 & 3043@xmath75 & 3055 & 3053 & 3054 & 3054 + & 2 & 1612 & 1567 & 1580 & 1579 & 1574 & 1571 + & 3 & 1251 & 1233 & 1233 & 1235 & 1230 & 1240 + & 4 & 1037 & 1020 & 1016 & 1015 & 1015 & 1015 + & 5 & 607 & 598 & 602 & 601 & 596 & 602 + b@xmath79 & 11 & 3147 & 3024@xmath75 & 3040 & 3062 & 3040 & 3040 + & 12 & 1576 & 1533 & 1525 & 1522 & 1529 & 1529 + & 13 & 1372 & 1342 & 1346 & 1353 & 1343 & 1359 + & 14 & 715 & 706 & 704 & 698 & 701 & 698 + b@xmath74 & 15 & 3148 & 2990@xmath75 & 3012 & 3017 & 3018 & 3018 + & 16 & 1513 & 1481 & 1483 & 1484 & 1482 & 1485 + & 17 & 1163 & 1136 & 1130 & 1135 & 1135 & 1062 + & 18 & 1030 & 1013 & 1018 & 1020 & 1020 & 1007 + b@xmath73 & 19 & 3162 & 3055@xmath75 & 3069 & 3069 & 3069 & 3069 + & 20 & 1440 & 1413 & 1411 & 1413 & 1416 & 1427 + & 21 & 1214 & 1174 & 1149 & 1337 & 1146 & 1146 + & 22 & 1084 & 1061 & 1063 & 1063 & 1062 & 1131 +   + a@xmath80 & 6 & 1002 & 982 & [ 960 ] & & [ 997 ] & 1065 + & 7 & 347 & 338 & [ 350 ] & & [ 422 ] & 350 + b@xmath71 & 8 & 946 & 927 & 927 & 925 & 919 & 940 + b@xmath72 & 9 & 987 & 968 & 983 & 976 & 983 & 983 + & 10 & 775 & 762 & 756 & 755 & 754 & 754 + b@xmath63 & 23 & 804 & 791 & 785 & 785 & 786 & 818 + & 24 & 429 & 421 & 418 & 417 & 420 & 414 + & & 6.9 ( 5.3@xmath68 ) & 7.1 ( 6.7@xmath68 ) & 6.5 ( 5.5@xmath68 ) & +   +   +     & & & & ins + & mode & harmonic&fundamental & ref.@xcite & irm@xcite & ref.@xcite +   + a@xmath81 & 1 & 3164 & 3025@xmath75 & 3082 & 3042 & 3042 + & 2 & 1156 & 1137 & 1130 & 1137 & 1125 + & 3 & 1005 & 989 & 989 & 989 & 991 + a@xmath82 & 4 & 1400 & 1375 & [ 1556 ] & & 1375 + & 5 & 1168 & 1124 & [ 1381 ] & & 1000 + e & 8 & 3159 & 3034@xmath75 & 3081 & 3059 & 3056 + & 9 & 1594 & 1545 & 1560 & 1556 & 1555 + & 10 & 1439 & 1405 & 1404 & 1410 & 1414 + & 11 & 1194 & 1167 & 1168 & 1173 & 1165 + & 12 & 688 & 675 & 666 & 676 & 675 +   + a@xmath83 & 6 & 947 & 932 & & 925 & 940 + & 7 & 755 & 740 & & 737 & 732 + e  & 13 & 1043 & 1017 & & 1034 & 1030 + & 14 & 345 & 338 & & 339 & 333 + & & & 7.5 ( 5.0@xmath68 ) & +   +   +     & & & ir & raman & ir reassigned + & mode & harmonic&fundamental & ref.@xcite & ref.@xcite & +   + a@xmath76 & 1 & 3195 & 3032@xmath75 & 3107 & 3110 & 3046 + & 2 & 3165 & 3028 & & 3045 & + & 3 & 1589 & 1548 & 1597 & 1594 & 1546 + & 4 & 1384 & 1347 & 1336 & 1329 & 1336 + & 5 & 1136 & 1122 & 1080 & 1088 & 1124 + & 6 & 1092 & 1081@xmath75 & 1069 & 1064 & 1080 + & 7 & 1006 & 985 & 979 & 977 & 979 + & 8 & 675 & 667 & 660 & 660 & 660 + b@xmath8 & 14 & 3171 & 3042@xmath75 & 3046 & & 3046 + & 15 & 1584 & 1538 & 1545 & 1547 & 1545 + & 16 & 1438 & 1410 & 1410 & & 1410 + & 17 & 1219 & 1199 & 1195 & 1198 & 1195 + & 18 & 1101 & 1078 & 1124 & 1127 & 1080 + & 19 & 974 & 929 & 935 & & 935 + & 20 & 662 & 659@xmath75 & 653 & & 660 +   + b@xmath76 & 10 & 1017 & 998 & & & + & 11 & 824 & 811 & 819 & & 819 + & 12 & 783 & 769 & 769 & & 769 + & 13 & 302 & 296 & [ 318 ] & & + a@xmath8 & 8 & 992 & 973 & & & + & 9 & 358 & 350 & & 365 & + & & 22.1 & & 5.2 +   +     deperturbed & @xmath84 & @xmath17 + @xmath49 + & 1074.5 & + & 10.6 & 1061.4 +   + deperturbed & @xmath22 & @xmath41 + @xmath49 + & 655.1 & + & -6.3 & 647.6 +   +   +   +   + perturbed & @xmath84 & @xmath17 + @xmath49 + coefficient of & 1055.5 & 1080.7 + @xmath84 & -0.4804 & -0.8770 + @xmath17 + @xmath49 & 0.8770 & -0.4804 +   + perturbes & @xmath22 & @xmath59 + @xmath49 + coefficient of & 644.0 & 658.7 + @xmath59 & 0.4939 & -0.8695 + @xmath59 + @xmath49 & 0.8175 & 0.4939 +     & & & ir & hires - ir + & mode & harmonic&fundamental & ref.@xcite & refs.@xcite +   + a & 1 & 3178 & 3051@xmath75 & 3090 & + & 2 & 3175 & 3075@xmath75 & 3060 & + & 3 & 3150 & 3007@xmath75 & 3035 & + & 4 & 1596 & 1556 & 1560 & + & 5 & 1563 & 1524 & 1529 & + & 6 & 1470 & 1438 & 1435 & + & 7 & 1401 & 1376 & 1380 & + & 8 & 1314 & 1286 & 1295 & + & 9 & 1179 & 1171@xmath75 & 1163 & + & 10 & 1139 & 1107 & 1136 & + & 11 & 1096 & 1067@xmath75 & 1113 & + & 12 & 1063 & 1044 & 1050 & 1043.5 + & 13 & 1011 & 989 & 955 & + & 14 & 726 & 717 & 713 & + & 15 & 632 & 629@xmath75 & & +   + a  & 16 & 1014 & 1001 & & + & 17 & 990 & 976 & & + & 18 & 862 & 853 & 851 & + & 19 & 791 & 779 & 768 & 768.7 + & 20 & 379 & 372 & & 367.9 + & 21 & 316 & 308 & & 311.3 + & & 16.1(5.8@xmath68 ) & +   +   +     deperturbed & @xmath36 & @xmath42 + @xmath30 + & 1160.4 & + & -13.7 & 1153.4 +   + deperturbed & @xmath19&2 @xmath30 + & 622.8 & + & 8.4 & 617.7 +   +   +   +   + perturbed & @xmath36 & @xmath42 + @xmath30 + coefficient of & 1142.7 & 1171.0 + @xmath36 & 0.6136 & -0.7896 + @xmath42 + @xmath30 & 0.7896 & 0.6136 +   + perturbed & @xmath19&2 @xmath30 + coefficient of & 611.4 & 629.0 + @xmath59 & -0.5949 & 0.8038 + 2 @xmath30 & 0.8038 & -0.5949 +     & & cc - pvtz & ano-4321 & cc - pvtz & tz2p + & @xcite & full@xcite & & @xcite & +   + a@xmath85 & 3253 & 3223 & 3221 & 3223 & 3226 & 3196 & 3192 + & 1524 & 1456 & 1452 & 1455 & 1457 & 1472 & 1473 + & 1059 & 1018 & 1018 & 1025 & 1025 & 1050 & 1047 + & 758 & 741 & 740 & 747 & 745 & 756 & 751 + b@xmath74 & 3251 & 3222 & 3218 & 3222 & 3225 & 3199 & 3191 + & 1250 & 1212 & 1210 & 1218 & 1216 & 1232 & 1231 + & 1109 & 1088 & 1088 & 1096 & 1098 & 1099 & 1096 + b@xmath73 & 1505 & 1471 & 1465 & 1472 & 1470 & 1479 & 1478 + & 1167 & 1137 & 1134 & 1138 & 1137 & 1151 & 1147 + & 845 & 908 & 905 & 932 & 927 & 964 & 969 + b@xmath79 & 1611 & 1558 & 1552 & 1561 & 1561 & 1559 & 1560 + & 1345 & 1321 & 1319 & 1323 & 1323 & 1330 & 1325 + & 648 & 637 & 635 & 640 & 640 & 656 & 650 +   + a@xmath86 & 344 & 257 & 319 & 344 & 346 & 354 & 344 + b@xmath72 & 985 & 947 & 957 & 998 & 996 & 1011 & 999 + & 803 & 772 & 783 & 816 & 816 & 837 & 829 + b@xmath63 & 930 & 896 & 904 & 919 & 921 & 940 & 930 + & 297 & 228 & 270 & 276 & 270 & 258 & 250 +     & mode & b97 - 1 & ccsd(t)/b97 - 1 & irm@xcite & ref.@xcite & sp@xcite & fmi@xcite +   + a@xmath85 & 1 & 3073 & 3107 & 3010 & 3040 & 3089 & 3090 + & 2 & 1424 & 1408 & 1415 & 1489 & 1417@xmath87 & 1418@xmath87 + & 3 & 1024 & 1000 & 1009 & 990 & 1015 & 1017 + & 4 & 741 & 735 & 736 & 737 & 734 & 736 + b@xmath74 & 11 & 3070 & 3104 & 3086 & 3070 & 3086 & 3090 + & 12 & 1204 & 1189 & 1204 & 1200 & 1204 & 1200 + & 13 & 1079 & 1081 & 1093 & 1106 & 1109@xmath70 & 1103@xmath70 + b@xmath73 & 14 & 1445 & 1437 & 1448 & 1434 & 1448 & 1440 + & 15 & 1112@xmath75 & 1102 & 1108 & 1187 & 1151 & + & 16 & 927 & 885 & 883 & 1085 & 893@xmath70 & 881@xmath70 + b@xmath79 & 8 & 1511 & 1511 & 1525 & 1543 & 1523@xmath87 & 1521@xmath87 + & 9 & 1294 & 1293 & 1290 & 1278 & 1302 & 1303 + & 10 & 645 & 638 & 640 & 679 & 649 & 651 +   + a@xmath86 & 5 & 336 & 338 & 335 & 319 & 335 & 335 + b@xmath72 & 6 & 977 & 974 & 994 & 925 & 1015 & + & 7 & 814 & 801 & 801 & 775 & 799 & 800 + b@xmath63 & 17 & 911 & 902 & 929 & 890 & 929 & 904 + & 18 & 244 & 261 & 254 & 340 & 254 & 254 + & 15.2@xmath68(15.1,14.8@xmath69 ) & 54.8 ( 57.5 ) & 19.0 ( 19.4 ) & 16.0 ( 15.6 ) + & 12.2@xmath68(11.3,9.4@xmath69 ) & 63.5 ( 64.5 ) & 20.1 ( 20.3 ) & 11.0 ( 10.2 ) +   +   +   +   +   +