Open access peer-reviewed chapter

Structure-Property Relationships in Benzofurazan Derivatives: A Combined Experimental and DFT/TD-DFT Investigation

Written By

Hanen Raissi, Imen Chérif, Hajer Ayachi, Ayoub Haj Said, Fredj Hassen, Sahbi Ayachi and Taoufik Boubaker

Submitted: 31 May 2021 Reviewed: 05 July 2021 Published: 23 November 2021

DOI: 10.5772/intechopen.99246

From the Edited Volume

Density Functional Theory - Recent Advances, New Perspectives and Applications

Edited by Daniel Glossman-Mitnik

Chapter metrics overview

295 Chapter Downloads

View Full Metrics


In this work we seek to understand and to quantify the reactivity of benzofurazan derivatives toward secondary cyclic amines, like pyrrolidine, piperidine and morpholine, acting as nucleophile groups in SNAr reactions. For this aim, physico-chemical and structural descriptors were determined experimentally and theoretically using the DFT/B3LYP/6-31+ g (d,p) methodology. Thus, different 4-X-7-nitrobenzofurazans (X = OCH3, OC6H5 and Cl) and products corresponding to the electrophilic aromatic substitution by pyrrolidine, piperidine and morpholine, were investigated. Particularly, the HOMO and LUMO energy levels of the studied compounds, determined by Cyclic Voltammetry (CV) and DFT calculations, were used to evaluate the electrophilicity index (ω). The latter was exploited, according to Parr’s approach, to develop a relationship which rationalizes the kinetic data previously reported for the reactions of the 4-X-7-nitrobenzofurazans with nucleophiles cited above. Moreover, the Parr’s electrophilicity index (ω) of these benzofurazans determined in this work were combined with their electrophilicity parameters (E), reported in preceding papers, was found to predict the unknown electrophilicity parameters E of 4-piperidino, 4-morpholino and 4-pyrrolidino-7-nitrobenzofurazan. In addition, the relationship between the Parr’s electrophilicity index (ω) and Hammett constants σ, has been used as a good model to predict the electronic effect of the nucleophile groups. Finally, we will subsequently compare the electrophilicity index (ω) and the electrophilicity parameters (E) of these series of 7-X-4-nitrobenzofurazans with the calculated dipole moment (μ) in order to elucidate general relationships between E, ω and μ.


  • Benzofurazan
  • Structure–property correlations
  • Optical properties
  • Mayr’s approach
  • Electrophilicity
  • Parr’s approach
  • Voltammetry Cyclic
  • DFT
  • ICT
  • NLO properties

1. Introduction

4-Nitrobenzofurazan (NBD) is neutral 10π electron-deficient heteroaromatic substrates. Organic materials containing the NBD moiety have been tested in potential biomedical and bio-analytical applications [1, 2, 3, 4, 5, 6, 7]. They have been also used as an effective electron-withdrawing building block for organic solar cell materials [8], and have been discovered to be nonlinear optical (NLO) materials [9]. In addition, the strong electron withdrawing properties of nitrobenzofurazan (NBD) derivatives make it an exceptional electrophilic product of the nucleophilic aromatic substitution (SNAr).

In addition, the activated chlorine atom in 4-chloro-7-nitrobenzofurazan (NBD chloride, namely NBD-Cl) can undergo electrophilic aromatic substitution by phenoxide, primary or secondary cyclic amines [10, 11, 12, 13] and some of 4,7 di-substituted benzofurazan compounds to generate weakly or non-fluorescents products [14, 15, 16, 17, 18].

A previous kinetic study of amino-substituted NBD showed that methoxy and aryloxy substituent’s are better leaving groups than chloro and that a red-colored Meisenheimer complex immediate was formed during the reactions [19].

On the other hand, the frontier molecular orbital descriptors including HOMO and LUMO energy levels play major roles in governing many chemical reactions as function of reactivity and stability of related molecules [20, 21]. Alternatively, contacts from intra-molecular interactions have been identified as an important driving force on the stabilization of various planar molecular structures [22, 23]. However, these interactions lead to intra-molecular charge transfer of formed compounds and thus provides a non-radiative decay from the excited state [24, 25, 26].

Recently, we reported a kinetic study for SNAr reactions of NBD with secondary amines [27, 28, 29, 30]. The present chapter is focused on the elucidation of the combined experimental–theoretical investigation relationships between the optical (UV–Vis, PL and time resolved photoluminescence (TR-PL), electrochemical (CV) properties and the chemical structures of 4.7-di-substituted benzofurazans (NBDs), based on the kinetic reactivity process of compounds. Parr’s approach was employed to develop a relationship which rationalizes the kinetic data previously reported for the reactions of 4-X-7-nitrobenzofurazans with various nucleophiles [27, 28, 29].

Then, the substituent effects were particularly limited only to morpholine, piperidine and pyrrolidine groups. Thus, we seek to study the structure–property relationships for molecules containing 4-nitrobenzofurazan moiety for a better understanding of their optical behavior and their chemical reactivity. For this aim, theoretical calculations using quantum chemical calculations within density functional theory (DFT) and its extension TD-DFT in the 6-31 + g(d,p) basis set have been developed to facilitate the in-depth understanding of structure–property relationships.

The Molecular structures of the investigated compounds are illustrated in Figures 1 and 2.

Figure 1.

Chemical structures of the studied electrophile compounds.

Figure 2.

Reaction mechanism between NBD-Cl and three examples of nucleophiles (morpholine, piperidine and pyrrolidine) groups.

To the best of our knowledge, no similar study has been published on electrochemical and photo-physical properties of the studied compounds. Cyclic Voltammetry (CV) is used to estimate the band gaps, electron affinities (i.e. the energy of LUMO, εLUMO) and ionization potentials (i.e., the energy of HOMO, εHOMO) [31]. Besides, the computed results from DFT and TD-DFT calculations enable the access to the structure geometry, the Mullikan charge distribution and the dipole moment of the 4,7 di-substituted benzofurazan at the ground- (S0) and excited (S1)-states [32, 33].


2. Experimental details

2.1 Materials

The 7-chloro-4-nitrobenzofurazan (NBD-Cl) and the tetraethyl ammonium tetrafluoroborate (TEAF) were commercial products (Aldrich, Darmstadt-Germany). Samples of 7-methoxy-4-nitrobenzofurazan (NBD-OCH3) and 7-phenoxy-4-nitrobenzofurazan (NBD-OC6H5) were prepared according to the standard methods described in the literature [27, 28, 29]. The secondary cyclic amines (morpholine, piperidine and pyrrolidine) were commercial products (Aldrich, Darmstadt-Germany) and were redistilled before use whenever necessary.

2.2 Apparatus

Cyclic Voltammetry (CV) measurements were carried out under nitrogen gas with a Voltalab 10 apparatus from Radiometer driven by the Volta Master software at a potential scan rate of 0.1 V/s, at room temperature, in 25 mL of CH3CN solution containing TEAF (0.1 M) as a supporting electrolyte. The cell used three electrodes. The working electrode was a 2 mm diameter platinum disk (Tacussel type EDI) and a platinum (Pt) wire as a counter electrode. The concentrations of the substrates were 10−3 M. The potentials EOx and ERed give us the information on HOMO and LUMO energy levels of studied materials.

The UV–visible absorption spectra of compounds have carried out on a Shimadzu UV–visible model 1650 spectrometer. The samples were dissolved in acetonitrile medium. The optical properties of the sample are secondly investigated by photoluminescence spectroscopy. The 514.5 nm line of the continuous-wave argon Ar + laser is used as an excitation source. Spectral analysis of the PL was performed using a Jobin Yvon iHR320 monochromator through the lock-in amplifier technique.

For the Time-Resolved PL (TRPL) study, the sample was illuminated by a pulsed laser diode with photon energy of 1.58 eV with a repetition rate of 80 MHz, using the time correlated photo counting technique. The details experiment was reported elsewhere [34].

2.3 Computational methods

Ground-state (S0) optimized geometries of the studied compounds were performed by means of DFT at B3LYP level of theory [35, 36] in conjunction with the 6-31+ g(d,p) basis set. No constraints were used and all structures were free to optimize in an acetonitrile solution by applying the conductor polarizable continuum model (CPCM) with the Gaussian 09 program [37] and the output files were visualized with the Gauss View (5.0.8) molecular visualization program [38]. To predict the optical absorption properties of the studied compounds, UV–Vis absorption spectra were simulated at the time-dependent TD-DFT level with the same basis set. The electronic properties such as highest occupied molecular orbital (εHOMO) and lowest unoccupied molecular orbital (εLUMO) energy levels and dipole moment have been extracted from calculations.

The semi-empirical quantum-chemical ZINDO level was used to predict the optical emission spectra on S1 optimized structures. Electronic transitions assignment and oscillator strengths were also calculated using the same method of calculations.


3. Results and discussion

3.1 Electrochemical analysis

Electrochemical comparative studies were undertaken using Cyclic Voltammetry (CV) for the six compounds 4.7 di-substituted benzofurazan (See Figures 3 and 4). The CVs profiles were performed in the potential range of −2.0 to +3.0 V.

Figure 3.

Oxidation and reduction voltammograms of NBD-Cl, NBD-OCH3 and NBD-OC6H5 products (10−3 M) recorded in acetonitrile.

Figure 4.

Oxidation and reduction voltammograms of NBD-Pyrr, NBD-pip and NBD-Morph products (10−3 M) recorded in acetonitrile.

The potentials (EOX and ERed) peaks can be used to determine the HOMO and LUMO energy levels of the studied compounds. Indeed, the value of their energy gaps can be estimated from the potentials difference.

The highest occupied molecular orbital (HOMO) energy and the lowest unoccupied molecular orbital (LUMO) energy were estimated from the first oxidation and reduction potentials, respectively. The energy levels were calculated according to an empirical method [39] and by assuming that the energy level of the ferrocene/ferrocenium is 4.8 V below the vacuum level.


Here, EOx and ERed are the potential of the oxidation peak and reduction peak, respectively. These values were obtained from the CV curve for each compound; EFc/Fc+0 is 0.59 eV, the average of oxidation and reduction peak potentials of ferrocene measured under the same experimental condition; 4.8 eV is the energy difference between the energy level of ferrocene and the vacuum. Note that when the oxidation/reduction peak of the studied compound was reversible, the Eox/red was replaced in the equation by E°ox/red determined experimentally.

For the electrophile compounds (Figure 3 and Table 1), the peaks appeared at the reduction potentials of about: -0.440 V (NBD-Cl), −0.716 (NBD-OCH3) and − 0.610 V (NBD-OC6H5). They correspond to LUMO energy levels at −3.77 eV, −3.49 eV and − 3.60 eV, respectively. In the same way, the oxidation potentials located at: 2.90 V, 2.28 V and 2.40 V led to HOMO energy level values of −7.11 eV, −6.49 eV and − 6.61 eV for NBD-Cl, NBD-OCH3 and NBD-OC6H5, respectively. The electrochemical energy gaps deduced from these measurements were: 3.34 eV (NBD-Cl), 3.00 eV (NBD-OCH3) and 3.01 eV (NBD-OC6H5).

CompoundsCyclic voltammetry (CV)
EOX(V vs. FC/FC+)ERed(V vs. FC/FC+)ERed0VEnergy levelsChemical reactivity descriptors
εHOMOCV (eV)εLUMOCV (eV)μ (eV)η (eV)WCV (eV)

Table 1.

Calculated HOMO and LUMO energy levels and extracted chemical reactivity descriptors for studied compounds.

For the three first compounds used as electrophile groups, they were oxidized and reduced at relatively distinct potentials. Indeed, the NBD-Cl is oxidized and reduced at higher potentials than methoxy and phenoxy substituted NBD.

For the amino-substitued NBD (Figure 4 and Table 1), the oxidation potential peak values were detected nearly 1.66 V, 1.58 V and 1.55 V for morpholine, piperidine and pyrrolidine substituted NBD, respectively. They correspond to HOMO energy levels at −5.87 eV, −5.79 eV and − 5.76 eV, respectively. However, the reduction peaks of the studied compounds were reversible and thus the ERed was replaced in the equation by ERed0 determined experimentally. Accordingly, the estimated LUMO energy levels are of about: −3.62 eV (NBD-Morph), −3.55 eV (NBD-Pip) and − 3.41 eV (NBD-Pyrr).

3.2 Photo-physical properties

3.2.1 UV-Vis optical absorption and emission (PL) analysis

The optical absorption spectra of the four substituted benzofurazans (NBDs) compounds are illustrated in Figure 5. The NBD-Cl, as starting material, has two distinct absorption bands (262 nm and 337 nm). The spectral characteristics of NBD-Cl were maintained and appeared as the same general features for amino-NBD derivatives. For the later systems, a new optical band in the range 478-488 nm was assigned to the formation of amino-NBD derivatives by SNAr reactions. In fact, typical three distinct absorption bands are presented for NBD substituted with pyrrolidine, piperidine or morpholine, similar to that previously reported NBD-amino derivatives [1, 2, 3]. Besides, the maximum absorption wavelengths of amino-NBD derivatives shifted more bathochromically than that of NBD-Cl.

Figure 5.

Experimental UV–vis optical absorption spectra of studied compounds as well as experimental PL spectrum of NBD-Morph.

According to assignments reported previously by Heberer and coworkers for NBD-amino derivatives [40] the short-wavelength band is associated to the aromatic benzofurazan compounds. The middle band at around 340 nm was attributed to a π → π* electronic transition. However, the band with lower energy in the range of 478-487 nm has ascribed to ICT arranged between the electron-donor and the nitro (NO2) electron-withdrawing groups within the molecule) (NBD-Cl, free band). Importantly, the maxima for the charge transfer band in these systems are detected in the visible region of the spectra.

In addition, similar profiles of optical absorption spectra for amino-NBD were obtained. Thus, the excitation corresponding to the lowest energy band of title compounds was used in order to investigate the excitation process from the ground (S0) to the first excited state (S1).

Experimentally, the 4,7-disubstitued benzofurazan compound, NBD-Morph was excited at the maximum absorption wavelength and show two distinct fluorescence bands (See Figure 5). The PL spectrum shows of two dominant features peaking at approximately 2.21 eV (559 nm) and 2.05 eV (603 nm), thus, a broad green to orange emission band has been identified for the tested NBD-Morph compound.

3.2.2 Time-resolved photoluminescence (TR-PL) decay kinetics

The performed time-resolved photoluminescence (TR-PL) decay kinetics is analyzed and illustrated by experiment (Figure 6). A particular attention is devoted to the thorough analysis of non-exponential decay kinetics. However, the TR-PL decay recorded at longer wavelength emission was relatively affected. It should be mentioned that lifetimes are inaccessible from the decay process by the use experimental equipment and could be related to an efficient intra-molecular charge transfer occurred in NBD-Morph molecule. The result was supported by the ICT characteristic optical band for a compound that had a lower intensity compared to those of analogue compounds (See Figure 5). The flat TR-PL spectrum means that the photo-carriers have a long lifetime (more than the nanosecond regime), leading to an important diffusion length.

Figure 6.

Time-resolved decay curves recorded at emission energy of 2.21 eV (560 nm) and 2.05 eV (603 nm) for NBD-Morph film.

3.3 Quantum chemistry computation results

Complete geometry optimizations of studied compounds were firstly carried out at DFT using B3LYP functional and the polarized 6-31 g+(d,p) basis set. The optimized molecular structures are illustrated in Figure 7. Particularly, the title amino-NBD compounds contain mainly three types of chemical bonds, namely O---H, N---H and N---O which are associated to an intra-molecular non covalent interactions.

Figure 7.

Molecular optimized structures with intra-molecular short contacts.

The computed geometrical parameters related to the short contacts at ground (S0) and excited (S1) states are given in Figure 7 and tabulated in Table 2. The S0 calculated short contact O---H was 2.39 Å (versus 2.38 Å, S1) in all systems, considerably smaller than the sum of the van der Waals radii (2.72 Å). The N---O short contact is fond 2.78 Å versus 2.74 Å, S1) significantly too smaller than the sum of the van der Waals radii (3.07 Å). In addition, the third short contact of N---H (smaller than the sum of the van der Waals radii (2.75 Å)) is allowed to vary from 2.11 Å -2.47 Å (S0) to 2.22 Å-2.53 Å (S1). The later (N---H) contacts are sensitive to amino substituted groups. As a result, these contacts have been identified as an important driving force on the stabilization of the studied co-planar molecular structures [41]. The major difference in the ground and excited state is the change in torsional angle ϕ and in the dipole moment. Calculations reveal that the existence of the amino groups significantly modifies the properties of their excited states. While the ground state is predicted to be planar, the excited-state geometry is twisted by about 2-7°. Due to conformational changes, the singlet excited-state dipole moment was found to be greater than ground-state dipole moment, unexpected for NBD-Pip (16.84 D (S1) versus 16.53 D (S0)). Importantly, the C-NO2 bond lengths indexed as (1) decreases by 0.4 Å upon excitation, due the nitro strong electron-withdrawing group effect.

Intra-molecular Short ContactsCovalent contactsDihedral AngleDipole Moment (D)
D1 O---H (Å)D2 N---O (Å)D3 N---H (Å)dC-N (1) (Å)dC-N (2) (Å)ϕ (°)
NBD-Pyrr2.39 (2.38)2.78 (2.74)2.47 (2.53)1.41 (1.37)1.33 (1.32)3.47
NBD-Pip2.39 (2.38)2.78 (2.74)2.11 (2.25)1.41 (1.37)1.34 (1.33)9.66 (16.76)16.84
NBD-Morph2.39 (2.38)2.78 (2.74)2.11 (2.22)1.41 (1.37)1.35 (1.33)5.16 (12.34)13.57

Table 2.

Calculated non covalent interactions for O---H N---O and N---H participating in short intra-molecular contacts within molecular at ground (S0) and excited (S1) states. The values presented in parentheses are obtained at S1.

Note that the band gap values estimated using the electrochemical method may be different from those calculated by optical or theoretical methods (see Table 3). This result could be explained by the fact that the oxidation/reduction process at the electrode are reactions generating species in ground states, whereas the optical electron transition leads to the formation of excited states. Moreover, the electro-generation of a radical cation or radical anion in solution involves other thermodynamic (salvation…) and kinetic effects. Consequently, it is expected that the peak potential will be lowered due to fast chemical reactions following the primarily electron transfer generating chemically reactive radical ions.

CompoundsEnergy levelsChemical reactivity descriptorsIP(eV)EA(eV)
εHOMOeVεLUMOeVμ (eV)η (eV)W (eV)
X = Cl−7.648−3.929−5.7881.8599.0097.4814.195
X = OCH3−6.988−3.561−5.2741.7138.1176.8063.830
X = OC6H5−6.985−3.684−5.3341.6508.6206.6294.041

Table 3.

Calculated HOMO and LUMO energy levels and extracted chemical reactivity descriptors for studied compounds using DFT//B3LYP/6-31+ g(d,p) in acetonitrile.

To better understand the optical responses, we used the time-dependent density functional theory (TD-DFT) to simulate the optical absorption spectra (See Figure 8). Their corresponding electronic transition assignments are given in Table 4. The maximum absorption wavelengths and the molecular extinction coefficients (ε) of the chloro- and amino-substituted NBD, in acetonitrile, are also given.

Figure 8.

Simulated optical absorption spectra of Cl-NBD and its related amio-substituted NBD by means of DFT//B3LYP/6-31+ g(d,p) level of theory, in acetonitrile.

Compoundsλmaxabs(nm)Osc. strengthMajor contribsMinor contribs
NBD-Cl3410.2729H → L (99%)
2550.0709H-4 → L (79%)H → L + 2 (10%)
H → L + 1 (6%)
λmaxExp: 262 nm, 337 nm
NBD-Pyrr2270.0181H → L + 3 (94%)H → L + 5 (3%)
3410.2781H → L + 1 (87%)H-1 → L (6%)
H → L (5%)
4380.4065H → L (94%)H → L + 1 (5%)
λmaxExp: 268 nm, 349 nm, 484 nm
NBD-Pip2300.0012H-7 → L (78%)
H-4 → L + 1 (13%)
H-9 → L (3%)
H-4 → L(2%)
3430.2083H-1 → L (30%)
H1 → L + 1 (59%)
H-1 → L + 1 (7%)
H → L (2%)
4520.4129H → L (94%)H → L + 1 (5%)
λmaxExp: 272 nm, 347 nm, 487 nm
NBD-Morph2370.0859H → L + 2 (64%)
H-3 → L + 1 (15%)
H-6 → L (12%)
H-6 → L + 1 (2%)
3410.2608H → L + 1 (89%)H-1 → L (5%)
H → L (4%)
4510.4209H → L (95%)H → L + 1 (4%)
λmaxExp: 274 nm, 340 nm, 478 nm

Table 4.

The vertical excited energies (nm) and their oscillator strengths (f) for the ground (S0 → S1) states of studied compounds by TD//B3LYP/6-31+ g(d,p) level of theory.

It is worthwhile mentioning that the extensive DFT and TD-DFT calculations show a good correlation between observed absorption and theoretical vertical excitation.

For the emission properties, it should be noted that the ZINDO semi-empirical quantum chemistry [42] was used to predict the emission spectra (See Figure 9). Table 4 lists the emission optical bands of the amino substituted benzofurazan and their assignments. The emission properties are mainly from LUMO and LUMO+1 to HOMO.

Figure 9.

Simulated emission spectra of amino-substituted NBD compounds (a,b,c). Superposition of emission spectra with experimental data for NBD-Morph is given (Figure 9D).

It is found that compounds emit at appreciably higher wavelengths in the range 620-640 nm. It should be noted also that bathochromic effect is due to the NBD moiety and the amino substituted groups. It should be noted that NBD-Morph present absorption at wavelength maximum of 478 nm and emit orange light at 559-603 nm. This provides a large apparent Stokes shift of 81 nm attributed to the ICT process (Table 5). The combined experimental and computed results suggest that the new NBD-Morph compound is still considered fluorescent. Then, considerable efforts have been undertaken to test the other analogue compounds.

CompoundsλmaxEm(nm)Osc. strengthMajor contribsMinor contribs
NBD-Pyrr4180.0498L + 1 → H (95%)
6230.6989L → H (97%)
NBD-Pip4180.0564L + 1 → H (95%)
6380.6806L → H (97%)
NBD-Morph4180.0552L + 1 → H (95%)
6310.6747L → H (97%)
λmaxExp: 559 nm, 603 nm

Table 5.

The vertical excited energies (nm) and their oscillator strengths (f) for the S1 → S0 states of studied compounds obtained by ZINDO method. The wavelengths excitation was ranging from 424 to 435 nm.

4-amino substituted NBD is recognized among to the broad family of intra-molecular charge transfer (ICT) complexes. The ICT characteristics are identified by the presence of both electron donor and withdrawing acceptor substituent moieties within the same molecule. Thus, the photo-initiated electron transfer from the donor to the acceptor sites yields two kinds positive and negative charges within separated functional parts of the molecule. Herein, we have extracted the atomic Mulliken charges from both ground- and excited-geometry structures. The results are illustrated in Figure 10. In fact the 4-amino substituted NBD belongs to the broad family of ICT molecules, with the amino group acting as an electron donor upon photo-excitation, and the nitro (NO2) group as an electron acceptor.

Figure 10.

Mulliken charge distribution of ground- and excited (bold values) states optimized structures of studied compounds.

As illustrated in Figure 11, we have presented the frontier molecular orbitals (FMOs) to obtain insight into the molecular structure and optoelectronic changes from the ground (S0) to excited (S1) States. We can observe that for all systems, the nitro and the 4-amino substituent groups significantly contribute to the optoelectronic properties. Particularly, the NBD-Morph differs from the other two compounds, where the oxygen atom did not contribute exclusively in excited state. This could explain the electron transfer from the nitrogen group of the electron donor moiety.

Figure 11.

Frontier molecular orbitals (FMOS) obtained at the ground- (S0) and excited- (S1) states of the optimized geometries of studied compounds.

3.4 Relationships of chemical structure and reactivity properties

It is relevant to note that good accuracy of computational predictions of optical properties is already attainable. Accordingly, we utilize this approach for accurate prediction of the global reactivity descriptors including chemical potential (μ), Chemical hardness (η) and electrophilicity index (ω) (see Table 5).

Parr’s electrophilicity ω values are calculated according to Eq. (3) [43, 44] based on μ and η which can be evaluated using Eqs. (4) and (5).


3.4.1 Correlation analysis

In a recent study, we showed that the second-order rate constants for the reactions of 7-X-4-nitrobenzofurazans 1 (1a, X = Cl, 1b, X = OC6H5 and 1c: X = OCH3) with secondary cyclic amines 2 (2a, morpholine, 2b, piperidine and 2c, pyrrolidine) in acetonitrile at 20°C [27, 29] can be described by the linear free energy relationship log k = sN(N+ E) [45, 46, 47, 48], and we determined the E parameters of these electrophiles 1a-c which are collected in Figure 12. We have now used these parameters E in order to elucidate the relationship between Mayr’s electrophilicities E [45] and Parr’s global electrophilicity index w values [49, 50]. As will be seen, satisfactory correlations between experimentally electrophilic reactivities E and other properties of these series of para-substituted nitrobenzofurazans 1a-c are found and discussed. On the other hand, we discuss how the structure–property relationships can be used to evaluate electrophilicity parameter E and Hammett constant σ values which are not directly accessible.

Figure 12.

Electrophilicity parameters E of 7-X-4-nitrobenzofurazans 1a-c. [27, 29].

3.4.2 Correlation between electrophilicity parameter (E) and global electrophilicity index (w)

Based on the data listed in Table 6, the plot of electrophilicity parameters E for these series of benzofurazans 1a-c against their global electrophilicity index w, determined in this work, was constructed. As can be seen in Figure 13, the relationship between E and w values is linear (r2 > 0.9799) and results in the following equation:


Table 6.

Electrophilicity parameter (E) and global electrophilicity index (w).

The E values were taken from Ref. [29].

The E value was taken from Ref. [27].

The E values estimated in this work.

The E values were taken from Ref. [27].

The w values calculated in this work.

Figure 13.

Correlation between experimentally determined electrophilicities (E) and global electrophilicity index (ω) calculated at DFT//B3LYP/6-31+ g(d,p) in acetonitrile for the benzofurazans 1a-c.


It is interesting to note that the relationships parameters E versus w have also been reported by many authors [51, 52, 53, 54, 55, 56]. The data in Table 7 are illustrative in this regard. Interestingly, the plot of Eq. (6) can be used to estimate the unknown electrophilicity parameter E values of other para-substituted benzofurazans. Using the w values calculated in the present work by DFT//B3LYP/6-31+ g(d,p), the E of 1d (X = morpholine, E = −18.01), 1e (X = piperidine, E = −18.52), and 1f (X = pyrrolidine, E = -18.78) have been obtained (The structures of title compounds are given in Figure 2). The detailed results are already listed in Table 6.

Table 7.

Correlations of the electrophilicity parameters E of some representative electrophiles with their global electrophilicity w values.

It should be noted that these predicted values of electrophilicity parameters E are smaller by 1.5 units than those reported values by Raissi and co-workers [27]. In all cases, the agreement between the E values of benzofurazans 1d-f estimated using DFT//B3LYP/6-31+ g(d,p) and previously reported with the same method of calculation, but recorded in water medium [27] is remarkably good [57, 58, 59, 60, 61]. It is noteworthy to mention that the observed deviations between results give lower accuracy owing to their dependence on solvent polarity and to the added diffuse basis functions, that have also been observed in numerous systems. Chamorro and co-workers, have also observed a satisfactory correlations between the electrophilicity parameter E of various benzhydrylium cations and their global electrophilicity index ω calculated at B3LYP/6- 31 g(d) and HF/6-31 g(d) levels of theory Model I (E = −41.60+ 3.57 ω+), Model II (E = −40.42+ 2.08 ω+), Model III (E = −38.66+ 2.44 ω+) and Model IV (E = −47.2+ 20.0 ω+) [62].

3.5 Effect of substituents X on electrophilicity: E vs. σ and w vs. σ correlations

Examination of the data in Tables 6 and 8 show that the electrophilicity parameter (E) or global electrophilicity index (w) of benzofurazans 1a-c appears to be significantly dependent on the electronic nature of the substituents X, i.e. the E or w decrease regularly from X = Cl to X = OCH3 leading to a nice Hammett-type plots [63, 64] (Figure 14), which are defined by the following Equations

Table 8.

Hammett constants (σ) values and dipole moment (μ, D) values.

The σ values were taken from Ref. [63, 64].

The σ values calculated by using Eq. (7) with E from Table 6.

Calculated by using Eq. (8) with w from Table 6.

The dipole moment μ values calculated in this work.

Figure 14.

Correlations of the electrophilicity parameter E and global electrophilicity index w versus the Hammett’s substituted σp values. The values of E and w are given in Table 6. The σp values were taken from Ref. [63, 64].


Correlations between E and σ parameters or E and w of various system have also been reported by many authors [65, 66, 67, 68, 69, 70, 71, 72]. Zenz and Mayr, have established that the E parameters for a series of trans-β-nitrostyrenes correlate well with σp values of their substituents (E = −13.95+ 2.08 σp (r2 = 0.9847)) [69]. Recently, Rammah and co-workers, have also observed a linear correlation between the electrophilicity parameters E of a series of 2-N-(4′-X-pheny1)-4,6-dinitrobenzotriazole 1-oxides and the Hammett’s substituent constants σ (E = −11.15+ 2.28 σ (r2 = 0.9878)). [65]

In order to test the applicability of Eq. (7) for the prediction of the Hammett’s constants σp of other groups, we have calculated the global electrophilicity index (w) of benzofurazans 1 g-j (1 g: X = HNOH, 1 h: X = HNNH2, 1i: X = HNH and 1j: X = HNCH3). The results are listed in Table 9 which also includes the σp values previously published [63, 64]. As well be seen, Table 2 clearly shows that the calculated σp values are in good agreement with the experimental data.

Table 9.

Comparison between calculated and experimental reported Hammett’s constants σp for the HNOH, HNNH2, HNH and HNCH3 groups.

The w values calculated in this work.

The dipole moment μ values calculated in this work.

The σpCal values calculated by using Eq. (8) with w from Table 5.

The σpExp values were taken from Ref. [63, 64].

Most importantly, the satisfactory correlations shown in Figure 14 can be employed to obtain more information about the unknown Hammett’s substituent constants (σ) of morpholine, piperidine and pyrrolidine groups.

Substitution of the E and w values for benzofurazans 1d-f into the correlation Eqs. (7) and (8), the Hammett’s σp values are thus obtained for the three leaving-groups. As can be seen in Table 8, the σp values estimated in the present work using work using Eqs. (7) or (8) are similar but relatively large compared with those previously reported (−0.47 < σp < −0.24) [27]. The major reason for the overestimation can, at least in part, be understood in terms of solvent polarity and also the used on the basis set with diffuse functions.

3.6 Correlation between dipole moment (μ) and Hammett’s substituent constants (σ)

As shown in Table 8 and Figure 15, the dipole moments (μ) calculated at DFT//B3LYP/6-31+g(d,p) for these series of benzofurazans 1a-f are highly sensitive to the electronic nature of the substituent X, e.g., it decreases from μ = 17.08 D to μ = 7.23 D as the substituent X changes from a strong electron-withdrawing group (X = Cl) to a strong electron-donating group (X = morpholine). The linear relationship as expressed by Eq. (9) relates the dipole moment μ for 1a-f directly to the Hammett’s substituent constant σ.

Figure 15.

Correlation between calculated dipole moment (μ) for the benzofurazans 1a-f and Hammett’s substituent constant σ calculated at DFT//B3LYP/6-31+ g(d,p) in acetonitrile.


On the other hand, comparing the electrophilicity parameter E and global electrophilicity index w with the dipole moment (μ), we found that the E and values w correlate well with the μ as shown in Figure 16. In addition, the high dipole moment values being related to the strong effect exerted by the electron-donating group X clearly indicate that our benzofurazans 1b-j can offer certain potential for nonlinear optical (NLO) applications.

Figure 16.

Correlation E vs. μ and w vs. μ for the benzofurazans 1a-f calculated at DFT//B3LYP/6-31+ g(d,p) in acetonitrile.

3.7 Correlation between global electrophilicity index wCV and Hammett’s substituent constants (σp)

To further extend the generality of the structure–property relationships to the global electrophilicity index wCV calculated from the Cyclic Voltammetry (CV) listed in Table 1, we examined correlations wCV vs. σp, wCV vs. E and wCV vs. w. As can be seen, the plotting of the wCV values determined above (see Table 1) versus the known Hammett’s constant σp values of substituent X gives rise to two separate linear correlations corresponding to a different behavior of Cl, OCH3, OC6H5 groups and morpholinyl, piperidinyl and pyrrolidinyl groups (Figure 17).

Figure 17.

Correlations of the global electrophilicity index wVC of the benzofurazans 1a-f versus the Hammett’s substituted σp values. The values of wCV are given in Table 1. The σp values were taken from Ref. [63, 64].

Similar behaviors have also been obtained when the global electrophilicity index wCV values were plotted against the electrophilicity parameter E or the global electrophilicity index w (Figures 18 and 19). The reason for this apparent disparity is not clear at this stage. Clearly, more investigate is necessary to explicate this behavior.

Figure 18.

Correlations between global electrophilicity index wCV and experimentally determined electrophilicities E for the benzofurazans 1a-c. the values of wCV are given in Table 1 and values of E are given in Table 6.

Figure 19.

Correlations of the global electrophilicity index wCV and global electrophilicity index ω calculated at DFT//B3LYP/6-31 + g(d,p) in acetonitrile for the benzofurazans 1a-f. the values of wCV are given in Table 1 and values of w are given in Table 6.


4. Conclusion

The reaction between NBD-Cl and morpholine, piperidine or pyrrolidine, as nucleophile group, proceeds normally by a SNAr attack affording 4.7-di-substituted benzofurazans. They were thoroughly analyzed electrochemically and optically by means of Cyclic Voltammetry (CV) and UV–Vis and time resolved photoluminescence (TR-PL), respectively.

Complementary study based on density function theory (DFT) and its extent TD-DFT were also conducted to unravel the electronic structure, the simulated optical spectra and to further understand the structure property relationships of the investigated compounds. The used procedure based theoretical calculations allows deducing the electronic parameters related to the HOMO and LUMO energy levels and their energy differences. These levels are accessible by measuring the peak potentials of the oxidation and reduction reactions. It is found that the measured energy levels are in close agreement with the values computed from DFT method. Specific nucleophile-electrophile interactions in SNAr reactions were expected to be governed by three main types of non-covalent bonds (O---H, N---O and N---H). Importantly, the compounds absorb visible light at longer wavelengths. In addition, the SNAr attack leads to intra-molecular charge transfer of formed compounds and thus provides a remarkably slow non-radiative decay from the excited state. Investigation of photo-physical properties revealed importance of intra-molecular twisting on excited (S1) states.

As verified by experimental and calculated results, the obtained molecules, by SNAr reactions, were modulated by intra-molecular non-covalent interactions that force the co-planarity and the rigidity of the structures yielding an efficient intra-molecular charge transfer (ICT). Thus, the nucleophile substituents acting as electron-donating groups in SNAr reaction afford to obtain materials with interesting nonlinear optical (NLO) response. In addition, the photo-optical data corroborated with TD-DFT approach were discussed in correlation with the structural architecture of each compound. We have shown the high efficiency fluorescent NBD-Morph compound that emits green to orange color at wavelength maximum of 559-603 nm.

The theoretical optimized parameters have been found to be in good agreement with the corresponding experimental data and results in the literature. Satisfactory linear correlation has been demonstrated between the parameters E and the global electrophilicity index w for the electrophilic reactivity at the C-X position of nitrobenzofurazans 1a-d. Accordingly, structure–property relationships were found to be able to evaluate the unknown electrophilicity parameter E of 7-X-4-nitrobenzofurazans 1d-f. Alternatively, the validity of Eqs. (7) and (8) has been satisfactorily verified by comparison of calculated and previously reported, in the literature, of the Hammett constant σp values of HNOH, HNNH2, HNH and HNCH3 groups.


Conflict of interest

The authors declare no conflict of interest.


  1. 1. F. Qian, C. Zhang, Y. Zhang, W. He, X. Gao, P. Hu, Z. Guo, J. Am. Chem. Soc. 131 (2009) 1460-1468.
  2. 2. T. Yamaguchi, M. Asanuma, S. Nakanishi, Y. Saito, M. Okazaki, K. Dodo, M. Sodeoka, Chem. Sci. 5 (2014) 1021-1029.
  3. 3. P. B. Ghosh, M.W.Whitehouse. Biochem. J. 108 (1968) 155-156.
  4. 4. P. B. Ghosh, M.W.Whitehouse. J. Med. Chem. 11(1968) 305-311.
  5. 5. D. Castagnolo, M. Pagano, M. Bernardini, M. Botta, Tetrahed. Lett. 53 (2012) 5008-5011.
  6. 6. C. Toriumi and K. Imai, Anal. Chem., 74 (2002) 2321-2327.
  7. 7. T. Kajiro, Y. Nakajima, T. Fukushima and K. Imai, Anal. Chem., 74 (2002) 4519-4525.
  8. 8. H.C. Ting, Y.T. Yang, C.H. Chen, J.H. Lee, J. H. Chang, C.I. Wu, T. L. Chiu, C.F. Lin, C. L. Chung, K. T. Wong, ChemSusChem 9 (2016) 1433-1441.
  9. 9. J. Mori, T. Kaino, Phys. Lett. A 127 (1988) 259.
  10. 10. L. Johnson, S. Lagerkvist, P. Lindroth, M. Ahnoff, K. Martinsson, Anal. Chem., 54 (1982) 939-942.
  11. 11. D. Dal Monte, E. Sandri, L. Di Nunno, S. Florio, P. Todesco, Chimica e Industria, 53 (1971) 940-942.
  12. 12. R.A. Manderville and E. Buncel, J. Chem. Soc. Perk. T. 2, (1993)1887-1894.
  13. 13. S. Uchiyama, T. Santa, T. Fukushima, H. Homma, K. Imai, J. Chem. Soc. Perk. T. 2, (1998) 2165-2173.
  14. 14. M. Bem, M. T. Caproiu, D. Stoicescu, T. Constantinescu, and A. T. Balaban, Cent. Eur. J. Chem., 2003, 1, 260.
  15. 15. I. Jamaoui, T. Boubaker, R. Goumont, Int. J. Chem. Kinet. 45 (2013) 152-160.
  16. 16. W.J. Wang, J.M. Xia, X. Hai, M.L. Chen, J.H. Wang, Environ. Sci. Nano, 4 (2017) 1037-1044.
  17. 17. R. Poulain, D. Horvath, B. Bonnet, C. Eckhoff, B. Chapelain, M.-C. Bodinier, and B. Deprez, J. Med. Chem., 44 (2001) 3378.
  18. 18. N. Lavignac, C. J. Allender, and K. R. Brain, Tetrahedron Lett., 45 (2004) 3625.
  19. 19. M. Bem, M. Vasilescu, M. T. Caproiu, C. Draghici, A. Beteringhe, T. Constantinescu, M. D. Banciu, and A. T. Balaban,Cent. Eur. J. Chem., 2 (2004) 672-685.
  20. 20. Jun-ichi Aihara, Reduced HOMO-LUMO Gap as an Index of Kinetic Stability for Polycyclic Aromatic Hydrocarbons, J. Phys. Chem. A 103, 37 (1999) 7487–7495.
  21. 21. Jun-ichi Aihara, Phys. Chem. Chem. Phys., 2000, 2, 3121-3125.
  22. 22. R. Zaier, S. Ayachi, Optik, 239 (2021) 166787.
  23. 23. R. Zaier, S. Hajaji, M. Kozaki, S. Ayachi, Opt. Mat. 91C (2019) 108-114.
  24. 24. A. J. Boulton, A. R. Katrizky, A. M. Hamid, J. Chem. Soc., 1967 (2005) 2007.
  25. 25. J. C. Hallé, M. Mokhtari, P. Soulié, M. J. Pouet, Can. J. Chem 75 (1997) 1240.
  26. 26. M. R. Crampton, L. M. Pearce, L. C. Rabbitt, J. Chem. Soc., Perkin trans 2 (2002) 257.
  27. 27. H. Raissi, H. Ayachi, F. Mahdhaoui, S. Ayachi, T. Boubaker, J. Mol. Str. 1224 (2021) 128843.
  28. 28. H. Ayachi, H. Raissi, F. Mahdhaoui, T. Boubaker, Int. J. Chem. Kinet. 52 (2020) 655-668.
  29. 29. H. Raissi, I. Jamaoui, R. Goumont, T. Boubaker, Int. J. Chem. Kinet. 49 (2017) 835-846.
  30. 30. F. Mahdhaoui, R. Zaier, N. Dhahri, S. Ayachi, T. Boubaker, Int. J. Chem. Kinet. 51 (2019) 249-257.
  31. 31. I. Messaoudi, I. Aribi, Z. Zouhour S. Ayachi, M. Othman, A.H. Said, J. Mol. Str. 1231 (2021) 129810.
  32. 32. A. Abbotto, L. Beverina, N. manfredi, G. A. Pagani, G. Archetti, H. G. Kuball, C. Wittenburg, J. heck, J. Holtman, Chem. Eur. J. 15 (2009) 6175-6185.
  33. 33. J. A. Davies, A. Elangovan, P. A. Sullivan, B. C. Olbricht, D. H. Bale, T.R. Ewy, C.M. Isborn, B.E. Eichinger, B.H. Robinson, P.J. Reid, X. Li, L. R. Dalton, J. Am. Chem. Soc. 130 (2008) 10565-10575.
  34. 34. Z. Zaaboub, F. Hassen, H. Maaref, Sol. Stat. Comm. 314–315 (2020) 113913.
  35. 35. A. D. Becke, J. Chem. Phys. 98 (1993) 5648-5652.
  36. 36. C. Lee, W. Yang, R. G. Parr, Phys. Rev. 37 (1988), 785-789.
  37. 37. M. J. Frisch., et al. Gaussian 09, Revision A.01, Gaussian: Wallingford, CT (2009).
  38. 38. Dennington, R.D., Keith, T.A. and Millam, J.M. (2008) GaussView 5.0.8, Gaussian.
  39. 39. J. L. Bredas, R. Silbey, D. S. Boudreaux, and R. R. Chance, J. Am. Chem. Soc 105 (22) (1983) 6555-6559.
  40. 40. Heberer, H.; Matschiner, H., J. Prakt. Chem., 328 (1986) 261-274.
  41. 41. S. Ayachi, S. Ghomrasni, M. Bouachrine, M. Hamidi, K. Alimi, J. Mol. Str. 1036, (2013) 7-18.
  42. 42. S. Ayachi, S. Ghomrasni, K. Alimi, J. App. Polym.Sci. 123 (2012) 2684-2696.
  43. 43. R. G. Parr, L. V. Szentpaly, S. Liu, J Am Chem Soc. 1999, 121, 1922-1924.
  44. 44. Parr, R. G, Yang, W. Density Functional Theory of Atoms and Molecules; Oxford University Press: New York 1989.
  45. 45. H. Mayr, M. Patz, Angew. Chem. Int. Ed. Engl. 33 (1994) 938-957.(b) H. Mayr, B. Kempf, A. R. Ofial, Acc. Chem. Res. 36 (2003) 66-77.
  46. 46. S. Ben Salah, F. Necibi, R. Goumont, T. Boubaker. ChemistrySelect. 5 (2020) 7648-7657.
  47. 47. S. Souissi, W. Gabsi, A. Echaieb, J. C. Hierso, P. Fleurat-Lessard, T. Boubaker. RSC Advances. 10 (2020) 28635-28643.
  48. 48. Z. Li, R. J. Mayer, A. R. Ofial, and H. Mayr. J Am Chem Soc. 142 (2020) 18, 8383-8402.
  49. 49. R. G. Parr, R. G. Pearson, J. Am. Chem. Soc. 105 (1983) 7512-7516.
  50. 50. R. G. Parr, L. V. Szentpaly, S. Liu. J. Am. Chem. Soc. 121 (1999) 1922-1924.
  51. 51. D. S. Allgäuer, H. Jangra, H. Asahara, Z. Li, Q. Chen, H. Zipse, A. R. Ofial, H. Mayr. J. Am. Chem. Soc., 139 (2017) 13318–13329.
  52. 52. Z. Li, H. Jangra, Q. Chen, P. Mayer, A. R. Ofial, H. Zipse, H. Mayr. J. Am. Chem. Soc. 140 (2018) 5500−5515.
  53. 53. W. Gabsi, T. Boubaker, R. Goumont. Int. J. Chem. Kinet. 48 (2016) 266-273.
  54. 54. P. Pérez, A. Labbé. J. Org. Chem. 67 (2002) 4747-4752.
  55. 55. L. R. Domingo, P. Pérez, R. Contreras. L. R. Domingo, P. Pérez, R. Contreras. Tetrahed. 60 (2004) 6585-6591.
  56. 56. D. S. Allgäuer, H. Jangra, H. Asahara, Z. Li, Q. Chen, H. Zipse, A. R. Ofial, H. Mayr. J. Am. Chem. Soc. 139 (2017) 13318-13329.
  57. 57. Q. Chen, P. Mayer, H. Mayr. Angew Chem Int Ed 55 (2016) 12644-12667.
  58. 58. R. J. Mayer, P. W. A. Allihn, N. Hampel, P. Mayer, S. A. Sieber, A. R. Ofia. Chem. Sci. 12 (2021) 4850-4865.
  59. 59. S. Ben Salah, T. Boubaker, R. Goumont. Int J Chem Kinet 49 (2017) 576-583.
  60. 60. E. Follet, P. Mayer, H. Mayr, Eur J Org Chem. 2016, 4050-4058.
  61. 61. N. Dhahri, T. Boubaker, R. Goumont. J Phys Org Chem. , 27 (2014) 484.
  62. 62. E. Chamorro, M. Duque-Noreña, P. Pérez. J Mol Struct. 896 (2009) 73-79.
  63. 63. C. Hansch, A. Leo, R. W. Taft. Chem. Rev. 91 (1991) 185-195.
  64. 64. Hammett L. P. Reaction Rates, Equilibria and Mechanism. Mc Graw-Hill, New York, 2nd edn, 1970.
  65. 65. M. Rammah, F. Mahdhaoui, W. Gabssi, T. Boubaker. ChemistrySelect. 6 (2021) 4424-4431.
  66. 66. O. Kaumanns, H. Mayr. J Org Chem. 73 (2008) 2738-2745.
  67. 67. H. Asahara, H. Mayr. Chem Asian J. 7 (2012) 1401, 1407.
  68. 68. K. Azaiez, N. Dhahri, T. Boubaker. Int J Chem Kinet. 49 (2017) 847-858.
  69. 69. I. Zenz, H. Mayr. J Org Chem. 76 (2011) 9370-937.
  70. 70. J. Ammer, C. Nolte, H. Mayr. J Am Chem Soc. 134 (2012) 13902-13911.
  71. 71. K. Troshin, C. Schindele, H. Mayr. J Org Chem. 76 (2011) 9391-9408.
  72. 72. D. Richter, N. Hampel, T. Singer, A. R. Ofial, H. Mayr. Eur J Org Chem. 2009, 3203-3211.

Written By

Hanen Raissi, Imen Chérif, Hajer Ayachi, Ayoub Haj Said, Fredj Hassen, Sahbi Ayachi and Taoufik Boubaker

Submitted: 31 May 2021 Reviewed: 05 July 2021 Published: 23 November 2021