starting from Neuron 7.0, impedance affects simulation!
Posted: Tue Jan 26, 2010 10:12 am
Hi,
I want to monitor the impedance of the model neuron while simulation is running, let's say every 50 ms. So I made a script that advances 50 ms, pauses the simulations and recalculates an impedance object that is inserted. However, starting from Neuron 7.0 I found out that I can't do this any longer, as something (I don't know what yet) changes when the impedance is computed and the rest of the simulation is completely different!
I have a model of slow oscillation bursting. It is only one section, one segment, and has an impedance object inserted. If I do:
The last 2000 ms of the simulation will have a completely different pattern of oscillation. This is regardless of the presence of spiking (i.e. fast sodium and potassium) conductances. Looking at the oscillations and how they changed it is as if the leak (passive) conductance was decreased. However, it doesn't happen when the parameters are changed such that oscillations don't occur (i.e., the model is stable).
But the most important aspect of this problem is that it didn't happen with Neuron 6.2. What has changed since with the impedance class? I can still keep working with neuron 6.2 for a while, but I think this is something that should be addressed in some way... or not?
Best Regards
EDIT: You can reproduce it in a section with HH mechanism, and a IClamp that makes it fire continuously.
I want to monitor the impedance of the model neuron while simulation is running, let's say every 50 ms. So I made a script that advances 50 ms, pauses the simulations and recalculates an impedance object that is inserted. However, starting from Neuron 7.0 I found out that I can't do this any longer, as something (I don't know what yet) changes when the impedance is computed and the rest of the simulation is completely different!
I have a model of slow oscillation bursting. It is only one section, one segment, and has an impedance object inserted. If I do:
Code: Select all
** run for 2000 ms
imp.compute(0)
** run for 2000 ms
But the most important aspect of this problem is that it didn't happen with Neuron 6.2. What has changed since with the impedance class? I can still keep working with neuron 6.2 for a while, but I think this is something that should be addressed in some way... or not?
Best Regards
EDIT: You can reproduce it in a section with HH mechanism, and a IClamp that makes it fire continuously.