Different initial concentration and potential in different cells
Posted: Tue Feb 01, 2022 9:32 am
Hi all,
I'm building a synapse model with a neuron and an astrocyte.
I would need to initialise the initial concentrations differently, however, not for the extracellular space.
I defined the ionic species as follow:
without the astrocyte, I would normally initialise the concentrations and the potential with:
But now that I need the initial concentration in cyt and cyt_astro different I don't know how am I supposed to define those. I tried with the following for all the ions and sections but the code just ignore these commands and sets the default.
Also, about the initialization of the voltage, how should I do it?
Thanks in advance for any help! :)
I'm building a synapse model with a neuron and an astrocyte.
I would need to initialise the initial concentrations differently, however, not for the extracellular space.
I defined the ionic species as follow:
Code: Select all
rxd_na = rxd.Species([extracellular, cyt], name = 'na', charge = 1, d = 1)
rxd_k = rxd.Species([extracellular, cyt], name = 'k', charge = 1, d = 1)
rxd_ca = rxd.Species([extracellular, cyt], name = 'ca', charge = 2, d = 1)
na_a = rxd.Species([extracellular, cyt_astro], name='na_a', charge=1, d=1)
k_a = rxd.Species([extracellular, cyt_astro], name='k_a', charge=1, d=1)
ca_a = rxd.Species([extracellular, cyt_astro], name='ca_a', charge=2, d=1)
Code: Select all
h.cao0_ca_ion = 1.8
h.nao0_na_ion = 145
h.ko0_k_ion = 3
h.cai0_ca_ion = 0.006
h.nai0_na_ion = 5
h.ki0_k_ion = 130
h.finitialize(-65*mV)
Code: Select all
rxd_ca[extracellular].initial = 1.8
...
Thanks in advance for any help! :)