Modeling Pyramidal L2/3 Neurons by Compressing Adjoint BEM Using H-Matrix Approximation

Electromagnetic brain stimulation uses electrodes or coils to induce electric fields (E-fields) in the brain and affect its activity. Our understanding of the precise effects of the device-induced E-fields on neural activity is limited. In this work, we present a novel boundary integral equation that enables the modeling of fully coupled E-fields from both neurons and stimulation devices. This approach is accelerated using fast direct solvers to allow for the analysis of realistic scenarios. We present examples, indicating the ability of our solver to analyze rat L2/3 pyramidal neurons derived from the blue brain project.