<?xml version="1.0" encoding="UTF-8" standalone="yes"?><oembed><version><![CDATA[1.0]]></version><provider_name><![CDATA[Azimuth]]></provider_name><provider_url><![CDATA[https://johncarlosbaez.wordpress.com]]></provider_url><author_name><![CDATA[John Baez]]></author_name><author_url><![CDATA[https://johncarlosbaez.wordpress.com/author/johncarlosbaez/]]></author_url><title><![CDATA[Network Theory (Part&nbsp;14)]]></title><type><![CDATA[link]]></type><html><![CDATA[<p>We&#8217;ve been doing a lot of hard work lately.  Let&#8217;s take a break and think about a fun example from chemistry! </p>
<h4> The ethyl cation </h4>
<p>Suppose you start with a molecule of <a href="http://en.wikipedia.org/wiki/Ethane">ethane</a>, which has 2 carbons and 6 hydrogens arranged like this:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Ethane"><img border="none" src="https://i2.wp.com/math.ucr.edu/home/baez/networks/200px-Ethane-A-3D-balls.png" alt="" /></a>
</div>
<p>Then suppose you remove one hydrogen.  The result is a positively charged ion, or <a href="http://en.wikipedia.org/wiki/Ion">&#8216;cation&#8217;</a>.  When I was a kid, I thought the opposite of a cation should be called a &#8216;dogion&#8217;.  Alas, it&#8217;s <a href="http://en.wikipedia.org/wiki/Ion#Anions_and_cations">not</a>.</p>
<p>This particular cation, formed from removing one hydrogen from an ethane molecule, is called an &#8216;ethyl cation&#8217;.  People used to think it looked like this:</p>
<div align="center"><a href="http://comporgchem.com/blog/?p=60"><img border="none" src="https://i0.wp.com/math.ucr.edu/home/baez/networks/ethyl_cation_classical.jpg" alt="" /></a></div>
<p>They also thought a hydrogen could hop from the carbon with 3 hydrogens attached to it to the carbon with 2.  So, they drew a graph with a vertex for each way the hydrogens could be arranged, and an edge for each hop.  It looks really cool:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Desargues_graph"><img border="none" src="https://i0.wp.com/math.ucr.edu/home/baez/networks/300px-Desargues_graph_2COL.svg.png" alt="" /></a>
</div>
<p>The red vertices come from arrangements where the <i>first</i> carbon has 2 hydrogens attached to it, and the blue vertices come from those where the <i>second</i> carbon has 2 hydrogens attached to it.  So, each edge goes between a red vertex and a blue vertex.  </p>
<p>This graph has 20 vertices, which are arrangements or &#8216;states&#8217; of the ethyl cation.  It has 30 edges, which are hops or &#8216;transitions&#8217;.  Let&#8217;s see why those numbers are right.  </p>
<p>First I need to explain the rules of the game.  The rules say that the 2 carbon atoms are distinguishable: there&#8217;s a &#8216;first&#8217; one and a &#8216;second&#8217; one.  The 5 hydrogen atoms are also distinguishable.  But, all we care about is which carbon atom each hydrogen is bonded to: we don&#8217;t care about further details of its location.  And we require that 2 of the hydrogens are bonded to one carbon, and 3 to the other.</p>
<p>If you&#8217;re a physicist, you may wonder why the rules work this way: after all, at a fundamental level, identical particles <i>aren&#8217;t</i> really distinguishable.  I&#8217;m afraid I can&#8217;t give a fully convincing explanation right now: I&#8217;m just reporting the rules as they were told to me!</p>
<p>Given these rules, there are 2 choices of which carbon has two hydrogens attached to it. Then there are </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B+%5Cbinom%7B5%7D%7B2%7D+%3D+%5Cfrac%7B5+%5Ctimes+4%7D%7B2+%5Ctimes+1%7D+%3D+10%7D+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{ &#92;binom{5}{2} = &#92;frac{5 &#92;times 4}{2 &#92;times 1} = 10} ' title='&#92;displaystyle{ &#92;binom{5}{2} = &#92;frac{5 &#92;times 4}{2 &#92;times 1} = 10} ' class='latex' /> </p>
<p>choices of which two hydrogens are attached to it.  This gives a total of 2 &times; 10 = 20 states.  These are the vertices of our graph: 10 red and 10 blue.</p>
<p>The edges of the graph are transitions between states.  Any hydrogen in the group of 3 can hop over to the group of 2.  There are 3 choices for which hydrogen atom makes the jump.  So, starting from any vertex in the graph there are 3 edges.  This means there are <img src='https://s0.wp.com/latex.php?latex=3+%5Ctimes+20+%2F+2+%3D+30&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='3 &#92;times 20 / 2 = 30' title='3 &#92;times 20 / 2 = 30' class='latex' /> edges.  </p>
<p>Why divide by 2?  Because each edge touches 2 vertices.  We have to avoid double-counting them.</p>
<h4> The Desargues graph </h4>
<p>The idea of using this graph in chemistry goes back to this paper:</p>
<p>&bull; A. T. Balaban, D. Fǎrcaşiu and R. Bǎnicǎ, Graphs of multiple 1,2-shifts in carbonium ions and related systems, <i>Rev. Roum. Chim.</i> <b>11</b> (1966), 1205.</p>
<p>This paper is famous because it was the first to use graphs in chemistry to describe molecular transitions, as opposed to using them as pictures of molecules!</p>
<p>But this particular graph was already famous for other reasons.  It&#8217;s called the <b>Desargues-Levi graph</b>, or <b>Desargues graph</b> for short:</p>
<p>&bull; <a href="http://en.wikipedia.org/wiki/Desargues_graph">Desargues graph</a>, Wikipedia.</p>
<p>Later I&#8217;ll say why it&#8217;s called this.  </p>
<p>There are lots of nice ways to draw the Desargues graph.  For example:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Desargues_graph"><img border="none" src="https://i2.wp.com/math.ucr.edu/home/baez/networks/300px-DesarguesGraph.svg.png" alt="" /></a>
</div>
<p>The reason why we can draw such pretty pictures is that the Desargues graph is very symmetrical.  Clearly any permutation of the 5 hydrogens acts as a symmetry of the graph, and so does any permutation of the 2 carbons.  This gives a symmetry group <img src='https://s0.wp.com/latex.php?latex=S_5+%5Ctimes+S_2&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='S_5 &#92;times S_2' title='S_5 &#92;times S_2' class='latex' />, which has <img src='https://s0.wp.com/latex.php?latex=5%21+%5Ctimes+2%21+%3D+240&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='5! &#92;times 2! = 240' title='5! &#92;times 2! = 240' class='latex' /> elements.  And in fact this turns out to be the full symmetry group of the Desargues graph.</p>
<p>The Desargues graph, its symmetry group, and its applications to chemistry are discussed here:</p>
<p>&bull; Milan Randic, Symmetry properties of graphs of interest in chemistry: II: Desargues-Levi graph, <i><a href="http://onlinelibrary.wiley.com/doi/10.1002/qua.560150610/abstract;jsessionid=0FB39FF9D7D985C1B6454CF4E30D9E92.d03t02">Int. Jour. Quantum Chem.</a></i> <b>15</b> (1997), 663-682.</p>
<h4> The ethyl cation, revisited </h4>
<p>We can try to describe the ethyl cation using probability theory.  If at any moment its state corresponds to some vertex of the Desargues graph, and it hops randomly along edges as time goes by, it will trace out a random walk on the Desargues graph.  This is a nice example of a <a href="http://en.wikipedia.org/wiki/Continuous-time_Markov_process">Markov process</a>!  </p>
<p>We could also try to describe the ethyl cation using quantum mechanics.   Then, instead of having a <i>probability</i> of hopping along an edge, it has an <i>amplitude</i> of doing so.   But as we&#8217;ve seen, a lot of similar math will still apply.  </p>
<p>It should be fun to compare the two approaches.  But I bet you&#8217;re wondering which approach is correct.  This is a somewhat tricky question, at least for me.  The answer would seem to depend on how much the ethyl cation is interacting with its environment&mdash;for example, bouncing off other molecules.  When a system is interacting a lot with its environment, a probabilistic approach seems to be more appropriate.  The relevant buzzword is <a href="http://en.wikipedia.org/wiki/Quantum_decoherence">&#8216;environmentally induced decoherence&#8217;</a>.   </p>
<p>However, there&#8217;s something much more basic I have tell you about.</p>
<p>After the paper by Balaban, Fǎrcaşiu and Bǎnicǎ came out, people gradually realized that the ethyl cation <i>doesn&#8217;t really look like the drawing I showed you!</i>  It&#8217;s what chemists call &#8216;nonclassical&#8217; ion.  What they mean is this: its actual structure is not what you get by taking the traditional ball-and-stick model of an ethane molecule and ripping off a hydrogen.  The ethyl cation really looks like this:</p>
<div align="center"><a href="http://comporgchem.com/blog/?p=60"><img border="none" src="https://i0.wp.com/math.ucr.edu/home/baez/networks/ethyl_cation_nonclassical.jpg" alt="" /></a></div>
<p>For more details, and pictures that you can actually rotate, see:</p>
<p>&bull; Stephen Bacharach, Ethyl cation, <a href="http://comporgchem.com/blog/?p=60">Computational Organic Chemistry</a>.</p>
<p>So, if we stubbornly insist on applying the Desargues graph to realistic chemistry, we need to find some other molecule to apply it to.</p>
<h4> Trigonal bipyramidal molecules </h4>
<p>Luckily, there are lots of options!  They&#8217;re called <a href="http://en.wikipedia.org/wiki/Trigonal_bipyramidal_molecular_geometry">trigonal bipyramidal molecules</a>.  They look like this:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Trigonal_bipyramidal_molecular_geometry"><img border="none" src="https://i0.wp.com/math.ucr.edu/home/baez/networks/Trigonal-bipyramidal-3D-balls.png/200px-Trigonal-bipyramidal-3D-balls.png" alt="" /></a>
</div>
<p>The 5 balls on the outside are called &#8216;ligands&#8217;: they could be atoms or bunches of atoms.   In chemistry, &#8216;<a href="http://en.wikipedia.org/wiki/Ligand">ligand</a>&#8216; just means something that&#8217;s stuck onto a central thing.   For example, in <a href="http://en.wikipedia.org/wiki/Phosphorus_pentachloride">phosphorus pentachloride</a> the ligands are chlorine atoms, all attached to a central phosphorus atom:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Phosphorus_pentachloride"><img border="none" src="https://i1.wp.com/math.ucr.edu/home/baez/networks/200px-Phosphorus-pentachloride-3D-balls.png" alt="" /></a>
</div>
<p>It&#8217;s a colorless solid, but as you might expect, it&#8217;s pretty <a href="http://www.cbwinfo.com/Chemical/Precursors/p36.html">nasty stuff</a>: it&#8217;s not flammable, but it reacts with water or heat to produce toxic chemicals like hydrogen chloride.</p>
<p>Another example is <a href="http://en.wikipedia.org/wiki/Iron_pentacarbonyl">iron pentacarbonyl</a>, where 5 carbon-oxygen ligands are attached to a central iron atom:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Iron_pentacarbonyl"><img border="none" src="https://i2.wp.com/math.ucr.edu/home/baez/networks/200px-Iron-pentacarbonyl-from-xtal-3D-balls.png" alt="" /></a>
</div>
<p>You can make this stuff by letting powdered iron react with carbon monoxide.  It&#8217;s a straw-colored liquid with a pungent smell!</p>
<p>Whenever you&#8217;ve got a molecule of this shape, the ligands come in two kinds.  There are the 2 &#8216;axial&#8217; ones, and the 3 &#8216;equatorial&#8217; ones:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Trigonal_bipyramidal_molecular_geometry"><img border="none" src="https://i1.wp.com/math.ucr.edu/home/baez/networks/200px-Trigonal-bipyramidal-3D-balls-ax-eq.png" alt="" /></a>
</div>
<p>And the molecule has 20 states&#8230; but only if count the states a certain way.   We have to treat all 5 ligands as distinguishable, but think of two arrangements of them as the same if we can rotate one to get the other.   The trigonal bipyramid has a rotational symmetry group with 6 elements.  So, there are 5! / 6 = 20 states.  </p>
<p>The transitions between states are devilishly tricky.  They&#8217;re called <b>pseudorotations</b>, and they look like this:</p>
<div align="center"><a href="http://en.wikipedia.org/wiki/Berry_mechanism"><img border="none" width="450" src="https://i1.wp.com/math.ucr.edu/home/baez/networks/700px-Iron-pentacarbonyl-Berry-mechanism.png" alt="" /></a></div>
<p>If you look very carefully, you&#8217;ll see what&#8217;s going on.  First the 2 axial ligands move towards each other to become equatorial.<br />
Now the equatorial ones are no longer in the horizontal plane: they&#8217;re in the plane facing us!   Then 2 of the 3 equatorial ones swing out to become axial.    This fancy dance is called the <a href="http://en.wikipedia.org/wiki/Berry_mechanism">Berry pseudorotation mechanism</a>.</p>
<p>To get from one state to another this way, we have to pick 2 of the 3 equatorial ligands to swing out and become axial.  There are 3 choices here.   So, if we draw a graph with states as vertices and transitions as edges, it will have 20 vertices and 20 &times; 3 / 2 = 30 edges.  That sounds suspiciously like the Desargues graph!</p>
<p><b>Puzzle 1</b>.  Show that the graph with states of a trigonal bipyramidal molecule as vertices and pseudorotations as edges is indeed the Desargues graph.</p>
<p>I think this fact was first noticed here:</p>
<p>&bull; Paul C. Lauterbur and Fausto Ramirez, Pseudorotation in trigonal-bipyramidal molecules, <a href="http://pubs.acs.org/doi/abs/10.1021/ja01026a029"><i>J. Am. Chem. Soc.</i></a> <b>90</b> (1968), 6722–6726.</p>
<p>Okay, enough for now!  Next time I&#8217;ll say more about the Markov process or quantum process corresponding to a random walk on the Desargues graph.  But since the Berry pseudorotation mechanism is so hard to visualize, I&#8217;ll <i>pretend</i> that the ethyl cation looks like this:</p>
<div align="center"><a href="http://comporgchem.com/blog/?p=60"><img border="none" src="https://i0.wp.com/math.ucr.edu/home/baez/networks/ethyl_cation_classical.jpg" alt="" /></a></div>
<p>and I&#8217;ll use this picture to help us think about the Desargues graph.</p>
<p>That&#8217;s okay: everything we&#8217;ll figure out can easily be translated to apply to the real-world situation of a trigonal bipyramidal molecule.  The virtue of math is that when two situations are &#8216;mathematically the same&#8217;, or &#8216;isomorphic&#8217;, we can talk about either one, and the results automatically apply to the other.  This is true even if the one we talk about doesn&#8217;t actually exist in the real world!</p>
]]></html><thumbnail_url><![CDATA[https://i2.wp.com/math.ucr.edu/home/baez/networks/200px-Ethane-A-3D-balls.png?fit=440%2C330]]></thumbnail_url><thumbnail_height><![CDATA[154]]></thumbnail_height><thumbnail_width><![CDATA[200]]></thumbnail_width></oembed>