12.7 Second-Generation ALMO-EDA Method

12.7.6 ALMO-EDA with non-aufbau Electronic Configurations

(May 16, 2021)

Q-Chem also supports ALMO-EDA calculations with one or multiple fragments in non-aufbau electronic configurations. This method is particularly useful for cases where the energetically most stable electronic state of a given fragment changes upon the formation of intermolecular complex, which is common for systems involving open-shell species such as carbenes, transition metal cations, etc. For example, the lowest-energy electronic configuration of an isolated Ba+ radical cation is 6s1; however, when it is ligated with CO, the complex is more stable when the unpaired electron of Ba+ is promoted from 6s1 to 5d1, which is thus in a non-aufbau state relative to the isolated Ba+ cation. A sensible strategy to study the formation of such a complex is to perform an ALMO-EDA calculation with Ba+ in the 5d1 electronic configuration, and then evaluate the energy difference between Ba+(5d1) and Ba+(6s1) and interpret that as the monomer “electronic preparation” energy. To ensure that the system stays in a given non-aufbau electronic configuration throughout, the Maximum Overlap Method 352 Gilbert A. T. B., Besley N. A., Gill P. M. W.
J. Phys. Chem. A
(2008), 112, pp. 13164.
(MOM, see Sec. 4.5.11 for details) is applied to almost all stages in ALMO-EDA, including (i) isolated fragment calculations, (ii) calculation for the polarized wavefunction, and (iii) full SCF calculation for the whole system. Among them, (i) and (iii) are standard SCF calculations using MOM, while (ii) involves the application of MOM to an SCFMI calculation, 720 Mao Y., Montoya-Castillo A., Markland T. E.
J. Chem. Phys.
(2020), 153, pp. 244111.
which is compatible with the use of DIIS algorithm to solve locally projected SCF equations (with FRGM_METHOD = STOLL or GIA, see Section 12.4).

Specifically, an ALMO-EDA calculation for a complex within a non-aufbau electronic configuration (assuming only one of the fragments is excited) comprises the following steps:

  • Perform ground-state SCF calculations for all fragments

  • Calculate the non-aufbau electronic configuration for the fragment as specified in the input using MOM; construct the frozen state for the system with that one fragment in the non-aufbau state and the rest in ground state

  • Starting from the frozen state, perform an SCFMI calculation with MOM to obtain the polarized state within the non-aufbau electronic configuration

  • Perform a full SCF calculation with MOM, starting from the polarized wavefunction in the non-aufbau state

Such a calculation requires EDA2_MOM = TRUE in addition to the setup of standard EDA2 jobs. The non-aufbau electronic state is specified through the $scfmi_mom input section, which has the following format:

frag_idx1   X_1   Y_1   spin_1
frag_idx2   X_2   Y_2   spin_2

In each row, the first entry specifies the index of the non-aufbau fragment (starting from 1). The second and third entries specify the non-aufbau electronic configuration, which is prepared by promoting an electron from a given fragment’s HOMO-X to LUMO+Y once the ground-state SCF calculation of that fragment is finished. The last entry specifies the spin of the occupied and virtual orbitals that are being swapped (“a” for α electrons and “b” for β electrons, and α orbitals will be used by default if the fourth entry is left blank). For instance, for a HOMOLUMO excitation on the first fragment, this section would simply look like

¯1  0  0  a

Q-Chem 5 also allows one to perform SCFMI calculations for non-aufbau electronic configurations, 720 Mao Y., Montoya-Castillo A., Markland T. E.
J. Chem. Phys.
(2020), 153, pp. 244111.
without going through the entire ALMO-EDA procedure. To do that, one can simply add SCFMI_MOM = TRUE to the $rem setup of an SCFMI calculation (see Section 12.7.1). The non-aufbau electronic configuration can be specified through the $scfmi_mom section in the same way.

       Perform ALMO-EDA calculation with non-aufbau electronic configurations using MOM
       FALSE Standard ALMO-EDA calculation TRUE ALMO-EDA for non-aufbau states

Example 12.19  ALMO-EDA calculation for the [Ba(CO)]+ complex with Ba+ in the 5d1 electronic configuration at the B3LYP/def2-TZVP level. The $scfmi_mom section specifies that the unpaired α electron in the 6s orbital is promoted to one of the 5d orbitals.

1 2
1 2
Ba   0.0  0.0  0.0
0 1
C    0.0  0.0  2.734
O    0.0  0.0  3.876

   JOBTYPE           EDA    !by default EDA2 = 2 with CT analysis
   METHOD            B3LYP
   BASIS             DEF2-TZVP
   ECP               DEF2-ECP
   EDA2_MOM          TRUE

1 0 2

View output