Adaptive octree meshes for simulation of extracellular electrophysiology

Abstract

Objective. The interaction between neural tissues and artificial electrodes is crucial for understanding and advancing neuroscientific research and therapeutic applications. However, accurately modeling this space around the neurons rapidly increases the computational complexity of neural simulations. Approach. This study demonstrates a dynamically adaptive simulation method that greatly accelerates computation by adjusting spatial resolution of the simulation as needed. Use of an octree structure for the mesh, in combination with the admittance method for discretizing conductivity, provides both accurate approximation and ease of modification on-the-fly. Main results. In tests of both local field potential estimation and multi-electrode stimulation, dynamically adapted meshes achieve accuracy comparable to high-resolution static meshes in an order of magnitude less time. Significance. The proposed simulation pipeline improves model scalability, allowing greater detail with fewer computational resources. The implementation is available as an open-source Python module, providing flexibility and ease of reuse for the broader research community.