Quantum Dynamics Simulations of Excited State ... - ACS Publications

Sep 23, 2016 - Quantum Dynamics Simulations of Excited State Energy Transfer in a. Zinc−Free-Base Porphyrin Dyad. Judah S. High,. †. Luis G. C. Re...
0 downloads 0 Views 2MB Size
Subscriber access provided by Heriot-Watt | University Library

Article

Quantum Dynamics Simulations of Excited State Energy Transfer in a Zinc - Free-Base Porphyrin Dyad Judah Sterling High, Luis G.C. Rego, and Elena Jakubikova J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.6b05739 • Publication Date (Web): 23 Sep 2016 Downloaded from http://pubs.acs.org on September 24, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Quantum Dynamics Simulations of Excited State Energy Transfer in a Zinc - Free-Base Porphyrin Dyad Judah S. High,a Luis G.C. Rego,b Elena Jakubikovaa,* a

Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695,

United States b

Department of Physics, Federal University of Santa Catarina, Florianopolis, Brazil

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 26

Abstract

Rational design of artificial light-harvesting molecular architectures entails building systems that absorb strongly in the visible and near-IR region of the electromagnetic spectrum and also funnel excited state energy to a single site. Ability to model nonadiabatic processes, such as excited state energy transfer (EET), that occur in such systems on picosecond time scale can aid in the development of novel artificial light harvesting arrays. A combination of density functional theory (DFT), time-dependent DFT, tight-binding molecular dynamics, and quantum dynamics is employed here to simulate EET in the ZnFbΦ dyad, a model artificial light-harvesting array that undergoes EET with an experimentally measured rate constant of (3.5 ps)-1 upon excitation at 550 nm in toluene [Yang et al., J. Phys. Chem. B 1998, 102, 9426-9436]. We find that in order to successfully simulate the EET process, it is important to (1) include coupling between nuclear and electronic degrees of freedom in the QD simulation, (2) account for Coulomb coupling between the electron and hole wavepackets, and (3) parameterize the extended Hückel model Hamiltonian employed in the QD simulations with respect to the DFT.

ACS Paragon Plus Environment

2

Page 3 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1. Introduction Photosynthetic organisms employ light harvesting systems1-9 to transform light energy into chemical energy. Natural photosynthesis relies heavily on excitation energy transfer10-16 (EET) such that, upon photo-absorption, excited state energy is funneled toward a reaction center to ultimately be used in chemical reactions.17-21 In nature, bacteriochlorophylls (BChl) and chlorophylls (Chl) are arranged into large macrostructures, which undergo EET with near perfect quantum efficiency.22 Substantial effort has been dedicated to designing artificial systems that mimic the lightabsorption, energy transfer and charge separation capabilities of pigment assemblies in plants and bacteria.23-25 In particular, studies of porphyrin dyads connected through various linkers have shown that EET rates in such artificial systems depend on donor-acceptor distance as well as orientation,26-29 steric hindrance of the linker,30 linker-attachment site,31 array topology29,30,32 and the presence of electron donating or withdrawing groups.31,33 Yang and Kim have summarized important factors that affect EET rates in multiporphyrin arrays that include the extent of excitonic dipole-dipole interactions and extensive π-conjugation.34 Porphyrin dyads, such as the meso-meso phenylene-linked zinc-free-base porphyrin (ZnFbΦ, see Figure 1), represent simple model systems for studies of EET in multiporphyrin arrays. The ZnFbΦ dyad undergoes rapid EET from the zinc porphyrin (Zn) to the freebase porphyrin (Fb), (kEET = (3.5 ps)-1) upon photoexcitation at 550 nm.33 Experimental studies comparing absorption and redox properties of the ZnFbΦ dyad with those of the Zn and Fb monomers confirmed relatively weak electronic coupling between the monomers in the dyad.33 However, it was noted that the electronic coupling in ZnFbΦ is stronger than in porphyrin dyads with longer linkers and slower EET rates.33 Due to the abundance of experimental data, ZnFbΦ represents an ideal model system for the testing and development of new computational approaches to simulations of EET processes. Computational studies play an important role in elucidating the underlaying mechanism of charge and excited state energy transfer in molecular arrays as well as dye-nanoparticle assemblies. Approaches based on the reduced density matrix, such as the Redfield theory, are often used to model quantum dynamics in light-harvesting systems.35 For example, Redfield theory has been successfully employed in modeling energy transfer within the reaction center of photosystem II,36 and a modified Redfield theory was used to model energy transfer and trapping

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 26

kinetics in Photosystem I.36 A number of other quantum dynamics methodologies, both fully quantum mechanical and mixed quantum-classical (MQC),37,38 have been developed and applied to simulate photoinduced dynamics in molecular and condensed matter systems. Fully quantum mechanical methodologies include the time dependent hartree (TDH) methods.37-44 Some of the predominant MQC methods inculde the Ehrenfest,45 fewest switches surface hopping,38,46-51 density matrix propagation,52,53 quantum Liouville46 and Bohmian mechanics methodologies.54-56

Figure 1. Molecular structure of meso-meso p-phenylene linked zinc-freebase porphyrin dyad (ZnFbΦ).

In this work, we investigate the EET in a model ZnFbΦ dyad (see Figure 1), utilizing a mixed time-dependent quantum-classical dynamics method that propagates wavepackets by solving the time-dependent Schrödinger equation using an extended Hückel (EH) Hamiltonian.57 This approach has been previously employed in the simulations of excited state charge separation in a carotenoid-porphyrin-fullerene system,58 dynamics of metal-to-ligand charge transfer excited states in solvated tris-bipyridine ruthenium(II) complexes,59 and interfacial electron transfer in sensitized TiO2 anatase.60-65 Compared to conventional (ground-state) molecular dynamics, nonadiabatic molecular dynamics (NA-MD) is a much more demanding computational method, owing to the fact that it is carried out in the excited state and also because the time steps of simulation must be considerably smaller than in common ground-state MD. The most common NA-MD methods separate the electron and nuclear dynamics propagation, by treating the nuclei by classical mechanics in mixed quantum–classical approaches. Some high level of ab initio NA-MD methods are Newton-X,66 Molpro,67 CPMD,68 Octopus69 and JADE.70 Essentially, the theory used to calculate the potential energy surfaces determines the computational cost of NA-MD methods. Therefore, the large number of degrees of freedom hinders the use of first principle methods, especially if there is the necessity to describe excited state dynamics. The demand for

ACS Paragon Plus Environment

4

Page 5 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

computationally feasible NA-MD methods promoted the use of semiempirical methods for mixed quantum–classical dynamics based on the extended Hückel method, such as the one used in this work. Another approach to accelerate the simulation consists of calculating the electron dynamics upon a previously calculated nuclear trajectory. This is known as sequential MM-QM approach58,65 or, equivalently, the classical path approximation (CPA),71 as implemented in the PYXAID package.72 It is based on the assumption that the excited potential energy surface does not change remarkably from ground to excited state, or when it changes more slowly than the process studied. Since the ZnFbΦ dyad studied in this work does not undergo large structural changes upon excitation (i.e., isomerization, dissociation, etc., see ref. 33), this first approximation to nuclear-electron coupling should provide a good qualitative description of the process investigated. The sequential MM-QD approach is faster from the computational viewpoint than the direct coupling between the electronic and nuclear degrees of freedom, and provides a semi-quantitative description of dynamic phenomena in the excited state. Thus it allows us to model the quantum dynamics of energy transfer in a molecular system comprised of almost 200 atoms within the atomistic framework on a picosecond time scale, using only moderate computing resources. In addition to the simulations of EET in ZnFbΦ dyad, a new approach for construction of initial wavepacket based on the results of time-dependent density functional theory calculations is also presented. This work represents the first time this particular QD approach has been used to simulate EET in a molecular system with a semi-quantitative accuracy and thus opens up its application to a wider class chemical systems. 2. Computational Methodology Electronic Structure Calculations Structure optimization of ZnFbΦ was performed with the B3LYP exchange-correlation functional.73-75 The 6-31G* basis set76,77 was employed for C, H, N and O atoms, while the LANL0878,79 effective core potential and associated basis set were used for zinc. A frequency calculation verified that the optimized structure lies at a minimum of its potential energy surface. The Gaussian 09 software package80 was employed for all DFT calculations.

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 26

All UV-vis spectra were calculated using time-dependent DFT (TD-DFT),81-83 as implemented in Gaussian 09,84 at the B3LYP/LANL08(Zn),6-31G*(C,H,N) level of theory. Natural transition orbitals (NTOs) were employed to characterize all excitations (see Figures S2-S11 in the SI).85 Molecular Dynamics Simulations Using the optimized ZnFbΦ as the initial structure, ten independent MD simulations of ZnFbΦ were performed using the tight-binding density functional theory (DFTB) methodology employing the DFTB+ software package.50,86-92 The MD simulations consisted of a 20 ps equilibration period (NVT ensemble) at 298 K, followed by a production simulation of 10 ps in the NVE ensemble. Nuclear configurations used in quantum dynamics were those obtained in the NVE ensemble. A 5-chain Nosé-Hoover93-97 thermostat with a coupling frequency of 3600 cm-1 was employed during the 20 ps equilibration period. All simulations were run with a 0.1 fs time step. Due to computational limitations, no solvent model was used during the MD simulations. Extended Hückel Parametrization The EH Hamiltonian was reparametrized with respect to the optimized ZnFbΦ Kohn-Sham Molecular Orbital (KS MO) shapes, energetic ordering, and energy gaps. Reparametrization involved first tuning atomic orbital (AO) ionization potentials (IP). The Wolfsberg-Helmholtz98101

parameter, κ, was set to 2.05 (κdefault = 1.75) (see Table S1 in the SI) to assist in reproducing

the KS HOMO-LUMO gap as well as various MO energy gaps from the HOMO-3 to LUMO+3. All other IPs were left default as prescribed by Alvarez.102 The manual EH parameters were than optimized further using the genetic algorithm employed in DynEMol.57 The optimized EH parameters were utilized during all QD simulations. Excited State Energy Transfer Simulations The quantum dynamics methodology developed by Rego and Batista57-59,65,103,104 was utilized to model EET in ZnFbΦ. In this model, initial electron and hole wavepackets, Φ El (0) and

Φ Hl (0) , are constructed as linear combinations of the EH MOs of the donor zinc trimesitylporphyrin subunit. The EH molecular orbitals (MOs) of Zn selected to construct the wavepackets correspond to Kohn-Sham (KS) MOs describing Q-band excitations as calculated by the TD-DFT.82-84,105-108 The selected excitations correspond with the excitation (550 nm) employed in experimental investigation of the EET rate in ZnFbΦ.33,109,110

ACS Paragon Plus Environment

6

Page 7 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The donor KS MOs and weights used in the construction of Φ El (0) and Φ Hl (0) were calculated from the transition density eigenvectors of the Casida equations81,111 for selected excitations of the Zn subunit. The EH MOs that best represent the KS MOs were determined before construction of the initial linear combinations of EH MOs. The mathematical details for wavepacket construction are given in the Supporting Information. Ten quantum dynamics simulations of EET were done through time propagation of electron and hole wavepackets, Φ El (t) and Φ Hl (t) , by solving the time-dependent Schrödinger equation with the EH98,112-122 Hamiltonian.59 The wavepackets were propagated for 10 ps, with a 0.1 fs time step, on pre-computed trajectories obtained from molecular dynamics (MD) simulations. Electrostatic Coulomb coupling between Φ El (t) and Φ Hl (t) was calculated at each step of the QD simulation. At each step in the time-slice propagation procedure, Φ El (t) and Φ Hl (t) were projected onto the atomic basis of Zn, p-phenylene, and Fb to obtain survival probabilities P(t), i.e., partial occupations of Φ El (t) and Φ Hl (t) on the Zn, p-phenylene, or Fb subunits.59 Since the wavepacket at t = 0 is constructed to be fully localized on the Zn-porphyrin donor, P(t=0) is always one for the Zn-subunit and zero for the acceptor (Fb) or bridge (pphenylene). To describe EET in ZnFbΦ, P(t) for the electron and hole wavepackets on the Zn subunit must decrease to zero, while simultaneously increasing to unity for the Fb subunit. 3. Results and Discussion 3.1. Ground state structure The optimized structure of ZnFbΦ is such that the Zn and Fb porphyrin subunits are coplanar while the p-phenylene bridge plane lies at an angle of about 60° from the planes of each tetrapyrrole. The incorporation of phenylene linkers in multiporphyrin arrays is thought to preserve the electronic characteristics of the porphyrin subunits.33 This is supported by the shapes and energies of the frontier molecular orbitals of the optimized ZnFbΦ system (see Figure 2). The eight frontier Kohn-Sham (KS) MOs of ZnFbΦ spanning from the HOMO-3 to LUMO+3 are localized on either the Zn or Fb subunits. These localized MOs resemble the Zn and Fb subunit frontier MOs as shown in Figure 2. Furthermore, there is little splitting of the frontier KS

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 26

MO energy levels in going from a situation in which the chromophores are isolated to the situation in which they are covalently linked via the p-phenylene bridge (see Figure 2).

Figure 2. Frontier molecular orbitals of the Zn and Fb porphyrin subunits as well as the ZnFbΦ dyad calculated at the B3LYP level of theory.

Previous experimental work showed that the one-electron redox potentials of the ZnFbΦ system remain unchanged in comparison to the one-electron redox potentials of isolated Zn and Fb chromophores.33 Taking Koopman’s theorem as a first approximation to ground state redox potentials,123,124 the ground state electronic structure of ZnFbΦ obtained here is consistent with the experimental observations. This is reflected in the small changes to the HOMO and LUMO energies going from the isolated Zn and Fb subunits to ZnFbΦ. Since EET simulations are done at the extended Hückel level of theory, the EH Hamiltonian was reparametrized to reproduce the ZnFbΦ frontier KS MO shapes and energies. A comparison of the DFT and reparametrized EH MOs is shown in Figure 3. The optimized EH parameters used in all EET simulations are shown in the Table S1 of the SI. Note that the EH MOs are in very close agreement with the DFT KS MOs. As both the EH MO energies and shapes are used during the wavepacket propagation procedure in DynEMol,57-59,65,103,104 it is critical that the EH Hamiltonian reproduces these properties with respect to results obtained at higher level of theory, DFT in this case.

ACS Paragon Plus Environment

8

Page 9 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. A comparison of the MOs of ZnFbΦ obtained at the B3LYP level (left) and those calculated using optimal parameters at the EH level of theory (right).

3.2. Absorption properties and initial state construction Previous experimental investigation of ZnFbΦ revealed that its UV-vis absorption spectrum can be approximately described as the superposition of the Zn and Fb subunit’s absorption spectra, suggesting that the excitation at 550 nm is a Q-band excitation localized on Zn subunit.33 The excitation at 550 nm was also utilized as the initial excitation for experimental EET studies.33 The calculated electronic absorption spectra of ZnFbΦ and the constituent monomers shown in Figure 4 reproduce these experimental findings. The calculated absorption spectrum of ZnFbΦ displays two near-degenerate Zn-localized excitations at 542.25 nm and 540.35 nm, which were confirmed to be nearly indistinguishable from the Q-band excitations of the isolated Zn subunit at 531.60 nm and 531.49 nm by the NTO analysis (see Figure 5). The initial electron and hole wavepackets for the EET simulations were therefore constructed based on the results of TD-DFT calculations on the Zn subunit of the ZnFbΦ dyad (see Appendix).

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

Figure 4. B3LYP calculated electronic absorption spectrum of the Zn and Fb subunits as well as the ZnFbΦ dyad. The subunits exhibit the Q and Soret band and the ZnFbΦ spectrum resembles the superposition of the subunit spectra.

Figure 5. Comparison of NTOs for the second and third lowest energy calculated excitations of ZnFbΦ and the near-degenerate Q-band excitations of the isolated Zn subunit.

Since EET in ZnFbΦ occurs in the pico-second regime,33 electronic coupling to nuclear degrees of freedom is likely to play an important role in the EET process. Ten MD simulations of ZnFbΦ equilibrated at 298 K were carried out to obtain nuclear trajectories over which the

ACS Paragon Plus Environment

10

Page 11 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

electron and hole wavepackets propagate in time. The electronic absorption spectra for each of the initial structures of Zn subunit, taken from the 10 MD simulations (MDSIM1-10), are shown in Figure 6.

Figure 6. B3LYP calculated UV-vis spectra of Zn subunit structure after 20 ps of equilibration at 298 K of ZnFbΦ via MD. All spectra exhibit Soret- and Q-band excitations.

All absorption spectra of Zn subunits shown in Figure 6 exhibit a Soret- and a Q-band with excited state energies and intensities that are within 0.18 eV of the excitation energies and intensities of the optimized Zn subunit. The Q-band excitations found near 550 nm for each Zn subunit involve two contributions from one-electron transitions between the HOMO, HOMO-1 and the LUMO, LUMO+1 of each Zn structure. The contributions from the one-electron transitions to the Q-band excitations are used directly in the wavepacket construction procedure summarized in the Appendix. Ultimately, weighting coefficients are calculated for EH MOs, which are used in forming the electron and hole wavepackets as a linear combination of EH MOs. Care is taken during the wavepacket construction procedure to identify the EH MOs that match the KS MOs contributing to the excitation, which were calculated for the thermalized Zn subunit structures. Figure 7 illustrates the result of the wavepacket construction procedure by showing the Zn EH MOs contributing to the electron and hole wavepacket along with their weighting coefficients. The electron and hole wavepacket densities resulting from the linear combination of EH MOs are also shown in Figure 7.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 26

Figure 7. Initial EH electron and hole wavepacket densities for a representative QD simulation. Wavepackets were constructed as a linear combination of the EH MOs that best resemble the geometry optimized KS MOs involved with a given excitation as prescribed by TD-DFT. The coefficients are calculated using the procedure outlined in the Appendix.

3.3. Excited state energy transfer simulations Wavepacket dynamics was employed to perform EET simulations in ZnFbΦ dyad, with the initial hole and electron excitations set up based on TD-DFT calculations on the Zn subunit. The survival probabilities, P(t), for the electron and hole wavepackets on the Zn, p-phenylene and Fb subunits of the ZnFbΦ dyad are shown in Figure 8 for a series of QD simulations utilizing computational models of increasing complexity. Case A corresponds to the simplest surmise, in which QD simulations were performed upon the static optimized ZnFbΦ structure, with Coulomb coupling between the electron and hole wavepackets neglected. Both electron and hole wavepackets remain mainly localized on the Zn subunit throughout the simulation. Neither the addition of nuclear motion (model B) nor the addition of Coulomb coupling (model C) into the model alone results in the EET in ZnFbΦ dyad on 1 ps time scale. The complete EET transfer is only observed when both nuclear motion and Coulomb coupling are included in the model (model D). Reparametrization of EH Hamiltonian is equally important for the success of the EET simulations. Utilization of the default Alvarez parameters in the QD simulations leads to incorrect description of the electron and hole wavepacket dynamics set, even if nuclear motion and Coulomb coupling are included in the model (model E). It is clear that in order to describe EET in ZnFbΦ system, one must include the following in the model: (1) nuclear motion, (2) electron-hole Coulomb coupling, and (3) EH Hamiltonian parametrized with respect to the DFT results.

ACS Paragon Plus Environment

12

Page 13 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 8. Survival probabilities for the Zn, p-phenylene and Fb subunits during QDSIM1 in which combinations of having nuclear motion, Coulombic wavepacket coupling and optimal EH parameters were included.

Figure 9 displays survival probabilities of electron and hole wavepackets on Zn subunit from ten independent QD simulations of EET in ZnFbΦ. All simulations utilize parametrized extended Hückel Hamiltonian and include nuclear motion as well as Coulomb coupling between electron and hole wavepackets. The survival probability for both electron and hole wavepackets on Zn subunit becomes negligible within 1.5 ps upon initial excitation of Zn subunit, suggesting a complete excited state energy transfer. Complete EET can be confirmed by inspecting survival probabilities on Fb and p-phenylene subunits, shown in Figures S15 and S16 in the SI. To further evaluate the importance of nuclear motion for simulating EET in molecular systems, ten additional QD simulations were performed on ten different static structures. The static structures were chosen to be the same as the initial structures of the ten MD production runs. The QD simulations utilized optimized EH parameters and included electron-hole Coulomb coupling in the model. Interestingly, two of the ten QD simulations displayed complete EET transfer in the time frame of 1.5 ps, while the wavepacket remained localized on the Zn subunit in the rest of the simulations (see Figure S17 in the SI). This suggests that certain conformations of ZnFbΦ dyad will facilitate the EET even at the frozen geometry. The structures resulting in EET are only accessible through the exploration of phase space and it may prove difficult to consistently

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 26

predict which static structures will result in EET a priori. A reliable simulation of the EET in this system thus requires the inclusion of nuclear motion into the model.

Figure 9. Survival probability for the electron and hole on the Zn subunit of ZnFbΦ for QDSIM1-10. All simulations show EET by means of the survival probabilities decaying to zero.

Examination of the electron and hole wavepacket densities for each EET simulation, as shown in Figure 10 for one of the QD simulations, reveals how electronic excitation energy is transferred from the Zn to the Fb. It is evident from Figure 10 that by 1.0 fs, wavepacket density has begun to delocalize onto the p-phenylene bridge and by 10 fs, wavepacket density has delocalized across the entire dyad. By 1.0 ps, EET has been achieved and the wavepacket density is completely localized on the Fb subunit.

ACS Paragon Plus Environment

14

Page 15 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 10. Shown here are snapshots of the electron and hole wavepacket density of QDSIM1 at various times. Note that by 1.0 ps the wavepacket densities are fully localized on the Fb subunit.

Since the transient absorption experiments used to monitor EET rates in ZnFbΦ measure the average response of a collection of ZnFbΦ molecules to light at 550 nm, the survival probabilities from all 10 QD simulations were averaged (see Figure 11), with the goal of obtaining an average EET rate to be compared to that determined from experiment ((3.5 ps)-1).33 Survival probabilities for the electron and hole wavepackets on the Zn- and Fb- porphyrin were fit to bi-exponential decay functions to obtain rate constants and characteristic times for EET in ZnFbΦ (see Figure 11). No fitting scheme was applied to the survival probability of the pphenylene bridge as it stays close to zero throughout the simulation. The calculated average characteristic time for the EET from the Zn onto the Fb subunit is 0.45 ps. This characteristic time is one order of magnitude faster than the experimentally determined characteristic time of 3.5 ps, suggesting that while our model correctly describes the qualitative behavior of this system, it does not possess quantitative accuracy. We attribute this lack of quantitative accuracy to the neglect of the solvent effects in our model. Solvation dynamics have been previously

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 26

shown to play an important role in obtaining the correct description of excited state dynamics in molecular systems.59

Figure 11. Average survival probabilities of the electron and hole wavepackets averaged over QDSIM110. The average survival probabilities were fitted to bi-exponentials with characteristic times τ. No fit was applied to the p-phenylene electron and hole P(t).

4. Conclusions We have presented a successful simulation of non-adiabatic excited state energy transfer in a prototypical meso-meso phenylene linked zinc-free-base porphyrin dyad, ZnFbΦ, utilizing a mixed time-dependent quantum mechanics/molecular mechanics computational method.57,59,103A procedure was also developed for constructing the initial wavepackets as a linear combination of extended Hückel MOs whose weighting coefficients were determined through excitation analysis of the donor Zn subunit via TD-DFT. This was necessary, since porphyrin excitations cannot be approximated as single one-electron MO-to-MO transitions. The inclusion of nuclear motion, Coulomb coupling between electron and hole wavepackets, and a reparametrized Extended Hückel Hamiltonian was critical for achieving Zn to Fb excitation energy migration in the quantum dynamics simulations. Our simulations yielded EET with a time constant of 0.45 ps, which is about an order of magnitude faster than the EET time

ACS Paragon Plus Environment

16

Page 17 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

constant of 3.5 ps determined from transient absorption measurements.33 We attribute this discrepancy to the neglect of solvation dynamics in our computational model. Overall, we have shown that a mixed time-dependent quantum mechanics/molecular dynamics computational approach can be used to simulate EET in rigid molecular systems with a semi-quantitative accuracy. Associated Content Supporting Information Description of a procedure for construction of initial wavepackets based on TD-DFT calculations. Optimal extended Hückel parameters used in QD simulations; frontier B3LYP MOs of initial ZnFbΦ structures used in QD simulations; NTOs of Q band excitations in ZnFbΦ and Zn donor at all initial QD structures; initial electron and hole wavepacket densities for QDSIM110; QD simulation P(t) plots demonstrating importance of factors resulting in EET; P(t) associated with the p-phenylene bridge and Fb-porphyrin acceptor for QDSIM1-10; and P(t) for ten static QD simulations. Author Information Corresponding Author *E-mail: [email protected], Tel. 919-515-1808 Notes The authors declare no competing financial interest. Acknowledgments This work was supported by the Department of Chemistry at North Carolina State University (NCSU) and the High Performance Computing (HPC) services at NCSU.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 26

References 1. Deisenhofer, J.; Epp, O.; Miki, K.; Huber, R.; Michel, H. X-Ray Structure Analysis of a Membrane Protein Complex. J. Mol. Biol. 1984, 180, 385-398. 2. Oh-oka, H. Type 1 Reaction Center of Photosynthetic Heliobacteria. Photochem. Photobiol. 2007, 83, 177-186. 3. Camara-Artigas, A.; Blankenship, R. E.; Allen, J. P. The Structure of the Fmo Protein from Chlorobium Tepidum at 2.2 Å Resolution. Photosynth. Res. 2003, 75, 49-55. 4. Cogdell, R. J.; Gardiner, A. T.; Roszak, A. W.; Law, C. J.; Southall, J.; Isaacs, N. W. Rings, Ellipses and Horseshoes: How Purple Bacteria Harvest Solar Energy. Photosynth. Res. 2004, 81, 207-214. 5. Olson, J. M. The FMO Protein. In Discoveries in Photosynthesis. Govindjee; Beatty, J. T.; Gest, H.; Allen, J. F.; Eds. Advances in Photosynthesis and Respiration. Springer Netherlands: Dordrecht, Netherlands, 2005, 20, 421-427. 6. Deisenhofer, J.; Epp, O.; Sinning, I.; Michel, H. Crystallographic Refinement at 2.3 a Resolution and Refined Model of the Photosynthetic Reaction Centre from Rhodopseudomonas Viridis. J. Mol. Biol. 1995, 246, 429-457. 7. Koepke, J.; Hu, X.; Muenke, C.; Schulten, K.; Michel, H. The Crystal Structure of the Light-Harvesting Complex II (B800–850) from Rhodospirillum Molischianum. Structure 1996, 4, 581-597. 8. McDermott, G.; Prince, S. M.; Freer, A. A.; Hawthornthwaite-Lawless, A. M.; Papiz, M. Z.; Cogdell, R. J.; Isaacs, N. W. Crystal Structure of an Integral Membrane LightHarvesting Complex from Photosynthetic Bacteria. Nature 1995, 374, 517-521. 9. Sener, M. K.; Olsen, J. D.; Hunter, C. N.; Schulten, K. Atomic-Level Structural and Functional Model of a Bacterial Photosynthetic Membrane Vesicle. Proc. Natl. Acad. Sci. U. S. A. 2007, 104, 15723-15728. 10. Davydov, A. S. The Theory of Molecular Excitons. Sov. Phys. Uspekhi 1964, 7, 145-178. 11. Harcourt, R. D.; Scholes, G. D.; Ghiggino, K. P. Rate Expressions for Excitation Transfer. II. Electronic Considerations of Direct and through–Configuration Exciton Resonance Interactions. J. Chem. Phys. 1994, 101, 10521. 12. Franck, J.; Teller, E. Migration and Photochemical Action of Excitation Energy in Crystals. J. Chem. Phys. 1938, 6, 861. 13. Scholes, G. D. Long-Range Resonance Energy Transfer in Molecular Systems. Annu. Rev. Phys. Chem. 2003, 54, 57-87. 14. Kasha, M.; Rawls, H. R.; Ashraf El-Bayoumi, M. The Exciton Model in Molecular Spectroscopy. Pure. Appl. Chem. 1965, 11. 15. Caprasecca, S.; Mennucci, B. Excitation Energy Transfer in Donor-BridgeAcceptor Systems: A Combined Quantum-Mechanical/Classical Analysis of the Role of the Bridge and the Solvent. J. Phys. Chem. A 2014. 16. Albinsson, B.; Martensson, J. Excitation Energy Transfer in Donor-BridgeAcceptor Systems Phys. Chem. Chem. Phys. 2010, 12, 7338-7351.

ACS Paragon Plus Environment

18

Page 19 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

17. Hu, X.; Ritz, T.; Damjanovic, A.; Autenrieth, F.; Schulten, K. Photosynthetic Apparatus of Purple Bacteria. Q. Rev. Biophys. 2002, 35. 18. Karreman, G.; Steele, R. H. On the Possibility of Long Distance Energy Transfer by Resonance in Biology. Biochim. Biophys. Acta 1957, 25, 280-291. 19. Hsin, J.; Strumpfer, J.; Sener, M.; Qian, P.; Hunter, C. N.; Schulten, K. Energy Transfer Dynamics in an RC-LH1-PufX Tubular Photosynthetic Membrane. New J. Phys. 2010, 12. 20. van Grondelle, R.; Dekker, J. P.; Gillbro, T.; Sundstrom, V. Energy Transfer and Trapping in Photosynthesis. Biochim. Biophys. Acta, Bionenerg. 1994, 1187, 1-65. 21. Muh, F.; Madjet , M.-A.; Adolphs, J.; Abdurahman, A.; Rabenstein, B.; Ishikita, H.; Knapp, E. W.; Renger, T. Alpha-Helices Direct Excitation Energy Flow in the Fenna Matthews Olson Protein Proc. Natl. Acad. Sci. U. S. A. 2007, 104, 16862-16867. 22. Cogdell, R. J.; Gall, A.; Kohler, J. The Architecture and Function of the LightHarvesting Apparatus of Purple Bacteria: From Single Molecules to in Vivo Membranes. Q. Rev. Biophys. 2006, 39, 227-324. 23. Gust, D.; Moore, T. A. Photosynthetic Model Systems. In Photoinduced Electron Transfer III. Mattay, J.; Eds. Topics in Current Chemistry. Springer Berlin Heidelberg: Berlin, Heidelberg, Germany, 1991, 159, 103-151. 24. Konovalova, N. V.; Evstigneeva, R. P.; Luzgina, V. N. Synthetic Molecular Systems Based on Porphyrins as Models for the Study of Energy Transfer in Photosynthesis. Russ. Chem. Rev. 2001, 70, 939-969. 25. Yang, J.; Yoon, M. C.; Yoo, H.; Kim, P.; Kim, D. Excitation Energy Transfer in Multiporphyrin Arrays with Cyclic Architectures: Towards Artificial Light-Harvesting Antenna Complexes. Chem. Soc. Rev. 2012, 41, 4808-4826. 26. Strachan, J.-P.; Gentemann, S.; Seth, J.; Kalsbeck, W. A.; Lindsey, J. S.; Holten, D.; Bocian, D. F. Effects of Orbital Ordering on Electronic Communication in Multiporphyrin Arrays. J. Am. Chem. Soc. 1997, 119, 11191-11201. 27. Cho, H. S.; Jeong, D. H.; Yoon, M.-C.; Kim, Y. H.; Kim, Y.-R.; Kim, D.; Jeoung, S. C.; Kim, S. K.; Aratani, N.; Shinmori, H.; et al. Excited-State Energy Transfer Processes in Phenylene- and Biphenylene-Linked and Directly-Linked Zinc(II) and Free-Base Hybrid Diporphyrins. J. Phys. Chem. A 2001, 105, 4200-4210. 28. Maruyama, K.; Osuka, A.; Mataga, N. A Chemical Approach Towards Photosynthetic Reaction Center (2). Pure. Appl. Chem. 1994, 66. 29. Kim, D.; Osuka, A. Directly Linked Porphyrin Arrays with Tunable Excitonic Interactions. Acc. Chem. Res. 2004, 37, 735-745. 30. Song, H. E.; Kirmaier, C.; Schwartz, J. K.; Hindin, E.; Yu, L.; Bocian, D. F.; Lindsey, J. S.; Holten, D. Effects of Multiple Pathways on Excited-State Energy Flow in SelfAssembled Wheel-and-Spoke Light-Harvesting Architectures. J. Phys. Chem. B 2006, 110, 19131-19139. 31. Yang, S. I.; Seth, J.; Balasubramanian, T.; Kim, D.; Lindsey, J. S.; Holten, D.; Bocian, D. F. Interplay of Orbital Tuning and Linker Location in Controlling Electronic Communication in Porphyrin Arrays. J. Am. Chem. Soc. 1999, 121, 4008-4018.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 26

32. Yoon, M. C.; Cho, S.; Kim, P.; Hori, T.; Aratani, N.; Osuka, A.; Kim, D. Structural Dependence on Excitation Energy Migration Processes in Artificial Light Harvesting Cyclic Zinc(II) Porphyrin Arrays. J. Phys. Chem. B 2009, 113, 15074-15082. 33. Yang, S. I.; Lammi, R. K.; Seth, J.; Riggs, J. A.; Arai, T.; Kim, D.; Bocian, D. F.; Holten, D.; Lindsey, J. S. Excited-State Energy Transfer and Ground-State Hole/Electron Hopping Inp-Phenylene-Linked Porphyrin Dimers. J. Phys. Chem. B 1998, 102, 9426-9436. 34. Yang, J.; Kim, D. Excitation Energy Migration Processes in Various MultiPorphyrin Assemblies. Phil. Trans. R. Soc. A 2012, 370, 3802-3818. 35. Yang, M.; Fleming, G. R. Influence of Phonons on Exciton Transfer Dynamics: Comparison of the Redfield, Förster, and Modified Redfield Equations. Chem. Phys. 2002, 275, 355-372. 36. Leegwater, J. A.; Durrant, J. R.; Klug, D. R. Exciton Equilibration Induced by Phonons: Theory and Application to PS II Reaction Centers. J. Phys. Chem. B 1997, 101, 72057210. 37. Herman, M. F. Dynamics by Semiclassical Methods. Annu. Rev. Phys. Chem. 1994, 45, 83-111. 38. Tully, J. C. Perspective: Nonadiabatic Dynamics Theory. J. Chem. Phys. 2012, 137, 22A301. 39. Gerber, R. B. Time-Dependent Self-Consistent Field Approximation for Intramolecular Energy Transfer. I. Formulation and Application to Dissociation of Van Der Waals Molecules. J. Chem. Phys. 1982, 77, 3022. 40. Beck, M. The Multiconfiguration Time-Dependent Hartree (MCTDH) Method: A Highly Efficient Algorithm for Propagating Wavepackets. Phys. Rep. 2000, 324, 1-105. 41. Wang, H.; Thoss, M. Multilayer Formulation of the Multiconfiguration TimeDependent Hartree Theory. J. Chem. Phys. 2003, 119, 1289. 42. Li, J.; Kondov, I.; Wang, H.; Thoss, M. Quantum Dynamical Simulation of Photoinduced Electron Transfer Processes in Dye-Semiconductor Systems: Theory and Application to Coumarin 343 at TiO2. J. Phys. Condens. Matter 2015, 27, 134202. 43. Borrelli, R.; Thoss, M.; Wang, H.; Domcke, W. Quantum Dynamics of ElectronTransfer Reactions: Photoinduced Intermolecular Electron Transfer in a Porphyrin–Quinone Complex. Mol. Phys. 2012, 110, 751-763. 44. Wang, L.; Long, R.; Prezhdo, O. V. Time-Domain Ab Initio Modeling of Photoinduced Dynamics at Nanoscale Interfaces. Annu. Rev. Phys. Chem. 2015, 66, 549-579. 45. Sawada, S.-I.; Nitzan, A.; Metiu, H. Mean-Trajectory Approximation for Chargeand Energy-Transfer Processes at Surfaces. Phys. Rev. B 1985, 32, 851-867. 46. Wang, L.; Sifain, A. E.; Prezhdo, O. V. Fewest Switches Surface Hopping in Liouville Space. J. Phys. Chem. Lett. 2015, 6, 3827-3833. 47. Drukker, K. Basics of Surface Hopping in Mixed Quantum/Classical Simulations. J. Comput. Phys. 1999, 153, 225-272. 48. Saalmann, U.; Schmidt, R. Non-Adiabatic Quantum Molecular Dynamics: Basic Formalism and Case Study. Z. Phys. D 1996, 38, 153-163. 49. Neukirch, A. J.; Park, J.; Zobac, V.; Wang, H.; Jelinek, P.; Prezhdo, O. V.; Zhou, H. C.; Lewis, J. P. Calculated Photo-Isomerization Efficiencies of Functionalized Azobenzene

ACS Paragon Plus Environment

20

Page 21 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Derivatives in Solar Energy Materials: Azo-Functional Organic Linkers for Porous Coordinated Polymers. J. Phys. Condens. Matter 2015, 27, 134208. 50. Seifert, G. Tight-Binding Density Functional Theory: An Approximate KohnSham DFT Scheme. J. Phys. Chem. A 2007, 111, 5609-5613. 51. Gao, X.; Peng, Q.; Niu, Y.; Wang, D.; Shuai, Z. Theoretical Insight into the Aggregation Induced Emission Phenomena of Diphenyldibenzofulvene: A Nonadiabatic Molecular Dynamics Study. Phys. Chem. Chem. Phys. 2012, 14, 14207-14216. 52. Huo, P.; Coker, D. F. Consistent Schemes for Non-Adiabatic Dynamics Derived from Partial Linearized Density Matrix Propagation. J. Chem. Phys. 2012, 137, 22A535. 53. Micha, D. A. Density Matrix Treatment of Non-Adiabatic Photoinduced Electron Transfer at a Semiconductor Surface. J. Chem. Phys. 2012, 137, 22A521. 54. Burant, J. C.; Tully, J. C. Nonadiabatic Dynamics via the Classical Limit Schrödinger Equation. J. Chem. Phys. 2000, 112, 6097. 55. Zamstein, N.; Tannor, D. J. Non-Adiabatic Molecular Dynamics with Complex Quantum Trajectories. I. The Diabatic Representation. J. Chem. Phys. 2012, 137, 22A517. 56. Zamstein, N.; Tannor, D. J. Non-Adiabatic Molecular Dynamics with Complex Quantum Trajectories. II. The Adiabatic Representation. J. Chem. Phys. 2012, 137, 22A518. 57. Rego, L. G.; Batista, V. S. Quantum Dynamics Simulations of Interfacial Electron Transfer in Sensitized TiO2 Semiconductors. J. Am. Chem. Soc. 2003, 125, 7989-7997. 58. Rego, L. G. C.; Hames, B. C.; Mazon, K. T.; Joswig, J.-O. Intramolecular Polarization Induces Electron–Hole Charge Separation in Light-Harvesting Molecular Triads. J. Phys. Chem. C 2014, 118, 126-134. 59. Hoff, D. A.; Silva, R.; Rego, L. G. C. Subpicosecond Dynamics of Metal-toLigand Charge-Transfer Excited States in Solvated [Ru(bpy)3]2+ Complexes. J. Phys. Chem. C 2011, 115, 15617-15626. 60. Bowman, D. N.; Mukherjee, S.; Barnes, L. J.; Jakubikova, E. Linker Dependence of Interfacial Electron Transfer Rates in Fe(II)-Polypyridine Sensitized Solar Cells. J. Phys. Condens. Matter 2015, 27, 134205. 61. Jakubikova, E.; Bowman, D. N. Fe(II)-Polypyridines as Chromophores in DyeSensitized Solar Cells: A Computational Perspective. Acc. Chem. Res. 2015, 48, 1441-1449. 62. Bowman, D. N.; Blew, J. H.; Tsuchiya, T.; Jakubikova, E. Elucidating BandSelective Sensitization in Iron(II) Polypyridine-TiO2 Assemblies. Inorg. Chem. 2013, 52, 86218628. 63. Mukherjee, S.; Bowman, D. N.; Jakubikova, E. Cyclometalated Fe(II) Complexes as Sensitizers in Dye-Sensitized Solar Cells. Inorg. Chem. 2015, 54, 560-569. 64. Mara, M. W.; Bowman, D. N.; Buyukcakir, O.; Shelby, M. L.; Haldrup, K.; Huang, J.; Harpham, M. R.; Stickrath, A. B.; Zhang, X.; Stoddart, J. F.; et al. Electron Injection from Copper Diimine Sensitizers into TiO2: Structural Effects and Their Implications for Solar Energy Conversion Devices. J. Am. Chem. Soc. 2015, 137, 9670-9684. 65. Monti, A.; Negre, C. F.; Batista, V. S.; Rego, L. G.; de Groot, H. J.; Buda, F. Crucial Role of Nuclear Dynamics for Electron Injection in a Dye-Semiconductor Complex. J Phys. Chem. Lett. 2015, 6, 2393-2398.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 26

(66) Barbatti, M.; Granucci, G.; Persico, M.; Ruckenbauer, M.; Vazdar, M.; EckertMaksić, M.; Lischka, H. The On-the-Fly Surface-Hopping Program System Newton-X: Application to Ab Initio Simulation of the Nonadiabatic Photodynamics of Benchmark Systems. J. Phot. Phot. A: Chem. 2007, 190, 228-240. 67. Levine, B. G.; Coe, J. D.; Virshup, A. M.; Martínez, T. J. Implementation of Ab Initio Multiple Spawning in the Molpro Quantum Chemistry Package. Chem. Phys. 2008, 347, 316. 68. Doltsinis, N. L.; Marx, D. Nonadiabatic Car-Parrinello Molecular Dynamics. Phys. Rev. Lett. 2002, 88, 166402. 69. Castro, A.; Appel, H.; Oliveira, M.; Rozzi, C. A.; Andrade, X.; Lorenzen, F.; Marques, M. A. L.; Gross, E. K. U.; Rubio, A. Octopus: A Tool for the Application of TimeDependent Density Functional Theory. Phys. Status Solidi B 2006, 243, 2465-2488. 70. Du, L.; Lan, Z. An On-the-Fly Surface-Hopping Program JADE for Nonadiabatic Molecular Dynamics of Polyatomic Systems: Implementation and Applications. J. Chem. Theory Comput. 2015, 11, 1360-1374. 71. Tully, J. C. Nonadiabatic Molecular Dynamics. Int. J. Quantum Chem. 1991, 40, 299-309. 72. Akimov, A. V.; Prezhdo, O. V. The Pyxaid Program for Non-Adiabatic Molecular Dynamics in Condensed Matter Systems. J. Chem. Theory Comput. 2013, 9, 4959-4972. 73. Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648. 74. Parr, R. G.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, 1989. 75. Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields. J. Phys. Chem. 1994, 98, 11623-11627. 76. Hehre, W. J. Self-Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian-Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. J. Chem. Phys. 1972, 56, 2257. 77. Hariharan, P. C.; Pople, J. A. The Influence of Polarization Functions on Molecular Orbital Hydrogenation Energies. Theor. Chim. Acta 1973, 28, 213-222. 78. Roy, L. E.; Hay, P. J.; Martin, R. L. Revised Basis Sets for the LANL Effective Core Potentials. J. Chem. Theory Comput. 2008, 4, 1029-1031. 79. Hay, P. J.; Wadt, W. R. Ab Initio Effective Core Potentials for Molecular Calculations. Potentials for the Transition Metal Atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270. 80. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; et al. Gaussian, Inc.: Willingford CT, 2009. 81. Casida, M. E. Time-Dependent Density-Functional Theory for Molecules and Molecular Solids. J. Mol. Struct.: THEOCHEM 2009, 914, 3-18. 82. Casida, M. E.; Jamorski, C.; Casida, K. C.; Salahub, D. R. Molecular Excitation Energies to High-Lying Bound States from Time-Dependent Density-Functional Response 

ACS Paragon Plus Environment

22

Page 23 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Theory: Characterization and Correction of the Time-Dependent Local  Density Approximation Ionization Threshold. J. Chem. Phys. 1998, 108, 4439. 83. Runge, E.; Gross, E. K. U. Density-Functional Theory for Time-Dependent Systems. Phys. Rev. Lett. 1984, 52, 997. 84. Stratmann, R. E.; Scuseria, G. E.; Frisch, M. J. An Efficient Implementation of Time-Dependent Density-Functional Theory for the Calculation of Excitation Energies of Large Molecules. J. Chem. Phys. 1998, 109, 8218. 85. Martin, R. L. Natural Transition Orbitals. J. Chem. Phys. 2003, 118, 4775. 86. Slater, J. C.; Koster, G. F. Simplified LCAO Method for the Periodic Potential Problem. Physical Review 1954, 94, 1498-1524. 87. Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-Consistent-Charge Density-Functional Tight-Binding Method for Simulations of Complex Materials Properties. Phys. Rev. B 1998, 58, 7260-7268. 88. Elstner, M. SCC-DFTB: What Is the Proper Degree of Self-Consistency?. J. Phys. Chem. A 2007, 111, 5614-5621. 89. Otte, N.; Scholten, M.; Thiel, W. Looking at Self-Consistent-Charge Density Functional Tight Binding from a Semiempirical Perspective. J. Phys. Chem. A 2007, 111, 57515755. 90. Koskinen, P.; Mäkinen, V. Density-Functional Tight-Binding for Beginners. Comput. Mater. Sci. 2009, 47, 237-253. 91. Porezag, D.; Frauenheim, T.; Köhler, T.; Seifert, G.; Kaschner, R. Construction of Tight-Binding-Like Potentials on the Basis of Density-Functional Theory: Application to Carbon. Phys. Rev. B 1995, 51, 12947-12957. 92. Frauenheim, T.; Seifert, G.; Elstner, M.; Niehaus, T.; Köhler, C.; Amkreutz, M.; Sternberg, M.; Hajnal, Z.; Carlo, A. D.; Suhai, S. Atomistic Simulations of Complex Materials: Ground-State and Excited-State Properties. J. Phys.: Condens. Matter 2002, 14, 3015-3047. 93. Nosé, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511. 94. Hünenberger, P. H. Thermostat Algorithms for Molecular Dynamics Simulations. In Advanced Computer Simulation. Holm, C.; Kremer, K.; Eds. Advances in Polymer Science. Springer Berlin Heidelberg: Berlin, Heidelberg, Germany, 2005, 173, 105-149. 95. Martyna, G. J.; Klein, M. L.; Tuckerman, M. Nosé–Hoover Chains: The Canonical Ensemble via Continuous Dynamics. J. Chem. Phys. 1992, 97, 2635. 96. Martyna, G. J.; Tuckerman, M. E.; Tobias, D. J.; Klein, M. L. Explicit Reversible Integrators for Extended Systems Dynamics. Mol. Phys. 1996, 87, 1117-1157. 97. Hoover, W. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695-1697. 98. Wolfsberg, M.; Helmholz, L. The Spectra and Electronic Structure of the Tetrahedral Ions MnO4−, CrO4−−, and ClO4−. J. Chem. Phys. 1952, 20, 837. 99. Ammeter, J. H.; Buergi, H. B.; Thibeault, J. C.; Hoffmann, R. Counterintuitive Orbital Mixing in Semiempirical and Ab Initio Molecular Orbital Calculations. J. Am. Chem. Soc. 1978, 100, 3686-3692.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 26

100. Whangbo, M.-H.; Hoffmann, R. Counterintuitive Orbital Mixing. J. Chem. Phys. 1978, 68, 5498. 101. Jorge, F. E.; Giambiagi, M.; Degiambiagi, M. S. About Counterintuitive Orbital Mixing and Bond Populations. Theor. Chim. Acta 1983, 63, 529-539. 102. Alvarez, S. Table of Parameters for Extended Hueckel Calculations Universitat de Barcelona, 1989. 103. Hoff, D. A.; da Silva, R.; Rego, L. G. C. Coupled Electron–Hole Quantum Dynamics on D−π–A Dye-Sensitized TiO2 semiconductors. J. Phys. Chem. C 2012, 116, 2116921178. 104. Abuabara, S. G.; Rego, L. G.; Batista, V. S. Influence of Thermal Fluctuations on Interfacial Electron Transfer in Functionalized TiO2 Semiconductors. J. Am. Chem. Soc. 2005, 127, 18234-18242. 105. Bauernschmitt, R.; Ahlrichs, R. Treatment of Electronic Excitations within the Adiabatic Approximation of Time Dependent Density Functional Theory. Chem. Phys. Lett. 1996, 256, 454-464. 106. Scalmani, G.; Frisch, M. J.; Mennucci, B.; Tomasi, J.; Cammi, R.; Barone, V. Geometries and Properties of Excited States in the Gas Phase and in Solution: Theory and Application of a Time-Dependent Density Functional Theory Polarizable Continuum Model. J. Chem. Phys. 2006, 124, 94107. 107. Petersilka, M.; Grossmann, U. J.; Gross, E. K. U. Excitation Energies from TimeDependent Density-Functional Theory. Phys. Rev. Lett. 1996, 76, 1212. 108. Furche, F.; Ahlrichs, R. Adiabatic Time-Dependent Density Functional Methods for Excited State Properties. J. Chem. Phys. 2002, 117, 7433. 109. Song, H. E.; Taniguchi, M.; Diers, J. R.; Kirmaier, C.; Bocian, D. F.; Lindsey, J. S.; Holten, D. Linker Dependence of Energy and Hole Transfer in Neutral and Oxidized Multiporphyrin Arrays. J. Phys. Chem. B 2009, 113, 16483-16493. 110. Song, H. E.; Taniguchi, M.; Kirmaier, C.; Bocian, D. F.; Lindsey, J. S.; Holten, D. Probing Ground-State Hole Transfer between Equivalent, Electrochemically Inaccessible States in Multiporphyrin Arrays Using Time-Resolved Optical Spectroscopy. Photochem. Photobiol. 2009, 85, 693-704. 111. Casida, M. E. Time-Dependent Density Functional Response Theory for Molecules. In Recent Advances in Computational Chemistry; 1st ed.; Chong, D. E., Eds. World Scientific: Singapore, 1995, 1, p 155-192. 112. Hoffmann, R. An Extended Hückel Theory. I. Hydrocarbons. J. Chem. Phys. 1963, 39, 1397. 113. Hoffmann, R. Extended Hückel Theory. II. σ Orbitals in the Azines. J. Chem. Phys. 1964, 40, 2745. 114. Hoffmann, R. Extended Hückel Theory. III. Compounds of Boron and Nitrogen. J. Chem. Phys. 1964, 40, 2474. 115. Hoffmann, R. Extended Hückel Theory. IV. Carbonium Ions. J. Chem. Phys. 1964, 40, 2480. 116. Hoffmann, R. Extended Hückel Theory – V. Tetrahedron 1966, 22, 521-538. 117. Hoffmann, R. Extended Hückel Theory – VI. Tetrahedron 1966, 22, 539-545.

ACS Paragon Plus Environment

24

Page 25 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

118. Boer, F. P.; Newton, M. D.; Lipscomb, W. N. Extended Hückel Theory and Molecular Hartree-Fock SCF Theory. Proc. Natl. Acad. Sci. 1964, 52, 890-893. 119. Fitzpatrick, N. J.; Murphy, G. H. Double Zeta D Radial Wave Functions for Transition Elements. Inorg. Chim. Acta 1986, 111, 139-140. 120. Blyholder, G.; Coulson, C. A. Basis of Extended Hückel Formalism. Theor. Chim. Acta 1968, 10, 316-324. 121. Slater, J. Atomic Shielding Constants. Phys. Rev. 1930, 36, 57-64. 122. Zener, C. Analytic Atomic Wave Functions. Phys. Rev. 1930, 36, 51-56. 123. Phillips, H.; Zheng, Z.; Geva, E.; Dunietz, B. D. Orbital Gap Predictions for Rational Design of Organic Photovoltaic Materials. Org. Electron. 2014, 15, 1509-1520. 124. Stein, T.; Autschbach, J.; Govind, N.; Kronik, L.; Baer, R. Curvature and Frontier Orbital Energies in Density Functional Theory. J. Phys. Chem. Lett. 2012, 3, 3740-3744.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 26

TOC Figure

ACS Paragon Plus Environment

26