Dipole Derivative, Polarizability Derivative, And Vibrational Polarizability Contribution Output From Gaussian09 With IOp(7/33)

For those itching for polarizability derivative orientation information and wondering where it is when you ask for it… what’s included below is a combination of a few points in one, specifically pointing out that the IOp options are not just “another part” of the Gaussian input file (with the IOp Overlays currently linked HERE).

The problem I realized after an email from Gaussian HQ was that, as was the case for the KMLYP density functional call discussed in previous posts about [18]-annulene, “opt” and “freq” keyword combinations are seen as two distinct runs in Gaussian that don’t pass the IOp information along (and, admittedly, I should have remembered that). Specifically, the additional print-out for the polarizability info is called by IOp(7/33=3).

What I provide below is a two-in-one input file that saves you from having to run double-duty input files in the checkpoint file. This also serves as a template for those looking for examples of combining multi-step input files that include mixed basis sets (as many of the problems I’ve been emailed stem from carriage return issues more than anything else). Note that the input file is set to run Raman intensities and produce higher-precision (hpmodes) eigenvectors (so, if you just want to test this, remove the “raman”).

#p scf=tight opt=tight b3lyp/GEN pseudo=read

C4H5Cl_B3LYP_631Gdp_LanL2DZ_IR_Raman Opt

0 1
 C                 -1.74671095   -0.64168298    0.00000000
 H                 -1.53944096   -1.69141587    0.00000000
 C                 -0.73010315    0.25446188    0.00000000
 H                 -0.93737314    1.30419477    0.00000000
 C                  0.73010315   -0.25446188    0.00000000
 H                  0.93737314   -1.30419477    0.00000000
 C                  1.74671095    0.64168298    0.00000000
 H                  1.53944096    1.69141587    0.00000000
 H                 -3.73526840    0.03531673    0.00000000
 Cl                 3.73526840   -0.03531673    0.00000000

C H 0


#p Geom=Check Guess=Read freq(raman,hpmodes) iop(7/33=3)
C4H5Cl_B3LYP_631Gdp_LanL2DZ_IR_Raman Freq
0 1

Note the carriage return after the second “0 1”.

For the demo molecule above, additional print-out below.

 Dipole derivatives wrt mode   1:  3.96988D-14 -1.15747D-14 -1.96904D-01
 Polarizability derivatives wrt mode          1
                 1             2             3 
      1   0.000000D+00  0.000000D+00  0.206435D+00
      2   0.000000D+00  0.000000D+00  0.143916D-01
      3   0.206435D+00  0.143916D-01  0.000000D+00
 Vibrational polarizability contributions from mode   1       0.0000000       0.0000000       0.0257731
 IFr=  0 A012= 0.23D-23 0.77D+00 0.13D+00 Act= 0.90D+00 DepolP= 0.75D+00 DepolU= 0.86D+00

Alternately, keep track of the checkpoint file.

The Structure Of [18]-Annulene: Computed Raman Spectra, Zero-Point Level And Proton NMR Chemical Shifts

In press (DOI:10.1016/j.molstruc.2012.05.016) in the Journal Of Molecular Structure (Volume 1023, 12 September 2012, Pages 212–215) in the special issue: MOLECULAR VIBRATIONS AND STRUCTURES: THEORY AND EXPERIMENT — A collection of papers dedicated to Professor Jaan Laane on the occasion of his 70th birthday.

This paper on the “actual” geometry of [18]-annulene is part of several larger stories addressing a larger polyene (or larger-polyene) issue. First among these is the meaning of experimental results obtained by various spectroscopic methods (in this case, using previous X-ray, Raman (with the C2 (blue) and D6h (red) simulated spectra shown in the image above), IR, and NMR data that produce different results within the limitations of the methods to study the single molecule). Second is the quality of the theoretical method for reproducing certain types of spectroscopic data. In the case of the [N]-annulene series, the ever-present B3LYP density functional is found to produce the time-average geometry of [18]-annulene found in X-ray data, but another density functional (in this case, KMLYP), finds that bond-alternate minima exist. Third is the importance of the zero-point level in the treatment of systems for which bond-alternate geometries exist with transition-state barriers calculated to be below the zero-point level in the classical approximation of nuclear positions (the Born-Oppenheimer Approximation).

NOTE 1: The KMYLP density functional is called in Gaussian with the following keyword set:

BLYP iop(3/76=1000005570) iop(3/77=0000004430) iop(3/78=0448010000)

NOTE 2: Optimization and Frequency calculations must be performed as TWO SEPARATE CALCULATIONS. The iop-called density functional does not carry itself over between opt + freq (or other properties) in the same input file. If you opt + freq in the same input file, you will Opt with KMLYP but freq with BLYP. This will be obvious by the number of imaginary modes.

Bruce S. Hudson and Damian G. Allis

Abstract. [18]-annulene has been of great interest from the structural point of view of its bond alternation. High-level calculations based on structures selected for agreement with NMR spectra lead to a bond-alternate C2 form over a non-alternating planar D6h structure deduced from diffraction, infrared (IR) and electronic spectral studies. Here it is shown that computed Raman spectra for the D6h and C2 forms are expected to be very different. However, two equivalent non-D6h bond-alternate minima of D3h or C2 geometries are separated by only a small barrier along a motion that involves CC stretching and compression. It is shown here that the zero-point level is above the barrier for this species. In light of that fact, the NMR calculations are reconsidered with inclusion of zero-point level averaging.