The GUGADRT program generates distict row table (DRT) used in the GUGACI in Direct CI calculations . Only DRT in active space are generated because the hole-particle symmetry is used in GUGACI [93, 94]. These DRT are used to evaluated the coupling coefficients by the Graphical Unitary Group Approach [72, 73, 74], for wavefunctions with at most two electrons excited from a set of reference configurations. The reference configurations can be specified as a list, where the occupation numbers are given for each active orbital (see below) in each reference configuration, or as a Full CI within the space defined by the active orbitals. In the GUGADRT and GUGACI the orbitals are classified as follows: Frozen, Inactive, Active, Secondary, and Deleted orbitals. Within each symmetry type, they follow this order. For the GUGADRT program, only the active orbitals are relevant.
Inactive: Inactive orbitals are doubly occupied in all reference configurations, but excitations out of this orbital space are allowed in the final CI wavefunction, i.e., they are correlated but have two electrons in all reference configurations. Since only single and double excitations are allowed, there can be no more than two holes in the active orbitals.
Active: Active orbitals are those which may have different occupation in different reference configurations.
188.8.131.52.1. Input files¶
Transformed one-electron integrals from MOTRA. Orbital information such as frozen, deleted orbitals will be read from this file.
184.108.40.206.2. Output files¶
This file contains the DRT that is needed in subsequent CI calculations.
This section describes the input to the GUGADRT program in the Molcas program system, with the program name:
The first four characters of the keywords are decoded and the rest are ignored.
Formally, there are no compulsory keyword. Obviously, some input must be given for a meaningful calculation.
The lines following this keyword are treated as title lines, until another keyword is encountered.
The spin degeneracy number, i.e. 2S+1. The value is read from the line following the keyword, in free format. The default value is 1, meaning a singlet wave function.
The number of electrons to be correlated in the CI calculation. The value is read from the line following the keyword, in free format. Note that this number should include the nr of electrons in inactive orbitals. An alternative input specification is NACTEL. Default: Twice nr of inactive orbitals.
The number of electrons in active orbitals in the reference configurations. The value is read from the line following the keyword, in free format. Note that this number includes only the of electrons in active orbitals. An alternative input specification is ELECTRONS. Default: Zero.
The number of inactive orbitals, i.e. orbitals that have occupation numbers of 2 in all reference configurations. Specified for each of the symmetries. The values are read from the line following the keyword, in free format.
The number of active orbitals, i.e. orbitals that have varying occupation numbers in the reference configurations. Specified for each of the symmetries. The values are read from the line following the keyword, in free format.
At least one of the Inactive or Active keywords must be present for a meaningful calculation. If one of them is left out, the default is 0 in all symmetries.
Specify selected reference configurations. The additional input that is required usually spans more than one line. The first line after the keyword contains the number of reference configurations, and the total number of active orbitals, and these two numbers are read by free format. Thereafter the input has one line per reference configuration, specifying the occupation number for each of the active orbitals, read by 80I1 format. Note that Reference and CIall are mutually exclusive.
Specify the selected symmetry type (the irrep) of the wave function as a number between 1 and 8 (see SEWARD). Default is 1, which always denote the totally symmetric irrep.
Use a Full CI within the subspace of the active orbitals as reference configurations. The symmetry of the wavefunction must be specified. Note that CIall and Reference are mutually exclusive. One of these two alternatives must be chosen for a meaningful calculation.
Printlevel of the program. Default printlevel (0) produces very little output. Printlevel 5 gives some information that may be of interest. The value is read from the line following the keyword, in free format.
220.127.116.11.2. Input example¶
&GUGADRT Title = CH2 molecule. Electrons = 8 Spin = 1 Inactive = 1 0 0 0 Active = 2 2 2 0 Symmetry = 1 Ciall