The NEO method is a natural extension of the SCF method and it inherits many of its functionalities. Some keywords that are used in the SCF Job Control are used in the NEO methods with a few additional keywords. The following $rem variables must be specified in order to run NEO calculations:
Note: Some $rem variables, which we find mostly self-explanatory but not discussed below, that also influence NEO jobs in regards to geometry and symmetry considerations include INPUT_BOHR, SYM_IGNORE, and POINT_GROUP_SYMMETRY.
NEO
NEO
Enable a NEO-SCF calculation.
TYPE:
BOOLEAN
DEFAULT:
FALSE
OPTIONS:
TRUE
Enable a NEO-SCF calculation.
FALSE
Disable a NEO-SCF calculation.
RECOMMENDATION:
Set to TRUE if desired.
JOBTYPE
JOBTYPE
Specifies the type of calculation.
TYPE:
STRING
DEFAULT:
SP
OPTIONS:
SP
Single point energy.
FORCE
Analytical force calculation.
OPT
Geometry minimization.
TS
Transition structure search.
FREQ
Frequency calculation.
RECOMMENDATION:
Application-dependent. Always use POINT_GROUP_SYMMETRY = FALSE with geometry optimization.
METHOD
METHOD
Specifies the exchange-correlation functional.
TYPE:
STRING
DEFAULT:
No default
OPTIONS:
NAME
Use METHOD = NAME, where NAME is one of the following:
HF for Hartree-Fock theory;
one of the DFT methods listed in Section 5.3.5.;
RECOMMENDATION:
In general, consult the literature to guide your selection. Our recommendations for DFT are indicated
in bold in Section 5.3.5.
BASIS
BASIS
Specifies the electronic basis sets to be used.
TYPE:
STRING
DEFAULT:
No default basis set
OPTIONS:
General, Gen
User defined ($basis keyword required).
Symbol
Use standard basis sets as per Chapter 8.
Mixed
Use a mixture of basis sets (see Chapter 8).
RECOMMENDATION:
Consult literature and reviews to aid your selection.
The format for the user-defined basis section is as follows:
$basis 0 **** $end
where
Atomic symbol of the atom (atomic number not accepted) | |
Angular momentum symbol (S, P, SP, D, F, G) | |
Degree of contraction of the shell (integer) | |
Scaling to be applied to exponents (default is 1.00) | |
Gaussian primitive exponent (positive real number) | |
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 required for compound shells like SP. As with all Q-Chem input deck information, all input is case-insensitive.
Additionally, the NEO-SCF methods require definition of the nuclear basis sets.
Refer to Ref.
1462
J. Chem. Phys.
(2020),
152,
pp. 244123.
Link
for selection of the protonic basis sets. Note that the pure (spherical harmonic) and Cartesian angular forms for the protonic basis sets are found
in the main article and in the supported information, respectively. Also, see details about the NEO_PURECART variable described in the next section.
The format for the user-defined neo_basis section is as follows:
$neo_basis 0 **** $end
Note: Several of the post-NEO-HF methods utilizing density fitting also require the nuclear auxiliary basis set to be defined. The format for such user-defined $neo_aux_basis section is identical to the $basis and $neo_basis input sections. See examples in Section 13.5.6 for more details.
More $rem variables that can be used to customize the calculation. Some of these are NEO specific variables, while others are shared with other Q-Chem modules:
BASIS_LIN_DEP_THRESH
BASIS_LIN_DEP_THRESH
Sets the threshold for determining linear dependence in the basis set
TYPE:
INTEGER
DEFAULT:
6
Corresponding to a threshold of
OPTIONS:
Sets the threshold to
RECOMMENDATION:
Set to 5 or smaller if you have a poorly behaved SCF and you suspect linear
dependence in you basis set. Lower values (larger thresholds) may affect the
accuracy of the calculation.
NEO_BASIS_LIN_DEP_THRESH
NEO_BASIS_LIN_DEP_THRESH
This keyword is used to set the liner dependency threshold for nuclear basis sets. It is defined as .
TYPE:
DOUBLE
DEFAULT:
5.0
OPTIONS:
User-defined
RECOMMENDATION:
No recommendation.
PURECART
PURECART
INTEGER
TYPE:
Controls the use of pure (spherical harmonic) or Cartesian angular forms
DEFAULT:
1111
Pure functions
OPTIONS:
Use 1 for pure and 2 for Cartesian.
RECOMMENDATION:
This is pre-defined for all standard basis sets
NEO_PURECART
NEO_PURECART
This keyword is used to specify Cartesian or spherical Gaussians for nuclear basis functions.
TYPE:
INTEGER
DEFAULT:
2222
OPTIONS:
User-defined
RECOMMENDATION:
The default value corresponds to the use of Cartesian Gaussians for all angular momentum classes.
The value NEO_PURECART = 1111 would use spherical Gaussians instead, similar to the use of PURECART.
NEO_ISOTOPE
NEO_ISOTOPE
Enable calculations of different types of isotopes. Only one type of isotope is allowed at present.
TYPE:
INTEGER
DEFAULT:
1
Default is the proton isotope.
OPTIONS:
1
This NEO calculation is using proton isotope.
2
This NEO calculation is using deuterium isotope.
3
This NEO calculation is using tritium isotope.
RECOMMENDATION:
Refer to the NEO literature for the best performance on the isotope effects calculations.
NEO_VPP
NEO_VPP
This keyword is used to control whether to remove terms from the nuclear Fock matrix and the corresponding kernel terms for NEO excited state methods.
The purpose for the removal of these terms in the case of one quantum proton is to save on computational cost. Note that the ground state NEO-HF or NEO-DFT energy
is the same with or without removal of these nuclear terms. The nuclear virtual orbitals will however be different, which may affect the results for perturbative
methods (i.e., multicomponent excited-state methods and post-NEO-HF wavefunction methods).
TYPE:
LOGICAL/INTEGER
DEFAULT:
TRUE
OPTIONS:
TRUE (or 1)
Enable this option (include nuclear terms).
FALSE (or 0)
Disable this option (remove nuclear terms).
RECOMMENDATION:
Set NEO_VPP = 0 only in the case of one quantum hydrogen.
XC_GRID
XC_GRID
Specifies the type of grid to use for DFT calculations.
TYPE:
INTEGER
DEFAULT:
Functional-dependent; see Table 5.3.
OPTIONS:
0
Use SG-0 for H, C, N, and O; SG-1 for all other atoms.
Use SG- for all atoms, , or 3
A string of two six-digit integers and , where is the number of radial points
and is the number of angular points where possible numbers of Lebedev angular
points, which must be an allowed value from Table 5.2 in
Section 5.5.
Similar format for Gauss-Legendre grids, with the six-digit integer corresponding
to the number of radial points and the six-digit integer providing the number of
Gauss-Legendre angular points, .
RECOMMENDATION:
We strongly recommend using an unpruned, Euler-Maclaurin-Lebedev grid type for NEO-DFT calculations. Larger grids may be required for optimization and frequency calculations.
UNRESTRICTED
UNRESTRICTED
Controls the use of restricted or unrestricted orbitals.
TYPE:
LOGICAL
DEFAULT:
FALSE
Closed-shell systems.
TRUE
Open-shell systems.
OPTIONS:
FALSE
Constrain the spatial part of the alpha and beta orbitals to be the same.
TRUE
Do not Constrain the spatial part of the alpha and beta orbitals.
RECOMMENDATION:
The ROHF method is not available. Note that for unrestricted calculations on
systems with an even number of electrons it is usually necessary to break
/ symmetry in the initial guess, by using SCF_GUESS_MIX or
providing $occupied information (see Section 4.4 on initial guesses).
S2THRESH
S2THRESH
Cutoff for neglect of overlap integrals, defined via a two-electron shell-pair threshold of (S2THRESH
).
TYPE:
INTEGER
DEFAULT:
Same as THRESH.
OPTIONS:
for a threshold of .
RECOMMENDATION:
Increase the value of S2THRESH if the program finds negative eigenvalues for the overlap matrix.
THRESH
THRESH
Cutoff for neglect of two electron integrals. (THRESH
).
TYPE:
INTEGER
DEFAULT:
8
For single point energies.
10
For optimizations and frequency calculations.
14
For coupled-cluster calculations.
OPTIONS:
for a threshold of .
RECOMMENDATION:
The value should satisfy THRESH 3 + SCF_CONVERGENCE, although
THRESH = 4 + SCF_CONVERGENCE is the default (in most cases) since Q-Chem v. 6.0.
See Ref.
459
J. Phys. Chem. A
(2024),
128,
pp. 7739.
Link
for recommended values of THRESH in the presence of diffuse
basis functions, where tighter thresholds are often required and too-loose thresholds may lead to slower
convergence or convergence failure.
MEM_STATIC
MEM_STATIC
Sets the memory for AO-integral evaluations and their transformations in Q-Chem 4.1 or older versions.
TYPE:
INTEGER
DEFAULT:
192
corresponding to 192 MB.
OPTIONS:
User-defined number of megabytes.
RECOMMENDATION:
Use the default.
MEM_TOTAL
MEM_TOTAL
Sets the total memory available to Q-Chem, in megabytes.
TYPE:
INTEGER
DEFAULT:
2000
2 GB
OPTIONS:
User-defined number of megabytes.
RECOMMENDATION:
Use the default, or set to the physical memory of your machine.
The minimum requirement is .
The following $rem variables can be used to specify the tolerances for NEO-SCF convergence. Some of these $rem variables are inherited from conventional SCF calculations:
SCF_CONVERGENCE
SCF_CONVERGENCE
NEO-SCF is considered converged when the electronic wave function error is less that
. Adjust the value of THRESH at the same
time. (Starting with Q-Chem 3.0, the DIIS error is measured by the maximum error
rather than the RMS error as in earlier versions.)
TYPE:
INTEGER
DEFAULT:
5
For single point energy calculations.
8
For geometry optimizations.
OPTIONS:
User-defined
RECOMMENDATION:
None.
NEO_N_SCF_CONVERGENCE
NEO_N_SCF_CONVERGENCE
NEO-SCF is considered converged when the nuclear wave function error is less that
.
TYPE:
INTEGER
DEFAULT:
7
OPTIONS:
User-defined
RECOMMENDATION:
None.
NEO_E_CONV
NEO_E_CONV
Energy convergence criteria in the NEO-SCF calculations so that the difference in energy between consecutive iterations is less than
.
TYPE:
INTEGER
DEFAULT:
8
OPTIONS:
User-defined
RECOMMENDATION:
Tighter criteria for geometry optimization are recommended. Note that for stepwise calculations, this convergence criteria checks
the difference in energy between electronic and protonic microiterations blocks.
MAX_SCF_CYCLES
MAX_SCF_CYCLES
Controls the maximum number of SCF iterations permitted.
TYPE:
INTEGER
DEFAULT:
50
OPTIONS:
User-selected.
RECOMMENDATION:
Increase for slowly converging systems such as those containing transition
metals.
NEO_E_MAX_SCF_CYCLES
NEO_E_MAX_SCF_CYCLES
Controls the maximum number of electronic SCF micoiterations permitted. This variable is only used for
stepwise NEO-SCF optimizations.
TYPE:
INTEGER
DEFAULT:
150
OPTIONS:
User-selected.
RECOMMENDATION:
None.
NEO_N_MAX_SCF_CYCLES
NEO_N_MAX_SCF_CYCLES
Controls the maximum number of nuclear SCF micoiterations permitted. This variable is only used for
stepwise NEO-SCF optimizations.
TYPE:
INTEGER
DEFAULT:
150
OPTIONS:
User-selected.
RECOMMENDATION:
None.
The NEO-SCF iteration procedure has traditionally involved solving the
electronic and quantum nuclear Roothaan/Kohn-Sham equations in alternating fashion until convergence to the variational solution is
reached. This approach can become computationally expensive, especially when many SCF iterations are required to achieve self-consistency.
To accelerate this procedure, those equations may instead be solved simultaneously at every step of the NEO-SCF iteration procedure.
This is now the default scheme for solving the NEO-SCF equations. We support the simultaneous implementation of the Direct Inversion in the Iterative Subspace (DIIS) algorithm
809
J. Phys. Chem. A
(2022),
126,
pp. 7033.
Link
as well as the Geometric Direct Minimization (GDM)
1331
Mol. Phys.
(2002),
100,
pp. 1713.
Link
approach. The following $rem variables permit some customization of the NEO-SCF iterations:
NEO_SIMULTANEOUS_SCF
NEO_SIMULTANEOUS_SCF
Enables simultaneous optimization algorithm.
TYPE:
LOGICAL
DEFAULT:
TRUE
OPTIONS:
TRUE
FALSE
RECOMMENDATION:
Use the default unless issues with convergence are observed, in which cause set to FALSE to utilize the traditional stepwise approach.
SCF_ALGORITHM
SCF_ALGORITHM
Algorithm used for converging the SCF.
TYPE:
STRING
DEFAULT:
DIIS
Pulay DIIS.
OPTIONS:
DIIS
Pulay DIIS.
GDM
Geometric Direct Minimization.
DIIS_GDM
Use DIIS and then later switch to geometric direct minimization.
RECOMMENDATION:
In the NEO methods, the GDM procedure is recommended. Note that DIIS_GDM is only supported for simultaneous optimizations.
NEO_SCF_GUESS_N
NEO_SCF_GUESS_N
Specifies the initial guess procedure for the nuclear orbitals in a NEO-SCF calculation.
TYPE:
STRING
DEFAULT:
TIGHT
Default unless ghost atoms are present.
CORE
Default when there are ghost atoms.
OPTIONS:
CORE
Diagonalize the nuclear core Hamiltonian
TIGHT
Occupies the least diffuse nuclear MO
AUTOSAD
On-the-fly superposition of nuclear atomic densities from an atomic NEO calculation
READ
Read previous nuclear MOs from disk
RECOMMENDATION:
Use the default. Note that if NEO_SCF_GUESS_N is set to READ, only the nuclear MOs from a previous NEO-SCF
job are read from disk. The starting guess for the electronic MOs in a NEO calculation remains the converged result from the
conventional electronic calculation of the current job. To read in both the electronic and nuclear MOs from a previous NEO-SCF job,
NEO_SCF_GUESS_N must be set to READ in conjunction with SKIP_SCFMAN set to TRUE.
NEO_STEPWISE_SCF_STEPS
NEO_STEPWISE_SCF_STEPS
Specifies the number of NEO-SCF stepwise/alternating macro-iterations to perform before switching to simultaneous algorithm.
TYPE:
INTEGER
DEFAULT:
0
OPTIONS:
User-defined
RECOMMENDATION:
The rate of convergence of the NEO-SCF
procedure is dependent on the initial guess for the electronic and protonic orbitals. For especially difficult systems,
it may be helpful to do stepwise optimization before switching to simultaneous.
Note: The following variables are specific to SCF_ALGORITHM = DIIS.
DIIS_SUBSPACE_SIZE
DIIS_SUBSPACE_SIZE
Controls the size of the DIIS subspace during the SCF.
TYPE:
INTEGER
DEFAULT:
10
OPTIONS:
User-defined
RECOMMENDATION:
None.
DIIS_ERR_RMS
DIIS_ERR_RMS
Changes the DIIS convergence metric from the maximum to the RMS error.
TYPE:
LOGICAL
DEFAULT:
FALSE
OPTIONS:
TRUE
FALSE
RECOMMENDATION:
Use the default, the maximum error provides a more reliable criterion.
Note: The following variables are specific to SCF_ALGORITHM = DIIS_GDM.
MAX_DIIS_CYCLES
MAX_DIIS_CYCLES
The maximum number of DIIS iterations before switching to (geometric) direct
minimization.
TYPE:
INTEGER
DEFAULT:
50
OPTIONS:
1
Only a single Roothaan step before switching to (G)DM
DIIS iterations before switching to (G)DM.
RECOMMENDATION:
None
THRESH_DIIS_SWITCH
THRESH_DIIS_SWITCH
The threshold for switching between DIIS extrapolation and direct minimization
of the SCF energy is .
TYPE:
INTEGER
DEFAULT:
2
OPTIONS:
User-defined.
RECOMMENDATION:
None