Valence shell electron pair repulsion theory, or VSEPR theory (/ /, VESP-ər,:410 və-SEP-ər), 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.
The premise of VSEPR is that the valence electron pairs surrounding an atom tend to repel each other and will, therefore, adopt an arrangement that minimizes this repulsion. This in turn decreases the molecule's energy and increases its stability, which determines the molecular geometry. Gillespie has emphasized that the electron-electron repulsion due to the Pauli exclusion principle is more important in determining molecular geometry than the electrostatic repulsion.
The insights of VSEPR theory are derived from topological analysis of the electron density of molecules. Such quantum chemical topology (QCT) methods include the electron localization function (ELF) and the quantum theory of atoms in molecules (AIM or QTAIM). Hence, VSEPR is unrelated to wave function based methods such as orbital hybridisation in valence bond theory.
The idea of a correlation between molecular geometry and number of valence electron pairs (both shared and unshared pairs) was originally proposed in 1939 by Ryutaro Tsuchida in Japan, and was independently presented in a Bakerian Lecture in 1940 by Nevil Sidgwick and Herbert Powell of the University of Oxford. In 1957, Ronald Gillespie and Ronald Sydney Nyholm of University College London refined this concept into a more detailed theory, capable of choosing between various alternative geometries.
VSEPR theory is used to predict the arrangement of electron pairs around central atoms in molecules, especially simple and symmetric molecules. A central atom is defined in this theory as an atom which is bonded to two or more other atoms, while a terminal atom is bonded to only one other atom.:398 For example in the molecule methyl isocyanate (H3C-N=C=O), the two carbons and one nitrogen are central atoms, and the three hydrogens and one oxygen are terminal atoms.:416 The geometry of the central atoms and their non-bonding electron pairs in turn determine the geometry of the larger whole molecule.
The number of electron pairs in the valence shell of a central atom is determined after drawing the Lewis structure of the molecule, and expanding it to show all bonding groups and lone pairs of electrons.:410–417 In VSEPR theory, a double bond or triple bond is treated as a single bonding group. The sum of the number of atoms bonded to a central atom and the number of lone pairs formed by its nonbonding valence electrons is known as the central atom's steric number.
The electron pairs (or groups if multiple bonds are present) are assumed to lie on the surface of a sphere centered on the central atom and tend to occupy positions that minimize their mutual repulsions by maximizing the distance between them.:410–417 The number of electron pairs (or groups), therefore, determines the overall geometry that they will adopt. For example, when there are two electron pairs surrounding the central atom, their mutual repulsion is minimal when they lie at opposite poles of the sphere. Therefore, the central atom is predicted to adopt a linear geometry. If there are 3 electron pairs surrounding the central atom, their repulsion is minimized by placing them at the vertices of an equilateral triangle centered on the atom. Therefore, the predicted geometry is trigonal. Likewise, for 4 electron pairs, the optimal arrangement is tetrahedral.:410–417
As a tool in predicting the geometry adopted with a given number of electron pairs, an often used physical demonstration of the principle of minimal electron pair repulsion utilizes inflated balloons. Through handling, balloons acquire a slight surface electrostatic charge that results in the adoption of roughly the same geometries when they are tied together at their stems as the corresponding number of electron pairs. For example, five balloons tied together adopt the trigonal bipyramidal geometry, just as do the five bonding pairs of a PCl5 molecule.
The steric number of a central atom in a molecule is the number of atoms bonded to that central atom, called its coordination number, plus the number of lone pairs of valence electrons on the central atom. In the molecule SF4, for example, the central sulfur atom has four ligands; the coordination number of sulfur is four. In addition to the four ligands, sulfur also has one lone pair in this molecule. Thus, the steric number is 4 + 1 = 5.
Degree of repulsion
The overall geometry is further refined by distinguishing between bonding and nonbonding electron pairs. The bonding electron pair shared in a sigma bond with an adjacent atom lies further from the central atom than a nonbonding (lone) pair of that atom, which is held close to its positively charged nucleus. VSEPR theory therefore views repulsion by the lone pair to be greater than the repulsion by a bonding pair. As such, when a molecule has 2 interactions with different degrees of repulsion, VSEPR theory predicts the structure where lone pairs occupy positions that allow them to experience less repulsion. Lone pair–lone pair (lp–lp) repulsions are considered stronger than lone pair–bonding pair (lp–bp) repulsions, which in turn are considered stronger than bonding pair–bonding pair (bp–bp) repulsions, distinctions that then guide decisions about overall geometry when 2 or more non-equivalent positions are possible.:410–417 For instance, when 5 valence electron pairs surround a central atom, they adopt a trigonal bipyramidal molecular geometry with two collinear axial positions and three equatorial positions. An electron pair in an axial position has three close equatorial neighbors only 90° away and a fourth much farther at 180°, while an equatorial electron pair has only two adjacent pairs at 90° and two at 120°. The repulsion from the close neighbors at 90° is more important, so that the axial positions experience more repulsion than the equatorial positions; hence, when there are lone pairs, they tend to occupy equatorial positions as shown in the diagrams of the next section for steric number five.
The difference between lone pairs and bonding pairs may also be used to rationalize deviations from idealized geometries. For example, the H2O molecule has four electron pairs in its valence shell: two lone pairs and two bond pairs. The four electron pairs are spread so as to point roughly towards the apices of a tetrahedron. However, the bond angle between the two O–H bonds is only 104.5°, rather than the 109.5° of a regular tetrahedron, because the two lone pairs (whose density or probability envelopes lie closer to the oxygen nucleus) exert a greater mutual repulsion than the two bond pairs.:410–417
A bond of higher bond order also exerts greater repulsion since the pi bond electrons contribute. For example in isobutylene, (H3C)2C=CH2, the H3C−C=C angle (124°) is larger than the H3C−C−CH3 angle (111.5°). However, in the carbonate ion, CO2−
3, all three C−O bonds are equivalent with angles of 120° due to resonance.
The "AXE method" of electron counting is commonly used when applying the VSEPR theory. The electron pairs around a central atom are represented by a formula AXnEm, where A represents the central atom and always has an implied subscript one. Each X represents a ligand (an atom bonded to A). Each E represents a lone pair of electrons on the central atom.:410–417 The total number of X and E is known as the steric number. For example in a molecule AX3E2, the atom A has a steric number of 5.
When the substituent (X) atoms are not all the same, the geometry is still approximately valid, but the bond angles may be slightly different from the ones where all the outside atoms are the same. For example, the double-bond carbons in alkenes like C2H4 are AX3E0, but the bond angles are not all exactly 120°. Likewise, SOCl2 is AX3E1, but because the X substituents are not identical, the X–A–X angles are not all equal.
Based on the steric number and distribution of Xs and Es, VSEPR theory makes the predictions in the following tables.
For main-group elements, there are stereochemically active lone pairs E whose number can vary between 0 to 3. Note that the geometries are named according to the atomic positions only and not the electron arrangement. For example, the description of AX2E1 as a bent molecule means that the three atoms AX2 are not in one straight line, although the lone pair helps to determine the geometry.
0 lone pairs
1 lone pair
2 lone pairs
3 lone pairs
|3|| || |
|4|| || || |
|5|| || || || |
|6|| || || |
|7|| || || |
including lone pairs, shown in pale yellow
excluding lone pairs
2, SO2,:413–414 O3, CCl2
|AX3E0||Trigonal planar||BF3,:413–414 CO2−|
|AX3E1||Trigonal pyramidal||NH3,:413–414 PCl3:407|
|AX4E1||Seesaw or disphenoidal||SF4:413–414:45|
|AX5E1||Square pyramidal||ClF5,:481 BrF5,:413–414 XeOF4|
Transition metals (Kepert model)
The lone pairs on transition metal atoms are usually stereochemically inactive, meaning that their presence does not change the molecular geometry. For example, the hexaaquo complexes M(H2O)6 are all octahedral for M = V3+, Mn3+, Co3+, Ni2+ and Zn2+, despite the fact that the electronic configurations of the central metal ion are d2, d4, d6, d8 and d10 respectively.:542 The Kepert model ignores all lone pairs on transition metal atoms, so that the geometry around all such atoms corresponds to the VSEPR geometry for AXn with 0 lone pairs E.:542 This is often written MLn, where M = metal and L = ligand. The Kepert model predicts the following geometries for coordination numbers of 2 through 9:
|Capped trigonal prismatic||TaF2−|
|Bicapped trigonal prismatic||ZrF4−|
|AX9||Tricapped trigonal prismatic||ReH2−|
|Capped square antiprismatic|
The methane molecule (CH4) is tetrahedral because there are four pairs of electrons. The four hydrogen atoms are positioned at the vertices of a tetrahedron, and the bond angle is cos−1(−1⁄3) ≈ 109° 28′. This is referred to as an AX4 type of molecule. As mentioned above, A represents the central atom and X represents an outer atom.:410–417
The ammonia molecule (NH3) has three pairs of electrons involved in bonding, but there is a lone pair of electrons on the nitrogen atom.:392–393 It is not bonded with another atom; however, it influences the overall shape through repulsions. As in methane above, there are four regions of electron density. Therefore, the overall orientation of the regions of electron density is tetrahedral. On the other hand, there are only three outer atoms. This is referred to as an AX3E type molecule because the lone pair is represented by an E.:410–417 By definition, the molecular shape or geometry describes the geometric arrangement of the atomic nuclei only, which is trigonal-pyramidal for NH3.:410–417
Steric numbers of 7 or greater are possible, but are less common. The steric number of 7 occurs in iodine heptafluoride (IF7); the base geometry for a steric number of 7 is pentagonal bipyramidal. The most common geometry for a steric number of 8 is a square antiprismatic geometry.:1165 Examples of this include the octacyanomolybdate (Mo(CN)4−
8) and octafluorozirconate (ZrF4−
8) anions.:1165 The nonahydridorhenate ion (ReH2−
9) in potassium nonahydridorhenate is a rare example of a compound with a steric number of 9, which has a tricapped trigonal prismatic geometry.:254
Possible geometries for steric numbers of 10, 11, 12, or 14 are bicapped square antiprismatic (or bicapped dodecadeltahedral), octadecahedral, icosahedral, and bicapped hexagonal antiprismatic, respectively. No compounds with steric numbers this high involving monodentate ligands exist, and those involving multidentate ligands can often be analysed more simply as complexes with lower steric numbers when some multidentate ligands are treated as a unit.:1165,1721
There are groups of compounds where VSEPR fails to predict the correct geometry.
Some AX2E0 molecules
Some AX2E2 molecules
One example of the AX2E2 geometry is molecular lithium oxide, Li2O, a linear rather than bent structure, which is ascribed to its bonds being essentially ionic and the strong lithium-lithium repulsion that results. Another example is O(SiH3)2 with an Si–O–Si angle of 144.1°, which compares to the angles in Cl2O (110.9°), (CH3)2O (111.7°), and N(CH3)3 (110.9°). Gillespie and Robinson rationalize the Si–O–Si bond angle based on the observed ability of a ligand's lone pair to most greatly repel other electron pairs when the ligand electronegativity is greater than or equal to that of the central atom. In O(SiH3)2, the central atom is more electronegative, and the lone pairs are less localized and more weakly repulsive. The larger Si–O–Si bond angle results from this and strong ligand-ligand repulsion by the relatively large -SiH3 ligand. Burford et al showed through X-ray diffraction studies that Cl3Al–O–PCl3 has a linear Al–O–P bond angle and is therefore a non-VSEPR molecule.
Some AX6E1 and AX8E1 molecules
Some AX6E1 molecules, e.g. xenon hexafluoride (XeF6) and the Te(IV) and Bi(III) anions, TeCl2−
6 and BiI3−
6, are octahedra, rather than pentagonal pyramids, and the lone pair does not affect the geometry to the degree predicted by VSEPR. Similarly, the octafluoroxenate ion (XeF2−
8) in nitrosonium octafluoroxenate(VI):498 is a square antiprism and not a bicapped trigonal prism (as predicted by VSEPR theory for an AX8E1 molecule), despite having a lone pair. One rationalization is that steric crowding of the ligands allows little or no room for the non-bonding lone pair; another rationalization is the inert pair effect.:214
Square planar transition metal complexes
The Kepert model predicts that AX4 transition metal molecules are tetrahedral in shape, and it cannot explain the formation of square planar complexes.:542 The majority of such complexes exhibit a d8 configuration as for the tetrachloroplatinate (PtCl2−
4) ion. The explanation of the shape of square planar complexes involves electronic effects and requires the use of crystal field theory.:562–4
Complexes with strong d-contribution
Some transition metal complexes with low d electron count have unusual geometries, which can be ascribed to d subshell bonding interaction. Gillespie found that this interaction produces bonding pairs that also occupy the respective antipodal points (ligand opposed) of the sphere. This phenomenon is an electronic effect resulting from the bilobed shape of the underlying sdx hybrid orbitals. The repulsion of these bidirectional bonding pairs leads to a different set of shapes.
|AX6||C3v Trigonal prismatic||W(CH3)6|
The gas phase structures of the triatomic halides of the heavier members of group 2, (i.e., calcium, strontium and barium halides, MX2), are not linear as predicted but are bent, (approximate X–M–X angles: CaF2, 145°; SrF2, 120°; BaF2, 108°; SrCl2, 130°; BaCl2, 115°; BaBr2, 115°; BaI2, 105°). It has been proposed by Gillespie that this is also caused by bonding interaction of the ligands with the d subshell of the metal atom, thus influencing the molecular geometry.
Relativistic effects on the electron orbitals of superheavy elements is predicted to influence the molecular geometry of some compounds. For instance, the 6d5/2 electrons in nihonium play an unexpectedly strong role in bonding, so NhF3 should assume a T-shaped geometry, instead of a trigonal planar geometry like its lighter congener BF3. In contrast, the extra stability of the 7p1/2 electrons in tennessine are predicted to make TsF3 trigonal planar, unlike the T-shaped geometry observed for IF3 and predicted for AtF3; similarly, OgF4 should have a tetrahedral geometry, while XeF4 has a square planar geometry and RnF4 is predicted to have the same.
The VSEPR theory can be extended to molecules with an odd number of electrons by treating the unpaired electron as a "half electron pair" — for example, Gillespie and Nyholm:364–365 suggested that the decrease in the bond angle in the series NO+
2 (180°), NO2 (134°), NO−
2 (115°) indicates that a given set of bonding electron pairs exert a weaker repulsion on a single non-bonding electron than on a pair of non-bonding electrons. In effect, they considered nitrogen dioxide as an AX2E0.5 molecule, with a geometry intermediate between NO+
2 and NO−
2. Similarly, chlorine dioxide (ClO2) is an AX2E1.5 molecule, with a geometry intermediate between ClO+
2 and ClO−
Finally, the methyl radical (CH3) is predicted to be trigonal pyramidal like the methyl anion (CH−
3), but with a larger bond angle (as in the trigonal planar methyl cation (CH+
3)). However, in this case, the VSEPR prediction is not quite true, as CH3 is actually planar, although its distortion to a pyramidal geometry requires very little energy.
- Bent's rule (effect of ligand electronegativity)
- Linear combination of atomic orbitals
- Molecular geometry
- Molecular modelling
- Software for molecular modeling
- Thomson problem
- Valency interaction formula
- Petrucci, R. H.; W. S., Harwood; F. G., Herring (2002). General Chemistry: Principles and Modern Applications (8th ed.). Prentice-Hall. ISBN 978-0-13-014329-7.
- Stoker, H. Stephen (2009). General, Organic, and Biological Chemistry. Cengage Learning. p. 119. ISBN 978-0-547-15281-3.
- Jolly, W. L. (1984). Modern Inorganic Chemistry. McGraw-Hill. pp. 77–90. ISBN 978-0-07-032760-3.
- Gillespie, R. J. (2008). "Fifty years of the VSEPR model". Coord. Chem. Rev. 252 (12–14): 1315–1327. doi:10.1016/j.ccr.2007.07.007.
- Bader, Richard F. W.; Gillespie, Ronald J.; MacDougall, Preston J. (1988). "A physical basis for the VSEPR model of molecular geometry". J. Am. Chem. Soc. 110 (22): 7329–7336. doi:10.1021/ja00230a009.
- Gillespie, R. J. (2004), "Teaching molecular geometry with the VSEPR model", J. Chem. Educ., 81 (3): 298–304, Bibcode:2004JChEd..81..298G, doi:10.1021/ed081p298
- Tsuchida, Ryutarō (1939). "A New Simple Theory of Valency" 新簡易原子價論 [New simple valency theory]. Nippon Kagaku Kaishi (in Japanese). 60 (3): 245–256. doi:10.1246/nikkashi1921.60.245.
- Sidgwick, N. V.; Powell, H. M. (1940). "Bakerian Lecture. Stereochemical Types and Valency Groups". Proc. Roy. Soc. A. 176 (965): 153–180. Bibcode:1940RSPSA.176..153S. doi:10.1098/rspa.1940.0084.
- Gillespie, R. J.; Nyholm, R. S. (1957). "Inorganic stereochemistry". Q. Rev. Chem. Soc. 11 (4): 339. doi:10.1039/QR9571100339.
- Gillespie, R. J. (1970). "The electron-pair repulsion model for molecular geometry". J. Chem. Educ. 47 (1): 18. Bibcode:1970JChEd..47...18G. doi:10.1021/ed047p18.
- Miessler, G. L.; Tarr, D. A. (1999). Inorganic Chemistry (2nd ed.). Prentice-Hall. pp. 54–62. ISBN 978-0-13-841891-5.
- Miessler, G. L.; Tarr, D. A. (1999). Inorganic Chemistry (2nd ed.). Prentice-Hall. p. 55. ISBN 978-0-13-841891-5.
- Petrucci, R. H.; W. S., Harwood; F. G., Herring (2002). General Chemistry: Principles and Modern Applications (8th ed.). Prentice-Hall. pp. 413–414 (Table 11.1). ISBN 978-0-13-014329-7.
- Housecroft, C. E.; Sharpe, A. G. (2005). Inorganic Chemistry (2nd ed.). Pearson. ISBN 978-0-130-39913-7.
- Baran, E. (2000). "Mean amplitudes of vibration of the pentagonal pyramidal XeOF−
5 and IOF2−
5 anions". J. Fluorine Chem. 101: 61–63. doi:10.1016/S0022-1139(99)00194-3.
- Anderson, O. P. (1983). "Book reviews: Inorganic Stereochemistry (by David L. Kepert)" (PDF). Acta Crystallographica B. 39: 527–528. doi:10.1107/S0108768183002864. Retrieved 14 September 2020.
based on a systematic quantitative application of the common ideas regarding electron-pair repulsion
- Brittin, W. E. (1945). "Valence Angle of the Tetrahedral Carbon Atom". J. Chem. Educ. 22 (3): 145. Bibcode:1945JChEd..22..145B. doi:10.1021/ed022p145.
- "Angle Between 2 Legs of a Tetrahedron" Archived 2018-10-03 at the Wayback Machine – Maze5.net
- Wiberg, E.; Holleman, A. F. (2001). Inorganic Chemistry. Academic Press. ISBN 978-0-12-352651-9.
- Power, Philip P. (September 2003). "Silicon, germanium, tin and lead analogues of acetylenes". Chem. Commun. (17): 2091–2101. doi:10.1039/B212224C. PMID 13678155.
- Nagase, Shigeru; Kobayashi, Kaoru; Takagi, Nozomi (6 October 2000). "Triple bonds between heavier Group 14 elements. A theoretical approach". J. Organomet. Chem. 11 (1–2): 264–271. doi:10.1016/S0022-328X(00)00489-7.
- Sekiguchi, Akira; Kinjō, Rei; Ichinohe, Masaaki (September 2004). "A Stable Compound Containing a Silicon–Silicon Triple Bond" (PDF). Science. 305 (5691): 1755–1757. Bibcode:2004Sci...305.1755S. doi:10.1126/science.1102209. PMID 15375262. S2CID 24416825.
- Bellert, D.; Breckenridge, W. H. (2001). "A spectroscopic determination of the bond length of the LiOLi molecule: Strong ionic bonding". J. Chem. Phys. 114 (7): 2871. Bibcode:2001JChPh.114.2871B. doi:10.1063/1.1349424.
- Gillespie, R. J.; Robinson, E. A. (2005). "Models of molecular geometry". Chem. Soc. Rev. 34 (5): 396–407. doi:10.1039/b405359c. PMID 15852152.
- Wells, A. F. (1984). Structural Inorganic Chemistry (5th ed.). Oxford Science Publications. ISBN 978-0-19-855370-0.
- Peterson, W.; Holloway, H.; Coyle, A.; Williams, M. (Sep 1971). "Antiprismatic Coordination about Xenon: the Structure of Nitrosonium Octafluoroxenate(VI)". Science. 173 (4003): 1238–1239. Bibcode:1971Sci...173.1238P. doi:10.1126/science.173.4003.1238. ISSN 0036-8075. PMID 17775218. S2CID 22384146.
- Hanson, Robert M. (1995). Molecular origami: precision scale models from paper. University Science Books. ISBN 978-0-935702-30-9.
- Kaupp, Martin (2001). ""Non-VSEPR" Structures and Bonding in d0 Systems". Angew. Chem. Int. Ed. Engl. 40 (1): 3534–3565. doi:10.1002/1521-3773(20011001)40:19<3534::AID-ANIE3534>3.0.CO;2-#.
- Gillespie, Ronald J.; Noury, Stéphane; Pilmé, Julien; Silvi, Bernard (2004). "An Electron Localization Function Study of the Geometry of d0 Molecules of the Period 4 Metals Ca to Mn". Inorg. Chem. 43 (10): 3248–3256. doi:10.1021/ic0354015. PMID 15132634.
- Landis, C. R.; Cleveland, T.; Firman, T. K. (1995). "Making sense of the shapes of simple metal hydrides". J. Am. Chem. Soc. 117 (6): 1859–1860. doi:10.1021/ja00111a036.
- Landis, C. R.; Cleveland, T.; Firman, T. K. (1996). "Structure of W(CH3)6". Science. 272 (5259): 179–183. doi:10.1126/science.272.5259.179f.
- King, R. Bruce (2000). "Atomic orbitals, symmetry, and coordination polyhedra". Coord. Chem. Rev. 197: 141–168. doi:10.1016/s0010-8545(99)00226-x.
- Haalan, A.; Hammel, A.; Rydpal, K.; Volden, H. V. (1990). "The coordination geometry of gaseous hexamethyltungsten is not octahedral". J. Am. Chem. Soc. 112 (11): 4547–4549. doi:10.1021/ja00167a065.
- Greenwood, Norman N.; Earnshaw, Alan (1997). Chemistry of the Elements (2nd ed.). Butterworth-Heinemann. ISBN 978-0-08-037941-8.
- Seijo, Luis; Barandiarán, Zoila; Huzinaga, Sigeru (1991). "Ab initio model potential study of the equilibrium geometry of alkaline earth dihalides: MX2 (M=Mg, Ca, Sr, Ba; X=F, Cl, Br, I)" (PDF). J. Chem. Phys. 94 (5): 3762. Bibcode:1991JChPh..94.3762S. doi:10.1063/1.459748. hdl:10486/7315.
- Seth, Michael; Schwerdtfeger, Peter; Fægri, Knut (1999). "The chemistry of superheavy elements. III. Theoretical studies on element 113 compounds". Journal of Chemical Physics. 111 (14): 6422–6433. Bibcode:1999JChPh.111.6422S. doi:10.1063/1.480168. S2CID 41854842.
- Bae, Ch.; Han, Y.-K.; Lee, Yo. S. (18 January 2003). "Spin−Orbit and Relativistic Effects on Structures and Stabilities of Group 17 Fluorides EF3 (E = I, At, and Element 117): Relativity Induced Stability for the D3h Structure of (117)F3". The Journal of Physical Chemistry A. 107 (6): 852–858. Bibcode:2003JPCA..107..852B. doi:10.1021/jp026531m.
- Han, Young-Kyu; Lee, Yoon Sup (1999). "Structures of RgFn (Rg = Xe, Rn, and Element 118. n = 2, 4.) Calculated by Two-component Spin-Orbit Methods. A Spin-Orbit Induced Isomer of (118)F4". Journal of Physical Chemistry A. 103 (8): 1104–1108. Bibcode:1999JPCA..103.1104H. doi:10.1021/jp983665k.
- Anslyn, E. V.; Dougherty, D. A. (2006). Modern Physical Organic Chemistry. University Science Books. p. 57. ISBN 978-1891389313.
|The Wikibook A-level Chemistry/OCR (Salters) has a page on the topic of: Molecular geometry and lone pairs|
- VSEPR AR - 3D VSEPR Theory Visualization with Augmented Reality app
- 3D Chem – Chemistry, structures, and 3D Molecules
- IUMSC – Indiana University Molecular Structure Center