Search:

# Manual For GUINEA

## Overview:

GUINEA evaluates vibrational energies and properties using explicitly summed Rayleigh-Schroedinger type perturbation theory. Vibrational energies, transition intensities, dipole moment corrections and harmonic derivatives may all be calculated at the VPT2 (second-order) and VPT4 (fourth-order) levels of perturbation theory. The numerical sums are written in such a way as to minimize the number of unnecessary intermediate states generated. This formulation is quite efficient, even for large numbers of modes (up to 24 modes have been tested and such calculations may be run in real-time). GUINEA may also explicitly form and diagonalize an effective Hamiltonian, followed by a transformation of other quantities such as transition moments, which may be useful for resonant systems. GUINEA uses for defining the necessary operators the files generated by the xcubic program of CFOUR, e.g. "quadratic", "cubic", "rota", "coriolis", etc. The files defining the elecronic potential must be present, while the rotational files and dipole moment files are optional, as are the main job archive files JOBARC/JAINDX. GUINEA will optionally use the additional "didq" file for third- and fourth-order Coriolis terms.

## Input:

GUINEA accepts user input either from a file, such as with the syntax "guinea < input_file", or interactively from the terminal. By default, interactive sessions prefix input lines with a prompt, ">", for clarity, while input from a file does not. Keywords are given one to a line, with any extra input for that keyword on one or more additional lines. Blank lines are ignored.

Keywords are grouped into two categories: "state commands" and "action commands". State commands add to, remove from, or otherwise modify the groups of states, called "sets", upon which the action commands act. A set may have any number of states, although only 100 sets may exist at any one time. Each set has its own settings for diagonalization, while the reference state and temperature (used in the Boltzmann factor for hot bands) are shared among all sets. Action commands then run the various calculations, or specify which calculations to run later, on the states in the selected set.

## Units:

• All energies are in wavenumber units (cm-1).
• Intensities are in km/mol.
• Dipole moment corrections are in milli-Debye (mD).
• First harmonic derivatives are unitless, and second harmonic derivatives are in reciprocal wavenumber units (cm).

## Memory Management:

Unlike ACES2, GUINEA uses its own memory management in which segments are allocated dynamically rather than from a common block. Total memory usage is limited by a user-defined limit (default is unlimited). The current version may have a few memory leaks and/or inconsistencies in the reported memory usage. However, these should be slight enough as to not affect normal usage (unless perhaps a huge number of states are shuffled in and out or the number of modes is very very large).

## Keywords:

This is a summary of all of the available keywords. If a keyword accepts additional input, the syntax will be given. Note that when the summary for a calculation is printed, states are listed in order of increasing energy, not the order in which they were added to the set, and that states are sorted by symmetry as they are added to a set. For diagonalized sets, the state printed for each level is the dominant contribution the that level's eigenvector.

### I. Global commands

print
Sets the print verbosity level: 1 = display < 2 = interactive < 3 = full; interactive is the default for terminal input, and display is the default for input from a file.
    print
[1,2,3]

OR

print
[display,interactive,full]

?print
Lists the current print verbosity level.
relative
Sets that printed energies should be relative to ground (i.e ZPVE is not included).
!relative
Sets that printed energies should be absolute (i.e ZPVE is included).
?relative
Lists whether printed energies are relative to ground.
freqformat
Sets the format for printed frequencies.
    freqformat
width.precision

intformat
Sets the format for printed intensities.
    intformat
width.precision

dipformat
Sets the format for printed dipole moment corrections.
    dipformat
width.precision

inttype
Sets the type of transition intensities to print.
    inttype
[total,cartesian,both]

diptype
Sets the type of dipole moment corrections to print.
    diptype
[total,cartesian,both]

memory
Sets the maximum amount of memory to use, initially unlimited. The amount must be integer for byte amounts but may be floating point for larger units.
    memory
n    [in bytes, 0=unlimited]

OR

memory
n[KB,MB,GB]

?memory
Lists the current memory limit and use.
polyad
Define a polyad from a set of states given as a list of integers, one for each normal mode. Defined polyads can be included in subsequent definitions of states as if they were additional normal modes (listed in the order defined). A polyad definition may include previously defined polyads. Polyads may be used in any set.
    polyad
n
p_1,1 p_1,2 ... p_1,m
p_2,1 p_2,2 ... p_2,m
...
p_n,1 p_n,2 ... p_n,m

!polyad
Remove the specified polyad by index.
    !polyad
i

?polyad
List all currently defined polyads and their order.
reference
Set the reference for hot bands.
    reference
r_1 r_2 ... r_m

?reference
List the reference for hot bands.
temp
Sets the temperature for hot bands.
?temp
Prints the temperature for hot bands.
adjfreq
    adjfreq
i
frequency

load
Runs all of the commands in the given file.
    load
filename

save
Write a file containing all of the commands necessary to reproduce the current program state (excluding the output of any calculations or queries).
    save
filename

done
Bye-bye! (also "quit", "exit", "bye")

### II. State commands

states
Add all of the specified states to the current set. States are given as a list of integers, one for each normal mode (and optionally, each polyad).
    states
n
s_1,1 s_1,2 ... s_1,m
s_2,1 s_2,2 ... s_2,m
...
s_n,1 s_n,2 ... s_n,m

!states
Remove all of the specified states from the current set.
    !states
n
s_1,1 s_1,2 ... s_1,m
s_2,1 s_2,2 ... s_2,m
...
s_n,1 s_n,2 ... s_n,m

?states
List all of the states in the current set.
!indices
Remove the states with the given indices from the current set.
    !indices
n
i_1
i_2
...
i_n

combination
Add all of the states from the given combination to the current set. This is all states with total number of quanta m in modes q_1 q_2 ... q_n. n=0 selects states with total number of quanta m in any mode
    combination
m n
q_1 q_2 ... q_n

OR

combination
m 0

!combination
Remove all of the states from the given combination from the current set. Same syntax as 'combination'.
    !combination
m n
q_1 q_2 ... q_n

OR

!combination
m 0

copyset
Add all states from the given set to the current set.
    copyset
i

fermi
Enables Fermi resonance terms in the effective Hamiltonian for the current set.
!fermi
Disables Fermi resonance terms in the effective Hamiltonian for the current set.
?fermi
Lists whether the current set includes Fermi resonance terms or not.
dd
Enables Darling-Dennison resonance terms in the effective Hamiltonian for the current set.
!dd
Disables Darling-Dennison resonance terms in the effective Hamiltonian for the current set.
?dd
Lists whether the current set includes Darling-Dennison resonance terms or not.
diagonalize
Enables diagonalization of the effective Hamiltonian for the current.
!diagonalize
Disables diagonalization of the effective Hamiltonian for the current.
?diagonalize
List whether the current set is diagonalized.
set
Creates a new set.
setset
Sets the current set.
    setset
i

?set
Lists basic information about all of the sets (does not list states).
!set
Removes the given set, 0 for current set.
    !set
i

copyset
Copies all states from the specified set to the current set.
    copyset
i


### III. Action commands

vibration
Set the calculation level for vibrational energies (default vpt2).
    vibration
[harmonic,vpt2,vpt4]

?vibration
List the current calculation level for vibrational energies.
intensity
Set the calculation level for intensities (default off).
    intensity
[off,harmonic,vpt2]

?intensity
List the current calculation level for intensities.
dipole
Set the calculation level for dipole moment corrections (default off).
    dipole
[off,harmonic,vpt2]

?dipole
List the current calculation level for dipole moment corrections.
harmderiv
Runs a harmonic derivative calculation.
calc
Runs the specified calculations on the current set.
?calc
Lists which calculations are turned on.

## Examples

Example 1 Input Example 1 Output
(this example uses a quartic force field of water)