HHNeuron.h

Go to the documentation of this file.

#ifndef _HHNEURON_H_
#define _HHNEURON_H_

#include "ODESystemBasedSpikingNeuron.h"



/* with an additional M channel and OU Noise
 * 
 */
class HHNeuron: public ODESystemBasedSpikingNeuron
{
    SIMOBJECT(HHNeuron, AdvancePhase::One)

public:
    double Amembrane; 
    double TADJM;

    double gbar_Na;
    double E_Na;
    double gbar_K;
    double E_K;
    double gbar_M;
    double E_M;

    HHNeuron(double Rm=100e6, 
               double Cm=3e-10, 
               double Vresting=-60e-3, 
               double Vthresh=0.03,  
               double Vinit=-60e-3, 
               double Trefract=5e-3, 
               double Inoise=0.0, 
               double Iinject=0.0) 
       : ODESystemBasedSpikingNeuron(5)
    {
        Amembrane=34636e-12;
        TADJM=2.9529;

        this->Rm=Rm;
        this->Cm=Cm;
        this->Vresting=Vresting;
        this->Vthresh=Vthresh;
        this->Vinit=Vinit;
        this->Trefract=Trefract;
        this->Inoise=Inoise;
        this->Iinject=Iinject;

//        gbar_Na = 500.0 * Amembrane;
        gbar_Na = 4.3 * 120 * Amembrane;
        E_Na    = 50e-3;

        gbar_K  = 100.0 * Amembrane;
        E_K     = -90e-3;

        gbar_M  = 5.0 * Amembrane;
        E_M     = -90e-3;
    };

    double alpha_m(double V) { return 1000.0*0.32*(13.0 - V)/(exp((13.0 - V)/4.0) - 1.0); }
    double beta_m(double V)  { return 1000.0*0.28*(V - 40.0) / (exp((V - 40.0)/5.0) - 1.0); }
    double infty_m(double V) { return alpha_m(V)/(beta_m(V)+alpha_m(V)); }

    double alpha_h(double V) { return 1000.0*0.128*exp((17.0 - V - 10.0)/18.0); }
    double beta_h(double V)  { return 1000.0*4/(1.0 + exp((40.0 - V - 10.0)/5.0)); }
    double infty_h(double V) { return alpha_h(V)/(beta_h(V)+alpha_h(V)); }

    double alpha_n(double V) { return 1000.0*0.032*(15.0 - V)/(exp((15.0 - V)/5.0) - 1.0); }
    double beta_n(double V)  { return 1000.0*0.5*exp((10.0 - V)/40.0); }
    double infty_n(double V) { return alpha_n(V)/(beta_n(V)+alpha_n(V)); }

    double alpha_nM(double V) { return TADJM*1000.0*0.0001*(V+30.0)/(1.0 - exp(-(V+30.0)/9.0)); }       
    double beta_nM(double V)  { return TADJM*1000.0*-0.0001*(V+30.0)/(1.0 - exp((V+30.0)/9.0)); }
    double infty_nM(double V) { return alpha_nM(V)/(beta_nM(V)+alpha_nM(V)); }

    virtual void derivatives(double t, const double y[], double f[]);

    virtual int reset(double dt); 
};


class HHINeuron : public ODESystemBasedSpikingNeuron 
{
    SIMOBJECT(HHINeuron, AdvancePhase::One)

public:
    double Amembrane;

    double gbar_Na;
    double E_Na;

    double gbar_K;
    double E_K;

    HHINeuron(double Rm=100e6, 
                double Cm=3e-10, 
                double Vresting=-60e-3, 
                double Vthresh=0.03,  
                double Vinit=-60e-3, 
                double Trefract=5e-3, 
                double Inoise=0.0, 
                double Iinject=0.0) 
        : ODESystemBasedSpikingNeuron(4)
    {
        Amembrane=34636e-12;

        this->Rm=Rm;
        this->Cm=Cm;
        this->Vresting=Vresting;
        this->Vthresh=Vthresh;
        this->Vinit=Vinit;
        this->Trefract=Trefract;
        this->Inoise=Inoise;
        this->Iinject=Iinject;

//        gbar_Na = 500.0 * Amembrane;
        gbar_Na = 4.3 * 120 * Amembrane;
        E_Na    = 50e-3;

        gbar_K  = 100.0 * Amembrane;
        E_K     = -90e-3;
    };

    double alpha_m(double V) { return 1000.0*0.32*(13.0 - V)/(exp((13 - V)/4.0) - 1.0); }
    double beta_m(double V)  { return 1000.0*0.28*(V - 40.0) / (exp((V - 40)/5.0) - 1.0); }
    double infty_m(double V) { return alpha_m(V)/(beta_m(V)+alpha_m(V)); }

    double alpha_h(double V) { return 1000.0*0.128*exp((17.0 - V - 10.0)/18.0); }
    double beta_h(double V)  { return 1000.0*4.0/(1.0 + exp((40.0 - V - 10.0)/5.0)); }
    double infty_h(double V) { return alpha_h(V)/(beta_h(V)+alpha_h(V)); }

    double alpha_n(double V) { return 1000.0*0.032*(15.0 - V)/(exp((15.0 - V)/5.0) - 1.0); }
    double beta_n(double V)  { return 1000.0*0.5*exp((10.0 - V)/40.0); }
    double infty_n(double V) { return alpha_n(V)/(beta_n(V)+alpha_n(V)); }

    virtual void derivatives(double t, const double y[], double f[]);

    virtual int reset(double dt); 
};

#endif // _HHNEURON_H_

---