F. Aydın and N.B. Arslan
Journal of Molecular Structure 1232 (2021) 129976
part in chemical reaction and give information about the reactivity
or stability of specific regions of the molecule. The frontier molec-
ular orbitals of the title compound were calculated by using the
B3LYP/6–311G(d,p) level of theory.
The distribution of HOMO, HOMO-1 and LUMO, LUMO+1 ener-
gies for the compound are shown in Fig. 3.
The ionization potential (I) and electron affinity (A) can be ex-
pressed as follows in terms of HOMO and LUMO orbital energy
I = −EHOMO and A = −ELU MO
In addition, theoretical reactivities such as chemical hardness
(η), chemical softness (ζ), nucleophilicity (ε), chemical potential
(μ), electrophilicity index (ω), electronegativity (χ) and energy
band gap (ꢀE) were calculated by using the following defined for-
Scheme 3. 1H and 13 C NMR chemical shifts of the title compound.
Table 5
orbital energies, HOMO-LUMO energy gap, ionization potential (I),
Electron affinity (A), chemical potential (μ), chemical hardness (η),
absolute softness (ζ), absolute electronegativity (χ), chemical poten-
tial (μ), and electrophilicity index (ω) of title compound for ground
state geometries in the gas phase at the B3LYP/6–311G(d,p) level.
Based on the results in Table 6, the HOMO and LUMO energies
for the compound are −6.705 and −2.698 eV, respectively, and the
energy gap value of title compound is found to be 4.007 eV. A
small gap implies low stability and a large gap implies high sta-
bility. A molecule with a large HOMO-LUMO gap is less polariz-
ability and is generally associated with a low chemical reactivity.
The chemical hardness (ƞ), softness (σ), and chemical potential
(μ) of the title compound were calcalated to be 2.003, 0.499 and
−4.702 eV. It should be noted that the electro-negativity provides
information about the electron accepting capability of a molecule
while nucleophilicity is a measure of Lewis basicity [34]. Calculated
electro-negativity (χ) and nucleophilicity (ε) values for this com-
pound are 4.702 and 0.181 eV, respectively. Hence, these values in-
dicate that the title compound is a chemically hard, high kinetic
stable, and less reactive molecule.
Parameters B3LYP/6–311G(d,p)
(eV)
Kcal/mol
Electronic Energy (a.u.)
EHOMO
−2440,521
- 6.705
- 2.698
+6.705
+2.698
+4.007
2.003
154.61
62.21
154.61
62.21
92.40
46.18
11.50
4.17
ELUMO
Ionization energy, I= -EHOMO
Electron affinity, A= -ELUMO
Energy band gap, [ꢀE=EHOMO-ELUMO
Chemical hardness, η = (I-A) / 2
Chemical softness, ζ = 1/2η
Nucleophilicity, ε = 1/ω
Proton affinity Pı
]
0.499
0.181
−4.702
−4.007
5.518
108.42
92.40
127.24
108.42
Chemical potential, μ = - (I + A) / 2
Electrophilicity index, ω = μ²/2 η
Electronegativity, χ = (I + A) / 2
Dipol moment (Debye)
4.702
7.891
3.5. Molecular Electrostatic Potential (MEP)
The structure of the compound was confirmed by the presence
of eight peaks in the 13C NMR spectrum. In the 13C NMR spec-
trum, the methanetriyl carbon connected with diamine (C1) has
a sp3 hybrid, but it was located at δ 100.98 ppm because of the
geminal diamine and trichloromethyl group. The chemical shift of
carbon, C1 was calculated as δ 134.97 ppm (in gas phase), 135.18
and 135.41 ppm (in CDCl3 and DMSO–d6 solvents), respectively.
The 13C NMR signal of the trichloromethyl group (C2) appeared
at δ 73.93 ppm due to effect of the attached electron-withdrawing
chloride atoms. Likewise, the calculated chemical shift values for
C2 atom, known carbon of trichloromethyl were calculated at δ
149.39 ppm (in gas phase) and 77.12, 77.15 ppm (in CDCl3 and
DMSO–d6 solvents), respectively. The carbon atoms in groups are
shielded by their own three chlorine atoms and therefore they give
a resonance signal in high field [31].
The molecular electrostatic potential (MEP) surface analysis is
an important property that can be successfully used to explain the
electrostatic interactions. Furthermore, MEP is useful for mapping
sites for nucleophilic reactions and electrophilic attacks by show-
ing electron density [35]. To predict the reactive sites of the in-
vestigated title compound, MEP was calculated using DFT at the
B3LYP/6–311G(d,p) basis set in the gas phase, as shown in Fig. 4.
In the color scheme adopted, red indicates an electron-rich region
with a partially negative charge and blue an electron-deficient re-
gion with partially positive charge, light blue indicates a slightly
electron-deficient, yellow a slightly electron-rich region and green
a neutral region [36]. Negative regions of the MEP (red color)
are related to electrophilic reactivity, while positive regions (blue
color) explain the nucleophilic attacks.
Twelve aromatic carbons were observed in the region of
δ 142.13–114.14 ppm as expected. The highest chemical shifts
142.13 ppm were belonging to the amino-bound carbons C3 and C9
as doublets. These carbons were found as δ 134.97 (in gas phase),
135.18 and 135.41 ppm (in CDCl3 and DMSO–d6 solvents) by theo-
retical calculations, respectively.
In Fig. 4, it appears that the negative region is mainly around
the nitro groups (red) [37]. Amine sites, which could be easily at-
tributed to the lone pair electrons on nitrogen atoms are electron-
rich and hence could be reactive atomic sites. Positive regions are
mainly distributed around some hydrogen atoms (blue). Besides,
the most positive region (0.351 a.u.) is the carbon (C2) connected
with chloride atoms. The positive region (blue) is over the C1-H
group, while the delocalization of electrons from N1 and N2 to
the nitrophenyl rings is confirmed by the yellowish color. In addi-
tion to these, MEP is related to the electron density on atoms in a
molecule and is a very useful descriptor in understanding sites for
hydrogen bonding interactions [38]. Regions of positive and nega-
tive MEP on the surface of hydrogen bond donors and acceptors
are influenced by the formation of intramolecular contacts with
the molecules. In the compound, the C6 and C8 carbons are meta
position atoms with respect to nitro group, and ortho and para po-
sition atoms with respect to amino group. The C6 and C8 carbons
became less positive due to the electron withdrawing of the nitro
Other aromatic carbons were observed at 136.60, 133.77, 127.44,
118.48 and 114.14 ppm in the spectrum (Scheme 3 and Fig. S4).
These carbons were calculated similarly to the experimental val-
ues (Table 5). Besides the change in chemical shift values of car-
bon atoms depending on the substrate concentration and relax-
ation time, approaches in the theoretical calculation can also be
different, so the correlation is not sometimes perfect.
3.4. Frontier Molecular Orbitals (FMOs) analysis and global reactivity
The frontier orbitals, HOMO and LUMO are the most important
orbitals in a molecule. These orbitals are the main orbitals taking
6