{"id":137296,"date":"2026-06-24T09:30:17","date_gmt":"2026-06-24T13:30:17","guid":{"rendered":"https:\/\/www.simonsfoundation.org\/?p=137296"},"modified":"2026-06-21T17:13:32","modified_gmt":"2026-06-21T21:13:32","slug":"mathematicians-unleash-multifold-speed-boost-for-supercomputer-simulations-of-molecules","status":"publish","type":"post","link":"https:\/\/www.simonsfoundation.org\/2026\/06\/24\/mathematicians-unleash-multifold-speed-boost-for-supercomputer-simulations-of-molecules\/","title":{"rendered":"Mathematicians Unleash Multifold Speed Boost for Supercomputer Simulations of Molecules"},"content":{"rendered":"","protected":false},"excerpt":{"rendered":"","protected":false},"author":432,"featured_media":137298,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"_relevanssi_hide_post":"","_relevanssi_hide_content":"","_relevanssi_pin_for_all":"","_relevanssi_pin_keywords":"","_relevanssi_unpin_keywords":"","_relevanssi_related_keywords":"","_relevanssi_related_include_ids":"","_relevanssi_related_exclude_ids":"","_relevanssi_related_no_append":"","_relevanssi_related_not_related":"","_relevanssi_related_posts":"","_relevanssi_noindex_reason":"","footnotes":""},"categories":[451],"tags":[],"class_list":["post-137296","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-ccm","news_type-press-release"],"acf":{"block_editor":[{"acf_fc_layout":"video","embed_code":"https:\/\/vimeo.com\/1201803836?share=copy&fl=sv&fe=ci","title":"An atomistic molecular dynamics simulation of a dense ionic liquid made of LiTFSI \u2014 a key lithium salt used to study next-generation battery electrolytes. Each sphere represents an atom, with colors distinguishing lithium ions and atoms in the TFSI anions. The motion illustrates how ions and their local environments continually rearrange and interact within the electrolyte. The movie was made from a molecular dynamics simulation of a 1 million-atom system.","caption":"Jiuyang Liang\/Flatiron Institute","video_autoplay":true,"video_loop":true,"video_width":695,"video_height":390,"video_container_size":"video_large"},{"acf_fc_layout":"text","text":"More than 20 percent of the workload on the world\u2019s 500 fastest supercomputers is spent simulating how atoms and molecules move \u2014 with applications ranging from material design to identifying drug interactions to understanding protein folding.\r\n\r\nBy leveraging a classical mathematical function, researchers at the Simons Foundation\u2019s <a href=\"https:\/\/www.simonsfoundation.org\/flatiron\/\">Flatiron Institute<\/a> developed a new method that enables these simulations to run 2.5 to seven times as fast \u2014 all without sacrificing accuracy. For the most popular software package for molecular dynamics simulations, <a href=\"https:\/\/www.gromacs.org\/\">GROMACS<\/a>, the researchers achieved a fivefold speed increase in simulations run at high accuracy.\r\n\r\nThe new method can be rapidly and easily integrated into existing software workflows, meaning the field could soon see drastically reduced time and energy demands for the simulations. The researchers present their breakthrough in a paper <a href=\"https:\/\/www.nature.com\/articles\/s41467-026-73232-8\">published online May 21 in <em>Nature Communications<\/em><\/a>.\r\n\r\n\u201cThere are so many fields in science that rely on molecular simulations that can now take less energy and computing time,\u201d says the study\u2019s senior author, <a href=\"https:\/\/www.simonsfoundation.org\/people\/shidong-jiang\/\">Shidong Jiang<\/a>, a senior research scientist at the Flatiron Institute\u2019s Center for Computational Mathematics (CCM). \u201cWe think this work is going to have a broad impact.\u201d\r\n\r\nThe CCM\u2019s new work \u201cholds tremendous potential to meaningfully accelerate molecular dynamics workloads \u2014 a field that for the past decade has only seen incremental improvements in speed,\u201d says Anthony Costa, director of digital biology at Nvidia, who was not involved in the study. \u201cThe work is a testament to the importance of applied mathematical research and its enormous impact in multiple domains, including but not limited to life sciences and materials science.\u201d\r\n\r\nJiang worked on the study with lead author <a href=\"https:\/\/www.simonsfoundation.org\/people\/jiuyang-liang\/\">Jiuyang Liang<\/a> \u2014 a CCM affiliate research fellow and a researcher at Shanghai Jiao Tong University \u2014 as well as <a href=\"https:\/\/www.simonsfoundation.org\/people\/libin-lu\/\">Libin Lu<\/a>, a CCM software engineer; <a href=\"https:\/\/www.simonsfoundation.org\/people\/alex-barnett\/\">Alex Barnett<\/a>, a CCM project leader; and CCM Director <a href=\"https:\/\/www.simonsfoundation.org\/people\/leslie-greengard\/\">Leslie Greengard<\/a>.\r\n<h1>Dancing Molecules<\/h1>\r\nMolecular dynamics simulations model the behavior of millions of molecules in a box, where each molecule is composed of multiple atoms. \u201cWe try to simulate the evolution of the system \u2014 say a protein in a box of water \u2014 as a function of time,\u201d says <a href=\"https:\/\/www.simonsfoundation.org\/people\/pilar-cossio\/\">Pilar Cossio<\/a>, a Flatiron Institute senior research scientist who uses molecular dynamics simulations as part of her work.\r\n\r\nThe challenge comes with the time scales involved. Like a video game rendering an environment, molecular dynamics simulations chop time up into small slices. While a video game might run at 60 video frames per second of gameplay, capturing the vibrations of molecular bonds requires around 500 trillion time slices per second.\r\n\r\nThat quickly adds up, even when most of the frequently studied scientific processes take just a few microseconds or milliseconds. For a simulation to be meaningful, it can require as many as a trillion steps. \u201cEven with great hardware, we are limited to hundreds of nanoseconds a day, which gets us to microseconds in weeks,\u201d says <a href=\"https:\/\/www.simonsfoundation.org\/people\/sonya-hanson\/\">Sonya Hanson<\/a>, a Flatiron Institute research scientist who also works with molecular dynamics simulations.\r\n<h1>Speeding Things Up<\/h1>\r\nIt\u2019s not hard to imagine why such simulations are so computationally intensive. The simulations must account for the pushing and pulling of the charged particles that make up the molecules. These long-range \u2018electrostatic forces\u2019 mean that \u201cwe basically have to calculate the distances between all of the atoms in this big box,\u201d Cossio says.\r\n\r\n\u201cIf you do this naively, you end up with a number of operations equal to the number of atoms squared,\u201d says Jiang. There are shortcuts, though. Mathematicians have invented and adapted algorithms such as the fast Fourier transform and, later, the fast multipole method (co-invented by Greengard), both of which reduce the required operations. Even with those advancements, long-range force calculations in molecular dynamics simulations have continued to consume vast computational time and resources.\r\n\r\n\u201cWe saw potential to improve a tool that has been unquestioned for decades,\u201d Barnett says."},{"acf_fc_layout":"image","type":"image_medium","image":137302,"title":"","caption":"Lucy Reading-Ikkanda\/Simons Foundation","caption_full":"","less_margin":false},{"acf_fc_layout":"text","text":"The CCM researchers turned to a set of mathematical functions known as prolate spheroidal wave functions, which were invented in 1880 and later applied to signal processing at Bell Labs in the 1960s. The functions make two critical decisions in the simulations: how to split electrostatic interactions between atoms into short-range and long-range components and how to spread atomic charges onto a grid for the long-range force calculation.\r\n\r\nFor these tasks, the functions need to be both spatially localized (confined to a narrow region) and as smooth as possible (having the minimum allowed spatial frequency range, or \u2018bandlimit\u2019). These are seemingly conflicting demands that the prolate function satisfies much better than previously used functions.\r\n\r\n\u201cIt took knowing the mathematical literature about the right way to do splits and spreading,\u201d combined with knowing about the bottlenecks in molecular dynamics simulations, Barnett says.\r\n\r\n\u201cSomething I learned from Leslie [Greengard],\u201d Jiang says, \u201cis that as computational mathematicians, our job is to give the user the most accurate solution in the shortest time. Sometimes an advance in one field \u2014 in this case, molecular dynamics \u2014 comes from people who are outside the discipline.\u201d\r\n\r\nThe researchers ensured that the new code could be seamlessly integrated into the code of the most widely used software packages, such as LAMMPS, GROMACS and OpenMM.\r\n\r\n\u201cWe needed to convince people that our method can beat the best ones out there,\u201d says Jiang, adding that their code has officially been accepted by the LAMMPS developers.\r\n<h1>Record-Breaking Speed<\/h1>\r\nThe researchers tested the functions to simulate molecular dynamics in several different systems, including a collection of water molecules, a protein system important in immunity, and a lithium-ion solution used in batteries. In each case, the new method completed the tests 2.5 to seven times as fast.\r\n\r\n\u201cWe did have this lurking fear that even though this was clearly a mathematically better idea, what if we missed something in thinking about applying it? But it worked out as we anticipated,\u201d Greengard says.\r\n\r\nThe findings underscore the importance of computational mathematics and out-of-the-box thinking. Says Greengard, \u201cThe Flatiron is a place where we have both the freedom to carry out and the recognition for doing this kind of interdisciplinary work.\u201d"},{"acf_fc_layout":"subtitle","title":"Information for Press","is_chapter":false,"is_inline":false,"margin_bottom":""},{"acf_fc_layout":"text","text":"For more information, please contact <a href=\"mailto:press@simonsfoundation.org\" target=\"_blank\" rel=\"noopener\">press@simonsfoundation.org<\/a>.\r\n<ul>\r\n \t<li><a href=\"https:\/\/www.nature.com\/articles\/s41467-026-73232-8\">Link to scientific paper<\/a><\/li>\r\n \t<li><a href=\"https:\/\/sf-web-assets-prod.s3.amazonaws.com\/wp-content\/uploads\/2026\/06\/16081114\/FI_MolecularDynamics_c.Jiuyang-Liang_Flatiron-Institute.jpg\" rel=\"noopener\">Link to high-resolution illustration<\/a><\/li>\r\n \t<li><a href=\"https:\/\/sf-web-assets-prod.s3.amazonaws.com\/wp-content\/uploads\/2026\/06\/16085148\/Graphic_SpeedBumps_c.Lucy-Reading-Ikkanda-Simons-Foundation-scaled.jpg\" target=\"_blank\" rel=\"noopener\">Link to high-resolution infographic<\/a><\/li>\r\n<\/ul>"}]},"_links":{"self":[{"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/posts\/137296"}],"collection":[{"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/users\/432"}],"replies":[{"embeddable":true,"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/comments?post=137296"}],"version-history":[{"count":12,"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/posts\/137296\/revisions"}],"predecessor-version":[{"id":137733,"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/posts\/137296\/revisions\/137733"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/media\/137298"}],"wp:attachment":[{"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/media?parent=137296"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/categories?post=137296"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.simonsfoundation.org\/wp-json\/wp\/v2\/tags?post=137296"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}