{"id":23855,"date":"2021-06-02T12:18:00","date_gmt":"2021-06-02T11:18:00","guid":{"rendered":"https:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=23855"},"modified":"2021-06-02T12:18:00","modified_gmt":"2021-06-02T11:18:00","slug":"a-closer-look-at-that-fourth-bond-in-c2","status":"publish","type":"post","link":"https:\/\/www.rzepa.net\/blog\/?p=23855","title":{"rendered":"A closer look at that fourth bond in C2."},"content":{"rendered":"<div class=\"kcite-section\" kcite-section-id=\"23855\">\n<p>From the last few posts here, you might have noticed much discussion about how the element carbon might sustain a quadruple bond. The <a href=\"https:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=3065\" rel=\"noopener\" target=\"_blank\">original post on this topic<\/a> from some years ago showed the molecular orbitals of the species CN<sup>+<\/sup>, which included two bonding &pi;-types and a low lying nodeless bonding &sigma;-orbital, all with double occupancies and adding up to a triple bond. Discussing now C<sub>2<\/sub> itself, there are two remaining orbitals for consideration which we will for the purpose here call the highest occupied &sigma;-MO or HO&sigma;MO (&Sigma;<sub>u<\/sub>) and the lowest unoccupied &sigma;-MO or LU&sigma;MO (&Sigma;<sub>g<\/sub>) and which are more mysterious.<\/p>\n<p>The HO&sigma;MO itself has one node (the lowest unoccupied or LU&sigma;MO has a further second node) bisecting the centre of the C-C bond, which makes it anti-bonding. This is emphasised by squaring the orbital (below), which shows a clear void of electron density in the C-C region. For this reason, many text books illustrating the main group diatomic molecules represent C<sub>2<\/sub> with two bonds: 3-1 = <b>2<\/b>.<\/p>\n<table border=\"1\">\n<tbody>\n<tr>\n<th colspan=\"2\">CASSCF(8,6) &#8220;HOMO&#8221;\/&#8221;LUMO&#8221; orbitals and densities of C<sub>2<\/sub><\/th>\n<\/tr>\n<tr>\n<th>(HO&sigma;MO)@0.02au<\/th>\n<th>(LU&sigma;MO)@0.02au<\/th>\n<\/tr>\n<tr>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([200,200],'load wp-content\/uploads\/2021\/06\/C2_mo6.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_mo6.jvxl;zoom 100;spin 3;','c1');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_mo6.png?w=200&#038;ssl=1\" alt=\"\"  \/><\/td>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([200,200],'load wp-content\/uploads\/2021\/06\/C2_mo7.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_mo7.jvxl;zoom 100;spin 3;','c2');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_mo7.png?w=200&#038;ssl=1\" alt=\"\"  \/><\/td>\n<\/tr>\n<tr>\n<th>(HO&sigma;MO<sup>2<\/sup>)@0.0004au<\/th>\n<th>(LU&sigma;MO<sup>2<\/sup>)@0.0004au<\/th>\n<\/tr>\n<tr>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([200,200],'load wp-content\/uploads\/2021\/06\/C2_6-sqrd.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_6-sqrd.jvxl;zoom 100;spin 3;','c3');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_6-sqrd.png?w=200&#038;ssl=1\" alt=\"\"  \/><\/td>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([200,200],'load wp-content\/uploads\/2021\/06\/C2_7-sqrd.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_7-sqrd.jvxl;zoom 100;spin 3;','c4');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_7-sqrd.png?w=200&#038;ssl=1\" alt=\"\"  \/><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Now to a MCSCF(8,6)\/Def2-SVPD calculation (FAIR DOI: <a href=\"https:\/\/doi.org\/10.14469\/hpc\/8307\" rel=\"noopener\" target=\"_blank\">10.14469\/hpc\/8307<\/a>) which means a multi-configuration calculation. Regular <i>e.g.<\/i> DFT methods assume only a single electronic configuration in which one set of doubly-occupied orbitals is variationally optimized. Thankfully, for the vast majority of molecules, this is actually a very good approximation. However for some species, and C<sub>2<\/sub> is one such, this is no longer true. A CASSCF(8,6) calculation uses the 105 different electronic configurations generated by using eight electrons in an active space of six orbitals and variationally optimises them all for a self-consistent-field. The orbitals corresponding to the erstwhile single-configurational HO&sigma;MO and LU&sigma;MO are the ones shown above. Their squares are shown underneath, the latter as noted above relating to the electron density distribution in the molecule. <\/p>\n<p> The MCSCF calculation for C<sub>2<\/sub> shows that primarily <b>two<\/b> different electronic configurations contribute significantly to the total wavefunction, the one with two electrons in the original HO&sigma;MO (now of course a misnomer) having a weight of 1.573e and the configuration with two electrons promoted to the now similarly misnamed LU&sigma;MO orbital by virtue of having a weight of 0.427e.<sup>&Dagger;<\/sup> This shows that this final 2e really must be described by two electronic configurations rather than one (and which reminds that the terms HO&sigma;MO and LU&sigma;MO really only apply to single-configuration methods). What difference does that make to the picture? The scaled linear combination of the two orbitals deriving from the two dominant electronic configurations for C<sub>2<\/sub> below shows that each now has an &#8220;extrusion&#8221; of the original orbital creeping along the C-C axis.<\/p>\n<table border=\"1\">\n<tbody>\n<tr>\n<th>(1.573*HO&sigma;MO + 0.427*LU&sigma;MO)@0.02au<\/th>\n<th>(1.573*HO&sigma;MO &#8211; 0.427*LU&sigma;MO)@0.02au<\/th>\n<\/tr>\n<tr>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([200,200],'load wp-content\/uploads\/2021\/06\/C2_add.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_add.jvxl;zoom 100;spin 3;','c5');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_add.png?w=200&#038;ssl=1\" alt=\"\"  \/><\/td>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([200,200],'load wp-content\/uploads\/2021\/06\/C2_sub.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_sub.jvxl;zoom 100;spin 3;','c6');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_sub.png?w=200&#038;ssl=1\" alt=\"\"  \/><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Squaring and weighted adding shows us what is happening to the electron density. That void in the density along the C-C axis apparent in the HO&sigma;MO above has now been nicely filled with density from those extrusions deriving from the partially occupied &#8220;LU&sigma;MO&#8221;. As a result, the node in the density along the bond has now vanished. By elevating 0.427 of the electrons from the original anti-bonding HO&sigma;MO into the complementary LU&sigma;MO, a new weakly-bonding &#8220;baby&#8221; orbital with ~two-electron occupancy replaces the original antibonding HO&sigma;MO. There is however relatively little additional density placed into the C-C region because of only 0.427e transfer into the &#8220;LU&sigma;MO&#8221;. The weak bonding character also matches the &#8220;bond dissociation energy&#8221; of this fourth bond of ~17 kcal\/mol as inferred by experimental measurement of the energies of the two reactions HC&equiv;CH &rarr; HC&equiv;C<sup>\u2022<\/sup> + H<sup>\u2022<\/sup>; HC&equiv;C<sup>\u2022<\/sup> &rarr; CC + H<sup>\u2022<\/sup>.\n<\/p>\n<table border=\"1\">\n<tbody>\n<tr>\n<th>(1.573*(HO&sigma;MO)<sup>2<\/sup> + 0.427*(HO&sigma;MO)<sup>2<\/sup>)@0.0004au<\/th>\n<\/tr>\n<tr>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([440,350],'load wp-content\/uploads\/2021\/06\/C2_1.564620.43072.cub_.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_1.564620.43072.cub_.jvxl;zoom 70;spin 3;','c7');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_1.564620.43072.png?w=300&#038;ssl=1\" alt=\"\"  \/><\/td>\n<\/tr>\n<tr>\n<th>The (&Sigma;<sub>g<\/sub>)<sup>2<\/sup>@0.0004au conventional &sigma;-bond for comparison<\/th>\n<\/tr>\n<tr>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([440,350],'load wp-content\/uploads\/2021\/06\/C2_3-sqrd.xyz;isosurface wp-content\/uploads\/2021\/06\/C2_3-sqrd.jvxl;zoom 70;spin 3;','c8');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/C2_3-sqrd.png?w=190&#038;ssl=1\" alt=\"\"  \/><\/td>\n<\/tr>\n<tr>\n<th><small>(&Sigma;<sub>g<\/sub>)<sup>2<\/sup> &#8211; (1.573*(HO&sigma;MO)<sup>2<\/sup> + 0.427*(HO&sigma;MO)<sup>2<\/sup>)@0.0004au<\/small><\/th>\n<\/tr>\n<tr>\n<td><img data-recalc-dims=\"1\" decoding=\"async\" class=\"aligncenter size-medium wp-image-23596\" onclick=\"jmolApplet([440,350],'load wp-content\/uploads\/2021\/06\/CC-1212_3-67.xyz;isosurface wp-content\/uploads\/2021\/06\/CC-1212_3-67.jvxl;zoom 70;spin 3;','c9');\" src=\"https:\/\/i0.wp.com\/www.rzepa.net\/blog\/wp-content\/uploads\/2021\/06\/CC-1212_3-67.png?w=190&#038;ssl=1\" alt=\"\"  \/><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>So by combining the appropriate occupancies of the HO&sigma;MO and LU&sigma;MO in a multi-configurational approach to C<sub>2<\/sub>, a new weak bond emerges, which when added to the three existing bonds referred to above gives a representation of four rather than two bonds for this molecule.<\/p>\n<hr \/>\n<p>The DOI of this post is <a href=\"https:\/\/doi.org\/gf9s\">https:\/\/doi.org\/gf9s<\/a><\/p>\n<hr \/>\n<p><sup>&Dagger;<\/sup>The weights for a CASSCF(12,12) calculation with 427350 configurations are 1.568 and 0.426. Conversely, a  CASSCF(8,5) calculation on 15 configurations yields 1.533 and 0.467. The optimised geometries also show an interesting trend. Thus 8,5 = 1.198&Aring;, 8,6 = 1.230  and 12,12 = 1.262&Aring;. As the active space decreases, so the weight of the configuration with a populated &Sigma;<sub>g<\/sub> orbital increases by &#8220;concentration&#8221; into this orbital and hence the C-C bond length also decreases as the amount of density injected into the C-C region increases.<\/p>\n<!-- kcite active, but no citations found -->\n<\/div> <!-- kcite-section 23855 -->","protected":false},"excerpt":{"rendered":"<p>From the last few posts here, you might have noticed much discussion about how the element carbon might sustain a quadruple bond. The original post on this topic from some years ago showed the molecular orbitals of the species CN+, which included two bonding &pi;-types and a low lying nodeless bonding &sigma;-orbital, all with double [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":true,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[6],"tags":[1526],"class_list":["post-23855","post","type-post","status-publish","format-standard","hentry","category-interesting-chemistry","tag-interesting-chemistry"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/p1gPyz-6cL","jetpack_likes_enabled":true,"_links":{"self":[{"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=\/wp\/v2\/posts\/23855","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=23855"}],"version-history":[{"count":0,"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=\/wp\/v2\/posts\/23855\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=23855"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=23855"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.rzepa.net\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=23855"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}