Create multi-compartment models
This document describes how to create and work with multi compartment models in xolotl.
Introduction and Limitations¶
xolotl is limited in how it can work with multi-compartment models:
- The only type of multi-compartment models fully supported are linear cables, with the soma at one end and the neurite coming out of it.
- You can in principle create arbitrary topologies with branching, but you are limited in the integration methods you can use (you can only use Exponential Euler, with no support for implicit solvers or multi-step solvers).
- For linear cables, xolotl can use an implicit Crank-Nicholson scheme to solve for the voltages in different parts of the cable. This is efficient and more correct, and requires you to connect parts of the cable together using the
Axial
synapse. - You cannot use the Runge-Kutta 4 solver with multi-compartment models.
With this out of the way, let's understand how multi-compartment models can be specified in xolotl.
Making a cable using the Axial
synapse type¶
The Axial
synapse is a special type of electrical synapse that forces you to think of compartments as cylinders, that have a radius and a length. When two compartments are connected using the a Axial synapse, they are assumed to be part of a cable, and xolotl will attempt to use the Crank-Nicholson scheme to solve for the voltage in all compartments on the cable implicitly.
Thus, the general recipe for creating a multi-compartment cable is:
x = xolotl;
x.add('compartment','CellBody')
x.add('compartment','Neurite','radius',.01,'len',.35*5)
x.slice('Neurite',10)
x.connect('Neurite01','CellBody')
x.CellBody.tree_idx = 0; % mark this as the cell body
Warning
tree_idx
, and the presence of Axial
synapses, is how xolotl "knows" a model is a multi-compartment model vs a set of a individual neurons that are electrically coupled. You must specify some compartment as the cell body using tree_idx
, otherwise, xolotl will not consider it a multi-compartment model.
Note that we're using the slice
function to slice up a cylinder into shorter sections, which also automatically wires them up using Axial
synapses.
An example showing a multi-compartment model is found in demo_multi_compartment.m
Maximal conductance of Axial
synapses¶
The maximal conductance of an Axial
synapse from the presynaptic compartment to the postsynaptic compartment is computed using Equation 6.30 from Dayan & Abbott (2001).
where is the axial radius and the length. is the axial resistivity in the presynaptic compartment. Subscripts "pre" and "post" refer to the presynaptic and postsynaptic compartments.
Size and shape of compartments¶
In single-compartment models, or networks of these, it is acceptable to use the A
and V
properties of compartments to specify the surface area and volume respectively.
In multi-compartment models in which the compartments are connected by Axial
synapses, the
len
and radius
properties, for the length and axial radius must be specified instead.
The shape of the compartment is cylindrical, and the surface area is computed automatically,
according to the "body" of the cylinder (i.e. ).