Journal of Molecular Structure (Theochem) 468 (1999) 209–221

Density functional studies of aniline and substituted anilines

Mariana E. Vaschetto*, Bernardo A. Retamal, Andrew P. Monkman
OEM, Department of Physics, University of Durham, Durham DH1 3LE, UK
Received 22 October 1998; received in revised form 2 December 1998; accepted 4 December 1998

Abstract
In this work, we present an exhaustive study of the structural as well as the electronic properties of aniline obtained using
several DFT based methods. The geometry of aniline obtained using the hybrid DFT methods (B3LYP, B3P86 and B3PW91)
shows better agreement with the experimental values, than the methods containing the gradient corrected exchange functional
(Becke 88). The new feature of the hybrid methods is that the exact Hartree–Fock exchange functional was taken in account.
Negative highest occupied and lowest unoccupied Kohn–Sham orbitals (we shall call them as HOMO, LUMO) seem to provide
a fair estimation of the experimental IPs and EAs values (particularly those obtained by using the gradient correct exchange
functional Becke 88). IPs and EAs from DFT calculations present systematic errors are observed owing to the HOMO and
LUMO, which do not correspond formally to any molecular orbital.We, also, study in a systematic way the influence of a wide
range of s – p electron donor-withdrawing substituent groups on the structural and electronic properties of aniline. In this kind
of systems, the position of the substituent group in the benzene ring as well as its electron donor–acceptor capabilities play a
very important role on the molecular structural and electronic properties. We also find that the intramolecular hydrogen bonds
are responsible for many particular experimental features. q 1999 Elsevier Science B.V. All rights reserved.
Keywords: Aniline; Polyaniline and derivatives; Density functional calculations; DFT; Hydrogen bond

1. Introduction microwave spectroscopy [8,9] and in the solid state by

X-ray crystallography [10]. Moreover, several theore-
Aromatic amines are very important in biological tical studies including semiempirical [11–13] and ab
as well as in materials science [1,2]. Particularly, initio calculations [6,14–16] were carried out to deter-
aniline was studied because of its technological appli- mine the structural and electronic properties of
cation in a vast amount of industrial processes [3,4]. aniline. However, the conformation of the amine
In fact, it was used as a precursor for the synthesis of group and its degree of tilting related to the phenylene
dye and sensitizer molecules, for the synthesis of plane are still a matter of study. This large amount of
molecules having potential large non-linear optical previous data available on aniline provides an ideal
responses, etc [5,6]. basis for using this compound as a model compound
Aniline is the classic molecular model for aromatic to test the accuracy of the density functional theory
amines [7]. Its structure is known in the gas phase by (DFT) based-on methods. Subsequently, the most
precise algorithms can be applied to a family of
anilines in order to generalize the trends of the physi-
* Corresponding author. Tel.: 144-191-374-2159; fax: 144-
191-374-3749.
cochemical properties.
E-mail address: [email protected] (M.E. Vaschetto) The inclusion of a substituent group in the aniline
0166-1280/99/$ - see front matter q 1999 Elsevier Science B.V. All rights reserved.
PII: S0166-128 0(98)00624-1
210 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221

leads to the variation of the charge distribution in the the treatment of electronic excitations has emerged.
molecule, and consequently this greatly affects the Ahlrichs and Bauernschmitt [26] consider the time
structural as well as the electronic parameters. Parti- dependent density functional formalism (TDDFT)
cularly, it is expected that large modifications of the within the adiabatic approximation. They consider
ionization potential (IP), electronic affinity (EA) and that the excitation energies and the oscillator strengths
excitation energies should be found. Determination of are given as poles and residues of the response func-
the excitation energies of finite and infinite systems tion as a function of frequency (v ). Results given by
have been a subject of many research during the last this method show that, although there is a constant
50 years [17,18]. shifting of the values obtained for the excitation ener-
In principle, there are several ways to calculate the gies, the TDDFT is a very promising and new
excitation energies. The first, and the simplest one, approach for studying electronic excitations [27].
involves the difference between the highest occupied For a detailed description of this method the reader
molecular orbital (HOMO) and the lowest unoccupied is referred to Refs. [26] and [28].
molecular orbital (LUMO) of a neutral system. This In this work, we present the excitation energies of
form corresponds to the frozen orbital approximation, aniline and substitute anilines using the simplest and
as the ground state properties are used to calculate less computationally expensive methods. The most
excitation values. This method is very practical, parti- advanced way of calculating the excitation energies
cularly for calculating large systems; however, the use will be considered in future work.
of DFT one-electron energy differences to estimate This work is mainly divided in two parts. The first
excitation energies is still a matter of debate [19]. is a systematic analysis of the DFT methods, by apply-
Rigorously, the density functional methods (which ing them to analyze the geometry of aniline (model
are based on Hohemberg and Kohn theorem [20]) compound) and its electronic parameters, such as, the
are designed to yield total energies, for example, as IP, the EA, and the first electronic excitation. We shall
a function of structure. However, the orbitals in this compare the results obtained by using the gradient-
case (Kohn–Sham orbitals [21]) are not electronic corrected exchange and the hybrid exchange func-
orbitals but mathematical arrangements and their tional Becke three parameters. We will show that
eigenvalues are not directly related to the electronic the inclusion of Hartree–Fock exchange into the
excitation energies. Nevertheless, it is often a good DFT calculation, to give the so-called hybrid meth-
approximation to neglect these formal inconsistencies ods, leads to a significant improvement of the
and consider them as electronic orbitals and orbital predicted structures.
energies, as they are in the HF framework [22]. In the second part of this work, we will apply the
The second way to obtain the excitation energies DFT-hybrid scheme to calculate the structural para-
involves the difference between the IP and the EA. In meters as well as the electronic properties of aniline
this case, accurate IPs and EAs are obtained as the and aniline derivatives. In fact, we will present a
total energy difference between the system with N21 detailed study of anilines substituted with a wide
or N11 electrons. This form of calculation is very range of s – p electron-donor to electron-withdrawing
efficient, but in some cases can become very difficult groups (–NH2, –OH, –OCH3, –CH3, –CO2H, –Cl, –
because of the inclusion of open shell considerations CF3, –CN, –NO2) [29–32]. The purpose of this study
in the calculations. is to show in a systematic manner the variation of the
The third option is computationally more expensive geometrical parameters depending on the substituent
and involves the rigorous inclusion of electronic considered. Also, the electronic consequences of
correlation. This requires either extended multiple including a substituent on position ortho, meta or
reference configuration interaction (MR–CI)[23], or para in aniline will be presented. Moreover, we will
complete active space plus perturbation theory to try to connect our calculated results with the experi-
second order (CAS–PT2)[24] type calculations. mental parameters, such as pKb in order to determine
Equation of motion coupled cluster (EOM–CC) treat- the predictive capability of the theoretical method
ment achieves high accurate results as well [25]. employed.
Finally, in the last few years a new development for In Section 2, we summarize the theoretical methods
 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221 211

Table 1
Structural parameters of aniline calculated at several theoretical levels. The basis set used is 6-31G*. Reported bond lengths are in Å and the
angles are in degrees. See Fig. 1 for more details

RHF-AM1 B3LYP B3P86 B3PW91 BLYP BP86 BPW91 Expt. a

r(N1–C1) 1.400 1.400 1.392 1.394 1.411 1.404 1.403 1.402

r(C1–C2) 1.415 1.405 1.402 1.404 1.416 1.414 1.413 1.397
r(C2–C3) 1.390 1.393 1.390 1.391 1.403 1.401 1.399 1.394
r(C3–C4) 1.394 1.397 1.393 1.394 1.407 1.405 1.403 1.396
a (N1–C1–C2) 120.7 120.7 120.7 120.7 120.7 120.7 120.7 120.3
a (C6–C1–C2) 118.5 118.5 118.5 118.5 118.6 118.5 118.5 119.26
a (C1–C2–C3) 120.3 120.5 120.5 120.5 120.4 120.5 120.5 120.7
a (C2–C3–C4) 120.7 120.8 120.8 120.8 120.8 120.8 120.8 120.4
a (C3–C4–C5) 119.5 118.9 118.8 118.8 118.9 118.9 118.9 118.55
g 1b 25.8 26.5 25.7 25.6 27.1 26.6 26.5
g 2c 158.2 156.5 157.1 157.1 156.1 156.5 156.6
ud 4.6 2.9 2.7 2.7 3.3 3.2 3.1
de 45.4 40.3 40.2 41.0 40.2 40.0 37.2
a
From Ref. [8].
b
Torsion angle defined by H1–N1–C1–C6.
c
Torsion angle defined by H2–N1–C1–C6.
d
Angle defined by the plane of the benzene ring and the nitrogen.
e
Angle defined by the plane of the amine group and the plane of the benzene ring (see Fig. 1).

employed to carry out these calculations. In Section expressed as:

3.1 we present the theoretical results of the aniline,
meanwhile, Section 3.2, we introduce the theoretical A × ExSlater 1 12A† × ExHF 1B × ExBecke 1C
results of substituted anilines. Finally, conclusions are
provided in Section 4. × EcLYP 1EcVWN :

B3P86 and B3PW91 are defined analogously

2. Computational methods and differ only in the use of the respective corre-
lation functional Ec. The values of A, B, and C are
The geometry of aniline was completely optimized 0.80, 0.72 and 0.81 as determined by Becke. The
using the DFT schemes described later. No geometri- main difference from the other exchange func-
cal constraints were imposed to obtain the minimum tionals lies in the incorporation of the exact
energy structure of aniline. Optimized geometries Hartree–Fock exchange with a weight of 20%.
were also determined using the semiempirical The calculations were carried out using 6-31G*
Hartree–Fock AM1 [33] in order to compare the basis set, as well as using the Stephens–Basch–
results obtained using the DFT methods. Krauss ECPs with split valence basis set augmen-
The methods employed are based on DFT and we ted with polarization function (CEP-31G*). The
considered various combinations of exchange-corre- ‘fine grid’ option was used to avoid numerical
lation functionals. The exchange functionals used are: inaccuracies in the energy calculations. Spin
Becke 88 (B) [34], and the Becke three parameters restricted theory is used.
hybrid (B3) [35]. Meanwhile, the correlation func- Substituted anilines were fully optimized using
tionals used are: Lee–Yang–Parr (LYP) [36], Perdew the restricted version of B3P86 and B3PW91. In
(P86) [37], and Perdew and Wang (PW91) [38,39]. these cases, only the basis set 6-31G* were
B3LYP indicates that Becke’s three-parameter hybrid employed. All the calculations were performed
functional is used in combination with the LYP corre- using Gaussian 94 [40] and the Cerius 3 package
lation functional, etc. This functional can be [41].
212 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221

3. Results and discussion

3.1. Theoretical calculations of aniline

Table 1 depicts the calculated and experimental

structural parameters of aniline. In order to compare
the accuracy of several theoretical methods, we opti-
mized the geometry of aniline using manifold of DFT
methods. Earlier, Johnson et al. [42] and Mannfors et
al. [43] found that all functionals containing the gradi-
ent corrected exchange (Becke 88) consistently
produces bond lengths between heavy atoms (i.e. C,
N, etc) which are too long for about 0.01–0.02 Å. The
inclusion of a Hartree–Fock component in the hybrid
methods that consider the Becke three parameters Fig. 1. Schematic representation of aniline molecule.
exchange functionals give better agreement with the
experimentally determined bond lengths. This by the amino group and the plane of the ring was
improvement in the prediction of the bond lengths 37.58 ^ 28 in gas phase (see Fig. 1). Unfortunately,
using hybrid methods is observed in our aniline calcu- this estimation was made considering that the
lations (see Table 1). The explanation of this is nitrogen atom lies in the same plane as the ring.
obvious and is related to the fact that at the However, X-ray measurements [10], previous ab
Hartree–Fock level the bond lengths are calculated initio calculations [7,14–16] and results reported
too short, while at the DFT level they are overesti- here predict that the nitrogen atom is out of the
mated. Related to the employed basis set, when the ring plane for about 28–38. In fact, this value is in
CEP-31G* is used an important overestimation of the good agreement with that determined by Fukuyo
bond lengths (around 0.01–0.02 Å) is observed et al. [9] by X-ray crystallography. In this case,
compared to the experimental values. Therefore, the nitrogen atoms are out of the plane by about
here we shall present some selected results obtained 0.113 and 0.122 Å, which yields to a d angle of
with the 6-31G* basis set, which gives structural 378 ^ 48 and 388 ^ 48, very close to that deter-
values in quite good agreement with the experimental mined in the gas phase and also very close to the
ones [8–10]. value obtained by DFT calculations. Our predicted
One of the most remarkable features of these calcu- d angle, around 408, is very close to the experi-
lations is the accurate prediction of the degree of non- mental value. This angle does not show any
planarity of the amine group. The experimental remarkable variation with the used combination
results, obtained from microwave spectroscopy of exchange-correlation functionals. However, we
[8,9], show that the angle d between the plane defined should mention here that when the CEP-31G*

Table 2
Energy of the frontier orbitals of aniline calculated at several theoretical levels. The used basis set is 6-31G*. All values are in eV.

RHF-AM1 B3LYP B3P86 B3PW91 BLYP BP86 BPW91 Expt.

HOMO 28.523 25.391 26.025 25.467 24.360 24.661 24.572 7.72 a

LUMO 0.639 0.251 20.366 0.185 20.344 20.632 20.527
DE(LUMO-HOMO) 9.162 5.642 5.629 5.652 4.016 4.029 4.045 4.4 b
Eanion c2Ecation c 5.111 6.360 4.720 4.4
a
From Ref. [59].
b
From Ref. [60].
c
Total energy
 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221 213

basis set is used an increase in the d angle of 10% are found. One of the most probable explanations
about 18–28 is found, which is, actually, a devia- of this deviation of the excitation energies can be
tion from the experimental value. found in the Kohn–Sham theoretical formulation.
A selected group of calculated electronic para- Despite this, this form of calculation of the IP and
meters of aniline determined by using several differ- the excitation energy is very controversial owing to
ent methods are given in Table 2. Owing to the the consideration of the Kohn–Sham orbitals as elec-
analogy between the highest occupied and the lowest tronic orbitals, as was mentioned earlier, we can
unoccupied Kohn–Sham orbitals and the frontier consider it as a good first approximation. Table 2
molecular orbitals, we name them as HOMO and shows some selected values of excitation energies
LUMO respectively. According to the Koopmans calculated as the difference between the total energy
theorem, it is possible to obtain IP as the absolute of the cation and the anion of the aniline. Although
value of the HOMO [44]. It is well known, that there are some improvements to the values, the
Hartree–Fock methods are characterised by their computational cost is, at least, three times larger
accurate predictions of the IP [31], while methods than in the first case (direct difference between
based on the DFT scheme largely underestimate the LUMO and HOMO in the neutral state). The advan-
IP [17,18]. However, better values of IP could be tage of this method is, however that it accounted for
achieved depending on the exchange-correlation any structural changes suffered by the system when
functional combination used. In fact, the IP values the oxidation state changes.
obtained using the hybrid methods are better than In summary, we have shown that the DFT calcula-
the IPs achieved when using the gradient corrected tions that employ the Becke three parameters
exchange functional, Becke 88. Also, the correlation exchange functional reproduce accurately the geome-
functional plays an important role in the IP calcula- try of aniline compared to the methods based on the
tions. The B3P86 method shows the best IP value Becke 88 (non-hybrid) exchange functional. Also, the
compared to the experimental one than the B3LYP LYP correlation functional appears to give the most
and the B3PW91 methods. accurate N-ring distance. Therefore, the B3LYP
The EA of a system is related to its LUMO. combination of exchange correlation functionals
Hartree–Fock methods show a large overestimation seems to be most precise in order to estimate struc-
of the LUMO. However, the DFT methods give more tural properties of aniline. However, the evaluation of
accurate EA than the Hartree–Fock methods. Here, DFT based methods for predicting the electronic
again the value predicted using the Perdew 86 corre- properties is not so direct. Considering the adiabatic
lation functional gives the most accurate EA [22]. approximation, Becke three parameters exchange
Previous studies [18,45] have shown that the IP and functional gives the closest IP values, however the
EA calculated using the DFT have almost the same most accurate excitation energies are obtained by
error, which produces a shift of the frontier orbitals using the Becke 88 functional. In these cases, the
through less negative energies. These systematic Perdew 86 correlation functional gives better IP
errors of the IP and EA values compensate when values. This relative improvement of the excitation
calculating excitation energies, giving a good approx- energy values obtained by using BP86, could be asso-
imation to this parameter. In Table 2, it can be seen ciated to an unrealistic structure predicted by this
that the methods based on the DFT give excitation method (see Table 1). Therefore, the better electronic
energies that are in relatively good agreement with properties predicted by BP86 could be an artifact
the experimental value. Several authors [18,22,45] originating from the structural distortions arising
have claimed those hybrid methods, particularly from this method. Consequently, based on this
those that included P86 as the correlation functional previous analysis we have selected the B3LYP
give accurate IP and band gap values, while the LSDA method, which provides reliable structural properties
method give better EA’s. In our case, the BP86 func- that are critical when analyzing a family of
tional combination is superior to the rest of methods compounds, particularly conjugated systems. As our
for calculating excitation energies, while hybrid meth- main aim is to compare the physicochemical para-
ods overestimate it. In all cases, errors larger than meters of a family of anilines, and although the
214 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221

Table 3 B3LYP shows some weakness for predicting electro-

Selected structural parameters of substituted anilines, calculated nic properties, trends predicted by this method will be
using B3LYP (basis set 6-31G*). Distances are in Å, angles in
degrees
qualitatively correct. Following these results, we used
B3LYP in order to predict structural and electronic
Substituent Position r N1–C1 g1 g2 d u effects of the substituents covalentely added to an
aniline molecule.
NH2 ortho 1.411 39.4 165.3 43.3 3.8
meta 1.401 26.3 156.1 2.9
para 1.411 28.8 153.5 43.8 3.2 3.2. Theoretical calculations of aniline derivatives
OH ortho 1.398 22.4 153.0 41.1 2.8
meta 1.400 27.7 157.8 2.9 The interaction of a molecule of aniline with a
para 1.408 28.1 154.9 3.1 substituent group is of great importance in determin-
CH3 ortho 1.403 33.0 162.5 40.8 2.8 ing its structural and electronic properties [31,32,46–
meta 1.401 26.9 156.6 40.5 3.0
para 1.403 27.1 156.0 41.1 3.1
48]. As aniline and polyaniline are very relevant in the
OCH3 ortho 1.397 22.3 152.9 41.1 2.9 field of biological and materials science, we shall try
meta 1.399 25.3 155.5 40.2 2.9 to elucidate the influence of several substituent groups
para 1.408 27.2 154.1 42.9 3.1 on their physicochemical properties [49]. In this
Cl ortho 1.389 23.0 157.3 37.0 3.0 section, three aspects about substitution in aniline
meta 1.396 25.2 157.2 39.0 2.8
para 1.397 25.9 157.0 2.9
are investigated in detail. The first one is the influence
CN ortho 1.378 17.7 159.7 30.3 2.8 of these substituents on the structural conformation of
meta 1.393 24.0 157.6 37.1 2.8 aniline, including the flattening of the amine group.
para 1.385 22.6 160.0 33.8 2.6 Secondly, we focus our discussion on the electronic
CF3 ortho 1.388 26.8 162.1 36.1 2.5 properties of anilines, mainly we shall be concerned
meta 1.395 24.7 157.1 38.1 2.5
para 1.391 23.9 158.4 36.4 2.8
about the excitation energies and the ionization poten-
COOH ortho1 1.366 11.2 165.8 21.1 2.2 tial of these systems and their relation with the charge
ortho2 1.376 16.3 163.0 26.8 3.2 distribution in the molecule. In addition, we will
meta 1.396 24.1 155.9 2.8 concentrate our attention in those substituents that,
para 1.387 22.8 159.6 34.4 2.5 in the ortho position, could make an intramolecular
NO2 ortho 1.356 1.81 177.2 2.0 3.1
meta 1.391 23.4 157.7 36.5 2.8
hydrogen bond with the amine group.
para 1.380 20.8 161.6 31.0 2.4
3.2.1. Structural properties of substituted anilines
Table 3 shows some selected structural parameters
of aniline substituted in the ortho, meta or para posi-
tions. The investigated substituents cover a wide
range of s – p electron donor-withdrawing capability.
The torsion angles defined by H1–N1–C1–C2 and
H2–N1–C1–C2 (g 1 andg 2), the u angle, and the d
angle (see Fig. 1) represent the extent of the interac-
tion between the pz orbitals of nitrogen and the p
system of the benzenoid ring. It is particularly inter-
esting that additional planarity of the amine group is
observed when an electron-withdrawing substituent
group is added to the ring. In particular, the nitro
group shows the strongest effects on the aniline
geometry, with all the atoms of o-nitroaniline lying
(almost) in the same plane as the benzenoid ring (see
Fig. 2). When a carboxylic or cyano group is added to
the ring, a similar effect is observed but to a lesser
Fig. 2. Optimized structure of o-nitroaniline. extent. In fact, in these cases coplanarity is not
 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221 215

Fig. 3. Optimized structures of the two stable conformers of o-

aminobenzoic acid.

reached. This additional planarity was reported before

in some experimental and theoretical studies
[32,33,49–54]. This flattening of the amine group is
coupled with a shortening of the bond length between
the nitrogen and the benzene ring (yielding a quinoid
like structure). Fig. 4. Relative total energy of substituted aniline against the
Anilines substituted with nitro and carboxy groups substituent position. All cases were scaled by a constant factor in
in ortho position deserves special attention. As these order to plot them in one graph.
groups contain an oxygen atom that lies next to the
hydrogen of the amine group, an intramolecular
hydrogen bond is probably formed [49,55]. In fact,

Table 4
Mulliken charges of aniline and ortho substituted anilines (in me). Hydrogen charge is summed over the corresponding carbon atoms of the
benzene ring. Calculations were performed using the DFT-B3LYP method with the basis set 6-31G*

Aniline NH2 OCH3 Cl CF3 CO2H (I) CO2H (II) NO2

N1 2786.31 2783.16 2789.37 2791.55 2803.60 2822.20 2815.02 2816.81

H1 317.93 302.49 a 333.70 a 339.24 a 338.25 a 374.84 a 362.74 a 381.99 a
H2 317.93 315.80 315.21 324.57 325.49 325.64 326.03 336.19
C1 306.87 262.23 259.06 322.15 291.32 314.31 294.43 353.30
C2 258.42 262.23 347.60 2118.22 244.47 34.48 43.49 252.46
C3 210.41 274.05 275.92 10.92 229.87 217.41 27.37 22.06
C4 218.77 223.38 228.36 23.90 210.06 217.26 215.35 210.46
C5 210.41 223.38 216.75 21.70 9.83 8.60 9.53 27.07
C6 258.42 27405 269.23 242.76 252.24 258.75 258.02 259.89
N2 2783.16 365.38
C7 2210.94 794.60 555.53 546.61
Cl 238.75
O1 2534.57 2523.43 c 2482.52 2447.91 c
O2 2584.93 2617.21 c 2403.38
X1 b 302.49 c 165.12 2287.18 c 410.58 412.62
X2 b 315.80 153.73 2267.25 d
X3 b 150.72 2264.82
a
Adjacent to the substituent group.
b
X ˆ H or F.
c
Adjacent to amine hydrogen.
d
F adjacent to benzene hydrogen.
 216

Table 5
Total energy and some selected electronic properties of aniline and substituted anilines, obtained employing the DFT-B3LYP with the basis set 6-31G*

Substituent Position Total energy (atomic units) Charge NH2(umeu) HOMO (eV) LUMO (ev) DE(HOMO2LUMO) (eV) m (Debyes)

H 2287.60176 2150.41 25.391 0.251 5.642 1.713

NH2 ortho 2342.94974 2182.19 24.976 0.491 4.484 2.453
meta 2342.95488 2157.15 25.022 0.587 5.609 2.604
para 2342.95112 2162.71 24.645 0.291 4.936 0.0
OH ortho 2362.81737 2138.03 25.111 0.418 5.29 1.479
meta 2362.81813 2152.02 25.319 0.394 5.713 2.741
para 2362.81469 2157.98 24.990 0.123 4.8676 2.1188
CH3 ortho 2326.91964 2152.67 25.327 0.364 5.691 1.736
meta 2326.91974 2153.08 25.319 0.326 5.644 1.622
para 2326.91900 2153.53 25.233 0.266 5.499 1.551
OCH3 ortho 2402.12460 2140.44 25.043 0.509 5.552 1.705
meta 2402.12474 2152.48 25.325 0.416 5.741 1.168
para 2402.12142 2158.45 24.944 0.187 5.131 1.814
Cl ortho 2747.19900 2127.9 25.658 20.124 5.534 2.053
meta 2747.19869 2139.0 25.733 20.151 5.582 3.155
para 2747.19798 2142.6 25.573 20.188 5.384 3.593
CN ortho 2379.84981 2132.25 25.984 21.110 4.874 4.2047
meta 2379.84638 2136.00 26.061 21.184 4.877 5.621
para 2379.84900 2128.77 26.005 20.826 5.180 6.752
CF3 ortho 2624.64116 2139.86 25.811 20.467 5.344 2.541
meta 2624.64030 2139.38 25.853 20.533 5.320 3.705
para 2624.64173 2136.60 25.895 20.250 5.645 4.484
COOH ortho1 2476.18251 2121.65 25.555 21.063 4.492 1.137
ortho2 2476.17770 2126.20 25.638 21.036 4.601 3.130
meta 2476.17578 2142.00 25.722 21.118 4.603 2.384
para 2476.17899 2133.3 25.814 20.804 5.010 3.9176
M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221

NO2 ortho 2492.10951 2101.02 26.041 22.164 3.877 4.759

meta 2492.10497 2132.03 26.135 22.244 3.891 5.656
para 2492.10870 2122.67 26.253 –1.955 4.299 7.118
 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221 217

between a fluorine atom (1) and a hydrogen from

the amine group (2.161 Å) and between fluorine (2)
and the hydrogen bonded in meta position (2.332 Å).
The meta and para trifluoromethylaniline, also, show
a short distance (around 2.390 Å) between one of the
fluorine and the adjacent benzene hydrogen. All these
interactions (resembling hydrogen bonds) are electro-
static, therefore, we should expect large modifications
in the charge distribution of these systems, i.e.,
unequal charge on the hydrogens of the amine group.

3.2.2. Electronic properties of substituted anilines

The charge distribution of aniline and some ortho
substituted anilines are shown in Table 4. The para
and ortho carbon atoms of aniline have the largest
negative charges compared to the other carbon
atoms [11,13,51]. This is in agreement with the
ortho and para directing power of the amine group
for electrophilic substitution known from the experi-
Fig. 5. Charge of the amine group (me) versus pKb of aniline and ment. The ortho carbon atoms show the greatest nega-
substituted aniline (for additional details see Table 5). tive charge, attributed to a large mesomeric effect
(1M) owing to the non-bonding electrons of the
this hydrogen bond is observed in both cases, it amine group. This charge distribution is strongly
induces an additional planarity in the whole system, affected by substituents. There exists a compromise
and in principle, should induce an increase in the between the strength of the substituent and its position
electronic delocalization of the molecule. This inter- in the ring, i.e., considering that the meta position is
action originates a kind of six-member ring fused to deactivated. Considering the substituent position, the
the benzenoid ring, which should confer special elec- unequal charge distribution on the carbons of the ring
tronic properties to the whole system. Fig. 2 shows the will induce differences in the relative stabilization of
optimized structure of o-nitroaniline, which has a each particular system. Hence, anilines substituted
remarkable short distance between the hydrogen with an electron-withdrawing group should be more
(from the amine group) and the oxygen (from the stable when the substituent is placed in the ortho or
nitro group). Also, particular consideration deserves para positions, whereas electron donor substituted
the o-aminobenzoic acid. In this case, there are two anilines show a minimum total energy when the
possible conformations of the carboxylic substituent group is bonded at the para position. This effect is
as is shown in Fig. 3. Structure I depicts the strongest illustrated in Fig. 4 where the relative differences
interaction with the hydrogen of the amine group, this between the total energy of all anilines studied are
being the most stable conformation (see Table 4 and plotted against the substituent position in the ring. In
Fig. 4). Also, hydroxy, methoxy, trifluoromethyl and all cases, the tendency is clear, although some addi-
chlorine substituent groups interact with the hydrogen tional stabilization is observed for ortho substituted
of the amine group, although in these cases the inter- anilines, when the substituent are OCH3, OH, Cl or
action is weaker [32,33,49,56–58]. In o-methoxyani- COOH (I). Particularly, a stabilization of 3 kcal/mol is
line, the distance between the oxygen and the observed when comparing the two possible conforma-
hydrogen is 2.251 Å, and in o-aminophenol this tions (I and II) of the o-aminobenzoic acid. This
length is 2.288 Å, while in o-chloroaniline the chlor- energy difference indicates the extent of the interac-
ine–hydrogen distance is 2.593 Å, which is very short tion between the amine group and the oxygen of the
considering the radius of chlorine. Also, o-trifluoro- substituent adjacent to it.
methylaniline shows particularly short distances The lack of electrons in the benzene ring produced
218 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221

by electron-withdrawing substituents is the main

reason for the extra planarity of the amine group,
and it produces the shortening of the nitrogen–carbon
bond distance, which is characteristic of the charged
anilines. However, not only the strength of the substi-
tuent group is important, but also its position plays an
important role in the induced conjugation.
By analyzing the charge distribution of ortho
substituted anilines, it is possible to explain some of
the special behavior reported experimentally (ortho
effect) [49]. In fact, in all substituted anilines reported
in Table 5 the hydrogens of the amine group support
different amounts of charge. Having the hydrogen
next to the substituent the lowest electronic density.
Moreover, the oxygen atoms of o-nitroaniline also
exhibit an asymmetric charge distribution. These
facts, together with the short distance between oxygen
and hydrogen, support the assumption of intramole-
cular hydrogen bond formation, which induces addi-
tional stabilization to the molecule [59,60]. A similar
effect is reported in o-aminobenzoic acid. In this case,
the electronic density on the oxygen atoms depends on
the conformation of the carboxy group (I or II, see Fig.
3). There exists a more positive electrostatic interac-
tion between the hydrogen of the amine and the
oxygen from the carbonyl group (more than with the
hydroxyl group). In contrast, the fluorine atom of o-
trifluoromethylaniline next to the amine group shows
a more negative charge than the rest of the fluorine
atoms. Although, in this case the electrostatic effect is
less significant than in carboxy and nitroanilines, an
interaction between the substituent and the amine
group is observed. Similar behavior is found for o-
methoxyaniline and o-chloroaniline. Finally, the
opposite effect is observed in o-diaminobenzene
where the amine groups repel each other. As the
charge of the amine group is related to the basicity
of the aromatic amines, we plotted the charge of the
amine group against the experimental pKb of several
amines [49,61] in order to compare our results (see
Fig. 5). In fact it is possible to establish a quite linear
Fig. 6. (a) Plot of the energy of the HOMO of substituted anilines
against the substituent position in the benzene ring. The HOMO of
relationship between our calculated charges and the
aniline is denoted by a dotted line. (b) Plot of the energy of the experimentally determined pKb, and therefore it is
LUMO of substituted anilines against the substituent position in the possible to consider the intramolecular electrostatic
benzene ring. The LUMO of aniline is denoted by a dotted line. interaction as a very plausible explanation for the
ortho effect. It is known that DFT, particularly at
the LDA level, overestimate the hydrogen-bond
energy [62,63]. However, our calculations show that
 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221 219

related to the vertical IP by applying Koopmans theo-

rem. Although, as was explained earlier, IP values
calculated by using the adiabatic approximation
depict large error, we expect that the relative trends
of the experimental IP will be adequately described by
this approximation. A variation of the IP in agreement
with the electron donor-withdrawing capability of the
substituent group is observed [49].
Electron-donor groups, especially the amino group,
interact mainly with the HOMO orbital of aniline,
whereas their interaction with the LUMO is minor
(see Fig. 6 a,b). The HOMO perturbation is particu-
larly strong when the substituent is placed at position
ortho or para in the benzene ring, as is well known.
When considering electron-withdrawing groups, there
is a stabilization of the HOMO orbital but the stron-
gest perturbation is produced between LUMO of the
substituent and the LUMO of aniline. In most cases
the final IPs and EAs obtained agree with the electron
Fig. 7. Plot of the difference between the HOMO and the LUMO donor-withdrawing capability of the substituent
energies against the substituent position in the benzene ring. The group.
DE(HOMO2LUMO) of aniline is denoted by a dotted line. Large variations of the excitation energies are
induced when a substituent group is bonded to the
the results obtained are comparable to those obtained benzene ring, owing to a large perturbation between
by AM1 [31] (which is quite accurate at predicting orbitals. Fig. 7 shows the first excitation energy of a
hydrogen bonds). Further, the explanation of the family of anilines calculated as DE(HOMO-LUMO).
experimental ortho effect is one of the most powerful Even though it was demonstrated that there is a large
arguments in favour of hydrogen bonds predicted by error in the absolute value of the excitation energies
B3LYP. Moreover, as B3LYP is a hybrid method it is calculated in this way, it is possible to consider this
not so unexpected that an improvement in the predic- procedure as an appropriate first approximation. In
tion of the hydrogen bonds is found compared to non- general, most of the substituents diminish the first
hybrids methods. electronic transition. Particularly, anilines substituted
Although, the aniline orbitals were largely with strong electron-withdrawing groups (such as
described before [6,64–68], we will briefly mention nitro) depict the smallest excitation energies. More-
here some relevant aspects in terms of our discussion. over, anilines substituted in the position ortho or para
First, the HOMO shows a large contribution from the with electron-donor groups, show that the HOMO
amine group and the benzene ring, whereas the perturbation is so large that this produces a substantial
LUMO corresponds mainly to the benzene ring lowering of the first excitation energy.
having a very small contribution from the meta Finally, we shall mention that anilines substi-
carbons. Therefore, it is expected that the IP and EA tuted with electron withdrawing groups show a
of aniline will be strongly affected by substitution. large permanent dipolar moment (m ). Anilines
There is a very close relationship between the strength substituted in the para position have attracted
of the substituent, the charge of the amine group, the considerable interest because of their possible
IP, and the basicity of the system [32]. Consequently, application in electro-optic modulation and second
the IP, as well as the charge distribution of substituted and third harmonic generation (SHG, and THG).
anilines, should contain a lot of information relevant The dipolar moments of the studied systems are
to practical purposes. Table 5 shows the HOMO of included in Table 5, which shows p-nitroaniline
aniline and its derivatives, which can be directly and p-cyanoaniline as the most suitable
220 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221

compounds of this family to be utilised in such Acknowledgements

applications. In fact these molecules have already
been experimentally studied, showing them to This work was supported by EPSRC through the
have large m and b values grant GR/L76037. Also, we are very grateful for the
reviewer’s useful comments and suggestions.

4. Conclusions References

In this work, we have shown calculations of aniline [1] F. Kippert, D. Lloyd, FEMS Microbiol. Lett. 132 (1994) 215.
[2] L. Olmedo, P. Hourquebie, F. Jousse, in: H.S. Nalwa (Ed.),
performed using density functional methods. We have Handbook of Organic Conductive Molecules and Polymers, 3,
compared several kinds of exchange-correlation func- Chichester, UK, 1997, pp. 367.
tionals. The hybrid schemes (such as B3LYP, B3P86, [3] A.G. MacDiarmid, A.J. Epstein, in: S.A. Jenekhe, K.J. Wynne
and B3PW91) have proved to be very accurate in (Eds.), Photonic and Optoelectronic Polymers, American
predicting the structural properties. Electronic proper- Chemical Society, Washington, DC, 1997, pp. 395.
[4] A.J. Epstein, M.R.S. Bull. 22 (1997) 16.
ties calculated using the gradient-corrected exchange [5] N. E. Agbor, J. P. Cresswell, M. C. Petty, A. P. Monkman,
seem to be more accurate, however it could be a Sensors and Actuators B: Chemical 41 (1997) 137.
product of its poor structural prediction (change in [6] J.O. Morley, J. Phys. Chem. 99 (1995) 1923.
the structure ; change in the electronic properties). [7] Y. Wang, S. Saebø, C.U. Pittman Jr., J. Mol. Struct. (Theo-
The ionization potentials as well as the excitation chem) 281 (1993) 91.
[8] D.G. Lister, J.K. Tyler, Chem. Commun. 00 (1966) 152.
energies are estimated by applying the adiabatic [9] D.G. Lister, J.K. Tyler, J.H. Høg, N. Wessel Larsen, J. Mol.
approximation, which provides a very good first Struct. 23 (1974) 253.
approach to the experimental values, i.e. calculated [10] M. Fukuyo, K. Hirotsu, T. Higuchi, Acta Cryst. B38 (1982)
excitation energy 5.1–5.6 eV, experimental excitation 640.
energies 4.4 eV. [11] W.J. Hehre, L. Radom, J.A. Pople, J. Am. Chem. Soc. 94
(1972) 1496.
Substituents attached to the benzene ring of [12] W.J.E. Parr, R.E. Wasylishen, J. Mol. Struct. 38 (1977) 272.
aniline produce large structural as well as electro- [13] A.D. Gorse, M. Pesquer, J. Mol. Struct. (Theochem) 281
nic distortions. Particularly nitro groups, which (1993) 21.
exert the most radical effects over the charge [14] A. Wolf, U. Voets, H.-H. Schmidtke, Theoret. Chim. Acta 54
distribution and consequently over the IP, EA, (1980) 229.
[15] E. Von Nagy-Felsobuky, R. D. Topsom, S. Pollack, R.W.
and excitation energies. In general, substituents Taft, J. Mol. Struct. (Theochem) 88 (1982) 255.
in the ortho position show additional interactions [16] Z. Niu, J.E. Boggs, J. Mol. Struct. (Theochem) 109 (1984)
with the amine group. Nitro and carboxy groups, 381.
in particular produce strong intramolecular hydro- [17] T. Ziegler, Chem. Rev. 91 (1991) 651.
gen bonds with the amine group. Despite such [18] U. Salzner, P.G. Pickup, R.A. Poirier, J.B. Lagowski, J. Phys.
Chem. A102 (1998) 2572.
particular structural features, there is a linear rela- [19] L.A. Curtiss, P.C. Redfern, K. Raghavachari, J.A. Pople, J.
tionship between the charge on the amine group Chem. Phys. 109 (1998) 42.
and the experimental pKb, which allows us to [20] P. Hohenberg, W. Kohn, Phys. Rev. B136 (1964) 864.
conclude that these intramolecular hydrogen [21] W. Kohn, L. Sham, Phys. Rev. A 140 (1965) 1133.
bonds could be the origin of the experimentally [22] M. Springborg, in: M. Springborg (Ed.), Density-Functional
Methods in Chemistry and Materials Science, John Wiley and
well-known ‘‘ortho effect’’. Sons, Chichester, 1997, pp. 207.
In general, the electronic properties predicted [23] M. Peric, B. Engels, S.D. Peyerimhoff, in: S.R. Langhoff
here have a small error, which is a consequence (Ed.), Quantum Mechanical Structure Calculations with
of the well-known ‘overall shifting’ of the orbital Chemical Accuracy, Kluwer Academic Publisher, Dordrecht,
diagram observed in DFT calculations. We hope to 1995.
[24] B. O. Roos, M. Fülscher, P. Malmqvist, M. Merchán, L.
improve these results in the near future by using Serrano-Andrés, in: S.R. Langhoff (Ed.), Quantum Mechan-
more sophisticated methods such as MR–CI or ical Structure Calculations with Chemical Accuracy, Kluwer
TDDFT. Academic Publisher, Dordrecht, 1995.
 M.E. Vaschetto et al. / Journal of Molecular Structure (Theochem) 468 (1999) 209–221 221

[25] J.F. Stanton, R.J. Bartlett, J. Chem. Phys. 98 (1993) 7029. [46] G.L. Butt, R.D. Topsom, Spectrochim. Acta Part A 34 (1978)
[26] R. Bauernschmidtt, R. Ahlrichs, Chem. Phys. Lett. 256 (1996) 975.
454. [47] J.P. LaFemina, C.B. Duke, A. Paton, J. Chem. Phys. 89 (1988)
[27] R. Bauernschmidtt, R. Ahlrichs, F.H. Hennrich, M.M. Kappes, 2668.
J. Am. Chem. Soc. 120 (1998) 5052. [48] V.E. Borisenko, A.V. Baturin, M. Przeslawska, A. Koll, J.
[28] M. Casida, in: D.P. Chong (Ed.), Recent Advances in Density Mol. Struct. 407 (1997) 53.
Functional Methods, 1, World Scientific, Singapore, 1995. [49] R. Stewart, in: H.H. Wasserm (Ed.), The Proton: Applications
[29] R. Disselkamp, H.S. Im, E.R. Berstein, J. Chem. Phys. 97 to Organic Chemistry, 46, Academic Press, London, 1985.
(1992) 7889. [50] R. Nakagaki, I. Aoyama, K. Shimizu, M. Akagi, J. Phys. Org.
[30] V.E. Borisenko, E.I. Greseva, E.V. Dukhnova, I.V. Nacheva, Chem. 6 (1993) 261.
J. Mol. Struct. 324 (1994) 199. [51] L. Serrano-Andres, M. Merchan, B.J. Roos, R. Lindh, J. Am.
[31] M.E. Vaschetto, B.A. Retamal, J. Phys. Chem. A 101 (1997) Chem. Soc. 117 (1995) 3189.
6945. [52] M. Pierre, J.B. Ehooman, G.R. Neergard, J.L. Coutaz, I.
[32] M.E. Vaschetto, A.P. Monkman, M. Springborg. J. Mol. Ledoux, Optics. Commun. 98 (1993) 209.
Struct. (Theochem) 468 (1999) 181. [53] P.-C. Maria, I. Leito, J.-F. Gal, O. Exner, M. Decouzon, Bull.
[33] M.J.S. Dewar, E.G. Zoebisch, E.F. Healy, J.J.P. Stewart, J. Soc. Chim. Fr. 132 (1995) 394.
Am. Chem. Soc. 107 (1985) 3902. [54] R. Nakagaki, D.G. Frost, C.A. McDowell, J. Electron. Spec-
[34] A.D. Becke, Phys. Rev. A 38 (1988) 3098. trosc. Relat. Phenom. 22 (1981) 289.
[35] A.D. Becke, The original three parameters hybrid, J. Chem. [55] G.S. Denisov, L.A. Kuzina, Mol. J. Struct. 271 (1992) 9.
Phys. 98 (1993) 5648. The form implemented in Gaussian [56] M. Khalil, R.J. Woods, D.F. Weaver, V.H. Smith Jr., J.
packages is slightly modified. Comput. Chem. 12 (1991) 584.
[36] C. Lee, W. Yang, R.G. Parr, Phys. Rev. B 37 (1988) 785. [57] E. Borisenko, A.V. Morev, A. Koll, J. Mol. Struct. 444 (1998)
[37] J. P. Perdew, Phys. Rev. B 33 (1986) 8822. 183.
[38] J.P. Perdew, in: P. Ziesche, H. Eschrig (Eds.), Electronic Prop- [58] A. Simperler, H. Lampert, W. Mikenda, J. Mol. Struct. 448
erties of Solids ’91, Akademie Verlag, Berlin, 1991. (1998) 191.
[39] J.P. Perdew, Y. Wang, Phys. Rev. B 45 (1992) 13 244. [59] P.I. Nagy, G.J. Durant, D.A. Smith, in: D.A. Smith (Ed.),
[40] Gaussian 94, Revision E.2, M.J. Frisch, G.W. Trucks, H.B. Modeling the Hydrogen Bond, American Chemical Society,
Schlegel, P.M.W. Gill, B.G. Johnson, M.A. Robb, J.R. Washington, DC, 1994, pp. 60.
Cheeseman, T. Keith, G.A. Petersson, J.A. Montgomery, K. [60] R.J. Woods, W.A. Szarek, V. H. Smith Jr., J. Chem Soc.
Raghavachari, M.A. Al-Laham, V.G. Zakrzewski, J.V. Ortiz, Chem. Commun. 00 (1991) 334.
J.B. Foresman, J. Cioslowski, B.B. Stefanov, A. Nanayakkara, [61] Dictionary of Organic Compounds, 5th edition, Chapman and
M. Challacombe, C.Y. Peng, P.Y. Ayala, W. Chen, M.W. Hall, London, 1985.
Wong, J.L. Andres, E.S. Replogle, R. Gomperts, R.L. Martin, [62] F. Sim, A. St-Amant, I. Papai, D.R. Salahub, J. Am. Chem.
D.J. Fox, J.S. Binkley, D.J. Defrees, J. Baker, Washington, Soc. 114 (1992) 4391.
DC, American Chemical Society, J.P. Stewart, M. Head- [63] J.M. Pérez-Jordá, A.D. Becke, Chem. Phys. Lett. 233 (1995)
Gordon, C. Gonzalez, and J.A. Pople, Gaussian, Inc., 134.
Pittsburgh, PA, 1995. [64] B. Sjögren, S. Stafström, J. Chem. Phys. 88 (1988) 3840.
[41] Cerius 2 User Guide, Molecular Simulations Inc., San Diego. [65] M. Quack, M. Stockburger, J. Mol. Spectrosc. 43 (1972) 87.
[42] B.G. Johnson, P.M.W. Gill, J.A. Pople, J. Chem. Phys. 98 [66] A. R. Bacon, J. M. Hollas, Faraday Discuss. Chem. Soc. 86
(1993) 5612. (1988) 129.
[43] B. Mannfors, J.T. Koskinen, L.-O. Pietilä, L. Ahjopalo, J. [67] X. Song, M. Yang, E.R. Davidson, J.P. Reilly, J. Chem. Phys.
Mol. Struct. (Theochem) 393 (1997) 39. 99 (1993) 3224.
[44] T. Koopmans, Physica (The Hague) 10 (1933) 104. [68] K. Kimura, S. Nagakura, Mol. Phys. 9 (1965) 117.
[45] U. Salzner, J. B. Lagowski, P.G. Pickup, R.A. Poirier, J.
Comput. Chem. 18 (1997) 1943.