Combining NEURON with a mesoscale equation
Posted: Fri Jun 16, 2017 10:44 am
Hello,
I need advise on how to combine NEURON with a mesoscopic equation (ODE).
I have built a network of simple cells (soma + dendrite with synapses) in NEURON using mainly python. The cells are connected with excitatory synapses (ExpSyn). I have chosen the cell properties in such a way that each cell generates a continuous spike train.
In order to combine the equation and NEURON's simulation, my plan is to change the mod file of the ExpSyn synapse so that at EVERY timestep the voltage change for each cell is calculated for the next time step given the voltage values from the previous timestep of the cell and all connected cells (using the mesoscopic equation). Then I would set the membrane voltage at the soma according to my calculated value. To make sure this happens at each time point, I want to set the threshold of the NetCon pointer in such a way that it is always crossed.
Because I am a student with limited amount of time, I decided to post this before I run into major implementation issues and to avoid spending a lot of time trying to implement something that couldn't work.
I would very much appreciate if someone (with more experience of and knowlegde about NEURON) could give me feedback about my plan. Do you think the ExpSyn mechanism is a good starting point to implement this kind of connection between cells? Are there any aspects that could cause problems that I haven't thought about yet?
I need advise on how to combine NEURON with a mesoscopic equation (ODE).
I have built a network of simple cells (soma + dendrite with synapses) in NEURON using mainly python. The cells are connected with excitatory synapses (ExpSyn). I have chosen the cell properties in such a way that each cell generates a continuous spike train.
In order to combine the equation and NEURON's simulation, my plan is to change the mod file of the ExpSyn synapse so that at EVERY timestep the voltage change for each cell is calculated for the next time step given the voltage values from the previous timestep of the cell and all connected cells (using the mesoscopic equation). Then I would set the membrane voltage at the soma according to my calculated value. To make sure this happens at each time point, I want to set the threshold of the NetCon pointer in such a way that it is always crossed.
Because I am a student with limited amount of time, I decided to post this before I run into major implementation issues and to avoid spending a lot of time trying to implement something that couldn't work.
I would very much appreciate if someone (with more experience of and knowlegde about NEURON) could give me feedback about my plan. Do you think the ExpSyn mechanism is a good starting point to implement this kind of connection between cells? Are there any aspects that could cause problems that I haven't thought about yet?