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
- 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¶
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.add('compartment','Soma') x.add('compartment','Neurite') x.slice('Neurite',10) x.connect('Neurite01','Soma') x.Soma.tree_idx = 0; % mark this as the cell body
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
An example showing a multi-compartment model is found in
Maximal conductance of
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
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
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. ).
V compartment properties aren't used in this case, you can still set them to values.
If you want them to automatically update, you can set them to anonymous functions:
x.comp.A = @() x.comp.len * 2 * pi * x.comp.radius;