Researchers at the Simons Foundation Flatiron Institute have developed a novel computational method that accelerates molecular dynamics simulations by a factor of two to seven, delivering significant gains in processing speed and energy efficiency. The breakthrough, detailed in a paper published online on May 21 in Nature Communications, addresses a longstanding bottleneck in high-performance computing: the calculation of long-range electrostatic forces between atoms. Molecular dynamics simulations model the temporal evolution of millions of interacting atoms, a task essential for advancements in drug discovery, protein folding analysis, and next-generation battery research. Historically, these simulations require trillions of computational steps to capture atomic vibrations, with force calculations scaling quadratically relative to atom count. Even with established algorithms like the fast multipole method, electrostatic computations have consumed over twenty percent of the workload on the world’s top five hundred supercomputers. To overcome this limitation, the Flatiron Institute Center for Computational Mathematics team turned to an underutilized mathematical tool: prolate spheroidal wave functions, originally devised in the nineteenth century and later adapted for signal processing. Lead author Jiuyang Liang and senior author Shidong Jiang, working alongside software engineer Libin Lu, project leader Alex Barnett, and center director Leslie Greengard, applied these functions to optimize two critical simulation parameters. The method determines how to partition electrostatic interactions into short and long range components while simultaneously mapping atomic charges onto computational grids. Prolate spheroidal functions satisfy the competing mathematical demands of spatial localization and smoothness far more efficiently than traditional approaches. The implementation integrates directly into leading molecular dynamics platforms, including GROMACS, LAMMPS, and OpenMM. Performance benchmarks across diverse systems, including water networks, immune-related proteins, and lithium-ion battery electrolytes, demonstrate consistent speedups ranging from two and a half to seven times, with GROMACS achieving a fivefold increase at high precision. The code has already been officially incorporated into the LAMMPS distribution, signaling immediate adoption potential. Senior research scientist Pilar Cossio and computational chemist Sonya Hanson, both affiliated with the Flatiron Institute, emphasized that the optimization drastically reduces the temporal constraints that previously limited simulations to mere microseconds. By lowering computational overhead, the method promises to cut both energy consumption and processing time across multiple scientific disciplines. Anthony Costa of Nvidia noted that the advancement represents a rare decade-spanning leap in simulation performance, highlighting the practical value of foundational mathematical research. The Flatiron Institute cross-disciplinary approach underscores a growing trend in computational science, where abstract mathematical theory is leveraged to resolve hardware-scale bottlenecks. As high-performance computing continues to power materials design and biomedical research, this algorithmic refinement positions molecular dynamics to operate at unprecedented scales without requiring proportional increases in computational infrastructure.