{"id":28137,"date":"2025-03-17T19:57:21","date_gmt":"2025-03-17T19:57:21","guid":{"rendered":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137"},"modified":"2025-07-22T09:03:53","modified_gmt":"2025-07-22T08:03:53","slug":"crystallography-meets-dft-quantum-modelling","status":"publish","type":"post","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137","title":{"rendered":"Crystallography meets DFT Quantum modelling."},"content":{"rendered":"<div class=\"kcite-section\" kcite-section-id=\"28137\">\n<p>X-ray crystallography is the technique of using the diffraction of x-rays by the electrons in a molecule to determine the positions of all the atoms in that molecule. Quantum theory teaches us that the electrons are to be found in shells around the atomic nuclei. There are two broad types, the outermost shell (also called the valence shell) and all the inner or core shells. The density of the core electrons is much higher (more compact) than the more diffuse valence shell for all but the hydrogen atom, which only has valence electrons. How does this relate to x-ray diffraction by electrons? Well, core electrons, because of their relative compactness, diffract X-rays more strongly than the valence electrons. This compactness of the core also means that its electron density distribution can be well (but not exactly) approximated by a sphere, with the nucleus at the centre of that sphere. And from this it follows that the density for each atom can be treated independently, the so-called IAM or independent atom model. For example all the carbon atoms in a molecule are approximated as having the same value for the electron density of their core shell. But the IAM approximation is much less good for hydrogen atoms, especially when they are attached to very polar atoms (Li, O, F, etc) and even atoms such as carbon or oxygen have noticeable deviations as illustrated in \u00a0figure 1 below. <span id=\"cite_ITEM-28137-0\" name=\"citation\"><a href=\"#ITEM-28137-0\">[1]<\/a><\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-full wp-image-28356\" src=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/d0sc05526c-f1.jpg\" alt=\"\" width=\"530\" height=\"279\" \/><br \/>\n<strong>Figure 1<\/strong> from <span id=\"cite_ITEM-28137-0\" name=\"citation\"><a href=\"#ITEM-28137-0\">[1]<\/a><\/span> with caption: <i>Deformation Hirshfeld densities for the carbon (left) and oxygen (right) atoms in the carboxylate group of Gly-l-Ala, i.e. difference between the spherical atomic electron density used in the IAM and the non-spherical Hirshfeld atom density used in Hirshfeld atom refinement=HAR (IAM minus HAR). Red = negative, blue = positive. Isovalue = 0.17 e\u00c5<sup>\u22123<\/sup><\/i>.<sup>\u2021<\/sup><\/p>\n<p>X-ray crystallography is all about matching the electron density map of a model structure with the electron density map derived from the diffraction data. In \u201cconventional\u201d X-ray crystallography &#8211; i.e. that used by most crystallographers &#8211; \u00a0the electron density map of the model is calculated using the IAM approach, where no consideration is given to any distortion of the electron density distribution caused by things like bonds \u2013 each atom is treated independently (hence the name). This method especially struggles with hydrogens and hence the inferred position of the hydrogen nucleus at the centre of an assumed spherical distribution is often difficult to obtain accurately. Enter quantum crystallography, whereby a model of the electron density distribution in a molecule can be calculated by solving the Schrodinger equation, nowadays to a very reasonable approximation in a reasonable time (minutes) using so-called density functional theory, or DFT. The resulting electron density map for the model structure might be expected to more closely match reality than the IAM approach. Most obviously affected by this change is the handling of hydrogen atoms. If one considers a C\u2013H bond from an sp<sup>3<\/sup> carbon atom, using an IAM approach the hydrogen atom (i.e. its nucleus or proton) would be placed at the centre of maximum electron density, in the full knowledge that this is not actually where the hydrogen atom nucleus itself is. The direction of the C\u2013H vector would be correct, but the distance would be too short. In the quantum crystallography approach, the positions of e.g. hydrogen atom nuclei are not exactly coincident with the electron density maxima, amounting in effect to non-spherical atoms, thus avoiding the systematic errors seen in the IAM approach. Smaller, but possibly still significant such errors might be expected for e.g. the 2nd row elements and beyond.<\/p>\n<p>Getting reliable hydrogen atom positions has previously required a neutron diffraction study, which is difficult, expensive and time consuming. So the idea of using the non-spherical DFT densities rather than the spherical IAM approach to build a model using X-ray diffraction data is very appealing. But does it work? To test this, we decided to go back to some previously published structures that were handled using the IAM approach, and re-refining them using quantum crystallography. We do not have the corresponding neutron studies to check the answers against, but we can still see how well the structures themselves refine and what new problems this approach might throw up.<\/p>\n<p><strong>Method<\/strong><\/p>\n<p>The original published structures<span id=\"cite_ITEM-28137-1\" name=\"citation\"><a href=\"#ITEM-28137-1\">[2]<\/a><\/span> were refined with SHELX-2014<span id=\"cite_ITEM-28137-2\" name=\"citation\"><a href=\"#ITEM-28137-2\">[3]<\/a><\/span> which uses an independent atom model (IAM) approach. The results reported here employed <a href=\"https:\/\/www.iucr.org\/resources\/other-directories\/software\/nosphera2\">NoSpherA2<\/a><span id=\"cite_ITEM-28137-0\" name=\"citation\"><a href=\"#ITEM-28137-0\">[1]<\/a><\/span>, <span id=\"cite_ITEM-28137-3\" name=\"citation\"><a href=\"#ITEM-28137-3\">[4]<\/a><\/span> using <a href=\"https:\/\/journals.iucr.org\/m\/issues\/2014\/05\/00\/fc5002\/fc5002.pdf\">Hirshfeld atom refinement<\/a><span id=\"cite_ITEM-28137-4\" name=\"citation\"><a href=\"#ITEM-28137-4\">[5]<\/a><\/span> and selecting Def2-SVP as the (all-atom) basis set<sup>\u2020<\/sup> and \u03c9B97X-V as the DFT method (the results seem relatively insensitive to either), implemented in the ORCA program.<span id=\"cite_ITEM-28137-5\" name=\"citation\"><a href=\"#ITEM-28137-5\">[6]<\/a><\/span> For the first attempts no changes were made to the structures beyond the anisotropic refinement of the now unconstrained hydrogen atoms. For four of the structures a number of the hydrogen atoms went non-positive definite (i.e. one of the radii of a thermal ellipsoid refined to a negative length), which is physically nonsensical and would be a significant barrier to publication. (we don&#8217;t quite want to say \u201cunpublishable\u201d as there are almost always exceptions, but a non-positive definite thermal parameter is pretty close to being unacceptable.) For these cases, a second version was created (V2) where all of the hydrogen atoms were refined isotropically but with the distances and thermal parameters still allowed to refine. For AB1709 (18b), this still had the isotropic thermal parameter of one of the hydrogen atoms (H11) go non-positive definite, so for that one hydrogen atom the free isotropic thermal parameter was replaced with a riding one.<\/p>\n<p><strong>The results<\/strong><\/p>\n<p>We chose a set of seven structures published in 2017<span id=\"cite_ITEM-28137-6\" name=\"citation\"><a href=\"#ITEM-28137-6\">[7]<\/a><\/span> and refined as noted above using conventional methods. These seven also comprise one of the very first sets of crystal structures for which full diffraction data were made available,<span id=\"cite_ITEM-28137-1\" name=\"citation\"><a href=\"#ITEM-28137-1\">[2]<\/a><\/span> rather than just the refined structure in the form of a CIF file. The \u00a0new results have also been deposited<span id=\"cite_ITEM-28137-7\" name=\"citation\"><a href=\"#ITEM-28137-7\">[8]<\/a><\/span> to augment the record for these compounds. Spreadsheets corresponding to the images below can be obtained by clicking on the image.<\/p>\n<p><a href=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/IAMvsQM1.xlsx\"><img decoding=\"async\" src=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/Results1.jpg\" alt=\"\" width=\"850\" \/><\/a><\/p>\n<p><a href=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/IAMvsQM2.xlsx\"><img decoding=\"async\" src=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/Results2.jpg\" alt=\"\" width=\"850\" \/><\/a><\/p>\n<p><a href=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/IAMvsQM3.xlsx\"><img decoding=\"async\" src=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/Results3.jpg\" alt=\"\" width=\"850\" \/><\/a><\/p>\n<ol>\n<li>All seven structures saw a reduction in the final R-factor.<span id=\"cite_ITEM-28137-7\" name=\"citation\"><a href=\"#ITEM-28137-7\">[8]<\/a><\/span> However, all of the structures also saw a significant increase in the number of parameters (as the hydrogen atoms went from using zero parameters each in a fully riding model to nine parameters each in a fully free anisotropic model). However, all the QM refinements passed the <strong>Hamilton test<\/strong>, suggesting that the reduced R-factors <strong>do indeed reflect a better model<\/strong>, rather than just being a consequence of the significantly increased number of parameters.<\/li>\n<li>All four of the structures containing bromine atoms had a number of the hydrogen atoms go non-positive definite when refined anisotropically. It is not clear exactly why this happened \u2013 there does not appear to be any correlation with data quality or intensity (as crudely measured by R(int) and mean I\/\u03c3 respectively), and though the redundancy for these structures is fairly low (between 1.5 and 1.7), those for the non-bromine structures are not much better (1.5, 2.3 and 4.9). These data sets were the result of experiments designed to collect 98.5% of the symmetry unique data with no consideration for redundancy at all. However comparison of the initial and secondary versions of the refinements of these four structure does show that the substantial majority of the observed R-factor decrease can be achieved without using anisotropic hydrogen atoms.<\/li>\n<li>As regards the precision of the structures, using one C(sp2)\u2013C(sp3) bond as a proxy (the C7\u2013C8 bond) we can see that the estimated standard deviation is either the same or only slightly lower in all seven structures, suggesting that getting lower e.s.d.s would not be a motivating factor for using quantum crystallography.<\/li>\n<li>One of the more unexpected results was the variation in F(000). In X-ray crystallography (deliberate emphasis on X-ray, as neutron diffraction is different) F(000) is supposed to be the total number of electrons present in the unit cell, and is used as an overall scale factor for the electron density map. It is very much <b>not<\/b> supposed to be variable, and any discrepancy would indicate an error in the calculated or reported formula and should be corrected. We do not understand why the <strong>QM refinements give a different answer than the IAM ones<\/strong> (some up and some down \u2014 normalised to a per molecule basis the range is \u20131.1 to +2.2), though it seems likely to be associated with cut-offs (boundaries) in measuring the \u201csmeared out\u201d electron density in the QM models, The IAM models all give the expected \u201ccorrect\u201d values.<\/li>\n<li>Based on the checkCIF reports for the QM structures, if quantum crystallography catches on in a big way, then checkCIF will probably need to be updated, there now being a number of high level alerts for long X\u2013H bonds.<\/li>\n<li>One of the major areas of uncertainty with quantum crystallography is what\/how much data needs to be collected. Symmetry unique data to 0.84 \u00c5 seems insufficient, but what would be sufficient \u2014 full sphere, redundancy, higher resolution? Would the final results be worth the extra time investment? None of the above aspects are clear at this stage, but it will be interesting to see how the technique develops.<\/li>\n<\/ol>\n<p>These seven crystal structures also occupy an interesting position for posterity. Data for them has been made available spanning eight years which illustrates two significantly different refinement methods being used during this period, as well as having access to \u00a0the original complete diffraction image data to allow any completely new analysis to be made in the future. \u00a0Who knows, maybe in eight years time an even better method may become available for comparison with the results reported here.<\/p>\n<hr \/>\n<p><sup>\u2021<\/sup>To put this into context, 0.17 eA<sup>-3<\/sup> would generally be regarded as a pretty low level background noise, similar to the value of the maximum residual electron density a crystallographer might be happy with. <sup>\u2020<\/sup>The structure which showed the smallest change in R factor on using quantum crystallography, <i>i.e.<\/i> AB1608b, was re-run with the triple-\u03b6 Def2-TZVPP basis set. This did give lower R factors but by very little (3.38% to 3.36% aniso with npd; 3.39 to 3.38 iso).<\/p>\n<h2>References<\/h2>\n    <ol class=\"kcite-bibliography csl-bib-body\"><li id=\"ITEM-28137-0\">F. Kleemiss, O.V. Dolomanov, M. Bodensteiner, N. Peyerimhoff, L. Midgley, L.J. Bourhis, A. Genoni, L.A. Malaspina, D. Jayatilaka, J.L. Spencer, F. White, B. Grundk\u00f6tter-Stock, S. Steinhauer, D. Lentz, H. Puschmann, and S. Grabowsky, \"Accurate crystal structures and chemical properties from NoSpherA2\", <i>Chemical Science<\/i>, vol. 12, pp. 1675-1692, 2021. <a href=\"https:\/\/doi.org\/10.1039\/d0sc05526c\">https:\/\/doi.org\/10.1039\/d0sc05526c<\/a>\n\n<\/li>\n<li id=\"ITEM-28137-1\">J. Almond-Thynne, \"Crystal structure data for Synthesis and Reactions of Benzannulated Spiroaminals; Tetrahydrospirobiquinolines\", 2017. <a href=\"https:\/\/doi.org\/10.14469\/hpc\/2297\">https:\/\/doi.org\/10.14469\/hpc\/2297<\/a>\n\n<\/li>\n<li id=\"ITEM-28137-2\">G.M. Sheldrick, \"Crystal structure refinement with&lt;i&gt;SHELXL&lt;\/i&gt;\", <i>Acta Crystallographica Section C Structural Chemistry<\/i>, vol. 71, pp. 3-8, 2015. <a href=\"https:\/\/doi.org\/10.1107\/s2053229614024218\">https:\/\/doi.org\/10.1107\/s2053229614024218<\/a>\n\n<\/li>\n<li id=\"ITEM-28137-3\">O.V. Dolomanov, L.J. Bourhis, R.J. Gildea, J.A.K. Howard, and H. Puschmann, \"&lt;i&gt;OLEX2&lt;\/i&gt;: a complete structure solution, refinement and analysis program\", <i>Journal of Applied Crystallography<\/i>, vol. 42, pp. 339-341, 2009. <a href=\"https:\/\/doi.org\/10.1107\/s0021889808042726\">https:\/\/doi.org\/10.1107\/s0021889808042726<\/a>\n\n<\/li>\n<li id=\"ITEM-28137-4\">S.C. Capelli, H. B\u00fcrgi, B. Dittrich, S. Grabowsky, and D. Jayatilaka, \"Hirshfeld atom refinement\", <i>IUCrJ<\/i>, vol. 1, pp. 361-379, 2014. <a href=\"https:\/\/doi.org\/10.1107\/s2052252514014845\">https:\/\/doi.org\/10.1107\/s2052252514014845<\/a>\n\n<\/li>\n<li id=\"ITEM-28137-5\">F. Neese, \"The ORCA program system\", <i>WIREs Computational Molecular Science<\/i>, vol. 2, pp. 73-78, 2011. <a href=\"https:\/\/doi.org\/10.1002\/wcms.81\">https:\/\/doi.org\/10.1002\/wcms.81<\/a>\n\n<\/li>\n<li id=\"ITEM-28137-6\">J. Almond-Thynne, A.J.P. White, A. Polyzos, H.S. Rzepa, P.J. Parsons, and A.G.M. Barrett, \"Synthesis and Reactions of Benzannulated Spiroaminals: Tetrahydrospirobiquinolines\", <i>ACS Omega<\/i>, vol. 2, pp. 3241-3249, 2017. <a href=\"https:\/\/doi.org\/10.1021\/acsomega.7b00482\">https:\/\/doi.org\/10.1021\/acsomega.7b00482<\/a>\n\n<\/li>\n<li id=\"ITEM-28137-7\">H. Rzepa, \"Crystallography meets DFT Quantum modelling\", 2025. <a href=\"https:\/\/doi.org\/10.14469\/hpc\/15030\">https:\/\/doi.org\/10.14469\/hpc\/15030<\/a>\n\n<\/li>\n<\/ol>\n\n<\/div> <!-- kcite-section 28137 -->","protected":false},"excerpt":{"rendered":"<p>X-ray crystallography is the technique of using the diffraction of x-rays by the electrons in a molecule to determine the positions of all the atoms in that molecule. Quantum theory teaches us that the electrons are to be found in shells around the atomic nuclei. There are two broad types, the outermost shell (also called [&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,"activitypub_content_warning":"","activitypub_content_visibility":"","activitypub_max_image_attachments":5,"activitypub_interaction_policy_quote":"anyone","activitypub_status":"","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":[2],"tags":[],"ppma_author":[2661,2662],"class_list":["post-28137","post","type-post","status-publish","format-standard","hentry","category-chemical-it"],"yoast_head":"<!-- This site is optimized with the Yoast SEO plugin v27.5 - https:\/\/yoast.com\/product\/yoast-seo-wordpress\/ -->\n<title>Crystallography meets DFT Quantum modelling. - Henry Rzepa&#039;s Blog<\/title>\n<meta name=\"robots\" content=\"index, follow, max-snippet:-1, max-image-preview:large, max-video-preview:-1\" \/>\n<link rel=\"canonical\" href=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137\" \/>\n<meta property=\"og:locale\" content=\"en_GB\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Crystallography meets DFT Quantum modelling. - Henry Rzepa&#039;s Blog\" \/>\n<meta property=\"og:description\" content=\"X-ray crystallography is the technique of using the diffraction of x-rays by the electrons in a molecule to determine the positions of all the atoms in that molecule. Quantum theory teaches us that the electrons are to be found in shells around the atomic nuclei. There are two broad types, the outermost shell (also called [&hellip;]\" \/>\n<meta property=\"og:url\" content=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137\" \/>\n<meta property=\"og:site_name\" content=\"Henry Rzepa&#039;s Blog\" \/>\n<meta property=\"article:published_time\" content=\"2025-03-17T19:57:21+00:00\" \/>\n<meta property=\"article:modified_time\" content=\"2025-07-22T08:03:53+00:00\" \/>\n<meta name=\"author\" content=\"Henry Rzepa, Andrew White\" \/>\n<meta name=\"twitter:card\" content=\"summary_large_image\" \/>\n<meta name=\"twitter:label1\" content=\"Written by\" \/>\n\t<meta name=\"twitter:data1\" content=\"Henry Rzepa\" \/>\n\t<meta name=\"twitter:label2\" content=\"Estimated reading time\" \/>\n\t<meta name=\"twitter:data2\" content=\"9 minutes\" \/>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Crystallography meets DFT Quantum modelling. - Henry Rzepa&#039;s Blog","robots":{"index":"index","follow":"follow","max-snippet":"max-snippet:-1","max-image-preview":"max-image-preview:large","max-video-preview":"max-video-preview:-1"},"canonical":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137","og_locale":"en_GB","og_type":"article","og_title":"Crystallography meets DFT Quantum modelling. - Henry Rzepa&#039;s Blog","og_description":"X-ray crystallography is the technique of using the diffraction of x-rays by the electrons in a molecule to determine the positions of all the atoms in that molecule. Quantum theory teaches us that the electrons are to be found in shells around the atomic nuclei. There are two broad types, the outermost shell (also called [&hellip;]","og_url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137","og_site_name":"Henry Rzepa&#039;s Blog","article_published_time":"2025-03-17T19:57:21+00:00","article_modified_time":"2025-07-22T08:03:53+00:00","author":"Henry Rzepa, Andrew White","twitter_card":"summary_large_image","twitter_misc":{"Written by":"Henry Rzepa","Estimated reading time":"9 minutes"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"Article","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#article","isPartOf":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137"},"author":{"name":"Henry Rzepa","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#\/schema\/person\/2b40f7b9c872a4dc1547e040a11b6281"},"headline":"Crystallography meets DFT Quantum modelling.","datePublished":"2025-03-17T19:57:21+00:00","dateModified":"2025-07-22T08:03:53+00:00","mainEntityOfPage":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137"},"wordCount":1710,"commentCount":0,"image":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#primaryimage"},"thumbnailUrl":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/d0sc05526c-f1.jpg","articleSection":["Chemical IT"],"inLanguage":"en-GB","potentialAction":[{"@type":"CommentAction","name":"Comment","target":["https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#respond"]}]},{"@type":"WebPage","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137","url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137","name":"Crystallography meets DFT Quantum modelling. - Henry Rzepa&#039;s Blog","isPartOf":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#website"},"primaryImageOfPage":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#primaryimage"},"image":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#primaryimage"},"thumbnailUrl":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/d0sc05526c-f1.jpg","datePublished":"2025-03-17T19:57:21+00:00","dateModified":"2025-07-22T08:03:53+00:00","author":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#\/schema\/person\/2b40f7b9c872a4dc1547e040a11b6281"},"breadcrumb":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#breadcrumb"},"inLanguage":"en-GB","potentialAction":[{"@type":"ReadAction","target":["https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137"]}]},{"@type":"ImageObject","inLanguage":"en-GB","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#primaryimage","url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/d0sc05526c-f1.jpg","contentUrl":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/d0sc05526c-f1.jpg","width":530,"height":279},{"@type":"BreadcrumbList","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=28137#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Home","item":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog"},{"@type":"ListItem","position":2,"name":"Crystallography meets DFT Quantum modelling."}]},{"@type":"WebSite","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#website","url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/","name":"Henry Rzepa&#039;s Blog","description":"Chemistry with a twist","potentialAction":[{"@type":"SearchAction","target":{"@type":"EntryPoint","urlTemplate":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?s={search_term_string}"},"query-input":{"@type":"PropertyValueSpecification","valueRequired":true,"valueName":"search_term_string"}}],"inLanguage":"en-GB"},{"@type":"Person","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#\/schema\/person\/2b40f7b9c872a4dc1547e040a11b6281","name":"Henry Rzepa","image":{"@type":"ImageObject","inLanguage":"en-GB","@id":"https:\/\/secure.gravatar.com\/avatar\/897b6740f7f599bca7942cdf7d7914af5988937ae0e3869ab09aebb87f26a731?s=96&d=blank&r=g370be3a7397865e4fd161aefeb0a5a85","url":"https:\/\/secure.gravatar.com\/avatar\/897b6740f7f599bca7942cdf7d7914af5988937ae0e3869ab09aebb87f26a731?s=96&d=blank&r=g","contentUrl":"https:\/\/secure.gravatar.com\/avatar\/897b6740f7f599bca7942cdf7d7914af5988937ae0e3869ab09aebb87f26a731?s=96&d=blank&r=g","caption":"Henry Rzepa"},"description":"Henry Rzepa is Emeritus Professor of Computational Chemistry at Imperial College London.","sameAs":["https:\/\/orcid.org\/0000-0002-8635-8390"],"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?author=1"}]}},"jetpack_publicize_connections":[],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/pDef7-7jP","jetpack-related-posts":[{"id":2599,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=2599","url_meta":{"origin":28137,"position":0},"title":"Hypervalency: Is it real?","author":"Henry Rzepa","date":"October 16, 2010","format":false,"excerpt":"The Wikipedia page on hypervalent compounds reveals that the concept is almost as old as that of normally valent compounds. The definition there, \u00a0is \"a molecule that contains one or more\u00a0main group elements formally bearing more than eight\u00a0electrons in their\u00a0valence shells\" (although it could equally apply to e.g. transition elements\u2026","rel":"","context":"In &quot;Hypervalency&quot;","block_context":{"text":"Hypervalency","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=7"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2010\/10\/IH7.jpg?resize=350%2C200","width":350,"height":200},"classes":[]},{"id":2619,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=2619","url_meta":{"origin":28137,"position":1},"title":"Hypervalency: I(CN)7 is not hypervalent!","author":"Henry Rzepa","date":"October 17, 2010","format":false,"excerpt":"In the last post, IH7 was examined to see if it might exhibit true hypervalency. The iodine, despite its high coordination, turned out not to be hypervalent, with its (s\/p) valence shell not exceeding eight electrons (and its d-shell still with 10, and the 6s\/6p shells largely unoccupied). Instead, the\u2026","rel":"","context":"In &quot;Hypervalency&quot;","block_context":{"text":"Hypervalency","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=7"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2010\/10\/ICN7.jpg?resize=350%2C200","width":350,"height":200},"classes":[]},{"id":18993,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=18993","url_meta":{"origin":28137,"position":2},"title":"VSEPR Theory: Octet-busting or not with trimethyl chlorine, ClMe3.","author":"Henry Rzepa","date":"November 12, 2017","format":false,"excerpt":"A few years back, I took a look at the valence-shell electron pair repulsion approach to the geometry of chlorine trifluoride, ClF3 using so-called ELF basins to locate centroids for both the covalent F-Cl bond electrons and the chlorine lone-pair electrons. Whereas the original VSEPR theory talks about five \"electron\u2026","rel":"","context":"In &quot;Chemical IT&quot;","block_context":{"text":"Chemical IT","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=2"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":3956,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=3956","url_meta":{"origin":28137,"position":3},"title":"Beryllocene and Uranocene: The 8, 18 and 32-electron rules.","author":"Henry Rzepa","date":"April 25, 2011","format":false,"excerpt":"In discussing ferrocene in the previous post, I mentioned Irving Langmuir's 1921 postulate that filled valence shells in what he called\u00a0complete molecules would have magic numbers of 2, 8, 18 or 32 electrons (deriving from the sum of terms in 2[1+3+5+7]). The first two dominate organic chemistry of course, whilst\u2026","rel":"","context":"In &quot;Interesting chemistry&quot;","block_context":{"text":"Interesting chemistry","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=4"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/beryllocene-elf.jpg?resize=350%2C200","width":350,"height":200},"classes":[]},{"id":19102,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=19102","url_meta":{"origin":28137,"position":4},"title":"Hypervalence and octet-expansion in trimethylene-\u03bb6-sulfane and related species.","author":"Henry Rzepa","date":"November 27, 2017","format":false,"excerpt":"Previously: \"Non-polar\" species such as SeMe6, SMe6, ClMe3, ClMe5 all revealed interesting properties for the Se-C, S-C or Cl-C \"single\" bonds. The latter two examples in particular hinted at internal structures for these single bonds, as manifested by two ELF basins for some of the bonds. Here I take a\u2026","rel":"","context":"In &quot;Hypervalency&quot;","block_context":{"text":"Hypervalency","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=7"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":16758,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=16758","url_meta":{"origin":28137,"position":5},"title":"What&#8217;s in a name? Carbenes:  a reality check.","author":"Henry Rzepa","date":"September 11, 2016","format":false,"excerpt":"To quote from Wikipedia: in chemistry, a carbene is a molecule containing a neutral carbon atom with a valence of two and two unshared valence electrons. The most ubiquitous type of carbene of recent times is the one shown below as 1, often referred to as a resonance stabilised or\u2026","rel":"","context":"In &quot;crystal_structure_mining&quot;","block_context":{"text":"crystal_structure_mining","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=1745"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]}],"jetpack_likes_enabled":false,"authors":[{"term_id":2661,"user_id":1,"is_guest":0,"slug":"admin","display_name":"Henry Rzepa","avatar_url":"https:\/\/secure.gravatar.com\/avatar\/897b6740f7f599bca7942cdf7d7914af5988937ae0e3869ab09aebb87f26a731?s=96&d=blank&r=g","0":null,"1":"","2":"","3":"","4":"","5":"","6":"","7":"","8":""},{"term_id":2662,"user_id":8,"is_guest":0,"slug":"andrew-white","display_name":"Andrew White","avatar_url":{"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/AJPW-photo.jpg","url2x":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2025\/03\/AJPW-photo.jpg"},"0":null,"1":"","2":"","3":"","4":"","5":"","6":"","7":"","8":""}],"_links":{"self":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/28137","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=28137"}],"version-history":[{"count":33,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/28137\/revisions"}],"predecessor-version":[{"id":28365,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/28137\/revisions\/28365"}],"wp:attachment":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=28137"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=28137"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=28137"},{"taxonomy":"author","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fppma_author&post=28137"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}