Calculation of electronic coupling matrix elements for ground and excited state electron transfer reactions: Comparison of the generalized Mulliken-Hush and block diagonalization methods

Journal of Chemical Physics

Volumn 106, Issue 22, 1997, Pages 9213-9226

  Cave, Robert J ^a   Newton, Marshall D ^b  

^a HARVEY MUDD COLLEGE   (United States)

^b BROOKHAVEN NATIONAL LABORATORY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

EID: 4243234418     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.474023     Document Type: Article

References (87)

1. (a) M. D. Newton, Chem. Rev. 91, 767 (1991);
2. (b) M. D. Newton and R. J. Cave, in Molecular Electronics, edited by J. Jortner and M. A. Ratner (Blackwell Science, in press).
3. R. A. Marcus and N. Sutin, Biochim. Biophys. Acta 811, 265 (1985).
4. I. R. Gould, R. H. Young, L. J. Mueller, A. C. Albrecht, and S. Farid, J. Am. Chem. Soc. 116, 3147, 8188 (1994).
5. (a) J. Halpern and L. E. Orgel, Discuss. Faraday Soc. 29, 32 (1960);
6. (b) H. M. McConnell, J. Chem. Phys. 35, 508 (1961).
7. S. Larsson, J. Am. Chem. Soc. 103, 4034 (1981).
8. (a) B. Källebring and S. Larsson, Chem. Phys. Lett. 138, 76 (1987);
9. (b) S. F. Fischer and P. O. Scherer, Chem. Phys. 115, 151 (1987);
10. (c) P. O. Scherer and S. F. Fischer, 131, 115 (1989);
11. (d) M. Plato, K. Möbius, M. E. Michel-Beyerle, M. E. Bixon, and J. Jortner, J. Am. Chem. Soc. 110, 7279 (1988);
12. (e) A. Warshel, S. Creighton, and W. W. Parson, J. Phys. Chem. 92, 2696 (1988).
13. R. J. Cave, D. V. Baxter, W. A. Goddard III, and J. D. Baldeschwieler, J. Chem. Phys. 87, 926 (1987).
14. (a) M. D. Newton, J. Phys. Chem. 95, 30 (1991), and references cited therein;
15. (b) 92, 3049 (1988);
16. (c) M. D. Newton in Cluster Models for Surface and Bulk Phenomena, edited by G. Pacchioni et al. (Plenum, New York, 1992);
17. (d) M. D. Newton, J. Phys. Chem. 95, 30 (1991);
18. (e) A. Broo and S. Larsson, Chem. Phys. 161, 363 (1992).
19. M. D. Newton, K. Ohta, and E. Zhong, J. Phys. Chem. 95, 2317 (1991).
20. (a) K. D. Jordan and M. N. Paddon-Row, J. Phys. Chem. 96, 1188 (1992);
21. (b) K. Kim, K. D. Jordan, and M. N. Paddon-Row, ibid. 98, 11 053 (1994).
22. (a) L. A. Curtiss, C. A. Naleway, and J. R. Miller, J. Phys. Chem. 97, 4050 (1993);
23. (b) 99, 1182 (1995).
24. (a) C. Liang and M. D. Newton, J. Phys. Chem. 96, 2855 (1992);
25. (b) 97, 3199 (1993).
26. (a) M. Braga and S. Larsson, Chem. Phys. Lett. 213, 217 (1993);
27. (b) M. Braga, A. Broo, and S. Larsson, Chem. Phys. 156, 1 (1991).
28. P. Siddarth and R. A. Marcus, J. Phys. Chem. 96, 3213 (1992).
29. (a) M. A. Ratner, J. Phys. Chem. 94, 4877 (1990);
30. (b) see, for example, J. J. Regan, S. M. Risser, D. N. Beratan, and J. N. Onuchic, ibid. 97, 13 083 (1993), and references cited therein.
31. R. J. Cave, M. D. Newton, K. Kumar, and M. B. Zimmt, J. Phys. Chem. 99, 17501 (1995).
32. (a) I. Shavitt, in Modern Theoretical Chemistry, Vol. 3 Methods of Electronic Structure Theory, edited by H. F. Schaefer, III (Plenum, New York, 1977), p. 189;
33. (b) D. K. W. Mok, R. Neumann, and N. C. Handy, J. Phys. Chem. 100, 6225 (1996).
34. (a) R. J. Cave and M. D. Newton, Chem. Phys. Lett. 249, 15 (1996); (b) in footnote 5 of Ref. 18(a), "2√β" should be replaced with the correct expression, "√2β"; the GMH model is an extension of the earlier two-state Mulliken-Hush model (Refs. 18(c)-18(f));
35. (c) R. S. Mulliken, J. Am. Chem. Soc. 64, 811 (1952);
36. (d) R. S. Mulliken and W. B. Person, Molecular Complexes (Wiley, New York, 1969);
37. (e) N. S. Hush, Prog. Inorg. Chem. 8, 391 (1967);
38. (f) N. S. Hush, Electrochim. Acta 13, 1005 (1968);
39. (g) J. R. Reimers and N. S. Hush, J. Phys. Chem. 95, 9773 (1991);
40. (h) C. Creutz, M. D. Newton, and N. Sutin, J. Photochem. Photobiol. A 82, 47 (1994).
41. (a) T. Pacher, L. S. Cederbaum, and H. Köppel, J. Chem. Phys. 89, 7367 (1988);
42. (b) L. S. Cederbaum, J. Schirmer, and H.-D. Meyer, J. Phys. A 22, 2427 (1989);
43. (c) T. Pacher, H. Köppel, and L. S. Cederbaum, J. Chem. Phys. 95, 6668 (1991);
44. (d) T. Pacher, L. S. Cederbaum, and H. Köppel, Adv. Chem. Phys. 84, 293 (1993); see also V. Sidis, ibid. 82, 73 (1992).
45. (d) T. Pacher, L. S. Cederbaum, and H. Köppel, Adv. Chem. Phys. 84, 293 (1993); see also V. Sidis, ibid. 82, 73 (1992).
46. (a) W. Domcke and C. Woywood, Chem. Phys. Lett. 216, 362 (1993);
47. (b) W. Domcke, C. Woywood, and M. Stengle, ibid. 226, 257 (1994).
48. W. Jarzeba, K. Thakur, A. Hermann, and P. F. Barbara, J. Phys. Chem. 99, 2016 (1995).
49. Additional comparative studies are currently underway for systems which include through-bond (TB) coupling of D-A sites via superexchange mechanisms (Ref. 1).
50. (a) H.-J. Werner and W. Meyer, J. Chem. Phys. 74, 5802 (1981);
51. (b) A. Macias and A. Riera, J. Phys. B 11, L489 (1978);
52. (c) S. Kato and Y. Amatatsu, J. Chem. Phys. 92, 7241 (1990);
53. (d) H. J. Kim, R. Bianco, B. J. Gertner, and J. T. Hynes, J. Phys. Chem. 97, 1723 (1993);
54. (e) F. A. Gadea and M. Pelissier, J. Chem. Phys. 93, 545 (1990).
55. M. D. Newton and N. Sutin, Annu. Rev. Phys. Chem. 34, 437 (1984).
56. B. O. Roos, Adv. Chem. Phys. 69, 399 (1987).
57. An alternative method for the formulation of diabatic orbitals entails use of the GMH method. Given a set of adiabatic states spanning the localized (diabatic) states of importance for the processes in question, one first performs a GMH analysis to obtain the GMH diabatic states. Formation of the natural orbitals (based on the one-electron density) for each of the GMH diabatic states yields, in the CASSCF active space, localized natural orbitals. The localized natural orbitals from one of the states may be used directly as reference diabatic orbitals in the BD procedure, or it may be that the active spaces from two or more GMH diabatic states need to be combined in order to span the set of localized orbitals required to describe the system.
58. ZnZn. Projection of these multiconfigurational zeroth-order states onto the adiabatic CI vectors at a given geometry is performed analogously to the single-CSF case; (b) For doublets involving three half-filled orbitals there are two linearly-independent spin eigenfunctions, and hence two CSFs which are included in the projector.
59. P. O. Löwdin, Adv. Quantum Chem. 5, 185 (1970).
60. Since the adiabatic states are obtained as eigenvectors of a CI matrix, the P and Q spaces are automatically orthogonal, independent of any rotations done within the P space.
61. C. F. Bender and H. F. Schaefer III, J. Chem. Phys. 55, 4798 (1977).
62. In cases where the active space canonical molecular orbitals obtained from the nSA/CASSCF procedure (at points of high symmetry or large D-A separation) are adequately localized on the sites of interest, they could be substituted for the ANO's and serve equally well as reference diabatic orbitals.
63. MOLCAS versions 2 and 3, K. Andersson, M. P. Fülscher, R. Lindh, P.-A. Malmqvist, J. Olsen, B. O. Roos, and A. J. Sadlej, University of Lund, Sweden, R. A. Blomberg and P. E. M. Siegbahn, University of Stockholm, Sweden, V. Kello, J. Noga, and M. Urban, Commenus University, Slovakia and P.-O. Widmark, IBM, Sweden, 1991.
64. A. J. H. Wachters, J. Chem. Phys. 52, 1033 (1970).
65. R. C. Raffenetti, J. Chem. Phys. 58, 4452 (1973).
66. GAMESS is a general purpose electronic structure program. The original version of the program was assembled by M. Dupuis, D. Spangler, and J. J. Wendoloski. The current version is described in M. W. Schmidt, K. K. Baldridge, J. A. Boatz, S. T. Elbert, M. S. Gordon, J. H. Jensen, S. Koseki, N. Matsunaga, K. A. Nguyen, S. J. Shu, T. L. Windus, M. Dupuis, J. A. Montgomery, J. Comput. Chem. 14, 1347 (1993).
67. W. S. Benedict, N. Gailer, and E. K. Plyler, J. Chem. Phys. 24, 1139 (1956).
68. (a) P.-O. Widmark, P.-Å. Malmqvist, and B. O. Roos, Theor. Chim. Acta. 77, 291 (1990);
69. (b) P.-O. Widmark, B. J. Persson, and B. O. Roos, ibid. 79, 419 (1991).
70. 2.
71. ab) and other energetics, will be reported elsewhere, R. J. Cave, Y.-P. Liu, and M. D. Newton (to be published).
72. C. E. Moore, Atomic Energy Levels, VoI II, Natl. Bur. Stand. (U.S.) Circ 467 (1949).
73. p′ (in each case the members of a given pair differed by <10%).
74. J. R. Reimers and N. S. Hush, Inorg. Chim. Acta 226, 33 (1994).
75. M. E. Sigman and G. L. Closs, J. Phys. Chem. 95, 5012 (1991).
76. (a) Other formulations similar in spirit to the BD method employed here (Refs. 19 and 20) have been reported by Ruedenberg and co-workers (Refs. 44(b)-44(d)). We thank Professor Ruedenberg for supplying us with a copy of Ref. 43(d) prior to publication, and also for a number of helpful discussions;
77. (b) K. Ruedenberg and G. Atchity, J. Chem. Phys. 99, 3799 (1993);
78. (c) G. Atchity and K. Ruedenberg, ibid. 99, 3790 (1993);
79. (d) K. Ruedenberg and G. Atchity, preprint.
80. L. E. Nitzsche and E. R. Davidson, J. Chem. Phys. 68, 3103 (1978).
81. P. M. Kozlowski and E. R. Davidson, J. Chem. Phys. 100, 3672 (1994).
82. M. R. Hoffmann, Chem. Phys. Lett. 210, 193 (1993).
83. (a) D. Mukherjee and S. Pal, Adv. Quantum Chem. 20, 292 (1989);
84. (b) K. Jankowski, J. Paldus, I. Grabowski, and K. Kowalski, J. Chem. Phys. 97, 7600 (1992);
85. (c) L. Meissner and R. J. Bartlett, ibid. 91, 4800 (1989);
86. L. Meissner, S. Kucharski, and R. J. Bartlett, ibid. 93, 1847 (1990).
87. Use of the BD transformation of the adiabatic nSA/CASSCF Hamiltonian matrix yields diabatic states as in the present application. Once one has obtained the CASSCF diabatic Hamiltonian matrix, the diagonal Hamiltonian matrix elements can be augmented with correlation energies obtained at a geometry where the diabatic states are noninteracting (e.g., a point of high symmetry or large distance). Rediagonalization of the augmented Hamiltonian matrix yields new adiabatic states based on the above diagonal dressing of the diabatic states. This procedure has been used to study the interaction of the low-lying triplet states of pyrazine as a function of bending, R. J. Cave (to be published).