TDDFT Study on the Sensing Mechanism of a ... - ACS Publications

Jun 26, 2017 - Wang , R. Two's Company, Three's a Crowd: Can H2S be the Third Endogenous Gaseous Transmitter? FASEB J. 2002, 16, 1792– 1798 DOI: ...
1 downloads 0 Views 3MB Size
Article pubs.acs.org/JPCA

DFT/TDDFT Study on the Sensing Mechanism of a Fluorescent Probe for Hydrogen Sulfide: Excited State Intramolecular Proton Transfer Coupled Twisted Intramolecular Charge Transfer Yang Li†,§ and Tian-Shu Chu*,†,‡ †

State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China ‡ Institute for Computational Sciences and Engineering, Laboratory of New Fiber Material and Modern Textile, the Growing Base for State Key Laboratory, Qingdao University, Qingdao 266071, P. R. China § University of Chinese Academy of Sciences, Beijing 100049, P. R. China ABSTRACT: By using density functional theory (DFT) and time-dependent density functional theory (TDDFT) methods, the sensing mechanism of a fluorescent probe 2-(2hydroxyphenyl) benzothiazole (HBT) derivative HBTPP-S for hydrogen sulfide has been thoroughly studied. The thiolysis reaction has a moderate reaction barrier of 18.40 kcal mol−1, which indicates that the hydrogen sulfide sensing process has a favorable response speed. Because of the nonradiative donor-excited photoinduced electron transfer (d-PET, fluorophore as the electron donor) from the excited HBTPP group to the electron-withdrawing 2,4-dinitrophenyl group, as well as the inhibition of the proton transfer (PT) and the excited state intramolecular proton transfer (ESIPT) process by 2,4dinitrophenyl group, the probe HBTPP-S is essentially nonfluorescent. On the other hand, the added hydrogen sulfide induces the thiolysis of the 2,4-dinitrophenyl ether bond, and then the thiolysis product HBTPP comes into existence. The theoretically simulated potential energy surface demonstrates that without the electron-withdrawing 2,4-dinitrophenyl group, the thiolysis product HBTPP undergoes the excited state intramolecular proton transfer (ESIPT) coupled twisted intramolecular charge transfer (TICT) processes in the first excited state. The absence of the d-PET and the process mentioned above may explain the significant fluorescent turn-on response and large Stokes shift of the thiolysis product HBTPP.

1. INTRODUCTION In spite of its unpleasant smell and poisonous attribute, hydrogen sulfide (H2S) has been deeply researched due to its functions as the third gasotransmitter following nitric oxide (NO) and carbon monoxide (CO)1−4 in recent years. As an important endogenous antioxidant, H2S is of great significance to various pathophysiological processes, including regulation of cell growth, mediation of neurotransmission, and inhibition of insulin signaling.5−7 However, deregulation of H2S in cells will cause Alzheimer’s and Down’s syndrome, diabetes, and liver cirrhosis.8−11 Therefore, efficient detection of H2S in living system has become an important subject in chemical and biological field. Lately, scientists have reported several sorts of fluorescent probes, which are based on the H2S-triggered specific reactions, such as Cu2+ quencher removal,12−14 reduction of azide15−17 or nitro,18,19 nucleophilic reaction,20−22 thiolysis of dinitrophenyl ether23−25 and so on,26−28 to achieve fluorescent turn-on or ratiometric response. Among these fluorescent probes, most fluorophores show small Stokes shift, which restrains the potential of their biological applications. Considering this factor, it is highly desirable to design a fluorophore with large Stokes shift. © 2017 American Chemical Society

With this, Zeng and co-workers have synthesized 2-(2-hydroxyphenyl) benzothiazole (HBT) derivative HBTPP-S as a novel H2S probe based on thiolysis of 2,4-dinitrophenyl ether and the “turn-on” excited state intramolecular proton transfer (ESIPT) coupled intramolecular charge transfer (ICT) mechanism (see Scheme 1).29 The probe HBTPP-S is nonfluorescent and the thiolysis product HBTPP shows remarkably red keto emission (620 nm) and large Stokes shift (240 nm), which could be used for biological imaging in living cells.30,31 However, we are not clear about the mechanism of the fluorescence quenching of the probe HBTPP-S, and whether the proton transfer (PT) process would occur in the ground state of the thiolysis product HBTPP? Moreover, depending on the reaction time domain, the ESIPT and ICT can occur separately or concomitantly and can be classified into three categories (PT* and CT* are the excited state proton transfer form and the excited state charge transfer form, respectively): (1) PT* → CT* (the ESIPT takes prior to ICT). (2) CT* → PT* (the rate of ICT is faster than that of ESIPT), Received: March 20, 2017 Revised: June 20, 2017 Published: June 26, 2017 5245

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A Scheme 1. Proposed Reaction Mechanism of the Probe HBTPP-S with H2S

Figure 1. Optimized geometries for the probe HBTPP-S and its thiolysis product HBTPP in acetonitrile.

is appended in the ortho-position of phenol through a carbon− carbon double bond, therefore, some of the issues worth studying include whether the structures of HBTPP-S and HBTPP are twisted by the rotating carbon−carbon double bond or not. Does the twisting motion affect the fluorescence emission or not? To answer the above questions to provide a clear idea about the sensing mechanism of this H2S probe, in this work, we carried out calculations based on the density functional theory (DFT)57 and time-dependent density functional theory (TDDFT).58 First, we studied the properties of the ground state and the first excited state of HBTPP-S, HBTPP, and other molecular structures involved in the sensing process. Then, based on these optimized structures, we theoretically simulated and reproduced the experimental absorption and emission spectra. Finally, we calculated the potential energy curves in the S0 and S1 states of HBTPP to explore the PT, ESIPT, and ICT processes and also analyzed the frontier molecular orbitals to explain the sensing mechanism. The DFT/TDDFT methods are suitable to study the mechanism involving PT, ESIPT, and ICT process.59−74 For instance, by using these calculation methods, our group has successfully proposed the mechanism for thiazolidinedione derivative in dimethylformamide solution that TICT could facilitate ESIPT initiated by intermolecular hydrogen-bond strengthening in the S1 state.62 Zhao’s group clearly demonstrated that for N,N-dimethylanilino-1,3-diketone (DMADK) chromophore, the ESIPT and TICT process were coupled in the S1 state rather than sequential ESIPT and TICT process,71 and these calculated results are consistent with those reported by Ciofini’s group.75

and (3) the rates of ESIPT and ICT events are about the same.32−35 Recently, the hydroxyphenylbenzazole family has been studied due to their ability to lead to ESIPT process, such as 2-(2′-hydroxyphenyl)benzoxazole (HBO)36−40 and 2-(2′-hydroxyphenyl)benzothiazole (HBT).41−49 Researchers have found that with the presence of an electron acceptor or electron donor substituent in the molecular framework, HBO and HBT show quite different mechanisms for ESIPT and ICT in the S1 state. When HBO and HBT have an electron acceptor substituent (−NO2, −COOH, −COOR or −RCN), the ESIPT reaction occurs prior to the ICT process.50,51 Nevertheless, when HBO and HBT have an electron donor group (−NH2, NR2, or −CH3), the ICT process is followed by the ESIPT reaction.52 Here, for the thiolysis product HBTPP, it is interesting to investigate whether the ESIPT and ICT processes occur separately or simultaneously? If the former, which comes first? Furthermore, among different ICT models, twisted ICT (TICT) and planar ICT (PICT) have attracted intensive interest and have been extensively reported.53−55 In the TICT state, the molecular configuration is twisted by rotating the acceptor group or donor group, whereas the acceptor and donor group are coplanar and conjugated in the PICT state. The molecule with donor and acceptor group linked by a π-conjugated bridging group may bear torsional relaxation in the excited ICT state, which could cause different degrees of charge separation depending on the torsional angle. More importantly, the fluorescence emission of such system is always influenced by the torsional relaxation process.56 Since, in the probe HBTPP-S and the product HBTPP designed by Zeng and co-workers,29 an electron acceptor group (methylpyridine) 5246

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A

2. METHODS All calculations were carried out by using the Gaussian 09 program package.76 By applying DFT/TDDFT methods, optimal geometries were computed without any symmetry constraint of HBTPP-S and HBTPP in the ground electronic state and in the first excited electronic state. A number of functionals, including B3LYP,77−79 B3PW91,80 PBEPBE,81 MPW1PW91,82 and CAM-B3LYP83 were pretested before we carried out this study, and we found that the long-rangecorrected (LRC) CAM-B3LYP functional shows the most satisfactory agreement with the experimental results. Therefore, the CAM-B3LYP functional was chosen to study the properties of the ground state and the first excited state in the DFT/TDDFT method. The TZVP84,85 was also chosen as the basis set during both the ground state and the first excited state optimizations, and the subsequent frequency analysis ensured that our obtained structures had no imaginary vibration and were located at the local minima in the potential energy surfaces. The intermediate and the transition states were found by scanning the potential energy surface. The frequency analysis confirmed that the intermediate has no imaginary frequency and the transition state has only one imaginary frequency, and the intrinsic reaction coordinate (IRC) calculations make sure that the transition state connects the intermediates and the proper products. The free energy profile for the thiolysis reaction was computed with Gibbs free energy correction. Note that all Gibbs free energies computed were relative energies using the reagent energy as a benchmark. In all calculations, unless there is a special declaration, the solvent effects are included using the integral equation formalism (IEF)86,87 version of polarizable continuum model (PCM)88 and the dielectric constant of acetonitrile (ε = 35.69).

Figure 2. Potential energy curves of S1 state (red line) for the thiolysis product HBTPP calculated at the TDDFT/CAM-B3LYP/TZVP level as functions of the O−H bond length (ESIPT coordinate). The energies of S0 state (black line) were calculated under the geometries of the corresponding S1 state.

O−H bond length has been calculated to aid the study of the PT process. As shown in Figure 2, the transfer of the proton H from O to N in the S0 state exhibits a barrier of 4.24 kcal mol−1 with the O−H bond length increasing and the PT process of HBTPP is endothermic. This proves that, in the ground state, the proton of the thiolysis product HBTPP would not transfer from the proton donor (OH group) to the proton acceptor (neighboring nitrogen group). Considering that the absorption and fluorescence titration experiments were done in a mixed solution of NaSH (a commonly employed H2S donor) and CH3CN/PBS,29 the thiolysis reaction should take place between the probe HBTPP-S and HS−. For further comprehension of this thiolysis dynamics process, we computed the ground state potential energy curve of the probe HBTPP-S and HS− as a function of the S···C distance and displayed the free energy profile for the thiolysis reaction in acetonitrile solvent. As shown in Figure 3, the formation of the intermediate (a complex formed by HBTPP-S and HS−) is endergonic by 2.86 kcal mol−1. The reaction barrier for the S−C bond formation and the C−O bond cleavage is 18.40 kcal mol−1. The moderate reaction barrier (18.40 kcal mol−1) implies that the thiolysis reaction between HBTPP-S and HS− would have a favorable reaction rate, hence, after the formation of HBTPP−, HBTPP− would continue to capture one free H+ from the mixed solution to form the end product HBTPP of the thiolysis reaction.63 On the basis of our optimized ground-state structures, we calculated the low-lying singlet electronic transition energies and the corresponding oscillator strengths of HBTPP-S and HBTPP(E). Our calculated results are demonstrated in Table 1 and Figure 4b. For the probe HBTPP-S, the calculated absorption peak is located at 316 nm and an intense S0 → S1 transition for the thiolysis product HBTPP(E) is located at 347 nm. The calculated bathochromic shift (31 nm) agrees with the experimental results where upon addition of NaHS, the red-shifted new absorption peak is observed at 380 nm (in Figure 4a) and the initial absorption peak at 320 nm decreases.29 We also calculated the UV−vis absorption spectra of HBTPP−. The result displayed that the first singlet transition

3. RESULTS AND DISCUSSION 3.1. The Ground State Structures and Absorption Spectra. The equilibrium geometric structures of the probe HBTPP-S and the thiolysis product HBTPP were displayed in Figure 1. For HBTPP-S and HBTPP, the dihedral angle C1−C2−C3−C4 is respectively 136.21° and 140.88° in the S0 state, which indicates that both HBTPP-S and HBTPP have the twisted geometry by rotating the ethylenic double bond and such torsional motion is usually accompanied by the ICT process. As illustrated, for the probe HBTPP-S, the phenolic hydroxyl group is masked by 2,4-dinitrophenyl group and is converted to the ether group, which induce that the proton transfer process cannot proceed in the ground state and the first excited state. After formation of the product HBTPP through thiolysis of 2,4-dinitrophenyl ether bond, the O−H bond lengths of the phenolic hydroxyl group and the intramolecular hydrogen bond H···N in the ground state are 0.99 and 1.68 Å. To examine whether the thiolysis product HBTPP involves the PT process (via intramolecular hydrogen bond from O to N) or not, we further optimized the proton transferred structure and compared the energy between the enol form (here referred to as HBTPP(E)) and keto form (here referred to as HBTPP(K)). In HBTPP(K), the dihedral angle C1−C2− C3−C4 is 146.06° in the S0 state and the O···H bond is lengthened to 1.64 Å, while the H−N bond is shortened to 1.05 Å. Correspondingly, the optimized SCF energy increases from −1393.59 a.u. to −1393.04 a.u., implying that the enol form structure is a more stable configuration. Meanwhile, the S0 state potential energy curve of HBTPP as a function of the 5247

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A

Figure 3. Free energy profiles for the thiolysis reaction. R reactant, IM intermediate, TS transition state, P product. All free energies are in kcal mol−1; all bond lengths are in Å. The inset shows the ground state potential energy curve of the probe HBTPP-S and HS− as a function of the S···C distance.

Table 1. Comparison of Experimental and Calculated UV−Vis Absorption and Fluorescence Emission for the Probe HBTPP-S and Its Thiolysis Product HBTPP in Acetonitrile electronic transition HBTPP-S HBTPP

absorption emission absorption emission

S0 S1 S0 S1

→ → → →

S1 S0 S1 S0

energy (nm (eV)) 316 691 347 571

(3.91) (1.79) (3.57) (2.17)

fa

contribb

expt (nm)c

0.0175 0.0007 0.9940 1.2880

H→L+1 L→H−5 H→L L→H

320 380 620

a

Oscillator strength. bH, HOMO (highest occupied molecular orbital) and L, LUMO (lowest unoccupied molecular orbital). cThe experimental absorbance and emission band taken from ref29.

Figure 4. Experimental absorption/fluorescence spectra (a) (the data were taken from ref 29) and the calculated absorption/fluorescence spectra (b) of the thiolysis product HBTPP in acetonitrile.

(S0 → S1) is located at 489 nm, which is not in agreement with the experimental absorption spectra of the product (380 nm). Hence, we further make sure that the end product of the thiolysis reaction should be HBTPP.

3.2. The Frontier Molecular Orbitals of HBTPP-S and HBTPP(E) and the Two Indexes of DCT and ξCT. The frontier molecular orbitals of HBTPP-S and HBTPP(E) are shown in Figure 5. The first singlet transition for HBTPP-S is 5248

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A

Figure 5. Frontier molecular orbitals (MOs) of the probe HBTPP-S (left column) and its thiolysis product HBTPP (right column) in acetonitrile. E and K are the enol and keto form, respectively. The asterisk indicates the first excited state. d-PET stands for donor-excited photoinduced electron transfer, ESIPT stands for excited state intramolecular proton transfer, and TICT stands for twisted intramolecular charge transfer.

Figure 6. Divided part 1 and part 2 of the probe HBTPP-S and its thiolysis product HBTPP.

Table 2. Contribution of Part 1 to Frontier Molecular Orbitals for HBTPP-S and HBTPP Using the Ground State and the First Excited State Orbitals by C-Squared Population Analysisa first excited state orbitals

ground state orbitals HBTPP-S (%) HOMO LUMO+1

7.00 12.44

HBTPP(E) (%) HOMO LUMO

9.82 36.64

HBTPP(K*) (%) HOMO LUMO

14.74 30.49

a

E and K are the enol and keto form of HBTPP, respectively. The asterisk (∗) indicates the first excited state.

dominated by the HOMO → LUMO+1 transition, while for HBTPP(E) the first singlet transition corresponds to the transition of HOMO → LUMO. For HBTPP-S or HBTPP(E), the electron density of methylpyridine moiety increases after the transition from HOMO to LUMO+1 or to LUMO, suggesting that the S1 states of HBTPP-S and HBTPP(E) are both TICT states. To describe the TICT character more accurately, we further divided HBTPP-S and HBTPP(E) into two parts and calculated the contribution of part 1 to the HOMO and LUMO+1/LUMO using the ground state orbitals of HBTPP-S and HBTPP(E) (see Figure 6 and Table 2) by C-squared population analysis.89 For HBTPP-S, there is 7.00% of the electron density residing on the part 1 for HOMO and 12.44% for LUMO+1, while for HBTPP(E), 9.82% of the electron

Figure 7. Constructed two-dimensional potential energy surface of the thiolysis product HBTPP in the S1 state as functions of the methylpyridine twisting angle ranging from 140.88° to 175° and the O−H bond length ranging from 0.95 to1.9 Å, calculated at the TDDFT/ CAM-B3LYP/TZVP level. The white arrow indicates the direction of both ESIPT and twisting process in the S1 state. The structures corresponding to HBTPP(E) and HBTPP(K*) are also shown. The asterisk (∗) indicates the first excited state.

density resides on the part 1 for HOMO and 36.44% for LUMO. Therefore, the S1 state of HBTPP(E) shows more obvious TICT character than HBTPP-S, which may give rise to the red shift in the absorbance spectra that is observed in the experimental UV−vis absorbance.29 5249

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A Table 3. Comparison of the Calculated S0 and S1 States Geometric Parameters of HBTPP(E) and HBTPP(K*) in Gas Phase and in Acetonitrile Solution ground state

first excited state

molecular parameters

gas phase

acetonitrile

gas phase

acetonitrile

O−H bond (Å) C1−C2−C3−C4 (deg)

1.00 145.95

0.99 140.88

1.80 159.14

1.85 169.35

Figure 10. Values of the DCT index calculated for the first excited state on the S1 relaxed geometries along the O−H bond length of HBTPP at the TDDFT/CAM-B3LYP/TZVP level. The energies of the S0 state were calculated under the geometries of the corresponding S1 state.

depletion and increase regions produced upon excitation, and ξCT means the dihedral angle formed between the planes defined by the two barycenters of charges and two carbon atoms belonging to the ethylene bridge.75,90−93 We also calculated the DCT and ξCT indexes using the Multiwfn software,94 and we got the DCT value of 2.44 Å (0.47 Å) and the ξCT value of 17.75° (1.17°) for HBTPP(E) (HBTPP-S), which demonstrated more obvious TICT character in HBTPP(E), in accordance with the picture of the frontier molecule orbitals. 3.3. The First Excited State Structures, the S1-State Potential Energy Surface/Curves, and the Mechanism of Coupled ESIPT and TICT in HBTPP. The geometries of the S1 state of HBTPP-S and the enol form (HBTPP(E*), the asterisk indicates the first excited state) as well as the keto form (HBTPP(K*)) of HBTPP were also presented in Figure 1.

Figure 8. Potential energy curves of the thiolysis product HBTPP(K), as a function of the twisted dihedral angle (C1−C2−C3−C4) in the S0 and S1 states (g, in the gas phase; s, in acetonitrile solution), calculated at the TDDFT/CAM-B3LYP/TZVP level. The small circles mean the structures of HBTPP(K*) before and after the TICT process. The asterisk (∗) indicates the first excited state.

The DCT and ξCT indexes are also suitable for describing the degree of the TICT character (that is, the larger the DCT and ξCT values, the larger the TICT character). Here, DCT means the distance between the barycenter of charges of the density

Figure 9. Schematic illustration of d-PET from the excited HBTPP to 2,4-dinitrophenyl group in the unreacted probe HBTPP-S, but d-PET lacks in the keto form of the thiolysis product HBTPP. 5250

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A Scheme 2. Four-Level Model for the ESIPT Process of HBTPPa

a

E and K are the enol and keto forms, respectively. The asterisk indicates the first excited state.

In the first excited state, the dihedral angle C1−C2−C3−C4 of HBTPP-S is 134.57° and HBTPP-S has a similar structure to that in the ground state. In HBTPP(E*) (see Figure 1), the bond length of O−H is lengthened from 0.99 Å in the S0 state to 1.03 Å in the S1 state and the intramolecular hydrogen bond H···N is shortened from 1.68 Å in the S0 state to 1.56 Å in the S1 state. The strengthening of the S1 state intramolecular hydrogen bond implies that the ESIPT process would take place via intramolecular hydrogen bond from O to N. We also note that the energy of HBTPP(K*) is lower than that of HBTPP(E*), thus, HBTPP(K*) is a more stable structure than HBTPP(E*). For HBTPP, the dihedral angle C1−C2− C3−C4 changes from 140.88° in the S0 state (HBTPP(E)) to 169.35° in the S1 state (HBTPP(K*)). After excitation to the S1 state, HBTPP undergoes the ESIPT process and the torsional relaxation by twisting the methylpyridine group. The torsional angle (C1−C2−C3−C4) changes from 140.88° in the S0 state to 169.35° in the S1 state, suggesting that the product HBTPP is more coplanar in the S1 state than that in the S0 state. Meanwhile, the dipole moment for HBTPP decreases from 18.88 D in the S0 state to 11.29 D in the S1

state, and the larger dipole moment of HBTPP(E) may be owing to greater charge transfer character. Figure 2 shows the S1-state potential energy curve of HBTPP as a function of the O−H bond length (the S0-state curve of HBTPP is also presented in the figure). To obtain this curve, structure optimizations were done by fixing the O−H bond while relaxing other parameters. The S1-state curve exhibits a very low barrier of 0.65 kcal mol−1, after overcoming the barrier, the energy of S1 state decreases along the ESIPT coordinate. The energy of the keto form in the S1 state is lower than that of enol form, proving that the ESIPT process of HBTPP is exothermic and formation of the keto form is thermodynamically allowed. Consequently, the transfer of the proton H from O to N via the intramolecular hydrogen bond of HBTPP in the S1 state is energetically favorable and ESIPT process can proceed to form the stable keto form structure. To sum up, in the ground state, the ESIPT process is inhibited because phenolic hydroxyl group is masked by 2,4-dinitrophenyl group, however, after thiolysis of 2,4-dinitrophenyl ether by H2S and upon photoexcitation, the proton of phenolic hydroxyl group can transfer to the 5251

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A neighboring nitrogen group and HBTPP(E) is converted into the stable HBTPP(K*). On the other hand, after excitation, the product HBTPP may undergo TICT accompanied by the internal rotational motion of the methylpyridine group. To investigate such a mechanism of ESIPT and TICT in the S1 state, we constructed the two-dimensional S1-state potential energy surface (PES) of the thiolysis product HBTPP, as a function of methylpyridine twisting angle C1−C2−C3−C4 (ranging from 140.88° to 175°) and the distance between O and H (ranging from 0.95 to 1.9 Å), by going through a time-consuming process. It is necessary to explain that the bond length at 0.99 and 1.85 Å corresponds to the O−H bond length of HBTPP(E) and HBTPP(K*), and the angles of 140.88° and 169.35° correspond to the dihedral angle of HBTPP(E) and HBTPP(K*). The two-dimensional PES are shown in Figure 7, together with the two critical points of HBTPP(E) and HBTPP(K*). From the PES, it can be seen that HBTPP first undergoes the TICT process by rotating the methylpyridine moiety upon photoexcitation. Then, the concomitant action of ESIPT and methylpyridine rotational motion comes into existence when the twisting angle of the enol form is larger than ∼155° (the coupling of ESIPT and TICT processes appears in the region of the twisting angle ranging from 155° to 169.35°). Afterward, the energy drops markedly with the increase in the O−H bond and reaches a stable point at the O−H bond length of ∼1.85 Å. Here the ESIPT process is very likely to play the main part in the sensing process, due to only after undergoing the coupled ESIPT and TICT processes, the stable near-planar HBTPP(K*) can be formed eventually. To reveal the solvent effect on the torsional motion of HBTPP after ESIPT, Table 3 compares the geometric parameters of HBTPP(E) and HBTPP(K*) in gas phase and in the acetonitrile solution. In the gas phase, the torsional angle (C1−C2−C3−C4) changes from 145.95° (S0) to 159.14° (S1). Figure 8 compares the S1-state potential energy curves of HBTPP(K*) in the gas phase and in acetonitrile solution, as a function of the twisted dihedral angle C1−C2−C3−C4 and at the fixed O−H bond lengths of 1.80 and 1.85 Å, respectively. The energies of S0 state calculated under the geometries of the corresponding S1 state are also shown in the figure. The two barrierless excited-state potential energy curves might demonstrate that TICT occurs easily in both gas phase and acetonitrile solution. The whole picture of the potential curves shows that, after the ESIPT, there probably is no intersystem crossing or conical intersection (CI) between the S0 and S1 states of HBTPP, which might exclude the competitive nonradiative channel of CI for HBTPP. Additionally, the energy difference before and after the TICT process in the S1 state is rather small in both cases, suggesting that the rotational motion lowers little energy for HBTPP in the S1 state in both gas phase and acetonitrile solution. So the internal rotational motion of methylpyridine group seems to play the minor part in the sensing process of the thiolysis product HBTPP. 3.4. Fluorescence Spectra, d-PET Process, ESIPTCoupled TICT, and Sensing Mechanism. Table 1 shows the calculated S1 → S0 transition for HBTPP-S* and HBTPP(K*). It is found that for the probe HBTPP-S the oscillator strength (f = 0.0007) of this transition approaches zero, thus, the S1 state of HBTPP-S is a dark state which decays to the ground state via the nonradiative channels, causing the quenched fluorescence of HBTPP-S. On the other hand, after the thiolysis of 2,4-dinitrophenyl ether bond, the S1 → S0

vertical excitation energy for HBTPP(K*) is predicted at 571 nm with the oscillator strength of 1.2880. Here, the large Stokes shift (224 nm) could be ascribed to the ESIPT process. All these theoretical results agree well with the previous experimental observations that the probe HBTPP-S is nonfluorescence and the thiolysis product HBTPP shows remarkably red keto emission (620 nm) and large Stokes shift (240 nm)29 (see Figure 4). On the basis of the molecular orbitals in Figure 5, for HBTPP-S, the π electrons on the LUMO+1 are primarily resided on the HBTPP group while the electrons on LUMO is delocalized over the electron-withdrawing 2,4-dinitrophenyl moiety. Therefore, upon photoexcitation, the nonradiative donor-excited photoinduced electron transfer (d-PET) can take place from the excited HBTPP group to the electronwithdrawing 2,4-dinitrophenyl group. Because of such efficient d-PET process (see Figure 9) and the notorious fluorescence quenching effect of 2,4-dinitrophenyl group,23−25,95,96 the probe HBTPP-S is essentially nonfluorescent. In addition to the d-PET process, other decay paths such as the nonradiative decay via nonadiabatic CI could also contribute to the fluorescence quenching of HBTPP-S. Because of the present TDDFT method is unable to locate the CI in HBTPP-S, the search of these other decay paths is temporally beyond the scope of this study. Thiolysis of the 2,4-dinitrophenyl ether bond by H2S has converted the nonfluorescent HBTPP-S into the fluorescent HBTPP, which lacks the d-PET process (see Figure 9) and undergoes the ESIPT coupled the TICT at the first excited state leading to the overlap of electrons on the transition orbitals of HBTPP(K*) (Figure 5). Further evidence can be seen from the charge transfer characters of HBTPP in Table 2. The values in this table were calculated by C-squared population analysis using the ground state orbital of HBTPP(E) and the first excited state orbital of HBTPP(K). Seen from the table, there is 9.82% of the electron density residing on part 1 for HOMO and 36.64% for LUMO using the ground state orbitals in the calculation, while there is 14.74% of the electron density residing on the part 1 for HOMO and 30.49% for LUMO using the first excited state orbitals in the calculation, demonstrating the ESIPT process (coupled to the TICT process) leads to the more overlap of electrons on the transition orbitals. The charge redistribution upon photoexcitation gives rise to that HBTPP(K*) possesses relatively weak ICT character accounting also for its fluorescent nature. Figure 10 shows our calculated DCT index values along the ESIPT coordinate of the O−H distance where the highest DCT value is found at the O−H distance of 1.03 Å corresponding to the structure of HBTPP(E*). The value decreases to the end of the O−H bond length at 1.85 Å corresponding to the structure of HBTPP(K*). Our analysis of the DCT index further confirmed the weak ICT character of HBTPP(K*). Some of these calculated results agree well with the previous studies about the fluorescent turn-on probe based on the H2S-promoted thiolysis of dinitrophenyl ether.23−25,95,96 According to previous literature work,97−99 ESIPT can be viewed as a very fast photo tautomerism process taking place along an intramolecular hydrogen bond, based on a four-level model (E → E* → K* → K → E). Scheme 2 presents the four-level model deduced from our calculation results for the ESIPT process in HBTPP. As illustrated, the energy level of HBTPP(K*) is lower than that of HBTPP(E*) while the energy level of HBTPP(K) is higher than that of HBTPP(E). 5252

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

The Journal of Physical Chemistry A



Therefore, the energy gap between HBTPP(K*) and HBTPP(K) is smaller than that between HBTPP(E*) and HBTPP(E), demonstrating that HBTPP(K*) emits long-wavelength fluorescence. Therefore, the large Stokes shift observed in HBTPP is due to the ESIPT after the UV-excitation. Overall, our theoretical results confirmed that the nonradiative d-PET process can lead to the fluorescence quenching of the probe HBTPP-S, while thiolysis of the ether bond by H2S, the lack of d-PET process, and the coupled ESIPT and TICT processes are all responsible for the strong fluorescence emission and large Stokes shift of the thiolysis product HBTPP. As a results, an intensive fluorescence emitted by HBTPP(K*) can be easily observed by the naked eye and the probe HBTPP-S can serve as an efficient fluorescent turn-on probe for the detection of the hydrogen sulfide.

REFERENCES

(1) Wang, R. Two’s Company, Three’s a Crowd: Can H2S be the Third Endogenous Gaseous Transmitter? FASEB J. 2002, 16, 1792− 1798. (2) Wang, R. The Gasotransmitter Role of Hydrogen Sulfide. Antioxid. Redox Signaling 2003, 5, 493−501. (3) Yang, G. D.; Sun, X. F.; Wang, R. Hydrogen Sulfide-Induced Apoptosis of Human Aorta Smooth Muscle Cells via the Activation of Mitogen-Activated Protein Kinases and Caspase-3. FASEB J. 2004, 18, 1782−1784. (4) Gadalla, M. M.; Snyder, S. H. Hydrogen Sulfide as a Gasotransmitter. J. Neurochem. 2010, 113, 14−26. (5) Yang, G.; Wu, L.; Jiang, B.; Yang, W.; Qi, J.; Cao, K.; Meng, Q.; Mustafa, A. K.; Mu, W.; Zhang, S.; et al. H2S as a Physiologic Vasorelaxant: Hypertension in Mice with Deletion of Cystathionine Gamma-Lyase. Science 2008, 322, 587−590. (6) Abe, K.; Kimura, H. The Possible Role of Hydrogen Sulfide as an Endogenous Neuromodulator. J. Neurosci. 1996, 16, 1066−1071. (7) Hosoki, R.; Matsuki, N.; Kimura, H. The Possible Role of Hydrogen Sulfide as an Endogenous Smooth Muscle Relaxant in Synergy with Nitric Oxide. Biochem. Biophys. Res. Commun. 1997, 237, 527−531. (8) Eto, K.; Asada, T.; Arima, K.; Makifuchi, T.; Kimura, H. Brain Hydrogen Sulfide is Severely Decreased in Alzheimer’s Disease. Biochem. Biophys. Res. Commun. 2002, 293, 1485−1488. (9) Kamoun, P.; Belardinelli, M. C.; Chabli, A.; Lallouchi, K.; Chadefaux-Vekemans, B. Endogenous Hydrogen Sulfide Overproduction in Down Syndrome. Am. J. Med. Genet. A 2003, 116a, 310−311. (10) Yang, W.; Yang, G. D.; Jia, X. M.; Wu, L. Y.; Wang, R. Activation of K-ATP Channels by H2S in Rat Insulin-Secreting Cells and the Underlying Mechanisms. J. Physiol. 2005, 569, 519−531. (11) Nagy, P.; Winterbourn, C. C. Reaction of Hydrogen Sulfide with the Neutrophil Oxidant Hypochlorous Acid to Generate Polysulfides. Chem. Res. Toxicol. 2010, 23, 1541−1543. (12) Sasakura, K.; Hanaoka, K.; Shibuya, N.; Mikami, Y.; Kimura, Y.; Komatsu, T.; Ueno, T.; Terai, T.; Kimura, H.; Nagano, T. Development of a Highly Selective Fluorescence Probe for Hydrogen Sulfide. J. Am. Chem. Soc. 2011, 133, 18003−18005. (13) Wang, J.; Long, L.; Xie, D.; Zhan, Y. Highly Selective Fluorescence Turn-on Sensor for Hydrogen Sulfide and Imaging in Living Cells. J. Lumin. 2013, 139, 40−46. (14) Hou, F.; Huang, L.; Xi, P.; Cheng, J.; Zhao, X.; Xie, G.; Shi, Y.; Cheng, F.; Yao, X.; Bai, D.; et al. A Retrievable and Highly Selective Fluorescent Probe for Monitoring Sulfide and Imaging in Living Cells. Inorg. Chem. 2012, 51, 2454−2460. (15) Lippert, A. R.; New, E. J.; Chang, C. J. Reaction-Based Fluorescent Probes for Selective Imaging of Hydrogen Sulfide in Living Cells. J. Am. Chem. Soc. 2011, 133, 10078−10080. (16) Peng, H.; Cheng, Y.; Dai, C.; King, A. L.; Predmore, B. L.; Lefer, D. J.; Wang, B. A Fluorescent Probe for Fast and Quantitative Detection of Hydrogen Sulfide in Blood. Angew. Chem., Int. Ed. 2011, 50, 9672−9675. (17) Wu, Z.; Li, Z.; Yang, L.; Han, J.; Han, S. Fluorogenic Detection of Hydrogen Sulfide via Reductive Unmasking of O-Azidomethylbenzoyl-Coumarin Conjugate. Chem. Commun. 2012, 48, 10120−10122. (18) Wang, R.; Yu, F.; Chen, L.; Chen, H.; Wang, L.; Zhang, W. A Highly Selective Turn-on Near-Infrared Fluorescent Probe for Hydrogen Sulfide Detection and Imaging in Living Cells. Chem. Commun. 2012, 48, 11757−11759. (19) Wu, M. Y.; Li, K.; Hou, J. T.; Huang, Z.; Yu, X. Q. A Selective Colorimetric and Ratiometric Fluorescent Probe for Hydrogen Sulfide. Org. Biomol. Chem. 2012, 10, 8342−8347. (20) Li, X.; Zhang, S.; Cao, J.; Xie, N.; Liu, T.; Yang, B.; He, Q.; Hu, Y. An ICT-Based Fluorescent Switch-on Probe for Hydrogen Sulfide in Living Cells. Chem. Commun. 2013, 49, 8656−8658. (21) Chen, Y.; Zhu, C.; Yang, Z.; Chen, J.; He, Y.; Jiao, Y.; He, W.; Qiu, L.; Cen, J.; Guo, Z. A Ratiometric Fluorescent Probe for Rapid Detection of Hydrogen Sulfide in Mitochondria. Angew. Chem., Int. Ed. 2013, 52, 1688−1691.

4. CONCLUSIONS DFT and TDDFT methods have been applied to elucidate the sensing mechanism of the hydrogen sulfide probe HBTPP-S based on the thiolysis of 2,4-dinitrophenyl ether bond. The free energy profile displays that the thiolysis reaction has a moderate barrier (18.40 kcal mol−1) that can lead to a favorable response speed of the probe HBTPP-S for HS−. The red shift in the absorbance spectra, the strong fluorescence emission and the large Stokes shift, all coincide with the experimental observations. The 2,4-dinitrophenyl group in the probe HBTPP-S efficiently quenches the fluorescence of the HBTPP fragment through the d-PET process, and it also masks the phenolic hydroxyl group to inhibit the PT and ESIPT processes. Accordingly, the probe HBTPP-S is essentially nonfluorescent. After removing the masking 2,4-dinitrophenyl group through the H2S-promoted thiolysis, the twisted enol form of thiolysis product HBTPP(E) is formed. Upon photoexcitation, HBTPP(E*), first undergoes the TICT process accompanied by the torsional relaxation, and meanwhile, the phenol proton of HBTPP(E*) is transferred to the neighboring nitrogen group through the ESIPT process. In the end, the HBTPP(E*) is converted into the stable near-planar keto form HBTPP(K*) which can emit strong fluorescent and show large Stokes shift. Here, the H2S-promoted thiolysis, the absence of d-PET process and the ESIPT-coupled TICT processes all account for the dramatic change in fluorescent properties of HBTPP such as strong fluorescent emission and large Stokes shift. In conclusion, the thiolysis of the ether bond, the nonradiative d-PET process and the coupling of ESIPT with TICT processes play the decisive roles during the sensing processes of the hydrogensulfide probe HBTPP-S.



Article

AUTHOR INFORMATION

Corresponding Author

*(T.-S.C.) E-mail: [email protected]; [email protected]. ORCID

Tian-Shu Chu: 0000-0002-3519-8737 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work is supported by the National Natural Science Foundation of China (Grant No. 21273234), the Shandong Provincial Natural Science Foundation of China (Grant No. ZR2014AM025), and the National Basic Research Program of China (Grant No. 2013CB834604). 5253

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A (22) Feng, X.; Zhang, T.; Liu, J. T.; Miao, J. Y.; Zhao, B. X. A New Ratiometric Fluorescent Probe for Rapid, Sensitive and Selective Detection of Endogenous Hydrogen Sulfide in Mitochondria. Chem. Commun. 2016, 52, 3131−3134. (23) Cao, X.; Lin, W.; Zheng, K.; He, L. A Near-Infrared Fluorescent Turn-on Probe for Fluorescence Imaging of Hydrogen Sulfide in Living Cells Based on Thiolysis of Dinitrophenyl Ether. Chem. Commun. 2012, 48, 10529−10531. (24) Zheng, K.; Lin, W.; Tan, L.; Cheng, D. A Two-Photon Fluorescent Probe with a Large Turn-on Signal for Imaging Hydrogen Sulfide in Living Tissues. Anal. Chim. Acta 2015, 853, 548−554. (25) Yuan, L.; Zuo, Q. P. Reaction-Based Fluorescent Probe for Hydrogen Sulfide with Large Signal-to-Noise Ratio in Living Cells and Tissues. Sens. Actuators, B 2014, 196, 151−155. (26) Wang, B. H.; Li, P.; Yu, F. B.; Chen, J. S.; Qu, Z. J.; Han, K. L. A Near-Infrared Reversible and Ratiometric Fluorescent Probe Based on Se-BODIPY for the Redox Cycle Mediated by Hypobromous Acid and Hydrogen Sulfide in Living Cells. Chem. Commun. 2013, 49, 5790− 5792. (27) Wang, B. H.; Li, P.; Yu, F. B.; Song, P.; Sun, X. F.; Yang, S. Q.; Lou, Z. R.; Han, K. L. A Reversible Fluorescence Probe Based on SeBODIPY for the Redox Cycle Between HClO Oxidative Stress and H2S Repair in Living Cells. Chem. Commun. 2013, 49, 1014−1016. (28) Yu, F. B.; Li, P.; Song, P.; Wang, B. H.; Zhao, J. Z.; Han, K. L. An ICT-Based Strategy to a Colorimetric and Ratiometric Fluorescence Probe for Hydrogen Sulfide in Living Cells. Chem. Commun. 2012, 48, 2852−2854. (29) Xu, P. F.; Liu, M. H.; Gao, T.; Zhang, H. L.; Li, Z. W.; Huang, X. Y.; Zeng, W. B. An ESIPT-Based Highly Selective and Sensitive Probe for the Detection of Hydrogen Sulfide. Tetrahedron Lett. 2015, 56, 4007−4010. (30) Yu, F. B.; Li, P.; Wang, B. H.; Han, K. L. Reversible NearInfrared Fluorescent Probe Introducing Tellurium to Mimetic Glutathione Peroxidase for Monitoring the Redox Cycles Between Peroxynitrite and Glutathione in Vivo. J. Am. Chem. Soc. 2013, 135, 7674−7680. (31) Lou, Z. R.; Li, P.; Han, K. L. Redox-Responsive Fluorescent Probes with Different Design Strategies. Acc. Chem. Res. 2015, 48, 1358−1368. (32) Hsieh, C. C.; Cheng, Y. M.; Hsu, C. J.; Chen, K. Y.; Chou, P. T. Spectroscopy and Femtosecond Dynamics of Excited-State Proton Transfer Induced Charge Transfer Reaction. J. Phys. Chem. A 2008, 112, 8323−8332. (33) Gutierrez, M.; Alarcos, N.; Liras, M.; Sanchez, F.; Douhal, A. Switching to a Reversible Proton Motion in a Charge-Transferred Dye. J. Phys. Chem. B 2015, 119, 552−562. (34) Alarcos, N.; Gutierrez, M.; Liras, M.; Sanchez, F.; Douhal, A. An Abnormally Slow Proton Transfer Reaction in a Simple HBO Derivative due to Ultrafast Intramolecular-Charge Transfer Events. Phys. Chem. Chem. Phys. 2015, 17, 16257−16269. (35) Alarcos, N.; Gutierrez, M.; Liras, M.; Sanchez, F.; Douhal, A. From Intra- to Inter-Molecular Hydrogen Bonds with the Surroundings: Steady-State and Time-Resolved Behaviours. Photoch. Photob. Sci. 2015, 14, 1306−1318. (36) Abeywickrama, C. S.; Pang, Y. Synthesis of Fused 2-(2′Hydroxyphenyl)Benzoxazole Derivatives: the Impact of Meta-/ParaSubstitution on Fluorescence and Zinc Binding. Tetrahedron Lett. 2016, 57, 3518−3522. (37) Honda, T.; Ishida, Y.; Arai, T. Substituent Effect of Formyl Group on Photochemical Properties of 2-(2′-Hydroxyphenyl)Benzazole Derivative. Bull. Chem. Soc. Jpn. 2016, 89, 876−878. (38) Li, C. Z.; Yang, Y. G.; Ma, C.; Liu, Y. F. Effect of Amino Group on the Excited-State Intramolecular Proton Transfer (ESIPT) Mechanisms of 2-(2′-Hydroxyphenyl)Benzoxazole and its Amino Derivatives. RSC Adv. 2016, 6, 5134−5140. (39) Le Guennic, B.; Chibani, S.; Charaf-Eddin, A.; Massue, J.; Ziessel, R.; Ulrich, G.; Jacquemin, D. The NBO Pattern in Luminescent Chromophores: Unravelling Excited-State Features Using TD-DFT. Phys. Chem. Chem. Phys. 2013, 15, 7534−7540.

(40) Houari, Y.; Charaf-Eddin, A.; Laurent, A. D.; Massue, J.; Ziessel, R.; Ulrich, G.; Jacquemin, D. Modeling Optical Signatures and Excited-State Reactivities of Substituted Hydroxyphenylbenzoxazole (HBO) ESIPT Dyes. Phys. Chem. Chem. Phys. 2014, 16, 1319−1321. (41) Aly, S. M.; Usman, A.; Alzayer, M.; Hamdi, G. A.; Alarousu, E.; Mohammed, O. F. Solvent-Dependent Excited-State Hydrogen Transfer and Intersystem Crossing in 2-(2′-Hydroxyphenyl)-Benzothiazole. J. Phys. Chem. B 2015, 119, 2596−2603. (42) Geng, L. H.; Yang, X. F.; Zhong, Y. G.; Li, Z.; Li, H. ″QuinonePhenol″ Transduction Activated Excited-State Intramolecular Proton Transfer: a New Strategy toward Ratiometric Fluorescent Probe for Sulfite in Living Cells. Dyes Pigm. 2015, 120, 213−219. (43) Raoui, M.; Massue, J.; Azarias, C.; Jacquemin, D.; Ulrich, G. Highly Fluorescent Extended 2-(2′-Hydroxyphenyl)Benzazole Dyes: Synthesis, Optical Properties and First-Principle Calculations. Chem. Commun. 2016, 52, 9216−9219. (44) Cheng, J. L.; Liu, D.; Li, W.; Bao, L. J.; Han, K. L. Comprehensive Studies on Excited-State Proton Transfer of a Series of 2-(2 ′-Hydroxyphenyl)Benzothiazole Derivatives: Synthesis, Optical Properties, and Theoretical Calculations. J. Phys. Chem. C 2015, 119, 4242−4251. (45) Sahana, S.; Mishra, G.; Sivakumar, S.; Bharadwaj, P. K. A 2-(2′Hydroxyphenyl) Benzothiazole (HBT)-Quinoline Conjugate: a Highly Specific Fluorescent Probe For Hg2+ Based on ESIPT and its Application in Bioimaging. Dalton. T 2015, 44, 20139−20146. (46) Majumdar, P.; Zhao, J. Z. 2-(2-Hydroxyphenyi)-Benzothiazole (HBT)-Rhodamine Dyad: Acidswitchable Absorption and Fluorescence of Excited-State Intramolecular Proton Transfer (ESIPT). J. Phys. Chem. B 2015, 119, 2384−2394. (47) Barman, S.; Mukhopadhyay, S. K.; Biswas, S.; Nandi, S.; Gangopadhyay, M.; Dey, S.; Anoop, A.; Pradeep Singh, N. D. A pHydroxyphenacyl-Benzothiazole-Chlorambucil Conjugate as a RealTime-Monitoring Drug-Delivery System Assisted by Excited-State Intramolecular Proton Transfer. Angew. Chem., Int. Ed. 2016, 55, 4194−4198. (48) Purkayastha, P.; Chattopadhyay, N. Role of Rotamerisation and Excited State Intramolecular Proton Transfer in the Photophysics of 2(2′-hydroxyphenyl)benzoxazole, 2-(2′-hydroxyphenyl)benzimidazole and 2-(2′-hydroxyphenyl)benzothiazole: a Theoretical Study. Phys. Chem. Chem. Phys. 2000, 2, 203−210. (49) Barbatti, M.; Aquino, A. J. A.; Lischka, H.; Schriever, C.; Lochbrunner, S.; Riedle, E. Ultrafast Internal Cconversion Pathway and Mechanism in 2-(2′-hydroxyphenyl)benzothiazole: a Case Study for Excited-State Intramolecular Proton Transfer Systems. Phys. Chem. Chem. Phys. 2009, 11, 1406−1415. (50) Seo, J.; Kim, S.; Park, S. Y. Strong Solvatochromic Fluorescence from the Intramolecular Charge-Transfer State Created by ExcitedState Intramolecular Proton Transfer. J. Am. Chem. Soc. 2004, 126, 11154−11155. (51) Kim, C. H.; Park, J.; Seo, J.; Park, S. Y.; Joo, T. Excited State Intramolecular Proton Transfer and Charge Transfer Dynamics of a 2(2′-Hydroxyphenyl)Benzoxazole Derivative in Solution. J. Phys. Chem. A 2010, 114, 5618−5629. (52) Cheng, Y. M.; Pu, S. C.; Hsu, C. J.; Lai, C. H.; Chou, P. T. Femtosecond Dynamics on 2-(2′-Hydroxy-4′-Diethylaminophenyl)Benzothiazole: Solvent Polarity in the Excited-State Proton Transfer. ChemPhysChem 2006, 7, 1372−1381. (53) Zhao, G. J.; Han, K. L. Role of Intramolecular and Intermolecular Hydrogen Bonding in Both Singlet and Triplet Excited States of Aminofluorenones on Internal Conversion, Intersystem Crossing, and Twisted Intramolecular Charge Transfer. J. Phys. Chem. A 2009, 113, 14329−14335. (54) Chai, S.; Cong, S. L. Intermolecular Hydrogen Bond and Twisted Intramolecular Charge Transfer in Excited State of Fast Violet B (FVB) in Methanol Solution. Spectrochim. Acta, Part A 2014, 125, 67−72. (55) Grabowski, Z. R.; Rotkiewicz, K.; Rettig, W. Structural Changes Accompanying Intramolecular Electron Transfer: Focus on Twisted 5254

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A Intramolecular Charge-Transfer States and Structures. Chem. Rev. 2003, 103, 3899−4031. (56) Ghosh, R.; Nandi, A.; Palit, D. K. Solvent Sensitive Intramolecular Charge Transfer Dynamics in the Excited States of 4N,N-Dimethylamino-4′-Nitrobiphenyl. Phys. Chem. Chem. Phys. 2016, 18, 7661−7671. (57) Parr, R. G.; Yang, W. Density Functional Theory of Atoms and Molecules; Oxford University Press: Oxford, U.K., 1989. (58) Runge, E.; Gross, E. K. U. Density-Functional Theory for TimeDependent Systems. Phys. Rev. Lett. 1984, 52, 997−1000. (59) 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. (60) Li, G. Y.; Zhao, G. J.; Han, K. L.; He, G. Z. A TD-DFT Study on the Cyanide-Chemosensing Mechanism of 8-Formyl-7-Hydroxycoumarin. J. Comput. Chem. 2011, 32, 668−674. (61) 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. (62) Wang, D. D.; Lü, R.; Yuan, M. H.; Chen, J. S.; Feng, L. Q.; Fu, A. P.; Tian, F. H.; Varandas, A. J. C.; Chu, T. S. Photoinduced Coupled Twisted Intramolecular Charge Transfer and Excited-State Proton Transfer via Intermolecular Hydrogen Bonding: a DFT/TDDFT Study. Chem. Phys. Lett. 2014, 610-611, 179−185. (63) Chen, J. S.; Zhou, P. W.; Zhao, L.; Chu, T. S. A DFT/TDDFT Study of the Excited State Intramolecular Proton Transfer Based Sensing Mechanism for the Aqueous Fluoride Chemosensor BTTPB. RSC Adv. 2014, 4, 254−259. (64) Zhao, G. J.; Han, K. L. Hydrogen Bonding in the Electronic Excited State. Acc. Chem. Res. 2012, 45, 404−413. (65) 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. (66) Zhao, G. J.; Han, K. L. Site-Specific Solvation of the Photoexcited Protochlorophyllide a in Methanol: Formation of the Hydrogen-Bonded Intermediate State Induced by Hydrogen-Bond Strengthening. Biophys. J. 2008, 94, 38−46. (67) Chen, J. S.; Zhou, P. W.; Li, G. Y.; Chu, T. S.; He, G. Z. Fluoride Anion Sensing Mechanism of 2-Ureido-4[1H]-Pyrimidinone Quadruple Hydrogen-Bonded Supramolecular Assembly: Photoinduced Electron Transfer and Partial Configuration Change. J. Phys. Chem. B 2013, 117, 5212−5221. (68) Chen, J. S.; Zhou, P. W.; Yang, S. Q.; Fu, A. P.; Chu, T. S. Sensing Mechanism for a Fluoride Chemosensor: Invalidity of ExcitedState Proton Transfer Mechanism. Phys. Chem. Chem. Phys. 2013, 15, 16183−16189. (69) Chen, J. S.; Liu, R. Z.; Yang, Y.; Chu, T. S. Intramolecular Charge Transfer and Sensing Mechanism for a Colorimetric Fluoride Sensor Based on 1,8-Naphthalimide Derivatives. Theor. Chem. Acc. 2014, 133, 1−7. (70) Chen, J. S.; Yuan, M. H.; Wang, J. P.; Yang, Y.; Chu, T. S. Sensing Mechanism for Biothiols Chemosensor DCO: Roles of Excited-State Hydrogen-Bonding and Intramolecular Charge Transfer. J. Phys. Chem. A 2014, 118, 8986−8995. (71) Zhao, J. F.; Yang, Y. Excited State Proton Transfer Coupled with Twisted Intermolecular Charge Transfer for N,N-Dimethylanilino-1,3Diketone in High Polar Acetonitrile Solvent. J. Mol. Liq. 2016, 220, 735−741. (72) Chu, T. S.; Xu, J. M. Photoinduced Hydrogen-Bonding Dynamics. J. Mol. Model. 2016, 22, 200−219. (73) 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. (74) 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. (75) Savarese, M.; Brémond, É.; Adamo, C.; Rega, N.; Ciofini, I. Excited-State Proton Transfer and Intramolecular Charge Transfer in 1,3-Diketone Molecules. ChemPhysChem 2016, 17, 1530−1538. (76) 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 A.1; Gaussian, Inc.: Wallingford, CT, 2009. (77) Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (78) 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: Condens. Matter Mater. Phys. 1988, 37, 785−789. (79) Miehlich, B.; Savin, A.; Stoll, H.; Preuss, H. Results Obtained with the Correlation Energy Density Functionals of Becke and Lee, Yang and Parr. Chem. Phys. Lett. 1989, 157, 200−206. (80) Curtiss, L. A.; Redfern, P. C.; Raghavachari, K.; Pople, J. A. Assessment of Gaussian-2 and Density Functional Theories for the Computation of Ionization Potentials and Electron Affinities. J. Chem. Phys. 1998, 109, 42−55. (81) Perdew, J. P.; Ernzerhof, M.; Burke, K. Rationale for Mixing Exact Exchange with Density Functional Approximations. J. Chem. Phys. 1996, 105, 9982−9985. (82) Lynch, B. J.; Fast, P. L.; Harris, M.; Truhlar, D. G. Adiabatic Connection for Kinetics. J. Phys. Chem. A 2000, 104, 4811−4815. (83) Yanai, T.; Tew, D. P.; Handy, N. C. A New Hybrid ExchangeCorrelation Functional Using the Coulomb-Attenuating Method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51−57. (84) Schafer, A.; Horn, H.; Ahlrichs, R. Fully Optimized Contracted Gaussian-Basis Sets for Atoms Li to Kr. J. Chem. Phys. 1992, 97, 2571− 2577. (85) Schafer, 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. (86) Cancès, E.; Mennucci, B.; Tomasi, J. A New Integral Equation Formalism for the Polarizable Continuum Model: Theoretical Background and Applications to Isotropic and Anisotropic Dielectrics. J. Chem. Phys. 1997, 107, 3032−3041. (87) Mennucci, B.; Cances, E.; Tomasi, J. Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation Method: Theoretical Bases, Computational Implementation, and Numerical Applications. J. Phys. Chem. B 1997, 101, 10506−10517. (88) Cammi, R.; Tomasi, J. Remarks on the Use of the Apparent Surface-Charges (ASC) Methods in Solvation Problems-Iterative Versus Matrix-Inversion Procedures and the Renormalization of the Apparent Charges. J. Comput. Chem. 1995, 16, 1449−1458. (89) Ros, P.; Schuit, G. C. A. Molecular Orbital Calculations on Copper Chloride Complexes. Theor. Chim. Acta 1966, 4, 1−12. (90) Le Bahers, T.; Adamo, C.; Ciofini, I. A Qualitative Index of Spatial Extent in Charge Transfer Excitations. J. Chem. Theory Comput. 2011, 7, 2498−2506. (91) Savarese, M.; Netti, P. A.; Adamo, C.; Rega, N.; Ciofini, I. Exploring the Metric of Excited State Proton Transfer Reactions. J. Phys. Chem. B 2013, 117, 16165−16173. (92) Savarese, M.; Raucci, U.; Adamo, C.; Netti, P. A.; Ciofini, I.; Rega, N. Non-radiative Decay Paths in Rhodamines: New Theoretical Insights. Phys. Chem. Chem. Phys. 2014, 16, 20681−20688. (93) Savarese, M.; Raucci, U.; Netti, P. A.; Adamo, C.; Rega, N.; Ciofini, I. A Qualitative Model to Identify Non-radiative Decay Channels: the Spiropyran as Case Study. Theor. Chem. Acc. 2016, 135, 211−218. (94) Lu, T.; Chen, F. W. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem. 2012, 33, 580−592. (95) Wang, J. L.; Lin, W. Y.; Li, W. L. Three-Channel Fluorescent Sensing via Organic White Light-Emitting Dyes for Detection of Hydrogen Sulfide in Living Cells. Biomaterials 2013, 34, 7429−7436. 5255

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256

Article

The Journal of Physical Chemistry A (96) Liu, Y.; Feng, G. Q. A Visible Light Excitable Colorimetric and Fluorescent ESIPT Probe for Rapid and Selective Detection of Hydrogen Sulfide. Org. Biomol. Chem. 2014, 12, 438−445. (97) Mutai, T.; Ohkawa, T.; Shono, H.; Araki, K. The Development of Aryl-Substituted 2-phenylimidazo[1,2-a]pyridines (PIP) with Various Colors of Excited-State Intramolecular Proton Transfer (ESIPT) Luminescence in the Solid State. J. Mater. Chem. C 2016, 4, 3599−3606. (98) 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. Photochem. Photobiol., A 2016, 317, 175−185. (99) Ali, R.; Saleh, S. M.; Elshaarawy, R. F. M. Turn-on pH NanoFluorosensor Based on Imidazolium Salicylaldehyde Ionic LiquidLabeled Silica Nanoparticles. RSC Adv. 2016, 6, 86965−86975.

5256

DOI: 10.1021/acs.jpca.7b02606 J. Phys. Chem. A 2017, 121, 5245−5256