Understanding and Interpreting Molecular Electron Density Distributions
Total Page:16
File Type:pdf, Size:1020Kb
Research: Science and Education Understanding and Interpreting Molecular Electron Density Distributions C. F. Matta and R. J. Gillespie* Department of Chemistry, McMaster University, Hamilton, ON L8S 4M1, Canada; *[email protected] Advances in the power and speed of computers have made concerned in this paper. A more formal discussion of electron ab initio and density functional theory (DFT) calculations density is presented in Appendix 2. an almost routine procedure on a PC (1). From these calcu- The electron density is key to the bonding and geometry lations the equilibrium geometry, the energy, and the wave of a molecule because the forces holding the nuclei together function of a molecule can be determined. From the wave in a molecule are the attractive forces between the electrons and function one can obtain all the properties of the molecule, the nuclei. These attractive forces are opposed by the repulsions including the distribution of electronic charge, or electron between the electrons and the repulsions between the nuclei. density. However, the electron density is often not calculated In the equilibrium geometry of a molecule these electrostatic or discussed, perhaps because it is not widely realized that forces just balance. The fundamentally important Hellman– very useful information on bonding and geometry can be Feynman theorem (4–7) states that the force on a nucleus in a obtained from it. It seems particularly important to discuss molecule is the sum of the Coulombic forces exerted by the other electron densities in introductory chemistry courses because nuclei and by the electron density distribution ρ. This means students can grasp the concept of electron density much more that the energy of interaction of the electrons with the nuclei readily than the abstract mathematical concept of an orbital. It can be found by a consideration of the classical electrostatic is also not as widely understood as it should be that orbitals forces between the nuclei and the electronic charge cloud. There are not physical observables but only mathematical constructs are no mysterious quantum mechanical forces, and no other that cannot be determined by experiment (2). In contrast, the force, such as the gravitational force, is of any importance in electron density distribution in a molecule or crystal can be holding the atoms in a molecule together. The atoms are held observed by electron diffraction and X-ray crystallography together by the electrostatic force exerted by the electronic charge (3); and it can also, and often more readily, be obtained from on the nuclei. But it is quantum mechanics, and particularly ab initio and density functional theory calculations. the Pauli principle, that determines the distribution of elec- This article gives a simple introduction to the electron tronic charge, as we shall see. densities of molecules and how they can be analyzed to ob- tain information on bonding and geometry. More detailed The Representation of the Electron Density discussions can be found in the books by Bader (4), Popelier (5), and Gillespie and Popelier (6 ). Computational details The electron density (ρ) varies in three dimensions (i.e., to reproduce the results presented in this paper are presented it is a function of the three spatial coordinates [x,y,z]), so a in Appendix 1. full description of how ρ varies with position requires a fourth dimension. A common solution to this problem is to show how ρ The Electron Density varies in one or more particular planes of the molecule. Figure 1a shows a relief map of the electron density, ρ, of the σ Quantum mechanics allows the determination of the SCl2 molecule in the v(xz) plane. The most striking features probability of finding an electron in an infinitesimal volume of this figure are that ρ is very large in an almost spherical surrounding any particular point in space (x,y,z); that is, the region around each nucleus while assuming relatively very probability density at this point. Since we can assign a prob- small values, and at first sight featureless topology, between ability density to any point in space, the probability density these nuclear regions. The high electron density in the nearly defines a scalar field, which is known as the probability density spherical region around each nucleus arises from the tightly distribution. When the probability density distribution is held core electrons; the relatively very small and more diffuse multiplied by the total number of electrons in the molecule, density between these regions arises from the more weakly N, it becomes what is known as the electron density distribution held bonding electrons. In fact, it was necessary to truncate or simply the electron density and is given the symbol ρ(x,y,z). the very high maxima in Figure 1a (at ρ = 2.00 au)1 to make It represents the probability of finding any one of the N it possible to show the features of the electron density distri- electrons in an infinitesimal volume of space surrounding bution between the nuclei. In particular, there is a ridge of the point (x,y,z), and therefore it yields the total number of increased electron density between the sulfur atom and each electrons when integrated over all space. The electron density of the chlorine atoms. The electron density has values of can be conveniently thought of as a cloud or gas of negative 3.123 × 103 and 2.589 × 103 au at the S and Cl nuclei, re- × ᎑1 charge that varies in density throughout the molecule. Such a spectively, but a value of only 1.662 10 au at the minimum charge cloud, or an approximate representation of it, is often of the ridge between the peak around the sulfur nucleus and used in introductory texts to represent the electron density each of the chlorine nuclei. This ridge of increased electron ψ2 of an atomic orbital. It is also often used incorrectly to density between the S atom and each of the Cl atoms, small as depict the orbital ψ itself. In a multielectron atom or molecule it is, is the density in the bonding region that is responsible for only the total electron density can be experimentally observed pulling the nuclei together. Along a line at the top of this ridge or calculated, and it is this total density with which we are the electron density is locally greater than in any direction JChemEd.chem.wisc.edu • Vol. 79 No. 9 September 2002 • Journal of Chemical Education 1141 Research: Science and Education a Figure 1. Relief maps of the electron density of (a) SCl2 and (b) H2O in the plane of the nuclei (density and distances from the origin of the coordinate system in au). Isodensity contour lines are shown in the order 0.001, 0.002, 0.004, 0.008 (four outermost contours); 0.02, 0.04, 0.08 (next three); 0.2, 0.4, 0.8 (next three). The density is truncated at 2.00 au (innermost contour). These contours are b shown in blue, violet, magenta, and green, respectively, on the Figure 2. Contour maps of the electron density of (a) SCl2 and (b) figure in the table of contents (p 1028). H2O. The density increases from the outermost 0.001 au isodensity contour in steps of 2 × 10n, 4 × 10n, and 8 × 10n au with n start- ing at ᎑3 and increasing in steps of unity. The lines connecting the away from the line. This line coincides with the bond between nuclei are the bond paths, and the lines delimiting each atom are the S and Cl atoms as it is normally drawn and is called a the intersection of the respective interatomic surface with the plane bond path (4–6, 10). The point of minimum electron density of the drawing. The same values for the contours apply to subse- along the bond path is called the bond critical point. quent contour plots in this paper. Figure 1b shows a relief map of the electron density of the water molecule. Its features are similar to those of the density map for SCl2, but the electron density around the hydrogen of the electron population. Thus this outer contour shows atoms is much smaller than around the oxygen atom, as we the shape of the molecule in the chosen plane. In a condensed would expect. The electron density at the maximum at the phase the effective size of a molecule is a little smaller. We oxygen nucleus has a value of 2.947 × 102 au, whereas that see more clearly here that the bond paths (the lines along at the hydrogen nucleus is only 4.341 × 10᎑1 au, which is only the top of the density ridges between the nuclei) coincide slightly greater than the value of 3.963 × 10᎑1 au at the with the bonds as they are normally drawn. minimum at the bond critical point. The very small electron Figure 3 shows the electron density contour maps for density surrounding the hydrogen nucleus is due to less than the period 2 fluorides LiF, BF3, CF4, OF2, and for the iso- one electron because the more electronegative oxygen atom lated B atom. In LiF each atom is almost spherical, consis- attracts electron density away from the hydrogen atom so that tent with the usual model of this molecule as consisting of it has a positive charge. the ions Li+ and F᎑. The volume of the lithium atom is much Another common way to represent the electron density smaller than that of the F atom, again consistent with the ionic distribution is as a contour map, analogous to a topographic model. We will see later that we can also obtain the atomic contour map representing the relief of a part of the earth’s charges from the electron density and that the charges on the surface.