Numerical calculation of the light propagation in tapered optical fibers for optical neural interfaces

Abstract As implantable optical systems recently enabled new approaches to study the brain with optical radiations, tapered optical fibers emerged as promising implantable waveguides to deliver and collect light from sub-cortical structures of the mouse brain. They rely on a specific feature of multimodal fiber optics: as the waveguide narrows, the number of guided modes decreases and the radiation can gradually couple with the environment. This happens along a taper segment whose length can be tailored to match with the depth of functional structures of the mouse brain, and can extend for a few millimeters. This anatomical requirement results in optical systems with an active area very long compared to the wavelength of the light they guide and their behaviour is typically estimated by ray tracing simulations, being finite-elements methods computationally too heavy. Here we present a computational technique that exploits the beam-envelope method and the cylindrical symmetry of the fibers to provide an efficient and exact calculation of the electric field along the fibers, which may enable the design of neural interfaces optimized to meet different goals.