Thermal

Subscriber access provided by READING UNIV

C: Plasmonics; Optical, Magnetic, and Hybrid Materials

Mechanisms of the Aggregation-Induced Emission and Photo/Thermal E/Z Isomerization of a Cyanostilbene Derivative: Theoretical Insights Norifumi Yamamoto J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b02147 • Publication Date (Web): 29 May 2018 Downloaded from http://pubs.acs.org on May 30, 2018

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 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 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.

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 22 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

Mechanisms of the Aggregation-Induced Emission and Photo/Thermal E/Z Isomerization of a Cyanostilbene Derivative: Theoretical Insights Norifumi Yamamoto* Department of Applied Chemistry, Faculty of Engineering, Chiba Institute of Technology, 2-17-1 Tsudanuma, Narashino, Chiba 275-0016, Japan

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 22

ABSTRACT

In this study, the aggregation-induced emission (AIE) of a cyanostilbene derivative, called CN-MBE, was investigated using molecular simulations. The E-form of CN-MBE has been known to exhibit AIE, which is non-emissive in dilute solutions but becomes highly emissive in aggregated states, while its Zform is non-emissive even in its crystalline form. In addition, upon UV light irradiation, the Z-form exhibits a Z-to-E isomerization, while the E-form keeps its conformation at room temperature and undergoes a nonradiative E-to-Z isomerization only at a high temperature. The results from the electronic structure calculations employed in this work showed that the potential energies of CN-MBE for the electronic ground (S0) and first excited (S1) states were degenerate at a twisted conformation around the ethylenic C=C π-bond, which led to the fluorescence quenching of the molecule. Molecular dynamics simulations and free-energy analyses revealed that the E-form molecules assembled closely, with the C=C bond rotation markedly restricted. This, in turn, prevented the fluorescence quenching via the S0/S1 conical intersection. In contrast, the Z-form molecules aggregated relatively sparsely, allowing for the nonradiative Z-to-E isomerization to proceed. All in all, the theoretical insights presented herein give a clear picture on the AIE and photo/thermal isomerization mechanisms of CNMBE.


2

Page 3 of 22 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


INTRODUCTION Fluorescent molecules with aggregation-induced emission (AIE), which are non-emissive in dilute solutions but become highly emissive in solid or aggregate phases, have come to attract a lot of attention as novel luminogenic materials for various optical devices,1,2 as well as chemical and biological sensing.3,4 The cyanostilbene derivative 2,3-bis(4′-methylbiphenyl-4-ly)acrylonitrile (CNMBE) is a prototypical AIE luminogen (AIE-gen),5,6 which has a small fluorescence quantum yield (Φ < 0.01) in dilute solutions that increases remarkably (Φ = 0.69) after nanoparticle formation.5 CN-MBE has two different geometric isomers called E- and Z-forms (Figure 1). While (E)-CN-MBE has been known to exhibit AIE, (Z)-CN-MBE is non-emissive even in crystalline phase.7 In addition, the Z-form is nonfluorescent but exhibits a Z-to-E isomerization upon UV light irradiation, while the E-form keeps its conformation upon UV light illumination at room temperature and undergoes a nonradiative E-to-Z isomerization only after UV light irradiation at a high temperature of about 500 K.7 Therefore, CNMBE has bright prospects as a novel functional material that exhibits AIE and photo/thermal E/Z isomerization beyond conventional molecular crystals and/or soft matters. The AIE of CN-MBE has been attributed to the combined effects of aggregation-induced planarization and J-aggregate formation in the solid phase.5,6,8,9 However, no definite view on the nature of AIE in CN-MBE has been given hitherto.


3


Page 4 of 22

Figure 1. Chemical structures of (E)- and (Z)-CN-MBE.

Apart from CN-MBE, tetraphenylethylene (TPE) is another archetypal AIE-gen.10 Because of its simple structure and easy synthesis, the AIE mechanism of TPE has been more widely investigated, both experimentally and theoretically, than that of CN-MBE.11-16 The main cause of AIE in TPE and its derivatives has been attributed to the restriction of intramolecular rotation (RIR) of their phenyl rings that are connected to the central ethylenic moiety via C−C σ-bonds.11,12 The RIR hypothesis suggests that in dilute solutions TPE undergoes fluorescence quenching through rotations of the phenyl rings by friction with solvent molecules. Upon aggregate formation, these rotations become restricted from the surrounding molecules, and a strong light emission is induced. Recently, however, it has been recognized that the ethylenic C=C π-bond rotation in TPE plays a crucial role in fluorescence quenching.14-16 Several computational studies based on electronic structure calculations revealed that the ground- and excited-states of TPE have the same energy with a twisted geometry about the ethylenic C=C bond site, and intersecting each other to form conical intersections (CIs).14,15 Even though this demonstrates that the fluorescence quenching in TPE can be caused by a rapid internal


4

Page 5 of 22 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


conversion via the CIs in a similar manner as in stilbene,17 the AIE mechanism is still a controversy.12,13 Considering that CN-MBE is also a substituted stilbene like TPE, the key to understanding the AIE mechanism of CN-MBE would be to examine whether and how the fluorescence quenching occurs, and if it involves an ethylenic C=C bond rotation. However, no study has elucidated this problem so far. Moreover, even though the steric hindrance from the surrounding molecules is essential to the AIE process, computational studies on the AIE-gens have been limited to single-point potential energy (PE) calculations of the molecules that ignore the thermodynamic influence from the environment. In this study, the mechanisms of AIE and photo/thermal E/Z isomerization of CN-MBE were investigated based on the free-energy (FE) profiles of the molecule in condensed phases as functions of the ethylenic C=C bond rotation.

COMPUTATIONAL DETAILS Electronic Structure Calculations. In order to elucidate the characteristics of the PE profiles for the photochemical processes of CN-MBE, the minimum energy paths (MEPs) for the E/Z isomerization reaction of CN-MBE in an isolated phase were investigated using electronic structure calculations. Constrained geometry optimizations were performed by fixing the torsional angle of the ethylenic C=C bond in CN-MBE, φ, (Figure 1) over a range of 0° to 180° for the electronic ground (S0) and first excited (S1) states. The spin-flip approach18 within the time-dependent density functional theory (SFTDDFT) method19 was used to compute the PEs and analytical gradients for the S0 and S1 states of the molecule. Since the SF-TDDFT method is known to suffer from spin-contamination, the lowest two electronic states with the values of 〈S S 〉 smaller than 0.5 were classified as the S0 and S1 states. The BHHLYP hybrid functional (50% Hartree-Fock plus 50% Becke exchange,20 the Lee-Yang-Parr correlation21) and 6-31G(d) basis set were employed in the SF-TDDFT calculations. In previous studies, the SF-TDDFT method was utilized to examine the photochemical reactions of ethylene22 and ACS Paragon Plus Environment

5


Page 6 of 22

stilbene.23 The results from these investigations were comparable to those obtained using the multiconfigurational ab initio molecular orbital methods, e.g., at the complete active space selfconsistent field (CASSCF)24 and second-order multireference perturbation theory (CASPT2)25 levels. Therefore, taking into account that CN-MBE is considered to be a substituted stilbene, the SF-TDDFT method would be a promising approach in this study. The ground- and excited-state PE surfaces of ethylene and stilbene are known to be nearly degenerate with twisted geometries about the central C=C bond site, and intersecting with each other to form CIs.17 Therefore, in this work, we investigated the minimum energy point of the conical intersection (MECI) between the S0 and S1 states of CN-MBE at the SF-TDDFT level using the penalty-constrained optimization method.26,27 All electronic structure calculations of CN-MBE presented herein were performed using the GAMESS program.28 Molecular Dynamics Simulations. The FE profile of the E/Z isomerization of CN-MBE in condensed phase was computed using molecular dynamics (MD) simulations based on an empirical force field representation. All MD simulations of CN-MBE presented herein were carried out using the GROMACS 2016.3 program with a GPU acceleration.29 Following the energy minimization and equilibration processes, production runs were performed under the isothermal-isobaric (NpT) ensemble. The temperature was maintained either at room temperature (300 K) or higher (500 K) using a velocityrescaling thermostat.30 The pressure was kept constant at 1 bar using a Parrinello−Rahman barostat.31 The time step used in all MD simulations was 2 fs and the list of nonbonded interaction pairs was updated every 10 steps with a cutoff radius of 1 nm. The electrostatic interactions were computed using the particle mesh Ewald algorithm32 with a cutoff radius of 1 nm, in which the reciprocal part was performed by fast Fourier transformations. All bond lengths were constrained using the linear constraint solver (LINCS) algorithm.33 Fitting Force Field Parameters. As demonstrated by previous studies on TPE-based AIE-gens,14-16 the rotation around the ethylenic C=C bond can play an important role in the photochemical processes of CN-MBE. Therefore, the empirical force fields used in the MD simulations were specifically ACS Paragon Plus Environment

6

Page 7 of 22 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


designed to comply with the SF-TDDFT calculations on the ethylenic C=C bond rotation by varying the parameters of the energy functions. In this study, an analytic linear least-square fitting method34 was employed for the optimization of the force field parameters. These parameters for the rotation around the ethylenic C=C bond of CN-MBE, φ, (Figure 1) were derived separately for the S0 and S1 states. Moreover, since the four C−C σ-bonds that connect the four phenyl rings in CN-MBE are rotatable, a set of corresponding dihedral terms were reparametrized based on the electronic structure calculations for model systems at the MP2/6-31G(d)//HF/6-31G(d) level. Computational details on the optimization of the force field parameters are included in the Supporting Information. In addition, other force field parameters, apart from the aforementioned dihedral terms, i.e., bond, angle, and nonbonding parameters, were determined based on the OPLS-UA force field.35 Free-Energy Calculations. The structural transition between the E- and Z-forms of CN-MBE was examined and the corresponding FE changes were computed with a well-tempered (WT) variant36 of the metadynamics (MetaD) algorithm37,38 by biasing the MD simulations using the PLUMED 2.4.0 plugin.39 While conventional MD simulations tend to do oversampling of the same conformations, MetaD avoids the oversampling problem by constructing a history-dependent bias potential, VG(s, t), as a sum of Gaussians centered along the MD trajectory followed by collective variables (CVs), s, up to time t.37,38 The dihedral angle of the ethylenic C=C bond rotation (φ) was used as a CV in the process. A bias Gaussian with a width of 10° was deposited every 250 timesteps of the simulations and the height of the bias was set to 1 kJ mol−1 with a bias factor of 40. The WT-MetaD simulations were parallelized using the multiple walker method,40 in which four 100-ns simulations, starting from different initial coordinates, were run in parallel. The FE profiles along the CVs, F(s), were obtained from the accumulated bias potential, VG, according to the equation F(s) = −VG(s, t → ∞) + C, where C is an irrelevant additive constant. Statistical errors in the FE profile determination were estimated through a block-averaging analysis.38


7


Page 8 of 22

Predicting Aggregated Structures. The aggregated structures of CN-MBE were predicted using the CONFLEX program.41 Although the crystal structures of other cyanostilbene derivatives similar to CNMBE have been determined,42,43 those of the (E)- and (Z)-CN-MBE molecules have not yet been clarified. Therefore, in this study, we explored possible crystal structures of these two molecules by randomly generating trial structures and subsequently optimizing the lattice parameters. In total, 3,000 trial molecular orientations were generated by rotating the rigid molecule in the unit cells. The spacegroup symmetries were constrained in 1 and P21/c during the lattice optimizations because the previously reported crystal structures of CN-MBE-like molecules were contained in one of these space groups.42,43 The crystal energies of all trial crystal structures were estimated based on the Merck Molecular Force Filed 94 (MMFF94)44 and the most stable ones were utilized in the MD simulations.

RESULTS AND DISCUSSION PE Profile of CN-MBE. Figure 2A shows the PE profiles along the S0- and S1-MEPs of CN-MBE in isolated phase. The change in PE of the S0 state along the S0-MEP (denoted as S0//S0-MEP) demonstrated that the isomerization between the E- (φ = 180°) and Z-forms (φ = 0°) of the molecule proceeded very slowly in the S0 state due to the large energy barrier (34 kcal mol−1) at φ = 90°. In contrast, the S1 PE change along the S1-MEP (denoted as S1//S1-MEP) indicated that the torsional motion about the ethylenic C=C bond reduced the energy from the vertically excited Franck-Condon (FC) points of the E- (φ = 180°) or Z-form (φ = 0°) to the minimum energy point (φ = 90°) with no energy barrier. Figure 2B shows a geometry corresponding to the S1-MEP minimum at φ = 90°. As the ethylenic C=C bond rotated along the S1-MEP, the PE gap between the S0 and S1 states decreased, with values within 1.23 eV in the region of φ = 60°−130° and a minimum energy difference of 0.84 eV at

φ = 120°. This result demonstrated that the S1-MEP passed the region where the S0 and S1 PE surfaces approach each other but do not intersect. Moreover, the S0/S1 MECI, in which the PEs of the S0 and S1


8

Page 9 of 22 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


states were almost identical (within 1 kcal mol−1), was discovered. Figure 2C shows a geometry corresponding to the S0/S1 MECI of CN-MBE, which was characterized not only by the twisting (φ = 75°) but also the pyramidalization of the ethylenic C=C bond, thereby indicating that the S0/S1 MECI did not lie on the S1-MEP only with the ethylenic C=C bond rotation. This phenomenon is similar to the one occurring in the ethylene and stilbene molecules.17 All in all, these results indicated that the rotation around the ethylenic C=C bond of CN-MBE was one of the principal coordinates that led to the S0/S1 CIs, where the molecule dissipates the excitation energy nonradiatively via a rapid S0/S1 internal conversion to cause fluorescence quenching.

Figure 2. (A) Potential energy profiles of the S0 and S1 states along the S0- and S1-MEPs with a change in the torsional angle (φ) of the ethylenic C=C bond of CN-MBE in an isolated phase. Geometries corresponding to the S1-MEP minimum at φ = 90° (B) and the S0/S1 MECI (C).

Notably, the 〈S S 〉 values of the lowest two states of geometries on the S0-MEP were almost constant, with average 〈S S 〉 values of 0.17 and 0.38, and were thus successfully classified as the S0 and S1 states, respectively. In contrast, along the S1-MEP, the 〈S S 〉 values for the lowest state, which was supposed to be the S0 state, increased with the twisting about the ethylenic C=C bond. As a result, the value of 〈S S 〉 was 1.35 at the twisted geometry (φ = 90°) and around 0.2 at the E- (φ = 180°) and Z-form (φ = 0°) ACS Paragon Plus Environment

9


Page 10 of 22

structures. This indicates that the lowest singlet state was strongly mixed with the triplet state at the twisted geometry of the S1-MEP, which was due to the fact that, at the 90°-twisted geometry, the ethylenic C=C π-bond is broken to produce two degenerate nonbonding p-orbitals containing two electrons, called diradicaloid.45

FE Profile of CN-MBE in THF Solution. Figure 3A shows the FE changes accompanying the rotation around the ethylenic C=C bond of CN-MBE in THF solution, with the overall errors in the FE profiles estimated to be smaller than 1 kJ mol−1. These were obtained via the WT-MetaD simulations of a system in which one CN-MBE molecule was solvated in 434 THF molecules enclosed in a cubic box of 62.8 nm3 under periodic boundary conditions (Figure 3B). According to the obtained results, no energy barrier existed in the FE profile of the S1 state from the vertically excited FC geometries of the E- and Z-forms (φ = 180° and 0°) to the twisted geometry (φ = 90°) (Figure 3A). Importantly, neither the S0/S1 CI nor the MECI of CN-MBE was determined from the FE profiles in THF solution and aggregated state. Instead, the MECI was obtained in an isolated phase, which was characterized by the twisting and pyramidalization of the ethylenic C=C bond. Presumably, the intermolecular interaction with the surrounding solvent or aggregated molecules had an important contribution to the entropy term of free-energy through steric hindrance but did not significantly affect the electronic structures. In other words, both in THF solution and as aggregates, the essential characteristics of the PE surfaces observed in isolated phase might be maintained. Taking this into account, we assumed that the twisting about the ethylenic C=C bond could be one of the principal coordinates that leads to the S0/S1 CIs both in THF solution and aggregated state. Hence, the FE profile shown in Figure 3A indicates that the corresponding structural change, i.e., the torsional motion of the ethylenic C=C bond site, occurred spontaneously even in THF solution, which, in turn, induced the nonradiative de-excitation transition from the S1 to S0 states in the vicinity of the 90°-twisted geometry. Therefore, it was concluded that


10

Page 11 of 22 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


since the S0/S1 CIs could be reached efficiently and facilitate fluorescence quenching after CN-MBE photoexcitation, the molecule would show no emission when dispersed in dilute solutions.

Figure 3. (A) Free-energy profile of the changes in torsional angle (φ) of the ethylenic C=C bond site of CN-MBE in THF solution. The point at which the torsion angle becomes φ = 90° on the S1-FE profile is indicated by a filled circle, while local minima on the S0-FE profile and FC points on the S1FE profile at the E- (φ = 180°) and Z-form (φ = 0°) geometries are indicated by crosses. (B) CN-MBE in THF obtained from MD simulations.

Structures of the CN-MBE Aggregates. Figure 4 shows the crystal structures of the E- and Z-forms of CN-MBE predicted in this study. The most stable structures of the E- and Z-form aggregates were found to be monoclinic systems with a P21/c space group, containing four molecules in the unit cell. Structural details of the CN-MBE aggregates are included in the Supporting Information. MD simulations of the aggregated states of CN-MBE were performed for the systems consisting of 196 Eand 192 Z-form molecules in monoclinic 1×7×7 and 6×2×4 supercells with sizes of 108.3 and 109.1 nm3, respectively, under periodic boundary conditions. Figures 5A and 6A present snapshot pictures of the (E)- and (Z)-CN-MBE aggregates, respectively, which were obtained from these simulations. The fractional free volumes (FFVs)46 were analyzed from the results of the MD simulations performed at ACS Paragon Plus Environment

11


Page 12 of 22

300 K. The FFVs of the E- and Z-forms of CN-MBE were found to be 22.1% and 24.2%, respectively, which indicated that the E-forms were more densely packed than the Z-forms in the aggregated phase.

Figure 4. Crystal structures of the E- (A) and Z-forms (B) of CN-MBE predicted in this study.

FE Profile of the (E)-CN-MBE Aggregates at 300 K. Figure 5B shows the FE profile of (E)-CNMBE in its aggregate phase obtained from the WT-metaD simulations performed at 300 K. As can be observed, the resulting FE profile is significantly different from that of CN-MBE in THF (Figure 3A). For the E-form aggregate, the conformational change leading to the twisted geometry (φ = 90°) from the vertically excited FC geometry (φ = 180°) on the S1 energy profile was revealed to be energetically unfavorable and experienced an increase in the FE of 57 kJ mol−1. Therefore, it was demonstrated that when (E)-CN-MBE was photoexcited in the aggregate phase, the rotation around the ethylenic C=C bond required to reach the S0/S1 CIs was energetically demanding because of the steric hindrance from the firmly stacked surrounding molecules. As a result, it was concluded that the sizeable steric hindrance that prevented the nonradiative transition via the S0/S1 CIs at the twisted conformations was responsible for the strong emission of (E)-CN-MBE in the aggregated phase. ACS Paragon Plus Environment

12

Page 13 of 22 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 5. (A) The (E)-CN-MBE aggregated structure obtained from MD simulations. Free-energy changes accompanying the rotation around the ethylenic C=C bond site of (E)-CN-MBE in the aggregated structure at 300 K (B) and 500 K (C).

FE Profile of the (Z)-CN-MBE Aggregates at 300 K. Figure 6B shows the FE profile of the Z-form of CN-MBE in its aggregate phase obtained from the WT-metaD simulations at 300 K. In contrast to the (E)-CN-MBE profile (Figure 5B), the FE profile of the Z-form aggregates resembled that of CNMBE in THF (Figure 3A). Notably, no energy barrier was detected in the FE profile of the S1 state from the vertically excited FC geometry of the Z-form (φ = 0°) to the twisted geometry (φ = 90°). This indicated that after the photoexcitation, (Z)-CN-MBE could relax from the FC point to the S0/S1 CIs on the S1 energy surface in a barrierless manner, which was similar to that of CN-MBE in THF but different from that of the (E)-CN-MBE aggregate. Therefore, it was concluded that the Z-form of CNMBE did not exhibit AIE due to the weakly packed structure that allowed for the rotation around the ethylenic C=C bond to reach the S0/S1 CIs. In this process, the molecules that reached the S0/S1 CIs separated in two paths - one that returned to the Z-form and one that changed to the E-from. Thus, it was evident that the Z-to-E isomerization reaction of CN-MBE occurred readily even at room temperature after photoexcitation.


13


Page 14 of 22

Figure 6. (A) The (Z)-CN-MBE aggregated structure obtained from MD simulations. Free-energy changes accompanying the rotation around the ethylenic C=C bond site of (Z)-CN-MBE in the aggregated structure at 300 K (B) and 500 K (C).

FE Profiles of the CN-MBE Aggregates at 500 K. Finally, MD simulations of the (E)- and (Z)-CNMBE aggregates were performed at a higher temperature (500 K). As a natural consequence, the aggregated structures became loose due to the dynamic motion of the molecules induced by the heat. The resulting FFVs at 500 K were determined to be 33.7% and 29.4% for the E- and Z-form aggregates, respectively. As can be observed from Figure 6B and C, the FE profile of the Z-form aggregate obtained from the WT-metaD simulations at 500 K was almost identical to that at 300 K. In contrast, the FE profile of the E-form aggregate at 500 K appeared very different from that at 300 K (Figure 5C and B, respectively). These results indicated that the photoexcited molecule could reach the S0/S1 CIs and cause fluorescence quenching due to the marked reduction in the energy barrier of the ethylenic C=C bond torsion. Therefore, it was concluded that at the higher temperature (500 K) the rotation around the ethylenic C=C bond was allowed even in the E-form aggregate, which led to a nonradiative E-to-Z isomerization.


14

Page 15 of 22 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


CONCLUSIONS In this study, we investigated the mechanisms of AIE and photo/thermal E/Z isomerization of CNMBE using molecular simulations. The PE profiles showed that the S0 and S1 states were degenerate at the twisted conformation around the ethylenic C=C bond, which led to the fluorescence quenching of the molecule via the S0/S1 CIs. The FE profiles showed that the rotation around the ethylenic C=C bond of the (E)-CN-MBE aggregate was markedly restricted at room temperature, thus preventing the quenching through the S0/S1 CIs at its twisted conformation and leading to strong fluorescence. In contrast, the ethylenic C=C bond in the Z-form aggregate could easily rotate to reach the S0/S1 CIs, which led to the barrierless nonradiative transition and subsequent Z-to-E isomerization reaction. However, the rotation of the ethylenic C=C bond in the E-form aggregate was allowed when exposed to high temperatures, thereby creating a nonradiative E-to-Z isomerization. All in all, the method employed in this study combines electronic structure calculations and MD simulations to compute the PE and FE profiles of the photochemical processes in CN-MBE, which can not only be used for analysis of various other AIE-gens but also for rational design of new photochemical materials.

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website. Computational details about the optimization of the force field parameters, structural details of (E)and (Z)-CN-MBE aggregates (PDF)

AUTHOR INFORMATION Corresponding Author ACS Paragon Plus Environment

15


Page 16 of 22

*N.Y.: E-mail: [email protected] ORCID Norifumi Yamamoto: 0000-0001-9900-8288 Notes The authors declare no competing financial interests.

ACKNOWLEDGMENT This work was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI (Grantin-Aid for Scientific Research (C) 17K07315), and the JGC-S Scholarship Foundation. The molecular simulations were carried out on the TSUBAME 3.0 supercomputer in Tokyo Institute of Technology, and the Research Center for Computational Science in the National Institutes of Natural Sciences (NINS).


16

Page 17 of 22 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


REFERENCES (1)

Hong, Y.; Lam, J. W. Y.; Tang, B. Z. Aggregation-Induced Emission. Chem. Soc. Rev. 2011, 40, 5361–5388.

(2)

Mei, J.; Leung, N. L. C.; Kwok, R. T. K.; Lam, J. W. Y.; Tang, B. Z. Aggregation-Induced Emission: Together We Shine, United We Soar! Chem. Rev. 2015, 115, 11718–11940.

(3)

Ding, D.; Li, K.; Liu, B.; Tang, B. Z. Bioprobes Based on AIE Fluorogens. Acc. Chem. Res. 2013, 46, 2441–2453.

(4)

Gao, M.; Tang, B. Z. Fluorescent Sensors Based on Aggregation-Induced Emission: Recent Advances and Perspectives. ACS Sens. 2017, 2, 1382–1399.

(5)

An, B.-K.; Kwon, S.-K.; Jung, S.-D.; Park, S. Y. Enhanced Emission and Its Switching in Fluorescent Organic Nanoparticles. J. Am. Chem. Soc. 2002, 124, 14410–14415.

(6)

An, B.-K.; Gierschner, J.; Park, S. Y. π-Conjugated Cyanostilbene Derivatives: A Unique SelfAssembly Motif for Molecular Nanostructures with Enhanced Emission and Transport. Acc. Chem. Res. 2012, 45, 544–554.

(7)

Chung, J. W.; Yoon, S.-J.; An, B.-K.; Park, S. Y. High-Contrast On/Off Fluorescence Switching via Reversible E–Z Isomerization of Diphenylstilbene Containing the αCyanostilbenic Moiety. J. Phys. Chem. C 2013, 117, 11285–11291.

(8)

Ito, F.; Fujimori, J. Fluorescence Visualization of the Molecular Assembly Processes During Solvent Evaporation via Aggregation-Induced Emission in a Cyanostilbene Derivative. CrystEngComm 2014, 16, 9779–9782.

(9)

Ito, F.; Fujimori, J.; Oka, N.; Sliwa, M.; Ruckebusch, C.; Ito, S.; Miyasaka, H. AIE Phenomena of a Cyanostilbene Derivative as a Probe of Molecular Assembly Processes. Faraday Discuss. 2017, 196, 231–243.

(10)

Zhao, Z.; Lam, J. W. Y.; Tang, B. Z. Tetraphenylethene: A Versatile AIE Building Block for


17


Page 18 of 22

the Construction of Efficient Luminescent Materials for Organic Light-Emitting Diodes. J. Mater. Chem. 2012, 22, 23726–23740. (11)

Shi, J.; Chang, N.; Li, C.; Mei, J.; Deng, C.; Luo, X.; Liu, Z.; Bo, Z.; Dong, Y. Q.; Tang, B. Z. Locking the Phenyl Rings of Tetraphenylethene Step by Step: Understanding the Mechanism of Aggregation-Induced Emission. Chem. Commun. 2012, 48, 10675–10677.

(12)

Tseng, N.-W.; Liu, J.; Ng, J. C. Y.; Lam, J. W. Y.; Sung, H. H. Y.; Williams, I. D.; Tang, B. Z. Deciphering Mechanism of Aggregation-Induced Emission (AIE): Is E–Z Isomerisation Involved in an AIE Process? Chem. Sci. 2012, 3, 493–497.

(13)

Yang, Z.; Qin, W.; Leung, N. L. C.; Arseneault, M.; Lam, J. W. Y.; Liang, G.; Sung, H. H. Y.; Williams, I. D.; Tang, B. Z. A Mechanistic Study of AIE Processes of TPE Luminogens: Intramolecular Rotation vs. Configurational Isomerization. J. Mater. Chem. C 2015, 4, 99–107.

(14)

Gao, Y.-J.; Chang, X.-P.; Liu, X.-Y.; Li, Q.-S.; Cui, G.; Thiel, W. Excited-State Decay Paths in Tetraphenylethene Derivatives. J. Phys. Chem. A 2017, 121, 2572–2579.

(15)

Kokado, K.; Machida, T.; Iwasa, T.; Taketsugu, T.; Sada, K. Twist of C=C Bond Plays a Crucial Role in the Quenching of AIE-Active Tetraphenylethene Derivatives in Solution. J. Phys. Chem. C 2017, 122, 245–251.

(16)

Xiong, J.-B.; Yuan, Y.-X.; Wang, L.; Sun, J.-P.; Qiao, W.-G.; Zhang, H.-C.; Duan, M.; Han, H.; Zhang, S.; Zheng, Y.-S. Evidence for Aggregation-Induced Emission from Free Rotation Restriction of Double Bond at Excited State. Org. Lett. 2018, 20, 373–376.

(17)

Levine, B. G.; Martı́nez, T. J. Isomerization Through Conical Intersections. Annu. Rev. Phys. Chem. 2007, 58, 613–634.

(18)

Krylov, A. I. Size-Consistent Wave Functions for Bond-Breaking: The Equation-of-Motion Spin-Flip Model. Chem. Phys. Lett. 2001, 338, 375–384.

(19)

Shao, Y.; Head-Gordon, M.; Krylov, A. I. The Spin-Flip Approach Within Time-Dependent Density Functional Theory: Theory and Applications to Diradicals. J. Chem. Phys. 2003, 118, ACS Paragon Plus Environment

18

Page 19 of 22 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


4807–4818. (20)

Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098–3100.

(21)

Lee, C.; Yang, W.; 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.

(22)

Minezawa, N.; Gordon, M. S. Optimizing Conical Intersections by Spin-Flip Density Functional Theory: Application to Ethylene. J. Phys. Chem. A 2009, 113, 12749–12753.

(23)

Minezawa, N.; Gordon, M. S. Photoisomerization of Stilbene: A Spin-Flip Density Functional Theory Approach. J. Phys. Chem. A 2011, 115, 7901–7911.

(24)

Roos, B. O. The Complete Active Space Self‐Consistent Field Method and Its Applications in Electronic Structure Calculations. Adv. Chem. Phys. 1987, 69, 399–445.

(25)

Roos, B. O. Theoretical Studies of Electronically Excited States of Molecular Systems Using Multiconfigurational Perturbation Theory. Acc. Chem. Res. 1999, 32, 137–144.

(26)

Levine, B. G.; Ko, C.; Quenneville, J.; Martı́nez, T. J. Conical Intersections and Double Excitations in Time-Dependent Density Functional Theory. Mol. Phys. 2006, 104, 1039–1051.

(27)

Levine, B. G.; Coe, J. D.; Martı́nez, T. J. Optimizing Conical Intersections Without Derivative Coupling Vectors: Application to Multistate Multireference Second-Order Perturbation Theory (MS-CASPT2). J. Phys. Chem. B 2008, 112, 405–413.

(28)

Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.; Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S.; et al. General Atomic and Molecular Electronic Structure System. J. Comput. Chem. 1993, 14, 1347–1363.

(29)

Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS: High Performance Molecular Simulations Through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1-2, 19–25.

(30)

Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling Through Velocity Rescaling. J. ACS Paragon Plus Environment

19


Page 20 of 22

Chem. Phys. 2007, 126, 014101. (31)

Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52, 7182–7190.

(32)

Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N⋅log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089–10092.

(33)

Hess, B.; Bekker, H.; Berendsen, H. J. C.; Fraaije, J. G. E. M. LINCS: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463–1472.

(34)

Hopkins, C. W.; Roitberg, A. E. Fitting of Dihedral Terms in Classical Force Fields as an Analytic Linear Least-Squares Problem. J. Chem. Inf. Model. 2014, 54, 1978–1986.

(35)

Jorgensen, W. L.; Tirado-Rives, J. The OPLS Potential Functions for Proteins. Energy Minimizations for Crystals of Cyclic Peptides and Crambin. J. Am. Chem. Soc. 1988, 110, 1657–1666.

(36)

Barducci, A.; Bussi, G.; Parrinello, M. Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method. Phys. Rev. Lett. 2008, 100, 020603.

(37)

Laio, A.; Parrinello, M. Escaping Free-Energy Minima. PNAS 2002, 99, 12562–12566.

(38)

Laio, A.; Gervasio, F. L. Metadynamics: A Method to Simulate Rare Events and Reconstruct the Free Energy in Biophysics, Chemistry and Material Science. Rep. Prog. Phys. 2008, 71, 126601.

(39)

Tribello, G. A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G. PLUMED 2: New Feathers for an Old Bird. Comput. Phys. Comm. 2014, 185, 604–613.

(40)

Raiteri, P.; Laio, A.; Gervasio, F. L.; Micheletti, C.; Parrinello, M. Efficient Reconstruction of Complex Free Energy Landscapes by Multiple Walkers Metadynamics. J. Phys. Chem. B 2006, 110, 3533–3539.

(41)

Obata, S.; Goto, H. High-Speed Prediction of Crystal Structures for Organic Molecules. AIP Conf. Proc. 2015, 1649, 130–134. ACS Paragon Plus Environment

20

Page 21 of 22 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

(42)


Chung, J. W.; You, Y.; Huh, H. S.; An, B.-K.; Yoon, S.-J.; Kim, S. H.; Lee, S. W.; Park, S. Y. Shear- and UV-Induced Fluorescence Switching in Stilbenic π-Dimer Crystals Powered by Reversible [2 + 2] Cycloaddition. J. Am. Chem. Soc. 2009, 131, 8163–8172.

(43)

An, B.-K.; Gihm, S. H.; Chung, J. W.; Park, C. R.; Kwon, S.-K.; Park, S. Y. Color-Tuned Highly Fluorescent Organic Nanowires/Nanofabrics: Easy Massive Fabrication and Molecular Structural Origin. J. Am. Chem. Soc. 2009, 131, 3950–3957.

(44)

Halgren, T. A. Merck Molecular Force Field. I. Basis, Form, Scope, Parameterization, and Performance of MMFF94. J. Comput. Chem. 1996, 17, 490–519.

(45)

Turro, N. J.; Ramamurthy, V.; Scaiano, J. C. Principles of Molecular Photochemistry: An Introduction; University Science Books: Sausalito, California, 2009.

(46)

Lourenço, T. C.; Coelho, M. F. C.; Ramalho, T. C.; van der Spoel, D.; Costa, L. T. Insights on the Solubility of CO2 in 1-Ethyl-3-Methylimidazolium Bis(Trifluoromethylsulfonyl)Imide from the Microscopic Point of View. Environ. Sci. Technol. 2013, 47, 7421–7429.


21


Page 22 of 22

TOC Graphic


22

Thermal

Recommend Documents