Document Type


Publication Date



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.


This is a pre-copy-editing, author-produced PDF of an article accepted for publication in Journal of Neural Engineering, volume 20, in 2023 following peer review. This article may not exactly replicate the final published version. The definitive publisher-authenticated version is available online at

The Creative Commons license below applies only to this version of the article.


The authors

Creative Commons License

Creative Commons License
This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 4.0 License.

Available for download on Wednesday, September 18, 2024