A charge transfer involves a change in the electron numbers in a pair of molecular fragments. As an example, we will use the following reaction when necessary, and a generalization to other cases is straightforward:
where an extra electron is localized to the donor (D) initially, and it becomes localized to the acceptor (A) in the final state. The two-state secular equation for the initial and final electronic states can be written as
This is very close to an eigenvalue problem except for the non-orthogonality between the initial and final states. A standard eigenvalue form for Eq. (10.77) can be obtained by using the Löwdin transformation:
where the off-diagonal element of the effective Hamiltonian matrix represents the electronic coupling for the reaction, and it is defined by
In a general case where the initial and final states are not normalized, the electronic coupling is written as
Thus, in principle, can be obtained when the matrix elements for the Hamiltonian and the overlap matrix are calculated.
The direct coupling (DC) scheme calculates the electronic coupling values via Eq. (10.80), and it is widely used to calculate charge transfer coupling.Ohta:1986, Broo:1990, Farazdel:1990, Zhang:1997 In the DC scheme, the coupling matrix element is calculated directly using charge-localized determinants (the “diabatic states” in electron transfer literature). In electron transfer systems, it has been shown that such charge-localized states can be approximated by symmetry-broken unrestricted Hartree-Fock (UHF) solutions.Ohta:1986, Broo:1990, Newton:1991 The adiabatic eigenstates are assumed to be the symmetric and anti-symmetric linear combinations of the two symmetry-broken UHF solutions in a DC calculation. Therefore, DC couplings can be viewed as a result of two-configuration solutions that may recover the non-dynamical correlation.
The core of the DC method is based on the corresponding orbital transformationKing:1967 and a calculation for Slater’s determinants in and .Farazdel:1990, Zhang:1997 Unfortunately, the calculation of is not available for DFT method because a functional of the two densities and is unknown and there are no existing approximate forms for .Wu:2006c To calculate charge transfer coupling with DFT, we can use the CDFT-CI method (Section 5.13.3), the frontier molecular orbital (FMO) approachValeev:2006, You:2015a (Section 10.15.2.5) or a hybrid scheme – DC with CDFT wave functions (Section 10.15.2.4).
Let and be two single Slater-determinant wave functions for the initial and final states, and and be the spin-orbital sets, respectively:
Since the two sets of spin-orbitals are not orthogonal, the overlap matrix can be defined as:
We note that is not Hermitian in general since the molecular orbitals of the initial and final states are separately determined. To calculate the matrix elements and , two sets of new orthogonal spin-orbitals can be used by the corresponding orbital transformation.King:1967 In this approach, each set of spin-orbitals and are linearly transformed,
where and are the left-singular and right-singular matrices, respectively, in the singular value decomposition (SVD) of :
The overlap matrix in the new basis is now diagonal
The Hamiltonian for electrons in molecules are a sum of one-electron and two-electron operators. In the following, we derive the expressions for the one-electron operator and two-electron operator ,
where and , for the molecular Hamiltonian, are
The evaluation of matrix elements can now proceed:
In an atomic orbital basis set, , we can expand the molecular spin orbitals and ,
The one-electron terms, Eq. (10.92), can be expressed as
where and define a generalized density matrix, :
Similarly, the two-electron terms, Eq. (10.94), are
where and are generalized density matrices as defined in Eq. (10.99) except in is replaced by .
The - and -spin orbitals are treated explicitly. In terms of the spatial orbitals, the one- and two-electron contributions can be reduced to
It is important to obtain proper charge-localized initial and final states for the DC scheme, and this step determines the quality of the coupling values. Q-Chem provides three approaches to construct charge-localized states:
The “1+1” approach
Since the system consists of donor and acceptor molecules or fragments, with a charge being localized either donor or acceptor, it is intuitive to combine wave functions of individual donor and acceptor fragments to form a charge-localized wave function. We call this approach “1+1” since the zeroth order wave functions are composed of the HF wave functions of the two fragments.
For example, for the case shown in Example (10.76), we can use Q-Chem to calculate two HF wave functions: those of anionic donor and of neutral acceptor and they jointly form the initial state. For the final state, wave functions of neutral donor and anionic acceptor are used. Then the coupling value is calculated via Eq. (10.80).
$molecule -1 2 -- -1 2, 0 1 C 0.662489 0.000000 0.000000 H 1.227637 0.917083 0.000000 H 1.227637 -0.917083 0.000000 C -0.662489 0.000000 0.000000 H -1.227637 -0.917083 0.000000 H -1.227637 0.917083 0.000000 -- 0 1, -1 2 C 0.720595 0.000000 4.5 H 1.288664 0.921368 4.5 H 1.288664 -0.921368 4.5 C -0.720595 0.000000 4.5 H -1.288664 -0.921368 4.5 H -1.288664 0.921368 4.5 $end $rem JOBTYPE SP METHOD HF BASIS 6-31G(d) SCF_PRINT_FRGM FALSE SYM_IGNORE TRUE SCF_GUESS FRAGMO STS_DC TRUE $end
In the $molecule subsection, the first line is for the charge and multiplicity of the whole system. The following blocks are two inputs for the two molecular fragments (donor and acceptor). In each block the first line consists of the charge and spin multiplicity in the initial state of the corresponding fragment, a comma, then the charge and multiplicity in the final state. Next lines are nuclear species and their positions of the fragment. For example, in the above example, the first block indicates that the electron donor is a doublet ethylene anion initially, and it becomes a singlet neutral species in the final state. The second block is for another ethylene going from a singlet neutral molecule to a doublet anion.
Note that the last three $rem variables in this example, SYM_IGNORE, SCF_GUESS and STS_DC must be set to be the values as in the example in order to perform DC calculation with “1+1” charge-localized states. An additional $rem variable, SCF_PRINT_FRGM is included. When it is TRUE a detailed output for the fragment HF self-consistent field calculation is given.
The “relaxed” approach
In “1+1” approach, the intermolecular interaction is neglected in the initial and final states, and so the final electronic coupling can be underestimated. As a second approach, Q-Chem can use “1+1” wave function as an initial guess to look for the charge-localized wave function by further HF self-consistent field calculation. This approach would ‘relax’ the wave function constructed by “1+1” method and include the intermolecular interaction effects in the initial and final wave functions. However, this method may sometimes fail, leading to either convergence problems or a resulting HF wave function that cannot represent the desired charge-localized states. This is more likely to be a problem when calculations are performed with diffusive basis functions, or when the donor and acceptor molecules are very close to each other.
To perform relaxed DC calculation, set STS_DC = RELAX.
A hybrid scheme – constrained DFT charge-localized states
Constrained DFT (Section 5.13) can be used to obtain charge-localized states. It is recommended to set both charge and spin constraints in order to generate proper charge localization. To perform DC calculation with CDFT states, set SAVE_SUBSYSTEM = 10 and SAVE_SUBSYSTEM = 20 to save CDFT molecular orbitals in the first two jobs of a batch jobs, and then in the third job of the batch job, set SCF_GUESS = READ and STS_DC = TRUE to compute electronic coupling values.
$molecule -1 2 C 0.662489 0.000000 0.000000 H 1.227637 0.917083 0.000000 H 1.227637 -0.917083 0.000000 C -0.662489 0.000000 0.000000 H -1.227637 -0.917083 0.000000 H -1.227637 0.917083 0.000000 C 0.720595 0.000000 4.500000 H 1.288664 0.921368 4.500000 H 1.288664 -0.921368 4.500000 C -0.720595 0.000000 4.500000 H -1.288664 -0.921368 4.500000 H -1.288664 0.921368 4.500000 $end $rem METHOD = wb97xd BASIS = cc-pvdz CDFT = true SYM_IGNORE = true SAVE_SUBSYSTEM = 10 $end $cdft 1.0 1.0 1 6 1.0 1.0 1 6 s $end @@@ $molecule read $end $rem method = wb97xd basis = cc-pvdz cdft = true sym_ignore = true save_subsystem = 20 $end $cdft 1.0 1.0 7 12 1.0 1.0 7 12 s $end @@@ $molecule read $end $rem method = hf basis = cc-pvdz sym_ignore = true scf_guess = read sts_dc = true $end
The frontier molecular orbital (FMO) approach is often used with DFT to calculate ET coupling.Valeev:2006, You:2015a FMO coupling value is essentially an off-diagonal Kohn–Sham matrix element with the overlap effect accounted
where , with being the Kohn–Sham operator of the donor-acceptor system. is the Kohn–Sham frontier molecular orbital for the donor (acceptor) fragment, which represents one-particle scheme of a charge transfer process.
In this approach, computations are often performed separately in the two fragments, and the off-diagonal Kohn–Sham operator (and the overlap matrix) in the FMOs is subsequently calculated. To compute FMO couplings, Q-Chem has a setup that is similar to the “1+1” approach
$molecule 0 1 -- 0 1 C 0.662489 0.000000 0.000000 H 1.227637 0.917083 0.000000 H 1.227637 -0.917083 0.000000 C -0.662489 0.000000 0.000000 H -1.227637 -0.917083 0.000000 H -1.227637 0.917083 0.000000 -- 0 1 C 0.720595 0.000000 4.5 H 1.288664 0.921368 4.5 H 1.288664 -0.921368 4.5 C -0.720595 0.000000 4.5 H -1.288664 -0.921368 4.5 H -1.288664 0.921368 4.5 $end $rem method = lrcwpbe omega = 370 basis = dz* scf_print_frgm = true sym_ignore = true scf_guess = fragmo sts_dc = fock sts_trans_donor = 2-3 ! use HOMO, HOMO-1 and LUMO, LUMO+1, LUMO+2 of donor sts_trans_acceptor = 1-2 ! use HOMO and LUMO, LUMO+1 of acceptor $end $rem_frgm print_orbitals = 5 $end
Note: The FMOs are not always HOMO or LUMO of fragments. We can use STS_TRANS_DONOR (and STS_TRANS_ACCEPTOR) to select a range of occupied and virtual orbitals for FMO coupling calculations.