Theoretical Investigation on the “ON-OFF” Mechanism of a Fluorescent Probe for Thiophenols: Photoinduced Electron Transfer and Intramolecular Charge Transfer

In this study, the sensing mechanism of (2E,4E)-5-(4-(dimethylamino)phenyl)-1-(2-(2,4dinitrophenoxy)phenyl)penta-2,4-dien-1-one (DAPH-DNP) towards thiophenols was investigated by density functional theory (DFT) and time-dependent DFT (TD-DFT). The DNP group plays an important role in charge transfer excitation. Due to the typical donor-excited photo-induced electron transfer (d-PET) process, DAPH-DNP has fluorescence quenching behavior. After the thiolysis reaction between DAPH-DNP and thiophenol, the hydroxyl group is released, and DAPH is generated with the reaction showing strong fluorescence. The fluorescence enhancement of DAPH is not caused by an excited-state intramolecular proton transfer (ESIPT) process. The potential energy curves (PECs) show that DAPH-keto is less stable than DAPH-enol. The frontier molecular orbitals (FMOs) of DAPH show that the excitation process is accompanied by intramolecular charger transfer (ICT), and the corresponding character of DAPH was further confirmed by hole-electron and interfragment charge transfer (IFCT) analysis methods. Above all, the sensing mechanism of the turn-on type probe DAPH-DNP towards thiophenol is based on the PET mechanism.


Introduction
Thiophenols (PhSH) are colorless liquids with a foul odor.Due to the characteristics of the S-H groups on the benzene rings [1,2], they are often utilized as intermediates in the manufacture of fine chemicals such as colors, pharmaceuticals, insecticides, polymer inhibitors, and antioxidants.Thiophenols are highly toxic and corrosive, when inhaled or ingested by the human body, they can lead to larynx, bronchospasm, edema, chemical pneumonia, and even death.The Environmental Protection Agency has identified thiophenols as priority contaminants.In recent years, the detection of thiophenols in biological and environmental samples has aroused great interest among scholars [3][4][5][6][7].
Traditional methods for detecting thiophenols include inductively coupled plasma mass spectrometry (ICP-MS), electrochemistry, and high-performance liquid chromatography (HPLC) [8,9].Due to limitations such as poor coordination and low sensitivity, it is difficult to detect appropriately the concentration of thiophenols in vitro and in vivo.Recently, near-infrared (NIR) fluorescent probes for the detection of thiophenols in water samples and living cells have been reported, which have simple synthesis routes, high selectivity, and high sensitivity [10][11][12][13].Fluorescent probes can detect trace components by observable wavelength changes and fluorescence intensities.Compared to other methods, fluorescent probes can penetrate deep into the tissue, reducing the interference of the autofluorescence of samples, and thus realizing the detection of biological samples.Fluorescent probes other methods, fluorescent probes can penetrate deep into the tissue, reducing the interference of the autofluorescence of samples, and thus realizing the detection of biological samples.Fluorescent probes have been synthesized primarily by relying on a diverse array of sensing mechanisms, including excited-state intramolecular proton transfer (ESIPT) [14][15][16], chelation-enhanced fluorescence (CHEF) [17][18][19], photoinduced electron transfer (PET) [20][21][22], intramolecular charge transfer (ICT) [23,24], and aggregation-induced enhanced emission (AIEE) [25,26].ESIPT based fluorescent probes exhibit unique excited-state photophysical properties of high fluorescence quantum yield and large Stokes shift.Fluorescence emission peaks of ESIPT molecules are often red-shifted by tens to hundreds of nanometers compared to their absorption peaks, and their fluorescence spectra have narrow peak widths, resulting in high optical resolution and detection sensitivity.ESIPT is a completely reversible cyclic process.Upon photoexcitation, the H atom of the hydrogen bond H-X...Y is transferred from the electronegative atom X to another electronegative atom Y, accompanied by the generation of a new hydrogen bond (H...X-Y), which is a typical enol-keto tautomerization.As a result, they have broad applications in fields such as biological imaging and drug screening [27][28][29].
Recently, a simple NIR-emitting fluorescent probe (DAPH-DNP), consisting of a thiophenol recognition unit (2,4-dinitrophenyl (DNP) group) and a fluorophore (2E,4E)5-(4-(dimethylamino)phenyl)-1-(2-hydroxyphenyl)penta-2,4dien-1-one (DAPH), with large Stokes shift based on the ESIPT mechanism was presented by Li et al. [30], which could be applied to detect effectively thiophenol in water samples and living cells with a remarkable recovery rate (see Scheme 1).Because the DNP group is highly selective to thiophene, the fluorescence of the probe DAPH-DNP is quenched, which hinders the ESIPT process.The nucleophilic substitution process of dinitrophenyl ethers involved in thiophenol release the hydroxyl group and thus restore the ESIPT process, resulting in fluorescence emission (Scheme Ⅰ).However, the fluorescence quenching mechanism of DAPH-DNP and the reason for the fluorescence emission of DAPH were not discussed in detail.In this paper, a theoretical calculation approach was employed to provide a comprehensive account of the fluorescence quenching process of the probe.The potential energy curves were established in order to analyze the ESIPT process.To further investigate the photophysical characteristics and sensing mechanism of DAPH-DNP and DAPH, the charge transfer process was studied utilizing hole-electron [31,32] and IFCT [33] analysis methods.

DAPH-DNP
We first computed the spectral data with different functionals, and the results are listed in Table S1, which indicate that the calculated absorption and emission energies of DAPH-Enol at the B3LYP-D3(BJ) [34,35] level (2.49 eV/2.14 eV) are more comparable to the experimental values (2.68 eV/1.90 eV).However, it is widely recognized that the conventional functional B3LYP might tend to underestimate the energy of the charge transfer (CT) state.On the other hand, the inclusion of solvent effects within the linearresponse (LR) scheme may lead to an overestimation of the CT state energy.The interplay between these factors often results in an error-cancellation phenomenon in LR-B3LYP, effectively describing the energy of the CT state [36,37].By combining a range-separated functional with a corrected LR (cLR) solvation scheme, it is possible to achieve a genuinely accurate depiction of the CT state's energy.Therefore, we performed the geometry optimizations at the LR-B3LYP level, and then recalculated the energy landscape using the cLR-ωB97XD approach.
Figure 1a shows the optimized structure of DAPH-DNP in the ground (S 0 -min) state, and the major atoms are marked with numbers.The key geometric parameters of DAPH-DNP are also marked (Figure 1).In the S 0 state, the DNP group on the benzene ring is not parallel to the fluorophore (DAPH).The plane of the DNP group is approximately 145.5 • angle from the plane of the DAPH.As shown in Table 1, the excitation process is dominated by the S 0 → S 1 , corresponding to the electron transition from the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular orbital+2 (LUMO+2), and the oscillator strength is 1.8031, which indicates that the first (S 1 ) excited state is a bright state.The near-zero oscillator strength in S 0 → S 2 implies that the second (S 2 ) excited state of DAPH-DNP is a dark state, which corresponds to the electron transition from HOMO-4 → LUMO+2.The key bond lengths and dihedral angles of DAPH-DNP in the S 1 (LE-min) and S 2 (CT-min) states are listed in Figure 1b,c   The plotted FMOs (shown in Figure 2) show that no significant change in electron distribution is observed during the HOMO → LUMO+2 while a substantial charge separation occurs in the HOMO-4 → LUMO+2 transition.Thus, it can be preliminarily judged that the S1 and S2 states are the LE state and the CT state, respectively.To describe more accurately the characteristics of electron excitation, hole-electron analysis is used.For S0 → S1, there is incomplete separation of the distribution regions of holes and electrons (corresponding to the location of center of the blue and green iso-surfaces), the centers of the blue and green iso-surfaces are slightly farther apart, the D-index is 2.979 Å, and the t-index is negative (−0.122), indicating a small degree of charge transfer within the fluorophore.For S0 → S2, there is an obvious charge separation in which holes and electrons are mainly distributed in the o-hydroxyacetophenone group and DAPH group, respectively.Moreover, the centers of the hole and electron are relatively far apart, indicating that a significant charge transfer has occurred during the excitation process.As listed in Table 2, the D index of S0 → S2 excitation is as large as 2.225, the Sr index is very small (0.384), and the t index is 0.201, significantly larger than 0. Therefore, the S0 → S2 The plotted FMOs (shown in Figure 2) show that no significant change in electron distribution is observed during the HOMO → LUMO+2 while a substantial charge separation occurs in the HOMO-4 → LUMO+2 transition.Thus, it can be preliminarily judged that the S 1 and S 2 states are the LE state and the CT state, respectively.To describe more accurately the characteristics of electron excitation, hole-electron analysis is used.For S 0 → S 1 , there is incomplete separation of the distribution regions of holes and electrons (corresponding to the location of center of the blue and green iso-surfaces), the centers of the blue and green iso-surfaces are slightly farther apart, the D-index is 2.979 Å, and the t-index is negative (−0.122), indicating a small degree of charge transfer within the fluorophore.For S 0 → S 2 , there is an obvious charge separation in which holes and electrons are mainly distributed in the o-hydroxyacetophenone group and DAPH group, respectively.Moreover, the centers of the hole and electron are relatively far apart, indicating that a significant charge transfer has occurred during the excitation process.As listed in Table 2, the D index of S 0 → S 2 excitation is as large as 2.225, the Sr index is very small (0.384), and the t index is 0.201, significantly larger than 0. Therefore, the S 0 → S 2 are charge transfer excitations in the o-hydroxyacetophenone group in the N,N-dimethylamino direction.Meanwhile, we also used the IFCT method to quantitatively describe the contribution of each segment of DAPH-DNP to charge transfer.The detailed information is shown in Table S2 and Figure 3. DAPH-DNP is divided into two fragments (Figure 3a).The calculated results show that the contribution of fragment 2 in the S 1 state to electrons is as high as 96.87%, and differences in the electron transfer of fragment 1 (−12.84%)and fragment 2 (12.84%) indicate that only a small amount of electron transfer from fragment 2 to fragment 1 during S 0 → S 1 excitation occurs.Figure 3c shows that the fragment in the S 2 state has the highest contribution to holes (96.87%), which can be seen from the difference of charge transfer (Table S2) in that the electrons of fragment 1 are significantly reduced (−29.55%),while the electrons of fragment 2 are significantly increased.This means that a significant intramolecular charge transfer during S 0 → S 2 excitation is shown.The calculated vertical excitation energy (VEE) in the S 1 state is 3.14 eV (395 nm) (Table 1), which slightly overestimates the experimental date (2.77 eV, 447 nm).
shows that the fragment in the S2 state has the highest contribution to holes (96.87%), which can be seen from the difference of charge transfer (Table S2) in that the electrons of fragment 1 are significantly reduced (−29.55%),while the electrons of fragment 2 are significantly increased.This means that a significant intramolecular charge transfer during S0 → S2 excitation is shown.The calculated vertical excitation energy (VEE) in the S1 state is 3.14 eV (395 nm) (Table 1), which slightly overestimates the experimental date (2.77 eV, 447 nm).To further explain the fluorescence quenching mechanism of the probe, the calculated energy diagrams involved in the excitation and emission processes are plotted.As indicated in Figure 4, at both of the FC point and LE-min, the S1 and S2 states of DAPH-DNP are in LE and CT states, respectively, and the energy gaps between these two states (CT and LE state) are 0.62 eV and 0.58 eV, respectively.However, at the CT-min, the order of the CT and LE states is reversed, with the S1 state transitioning into a CT state.This suggests that the fluorescence quenching mechanism at the cLR-ωB97XD level involves a transformation from the LE to the CT state through a nearly barrierless minimal energy conical intersection (MECI).This mechanism, as reported in previous studies [38,39], aligns well with experimental observations.According to Kasha's rule [40], the fluorescence emitted by a molecule can only be excited from the S1 state.Since the S1 state has become a dark CT state at CT-min, this can only return to the ground state by non-radiative decay and is in accordance with the experimentally observed fluorescence quenching behavior of the probe.Briefly, the significant charge transfer from the DAPH unit to the DNP group in relaxation is well in line with the definition of the To further explain the fluorescence quenching mechanism of the probe, the calculated energy diagrams involved in the excitation and emission processes are plotted.As indicated in Figure 4, at both of the FC point and LE-min, the S 1 and S 2 states of DAPH-DNP are in LE and CT states, respectively, and the energy gaps between these two states (CT and LE state) are 0.62 eV and 0.58 eV, respectively.However, at the CT-min, the order of the CT and LE states is reversed, with the S 1 state transitioning into a CT state.This suggests that the fluorescence quenching mechanism at the cLR-ωB97XD level involves a transformation from the LE to the CT state through a nearly barrierless minimal energy conical intersection (MECI).This mechanism, as reported in previous studies [38,39], aligns well with experimental observations.According to Kasha's rule [40], the fluorescence emitted by a molecule can only be excited from the S 1 state.Since the S 1 state has become a dark CT state at CT-min, this can only return to the ground state by non-radiative decay and is in accordance with the experimentally observed fluorescence quenching behavior of the probe.Briefly, the significant charge transfer from the DAPH unit to the DNP group in relaxation is well in line with the definition of the typical d-PET (donor-excited PET) [41,42].

DAPH
The optimized geometries of DAPH and the corresponding parameters in the S0 (S0-min) and S1 (S1-min) states are shown in Figure 5.All structures are ensured to be located at the minimum energy point, without virtual frequency.Compared with the structure in the S0 state, the dihedral angles of O3-C4-C5-C6 and C7-C8-C9-C10 in the S1 state change from −179.924° to 179.980° and −0.029° to −0.004°, respectively, indicating that the geometry of DAPH tends to be flattened after photoexcitation.In addition, after the reaction with thiophenol, the DNP group of DAPH-DNP is eliminated, releasing the phenol group, which can form an intramolecular hydrogen bond with the adjacent O atom.The O1-H2 bond length of DAPH is increased from 1.003 Å (S0-min) to 1.025 Å (S1-min), and the H2…O3 bond length is reduced from 1.582 Å (S0-min) to 1.499 Å (S1-min).Moreover, the bond angles of O1-H2…O3 increased by 4.0° from 150.9° (S0-min) to (S1-min).After photoexcitation, the strength of the hydrogen bond is increased, thereby promoting the ESIPT process.To further verify whether the hydrogen bonding in the S1 state is enhanced, the calculated IR vibrational spectra of the relevant hydrogen bond are presented.The results depicted in Figure 6 indicate a significant red-shift of 403 cm −1 in the IR vibrational frequency of O1-H2, from 3035 cm −1 (S0-min) to 2632 cm −1 (S1-min), confirming the enhancement of the hydrogen bonding phenomenon.This observation is consistent with the structural analysis view.Additionally, we conducted single point energy calculations using the COSMO-ADC(2) method within the perturbation theory on the energy and density (PTED) reaction field scheme as implemented in Turbomole [43].The COSMO-ADC(2)/PTED results (Table S3) exhibit a similar energy landscape to that of the cLR-ωB97XD.Moreover, the computed vertical excitation energy is in closer agreement with the experimentally measured absorption maximum compared to the cLR-ωB97XD.However, at the CT-min, the computed energy gap between the ground state and the CT state is only 0.08 eV, suggesting that the COSMO-ADC(2)/PTED might underestimate the energy of the CT state, in line with previous benchmark studies [44].

DAPH
The optimized geometries of DAPH and the corresponding parameters in the S 0 (S 0 -min) and S 1 (S 1 -min) states are shown in Figure 5.All structures are ensured to be located at the minimum energy point, without virtual frequency.Compared with the structure in the S 0 state, the dihedral angles of O 3 -C 4 -C 5 -C 6 and C 7 -C 8 -C 9 -C 10 in the S 1 state change from −179.924 • to 179.980 • and −0.029 • to −0.004 • , respectively, indicating that the geometry of DAPH tends to be flattened after photoexcitation.In addition, after the reaction with thiophenol, the DNP group of DAPH-DNP is eliminated, releasing the phenol group, which can form an intramolecular hydrogen bond with the adjacent O atom.The O 1 -H 2 bond length of DAPH is increased from 1.003 Å (S 0 -min) to 1.025 Å (S 1 -min), and the H 2 . . .O 3 bond length is reduced from 1.582 Å (S 0 -min) to 1.499 Å (S 1 -min).Moreover, the bond angles of O 1 -H 2 . . .O 3 increased by 4.0 • from 150.9 • (S 0 -min) to 154.9 • (S 1 -min).After photoexcitation, the strength of the hydrogen bond is increased, thereby promoting the ESIPT process.To further verify whether the hydrogen bonding in the S 1 state is enhanced, the calculated IR vibrational spectra of the relevant hydrogen bond are presented.The results depicted in Figure 6 indicate a significant red-shift of 403 cm −1 in the IR vibrational frequency of O 1 -H 2 , from 3035 cm −1 (S 0 -min) to 2632 cm −1 (S 1 -min), confirming the enhancement of the hydrogen bonding phenomenon.This observation is consistent with the structural analysis view.
thereby promoting the ESIPT process.To further verify whether the hydrogen bonding in the S1 state is enhanced, the calculated IR vibrational spectra of the relevant hydrogen bond are presented.The results depicted in 6 indicate a significant red-shift of 403 cm −1 in the IR vibrational frequency of O1-H2, from 3035 cm −1 (S0-min) to 2632 cm −1 (S1-min), confirming the enhancement of the hydrogen bonding phenomenon.This observation is consistent with the structural analysis view.PT process entails a rearrangement of the electron structure within the π system, which subsequently impacts the molecule's geometry.To gain a deeper understanding of the PT process, we generated PECs for DAPH in both the S0 and S1 states by conducting a relaxed potential energy scan along the O1-H2 bond (as illustrated in Figure 7).The results show that in the S0 state, the energy of DAPH system continues to rise along the direction of increasing bond length, and there is no low energy point, which means that DAPH-keto cannot exist stably.The PEC of the S1 state shows that DAPH-enol can be transformed into DAPH-keto at a low energy point after overcoming an energy barrier of 2.36 kcal/mol at the cLR-ωB97XD/TZVP level.The forward reaction exhibits a higher energy barrier compared to the reverse process (0.65 kcal/mol), suggesting that the conversion of DAPH-keto to DAPH-enol can occur with ease.The above results illustrate that DAPH is not prone to the GSIPT or ESIPT process; this is the reason why DAPH-keto could not be optimized in the S0 and S1 states.Therefore, the onset of the fluorescence emission of DAPH may not be caused by the ESIPT mechanism proposed in the experiment.It is important to highlight that a relaxed potential energy scan may yield energy profiles that differ significantly from those obtained through a rigid potential energy scan.To verify that the relaxation of molecular geometry does not introduce qualitative changes or alter our conclusions, we also presented the results of a rigid potential energy scan (see Figure S7), which supports the same conclusion.PT process entails a rearrangement of the electron structure within the π system, which subsequently impacts the molecule's geometry.To gain a deeper understanding of the PT process, we generated PECs for DAPH in both the S 0 and S 1 states by conducting a relaxed potential energy scan along the O 1 -H 2 bond (as illustrated in Figure 7).The results show that in the S 0 state, the energy of DAPH system continues to rise along the direction of increasing bond length, and there is no low energy point, which means that DAPH-keto cannot exist stably.The PEC of the S 1 state shows that DAPH-enol can be transformed into DAPH-keto at a low energy point after overcoming an energy barrier of 2.36 kcal/mol at the cLR-ωB97XD/TZVP level.The forward reaction exhibits a higher energy barrier compared to the reverse process (0.65 kcal/mol), suggesting that the conversion of DAPH-keto to DAPH-enol can occur with ease.The above results illustrate that DAPH is not prone to the GSIPT or ESIPT process; this is the reason why DAPH-keto could not be optimized in the S 0 and S 1 states.Therefore, the onset of the fluorescence emission of DAPH may not be caused by the ESIPT mechanism proposed in the experiment.It is important to highlight that a relaxed potential energy scan may yield energy profiles that differ significantly from those obtained through a rigid potential energy scan.To verify that the relaxation of molecular geometry does not introduce qualitative changes or alter our conclusions, we also presented the results of a rigid potential energy scan (see Figure S7), which supports the same conclusion.To study the fluorescence enhancement mechanism, we calculated the absorption and emission spectra of DAPH at the ωB97XD/TZVP level.Moreover, the vertical excitation and emission energies of DAPH involved in orbital transitions, as well as the oscillator strengths, are listed in Table 1.The calculated results show that the S0 → S1 is the dominant excitation pathway characterized by the transition from HOMO → LUMO, which blocks the occurrence of the PET process.And the excitation energy of DAPH is 2.96 eV (419 nm), which is very consistent with the experiment value (2.68 eV, 462 nm).Compared with DAPH-DNP, the maximum absorption peak of DAPH is shifted in the direction of long wavelength, the 24 nm red-shift indicates that the n-π conjugation effect of the DNP group on the probe is eliminated after the thiolysis reaction.The calculated emission data show that the oscillator strength of DAPH-enol during emission is 1.8406, which means that the decay of DAPH-enol from the S1 state to the S0 state occurs via radiative transition.In addition, the emission energy of DAPH-enol is also calculated to be 2.74 eV (453 nm), which is significantly overestimated by the ωB97XD functional.This is not surprising because the range-separated functional typically overestimates the energy of the LE state.We then computed the computed emission energy at the LR-B3LYP level, which is 2.14 eV and agrees well with the experimental value (1.90 eV).That is to say, the observed fluorescence should be emitted from DPAH-enol, which agrees with the conclusion we discussed in the PEC section.
It is worth mentioning that, according to our conclusion, the significant Stokes shift observed is not attributed to the ESIPT mechanism and requires further investigation to elucidate the underlying factors contributing to the observed phenomenon.As shown in Figure 8, the contribution of the phenol group to HOMO is not obvious, while the contribution to LUMO is significantly increased; the electron distribution of LUMO delocalizes over the DAPH conjugate system after photoexcitation, which indicates that DAPH in the S1 state exhibits ICT character.Moreover, to describe the characters of electron excitation more accurately, the hole-electron distribution of DAPH was analyzed.Compared to the hole, the centroid of the electrons shifts in the direction of the phenol group.As listed in Table 2, the D index is 2.159 Å, and the t-index has a negative value (−1.399Å), significantly less than 0, indicating an obvious separation between holes and elec- To study the fluorescence enhancement mechanism, we calculated the absorption and emission spectra of DAPH at the ωB97XD/TZVP level.Moreover, the vertical excitation and emission energies of DAPH involved in orbital transitions, as well as the oscillator strengths, are listed in Table 1.The calculated results show that the S 0 → S 1 is the dominant excitation pathway characterized by the transition from HOMO → LUMO, which blocks the occurrence of the PET process.And the excitation energy of DAPH is 2.96 eV (419 nm), which is very consistent with the experiment value (2.68 eV, 462 nm).Compared with DAPH-DNP, the maximum absorption peak of DAPH is shifted in the direction of long wavelength, the 24 nm red-shift indicates that the n-π conjugation effect of the DNP group on the probe is eliminated after the thiolysis reaction.The calculated emission data show that the oscillator strength of DAPH-enol during emission is 1.8406, which means that the decay of DAPH-enol from the S 1 state to the S 0 state occurs via radiative transition.In addition, the emission energy of DAPH-enol is also calculated to be 2.74 eV (453 nm), which is significantly overestimated by the ωB97XD functional.This is not surprising because the range-separated functional typically overestimates the energy of the LE state.We then computed the computed emission energy at the LR-B3LYP level, which is 2.14 eV and agrees well with the experimental value (1.90 eV).That is to say, the observed fluorescence should be emitted from DPAH-enol, which agrees with the conclusion we discussed in the PEC section.
It is worth mentioning that, according to our conclusion, the significant Stokes shift observed is not attributed to the ESIPT mechanism and requires further investigation to elucidate the underlying factors contributing to the observed phenomenon.As shown in Figure 8, the contribution of the phenol group to HOMO is not obvious, while the contribution to LUMO is significantly increased; the electron distribution of LUMO delocalizes over the DAPH conjugate system after photoexcitation, which indicates that DAPH in the S 1 state exhibits ICT character.Moreover, to describe the characters of electron excitation more accurately, the hole-electron distribution of DAPH was analyzed.Compared to the hole, the centroid of the electrons shifts in the direction of the phenol group.As listed in Table 2, the D index is 2.159 Å, and the t-index has a negative value (−1.399Å), significantly less than 0, indicating an obvious separation between holes and electrons.Meanwhile, the IFCT analysis (Figure 9) is provided to describe quantitatively the charge transfer characteristics of DAPH, and is divided into three fragments.Table S4 shows that the contributions of fragment 2 and fragment 3 to electrons and holes are 70.68% and 50.54%, respectively, and the charge transfer differences of fragment 1, fragment 2, and fragment 3 in the excitation of S 0 → S 1 are 2.31%, 27.79%, and −30.10%, respectively, indicating that the dominant contribution of electron transfer is from fragment 3 to fragment 2 during excitation.Thus, it further confirms the ICT character of the S 1 state, and the large Stokes shift observed in the experiment could be caused by the ICT process of DAPH.
Molecules 2023, 28, x FOR PEER REVIEW 10 of 15 trons.Meanwhile, the IFCT analysis (Figure 9) is provided to describe quantitatively the charge transfer characteristics of DAPH, and is divided into three fragments.Table S4 shows that the contributions of fragment 2 and fragment 3 to electrons and holes are 70.68% and 50.54%, respectively, and the charge transfer differences of fragment 1, fragment 2, and fragment 3 in the excitation of S0 → S1 are 2.31%, 27.79%, and −30.10%, respectively, indicating that the dominant contribution of electron transfer is from fragment 3 to fragment 2 during excitation.Thus, it further confirms the ICT character of the S1 state, and the large Stokes shift observed in the experiment could be caused by the ICT process of DAPH.trons.Meanwhile, the IFCT analysis (Figure 9) is provided to describe quantitatively the charge transfer characteristics of DAPH, and is divided into three fragments.Table S4 shows that the contributions of fragment 2 and fragment 3 to electrons and holes are 70.68% and 50.54%, respectively, and the charge transfer differences of fragment 1, fragment 2, and fragment 3 in the excitation of S0 → S1 are 2.31%, 27.79%, and −30.10%, respectively, indicating that the dominant contribution of electron transfer is from fragment 3 to fragment 2 during excitation.Thus, it further confirms the ICT character of the S1 state, and the large Stokes shift observed in the experiment could be caused by the ICT process of DAPH.

Computational Details
Gaussian 16 program package [45] was utilized for all theoretical calculations con- structure information and infrared vibration spectra data show that the intramolecular hydrogen bonds of DAPH are enhanced in the S 1 state.The calculated PECs indicate that the ESIPT process could not occur in DAPH due to the unstable nature of the keto structure; it quickly reverted back to the enol structure.ICT is the main cause of the large Stokes shift, which is inconsistent with the conclusion proposed in the experiment.This further suggests that the fluorescence enhancement of DAPH is not based on the ESIPT process.By the analysis of the FMOs of DAPH, the process of electronic excitation was found to be accompanied by the ICT process, indicating that the large Stokes shift observed in the experiment can be attributed to the ICT characteristic in the S 1 state of DAPH-enol.Our theoretical study not only explains the sensing mechanism but also confirms the inefficiency of the ESIPT process in the probe DAPH-DNP, which is of great significance for the synthesis of novel probes in the future.

Scheme 1 . 1 .
Scheme 1. Sensing mechanism of probe DAPH-DNP towards thiophenols.In this paper, a theoretical calculation approach was employed to provide a comprehensive account of the fluorescence quenching process of the probe.The potential energy curves were established in order to analyze the ESIPT process.To further investigate the photophysical characteristics and sensing mechanism of DAPH-DNP and

Figure 1 .
Figure 1.The optimized structures of DAPH-DNP in the S 0 (a), LE (S 1 ) (b), and CT (S 2 ) states (c).The labeling of atomic color: O: red; C: blue; H: white; N: yellow, the bond length in Å, the dihedral angle in degrees.

Figure 2 .
Figure 2. Excitation processes of DAPH-DNP.Molecular orbitals are given in blue and red iso-surfaces, holes and electrons are given in blue and green iso-surfaces, respectively (ωB97XD/TZVP/IEFPCM).

Figure 2 .
Figure 2. Excitation processes of DAPH-DNP.Molecular orbitals are given in blue and red iso-surfaces, holes and electrons are given in blue and green iso-surfaces, respectively (ωB97XD/TZVP/IEFPCM).

Figure 3 .
Figure 3.The IFCT analyzing the electron excitation process of DAPH-DNP molecular fragments.(a) The amount of electron transfer between fragments from S 0 to S 1 (b) and S 2 states (c) (ωB97XD/TZVP/IEFPCM).

Figure 5 .
Figure 5.The optimized geometries of DAPH in the S0 and S1 states.The labeling of atomic color: O: red; C: blue; H: white; N: yellow, the bond length in Å, the angles are in degrees.

Figure 5 .
Figure 5.The optimized geometries of DAPH in the S0 and S1 states.The labeling of atomic color: O: red; C: blue; H: white; N: yellow, the bond length in Å, the angles are in degrees.

Figure 5 . 15 Figure 6 .
Figure 5.The optimized geometries of DAPH in the S 0 and S 1 states.The labeling of atomic color: O: red; C: blue; H: white; N: yellow, the bond length in Å, the angles are in degrees.Molecules 2023, 28, x FOR PEER REVIEW 8 of 15

Figure 6 .
Figure 6.The calculated IR spectra of the DAPH in the spectral region of the O 1 -H 2 stretching band calculated.

Figure 7 .
Figure 7.The PECs of the S 0 and S 1 states for DAPH along with the O 1 -H 2 bond length (ωB97XD/TZVP/IEFPCM).

Figure 8 .
Figure 8. Excitation processes of DAPH.Molecular orbitals are given in blue and red iso-surfaces, holes and electrons are given in blue and green iso-surfaces, respectively (ωB97XD/TZVP/IEFPCM).

Figure 9 .
Figure 9.The DAPH molecular fragments electron excitation process analyzed by the IFCT (a) and the amount of inter-fragment electron transfer from the S0 to S1 States (b) (ωB97XD/TZVP/IEFPCM).As Figure10shows, at the FC point, the S1 state is a bright ICT state.By calculating the MOs, the LUMO transformation of DAPH includes intramolecular CT (ICT) from o-hydroxyacetophenone to N,N-dimethylamino.Thus, under photoexcitation, DAPH relaxes directly from the FC point to S1-min via vibration relaxation, then returns to the ground state by fluorescence emission.

Figure 8 .
Figure 8. Excitation processes of DAPH.Molecular orbitals are given in blue and red iso-surfaces, holes and electrons are given in blue and green iso-surfaces, respectively (ωB97XD/TZVP/IEFPCM).

Figure 8 .
Figure 8. Excitation processes of DAPH.Molecular orbitals are given in blue and red iso-surfaces, holes and electrons are given in blue and green iso-surfaces, respectively (ωB97XD/TZVP/IEFPCM).

Figure 9 .
Figure 9.The DAPH molecular fragments electron excitation process analyzed by the IFCT (a) and the amount of inter-fragment electron transfer from the S0 to S1 States (b) (ωB97XD/TZVP/IEFPCM).As Figure10shows, at the FC point, the S1 state is a bright ICT state.By calculating the MOs, the LUMO transformation of DAPH includes intramolecular CT (ICT) from o-hydroxyacetophenone to N,N-dimethylamino.Thus, under photoexcitation, DAPH relaxes directly from the FC point to S1-min via vibration relaxation, then returns to the ground state by fluorescence emission.

Figure 9 .
Figure 9.The DAPH molecular fragments electron excitation process analyzed by the IFCT (a) and the amount of inter-fragment electron transfer from the S 0 to S 1 States (b) (ωB97XD/TZVP/IEFPCM).As Figure10shows, at the FC point, the S 1 state is a bright ICT state.By calculating the MOs, the LUMO transformation of DAPH includes intramolecular CT (ICT) from ohydroxyacetophenone to N,N-dimethylamino.Thus, under photoexcitation, DAPH relaxes directly from the FC point to S 1 -min via vibration relaxation, then returns to the ground state by fluorescence emission.
: The B3LYP/TZVP/IEFPCM calculated energies of DAPH-DNP showing the PET mechanism; Figure S2: Excitation processes of DAPH-DNP (molecular orbitals are given in blue and red iso-surfaces, holes and electrons are given in blue and green iso-surfaces, respectively); Figure S3: The IFCT analyzing the electron excitation process of DAPH-DNP molecular fragments (a), the amount of electron transfer between fragments from the S 0 to S 1 (b) and S 3 states (c); Figure S4: Frontier molecular orbitals showing the excitation process of DAPH; Figure S5: The B3LYP/TZVP/IEFPCM calculated energies of DAPH showing the ICT mechanism; Figure S6: The PECs of the S 0 and S 1 states for DAPH along with the O 1 -H 2 bond length; Figure S7: The PECs of the S 0 and S 1 states for DAPH along with the O 1 -H 2 bond length (rigid scan).Author Contributions: Y.W. (Yuxi Wang): Investigation, Writing-original draft.M.Z.: Investigation, Data curation, Formal analysis.W.L.: Data curation.Y.W. (Yi Wang): Conceptualization, Writingreview and editing.P.Z.: Writing-review and editing, supervision, funding acquisition.All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by the Natural Science Foundation of Shandong Province, China (Grant No: ZR2023MB146).Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.

This indicates that the fluorescence quenching behavior of DAPH-DNP may be affected by the DNP group on the benzene ring. Table 1. Detailed theoretical and experimental spectral data for DAPH-DNP and
. The parameters of the optimized S 1 -state geometry are similar to those in the S 0 state.The main difference is that the dihedral angle of ∠C 4 C 6 C 7 O 8 is reduced from 46.7 • to 34.8 • , which means that the flatness of the DAPH-DNP increases after reaching the S 1 state.Compared with the S 0 state, the structural changes in the S 1 state are mainly manifested by the slight torsion of the DNP group.Unlike the S 1 state, the dihedral angle of ∠C 1 C 2 O 3 C 4 in the S 2 state changes from 132.7 • to 119.3 • , and the C 2 -O 3 bond length is increased from 1.353 Å to 1.377 Å. DPAH (ωB97XD/TZVP/IEFPCM).

Table 2 .
The calculated results of the excited states for DAPH-DNP and DAPH, including the centroid distance (D), the degree of overlap (Sr), the width distribution (H), degree of separation (t),

Table 2 .
The calculated results of the excited states for DAPH-DNP and DAPH, including the centroid distance (D), the degree of overlap (Sr), the width distribution (H), degree of separation (t), hole delocalization index (HDI), and electron delocalization index (EDI) (ωB97XD/TZVP/IEFPCM).