The SelectNormalModes and SelectAnharmonicModes options require additional input. The modes to select are specified in a separate blankline terminated input section. The initial mode list is always empty.
Integers and integer ranges without a keyword are interpreted as mode numbers; although,This option is only available for methods with the [not]mode keywords may be used. The keywords atoms and notatoms can be used to define an atom list whose modes should be included/excluded (respectively). Atoms can also be specified by ONIOM layer via the [not]layer keywords, which accept these values: real for the real system, model for the model system in a 2layer ONIOM, middle for the middle layer in a 3layer ONIOM, and small for the model layer of a 3layer ONIOM. Atoms may be similarly included/excluded by residue with residue and notresidue, which accept lists of residue names or numbers. Both keyword sets function as shorthand forms for atom lists.
Here are some examples:
Retrieving Force Constants
ReadFC
Requests that the force constants from a previous frequency calculation be read from the checkpoint file, and the mode and thermochemical analysis be repeated, presumably using a different temperature, pressure, or isotopes, at minimal computational cost. Note that since the basis set is read from the checkpoint file, no general basis should be input. If the Raman option was specified in the previous job, then do not specify it again when using this option.
Requesting Specific Spectra
Raman
Compute Raman intensities in addition to IR intensities. This is the default for HartreeFock. It may be specified for DFT and MP2 calculations. For MP2, Raman intensities are produced by numerical differentiation of dipole derivatives with respect to the electric field (span class=”kwit”>Raman is equivalent to NRaman for this method).
NRaman
Compute polarizability derivatives by numerically differentiating the analytic dipole derivatives with respect to an electric field. This is the default for MP2 if Freq=Raman.
NNRaman
Compute polarizability derivatives by numerically differentiating the analytic polarizability with respect to nuclear coordinates.
NoRaman
Skips the extra steps required to compute the Raman intensities during HartreeFock analytic frequency calculations, saving 1030% in CPU time.
VCD
Compute the vibrational circular dichroism (VCD) intensities in addition to the normal frequency analysis. This option is valid for HartreeFock and DFT methods. This option also computes optical rotations (see Polar=OptRot).
ROA
Compute dynamic analytic Raman optical activity intensities using GIAOs [Cheeseman11a]. This procedure requires one or more incident light frequencies to be supplied in the input to be used in the electromagnetic perturbations (CPHF=RdFreq is the default with Freq=ROA). This option is valid for HartreeFock and DFT methods. NNROA says to use the numerical ROA method from Gaussian 03; this is useful only for reproducing the results of prior calculations.
Anharmonic Frequency Analysis
Anharmonic
Do numerical differentiation along modes to compute zeropoint energies, anharmonic frequencies, and anharmonic vibrationalrotational couplings if VibRot is also specified. This option is only available for methods with analytic second derivatives: HartreeFock, DFT, CIS, and MP2.
ReadAnharm
Read an input section with additional parameters for the vibrationalrotational coupling and/or anharmonic vibrational analysis (VibRot or Anharmonic options). Available input options are documented below following the examples.
ReadHarmonic
Read the central point force constants and normal modes from a previous harmonic frequency calculation and avoid repeating the calculation at the central pointn.
ReadDifferentNtharmonic
Read the central point energy, forces, and force constants from a previous calculation and then compute 3^{rd} and 4^{th} derivatives at the current (presumably lower) level of theory for anharmonic spectra.
SelectAnharmonicModes
Read an input section selecting which modes are used for differentiation in anharmonic analysis. The format of this input section is discussed above. SelAnharmonicModes is a synonym for this option.
Vibronic Spectra: FranckCondon, HerzbergTeller and FCHT
The following options perform an analysis for an electronic excitation using the corresponding method; these jobs use vibrational analysis calculations for the ground state and the excited state to compute the amplitudes for electronic transitions between the two states. The vibrational information for the ground state is taken from the current job (Freq or Freq=ReadFC), and the vibrational information for the excited state is taken from a checkpoint file, whose name is provided in a separate input section (enclose the path in quotes if it contains internal spaces). The latter will be from a CISingles or TDDFT Freq=SaveNormalModes calculation.
The ReadFCHT option can be added to cause additional input to be read to control these calculations (see below). In the latter case, the excited state checkpoint file would typically have been generated with Freq=(SelectNormalModes, SaveNormalModes) with the same modes selected.
If CIS frequencies are to be used with the HerzbergTeller or FranckCondonHerzbergTeller analysis, the CIS frequencies must be done numerically (Freq=Numer rather than Freq). This is because the transition dipole derivatives are not computed during the analytic force constant evaluation. The corresponding HF frequency calculation on the ground state, which is also required, can be done analytically as usual.
FranckCondon
Use the FranckCondon method [Sharp64, Doktorov77, Kupka86, Zhixing89, Berger97, Peluso97, Berger98, Borrelli03, Weber03, Coutsias04, Dierksen04, Lami04, Dierksen04a, Dierksen05, Liang05, Jankowiak07, Santoro07, Santoro07a, Barone09] (the implementation is described in [Santoro07, Santoro07a, Santoro08, Barone09]). FC is a synonym for this option. Transitions for ionizations can be analyzed instead of excitations. In this case, the molecule specification corresponds to the neutral form, and the additional checkpoint file named in the input section corresponds to the cation.
HerzbergTeller
Use the HerzbergTeller method [Herzberg33, Sharp64, Small71, Orlandi73, Lin74, Santoro08] (the implementation is described in [Santoro08]). HT is a synonym for this option.
FCHT
Use the FranckCondonHerzbergTeller method [Santoro08].
Emission
Indicates that emission rather than absorption should be simulated for a FranckCondon and/or HerzbergTeller analysis. In this case, within the computation, the initial state is the excited state, and the final state is the ground state (although,This option allows you to specify alternatives to the default temperature, pressure, frequency scale factor the sources of frequency data for the ground and excited state are as described above: current job=ground state, second checkpoint file=excited state).
ReadFCHT
Read an input section containing parameters for the calculation. Available input options are documented below following the examples. This input section precedes that for ReadAnharmon if both are present.
Other Calculation Variations and Properties
VibRot
Analyze vibrationalrotational coupling.
Projected
For a point on a massweighted reaction path (IRC), compute the projected frequencies for vibrations perpendicular to the path. For the projection, the gradient is used to compute the tangent to the path. Note that this computation is very sensitive to the accuracy of the structure and the path [Baboul97]. Accordingly, the geometry should be specified to at least 5 significant digits. This computation is not meaningful at a minimum.
TProjected
Perform a projected harmonic frequency analysis if the RMS force is ≥ 1.d3 Hartree/Bohr and perform regular harmonic analysis if the RMS force is smaller.
HinderedRotor
Requests the identification of internal rotation modes during the harmonic vibrational analysis [McClurg97, Ayala98, McClurg99]. If any modes are identified as internal rotation, hindered or free, the thermodynamic functions are corrected. The identification of the rotating groups is made possible by the use of redundant internal coordinates. Because some structures, such as transition states, may have a specific bonding pattern not automatically recognized, the set of redundant internal coordinates may need to be altered via the Geom=Modify keyword. Rotations involving metals require additional input via the ReadHinderedRotor option (see below).
If the force constants are available on a previously generated checkpoint file, additional vibrational/internal rotation analyses may be performed by specifying Freq=(ReadFC, HinderedRotor). Since Opt=CalcAll automatically performs a vibrational analysis on the optimized structure, Opt=(CalcAll, HinderedRotor) may also be used.
ReadHinderedRotor
Causes an additional input section to be read containing the rotational barrier cutoff height (in kcal/mol) and optionally the periodicity, symmetry number and multiplicity for rotational modes. Rotations with barrier heights larger than the cutoff value will be automatically frozen. If the periodicity value is negative, then the corresponding rotor is also frozen. You must provide the periodicity, symmetry and spin multiplicity for all rotatable bonds contain metals. The input section is terminated with a blank line, and has the following format:
25 
Includes modes 2 through 5. 
atoms=O 
Includes modes involving oxygen atoms. 
120 atoms=Fe 
Includes modes 1 through 20 and any modes involving iron atoms. 
layer=real notatoms=H 
Includes modes for heavy atoms in low layer (subject to default threshold). 
VMaxvalue 

Atom1 Atom2 periodicity symmetry spin 
Repeated as necessary. 
… 

Normal Modes
HPModes
Include the high precision format (to five figures) vibrational frequency eigenvectors in the frequency output in addition to the normal threefigure output.
InternalModes
Print modes as displacements in redundant internal coordinates. IntModes is a synonym for this option.
SaveNormalModes
Save all modes in the checkpoint file. SaveNM is a synonym for this option. NoSaveNormalModes, or NoSaveNM, is the default.
ReadNormalModes
Read saved modes from the checkpoint file. ReadNM is a synonym for this option. NoReadNormalModes, or NoReadNM, is the default.
SelectNormalModes
Read input selecting the particular modes to display. SelectNM is a synonym for this option. NoSelectNormalModes, or NoSelectNM, is the default. AllModes says to include all modes in the output. The format of this input section is discussed above. Note that this option does not affect the functioning of SaveNormalModes, which always saves all modes in the checkpoint file.
SortModes
Sort modes by ONIOM layer in the output.
ModelModes
Display only modes involving the smallest model system in an ONIOM calculation.
MiddleModes
Display only modes involving the two model systems in a 3layer ONIOM.
PrintDerivatives
Print normal mode derivatives of the dipole moment, polarizability, and so on.
PrintFrozenAtoms
By default, the zero displacements for frozen atoms are not printed in the mode output. This option requests that all atoms be listed.
NoPrintNM
Used to suppress printing of the normal mode components during a frequency calculation. The frequencies and intensities are still reported for each mode.
GeometryRelated Options
ModRedundant
Readin modifications to redundant internal coordinates (i.e., for use with InternalModes). Note that the same coordinates are used for both optimization and mode analysis in an Opt Freq, for which this is the same as Opt=ModRedundant. See the discussion of the Opt keyword for details on the input format.
ReadIsotopes
This option allows you to specify alternatives to the default temperature, pressure, frequency scale factor and/or isotopes—298.15 K, 1 atmosphere, no scaling, and the most abundant isotopes (respectively). It is useful when you want to rerun an analysis using different parameters from the data in a checkpoint file.
Be aware, however, that all of these can be specified in the route section (Temperature, Pressure and Scale keywords) and molecule specification (the Iso parameter), as in this example:
#T Method/631G(d) JobType Temperature=300.0 …
…
0 1
C(Iso=13)
…
ReadIsotopes input has the following format:
temp pressure [scale] 
Values must be real numbers. 
isotope mass for atom 1


isotope mass for atom 2


…


isotope mass for atom n


Where temp, pressure, and scale are the desired temperature, pressure, and an optional scale factor for frequency data when used for thermochemical analysis (the default is unscaled). The remaining lines hold the isotope masses for the various atoms in the molecule, arranged in the same order as they appeared in the molecule specification section. If integers are used to specify the atomic masses, the program will automatically use the corresponding actual exact isotopic mass (e.g., 18 specifies ^{18}O, and Gaussian uses the value 17.99916).
Analytic frequencies are available for the AM1, PM3, PM3MM, PM6, PDDG, DFTB, DFTBA, HF, DFT, MP2, CIS, TD and CASSCF methods.
Numerical frequencies are available for MP3, MP4(SDQ), CID, CISD, CCD, CCSD, EOMCCSD and QCISD.
Raman is available for the HF, DFT and MP2 methods.
VCD and ROA are available for HF and DFT methods.
Anharmonic is available for HF, DFT, MP2 and CIS methods.
Freq and NMR can both be on the same route for HF and DFT.
Frequency Output. The basic components of the output from a frequency calculation are discussed in detail in chapter 4 of Exploring Chemistry with Electronic Structure Methods [Foresman15].
New Gaussian users are often surprised to see that the final part frequency calculation output that looks that of a geometry optimization at the beginning of a frequency job:
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Initialization pass.
Link 103, which performs geometry optimizations, is executed at the beginning and end of all frequency calculations. This is done so that the quadratic optimization step can be computed using the correct second derivatives. Occasionally an optimization will complete according to the normal criterion using the approximate Hessian matrix, but the step size is actually larger than the convergence criterion when the correct second derivatives are used. The next step is printed at the end of a frequency calculation so that such problems can be identified. If you think this concern is applicable, use Opt=CalcAll instead of Freq in the route section of the job, which will complete the optimization if the geometry is determined not to have fully converged (usually, given the full second derivative matrix near a stationary point, only one additional optimization step is needed), and will automatically perform a frequency analysis at the final structure.
Specifying #P in the route section produces some additional output for frequency calculations. Of most importance are the polarizability and hyperpolarizability tensors (the latter in Raman calculations only); although, they still may be found in the archive entry in normal printlevel jobs. They are presented in lower triangular and lower tetrahedral order, respectively (i.e., α_{xx}, α_{xy}, α_{yy}, α_{xz}, α_{yz}, α_{zz} and β_{xxx}, β_{xxy}, β_{xyy}, β_{yyy}, β_{xxz}, β_{xyz}, β_{yyz}, β_{xzz}, β_{yzz}, β_{zzz}), in the standard orientation:
Dipole = 2.37312183D16 6.66133815D16 9.39281319D01
Polarizability= 7.83427191D01 1.60008472D15 6.80285860D+00
3.11369582D17 2.72397709D16 3.62729494D+00
HyperPolar = 3.08796953D16 6.27350412D14 4.17080415D16
5.55019858D14 7.26773439D01 1.09052038D14
2.07727337D+01 4.49920497D16 1.40402516D13
1.10991697D+01
#P also produces a bargraph of the simulated spectra for small cases.
Thermochemistry analysis follows the frequency and normal mode data:
Zeropoint correction= .023261 (Hartree/Particle) 

Thermal correction to Energy= .026094 

Thermal correction to Enthalpy= .027038 

Thermal correction to Gibbs Free Energy= .052698 

Sum of electronic and zeropoint Energies=527.492585 
E_{0}=E_{elec}+ZPE 
Sum of electronic and thermal Energies= 527.489751 
E= E_{0}+ E_{vib}+ E_{rot}+E_{trans} 
Sum of electronic and thermal Enthalpies=527.488807 
H=E+RT 
Sum of electronic and thermal Free Energies=527.463147 
G=HTS 
The raw zeropoint energy correction and the thermal corrections to the total energy, enthalpy, and Gibbs free energy (all of which include the zeropoint energy) are listed, followed by the corresponding corrected energy. The analysis uses the standard expressions for an ideal gas in the canonical ensemble. Details can be found in McQuarrie [McQuarrie73] and other standard statistical mechanics texts. In the output, the various quantities are labeled as follows:
E (Thermal) 
Contributions to the thermal energy correction 
CV 
Constant volume molar heat capacity 
S 
Entropy 
Q 
Partition function 
The thermochemistry analysis treats all modes other than the free rotations and translations as harmonic vibrations. For molecules having hindered internal rotations, this can produce slight errors in the energy and heat capacity at room temperatures and can have a significant effect on the entropy. The contributions of any very low frequency vibrational modes are listed separately so that their harmonic contributions can be subtracted from the totals and their correctly computed contributions included should they be group rotations and high accuracy is required. Expressions for hindered rotational contributions to these terms can be found in Benson [Benson68]. The partition functions are also computed, with both the bottom of the vibrational well and the lowest (zeropoint) vibrational state as reference.
Preresonance Raman. This calculation type is requested with one of the Raman options in combination with CPHF=RdFreq. The frequency specified for the latter should be chosen as follows:

Determine the difference in frequency between the peak of interest in the UV/visible absorption spectrum and the incident light used in the Raman experiment.
 Perform a TD calculation using a DFT method in order to determine the predicted location of the same peak.
 Specify a frequency for CPHF=RdFreq which is shifted from the predicted peak by the same amount as the incident light differs from the observed peak.
Preresonance Raman results are reported as additional rows within the normal frequency tables:
Harmonic frequencies (cm**1), IR intensities (KM/Mole), Raman
scattering activities (A**4/AMU), depolarization ratios for plane
and unpolarized incident light, reduced masses (AMU), force constants
(mDyne/A), and normal coordinates:
1
B1
Frequencies  1315.8011
Red. masses  1.3435
Frc consts  1.3704
IR Inten  7.6649
Raman Activ  0.0260
Depolar (P)  0.7500
Depolar (U)  0.8571
RamAct Fr= 1 0.0260 Additional output lines begin here.
DepP Fr= 1 0.7500
DepU Fr= 1 0.8571
RamAct Fr= 2 0.0023
DepP Fr= 2 0.7500
DepU Fr= 2 0.8571
VibrationRotation Coupling Output. If the VibRot option is specified, then the harmonic vibrationalrotational analysis appears immediately after the normal thermochemistry analysis in the output, introduced by this header:
Harmonic VibroRotational Analysis
If anharmonic analysis is requested as well (i.e., VibRot and Anharmonic are both specified), then the anharmonic vibrationalrotational analysis results follow the harmonic ones, introduced by the following header:
Secondorder Perturbative Anharmonic Analysis
Anharmonic Frequency Calculations. Freq=Anharmonic jobs produce additional output following the normal frequency output. (It follows the vibrationalrotational coupling output if this was specified as well.) We will briefly consider the most important items.
The output displays the equilibrium geometry (i.e., the minimum on the potential energy surface), followed by the anharmonic vibrationally averaged structure at 0 K:
Internal coordinates for the Equilibrium structure (Se)
Interatomic distances:
1 2 3 4
1 C 0.000000
2 O 1.206908 0.000000
3 H 1.083243 2.008999 0.000000
4 H 1.083243 2.008999 1.826598 0.000000
Interatomic angles:
O2C1H3=122.5294 O2C1H4=122.5294 H3C1H4=114.9412
O2H3H4= 62.9605
Dihedral angles:
H4C1H3O2= 180.
Internal coordinates for the vibrationally average structure at 0K (Sz)
Interatomic distances:
1 2 3 4
1 C 0.000000
2 O 1.210431 0.000000
3 H 1.097064 2.024452 0.000000
4 H 1.097064 2.024452 1.849067 0.000000
Interatomic angles:
O2C1H3=122.57 O2C1H4=122.57 H3C1H4=114.8601
O2H4H3= 62.8267
Dihedral angles:
H4C1H3O2= 180.
Note that the bond lengths are slightly longer in the latter structure. The predicted coordinates at STP follow in the output.
The anharmonic zero point energy is given shortly thereafter in the output:
Anharmonic Zero Point Energy

Harmonic : cm1 = 5008.40626 ; Kcal/mol = 14.320 ; KJ/mol = 59.914
Anharmonic Pot.: cm1 = 53.31902 ; Kcal/mol = 0.152 ; KJ/mol = 0.638
Watson+Coriolis: cm1 = 12.83227 ; Kcal/mol = 0.037 ; KJ/mol = 0.154
Total Anharm : cm1 = 4942.25496 ; Kcal/mol = 14.131 ; KJ/mol = 59.122
The anharmonic frequencies themselves appear just a bit later in this table, in the column labeled E(anharm):
==================================================
Anharmonic Infrared Spectroscopy
==================================================
Units: Transition energies (E) in cm^1
Integrated intensity (I) in km.mol^1
Fundamental Bands

Mode(n) E(harm) E(anharm) I(harm) I(anharm)
1(1) 2938.531 2788.983 55.17567187 55.41312200
2(1) 1888.862 1864.231 101.42877427 104.63741421
...
Overtones

Mode(n) E(harm) E(anharm) I(anharm)
1(2) 5877.061 5517.149 0.00211652
2(2) 3777.724 3710.383 3.68324904
...
Combination Bands

Mode(n) Mode(n) E(harm) E(anharm) I(anharm)
2(1) 1(1) 4827.393 4654.114 1.74785224
3(1) 1(1) 4490.139 4271.343 0.04557003
...
The harmonic frequencies are also listed for convenience.
Vibronic Analysis. The following input file predicts the vibronic spectrum:
%OldChk=excited Excited state calculation.
%Chk=fcht
# Freq=(ReadFC,FCHT,ReadFCHT) Geom=AllCheck …
TimeIndependent ReadFCHT additional input.
Output=Matrix=JK Output Duschinsky matrix and shift vector.
final blank line
The molecule specification is taken from the checkpoint file from the excited state, as are the force constants for the excited states.
FCHT analysis produces many results. The final Duschinsky (state overlap) matrix appears as follows:
Final Duschinsky matrix

Note: The normal coordinates of the final state (columns) are expressed
in the basis set of the normal coordinates of the initial state (rows)
1 2 3 4 5
1 0.539484D+00 0.839747D+00 0.139916D01 0.147815D01 0.167387D02
2 0.594185D+00 0.373849D+00 0.647845D+00 0.757424D01 0.627709D02
3 0.303582D01 0.276954D01 0.572527D02 0.354162D+00 0.933518D+00
…
Note that this output reports the value of J_{ij} for each pair of states. Generally, what is plotted is J^{2}.
The locations and intensities of the predicted bands are reported as follows:
==================================================
Information on Transitions
==================================================
Energy of the 00 transition: 31327.1976 cm^(1)
NOTE: The energy (transition energy) refers to the relative energy,
with respect to the 00 transition energy.
The intensity is the line intensity.
DipStr is the dipole strength.
Energy = 0.0000 cm^1: 0> > 0> Frequency and transition (states).
> Intensity = 7003. (DipStr = 0.9135E01)
Energy = 457.9310 cm^1: 0> > 9^1> Location is ~31875 cm^{1}.
> Intensity = 650.2 (DipStr = 0.8360E02) Intensity is in dm^{3}cm^{1}mol^{1}; dipole strength in au.
…
The final predicted spectrum follows in a form suitable for plotting:
==================================================
Final Spectrum
==================================================
Band broadening simulated by mean of Gaussian functions with
HalfWidths at HalfMaximum of 135.00 cm^(1)
Legend:

1st col.: Energy (in cm^1)
2nd col.: Intensity at T=0K
Intensity: Molar absorption coefficient (in dm^3.mol^1.cm^1)

30327.1976 0.000000D+00
…
31319.1976 0.699549D+04
31327.1976 0.701428D+04
31335.1976 0.699927D+04
…
Resonance Raman Spectra. The following input file computes the resonance Raman intensities from two previously run frequency calculations.
%Chk=S0_freq Ground state checkpoint file.
# Freq=(FC,ReadFC,ReadFCHT) Geom=AllCheck …
TimeIndependent ReadFCHT additional input.
Spectroscopy=ResonanceRaman Predict resonance Raman spectrum.
Spectrum=(Lower=800.,Upper=2800.,Broadening=Stick) Spectrum specifications.
Intermediate=Source=Chk Get second state data from checkpoint file (named below).
RR=(OmegaMin=55000,OmegaMax=56000,OmegaStep=100) RR analysis parameters: ω range and step size.
S2_freq.chk Excited state checkpoint file.
See the section on Freq=ReadFCHT for details about the additional input. For each of the Raman modes, the following output appears for each point in the specified range of incident energies (omega):
==================================================
Information on Transitions
==================================================
Energy of the 00 transition: 54854.2397 cm^(1)
Alp2: alpha^2, BsAl: beta_s(alpha)^2, BaAl: beta_a(alpha)^2
Energy = 0.0000 cm^1: 0> > 0> Relative energy and involved states.
> Omega = 55000.0 cm^1, Sigma = 1.1332
Alp2 = 0.33009E+02, BsAl = 0.29859E+03, BaAl = 0.00000E+00
Following this output, the same data is presented in a tabular form:
==================================================
Final Spectrum
==================================================
No band broadening applied (stick spectrum)
Legend:

1st col.: Raman shift (in cm^1)
2nd col.: Intensity at T=0K for incident energy: 55000.00 cm^1
3rd col.: Intensity at T=0K for incident energy: 55100.00 cm^1
4th col.: Intensity at T=0K for incident energy: 55200.00 cm^1
5th col.: Intensity at T=0K for incident energy: 55300.00 cm^1
Raman scattering intensity in cm^3.mol^1.sr^1

…
1188.0000 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00
1190.0000 0.134622D21 0.213038D21 0.358179D21 0.644832D21
1192.0000 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00
…
Since no spectral broadening was requested here (Spectrum=Broadening=Stick, the only rows with nonzero intensities correspond to the Raman active frequencies.
Examining LowLying Frequencies. The output from the full force constant matrix diagonalization (the default Freq=DiagFull), in which the rotational and translational degrees of freedom are retained, appears as following in the output:
Low frequencies  19.9673 0.0011 0.0010 0.0010 14.2959
Low frequencies  25.6133 385.4672 988.9028 1083.0692
This output is from an Opt Freq calculation on methanol. Ignoring sign, there are 3 lowlying modes, located at around 14, 19, and 25 wavenumbers (in addition to the three that are ~0). However, if we rerun the calculation using tight optimization criteria (Opt=Tight) and a larger integration grid, the lowest modes become:
Low frequencies  7.4956 5.4813 2.6908 0.0003 0.0007
Low frequencies  0.0011 380.1699 988.1436 1081.9083
The lowlying modes are now quite small, and the lowest frequencies have moved slightly as a result.
This analysis is especially important for molecular systems having frequencies at small wavenumbers. For example, if the lowest reported frequency is around 30 and there is a lowlying mode around 25 as above, then the former value is in considerable doubt (as is whether the molecular structure is even a minimum).
Rerunning a Frequency Calculation with Different Thermochemistry Parameters. The following twostep job contains an initial frequency calculation followed by a second thermochemistry analysis using a different temperature, pressure, and selection of isotopes:
%Chk=freq
# B3LYP/6311+G(2d,p) Freq
Frequencies at STP
molecule specification
Link1
%Chk=freq
%NoSave
# B3LYP/6311+G(2d,p) Freq(ReadIso,ReadFC) Geom=Check
Repeat at 300 K
0,1
300.0 1.0
16
2
3
…
Note also that the freqchk utility may be used to rerun the thermochemical analysis from the frequency data stored in a Gaussian checkpoint file.
The following keywords for specifying various aspects of Freq=Anharm calculations are included as additional input within the Gaussian input file. They control various aspects of anharmonic frequency analyses. Note that these keywords are completely different from those supported in Gaussian 09 (a few of these changes were introduced in Gaussian 09 revision D.01).
Data Sources and Format
The DataSrc, DataAdd and DataMod input items locate the various data required by the anharmonic frequency analysis. They each take a list of parameters and associated values which specify locations from which to retrieve different data items. In general, parameters specify what data is to be read and their values specify the location of that data. The available options for the latter are listed below. Generally, they may be optionally followed by a format suffix.
Source keywords are used to specify where the data is located:
 Src: Use data from the RWF file for the current job. RWF is a synonym.
 Chk: Retrieve data from the current checkpoint file (as defined with %Chk or %OldChk).
 In: Read data from the input stream.
 InChkn: Retrieve data from the nth file in the file list (see the discussion of additional input sections below). Valid values of n run from 1 to 6.
Format Suffixes are appended directly to the source item, and they specify a nondefault format for various readin data. For example, InQMW says to read derivative data from the input stream in massweighted normal coordinates. The following suffixes are available:
 QMW: Derivatives are with respect to normal modes in massweighted normal coordinates. Q is a synonym for this format suffix.
 QMWX: Harmonic derivatives are with respect to normal modes in Cartesian coordinates, and anharmonic derivatives are with respect to normal modes in massweighted normal coordinates. X is a synonym for this format suffix.
 QRedX: Harmonic derivatives are with respect to normal modes in Cartesian coordinates, and anharmonic derivatives are with respect to normal modes in dimensionless normal coordinates.
By default, QMWX is tried first, followed by QMW.
DataSrc=param: Specify the source(s) of various readin data. The parameter consists of a keyword indicating the data to which it applies and a source keyword indicating its location (and possibly its format). The available parameters are:
 source: Sets the data source for all data.
 Harm=source: Sets the data source for harmonic data. The default is taken from the source file.
 Anharm=source: Sets the data source for harmonic data. The default is taken from the source file.
 Coriolis=source: Sets the data source for the Coriolis couplings. At present, the only supported items are Src and In, and format suffixes may not be used.
 NMOrder=ordering: Specifies the order normal modes are stored in the input source, selected from the following list. This item may be specified in addition to a source item.
 AscNoIrrep: Ascending order. Do not sort by irreducible representation. This is the default.
 Asc: Ascending order. Sort by irreducible representation if possible.
 Desc: Descending order. Sort by irreducible representation if possible.
 DescNoIrrep: Descending order. Do not sort by irreducible representation.
 Print: Use same order as for printing.
The following DataSrc items are deprecated, and are included only for backward similarity to Gaussian 09 (where they functioned as toplevel additional input items).
 InDerAU: Use data from the input stream in atomic units.
 InDerAJ: Use data from the input stream in attoJoules
 InDerRed: Use data from the input stream in reduced form. Reduced is an alternate name for this item.
 InDerGau: Use data from the input stream with the layout of the Gaussian output. InGauDer is an alternate name for this item.
DataAdd=params: Read alternate data to replace or complete the original data. Using this option will replace the already existing information in the original data with the data specified here.
 Freq: Replace harmonic frequencies with values given in the input stream (in cm^{1}). A data source may also be specified as a parameter: Freq=source, but format suffixes are not allowed.
 PESFull=sources: Read force constraints from specified specified source.
 PESHarm=sources: Read harmonic force constants from specified specified source.
 PESAnh=sources: Read anharmonic force constants from specified specified source.
 EDipFull=sources: Read the electric dipole from the specified specified source.
 EDipHarm=sources: Read the harmonic components of the electric dipole from the specified specified source.
 EDipAnh=sources: Read the anharmonic components of the electric dipole from the specified specified source.
 MDipFull=sources: Read the magnetic dipole from the specified source.
 MDipHarm=sources: Read the harmonic components of the magnetic dipole from the specified source.
 MDipAnh=sources: Read the anharmonic components of the magnetic dipole from the specified source.
 PolFull=sources: Read the polarizability tensor from the specified source.
 PolHarm=sources: Read the harmonic components of the polarizability tensor from the specified source.
 PolAnh=sources: Read the anharmonic components of the polarizability tensor from the specified source.
 MagFFull=sources: Read the magneticfield properties from the specified source.
 MagFHarm=sources: Read the harmonic components of the magneticfield properties from the specified source.
 MagFAnh=sources: Read the anharmonic components of the magneticfield properties from the specified source.
 FreqDepPFull=sources: Read the frequencydependent properties from the specified source.
 FreqDepPHarm=sources: Read the harmonic components of the frequencydependent properties from the specified source.
 FreqDepPAnh=sources: Read the anharmonic components of the frequencydependent properties from the specified source.
DataMod=params: Modify the data in various manners.
 ScHarm=value: Scales harmonic frequencies with a constant scaling factor (default is 1.0).
 NoCor: Discards Coriolis couplings in calculations. By default, all couplings will be retained.
 DerOrder=N: Selects the derivatives order to keep. E.g. DerOrder=123 discards all quartic force constants.
 DerIndex=N: Sets the maximum number of independent indexes for a derivative. E.g. DerIndex=2 keeps k_{iij} but discards k_{ijk}.
 SkipPT2=what: Selectively removes derivatives based on the parameter, whose possible values are listed below:
 No: Do not remove the data. This is the default option.
 Modes: Removes the derivatives with respect to any of the normal modes given in the input stream.
 Constants: Removes the derivatives based on the indexes given in the input stream. The input explicitly specifies the force constants (energy derivatives) to be removed. Each line specifies the involved normal modes, with the derivative order implied by the number of indexes. For example, to remove the third derivatives with respect to normal coordinates Q_{1} Q_{2} Q_{5}, the input line would be:
1 2 5
 OptModes: Modify derivatives based on additional instructions in the input stream (see input ordering section below).
Tolerances=data: Modify the tolerance threshold to include/discard derivative data.
 Gradient=value: Threshold for the energy first derivatives (default is 3.7074×10^{−3}).
 Hessian=value: Threshold for the energy second derivatives (default is 3.7074×10^{−5}).
 Cubic=value: Threshold for the energy third derivatives (default is 3.7074×10^{−5}).
 Quartic=value: Threshold for the energy fourth derivatives (default is 3.7074×10^{−5}).
 Coriolis=value: Threshold for the Coriolis couplings (default is 1.0×10^{−3}).
 Inertia=value: Threshold for the principal moments of inertia (default is 1.0×10^{−4}Å^{2}).
 Symm=value: Tolerance for anharmonic data with respect to symmetry rules (default is 2%).
Output Control
This section specifies the contents and destination of the calculation output.
Print=items: Include items in the output file. Available items are the following:
 InDataX: Include data compatible with DataSrc=InQMWX. The form Print=InDataX=Ext writes the data to the external file input_data.dat.
 InDataNM: Include data compatible with DataSrc=InQMW. The form Print=InDataNM=Ext writes the data to the external file input_data.dat
 YMatrix: Include the Y matrix (a variant of the χ matrix).
 Verbosity=n: Specify the verbosity level. The default is 0.
 ITop=rep: Selects the representation used for rotational spectroscopy. By default, it is defined automatically by Gaussian from the principal moments of inertia. Available representations are:
 Ir: Ir Representation: I_{z} < I_{x} < I_{y}
 IIr: IIr Representation: I_{y} < I_{z} < I_{x}
 IIIr: IIIr Representation: I_{x} < I_{y} < I_{z}
 Il: Il Representation: I_{z} < I_{y} < I_{x}
 IIl: IIl Representation: I_{x} < I_{z} < I_{y}
 IIIl: IIIl Representation: I_{y} < I_{x} < I_{z}
 ZAxisSymm=axis: Sets the Eckart axis to be used as Z for the definition of the reduced Hamiltonians for the vibrorotational analysis. Available choices are:
 X: Z collinear with X.
 Y: Z collinear with Y.
 Z: Z collinear with Z.
 NMOrder=ordering: Specifies the order in which normal modes are listed:
 Asc: Ascending order. Sort by irreducible representation if possible.
 Desc: Descending order. Sort by irreducible representation if possible. This is the default.
 AscNoIrrep: Ascending order. Do not sort by irreducible representation.
 DescNoIrrep: Descending order. Do not sort by irreducible representation.
 PT2VarEVec: Include the eigenvector matrix from the diagonalization of the variational matrix.
 PT2VarStates: Include the projection of the variational states on the deperturbed ones.
 PT2VarProj: Include the projection of the DVPT2 states on the new variational states.
 InDataAU: Write data compatible with DataSrc=InDerAU (deprecated).
 Polymode: Write data to use in the Polymode program.
ReducedDimensionality Schemes
RedDim=items: Specifies which normal modes are active in the analysis. Items are:
 Active=n: Activate the n modes specified in the input stream. By default, all modes are active.
 Inactive=n: Read list of n inactive modes from the input stream.
 Frozen=n: Read n modes to be frozen from the input stream.
 MinFreqAc=freq: Sets the normal modes with a frequency above the specified value to be active (default is 0). Only valid if MaxFreqAc>MinFreqAc.
 MaxFreqAc=freq: Sets the normal modes with a frequency below the specified value to be active (default is infinity). Only valid if MaxFreqAc>MinFreqAc.
 MinFreqIn=freq: Sets the normal modes with a frequency above the specified value to be inactive (default is 0). Only valid if MaxFreqIn>MinFreqIn.
 MaxFreqIn=freq: Sets the normal modes with a frequency below the specified value to be active (default is infinity). Only valid if MaxFreqIn>MinFreqIn.
 MinFreqFr=freq: Sets the normal modes with a frequency above the specified value to be frozen (default is 0). Only valid if MaxFreqFr>MinFreqFr.
 MaxFreqFr=freq: Sets the normal modes with a frequency below the specified value to be frozen (default is infinity).Only valid if MaxFreqFr>MinFreqFr.
SecondOrder Vibrational Perturbation Theory (VPT2)
PT2Model=data: Sets the VPT2 model to use.
 HDCPT2: Use the Hybrid DegeneracyCorrected 2ndorder Perturbation Theory. This is the default.
 VPT2: Use the original 2ndorder Vibrational Perturbation Theory. Vibrational spectroscopy intensities are available for this model.
 DVPT2: Use the Deperturbed 2ndorder Vibrational Perturbation Theory. Vibrational spectroscopy intensities are available for this model. The form DVPT2=all selects all possibly resonant terms as Fermi resonances, and it is equivalent to Resonances=(DFreqFrm=∞,DPT2Var=0.)
 GVPT2: Use the Generalized 2ndorder Vibrational Perturbation Theory. Equivalent to DVPT2, but the removed terms are treated variationally in a second step. Vibrational spectroscopy intensities are available for this model. The form GVPT2=all selects all possibly resonant terms as Fermi resonances, and it is equivalent to Resonances=(DFreq12=∞,K12Min=0.)
 DCPT2: Use the DegeneracyCorrected 2ndorder Perturbation Theory.
HDCPT2=params: Set the parameters for the model with the Alpha and Beta options (i.e., HDCPT2=Alpha=value), which specify values for the corresponding variables in the expression for Λ:
Resonances=params: Set resonance thresholds and parameters for DVPT2 (Fermirelated items only) and GVPT2 calculations.
Spectroscopy
Spectro=MaxQuanta=quanta: Compute transition integrals to states with up to the specified quanta. The default is 2.
ROA=params: Options related to Raman optical activity. If the keyword is specified, then only those scatterings explicitly requested will be computed. If the ROA is not specified, then intensities are computed for all supported scatterings.
 ICP0: Compute ROA intensity for ICP forward scattering.
 ICP90x: Compute ROA intensity for incident circular polarization (ICP) rightangle scattering (polarized).
 ICP90z: Compute ROA intensity for ICP rightangle scattering (depolarized).
 ICP90*: Compute ROA intensity for ICP rightangle scattering (magic angle).
 ICP180: Compute ROA intensity for ICP backward scattering.
 SCP0: Compute ROA intensity for scattered circular polarization(SCP) forward scattering.
 SCP90x: Compute ROA intensity for SCP rightangle scattering (polarized).
 SCP90z: Compute ROA intensity for SCP rightangle scattering (depolarized).
 SCP90*: Compute ROA intensity for SCP rightangle scattering (magic angle).
 SCP180: Compute ROA intensity for SCP backward scattering.
 DCP180: Compute ROA intensity for double circular polarization (DCP) backward scattering.
 All: Compute the ROA intensity for all supported scatterings.
Freq=ReadAnharm Additional Input Ordering
The potential input sections for the various Freq=ReadAnharm additional input items should follow the keyword list section in the following order. Blank lines separate input sections, but each section and its terminating blank line should be included only when the corresponding keyword is specified.
Freq=ReadAnharm input keywords
blank line
checkpoint file list for DataSrc=InChkn or DataAdd=…=InChkn
blank line
data for DataSrc=In or DataAdd=…=In (harmonic followed immediately by anharmonic)
blank line
data for DataAdd=Freq
blank line
n modes for RedDim=Active
blank line
n modes for RedDim=Inactive
blank line
n modes for RedDim=Frozen
blank line
keyword for DataMod=SkipPT2=OptModes (see below)
modes for DataMod=SkipPT2=Modes or Constants or DataMod=SkipPT2=OptModes
blank line
data for Resonances=List
blank line
DataMod=SkipPT2=OptModes Keywords: Removal of Contributions from Selected Normal Modes
One or more options are specified on a single line with no blank line following. The options available are:
 MinInd=n: Controls the minimum number of times a selected normal mode must appear to discard the derivative. The default is 1. For example, a value of 2 means that k_{ijk} is kept, but k_{iij} and k_{iii} are removed.
 EnOrd=mask: Controls the energy derivative orders to consider. The default is 1234 which says to include all energy derivatives. The value 14 says to only treat the first and fourth energy derivatives and to ignore the second and third energy derivatives.
Here is an example calculation using DataMod=SkipPT2=OptModes:
# Freq=(Anharm,ReadAnharm) …
Formaldehyde
0 1
C 0.6067825443565 0.0000000216230 0.0000000000000
O 0.6033290944914 0.0000000215000 0.0000000000000
H 1.1752074085613 0.9201232113261 0.0000000000000
H 1.1752073429832 0.9201232950844 0.0000000000000
PT2Model=GVPT2 Resonances=(NoDarDen,NoFermi)
DataMod=SkipPT2=OptModes
MinInd=2 EnOrd=34 Control keywords: keep only k_{ijk} for third and fourth derivatives
5 6 Modes for which to remove derivatives
terminating blank line
Specifying Resonances
For Resonances=Add, Resonances=List=Replace and Resonances=Delete, each line contains the type of resonance and the indexes of the normal modes involved in the resonances. Supported types are:
 12 or 12: 12 Fermi resonance. 3 indexes needed, given in the order: ω_{1} ≈ ω_{2} + ω_{3}
 11 or 11: 11 DarlingDennison resonance. 2 indexes needed, given in the order: ω_{1} ≈ ω_{2}
 13 or 13: 13 DarlingDennison resonance. 4 indexes needed, given in the order: ω_{1} ≈ ω_{2} + ω_{3} + ω_{4}
 22 or 22: 22 DarlingDennison resonance. 4 indexes needed, given in the order: ω_{1} + ω_{2} ≈ ω_{3} + ω_{4}
For Resonances=List=Modify, an action must be specified at the beginning of the line, before the resonance type and the normal modes. The available actions are ADD (add resonance) and DEL (remove resonance if previously identified).
Examples
Example 1: Frequency data calculated in the current job:
%Chk=example1
# B3LYP/6311+G(d,p) Freq=(Anharmonic,ReadAnharm)
Anharmonic frequencies example 1
molecule specification
DataMod=SkipPT2=Modes Freq=ReadAnharm additional input
RedDim=Inactive=1 # modes to make inactive
PT2Model=GVPT2
Resonances=Add
Print=InDataX=Ext Write data to an external file
blank line
6 Modes to make inactive: RedDim=Inactive
blank line
4 5 Modes for which to discard derivatives: DataMod=SkipPT2=Modes
blank line
12 2 3 3 List of additional resonances: Resonances=Add
blank line
Example 2: Frequency data read from input stream. It uses the checkpoint file from Example 1 to retrieve the molecular geometry and Hessian:
%OldChk=example1
%Chk=example2
# B3LYP/6311+G(d,p) Freq=(ReadFC,Anharmonic,ReadAnharm) Geom=Check
Anharmonic frequencies example 2
0 1
PT2Model=GVPT2 Freq=ReadAnharm additional input
DataSrc=(InQMWX,NMOrder=Desc) Read harmonic/anharmonic data from input
DataAdd=Freq Replace harmonic frequencies with input values
RedDim=Inactive=1
blank line
@input_data.dat Data file from previous job with Print=InDataX=Ext
blank line
1198.351 List of harmonic frequencies: DataAdd=Freq
1259.285
1530.636
1814.590
2884.164
2941.556
blank line
6
blank line
Example 3: The following example specifies alternate values for several parameters:
# B3LYP/631G(d) Freq=(Anharm,ReadAnharm)
Anharmonic frequencies example 3
molecule specification
Tolerances=Coriolis=0.25 Freq=ReadAnharm additional input
Resonances=(DFreq12=220.0,K12Min=1.1)
DataMod=SCHarm=0.98
blank line
Example 4: Anharmonic VCD and ROA spectra calculation:
%Chk=project4
# Freq=(ROA,VCD,Anharm,ReadAnharm) CPHF=RdFreq …
Anharmonic VCD and ROA spectrum
0 1
molecule specification
589.3nm incident light frequency
PT2Model=GVPT2 ReadAnharm input section
Spectro=MaxQuanta=3
blank line
The following keywords for specifying various aspects of Freq=FC, HT or FCHT calculations are included as additional input within the Gaussian input file. They control various aspects of these analyses. Note that these keywords are different from those supported in Gaussian 09.
Specifying Calculation Options and Parameters
Transition=type: Specify the type of electronic transition.
 Absorption: Absorption. This is the default.
 Emission: Emission.
Spectroscopy=type(s): Spectroscopy to simulate. The default is Spectroscopy=OnePhoton. The following items are available (names may be abbreviated to just the capital letters below: e.g., RR for ResonanceRaman):
 OnePhoton: Onephoton process. This is the default.
 CircularDichroism: Circular Dichroism: electronic circular dichroism (ECD) for absorption and circularly polarized luminescence (CPL) for emission.
 OnePhotonAbsorption: Onephoton absorption. Implies Transition=Absorption.
 OnePhotonEmission: Onephoton emission. Implies Transition=Emission.
 ElectronicCircularDichroism: Electronic circular dichroism. Implies Transition=Absorption.
 CircularlyPolarizedLuminescence: Circularly polarized luminescence. Implies Transition=Emission.
 ResonanceRaman: Vibrational Resonance Raman. Options are specified via a separate item (see below).
ResonanceRaman=options: The following RRspecific options are available:
 CombOnly: Compute only combination bands.
 NoComb: Skip computation of combination bands.
 aMean=coeff: Set the coefficient of the mean polarizability. The default is 45.
 Damping=value: Damping constant or lifetime of intermediate states (in cm^{−1}). The default is 500 cm^{1}.
 dAnti=coeff: Set the coefficient of the antisymmetric anisotropy. The default is 5.
 gSymm=coeff: Set the coefficient of the symmetric anisotropy. The default is 7.
 Omega=value: Incident energy (in cm^{−1}). By default, it is calculated as the difference in energy between the vibrational fundamental states of the initial and intermediate states.
 OmegaList: Reads a list of incident energies (in cm^{−1}) in the input stream to build an energy profile.
 OmegaMax=value: Maximum energy (in cm^{−1}) for the incident energy profile.
 OmegaMin=value: Minimum energy (in cm^{−1}) for the incident energy profile.
 OmegaNum=n: Number of energies in the incident energy profile.
 OmegaStep=value: Energy interval (in cm^{−1}) between two energies for the incident energy profile.
 Scattering=params: Scattering geometry for the Resonance Raman simulation. Available items are (note that the asterisks below are part of the parameter names):
 ICP0: Compute RR intensity for Incident Circular Polarization (ICP) forward scattering.
 ICP90x: Compute RR intensity for ICP rightangle scattering (polarized).
 ICP90z: Compute RR intensity for ICP rightangle scattering (depolarized).
 ICP90*: Compute RR intensity for ICP rightangle scattering (magic angle).
 ICP180: Compute RR intensity for ICP backward scattering.
 SCP0: Compute RR intensity for Scattered Circular Polarization (SCP) forward scattering.
 SCP90x: Compute RR intensity for SCP rightangle scattering (polarized).
 SCP90z: Compute RR intensity for SCP rightangle scattering (depolarized).
 SCP90*: Compute RR intensity for SCP rightangle scattering (magic angle).
 SCP180: Compute RR intensity for SCP backward scattering.
 DCP180: Compute RR intensity for Double Circular Polarization (DCP) forward scattering.
 All: Compute the RR intensity for all supported scatterings.
Temperature: Specify temperature and related parameters:
 Value=temp: Temperature of simulation in K. The default is the value specified to Gaussian with the Temp keyword or via additional input; the Gaussian default is 298.15 K.
 MinPop=ratio: Minimum ratio between the Boltzmann populations of any vibrational initial state and the fundamental state for the former to be considered in the calculations. In other words, this item specifies the fraction of a vibrational state that must be populated in order for it to be treated as the starting point of a transition. The default value is 0.1 (10%).
TransProp=definition: Definition of the transition dipole moment(s). The first three items effectively select among FranckCondon, HerzbergTeller and FCHT analyses. However, the corresponding options to the Freq keyword are preferable. The following items are available:
 FC: The dipole moment is assumed constant during the electronic transition. This selection describes dipoleallowed transitions well. It is the default.
 FCHT: Computes zeroth and firstorder terms of the Taylor expansion of the transition dipole moment about the equilibrium geometry. It is needed to correctly treat weaklyallowed electronic transitions or CD spectra.
 HT: Implements linear variation of the dipole moment with respect to the normal mode. This item corresponds to the firstorder term of the Taylor expansion of the transition dipole moment about the equilibrium geometry.
 DipA=source: Explicit definition of the transition dipole moment d_{A}. The following options for the data source are available:
 Auto: Gaussian will choose the definition depending on the simulation parameters. This is the default.
 Read: Read data from main input source (see “Data Sources” below).
 Input: Read data from input stream.
 DipB=source: Definition of the transition dipole moment d_{B}. This item accepts the same data source options as DipA.
 EDip=source: Definition of the transition electric dipole moment. This item accepts the same data source options as DipA.
 MDip=source: Definition of the transition magnetic dipole moment. This item accepts the same data source options as DipA.
 NoUse: Discard the transition dipole read in the data source and replace it with a unitary one. This is the default behavior for electronic transition with a change in multiplicity or charge. Keyword is not supported for HT and FCHT calculations.
Method=method: Selects the representation model for the electronic transition. The nondefault options specify methods for approximating the excited state frequencies based on the ground state. These are not generally preferable to modeling the excited state explicitly. The following methods are available:
 AdiabaticHessian: Both PESs—ground state and excited state—are calculated at the harmonic level about their respective minimum. This is the default.
 AdiabaticShift: Both PESs are calculated at the harmonic level about their respective minimum, but the PES of the final state is assumed to be the same as the initial state. Only the equilibrium geometry of the final state is calculated.
 VerticalHessian: The PES of the final state is evaluated about the equilibrium geometry of the initial state.
 VerticalGradient: The PES of the final state is evaluated about the equilibrium geometry of the initial state, but the PES of the final state is assumed to be the same as the initial state. Only the energy gradient of the final state is calculated about the equilibrium geometry of the initial state. (This method is also known as the linear coupling model.)
For emission spectra, the form VerticalGradient=Abs causes Gaussian to compute the 00 transition energy in the same way as absorption spectra in order to get the correct 00 transition energy. It is needed when the emission spectrum is incorrectly computed as an absorption spectrum (equilibrium geometry of the ground state, frequencies of the ground state and forces of the excited state).
Prescreening=params: Sets the prescreening criteria for choosing the most intense transitions. Only available for the timeindependent framework. Available parameters:
 MaxC1=n: Maximum quantum number reached in C_{1} (C_{1}^{max}). The default is 20.
 MaxC2=n: Maximum quantum number reached by both modes involved in the combination bands in C_{2} (C_{2}^{max}). The default is 13.
 MaxInt=millions: Maximum number of integrals to compute for each class above C_{2} (N_{I}^{max}), in units of one million. The default value is 100 (100,000,000).
TimeIndependent: Use the timeindependent framework. This is the default for onephoton spectroscopy. TI is a synonym for this item.
TimeDependent=params: Compute the band shape using the pathintegral approach instead of the sumoverstates approach. This is the default for Resonance Raman spectroscopy. TD is a synonym for this item. Available parameters are:
 2NStep=m: Use 2^{m} steps for the integration. m defaults to 18 for OPA, OPE, ECD and CPL and 12 for RR.
 2NStepWin=k: Set the number of steps in which the correlation function Χ(t) is actually computed to as 2^{k}. The function is 0 outside this range. Note that k must be ≤ m from 2NStep above. By default, k is set equal to m.
 GauHWHM=n: Inhomogenous broadening, applied as a dephasing (in cm^{−1}). The default is 135 cm^{−1}.
 LorHWHM=n: Homogenous broadening, applied as a dephasing (in cm^{−1}). The default is 0 cm^{−1}.
 Time=seconds: Time interval Δt in seconds. The default value is 2^{m}×10^{17}, where m is from 2NStep above.
Termination=DeltaSP=value: Sets the termination criteria. Set the minimum difference between two consecutive classes of the final state to continue the calculation. The default is 0.0 (always continue).
Data Sources
The items in this section specify the locations and methods for obtaining various data used by the FCHT analysis.
Initial=items: Data source(s) and/or parameters for the initial state.
Final=items: Data source(s) and/or parameters for the final state.
Intermediate=items: Data sources(s) and/or parameters related to the intermediate state. Only valid for Resonance Raman.
The following parameters are available for the three preceding items:
 Source=source: Data source for the state:
 Calc: Current calculation. This is the default for Absorption spectra.
 Chk: Checkpoint file (the filename is given in the input stream). This is the default for Emission spectra.
 LogFile: Gaussian output filename (the filename is given in the input stream).
 Freq=params: Source and handling of vibrational frequency data:
 Read: Read frequencies from the main data source. This is the default.
 Input: Read frequencies from the input stream.
 Scale: Scale frequencies using a modespecific extrapolated scaling factor based on the Duschinsky transformation. The reference state used is the other electronic state. For example, Final=Freq=Scale uses the input initial state frequencies to scale those of the final state. The required frequencies are taken from the input stream in their usual position (see below).
 MaxBands=state: Set the highest class state to consider. The default is 3 for Initial and 7 otherwise.
 MaxStates=n: Maximum number of initial vibrational states actually considered in the calculations. Only valid with Initial. Note that this value is not the number of configurations printed by Gaussian. The default is 50.
 ExcState=state: Excited electronic state actually treated in the Gaussian output file. Only used if the data source for the excited state is a Gaussian output file.
DataAdd=DeltaE=value: Difference in energy between the electronic states (in Hartrees). By default, it is calculated from the data sources used for the initial and final states.
DataMod=Duschinsky=params: Duschinsky matrix to use in the calculation, By default, the true Duschinsky matrix is used. Note that the definition of the Duschinsky matrix depends on the model used to describe the transition. Other options are:
 Identity: Use the identity matrix as the Duschinsky matrix.
 Diagonal: Swap the columns of the correct Duschinsky matrix to be as diagonal as possible and replace it by the identity matrix. If this cannot be done because the matrix is not diagonal enough, an error occurs.
Superpose=params: Specifty the superposition procedure.
 Reorient: Reorient molecules to maximize overlap between the two structures. This is the default.
 NoReorient: Do not try to superpose the structures.
 Rotation=n: Rotation algorithm used for the superposition:
 0: Alternate usage. First quaternions, then rotation angles. This is the default.
 1: Use the program default (be aware that this may change between Gaussian versions).
 2: Use quaternions only.
 3: Use rotation angles only.
 RotNIter=xxxyyy: Maximum number of iterations for each superposition algorithm, where xxx is the number of iterations for the algorithm based on quaternions, and yyy is the number of iterations for the anglebased algorithm. The default parameter is 030100 or 30 for quaternions and 100 for rotation angles.
Output Selection
Print=params: This item controls the information included in the output.
 Spectra=params: Spectra to be included (the nondefault items are not valid in the timedependent framework nor for Resonance Raman):
 Final: Print only the final spectrum. This is the default.
 All: Print the maximum number of classspecific spectra.
 ClassI: Print one spectrum for each class of the initial state.
 ClassF: Print one spectrum for each class of the final state.
 Matrix=list: Matrices to be displayed: J, K, A, B, C, D, E. The string is composed of the matrices of interest, e.g., JK. A, B, C, D and E are the Sharp and Rosenstock matrices, J is the Duschinsky matrix, and K is the shift vector.
 HuangRhys: Print HuangRhys factors [Huang50].
 AssignThresh=n: Include data (assignment, energy, intensity) of the transitions, which contribute at least to 100n% to the total intensity of the spectrum obtained for each initial vibrational state. The default is 0.01 (1%).
 TDAutoCorr=n: Print n points of the timedependent autocorrelation function. The default is 0 (output disabled).
 Tensors: Print the transition tensors instead of the invariants. Only available for Resonance Raman.
 Conver: Print convergence data. Only available for Resonance Raman.
 Color=params: Specify RGB standard for color output:
 sRGB: Use the sRGB standard (white reference: D65, gamma correction based on the IEC 6196621 standard). This is the default.
 None: Do not print the color in an RGB format.
 CIE: Use the CIE RGB standard (white reference: D65, no gamma correction).
 HDTV: Use the HDTV RGB standard (white reference: D65, gamma correction based on the ITUR BT.7093 standard).
Spectrum=params: Control spectrum layout.
 Lower=value: Energy of the lower bound of the spectrum (in cm^{−1}). For absorption, the default is 1000 cm^{−1}. For emission, the default is 8000 cm^{−1}. The bound is defined with respect to the 00 transition.
 Upper=value: Energy of the upper bound of the spectrum (in cm^{−1}). For absorption, the default is +8000 cm^{−1}. For emission, the default is +1000 cm^{−1}. The bound is defined with respect to the 00 transition.
 AbsBound: Deactivate the default behavior of Gaussian to define the lower and upper bounds of the spectrum with respect to the 00 transition.
 Grain=value: Energy difference between two adjacent points for the discretization of the spectrum (in cm^{−1}). The default is 8cm^{1}.
 Broadening=params: Discribution function used to compute the bandshape:
 Gaussian: Use normal distribution functions to simulate the inhomogeneous broadening. This is the default.
 Lorentzian: Use Cauchy distribution functions to simulate the homogeneous broadening.
 Stick: Do not simulate the band broadening. Print the bands as sticks.
 HWHM=value: Halfwidth at halfmaximum of the distribution function used for the convolution (in cm^{−1}). The default is 135 cm^{1}.
ReducedDimensionality Schemes
RedDim=params: Activates and sets a reduceddimensionality scheme.
 Block: Use the correct Duschinsky matrix to find the projection of a list of modes given in the input stream on the other state to reduce consistently the dimension of the problem. By default, the state of reference is the lower one.
 ClearLowFreq=n: Remove all normal modes with a frequency below n cm^{−1} in absolute value. If this parameter is specified but n is omitted, it defaults to 50.
 ClearImFreq: Remove at most one normal mode per state with an imaginary frequency. If there is one mode with an imaginary frequency in each state, Gaussian checks that they are uncoupled to other modes and project themselves on each other in the electronic transition.
 BlockThresh=value: Threshold for the block(s) construction (between 0 and 1). A high value ensures that the selected modes are uncoupled from the remaining ones. The default is 0.9.
 BlockTol=n: Specifies the maximum ratio allowed between the final number of selected modes and the initial set: if more than n×#modes chosen by the user are selected, Gaussian stops the calculations. The default is 1.6.
Freq=ReadFCHT Additional Input Ordering
The potential input sections for the various Freq=ReadFCHT additional input items should follow the keyword list section in the following order. Blank lines separate input sections, but each section and its terminating blank line should be included only when the corresponding item is specified.
Freq=ReadFCHT input keywords
blank line
filename for the Chk or LogFile parameters to Input=Source
filename for the Chk or LogFile parameters to Final=Source (OPA, OPE, ECD, CPL) or Intermediate=Source (RR)
blank line
initial state frequency list: Initial=Freq=Input and/or Final=Freq=Scale
blank line
final/intermediate state frequency list: the Freq=Input option to Final (OPA, OPE, ECD, CPL) or Intermediate (RR) and/or Initial=Freq=Scale
blank line
transition dipole moment data for TransProp=DipA or TransProp=EDip
blank line
transition dipole moment data for TransProp=DipB or TransProp=MDip
blank line
list of modes for RedDim=Block
blank line
incident energies for ResonanceRaman=OmegaList
blank line
Examples
The following input file illustrates the use of additional input for FC analyses.
Example 1. The following calculation performs a FranckCondon frequency analysis for phenoxyl:
%Chk=phenoxyls0.chk
# Freq=(FC,ReadFC,ReadFCHT) Geom=Check …
Phenoxyl FrankCondon analysis
0 2
Final=Source=Chk ReadFCHT input: specify source for final state.
Print=(Spectra=All,Matrix=JK) Output all spectra, the Duschinsky matrix and the shift vector.
phenoxyls1.chk Checkpoint file for the final state.
terminal blank line
Example 2. The following job predicts the ECD spectrum (selected by the item in the additional input section). Note that the second checkpoint file for the final section is specified in the subsequent input section. It is read even though there is no explicit item in the additional input section since Final=Source=Chk is the default for emission spectra.
%Chk=initial
# Freq=(FCHT,ReadFC,ReadFCHT) Geom=AllCheck …
Transition=Absorption ReadFCHT input.
Spectroscopy=CD Select spectrum to predict: ECD.
final.chk Checkpoint file for the final state.
terminal blank line
Example 3. The following job performs the analysis at 500 K, using the timedependent framework.
%Chk=temp500init
# Freq=(FC,ReadFC,ReadFCHT,SaveNM) Geom=AllCheck …
Temperature=Value=500.0
temp500final.chk
terminal blank line
Example 4. The following job computes the resonance Raman spectrum:
%Chk=S0_freq
# Freq=(FC,ReadFC,ReadFCHT) Geom=Check …
RR spectrum
0 1
Spectroscopy=RR ReadFCHT input: select spectrum to predict.
TransProp=EDip=Input
RR=OmegaList
TD=(2NSTEP=12,2NSTEPWIN=12,Time=1.0d12,GauHWHM=100)
Print=(Tensors,Matrix=JK)
S2_freq.chk Checkpoint file for final state.
1.000D0 1.000D0 0.00D0 Transition dipole for TransProp=EDip.
55000 55500 57000 List of incident energies for RR=OmegaList.
terminal blank line
Example 5. The final job step below computes the vibrational envelope of a photoionization spectrum:
%chk=neutral Calculation on the neutral form (initial state).
# B3LYP/631+G(d,p) Freq=SaveNM
neutral form
0 1
molecule specification
Link1
%chk=cation Calculation on the cation (final state).
# B3LYP/631+G(d,p) Freq=SaveNM
cation
1 2
molecule specification
Link1
%oldchk=neutral Data for the neutral form.
%chk=fc
# B3LYP/631+G(d,p) Freq=(FC,ReadFCHT) Geom=Check
photoionization
0,1
Initial=Source=Calc Final=Source=Chk Retrieve final state from a checkpoint file.
Spectrum=(Lower=500.,Upper=+5000.,Grain=1.,HWHM=50.) Parameters for computed spectrum.
Prescreening=(MaxC1=30.,MaxC2=20.) Parameters to select transitions.
Print=Matrix=JK Output Duschinsky matrix and shift vector.
cation.chk Checkpoint filename for final state.
terminal blank line