Figures and data
Variability in grid cell properties within a module leads to enhanced encoding of local space.
When the activity of three idealized grid cells, all with the same grid spacing and orientation, are considered, the periodicity of the responses limits the amount of information conveyed about local space (Left column – “Fixed grid properties within module”). That is, there are multiple locations in physical space with identical population level activity. However, when three grid cells with variable grid spacing and orientation (in the realm of what is measured within individual grid modules – see Results), their joint activity contains considerably more information (Right column – “Variable grid properties within module”). This benefit of spatial inhomogeneity is expected to increase with larger populations of grid cells. Dashed squares in the joint activity map are enlarged below.
Grid properties are variable within a single grid module (recording ID R12).
(A) Overview of the standard procedure used to calculate the grid spacing and orientation of a given grid cell. First, spike maps are computed by identifying the location of the animal at the time of each spike. Gray line denotes the trajectory of the rat, red dots denote locations of spikes. A rate map is constructed by binning space and normalizing by the amount of time the rat spent in each spatial bin. A spatial autocorrelogram (SAC) is computed and, after the center peak is masked out (white pixels in the center of the spatial autocorrelogram – leading to change in color scale), the grid properties are fit by measuring the length and angle of the three peaks closest to 0°. (B) Example grid cells from the same module (recording ID R12), with estimated grid score, orientation (θ), and spacing (λ). (C)–(D) SAC overlaid for two pairs of grid cells [from (B)]; one pair with different θ and similar λ (C) and the other with similar θ and different λ (D). (E)–(F) Distribution of θ (E) and λ (F) across all grid cells with grid score > 0.85 (N = 74).
Variability of grid properties is a robust feature of individual grid module (recording ID R12).
Schematic overview of approach used to compute the between- and within-cell variability of grid orientation and spacing. (B)–(C) Distribution of within- and between-cell variability of θ and λ, respectively. Note that the distribution is across all 100 random shuffles of the data into two halves. (D) Average within-cell variability of grid orientation 
Variability of grid properties, restricted to the same axis, is a robust feature of individual grid module (recording ID R12).
Same analysis as in Fig. 3 (B) – (E), but for variability measured on each axis independently. For visualization, we exclude a small number of cells that were outside the axes limits, including 2, 5, and 10 cells for Axes 1–3, respectively (C); and 3, 4, and 3 cells for Axes 1–3, respectively (D); these cells were included in non-parametric statistical analyses.
Within module grid property variability is a robust feature across modules.
(A) Average within-cell variability of grid orientation 
Variability in grid properties enables improved decoding of local space from the activity of grid cells within a single module.
(A) Example noisy grid cell rate maps generated from a Poisson process. The size of the square arena is set to 1.5 m × 1.5 m to be consistent with what was used in the experimental set-up analyzed [34]. (B)–(C) Distribution of sampled grid spacing and orientation from synthetic population, when using 
Accepted and rejected cells, across all grid modules.
(Left) The percent of all cells, across all modules, that were rejected by our inclusion criteria. Cells that were rejected as not having SACs, computed with all data, that were well described by hexagonal structure (“Rejected: Poor grid fits”) are shown in red. Cells that were rejected as not reliably having SACs, computed from splits of the data, that were well described by hexagonal structure (“Rejected: Unreliable”) are shown in yellow. Cells that met these criteria (“Accepted”) are shown in blue. (Middle) Grid score of all cells, with coloring denoting whether they were accepted or rejected. Dashed gray line denotes population mean. (Right) The number of splits of the data (out of 100) that cells had SAC’s with poor grid fits, as a function of each cell’s grid score.
Variability in grid spacing, within a single module, exists when computing λ directly from the rate maps.
(A) (Left) Example grid cell rate maps from the same module (recording ID R12) with overlaid triangles, corresponding to the spacing between each of the three most central peaks. Grid spacing is computed as the average of the three lengths. (Right) To aid comparison, the triangles are enlarged (with their relative size fixed) and overlaid. Note that these cells are the same ones plotted in Fig. 2B. (B) Distribution of grid spacing computed using the SAC and the rate maps. (C) The grid spacing of all grid cells, from recording R12, computed using the SAC and the rate map. Red dashed line is the linear regression fit with R2 and p-value reported above.
Bin length does not affect the percent of cells, with good grid fits, that are accepted for analysis.
The percent of cells with good grid fits (i.e., those cells that do not get rejected for having “poor grid fits”) that are accepted by not being deemed unreliable, as a function of the size of the bins used in the shuffle analysis. All modules are plotted (each line is colored based on its recording ID).
Variability in grid properties improves decoding of local space for grid modules with different grid spacings.
(A)–(B) Same as Fig. 6F, for 
Variability in grid properties improves decoding of local space for multiple modules, when the modules have integer multiple spacing.
(A)–(B) Same as Fig. 6D, when decoding from multiple modules. (A) 
