B. Nakwanich, A. Koonwong, A. Suramitr et al.
Journal of Molecular Structure 1245 (2021) 131132
Scheme 1. Potential mechanism for F− detection of all sensors.
3.2. Theoretical calculation
SA1keto and SA1keto-F− complexes. For the SA1-enol and SA1-
keto, it shown that the O-H stretching vibrations of the SA1-enol
appeared at 3760 cm−1 and 3322 cm−1, while SA1keto’s O–H and
To investigate the active site of SA1, SA2, and SA3 molecules for
F− detection, theoretical calculations were conducted. Electrostatic
potential (ESP) analysis was undertaken using the Gauss View 5.0
program. The ESP molecular surfaces of the sensor molecules are
presented in Fig. 5. For all sensing compounds, it is revealed that
the H atom of the OH moieties gave a positive ESP surface, while
the aromatic unit showed a negative ESP surface. This resulted in
the preference of F− anions to close to the positively charged H
atoms rather than to the negatively charged aromatic units. In con-
clusion, detection of the sensor compounds can be performed via
interaction between F− anion and H atom of the OH moieties of
the sensors.
Structures of the SAenol, SAketo, and SAketo-F− intermediate and
the deprotonated form were optimized separately (Fig. 6) in DFT
formalism with 6-311 G(d,p) basis set and the CAM-B3LYP func-
tional. At each optimized structure, lowest few electronic transi-
tions of the molecules (SAenol, SAketo, and SAketo-F−) were exam-
ined with TDDFT formalism. All calculations were carried out in
acetonitrile as solvent, treated in Onsager’s SCRF formalism, with
polarisable continuum (PCM) approximation. A single imaginary
frequency confirmed the obtained structure as the desired transi-
tion state.
N–H stretching vibrations appeared at 3777 cm−1 and 2493 cm−1
,
respectively. Compared to SA-keto+F− complex, the N–H stretching
vibration of complex appeared at 3260 cm−1, indicating large red-
shifts. All calculation in Fig. 8 indicated that there were strong in-
teractions in the O•••H–X (X = N and F) hydrogen bond and there
were large red-shifts for H–X stretching vibrational frequencies in
the SA-enol, SA-keto and SA-keto+F− complexes.
The interaction energies with BSSE corrections for com-
plexes (SA1-keto+F−, SA2-keto+F− and SA3-keto+F−) are listed in
Table 3. All the values of the adsorption energies being negative, it
is implied that the adsorption structures are stable, and the pro-
cess is exothermic. It is clear that the BSSE-corrected interaction
energy (ꢀEBSSE) of complex SA1-keto+F−, SA2-keto+F− and SA3-
keto+F− are favorable for the distinct selectivity of F−. These calcu-
lation results are in good agreement with the reported experimen-
tal observations that intermolecular proton transfer (IPT) between
chemosensor substrate SA1, SA2 and SA3. Thus, one can conclude
that the host chemosensors have a much stronger affinity to F− ion
through intermolecular proton transfer, which leads to the forma-
tion of chemosensor anions by F−.
To simulate the absorption spectra of the molecules, vertical ex-
citation energies of the SAenol, SAketo, and SAketo-F− for SA1, SA2,
and SA3 were calculated based on the time-dependent density
functional theory (TD-DFT) method using their ground-state opti-
mized structures. The calculated electronic transitions energies and
corresponding transition oscillator strengths for transition from the
ground state (S0) to the excited state (S1) of the sensor molecules
as well as the experimental UV-vis absorption data are summa-
The optimized structure of the SAenol, SAketo, and SAketo-F− form
confirmed the mechanism of fluoride binding (Scheme 1) through
hydrogen bond formation between the O–H and N–H with F−.
Some geometrical parameters were shown in Fig. 6. The lengths
of O–H of SA1enol, SA1keto, and SA1keto-F− were increased to 0.991
˚
˚
˚
A, 1.561 A and 1.998 A, respectively, while the N–H lengths were
˚
˚
˚
decreased to 1.741 A, 1.063 A and 1.034 A, respectively. Moreover,
this study revealed that the NH–O and N–HO hydrogen bonds in
Figs. 6 and S2 are involved in NH and OH tautomers from which
it could be inferred that F− binding occurs through the SAketo-F−
form, followed by deprotonation of the ligand, formation of hydro-
gen bond between deprotonated oxygen, and -NH and release of
HF. In addition to, the frequency calculations were performed on
each optimized structure. Similar phenomenon was also observed
for SA2 and SA3 (See in Supporting Information, Fig. S3).
From the DFT/TDDFT studies, it was observed that the theoreti-
cal absorption peak (λabs) of SA1, SA2, and SA3 at 321, 327, and 329
nm of the enol form were caused due to H→L transition which
was compared with experimental values (λexpt) at 334, 343, and
349 nm, respectively. The theoretical absorption peaks of the keto
form, at 368, 382, and 383 nm may be compared with the exper-
imental values at 380, 398, and 400 nm, respectively. It is found
that all the calculated wavelengths (λcal) showed a slight blue shift
compared to that obtained from the experiments with an under-
estimated error of ~3-7%, which was in an acceptable range. The
absorption of the probable intermediate (SAketo-F−) were obtained
theoretically at 378, 386, and 395 nm which is very close to the ex-
perimental value of keto at 380, 398, and 400 nm, respectively. The
Frequency calculations were carried out on each optimized
structure, and their IR spectra were discussed. Fig. 7 showed the
simulated infrared (IR) spectra for the SA-enol, SA-keto and SA-
keto+F−, where the intensity was plotted against the harmonic vi-
brational frequencies. For the complexity of vibrational modes, it
was difficult to attribute all bands, so we have only analyzed some
H–X (X = O and N) vibrational frequencies [36-38] in the SA1enol,
7