Q-Chem 5.0 User’s Manual

C.1 Q-Chem Text Input Summary

C.1.1 Keyword: $molecule

Four methods are available for inputing geometry information:

C.1.2 Keyword: $rem

See also the list of $rem variables at the end of this Appendix. The general format is:

$rem
   REM_VARIABLE   VALUE   [optional comment]
$end

C.1.3 Keyword: $basis

The format for the user–defined basis section is as follows:

$basis

         
 

$X$

0

     
 

$L$

$K$

$scale$

   
 

$\alpha _1$

$C_1^{L_{min}} $

$C_1^{L_{min}+1}$

$\ldots $

$C_1^{L_{max}}$

 

$\alpha _2$

$C_2^{L_{min}} $

$C_2^{L_{min}+1}$

$\ldots $

$C_2^{L_{max}}$

 

$\vdots $

$\vdots $

$\vdots $

$\ddots $

$\vdots $

 

$\alpha _ K$

$C_ K^{L_{min}} $

$C_ K^{L_{min}+1}$

$\ldots $

$C_ K^{L_{max}}$

****

         

$end

         
where

$X$

Atomic symbol of the atom (atomic number not accepted)

$L$

Angular momentum symbol (S, P, SP, D, F, G)

$K$

Degree of contraction of the shell (integer)

$scale$

Scaling to be applied to exponents (default is 1.00)

$a_ i$

Gaussian primitive exponent (positive real number)

$C_ i^ L$

Contraction coefficient for each angular momentum (non–zero real numbers).

Atoms are terminated with **** and the complete basis set is terminated with the $end keyword terminator. No blank lines can be incorporated within the general basis set input. Note that more than one contraction coefficient per line is one required for compound shells like SP. As with all Q-Chem input deck information, all input is case–insensitive.

C.1.4 Keyword: $comment

Note that the entire input deck is echoed to the output file, thus making the $comment keyword largely redundant.

$comment
   User comments - copied to output file
$end

C.1.5 Keyword: $ecp

$ecp
For each atom that will bear an ECP
Chemical symbol for the atom
ECP name; the $L$ value for the ECP; number of core electrons removed
For each ECP component (in the order unprojected, $\hat{P}_0 $, $\hat{P}_1 $, , $\hat{P}_{L-1} $
The component name
The number of Gaussians in the component
For each Gaussian in the component
The power of $r$; the exponent; the contraction coefficient
****
$end

Note: (1) All of the information in the $ecp block is case–insensitive.
(2) The $L$ value may not exceed 4. That is, nothing beyond $G$ projectors is allowed.
(3) The power of $r$ (which includes the Jacobian $r^2$ factor) must be 0, 1 or 2.

C.1.6 Keyword: $empirical_dispersion

$empirical_dispersion
S6 S6_value
D D_value
C6 element_1 C6_value_for_element_1 element_2 C6_value_for_element_2
VDW_RADII element_1 radii_for_element_1 element_2 radii_for_element_2
$end

Note: This section is only for values that the user wants to change from the default values recommended by Grimme.

C.1.7 Keyword: $external_charges

All input should be given in atomic units.

Update: While charges should indeed be listed in atomic units, the units for distances depend on the user input. If the structure is specified in ngstroms (the default), the coordinates for external charges should also be in ngstroms. If the structure is specified in atomic units, the coordinates for external charges should also be in atomic units. (See INPUT_BOHR.)

$external_charges
   x-coord1   y-coord1   z-coord1   charge1
   x-coord2   y-coord2   z-coord2   charge2
$end

C.1.8 Keyword: $intracule

$intracule

int_type

0

Compute $P(u)$ only

 

1

Compute $M(v)$ only

 

2

Compute $W(u,v)$ only

 

3

Compute $P(u)$, $M(v)$ and $W(u,v)$

 

4

Compute $P(u)$ and $M(v)$

 

5

Compute $P(u)$ and $W(u,v)$

 

6

Compute $M(v)$ and $W(u,v)$

u_points

 

Number of points, start, end.

v_points

 

Number of points, start, end.

moments

0–4

Order of moments to be computed ($P(u)$ only).

derivs

0–4

order of derivatives to be computed ($P(u)$ only).

accuracy

$n$

($10^{-n})$ specify accuracy of intracule interpolation table ($P(u)$ only).

$end

C.1.9 Keyword: $isotopes

Note that masses should be given in atomic units.

$isotopes
   number_extra_loops  tp_flag
      number_of_atoms  [temp pressure]
      atom_number1 mass1
      atom_number2 mass2
   ...
$end

C.1.10 Keyword: $multipole_field

Multipole fields are all in atomic units.

$multipole_field
   field_component1   value1
   field_component2   value2
   ...
$end

C.1.11 Keyword: $nbo

Refer to Chapter 10 and the NBO manual for further information. Note that the NBO $rem variable must be set to ON to initiate the NBO package.

$nbo
   [ NBO options ]
$end

C.1.12 Keyword: $occupied

$occupied
   1  2  3  4 ...  nalpha
   1  2  3  4 ...  nbeta  
$end  

C.1.13 Keyword: $opt

Note that units are in ngstroms and degrees. Also see the summary in the next section of this Appendix.

$opt
CONSTRAINT
stre  atom1  atom2  value
...
bend  atom1  atom2  atom3  value
...
outp  atom1  atom2  atom3  atom4  value
...
tors  atom1  atom2  atom3  atom4  value
...
linc  atom1  atom2  atom3  atom4  value
...
linp  atom1  atom2  atom3  atom4  value
...
ENDCONSTRAINT

FIXED
atom   coordinate_reference
...
ENDFIXED

DUMMY
idum   type   list_length   defining_list
...
ENDDUMMY

CONNECT
atom   list_length   list
...
ENDCONNECT
$end

C.1.14 Keyword: $svp

$svp
<KEYWORD>=<VALUE>, <KEYWORD>=<VALUE>,...
<KEYWORD>=<VALUE>
$end

For example, the section may look like this:

$svp
    RHOISO=0.001, DIELST=78.39, NPTLEB=110
$end 

C.1.15 Keyword: $svpirf

 
$svpirf
    <# point> <x point> <y point> <z point> <charge> <grid weight>
    <# point> <x normal> <y normal> <z normal>
$end

C.1.16 Keyword: $plots

$plots
One comment line
Specification of the 3–D mesh of points on 3 lines:
$N_ x\quad x_{\ensuremath{\mathrm{min}}}\quad x_{\ensuremath{\mathrm{max}}}$
$N_ y\quad y_{\ensuremath{\mathrm{min}}}\quad y_{\ensuremath{\mathrm{max}}}$
$N_ z\quad z_{\ensuremath{\mathrm{min}}}\quad z_{\ensuremath{\mathrm{max}}}$
A line with 4 integers indicating how many things to plot:
$N_{\ensuremath{\mathrm{MO}}}\quad N_{\ensuremath{\mathrm{Rho}}}\quad N_{\ensuremath{\mathrm{Trans}}}\quad N_{\ensuremath{\mathrm{DA}}}$
An optional line with the integer list of MO’s to evaluate (only if $N_{\ensuremath{\mathrm{MO}}} > 0$)
MO(1)  MO(2) $\ldots $ MO($N_{\ensuremath{\mathrm{MO}}}$)
An optional line with the integer list of densities to evaluate (only if $N_{\ensuremath{\mathrm{Rho}}} > 0$)
Rho(1)  Rho(2) $\ldots $ Rho($N_{\ensuremath{\mathrm{Rho}}}$)
An optional line with the integer list of transition densities (only if $N_{\ensuremath{\mathrm{Trans}}}> 0$)
Trans(1)  Trans(2) $\ldots $ Trans($N_{\ensuremath{\mathrm{Trans}}}$)
An optional line with states for detachment/attachment densities (if $N_{\ensuremath{\mathrm{DA}}} > 0$)
DA(1)  DA(2) $\ldots $ DA($N_{\ensuremath{\mathrm{DA}}}$)
$end

C.1.17 Keyword: $localized_diabatization

$plots
One comment line.
One line with an an array of adiabatic states to mix together.
$<adiabat1> \;  \;  \;  \;  \;  \;  <adiabat2> \;  \;  \;  \;  \;  \;  \;  <adiabat3>$ $\ldots $
$end

Note:  We count adiabatic states such that the first excited state is $<adiabat> = 1$, the fifth is $<adiabat> = 5$, and so forth.

C.1.18 Keyword $van_der_waals

Note: All radii are given in ngstroms.

$van_der_waals
   1
   atomic_number   VdW_radius
$end

(alternative format)

$van_der_waals
   2
   sequential_atom_number   VdW_radius 
$end

C.1.19 Keyword: $xc_functional

$xc_functional
   X  exchange_symbol  coefficient
   X  exchange_symbol  coefficient
   ...
   C  correlation_symbol  coefficient
   C  correlation_symbol  coefficient
   ...
   K  coefficient
$end