{"id":23360,"date":"2021-02-07T10:12:23","date_gmt":"2021-02-07T10:12:23","guid":{"rendered":"https:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=23360"},"modified":"2021-03-03T09:20:40","modified_gmt":"2021-03-03T09:20:40","slug":"non-covalent-interaction-nci-surfaces-for-two-large-annulenes-revisited","status":"publish","type":"post","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360","title":{"rendered":"Non-covalent-interaction (NCI) surfaces for two large annulenes (revisited)."},"content":{"rendered":"
\n

The last post<\/a> addressed the concept of “steric clashes” in a pericyclic reaction transition state as an extension of the time honoured<\/a> practice of building\u00a0molecular models to analyse reaction outcomes. A modern computer generated model might express this in terms of a NCI (non-covalent-interaction) surface. A few posts ago, I had looked at some “molecules of the year” for 2020, one of which was a “figure-eight” twisted dodecaporphyrin<\/a> in which an aspect of\u00a0the reported[1]<\/a><\/span> geometry had struck me as potentially lacking features due to the so-called non-covalent dispersion or\u00a0van der Waals<\/a> attractions. So I am revisiting here by adding the NCI surface for this molecule and one other.<\/p>\n

The molecule in question has 720 atoms and can be regarded as a [162]-annulene (4n+2, n=40) with a linking number Lk =2\u03c0.[2]<\/a><\/span> The NCI surface is ideally computed from a “self-consistent-field or\u00a0SCF density” and so in this instance I used the PM7 SCF-density, which is derived from the valence shell only and does not include the core shells. That hardly matters since the non-covalent NCI surface does not use the core shells!<\/p>\n

\"\"

Click for 3D model<\/p><\/div>\n

You can see from the above that the porphyrin-stacking region has a very dense NCI green surface (arrow), indicating a lot of stabilisation is originating there; something lacking in the original proposed structure. There are lots of other features and\u00a0so I do encourage you to explore the 3D model.<\/p>\n

The second (hypothetical) molecule is a simpler CH-based [144]-annulene,[3]<\/a><\/span> comprising a twisted coil of 144 CH= units<\/a> with a linking number Lk = 18\u03c0 (the largest such ever proposed for a molecule!). The SCF-NCI surface (derived from an \u03c9B97XD\/6-31G(d,p) calculation) is contiguous all the way around the circuit and must be the ultimate \u03c0-\u03c0 stacked molecule!<\/p>\n

\"\"

144-annulene. Click for 3D<\/p><\/div>\n

\"\"

153-annulene. Click for 3D<\/p><\/div>\n

I should end with a brief tutorial on how to generate these surfaces. You need a density matrix (e.g.<\/em> DOI: 10.14469\/ch\/16967<\/a>). In programs such as Gaussian 16, this can be obtained from the checkpoint file, which contains it. A progam called cubgen is used by (e.g.<\/em>) Gaussian to create a 3D cube of electron density values (as well as other interesting properties). To get good resolution (~ 0.044\u00c5) the file will be between 500 – 800 Mbyte in size. If a resolution of ~ 0.088\u00c5 is used it will be eight times smaller. For cube files less than ~105 Mbyte in size, you can use this\u00a0Web-based tool (DOI: 10.14469\/hpc\/7864<\/a>) to get the NCI surface. For the larger files you will need the Jmol application<\/a> which can sustain files up to ~1 Gbyte (or larger, but\u00a0I have not tested) and where you will run the following script:\u00a0<\/p>\n

load density.cub;isosurface parameters [0.5 1 0.0005 0.05 0.95 1.00] NCI \"\";color isosurface \"bgyor\" range -0.04 0.04;write density.xyz;write density.jvxl;<\/tt><\/span><\/p>\n

\n(where this script assumes that the file density.cub<\/tt> file is in the same folder as the Jmol.java<\/tt> application).<\/p>\n

\n

Postscript<\/b>: The NCI analysis is based on computing the total density of the molecule. Close inspection of the top molecule as computed using the semi-empirical method PM7 reveals some interesting features extending beyond the C-H bonds. Analysis of this reveals it to be an artefact of the computed density, itself traced back to differences in how overlaps are handled in computing the density for this particular method. This error is not present for the MNDO semi-empirical method. When evaluated using MNDO, but at the geometry computed by PM7, these artefacts are removed. The NCI feature in the π-π stacking shown above however remains and hence is not an artefact.<\/p>\n

\"\"<\/p>\n

References<\/h2>\n
  1. M. Rickhaus, M. Jirasek, L. Tejerina, H. Gotfredsen, M.D. Peeks, R. Haver, H. Jiang, T.D.W. Claridge, and H.L. Anderson, \"Global aromaticity at the nanoscale\", Nature Chemistry<\/i>, vol. 12, pp. 236-241, 2020. https:\/\/doi.org\/10.1038\/s41557-019-0398-3<\/a>\n\n<\/li>\n
  2. S.M. Rappaport, and H.S. Rzepa, \"Intrinsically Chiral Aromaticity. Rules Incorporating Linking Number, Twist, and Writhe for Higher-Twist M\u00f6bius Annulenes\", Journal of the American Chemical Society<\/i>, vol. 130, pp. 7613-7619, 2008. https:\/\/doi.org\/10.1021\/ja710438j<\/a>\n\n<\/li>\n
  3. R.J.F. Berger, \"Prediction of a Cyclic Helical Oligoacetylene Showing Anapolar Ring Currents in the Magnetic Field\", Zeitschrift f\u00fcr Naturforschung B<\/i>, vol. 67, pp. 1127-1131, 2012. https:\/\/doi.org\/10.5560\/znb.2012-0189<\/a>\n\n<\/li>\n<\/ol>\n\n<\/div> ","protected":false},"excerpt":{"rendered":"

    The last post addressed the concept of “steric clashes” in a pericyclic reaction transition state as an extension of the time honoured practice of building\u00a0molecular models to analyse reaction outcomes. A modern computer generated model might express this in terms of a NCI (non-covalent-interaction) surface. A few posts ago, I had looked at some “molecules […]<\/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":[4],"tags":[],"ppma_author":[2661],"class_list":["post-23360","post","type-post","status-publish","format-standard","hentry","category-interesting-chemistry"],"yoast_head":"\nNon-covalent-interaction (NCI) surfaces for two large annulenes (revisited). - Henry Rzepa'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=23360\" \/>\n<meta property=\"og:locale\" content=\"en_GB\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Non-covalent-interaction (NCI) surfaces for two large annulenes (revisited). - Henry Rzepa's Blog\" \/>\n<meta property=\"og:description\" content=\"The last post addressed the concept of “steric clashes” in a pericyclic reaction transition state as an extension of the time honoured practice of building\u00a0molecular models to analyse reaction outcomes. A modern computer generated model might express this in terms of a NCI (non-covalent-interaction) surface. A few posts ago, I had looked at some “molecules […]\" \/>\n<meta property=\"og:url\" content=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360\" \/>\n<meta property=\"og:site_name\" content=\"Henry Rzepa's Blog\" \/>\n<meta property=\"article:published_time\" content=\"2021-02-07T10:12:23+00:00\" \/>\n<meta property=\"article:modified_time\" content=\"2021-03-03T09:20:40+00:00\" \/>\n<meta property=\"og:image\" content=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2021\/02\/new.jpg\" \/>\n<meta name=\"author\" content=\"Henry Rzepa\" \/>\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=\"3 minutes\" \/>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Non-covalent-interaction (NCI) surfaces for two large annulenes (revisited). - Henry Rzepa'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=23360","og_locale":"en_GB","og_type":"article","og_title":"Non-covalent-interaction (NCI) surfaces for two large annulenes (revisited). - Henry Rzepa's Blog","og_description":"The last post addressed the concept of “steric clashes” in a pericyclic reaction transition state as an extension of the time honoured practice of building\u00a0molecular models to analyse reaction outcomes. A modern computer generated model might express this in terms of a NCI (non-covalent-interaction) surface. A few posts ago, I had looked at some “molecules […]","og_url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360","og_site_name":"Henry Rzepa's Blog","article_published_time":"2021-02-07T10:12:23+00:00","article_modified_time":"2021-03-03T09:20:40+00:00","og_image":[{"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2021\/02\/new.jpg","type":"","width":"","height":""}],"author":"Henry Rzepa","twitter_card":"summary_large_image","twitter_misc":{"Written by":"Henry Rzepa","Estimated reading time":"3 minutes"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"Article","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360#article","isPartOf":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360"},"author":{"name":"Henry Rzepa","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#\/schema\/person\/2b40f7b9c872a4dc1547e040a11b6281"},"headline":"Non-covalent-interaction (NCI) surfaces for two large annulenes (revisited).","datePublished":"2021-02-07T10:12:23+00:00","dateModified":"2021-03-03T09:20:40+00:00","mainEntityOfPage":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360"},"wordCount":633,"commentCount":0,"image":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360#primaryimage"},"thumbnailUrl":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2021\/02\/new.jpg","articleSection":["Interesting chemistry"],"inLanguage":"en-GB","potentialAction":[{"@type":"CommentAction","name":"Comment","target":["https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360#respond"]}]},{"@type":"WebPage","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360","url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360","name":"Non-covalent-interaction (NCI) surfaces for two large annulenes (revisited). - Henry Rzepa's Blog","isPartOf":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#website"},"primaryImageOfPage":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360#primaryimage"},"image":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360#primaryimage"},"thumbnailUrl":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2021\/02\/new.jpg","datePublished":"2021-02-07T10:12:23+00:00","dateModified":"2021-03-03T09:20:40+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=23360#breadcrumb"},"inLanguage":"en-GB","potentialAction":[{"@type":"ReadAction","target":["https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360"]}]},{"@type":"ImageObject","inLanguage":"en-GB","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360#primaryimage","url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2021\/02\/new.jpg","contentUrl":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2021\/02\/new.jpg","width":601,"height":612,"caption":"Click for 3D model"},{"@type":"BreadcrumbList","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=23360#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Home","item":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog"},{"@type":"ListItem","position":2,"name":"Non-covalent-interaction (NCI) surfaces for two large annulenes (revisited)."}]},{"@type":"WebSite","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#website","url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/","name":"Henry Rzepa'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-64M","jetpack-related-posts":[],"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":""}],"_links":{"self":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/23360","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=23360"}],"version-history":[{"count":22,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/23360\/revisions"}],"predecessor-version":[{"id":23408,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/23360\/revisions\/23408"}],"wp:attachment":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=23360"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=23360"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=23360"},{"taxonomy":"author","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fppma_author&post=23360"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}