source: project/release/4/9ML-toolkit/trunk/examples/Brunel00/Brunel_network_alpha.xml @ 30932

<?xml version='1.0' encoding='UTF-8'?>
<NineML xmlns="http://nineml.incf.org/9ML/0.3"
    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
    xsi:schemaLocation="http://nineml.incf.org/9ML/0.3 ../NineML_v0.3.xsd"
    name="Brunel (2000) Model A Network with alpha synapses">
  
  <Property><!-- Scales the size of the network (total 5*order neurons) -->
    <label>order</label>
    <value>2500</value>
  </Property>

  <Property><!-- Number of excitatory neurons -->
    <label>Ne</label>
    <value><MathInline>4 * order</MathInline></value>
  </Property>

  <Property><!-- Number of inhibitory neurons -->
    <label>Ni</label>
    <value><MathInline>1 * order</MathInline></value>
  </Property>

  <Property><!-- Connectivity probability -->
    <label>epsilon</label>
    <value>0.1</value>
  </Property>

  <Property><!-- Total number of excitatory synapses -->
    <label>Ce</label>
    <value><MathInline>epsilon * Ne</MathInline></value>
  </Property>

  <Property><!-- Total number of inhibitory synapses -->
    <label>Ci</label>
    <value><MathInline>epsilon * Ni</MathInline></value>
  </Property>

  <Property><!-- Total number of external synapses -->
    <label>Cext</label>
    <value><MathInline>Ce</MathInline></value>
  </Property>

  <Property><!-- Global delay for all neurons in the group -->
    <label>delay</label>
    <value>1.5</value>
    <unit>ms</unit>
  </Property>
  
  <Property> <!-- Synaptic weight -->
    <label>J</label>
    <value>0.1</value> 
    <unit>mV</unit> 
  </Property>

  <Property> 
    <label>J_max_unit</label>
    <value><MathInline>1.209</MathInline></value>
  </Property>

  <Property> 
    <label>g</label> <!-- Relative strength of inhibitory synapses -->
    <value>3.0</value> 
  </Property>

  <Property> <!-- nu_ext / nu_thr -->
    <label>eta</label>
    <value>2.0</value> 
  </Property>

  <Property> <!-- Excitatory weights -->
    <label>Je</label>
    <value><MathInline>J / J_max_unit</MathInline></value> 
  </Property>

  <Property> <!-- Inhibitory weights -->
   <label>Ji</label>
    <value><MathInline>-g * Je</MathInline></value> 
  </Property>

  <Property> <!-- External weights -->
    <label>Jext</label>
    <value><MathInline>Je</MathInline></value> 
  </Property>

  <Property> <!-- Firing threshold -->
    <label>theta</label>
    <value>20.0</value>
    <unit>mV</unit>
  </Property>
  
  <Property><!-- Membrane time constant -->
    <label>tau</label>
    <value>20.0</value>
    <unit>ms</unit>
  </Property>

  <Property><!-- Synapse time constant -->
    <label>tau_syn</label>
    <value>0.5</value>
  </Property>

  <Property> <!-- Threshold rate -->
    <label>nu_thresh</label>
    <value><MathInline>theta / (Je * Ce * tau * exp (1.0) * tau_syn)</MathInline></value> 
  </Property>
  
  <Property><!-- External rate per synapse --> 
    <label>nu_ext</label>
    <value><MathInline>eta * nu_thresh</MathInline></value>
  </Property>

  <Property> <!-- Mean input spiking rate -->
    <label>input_rate</label>
    <value><MathInline>1000.0 * nu_ext</MathInline></value>
    <unit>Hz</unit>
  </Property>
  
  <Component name="E">

    <Definition>
      <link> BrunelIaF.xml </link>
    </Definition>

    <properties>

      <Property>
        <label>tau</label>
        <value>20.0</value>
      </Property>
      
      <Property>
        <label>theta</label>
        <value><MathInline>theta</MathInline></value>
      </Property>

      <Property>
        <label>tau_rp</label>
        <value>2.0</value>
      </Property>
      
      <Property>
        <label>Vreset</label>
        <value>10.0</value>
      </Property>
      
      <Property>
        <label>R</label>
        <value>1.0</value>
      </Property>
      
      <Initial>
        <label>Isyn</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>V</label>
        <value><MathInline>randomUniform ()</MathInline></value>
      </Initial>

      <Initial>
        <label>t_rpend</label>
        <value>0.0</value>
      </Initial>
      
      <Initial>
        <label>spikeOutput</label>
        <value>false</value>
      </Initial>
      
      <Initial>
        <label>refractoryEnd</label>
        <value>false</value>
      </Initial>

    </properties>

  </Component>

  <Component name="I">

    <Definition>
      <link> BrunelIaF.xml </link>
    </Definition>

    <properties>

      <Property>
        <label>tau</label>
        <value>20.0</value>
      </Property>
      
      <Property>
        <label>theta</label>
        <value><MathInline>theta</MathInline></value>
      </Property>

      <Property>
        <label>tau_rp</label>
        <value>2.0</value>
      </Property>
      
      <Property>
        <label>Vreset</label>
        <value>10.0</value>
      </Property>

      <Property>
        <label>R</label>
        <value>1.0</value>
      </Property>

      <Initial>
        <label>Isyn</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>V</label>
        <value><MathInline>randomUniform ()</MathInline></value>
      </Initial>
      
      <Initial>
        <label>t_rpend</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>spikeOutput</label>
        <value>false</value>
      </Initial>
      
      <Initial>
        <label>refractoryEnd</label>
        <value>false</value>
      </Initial>
      
      
    </properties>

  </Component>


  <Component name="Ext">
    
    <Definition>
      <link> http://www.NineML.org/stim/Poisson.9ml </link>
    </Definition>
    
    <properties>

      <Property> 
        <label>rate</label>
        <value><MathInline>input_rate</MathInline></value>
      </Property>
      
    </properties>
    
  </Component>

  <Component name="ISyn">
    
    <Definition>
      <link> AlphaPSR.xml </link>
    </Definition>
    
    <properties>

      <Property> 
        <label>tau_syn</label>
        <value><MathInline>tau_syn</MathInline></value>
      </Property>

      <Initial>
        <label>A</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>B</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>spike</label>
        <value>false</value>
      </Initial>

      <Initial>
        <label>tspike</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>Ispike</label>
        <value>0.0</value>
      </Initial>
      
      <Initial>
        <label>Isyn</label>
        <value>0.0</value>
      </Initial>
      
    </properties>
    
  </Component>

  <Component name="ESyn">
    
    <Definition>
      <link> AlphaPSR.xml </link>
    </Definition>
    
    <properties>

      <Property> 
        <label>tau_syn</label>
        <value><MathInline>tau_syn</MathInline></value>
      </Property>

      <Initial>
        <label>A</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>B</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>spike</label>
        <value>false</value>
      </Initial>

      <Initial>
        <label>tspike</label>
        <value>0.0</value>
      </Initial>

      <Initial>
        <label>Ispike</label>
        <value>0.0</value>
      </Initial>
      
      <Initial>
        <label>Isyn</label>
        <value>0.0</value>
      </Initial>
      
    </properties>
    
  </Component>


  <Component name="ExternalPlasticity">
    
    <Definition>
      <link> ConstantPlasticity.xml </link>
    </Definition>
    
    <properties>
      
      <Initial>
        <label>weight</label>
        <value><MathInline>Jext</MathInline></value>
      </Initial>
      
    </properties>
    
  </Component>

  <Component name="ExcitatoryPlasticity">
    
    <Definition>
      <link> ConstantPlasticity.xml </link>
    </Definition>
    
    <properties>
      
      <Initial>
        <label>weight</label>
        <value><MathInline>Je</MathInline></value>
      </Initial>
      
    </properties>
    
  </Component>

  <Component name="InhibitoryPlasticity">
    
    <Definition>
      <link> ConstantPlasticity.xml </link>
    </Definition>
    
    <properties>
      
      <Initial>
        <label>weight</label>
        <value><MathInline>Ji</MathInline></value>
      </Initial>
      
    </properties>
    
  </Component>


  <Component name="RandomUniform">
    
    <Definition>
      <link> RandomUniformConnection.xml </link>
    </Definition>
    
    <properties>
      <Property>
        <label>epsilon</label>
        <value><MathInline>epsilon</MathInline></value>
      </Property>
    </properties>

  </Component>

  <Group name="BrunelCaseAalpha">

    <Population name="Exc"> <!-- Excitatory population definition -->  
      <prototype><reference>E</reference></prototype>
      <number><MathInline>Ne</MathInline></number>
    </Population>

    <Population name="Inh"> <!-- Inhibitory population definition -->  
      <prototype><reference>I</reference></prototype>
      <number><MathInline>Ni</MathInline></number>
    </Population>

    <Population name="Ext"> <!-- External sources definition -->  
      <prototype><reference>Ext</reference></prototype>
      <number><MathInline>Cext</MathInline></number>
    </Population>

    <Set name="All neurons"> <!-- The set of all neurons -->  
      <union>
        <left><singleton>Exc</singleton></left> 
        <right><singleton>Inh</singleton></right>
      </union>
    </Set>

    <!-- Projection from external sources to all neurons -->  
    <Projection name="External" 
                input_event_port="spike" 
                input_analog_port="Isyn"
                output_event_port="spikeOutput"
                input_response_port="Ispike"> 
      <source><reference>Ext</reference></source>
      <target><reference>All neurons</reference></target>
      <rule operator="all-to-all"/>
      <response><reference>ESyn</reference></response>
      <plasticity><reference>ExternalPlasticity</reference></plasticity>
    </Projection>

    <!-- Projection from excitatory neurons to all neurons -->  
    <Projection name="Excitation" 
                input_event_port="spike" 
                input_analog_port="Isyn"
                output_event_port="spikeOutput"
                input_response_port="Ispike"> 
      <source><reference>Exc</reference></source>
      <target><reference>All neurons</reference></target>
      <rule operator="for-each">
        <reference>RandomUniform</reference>
        <state>connection</state>
      </rule>
      <response><reference>ESyn</reference></response>
      <plasticity><reference>ExcitatoryPlasticity</reference></plasticity>
    </Projection>

    <!-- Projection from inhibitory neurons to all neurons -->  
    <Projection name="Inhibition"
                input_event_port="spike" 
                input_analog_port="Isyn"
                output_event_port="spikeOutput"
                input_response_port="Ispike"> 
      <source><reference>Inh</reference></source>
      <target><reference>All neurons</reference></target>
      <rule operator="for-each">
        <reference>RandomUniform</reference>
        <state>connection</state>
      </rule>
      <response><reference>ISyn</reference></response>
      <plasticity><reference>InhibitoryPlasticity</reference></plasticity>
    </Projection>

    <!-- Record spikes from the given set -->
    <SpikeRecording set="All neurons"/>

  </Group>

</NineML>