Reaction Mechanism of Photodeamination ... - ACS Publications

Jun 4, 2018 - Excited-state intramolecular proton transfer (ESIPT) acts as the most significant .... the total spin density matrix P = 0, and eq 1 can...
7 downloads 0 Views 3MB Size
Subscriber access provided by Miami University Libraries

A: Spectroscopy, Molecular Structure, and Quantum Chemistry

Reaction Mechanism of Photodeamination Induced by ExcitedState Intramolecular Proton Transfer of Anthrols Molecule Yunfan Yang, Yong Ding, Yu Zhao, Feng-Cai Ma, and Yong-Qing Li J. Phys. Chem. A, Just Accepted Manuscript • Publication Date (Web): 04 Jun 2018 Downloaded from http://pubs.acs.org on June 4, 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 25 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

Reaction Mechanism of Photodeamination Induced by Excited-State Intramolecular Proton Transfer of Anthrols Molecule Yunfan Yang†, Yong Ding†, Yu Zhao,Fengcai Ma* and Yongqing Li* Department of Physics, Liaoning University, Shenyang 110036, P. R. China.

Abstract Photodeamination reaction of Anthrols molecule generating the high-activity quinone methides intermediate had been investigated (J. Org. Chem. 2017, 82, 6006-6021), while lacking careful explanation for its reaction mechanism. In our research, the mechanism of Anthrols molecule photodeamination induced by excited state intramolecular proton transfer was reported for the first time. Absorption and emission spectra calculated for the work presented here agreed well with experimental results. To propose a molecular-level explanation of photodeamination reaction, we calculated bond parameters, corresponding infrared vibrational frequencies, Mayer bond orders and visualized frontier molecular orbitals of different molecules, and the hydrogen bond strengthening behavior in excited states was confirmed, enhancing the excited state intramolecular proton transfer of Anthrols molecule. Moreover, we concluded that photoisomerization weakened the bond strength between nitrogen and carbon atoms, which promoted photodeamination reaction. Finally, for visually and quantificationally revealing the photodeamination mechanism, we have established the three-dimensional potential energy surfaces for deamination reaction in different electronic states and calculated the corresponding reaction Gibbs free energies, it can be confirmed that the photodeamination reaction of Anthrols molecule must be induced by proton transfer reaction in excited state.

† These authors contributed equally to this work. *E-mail: [email protected] and [email protected].

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

Introduction Excited state intramolecular proton transfer (ESIPT) acts as the most significant photochemical switching processes of excited-state hydrogen bonding dynamics.1-8 Hydrogen bond originating from electrostatic interaction offers driving force for proton transfer reactions, which depend on the changes of hydrogen-bond properties in different electronic states.9 Usually, the proton transfer reaction is an endothermic process in ground (S0) state, the keto form of molecule cannot be obtained from enol form, since keto form is thermodynamically unstable in S0 state. However, upon the photoexcitation process, the electron population will present distinct changes, its redistribution results in the stronger acidity in proton donor group and the stronger basicity in proton acceptor group, strengthening hydrogen bond interaction facilitates hydrogen proton to transfer to acceptor group.10 The exothermic photoisomerization process generates stable keto form in the first excited (S1) state. On the basis of Kasha’s rule, if a molecule emits fluorescence, in general it will be from the lowest excited state to ground state,11 while comparing with excitation energy, the emission of keto form is of a large Stokes’ shift with 8000-10000 cm-1,12 therefore the distinct dual fluorescent phenomenon originating enol and keto in S1 state will be observed, which is an important reason that the ESIPT reactions play a key role in lots of fields. For example, white organic light-emitting diodes (OLED)13,14, newer arenas of biotechnology15,16, fluorescence imaging17,18, fluorescence probing19-21, optoelectronic devices22, laser dyes23,24 and so on. Since the ESIPT reaction in salicylate derivatives was reported firstly in experiment by Weller,25 during the past five decades, the ESIPT reactions have been widely studied in theory and experiment by researchers.26-36 Recently, for example, Zhou et al. investigated the effect of different conditions and intermolecular hydrogen bond sites on ESIPT reaction for methyl salicylate in 2015.37 Zhao et al. proposed theoretically the new mechanism of the simultaneous excited state intramolecular double proton transfer in 2015.38 Li et al. summarized an advanced review of the fluorescence-based sensing mechanisms involving in ESIPT reactions in 2017.39 ACS Paragon Plus Environment

Page 2 of 25

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

Goswami et al. reported the photophysical phenomenon of near infra-red emission induced by ESIPT reaction for acetate in 2013.40 Wang et al. discussed the mechanism of ESIPT-induced photocyclization reactions in 2011.41 Škalamera et al.42 (in experiment) and our group43 (in theory) reported the mechanism of photodeamination induced by ESIPT reaction for p-cresol derivatives in 2015 and 2017, respectively. In 2017, Škalamera et al. have synthesized series of Anthrols derivative, and reported that the high-activity quinone methides (QMs) were observed by laser flash photolysis method in experiment. They referred to the deamination of Anthrols in CH3OH solvent might be induced by ESIPT reaction,44 while this conclusion lacks a carefully theoretical interpretation. To quantitatively and qualitatively illustrate the mechanism of photodeamination induced by ESIPT reaction, the theoretical calculation methods are carried out in our research work. The molecular-level explanations are presented by analyzing bond parameters and corresponding infrared (IR) vibrational frequencies, visualizing frontier molecular orbital (FMOs), calculating Mayer bond orders, constructing three-dimensional potential energy surfaces (3DPES) and calculating reactive Gibbs free energies.

Computational Method In order to obtain correct molecular constructions, the implicit solvation model based on density (SMD)45 is considered with density functional theory (DFT) and its Time-Dependent (TD) method. The TD-DFT functional method has been confirmed by contemporary researchers that it is efficient and accurate measure to investigate photochemical reactions.46-57 Dielectric constant (ε) has been set as 38.8 to simulate the acetonitrile solution environment in experiment. All geometric constructions are fully optimized by Gaussian 09 software58 carried out B3LYP-D3/6-31G(d,p) calculation level, the transition-state constructions are calculated by Berny arithmetic method.59 Based on the same calculation level as optimization, frequency analyses are conducted to ensure that each of the stable geometries is a real minimum, free from imaginary, while the transition-state constructions are confirmed by sole imaginary

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 4 of 25

frequency of specific vibrational mode. The veracity of transition-state constructions can be further confirmed by calculating intrinsic reaction coordinate (IRC).60 The geometries of Anthrols, isomer Anthrols molecule (Anthrols-T) and final product QM molecules are visualized with VMD program.61 Frontier molecular orbitals (FMOs) isosurfaces are visualized with Chemcraft program.62 Mayer bond orders63 of molecule are calculated and analyzed by Multiwfn program.64 The Mayer orders between N atom and C atom are defined as α β I NC = I NC + I NC = 2 ∑∑ [( P α S )ba ( P α S ) ab + ( P β S )ba ( P β S ) ab ]

(1)

a∈ N b∈C

Where P α and P β represents alpha and beta density matrix, respectively, S is overlap matrix, total spin density matrix is P = P α + P β . For restricted close-shell circumstance, total spin density matrix P = 0 , the equation (1) can be written as I NC =

∑∑ ( PS )

ba

( PS ) ab

(2)

a∈ N b∈C

The total valence of N atom in Multiwfn64 is defined as

VN = 2 ∑ ( PS ) aa − ∑ ∑ ( PS )ba ( PS ) ab a∈ N

(3)

a∈ N b∈ N

The free valence is defined as

FN = VN −

∑I

NC

C≠N

= ∑ ∑ ( P S S )ba ( P S S ) ab

(4)

a∈ N b∈ N

For our restrained closed shell molecular system, spin density matrix P S = 0 , the free valence FN of N atom is zero, its total valence is sum of the related Mayer bond orders:

VN =

∑I

NC

(5)

C≠N

Therefore, the bonding capacity of atoms can denote as the Mayer bond order. The three-dimensional potential energy surfaces (3DPES) in different electronic states are built by constraining optimizations method as function of bond O1-H2 and N3-C4, their step size is fixed at 0.1 Å to better clarify the photodeamination mechanism induced by ESIPT reaction. The Gibbs free energy correction values (∆Gcorr) of all stable and transition-state structures are calculated by B3LYP-D3/6-31G(d,p) calculation level upon 1 atm and 298.15 K. In addition, the larger atomic basis set

ACS Paragon Plus Environment

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

def2-tzvp65 is used for eliminating the impact of ampliative basis set and obtaining the more accurate single-point energy (SPE), the formula calculating to Gibbs free energies (G) as followes: G=SPE+∆Gcorr All theoretical calculation DFT/TDDFT methods consider the dispersion correction with ultrafine integration grid (99, 590 points).

Results and discussion Analyses of geometric constructions

To propose the molecular-level explanation about photodeamination mechanism induced by ESIPT reaction, the molecular configurations of Anthrols, isomer Anthrols-T and final product QM have been fully optimized as shown in Figure 1. To conveniently illustrate the changes of bond parameters during photochemical reactions, key atoms are numbered 1-4. The violet imaginary lines represent hydrogen bonds, proton transfer reactions can occur along them. Deamination reactions occur in bond N3-C4 moiety in different electronic states. The corresponding bond parameters are presented in Table 1. For Anthrols molecule, the bond length of O1-H2 is 1.009 Å in S0 state, which increases to 1.029 Å in S1 state, while the bond length of N3-H2 reduces from 1.716 Å in S0 state to 1.644 Å in S1 state. The hydrogen bond’s angle of δ (O1-H2-N3) is also exhibited, which changes from 152.1° to 152.9° from S0 to S1 state. It is concluded that the hydrogen bond O1-H2···N3 is strengthened upon photo-induced process. ESIPT reaction can be facilitated, since the hydrogen bond interaction provides driving force for it. Note that bond length of N3-C4 is 1.481 Å in S0 state, while the length is 1.487 Å in S1 state. The little change of bond N3-C4 is ∆r=0.00551 Å, which indicates that the intensity of bond N3-C4 hardly changes under the photoexcitation process, so photo-excitation

will not affect the deamination

reaction of Anthrols molecule. For isomer Anthrols-T, its instability of structure in S0 state has been verified via theoretical calculation method. Upon photo-induced process the isomer Anthrols-T can be obtained by ESIPT reaction, bond parameters of

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 6 of 25

isomer form occur notable change compared with that of normal Anthrols form in Table 1. Table 1 Bond Lengths (Å) and Angles (°) for the Anthrols, isomer Anthrols-T and final product QM molecules in S0 and S1 state. Anthrols

Anthrols-T

QM

S0

S1

S0

S1

S0

S1

O1-H2

1.00861

1.02883

-

1.85902

-

3.48437

N3-H2

1.71643

1.64426

-

1.04170

-

1.01731

δ(O1-H2-N3)

152.121

152.855

-

139.687

-

105.213

N3-C4

1.48109

1.48660

-

1.52971

-

2.30835

Note that the bond length of N3-C4 on Anthrols-T molecule changes to 1.530 Å from 1.487 Å on Anthrols molecule in S1 state. We can conclude that the intensity of bond N3-C4

became

weaker

on

Anthrols-T

molecule,

which

indicates

that

photoisomerization will promote the photodeamination reaction. For final product QM molecule, we only obtain its construction in S1 state, thereinto, the bond length of O1-H2 and N3-C4 are 3.484 Å and 2.308 Å separately, which indicates that hydrogen bond O1···H2-N3 and covalent bond N3-C4 are interrupted.

Fig. 1 Optimized geometrical configurations of Anthrols, isomer Anthrols-T and final product QM molecules. Key atoms are numbered 1−4. Atom color coding: O, red; N, blue; C, cyan; H, gray.

ACS Paragon Plus Environment

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

Analyses of infrared (IR) vibrational spectral The IR vibrational frequencies of Anthrols, Anthrols-T and final product QM molecules were also calculated at the same theory level as the optimization, and no imaginary frequencies were found. Thus, these forms were confirmed as the real, most stable structures. The shift in the IR vibrational frequency of a chemical bond can usually be taken as an indicator of bond strength.66 As shown in Figure 2 (a), the IR vibrational frequencies of hydrogen-bond O1-H2 moiety are exhibited in S0 and S1 state for normal Anthrols form, which are 2945 cm-1 and 2602 cm-1 in S0 and S1 state, respectively. After the photoexcitation process, the IR vibrational frequency of O1-H2 is of a large red shift with 343 cm-1. It indicates that the intensity of bond O1-H2 weaken in S1 state, leading to the strengthening in hydrogen bond O1-H2···N3. The IR vibrational frequencies of bond N3-C4 of Anthrols and Anthrols-T forms in S1 state are shown in Figure 2 (a), corresponding values of Anthrols and Anthrols-T are 1055.3 cm-1 and 945.4 cm-1, respectively. The IR vibrational frequency of bond N3-C4 has a red shift with 109.9 cm-1, which confirms that the ESIPT reaction causes the weakening of bond N3-C4, and then it is helpful to conduct the deamination.

Fig. 2 Infrared (IR) vibrational frequencies of (a) bond O1−H1 in the S0 and S1 states for Anthrols; (b) bond N3-C4 for Anthrols and isomer Anthrols-T in the S1 state.

Calculated frontier molecular orbitals and Mayer orders

The changes of bond length or IR vibrational frequency in most case serve as the

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 8 of 25

indicators of bond strength. To essentially explain the changes of bond strength for hydrogen bond and N3-C4, the frontier molecular orbitals (FMOs) and Mayer orders of molecule have been calculated. Upon photo-induced process, electron population will redistribute in molecule, further influencing the hydrogen-bond strength originated from electrostatic interaction. In Table 2, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) occupy predominant transition component with 98 % from S0 to S1 state. The relevant oscillator strength (f) is 0.0658, the excitation energy is 411 nm.

Table 2 Transition types with different excitation energies λ (nm) of Anthrols and their oscillator strength (f) and their transition compositions. process

excitation

species

Anthrols

transition

λ (nm)

f

compositions

S0→S1

411

0.0658

H→L (98%)

S0→S2

336

0.0299

H-1→L (60.7%) H→L+1 (36.7%) S0→S3

310

0.0033

H-2→L (97%)

The isosurfaces of FMOs are present in Figure 3, the electron population redistribution is characterized by ππ* in electronic transition process. Note that the electron density of proton donor group decreases locating in region in green ellipse from HOMO to LUMO process, which facilitates hydrogen proton to move to proton acceptor (dimethylamino nitrogen atom) because of the stronger acidity of proton donor group (carbonyl oxygen atom), therefore hydrogen bond O1-H2···N3 is strengthened in photoexcitation process. To investigating effect of ESIPT reaction on intensity of bond N3-C4, we calculate the Mayer orders of Anthrols and Anthrols-T molecules. The value of Mayer order of bond N3-C4 is 0.903 on normal Anthrols molecule, while it decreases to 0.829 on isomer Anthrols-T molecule. Obvious reduction of Mayer orders indicates the intensity of bond N3-C4 weakening through ESIPT reaction, which further facilitates the occurrence of deamination reaction.

ACS Paragon Plus Environment

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

Fig. 3 Frontier molecular orbitals HOMO and LUMO of Anthrols molecule (isovalue = 0.03 a. u.). The green elliptical regions are labeled with hydrogen bond moiety.

Analyses of electron spectra

The electronic transition energies can be calculated by using B3LYP-D3 method and def2-tzvp65 atomic basis set to reproduce experimental photophysical phenomena.44 According to Frank-Condon principle67, the approximation that an electronic transition is most likely to occur without changes in the positions of the nuclei in the molecular entity and its environment. The target state is named as Frank-Condon state, the transition belongs to vertical transition, which means that the calculation of absorption and emission energy based on ground-state geometry and excited-state geometry, respectively.68 On the basis of Kasha’s rule, the processes of fluorescent radiation are only presented from the lowest excited state to the ground state.11 To accurately calculate absorption and fluorescence between S0 and S1 state, we solve three excited states (nstates=3). In Figure 4, the absorption and fluorescence spectra are depicted and broadened by GaussView program based on its default Gaussian and Lorentz function.69 Herein, X and Y axis represent wave length and molar absorption coeffcient, respectively. The value of absorption peak locates in 411 nm, coinciding with value 420 nm of absorption in experiment. Note that emission spectra exhibit a distinct dual-fluorescent phenomenon in Figure 4, the large Stocks’ shift is about 147

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 25

nm that is generated by the ESIPT reaction, some contemporary researchers have attributed this photophysical phenomenon to photo-isomerization processes.70-76 The fluorescence signals of 478 nm and 558 nm are attributed to Anthrols and Anthrols-T molecules, respectively, which are agreement with the corresponding values in experiment. In Table 3, the obtained values of absorption and emission with theoretical computation approaches are compared with corresponding values in experiment. All results within the margin of error of 0.2 eV. It concludes that the theoretical methods that we have carried out are reliable to investigate relevant photochemical reactive mechanisms of Anthrols molecule.

Fig. 4 Absorption and fluorescence spectra of Anthrols and isomer Anthrols-T molecules employing the B3LYP-D3/def2-tzvp calculation level. The large Stokes’ shift is labeled as the difference of spectra between enol absorption and keto emission.

Table 3 Absorption and emission spectra peak values (nm) based on B3LYP-D3/def2-tzvp calculation level and corresponding values (nm) in experiment, the difference between them is represented as ∆λ (nm), E represents electron volt (eV) in parenthesis. B3LYP-D3/def2-tzvp

Theor.

Exp.

∆λ(E)

Anthrols(abs.)

411

420

18(0.13)

Anthrols(flu.)

478

450

28(0.16)

Anthrols-T(flu.)

558

575

17(0.07)

ACS Paragon Plus Environment

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

Analysis of the reaction pathways and mechanisms

To better explain the mechanisms of deamination, the reactive three-dimensional potential energy surfaces (3DPES) as function of bond O1-H2 and N3-C4 are built by restricted optimization method in Figure 5 and Figure 6. Herein, the X and Y axis denote the linear changes of bond O1-H2 and N3-C4, respectively. The Z axis represents single-point energies with the changing of configuration.

Fig. 5 (a) 3D potential-energy surface as functions of bond lengths O1−H2 (from 1.00 to 2.20 Å) and N3−C4 (from 1.47 to 2.67 Å) of Anthrols in the S0 state. The white imaginary arrow represents the infeasible reactions, (b) the corresponding contour map.

Fig. 6 3D potential-energy surface as functions of bond lengths O1−H2 (from 1.00 to 2.20 Å) and N3−C4 (from 1.47 to 2.22 Å) of Anthrols in the S1 state. The white imaginary arrow represents the infeasible reaction; the white solid arrows represent workable reactions, (b) the corresponding contour map.

In Figure 5 (a) the 3DPES indicates that both direct and induced by proton transfer 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

deamination reaction are infeasible in S0 state. The corresponding contour map has been exhibited on Figure 5 (b), note that Anthrols molecule (a) is incapable of generating thermodynamically stable Anthrols-T (b) and final QM (c) in S0 state. Figure 6 (a) shows the 3DPES of photodeamination reaction in S1 state, which indicates that the direct photodeamination reaction of Anthrols molecule (A) that plotted by imaginary line still is infeasible, while it is workable induced by ESIPT reaction. In Figure 6 (b) the isomer Anthrols-T (B) and final product QM (C) are locally most stable points in contour map, the reaction pathway of photodeamination is A→B→C. To quantificationally illustrate occurrence of photodeamination reaction, the Gibbs free energies have been calculated to approximate reactive activation energies. Herein, the transition-state constructions are sought with Berny arithmetic method59, their accuracy is confirmed via analyses of IR vibrational frequencies, that is, there is only one imaginary frequency and its vibrational mode is consistent with the reactive orientation. In addition, we calculate intrinsic reaction coordinate (IRC) to verify the transition-state constructions as shown in Figure 7. In Figure 7 (a) the summit TS1 represents transition-state construction of ESIPT reaction; the IRC connects correct reactant (Anthrols) and product (Anthrols-T). In Figure 7 (b) the summit TS2 represents transition-state construction of photodeamination reaction; the IRC connects correct reactant (Anthrols-T) and product (QM). Gibbs free energy corrected values (∆Gcorr) and single-point energies (SPE) of these constructions are calculated, the Gibbs free energies (G) are calculated with equation G=SPE+∆Gcorr as shown in Table 4. To obviously present the activation energies of ESIPT and photodeamination reactions, relevant schematic diagram is shown in Figure 8.

ACS Paragon Plus Environment

Page 12 of 25

Page 13 of 25 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

Fig. 7 Calculated Intrinsic Reaction Coordinate (IRC) of transition states, (a) ESIPT reaction, (b) photodeamination reaction.

Table 4 Single-point energies (SPE), Gibbs free energy corrected values (∆Gcorr) and Gibbs free energies (G) of minimum points and transition-state constructions. The all numerical values are of atomic unit (a. u.) Anthrols

TS1

(A)

Anthrols-T

TS2

(B)

QM (C)

SPE

-866.94240

-866.94130

-866.95290

-866.94330

-866.94690

∆Gcorr

0.307643

0.303453

0.309645

0.306546

0.302740

G=SPE+∆Gcorr

-866.63476

-866.63785

-866.64326

-866.63676

-866.64420

The free energy barrier ∆G1 of ESIPT reaction (A→B) is -1.95 Kcal/mol, it is a common phenomenon that obtained negative value is due to use the Gibbs free energy correction.77-80 The potential barrier of ESIPT reaction with 0.69 Kcal/mol is obtained by calculating single-point energies. The free energy barrier ∆G1r of reversed ESIPT

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

reaction (B→A) is 3.40 Kcal/mol. For photodeamination reaction (B→C), free energy barrier ∆G2 is 4.13 Kcal/mol, its reversed (C→B) free energy barrier ∆G2r is 4.69 Kcal/mol. Therefore, it concludes that the photodeamination reaction induced by ESIPT reaction occurs only in S1 state.

Fig. 8 Schematic diagram of ESIPT and photodeamination reaction, relevant reactive free energy barriers are calculated.

Conclusion In this work, the mechanism of ESIPT-induced photodeamination reaction is given a molecular-level explanation based on theoretical research. Absorption and emission spectra coincide with observed results in experiment. The analyses of bond parameters and IR vibrational frequencies of different constructions confirm that the intramolecular hydrogen bond O1-H2···N3 is strengthened upon photoexcitation, facilitating the occurrence of ESIPT reaction, this conclusion also confirms that the photoisomerization further weaken the strength of bond N3-C4, facilitating its

ACS Paragon Plus Environment

Page 14 of 25

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

breakage. To essentially explain changes of strength of hydrogen bond O1-H2···N3 and bond N3-C4, the FMOs and Mayer bond orders of Anthrols molecule are calculated. Because of electron density redistribution upon photoexcitation process, enhancing acidity of donor group facilitate hydrogen proton to transfer to acceptor group, which represents as hydrogen bond O1-H2···N3 strengthening. The weaker bonding capacity of N3 and C4 atoms is exhibited in Anthrols-T molecule, since Mayer bond orders of bond N3-C4 is decreasing compared with that in Anthrols molecule, which indicates that the strength of bond N3-C4 is weaken by ESIPT reaction. According to the constructed 3DPES and calculated reactive Gibbs free energies barriers, we can clearly indicate that deamination or photodeamination reactions cannot direct occur wherever in S0 or S1 state, while the photodeamination reaction can be induced only by ESIPT reaction.

Acknowledgements This work was supported by the National Natural Science Foundation of China (Grant No. 11474141 and 11274149), the Natural Science Foundation of Liaoning Province (Grant No. 20170540408), and the Science and Technology Plan Project of Shenyang City (Grant No.17-231-1-06).

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

References (1) Zhao G. J.; Han, K. L. Hydrogen Bonding in the Electronic Excited State. Acc. Chem. Res. 2011, 45, 404-413. (2) Santos, F. S.; Ramasamy, E.; Ramamurthy, V.; Rodembusch, F. S. Excited State Behavior of Benzoxazole Derivatives in a Confined Environment Afforded by a Water Soluble Octaacid Capsule. J Photoch. Photobio. A 2016, 317, 175-185. (3) Vazquez, S. R.; Rodriguez, M. C. R.; Mosquera, M.; Rodriguez-Prieto, F. Rotamerism, Tautomerism, and Excited-State Intramolecular Proton Transfer in 2-(4‘-N,N-Diethylamino-2‘-hydroxyphenyl)benzimidazoles:  Novel Benzimidazoles Undergoing Excited-State Intramolecular Coupled Proton and Charge Transfer. J. Phys. Chem. A 2008, 112, 376-387. (4) Serdiuk, I. E. White Light from a Single Fluorophore: A Strategy Utilizing Excited-State Intramolecular Proton-Transfer Phenomenon and Its Verification. J. Phys. Chem. C 2017, 121, 5277-5286. (5) Lan, R. F.; Yang, Y. F.; Ma, Y. Z.; Li, Y. Q. The Theoretical Study of Excited-State Intramolecular Proton Transfer of 2, 5-bis (Benzoxazol-2-yl) Thiophene-3, 4-Diol. Spectrochim. Acta A 2017, 183, 37-44. (6) Yang, Y. F.; Zhao. J. F.; Li, Y. Q. Theoretical Study of the ESIPT Process for a New Natural Product Quercetin Sci. Rep. 2016, 6, 32152. (7) Chu, T. S.; Liu, B. T. Establishing New Mechanisms with Triplet and Singlet Excited-State Hydrogen Bonding Roles in Photoinduced Liquid Dynamics. Int. Rev. Phys. Chem. 2016, 35, 187-208. (8) Manna, A.; Sarkar, D.; Goswami, S.; Quah. C. K.; Fun, H. K. Single Excited State Intramolecular Proton Transfer (ESIPT) Chemodosimeter Based on Rhodol for both Hg2+ and OCl−: Ratiometric Detection with Live-Cell Imaging. RSC Adv. 2016, 6, 57417-57423. (9) Shvadchak, V. V.; Falomir-Lockhart, L. J.; Yushchenko D. A.; Jovin, T. M. Specificity and Kinetics of α-Synuclein Binding to Model Membranes Determined with Fluorescent Excited State Intramolecular Proton Transfer (ESIPT) Probe. J. Biol. ACS Paragon Plus Environment

Page 16 of 25

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

Chem. 2011, 286, 13023-13032. (10) Stasyuk, A. J.; Cyranski, M. K.; Gryko, D. T.; Sola, M. Acidic C–H Bond as a Proton Donor in Excited State Intramolecular Proton Transfer Reactions. J. Chem. Theory Comput. 2015, 11, 1046-1054. (11) Kasha, M.; Rawls, H. R.; Elbayoumi, M. A.; The Exciton Model in Molecular Spectroscopy. Pure Appl. Chem. 1965, 11, 371-392. (12) Lakowicz, J. R. Principles of fluorescence spectroscopy, 3rd edn; Springer: New York Google Scholar, 2006. (13) Mamada, M.; Inada, K.; Komino, T.; Jr, W. J. P.; Nakanotani, H.; Adachi, C. Highly Efficient Thermally Activated Delayed Fluorescencefrom an Excited-State Intramolecular Proton Transfer System. Acs Central Sci. 2017, 3, 769-777. (14) Padalkar, V. S.; Seki, S. Excited-State Intramolecular Proton-Transfer (ESIPT)-Inspired Solid State Emitters. Chem. Soc. Rev. 2015, 45, 169-202. (15) Wu, J. S.; Liu, W. M.; Ge, J. C.; Zhang, H. Y.; Wang, P. F. New Sensing Mechanisms for Design of Fluorescent Chemosensors Emerging in Recent Years. Chem. Soc. Rev. 2011, 40, 3483-3495. (16) Bera, M. K.; Chakraborty, C.; Singh, P. K.; Sahu, C.; Sen, K.; Maji, S.; Das, A. K.; Malik, S. Fluorene-Based Chemodosimeter for “Turn-on” Sensing of Cyanide by Hampering ESIPT and Live Cell Imaging. J. Mater. Chem. B 2014, 2, 4733-4739. (17) Geng, L. H.; Yang, X. F.; Zhong, Y. G.; Li, Z.; Li, H. “Quinone-Phenol” Transduction Activated Excited-State Intramolecular Proton Transfer: A New Strategy toward Ratiometric Fluorescent Probe for Sulfite in Living Cells. Dyes Pigments.

2015, 120, 213-219. (18) Lv, H. M.; Chen, Y.; Lei, J.; Au, C. T.; Yin, S. F. An ESIPT-Based Ratiometric Fluorescent Probe for the Imaging of Nitroxyl in Living Cells. Anal. Methods. 2015, 7, 3883-3887. (19) Budzak, S.; Jacquemin, D. Mechanism of Fluorescence Switching in One ESIPT-Based Al3+ Probe. J. Phys. Chem. B 2016, 120, 6730-6738. (20) Escudero, D. Revising Intramolecular Photoinduced Electron Transfer (PET) from First-Principles. Acc. Chem. Res. 2016, 49, 1816-1824. 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

(21) Li, G. Y.; Zhao, G. J.; Liu, Y. H.; Han, K. L.; He, G. Z. TD-DFT Study on the Sensing Mechanism of a Fluorescent Chemosensor for Fluoride: Excited-State Proton Transfer. J. Comput. Chem. 2010, 31, 1759-1765. (22) Sakai, K.; Ishikawa, T.; Akutagawa, T. A Blue-White-Yellow Color-Tunable Excited State Intramolecular Proton Transfer (ESIPT) Fluorophore: Sensitivity to Polar-Nonpolar Solvent Ratios. J.Mater. Chem. C 2013, 1, 7866-7871. (23) Sakai, K. I.; Ichikawa, M.; Taniguchi, Y. Photoluminescent Mechanism of a Proton-Transfer Laser Dye in Highly Doped Polymer Films. Chem. Phys. Lett. 2006, 420, 405-409. (24) Sakai, K. I.; Tsuzuki, T.; Itoh, Y.; Ichikawa, M.; Taniguchi, Y. Using proton-transfer laser dyes for organic laser diodes. Appl. Phys. Lett. 2005, 86, 081103. (25) Weller, A. Using Proton-Transfer Laser Dyes for Organic Laser Diodes. Naturwiss. 1955, 42, 175-176. (26) Fischer, M.; Wan, P. Nonlinear Solvent Water Effects in the Excited-State (Formal) Intramolecular Proton Transfer (ESIPT) in m-Hydroxy-1,1-diaryl Alkenes:  Efficient Formation of m-Quinone Methides. J. Am. Chem. Soc. 1999, 121, 4555-4562. (27) Ameer-Beg, S.; Ormson, S. M.; Brown, R. G.; Matousek, P.; Towrie, M.; Nibbering, E. T. J.; Foggi, P.; Neuwahl, F. V. R. Ultrafast Measurements of Excited State Intramolecular Proton Transfer (ESIPT) in Room Temperature Solutions of 3-Hydroxyflavone and Derivatives. J. Phys. Chem. A 2001, 105, 3709-3718. (28) Doroshenko, A. O.; Posokhov, E. A.; Verezubova, A. A.; Ptyagina, L. M.; Skripkina, V. T.; Shershukov, V. M. Radiationless Deactivation of the Excited Phototautomer Form and Molecular Structure of ESIPT-Compounds. Photochem. Photobiol. Sci. 2002, 1, 92-99. (29) Lim, S. J.; Seo J.; Park, S. Y. Photochromic Switching of Excited-State Intramolecular Proton-Transfer (ESIPT) Fluorescence:  A Unique Route to High-Contrast Memory Switching and Nondestructive Readout. J. Am. Chem. Soc.

2006, 128, 14542-14547. (30) Cardoso, M. B.; Samios, D.; da Silveira, N. P.; Rodembusch, F. S.; Stefani, V. ACS Paragon Plus Environment

Page 18 of 25

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

ESIPT-Exhibiting Protein Probes: A Sensitive Method for Rice Proteins Detection during Starch Extraction. Photochem. Photobiol. Sci. 2007, 6, 99-102. (31) Li, G. Y.; Chu, T. S. TD-DFT Study on Fluoride-Sensing Mechanism of 2-(2′ -Phenylureaphenyl) Benzoxazole: The Way to Inhibit the ESIPT Process. Phys. Chem. Chem. Phys. 2011, 13, 20766-20771. (32) Chevalier, K.; Wolf, M. M.; Funk, A.; Andres, M.; Gerhards, M.; Diller, R. Transient IR Spectroscopy and ab Initio Calculations on ESIPT in 3-Hydroxyflavone Solvated in Acetonitrile. Phys. Chem. Chem. Phys. 2012, 14, 15007-15020. (33) Ciuciu, A. I.; Flamigni, L.; Skonieczny, K.; Gryko, D. T. Blue-Green Emitting Sulphonamido-Imidazole Derivatives: ESIPT Based Excited State Dynamics. Phys. Chem. Chem. Phys. 2013, 15, 16907-16916. (34) Goswami, S.; Manna, A.; Paul, S.; Das, A. K.; Nandi, P. K.; Maity, A. K.; Saha, P. A Turn on ESIPT Probe for Rapid and Ratiometric Fluorogenic Detection of Homocysteine and Cysteine in Water with Live Cell-Imaging. Tetrahedron Lett. 2014, 55, 490-494. (35) Zhao, J. F.; Song, P.; Ma, F. C. Theoretical Investigation on ESIPT Process of 3-(1,

3-Benzothiazol-2-yl)-2-Hydroxynaphthalene-1-Carbaldehyde

Chemosensor

between Polar and Non-Polar Solvent. Commun. Comput. Chem. 2015, 3, 44-54. (36) Chen, L.; Zhang, H.; Pan, M.; Wei, Z. W.; Wang, H. P.; Fan, Y. N.; Su, C. Y. An Efficient Visible and Near-Infrared (NIR) Emitting SmIII Metal-Organic Framework (Sm-MOF) Sensitized by Excited-State Intramolecular Proton Transfer (ESIPT) Ligand. Chem. Asian J. 2016, 11, 1765-1769. (37) Zhou, P. W.; Hoffmann, M. R.; Han, K. L.; He, G. Z. New Insights into the Dual Fluorescence of Methyl Salicylate: Effects of Intermolecular Hydrogen Bonding and Solvation. J. Phys. Chem. B 2015, 119, 2125-2131. (38) Zhao, J. F.; Chen, J. S.; Liu, J. Y.; Hoffmann, M. R. Competitive Excited-State Single

or

Double

Proton

Transfer

Mechanisms

for

bis-2,

5-(2-Benzoxazolyl)-Hydroquinone and Its Derivatives. Phys. Chem. Chem. Phys.

2015, 17, 11990-11999. (39) Li, G. Y.; Han, K. L. The Sensing Mechanism Studies of the Fluorescent Probes 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 20 of 25

with Electronically Excited State Calculations. WIREs Comput. Mol. Sci. 2017, e1351. (40) Goswami, S.; Maity, S.; Das, A. K.; Maity, A. C.; Mandal, T. K.; Samanta, S. Remarkable

ESIPT

Induced

NIR

Emission

by

a

Selective

Colorimetric

Dibenzimidazolo Diimine Sensor for Acetate. Tetrahedron Lett. 2013, 54, 5232-5235. (41) Wang, Y. H.; Wan, P. Excited State Intramolecular Proton Transfer (ESIPT) in Dihydroxyphenyl Anthracenes. Photochem. Photobiol. Sci. 2011, 10, 1934-1944. (42) Škalamera, D.; Bohne, C.; Landgraf, S.; Basaric, N. Photodeamination Reaction Mechanism in Aminomethyl p-Cresol Derivatives: Different Reactivity of Amines and Ammonium Salts. J. Org. Chem. 2015, 80, 10817-10828. (43) Yang, Y. F.; Ma, Y. Z.; Zhao, Y.; Zhao Y. L.; Li, Y. Q. Theoretical Investigation of the

Reaction

Mechanism

of

Photodeamination

Induced

by

Excited-State

Intramolecular Proton Transfer of Cresol Derivatives. J. Phys. Chem. A 2017, 122, 1011-1018. (44) Škalamera, D.; Mlinarićmajerski, K.; Martin, K. I.; Kralj, M.; Oake, J.; Wan, P.; Bohne, C.; Basaric, N. Photochemical Formation of Anthracene Quinone Methide Derivatives. J. Org. Chem. 2017, 82, 6006-6021. (45) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378-6396. (46) Pal, S.; Joy, S.; Paul, H.; Banerjee, S.; Maji, A.; Zangrando. E.; Chattopadhyay, P. Understanding the Difference in Photophysical Properties of Cyclometalated Iridium(III) and Rhodium(III) Complexes by Detailed Time-Dependent Density Functional Theory and Frontier Molecular Orbital Supports. J. Phys. Chem. C 2017, 121, 11632-11642. (47) Li, H.; Shi, Y.; Yin, H.; Wang, Y.; Cong, L.; Jin, M. X.; Ding, D. J. New Insights into the Solvent-Assisted Excited-State Double Proton Transfer of 2-(1H-Pyrazol-5-yl) Pyridine with Alcoholic Partners: A TDDFT Investigation. Spectrochim. Acta A 2015, 141, 211-215. ACS Paragon Plus Environment

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

(48) Liu, Y. H.; Mehata, M. S.; Lan, S. C. TDDFT Study of the Polarity Controlled Ion-Pair Separation in an Excited-State Proton Transfer Reaction. Spectrochim. Acta A

2014, 128, 280-284. (49) Fan, G. H.; Han, K. L.; He, G. Z. Time-Dependent Density Functional-Based Tight-bind Method Efficiently Implemented with OpenMP Parallel and GPU Acceleration. Chinese J. Chem. Phys. 2013, 26, 635-645. (50) Zhao, X. H.; Chen, M. D. A TDDFT Study on the Singlet and Triplet Excited-State Hydrogen Bonding and Proton Transfer of 10-Hydroxybenzo[h] Quinoline (HBQ) and 7,9-Diiodo-10-Hydroxybenzo[h] Quinoline (DIHBQ). Chem. Phys. Lett. 2011, 512, 35-39. (51) Zhao, G. J.; Chen, R. K.; Sun, M. T.; Liu, J. Y.; Li, G. Y.; Gao, Y. L.; Han, K. L.; Yang, X. C.; Sun, L. Photoinduced Intramolecular Charge Transfer and S2 Fluorescence in Thiophene-π-Conjugated Donor-Acceptor Systems: Experimental and TDDFT Studies. Chem-Eur. J 2008, 14, 6935-6947. (52) Zhao, J. Z.; Ji, S. M.; Chen, Y. H.; Guo, H. M.; Yang, P. Excited State Intramolecular Proton Transfer (ESIPT): From Principal Photophysics to the Development of New Chromophores and Applications in Fluorescent Molecular Probes and Luminescent Materials. Phys. Chem. Chem. Phys. 2012, 14, 8803-8817. (53) Li, Y. Q.; Ma, Y. Z.; Yang, Y. F.; Shi, W.; Lan, R. F.; Guo, Q. Effects of Different Substituents of Methyl 5-R-Salicylates on the Excited State Intramolecular Proton Transfer Process. Phys. Chem. Chem. Phys., 2018, 20, 4208-4215. (54) Zhao, G. J.; Han, K. L. Novel Infrared Spectra for Intermolecular Dihydrogen Bonding of the Phenol-Borane-Trimethylamine Complex in Electronically Excited State. J. Chem. Phys. 2007, 127, 024306. (55) Zhao, G. J.; Han, K. L. Ultrafast Hydrogen Bond Strengthening of the Photoexcited Fluorenone in Alcohols for Facilitating the Fluorescence Quenching. J. Phys. Chem. A 2007, 111, 9218-9223. (56) Zhao, G. J.; Han, K. L. Early Time Hydrogen-Bonding Dynamics of Photoexcited Coumarin 102 in Hydrogen-Donating Solvents: Theoretical Study. J. Phys. Chem. A

2007, 111, 2469-2474. 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

(57) Ma, Y. Z.; Yang, Y. F.; Lan, R. F.; Li, Y. Q. Effect of Different Substituted Groups on Excited-State Intramolecular Proton Transfer of 1-(Acylamino)-Anthraquinons. J. Phys. Chem. C 2017, 121, 14779-14786. (58) 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 D. 1; Gaussian, Inc.: Wallingford, CT, 2010. (59) Schlegel, H. B. Geometry Optimization on Potential Energy Surfaces. J. Comput. Chem. 1982, 3, 214-218. (60) Fukui, K. The Path of Chemical Reactions-the IRC Approach. Accounts Chem. Res. 1981, 14, 363-368. (61) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J Mol. Graph. Model. 1996, 14, 33-38. (62) Andrienko, G. A. Chemcraft. http://www.chemcraftprog.com. (accessed April 7, 2018). (63) Mayer, I. Charge, Bond Order and Valence in the AB Initio SCF Theory. Chem. Phys. Lett. 1985, 97, 270-274. (64) Lu, T. Chen, F. W. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem. 2012, 33, 580-592. (65) Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297-3305. (66) Karas, L. J.; Batista, P. R.; Viesser, R. V.; Tormena, C. F.; Rittner, R.; de Oliveira, P. R. Trends of Intramolecular Hydrogen Bonding in Substituted Alcohols: A Deeper Investigation. Phys. Chem. Chem. Phys. 2017, 19, 16904-16913. (67) Coolidge, A. S.; James, H. M.; Present, R. D. A Study of the Franck-Condon Principle. J. Chem. Phys. 1936, 4, 193-211. (68) Adamo, C.; Jacquemin, D. The Calculations of Excited-State Properties with Time-Dependent Density Functional Theory. Chem. Soc. Rev. 2013, 42, 845-856. (69) Dennington, R.; Keith, T. A.; Millam, J. M. GaussView, Version 6; Semichem Inc.: Shawnee Mission, KS, 2016. ACS Paragon Plus Environment

Page 22 of 25

Page 23 of 25 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

(70) Zhao, J. F.; Chen, J. S.; Cui, Y. L.; Wang, J.; Xia, L. X.; Dai, Y. M.; Song, P.; Ma, F. C. A Questionable Excited-State Double-Proton Transfer Mechanism for 3-Hydroxyisoquinoline. Phys. Chem. Chem. Phys. 2015, 17, 1142-1150. (71) Zhao, G. J.; Liu, J. Y.; Zhou, L. C.; Han, K. L. Site-Selective Photoinduced Electron Transfer from Alcoholic Solvents to the Chromophore Facilitated by Hydrogen Bonding:  A New Fluorescence Quenching Mechanism. J. Phys. Chem. B

2007, 111, 8940-8945. (72) Zhao, J. F.; Yao, H. B.; Liu, J. Y.; Hoffmann, M. R. New Excited-State Proton Transfer Mechanisms for 1,8-Dihydroxydibenzo [a,h] phenazine. J. Phys. Chem. A

2015, 119, 681-688. (73) Zhao, Y.; Yang, Y. F.; Ma, Y. Z.; Li, Y. Q. Stimuli-Responsive Luminescent Coumarin Thiazole Hybrid Dye: Mechanism of Excited-State Intramolecular Double Proton Transfer. J. Lumin. 2018, 201, 189-195. (74) Zhao, J. F.; Liu, X. Y.; Zheng, Y. J. Controlling Excited State Single versus Double Proton Transfer for 2,2′-Bipyridyl-3,3′-diol: Solvent Effect. J. Phys. Chem. A

2017, 121, 4002-4008. (75) Zhao, J. F.; Dong, H.; Zheng, Y. J. Theoretical Insights into the Excited State Double Proton Transfer Mechanism of Deep Red Pigment Alkannin. J. Phys. Chem. A

2018, 122, 1200-1208. (76) Ma, Y. Z.; Yang, Y. F.; Shi, W.; Song, Y. Z.; Li, Y. Q. The Order of Multiple Excited State Proton Transfer in Ternary Complex of Norharmane and Acetic Acids. Spectrochim. Acta A 2018, 202, 30-35. (77) Zhao, Y. X.; Wu, X. N.; Ma, J. B.; He, S. G.; Ding, X. L. Experimental and Theoretical Study of the Reactions between Vanadium-Silicon Heteronuclear Oxide Cluster Anions with n-Butane. J. Phys. Chem. C 2010, 114, 12271-12279. (78) Feyel, S.; Döbler, J.; Höckendorf, R.; Beyer, M. K.; Sauer, J.; Schwarz, H. Activation of Methane by Oligomeric (Al2O3)x+ (x= 3, 4, 5): The Role of Oxygen-Centered Radicals in Thermal Hydrogen-Atom Abstraction. Angew. Chem.

2008, 47, 1946-1950. (79) Ma, J. B.; Wu, X. N.; Zhao, X. X.; Ding, X. L.; He, S. G. Methane Activation by 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

V 3 PO 10˙+ and V 4 O 10˙+ Clusters: A Comparative Study. Phys. Chem. Chem. Phys. 2010, 12, 12223-12228. (80) Wu, X. N.; Zhao, Y. X.; Xue, W.; Wang, Z. C.; He, S. G.; Ding, X. L. Active Sites of Stoichiometric Cerium Oxide Cations (Ce m O 2m+) Probed by Reactions with Carbon Monoxide and Small Hydrocarbon Molecules. Phys. Chem. Chem. Phys.

2010, 12, 3984-3997.

ACS Paragon Plus Environment

Page 24 of 25

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

Table of Contents (TOC) Graphic

ACS Paragon Plus Environment