{"id":5290,"date":"2011-10-28T15:40:05","date_gmt":"2011-10-28T14:40:05","guid":{"rendered":"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=5290"},"modified":"2011-10-29T09:51:35","modified_gmt":"2011-10-29T08:51:35","slug":"computers-1967-2011-a-personal-perspective-part-4-moores-law-and-molecules","status":"publish","type":"post","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290","title":{"rendered":"Computers 1967-2011: a personal perspective. Part 4. Moore’s Law and Molecules."},"content":{"rendered":"
\n

Moore’s law<\/a><\/strong>\u00a0describes a long-term trend in the evolution of computing hardware, and it is often interpreted in terms of processing speed. Here I chart this rise in terms of the size of computable molecules. By computable I mean specifically how long it takes to predict the geometry of a given molecule using a quantum mechanical procedure.<\/p>\n

\"\"<\/a>

LSD, the 1975 benchmark for computable molecules.<\/p><\/div>\n

The geometry (shape) of a molecule is defined by 3N-6 variables, where N is the number of atoms it contains. Optimising the value of variables in order to obtain the minimum value of a function was first conducted by chemical engineers, who needed to improve the function of chemical reactor plants. The mathematical techniques they developed were adopted to molecules in the 1970s, and in 1975<\/strong> a milestone was reached with the molecule above. Here, N=49, and 3N-6=141. The function used was one describing its computed enthalpy of formation, using a quantum mechanical procedure known as MINDO\/3. The computer used was what passed then for a supercomputer, a CDC 6600 (of which a large well endowed university could probably afford one of). It was almost impossible to get exclusive access to such a beast (its computing power was shared amongst the entire university, in this case of about 50,000 people), but during a slack period over a long weekend, the optimised geometry of LSD was obtained (it’s difficult to know how many hours the CDC 6600 took to perform this feat, but I suspect it might have been around 72). The result was announced by Paul Weiner to the group I was then part of (the Dewar research group), and Michael immediately announced that this deserved an unusual Monday night sojourn to the Texas Tavern, where double pitchers of beer would be available. You might be tempted to ask what the reason for the celebration was. Well, LSD was a “real molecule” (and not a hallucination). It meant one could predict for the first time the geometry of realistic molecules such as drugs and hence be taken seriously by people who dealt with molecules of this size for a living. And if you could predict the energy of its equilibrium geometry, you could then quickly move on to predicting the barriers to its reaction. A clear tipping point had been reached in computational simulation.<\/p>\n

In 1975, MINDO\/3 was thought to compute an energy function around 1000 to 10,000 faster than the supposedly more accurate ab initio<\/em> codes then available (in fact you could not then routinely optimise geometries with the common codes of this type). With this in mind, one can subject the same molecule to a modern \u03c9B97XD\/6-311G(d,p) optimisation. This level of theory is probably closer to 104<\/sup> to 105<\/sup> times slower to compute than MINDO\/3. On a modest “high performance” resource (which nowadays runs in parallel, in fact on 32 cores in this case), the calculation<\/a> takes about an hour (starting from a 1973 X-ray structure, which turns out to be quite a poor place to start from, and almost certainly poorer than the 1975 point). In (very) round numbers, the modern calculation is about a million times faster. Which (coincidentally) is approximately the factor predicted by Moore’s law.<\/p>\n

I will give one more example, this time for an example dating from around 2003<\/strong>, 28 years on from the original benchmark.<\/p>\n

\"\"<\/a>

Transition state for lactide polymerisation.<\/p><\/div>\n

This example has 114 atoms, and hence 3N-6 =336, or 2.42 times the 1975 size. It is a transition state, which is a far slower calculation then an equilibrium geometry. It is also typical of the polymerisation chemistry of the naughties. Each run on the computer (B3LYP\/6-31G(d), with the alkyl groups treated at STO-3G) now took about 8-10 days (on a machine with 4 cores), and probably 2-4 runs in total would have been required per system (of which four needed to be studied to derive meaningful conclusions). Let us say 1000 hours per transition state. Together with false starts etc, the project took about 18 months to complete. Move on to 2010<\/strong>; added to the model was a significantly better (= slower) basis set and a solvation correction, and a single calculation now took 67 hours<\/a>. In 2011<\/strong>, it would be reduced to ~10 hours (by now we are up to 64-core computers).<\/p>\n

In 2011, calculations involving ~250 atoms are now regarded as almost routine, and molecules with up to this number of atoms cover most of the discrete (i.e.<\/em> non repeating) molecular systems of interest nowadays. But the 1975 LSD calculation still stands as the day that realistic computational chemistry came of age<\/a>.<\/p>\n\n<\/div> ","protected":false},"excerpt":{"rendered":"

Moore’s law\u00a0describes a long-term trend in the evolution of computing hardware, and it is often interpreted in terms of processing speed. Here I chart this rise in terms of the size of computable molecules. By computable I mean specifically how long it takes to predict the geometry of a given molecule using a quantum mechanical […]<\/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":false,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[4],"tags":[183,717,720,34,24,718,716,2651,712,719,590,721,723,722,42],"ppma_author":[2661],"class_list":["post-5290","post","type-post","status-publish","format-standard","hentry","category-interesting-chemistry","tag-3g","tag-chemical-engineers","tag-chemical-reactor-plants","tag-computational-chemistry","tag-energy","tag-energy-function","tag-hallucination","tag-historical","tag-lsd","tag-molecular-systems","tag-paul-weiner","tag-simulation","tag-sojourn","tag-texas-tavern","tag-x-ray"],"yoast_head":"\nComputers 1967-2011: a personal perspective. Part 4. Moore's Law and Molecules. - 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=5290\" \/>\n<meta property=\"og:locale\" content=\"en_GB\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Computers 1967-2011: a personal perspective. Part 4. Moore's Law and Molecules. - Henry Rzepa's Blog\" \/>\n<meta property=\"og:description\" content=\"Moore’s law\u00a0describes a long-term trend in the evolution of computing hardware, and it is often interpreted in terms of processing speed. Here I chart this rise in terms of the size of computable molecules. By computable I mean specifically how long it takes to predict the geometry of a given molecule using a quantum mechanical […]\" \/>\n<meta property=\"og:url\" content=\"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290\" \/>\n<meta property=\"og:site_name\" content=\"Henry Rzepa's Blog\" \/>\n<meta property=\"article:published_time\" content=\"2011-10-28T14:40:05+00:00\" \/>\n<meta property=\"article:modified_time\" content=\"2011-10-29T08:51:35+00:00\" \/>\n<meta property=\"og:image\" content=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/10\/LSD.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=\"4 minutes\" \/>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Computers 1967-2011: a personal perspective. Part 4. Moore's Law and Molecules. - 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=5290","og_locale":"en_GB","og_type":"article","og_title":"Computers 1967-2011: a personal perspective. Part 4. Moore's Law and Molecules. - Henry Rzepa's Blog","og_description":"Moore’s law\u00a0describes a long-term trend in the evolution of computing hardware, and it is often interpreted in terms of processing speed. Here I chart this rise in terms of the size of computable molecules. By computable I mean specifically how long it takes to predict the geometry of a given molecule using a quantum mechanical […]","og_url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290","og_site_name":"Henry Rzepa's Blog","article_published_time":"2011-10-28T14:40:05+00:00","article_modified_time":"2011-10-29T08:51:35+00:00","og_image":[{"url":"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/10\/LSD.jpg","type":"","width":"","height":""}],"author":"Henry Rzepa","twitter_card":"summary_large_image","twitter_misc":{"Written by":"Henry Rzepa","Estimated reading time":"4 minutes"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"Article","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290#article","isPartOf":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290"},"author":{"name":"Henry Rzepa","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/#\/schema\/person\/2b40f7b9c872a4dc1547e040a11b6281"},"headline":"Computers 1967-2011: a personal perspective. Part 4. Moore’s Law and Molecules.","datePublished":"2011-10-28T14:40:05+00:00","dateModified":"2011-10-29T08:51:35+00:00","mainEntityOfPage":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290"},"wordCount":777,"commentCount":0,"image":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290#primaryimage"},"thumbnailUrl":"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/10\/LSD.jpg","keywords":["3g","chemical engineers","chemical reactor plants","computational chemistry","energy","energy function","hallucination","Historical","LSD","molecular systems","Paul Weiner","simulation","sojourn","Texas Tavern","X-ray"],"articleSection":["Interesting chemistry"],"inLanguage":"en-GB","potentialAction":[{"@type":"CommentAction","name":"Comment","target":["https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290#respond"]}]},{"@type":"WebPage","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290","url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290","name":"Computers 1967-2011: a personal perspective. Part 4. Moore's Law and Molecules. - 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=5290#primaryimage"},"image":{"@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290#primaryimage"},"thumbnailUrl":"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/10\/LSD.jpg","datePublished":"2011-10-28T14:40:05+00:00","dateModified":"2011-10-29T08:51:35+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=5290#breadcrumb"},"inLanguage":"en-GB","potentialAction":[{"@type":"ReadAction","target":["https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290"]}]},{"@type":"ImageObject","inLanguage":"en-GB","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290#primaryimage","url":"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/10\/LSD.jpg","contentUrl":"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/10\/LSD.jpg"},{"@type":"BreadcrumbList","@id":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=5290#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Home","item":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog"},{"@type":"ListItem","position":2,"name":"Computers 1967-2011: a personal perspective. Part 4. Moore’s Law and Molecules."}]},{"@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-1nk","jetpack-related-posts":[{"id":20601,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=20601","url_meta":{"origin":5290,"position":0},"title":"Impossible molecules.","author":"Henry Rzepa","date":"April 1, 2019","format":false,"excerpt":"Members of the chemical FAIR data community have just met in Orlando (with help from the NSF, the American National Science Foundation)\u00a0to discuss how such data is progressing in chemistry. There are a lot of themes converging at the moment. Thus this article extolls the virtues of having raw NMR\u2026","rel":"","context":"In "Chemical IT"","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":19550,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=19550","url_meta":{"origin":5290,"position":1},"title":"A record polarity for a neutral compound?","author":"Henry Rzepa","date":"April 13, 2018","format":false,"excerpt":"In several posts a year or so ago I considered various suggestions for the most polar neutral molecules, as measured by the dipole moment. A record had been claimed for a synthesized molecule of ~14.1\u00b10.7D. I pushed this to a calculated 21.7D for an admittedly hypothetical and unsynthesized molecule. Here\u2026","rel":"","context":"In "Interesting chemistry"","block_context":{"text":"Interesting chemistry","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=4"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":221,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=221","url_meta":{"origin":5290,"position":2},"title":"How do molecules interact with each other?","author":"Henry Rzepa","date":"April 12, 2009","format":false,"excerpt":"Understanding how molecules interact (bind) with each other when in close proximity is essential in all areas of chemistry. One specific example of this need is for the molecule shown below. This is the so-called Pirkle Reagent and is much used to help resolve the two enantiomers of a racemic\u2026","rel":"","context":"In "Interesting chemistry"","block_context":{"text":"Interesting chemistry","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=4"},"img":{"alt_text":"The Pirkle reagent","src":"https:\/\/i0.wp.com\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2009\/04\/pirkle.jpg?resize=350%2C200","width":350,"height":200},"classes":[]},{"id":20778,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=20778","url_meta":{"origin":5290,"position":3},"title":"Imaging normal vibrational modes of a single molecule of CoTPP: a mystery about the nature of the imaged species.","author":"Henry Rzepa","date":"April 25, 2019","format":false,"excerpt":"Previously, I explored (computationally) the normal vibrational modes of Co(II)-tetraphenylporphyrin (CoTPP) as a \"flattened\" species on copper or gold surfaces for comparison with those recently imaged. The initial intent was to estimate the \"flattening\" energy. There are six electronic possibilities for this molecule on a metal surface. Respectively positively, or\u2026","rel":"","context":"In "Interesting chemistry"","block_context":{"text":"Interesting chemistry","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=4"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":1135,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=1135","url_meta":{"origin":5290,"position":4},"title":"The SN1 Reaction- revisited","author":"Henry Rzepa","date":"November 11, 2009","format":false,"excerpt":"In an earlier post I wrote about the iconic SN1 solvolysis reaction, and presented a model for the transition state involving 13 water molecules. Here, I follow this up with an improved molecule containing 16 water molecules, and how the barrier for this model compares with experiment. This latter is\u2026","rel":"","context":"In "Interesting chemistry"","block_context":{"text":"Interesting chemistry","link":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?cat=4"},"img":{"alt_text":"Transition state for Sn1 solvolysis of tert-butyl chloride","src":"https:\/\/i0.wp.com\/www.ch.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2009\/11\/sn1.jpg?resize=350%2C200","width":350,"height":200},"classes":[]},{"id":20679,"url":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/?p=20679","url_meta":{"origin":5290,"position":5},"title":"Imaging vibrational normal modes of a single molecule.","author":"Henry Rzepa","date":"April 18, 2019","format":false,"excerpt":"The topic of this post originates from a recent article which is attracting much attention. The technique uses confined light to both increase the spatial resolution by around three orders of magnitude and also to amplify the signal from individual molecules to the point it can be recorded. To me,\u2026","rel":"","context":"In "Interesting chemistry"","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.ic.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2019\/04\/Figure3a-1024x529.jpg?resize=350%2C200&ssl=1","width":350,"height":200},"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":""}],"_links":{"self":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/5290","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=5290"}],"version-history":[{"count":0,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=\/wp\/v2\/posts\/5290\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=5290"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=5290"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=5290"},{"taxonomy":"author","embeddable":true,"href":"https:\/\/www.ch.ic.ac.uk\/rzepa\/blog\/index.php?rest_route=%2Fwp%2Fv2%2Fppma_author&post=5290"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}