A QM/MM job is requested by setting the $rem variables QM_MM_INTERFACE and FORCE_FIELD. Also required are a $qm_atoms input section and appropriate modifications to the $molecule section, as described above. Additional job control variables are detailed here.
QM_MM_INTERFACE
QM_MM_INTERFACE
      
Enables internal QM/MM calculations.
TYPE:
      
STRING
DEFAULT:
      
NONE
OPTIONS:
      
MM
Molecular mechanics calculation (i.e., no QM region)
ONIOM
QM/MM calculation using two-layer mechanical embedding
JANUS
QM/MM calculation using electronic embedding
RECOMMENDATION:
      
The ONIOM model and Janus models are described above. Choosing MM
leads to no electronic structure calculation. However, when using MM, one
still needs to define the $rem variables BASIS and EXCHANGE
in order for Q-Chem to proceed smoothly.
FORCE_FIELD
FORCE_FIELD
      
Specifies the force field for MM energies in QM/MM calculations.
TYPE:
      
STRING
DEFAULT:
      
NONE
OPTIONS:
      
AMBER99
AMBER99 force field
CHARMM27
CHARMM27 force field
OPLSAA
OPLSAA force field
RECOMMENDATION:
      
None.
CHARGE_CHARGE_REPULSION
CHARGE_CHARGE_REPULSION
      
The repulsive Coulomb interaction parameter for YinYang atoms.
TYPE:
      
INTEGER
DEFAULT:
      
550
OPTIONS:
      
Use Q = 
RECOMMENDATION:
      
The repulsive Coulomb potential maintains bond lengths involving YinYang
atoms with the potential . The default is parameterized for carbon atoms.
FORCEMAN_PRINT
FORCEMAN_PRINT
      
Controls printing level for MM (and thus QM/MM) jobs.
TYPE:
      
INTEGER
DEFAULT:
      
0
OPTIONS:
      
0
Minimal printing, as befits a large MM job that runs for many steps.
1
Additional information including MM gradient.
2
Print individual terms in the MM energy expression.
3
Print connectivity information.
4
Print individual terms in the MM gradient.
RECOMMENDATION:
      
Use the default unless trying to diagnose a problem.
GAUSSIAN_BLUR
GAUSSIAN_BLUR
      
Enables the use of Gaussian-delocalized external charges in a QM/MM calculation.
TYPE:
      
LOGICAL
DEFAULT:
      
FALSE
OPTIONS:
      
TRUE
Delocalizes external charges with Gaussian functions.
FALSE
Point charges
RECOMMENDATION:
      
None
GAUSS_BLUR_WIDTH
GAUSS_BLUR_WIDTH
      
Delocalization width for external MM Gaussian charges in a Janus calculations.
TYPE:
      
INTEGER
DEFAULT:
      
NONE
OPTIONS:
      
Use a width of  Å.
RECOMMENDATION:
      
Blur all MM external charges in a QM/MM calculation with the specified width.
Gaussian blurring is currently incompatible with PCM calculations. Values of
1.0–2.0 Å are recommended in Ref. 
      292
      
         
            
           J. Chem. Phys.
 
           (2002), 
           117,
           pp. 10534.
        
        
            
               Link
            
        
     
   
.
MM_SUBTRACTIVE
MM_SUBTRACTIVE
      
Specifies whether a subtractive scheme is used in the , Eq. (11.50), portion of the calculation.
TYPE:
      
LOGICAL
DEFAULT:
      
FALSE
OPTIONS:
      
FALSE
Only pairs that are not 1-2, 1-3, or 1-4 pairs are used.
TRUE
All pairs are calculated, and then the pairs that are double counted (1-2, 1-3, and 1-4) are subtracted out.
RECOMMENDATION:
      
When running QM/MM or MM calculations there is not recommendation. When running a QM/MM-Ewald calculation
the value must be set to TRUE.
MODEL_SYSTEM_CHARGE
MODEL_SYSTEM_CHARGE
      
Specifies the QM subsystem charge if different from the $molecule section.
TYPE:
      
INTEGER
DEFAULT:
      
NONE
OPTIONS:
      
The charge of the QM subsystem.
RECOMMENDATION:
      
This option only needs to be used if the QM subsystem (model system)
has a charge that is different from the total system charge.
MODEL_SYSTEM_MULT
MODEL_SYSTEM_MULT
      
Specifies the QM subsystem multiplicity if different from the $molecule section.
TYPE:
      
INTEGER
DEFAULT:
      
NONE
OPTIONS:
      
The multiplicity of the QM subsystem.
RECOMMENDATION:
      
This option only needs to be used if the QM subsystem (model system)
has a multiplicity that is different from the total system multiplicity.
ONIOM calculations must be closed shell.
PRINT_INPUT
PRINT_INPUT
      
Specifies whether to echo the input file in the output.
TYPE:
      
INTEGER
DEFAULT:
      
TRUE
OPTIONS:
      
TRUE
Print the Q-Chem input file as part of the Q-Chem output file.
FALSE
Do not repeat the input file.
RECOMMENDATION:
      
This is useful for record-keeping and defaults to TRUE except for QM/MM jobs,
where the default value is FALSE due to the size of some input files.
USER_CONNECT
USER_CONNECT
      
Enables explicitly defined bonds.
TYPE:
      
STRING
DEFAULT:
      
FALSE
OPTIONS:
      
TRUE
Bond connectivity is read from the $molecule section
FALSE
Bond connectivity is determined by atom proximity
RECOMMENDATION:
      
Set to TRUE if bond connectivity is known, in which case this connectivity must be
specified in the $molecule section. This greatly accelerates MM calculations.