www.neuron.yale.edu

My colleague and I have dissected the method and code of the paper by Pashut et al. and identified the majority of problems related to the coupling of the electric field to the neuron model. The details are posted in an online comment to the publication: https://journals.plos.org/ploscompbiol/articl...

Sorry for not being clear. Since I was using the default HH channels in NEURON, I indeed meant the HH.mod file distributed with NEURON, which I located under the folder: \src\nrnoc. Thanks for clarifying the vtrap function and continuity requirement. I thought that the 1e-6 band was narrow enough th...

Ted, thank you for the suggestion!! I have doubted the HH.mod file as well, but wasn't very confident to change it. After all, this is the first time I worked with hoc code. I have adjusted the .mod file and the problem has been resolved even without using LinearMechanism. Thresholds could be obtain...

I implemented the model with LinearMechanism, and unfortunately the S-D curve was pretty much identical compared to before. Therefore, I don't think that the very low axial resistance was the problem, but something else contributed to the behavior with short pulse duration and small time steps.

I have actually tried this "spoke-hub" topology as well. The results were almost identical as the other topology, with the same issue for short pulses. Although, I have not experimented so much with changing Ra (set only to 1e-3) or the section's size (set to cylinders of equal length and ...

Hello Ted, thanks for your clarification and suggestion!

I'll explore both the Linear Circuit Builder and Linear Mechanism and see whether it works.

Hello all! My questions are: 1. How does NEURON handle axial/intracellular resistance in its internal calculation, and how would setting intracellular resistivity/resistance to very small values (Ra<<default) affect the accuracy and stability of the simulation? Could very short pulse duration (<10 u...