From NEST

[ Example networks ]

Example for the adapting exponential integrate and fire model according to Brette and Gerstner (2005) J. Neurophysiology.

The first script reproduces figure 2.C of the paper, the second script reproduces figure 3d.

Note that Brette&Gerstner give the value for b in nA. To be consistent with the other parameters in the equations, b must be converted to pA (pico Ampere).

   #! /usr/bin/env python
   
   # Test for the adapting exponential integrate and fire model according to
   # Brette and Gerstner (2005) J. Neurophysiology.
   # This script reproduces figure 2.C of the paper.
   # Note that Brette&Gerstner give the value for b in nA.
   # To be consistent with the other parameters in the equations, b must be
   # converted to pA (pico Ampere).
   
   import nest
   import nest.voltage_trace
   
   nest.ResetKernel()
   res=0.1
   nest.SetStatus([0],[{"resolution": res}])
   neuron=nest.Create("aeif_cond_alpha")
   nest.SetStatus(neuron,[{"V_peak": 0.0, "a": 4.0, "b":80.5}])
   dc=nest.Create("dc_generator",2)
   
   nest.SetStatus(dc,[{"amplitude":500.0, 
                       "start":0.0, 
                       "stop":200.0},
                      {"amplitude":800.0,
                       "start":500.0,
                       "stop":1000.0}])
   
   nest.ConvergentConnect(dc,neuron)
   
   voltmeter= nest.Create("voltmeter")
   nest.SetStatus(voltmeter,[{"to_file": True, "withtime": True}])
   nest.Connect(voltmeter,neuron)
   
   nest.Simulate(1000.0)
   
   nest.voltage_trace.from_device(voltmeter)

---

   #! /usr/bin/env python
   
   # Test for the adapting exponential integrate and fire model according to
   # Brette and Gerstner (2005) J. Neurophysiology.
   # This script reproduces figure 3.d of the paper.
   # Note that Brette&Gerstner give the value for b in nA.
   # To be consistent with the other parameters in the equations, b must be
   # converted to pA (pico Ampere).
   
   import nest
   import nest.voltage_trace
   
   nest.ResetKernel()
   res=0.1
   nest.SetStatus([0],[{"resolution": res}])
   neuron=nest.Create("aeif_cond_alpha")
   nest.SetStatus(neuron,[{"V_peak": 0.0, "E_L":-60.0, "a":80.0, "b":80.5, "tau_w": 720.0}])
   dc=nest.Create("dc_generator")
   
   nest.SetStatus(dc,[{"amplitude":-800.0, "start":0.0, "stop":400.0}])
   
   nest.ConvergentConnect(dc,neuron)
   
   voltmeter= nest.Create("voltmeter")
   nest.SetStatus(voltmeter,[{"to_file": True, "withtime": True}])
   nest.Connect(voltmeter,neuron)
   
   nest.Simulate(1000.0)
   
   nest.voltage_trace.from_device(voltmeter)