Photochromic Mechanism of a Bridged Diarylethene: Combined

Jan 4, 2017 - *Phone: 86-010-58806770. ... Intramolecularly bridged diarylethenes exhibit improved photocyclization quantum yields because the anti-sy...
0 downloads 0 Views 5MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Photochromic Mechanism of a Bridged Diarylethene: Combined Electronic Structure Calculations and Nonadiabatic Dynamics Simulations Yating Wang, Yuan-Jun Gao, Qian Wang, and Ganglong Cui J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.6b11682 • Publication Date (Web): 04 Jan 2017 Downloaded from http://pubs.acs.org on January 5, 2017

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 32

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

Photochromic Mechanism of a Bridged Diarylethene: Combined Electronic Structure Calculations and Nonadiabatic Dynamics Simulations Ya-Ting Wang, Yuan-Jun Gao, Qian Wang, and Ganglong Cui∗ Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China E-mail: [email protected]

Abstract Intramolecularly bridged diarylethenes exhibit improved photocyclization quantum yields because the anti-syn isomerization that originally suppresses photocyclization in classical diarylethenes is blocked. Experimentally, three possible channels have been proposed to interpret experimental observation, but many details of photochromic mechanism remain ambiguous. In this work we have employed a series of electronic structure methods (OM2/MRCI, DFT, TDDFT, RI-CC2, DFT/MRCI, and CASPT2) to comprehensively study excited state properties, photocyclization, and photoreversion dynamics of 1,2-dicyano[2,2]metacyclophan-1-ene. On the basis of optimized stationary points and minimum-energy conical intersections, we have refined experimentally proposed photochromic mechanism. Only an S1 /S0 minimum∗ To

whom correspondence should be addressed

1 ACS Paragon Plus Environment

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

energy conical intersection is located; thus, we can exclude the third channel experimentally proposed. In addition, we find that both photocyclization and photoreversion processes use the same S1 /S0 conical intersection to decay the S1 system to the S0 state, so we can unify the remaining two channels into one. These new insights are verified by our OM2/MRCI nonadiabatic dynamics simulations. The S1 excited-state lifetimes of photocyclization and photoreversion are estimated to be 349 and 453 fs, respectively, which are close to experimentally measured values: 240±60 fs and 250 fs in acetonitrile solution. The present study not only interprets experimental observations and refines previously proposed mechanism, but also provides new physical insights that are valuable for future experiments.

Introduction Photochromism is a light-induced reversible chemical transformation between two stable forms of a compound with different absorption wavelengths. 1 This kind of reactions are of great significance due to their actual implications for the development of photomemories, photoswitches, and data storage materials. 2,3 Among available photochromic families, diarylethenes are a class of the most efficient photoswitch compounds because of their high thermodynamic stability, high fatigue resistance, and good photoconversion yield. 4,5 For this reason numerous advances on diarylethenes in various environments have been reported, including chemical synthesis of various derivatives, multi-photon absorption, luminescence, Förster resonance energy transfer, magnetic properties, electronic transport properties, laser flash photolysis, photocyclization dynamics, etc. 6–41 However, photocyclization yields of traditional diarylethenes are not satisfactory because of the co-existence of open-form anti and syn configurations (only the latter is active for photocyclization, see Eq. 1). In addition, both isomerization and cyclization reactions can take place from the syn configuration; thus, photocyclization is to certain extent suppressed by the competing photoisomerization. 1,4 2 ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

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 1: Photochemical reactions in classical diarylethenes: isomerization between synand anti-open configurations, and photocyclization and photoreversion (ring-opening) between syn-open and syn-closed configurations.

Figure 2: Bridged photochromic diarylethene (meta-cyclophan-1-ene) synthesized by Takeshita et al. 42–44 In this series of compounds, two reacting aryl groups are maintained in an antiparallel conformation by the presence of an additional alkyl bridge of variable length n.

3 ACS Paragon Plus Environment

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

In order to improve their photocyclization yields, Takeshita et al. successfully synthesized a series of meta-cyclophan-1-ene photochromic compounds (Fig. 2) in which the two reacting aryl groups are maintained in an antiparallel conformation by an additional alkyl bridge. 42–44 In their early work, only a marginal improvement of photocyclization yield was reported. 45 To understand the reason for this moderate improvement and to figure out parameters controlling photoswitch efficiency, Aloïse et al. re-investigated photochemical reactivities of some bridged diarylethenes in pure acetonitrile and mixed hexane-ethanol (12:1) solutions using stationary and femtosecond transient absorption spectroscopy. 27 Many spectroscopic properties were obtained, as well as photochromic quantum yields at different wavelengths and solutions. Interestingly, the photocyclization quantum yield of a bridged diarylethene is close to unity when the wavelength is longer than 400 nm, but decreases to 0.69 at 313 nm in acetonitrile solution. Furthermore, in terms of time-resolved absorption spectra, time-dependent density functional theory (TD-DFT), and semiempirical configuration interaction with single excitation (AM1/CIS) calculations, they proposed three possible photocyclization pathways to interpret their experimental observation. In the first channel, Eq. 1, the S1 state is directly populated in the Franck-Condon region, followed by an ultrafast relaxation (less than 100 fs) to the open-form S1 minimum [S∗1 (OF)]; then the system decays to the S0 state via a conical intersection [CI1 (S1 /S0 )] within less than 120 fs, eventually arriving at the closed-form [S0 (CF)]. In the second channel, Eq. 2, the higher excited state [SFC n (OF)] is first populated and then this state quickly decays to the open-form hot S1 minimum [S∗1 (OF)]. After overcoming a transition state [S1 (TS)], the system arrives at a closed-form S1 minimum [S1 (CF)], followed by internal conversion to the S0 state [S0 (CF)]. Finally, they proposed the third channel (Eq. 3) to be as an alternative to the second channel because there is no clear isosbestic point involved in the closed-form visible band observed in their transient spectra. In such channel, another S1 /S0 conical intersection [CI2 (S1 /S0 )] directly funnels the open-form hot S1 minimum [S∗1 (OF)] to the closed-form S0 minimum [S0 (CF)]. Nevertheless, one can easily find that

4 ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

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

many mechanistic details at the atomistic level are elusive in these three pathways, for example, what do these two different S1 /S0 conical intersections look like? etc.

∗ SFC 1 (OF) → S1 (OF) → CI1 → S0 (CF) → S0 (CF)

(1)

∗ SFC n (OF) → S1 (OF) → S1 (TS) → S1 (CF) → S0 (CF)

(2)

∗ SFC n (OF) → S1 (OF) → CI2 → S0 (CF)

(3)

Computer simulations are complementary to experiments. They can provide invaluable insights and give more mechanistic details that cannot be gained by experiments alone. Theoretical computations on photocyclization and photoreversion reactions of bridged diarylethenes are rarely reported. To our best knowledge, there are merely a few AM1/CIS and TD-DFT electronic structure calculations by Aloïse et al., 27 who aimed at interpreting experimentally observed time-resolved transient absorption spectra. Even though the AM1/CIS computations already gave insights into the photochromic mechanism of bridged diarylethenes, many atomistic details remain ambiguous. Hence further electronic structure calculations and nonadiabatic dynamics simulations are very meaningful. In this work, we employed a series of electronic structure methods i.e. TDDFT, RI-CC2, DFT/MRCI, and CASPT2 to re-investigate the S1 photocyclization and photoreversion processes of one bridged diarylethene studied experimentally. 27 Subsequently, we carried out trajectory-based nonadiabatic surface-hopping dynamics simulations to verify our proposed mechanism.

5 ACS Paragon Plus Environment

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

Figure 3: Photocyclization and photoreversion (ring-opening) reactions of the bridged diarylethene proposed by Aloïse et al. 27 The corresponding π orbitals are also shown.

Simulation Details Ab Initio and Density Functional Methods Stationary points of the S0 state were optimized using the B3LYP method. 46–48 Vertical excitation energies were computed using the TD-DFT method with different exchangecorrelation functionals (B3LYP, 46–48 CAM-B3LYP, 49 ω B97XD, 50 M06-2X 51 ), the coupled cluster theory (RI-CC2) method, 52–54 the complete active space second-order perturbation theory (CASPT2) method, 55 and the DFT-based multireference configuration interaction (DFT/MRCI) method. 56 In the DFT/MRCI calculations, the BHLYP functional 57 was used in order to match with the original parametrization. To achieve more accurate potential energy profiles, the CASPT2 method was employed to refine single-point energies of all optimized structures. In the CASPT2//CASSCF calculations, 10 π electrons were placed into the active space that is comprised of 10 π and π ∗ orbitals; three roots with equal weights were used as well as the cholesky decomposition technique 58 and the imaginary shift technique (0.1 a.u. in our calculations). 59 The following basis sets and packages were used: cc-pVDZ for CASPT2//CASSCF (MOLCAS7.6), 60,61 SVP for RI-CC2 (TURBOMOLE6.3), 62,63 TZVP for DFT and TDDFT (GAUSSIAN09), 64,65 and TZVP for DFT/MRCI (TURBOMOLE5.7.1). 56,64,66

6 ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32

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

Semiempirical Methods All semiempirical calculations (geometric optimizations and nonadiabatic dynamics) were performed using the OM2/MRCI method implemented in the MNDO99 program. 67–72 The restricted open-shell Hartree-Fock formalism was applied in the self-consistent field (SCF) treatment (i.e., the orbitals were optimized for the leading electronic configuration of the S1 state with two singly occupied π orbitals). The active space in the MRCI calculations included six electrons in six orbitals: in terms of the SCF configuration, it comprised the two highest doubly occupied π orbitals, the two singly occupied π orbitals, and the two lowest unoccupied π ∗ orbitals. For the MRCI treatment, three configuration state functions were chosen as references, namely the SCF configuration and the other two closed-shell configurations derived therefrom (i.e., all singlet configurations that can be generated from HOMO and LUMO of the closed-shell ground state). The MRCI wavefunction was built by allowing all single and double excitations from these three references. During OM2/MRCI geometry optimizations and dynamics simulations, all required gradients and nonadiabatic coupling elements were computed analytically. 67–69

Nonadiabatic Dynamics Simulations Nonadiabatic dynamics simulations with Tully’s fewest-switches surface hopping (FSSH) scheme 73,74 were carried out at the semiempirical OM2/MRCI level. 67–69,75,76 The initial atomic coordinates and velocities were randomly selected from 5 ps trajectories of ground-state Born-Oppenheimer molecular dynamics. The number of excited-state dynamics runs for each chosen snapshot was chosen according to the computed S0 → S1 transition probabilities. A total of 272 (142 for photocyclization and 130 for photoreversion) surface-hopping trajectories were successfully finished and thus evaluated in the present nonadiabatic OM2/MRCI simulations, with all relevant energies, gradients, and nonadiabatic coupling elements computed “on the fly” when needed. The fewest-switches

7 ACS Paragon Plus Environment

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

criterion of Tully et al. was applied to decide when to hop and how to adjust the related velocities. 73,74 Time step was chosen to be 0.2 fs for nuclear motion and 0.002 fs for electronic propagation. The unitary propagator evaluated at middle-point was used to propagate electronic motion. 75,76 Translation and rotation motions were removed in each step. The empirical decoherence correction (0.1 in our simulations) proposed by Granucci et al. was employed. 77 Further technical details were given in previous publications. 78–92

Results and Discussion Open- and Closed-Form S0 Minima The bridged diarylethene has two stable conformers i.e. open and closed forms as shown in Fig. 3. Its two reacting phenyl groups are linked by the −C2 H4 − group. This chemical modification removes the anti configuration that originally exists in nonbridged diarylethenes (e.g. stilbene in Fig. 1). As a result, It blocks the anti-syn isomerization, in such way enhancing the quantum yield of photocyclization. Fig. 4 shows two stable conformers optimized by the OM2/MRCI method, which are referred to as S0-OF and S0-CF. In S0-OF, the C1 and C2 atoms are not bonded (2.457 Å); whereas, this bond is formed in S0-CF (1.530 Å). At the CASPT2/cc-pVDZ and OM2/MRCI levels, the potential energy of S0-CF is predicted to be 8.0 and 14.5 kcal/mol higher than that of S0-OF, respectively.

Vertical Excitation Energies The computed S0 → S1 vertical excitation energies of the open- and closed-form bridged diarylethene in vacuo are compiled in Table 1. At S0-OF, TD-B3LYP (2.96 eV), DFT/MRCI (3.17 eV), and AM1/CIS (3.20 eV) to different extent underestimate the S0 → S1 vertical excitation energy compared with CASPT2 (3.59 eV), CC2 (3.49 eV), TD-M062X (3.48 eV), TD-CAM-B3LYP (3.43 eV), and OM2/MRCI (3.87 eV). This implies that the S1 state 8 ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32

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

4

5

3

6

2 457

3 1.5

. 12

S0-OF

1.66

0

S0-CF

28

7

1.7

S1-CF

S1S0-CF

Figure 4: S0 and S1 minima and S1 /S0 conical intersection optimized at the OM2/MRCI level (OF: open form; CF: closed form). Also shown are the C1-C2 bond lengths. See supporting information for Cartesian coordinates.

9 ACS Paragon Plus Environment

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 32

has partial charge-transfer character because the conventional and non-range-separated functionals fail to give good prediction for charge-transfer states. At S0-CF, TD-B3LYP (1.97 eV) and DFT/MRCI (2.17 eV) still underestimate in comparison with CASPT2 (2.51 eV); but, AM1/CIS (2.33 eV) is close to OM2/MRCI (2.28 eV), TD-CAM-B3LYP (2.26 eV) and CC2 (2.28 eV). Table 1: Vertical S0 → S1 Excitation Energies (eV) of the Open- (first row) and ClosedForm (second row) Bridged Diarylethene in Vacuo Calculated by Various Electronic Structure Methods† OM2/MRCI CAM-B3LYP M062X CC2a DFT/MRCIa AM1/CISb 3.87 3.43 3.48 3.49 3.17 3.20 2.28 2.26 2.24 2.28 2.17 2.33 † TDDFT computations use the B3LYP/TZVP optimized structures; a TZVP basis set with CC2/SVP optimized structures; b from the work of Aloïse et al. (AM1-CIS//AM1); c CASPT2(10,10)/cc-pVDZ with the OM2/MRCI optimized structures.

CASPT2c 3.59 2.51

Table 2: Single-Point Energies Calculated by OM2/MRCI (first row) and SA3CASPT2(10,10)/cc-pVDZ//OM2/MRCI (second row) Methods S0-OF S1-OF-FC 0.0 89.2 0.0 82.8

S0-CF 14.5 8.0

S1-CF-FC 64.7 65.0

S1-CF 49.3 53.3

S1S0-CF(S0/S1) 59.0/59.0 55.2/61.1

Table 3: CASPT2/cc-pVDZ Calculated Dipole Moments (Debye) at the Franck-Condon Points of Open- (column 2-3) and Closed-Form (column 4-5) (i.e. S0-OF and S0-CF) structure S0-OF S0-CF electronic state S0 S1 S0 S1 in vacuoa 8.1 13.1 9.5 8.7 a OM2/MRCI optimized structure; To understand electronic structures of the S1 state at the open- and closed-form Franck-Condon points i.e. S0-OF and S0-CF, we analyzed their gas-phase S1 electronic wavefunctions at the CASPT2 level, which can be expressed:

ψS1 (OF) = −0.903719ΦHOMO→LUMO + 0.134895ΦHOMO→LUMO+2 + ... 10 ACS Paragon Plus Environment

(4)

Page 11 of 32

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

ψS1 (CF) = −0.887181ΦHOMO→LUMO + 0.073979ΦHOMO→LUMO+1 + ...

(5)

where Φ represents an electronic configuration, e.g. ΦHOMO→LUMO corresponds to a single-excitation configuration due to an electron excited from HOMO to LUMO. It can be seen that the S1 electronic state is mainly comprised of the HOMO→LUMO electronic configuration at both structures. In acetonitrile solution, the ΦHOMO→LUMO electronic configuration still contributes overwhelmingly to the S1 state. In addition, we have plotted HOMO and LUMO (Fig. 5). At S0-OF, the HOMO→ LUMO electron transition causes partial charge transfer. In HOMO, the center of electronic cloud is located around the C1 and C2 atoms, whereas in LUMO, it moves toward the C4 and C5 atoms. This electron redistribution at S0-OF is also reflected in the CASPT2 computed dipole moments in vacuo/acetonitrile solution: 8.1/10.7 Debye [S0 ] vs. 13.1/16.4 Debye [S1 ]. By contrast, there just exists little charge transfer in the S1 state at S0-CF (see Table 3).

S1 Minimum and S1 /S0 Conical Intersection We tried to optimize the open-form S1 minimum from the open-form S0 minimum at the OM2/MRCI level. However, all attempts converge toward a closed-form minimum with the C1-C2 distance of 1.667 Å. This structure is actually the closed-form S1 minimum, which is referred to as S1-CF in Fig. 4. This implies that the relaxation process from the opento closed-form is essentially barrierless in the S1 state. At the OM2/MRCI level, we have also optimized an S1 /S0 minimum-energy conical intersection, which is referred to as S1S0-CF in Fig. 4. In terms of the C1-C2 reaction coordinate, one can easily find that S1S0-CF (1.728 Å) is structurally inbetween the openform Franck-Condon point (2.457 Å, S0-OF) and the closed-form S1 minimum (1.667 Å, S1-CF). In other words, it is in the middle of the S1 relaxation process. In addition, S1S0-CF is energetically close to S1-CF: 65.0 [64.7] kcal/mol vs. 53.3 [49.3] kcal/mol

11 ACS Paragon Plus Environment

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 32

OF

HOMO

LUMO

CF

Figure 5: HOMO and LUMO molecular orbitals (isosurface: 0.02) responsible for the S0 → S1 electronic transition at the open- and closed-form Franck-Condon points, respectively (i.e. S0-OF and S0-CF).

Figure 6: Time-dependent population of the S1 and S0 states for the open-form→ closedform photocyclization and closed-form→ open-form photoreversion processes of the bridged diarylethene. See text for discussion.

12 ACS Paragon Plus Environment

Page 13 of 32

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

at the CASPT2/cc-pVDZ [OM2/MRCI] level. Therefore, this S1 /S0 conical intersection is accessible both structurally and energetically and should play an important role for the photocyclization and photoreversion of the bridged diarylethene (see below).

Photocyclization and Photoreversion Dynamics To explore the S1 excited-state lifetime of the open-form bridged diarylethene, we have carried out 1 ps OM2/MRCI nonadiabatic dynamics simulations (142 successful runs). In the end of 1 ps simulations, 96.5% trajectories decay to the S0 state; 3.5% ones still survive in the S1 state. In terms of our trajectories, the photocyclization quantum yield is estimated to be 0.47, which is close to the experimentally measured value of 0.49 in methanol upon irradiation of 313 nm. 27 The left panel of Fig. 6 illustrates the time-dependent S1 and S0 populations. Two curves intersect at 242 fs, based on which the S1 lifetime is estimated to be 349 fs (see the fitting model in our previous work 78 ). Our in-vacuo prediction is a little longer than experimentally measured 240±60 fs in acetonitrile solution. 27 This difference could be due to the lacking of solvent effects in our dynamics simulations.

Figure 7: Distribution of the C1-C2 distance at all S1 →S0 hopping points of both (left) photocyclization and (right) photoreversion processes. The distribution of the C1-C2 distance (reaction coordinate of photocyclization) at the 13 ACS Paragon Plus Environment

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

S1 → S0 hopping points is shown in the left panel of Fig. 7. Most of the S1 → S0 hoppings take place around the S1 /S0 minimum-energy conical intersection S1S0-CF (1.728 Å), which proves that only a hopping region that centers around S1S0-CF is responsible for the S1 decay of the trajectories starting from the open-form bridged diarylethene. This dynamical feature is well consistent with the topological structure of the S0 and S1 potential energy surfaces (see above).

Figure 8: Time-dependent C1-C2 distance in four typical trajectories starting from the open form. Also shown are the hopping times (numbers, in fs). Four typical trajectories were seen in our dynamics simulations (see Fig. 8). In Trajectory_33, the system quickly relaxes to the S1 minimum (S1-CF) from the open-form Franck-Condon point S0-OF within about 50 fs; then, it oscillates around S1-CF until 146 fs, followed by an S1 → S0 de-excitation; finally, it reversely goes to the open-form S0 minimum (S0-OF). Trajectory_39 is similar to Trajectory_33 except that the system finally goes 14 ACS Paragon Plus Environment

Page 14 of 32

Page 15 of 32

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

to the closed-form S0 minimum (S0-CF). However, in Trajectory_16 and Trajectory_173, the system oscillates around S1-CF for a longer time before it decays to the S0 state (762 fs in Trajectory_16 and 628 fs in Trajectory_173). In addition to photocyclization, the photoreversion process from the closed-form bridged diarylethene is also simulated using the OM2/MRCI nonadiabatic dynamics simulations (130 runs with 1 ps time). The S1 excited-state lifetime, based on the time-dependent S1 and S0 populations in the right panel of Fig. 6, is estimated to be around 453 fs, which is a little bigger than the experimentally estimated 250 fs in acetonitrile solution. 27 Different from the photocyclization situation, 80% trajectories decay to the S0 state in the end of 1 ps simulation time. It can be seen that the S1 excited-state lifetime for the photoreversion reaction is a little longer than that for the photocyclization reaction. This difference could be interrelated with the energy of the S1 /S0 minimum-energy conical intersection (S1S0-CF) relative to their respective starting structures in the reaction pathways. According to our computations, S1S0-CF is 21.7 and 3.9 kcal/mol lower than S1-OF-FC and S1-CF-FC (see Table 2). With this in mind, it is rational that accessing to S1S0-CF from S1-OF-FC is faster than that from S1-CF-FC. The distribution of the C1-C2 distance at all S1 → S0 hopping points is shown in the right panel of Fig. 7. Similarly, the S1 → S0 radiationless transition mainly takes place near the same S1 /S0 minimum-energy conical intersection S1S0-CF. There exists as well four typical trajectories for the photoreversion process (see Fig. 9). Trajectory_21 and Trajectory_68 relax quickly to the S1 minimum (S1-CF) within 20 fs. After several periods of oscillation around S1-CF, the system decays to the S0 state (102 fs for Trajectory_21 and 182 fs for Trajectory_68), finally decaying to the open- and closed-form S0 minima S0-OF and S0-CF, respectively. In comparison, Trajectory_8 and Trajectory_58 oscillate around S1-CF much longer before they decay to the S0 state.

15 ACS Paragon Plus Environment

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

Figure 9: Time-dependent C1-C2 distance in four typical trajectories starting from the closed form. Also shown are the hopping times (numbers, in fs).

16 ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

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

Mechanism In terms of our present static electronic structure calculations and nonadiabatic dynamics simulations, we suggest the following S1 relaxation mechanism for the photocyclization and photoreversion processes of the bridged diarylethene (Fig. 10). It should be noted that both processes adopt a similar photochromic mechanism in the S1 state except starting from different Franck-Condon regions. Specifically, the photocyclization starts from the open-form conformer, while the photoreversion does from the closed-form one. For the photocyclization or photoreversion process, the S1 state is directly populated at S0OF or S0-CF (of course, the S1 state can be indirectly populated via internal conversion from higher excited states; not considered in our work); then, the system quickly relaxes to S1-CF; finally, it decays to the S0 state around S1S0-CF. Once decaying to the S0 state at this S1 /S0 conical intersection, the system bifurcates into either S0-OF or S0-CF.

S1-OF-FC (82.8) S1S0-CF (61.1)

S1-CF-FC (65.0) S1-CF (53.3)

S0-CF(8.0) S0-OF(0.0) Figure 10: Proposed photocyclization and photoreversion mechanism for the bridged diarylethene system. Energies are refined by the SA3-CASPT2(10,10)/cc-pVDZ method. The blue lines represent the S1 and S0 adiabatic potential energy surfaces, while the black one represents the S1 /S0 crossing seam. Our suggested S1 photocyclization and photoreversion mechanism overall matches with experimentally proposed one (Eq. 1). However, we did not find an open-form S1 min17 ACS Paragon Plus Environment

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

imum at the OM2/MRCI level and all optimizations starting from the open-form structure converge toward the closed-form S1 minimum (S1-CF). This essentially barrierless S1 relaxation from the open- to closed-form region is also confirmed by our OM2/MRCI-based nonadiabatic dynamic simulations (Fig. 8). In all four typical photocyclization trajectories, the C1-C2 distance rapidly decreases to ca. 1.6 Å from the initial ca. 2.5 Å within 50 fs. Actually, in the work of Aloïse and coworkers, they have merely estimated a negligible barrier of 0.7 kcal/mol at the AM1/CIS level. 27 It is also helpful to note that in many similar 6π electrocyclic systems, this kind of photocyclization reaction is always barrierless. Moreover, Aloïse et al. also proposed another two indirect S1 photocyclization pathways where the S1 state is populated via internal conversion from higher excited states. 27 In the former, Eq. 2, a transition state was suggested. However, this viewpoint is not supported by our present work: in static electronic structure calculations, we do not locate such kind of transition state; in dynamics simulations, such region is also not encountered. The latter, Eq. 3, can also be excluded because we only find a structurally and energetically allowed S1 /S0 conical intersection at the OM2/MRCI level. The most convincing evidence is that in our photocyclization and photoreversion dynamics simulations, only a hopping region that centers around the S1 /S0 conical intersection S1S0-CF is observed (see discussion above and Fig. 7). Taken these reasons, it is clear that these two indirect S1 relaxation channels are in essence the same as the S1 photocyclization and photoreversion mechanism in Fig. 10. In comparison with previous experiments and computations, 27 our current work not only refines previous mechanism but also provides new insights. Due to the impurity of initial samples and the nonnegligible contribution from the open-form bridged diarylethene, it was experimentally difficult to precisely assign the characteristic time for the S1 (CF)→ S0 (CF) internal conversion. 27 Instead, experiments only estimated an upper limit of 250 fs, which is close to our estimated excited-state lifetime of 453 fs. In addition, we have demonstrated that both photocyclization and photoreversion processes share an ener-

18 ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32

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

getically easily accessible conical intersection namely S1S0-CF to funnel the S1 excited state. In combination with that the S1 relaxation process from both S0-OF and S0-CF is barrierless, the photocyclization and photoreversion virtually adopt a similar mechanism, just with different initial conditions (see Fig. 10).

Conclusions By using electronic structure calculations and nonadiabatic dynamics simulations we have for the first time explored the S1 photocyclization and photoreversion dynamics of a recently studied bridged diarylethene. 27 In the first phase, we have characterized the S0 and S1 minima and the S1 /S0 minimum-energy conical intersection, and mapped the relevant S1 and S0 potential energy surfaces, based on which the photocyclization and photoreversion mechanism was suggested. In the second phase, we verified our proposed mechanism and predicted the S1 decay times with resort to trajectory-based nonadiabatic dynamics simulations. The estimated excited-state lifetimes are in good agreement with experimentally measured values. The present simulations interpret experimental observations and refine previously proposed pathways. Furthermore, although present electronic structure calculations and dynamics simulations provide new important insights into the mechanistic photochromism of bridged diarylethenes, further experimental and theoretical studies are still needed to understand wavelength dependent photocyclization quantum yields and solvent effects.

Supporting Information Product analysis and Cartesian coordinates of all optimized structures.

19 ACS Paragon Plus Environment

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

Acknowledgments This work was supported by the National Natural Science Foundation of China (21522302, 21520102005, and 21421003); G.C. is also grateful for financial support from the "Recruitment Program of Global Youth Experts" and "Fundamental Research Funds for Central Universities".

Corresponding Author Phone number: 86-010-58806770. Email: [email protected] (G.C.)

References (1) Dürr, H.; Bouas-Laurent, H. Photochromism: Molecules and Systems; Elsevier, 2003. (2) Berkovic, G.; Krongauz, V.; Weiss, V. Spiropyrans and Spirooxazines for Memories and Switches. Chem. Rev. 2000, 100, 1741–1754. (3) Hampp, N. Bacteriorhodopsin as a Photochromic Retinal Protein for Optical Memories. Chem. Rev. 2000, 100, 1755–1776. (4) Irie, M. Photochromism: Memories and Switches-Introduction. Chem. Rev. 2000, 100, 1683–1684. (5) Irie, M.; Kobatake, S.; Horichi, M. Reversible Surface Morphology Changes of a Photochromic Diarylethene Single Crystal by Photoirradiation. Science 2001, 291, 1769–1772. (6) Matsuda, K.; Irie, M. A Diarylethene with Two Nitronyl Nitroxides: Photoswitching of Intramolecular Magnetic Interaction. J. Am. Chem. Soc. 2000, 122, 7195–7201.

20 ACS Paragon Plus Environment

Page 20 of 32

Page 21 of 32

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

(7) Yagi, K.; Soong, C. F.; Irie, M. Synthesis of Fluorescent Diarylethenes Having a 2,4,5-Triphenylimidazole Chromophore. J. Org. Chem. 2001, 66, 5419–5423. (8) Shibata, K.; Muto, K.; Kobatake, S.; Irie, M. Photocyclization/Cycloreversion Quantum Yields of Diarylethenes in Single Crystals. J. Phys. Chem. A 2002, 106, 209– 214. (9) Morimoto, M.; Kobatake, S.; Irie, M. Multicolor Photochromism of Two- and ThreeComponent Diarylethene Crystals. J. Am. Chem. Soc. 2003, 125, 11080–11087. (10) Boggio-Pasqua, M.; Ravaglia, M.; Bearpark, M. J.; Garavelli, M.; Robb, M. A. Can Diarylethene Photochromism Be Explained by a Reaction Path Alone? A CASSCF Study with Model MMVB Dynamics. J. Phys. Chem. A 2003, 107, 11139–11152. (11) Yam, V. W.-W.; Ko, C.-C.; Zhu, N. Y. Photochromic and Luminescence Switching Properties of a Versatile Diarylethene-Containing 1,10-Phenanthroline Ligand and Its Rhenium (I) Complex. J. Am. Chem. Soc. 2004, 126, 12734–12735. (12) Murakami, M.; Miyasaka, H.; Okada, T.; Kobatake, S.; Irie, M. Dynamics and Mechanisms of the Multiphoton Gated Photochromic Reaction of Diarylethene Derivatives. J. Am. Chem. Soc. 2004, 126, 14764–14772. (13) Frigoli, M.; Welch, C.; Mehl, G. H. Design of Mesomorphic Diarylethene-Based Photochromes. J. Am. Chem. Soc. 2004, 126, 15382–15383. (14) Tanifuji, N.; Irie, M.; Matsuda, K. New Photoswitching Unit for Magnetic Interaction: Diarylethene with 2,5-Bis(arylethynyl)-3-thienyl Group. J. Am. Chem. Soc. 2005, 127, 13344–13353. (15) Moriyama, Y.; Matsuda, K.; Tanifuji, N.; Irie, S.; Irie, M. Electrochemical Cyclization/Cycloreversion Reactions of Diarylethenes. Org. Lett. 2005, 7, 3315–3318.

21 ACS Paragon Plus Environment

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

(16) Saita, S.; Yamaguchi, T.; Kawai, T.; Irie, M. Two-Photon Photochromism of Diarylethene Dimer Derivatives. ChemPhysChem 2005, 6, 2300–2306. (17) Lee, P. H.-M.; Ko, C.-C.; Zhu, N. Y.; Yam, V. W.-W. Metal Coordination-Assisted NearInfrared Photochromic Behavior: A Large Perturbation on Absorption Wavelength Properties of N,N-Donor Ligands Containing Diarylethene Derivatives by Coordination to the Rhenium (I) Metal Center. J. Am. Chem. Soc. 2007, 129, 6058–6059. (18) Pu, S. Z.; Liu, G.; Shen, L.; Xu, J. K. Efficient Synthesis and Properties of Isomeric Photochromic Diarylethenes Having a Pyrrole Unit. Org. Lett. 2007, 9, 2139–2142. (19) Katsonis, N.; Minoia, A.; Kudernac, T.; Mutai, T.; Xu, H.; Uji-i, H.; Lazzaroni, R.; De Feyter, S.; Feringa, B. L. Locking of Helicity and Shape Complementarity in Diarylethene Dimers on Graphite. J. Am. Chem. Soc. 2008, 130, 386–387. (20) Indelli, M. T.; Carli, S.; Ghirotti, M.; Chiorboli, C.; Ravaglia, M.; Garavelli, M.; Scandola, F. Triplet Pathways in Diarylethene Photochromism: Photophysical and Computational Study of Dyads Containing Ruthenium (II) Polypyridine and 1,2-Bis(2methylbenzothiophene-3-yl)maleimide Units. J. Am. Chem. Soc. 2008, 130, 7286– 7299. (21) Yam, V. W.-W.; Lee, J. K.-W.; Ko, C.-C.; Zhu, N. Y. Photochromic DiaryletheneContaining Ionic Liquids and N-Heterocyclic Carbenes. J. Am. Chem. Soc. 2009, 131, 912–913. (22) Al-Atar, U.; Fernandes, R.; Johnsen, B.; Baillie, D.; Branda, N. R. A Photocontrolled Molecular Switch Regulates Paralysis in a Living Organism. J. Am. Chem. Soc. 2009, 131, 15966–15967. (23) Yokoyama, Y.; Shiozawa, T.; Tani, Y.; Ubukata, T. A Unified Strategy for Exceptionally High Diastereoselectivity in the Photochemical Ring Closure of Chiral Diarylethenes. Angew. Chem. Int. Ed. 2009, 48, 4521–4523. 22 ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32

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

(24) Gallardo, I.; Guirado, G.; Prats, G.; Takeshita, M. Electrochemically Triggered Conversion between Metacyclophan-1-ene and Dihydropyrene Molecular Switching Systems. Phys. Chem. Chem. Phys. 2009, 11, 1502–1507. (25) Liu, H.-H.; Chen, Y. The Photochromism and Fluorescence of Diarylethenes with a Imidazole Bridge Unit: A Strategy for the Design of Turn-on Fluorescent Diarylethene System. J. Phys. Chem. A 2009, 113, 5550–5553. (26) Uchida, K.; Nishikawa, N.; Izumi, N.; Yamazoe, S.; Mayama, H.; Kojima, Y.; Yokojima, S.; Nakamura, S.; Tsujii, K.; Irie, M. Phototunable Diarylethene Microcrystalline Surfaces: Lotus and Petal Effects upon Wetting. Angew. Chem. Int. Ed. 2010, 49, 5942–5944. (27) Aloïse, S.; Sliwa, M.; Pawlowska, Z.; Réhault, J.; Dubois, J.; Poizat, O.; Buntinx, G.; Perrier, A.; Maurel, F.; Yamaguchi, S. et al. Bridged Photochromic Diarylethenes Investigated by Ultrafast Absorption Spectroscopy: Evidence for Two Distinct Photocyclization Pathways. J. Am. Chem. Soc. 2010, 132, 7379–7390. (28) Singer, M.; Jäschke, A. Reversibly Photoswitchable Nucleosides: Synthesis and Photochromic Properties of Diarylethene-Functionalized 7-Deazaadenosine Derivatives. J. Am. Chem. Soc. 2010, 132, 8372–8377. (29) Maafi, M. The Potential of AB (1Φ) Systems for Direct Actinometry. Diarylethenes as Successful Actinometers for the Visible Range. Phys. Chem. Chem. Phys. 2010, 12, 13248–13254. (30) Mori, K.; Ishibashi, Y.; Matsuda, H.; Ito, S.; Nagasawa, Y.; Nakagawa, H.; Uchida, K.; Yokojima, S.; Nakamura, S.; Irie, M. et al. One-Color Reversible Control of Photochromic Reactions in a Diarylethene Derivative: Three-Photon Cyclization and Two-Photon Cycloreversion by a Near-Infrared Femtosecond Laser Pulse at 1.28

µ m. J. Am. Chem. Soc. 2011, 133, 2621–2625. 23 ACS Paragon Plus Environment

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

(31) Fukaminato, T.; Doi, T.; Tamaoki, N.; Okuno, K.; Ishibashi, Y.; Miyasaka, H.; Irie, M. Single-Molecule Fluorescence Photoswitching of a Diarylethene-Perylenebisimide Dyad: Non-Destructive Fluorescence Readout. J. Am. Chem. Soc. 2011, 133, 4984– 4990. (32) Uchida, K.; Yamanoi, Y.; Yonezawa, T.; Nishihara, H. Reversible On/Off Conductance Switching of Single Diarylethene Immobilized on a Silicon Surface. J. Am. Chem. Soc. 2011, 133, 9239–9241. (33) Nishi, H.; Asahi, T.; Kobatake, S. Enhanced One-Photon Cycloreversion Reaction of Diarylethenes near Individual Gold Nanoparticles. J. Phys. Chem. C 2011, 115, 4564–4570. (34) Chan, J. C.-H.; Lam, W. H.; Wong, H.-L.; Zhu, N. Y.; Wong, W.-T.; Yam, V. W.W. Diarylethene-Containing Cyclometalated Platinum (II) Complexes: Tunable Photochromism via Metal Coordination and Rational Ligand Design. J. Am. Chem. Soc. 2011, 133, 12690–12705. (35) Perissinotto, S.; Garbugli, M.; Fazzi, D.; Bertarelli, C.; Carvelli, M.; Srimath Kandada, A. R.; Yue, Z. N.; Wong, K. S.; Lanzani, G. Optical Modulation of Amplified Emission in a Polyfluorene-Diarylethene Blend. ChemPhysChem 2011, 12, 3619– 3623. (36) Perrier, A.; Maurel, F.; Jacquemin, D. Interplay between Electronic and Steric Effects in Multiphotochromic Diarylethenes. J. Phys. Chem. C 2011, 115, 9193–9203. (37) Perrier, A.; Maurel, F.; Jacquemin, D. Diarylethene-Dihydroazulene Multimode Photochrome: A Theoretical Spectroscopic Investigation. Phys. Chem. Chem. Phys. 2011, 13, 13791–13799. (38) Yokojima, S.; Kobayashi, T.; Shinoda, K.; Matsuda, K.; Higashiguchi, K.; Naka-

24 ACS Paragon Plus Environment

Page 24 of 32

Page 25 of 32

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

mura, S. π -Conjugation of Two Nitronyl Nitroxides-Attached Diarylethenes. J. Phys. Chem. B 2011, 115, 5685–5692. (39) Liu, S. J.; Sun, S.; Wang, C. M.; Zhao, Q.; Sun, H. B.; Li, F. Y.; Fan, Q. L.; Huang, W. DFT/TDDFT Investigation of the Modulation of Photochromic Properties in an Organoboron-Based Diarylethene by Fluoride Ions. ChemPhysChem 2011, 12, 313–321. (40) Nishi, H.; Asahi, T.; Kobatake, S. Plasmonic Enhancement of Gold Nanoparticles on Photocycloreversion Reaction of Diarylethene Derivatives Depending on Particle Size, Distance from the Particle Surface, and Irradiation Wavelength. Phys. Chem. Chem. Phys. 2012, 14, 4898–4905. (41) Matis, J. R.; Schönborn, J. B.; Saalfrank, P. A Multi-Reference Study of the Byproduct Formation for a Ring-Closed Dithienylethene Photoswitch. Phys. Chem. Chem. Phys. 2015, 17, 14088–14095. (42) Takeshita, M.; Ogawa, M.; Miyata, K.; Yamato, T. Synthesis and Photochromic Reaction of 1,2-Diphenylperfluorocyclopentenes. J. Phys. Org. Chem. 2003, 16, 148–151. (43) Takeshita, M.; Jin-nouchi, H.; Ishikawa, J. Syntheses and Photochromic Properties of Metacyclothiophenophan-1-enes. Chem. Lett. 2009, 38, 982–983. (44) Takeshita, M.; Inoue, M.; Hisasue, H.; Maekawa, S.; Nakamura, T. Synthesis and Photochromism of 1,2-Dicyano[2.n]metacyclophan-1-enes. J. Phys. Org. Chem. 2007, 20, 830–837. (45) Takeshita, M.;

Yamato, T. Synthesis and Photochromic Properties of 1,2-

Dicyano[2.n]metacyclophan-1-enes. Tetrahedron Lett. 2001, 42, 4345–4347. (46) Vosko, S. H.; Wilk, L.; Nusair, M. Accurate Spin-Dependent Electron Liquid Correla-

25 ACS Paragon Plus Environment

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

tion Energies for Local Spin Density Calculations: A Critical Analysis. Can. J. Phys. 1980, 58, 1200–1211. (47) Lee, C.; Yang, W. T.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. (48) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098–3100. (49) Yanai, T.; Tew, D. P.; Handy, N. C. A New Hybrid Exchange-Correlation Functional Using the Coulomb-Attenuating Method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51–57. (50) Chai, J.-D.; Head-Gordon, M. Systematic Optimization of Long-Range Corrected Hybrid Density Functionals. J. Chem. Phys. 2008, 128, 084106. (51) Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215– 241. (52) Christiansen, O.; Koch, H.; Jørgensen, P. The Second-Order Approximate Coupled Cluster Singles and Doubles Model CC2. Chem. Phys. Lett. 1995, 243, 409–418. (53) Hättig, C.; Köhn, A. Transition Moments and Excited-State First-Order Properties in the Coupled-Cluster Model CC2 Using the Resolution-of-the-Identity Approximation. J. Chem. Phys. 2002, 117, 6939–6951. (54) Hättig, C. Geometry Optimizations with the Coupled-Cluster Model CC2 Using the Resolution-of-the-Identity Approximation. J. Chem. Phys. 2003, 118, 7751–7761.

26 ACS Paragon Plus Environment

Page 26 of 32

Page 27 of 32

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

(55) Andersson, K.; Malmqvist, P.-Å.; Roos, B. O. Second-Order Perturbation Theory with a Complete Active Space Self-Consistent Field Reference Function. J. Chem. Phys. 1992, 96, 1218–1226. (56) Grimme, S.; Waletzke, M. A Combination of Kohn-Sham Density Functional Theory and Multi-Reference Configuration Interaction Methods. J. Chem. Phys. 1999, 111, 5645–5655. (57) Becke, A. D. A New Mixing of Hartree-Fock and Local Density-Functional Theories. J. Chem. Phys. 1993, 98, 1372–1377. (58) Boström, J.; Delcey, M. G.; Aquilante, F.; Serrano-Andrés, L.; Pedersen, T. B.; Lindh, R. Calibration of Cholesky Auxiliary Basis Sets for Multiconfigurational Perturbation Theory Calculations of Excitation Energies. J. Chem. Theory Comput. 2010, 6, 747–754. (59) Forsberg, N.; Malmqvist, P.-Å. Multiconfiguration Perturbation Theory with Imaginary Level Shift. Chem. Phys. Lett. 1997, 274, 196–204. (60) Dunning Jr., T. H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007– 1023. (61) Aquilante, F.; De Vico, L.; Ferré, N.; Ghigo, G.; Malmqvist, P.-Å.; Neogrády, P.; Pedˇ ersen, T. B.; Pitonák, M.; Reiher, M.; Roos, B. O. et al. MOLCAS 7: The Next Generation. J. Comput. Chem. 2010, 31, 224–247. (62) Schäfer, A.; Horn, H.; Ahlrichs, R. Fully Optimized Contracted Gaussian Basis Sets for Atoms Li to Kr. J. Chem. Phys. 1992, 97, 2571–2577. (63) TURBOMOLE V6.3 2011, a Development of University of Karlsruhe and

27 ACS Paragon Plus Environment

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

Forschungszentrum Karlsruhe GmbH, 1989-2007, TURBOMOLE GmbH, since 2007; available from http://www.turbomole.com. (64) Schäfer, A.; Huber, C.; Ahlrichs, R. Fully Optimized Contracted Gaussian Basis Sets of Triple Zeta Valence Quality for Atoms Li to Kr. J. Chem. Phys. 1994, 100, 5829– 5835. (65) 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 09, Revision B.01. Gaussian, Inc., Wallingford CT, 2010. (66) TURBOMOLE V5.7.1 2004, a Development of University of Karlsruhe and Forschungszentrum Karlsruhe GmbH, 1989-2007, TURBOMOLE GmbH, since 2007; available from http://www.turbomole.com. (67) Weber, W. Ph.D. thesis, University of Zürich, 1996. (68) Weber, W.; Thiel, W. Orthogonalization Corrections for Semiempirical Methods. Theor. Chem. Acc. 2000, 103, 495–506. (69) Koslowski, A.; Beck, M. E.; Thiel, W. Implementation of a General Multireference Configuration Interaction Procedure with Analytic Gradients in a Semiempirical Context Using the Graphical Unitary Group Approach. J. Comput. Chem. 2003, 24, 714– 726. (70) Thiel, W. MNDO99 Program, Version 6.1. Max-Planck-Institut für Kohlenforschung: Mülheim, Germany, 2007. (71) Tuna, D.; Lu, Y.; Koslowski, A.; Thiel, W. Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks of Electronically Excited States. J. Chem. Theory Comput. 2016, 12, 4400–4422.

28 ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

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

(72) Dral, P. O.; Wu, X.; Spörkel, L.; Koslowski, A.; Weber, W.; Steiger, R.; Scholten, M.; Thiel, W. Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Theory, Implementation, and Parameters. J. Chem. Theory Comput. 2016, 12, 1082–1096. (73) Tully, J. C.; Preston, R. K. Trajectory Surface Hopping Approach to Nonadiabatic Molecular Collisions: The Reaction of H+ with D2 . J. Chem. Phys. 1971, 55, 562– 572. (74) Hammes-Schiffer, S.; Tully, J. C. Proton Transfer in Solution: Molecular Dynamics with Quantum Transitions. J. Chem. Phys. 1994, 101, 4657–4667. (75) Fabiano, E.; Keal, T. W.; Thiel, W. Implementation of Surface Hopping Molecular Dynamics Using Semiempirical Methods. Chem. Phys. 2008, 349, 334–347. (76) Keal, T. W.; Wanko, M.; Thiel, W. Assessment of Semiempirical Methods for the Photoisomerisation of a Protonated Schiff Base. Theor. Chem. Acc. 2009, 123, 145– 156. (77) Granucci, G.; Persico, M.; Zoccante, A. Including Quantum Decoherence in Surface Hopping. J. Chem. Phys. 2010, 133, 134111. (78) Cui, G. L.; Lan, Z. G.; Thiel, W. Intramolecular Hydrogen Bonding Plays a Crucial Role in the Photophysics and Photochemistry of the GFP Chromophore. J. Am. Chem. Soc. 2012, 134, 1662–1672. (79) Cui, G. L.; Thiel, W. Nonadiabatic Dynamics of a Truncated Indigo Model. Phys. Chem. Chem. Phys. 2012, 14, 12378–12384. (80) Cui, G. L.; Thiel, W. Photoinduced Ultrafast Wolff Rearrangement: A Non-Adiabatic Dynamics Perspective. Angew. Chem. Int. Ed. 2013, 52, 433–436.

29 ACS Paragon Plus Environment

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

(81) Spörkel, L.; Cui, G. L.; Thiel, W. Photodynamics of Schiff Base Salicylideneaniline: Trajectory Surface-Hopping Simulations. J. Phys. Chem. A 2013, 117, 4574–4583. (82) Xia, S.-H.; Xie, B.-B.; Fang, Q.; Cui, G. L.; Thiel, W. Excited-State Intramolecular Proton Transfer to Carbon Atoms: Nonadiabatic Surface-Hopping Dynamics Simulations. Phys. Chem. Chem. Phys. 2015, 17, 9687–9697. (83) Xie, B.-B.; Xia, S.-H.; Chang, X.-P.; Cui, G. L. Photophysics of Auramine-O: Electronic Structure Calculations and Nonadiabatic Dynamics Simulations. Phys. Chem. Chem. Phys. 2016, 18, 403–413. (84) Liu, X.-Y.; Chang, X.-P.; Xia, S.-H.; Cui, G. L.; Thiel, W. Excited-State ProtonTransfer-Induced Trapping Enhances the Fluorescence Emission of a Locked GFP Chromophore. J. Chem. Theory Comput. 2016, 12, 753–764. (85) Xia, S.-H.; Cui, G. L.; Fang, W.-H.; Thiel, W. How Photoisomerization Drives Peptide Folding and Unfolding: Insights from QM/MM and MM Dynamics Simulations. Angew. Chem. Int. Ed. 2016, 55, 2107–2112. (86) Wang, Y.-T.; Liu, X.-Y.; Cui, G. L.; Fang, W.-H.; Thiel, W. Photoisomerization of Arylazopyrazole Photoswitches: Stereospecific Excited-State Relaxation. Angew. Chem. Int. Ed. 2016, 55, 14009–14013. (87) Kazaryan, A.; Lan, Z. G.; Schäfer, L. V.; Thiel, W.; Filatov, M. Surface Hopping Excited-State Dynamics Study of the Photoisomerization of a Light-Driven Fluorene Molecular Rotary Motor. J. Chem. Theory Comput. 2011, 7, 2189–2199. (88) Weingart, O.; Lan, Z. G.; Koslowski, A.; Thiel, W. Chiral Pathways and Periodic Decay in cis-Azobenzene Photodynamics. J. Phys. Chem. Lett. 2011, 2, 1506–1509. (89) Schönborn, J. B.; Koslowski, A.; Thiel, W.; Hartke, B. Photochemical Dynamics of E-iPr-Furylfulgide. Phys. Chem. Chem. Phys. 2012, 14, 12193–12201. 30 ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32

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

(90) Schönborn, J. B.; Hartke, B. Photochemical Dynamics of E-MethylfurylfulgideKinematic Effects on Photorelaxation Dynamics of Furylfulgides. Phys. Chem. Chem. Phys. 2014, 16, 2483–2490. (91) Shemesh, D.; Blair, S. L.; Nizkorodov, S. A.; Gerber, R. B. Photochemistry of Aldehyde Clusters: Cross-Molecular versus Unimolecular Reaction Dynamics. Phys. Chem. Chem. Phys. 2014, 16, 23861–23868. (92) Gerber, R. B.; Shemesh, D.; Varner, M. E.; Kalinowski, J.; Hirshberg, B. Ab Initio and Semi-Empirical Molecular Dynamics Simulations of Chemical Reactions in Isolated Molecules and in Clusters. Phys. Chem. Chem. Phys. 2014, 16, 9760–9775.

31 ACS Paragon Plus Environment

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

32 ACS Paragon Plus Environment

Page 32 of 32