The HSAB concept can be developed into an universal scheme based on quantum chemistry

H. Bögel*, U. Laube, S. Hansen

Martin-Luther-University of Halle, Institute of Physical Chemistry Geusaerstr., D-06217 MERSEBURG (Germany)


1. Introduction

The electronic structure of chemical compounds is the fundamental quantity of that, what we call "STRUCTURE". In this sense the electronic structure is responsible for the strength of chemical bonds, for the spatial structure, the chemical behaviour, and all the physico-chemical properties. In the language of quantum mechanics, the properties (observables) are expectation values of the corresponding operators applied to the molecular wavefunctions. This should be kept in mind for all Structure-Property- Relations, we want to derive in our investigations or exploit in predicting new properties or activities for new structures.

But very often these structure activity relations (SAR) are based on simple quantities and indices, which were not derived from quantum mechanics. Much of our actual knowledge in chemistry is composed of empirical or semi-empirical rules and concepts. Some of the chemical terms in use have or had no strong definition, but they are very popular and were used to explain a lot of facts in chemistry in qualitative accuracy. This was the case for the well-known HARD and SOFT / ACIDS and BASES concept and related terms, like:

On the other hand such an empirical way of thinking has established in the chemistry - as thoroughly useful, which works with ideas, which do not originate from a ´strong´ theory, but clearly and vividly describes a series of relations among chemical data. Because of their simplicity they accompany our structural - thinking and the technical terms in our language in chemistry.

Tab.1 Some characteristics of the general term ´STRUCTURE´
spatialelectronicenergetic
position of the nucleielectron distributioninteractions
atomic coordinateselectron densitiesforce constants
bond length, bond angleatomic chargesmolecular orbitals
molecular volumn
CPK- or Ball-Stick-model
molecular electrostatic potentialtotal energy
heat of formation

In the interpretation of intermolecular interactions and chemical reactivity we have to take into account that the reactivity of the molecule as a whole may be divided into regions of different behaviour (regio-selectivity). So we have to replace some global characteristics by local ones. Visualisation of the extent of molecular properties in different areas of the molecule is a helpful and important technique to support structural concepts by computer graphics, to relate the expected or observed physico-chemical properties and the chemical behaviour of molecules to the calculated electronic structure.

2. The concept of Hard and Soft Acids and Bases

Hard- and softness in the sense of the HSAB (Hard and Soft, Acids and Bases) concepts [1] were developed on the basis of empirical data. Within the 30 years of development of the HSAB concept there is some remarkable progress, which can should be summarised with the following steps:

This contribution shows, how local atomic quantities of hardness and softness, generated from semi-empirical and ab initio MO-calculations followed by CSA, can be visualized.

In table 2 we summarise some conventional terms of the electronic structure compared with those used in the HSAB concept. Similar to Structure & Dynamics (Tab.2) we have to take into account both of these aspects of an electronic structure (Tab. 3).

Tab. 2 structure and dynamics in relation to the HSAB concept
STRUCTUREDYNAMICS
geometrical structurespatial arrangement of atomsfexibility
exp. observationX-ray scatteringNMR relaxation
electronic structureelectron densityelectron moveability
exp. observationpolarity, dipole moment
my (vector)
polarisability
alpha (tensor)
HSAB concepthardnesssoftness

Tab. 3: Convenient terms to describe the electronic structure
localglobalHSAB
staticatomic chargepolarityhardness
dynamicmoveabilitypolarizabilitysoftness

3. Charge Sensitive Analysis

Within the quantum mechanical Charge Sensitive Analysis (CSA) [6]) we are able to derive local atomic quantities for electronegativity, hardness and softness besides the global quantities from calculated geometrical and electronic structure data. These data can be converted from any quantum mechanical calculation. We use it for semiempirical MO programs (MOPAC, ZINDO) and some ab initio programs (GAUSSIAN94, TURBOMOLE, GAMESS, ADF, DGauss, deMon) successfully.

At the moment we can calculate the following CSA-values [8]:

For the visualisation of these local atomic quantities on a molecular surface we developed a program, written in C running under AIX operating system in a X-Windows environment on a IBM-RS/6000-Graphics-Workstation. The generated molecular surface (Connolly Surface) is coloured with the actual atomic quantities in two kinds:

For the visualisation one can chose from three different molecular representations:

The next two figures show some examples of application of the CSA_VIEW program. The calculated CSA atomic hardness of the aniline molecule are calculated from atomic MNDO-charges (MOPAC program). These pictures show the two possible colouring models applied to the same molecule.

Fig. 1: Aniline atomic hardness coloured with the "rain bow" model Fig. 2: Aniline atomic hardness coloured with the "intensity" model
Fig. 1: Aniline atomic hardness coloured with the "rain bow" model Fig. 2: Aniline atomic hardness coloured with the "intensity" model

It can be seen that the picture 1 gives a better survey about the different elements and their quantities than "intensity" colouring of picture 2. Both pictures are generated from the same set of CSA data.

Figure 1 gives a good impression of the equivalent atomic positions and show the electronic effect in ortho-, meta, and para-position to the amino group in the benzene ring.

The next example (Fig. 3) is the nitrobenzene molecule showing the withdrawing substituent effect of the nitro group.

Fig. 3: Atomic hardness in nitrobenzene
Fig. 3: Atomic hardness in nitrobenzene

The following three examples wants to show the electronic effect of hetero atoms in aromatic rings.This pictures show the CSA values of the electron distribution of the valence electrons, but there is no possibility for a separate representation of electrons belonging to the sigma- or to the pi-symmetry (sometimes do they have quite different effects, -I and +M).

Fig. 4: Atomic hardness of pyrrole Fig. 5: Atomic hardness of furane
Fig. 4: Atomic hardness of pyrrole Fig. 5: Atomic hardness of furane

Fig. 6 Atomic hardness of pyridine molecule
Fig. 6 Atomic hardness of pyridine molecule

Further development of the visualisation of the two important properties of an electronic structure, the hard- and softness will be the interpretation of the calculated Fukui function and a distance dependent contribution from neighbouring atoms to the value of an atomic surface area.

4. Summary

With the 30 years of history, the HSAB concept made an enormous progress. It has developed has been developed from an empirical qualitative to a quantitative concept on the basis of modern quantum mechanics. The global molecular quantities like electronegativity, hardness, and softness can be calculated. Within the Charge Sensitive Analysis, it is possible to go to the essential regional values to characterise the chemical behaviour and chemical reactivity of the alternative centres of a molecule.

We have programmed recently a pictorial representation of hardness and softness on the molecular surface of molecules. The derived numerical values can be visualised with the CSA_VIEW program [9] running with X-Windows on UNIX graphic workstations. The colouring of the atomic spheres is similar to those well-known pictures of the Molecular Electrostatic Potential (MEP) or the Molecular Lipophilic Potential (MLP) which were successfully used in molecular modelling.

We would like at this point cordially thank for the friendly support through the colleagues K. Jug (Hannover), H. Preuss (Stuttgart), and J. Hinze (Bielefeld). Most this work presented here had been supported by a project of the Deutsche Forschungsgemeinschaft (DFG).

5. References

[1] Pearson, R. G.: J. Am. Chem. Soc. 85 (1963) 3533
[2] Klopman, G.: J. Am. Chem. Soc. 90 (1968) 223
[3] Parr, R. G., Pearson, R. G.: J. Am. Chem. Soc. 105 (1983) 7512
[4] Yang, W., Parr, R.G.: Proc. Natl. Acad. Sci. USA 82 (1985) 6723
[5] Pearson, R. G.: J. Am. Chem. Soc. 107 (1985) 6301
[6] Pearson, R. G.: Coord. Chem. Rev. 100 (1990) 403
[7] Nalewajski, R. F.: Struct. Bonding 80 (1993) 115
[8] Laube, U., Boegel, H., Dettmann, J.: Software Developments in Chemistry (GDCh-Frankfurt) Vol. 8 (1994) 251 - 264
[9] CSA_VIEW, X-Windows program written in C, Merseburg (1995)

Arbeitsgruppe Computerchemie am Institut für physikalische Chemie der Martin-Luther-Universität Halle-Wittenberg