Atoms in molecules

Last updated

In quantum chemistry, the quantum theory of atoms in molecules (QTAIM), sometimes referred to as atoms in molecules (AIM), is a model of molecular and condensed matter electronic systems (such as crystals) in which the principal objects of molecular structure - atoms and bonds - are natural expressions of a system's observable electron density distribution function. An electron density distribution of a molecule is a probability distribution that describes the average manner in which the electronic charge is distributed throughout real space in the attractive field exerted by the nuclei. According to QTAIM, molecular structure is revealed by the stationary points of the electron density together with the gradient paths of the electron density that originate and terminate at these points.

Contents

QTAIM was primarily developed by Professor Richard Bader and his research group at McMaster University over the course of decades, beginning with analyses of theoretically calculated electron densities of simple molecules in the early 1960s and culminating with analyses of both theoretically and experimentally measured electron densities of crystals in the 90s. The development of QTAIM was driven by the assumption that, since the concepts of atoms and bonds have been and continue to be so ubiquitously useful in interpreting, classifying, predicting and communicating chemistry, they should have a well-defined physical basis.

QTAIM recovers the central operational concepts of the molecular structure hypothesis, that of a functional grouping of atoms with an additive and characteristic set of properties, together with a definition of the bonds that link the atoms and impart the structure. QTAIM defines chemical bonding and structure of a chemical system based on the topology of the electron density. In addition to bonding, QTAIM allows the calculation of certain physical properties on a per-atom basis, by dividing space up into atomic volumes containing exactly one nucleus, which acts as a local attractor of the electron density. In QTAIM an atom is defined as a proper open system, i.e. a system that can share energy and electron density which is localized in the 3D space. The mathematical study of these features is usually referred to in the literature as charge density topology.

QTAIM rests on the fact that the dominant topological property of the vast majority of electron density distributions is the presence of strong maxima that occur exclusively at the nuclei, certain pairs of which are linked together by ridges of electron density. In terms of an electron density distribution's gradient vector field, this corresponds to a complete, non-overlapping partitioning of a molecule into three-dimensional basins (atoms) that are linked together by shared two-dimensional separatrices (interatomic surfaces). Within each interatomic surface, the electron density is a maximum at the corresponding internuclear saddle point, which also lies at the minimum of the ridge between corresponding pair of nuclei, the ridge being defined by the pair of gradient trajectories (bond path) originating at the saddle point and terminating at the nuclei. Because QTAIM atoms are always bounded by surfaces having zero flux in the gradient vector field of the electron density, they have some unique quantum mechanical properties compared to other subsystem definitions, including unique electronic kinetic energy the satisfaction of an electronic virial theorem analogous to the molecular electronic virial theorem and some interesting variational properties. QTAIM has gradually become a method for addressing possible questions regarding chemical systems, in a variety of situations hardly handled before by any other model or theory in chemistry. [1] [2] [3] [4]

Applications

QTAIM is applied to the description of certain organic crystals with unusually short distances between neighboring molecules as observed by X-ray diffraction. For example, in the crystal structure of molecular chlorine the experimental Cl...Cl distance between two molecules is 327 picometres which is less than the sum of the van der Waals radii of 350 picometres. In one QTAIM result 12 bond paths start from each chlorine atom to other chlorine atoms including the other chlorine atom in the molecule. The theory also aims to explain the metallic properties of metallic hydrogen in much the same way.

The theory is also applied to so-called hydrogen-hydrogen bonds [5] as they occur in molecules such as phenanthrene and chrysene. In these compounds the distance between two ortho hydrogen atoms again is shorter than their van der Waals radii and according to in silico experiments based on this theory, a bond path is identified between them. Both hydrogen atoms have identical electron density and are closed shell and therefore they are very different from the so-called dihydrogen bonds which are postulated for compounds such as (CH3)2NHBH3 and also different from so-called agostic interactions.

In mainstream chemistry close proximity of two nonbonding atoms leads to destabilizing steric repulsion but in QTAIM the observed hydrogen hydrogen interactions are in fact stabilizing. It is well known that both kinked phenanthrene and chrysene are around 6 kcal/mol (25 kJ/mol) more stable than their linear isomers anthracene and tetracene. One traditional explanation is given by Clar's rule. QTAIM shows that a calculated stabilization for phenanthrene by 8 kcal/mol (33 kJ/mol) is the result of destabilization of the compound by 8 kcal/mol (33 kJ/mol) originating from electron transfer from carbon to hydrogen, offset by 12.1 kcal (51 kJ/mol) of stabilization due to a H..H bond path. The electron density at the critical point between the two hydrogen atoms is low, 0.012 e for phenanthrene. Another property of the bond path is its curvature.

Biphenyl, phenanthrene and anthracene HydrogenhydrogenbondingE.svg
Biphenyl, phenanthrene and anthracene

Another molecule studied in QTAIM is biphenyl. Its two phenyl rings are oriented at a 38° angle with respect to each other with the planar molecular geometry (encountered in a rotation around the central C-C bond) destabilized by 2.1 kcal/mol (8.8 kJ/mol) and the perpendicular one destabilized by 2.5 kcal/mol (10.5 kJ/mol). The classic explanations for this rotation barrier are steric repulsion between the ortho-hydrogen atoms (planar) and breaking of delocalization of pi density over both rings (perpendicular).

QTAIM has also been applied to study the electron topology of solvated post-translational modifications to protein. For example, covalently bonded force constants in a set of lysine-arginine derived advanced glycation end-products were derived using the electronic structure calculations and then bond paths were used to illustrate differences in each of the applied computational chemistry functionals. [6] Furthermore, QTAIM had been used to identify a bond path network of hydrogen bonds between glucosepane and nearby water molecules. [7]

In QTAIM the energy increase on decreasing the dihedral angle from 38° to 0° is a summation of several factors. Destabilizing factors are the increase in bond length between the connecting carbon atoms (because they have to accommodate the approaching hydrogen atoms) and transfer of electronic charge from carbon to hydrogen. Stabilizing factors are increased delocalization of pi-electrons from one ring to the other and the one that tips the balance is a hydrogen - hydrogen bond between the ortho hydrogens.

The hydrogen bond is not without its critics. According to one, the relative stability of phenanthrene compared to its isomers can be adequately explained by comparing resonance stabilizations. [8] Another critic [9] argues that the stability of phenanthrene can be attributed to more effective pi-pi overlap in the central double bond; the existence of bond paths are not questioned but the stabilizing energy derived from it is.

See also

Related Research Articles

<span class="mw-page-title-main">Chemical bond</span> Lasting attraction between atoms that enables the formation of chemical compounds

A chemical bond is a lasting attraction between atoms or ions that enables the formation of molecules, crystals, and other structures. The bond may result from the electrostatic force between oppositely charged ions as in ionic bonds, or through the sharing of electrons as in covalent bonds. The strength of chemical bonds varies considerably; there are "strong bonds" or "primary bonds" such as covalent, ionic and metallic bonds, and "weak bonds" or "secondary bonds" such as dipole–dipole interactions, the London dispersion force, and hydrogen bonding.

<span class="mw-page-title-main">Covalent bond</span> Chemical bond by sharing of electron pairs

A covalent bond is a chemical bond that involves the sharing of electrons to form electron pairs between atoms. These electron pairs are known as shared pairs or bonding pairs. The stable balance of attractive and repulsive forces between atoms, when they share electrons, is known as covalent bonding. For many molecules, the sharing of electrons allows each atom to attain the equivalent of a full valence shell, corresponding to a stable electronic configuration. In organic chemistry, covalent bonding is much more common than ionic bonding.

<span class="mw-page-title-main">Hydrogen bond</span> Intermolecular attraction between a hydrogen-donor pair and an acceptor

In chemistry, a hydrogen bond is primarily an electrostatic force of attraction between a hydrogen (H) atom which is covalently bonded to a more electronegative "donor" atom or group (Dn), and another electronegative atom bearing a lone pair of electrons—the hydrogen bond acceptor (Ac). Such an interacting system is generally denoted Dn−H···Ac, where the solid line denotes a polar covalent bond, and the dotted or dashed line indicates the hydrogen bond. The most frequent donor and acceptor atoms are the period 2 elements nitrogen (N), oxygen (O), and fluorine (F).

<span class="mw-page-title-main">Ionic bonding</span> Chemical bonding involving attraction between ions

Ionic bonding is a type of chemical bonding that involves the electrostatic attraction between oppositely charged ions, or between two atoms with sharply different electronegativities, and is the primary interaction occurring in ionic compounds. It is one of the main types of bonding, along with covalent bonding and metallic bonding. Ions are atoms with an electrostatic charge. Atoms that gain electrons make negatively charged ions. Atoms that lose electrons make positively charged ions. This transfer of electrons is known as electrovalence in contrast to covalence. In the simplest case, the cation is a metal atom and the anion is a nonmetal atom, but these ions can be more complex, e.g. molecular ions like NH+
4 or SO2−
4. In simpler words, an ionic bond results from the transfer of electrons from a metal to a non-metal to obtain a full valence shell for both atoms.

Quantum chemistry, also called molecular quantum mechanics, is a branch of physical chemistry focused on the application of quantum mechanics to chemical systems, particularly towards the quantum-mechanical calculation of electronic contributions to physical and chemical properties of molecules, materials, and solutions at the atomic level. These calculations include systematically applied approximations intended to make calculations computationally feasible while still capturing as much information about important contributions to the computed wave functions as well as to observable properties such as structures, spectra, and thermodynamic properties. Quantum chemistry is also concerned with the computation of quantum effects on molecular dynamics and chemical kinetics.

<span class="mw-page-title-main">Theoretical chemistry</span> Branch of chemistry

Theoretical chemistry is the branch of chemistry which develops theoretical generalizations that are part of the theoretical arsenal of modern chemistry: for example, the concepts of chemical bonding, chemical reaction, valence, the surface of potential energy, molecular orbitals, orbital interactions, and molecule activation.

<span class="mw-page-title-main">Conjugated system</span> System of connected p-orbitals with delocalized electrons in a molecule

In theoretical chemistry, a conjugated system is a system of connected p-orbitals with delocalized electrons in a molecule, which in general lowers the overall energy of the molecule and increases stability. It is conventionally represented as having alternating single and multiple bonds. Lone pairs, radicals or carbenium ions may be part of the system, which may be cyclic, acyclic, linear or mixed. The term "conjugated" was coined in 1899 by the German chemist Johannes Thiele.

<span class="mw-page-title-main">Octet rule</span> Chemical rule of thumb

The octet rule is a chemical rule of thumb that reflects the theory that main-group elements tend to bond in such a way that each atom has eight electrons in its valence shell, giving it the same electronic configuration as a noble gas. The rule is especially applicable to carbon, nitrogen, oxygen, and the halogens; although more generally the rule is applicable for the s-block and p-block of the periodic table. Other rules exist for other elements, such as the duplet rule for hydrogen and helium, and the 18-electron rule for transition metals.

In chemistry, resonance, also called mesomerism, is a way of describing bonding in certain molecules or polyatomic ions by the combination of several contributing structures into a resonance hybrid in valence bond theory. It has particular value for analyzing delocalized electrons where the bonding cannot be expressed by one single Lewis structure. It is considered the accurate structure for a compound.

In chemistry, molecular orbital theory is a method for describing the electronic structure of molecules using quantum mechanics. It was proposed early in the 20th century.

In chemistry, valence bond (VB) theory is one of the two basic theories, along with molecular orbital (MO) theory, that were developed to use the methods of quantum mechanics to explain chemical bonding. It focuses on how the atomic orbitals of the dissociated atoms combine to give individual chemical bonds when a molecule is formed. In contrast, molecular orbital theory has orbitals that cover the whole molecule.

<span class="mw-page-title-main">Molecular geometry</span> Study of the 3D shapes of molecules

Molecular geometry is the three-dimensional arrangement of the atoms that constitute a molecule. It includes the general shape of the molecule as well as bond lengths, bond angles, torsional angles and any other geometrical parameters that determine the position of each atom.

<span class="mw-page-title-main">VSEPR theory</span> Model for predicting molecular geometry

Valence shell electron pair repulsion (VSEPR) theory, is a model used in chemistry to predict the geometry of individual molecules from the number of electron pairs surrounding their central atoms. It is also named the Gillespie-Nyholm theory after its two main developers, Ronald Gillespie and Ronald Nyholm.

<span class="mw-page-title-main">Hyperconjugation</span> Concept in organic chemistry

In organic chemistry, hyperconjugation refers to the delocalization of electrons with the participation of bonds of primarily σ-character. Usually, hyperconjugation involves the interaction of the electrons in a sigma (σ) orbital with an adjacent unpopulated non-bonding p or antibonding σ* or π* orbitals to give a pair of extended molecular orbitals. However, sometimes, low-lying antibonding σ* orbitals may also interact with filled orbitals of lone pair character (n) in what is termed negative hyperconjugation. Increased electron delocalization associated with hyperconjugation increases the stability of the system. In particular, the new orbital with bonding character is stabilized, resulting in an overall stabilization of the molecule. Only electrons in bonds that are in the β position can have this sort of direct stabilizing effect — donating from a sigma bond on an atom to an orbital in another atom directly attached to it. However, extended versions of hyperconjugation can be important as well. The Baker–Nathan effect, sometimes used synonymously for hyperconjugation, is a specific application of it to certain chemical reactions or types of structures.

A molecular orbital diagram, or MO diagram, is a qualitative descriptive tool explaining chemical bonding in molecules in terms of molecular orbital theory in general and the linear combination of atomic orbitals (LCAO) method in particular. A fundamental principle of these theories is that as atoms bond to form molecules, a certain number of atomic orbitals combine to form the same number of molecular orbitals, although the electrons involved may be redistributed among the orbitals. This tool is very well suited for simple diatomic molecules such as dihydrogen, dioxygen, and carbon monoxide but becomes more complex when discussing even comparatively simple polyatomic molecules, such as methane. MO diagrams can explain why some molecules exist and others do not. They can also predict bond strength, as well as the electronic transitions that can take place.

In chemistry, a halogen bond occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a halogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. Like a hydrogen bond, the result is not a formal chemical bond, but rather a strong electrostatic attraction. Mathematically, the interaction can be decomposed in two terms: one describing an electrostatic, orbital-mixing charge-transfer and another describing electron-cloud dispersion. Halogen bonds find application in supramolecular chemistry; drug design and biochemistry; crystal engineering and liquid crystals; and organic catalysis.

<span class="mw-page-title-main">Multipole density formalism</span>

The Multipole Density Formalism is an X-ray crystallography method of electron density modelling proposed by Niels K. Hansen and Philip Coppens in 1978. Unlike the commonly used Independent Atom Model, the Hansen-Coppens Formalism presents an aspherical approach, allowing one to model the electron distribution around a nucleus separately in different directions and therefore describe numerous chemical features of a molecule inside the unit cell of an examined crystal in detail.

In theoretical chemistry, the bonding orbital is used in molecular orbital (MO) theory to describe the attractive interactions between the atomic orbitals of two or more atoms in a molecule. In MO theory, electrons are portrayed to move in waves. When more than one of these waves come close together, the in-phase combination of these waves produces an interaction that leads to a species that is greatly stabilized. The result of the waves’ constructive interference causes the density of the electrons to be found within the binding region, creating a stable bond between the two species.

<span class="mw-page-title-main">Linnett double-quartet theory</span>

Linnett double-quartet theory (LDQ) is a method of describing the bonding in molecules which involves separating the electrons depending on their spin, placing them into separate 'spin tetrahedra' to minimise the Pauli repulsions between electrons of the same spin. Introduced by J. W. Linnett in his 1961 monograph and 1964 book, this method expands on the electron dot structures pioneered by G. N. Lewis. While the theory retains the requirement for fulfilling the octet rule, it dispenses with the need to force electrons into coincident pairs. Instead, the theory stipulates that the four electrons of a given spin should maximise the distances between each other, resulting in a net tetrahedral electronic arrangement that is the fundamental molecular building block of the theory.

Dispersion stabilized molecules are molecules where the London dispersion force (LDF), a non-covalent attractive force between atoms and molecules, plays a significant role in promoting the molecule's stability. Distinct from steric hindrance, dispersion stabilization has only recently been considered in depth by organic and inorganic chemists after earlier gaining prominence in protein science and supramolecular chemistry. Although usually weaker than covalent bonding and other forms of non-covalent interactions like hydrogen bonding, dispersion forces are known to be a significant if not dominating stabilizing force in certain organic, inorganic, and main group molecules, stabilizing otherwise reactive moieties and exotic bonding.

References

  1. ↑ Bader, Richard (1994). Atoms in Molecules: A Quantum Theory. USA: Oxford University Press. ISBN   978-0-19-855865-1.
  2. ↑ Bader, R. (1991). "A quantum theory of molecular structure and its applications". Chemical Reviews. 91 (5): 893–928. doi:10.1021/cr00005a013.
  3. ↑ Bader, R.F.W. (2005). "The Quantum Mechanical Basis for Conceptual Chemistry". Monatshefte für Chemie. 136 (6): 819–854. doi:10.1007/s00706-005-0307-x. S2CID   121874327.
  4. ↑ Bader, R.F.W. (1998). "Atoms in Molecules". Encyclopedia of Computational Chemistry. 1: 64–86.
  5. ↑ Matta, Chérif F.; Hernández-Trujillo, Jesús; Tang, Ting-Hua; Bader, Richard F. W. (2003). "Hydrogen–Hydrogen Bonding: A Stabilizing Interaction in Molecules and Crystals". Chemistry - A European Journal. 9 (9): 1940–1951. doi:10.1002/chem.200204626. PMID   12740840.
  6. ↑ Nash, A., de Leeuw, N. H., Birch, H. L. (2018). "Bonded Force Constant Derivation of Lysine-Arginine Cross-linked Advanced Glycation End-Products". ChemRxiv.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  7. ↑ Nash, Anthony; Saßmannshausen, Jörg; Bozec, Laurent; Birch, Helen L.; De Leeuw, Nora H. (2017). "Computational study of glucosepane-water hydrogen bond formation: an electron topology and orbital analysis". Journal of Biomolecular Structure and Dynamics. 35 (5): 1127–1137. doi: 10.1080/07391102.2016.1172026 . PMID   27092586.
  8. ↑ Dunitz, Jack D.; Gavezzotti, Angelo (2005). "Molecular Recognition in Organic Crystals: Directed Intermolecular Bonds or Nonlocalized Bonding?". Angewandte Chemie International Edition. 44 (12): 1766–1787. doi:10.1002/anie.200460157. PMID   15685679.
  9. ↑ Poater, Jordi; Visser, Ruud; Solà, Miquel; Bickelhaupt, F. Matthias (2007). "Polycyclic Benzenoids: Why Kinked is More Stable than Straight". The Journal of Organic Chemistry. 72 (4): 1134–1142. Bibcode:2007JOCh...72.1134P. doi:10.1021/jo061637p. PMID   17288368.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.