4.3 Basic SCF Job Control

4.3.4 Symmetry

Symmetry is a powerful branch of mathematics and is often exploited in quantum chemistry, both to reduce the computational workload and to classify the final results obtained.908, 242, 243 Q-Chem is able to determine the point group symmetry of the molecular nuclei and, on completion of the SCF procedure, classify the symmetry of molecular orbitals, and provide symmetry decomposition of kinetic and nuclear attraction energy (see Chapter 11).

Molecular systems possessing point group symmetry offer the possibility of large savings of computational time, by avoiding calculations of integrals which are equivalent i.e., those integrals which can be mapped on to one another under one of the symmetry operations of the molecular point group. The Q-Chem default is to use symmetry to reduce computational time, when possible.

There are several keywords that are related to symmetry, which causes frequent confusion. SYM_IGNORE controls symmetry throughout all modules. The default is FALSE. In some cases it may be desirable to turn off symmetry altogether, for example if you do not want Q-Chem to reorient the molecule into the standard nuclear orientation, or if you want to turn it off for finite difference calculations. If the SYM_IGNORE keyword is set to TRUE then the coordinates will not be altered from the input, and the point group will be set to C1.

The SYMMETRY keyword controls symmetry in some integral routines. It is set to FALSE by default. Note that setting it to FALSE does not turn point group symmetry off, and does not disable symmetry in the coupled-cluster suite (CCMAN and CCMAN2), which is controlled by CC_SYMMETRY (see Chapters 6 and 7), although we noticed that sometimes it may interfere with the determination of orbital symmetries, possibly due to numerical noise. In some cases, SYMMETRY = TRUE can cause problems (poor convergence and wildly incorrect SCF energies) and turning it off can avoid these problems.

Note:  The user should be aware about different conventions for defining symmetry elements. The arbitrariness affects, for example, C2v point group. The specific choice affects how the irreducible representations in the affected groups are labeled. For example, b1 and b2 irreducible representations in C2v are flipped when using different conventions. Q-Chem uses non-Mulliken symmetry convention. See http://iopenshell.usc.edu/howto/symmetry for detailed explanations.

SYMMETRY
       Controls the efficiency through the use of point group symmetry for calculating integrals.
TYPE:
       LOGICAL
DEFAULT:
       TRUE Use symmetry for computing integrals.
OPTIONS:
       TRUE Use symmetry when available. FALSE Do not use symmetry. This is always the case for RIMP2 jobs
RECOMMENDATION:
       Use the default unless benchmarking. Note that symmetry usage is disabled for RIMP2, FFT, and QM/MM jobs.

SYM_IGNORE
       Controls whether or not Q-Chem determines the point group of the molecule and reorients the molecule to the standard orientation.
TYPE:
       LOGICAL
DEFAULT:
       FALSE Do determine the point group (disabled for RIMP2 jobs).
OPTIONS:
       TRUE/FALSE
RECOMMENDATION:
       Use the default unless you do not want the molecule to be reoriented. Note that symmetry usage is disabled for RIMP2 jobs.

SYM_TOL
       Controls the tolerance for determining point group symmetry. Differences in atom locations less than 10-SYM_TOL are treated as zero.
TYPE:
       INTEGER
DEFAULT:
       5 Corresponding to 10-5.
OPTIONS:
       User defined.
RECOMMENDATION:
       Use the default unless the molecule has high symmetry which is not being correctly identified. Note that relaxing this tolerance too much may introduce errors into the calculation.