CCL:G: Anharmonic Infrared Spectroscopy



Is it plausible that there are degeneracies in the vibrational/rotational spectrum, leading to the divergences you observe? This would be the most obvious explanation for why the numbers blow up for any finite-order perturbation theory. Perhaps you have coupled methyl groups?


With smaller basis sets, oftentimes the divergences in finite order perturbation theory are "masked" (I've seen this with cc-pvdz vs. aug-cc-pVDZ), simply because you have not yet encountered the denominator problem. It's a "better" answer for wrong reasons.


On 6/4/17 4:38 AM, Jens Spanget-Larsen spanget**ruc.dk wrote:

Dear CCL!
When I perform Gaussian calculations with the option freq=anharm, I frequently get strange results (Gaussian 09, Revision D.01). Below is a portion of the results of a b3lyp/6-311++G(d,p) calculation on a planar, aromatic, C2v symmetrical hydrocarbon. The anharmonic results are obviously ridiculous. With smaller basis sets, the results tend to be more realistic. Any explanations?

Yours, Jens >--<

 

     ==================================================
              Anharmonic Infrared Spectroscopy
     ==================================================

 

 Fundamental Bands
 -----------------
    Mode(Quanta)       E(harm)   E(anharm)       I(harm)        I(anharm)
     1(1)              3476.017   3200.902     43.37436356    538.05352307
     2(1)              3202.964   3109.532      0.01728522    255.35137878
     3(1)              3200.221   3033.726      8.04709323    191.91610101
     4(1)              3175.408   3220.259      4.37202654   2730.70615470
     5(1)              2205.565   2153.406      1.97473796     11.71516009
     6(1)              1604.309   1584.143      6.91156654    634.57714516
     7(1)              1432.153   1374.467      3.54182113  12790.38891389
     8(1)              1263.675   1263.359      0.79105395    139.40836060
     9(1)              1112.643   1056.598      5.36740074   1394.13686552
    10(1)              1010.305   1031.530      0.00269831     47.10463873
    11(1)               717.365    793.822      1.28159522 103894.31384900
    12(1)               697.978   1794.196     67.17706473   1444.46852987
    13(1)               505.134   1553.018      6.40712425  16868.05994564
    14(1)               465.540    608.353      2.04849136  24943.02416835
    15(1)               115.463    594.923      0.77231918   6377.20573512
    16(1)               930.354  -1253.748      0.00000000      0.00000039
    17(1)               653.549  -2649.903      0.00000000      0.00000067
    18(1)               621.501 -20992.872      0.00000000      0.01789241
    19(1)               355.466  -1676.263      0.00000000      0.00001474
    20(1)               156.168  -1653.479      0.00000000      0.00002682
    21(1)               990.196   7085.310      0.11097137 176674.25100856
    22(1)               925.704  -2638.620     16.28434626 228593.72377779
    23(1)               811.788   1939.226     41.57436680 251187.21629870
    24(1)               701.424   1649.690     14.23131322 ***************
    25(1)               650.217   -656.428    106.61652268 433910.83299326
    26(1)               493.743  16521.965      0.05960240 ***************
    27(1)               378.821  -1503.535      5.22048529 779479.83775810
    28(1)               121.625  -1543.939      2.61491704  76301.87107925
    .......

    .......

 

  ------------------------------------------------------

  JENS SPANGET-LARSEN        Office:       +45 4674 2710

  Dept. of Science           Fax:          +45 4674 3011

  Roskilde University        Mobile:       +45 2320 6246

  P.O.Box 260                E-Mail:      spanget||ruc.dk

  DK-4000 Roskilde, Denmark  Web: thiele.ruc.dk/~spanget

  ------------------------------------------------------