Part 2. Three Primary Areas of Theoretical Chemistry Chapter 5. An Overview of Theoretical Chemistry In this Chapter, many of the basic concepts and tools of theoretical chemistry are discussed only at an introductory level and without providing much of the background needed to fully comprehend them. Most of these topics are covered again in considerably more detail in Chapters 6-8, which focus on the three primary sub-disciplines of the field. The purpose of the present Chapter is to give you an overview of the field that you will learn the details of in these later Chapters . What is Theoretical Chemistry about The science of chemistry deals with molecules including the radicals, cations, and anions they produce when fragmented or ionized. Chemists study isolated molecules (e.g., as occur in the atmosphere and in astronomical environments), solutions of molecules or ions dissolved in solvents, as well as solid, liquid and plastic materials mprised of molecules. All such forms of molecular matter are what chemistry is about Chemical science includes how to make molecules(synthesis), how to detect and quantitate them(analysis ), how to probe their properties and how they undergo or change as reactions occur(physical)
1 Part 2. Three Primary Areas of Theoretical Chemistry Chapter 5. An Overview of Theoretical Chemistry In this Chapter, many of the basic concepts and tools of theoretical chemistry are discussed only at an introductory level and without providing much of the background needed to fully comprehend them. Most of these topics are covered again in considerably more detail in Chapters 6-8, which focus on the three primary sub-disciplines of the field. The purpose of the present Chapter is to give you an overview of the field that you will learn the details of in these later Chapters. I. What is Theoretical Chemistry About? The science of chemistry deals with molecules including the radicals, cations, and anions they produce when fragmented or ionized. Chemists study isolated molecules (e.g., as occur in the atmosphere and in astronomical environments), solutions of molecules or ions dissolved in solvents, as well as solid, liquid, and plastic materials comprised of molecules. All such forms of molecular matter are what chemistry is about. Chemical science includes how to make molecules (synthesis), how to detect and quantitate them (analysis), how to probe their properties and how they undergo or change as reactions occur (physical)
A. Molecular Structure-bonding, shapes, electronic structures One of the more fundamental issues chemistry addresses is molecular structure, which means how the molecule's atoms are linked together by bonds and what the interatomic distances and angles are. Another component of structure analysis relates to what the electrons are doing in the molecule that is how the molecule' s orbitals are ocupied and in which electronic state the molecule exists. For example, in the arginine molecule shown in Fig. 5. 1, a HooC-carboxylic acid group(its oxygen atoms are shown in red) is linked to an adjacent carbon atom (yellow) which itself is bonded to an -NH amino group (whose nitrogen atom is blue). Also connected to the a-carbon atom are a chain of three methylene-CH2-groups, a-NH-group, then a carbon atom attached both by a double bond to an imine -NH group and to an amino-NH, group
2 A. Molecular Structure- bonding, shapes, electronic structures One of the more fundamental issues chemistry addresses is molecular structure, which means how the molecule’s atoms are linked together by bonds and what the interatomic distances and angles are. Another component of structure analysis relates to what the electrons are doing in the molecule; that is, how the molecule’s orbitals are ocupied and in which electronic state the molecule exists. For example, in the arginine molecule shown in Fig. 5.1, a HOOC- carboxylic acid group (its oxygen atoms are shown in red) is linked to an adjacent carbon atom (yellow) which itself is bonded to an –NH2 amino group (whose nitrogen atom is blue). Also connected to the a-carbon atom are a chain of three methylene –CH2 - groups, a –NH- group, then a carbon atom attached both by a double bond to an imine –NH group and to an amino –NH2 group. 1.987
Figure 5. 1 The arginine molecule in its non-zwitterion form with dotted hydrogen bond The connectivity among the atoms in arginine is dictated by the well known valence preferences displayed by H, C, O, and N atoms. The internal bond angles are, to a large extent,also determined by the valences of the constituent atoms (i.e, the sp or sp nature of the bonding orbitals ). However, there are other interactions among the several functional groups in arginine that also contribute to its ultimate structure. In particular, the hydrogen bond linking the a-amino group's nitrogen atom to the -- group's hydrogen atom causes this molecule to fold into a less extended structure than it What does theory have to do with issues of molecular structure and why is knowledge of structure so important? It is important because the structure of a molecule has a very important role in determining the kinds of reactions that molecule will undergo, what kind of radiation it will absorb and emit, and to what active sites"in neighboring molecules or nearby materials it will bind. A molecules shape(e.g, rod like, flat, globular, etc. )is one of the first things a chemist thinks of when trying to predict where, anotner molecule or on a surface or a cell the molecule will fit""and be able to bind and perhaps react. The presence of lone pairs of electrons(which act as Lewis base sites), of T orbitals(which can act as electron donor and electron acceptor sites), and of highl polar or ionic groups guide the chemist further in determining where on the molecule's framework various reactant species(e.g, electrophylic or nucleophilic or radical)will be most strongly attracted Clearly, molecular structure is a crucial aspect of the chemists toolbox
3 Figure 5.1 The arginine molecule in its non-zwitterion form with dotted hydrogen bond. The connectivity among the atoms in arginine is dictated by the well known valence preferences displayed by H, C, O, and N atoms. The internal bond angles are, to a large extent, also determined by the valences of the constituent atoms (i.e., the sp3 or sp2 nature of the bonding orbitals). However, there are other interactions among the several functional groups in arginine that also contribute to its ultimate structure. In particular, the hydrogen bond linking the a-amino group’s nitrogen atom to the –NH- group’s hydrogen atom causes this molecule to fold into a less extended structure than it otherwise might. What does theory have to do with issues of molecular structure and why is knowledge of structure so important? It is important because the structure of a molecule has a very important role in determining the kinds of reactions that molecule will undergo, what kind of radiation it will absorb and emit, and to what “active sites” in neighboring molecules or nearby materials it will bind. A molecule’s shape (e.g., rod like, flat, globular, etc.) is one of the first things a chemist thinks of when trying to predict where, at another molecule or on a surface or a cell, the molecule will “fit” and be able to bind and perhaps react. The presence of lone pairs of electrons (which act as Lewis base sites), of p orbitals (which can act as electron donor and electron acceptor sites), and of highly polar or ionic groups guide the chemist further in determining where on the molecule’s framework various reactant species (e.g., electrophylic or nucleophilic or radical) will be most strongly attracted. Clearly, molecular structure is a crucial aspect of the chemists’ toolbox
How does theory relate to molecular structure? As we discussed in the Background Material, the Born-Oppenheimer approximation leads us to use quantum mechanics to predict the energy E of a molecule for any positions (Ra)of its nuclei given the number of electrons N in the molecule (or ion ). This means, for example, that the energy of the arginine molecule in its lowest electronic state (i.e, with the electrons occupying the lowest energy orbitals)can be determined for any location of the nuclei if the Schrodinger equation governing the movements of the electrons can be solved Ifyou have not had a good class on how quantum mechanics is used within chemistry, I urge you to take the time needed to master the background Material. In those pages, I introduce the central concepts of quantum mechanics and I show how they apply to several very important cases including 1. electrons moving in 1, 2, and 3 dimensions and how these models relate to electronic structures of polyenes and to electronic bands in solids 2. the classical and quantum probability densities and how they differ, 3. time propagation of quantum wave functions, 4. the Huckel or tight-binding model of chemical bonding among atomic orbitals, 5. harmonic vibrations 6. molecular rotations 7. electron tunneling 8. atomic orbitals' angular and radial characteristics, 9. and point group symmetry and how it is used to label orbitals and vibrations
4 How does theory relate to molecular structure? As we discussed in the Background Material, the Born-Oppenheimer approximation leads us to use quantum mechanics to predict the energy E of a molecule for any positions ({Ra}) of its nuclei given the number of electrons Ne in the molecule (or ion). This means, for example, that the energy of the arginine molecule in its lowest electronic state (i.e., with the electrons occupying the lowest energy orbitals) can be determined for any location of the nuclei if the Schrödinger equation governing the movements of the electrons can be solved. If you have not had a good class on how quantum mechanics is used within chemistry, I urge you to take the time needed to master the Background Material. In those pages, I introduce the central concepts of quantum mechanics and I show how they apply to several very important cases including 1. electrons moving in 1, 2, and 3 dimensions and how these models relate to electronic structures of polyenes and to electronic bands in solids 2. the classical and quantum probability densities and how they differ, 3. time propagation of quantum wave functions, 4. the Hückel or tight-binding model of chemical bonding among atomic orbitals, 5. harmonic vibrations, 6. molecular rotations, 7. electron tunneling, 8. atomic orbitals’ angular and radial characteristics, 9. and point group symmetry and how it is used to label orbitals and vibrations
You need to know this material if you wish to understand most of what this text offers,so I urge you to read the background Material if your education to date has not yet d to it Let us now return to the discussion of how theory deals with molecular structure. We assume that we know the energy e(R D at various locations of the nuclei. In some cases, we denote this energy V(R )and in others we use E(R,) because, within the Born Oppenheimer approximation, the electronic energy e serves as the potential v for the molecule's vibrational motions. As discussed in the backgound material. one can then perform a search for the lowest energy structure(. g. by finding where the gradient vector vanishes aE/aR =0 and where the second derivative or Hessian matrix (aE/ar aR,has no negative eigenvalues). By finding such a local-minimum in the energy landscape, theory is able to determine a stable structure of such a molecule. The word stable is used to describe these structures not because they are lower in energy than all other possible arrangements of the atoms but because the curvatures, as given in terms of eigenvalues of the Hessian matrix(OE/OR aR), are positive at this particular geometry. The procedures by which minima on the energy landscape are found may involve simply testing whether the energy decreases or increases as each geometrical coordinate is varied by a small amount. Alternatively, if the gradients aE/aR are known at a particular geometry, one can perform searches directed"downhill"along the negative of the gradient itself. By taking a small"step"along such a direction, one can move to a new geometry that is lower in energy. If not only the gradients aE/aR, but also the second derivatives(OE/aR,aR ) are known at some geometry, one can make a more"intelligent
5 You need to know this material if you wish to understand most of what this text offers, so I urge you to read the Background Material if your education to date has not yet adequately been exposed to it. Let us now return to the discussion of how theory deals with molecular structure. We assume that we know the energy E({Ra}) at various locations {Ra} of the nuclei. In some cases, we denote this energy V(Ra ) and in others we use E(Ra ) because, within the BornOppenheimer approximation, the electronic energy E serves as the potential V for the molecule’s vibrational motions. As discussed in the Backgound Material, one can then perform a search for the lowest energy structure (e.g., by finding where the gradient vector vanishes ¶E/¶Ra = 0 and where the second derivative or Hessian matrix (¶ 2E/¶Ra¶Rb ) has no negative eigenvalues). By finding such a local-minimum in the energy landscape, theory is able to determine a stable structure of such a molecule. The word stable is used to describe these structures not because they are lower in energy than all other possible arrangements of the atoms but because the curvatures, as given in terms of eigenvalues of the Hessian matrix (¶ 2E/¶Ra¶Ra ), are positive at this particular geometry. The procedures by which minima on the energy landscape are found may involve simply testing whether the energy decreases or increases as each geometrical coordinate is varied by a small amount. Alternatively, if the gradients ¶E/¶Ra are known at a particular geometry, one can perform searches directed “downhill” along the negative of the gradient itself. By taking a small “step” along such a direction, one can move to a new geometry that is lower in energy. If not only the gradients ¶E/¶Ra but also the second derivatives (¶ 2E/¶Ra¶Ra ) are known at some geometry, one can make a more “intelligent
step toward a geometry of lower energy For additional details about how such geometry optimization searches are performed within modern computational chemistry software see the background Material where this subject was treated in greater detail It often turns out that a molecule has more than one stable structure(isomer) for a given electronic state. Moreover, the geometries that pertain to stable structures of excited electronic state are different than those obtained for the ground state(because the orbital occupancy and thus the nature of the bonding is different). again using arginine an example, its ground electronic state also has the structure shown in Fig. 5. 2 as a stable isomer. Notice that this isomer and that shown earlier have the atoms linked together in identical manners, but in the second structure the a-amino group is involved in two hydrogen bonds while it is involved in only one in the former. In principle, the relative energies of these two geometrical isomers can be determined by solving the electronic Schrodinger equation while placing the constituent nuclei in the locations described in the two figures 6
6 step toward a geometry of lower energy. For additional details about how such geometry optimization searches are performed within modern computational chemistry software, see the Background Material where this subject was treated in greater detail. It often turns out that a molecule has more than one stable structure (isomer) for a given electronic state. Moreover, the geometries that pertain to stable structures of excited electronic state are different than those obtained for the ground state (because the orbital occupancy and thus the nature of the bonding is different). Again using arginine as an example, its ground electronic state also has the structure shown in Fig. 5.2 as a stable isomer. Notice that this isomer and that shown earlier have the atoms linked together in identical manners, but in the second structure the a-amino group is involved in two hydrogen bonds while it is involved in only one in the former. In principle, the relative energies of these two geometrical isomers can be determined by solving the electronic Schrödinger equation while placing the constituent nuclei in the locations described in the two figures
1.916 2 2144 Figure 5.2 Another stable structure for the arginine molecule If the arginine molecule is excited to another electronic state, for example, by promoting a non-bonding electron on its C=O oxygen atom into the neighboring C-O T' orbital. its stable structures will not be the same as in the ground electronic state. In particular, the corresponding C-o distance will be longer than in the ground state, but other internal geometrical parameters may also be modified (albeit probably less so than the C-o distance). Moreover, the chemical reactivity of this excited state of arginine will be different than that of the ground state because the two states have different orbitals available to react with attacking reagents
7 1.916 2.144 Figure 5.2 Another stable structure for the arginine molecule. If the arginine molecule is excited to another electronic state, for example, by promoting a non-bonding electron on its C=O oxygen atom into the neighboring C-O p* orbital, its stable structures will not be the same as in the ground electronic state. In particular, the corresponding C-O distance will be longer than in the ground state, but other internal geometrical parameters may also be modified (albeit probably less so than the C-O distance). Moreover, the chemical reactivity of this excited state of arginine will be different than that of the ground state because the two states have different orbitals available to react with attacking reagents
In summary, by solving the electronic Schrodinger equation at a variety of geometries and searching for geometries where the gradient vanishes and the Hessian matrix has all positive eigenvalues, one can find stable structures of molecules(and ions). The Schrodinger equation is a necessary aspect of this process because the movement of the electrons is governed by this equation rather than by Newtonian classical equations. The information gained after carrying out such a geometry optimization process include(1) all of the interatomic distances and internal angles needed to specify the equilibrium geometry(Rae) and(2)the total electronic energy E at this particular geometry It is also possible to extract much more information from these calculations. For example, by multiplying elements of the Hessian matrix(oE/OR, aR, by the inverse quare roots of the atomic masses of the atoms labeled a and b, one forms the mass- eighted Hessian(m, m)(aE/aR, OR,)whose non-zero eigenvalues give the harmonic vibrational frequencies(o, of the molecule. The eigenvectors (Rka of the mass wieghted Hessian mantrix give the relative displacements in coordinates r a that accompany vibration in the k normal mode (i.e, they describe the normal mode motions). Details about how these harmonic vibrational frequencies and normal modes are obtained were discussed earlier in the Background Material B. Molecular Change-reactions, isomerization, interactions Changes in bonding 8
8 In summary, by solving the electronic Schrödinger equation at a variety of geometries and searching for geometries where the gradient vanishes and the Hessian matrix has all positive eigenvalues, one can find stable structures of molecules (and ions). The Schrödinger equation is a necessary aspect of this process because the movement of the electrons is governed by this equation rather than by Newtonian classical equations. The information gained after carrying out such a geometry optimization process include (1) all of the interatomic distances and internal angles needed to specify the equilibrium geometry {Raeq} and (2) the total electronic energy E at this particular geometry. It is also possible to extract much more information from these calculations. For example, by multiplying elements of the Hessian matrix (¶ 2E/¶Ra¶Rb ) by the inverse square roots of the atomic masses of the atoms labeled a and b, one forms the massweighted Hessian (ma mb ) -1/2 (¶ 2E/¶Ra¶Rb ) whose non-zero eigenvalues give the harmonic vibrational frequencies {wk} of the molecule. The eigenvectors {Rk,a} of the masswieghted Hessian mantrix give the relative displacements in coordinates Rka that accompany vibration in the kth normal mode (i.e., they describe the normal mode motions). Details about how these harmonic vibrational frequencies and normal modes are obtained were discussed earlier in the Background Material. B. Molecular Change- reactions, isomerization, interactions 1. Changes in bonding
Chemistry also deals with transformations of matter including changes that occur when cues react are ex cited(electronically, vibrationally, or rotationally ) or undergo geometrical rearrangements. Again, theory forms the cornerstone that allows experimental probes of chemical change to be connected to the molecular level and that allows simulations of such changes Molecular excitation may or may not involve altering the electronic structure of the molecule. vibrational and rotational excitation do not but electronic excitation ionization, and electron attachment do. as illustrated in Fig. 5. 3 where a bi-molecular reaction is displayed, chemical reactions involve breaking some bonds and forming others and thus involve rearrangement of the electrons among various molecular orbitals Figure 5.3 Two bimolecular reactions a and b show an atom combining with a diatomic c and d show an atom abstracting an atom from a diatomic
9 Chemistry also deals with transformations of matter including changes that occur when molecules react, are excited (electronically, vibrationally, or rotationally), or undergo geometrical rearrangements. Again, theory forms the cornerstone that allows experimental probes of chemical change to be connected to the molecular level and that allows simulations of such changes. Molecular excitation may or may not involve altering the electronic structure of the molecule; vibrational and rotational excitation do not, but electronic excitation, ionization, and electron attachment do. As illustrated in Fig. 5.3 where a bi-molecular reaction is displayed, chemical reactions involve breaking some bonds and forming others, and thus involve rearrangement of the electrons among various molecular orbitals. Figure 5.3 Two bimolecular reactions; a and b show an atom combining with a diatomic; c and d show an atom abstracting an atom from a diatomic
In this example, in part(a) the green atom collides with the brown diatomic molecule and forms the bound triatomic(b). Alternatively, in(c)and(d), a pink atom collides with a green diatomic to break the bond between the two green atoms and form a new bond between the pink and green atoms. Both such reactions are termed bi-molecular because the basic step in which the reaction takes place requires a collision between to independent species(i.e, the atom and the diatomic) A simple example of a unimolecular chemical reaction is offered by the arginine molecule considered above. In the first structure shown for arginine, the carboxylic acid group retains its HOOC-bonding. However, in the zwitterion structure of this same molecule, shown in Fig. 5.4, the HOOC- group has been deprotonated to produce a carboxy late anion group-Coo, with the H ion now bonded to the terminal imine group, thus converting it to an amino group and placing the net positive charge on the adjacer carbon atom The unimolecular tautomerization reaction in which the two forms of arginine are interconverted involves breaking an O-H bond, forming a N-H bond, and changing a carbon-nitrogen double bond into a carbon-nitrogen single bond. In such a process, the electronic structure is significantly altered, and, as a result, the two isomers can display very different chemical reactivities toward other reagents. Notice that, once again, the ultimate structure of the zwitterion tautomer of arganine is determined by the valence preferences of its constitutent atoms as well as by hydrogen bonds formed among various functional groups(the carboxylate group and one amino group and one-NH
10 In this example, in part (a) the green atom collides with the brown diatomic molecule and forms the bound triatomic (b). Alternatively, in (c) and (d), a pink atom collides with a green diatomic to break the bond between the two green atoms and form a new bond between the pink and green atoms. Both such reactions are termed bi-molecular because the basic step in which the reaction takes place requires a collision between to independent species (i.e., the atom and the diatomic). A simple example of a unimolecular chemical reaction is offered by the arginine molecule considered above. In the first structure shown for arginine, the carboxylic acid group retains its HOOC- bonding. However, in the zwitterion structure of this same molecule, shown in Fig. 5.4, the HOOC- group has been deprotonated to produce a carboxylate anion group –COO- , with the H+ ion now bonded to the terminal imine group, thus converting it to an amino group and placing the net positive charge on the adjacent carbon atom. The unimolecular tautomerization reaction in which the two forms of arginine are interconverted involves breaking an O-H bond, forming a N-H bond, and changing a carbon-nitrogen double bond into a carbon-nitrogen single bond. In such a process, the electronic structure is significantly altered, and, as a result, the two isomers can display very different chemical reactivities toward other reagents. Notice that, once again, the ultimate structure of the zwitterion tautomer of arganine is determined by the valence preferences of its constitutent atoms as well as by hydrogen bonds formed among various functional groups (the carboxylate group and one amino group and one –NHgroup)