# Realization the model of network in Brain 2

**URL:** <https://brian.discourse.group/t/realization-the-model-of-network-in-brain-2/444>\
**Category:** Science, Projects, and Showcases\
**Created:** [5 August 2021 12:30 UTC](https://brian.discourse.group/t/realization-the-model-of-network-in-brain-2/444 "2021-08-05T12:30:57Z")\
**Posts on this page:** 2\
**Page:** 2

<div class="post-metadata">

**Author:** ![adam](https://yyz2.discourse-cdn.com/free1/user_avatar/brian.discourse.group/adam/32/197_2.png) [@adam](https://brian.discourse.group/u/adam)\
**Post date:** [17 January 2022 18:36 UTC](https://brian.discourse.group/t/realization-the-model-of-network-in-brain-2/444/21 "2022-01-17T18:36:36Z")

</div>

For a [different thread on Izhikevich neurons](https://brian.discourse.group/t/run-error-with-modeling-izhikevich-neuron/579/4) I worked through an example script with adds units to the (sort-of) dimensionless way the original equations are specified.

In doing so, I added some useful features to the script that may be worth taking a look at, including a dictionary of neuron types:

```auto
# a,b,c,d define the type of neuron
# -------------------------------------- a b c d
abcd_neurons = {'regular spiking': [0.02/ms, 0.2, -65*mV, 8*mV],
                'intrinsically bursting':[0.02/ms, 0.2, -55*mV, 4*mV],
                'chattering': [0.02/ms, 0.2, -50*mV, 2*mV],
                
                'fast spiking': [0.10/ms, 0.2, -65*mV, 2*mV],
                'low-threshold spiking': [0.02/ms, 0.25, -65*mV, 2*mV],
                
                'thalamo-cortical': [0.02/ms, 0.25, -65*mV, 0.05*mV],
                'resonator': [0.10/ms, 0.25, -65*mV, 8*mV],
                }
a,b,c,d = abcd_neurons[neuron_type]

```

of particular interest to @serge25, I was able to get intrinsic bursting:  
 ![izhi_intrinsically bursting](https://global.discourse-cdn.com/free1/uploads/brian/original/1X/554a4fb6e26d265e5cb8aeda218832439db07141.png)  
as well as thalamo-cortical post-inhibitory rebound (bursting)  
 ![izhi_thalamo-cortical_burst](https://global.discourse-cdn.com/free1/uploads/brian/original/1X/7cfd69da3f6996fd7b07c5e2744cf6b3a6d053d4.png)

[see the script here](https://brian.discourse.group/t/run-error-with-modeling-izhikevich-neuron/579/4)  
@mstimberg is this worth turning into an example for the documentation?

---

<div class="post-metadata">

**Author:** ![mstimberg](https://yyz2.discourse-cdn.com/free1/user_avatar/brian.discourse.group/mstimberg/32/11_2.png) [@mstimberg](https://brian.discourse.group/u/mstimberg)\
**Post date:** [20 January 2022 16:00 UTC](https://brian.discourse.group/t/realization-the-model-of-network-in-brain-2/444/22 "2022-01-20T16:00:33Z")

</div>

> [@adam](#):
>
> @mstimberg is this worth turning into an example for the documentation?

Would be great to have this in the example section! I wonder whether it would make sense to have the constants `a` etc. instead as constant parameters in the equations (i.e. `a : 1/second (constant)` etc.). This way, one could do a single simulation for all classes and plot them all at once. On the other hand, this would make the code a bit more complex to read, I guess. Having to choose one parametrization in the script would be perfectly fine as well, after all, e.g. [Example: Brette\_Gerstner\_2005 — Brian 2 2.5.4 documentation](https://brian2.readthedocs.io/en/stable/examples/frompapers.Brette_Gerstner_2005.html) does this as well.

[Previous page](https://brian.discourse.group/t/realization-the-model-of-network-in-brain-2/444.md?page=1)