X

Search Results

Searching....

4.3 Basic SCF Job Control

4.3.5 Symmetry

(April 13, 2024)

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. 1205 Takada T., Dupuis M., King H. F.
J. Chem. Phys.
(1981), 75, pp. 332.
Link
, 318 Dupuis M., King H. F.
Int. J. Quantum Chem.
(1977), 11, pp. 613.
Link
, 319 Dupuis M., King H. F.
J. Chem. Phys.
(1978), 68, pp. 3998.
Link
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 10).

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. POINT_GROUP_SYMMETRY controls symmetry throughout all modules. The default is TRUE. 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 POINT_GROUP_SYMMETRY keyword is set to FALSE then the coordinates will not be altered from the input, and the point group will be set to C1.

The INTEGRAL_SYMMETRY keyword controls symmetry in some integral routines. It is set to TRUE 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, INTEGRAL_SYMMETRY = TRUE can cause problems (poor convergence and wildly incorrect SCF energies) and turning it off can avoid these problems.

The symmetry is turned off when ghost atoms are used because their placement may break the symmetry of the molecular structure. However, the user can force calculation to use symmetry by FORCE_SYMMETRY_ON keyword; this is illustrated in example 7.11.8.2.

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.

INTEGRAL_SYMMETRY

INTEGRAL_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.

POINT_GROUP_SYMMETRY

POINT_GROUP_SYMMETRY
       Controls whether or not Q-Chem determines the point group of the molecule and reorients the molecule to the standard orientation.
TYPE:
       LOGICAL
DEFAULT:
       TRUE 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

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. Misidentified point group symmetry can affect rotational symmetry numbers in thermochemical calculations.

FORCE_SYMMETRY_ON

FORCE_SYMMETRY_ON
       Overrides turning off symmetry in calculations using ghost atoms.
TYPE:
       LOGICAL
DEFAULT:
       FALSE Turn symmetry off when using ghost atoms.
OPTIONS:
       TRUE Force symmetry. FALSE Do not use symmetry.
RECOMMENDATION:
       Use the default unless you know what you are doing.