- Description
- Options
- Restrictions
- Related Keywords
- Examples
Options
Algorithm-Related Options
Harris
Diagonalize the Harris functional [Harris85] for the initial guess. This is the default for all HF and DFT calculations.
Huckel
Requests that a Huckel guess be generated. This is the default for CNDO, INDO, MNDO, and MINDO3. The Huckel guess should be considered for PM6 calculations involving many second-row atoms.
RdScale
Read in the scale factor on atomic hardnesses used in iterative extended Huckel. The default is 7.0 times the QEq value.
OldHuckel
Use the old Huckel guess (pre-Gaussian 03) instead of CNDO or the updated Huckel.
INDO
Use the Gaussian 16 default guess: INDO for first-row systems and Huckel for second-row and beyond.
AM1
Do an AM1 calculation for the initial guess (currently only works with sparse matrix code). Guess=(AM1,Always) causes later steps in a geometry optimization to generate a new guess at each point and compare the energies with the density from the old point and the new guess and take the best one.
Core
Requests that the core Hamiltonian be diagonalized to form the initial guess. This is the default for AM1, PM3, PM3MM, PM6, and PDDG.
Orbital-Related Options
Permute
Read in a permutation of orbitals in the initial guess. The numbers of the generated guess orbitals are given in the order in which they should be used in the SCF. Ranges (e.g. 7–12) can be used, and all orbitals not listed are put in after the listed orbitals in their original order. Separate permutation lists for α and β orbitals must be specified (each list separated by a blank line) for open shell systems.
Alter
Indicates that the orbitals selected for occupation in the Hartree-Fock wavefunction should not be those of lowest energy. Normally, the occupied orbitals are selected as those with lowest eigenvalues for the one-electron Hamiltonian used in the initial guess programs. The alteration sections consist of a set of transpositions indicating that one of these occupied orbitals is to be replaced by one of the other (virtual) orbitals. Each such transposition is on a separate line and has two integers N1 and N2 (free format, separated by spaces or a comma as usual) indicating that orbital N1 is to be swapped with orbital N2. The list of orbital transpositions is terminated by the blank line at the end of the input section.
For UHF calculations, two such orbital alteration sections are required, the first specifying transpositions of α orbitals and the second specifying transpositions of β orbitals. Both sections are always required. Thus, even if only α transpositions are needed, the β section is required even though it is empty (and vice-versa). The second blank line to indicate an empty β section must be included.
Mix
Requests that the HOMO and LUMO be mixed so as to destroy α-β and spatial symmetries. This is useful in producing UHF wavefunctions for singlet states. Mixing orbitals is only done by default when generating a complex initial guess. NoMix says to not mix orbitals.
DensityMix[=N]
Whether to mix occupied and virtual orbital contributions in forming the initial guess density. N defaults to -3 (use Huckel eigenvalues to decide which orbitals to mix). Pre-defined values of N are 0 (no mixing), -1 (mix all orbitals equally), -2 (mix valence orbitals and an equal number of virtual orbitals), and -3 (mix valence occupied orbitals with 0.05 au (according to ZDO) of the HOMO and virtual orbitals within 0.15 au); a positive value for N says to use equal occupations of the lowest N virtual orbitals and highest N occupied orbitals.
CopyChk
Copy the MOs and density unedited, unprojected and unchecked from the checkpoint file for the guess. Useful when importing from matrix element files.
Biorthogonalize
For an unrestricted guess, biorthogonalize the α and β MOs to maximally pair electrons of opposite spin. Done automatically when a UHF wavefunction is read-in as a guess for ROHF. This option is also useful in the combination of Guess=(BiOrth,Read,Only,Save) to replace the canonical UHF orbitals with ones which have alpha and beta orbitals matched up as much as possible for viewing in GaussView or other visualization packages.
NaturalOrbitals
Include natural orbitals in the checkpoint file. This must be accomplished via a separate job step specifying this option as well as Only and Save. See the discussion of the Population keyword for details.
Procedural Options
Only
Guess=Only functions as a calculation type keyword and requests that the calculation terminate once the initial guess is computed and printed. Note that the amount of orbital information that is printed is controlled by the Pop keyword. Guess=Only may not be used with the MOPAC-based semi-empirical methods (INDO, CNDO, MNDO, MINDO3).
This option is useful in preliminary runs to check if configuration alteration is necessary. For example, Guess=Only may be specified with CASSCF in order to obtain information on the number of CI configurations in the CAS active space (as well as the initial orbitals).
Guess=(Only,Read) may also be used to produce population and other post-calculation analyses from the data in a checkpoint file. For example, these options alone will produce a population analysis using the wavefunction in the checkpoint file. Guess(Only,Read) Prop will cause electrostatic properties to be calculated using the wavefunction in the checkpoint file.
Always
Requests that a new initial guess be generated at each point of an optimization. By default, the SCF results from the last point are used for the guess at the next point.
Fragment=N
Generate a guess built from fragment guesses or SCF solutions. The guess is saved on the checkpoint file regardless of if it is used in the current calculation, allowing it to be used in subsequent ones. This option is usually combined with Guess=Only to generate a guess from fragment guess orbitals (otherwise a full SCF calculation is performed for each fragment). Assignment of atoms to fragments and fragment charges and multiplicities are specified as described in Molecule Specifications, except that a negative spin multiplicity means that the unpaired orbitals for the fragment are to become β spin orbitals in the combined set.
Local
Requests that orbitals be localized using the Boys method [Boys60]. Occupied and virtual orbitals are localized separately, and the irreducible representations (after possible merging using LowSymm or NoSymm) are not mixed. Localized orbital analysis of a converged SCF wavefunction may then be done using a second job step, which includes Guess(Read,Local,Only) and Pop=Full in its route section.
Sparse
Perform a sparse semi-empirical calculation for the initial guess. This option may be useful for very, very large HF or DFT calculations using the sparse matrix facility.
Extra
Do an extra, new initial guess when reading orbitals from the RWF (i.e., during geometry optimizations). By default, this is done if the default Harris guess is allowed, no alteration of configuration was requested, and the optimization did not take a small step as flagged by variable 4 in ILSW (interlink status word). Use NoExtra to disable this feature.
Fock
Reuse Fock matrices rather than orbitals when reading from previous results from the read-write or checkpoint files. This is the default for periodic boundary conditions calculations if Guess=Alter is not specified. NoFock disables this behavior, and it is the default for non-PBC calculations.
Guess Reading/Saving Options
Read
Requests that the initial guess be read from the checkpoint file (Guess=Read is often specified along with Geom=Checkpoint). This option may be combined with Alter, in which case the orbitals are read from the checkpoint file, projected onto the current basis set, and then the specified alterations are made. Checkpoint is a synonym for Read.
TCheck
Read SCF restart information from the checkpoint file if possible; otherwise, read the initial guess. Designed for multistep jobs so that the SCF can be restarted in a later step when possible or it can be run with Guess=Read otherwise. This option implies CIS=Read or TD=Read when applicable/appropriate.
Alpha
Use alpha orbitals for both alpha and beta guess during Guess=Read.
Translate
Translate requests that the coordinates of the atoms used to produce a guess, which is read in, be translated to the current atomic coordinates. This is the default. It may fail in unusual cases, such as when a wavefunction is used as a guess for a system with a different stoichiometry, in which case Guess=NoTranslate should be specified.
Cards
Specifies that after the initial guess is generated, some or all of the orbitals will be replaced with ones read from the input stream. This option can be used to read a complete initial guess from the input stream by replacing every orbital. The replacement orbitals are placed in the input section following the guess alteration commands, if any. For UHF, there are separate α and β replacement orbital input sections.
The replacement orbitals input section (the α replacement orbitals section for UHF) begins with a line specifying the Fortran format with which to read the replacement orbital input, enclosed in parentheses. For example: (4E20.8). The remainder of the section contains one or more instances of the following:
IVec Orbital to replace (0=end, -1=replace all orbitals in order). (A(I,IVec),I=1,N) New orbital in the format specified in the first line.
The format for the line containing IVec is Fortran I5. The β orbital replacement section for UHF calculations differs only in that it omits the initial format specification line. See the examples section for sample replacement orbital input.
Input
Read a line from the input file containing the name of a checkpoint file. Instead of a filename, you can specify one of the following keywords:
- generate says to generate a guess as usual.
- read and chk say to retrieve data from the job’s checkpoint file. They are thus equivalent to Guess=Read, and use the checkpoint file specified to %Chk or %OldChk as usual.
- none says to skip generating a guess at all.
- dummy generates a dummy guess that does not require any work.
See the examples section for information about using this option with ONIOM.
Restart
Start a new SCF calculation using the restart data saved on the checkpoint file from a job with a different geometry and/or basis set. (The option name is a bit of a misnomer.)
Save
Save the generated initial guess back into the checkpoint file at the conclusion of a Guess=Only run. This option is useful for saving localized orbitals.
Print the initial guess.
Symmetry-Related Options
LowSymm
Requests that irreducible representations of the molecular point group be combined in the symmetry information used in the N3 steps in the SCF to allow lowered symmetry of the wavefunction. This enables the orbitals (and possibly but not necessarily the total wavefunction) to have lower symmetry than the full molecular point group. This option is available only for GVB calculations, where it is often necessary for calculations on symmetric systems (see the discussion of the GVB keyword for an example using this option).
The option expects a single line of input (in the format 16I2) giving the numbers of the irreducible representations to combine, with the new groups separated by 0; the list itself must be terminated by a 9. The numbers correspond to the order in which the representations are listed by Link 301 in the output file (see the Examples).
Since this input section is always exactly one line long, it is not terminated by a blank line. Note that irreducible representations are combined before orbital localization is done and that localized orbitals retain whatever symmetry is kept. Guess=NoSymm removes all orbital symmetry constraints without reading any input.
NoSymm
Requests that all orbital symmetry constraints be lifted. Synonymous with SCF=NoSymm and Symm=NoSCF.
ForceAbelianSymmetry
Force the initial guess orbitals to transform according to irreducible representations of the Abelian point group. NoForceAbelianSymmetry is the default.
Valid Option Combinations
Only reasonable combinations are valid. For example Guess=(Always,Alter) and Guess=(Read,Alter) work as expected (in the former case, alterations are read once and the same interchanges are applied at each geometry). Conversely, Guess=(Always,Read) is contradictory and will lead to unpredictable results. Refer to the input sections order table to determine the ordering of the input sections for combinations of options like Guess=(Cards,Alter).
Examples
Transposing Two Orbitals with Guess=Alter
This example finds the UHF/STO-3G structure of the 2A1 excited state of the amino radical. First, a Guess=Only calculation is run to determine whether any alter (reordering) instructions are needed to obtain the desired electronic state. The HF/STO-3G theoretical model is used by default:
# Guess=Only Test Amino radical test of initial guess 0 2 n h 1 nh h 1 nh 2 hnh nh 1.03 hnh 120.0
Here is the orbital symmetry summary output from the job, which comes immediately before the population analysis in the output:
Initial guess orbital symmetries:
Alpha Orbitals:
Occupied (A1) (A1) (B2) (A1) (B1)
Virtual (A1) (B2)
Beta Orbitals:
Occupied (A1) (A1) (B2) (A1)
Virtual (B1) (A1) (B2)
The electronic state of the initial guess is 2-B1.
Initial guess <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7500 S= 0.5000
Since a doublet state is involved, α and β orbitals are given separately. From the orbital symmetries, the electron configuration in the initial guess is 2a12a12b22a11b1, yielding a 2B1 wavefunction. This is indeed the ground state of NH2. The expectation value of S2 for the unrestricted initial guess is printed, in this case, that corresponding to a pure doublet: 0.75.
Since we want to model the 2A1 excited state of the amino radical, we will need to alter this initial orbital configuration: a β electron must be moved from orbital 4 to orbital 5 (the electron configuration is then 2a12a12b22b11a1).
Guess=Alter may be used to accomplish this. Here is the input for the geometry optimization:
# UHF/6-31G(d) Opt Guess=Alter Amino radical: HF/6-31G(d) Structure of 2-A1 state 0 2 N H 1 nh H 1 nh 2 hnh Variables: nh 1.03 hnh 120.0 Blank line ends the molecule specification section. Blank line ends the alpha section (empty in this case). 4 5 Transpose beta orbitals 4 and 5. End of the beta alteration section.
Note that an extra blank line—line 14—is necessary to indicate an empty α alteration section. The final two lines then constitute the β alteration section.
The initial guess program prints a list of orbitals that were interchanged as a result of the Alter option:
Harris functional with IExCor= 205 diagonalized for initial guess.
…
No Alpha orbitals switched.
Pairs of Beta orbitals switched:
4 5
Initial guess orbital symmetries:
Alpha Orbitals:
Occupied (A1) (A1) (B2) (A1) (B1)
Virtual (A1) (B2) (B1) (A1) (B2) (A1) (B2) (A1) (A2) (A1)
(B1) (A1) (B2) (A1)
Beta Orbitals:
Occupied (A1) (A1) (B2) (B1)
Virtual (A1) (A1) (B2) (B1) (A1) (B2) (A1) (B2) (A1) (A2)
(A1) (B1) (A1) (B2) (A1)
The electronic state of the initial guess is 2-A1.
Initial guess <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7500 S= 0.5000
The SCF calculation follows, and the energy and S2 eigenvalue for the UHF wavefunction are printed. The S2 value which results if contamination of the wavefunction from the next possible spin multiplicity (quartets for doublets, quintets for triplets, etc.) is removed is also printed:
SCF Done: E(UHF) = -55.4915172451 A.U. after 12 cycles
Convg = 0.2693D-08 -V/T = 2.0038
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7534 S= 0.5017
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation 0.7534, after 0.7500
Although this calculation does in fact converge correctly to 2A1 state, it sometimes happens that the order of orbital symmetries switches during the course of the SCF iterations. If the orbital symmetries of the final wavefunction are different from those in the initial guess (whether or not you are using Guess=Alter), we recommend using the direct minimization routine, specified with the SCF=QC or SCF=DM keywords, which usually holds symmetry from one iteration to the next.
Reordering Orbitals with Guess=Permute
This option is often is the easiest way to perform a complex modification of the initial guess, as in this example:
# CASSCF/6-31G(d,p) Opt Guess=Permute Pop=Reg CAS job 0 1 molecule specification 1-60 65 63 64 66 68 67 61-62 69 Specify new orbital ordering.
Here, we have rearranged orbitals 61–68. Listing the final orbital (69) is not really necessary, but it helps to make the input easier to understand for humans.
Specifying Calculation-Specific Guesses for ONIOM with Guess=Input
The Guess=Input is useful for ONIOM jobs where the SCF is hard to converge for some of the sub-calculations at the initial geometry. Convergence can be accomplished using a procedure like the following:
Step 1. Use ONIOM=OnlyInput to print out input files for the individual calculations.
Step 2. Converge the SCF for the problem step separately (for example, using Stable=Opt).
Step 3. After verifying that the correct state has been found, run an ONIOM calculation with Guess=Input to retrieve the correct initial guesses for the individual states. In the ONIOM input file, you will then have a line for the guess for each calculation (i.e., 3 lines for a 2-layer MO:MO), where each line can be a checkpoint file name or read (for Guess=Read from the regular checkpoint file) or generate (generate a new initial guess), as appropriate. For example, if only the high level calculation on the model system was problematic, the input would be:
| generate | Low level calculation on the real system. |
| Blank line. | |
| guess_for_high.chk | High level calculation on the model system. |
| Blank line. | |
| generate | Low level calculation on the model system. |
Reading in Orbitals with Guess=Cards
Some or all of the orbitals may be replaced, after the initial guess is generated, using Guess=Cards. Here is some sample input for this option, which replaces orbitals 1 and 4 (note that the format for the third and following lines is specified in line 1):
(3E20.8)
1
0.5809834509E+00 0.4612416518E+00 -0.6437319952E-04
0.1724432549E-02 0.1282235396E-14 0.5417658499E-13
0.1639966912E-02 -0.9146282229E-15 -0.6407549694E-13
-0.4538843604E-03 0.6038992958E-04 -0.1131035485E-03
0.6038992969E-04 -0.1131035471E-03
4
0.7700779642E-13 0.1240395916E-12 -0.3110890228E-12
-0.4479190461E-12 -0.1478805861E-13 0.5807753928E+00
0.6441113412E-12 -0.3119296374E-14 0.1554735923E+00
-0.1190754528E-11 0.2567325943E+00 0.1459733219E+00
-0.2567325943E+00 -0.1459733219E+00
0
An orbital number of zero ends the replacement orbital input.
Antiferromagnetic Coupling
Here is an example of a fragment guess job. The first step generates guesses for each one of the fragments and then combines them together into a guess for the SCF on the full molecule. The second step reads in that guess for a subsequent calculation. The charge and multiplicity line has the overall molecular charge and multiplicity first, followed by those for each fragment.
This example is Fe2S2 + 4 S-R ligands, where R is a phenyl group. It is a singlet with an overall charge of -2. The guess does each bare sulfur as S(2-) closed-shell (fragments 2 and 4), the two irons as Fe(3+) sextets antiferromagnetically coupled (fragments 1 and 3, with 1 being alpha-spin and 3 beta-spin), and each of the S-R fragments (5-8) as closed-shell singlet anions.
Fragment Assignments
%chk=FragGuess %mem=64mw #P UBP86/6-311G(d) Guess=(Fragment=8,Only) Pop=None Fe2S2 cluster with phenylthiolates. Step 1: Generate fragment guess -2,1 3,6 -2,1 3,-6 -2,1 -1,1 -1,1 -1,1 -1,1 H(Fragment=7) 23.5010 2.2873 8.5744 S(Fragment=2) 14.8495 1.1490 7.0431 Fe(Fragment=3) 17.0430 1.0091 7.0068 S(Fragment=4) 17.4565 -1.1490 7.0431 S(Fragment=5) 14.3762 -2.1581 8.7983 C(Fragment=5) 12.5993 -2.1848 8.6878 H(Fragment=5) 12.3743 -3.6513 10.1678 C(Fragment=5) 10.4994 -3.1122 9.4309 H(Fragment=5) 9.9929 -3.7579 10.0022 C(Fragment=5) 9.8049 -2.2791 8.5639 H(Fragment=5) 8.8050 -2.2873 8.5744 C(Fragment=5) 10.4833 -1.4146 7.6615 H(Fragment=5) 9.9730 -0.8525 7.0106 S(Fragment=8) 14.3794 -1.8091 5.0446 C(Fragment=5) 11.9048 -1.3675 7.7057 H(Fragment=5) 12.4158 -0.7843 7.0743 C(Fragment=6) 17.2999 3.4265 4.6624 C(Fragment=6) 16.6376 4.1967 5.6090 H(Fragment=6) 16.5022 3.8494 6.5369 C(Fragment=6) 16.1530 5.4856 5.2463 H(Fragment=6) 15.6665 6.0472 5.9155 C(Fragment=6) 16.3468 5.9257 4.0431 H(Fragment=6) 16.0236 6.8408 3.8020 C(Fragment=6) 17.0091 5.1398 3.0522 H(Fragment=6) 17.1330 5.4944 2.1254 C(Fragment=6) 17.4775 3.8823 3.3884 H(Fragment=6) 17.9400 3.3149 2.7071 S(Fragment=7) 17.9298 2.1581 8.7983 C(Fragment=7) 19.7067 2.1848 8.6878 C(Fragment=7) 20.4174 3.0650 9.5194 H(Fragment=7) 19.9317 3.6513 10.1678 C(Fragment=7) 21.8066 3.1122 9.4309 H(Fragment=7) 22.3132 3.7579 10.0022 C(Fragment=7) 22.5011 2.2791 8.5639 C(Fragment=7) 21.8227 1.4146 7.6615 H(Fragment=7) 22.3330 0.8525 7.0106 C(Fragment=7) 20.4012 1.3675 7.7057 H(Fragment=7) 19.8902 0.7843 7.0743 C(Fragment=8) 15.0061 -3.4265 4.6624 Fe(Fragment=1) 15.2630 -1.0091 7.0068 C(Fragment=8) 15.6684 -4.1967 5.6090 H(Fragment=8) 15.8038 -3.8494 6.5369 C(Fragment=8) 16.1530 -5.4856 5.2463 H(Fragment=8) 16.6395 -6.0472 5.9155 C(Fragment=5) 11.8886 -3.0650 9.5194 C(Fragment=8) 15.9592 -5.9257 4.0431 H(Fragment=8) 16.2824 -6.8408 3.8020 C(Fragment=8) 15.2969 -5.1398 3.0522 H(Fragment=8) 15.1730 -5.4944 2.1254 C(Fragment=8) 14.8285 -3.8823 3.3884 H(Fragment=8) 14.3660 -3.3149 2.7071 S(Fragment=6) 17.9266 1.8091 5.0446 --Link1-- %chk=FragGuess %mem=64mw #P UBP86/6-311G*/Auto Guess=Read Geom=AllCheck …
- Description
- Options
- Restrictions
- Related Keywords
- Examples
Algorithm-Related Options
Harris
Diagonalize the Harris functional [Harris85] for the initial guess. This is the default for all HF and DFT calculations.
Huckel
Requests that a Huckel guess be generated. This is the default for CNDO, INDO, MNDO, and MINDO3. The Huckel guess should be considered for PM6 calculations involving many second-row atoms.
RdScale
Read in the scale factor on atomic hardnesses used in iterative extended Huckel. The default is 7.0 times the QEq value.
OldHuckel
Use the old Huckel guess (pre-Gaussian 03) instead of CNDO or the updated Huckel.
INDO
Use the Gaussian 16 default guess: INDO for first-row systems and Huckel for second-row and beyond.
AM1
Do an AM1 calculation for the initial guess (currently only works with sparse matrix code). Guess=(AM1,Always) causes later steps in a geometry optimization to generate a new guess at each point and compare the energies with the density from the old point and the new guess and take the best one.
Core
Requests that the core Hamiltonian be diagonalized to form the initial guess. This is the default for AM1, PM3, PM3MM, PM6, and PDDG.
Orbital-Related Options
Permute
Read in a permutation of orbitals in the initial guess. The numbers of the generated guess orbitals are given in the order in which they should be used in the SCF. Ranges (e.g. 7–12) can be used, and all orbitals not listed are put in after the listed orbitals in their original order. Separate permutation lists for α and β orbitals must be specified (each list separated by a blank line) for open shell systems.
Alter
Indicates that the orbitals selected for occupation in the Hartree-Fock wavefunction should not be those of lowest energy. Normally, the occupied orbitals are selected as those with lowest eigenvalues for the one-electron Hamiltonian used in the initial guess programs. The alteration sections consist of a set of transpositions indicating that one of these occupied orbitals is to be replaced by one of the other (virtual) orbitals. Each such transposition is on a separate line and has two integers N1 and N2 (free format, separated by spaces or a comma as usual) indicating that orbital N1 is to be swapped with orbital N2. The list of orbital transpositions is terminated by the blank line at the end of the input section.
For UHF calculations, two such orbital alteration sections are required, the first specifying transpositions of α orbitals and the second specifying transpositions of β orbitals. Both sections are always required. Thus, even if only α transpositions are needed, the β section is required even though it is empty (and vice-versa). The second blank line to indicate an empty β section must be included.
Mix
Requests that the HOMO and LUMO be mixed so as to destroy α-β and spatial symmetries. This is useful in producing UHF wavefunctions for singlet states. Mixing orbitals is only done by default when generating a complex initial guess. NoMix says to not mix orbitals.
DensityMix[=N]
Whether to mix occupied and virtual orbital contributions in forming the initial guess density. N defaults to -3 (use Huckel eigenvalues to decide which orbitals to mix).
CopyChk
Copy the MOs and density unedited, unprojected and unchecked from the checkpoint file for the guess. Useful when importing from matrix element files.
Biorthogonalize
For an unrestricted guess, biorthogonalize the α and β MOs to maximally pair electrons of opposite spin. Done automatically when a UHF wavefunction is read-in as a guess for ROHF. This option is also useful in the combination of Guess=(BiOrth,Read,Only,Save) to replace the canonical UHF orbitals with ones which have alpha and beta orbitals matched up as much as possible for viewing in GaussView or other visualization packages.
NaturalOrbitals
Include natural orbitals in the checkpoint file. This must be accomplished via a separate job step specifying this option as well as Only and Save. See the discussion of the Population keyword for details.
Procedural Options
Only
Guess=Only functions as a calculation type keyword and requests that the calculation terminate once the initial guess is computed and printed. Note that the amount of orbital information that is printed is controlled by the Pop keyword. Guess=Only may not be used with the MOPAC-based semi-empirical methods (INDO, CNDO, MNDO, MINDO3).
This option is useful in preliminary runs to check if configuration alteration is necessary. For example, Guess=Only may be specified with CASSCF in order to obtain information on the number of CI configurations in the CAS active space (as well as the initial orbitals).
Guess=(Only,Read) may also be used to produce population and other post-calculation analyses from the data in a checkpoint file. For example, these options alone will produce a population analysis using the wavefunction in the checkpoint file. Guess(Only,Read) Prop will cause electrostatic properties to be calculated using the wavefunction in the checkpoint file.
Always
Requests that a new initial guess be generated at each point of an optimization. By default, the SCF results from the last point are used for the guess at the next point.
Fragment=N
Generate a guess built from fragment guesses or SCF solutions. The guess is saved on the checkpoint file regardless of if it is used in the current calculation, allowing it to be used in subsequent ones. This option is usually combined with Guess=Only to generate a guess from fragment guess orbitals (otherwise a full SCF calculation is performed for each fragment). Assignment of atoms to fragments and fragment charges and multiplicities are specified as described in Molecule Specifications, except that a negative spin multiplicity means that the unpaired orbitals for the fragment are to become β spin orbitals in the combined set.
Local
Requests that orbitals be localized using the Boys method [Boys60]. Occupied and virtual orbitals are localized separately, and the irreducible representations (after possible merging using LowSymm or NoSymm) are not mixed. Localized orbital analysis of a converged SCF wavefunction may then be done using a second job step, which includes Guess(Read,Local,Only) and Pop=Full in its route section.
Sparse
Perform a sparse SE calculation for the initial guess. This option may be useful for very, very large HF or DFT calculations using the sparse matrix facility.
Extra
Do an extra, new initial guess when reading orbitals from the RWF (i.e., during geometry optimizations). By default, this is done if the default Harris guess is allowed, no alteration of configuration was requested, and the optimization did not take a small step as flagged by variable 4 in ILSW. Use NoExtra to disable this feature.
Fock
Reuse Fock matrices rather than orbitals when reading from previous results from the read-write or checkpoint files. This is the default for periodic boundary conditions calculations if Guess=Alter is not specified. NoFock disables this behavior, and it is the default for non-PBC calculations.
Guess Reading/Saving Options
Read
Requests that the initial guess be read from the checkpoint file (Guess=Read is often specified along with Geom=Checkpoint). This option may be combined with Alter, in which case the orbitals are read from the checkpoint file, projected onto the current basis set, and then the specified alterations are made. Checkpoint is a synonym for Read.
TCheck
Read SCF restart information from the checkpoint file if possible; otherwise, read the initial guess. Designed for multistep jobs so that the SCF can be restarted in a later step when possible or it can be run with Guess=Read otherwise. This option implies CIS=Read or TD=Read when applicable/appropriate.Alpha
Use alpha orbitals for both alpha and beta guess during Guess=Read.
Translate
Translate requests that the coordinates of the atoms used to produce a guess, which is read in, be translated to the current atomic coordinates. This is the default. It may fail in unusual cases, such as when a wavefunction is used as a guess for a system with a different stoichiometry, in which case Guess=NoTranslate should be specified.
Cards
Specifies that after the initial guess is generated, some or all of the orbitals will be replaced with ones read from the input stream. This option can be used to read a complete initial guess from the input stream by replacing every orbital. The replacement orbitals are placed in the input section following the guess alteration commands, if any. For UHF, there are separate α and β replacement orbital input sections.
The replacement orbitals input section (the α replacement orbitals section for UHF) begins with a line specifying the Fortran format with which to read the replacement orbital input, enclosed in parentheses. For example: (4E20.8). The remainder of the section contains one or more instances of the following:
IVec Orbital to replace (0=end, -1=replace all orbitals in order). (A(I,IVec),I=1,N) New orbital in the format specified in the first line.
The format for the line containing IVec is Fortran I5. The β orbital replacement section for UHF calculations differs only in that it omits the initial format specification line. See the examples section for sample replacement orbital input.
Input
Read a line from the input file containing the name of a checkpoint file. Instead of a filename, you can specify one of the following keywords:
- generate says to generate a guess as usual.
- read and chk say to retrieve data from the job’s checkpoint file. They are thus equivalent to Guess=Read, and use the checkpoint file specified to %Chk or %OldChk as usual.
- none says to skip generating a guess at all.
- dummy generates a dummy guess that does not require any work.
See the examples section for information about using this option with ONIOM.
Restart
Start a new SCF calculation using the restart data saved on the checkpoint file from a job with a different geometry and/or basis set. (The option name is a bit of a misnomer.)
Save
Save the generated initial guess back into the checkpoint file at the conclusion of a Guess=Only run. This option is useful for saving localized orbitals.
Print the initial guess.
Symmetry-Related Options
LowSymm
Requests that irreducible representations of the molecular point group be combined in the symmetry information used in the N3 steps in the SCF to allow lowered symmetry of the wavefunction. This enables the orbitals (and possibly but not necessarily the total wavefunction) to have lower symmetry than the full molecular point group. This option is available only for GVB calculations, where it is often necessary for calculations on symmetric systems (see the discussion of the GVB keyword below for an example using this option).
The option expects a single line of input (in the format 16I2) giving the numbers of the irreducible representations to combine, with the new groups separated by 0; the list itself must be terminated by a 9. The numbers correspond to the order in which the representations are listed by Link 301 in the output file (see the examples subsection below).
Since this input section is always exactly one line long, it is not terminated by a blank line. Note that irreducible representations are combined before orbital localization is done and that localized orbitals retain whatever symmetry is kept. Guess=NoSymm removes all orbital symmetry constraints without reading any input.
NoSymm
Requests that all orbital symmetry constraints be lifted. Synonymous with SCF=NoSymm and Symm=NoSCF.
ForceAbelianSymmetry
Force the initial guess orbitals to transform according to irreducible representations of the Abelian point group. NoForceAbelianSymmetry is the default.
Valid Option Combinations
Only reasonable combinations are valid. For example Guess=(Always,Alter) and Guess=(Read,Alter) work as expected (in the former case, alterations are read once and the same interchanges are applied at each geometry). Conversely, Guess=(Always,Read) is contradictory and will lead to unpredictable results. Refer to the input sections order table at the beginning of this chapter to determine the ordering of the input sections for combinations of options like Guess=(Cards,Alter).
Transposing Two Orbitals with Guess=Alter
This example finds the UHF/STO-3G structure of the 2A1 excited state of the amino radical. First, a Guess=Only calculation is run to determine whether any alter (reordering) instructions are needed to obtain the desired electronic state. The HF/STO-3G theoretical model is used by default:
# Guess=Only Test Amino radical test of initial guess 0 2 n h 1 nh h 1 nh 2 hnh nh 1.03 hnh 120.0
Here is the orbital symmetry summary output from the job, which comes immediately before the population analysis in the output:
Initial guess orbital symmetries:
Alpha Orbitals:
Occupied (A1) (A1) (B2) (A1) (B1)
Virtual (A1) (B2)
Beta Orbitals:
Occupied (A1) (A1) (B2) (A1)
Virtual (B1) (A1) (B2)
The electronic state of the initial guess is 2-B1.
Initial guess = 0.0000 = 0.0000 = 0.5000 = 0.7500 S= 0.5000
Since a doublet state is involved, α and β orbitals are given separately. From the orbital symmetries, the electron configuration in the initial guess is 2a12a12b22a11b1, yielding a 2B1 wavefunction. This is indeed the ground state of NH2. The expectation value of S2 for the unrestricted initial guess is printed, in this case, that corresponding to a pure doublet: 0.75.
Since we want to model the 2A1 excited state of the amino radical, we will need to alter this initial orbital configuration: a β electron must be moved from orbital 4 to orbital 5 (the electron configuration is then 2a12a12b22b11a1).
Guess=Alter may be used to accomplish this. Here is the input for the geometry optimization:
# UHF/6-31G(d) Opt Guess=Alter Amino radical: HF/6-31G(d) Structure of 2-A1 state 0 2 N H 1 nh H 1 nh 2 hnh Variables: nh 1.03 hnh 120.0 Blank line ends the molecule specification section. Blank line ends the alpha section (empty in this case). 4 5 Transpose beta orbitals 4 and 5. End of the beta alteration section.
Note that an extra blank line—line 14—is necessary to indicate an empty α alteration section. The final two lines then constitute the β alteration section.
The initial guess program prints a list of orbitals that were interchanged as a result of the Alter option:
Harris functional with IExCor= 205 diagonalized for initial guess.
…
No Alpha orbitals switched.
Pairs of Beta orbitals switched:
4 5
Initial guess orbital symmetries:
Alpha Orbitals:
Occupied (A1) (A1) (B2) (A1) (B1)
Virtual (A1) (B2) (B1) (A1) (B2) (A1) (B2) (A1) (A2) (A1)
(B1) (A1) (B2) (A1)
Beta Orbitals:
Occupied (A1) (A1) (B2) (B1)
Virtual (A1) (A1) (B2) (B1) (A1) (B2) (A1) (B2) (A1) (A2)
(A1) (B1) (A1) (B2) (A1)
The electronic state of the initial guess is 2-A1.
Initial guess = 0.0000 = 0.0000 = 0.5000 = 0.7500 S= 0.5000
The SCF calculation follows, and the energy and S2 eigenvalue for the UHF wavefunction are printed. The S2 value which results if contamination of the wavefunction from the next possible spin multiplicity (quartets for doublets, quintets for triplets, etc.) is removed is also printed:
SCF Done: E(UHF) = -55.4915172451 A.U. after 12 cycles
Convg = 0.2693D-08 -V/T = 2.0038
= 0.0000 = 0.0000 = 0.5000 = 0.7534 S= 0.5017
= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation 0.7534, after 0.7500
Although this calculation does in fact converge correctly to 2A1 state, it sometimes happens that the order of orbital symmetries switches during the course of the SCF iterations. If the orbital symmetries of the final wavefunction are different from those in the initial guess (whether or not you are using Guess=Alter), we recommend using the direct minimization routine, specified with the SCF=QC or SCF=DM keywords, which usually holds symmetry from one iteration to the next.
Reordering Orbitals with Guess=Permute
This option is often is the easiest way to perform a complex modification of the initial guess, as in this example:
# CASSCF/6-31G(d,p) Opt Guess=Permute Pop=Reg CAS job 0 1 molecule specification 1-60 65 63 64 66 68 67 61-62 69 Specify new orbital ordering.
Here, we have rearranged orbitals 61–68. Listing the final orbital (69) is not really necessary, but it helps to make the input easier to understand for humans.
Specifying Calculation-Specific Guesses for ONIOM with Guess=Input
The Guess=Input is useful for ONIOM jobs where the SCF is hard to converge for some of the sub-calculations at the initial geometry. Convergence can be accomplished using a procedure like the following:
Step 1. Use ONIOM=OnlyInput to print out input files for the individual calculations.
Step 2. Converge the SCF for the problem step separately (for example, using Stable=Opt).
Step 3. After verifying that the correct state has been found, run an ONIOM calculation with Guess=Input to retrieve the correct initial guesses for the individual states. In the ONIOM input file, you will then have a line for the guess for each calculation (i.e., 3 lines for a 2-layer MO:MO), where each line can be a checkpoint file name or read (for Guess=Read from the regular checkpoint file) or generate (generate a new initial guess), as appropriate. For example, if only the high level calculation on the model system was problematic, the input would be:
| generate | Low level calculation on the real system. |
| Blank line. | |
| guess_for_high.chk | High level calculation on the model system. |
| Blank line. | |
| generate | Low level calculation on the model system. |
Reading in Orbitals with Guess=Cards
Some or all of the orbitals may be replaced, after the initial guess is generated, using Guess=Cards. Here is some sample input for this option, which replaces orbitals 1 and 4 (note that the format for the third and following lines is specified in line 1):
(3E20.8)
1
0.5809834509E+00 0.4612416518E+00 -0.6437319952E-04
0.1724432549E-02 0.1282235396E-14 0.5417658499E-13
0.1639966912E-02 -0.9146282229E-15 -0.6407549694E-13
-0.4538843604E-03 0.6038992958E-04 -0.1131035485E-03
0.6038992969E-04 -0.1131035471E-03
4
0.7700779642E-13 0.1240395916E-12 -0.3110890228E-12
-0.4479190461E-12 -0.1478805861E-13 0.5807753928E+00
0.6441113412E-12 -0.3119296374E-14 0.1554735923E+00
-0.1190754528E-11 0.2567325943E+00 0.1459733219E+00
-0.2567325943E+00 -0.1459733219E+00
0
An orbital number of zero ends the replacement orbital input.
Antiferromagnetic Coupling
Here is an example of a fragment guess job. The first step generates guesses for each one of the fragments and then combines them together into a guess for the SCF on the full molecule. The second step reads in that guess for a subsequent calculation. The charge and multiplicity line has the overall molecular charge and multiplicity first, followed by those for each fragment.
This example is Fe2S2 + 4 S-R ligands, where R is a phenyl group. It is a singlet with an overall charge of -2. The guess does each bare sulfur as S(2-) closed-shell (fragments 2 and 4), the two irons as Fe(3+) sextets antiferromagnetically coupled (fragments 1 and 3, with 1 being alpha-spin and 3 beta-spin), and each of the S-R fragments (5-8) as closed-shell singlet anions.
Fragment Assignments
%chk=FragGuess %mem=64mw #P UBP86/6-311G(d) Guess=(Fragment=8,Only) Pop=None Fe2S2 cluster with phenylthiolates. Step 1: Generate fragment guess -2,1 3,6 -2,1 3,-6 -2,1 -1,1 -1,1 -1,1 -1,1 H(Fragment=7) 23.5010 2.2873 8.5744 S(Fragment=2) 14.8495 1.1490 7.0431 Fe(Fragment=3) 17.0430 1.0091 7.0068 S(Fragment=4) 17.4565 -1.1490 7.0431 S(Fragment=5) 14.3762 -2.1581 8.7983 C(Fragment=5) 12.5993 -2.1848 8.6878 H(Fragment=5) 12.3743 -3.6513 10.1678 C(Fragment=5) 10.4994 -3.1122 9.4309 H(Fragment=5) 9.9929 -3.7579 10.0022 C(Fragment=5) 9.8049 -2.2791 8.5639 H(Fragment=5) 8.8050 -2.2873 8.5744 C(Fragment=5) 10.4833 -1.4146 7.6615 H(Fragment=5) 9.9730 -0.8525 7.0106 S(Fragment=8) 14.3794 -1.8091 5.0446 C(Fragment=5) 11.9048 -1.3675 7.7057 H(Fragment=5) 12.4158 -0.7843 7.0743 C(Fragment=6) 17.2999 3.4265 4.6624 C(Fragment=6) 16.6376 4.1967 5.6090 H(Fragment=6) 16.5022 3.8494 6.5369 C(Fragment=6) 16.1530 5.4856 5.2463 H(Fragment=6) 15.6665 6.0472 5.9155 C(Fragment=6) 16.3468 5.9257 4.0431 H(Fragment=6) 16.0236 6.8408 3.8020 C(Fragment=6) 17.0091 5.1398 3.0522 H(Fragment=6) 17.1330 5.4944 2.1254 C(Fragment=6) 17.4775 3.8823 3.3884 H(Fragment=6) 17.9400 3.3149 2.7071 S(Fragment=7) 17.9298 2.1581 8.7983 C(Fragment=7) 19.7067 2.1848 8.6878 C(Fragment=7) 20.4174 3.0650 9.5194 H(Fragment=7) 19.9317 3.6513 10.1678 C(Fragment=7) 21.8066 3.1122 9.4309 H(Fragment=7) 22.3132 3.7579 10.0022 C(Fragment=7) 22.5011 2.2791 8.5639 C(Fragment=7) 21.8227 1.4146 7.6615 H(Fragment=7) 22.3330 0.8525 7.0106 C(Fragment=7) 20.4012 1.3675 7.7057 H(Fragment=7) 19.8902 0.7843 7.0743 C(Fragment=8) 15.0061 -3.4265 4.6624 Fe(Fragment=1) 15.2630 -1.0091 7.0068 C(Fragment=8) 15.6684 -4.1967 5.6090 H(Fragment=8) 15.8038 -3.8494 6.5369 C(Fragment=8) 16.1530 -5.4856 5.2463 H(Fragment=8) 16.6395 -6.0472 5.9155 C(Fragment=5) 11.8886 -3.0650 9.5194 C(Fragment=8) 15.9592 -5.9257 4.0431 H(Fragment=8) 16.2824 -6.8408 3.8020 C(Fragment=8) 15.2969 -5.1398 3.0522 H(Fragment=8) 15.1730 -5.4944 2.1254 C(Fragment=8) 14.8285 -3.8823 3.3884 H(Fragment=8) 14.3660 -3.3149 2.7071 S(Fragment=6) 17.9266 1.8091 5.0446 --Link1-- %chk=FragGuess %mem=64mw #P UBP86/6-311G*/Auto Guess=Read Geom=AllCheck …
Last updated on: 11 December 2024. [G16 Rev. C.01]
