Searching....

# 10.15.3 Fragment-Based Methods for Electronic Coupling

(July 14, 2022)

## 10.15.3.1 Approach based on absolutely localized molecular orbitals

One can use absolutely localized molecular orbitals (ALMOs, see Chapter 12) to construct charge-localized diabatic states directly from DFT calculations. The ALMOs on each fragment are expanded by the AO basis functions belonging to the same fragment alone, whose corresponding MO coefficient matrix is fragment block-diagonal. 599 Khaliullin R. Z., Head-Gordon M., Bell A. T.
J. Chem. Phys.
(2006), 124, pp. 204105.
In energy decomposition analysis methods, 598 Khaliullin R. Z. et al.
J. Phys. Chem. A
(2007), 111, pp. 8753.
, 515 Horn P. R., Mao Y., Head-Gordon M.
Phys. Chem. Chem. Phys.
(2016), 18, pp. 23067.
ALMOs are utilized to separate the effects of polarization and charge transfer in intermolecular binding, because they have the useful property that they do not allow for charge transfer between fragments under the Mulliken definition of charge population. Making use of this property, one can construct charge-localized diabats for hole and electron transfer. For example, considering the initial and final states of a hole transfer process, $|D^{+}A\rangle$ and $|DA^{+}\rangle$, the two diabats can be represented in the following form:

 $\displaystyle|{\psi_{a}}\rangle$ $\displaystyle=\frac{1}{\sqrt{(N-1)!}}\mathrm{det}\left\{\phi^{(a)}_{D1},\phi^{% (a)}_{D2},\dots,\phi^{(a)}_{Dn_{D}\mathrm{-}1}\phi^{(a)}_{A1},\phi^{(a)}_{A2},% \dots,\phi^{(a)}_{An_{A}}\right\}$ (10.135a) $\displaystyle|{\psi_{b}}\rangle$ $\displaystyle=\frac{1}{\sqrt{(N-1)!}}\mathrm{det}\left\{\phi^{(b)}_{D1},\phi^{% (b)}_{D2},\dots,\phi^{(b)}_{Dn_{D}}\phi^{(b)}_{A1},\phi^{(b)}_{A2},\dots,\phi^% {(b)}_{An_{A}\mathrm{-}1}\right\}$ (10.135b)

For systems where the donor and acceptor moieties are well-separated, one can construct the ALMO-based diabats by simply concatenating orbitals obtained from isolated fragment calculations: $D^{+}$ and $A$ for one diabat, and $D$ and $A^{+}$ for the other. The energy of each ALMO diabat can then be variationally optimized with respect to orbital rotations on fragment, using the SCFMI technique (see Section 12.4). 1151 Stoll H., Wagenblast G., Preuss H.
Theor. Chem. Acc.
(1980), 57, pp. 169.
, 379 Gianinetti E., Raimondi M., Tornaghi E.
Int. J. Quantum Chem.
(1996), 60, pp. 157.
, 599 Khaliullin R. Z., Head-Gordon M., Bell A. T.
J. Chem. Phys.
(2006), 124, pp. 204105.
These ALMO-based diabatic states are variationally optimized such that the associated nuclear forces can be easily computed. 778 Mao Y., Horn P. R., Head-Gordon M.
Phys. Chem. Chem. Phys.
(2017), 19, pp. 5944.
The mutual polarization of donor and acceptor moieties in the presence of each other is also taken into account.

To calculate the electronic coupling between two ALMO diabats, one should first construct the diabatic Hamiltonian in the ALMO state basis

 $\mathbf{H}^{\prime}=\begin{pmatrix}H^{\prime}_{aa}&H^{\prime}_{ab}\\ H^{\prime}_{ba}&H^{\prime}_{bb}\end{pmatrix}$ (10.136)

and then transform that into the Löwdin-orthogonalized basis

 $\mathbf{H}=\mathbf{S}^{-1/2}\mathbf{H}^{\prime}\mathbf{S}^{-1/2}$ (10.137)

whose off-diagonal element, $H_{ab}$, corresponds to the diabatic coupling to be evaluated. In the 2-state case, we have

 $H_{ab}=\frac{1}{1-S_{ab}^{2}}\left|H^{\prime}_{ab}-\frac{H^{\prime}_{aa}+H^{% \prime}_{bb}}{2}S_{ab}\right|$ (10.138)

which requires the overlap between two ALMO diabats and the diagonal and off-diagonal elements of $\mathbf{H}^{\prime}$. The interstate overlap is given by

 $S_{ab}=\langle\psi_{a}|\psi_{b}\rangle=\mathrm{det}[(\mathbf{C}^{(a)}_{\mathrm% {o}})^{\dagger}\mathbf{S}\mathbf{C}^{(b)}_{\mathrm{o}}].$ (10.139)

where $\mathbf{C}^{(a)}_{\mathrm{o}}$ and $\mathbf{C}^{(b)}_{\mathrm{o}}$ are MO coefficients for the occupied orbitals in diabats $|\psi_{a}\rangle$ and $|\psi_{b}\rangle$, respectively, and $\mathbf{S}$ is the AO overlap matrix.

The elements of the diabatic Hamiltonian matrix can be evaluated using the multi-state DFT (MSDFT) approach. 187 Cembran A. et al.
J. Chem. Theory Comput.
(2009), 5, pp. 2702.
, 1014 Ren H. et al.
J. Phys. Chem. Lett.
(2016), 7, pp. 2286.
, 782 Mao Y., Montoya-Castillo A., Markland T. E.
J. Chem. Phys.
(2019), 151, pp. 164114.
For the diagonal elements, it is straightforward to employ the KS energies of the two diabats:

 $H^{\prime}_{aa}=E_{a}^{\mathrm{KS}}[\mathbf{P}^{(a)}],\quad H^{\prime}_{bb}=E_% {b}^{\mathrm{KS}}[\mathbf{P}^{(b)}]$ (10.140)

where $\mathbf{P}^{(a)}$ and $\mathbf{P}^{(b)}$ are the one-electron density matrices associated with two ALMO states $|\psi_{a}\rangle$ and $|\psi_{b}\rangle$, respectively. The approximation for the off-diagonal element is theoretically more challenging. In the original MSDFT scheme, 187 Cembran A. et al.
J. Chem. Theory Comput.
(2009), 5, pp. 2702.
, 1014 Ren H. et al.
J. Phys. Chem. Lett.
(2016), 7, pp. 2286.

 $H^{\prime}_{ab}=S_{ab}\left[V_{\mathrm{nn}}+\mathbf{P}_{ab}\cdot\mathbf{h}+% \frac{1}{2}\mathbf{P}_{ab}\cdot\mathbf{II}\cdot\mathbf{P}_{ab}+\frac{1}{2}(% \Delta E_{a}^{\mathrm{c}}+\Delta E_{b}^{\mathrm{c}})\right]$ (10.141)

where $\mathbf{P}_{ab}$ is the one-particle transition density matrix between two ALMO states

 $\mathbf{P}_{ab}=\mathbf{C}_{\mathrm{o}}^{(a)}\left[(\mathbf{C}_{\mathrm{o}}^{(% b)})^{\dagger}\mathbf{S}\mathbf{C}_{\mathrm{o}}^{(a)}\right]^{-1}(\mathbf{C}_{% \mathrm{o}}^{(b)})^{\dagger}$ (10.142)

The first three terms on the right-hand side of Eq. (10.141) correspond to the contributions from nuclear repulsion, one-electron Hamiltonian (kinetic energy and nuclei-electron attraction), and full two-electron integrals (Coulomb and full HF exchange), which can be derived as in non-orthogonal CI. 1187 Thom A. J. W., Head-Gordon M.
J. Chem. Phys.
(2009), 131, pp. 124113.
The last term accounts for the contribution from exchange-correlation (XC) functional as a correction to the HF coupling, which is given by the average of the difference between the KS and HF energies calculated from the same one-electron density matrix for each diabat:

 $\displaystyle\Delta E_{a}^{\mathrm{c}}$ $\displaystyle=E^{\mathrm{KS}}_{a}[\mathbf{P}^{(a)}]-E^{\mathrm{HF}}_{a}[% \mathbf{P}^{(a)}]$ (10.143a) $\displaystyle\Delta E_{b}^{\mathrm{c}}$ $\displaystyle=E^{\mathrm{KS}}_{b}[\mathbf{P}^{(b)}]-E^{\mathrm{HF}}_{b}[% \mathbf{P}^{(b)}]\;.$ (10.143b)

This approach was denoted as ALMO(MSDFT) in Ref.  782 Mao Y., Montoya-Castillo A., Markland T. E.
J. Chem. Phys.
(2019), 151, pp. 164114.
and it was found to overestimate the electronic couplings for the tested hole and electron transfer systems. A modified approach, denoted as ALMO(MSDFT2), was proposed in Ref.  782 Mao Y., Montoya-Castillo A., Markland T. E.
J. Chem. Phys.
(2019), 151, pp. 164114.
, which evaluates the XC contribution using the XC energy of the symmetrized transition density matrix

 $H^{\prime}_{ab}=S_{ab}\left[V_{\mathrm{nn}}+\mathbf{P}_{ab}\cdot\mathbf{h}+% \frac{1}{2}\mathbf{P}_{ab}\cdot\mathbf{II}\cdot\mathbf{P}_{ab}+E_{\mathrm{xc}}% [\tilde{\mathbf{P}}_{ab}]\right]$ (10.144)

where

 $\tilde{\mathbf{P}}_{ab}=\frac{1}{2}(\mathbf{P}_{ab}+\mathbf{P}_{ba})\;.$ (10.145)

Note that in Eq. (10.144), $\mathbf{II}$ includes only Coulomb integrals and a fraction of exact exchange if hybrid functionals are employed.

According to the benchmark results in Ref.  782 Mao Y., Montoya-Castillo A., Markland T. E.
J. Chem. Phys.
(2019), 151, pp. 164114.
, ALMO(MSDFT2) shows better accuracy than the original MSDFT method for hole and electron transfer, and thus it is implemented as the default approach to compute electronic couplings between ALMO diabats in Q-Chem. We note that the results given by Eq. (10.144) may become inaccurate when the overlap between two states becomes near-singular, as

 $\boldsymbol{\sigma}_{ba}=(\mathbf{C}_{\mathrm{o}}^{(b)})^{\dagger}\mathbf{S}% \mathbf{C}_{\mathrm{o}}^{(a)}$ (10.146)

is inverted when constructing the transition density [Eq. (10.142)]. To circumvent this numerical issue, one can replace the inverse in Eq. (10.142) with the Penrose pseudo-inverse, which was suggested for a similar objective in Ref.  912 Pavanello M. et al.
J. Chem. Phys.
(2013), 138, pp. 054101.
.

## 10.15.3.2 Projection operator Diabatization (POD)

Besides ALMO-based diabatization method, other fragment-based diabatization methods are available in Q-Chem. The projection operator diabatization (POD) method 618 Kondov I. et al.
J. Phys. Chem. C
(2007), 111, pp. 11970.
starts from a standard KS-DFT calculation of the system and post-processes the converged Fock matrix. It first transforms the Fock matrix into the Löwdin-orthogonalized AO basis and then partitions that into the donor and acceptor blocks, assuming that these orthogonalized AO basis functions still retain their original fragment tags:

 $\tilde{\mathbf{F}}=\mathbf{S}^{-1/2}\mathbf{F}\mathbf{S}^{-1/2}=\begin{pmatrix% }\tilde{\mathbf{F}}_{dd}&\tilde{\mathbf{F}}_{da}\\ \tilde{\mathbf{F}}_{ad}&\tilde{\mathbf{F}}_{aa}\end{pmatrix}$ (10.147)

One then diagonalizes $\tilde{\mathbf{F}}_{dd}$ and $\tilde{\mathbf{F}}_{aa}$ separately

 $\boldsymbol{\epsilon}_{d}=\mathbf{D}_{d}^{\dagger}\tilde{\mathbf{F}}_{dd}% \mathbf{D}_{d},\ \boldsymbol{\epsilon}_{a}=\mathbf{D}_{a}^{\dagger}\tilde{% \mathbf{F}}_{aa}\mathbf{D}_{a},$ (10.148)

where the eigenvectors $\mathbf{D}_{d}$ and $\mathbf{D}_{a}$ define the single-particle “diabatic states”:

 \displaystyle\begin{aligned} \displaystyle\left|\bar{\varphi}^{(d)}_{p}\right% \rangle&\displaystyle=\sum_{\mu}\left|\tilde{\chi}^{(d)}_{\mu}\right\rangle(D_% {d})^{\mu}_{\ p}\\ \displaystyle\left|\bar{\varphi}^{(a)}_{p}\right\rangle&\displaystyle=\sum_{% \mu}\left|\tilde{\chi}^{(a)}_{\mu}\right\rangle(D_{a})^{\mu}_{\ p},\\ \end{aligned} (10.149)

and transforms the off-diagonal block of the Fock matrix into this diabatic basis

 $\bar{\mathbf{F}}_{da}=\mathbf{D}_{d}^{\dagger}\tilde{\mathbf{F}}_{da}\mathbf{D% }_{a}$ (10.150)

yielding

 $\bar{\mathbf{F}}=\begin{pmatrix}\boldsymbol{\epsilon}_{d}&\bar{\mathbf{F}}_{da% }\\ \bar{\mathbf{F}}_{ad}&\boldsymbol{\epsilon}_{a}\end{pmatrix}$ (10.151)

The couplings between these single-particle diabatic orbitals can then be directly read off from the elements of $\bar{\mathbf{F}}_{da}$.

The Q-Chem implementation of the POD method follows the description in Refs.  618 Kondov I. et al.
J. Phys. Chem. C
(2007), 111, pp. 11970.
and 1318 Yang C.-H., Yam C. Y., Wang H.
Phys. Chem. Chem. Phys.
(2018), 20, pp. 2571.
, where a closed-shell reference system is used to generate the Fock matrix to be processed, i.e., $\mathbf{F}$ in Eq. (10.147). By default, only the D(HOMO) - A(HOMO) coupling is calculated for the hole transfer cases, and the D(LUMO) - A(LUMO) coupling for the electron transfer cases. To calculate the couplings between multiple pairs of donor and acceptor orbitals, the user can set $rem variable POD_MULTI_PAIRS to TRUE and control the number of orbitals pairs through POD_WINDOW. See the instruction in Sec. 10.15.3.4. Because of the use of globally Löwdin-orthogonalized orbitals in Eq. (10.147), the diabatic orbitals created by POD cannot be strictly localized on fragments. This renders the POD results unstable with the change of employed AO basis sets: when larger basis sets are used, the mixing between AO basis functions on different fragments becomes stronger, and the resulting $H_{ab}$ decreases. To alleviate this problem, a revised POD method, which was named as “POD2”, was proposed by Ghan et al.. 376 Ghan S. et al. J. Chem. Theory Comput. (2020), 16, pp. 7431. It avoids the global Löwdin-orthogonalization of the AO basis; instead, it separately diagonalizes the the donor and acceptor blocks of the Fock matrix (in the original AO basis):  $\mathbf{F}_{dd}\mathbf{C}_{d}=\mathbf{S}_{dd}\mathbf{C}_{d}\boldsymbol{% \epsilon}_{d},\quad\mathbf{F}_{aa}\mathbf{C}_{a}=\mathbf{S}_{aa}\mathbf{C}_{a}% \boldsymbol{\epsilon}_{a}$ (10.152) The obtained diabatic MO coefficient matrix is fragment-block-diagonal in the AO basis:  $\mathbf{C}_{\text{diab}}=\begin{pmatrix}\mathbf{C}_{d}&\mathbf{0}\\ \mathbf{0}&\mathbf{C}_{a}\\ \end{pmatrix}$ (10.153) Transforming the AO Fock matrix into this diabatic MO basis, the D-D and A-A blocks of the resulting matrix are diagonal matrices:  $\bar{\mathbf{F}}_{\text{diab}}=\mathbf{C}_{\text{diab}}^{T}\mathbf{F}\mathbf{C% }_{\text{diab}}=\begin{pmatrix}\boldsymbol{\epsilon}_{d}&\bar{\mathbf{F}}_{da}% \\ \bar{\mathbf{F}}_{ad}&\boldsymbol{\epsilon}_{a}\\ \end{pmatrix}$ (10.154) Using the matrix elements in the off-diagonal block ($\bar{\mathbf{F}}_{da}$) directly would yield overestimated couplings since the diabatic MOs $\mathbf{C}_{d}$ and $\mathbf{C}_{a}$ are not orthogonal to each other. Therefore, a final orthogonalization step is required to obtain the diabatic coupling between a pair of orbitals that are located on the donor and acceptor, respectively. Denoting this pair of orbital as $\phi_{d}$ and $\phi_{a}$, one can construct the $2\times 2$ Hamiltonian and overlap matrices:  $\mathbf{H}=\begin{pmatrix}\epsilon_{d}&\bar{F}_{da}\\ \bar{F}_{ad}&\epsilon_{a}\\ \end{pmatrix},\quad\mathbf{S}=\begin{pmatrix}1&S_{da}\\ S_{ad}&1\\ \end{pmatrix}$ (10.155) Two orthogonalization schemes have been investigated by Ghan et al.. 376 Ghan S. et al. J. Chem. Theory Comput. (2020), 16, pp. 7431. The first approach performs a Löwdin orthgonalization on $\phi_{d}$ and $\phi_{a}$, which is denoted as POD2L. The resulting coupling between the orthogonalized diabatic orbitals are  $H_{da}^{\text{eff}}=\frac{1}{1-S_{da}^{2}}\left|\bar{F}_{da}-\frac{1}{2}(% \epsilon_{d}+\epsilon_{a})S_{da}\right|$ (10.156) The second approach employs the Gram-Schmidt orthogonalization, which keeps one of the two orbitals ($\phi_{d}$ or $\phi_{a}$) intact while ensures that the other is strictly orthogonal to it. This approach is denoted as POD2GS, and it might be better choice for asymmetric cases (e.g. surface and adsorbates) where one can choose to retain the orbital on the less sizable fragment. These two POD2 variants afford significantly improved accuracy over the original POD method, especially in terms of the robustness with regard to the use of extensive basis sets. ## 10.15.3.3 Fragment Orbital DFT (FODFT) Fragment orbital DFT (FODFT) 875 Oberhofer H., Blumberger J. Phys. Chem. Chem. Phys. (2012), 14, pp. 13846. , 1081 Senthilkumar K. et al. J. Chem. Phys. (2003), 119, pp. 9809. , 1067 Schober C., Reuter K., Oberhofer H. J. Chem. Phys. (2016), 144, pp. 054103. is an approach to compute the diabatic couplings for hole and electron transfer between fragments. There have been several different flavors of FODFT approaches developed in literature, and here we introduce the most recent variant by Schober et al. 1067 Schober C., Reuter K., Oberhofer H. J. Chem. Phys. (2016), 144, pp. 054103. Considering a hole transfer process $D^{+}+A\rightarrow D+A^{+}$ or an electron transfer process $D^{-}+A\rightarrow D+A^{-}$, where the donor ($D$) and acceptor ($A$) fragments have $n_{D}$ and $n_{A}$ electrons, respectively, the procedure is as follows: • Perform KS-DFT calculations for isolated donor and acceptor fragments; collect the converged fragment orbitals: $\{\phi_{D1},\phi_{D2},\dots,\phi_{Dn_{D}\pm 1}\}$ and $\{\phi_{A1},\phi_{A2},\dots,\phi_{An_{A}}\}$ • Löwdin-orthogonalize the occupied orbitals on two fragments. The reactant diabat ($D^{+}A$ or $D^{-}A$) can be represented as  $|\bar{\psi}_{a}\rangle=\frac{1}{\sqrt{(N-1)!}}\mathrm{det}\left\{\bar{\phi}_{D% 1},\bar{\phi}_{D2},\dots,\bar{\phi}_{Dn_{D}\pm 1}\bar{\phi}_{A1},\bar{\phi}_{A% 2},\dots,\bar{\phi}_{An_{A}}\right\}$ (10.157) where “$\bar{\phi}$” denotes Löwdin-orthogonalized orbitals, and $N=n_{D}+n_{A}$. Note that the lowest unoccupied orbital where the electron is transferring to, $\phi_{Dn_{D}}$ in the case of HT or $\phi_{An_{A}+1}$ in the case of ET, also needs to be made orthogonal to the space spanned by all occupied orbitals. • Construct the product diabat ($DA^{+}$ or $DA^{-}$), simply by moving the hole from $\bar{\phi}_{Dn_{D}}$ to $\bar{\phi}_{An_{A}}$ (HT), or the excess electron from $\bar{\phi}_{Dn_{D}+1}$ to $\bar{\phi}_{An_{A}+1}$ (ET)  $|\bar{\psi}_{b}\rangle=\frac{1}{\sqrt{(N-1)!}}\mathrm{det}\left\{\bar{\phi}_{D% 1},\bar{\phi}_{D2},\dots,\bar{\phi}_{Dn_{D}}\bar{\phi}_{A1},\bar{\phi}_{A2},% \dots,\bar{\phi}_{An_{A}\pm 1}\right\}$ (10.158) • Compute the electronic coupling between $|\bar{\psi}_{a}\rangle$ and $|\bar{\psi}_{b}\rangle$, which is approximated by the coupling of the orthogonalized fragment orbitals through the Kohn-Sham Fock operator (built from the reactant diabat)  $\displaystyle\langle\bar{\psi}_{a}|\hat{H}|\bar{\psi}_{b}\rangle\approx\begin{% cases}\langle\phi_{Dn_{D}}|\hat{f}_{\mathrm{KS}}|\phi_{An_{A}}\rangle,\ % \mathrm{HT}\\ \langle\phi_{Dn_{D}+1}|\hat{f}_{\mathrm{KS}}|\phi_{An_{A}+1}\rangle,\ \mathrm{% ET}\end{cases}$ (10.159) The approach described above is denoted as FODFT($\mathrm{2n-1}$)@$D^{+}A$ (HT) / FODFT($\mathrm{2n+1}$)@$D^{-}A$ (ET) 1067 Schober C., Reuter K., Oberhofer H. J. Chem. Phys. (2016), 144, pp. 054103. as the charged fragment is explicitly taken into account when preparing the fragment orbitals and the KS Fock matrix is built from $2n\mp 1$ occupied orbitals. Besides this, there are two other variants of FODFT: 1. 1. FODFT($2n$)@$DA$: 1081 Senthilkumar K. et al. J. Chem. Phys. (2003), 119, pp. 9809. fragment orbitals prepared with $D$ and $A$ both closed-shell; KS Fock operator constructed from $2n$ occupied orbitals 2. 2. FODFT($\mathrm{2n-1}$)@$DA$ (HT) / FODFT($\mathrm{2n+1}$)@$D^{-}A^{-}$ (ET): 875 Oberhofer H., Blumberger J. Phys. Chem. Chem. Phys. (2012), 14, pp. 13846. fragment orbitals prepared with the system having one excess electron ($DA$ for HT and $D^{-}A^{-}$ for ET), while one occupied orbital is removed when building the KS Fock operator According to the benchmark results, 1067 Schober C., Reuter K., Oberhofer H. J. Chem. Phys. (2016), 144, pp. 054103. , 782 Mao Y., Montoya-Castillo A., Markland T. E. J. Chem. Phys. (2019), 151, pp. 164114. FODFT($\mathrm{2n-1}$)@$D^{+}A$ (HT) / FODFT($\mathrm{2n+1}$)@$D^{-}A$ (ET) is the best-performing method, possibly because of its explicit account for charged fragments and consistent electron count in the preparation of fragment orbitals and in the construction of Fock matrix. One issue associated with the FODFT methods is that for asymmetric systems, the results would depend on how one chooses the initial and final states for an electron or hole transfer process (e.g. $D^{+}A$ vs. $DA^{+}$), especially for the two variants that build the Fock matrix with $2n\pm 1$ occupied orbitals. 782 Mao Y., Montoya-Castillo A., Markland T. E. J. Chem. Phys. (2019), 151, pp. 164114. The Q-Chem implementation of FODFT($\mathrm{2n-1}$)@$DA$ / FODFT($\mathrm{2n+1}$)@$D^{-}A^{-}$ automatically computes $H_{ab}$ in both ways and then reports the average, as it only requires an extra Fock matrix build. This, however, is not automatically done for FODFT($\mathrm{2n-1}$)@$D^{+}A$ / FODFT($\mathrm{2n+1}$)@$D^{-}A$. ## 10.15.3.4 Job control of fragment based diabatization methods POD, FODFT, and ALMO(MSDFT) calculations in Q-Chem require specification of fragments in the$molecule section (see Sec. 12.2). For ALMO(MSDFT) calculations, one also needs to specify the charge and multiplicity of each fragment in each diabatic state in the $almo_coupling section, where two hyphens indicate the separation of different diabats: $almo_coupling
charge_frag_1     mult_frag_1        !diabat 1
charge_frag_2     mult_frag_2
--
charge_frag_1     mult_frag_1        !diabat 2
charge_frag_2     mult_frag_2
$end  The current implementation of FODFT is limited to hole transfer between the HOMOs of two fragments or electron transfer between the LUMOs, and the current simplementation of ALMO(MSDFT) is limited to ground state electron or hole transfer involving two states. FRAG_DIABAT_METHOD FRAG_DIABAT_METHOD Specify fragment based diabatization method TYPE: STRING DEFAULT: NONE OPTIONS: ALMO_MSDFT Perform ALMO(MSDFT) diabatization POD Perform projection operator diabatization (the original method) POD2_L Perform POD2 with Löwdin orthogonalization POD2_GS Perform POD2 with Grad-Schmidt orthogonalization ESID The energy-split-in-dimer method, 1210 Valeev E. F. et al. J. Am. Chem. Soc. (2006), 128, pp. 9882. which is equivalent to the FMO approach introduced in Section 10.15.2.5 FODFT Calculate electronic coupling using fragment orbital DFT RECOMMENDATION: NONE FRAG_DIABAT_DOHT FRAG_DIABAT_DOHT Specify whether hole or electron transfer is considered TYPE: BOOLEAN DEFAULT: TRUE OPTIONS: TRUE Do hole transfer FALSE Do electron transfer RECOMMENDATION: Need to be specified for POD and FODFT calculations FRAG_DIABAT_PRINT FRAG_DIABAT_PRINT Specify the print level for fragment based diabatization calculations TYPE: INTEGER DEFAULT: 0 OPTIONS: 0 No additional prints $\geq 1$ Print additional details RECOMMENDATION: Use the default unless debug information is needed MSDFT_METHOD MSDFT_METHOD Specify the scheme for ALMO(MSDFT) TYPE: INTEGER DEFAULT: 2 OPTIONS: 1 The original MSDFT scheme [Eq. (10.141)] 2 The ALMO(MSDFT2) approach [Eq. (10.144)] RECOMMENDATION: Use the default method. Note that the method will be automatically reset to 1 if a meta-GGA functional is requested. MSDFT_PINV_THRESH MSDFT_PINV_THRESH Set the threshold for pseudo-inverse of the interstate overlap TYPE: INTEGER DEFAULT: 4 OPTIONS: $n$ Set the threshold to 10${}^{-n}$ RECOMMENDATION: Use the default value POD_MULTI_PAIRS POD_MULTI_PAIRS Calculate the couplings between multiple pairs of donor and acceptor orbitals in POD TYPE: BOOLEAN DEFAULT: FALSE OPTIONS: TRUE Calculate the couplings between multiple pairs of orbitals FALSE Only calculate the D(HOMO)–A(HOMO) coupling (for HT) or D(LUMO)–A(LUMO) coupling (for ET) RECOMMENDATION: None POD_WINDOW POD_WINDOW Specify the number of donor and acceptor orbitals when couplings between multiple pairs are requested TYPE: INTEGER DEFAULT: 5 OPTIONS: $n$ Including $n$ frontier occupied orbitals (from $\mathrm{HOMO}-n+1$ to HOMO) and $n$ frontier virtual orbitals (from LUMO to $\mathrm{LUMO}+n-1$) for both donor and acceptor RECOMMENDATION: None FODFT_METHOD FODFT_METHOD Specify the flavor of FODFT method TYPE: INTEGER DEFAULT: 1 OPTIONS: 1 FODFT($\mathrm{2n-1}$)@$D^{+}A$ (HT) / FODFT($\mathrm{2n+1}$)@$D^{-}A$ (ET) 2 FODFT($\mathrm{2n}$)@$DA$ 3 FODFT($\mathrm{2n-1}$)@$DA$ (HT) / FODFT($\mathrm{2n+1}$)@$D^{-}A^{-}$ (ET) RECOMMENDATION: The default approach shows the best overall performance FODFT_DONOR FODFT_DONOR Specify the donor fragment in FODFT calculation TYPE: INTEGER DEFAULT: 1 OPTIONS: 1 First fragment as donor 2 Second fragment as donor RECOMMENDATION: With FODFT_METHOD = 1, the charged fragment needs to be the donor fragment Example 10.55 ALMO(MSDFT2) calculation for hole transfer in ethylene dimer $molecule
1 2
--
1 2
C      0.000000    0.000000    0.000000
C      1.332000    0.000000    0.000000
H     -0.574301    0.000000   -0.928785
H     -0.574301    0.000000    0.928785
H      1.906301    0.000000    0.928785
H      1.906301    0.000000   -0.928785
--
0 1
C     -0.000000    4.000000    0.000000
C      1.332000    4.000000   -0.000000
H     -0.574301    4.000000    0.928785
H     -0.574301    4.000000   -0.928785
H      1.906301    4.000000   -0.928785
H      1.906301    4.000000    0.928785
$end$rem
METHOD               PBE0
BASIS                6-31+G(D)
UNRESTRICTED         TRUE
THRESH               14
SCF_CONVERGENCE      8
SYMMETRY             FALSE
SYM_IGNORE           TRUE
SCFMI_MODE           1
FRGM_METHOD          STOLL
FRAG_DIABAT_METHOD   ALMO_MSDFT
$end$almo_coupling
1  2
0  1
--
0  1
1  2
$end  Example 10.56 POD diabatization method for hole transfer in ethylene dimer. FRAG_DIABAT_METHOD can be set to POD2_L or POD2_GS for POD2 diabatization methods. $molecule
0 1
--
0 1
C      0.000000    0.000000    0.000000
C      1.332000    0.000000    0.000000
H     -0.574301    0.000000   -0.928785
H     -0.574301    0.000000    0.928785
H      1.906301    0.000000    0.928785
H      1.906301    0.000000   -0.928785
--
0 1
C     -0.000000    4.000000    0.000000
C      1.332000    4.000000   -0.000000
H     -0.574301    4.000000    0.928785
H     -0.574301    4.000000   -0.928785
H      1.906301    4.000000   -0.928785
H      1.906301    4.000000    0.928785
$end$rem
METHOD               lrc-wpbeh
BASIS                6-31+g(d)
FRAG_DIABAT_METHOD   pod
SCF_CONVERGENCE      8
THRESH               14
SYMMETRY             false
SYM_IGNORE           true
$end  Example 10.57 FODFT($\mathrm{2n-1}$)@$D^{+}A$ calculation for hole transfer in ethylene dimer $molecule
1 2
--
1 2
C      0.000000    0.000000    0.000000
C      1.332000    0.000000    0.000000
H     -0.574301    0.000000   -0.928785
H     -0.574301    0.000000    0.928785
H      1.906301    0.000000    0.928785
H      1.906301    0.000000   -0.928785
--
0 1
C     -0.000000    4.000000    0.000000
C      1.332000    4.000000   -0.000000
H     -0.574301    4.000000    0.928785
H     -0.574301    4.000000   -0.928785
H      1.906301    4.000000   -0.928785
H      1.906301    4.000000    0.928785
$end$rem
METHOD               wb97x-d
BASIS                6-31+g(d)
UNRESTRICTED         true
SCF_CONVERGENCE      8
THRESH               14
SYMMETRY             false
SYM_IGNORE           true
FRAG_DIABAT_METHOD   fodft
FODFT_METHOD         1
\$end