Alpha rhythm in vitro visual cortex (Traub et al 2020)

 Download zip file 
Help downloading and running models
Accession:263703
The paper describes an experimental model of the alpha rhythm generated by layer 4 pyramidal neurons in a visual cortex slice. The simulation model is derived from that of Traub et al. (2005) J Neurophysiol, developed for thalamocortical oscillations.
Reference:
1 . Traub RD, Hawkins K, Adams NE, Hall SP, Simon A, Whittington MA (2020) Layer 4 pyramidal neuron dendritic bursting underlies a post-stimulus visual cortical alpha rhythm Nature Communications Biology, in press
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Thalamus; Neocortex;
Cell Type(s): Thalamus geniculate nucleus/lateral principal GLU cell; Thalamus reticular nucleus GABA cell; Neocortex U1 L6 pyramidal corticalthalamic GLU cell; Neocortex U1 L2/6 pyramidal intratelencephalic GLU cell; Neocortex layer 4 pyramidal cell; Neocortex fast spiking (FS) interneuron; Neocortex spiking regular (RS) neuron; Neocortex spiking low threshold (LTS) neuron; Neocortex spiny stellate cell;
Channel(s): I Na,p; I Na,t; I L high threshold; I A; I K; I M; I h; I K,Ca; I Calcium; I A, slow;
Gap Junctions: Gap junctions;
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: FORTRAN;
Model Concept(s): Brain Rhythms; Dendritic Action Potentials;
Implementer(s): Traub, Roger D [rtraub at us.ibm.com];
Search NeuronDB for information about:  Thalamus geniculate nucleus/lateral principal GLU cell; Thalamus reticular nucleus GABA cell; Neocortex U1 L2/6 pyramidal intratelencephalic GLU cell; Neocortex U1 L6 pyramidal corticalthalamic GLU cell; I Na,p; I Na,t; I L high threshold; I A; I K; I M; I h; I K,Ca; I Calcium; I A, slow;
/
alpha_rhythm_code_RT
readme.txt
alphaY33.pdf
alphaY67.f
alphaY67.pdf
dexptablebig_setup.f *
dexptablesmall_setup.f *
durand.f *
fnmda.f *
gettime.c *
groucho_gapbld.f *
groucho_gapbld_mix.f *
integrate_deepaxaxx.f *
integrate_deepbaskx.f *
integrate_deepLTSx.f *
integrate_deepng.f *
integrate_nontuftRSXXB.f
integrate_nrtxB.f *
integrate_spinstelldiegoxB.f *
integrate_supaxaxx.f *
integrate_supbaskx.f *
integrate_supLTSX.f *
integrate_supng.f *
integrate_suppyrFRBxPB.f *
integrate_suppyrRSXPB.f
integrate_tcrxB.f *
integrate_tuftIBVx3B.f
integrate_tuftRSXXB.f *
makefile
otis.f *
otis_table_setup.f *
synaptic_compmap_construct.f *
synaptic_map_construct.f *
                            
      FUNCTION ran1(idum)
      INTEGER idum,IA,IM,IQ,IR,NTAB,NDIV
      DOUBLE PRECISION ran1,AM,EPS,RNMX
      PARAMETER (IA=16807,IM=2147483647,AM=1.d0/IM,IQ=127773,IR=2836,
     *NTAB=32,NDIV=1+(IM-1)/NTAB,EPS=3.d-16,RNMX=1.d0-EPS)
      INTEGER j,k,iv(NTAB),iy
      SAVE iv,iy
      DATA iv /NTAB*0/, iy /0/
      if (idum.le.0.or.iy.eq.0) then
        idum=max(-idum,1)
        do 11 j=NTAB+8,1,-1
          k=idum/IQ
          idum=IA*(idum-k*IQ)-IR*k
          if (idum.lt.0) idum=idum+IM
          if (j.le.NTAB) iv(j)=idum
11      continue
        iy=iv(1)
      endif
      k=idum/IQ
      idum=IA*(idum-k*IQ)-IR*k
      if (idum.lt.0) idum=idum+IM
      j=1+iy/NDIV
      iy=iv(j)
      iv(j)=idum
      ran1=min(AM*iy,RNMX)
      return
      END

      subroutine durand(seed, npts, x)
      implicit none
      integer npts, i, idum
      real*8  seed, ran1,  x(npts)
      if (seed .lt. 0.0d0  .or. seed .gt. 2147483648.0) then
         write(6,*) 'seed must be a positive integer < 2.1474*10**9'
         seed = 135791113.0
      end if
      idum = -int(seed)
      do i = 1, npts
           x(i) = ran1(idum)
      end do
      seed = dble(idum)
      end