<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Neurorobot.</journal-id>
<journal-title>Frontiers in Neurorobotics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Neurorobot.</abbrev-journal-title>
<issn pub-type="epub">1662-5218</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnbot.2017.00061</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Cognitive Mapping Based on Conjunctive Representations of Space and Movement</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Zeng</surname> <given-names>Taiping</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/320504"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Si</surname> <given-names>Bailu</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="cor1">&#x0002A;</xref>
<uri xlink:href="http://frontiersin.org/people/u/218657"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>State Key Laboratory of Robotics, Shenyang Institute of Automation, Chinese Academy of Sciences</institution>, <addr-line>Shenyang</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>University of Chinese Academy of Sciences</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Shuai Li, Hong Kong Polytechnic University, Hong Kong</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Yang Zhan, Shenzhen Institutes of Advanced Technology (CAS), China; Lin Xiao, Jishou University, China</p></fn>
<corresp content-type="corresp" id="cor1">&#x0002A;Correspondence: Bailu Si, <email>sibailu&#x00040;sia.ac.cn</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>22</day>
<month>11</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection">
<year>2017</year>
</pub-date>
<volume>11</volume>
<elocation-id>61</elocation-id>
<history>
<date date-type="received">
<day>05</day>
<month>01</month>
<year>2017</year>
</date>
<date date-type="accepted">
<day>19</day>
<month>10</month>
<year>2017</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2017 Zeng and Si.</copyright-statement>
<copyright-year>2017</copyright-year>
<copyright-holder>Zeng and Si</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>It is a challenge to build robust simultaneous localization and mapping (SLAM) system in dynamical large-scale environments. Inspired by recent findings in the entorhinal&#x02013;hippocampal neuronal circuits, we propose a cognitive mapping model that includes continuous attractor networks of head-direction cells and conjunctive grid cells to integrate velocity information by conjunctive encodings of space and movement. Visual inputs from the local view cells in the model provide feedback cues to correct drifting errors of the attractors caused by the noisy velocity inputs. We demonstrate the mapping performance of the proposed cognitive mapping model on an open-source dataset of 66&#x02009;km car journey in a 3&#x02009;km&#x02009;&#x000D7;&#x02009;1.6&#x02009;km urban area. Experimental results show that the proposed model is robust in building a coherent semi-metric topological map of the entire urban area using a monocular camera, even though the image inputs contain various changes caused by different light conditions and terrains. The results in this study could inspire both neuroscience and robotic research to better understand the neural computational mechanisms of spatial cognition and to build robust robotic navigation systems in large-scale environments.</p>
</abstract>
<kwd-group>
<kwd>medial entorhinal cortex</kwd>
<kwd>head-direction cells</kwd>
<kwd>conjunctive grid cells</kwd>
<kwd>cognitive map</kwd>
<kwd>attractor dynamics</kwd>
<kwd>path integration</kwd>
<kwd>monocular SLAM</kwd>
<kwd>bio-inspired robots</kwd>
</kwd-group>
<counts>
<fig-count count="10"/>
<table-count count="1"/>
<equation-count count="20"/>
<ref-count count="58"/>
<page-count count="16"/>
<word-count count="10617"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="introduction">
<label>1</label> <title>Introduction</title>
<p>One of the fundamental challenges to integrate robots into our daily life is to endow with the ability of spatial cognition. Over the past two decades, many efforts have been exerted to develop truly autonomous robot navigation systems that are able to simultaneously localize, map, and navigate without human interventions in domestic home and workplace. This line of research has been recognized as the problem of Simultaneous Localization and Mapping (SLAM) (Smith and Cheeseman, <xref ref-type="bibr" rid="B40">1986</xref>; Bailey and Durrant-Whyte, <xref ref-type="bibr" rid="B1">2006</xref>; Durrant-Whyte and Bailey, <xref ref-type="bibr" rid="B11">2006</xref>; Thrun and Leonard, <xref ref-type="bibr" rid="B48">2008</xref>) and constitutes a significant area of robotic research. The approaches to solving SLAM problem fall into two classes. The classical approach is the probabilistic one that uses the Extended Kalman Filters, the Particle Filters, etc., as the basis for robot SLAM (Thrun et al., <xref ref-type="bibr" rid="B47">2005</xref>). Numerous variants of filtering algorithms have been proposed to estimate the position and the heading direction of the robot in large environments. However, most of these algorithms require costly sensors, significant computational resources, and the assumption of static environment. Although filtering algorithms provide beautiful mathematical frameworks for SLAM, they are in lack of robustness, since bad data association (i.e., the matching between perception and internal representation) can often lead them into irreversible divergence, impairing their application to dynamic and large-scale environments. For a good review, one can refer to Thrun and Leonard (<xref ref-type="bibr" rid="B48">2008</xref>); Dudek and Jenkin (<xref ref-type="bibr" rid="B10">2010</xref>).</p>
<p>The other approach of SLAM research is sparkled by the inspirations from neuroscience and aims to build neural network models that work in analogy to the neural circuits of spatial cognition (Zeno et al., <xref ref-type="bibr" rid="B56">2016</xref>). Animals show amazing navigation ability (Klaassen et al., <xref ref-type="bibr" rid="B18">2011</xref>; Tsoar et al., <xref ref-type="bibr" rid="B51">2011</xref>), which has long been hypothesized to result from internal map-like representation in neural space, i.e., cognitive map of the environment (Tolman, <xref ref-type="bibr" rid="B49">1948</xref>). The idea of cognitive map finds embodiment, thanks to the discovery of place cells (PCs) in the hippocampus of mammalian brains (O&#x02019;Keefe and Dostrovsky, <xref ref-type="bibr" rid="B31">1971</xref>; O&#x02019;Keefe and Conway, <xref ref-type="bibr" rid="B30">1978</xref>). A random set of place cells are recruited to encode positions in an environment and fire selectively in some restricted portions of the environment. The firing locations of place cells are randomly redistributed across environments and do not conserve metric relation. Place cells are considered to form topological maps of environments. Head direction cells (HD cells), first found in the post-subiculum (Taube et al., <xref ref-type="bibr" rid="B46">1990</xref>), fire when the animal&#x02019;s head points to a specific direction relative to an allocentric reference frame. Head directions cells exist in various brain regions (Taube, <xref ref-type="bibr" rid="B45">2007</xref>) and provide directional information like a compass. The grid cells discovered in the dorsolateral band of the medial entorhinal cortex (dMEC) fire in multiple locations that collectively define triangular grids spanning the whole environment explored by the animal (Hafting et al., <xref ref-type="bibr" rid="B16">2005</xref>). The population codes of grid cells can accurately encode animal&#x02019;s spatial positions in an environment (Mathis et al., <xref ref-type="bibr" rid="B25">2012</xref>). The relative spatial phases of nearby grid cells in dMEC are kept fixed across environments, up to a coherent rotation and shift (Fyhn et al., <xref ref-type="bibr" rid="B14">2007</xref>). Grid cells are thought to be an important ingredient of the inner GPS in mammalian brains (Knierim, <xref ref-type="bibr" rid="B19">2015</xref>). Grid cells show a rich variety of response properties. Pure positional grid cells, abundant in layer II of MEC, do not show strong modulation by movements. Conjunctive grid cells, residing mainly in layer III, V, and VI of MEC, are selective also to running directions (Sargolini et al., <xref ref-type="bibr" rid="B34">2006</xref>). More recently, speed cells found in the MEC are characterized by a context-invariant positive, linear response to running speed, which is a key component of the dynamic representation of the animal&#x02019;s self location in MEC (Kropff et al., <xref ref-type="bibr" rid="B22">2015</xref>).</p>
<p>Sitting between sensory cortex and hippocampus, MEC is a converging area where grid cells, HD cells, and speed cells work together to integrate movement and sensory information. Computational models of grid cells reveal possible information processing mechanisms of MEC in forming spatial representations and provide building blocks for cognitive mapping systems (Zilli, <xref ref-type="bibr" rid="B58">2012</xref>). The large majority of grid cells can be divided into three classes: oscillatory interference (OI) models, continuous attractor network models (CAN) (Fuhs and Touretzky, <xref ref-type="bibr" rid="B13">2006</xref>; McNaughton et al., <xref ref-type="bibr" rid="B26">2006</xref>; Burak and Fiete, <xref ref-type="bibr" rid="B3">2009</xref>; Knierim and Zhang, <xref ref-type="bibr" rid="B20">2012</xref>; Si et al., <xref ref-type="bibr" rid="B37">2014</xref>), and self-organization (SO) models (Kropff and Treves, <xref ref-type="bibr" rid="B23">2008</xref>; Si and Treves, <xref ref-type="bibr" rid="B39">2013</xref>; Stella et al., <xref ref-type="bibr" rid="B42">2013</xref>; Stella and Treves, <xref ref-type="bibr" rid="B43">2015</xref>). OI models propose that the hexagonal grid pattern arises from the interference between velocity modulated oscillators and a theta oscillator (Burgess et al., <xref ref-type="bibr" rid="B4">2007</xref>; Welday et al., <xref ref-type="bibr" rid="B52">2011</xref>). The OI models are susceptible to drifting error and lack of robustness (Remme et al., <xref ref-type="bibr" rid="B32">2010</xref>) and may not be robust enough for robot navigation systems. The CAN models use structured connections between cells to generate hexagonal grid firing patterns in the population level (Fuhs and Touretzky, <xref ref-type="bibr" rid="B13">2006</xref>; McNaughton et al., <xref ref-type="bibr" rid="B26">2006</xref>; Burak and Fiete, <xref ref-type="bibr" rid="B3">2009</xref>). The grid firing pattern is a stable network state arising from the collective behavior of the neural population and lies in a continuous manifold of steady states. If coupled by velocity inputs, a continuous attractor network model of grid cells is able to propagate its network state along the manifold, perform accurate path integration, and, therefore produce single cell hexagonal grid firing patterns in the environment (Burak and Fiete, <xref ref-type="bibr" rid="B3">2009</xref>). The conjunctive grid cells observed in layer III and deeper layers of MEC are modeled by a CAN that encodes position and velocity explicitly (Si et al., <xref ref-type="bibr" rid="B37">2014</xref>). The conjunctive position-by-velocity representation produces intrinsically moving patterns that can be selected by velocity inputs. The movement of the pattern is not dependent on the amplitude of the velocity signal or on the firing properties of the cells and, therefore, is able to achieve robust path integration even if the connections are contaminated by noise. SO models of grid cells assume that grid cells receive broadly tuned spatial inputs, e.g., feedforward inputs from visual areas or the hippocampus. Single cell firing rate adaptation could result in triangular firing patterns (Kropff and Treves, <xref ref-type="bibr" rid="B23">2008</xref>). Since the feedforward inputs function as error correcting signals, the anchoring of grid maps in the environment is strongly influenced by the inputs experienced by the animal (Si et al., <xref ref-type="bibr" rid="B36">2012</xref>; Stensola et al., <xref ref-type="bibr" rid="B44">2015</xref>).</p>
<p>Although numerous computational models of grid cells have been proposed to account for the mechanisms of integrating movement and sensory information into spatial representations, so far, only a few robot navigation systems utilize computational mechanisms of spatial cognition. In the model of RatSLAM, pose cells were developed to encode positions in large-scale environment for long-term robotic navigation tasks (Milford and Wyeth, <xref ref-type="bibr" rid="B28">2008</xref>). Pose cell is a kind of abstract cell, which updates its activity by displacing a copy of the activity packet, rather than performing path integration according to the dynamics of the network. In Jauffret et al. (<xref ref-type="bibr" rid="B17">2015</xref>), a modulo projection of the path integration is used to generate the grid cell firing patterns. However, grid cells aid neither localization nor mapping in the mobile robot system. There is a sparse number of robot navigation systems, which use grid cell networks to do path integration and represent the pose of a robot (Sheynikhovich et al., <xref ref-type="bibr" rid="B35">2009</xref>; Yuan et al., <xref ref-type="bibr" rid="B55">2015</xref>; Mulas et al., <xref ref-type="bibr" rid="B29">2016</xref>). Sheynikhovich et al. (<xref ref-type="bibr" rid="B35">2009</xref>) built a biological plausible model including visual, medial entorhinal cortex, hippocampus, and striatum to perform navigation in a simulated environment. Their model successfully reproduced many experimental results, such as rescaling of firing fields of place cells and grid cells and solved water maze task by reinforcement learning based on spatial representations. Both (Yuan et al., <xref ref-type="bibr" rid="B55">2015</xref>) and (Mulas et al., <xref ref-type="bibr" rid="B29">2016</xref>) used a model of pure positional grid cells in layer II of MEC (Burak and Fiete, <xref ref-type="bibr" rid="B3">2009</xref>) to track the position of a robot and used visual inputs to anchor grids in the environment. They tested their system in indoor environments of dozens of meters in side. The robustness is yet to be verified in large-scale environments.</p>
<p>The theory of neural dynamics has been extensively adopted to model brain functions, such as memory (Treves, <xref ref-type="bibr" rid="B50">1990</xref>; Mi et al., <xref ref-type="bibr" rid="B27">2017</xref>), navigation (Gerstner and Abbott, <xref ref-type="bibr" rid="B15">1997</xref>), and sensory integration (Zhang et al., <xref ref-type="bibr" rid="B57">2016</xref>). The theory of neural dynamics has become popular in robotics to enhance the cognitive abilities of robots (Erlhagen and Bicho, <xref ref-type="bibr" rid="B12">2006</xref>), such as obstacle avoidance (Xiao et al., <xref ref-type="bibr" rid="B54">2017</xref>), coordinated path tracking (Xiao and Zhang, <xref ref-type="bibr" rid="B53">2016</xref>), motion tracking (Ding et al., <xref ref-type="bibr" rid="B8">2017</xref>), and grasping (Knips et al., <xref ref-type="bibr" rid="B21">2017</xref>).</p>
<p>In this work, we follow the research line of neural dynamics and developed a cognitive mapping model for mobile robots, taking advantage of the coding strategies of the spatial memory circuits in mammalian brains. The key components of the proposed model include HD cells, conjunctive grid cells, and local view cells. Both HD cells and conjunctive grid cells are modeled by continuous attractor networks that operate on the same principles. More specifically, HD cells in the model represent arbitrary conjunctions of head directions and rotations of the animals. Due to the asymmetric recurrent connections and the network dynamics, intrinsically rotating patterns spontaneously emerge in the network. Angular velocity inputs to the network pick up patterns rotating with appropriate velocities and track the head direction of the robot. Similarly, conjunctive grid cells in the model encode conjunctions of positions and translations in the two-dimensional environment. The inputs from HD cells and speed cells activate a subset of the conjunctive grid cells that produce triangular patterns moving intrinsically in the neural tissue with velocity proportional to the running velocity of the robot. Both HD cells and conjunctive grid cells get inputs from local view cells, which are activated once the robot moves to a familiar scene and which provide anchoring cues. The mapping performance of the proposed model is demonstrated on a 66&#x02009;km long urban car journey. Our cognitive mapping model generates a coherent map of the outdoor environment using single camera image data. The experimental results show that our model is stable and robust in large-scale outdoor environment of variable light conditions and terrains.</p>
<p>The major contributions in this paper are twofold. First, the neural network framework presented in this paper is based on recent experimental studies on the hippocampal&#x02013;entorhinal circuits and aims to model layer III and deep layers of the MEC and the hippocampus. In the system, conjunctive HD-by-velocity cells and conjunctive grid-by-velocity cells work hand in hand to integrate movement and sensory information and build a large-scale map. Second, the neural dynamics of the hippocampal&#x02013;entorhinal circuits is modeled. Provided with the inputs from local view cells, the neural dynamics of the system functions as a general mechanism for error correction and pattern completion.</p>
<p>This paper proceeds as follows. In Section <xref ref-type="sec" rid="S2">2</xref>, we describe the HD cell model, the conjunctive grid cell model, and other core components of the robot navigation model. The experimental setup and the resulting neural activities, rate maps, and the semi-metric topological map are presented in Section <xref ref-type="sec" rid="S3">3</xref>. Section <xref ref-type="sec" rid="S4">4</xref> discusses the results and points to future directions, with a brief conclusion in Section <xref ref-type="sec" rid="S5">5</xref>.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<label>2</label> <title>Materials and Methods</title>
<p>In this study, we propose a cognitive mapping model that has similar regions and cell types as in the entorhinal&#x02013;hippocampal circuits of mammalian brains. Localization is the key function of the proposed model and is performed by the HD cells and the conjunctive grid cells. The HD cells and the conjunctive grid cells in the model correspond to layer III, V, and VI of the MEC of mammalian brains. They integrate rotational or translational velocity and form directional or positional codes, respectively. The local view cells in the model correspond to the retrosplenial cortex or visual cortex and relay visual inputs to the MEC. Cognitive map in the system corresponds to the hippocampus.</p>
<sec id="S2-1">
<label>2.1</label> <title>Neural Network Architecture</title>
<p>The architecture of the proposed model is shown in Figure <xref ref-type="fig" rid="F1">1</xref>. The HD cells in the model form a ring attractor in the neural space. A bump of activity emerges in the network. The phase of the bump represents the robot&#x02019;s real angle of the head direction in the physical environment up to a constant shift. Tuned angular velocity inputs active a subset of the HD neurons, and the bump rotates intrinsically with the same velocity as the angular velocity of the robot.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>The neural network architecture of the model and the diagram of information flow. The HD-by-velocity cells, updated by angular velocity, represent head directions. The grid-by-velocity cells receiving translational velocity, converted from linear speed and the head direction representations, provide positional representation, which is in turn utilized to build a cognitive map. The acitivities of the HD-by-velocity cells and the grid-by-velocity cells are encoded by heat maps with red for high activity and blue for no activity.</p></caption>
<graphic xlink:href="fnbot-11-00061-g001.tif"/>
</fig>
<p>The conjunctive grid cells in the model form a torus attractor due to the periodic boundary conditions of the neural space. The phase of the grid pattern encodes the position of the robot in the physical environment up to modulo operation. The conjunctive grid cells are activated by the head direction cells and the speed cells, which are both available in the MEC of rodents. The grid pattern in the network moves with a velocity proportional to the movement of the animal in the physical environment. Thus, the HD cells and the conjunctive grid cells work hand in hand and form population codes of the head direction and the position of the robot in the environment.</p>
<p>The conjunctive representations of space and movement allow the networks to have stable states for all movement conditions. This coding property endows the network robustness regardless the movement conditions.</p>
</sec>
<sec id="S2-2">
<label>2.2</label> <title>HD-by-Velocity Cell Model</title>
<p>The HD-by-velocity cells integrate angular velocity and represent one-dimensional head direction. This is achieved by conjunctive encoding of head direction and rotation.</p>
<sec id="S2-2-1">
<label>2.2.1</label> <title>Neural Representation of Head Direction and Angular Velocity</title>
<p>Each unit in the network is labeled by its coordinate (<italic>&#x003B8;, &#x003BD;</italic>) on a two-dimensional manifold. <italic>&#x003B8;</italic> &#x02208; [0, 2&#x003C0;) is the internal representation of head directions in the environment. Not that <italic>&#x003B8;</italic> has periodic boundary condition. <italic>&#x003BD;</italic> &#x02208; [&#x02212;<italic>L<sub>r</sub>, L<sub>r</sub></italic>] encodes angular velocity. Both <italic>&#x003BD;</italic> and <italic>&#x003B8;</italic> are dimensionless quantities in the neural space. They reflect the real angular velocity and the head direction of the robot. We design the connection weights between units to sustain a single bump of activity both in the direction dimension <italic>&#x003B8;</italic> and in the rotation dimension <italic>&#x003BD;</italic>. The strength of the connection from a presynaptic unit (<italic>&#x003B8;&#x02032;</italic>, <italic>&#x003BD;&#x02032;</italic>) to a postsynaptic unit (<italic>&#x003B8;, &#x003BD;</italic>) is written as
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mi>J</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>&#x003BD;</mml:mi><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup><mml:mo class="MathClass-punc">,</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msub><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mtext>cos</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:mfenced><mml:mtext>cos</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mn>&#x003BB;</mml:mn><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003BD;</mml:mi><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow></mml:mfenced><mml:mspace width="0.3em" class="thinspace"/><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where <italic>J</italic><sub>0</sub>&#x02009;&#x0003C;&#x02009;0 is a uniform inhibition, <italic>J</italic><sub>1</sub>&#x02009;&#x0003E;&#x02009;0 describes the interaction strength, and &#x003BB; defines the spread of velocity tuning. The connection weight from the unit at <italic>&#x003B8;&#x02032;</italic> to the unit at <italic>&#x003B8;</italic> in the direction dimension is asymmetric, since the postsynaptic unit that is maximally activated by <italic>&#x003B8;&#x02032;</italic> is not centered at <italic>&#x003B8;&#x02032;</italic> but at <italic>&#x003B8;&#x02032;</italic>&#x02009;&#x0002B;&#x02009;<italic>&#x003BD;&#x02032;</italic>. The asymmetry causes the bump to move along the direction dimension with a velocity determined by <italic>&#x003BD;&#x02032;</italic>.</p>
</sec>
<sec id="S2-2-2">
<label>2.2.2</label> <title>Network Dynamics</title>
<p>We consider the activity of the network driven by velocity and sensory inputs. The firing rate <italic>m</italic>(<italic>&#x003B8;, &#x003BD;</italic>) of the unit at coordinate (<italic>&#x003B8;, &#x003BD;</italic>) is governed by the differential equation
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mn>&#x003C4;</mml:mn><mml:mover accent="true"><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x002D9;</mml:mo></mml:mover><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:mi>f</mml:mi><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-op">&#x0222C;</mml:mo><mml:mi>D</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mi>D</mml:mi><mml:mi>&#x003BD;</mml:mi><mml:mi>J</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>&#x003BD;</mml:mi><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup><mml:mo class="MathClass-punc">,</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup><mml:mo class="MathClass-punc">,</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="italic">view</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfenced><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where <italic>I<sub>&#x003BD;</sub></italic> and <italic>I<sub>view</sub></italic> are the velocity tuned input and the calibration current injected by local view cells, respectively, which we will explain in detail below. &#x003C4; is the time constant set to 10&#x02009;ms. <italic>f</italic> (<italic>x</italic>) is a threshold-linear function: <italic>f</italic> (<italic>x</italic>)&#x02009;&#x02261;&#x02009;[<italic>x</italic>]<sub>&#x0002B;</sub>&#x02009;&#x0003D;&#x02009;<italic>x</italic> when <italic>x</italic>&#x02009;&#x0003E;&#x02009;0 and 0 otherwise. The shorthand notations are used <inline-formula><mml:math id="M3"><mml:mo class="MathClass-op">&#x0222B;</mml:mo><mml:mi>D</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:mfrac><mml:msubsup><mml:mrow><mml:mo class="MathClass-op">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:msubsup><mml:mspace width="0.3em"/><mml:mi>d</mml:mi><mml:mi>&#x003B8;</mml:mi></mml:math></inline-formula>, and <inline-formula><mml:math id="M4"><mml:mo class="MathClass-op">&#x0222B;</mml:mo><mml:mi>D</mml:mi><mml:mi>&#x003BD;</mml:mi><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:msubsup><mml:mrow><mml:mo class="MathClass-op">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:mspace width="0.3em"/><mml:mi>d</mml:mi><mml:mi>&#x003BD;</mml:mi></mml:math></inline-formula>.</p>
</sec>
<sec id="S2-2-3">
<label>2.2.3</label> <title>Angular Velocity Inputs</title>
<p>In order to integrate angular velocity, the activity bump in the HD network should be placed at appropriate position in the velocity dimension of the neural space. The desired bump location in the <italic>&#x003BD;</italic> dimension corresponding to a given external angular velocity <italic>V</italic> is (Si et al., <xref ref-type="bibr" rid="B37">2014</xref>)
<disp-formula id="E3"><label>(3)</label><mml:math id="M5"><mml:mi>u</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mtext>arctan</mml:mtext><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mn>&#x003C4;</mml:mn><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>Here, &#x003C4; is the time constant defined in equation (<xref ref-type="disp-formula" rid="E2">2</xref>).</p>
<p>The velocity input to the HD-by-velocity units is simply modeled by a tuning function of Gaussian shape
<disp-formula id="E4"><label>(4)</label><mml:math id="M6"><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003BD;</mml:mi><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msub><mml:mfenced separators="" open="[" close="]"><mml:mrow><mml:mn>1</mml:mn><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mn>&#x003F5;</mml:mn><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:mn>&#x003F5;</mml:mn><mml:mspace width="0.3em"/><mml:mtext>exp</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003BD;</mml:mi><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mi>u</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:msubsup><mml:mrow><mml:mn>&#x003C3;</mml:mn></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced><mml:mspace width="0.3em" class="thinspace"/><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>Here, <italic>I<sub>r</sub></italic> is the amplitude of the rotational velocity input, &#x003F5; defines the strength of the velocity tuning, and &#x003C3;<italic><sub>r</sub></italic> is the sharpness of the velocity tuning.</p>
</sec>
<sec id="S2-2-4">
<label>2.2.4</label> <title>Estimation of the Head Direction</title>
<p>The center of the activity bump on the <italic>&#x003B8;</italic> axis and the <italic>&#x003BD;</italic> axis encodes the head direction and angular velocity of the robot. We define Fourier transformations to recover the head direction and the angular velocity of the robot from the network state
<disp-formula id="E5"><label>(5)</label><mml:math id="M7"><mml:mn>&#x003C8;</mml:mn><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mi>&#x02220;</mml:mi><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-op">&#x0222C;</mml:mo><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mspace width="0.3em"/><mml:mtext>exp</mml:mtext><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mi>D</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mi>D</mml:mi><mml:mi>&#x003BD;</mml:mi></mml:mrow></mml:mfenced><mml:mspace width="0.3em" class="thinspace"/><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
<disp-formula id="E6"><label>(6)</label><mml:math id="M8"><mml:mn>&#x003D5;</mml:mn><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mi>&#x02220;</mml:mi><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-op">&#x0222C;</mml:mo><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mspace width="0.3em"/><mml:mtext>exp</mml:mtext><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mn>&#x003BB;&#x003BD;</mml:mn></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mi>D</mml:mi><mml:mi>&#x003B8;</mml:mi><mml:mi>D</mml:mi><mml:mi>&#x003BD;</mml:mi></mml:mrow></mml:mfenced></mml:mrow><mml:mrow><mml:mn>&#x003BB;</mml:mn></mml:mrow></mml:mfrac><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>Here, <italic>i</italic> is the imaginary unit, and function &#x02220;(<italic>Z</italic>) takes the angle of a complex number <italic>Z</italic>. &#x003C8; &#x02208; [0, 2&#x003C0;) is the estimated phase of the bump in the direction axis of the neural space and corresponds to the head direction of the robot in the physical space. &#x003D5; &#x02208; (&#x02212;<italic>L<sub>r</sub>, L<sub>r</sub></italic>) is the estimated phase of the bump in the velocity axis of the neural space and can be mapped to the angular velocity of the robot in the physical space by inverting equation (<xref ref-type="disp-formula" rid="E3">3</xref>)
<disp-formula id="E7"><label>(7)</label><mml:math id="M9"><mml:mi>V</mml:mi><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mtext>tan</mml:mtext><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mn>&#x003D5;</mml:mn></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>&#x003C4;</mml:mn></mml:mrow></mml:mfrac><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>Note that <italic>L<sub>r</sub></italic> should be chosen large enough, so that the recovered velocity <italic>V</italic> is able to represent all possible angular velocities of the robot.</p>
</sec>
</sec>
<sec id="S2-3">
<label>2.3</label> <title>Grid-by-Velocity Cell Model</title>
<p>Now, we expand our HD-by-velocity cell model to do path integration in two-dimensional environments. To this end, grid-by-velocity cells need represent two-dimensional spatial locations and two-dimensional velocities. The units in the network are wired with appropriate connection profiles, so that hexagonal grid firing pattern is generated and translated in the spatial dimension of the neural space.</p>
<sec id="S2-3-5">
<label>2.3.1</label> <title>Neural Representation of Position and Velocity</title>
<p>Units in the grid-by-velocity network is labeled by coordinates <inline-formula><mml:math id="M10"><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> in a four-dimensional neural space. <inline-formula><mml:math id="M11"><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-punc">,</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> represents two-dimensional positions in the environment. <inline-formula><mml:math id="M12"><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-punc">,</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> encodes the velocity components in the environment. We assume <inline-formula><mml:math id="M51"><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:math></inline-formula> has periodic boundary conditions, i.e., <italic>&#x003B8;<sub>x</sub></italic>, <italic>&#x003B8;<sub>y</sub></italic>&#x02009;&#x02208;&#x02009;[0, 2&#x003C0;). <italic>&#x003BD;<sub>x</sub></italic> and <italic>&#x003BD;<sub>y</sub></italic> are chosen in [&#x02212;<italic>L<sub>t</sub>, L<sub>t</sub></italic>].</p>
<p>The connection weight from unit <inline-formula><mml:math id="M13"><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> to <inline-formula><mml:math id="M14"><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> is described as
<disp-formula id="E8"><label>(8)</label><mml:math id="M15"><mml:mi>J</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msub><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>J</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mtext>cos</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mi>k</mml:mi><mml:msqrt><mml:mrow><mml:mstyle displaystyle='true'><mml:munder><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo class="MathClass-rel">&#x02208;</mml:mo><mml:mrow><mml:mo class="MathClass-open">&#x0007B;</mml:mo><mml:mrow><mml:mi>x</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>y</mml:mi></mml:mrow><mml:mo class="MathClass-close">&#x0007D;</mml:mo></mml:mrow></mml:mrow></mml:munder></mml:mstyle><mml:mspace width="0.3em"/><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msub><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msub><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:msup><mml:mrow><mml:mo class="MathClass-rel">&#x0007C;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:msqrt></mml:mrow></mml:mfenced><mml:mtext>cos</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mn>&#x003BB;</mml:mn><mml:msqrt><mml:mrow><mml:mstyle displaystyle='true'><mml:munder><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo class="MathClass-rel">&#x02208;</mml:mo><mml:mrow><mml:mo class="MathClass-open">&#x0007B;</mml:mo><mml:mrow><mml:mi>x</mml:mi><mml:mo class="MathClass-punc">,</mml:mo><mml:mi>y</mml:mi></mml:mrow><mml:mo class="MathClass-close">&#x0007D;</mml:mo></mml:mrow></mml:mrow></mml:munder></mml:mstyle><mml:mspace width="0.3em"/><mml:msup><mml:mrow><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msub><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:msqrt></mml:mrow></mml:mfenced><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where the integer <italic>k</italic>&#x02009;&#x0003D;&#x02009;2, is chosen so that the network accommodates two bumps both in <italic>&#x003B8;<sub>x</sub></italic> axis and in <italic>&#x003B8;<sub>y</sub></italic> axis. There is only one bump in each of the velocity dimensions, however. &#x0007C;&#x0007C;<italic>d</italic>&#x0007C;&#x0007C; is the distance on a circle: &#x0007C;&#x0007C;<italic>d</italic>&#x0007C;&#x0007C;&#x02009;&#x0003D;&#x02009;mod(<italic>d</italic>&#x02009;&#x0002B;&#x02009;&#x003C0;, 2&#x003C0;)&#x02009;&#x02212;&#x02009;&#x003C0;, and mod(<italic>x, y</italic>) &#x02208; [0, <italic>y</italic>) gives x modulo y.</p>
</sec>
<sec id="S2-3-6">
<label>2.3.2</label> <title>Network Dynamics</title>
<p>Although grid-by-velocity units are organized with different manifold structure from the HD-by-velocity units, they share the same intrinsic dynamics as in equation (<xref ref-type="disp-formula" rid="E2">2</xref>)
<disp-formula id="E9"><label>(9)</label><mml:math id="M16"><mml:mn>&#x003C4;</mml:mn><mml:mover accent="true"><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x002D9;</mml:mo></mml:mover><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:mi>f</mml:mi><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-op">&#x0222C;</mml:mo><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mi>J</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mo class="MathClass-rel">&#x02032;</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="italic">view</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfenced><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>Note that <inline-formula><mml:math id="M17"><mml:mo class="MathClass-op">&#x0222B;</mml:mo><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mn>4</mml:mn><mml:msup><mml:mrow><mml:mi>&#x003C0;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:msubsup><mml:mrow><mml:mo class="MathClass-op">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:msubsup><mml:mspace width="0.3em"/><mml:msubsup><mml:mrow><mml:mo class="MathClass-op">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:msubsup><mml:mspace width="0.3em"/><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, and <inline-formula><mml:math id="M18"><mml:mo class="MathClass-op">&#x0222B;</mml:mo><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mn>4</mml:mn><mml:msubsup><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac><mml:msubsup><mml:mrow><mml:mo class="MathClass-op">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:mspace width="0.3em"/><mml:msubsup><mml:mrow><mml:mo class="MathClass-op">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msubsup><mml:mspace width="0.3em"/><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>.</p>
</sec>
<sec id="S2-3-7">
<label>2.3.3</label> <title>Translational Velocity Inputs</title>
<p>In order to perform accurate path integration, the velocity of the moving bumps in the neural space should always be kept proportional to the velocity of the robot in the physical space. The network requires tuned velocity input to pin the activity bumps at appropriate positions on the velocity axes in the neural space, so that the bumps move with the desired velocity.</p>
<p>The translational velocity <inline-formula><mml:math id="M19"><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-punc">,</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> of the robot is obtained by projecting the running speed to the axes of the reference frame using the HD estimated from the HD-by-velocity units (equation (<xref ref-type="disp-formula" rid="E5">5</xref>)). In the brain, running speed is encoded by the speed cells in the MEC. Given the translational velocity of the robot, the desired positions on the velocity axes in the neural space is given by <inline-formula><mml:math id="M20"><mml:mover accent="true"><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> (Si et al., <xref ref-type="bibr" rid="B37">2014</xref>)
<disp-formula id="E10"><label>(10)</label><mml:math id="M21"><mml:mover accent="true"><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:mfrac><mml:mtext>arctan</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:mn>&#x003C4;</mml:mn><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>S</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:mfenced><mml:mspace width="0.3em" class="thinspace"/><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where the function arctan operates on each dimension of V &#x02192;. S is a scaling factor between the external velocity of the robot in the physical environment and the velocity of the bumps in the neural space. S determines the spacing between the neighboring fields of a grid firing pattern in the environment.</p>
<p>The velocity-tuned inputs to the grid-by-velocity units take a Gaussian form for simplicity
<disp-formula id="E11"><label>(11)</label><mml:math id="M22"><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mfenced separators="" open="[" close="]"><mml:mrow><mml:mn>1</mml:mn><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mn>&#x003F5;</mml:mn><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:mn>&#x003F5;</mml:mn><mml:mspace width="0.3em"/><mml:mtext>exp</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mfrac><mml:mrow><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mo class="MathClass-rel">&#x0007C;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:msubsup><mml:mrow><mml:mn>&#x003C3;</mml:mn></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced><mml:mspace width="0.3em" class="thinspace"/><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where &#x0007C;&#x022C5;&#x0007C; is the Euclidean norm of a vector. <italic>I<sub>t</sub></italic> is the amplitude of the translational velocity input. &#x003F5; is the strength of velocity tuning as in equation (<xref ref-type="disp-formula" rid="E4">4</xref>), and &#x003C3;<italic><sub>t</sub></italic> the width of the translational velocity tuning.</p>
</sec>
<sec id="S2-3-8">
<label>2.3.4</label> <title>Estimation of the Bump Position and Localization</title>
<p>The center of the bumps on the neural manifold encodes spatial position and translational velocity. We use Fourier transformations to estimate the center of the bumps on the four-dimensional neural manifold.</p>
<p>First, we project the firing activity to the projection axes defined by the grid pattern in the spatial dimensions of the neural manifold and obtain the phases of the grid pattern along the projection axes
<disp-formula id="E12"><label>(12)</label><mml:math id="M26"><mml:msub><mml:mrow><mml:mn>&#x003C8;</mml:mn></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mi>&#x02220;</mml:mi><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-op">&#x0222C;</mml:mo><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mtext>exp</mml:mtext><mml:mspace width="0.3em" class="thinspace"/><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ik</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:mfrac><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where <italic>i</italic> is the imaginary unit, <inline-formula><mml:math id="M27"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> the unit direction of the <italic>j</italic>-th projection axis, namely <inline-formula><mml:math id="M28"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo class="MathClass-open">[</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo class="MathClass-punc">,</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo class="MathClass-close">]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>, <inline-formula><mml:math id="M29"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo class="MathClass-open">[</mml:mo><mml:mrow><mml:mtext>sin</mml:mtext><mml:mn>&#x003B1;</mml:mn><mml:mo class="MathClass-punc">,</mml:mo><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mtext>cos</mml:mtext><mml:mn>&#x003B1;</mml:mn></mml:mrow><mml:mo class="MathClass-close">]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>, <inline-formula><mml:math id="M30"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo class="MathClass-open">[</mml:mo><mml:mrow><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mtext>sin</mml:mtext><mml:mn>&#x003B1;</mml:mn><mml:mo class="MathClass-punc">,</mml:mo><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mtext>cos</mml:mtext><mml:mn>&#x003B1;</mml:mn></mml:mrow><mml:mo class="MathClass-close">]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> (see Figure <xref ref-type="fig" rid="F2">2</xref>). &#x003B1;, determined by the grid pattern in the torus, is the angle of the second grid axis. In our case, &#x003B1;&#x02009;&#x0003D;&#x02009;arctan(2) due to the square shape of the torus. <inline-formula><mml:math id="M31"><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> is the respective wavelength of the grid pattern along the projection axis <italic>j</italic>, i.e., <inline-formula><mml:math id="M32"><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>&#x02009;&#x0003D;&#x02009;1, <inline-formula><mml:math id="M33"><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>&#x02009;&#x0003D;&#x02009;<inline-formula><mml:math id="M34"><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>&#x02009;&#x0003D;&#x02009;sin &#x003B1;.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Estimation of bump position on the torus. The torus neural space of (<italic>&#x003B8;<sub>x</sub></italic>, <italic>&#x003B8;<sub>y</sub></italic>) is unfolded and shown as the black box. The grid pattern is projected by Fourier transforms along projection axes <inline-formula><mml:math id="M23"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula><mml:math id="M24"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>. The center of a bump (red dot) in the torus is represented by the phase in each projection axis. The position of the bump on the torus is recovered by mapping the phases in the projection axes to the neural space according to the trigonometric relation defined by &#x003B1;, which is the angle of the grid axis perpendicular to <inline-formula><mml:math id="M25"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>.</p></caption>
<graphic xlink:href="fnbot-11-00061-g002.tif"/>
</fig>
<p>Second, we estimate the center of the bumps in the spatial dimensions of the neural manifold according to the trigonometric relation between the projection axes and the spatial axes. As shown in Figure <xref ref-type="fig" rid="F2">2</xref>, the center of the bumps in the <italic>&#x003B8;<sub>y</sub></italic> dimension is simply
<disp-formula id="E13"><label>(13)</label><mml:math id="M35"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msub><mml:mrow><mml:mn>&#x003C8;</mml:mn></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>By mapping the scaled phase <inline-formula><mml:math id="M36"><mml:msub><mml:mrow><mml:mn>&#x003C8;</mml:mn></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula><mml:math id="M37"><mml:msub><mml:mrow><mml:mn>&#x003C8;</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> to the <italic>&#x003B8;<sub>x</sub></italic> direction, the center of the bumps in the <italic>&#x003B8;<sub>x</sub></italic> dimension is given by
<disp-formula id="E14"><label>(14)</label><mml:math id="M38"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mn>&#x003C8;</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mtext>sin</mml:mtext><mml:mi>&#x003B1;</mml:mi></mml:mrow></mml:mfrac><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mn>&#x003C8;</mml:mn></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mtext>tan</mml:mtext><mml:mi>&#x003B1;</mml:mi></mml:mrow></mml:mfrac><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>The phases estimated in equation (<xref ref-type="disp-formula" rid="E12">12</xref>) are constrained by the periodic boundary conditions. In order to estimate the true position of the robot in the environment, the phases should be unwrapped by considering the phases of two consecutive time steps. After unwrapping, the cumulative phases are used in equations (<xref ref-type="disp-formula" rid="E13">13</xref>) and (<xref ref-type="disp-formula" rid="E14">14</xref>) to estimate the position of the robot
<disp-formula id="E15"><label>(15)</label><mml:math id="M39"><mml:mover accent="true"><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:mo class="MathClass-bin">&#x02215;</mml:mo><mml:mi>k</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
<disp-formula id="E16"><label>(16)</label><mml:math id="M40"><mml:mover accent='true'><mml:mi>y</mml:mi><mml:mo>&#x0005E;</mml:mo></mml:mover><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:mo class="MathClass-bin">&#x02215;</mml:mo><mml:mi>k</mml:mi></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where <inline-formula><mml:math id="M41"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M42"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:math></inline-formula> are the cumulative phases after unwrapping at time t. The scaling factor <inline-formula><mml:math id="M43"><mml:mfrac><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:mo class="MathClass-bin">&#x02215;</mml:mo><mml:mi>k</mml:mi></mml:mrow></mml:mfrac></mml:math></inline-formula> is the ratio between the velocity of the robot in the physical environment and the velocity of bumps in the network.</p>
<p>To estimate the center of the bumps in the velocity dimensions, the same method as equation (<xref ref-type="disp-formula" rid="E6">6</xref>) is used in the following form
<disp-formula id="E17"><label>(17)</label><mml:math id="M44"><mml:msub><mml:mrow><mml:mn>&#x003D5;</mml:mn></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mi>&#x02220;</mml:mi><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-op">&#x0222C;</mml:mo><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mi>D</mml:mi><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mi>m</mml:mi><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-punc">,</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mspace width="0.3em"/><mml:mtext>exp</mml:mtext><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mn>&#x003BB;</mml:mn><mml:msub><mml:mrow><mml:mi>&#x003BD;</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow></mml:mfenced></mml:mrow><mml:mrow><mml:mn>&#x003BB;</mml:mn></mml:mrow></mml:mfrac><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where <italic>j</italic> &#x02208; &#x0007B;<italic>x, y</italic>&#x0007D;. &#x003D5;<italic><sub>j</sub></italic> &#x02208; (&#x02212;<italic>L<sub>t</sub>, L<sub>t</sub></italic>) is the position of the bump on the axes <italic>&#x003BD;<sub>j</sub></italic>.</p>
<p>According to the mapping between the velocity of the bumps in the neural space and the velocity of the robot, the velocity of the robot in the environment is given by
<disp-formula id="E18"><label>(18)</label><mml:math id="M45"><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mfrac><mml:mrow><mml:mi>S</mml:mi><mml:mtext>tan</mml:mtext><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:msub><mml:mrow><mml:mn>&#x003D5;</mml:mn></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:mn>&#x003C4;</mml:mn></mml:mrow></mml:mfrac><mml:mo class="MathClass-punc">.</mml:mo></mml:math></disp-formula></p>
<p>Note that <italic>L<sub>t</sub></italic> should be chosen large enough to cover the maximal translation velocity experienced by the robot.</p>
</sec>
</sec>
<sec id="S2-4">
<label>2.4</label> <title>Estimating Angular and Translational Velocity from Visual Inputs</title>
<p>The angular velocity and the translational speed of the robot are estimated by matching two consecutive frames from the camera, according to the method described in Milford and Wyeth (<xref ref-type="bibr" rid="B28">2008</xref>). The estimated angular velocity is utilized by the HD-by-velocity units to form HD representations (equation (<xref ref-type="disp-formula" rid="E4">4</xref>)). The translational velocity is obtained by combining the HD estimation from the HD-by-velocity network and the translational speed estimated from the image sequence, and is fed to the grid-by-velocity network to perform path integration (equation (<xref ref-type="disp-formula" rid="E11">11</xref>)).</p>
</sec>
<sec id="S2-5">
<label>2.5</label> <title>Calibration from Local View Cells</title>
<p>Since movement estimation is subject to noise, both the head direction cell network and the grid cell network need calibration by visual inputs. Following Milford and Wyeth (<xref ref-type="bibr" rid="B28">2008</xref>), we extract local view templates from images to encode the scenes observed from the camera. If a view is sufficiently different from previously observed views, a new view cell is added to the system, together with the view template and the estimated bump positions &#x003C8;, <inline-formula><mml:math id="M46"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula><mml:math id="M47"><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>. If the robot comes back to a previously visited location, and a familiar view reappears, the corresponding local view cell is activated and provides inputs to the networks. We model the current injected from a view cell to the HD cell network as
<disp-formula id="E19"><label>(19)</label><mml:math id="M48"><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="italic">view</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mtext>exp</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mfrac><mml:mrow><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mn>&#x003C8;</mml:mn><mml:mo class="MathClass-rel">&#x0007C;</mml:mo><mml:msup><mml:mrow><mml:mo class="MathClass-rel">&#x0007C;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:msubsup><mml:mrow><mml:mn>&#x003C3;</mml:mn></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow></mml:mfenced><mml:mspace width="0.3em" class="thinspace"/><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where <italic>I<sub>d</sub></italic> is the amplitude of the injected current to the HD cell network, &#x003C8; is the associated phase with the view template, &#x003C3;<italic><sub>d</sub></italic> is sharpness of the Gaussian tuning.</p>
<p>For the conjunctive grid cell network, the input from a view cell is given by
<disp-formula id="E20"><label>(20)</label><mml:math id="M49"><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="italic">view</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mfrac><mml:mstyle displaystyle='true'><mml:munderover><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo class="MathClass-rel">&#x0003D;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:munderover></mml:mstyle><mml:mspace width="0.3em"/><mml:mtext>cos</mml:mtext><mml:mfenced separators="" open="(" close=")"><mml:mrow><mml:mi>k</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo class="MathClass-open">(</mml:mo><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover><mml:mo class="MathClass-bin">&#x02212;</mml:mo><mml:mover accent="true"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-close">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:mfenced><mml:mo class="MathClass-bin">&#x0002B;</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:mfenced><mml:mspace width="0.3em" class="thinspace"/><mml:mo class="MathClass-punc">,</mml:mo></mml:math></disp-formula>
where <italic>I<sub>p</sub></italic> is the amplitude of the injected current, <italic>C</italic> is a constant to adjust the baseline, <inline-formula><mml:math id="M50"><mml:mover accent="true"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mo class="MathClass-op">&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mo class="MathClass-op">&#x0005E;</mml:mo></mml:mover></mml:math></inline-formula> is the phase associated with the local view template (see equations (<xref ref-type="disp-formula" rid="E13">13</xref>) and (<xref ref-type="disp-formula" rid="E14">14</xref>)), and other parameters are the same as in equation (<xref ref-type="disp-formula" rid="E12">12</xref>). Note that the injected current from a local view cell depends only on <italic>&#x003B8;</italic>&#x02192;, and is the same for different <italic>&#x003BD;</italic>&#x02192;.</p>
</sec>
<sec id="S2-6">
<label>2.6</label> <title>Neural Network Parameters</title>
<p>The HD-by-velocity units are regularly arranged on a rectangular grid in the two-dimensional manifold of direction and rotation, with 51 units in the direction dimension and 25 units in the rotation dimension. All together, there are 1,275 HD-by-velocity units. <italic>L<sub>r</sub></italic>&#x02009;&#x0003D;&#x02009;0.0095 is sufficient to represent the maximal angular velocity of the vehicle.</p>
<p>The grid-by-velocity units are regularly arranged on a rectangular grid in the four-dimensional neuronal manifold of position and velocity. There are 15 units in each of the positional dimensions and 7 units in each of the velocity dimensions, resulting in 11,025 grid-by-velocity units. To represent the maximal experienced translational velocity of the robot, <italic>L<sub>t</sub></italic> is chosen to be 0.3.</p>
<p>Table <xref ref-type="table" rid="T1">1</xref> summarizes the parameters used in the model.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Values of the parameters used in the system.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left">Parameter</th>
<th align="center"><italic>J</italic><sub>0</sub></th>
<th align="center"><italic>J</italic><sub>1</sub></th>
<th align="center"><italic>&#x003BB;</italic></th>
<th align="center"><italic>L<sub>r</sub></italic></th>
<th align="center"><italic>I<sub>r</sub></italic></th>
<th align="center"><italic>&#x003F5;</italic></th>
<th align="center"><italic>&#x003C3;<sub>r</sub></italic></th>
<th align="center"><italic>J<sub>k</sub></italic></th>
<th align="center"><italic>k</italic></th>
<th align="center"><italic>L<sub>t</sub></italic></th>
<th align="center"><italic>S</italic></th>
<th align="center"><italic>I<sub>t</sub></italic></th>
<th align="center"><italic>&#x003C3;<sub>t</sub></italic></th>
<th align="center"><italic>I<sub>d</sub></italic></th>
<th align="center"><italic>&#x003C3;<sub>d</sub></italic></th>
<th align="center"><italic>I<sub>p</sub></italic></th>
<th align="center">C</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left">Value</td>
<td align="center">&#x02212;60</td>
<td align="center">50</td>
<td align="center">0.8</td>
<td align="center">0.0095</td>
<td align="center">50</td>
<td align="center">0.8</td>
<td align="center">0.012</td>
<td align="center">50</td>
<td align="center">2</td>
<td align="center">0.3</td>
<td align="center">30</td>
<td align="center">60</td>
<td align="center">0.1</td>
<td align="center">60</td>
<td align="center">2.19</td>
<td align="center">200</td>
<td align="center">0.5</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="S2-7">
<label>2.7</label> <title>Map Representation</title>
<p>A map of the environment can be formed by reading out the spatial codes in the grid-by-velocity network. We could resort to place cell models (Solstad et al., <xref ref-type="bibr" rid="B41">2006</xref>; Si and Treves, <xref ref-type="bibr" rid="B38">2009</xref>; Cheng and Frank, <xref ref-type="bibr" rid="B6">2011</xref>) to build cognitive maps. In this study, we simply adopt the experience map representation (Milford and Wyeth, <xref ref-type="bibr" rid="B28">2008</xref>), which is a topological map storing the spatial positions and their transitions measured from visual odometry. A graph relaxation algorithm (Duckett et al., <xref ref-type="bibr" rid="B9">2002</xref>) is applied to optimize the topological map when the robot comes back to a familiar location and closes a loop (Milford and Wyeth, <xref ref-type="bibr" rid="B28">2008</xref>).</p>
</sec>
<sec id="S2-8">
<label>2.8</label> <title>Implementation of the Cognitive Mapping Model</title>
<p>The proposed cognitive mapping model is implemented in the C&#x0002B;&#x0002B; language and is run in the Robot Operating System (ROS) Indigo on Ubuntu 14.04 LTS (Trusty). We organize the software architecture of our cognitive mapping model into five nodes (Figure <xref ref-type="fig" rid="F3">3</xref>), similarly as that of the open RatSLAM system (Ball et al., <xref ref-type="bibr" rid="B2">2013</xref>), for the ease of reusing their visual processing algorithm and comparing with their results.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>The software architecture of the cognitive mapping model. The sensor/bagfile node provides camera images. Velocities are estimated by the visual odometry node. The local view cell node analyzes whether the current view is new or not. Path integration and decisions of creating links and nodes are performed by the spatial memory network node. The topological map is generated by the experience map node.</p></caption>
<graphic xlink:href="fnbot-11-00061-g003.tif"/>
</fig>
<p>The visual odometry node measures the angular velocity and the translational speed according to the changes of visual scenes. It receives frames from the ROS either from a camera or from stored data in a bagfile.</p>
<p>The local view cell node determines whether the current scene is new or is a familiar one. It provides calibration current to the networks in the spatial memory network node.</p>
<p>The spatial memory network node includes the head-direction cell network and the conjunctive grid cell network. This node receives two types of ROS messages as inputs: odometry and view templates. As shown in Section <xref ref-type="sec" rid="S2-2-3">2.2.3</xref>, <xref ref-type="sec" rid="S2-3-7">2.3.3</xref>, and <xref ref-type="sec" rid="S2-5">2.5</xref>, the HD cell network and the conjunctive grid cell network integrate velocity information and visual information to form neural codes of the pose of the robot. The spatial memory network node also makes decisions about the creation of nodes and links in the experience map and sends ROS messages of graph operations to the experience map node.</p>
<p>The experience map node reads out the neural codes of the conjunctive grid units and represents the key locations of the environment as the vertices in a topological graph. A vertex stores both the position estimated from the conjunctive grid units and the odometry relative to its previous vertex in the graph. On loop closure, a topological map relaxation method is used to find the minimum disagreement between the odometric transition information and the absolute position in the topological map space. When the position encoded in the conjunctive grid cell network is far enough from the position of the previous vertex, the spatial memory network node would inform the experience map to create a new vertex and a new edge linking to the previous vertex.</p>
<p>We write python scripts to visualize the live state of our robot navigation system. The key components of the system are shown, including the real time neural activity of the HD cells and the conjunctive grid cells (Figures <xref ref-type="fig" rid="F4">4</xref>A,B), the image of the scene and the local view templates (Figure <xref ref-type="fig" rid="F4">4</xref>C), as well as the current experience map (Figure <xref ref-type="fig" rid="F4">4</xref>D).</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Screenshots of the cognitive mapping model. <bold>(A)</bold> The neural activity of the HD-by-velocity units is depicted as a heat map, with warm colors for high activities and cool colors for acitivities diminishing to zero; <bold>(B)</bold> the activity of the conjunctive grid-by-velocity units is shown as heat maps, and each panel corresponds to a subspace of the neural manifold with the same velocity label. The same color code as in <bold>(A)</bold> is used; <bold>(C)</bold> the video frame (top) is the visual input to the robot. A visual feature map is constructed to represent the current scene (bottom). The best matched template is shown in the middle; <bold>(D)</bold> experience map showing paths with nodes and links.</p></caption>
<graphic xlink:href="fnbot-11-00061-g004.tif"/>
</fig>
</sec>
<sec id="S2-9">
<label>2.9</label> <title>St Lucia Suburb Dataset</title>
<p>The St Lucia Suburb dataset was first used in Milford and Wyeth (<xref ref-type="bibr" rid="B28">2008</xref>). The dataset was gathered in the suburb area of St Lucia in Brisbane, Australia by a vehicle. The vehicle was driven around the road network, visiting every street at least once. The trajectory is 66&#x02009;km long, finished within about 100&#x02009;min at typical driving speed up to 60&#x02009;km/h. The trajectory spans an area of 3&#x02009;km by 1.6&#x02009;km in East&#x02013;West and North&#x02013;South direction, respectively. The street views were recorded in 10&#x02009;frames/s into a video of resolution 640&#x02009;&#x000D7;&#x02009;480 pixels by a laptop&#x02019;s built-in web-cam mounted on the roof of the vehicle. The GPS information of the vehicle was not collected, however.</p>
<p>St Lucia is a very challenging environment to map due to its big variability and fairly large scale. There are 51 inner loops of varying size and shape, and more than 80 intersections in the road network (Figure <xref ref-type="fig" rid="F5">5</xref>A). The roads vary from busy multilane streets to small tracks, the terrains range from flat plains to steep hills, and the light conditions change a lot as well (Figure <xref ref-type="fig" rid="F5">5</xref>B).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>The complex environment of St Lucia. <bold>(A)</bold> The map of St Lucia from a bird&#x02019;s eye view. Black lines are the major road network in St Lucia urban area. The area is approximately 3.0&#x02009;km by 1.6&#x02009;km in size; <bold>(B)</bold> example scenes experienced by the robot in St Lucia. Panel <bold>(A)</bold> reused from <uri xlink:href="http://www.itee.uq.edu.au/think//filething/get/246/MichaelFig3.png">http://www.itee.uq.edu.au/think//filething/get/246/MichaelFig3.png</uri>.</p></caption>
<graphic xlink:href="fnbot-11-00061-g005.tif"/>
</fig>
</sec>
</sec>
<sec id="S3">
<label>3</label> <title>Results</title>
<p>We tested our cognitive mapping model on the St Lucia dataset (Ball et al., <xref ref-type="bibr" rid="B2">2013</xref>). Our cognitive mapping model is run on a personal computer with 3.4&#x02009;GHz six-core Intel i7 processor and 64 GB memory. The program is paralleled using OpenMP. Video <xref ref-type="supplementary-material" rid="SM1">S1</xref> in Supplementary Material shows the mapping process.</p>
<sec id="S3-10">
<label>3.1</label> <title>Neural Representation</title>
<p>Figures <xref ref-type="fig" rid="F6">6</xref>A,C show the activity of the HD-by-velocity units and the conjunctive grid-by-velocity units in the beginning of the experiment, when the robot is stationary. In the HD-by-velocity network, a single bump of activity emerges. Since the angular velocity is zero at this time, the bump is centered in the velocity dimension (red arrow in Figure <xref ref-type="fig" rid="F6">6</xref>A). The position of the bump in the direction dimension (black arrow in Figure <xref ref-type="fig" rid="F6">6</xref>A) encodes the direction of the robot in the environment. At this moment, the activity of the grid-by-velocity units is also centered in the center of the velocity dimensions (red arrows in Figure <xref ref-type="fig" rid="F6">6</xref>C) due to the zero translational velocity input. In spatial dimensions, the same grid pattern sustains, with peak activity decays gradually when the velocity labels of the units are away from the center. The phase of the grid pattern (short black arrows in Figure <xref ref-type="fig" rid="F6">6</xref>C) encodes the position of the robot in the environment.</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Population activities of the units in the cognitive mapping model. The activity of each unit is color-coded from red to blue for the maximal activity to zero activity. <bold>(A)</bold> The population activity of the HD-by-velocity units in the beginning of the experiment. The horizontal axis is the direction axis <italic>&#x003B8;</italic>, and the vertical axis is the velocity axis <italic>&#x003BD;</italic>. Positive (negative) <italic>&#x003BD;</italic> represent counterclockwise (clockwise) rotation. A localized bump pattern emerges in the network. At his moment, the center of the bump is at (0, 0), indicated by the black and red arrows in the direction and velocity dimension, respectively; <bold>(B)</bold> the activity of the HD-by-velocity network in the middle of the experiment. The center of the bump is at (3.25, 0.0069); <bold>(C)</bold> the population activity of the grid-by-velocity units in the beginning of the experiment. The four-dimensional neural activity is sliced in the <italic>&#x003BD;<sub>x</sub></italic> and <italic>&#x003BD;<sub>y</sub></italic> dimensions into 49 planes shown by each panel, with increasing <italic>&#x003BD;<sub>x</sub></italic> from left to right and increasing <italic>&#x003BD;<sub>y</sub></italic> from bottom to top. Each panel depicts the activity of the units with the same velocity label, with horizontal axis for <italic>&#x003B8;<sub>x</sub></italic> and vertical axis for <italic>&#x003B8;<sub>y</sub></italic>. The center of the bumps is marked by the short black arrows in the spatial dimensions and short red arrows in the velocity dimensions; <bold>(D)</bold> the activity of the grid-by-velocity network in the middle of the experiment. The center of the bumps is at (3.85, 3.91, &#x02212;0.0889, &#x02212;0.0045) in the 4D neural space (<italic>&#x003B8;<sub>x</sub></italic>, <italic>&#x003B8;<sub>y</sub></italic>, <italic>&#x003BD;<sub>x</sub></italic>, <italic>&#x003BD;<sub>y</sub></italic>).</p></caption>
<graphic xlink:href="fnbot-11-00061-g006.tif"/>
</fig>
<p>One example of the network states in the middle of the experiment is displayed in Figures <xref ref-type="fig" rid="F6">6</xref>B,D. The activity bump of the HD-by-velocity network is centered at (3.25, 0.0069) (Figure <xref ref-type="fig" rid="F6">6</xref>B), which means that the current HD angle of the robot is encoded as 186.21&#x000B0;, and meanwhile, the robot is rotating counterclockwise at 0.69&#x02009;rad/s in the physical environment. At this moment, the activity pattern of the grid-by-velocity network is centered at <italic>&#x003BD;<sub>x</sub></italic>&#x02009;&#x0003D;&#x02009;&#x02212;&#x02009;0.0889 and <italic>&#x003BD;<sub>y</sub></italic>&#x02009;&#x0003D;&#x02009;&#x02212;&#x02009;0.0045 in the velocity dimensions, corresponding to the movement of the robot in <italic>V<sub>x</sub></italic>&#x02009;&#x0003D;&#x02009;&#x02212;&#x02009;11.43 and <italic>V<sub>y</sub></italic>&#x02009;&#x0003D;&#x02009;&#x02212;&#x02009;0.57&#x02009;m/s.</p>
</sec>
<sec id="S3-11">
<label>3.2</label> <title>Pattern Completion with Visual Feedback</title>
<p>The network states are calibrated by the cues provided by the local view cells. The calibration is achieved by the pattern completion mechanism of the attractor networks. Figure <xref ref-type="fig" rid="F7">7</xref> and Video <xref ref-type="supplementary-material" rid="SM2">S2</xref> in Supplementary Materials show an example of pattern completion process during loop closure. In Figures <xref ref-type="fig" rid="F7">7</xref>A,B, the robot traverses a path for the first time, while the networks integrate velocity inputs, forming localized bumps that encode the head direction and the position of the robot. When the robot closes a loop and revisits familiar places (Figures <xref ref-type="fig" rid="F7">7</xref>C,D), local view cells get activated and provide strong cues to the networks. Weak bumps emerge at the phases associated to the local view cells and quickly become the dominating bumps in the network, retrieving phases similar to the phases when the robots first visited this place (Figures <xref ref-type="fig" rid="F7">7</xref>A vs. <xref ref-type="fig" rid="F7">7</xref>D). Once the robot moves to unfamiliar areas, local view cells become inactive, and the bumps revert to their localized shapes (Figures <xref ref-type="fig" rid="F7">7</xref>E,F).</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Calibration of the network states during loop closure by the help of visual feedback from the local view cells. There are three parts in each panel: the activity of the HD-by-velocity network (top), the activity of the grid-by-velocity network (middle), and the corresponding experience map at the time (bottom). <bold>(A)</bold> The robot drives along a path for the first time. The activity bump of the HD-by-velocity network is centered at (3.3763, 0.0023) and the bumps of the grid-by-velocity network are at (1.1781, 5.7298, &#x02212;0.1027, &#x02212;0.0191); <bold>(B)</bold> before the robot revisits a familiar place, the center of the bump of the HD-by-velocity network moves to (3.0973, 0.0021) and the center of the bumps of the grid-by-velocity network (bottom) is at (2.2176, 2.3562, &#x02212;0.1272, &#x02212;0.0158). <bold>(C)</bold> The robot revisits the same place as in <bold>(A)</bold> and starts to close a loop. The local view cells inject currents to the HD-by-velocity network and the grid-by-velocity network, and create weaker bumps at the same phases as the network states in <bold>(A)</bold>; <bold>(D)</bold> after multiple continuous current injections, the networks retrieve similar phases as those when the robot first visited the place; <bold>(E)</bold> once the robot enters a new area, the local view cells become inactive, and the network states stat to sharpen. <bold>(F)</bold> After loop closures, the network states are modulated by the velocity inputs. The bumps revert to their localized shapes.</p></caption>
<graphic xlink:href="fnbot-11-00061-g007.tif"/>
</fig>
</sec>
<sec id="S3-12">
<label>3.3</label> <title>Cognitive Map</title>
<p>The cognitive map created by the mapping system is shown in Figure <xref ref-type="fig" rid="F8">8</xref>. The thick green line is composed of the vertices of the topological graph. The positions of the vertices are determined by the conjunctive grid cells and graph relaxation. The link between two related vertices is shown by the fine blue line. Since the physical distance between two topological vertices is considered, the experience map is actually a semi-metric topological map, which can be compared visually with the ground truth map (Figure <xref ref-type="fig" rid="F5">5</xref>A) by naked eyes. The cognitive map conserves the overall layout of the road network of the environment. The cognitive map correctly represents all loop closures and intersections, although the orientation and length of the path are slightly different from the ground truth map. Overall, the map built by our method is consistent with the true map of the environment, and of similar quality as compared with that obtained in Ball et al. (<xref ref-type="bibr" rid="B2">2013</xref>).</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>The semi-metric cognitive map of St Lucia formed in the cognitive mapping model. The green thick line comprises topological graph vertices, and the blue thin line consists of links between connected vertices.</p></caption>
<graphic xlink:href="fnbot-11-00061-g008.tif"/>
</fig>
</sec>
<sec id="S3-13">
<label>3.4</label> <title>Firing Rate Maps</title>
<p>To see the responses of single units, we show the activities of example units on top of the cognitive map. Figure <xref ref-type="fig" rid="F9">9</xref>A depicts the firing rate map of the HD-by-velocity unit with label (0, 0) in the neural manifold. It fires strongly when the robot moves in broad directions close to southwest. The unit responses to a wide range of directions, since the width of the bump in the HD-by-velocity network is quite big, covering about 70&#x000B0; in head directions. The unit does not fire on many of the roads parallel to the southwest direction, since on those roads the robot moves northeast, i.e., opposite to the preferred firing direction of the unit. This can be confirmed by the firing rate map of a unit with opposite preferred head direction, shown in Figure <xref ref-type="fig" rid="F9">9</xref>B. This unit has the label (&#x003C0;, 0.0024). This unit fires in broad directions close to northeast, and its firing rate map does not overlap with that of the unit in Figure <xref ref-type="fig" rid="F9">9</xref>A, meaning that it keeps silent when the unit in Figure <xref ref-type="fig" rid="F9">9</xref>A fires. This unit prefers a counterclockwise rotational velocity of 0.24&#x02009;rad/s. Therefore, on bending paths, when the robot turns counterclockwise at about northeast directions, its firing rate always increases first and then decreases, which means that the head direction and the angular velocity of the robot become close to and away from the preferred conjunction of direction and rotational velocity of this unit. Figure <xref ref-type="fig" rid="F9">9</xref>C shows the firing rate map of the HD units with label 0 in the dimension <italic>&#x003B8;</italic>. The firing pattern in Figure <xref ref-type="fig" rid="F9">9</xref>C indicates these units encode quite large range of directions centered at southwest. In local regions, these cells conserve consistent preference for directions.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p>Firing rate maps of example HD-by-velocity units. In each panel, the firing rate is plotted at the locations in the experience map. Firing rate is color-coded by the same jet colormap, i.e., red for high firing rate and blue for zero firing rate. <bold>(A)</bold> The firing rate map of the unit at (0, 0) in the (<italic>&#x003B8;, &#x003BD;</italic>) manifold; <bold>(B)</bold> The firing rate map of the unit at (&#x003C0;, 0.0024), which has the opposite preference in head direction as the unit in <bold>(A)</bold>; <bold>(C)</bold> The total activity of all the units with the same label <italic>&#x003B8;</italic>&#x02009;&#x0003D;&#x02009;0, i.e., summing firing rates of all the units along the dimension <italic>&#x003BD;</italic> at <italic>&#x003B8;</italic>&#x02009;&#x0003D;&#x02009;0.</p></caption>
<graphic xlink:href="fnbot-11-00061-g009.tif"/>
</fig>
<p>Due to the periodic boundary conditions of the grid pattern in the conjunctive grid-by-velocity cell network, each grid unit fires at multiple distinct locations in the environment. Figures <xref ref-type="fig" rid="F10">10</xref>A,B give the firing rate maps of two example grid units in the network with the same spatial label but different velocity labels. In each of these firing fields, the firing rate always gradually increases when the robot moves closer to the center of the firing field and decreases bit by bit when the robot leaves the field center. The firing fields of the two example grid units overlap a lot, since they have the same spatial preference. The difference in firing field comes from the fact that the grid unit in Figure <xref ref-type="fig" rid="F10">10</xref>B prefers high translational velocity, however, the grid unit in Figure <xref ref-type="fig" rid="F10">10</xref>A prefers low translational velocity. Therefore, the grid units in Figure <xref ref-type="fig" rid="F10">10</xref>A fires in many locations close to the turns of the road, where the robot has to slow down significantly. The total activity of the grid units, which share the same spatial preference as the units in Figures <xref ref-type="fig" rid="F10">10</xref>A,B is shown in Figure <xref ref-type="fig" rid="F10">10</xref>C. The firing map is composed of many distributed firing fields, but they are lack of global grid structure. This is similar to the fragmented firing pattern observed in hairpin maze (Derdikman et al., <xref ref-type="bibr" rid="B7">2009</xref>), where the animal also runs on linear paths. Note that the distance between the grid fields of the units in the network is much larger than the typical spacing observed in animal experiments. This is due to the fact that the velocity input to the model is defined in the visual space of the input images and is different from the actual velocity of the robot.</p>
<fig id="F10" position="float">
<label>Figure 10</label>
<caption><p>Firing rate maps of example conjunctive grid-by-velocity units. In each panel, firing rate is plotted at the locations in the experience map. For all panels, the same colormap is used to code firing rate, red for high firing rate and gradually changing to blue for firing rate decreasing to zero. <bold>(A)</bold> The firing rate map of the unit at (&#x003C0;, &#x003C0;, 0, 0) in the (<italic>&#x003B8;<sub>x</sub></italic>, <italic>&#x003B8;<sub>y</sub></italic>, <italic>&#x003BD;<sub>x</sub></italic>, <italic>&#x003BD;<sub>y</sub></italic>) manifold; <bold>(B)</bold> the firing rate map of the unit at (&#x003C0;, &#x003C0;, 0.1, 0.1); <bold>(C)</bold> the map of the total activity of the units at (&#x003C0;, &#x003C0;) in the (<italic>&#x003B8;<sub>x</sub></italic>, <italic>&#x003B8;<sub>y</sub></italic>) manifold, which sums the firing rates of the units along the dimensions <italic>&#x003BD;<sub>x</sub></italic> and <italic>&#x003BD;<sub>y</sub></italic>.</p></caption>
<graphic xlink:href="fnbot-11-00061-g010.tif"/>
</fig>
</sec>
</sec>
<sec id="S4" sec-type="discussion">
<label>4</label> <title>Discussion</title>
<p>We proposed a robust cognitive mapping model, which integrates movement and sensory information and forms stable spatial representations of the environment. We demonstrated that our cognitive mapping model using only visual sensory information could successfully build a coherent topological map of large-scale outdoor environment from an open-source dataset of 66&#x02009;km urban car journey (see Video <xref ref-type="supplementary-material" rid="SM1">S1</xref> in Supplementary Materials).</p>
<p>There are two major contributions in this study. First, the network architecture of our cognitive mapping model is consistent with the synaptic organization of the MEC. Our model is built upon layer III and deep layers of the MEC. Layer III and deep layers of the MEC may be the main site to perform key computations for navigation, and the layer II of the MEC may merely relays information to the hippocampus and functions as the output of the MEC. Supporting experimental evidence shows that there exist strong projections from deep layers to superficial layers of the MEC. HD cells, conjunctive grid-by-head-direction cells, and speed cells are abundant in layer III and deep layers of the MEC, and encode direction and speed information needed in path integration. However, layer II of the MEC contains mainly pure positional grid cells. Second, our cognitive mapping model provides coherent robust mechanisms to perform path integration. Both the HD cells and the conjunctive grid cells in the model integrate motion and perception information using the same principle. The performance of the model is manifested by experiments in large-scale natural environment. Compared with the pose cell network in Ball et al. (<xref ref-type="bibr" rid="B2">2013</xref>), our model achieves similar mapping results as those shown in Ball et al. (<xref ref-type="bibr" rid="B2">2013</xref>). Although our model reuses some of the vision and map operation methods of Ball et al. (<xref ref-type="bibr" rid="B2">2013</xref>), the core components of the HD cell network and the conjunctive grid cell network in our system are fundamentally different from the pose cell network. In our model, path integration is performed by the intrinsic attractor dynamics driven by velocity inputs, rather than simply displacing the copy of the activity package. In the proposed model, the convergence process of local view calibration is also performed by the dynamics of the attractors, instead of global inhibition.</p>
<p>The grid cells in our model does not express grid-like firing pattern in large-scale environment (Figure <xref ref-type="fig" rid="F10">10</xref>). Meanwhile, the HD cells in the model fire at a broad range of directions (Figure <xref ref-type="fig" rid="F9">9</xref>). This is due to the fact that local view cells anchor grids to local cues and alters the firing patterns of the units. It has been shown that in environments with abundant local cues, such like corridors, grids cells show fragmented firing maps (Derdikman et al., <xref ref-type="bibr" rid="B7">2009</xref>). It is much harder for animals to form a globally consistent metric map in large than in small environments. Indeed, as shown in Carpenter et al. (<xref ref-type="bibr" rid="B5">2015</xref>), grid cells in rats acquire globally coherent representations after repeated exploration, while initially grid maps are only locally consistent.</p>
<p>Several potential limitations exist in the current study. First, the number of units in the system is large and requires substantial amount of computational resources. Second, place cell network is not included in our system. The function of place cells is not to provide a metric map. Experiments showed that in large-scale environments, place cell has multiple irregularly spaced place fields (Rich et al., <xref ref-type="bibr" rid="B33">2014</xref>; Liu et al., <xref ref-type="bibr" rid="B24">2015</xref>). A metric map is necessary if direct interpretation is preferred. For robotic applications, it would suffice to have distributed neural representations if the robots are able to distinguish different locations of the environment.</p>
<p>For future research, we plan to include a place cell network, which integrates multiple grid cell modules with diversity of spacings and orientations, and investigate its function in path planning and reward-based learning. We will also investigate neurobiologically inspired algorithms to improve the computational efficiency of the system for robot navigation.</p>
</sec>
<sec id="S5">
<label>5</label> <title>Conclusion</title>
<p>In summary, conjunctive space-by-movement attractor network models are proposed in this paper to achieve cognitive mapping in a large-scale natural environment. Head direction cells and conjunctive grid cells work on the same principle, and represent positions, head directions, and velocity at the same time. Our model facilitates to reach a better understanding of the neural mechanisms of spatial cognition and to develop more accurate models of the brain. Furthermore, it inspires innovative high performance cognitive mapping systems, which are able to function in dynamical natural environment with long-term autonomy.</p>
</sec>
<sec id="S6" sec-type="author-contributor">
<title>Author Contributions</title>
<p>TZ and BS conceived the model, performed experiments, analyzed the results, and wrote the paper.</p>
</sec>
<sec id="S9">
<title>Conflict of Interest Statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<fn-group>
<fn fn-type="financial-disclosure">
<p><bold>Funding.</bold> This work is supported by the Distinguished Young Scholar Project of the Thousand Talents Program of China (grant No. Y5A1370101) and the State Key Laboratory of Robotics of China (grant No. Y5A1203901).</p>
</fn>
</fn-group>
<sec id="S8" sec-type="supplementary-material">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at <uri xlink:href="http://www.frontiersin.org/article/10.3389/fnbot.2017.00061/full#supplementary-material">http://www.frontiersin.org/article/10.3389/fnbot.2017.00061/full&#x00023;supplementary-material</uri>.</p>
<supplementary-material xlink:href="Video_1.MP4" id="SM1" mimetype="applicationn/MP4" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Video S1</label>
<caption><p>Mapping process.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Video_2.MP4" id="SM2" mimetype="applicationn/MP4" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Video S2</label>
<caption><p>An example of pattern completion process.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Presentation_1.PDF" id="SM3" mimetype="applicationn/PDF" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bailey</surname> <given-names>T.</given-names></name> <name><surname>Durrant-Whyte</surname> <given-names>H.</given-names></name></person-group> (<year>2006</year>). <article-title>Simultaneous localization and mapping (SLAM): part II</article-title>. <source>IEEE Robot. Autom. Mag.</source> <volume>13</volume>, <fpage>108</fpage>&#x02013;<lpage>117</lpage>.<pub-id pub-id-type="doi">10.1109/MRA.2006.1678144</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ball</surname> <given-names>D.</given-names></name> <name><surname>Heath</surname> <given-names>S.</given-names></name> <name><surname>Wiles</surname> <given-names>J.</given-names></name> <name><surname>Wyeth</surname> <given-names>G.</given-names></name> <name><surname>Corke</surname> <given-names>P.</given-names></name> <name><surname>Milford</surname> <given-names>M.</given-names></name></person-group> (<year>2013</year>). <article-title>OpenRatSLAM: an open source brain-based SLAM system</article-title>. <source>Auton. Robots</source> <volume>34</volume>, <fpage>149</fpage>&#x02013;<lpage>176</lpage>.<pub-id pub-id-type="doi">10.1007/s10514-012-9317-9</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Burak</surname> <given-names>Y.</given-names></name> <name><surname>Fiete</surname> <given-names>I. R.</given-names></name></person-group> (<year>2009</year>). <article-title>Accurate path integration in continuous attractor network models of grid cells</article-title>. <source>PLoS Comput. Biol.</source> <volume>5</volume>:<fpage>e1000291</fpage>.<pub-id pub-id-type="doi">10.1371/journal.pcbi.1000291</pub-id><pub-id pub-id-type="pmid">19229307</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Burgess</surname> <given-names>N.</given-names></name> <name><surname>Barry</surname> <given-names>C.</given-names></name> <name><surname>O&#x02019;Keefe</surname> <given-names>J.</given-names></name></person-group> (<year>2007</year>). <article-title>An oscillatory interference model of grid cell firing</article-title>. <source>Hippocampus</source> <volume>17</volume>, <fpage>801</fpage>&#x02013;<lpage>812</lpage>.<pub-id pub-id-type="doi">10.1002/hipo.20327</pub-id><pub-id pub-id-type="pmid">17598147</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Carpenter</surname> <given-names>F.</given-names></name> <name><surname>Manson</surname> <given-names>D.</given-names></name> <name><surname>Jeffery</surname> <given-names>K.</given-names></name> <name><surname>Burgess</surname> <given-names>N.</given-names></name> <name><surname>Barry</surname> <given-names>C.</given-names></name></person-group> (<year>2015</year>). <article-title>Grid cells form a global representation of connected environments</article-title>. <source>Curr. Biol.</source> <volume>25</volume>, <fpage>1176</fpage>&#x02013;<lpage>1182</lpage>.<pub-id pub-id-type="doi">10.1016/j.cub.2015.02.037</pub-id><pub-id pub-id-type="pmid">25913404</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cheng</surname> <given-names>S.</given-names></name> <name><surname>Frank</surname> <given-names>L. M.</given-names></name></person-group> (<year>2011</year>). <article-title>The structure of networks that produce the transformation from grid cells to place cells</article-title>. <source>Neuroscience</source> <volume>197</volume>, <fpage>293</fpage>&#x02013;<lpage>306</lpage>.<pub-id pub-id-type="doi">10.1016/j.neuroscience.2011.09.002</pub-id><pub-id pub-id-type="pmid">21963867</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Derdikman</surname> <given-names>D.</given-names></name> <name><surname>Whitlock</surname> <given-names>J. R.</given-names></name> <name><surname>Tsao</surname> <given-names>A.</given-names></name> <name><surname>Fyhn</surname> <given-names>M.</given-names></name> <name><surname>Hafting</surname> <given-names>T.</given-names></name> <name><surname>Moser</surname> <given-names>M.-B.</given-names></name> <etal/></person-group> (<year>2009</year>). <article-title>Fragmentation of grid cell maps in a multicompartment environment</article-title>. <source>Nat. Neurosci.</source> <volume>12</volume>, <fpage>1325</fpage>&#x02013;<lpage>1332</lpage>.<pub-id pub-id-type="doi">10.1038/nn.2396</pub-id><pub-id pub-id-type="pmid">19749749</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ding</surname> <given-names>L.</given-names></name> <name><surname>Xiao</surname> <given-names>L.</given-names></name> <name><surname>Liao</surname> <given-names>B.</given-names></name> <name><surname>Lu</surname> <given-names>R.</given-names></name> <name><surname>Peng</surname> <given-names>H.</given-names></name></person-group> (<year>2017</year>). <article-title>An improved recurrent neural network for complex-valued systems of linear equation and its application to robotic motion tracking</article-title>. <source>Front. Neurorobot.</source> <volume>11</volume>:<fpage>45</fpage>.<pub-id pub-id-type="doi">10.3389/fnbot.2017.00045</pub-id><pub-id pub-id-type="pmid">28919855</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Duckett</surname> <given-names>T.</given-names></name> <name><surname>Marsland</surname> <given-names>S.</given-names></name> <name><surname>Shapiro</surname> <given-names>J.</given-names></name></person-group> (<year>2002</year>). <article-title>Fast, on-line learning of globally consistent maps</article-title>. <source>Auton. Robots</source> <volume>12</volume>, <fpage>287</fpage>&#x02013;<lpage>300</lpage>.<pub-id pub-id-type="doi">10.1023/A:1015269615729</pub-id></citation></ref>
<ref id="B10"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Dudek</surname> <given-names>G.</given-names></name> <name><surname>Jenkin</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <source>Computational Principles of Mobile Robotics</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>.</citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Durrant-Whyte</surname> <given-names>H.</given-names></name> <name><surname>Bailey</surname> <given-names>T.</given-names></name></person-group> (<year>2006</year>). <article-title>Simultaneous localization and mapping: part I</article-title>. <source>IEEE Robot. Autom. Mag.</source> <volume>13</volume>, <fpage>99</fpage>&#x02013;<lpage>110</lpage>.<pub-id pub-id-type="doi">10.1109/MRA.2006.1638022</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Erlhagen</surname> <given-names>W.</given-names></name> <name><surname>Bicho</surname> <given-names>E.</given-names></name></person-group> (<year>2006</year>). <article-title>The dynamic neural field approach to cognitive robotics</article-title>. <source>J. Neural Eng.</source> <volume>3</volume>, <fpage>R36</fpage>.<pub-id pub-id-type="doi">10.1088/1741-2560/3/3/R02</pub-id><pub-id pub-id-type="pmid">16921201</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fuhs</surname> <given-names>M. C.</given-names></name> <name><surname>Touretzky</surname> <given-names>D. S.</given-names></name></person-group> (<year>2006</year>). <article-title>A spin glass model of path integration in rat medial entorhinal cortex</article-title>. <source>J. Neurosci.</source> <volume>26</volume>, <fpage>4266</fpage>&#x02013;<lpage>4276</lpage>.<pub-id pub-id-type="doi">10.1523/JNEUROSCI.4353-05.2006</pub-id><pub-id pub-id-type="pmid">16624947</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fyhn</surname> <given-names>M.</given-names></name> <name><surname>Hafting</surname> <given-names>T.</given-names></name> <name><surname>Treves</surname> <given-names>A.</given-names></name> <name><surname>Moser</surname> <given-names>M.-B.</given-names></name> <name><surname>Moser</surname> <given-names>E.</given-names></name></person-group> (<year>2007</year>). <article-title>Hippocampal remapping and grid realignment in entorhinal cortex</article-title>. <source>Nature</source> <volume>446</volume>, <fpage>190</fpage>&#x02013;<lpage>194</lpage>.<pub-id pub-id-type="doi">10.1038/nature05601</pub-id><pub-id pub-id-type="pmid">17322902</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gerstner</surname> <given-names>W.</given-names></name> <name><surname>Abbott</surname> <given-names>L.</given-names></name></person-group> (<year>1997</year>). <article-title>Learning navigational maps through potentiation and modulation of hippocampal place cells</article-title>. <source>J. Comput. Neurosci.</source> <volume>4</volume>, <fpage>79</fpage>&#x02013;<lpage>94</lpage>.<pub-id pub-id-type="doi">10.1023/A:1008820728122</pub-id><pub-id pub-id-type="pmid">9046453</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hafting</surname> <given-names>T.</given-names></name> <name><surname>Fyhn</surname> <given-names>M.</given-names></name> <name><surname>Molden</surname> <given-names>S.</given-names></name> <name><surname>Moser</surname> <given-names>M.-B.</given-names></name> <name><surname>Moser</surname> <given-names>E. I.</given-names></name></person-group> (<year>2005</year>). <article-title>Microstructure of a spatial map in the entorhinal cortex</article-title>. <source>Nature</source> <volume>436</volume>, <fpage>801</fpage>&#x02013;<lpage>806</lpage>.<pub-id pub-id-type="doi">10.1038/nature03721</pub-id><pub-id pub-id-type="pmid">15965463</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jauffret</surname> <given-names>A.</given-names></name> <name><surname>Cuperlier</surname> <given-names>N.</given-names></name> <name><surname>Gaussier</surname> <given-names>P.</given-names></name></person-group> (<year>2015</year>). <article-title>From grid cells and visual place cells to multimodal place cell: a new robotic architecture</article-title>. <source>Front. Neurorobot.</source> <volume>9</volume>:<fpage>1</fpage>.<pub-id pub-id-type="doi">10.3389/fnbot.2015.00001</pub-id><pub-id pub-id-type="pmid">25904862</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Klaassen</surname> <given-names>R. H.</given-names></name> <name><surname>Alerstam</surname> <given-names>T.</given-names></name> <name><surname>Carlsson</surname> <given-names>P.</given-names></name> <name><surname>Fox</surname> <given-names>J. W.</given-names></name> <name><surname>Lindstr&#x000F6;m</surname> <given-names>&#x000C5;</given-names></name></person-group> (<year>2011</year>). <article-title>Great flights by great snipes: long and fast non-stop migration over benign habitats</article-title>. <source>Biol. Lett.</source> <volume>7</volume>, <fpage>833</fpage>&#x02013;<lpage>835</lpage>.<pub-id pub-id-type="doi">10.1098/rsbl.2011.0343</pub-id><pub-id pub-id-type="pmid">21613283</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Knierim</surname> <given-names>J. J.</given-names></name></person-group> (<year>2015</year>). <article-title>From the gps to hm: place cells, grid cells, and memory</article-title>. <source>Hippocampus</source> <volume>25</volume>, <fpage>719</fpage>&#x02013;<lpage>725</lpage>.<pub-id pub-id-type="doi">10.1002/hipo.22453</pub-id><pub-id pub-id-type="pmid">25788454</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Knierim</surname> <given-names>J. J.</given-names></name> <name><surname>Zhang</surname> <given-names>K.</given-names></name></person-group> (<year>2012</year>). <article-title>Attractor dynamics of spatially correlated neural activity in the limbic system</article-title>. <source>Annu. Rev. Neurosci.</source> <volume>35</volume>, <fpage>267</fpage>&#x02013;<lpage>285</lpage>.<pub-id pub-id-type="doi">10.1146/annurev-neuro-062111-150351</pub-id><pub-id pub-id-type="pmid">22462545</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Knips</surname> <given-names>G.</given-names></name> <name><surname>Zibner</surname> <given-names>S. K.</given-names></name> <name><surname>Reimann</surname> <given-names>H.</given-names></name> <name><surname>Sch&#x000F6;ner</surname> <given-names>G.</given-names></name></person-group> (<year>2017</year>). <article-title>A neural dynamic architecture for reaching and grasping integrates perception and movement generation and enables on-line updating</article-title>. <source>Front. Neurorobot.</source> <volume>11</volume>:<fpage>9</fpage>.<pub-id pub-id-type="doi">10.3389/fnbot.2017.00009</pub-id><pub-id pub-id-type="pmid">28303100</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kropff</surname> <given-names>E.</given-names></name> <name><surname>Carmichael</surname> <given-names>J. E.</given-names></name> <name><surname>Moser</surname> <given-names>M.-B.</given-names></name> <name><surname>Moser</surname> <given-names>E. I.</given-names></name></person-group> (<year>2015</year>). <article-title>Speed cells in the medial entorhinal cortex</article-title>. <source>Nature</source> <volume>523</volume>, <fpage>419</fpage>&#x02013;<lpage>424</lpage>.<pub-id pub-id-type="doi">10.1038/nature14622</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kropff</surname> <given-names>E.</given-names></name> <name><surname>Treves</surname> <given-names>A.</given-names></name></person-group> (<year>2008</year>). <article-title>The emergence of grid cells: intelligent design or just adaptation?</article-title> <source>Hippocampus</source> <volume>18</volume>, <fpage>1256</fpage>&#x02013;<lpage>1269</lpage>.<pub-id pub-id-type="doi">10.1002/hipo.20520</pub-id><pub-id pub-id-type="pmid">19021261</pub-id></citation></ref>
<ref id="B24"><citation citation-type="confproc"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>S.</given-names></name> <name><surname>Si</surname> <given-names>B.</given-names></name> <name><surname>Lin</surname> <given-names>Y.</given-names></name></person-group> (<year>2015</year>). <article-title>&#x0201C;Self-organization of hippocampal representations in large environments,&#x0201D;</article-title> in <conf-name>Proceedings of the 2015 International Joint Conference on Neural Networks</conf-name>, <conf-loc>Killarney</conf-loc>.</citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mathis</surname> <given-names>A.</given-names></name> <name><surname>Herz</surname> <given-names>A. V.</given-names></name> <name><surname>Stemmler</surname> <given-names>M.</given-names></name></person-group> (<year>2012</year>). <article-title>Optimal population codes for space: grid cells outperform place cells</article-title>. <source>Neural Comput.</source> <volume>24</volume>, <fpage>2280</fpage>&#x02013;<lpage>2317</lpage>.<pub-id pub-id-type="doi">10.1162/NECO_a_00319</pub-id><pub-id pub-id-type="pmid">22594833</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>McNaughton</surname> <given-names>B. L.</given-names></name> <name><surname>Battaglia</surname> <given-names>F. P.</given-names></name> <name><surname>Jensen</surname> <given-names>O.</given-names></name> <name><surname>Moser</surname> <given-names>E. I.</given-names></name> <name><surname>Moser</surname> <given-names>M.-B.</given-names></name></person-group> (<year>2006</year>). <article-title>Path integration and the neural basis of the &#x02018;cognitive map&#x02019;</article-title>. <source>Nat. Rev. Neurosci.</source> <volume>7</volume>, <fpage>663</fpage>&#x02013;<lpage>678</lpage>.<pub-id pub-id-type="doi">10.1038/nrn1932</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mi</surname> <given-names>Y.</given-names></name> <name><surname>Katkov</surname> <given-names>M.</given-names></name> <name><surname>Tsodyks</surname> <given-names>M.</given-names></name></person-group> (<year>2017</year>). <article-title>Synaptic correlates of working memory capacity</article-title>. <source>Neuron</source> <volume>93</volume>, <fpage>323</fpage>&#x02013;<lpage>330</lpage>.<pub-id pub-id-type="doi">10.1016/j.neuron.2016.12.004</pub-id><pub-id pub-id-type="pmid">28041884</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Milford</surname> <given-names>M. J.</given-names></name> <name><surname>Wyeth</surname> <given-names>G. F.</given-names></name></person-group> (<year>2008</year>). <article-title>Mapping a suburb with a single camera using a biologically inspired SLAM system</article-title>. <source>IEEE Trans. Robot.</source> <volume>24</volume>, <fpage>1038</fpage>&#x02013;<lpage>1053</lpage>.<pub-id pub-id-type="doi">10.1109/TRO.2008.2004520</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mulas</surname> <given-names>M.</given-names></name> <name><surname>Waniek</surname> <given-names>N.</given-names></name> <name><surname>Conradt</surname> <given-names>J.</given-names></name></person-group> (<year>2016</year>). <article-title>Hebbian plasticity realigns grid cell activity with external sensory cues in continuous attractor models</article-title>. <source>Front. Comput. Neurosci.</source> <volume>10</volume>:<fpage>13</fpage>.<pub-id pub-id-type="doi">10.3389/fncom.2016.00013</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>O&#x02019;Keefe</surname> <given-names>J.</given-names></name> <name><surname>Conway</surname> <given-names>D.</given-names></name></person-group> (<year>1978</year>). <article-title>Hippocampal place units in the freely moving rat: why they fire where they fire</article-title>. <source>Exp. Brain Res.</source> <volume>31</volume>, <fpage>573</fpage>&#x02013;<lpage>590</lpage>.<pub-id pub-id-type="pmid">658182</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>O&#x02019;Keefe</surname> <given-names>J.</given-names></name> <name><surname>Dostrovsky</surname> <given-names>J.</given-names></name></person-group> (<year>1971</year>). <article-title>The hippocampus as a spatial map. Preliminary evidence from unit activity in the freely-moving rat</article-title>. <source>Brain Res.</source> <volume>34</volume>, <fpage>171</fpage>&#x02013;<lpage>175</lpage>.<pub-id pub-id-type="doi">10.1016/0006-8993(71)90358-1</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Remme</surname> <given-names>M. W.</given-names></name> <name><surname>Lengyel</surname> <given-names>M.</given-names></name> <name><surname>Gutkin</surname> <given-names>B. S.</given-names></name></person-group> (<year>2010</year>). <article-title>Democracy-independence trade-off in oscillating dendrites and its implications for grid cells</article-title>. <source>Neuron</source> <volume>66</volume>, <fpage>429</fpage>&#x02013;<lpage>437</lpage>.<pub-id pub-id-type="doi">10.1016/j.neuron.2010.04.027</pub-id><pub-id pub-id-type="pmid">20471355</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rich</surname> <given-names>P. D.</given-names></name> <name><surname>Liaw</surname> <given-names>H.-P.</given-names></name> <name><surname>Lee</surname> <given-names>A. K.</given-names></name></person-group> (<year>2014</year>). <article-title>Large environments reveal the statistical structure governing hippocampal representations</article-title>. <source>Science</source> <volume>345</volume>, <fpage>814</fpage>&#x02013;<lpage>817</lpage>.<pub-id pub-id-type="doi">10.1126/science.1255635</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sargolini</surname> <given-names>F.</given-names></name> <name><surname>Fyhn</surname> <given-names>M.</given-names></name> <name><surname>Hafting</surname> <given-names>T.</given-names></name> <name><surname>McNaughton</surname> <given-names>B. L.</given-names></name> <name><surname>Witter</surname> <given-names>M. P.</given-names></name> <name><surname>Moser</surname> <given-names>M.-B.</given-names></name> <etal/></person-group> (<year>2006</year>). <article-title>Conjunctive representation of position, direction, and velocity in entorhinal cortex</article-title>. <source>Science</source> <volume>312</volume>, <fpage>758</fpage>&#x02013;<lpage>762</lpage>.<pub-id pub-id-type="doi">10.1126/science.1125572</pub-id><pub-id pub-id-type="pmid">16675704</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sheynikhovich</surname> <given-names>D.</given-names></name> <name><surname>Chavarriaga</surname> <given-names>R.</given-names></name> <name><surname>Str&#x000F6;sslin</surname> <given-names>T.</given-names></name> <name><surname>Arleo</surname> <given-names>A.</given-names></name> <name><surname>Gerstner</surname> <given-names>W.</given-names></name></person-group> (<year>2009</year>). <article-title>Is there a geometric module for spatial orientation? Insights from a rodent navigation model</article-title>. <source>Psychol. Rev.</source> <volume>116</volume>, <fpage>540</fpage>&#x02013;<lpage>566</lpage>.<pub-id pub-id-type="doi">10.1037/a0016170</pub-id><pub-id pub-id-type="pmid">19618986</pub-id></citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Si</surname> <given-names>B.</given-names></name> <name><surname>Kropff</surname> <given-names>E.</given-names></name> <name><surname>Treves</surname> <given-names>A.</given-names></name></person-group> (<year>2012</year>). <article-title>Grid alignment in entorhinal cortex</article-title>. <source>Biol. Cybern.</source> <volume>106</volume>, <fpage>483</fpage>&#x02013;<lpage>506</lpage>.<pub-id pub-id-type="doi">10.1007/s00422-012-0513-7</pub-id><pub-id pub-id-type="pmid">22892761</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Si</surname> <given-names>B.</given-names></name> <name><surname>Romani</surname> <given-names>S.</given-names></name> <name><surname>Tsodyks</surname> <given-names>M.</given-names></name></person-group> (<year>2014</year>). <article-title>Continuous attractor network model for conjunctive position-by-velocity tuning of grid cells</article-title>. <source>PLoS Comput. Biol.</source> <volume>10</volume>:<fpage>e1003558</fpage>.<pub-id pub-id-type="doi">10.1371/journal.pcbi.1003558</pub-id><pub-id pub-id-type="pmid">24743341</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Si</surname> <given-names>B.</given-names></name> <name><surname>Treves</surname> <given-names>A.</given-names></name></person-group> (<year>2009</year>). <article-title>The role of competitive learning in the generation of DG field s from EC inputs</article-title>. <source>Cogn. Neurodyn.</source> <volume>3</volume>, <fpage>177</fpage>&#x02013;<lpage>187</lpage>.<pub-id pub-id-type="doi">10.1007/s11571-009-9079-z</pub-id><pub-id pub-id-type="pmid">19301148</pub-id></citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Si</surname> <given-names>B.</given-names></name> <name><surname>Treves</surname> <given-names>A.</given-names></name></person-group> (<year>2013</year>). <article-title>A model for the differentiation between grid and conjunctive units in medial entorhinal cortex</article-title>. <source>Hippocampus</source> <volume>23</volume>, <fpage>1410</fpage>&#x02013;<lpage>1424</lpage>.<pub-id pub-id-type="doi">10.1002/hipo.22194</pub-id><pub-id pub-id-type="pmid">23966345</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>R. C.</given-names></name> <name><surname>Cheeseman</surname> <given-names>P.</given-names></name></person-group> (<year>1986</year>). <article-title>On the representation and estimation of spatial uncertainty</article-title>. <source>Int. J. Rob. Res.</source> <volume>5</volume>, <fpage>56</fpage>&#x02013;<lpage>68</lpage>.<pub-id pub-id-type="doi">10.1177/027836498600500404</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Solstad</surname> <given-names>T.</given-names></name> <name><surname>Moser</surname> <given-names>E. I.</given-names></name> <name><surname>Einevoll</surname> <given-names>G. T.</given-names></name></person-group> (<year>2006</year>). <article-title>From grid cells to place cells: a mathematical model</article-title>. <source>Hippocampus</source> <volume>16</volume>, <fpage>1026</fpage>&#x02013;<lpage>1031</lpage>.<pub-id pub-id-type="doi">10.1002/hipo.20244</pub-id><pub-id pub-id-type="pmid">17094145</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stella</surname> <given-names>F.</given-names></name> <name><surname>Si</surname> <given-names>B.</given-names></name> <name><surname>Kropff</surname> <given-names>E.</given-names></name> <name><surname>Treves</surname> <given-names>A.</given-names></name></person-group> (<year>2013</year>). <article-title>Grid maps for spaceflight, anyone? They are for free!</article-title> <source>Behav. Brain Sci.</source> <volume>36</volume>, <fpage>566</fpage>&#x02013;<lpage>567</lpage>.<pub-id pub-id-type="doi">10.1017/S0140525X13000575</pub-id><pub-id pub-id-type="pmid">24103622</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stella</surname> <given-names>F.</given-names></name> <name><surname>Treves</surname> <given-names>A.</given-names></name></person-group> (<year>2015</year>). <article-title>The self-organization of grid cells in 3d</article-title>. <source>Elife</source> <volume>4</volume>, <fpage>e05913</fpage>.<pub-id pub-id-type="doi">10.7554/eLife.05913</pub-id><pub-id pub-id-type="pmid">25821989</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stensola</surname> <given-names>T.</given-names></name> <name><surname>Stensola</surname> <given-names>H.</given-names></name> <name><surname>Moser</surname> <given-names>M.-B.</given-names></name> <name><surname>Moser</surname> <given-names>E. I.</given-names></name></person-group> (<year>2015</year>). <article-title>Shearing-induced asymmetry in entorhinal grid cells</article-title>. <source>Nature</source> <volume>518</volume>, <fpage>207</fpage>&#x02013;<lpage>212</lpage>.<pub-id pub-id-type="doi">10.1038/nature14151</pub-id><pub-id pub-id-type="pmid">25673414</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taube</surname> <given-names>J. S.</given-names></name></person-group> (<year>2007</year>). <article-title>The head direction signal: origins and sensory-motor integration</article-title>. <source>Annu. Rev. Neurosci.</source> <volume>30</volume>, <fpage>181</fpage>&#x02013;<lpage>207</lpage>.<pub-id pub-id-type="doi">10.1146/annurev.neuro.29.051605.112854</pub-id><pub-id pub-id-type="pmid">17341158</pub-id></citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taube</surname> <given-names>J. S.</given-names></name> <name><surname>Muller</surname> <given-names>R. U.</given-names></name> <name><surname>Ranck</surname> <given-names>J. B.</given-names></name></person-group> (<year>1990</year>). <article-title>Head-direction cells recorded from the postsubiculum in freely moving rats. I. Description and quantitative analysis</article-title>. <source>J. Neurosci.</source> <volume>10</volume>, <fpage>420</fpage>&#x02013;<lpage>435</lpage>.<pub-id pub-id-type="pmid">2303851</pub-id></citation></ref>
<ref id="B47"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Thrun</surname> <given-names>S.</given-names></name> <name><surname>Burgard</surname> <given-names>W.</given-names></name> <name><surname>Fox</surname> <given-names>D.</given-names></name></person-group> (<year>2005</year>). <source>Probabilistic Robotics (Intelligent Robotics and Autonomous Agents Series)</source>. <publisher-loc>Cambridge, MA</publisher-loc>: <publisher-name>The MIT Press</publisher-name>.</citation></ref>
<ref id="B48"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Thrun</surname> <given-names>S.</given-names></name> <name><surname>Leonard</surname> <given-names>J. J.</given-names></name></person-group> (<year>2008</year>). <article-title>&#x0201C;Simultaneous localization and mapping,&#x0201D;</article-title> in <source>Springer Handbook of Robotics</source>, eds <person-group person-group-type="editor"><name><surname>Jefferies</surname> <given-names>M. E.</given-names></name> <name><surname>Yeap</surname> <given-names>W.-K.</given-names></name></person-group> (<publisher-loc>Berlin, Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>871</fpage>&#x02013;<lpage>889</lpage>.</citation></ref>
<ref id="B49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tolman</surname> <given-names>E. C.</given-names></name></person-group> (<year>1948</year>). <article-title>Cognitive maps in rats and men</article-title>. <source>Psychol. Rev.</source> <volume>55</volume>, <fpage>189</fpage>.<pub-id pub-id-type="doi">10.1037/h0061626</pub-id></citation></ref>
<ref id="B50"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Treves</surname> <given-names>A.</given-names></name></person-group> (<year>1990</year>). <article-title>Graded-response neurons and information encodings in autoassociative memories</article-title>. <source>Phys. Rev. A</source> <volume>42</volume>, <fpage>2418</fpage>.<pub-id pub-id-type="doi">10.1103/PhysRevA.42.2418</pub-id></citation></ref>
<ref id="B51"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tsoar</surname> <given-names>A.</given-names></name> <name><surname>Nathan</surname> <given-names>R.</given-names></name> <name><surname>Bartan</surname> <given-names>Y.</given-names></name> <name><surname>Vyssotski</surname> <given-names>A.</given-names></name> <name><surname>DellOmo</surname> <given-names>G.</given-names></name> <name><surname>Ulanovsky</surname> <given-names>N.</given-names></name></person-group> (<year>2011</year>). <article-title>Large-scale navigational map in a mammal</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A.</source> <volume>108</volume>, <fpage>E718</fpage>&#x02013;<lpage>E724</lpage>.<pub-id pub-id-type="doi">10.1073/pnas.1107365108</pub-id><pub-id pub-id-type="pmid">21844350</pub-id></citation></ref>
<ref id="B52"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Welday</surname> <given-names>A. C.</given-names></name> <name><surname>Shlifer</surname> <given-names>I. G.</given-names></name> <name><surname>Bloom</surname> <given-names>M. L.</given-names></name> <name><surname>Zhang</surname> <given-names>K.</given-names></name> <name><surname>Blair</surname> <given-names>H. T.</given-names></name></person-group> (<year>2011</year>). <article-title>Cosine directional tuning of theta cell burst frequencies: evidence for spatial coding by oscillatory interference</article-title>. <source>J. Neurosci.</source> <volume>31</volume>, <fpage>16157</fpage>&#x02013;<lpage>16176</lpage>.<pub-id pub-id-type="doi">10.1523/JNEUROSCI.0712-11.2011</pub-id><pub-id pub-id-type="pmid">22072668</pub-id></citation></ref>
<ref id="B53"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiao</surname> <given-names>L.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name></person-group> (<year>2016</year>). <article-title>Dynamic design, numerical solution and effective verification of acceleration-level obstacle-avoidance scheme for robot manipulators</article-title>. <source>Int. J. Syst. Sci.</source> <volume>47</volume>, <fpage>932</fpage>&#x02013;<lpage>945</lpage>.<pub-id pub-id-type="doi">10.1080/00207721.2014.909971</pub-id></citation></ref>
<ref id="B54"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiao</surname> <given-names>L.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Liao</surname> <given-names>B.</given-names></name> <name><surname>Zhang</surname> <given-names>Z.</given-names></name> <name><surname>Ding</surname> <given-names>L.</given-names></name> <name><surname>Jin</surname> <given-names>L.</given-names></name></person-group> (<year>2017</year>). <article-title>A velocity-level bi-criteria optimization scheme for coordinated path tracking of dual robot manipulators using recurrent neural network</article-title>. <source>Front. Neurorobot.</source> <volume>11</volume>:<fpage>47</fpage>.<pub-id pub-id-type="doi">10.3389/fnbot.2017.00047</pub-id><pub-id pub-id-type="pmid">28928651</pub-id></citation></ref>
<ref id="B55"><citation citation-type="confproc"><person-group person-group-type="author"><name><surname>Yuan</surname> <given-names>M.</given-names></name> <name><surname>Tian</surname> <given-names>B.</given-names></name> <name><surname>Shim</surname> <given-names>V. A.</given-names></name> <name><surname>Tang</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>H.</given-names></name></person-group> (<year>2015</year>). <article-title>&#x0201C;An Entorhinal-Hippocampal Model for simultaneous cognitive map building,&#x0201D;</article-title> in <conf-name>Twenty-Ninth AAAI Conference on Artificial Intelligence</conf-name>, <conf-loc>Austin Texas</conf-loc>.</citation></ref>
<ref id="B56"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zeno</surname> <given-names>P. J.</given-names></name> <name><surname>Patel</surname> <given-names>S.</given-names></name> <name><surname>Sobh</surname> <given-names>T. M.</given-names></name></person-group> (<year>2016</year>). <article-title>Review of neurobiologically based mobile robot navigation system research performed since 2000</article-title>. <source>J. Rob.</source> <volume>2016</volume>:<fpage>8637251</fpage>.<pub-id pub-id-type="doi">10.1155/2016/8637251</pub-id></citation></ref>
<ref id="B57"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>W.-H.</given-names></name> <name><surname>Chen</surname> <given-names>A.</given-names></name> <name><surname>Rasch</surname> <given-names>M. J.</given-names></name> <name><surname>Wu</surname> <given-names>S.</given-names></name></person-group> (<year>2016</year>). <article-title>Decentralized multisensory information integration in neural systems</article-title>. <source>J. Neurosci.</source> <volume>36</volume>, <fpage>532</fpage>&#x02013;<lpage>547</lpage>.<pub-id pub-id-type="doi">10.1523/JNEUROSCI.0578-15.2016</pub-id><pub-id pub-id-type="pmid">26758843</pub-id></citation></ref>
<ref id="B58"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zilli</surname> <given-names>E. A.</given-names></name></person-group> (<year>2012</year>). <article-title>Models of grid cell spatial firing published 2005&#x02013;2011</article-title>. <source>Front. Neural Circuits</source> <volume>6</volume>:<fpage>16</fpage>.<pub-id pub-id-type="doi">10.3389/fncir.2012.00016</pub-id></citation></ref>
</ref-list>
</back>
</article>