Determinantal point process attention over grid cell code supports out of distribution generalization

eLife assessment

This important modeling work demonstrates out-of-distribution generalization using a grid cell coding scheme combined with an attentional mechanism that operates over these representations (Determinantal Point Process Attention). The simulations provide compelling evidence that the model can improve generalization performance for analogies, addition, and multiplication. The paper is significant in demonstrating how neural grid codes can support human-like generalization capabilities in analogy and arithmetic tasks, which has been a challenge for prior models.

https://doi.org/10.7554/eLife.89911.3.sa0

Significance of the findings:

Important: Findings that have theoretical or practical implications beyond a single subfield

Landmark
Fundamental
Important
Valuable
Useful

Strength of evidence:

Compelling: Evidence that features methods, data and analyses more rigorous than the current state-of-the-art

Exceptional
Compelling
Convincing
Solid
Incomplete
Inadequate

During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments

Abstract
Introduction
Methods
Results
Discussion
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

Deep neural networks have made tremendous gains in emulating human-like intelligence, and have been used increasingly as ways of understanding how the brain may solve the complex computational problems on which this relies. However, these still fall short of, and therefore fail to provide insight into how the brain supports strong forms of generalization of which humans are capable. One such case is out-of-distribution (OOD) generalization – successful performance on test examples that lie outside the distribution of the training set. Here, we identify properties of processing in the brain that may contribute to this ability. We describe a two-part algorithm that draws on speciﬁc features of neural computation to achieve OOD generalization, and provide a proof of concept by evaluating performance on two challenging cognitive tasks. First we draw on the fact that the mammalian brain represents metric spaces using grid cell code (e.g., in the entorhinal cortex): abstract representations of relational structure, organized in recurring motifs that cover the representational space. Second, we propose an attentional mechanism that operates over the grid cell code using determinantal point process (DPP), that we call DPP attention (DPP-A) – a transformation that ensures maximum sparseness in the coverage of that space. We show that a loss function that combines standard task-optimized error with DPP-A can exploit the recurring motifs in the grid cell code, and can be integrated with common architectures to achieve strong OOD generalization performance on analogy and arithmetic tasks. This provides both an interpretation of how the grid cell code in the mammalian brain may contribute to generalization performance, and at the same time a potential means for improving such capabilities in artiﬁcial neural networks.

Introduction

Deep neural networks now meet, or even exceed, human competency in many challenging task domains (He et al., 2016; Silver et al., 2017; Wu et al., 2016; He et al., 2017). Their success on these tasks, however, is generally limited to the narrow set of conditions under which they were trained, falling short of the capacity for strong forms of generalization that is central to human intelligence (Barrett et al., 2018; Lake and Baroni, 2018; Hill et al., 2019; Webb et al., 2020), and hence fail to provide insights into how our brain supports them. One such case is out-of-distribution (OOD) generalization where the test data lies outside the distribution of the training data. Here, we consider two challenging cognitive problems that often require a capacity for OOD generalization: (1) analogy and (2) arithmetic. What enables the human brain to successfully generalize on these tasks, and how might we better realize that ability in deep learning systems?

To address the problem, we focus on two properties of processing in the brain that we hypothesize are useful for OOD generalization: (1) the abstract representations of relational structure, in which relations are preserved across transformations like translation and scaling (such as observed for grid cells in mammalian medial entorhinal cortex Hafting et al., 2005); and (2) an attentional objective inspired from determinantal point processes (DPPs), which are probabilistic models of repulsion arising in quantum physics (Macchi, 1975), to attend to abstract representations that have maximum variance and minimum correlation among them, over the training data. We refer to this as DPP attention or DPP-A. The net eﬀect of these two properties is to normalize the representations of training and testing data in a way that preserves their relational structure, and allows the network to learn that structure in a form that can be applied well beyond the domain over which it was trained.

In previous work, it has been shown that such OOD generalization can be accomplished in a neural network by providing it with a mechanism for temporal context normalization (TCN) (Webb et al., 2020), a technique that allows neural networks to preserve the relational structure between the inputs in a local temporal context, while abstracting over the diﬀerences between contexts. [Temporal context normalization (Webb et al., 2020) is a normalization procedure proposed for use in training a neural network, similar to batch normalization (Ioﬀe and Szegedy, 2015), in which tensor normalization is applied over the temporal instead of the batch dimension, which is shown to help with OOD generalization. It is unrelated to the temporal context model (Howard et al., 2005), which is a computational model that proposes a role for temporal coding in the functions of the medialtemporal lobe in support of episodic recall, and spatial navigation.] Here, we test whether the same capabilities can be achieved using a well-established, biologically plausible embedding scheme – grid cell code – and an adaptive form of normalization that is based strictly on the statistics of the training data in the embedding space. We show that when deep neural networks are presented with data that exhibits such relational structure, grid cell code coupled with an error-minimizing/attentional objective promotes strong OOD generalization. We unpack each of these theoretical components in turn before describing the tasks, modeling architectures, and results.

Abstract representations of relational structure

The ﬁrst component of the proposed framework relies on the idea that a key element underlying human-like OOD generalization is the use of low-dimensional representations that emphasize the relational structure between data points. Empirical evidence suggests that, for spatial information, this is accomplished in the brain by encoding the organism’s spatial position using a periodic code consisting of diﬀerent frequencies and phases (akin to a Fourier transform of the space). Although grid cells were discovered for representations of space (Hafting et al., 2005; Sreenivasan and Fiete, 2011; Mathis et al., 2012) and used for guiding spatial behavior (Erdem and Hasselmo, 2014; Bush et al., 2015), they have since been identiﬁed in non-spatial domains, such as auditory tones (Aronov et al., 2017), odor (Bao et al., 2019), episodic memory (Chandra et al., 2023), and conceptual dimensions (Constantinescu et al., 2016). These ﬁndings suggest that the coding scheme used by grid cells may serve as a general representation of metric structure that may be exploited for reasoning about the abstract conceptual dimensions required for higher-level reasoning tasks, such as analogy and mathematics (McNamee et al., 2022). Of interest here, the periodic response function displayed by grid cells belonging to a particular frequency is invariant to translation by its period, and increasing the scale of a higher-frequency response gives a lower-frequency response and vice versa, making it invariant to scale across frequencies. This is particularly promising for prospects of OOD generalization: downstream systems that acquire parameters over a narrow training region may be able to successfully apply those parameters across transformations of translation or scale, given the shared structure (which can also be learned; Cueva and Wei, 2018; Banino et al., 2018; Whittington et al., 2020).

DPP-A

The second component of our proposed framework is a novel attentional objective that uses the statistics of the training data to sculpt the inﬂuence of grid cells on downstream computation. Despite the use of a relational encoding metric (i.e., grid cell code), generalization may also require identifying which aspects of this encoding that could potentially be shared across training and test distributions. Here, we implement this by identifying, and restricting further processing to those grid cell embeddings that exhibit the greatest variance, but are least redundant (i.e., pairwise uncorrelated) over the training data. Formally, this is captured by maximizing the determinant of the covariance matrix of the grid cell embeddings computed over the training data (Kulesza and Taskar, 2012). To avoid overﬁtting the training data, we attend to a subset of grid cell embeddings that maximize the volume in the representational space, diminishing the inﬂuence of low-variance codes (irrelevant), or codes with high similarity to other codes (redundant), which decrease the determinant of the covariance matrix.

DPP-A is inspired by mathematical work in statistical physics using DPPs that originated for modeling the distribution of fermions at thermal equilibrium (Macchi, 1975). DPPs have since been adopted in machine learning for applications in which diversity in a subset of selected items is desirable, such as recommender systems (Kulesza and Taskar, 2012). Recent work in computational cognitive science has shown DPPs naturally capture inductive biases in human inference, such as some word-learning and reasoning tasks (e.g., one noun should only refer to one object) while also serving as an eﬃcient memory code (Webb et al., 2020). In that context, the learner is biased to ﬁnd a set of possible word-meaning pairs whose representations exhibit the greatest variance and lowest covariance on a task-relevant dataset. DPPs also provide a formal objective for the type of orthogonal coding that has been proposed to be characteristic of representations in mammalian hippocampus, and integral for episodic memory (McClelland et al., 1995). Thus, using the DPP objective to govern attention over the grid cell code, known to be implemented in the entorhinal cortex (Hafting et al., 2005; Barry et al., 2007; Stensola et al., 2012; Giocomo et al., 2011; Brandon et al., 2011) (one synapse upstream of the hippocampus), aligns with the function and organization of cognitive and neural systems underlying the capability for abstraction.

Taken together, the representational and attention mechanisms outlined above deﬁne a two-component framework of neural computation for OOD generalization, by minimizing task-speciﬁc error subject to: (1) embeddings that encode relational structure among the data (grid cell code), and (2) attention to those embeddings that maximize the ‘volume’ of the representational space that is covered, while minimizing redundancy (DPP-A). Below, we demonstrate proof of concept by showing that these mechanisms allow artiﬁcial neural networks to learn representations that support OOD generalization on two challenging cognitive tasks and therefore serve as a reasonable starting point for examining the properties of interest in these networks.

Methods

Figure 1 illustrates the general framework. Task inputs, corresponding to points in a metric space, are represented as a set of grid cell embeddings $x_{t = 1.. T}$ , that are then passed to the inference module $R$ . The embedding of each input is represented by the pattern of activity of grid cells that respond selectively to diﬀerent combinations of phases and frequencies. Attention over these is a learned gating $g$ of the grid cells, the gated activations of which ( $x ⊙ g$ ) are passed to the inference module ( $R$ ). The parameterization of $g$ and $R$ are determined by backpropagation of the error signal obtained by two loss functions over the training set. Note that learning of parameter $g$ occurs only over the training space and is not further modiﬁed during testing (i.e., over the test spaces). The ﬁrst loss function, $L_{D P P}$ favors attentional gatings over the grid cells that maximize the DPP-A objective; that is, the ‘volume’ of the representational space covered by the attended grid cells. The second loss function, $L_{t a s k}$ is a standard task error term (e.g., the cross-entropy of targets $y$ and task outputs $\hat{y}$ over the training set). We describe each of these components in the following sections.

Figure 1

Download asset Open asset

Schematic of the overall framework.

Given a task (e.g., an analogy to solve), inputs (denoted as ${A, B, C, D}$ ) are represented by the grid cell code, consisting of units (grid cells) representing different combinations of frequencies and phases. Grid cell embeddings ( $x_{A}$ , $x_{B}$ , $x_{C}$ , $x_{D}$ ) are multiplied elementwise (represented as a Hadamard product $⊙$ ) by a set of learned attention gates $g$ , then passed to the inference module $R$ . The attention gates $g$ are optimized using $L_{D P P}$ , which encourages attention to grid cell embeddings that maximize the volume of the representational space. The inference module outputs a score for each candidate analogy (consisting of $A, B, C$ and a candidate answer choice $D$ ). The scores for all answer choices are passed through a softmax to generate an answer $\hat{y}$ , which is compared against the target $y$ to generate the task loss $L_{t a s k}$ .

Task setup

Analogy task

We constructed proportional analogy problems with four terms, of the form $A : B : : C : D$ , where the relation between $A$ and $B$ was the same as between $C$ and $D$ . Each of $A, B, C, D$ was a point in the integer space $ℤ^{2}$ , with each dimension sampled from the range $[0, M - 1]$ , where $M$ denotes the size of the training region. To form an analogy, two pairs of points ( $A, B$ ) and ( $C, D$ ) were chosen such that the vectors $A B$ and $C D$ were equal. Each analogy problem also contained a set of six foil items sampled in the range ${[0, M - 1]}^{2}$ excluding $D$ , such that they did not form an analogy with $A, B, C$ . The task was, given $A$ , $B$ , and $C$ , to select $D$ from a set of multiple choices consisting of $D$ and the six foil items. During training, the networks were exposed to sets of points sampled uniformly over locations in the training range, and with pairs of points forming vectors of varying length. The network was trained on 80% of all such sets of points in the training range, with 20% held out as the validation set.

To study OOD generalization, we created two cases of test data, that tested for OOD generalization in translation and scale. For the translation invariance case (Figure 2a), the constituents of the training analogies were translated along both dimensions by the same integer value (obtained by multiplying $K$ and $M$ , both of which are integer values) such that the test analogies were in the range ${[K M, (K + 1) M - 1]}^{2}$ after translation. [We transformed by the same amount along both dimensions so that the OOD generalization regimes are similar to Webb et al., 2020.] Non-overlapping test regions were generated for $K \in [1, 9]$ . Similar to the translation OOD generalization regime of Webb et al., 2020, this allowed the graded evaluation of OOD generalization to a series of increasingly remote test domains as the distance from the training region increased. For example a training analogy $A : B : : C : D$ after translation by $K M$ , would be $A + K M : B + K M : : C + K M : D + K M$ .

Figure 2

Download asset Open asset

Generation of test analogies from training analogies (region marked in blue) by: (a) translating both dimension values of $A, B, C, D$ by the same amount; and (b) scaling both dimension values of $A, B, C, D$ by the same amount.

Since both dimension values are transformed by the same amount, each input gets transformed along the diagonal.

Figure 2—source data 1 The zip file contains the data for the analogy task depicted in Figure 2.: https://cdn.elifesciences.org/articles/89911/elife-89911-fig2-data1-v1.zip
Download elife-89911-fig2-data1-v1.zip

For the scale invariance case (Figure 2b), we scaled each constituent of the training analogies by $K$ so that the test analogies after scaling were in the range ${[0, K M - 1]}^{2}$ . Thus, an analogy $A : B : : C : D$ after scaling by $K$ , would be $K A : K B : : K C : K D$ . By varying the value of $K$ from 1 to 9, we scaled the training analogies to occupy increasingly distant and larger regions of the test space. It is worth noting that while humans can exhibit complex and sophisticated forms of analogical reasoning (Holyoak, 2012; Lu et al., 2022; Webb et al., 2023), here we focused on a relatively simple form, that was inspired by Rumelhart’s parallelogram model of analogy (Mikolov et al., 2013; Rumelhart and Abrahamson, 1973) that has been used to explain traditional human verbal analogies (e.g., “king is to what as man is to woman?”). Our model, like that one, seeks to explain analogical reasoning in terms of the computation of simple Euclidean distances (i.e., $A - B = C - D$ , where $A, B, C, D$ are vectors in 2D space).

Arithmetic task

We tested two types of arithmetic operations, corresponding to the translation and scaling transformations used in the analogy tasks: elementwise addition and multiplication of two inputs $A$ and $B$ , each a point in $ℤ^{2}$ , for which $C$ was the point corresponding to the answer (i.e., $C = A + B$ or $C = A * B$ ). As with the analogy task, each arithmetic problem also contained a set of six foil items sampled in the range ${[0, M - 1]}^{2}$ , excluding $C$ . The task was to select $C$ from a set of choices consisting of $C$ and the six foil items. Similar to the analogy task, training data was constructed from a uniform distribution of points and vector lengths in the training range, with 20% held out as the validation set. To study OOD generalization, we created testing data corresponding to $K = 9$ non-overlapping regions, such that $C \in {[M, 2 M - 1]}^{2}, {[2 M, 3 M - 1]}^{2}, ... {[K M, (K + 1) M - 1]}^{2}$ .

Architecture

Grid cell code

As discussed above, the grid cell code is found in the mammalian neocortex, that support structured, low-dimensional representations of task-relevant information. For example, an organism’s location in 2D allocentric space (Hafting et al., 2005), the frequency of 1D auditory stimuli (Aronov et al., 2017), and conceptual knowledge in two continuous dimensions Doeller et al., 2010; Constantinescu et al., 2016 have all been shown to be represented as unique, similarity-preserving combinations of frequencies and phases. Here, these codes are of interest because the relational structure in the input is preserved in the code across translation and scale. This oﬀers a promising metric that can be used to learn structure relevant to the processing of analogies (Frankland et al., 2019) and arithmetic over a restricted range of stimulus values, and then used to generalize such processing to stimuli outside of the domain of task training.

To derive the grid cell code for stimuli, we follow the analytic approach described by Bicanski and Burgess, 2019 (https://github.com/bicanski/VisualGridsRecognitionMem; Bicanski, 2019). Speciﬁcally, the grid cell embedding for a particular stimulus location $A$ is given by:

x_{A} = max (0, \cos (z_{0}) + \cos (z_{1}) + \cos (z_{2}))

where,

z_{i} = b_{i} * (F A + A_{o f f s e t})

b_{0} = (\begin{matrix} \cos (0) \\ \sin (0) \end{matrix}), b_{1} = (\begin{matrix} \cos (\frac{π}{3}) \\ \sin (\frac{π}{3}) \end{matrix}), b_{2} = (\begin{matrix} \cos (\frac{2 π}{3}) \\ \sin (\frac{2 π}{3}) \end{matrix})

The spatial frequencies of grids ( $F$ ) begin at a value of $0.0028 * 2 π$ . Wei et al., 2015 have shown that, to minimize the number of variables needed to represent an integer domain of size $S$ , the ﬁring rate widths and frequencies should scale geometrically in a range $(\sqrt{2}, \sqrt{e})$ , closely matching empirically observed scaling in entorhinal cortex (Stensola et al., 2012). We choose a scaling factor of $\sqrt{2}$ to eﬃciently tile the space. One consequence of this eﬃciency is that the total number of discrete frequencies in the entorhinal cortex is expected to be small. Empirically, it has been estimated to be between 8 and 12 (Moser et al., 2015) [It seems likely that the use of grid cell code for abstraction in human cognition requires a considerably greater number of states S than that used by the rodent for sensory encoding. However, given exponential scaling, the total number of frequencies is expected to remain low, increasing as a logarithm of S.]. Here, we choose $N_{f} = 9$ (dimension of $F$ ) as the number of frequencies. $A$ refers to a particular location in a two dimensional space, and 100 oﬀsets ( $A_{o f f s e t}$ ) are used for each frequency to evenly cover a space of 1000 × 1000 locations using 900 grid cells. These oﬀsets represent diﬀerent phases within each frequency and since there are 100 of them, $N_{p} = 100$ . Hence $N_{p} \times N_{f} = 900$ , which denotes the number of grid cells. Each point from the set of 2D points for the stimuli in a task (described in Task setup), was represented using the ﬁring rate of 900 grid cells which constituted the grid cell embedding for that point to form the inputs to our model.

DPP-A

We hypothesize that the use of a relational encoding metric (i.e., grid cell code) is extremely useful, but not suﬃcient for a system to achieve strong generalization, which requires attending to particular aspects of the encoding that can capture the same relational structure across the training and test distributions. Toward this end, we propose an attentional objective that uses the statistics of the training data to attend to grid cell embeddings that can induce the inference module to achieve strong generalization. Our objective, which we describe in detail below, seeks to identify those grid cell embeddings that exhibit the greatest variance but are least redundant (pairwise uncorrelated over the training data). Formally, this is captured by maximizing the determinant of the covariance matrix of the grid cell embeddings computed over the training data (Kulesza and Taskar, 2012). Although in machine learning, DPPs have been particularly inﬂuential in work on recommender systems (Chen et al., 2018), summarization (Gong et al., 2014; Perez-Beltrachini and Lapata, 2021), and neural network pruning (Mariet and Sra, 2015), here, we propose to use maximization of the determinant of the covariance matrix as an attentional mechanism to limit the inﬂuence of grid cell embeddings with low-variance (which are less relevant) or with high similarity to other grid cell embeddings (which are redundant).

For the speciﬁc tasks that we study here, we have assumed the grid cell embeddings to be pre-learned to represent the entire space of possible test data points, and we are simply focused on the problem of how to determine which of these are most useful in enabling generalization for a task-optimized network trained on a small fraction of that space (Figure 2). That is, we look for a way to attend to a subset of grid cells frequencies whose embeddings capture recurring task-relevant relational structure. We ﬁnd that grid cell embeddings corresponding to the higher spatial frequency grid cells exhibit greater variance over the training data than the lower-frequency embeddings, while at the same time the correlations among those grid cell embeddings are lower than the correlations among the lower-frequency grid cell embeddings. The determinant of the covariance matrix of the grid cell embeddings is maximized when the variances of the grid cell embeddings are high (they are ‘expressive’) and the correlation among the grid cell embeddings is low (they ‘cover the representational space’). As a result, the higher-frequency grid cell embeddings more eﬃciently cover the representational space of the training data, allowing them to eﬃciently capture the same relational structure across training and test distributions which is required for OOD generalization.

Algorithm 1. Training with DPP-A
Parameters: inference module $R$ , attention gates $g$ Hyperparameters: number of frequencies $N_{f}$ , number of phases $N_{p}$ , number of epochs optimizing DPP attention $N_{E_{D P P}}$ , number of epochs optimizing task loss $N_{E_{t a s k}}$ , number of batches per epoch $N_{b}$ Inputs: covariance matrix $V$ , grid cell embeddings $x$ and targets $y$ for all training problems
Initialize $g$ , $R$ for $i = 1$ to $N_{E_{D P P}}$ do for $j = 1$ to $N_{b}$ do $\hat{F} (g, V, N_{f}, N_{p}) = \sum_{f = 1}^{N_{f}} \log det (d i a g (σ (g_{f})) (V_{f} - I) + I)$ $L_{D P P} = - \hat{F} (g, V, N_{f}, N_{p})$ Update $g$ end for end for ${\hat{F}}_{f \in [1, N_{f}]} = \log det (d i a g (σ (g_{f})) (V_{f} - I) + I)$ $f_{{m a x}_{D P P}} = {\arg max}_{f \in [1, N_{f}]} {\hat{F}}_{f}$ for $i = 1$ to $N_{E_{t a s k}}$ do for $j = 1$ to $N_{b}$ do $\hat{y} = R (x_{f = f_{{m a x}_{D P P}}})$ $L_{t a s k} = c r o s s - e n t r o p y (\hat{y}, y)$ Update $R$ end for end for

Algorithm 1. Training with DPP-A

Parameters: inference module

R

, attention gates

g

Hyperparameters: number of frequencies

N_{f}

, number of phases

N_{p}

, number of epochs optimizing DPP attention

N_{E_{D P P}}

, number of epochs optimizing task loss

N_{E_{t a s k}}

, number of batches per epoch

N_{b}

Inputs: covariance matrix

V

, grid cell embeddings

x

and targets

y

for all training problems

Initialize

g

R

for

i = 1

N_{E_{D P P}}

do
for

j = 1

N_{b}

\hat{F} (g, V, N_{f}, N_{p}) = \sum_{f = 1}^{N_{f}} \log det (d i a g (σ (g_{f})) (V_{f} - I) + I)

L_{D P P} = - \hat{F} (g, V, N_{f}, N_{p})

Update

g

end for
end for

{\hat{F}}_{f \in [1, N_{f}]} = \log det (d i a g (σ (g_{f})) (V_{f} - I) + I)

f_{{m a x}_{D P P}} = {\arg max}_{f \in [1, N_{f}]} {\hat{F}}_{f}

for

i = 1

N_{E_{t a s k}}

do
for

j = 1

N_{b}

\hat{y} = R (x_{f = f_{{m a x}_{D P P}}})

L_{t a s k} = c r o s s - e n t r o p y (\hat{y}, y)

Update

R

end for
end for

Formally, we treat obtaining $L_{D P P}$ as an approximation of a DPP. A DPP $P$ deﬁnes a probability measure on all subsets of a set of items $U = {1, 2, . . . N}$ . For every $u \subseteq U$ , $P (u) \propto det (V_{u})$ . Here, $V$ is a positive semideﬁnite covariance matrix and $V_{u} = {[V_{i j}]}_{i, j \in u}$ denotes the matrix $V$ restricted to the entries indexed by elements of $u$ . The maximum a posteriori (MAP) problem ${a r g m a x}_{u} det (V_{u})$ is NP-hard (Ko et al., 1995). However, $f (u) = \log (det (V_{u}))$ satisﬁes the property of a submodular function, and various algorithms exist for approximately maximizing them. One common way is to approximate this discrete optimization problem by replacing the discrete variables with continuous variables and extending the objective function to the continuous domain. Gillenwater et al., 2012 proposed a continuous extension that is eﬃciently computable and diﬀerentiable:

\hat{F} (w) = \log \sum_{u} \prod_{i \in u} w_{i} \prod_{i \in u} (1 - w_{i}) \exp (f (u)), w \in {[0, 1]}^{N} .

We use the following theorem from Gillenwater et al., 2012 to construct $L_{D P P}$ :

Theorem 2.1

For a positive semideﬁnite matrix $V$ and $w \in {[0, 1]}^{N}$ :

\sum_{u} \prod_{i \in u} w_{i} \prod_{i \in u} (1 - w_{i}) det (V_{u}) = det (d i a g (w) (V - I) + I)

We propose an attention mechanism that uses $L_{D P P}$ to attend to subsets of grid cell embeddings for further processing. Algorithm 1 describes the training procedure with DPP-A which consists of two steps, using $L_{D P P}$ as the ﬁrst step. This step maximizes the objective function:

\hat{F} (g, V, N_{f}, N_{p}) = \sum_{f = 1}^{N_{f}} \log det (d i a g (σ (g_{f})) (V_{f} - I) + I)

using stochastic gradient ascent for $N_{E_{D P P}}$ epochs, which is equivalent to minimizing $L_{D P P}$ , as $L_{D P P} = - \hat{F} (g, V, N_{f}, N_{p})$ . It involves attending to grid cell embeddings that exhibit the greatest within-frequency variance but are least redundant (i.e., that are least also pairwise uncorrelated) over the training data. This is achieved by maximizing the determinant of the covariance matrix over the within-frequency grid cell embeddings of the training data, and Equation 6 is obtained by applying log on both sides of the Theorem 2.1, and in our case $U$ refers to grid cells of a particular frequency. Here, $g$ are the attention gates corresponding to each grid cell, and $N_{f}$ is the number of distinct frequencies. The matrix $V$ captures a measure of the covariance of the grid cell embeddings over the training region. We used the $s y n t h_k e r n e l$ function (https://github.com/insuhan/fastdppmap/blob/db7a28c38ce654bdbfd5ab1128d3d5910b68df6b/test_greedy.m#L123; Han, 2017) to construct $V$ , where in our case $m$ are the variances of the grid cell embeddings $S$ computed over the training region $M$ , $N$ is the number of grid cells and $w_{m}, b$ are hyperparameters with values of 1 and 0.1, respectively. [ $S$ need not be a square matrix in our case, whose second dimension M is the size of the training region. L_kernel is same as $V$ .] The dimensionality of $V$ is $N_{f} N_{p} \times N_{f} N_{p} (900 \times 900)$ . $g_{f}$ are the gates of the grid cells belonging to the $f$ th frequency, so $g_{f} = g [f N_{p} : (f + 1) N_{p}]$ , where $N_{p}$ is the number of phases for each frequency. $V_{f} = V [f N_{p} : (f + 1) N_{p}, f N_{p} : (f + 1) N_{p}]$ is the restriction of $V$ to the grid cell embeddings for $f$ th frequency, so it captured the covariance of the grid cell embeddings belonging to the $f$ th frequency. $σ$ is sigmoid nonlinearity applied to $g_{f}$ , deﬁned as $σ (g_{f}) = 1 / (1 + e^{- g_{f}})$ , so that their values are between 0 and 1. $d i a g (σ (g_{f}))$ converts vector $σ (g_{f})$ into a matrix with $σ (g_{f})$ as the diagonal of the matrix and the rest entries are zero, which is multiplied with $V - I$ , which results in elementwise multiplication of $σ (g_{f})$ with the column vectors of $V - I$ . Equation 6 which involved summation of the logarithm of the determinant of the gated covariance matrix of grid cell embeddings within each frequency, over $N_{f}$ frequencies was used to compute the negative of $L_{D P P}$ . Maximizing $\hat{F}$ gave the approximate maximum within-frequency log determinant for each frequency $f \in [1, N_{f}]$ , which we denote for the $f$ th frequency as ${\hat{F}}_{f}$ . In the second step of the training procedure, we used the $f_{m a x_{D P P}}$ grid cell frequency, where $f_{m a x_{D P P}} = {a r g m a x}_{f \in [1, N_{f}]} {\hat{F}}_{f}$ . In other words, we used the grid cell embeddings for the frequency which had the maximum within-frequency log determinant at the end of the ﬁrst step, which we ﬁnd are best at capturing the relational structure across the training and testing data, thereby promoting OOD generalization. In this step, we trained the inference module minimizing $L_{t a s k}$ over $N_{E_{t a s k}}$ epochs. More details can be found in Appendix 1: ‘DPP-A attentional modulation’.

Inference module

We implemented the inference module $R$ in two forms, one using Long Short Term Memory (LSTM) (Hochreiter and Schmidhuber, 1997) and the other using a transformer (Vaswani et al., 2017) architecture. Separate networks were trained for the analogy and arithmetic tasks using each form of inference module. For each task, the attended grid cell embeddings of each stimulus obtained from the DPP-A process ( $x_{f = f_{m a x_{D P P}}}$ ), were provided to $R$ as its inputs, which we denote as $x_{R}$ for brevity. For the arithmetic task, we also concatenated to $x_{R}$ a one-hot tensor of dimension 2, before feeding to $R$ that speciﬁed which computation to perform (addition or multiplication). As proposed by Hill et al., 2019, we treated both the analogy and arithmetic tasks as scoring (i.e., multiple choice) problems. For each analogy, the inference module was presented with multiple problems, each consisting of three stimuli, $A, B, C$ , and a candidate completion from the set containing $D$ (the correct completion) and six foil completions. For each instance of the arithmetic task, it was presented with two stimuli, $A, B$ , and a candidate completion from the set containing $C$ (the correct completion) and six foil completions. Stimuli were presented sequentially for $R$ constructed using an LSTM, which consists of three gates, and computations using them deﬁned as below:

i n p u t_g a t e [t] = σ (W_{i i} x_{R} [t] + b_{i i} + W_{h i} h [t - 1] + b_{h i})

f o r g e t_g a t e [t] = σ (W_{i f} x_{R} [t] + b_{i f} + W_{h f} h [t - 1] + b_{h f})

o u t p u t_g a t e [t] = σ (W_{i o} x_{R} [t] + b_{i o} + W_{h o} h [t - 1] + b_{h o})

c [t] = f o r g e t_g a t e [t] ⊙ c [t - 1] + i n p u t_g a t e [t] ⊙ \tanh (W_{i c} x_{R} [t] + b_{i c} + W_{h c} h [t - 1] + b_{h c})

h [t] = o u t p u t_g a t e [t] ⊙ \tanh (c [t])

where $W_{i i}$ , $W_{h i}$ , $W_{i f}$ , $W_{h f}$ , $W_{i o}$ , $W_{h o}$ , $W_{i c}$ , $W_{h c}$ are weight matrices and $b_{i i}$ , $b_{h i}$ , $b_{i f}$ , $b_{h f}$ , $b_{i o}$ , $b_{h o}$ , $b_{i c}$ , $b_{h c}$ are bias vectors. $h [t]$ and $c [t]$ are the hidden state and the cell state at time $t$ , respectively. The hidden state for the last timestep was passed through a linear layer with a single output unit to generate a score for the candidate completions for each problem. We used a single layered LSTM of 512 hidden units, which corresponds to the size of the hidden state, cell state, bias vectors, and weight matrices. The hidden and cell state of the LSTM was reinitialized at the start of each sequence for each candidate completion.

For $R$ constructed using a transformer, we used the standard multi-head self attention ( $M H S A$ ) mechanism followed by a multi-layered perceptron ( $M L P$ ) – a feedforward neural network with one hidden layer, with layer normalization (Ba et al., 2016) ( $L N$ ) to transform $x_{R}$ at each layer of the transformer, deﬁned as below:

S A (Q, K, V) = s o f t m a x (\frac{{Q K}^{T}}{\sqrt{d_{k}}}) V

M H S A (x_{R}) = C o n c a t (S A (W_{1}^{Q} x_{R}, W_{1}^{K} x_{R}, W_{1}^{V} x_{R}), . . . ., S A (W_{H}^{Q} x_{R}, W_{H}^{K} x_{R}, W_{H}^{V} x_{R})) W^{O}

{\tilde{x}}_{R} = M L P (L N (M H S A (L N (x_{R})) + x_{R})) + M H S A (L N (x_{R})) + x_{R}

where $Q$ , $K$ , and $V$ are called the query, key, and value matrices, respectively, $d_{k}$ is the dimension of the matrices, $W^{Q}$ , $W^{K}$ , and $W^{V}$ are the corresponding projection matrices, $W^{O}$ is the output projection matrix which is applied to the concatenation of self attention ( $S A$ ) output for each head, and $H$ is the number of heads. The $s o f t m a x$ function is used to convert real valued vector inputs into a probability distribution, deﬁned as $s o f t m a x (z_{i}) = e^{z_{i}} / \sum_{i} e^{z_{i}}$ . We used a transformer with 6 layers, each of which had 8 heads, $d_{k} = 32$ , and $M L P$ hidden layer dimension of 512. The stimuli were presented together and projected into 128 dimensions to be more easily divisible by the number of heads, followed by layer normalization. Since a transformer is naturally invariant to the order of the stimuli, to make use of the order we added to $x_{R}$ a learnable positional encoding (Kazemnejad, 2019), which is the linear projection of one-hot tensors denoting the position of stimuli in the sequence. We then concatenated a learned CLS (short for ‘classiﬁcation’) token (analogous to the CLS token in Devlin et al., 2018) to $x_{R}$ , before transforming with Equation 14. We took the transformed value ( ${\tilde{x}}_{R}$ ) corresponding to the CLS token, and passed it to a linear layer with 1 output unit to generate a score for each candidate completion. This procedure was repeated for each candidate completion.

The seven scores (one for the correct completion and for six foil completions) were normalized using a softmax function, such that a higher score would correspond to a higher probability and vice versa, and the probabilities sum to 1. The inference module was trained using the task-speciﬁc cross-entropy loss ( $L_{t a s k} = c r o s s - e n t r o p y$ ) between the softmax-normalized scores and the index for the correct completion (target), deﬁned as $- \log (s o f t m a x (s c o r e s) [t a r g e t])$ . It is worth noting that the properties of equivariance hold, since the probability distribution after applying softmax remains the same when the transformation (translation or scaling) is applied to the scores for each of the answer choices obtained from the output of the inference module, and when the same transformation is applied to the stimuli for the task and all the answer choices before presenting as input to the inference module to obtain the scores. We also tried formulating the tasks as regression problems, the details of which can be found in Appendix 1: ‘Regression formulation’.

While our model is not construed to be as speciﬁc about the implementation of the $R$ module, we assume that – as a standard deep learning component – it is likely to map onto neocortical structures that interact with the entorhinal cortex and, in particular, regions of the prefrontal–posterior parietal network widely believed to be involved in abstract relational processes (Waltz et al., 1999; Christoff et al., 2001; Knowlton et al., 2012; Summerfield et al., 2020). In particular, the role of the prefrontal cortex in the encoding and active maintenance of abstract information needed for task performance (such as rules and relations) has often been modeled using gated recurrent networks, such as LSTMs (Frank et al., 2001; Braver and Cohen, 2000), and the posterior parietal cortex has long been known to support ‘maps’ that may provide an important substrate for computing complex relations (Summerfield et al., 2020).

Experiments

Experimental details

For each task, the sequence of stimuli for a given problem was encoded using grid cell code (see Grid cell code), that were then modulated by DPP-A (see DPP-A), and passed to the inference module $R$ (see Inference module). We trained three networks using each type of inference module. For networks using an LSTM in the inference module, we trained each network for number of epochs for optimizing DPP attention $N_{E_{D P P}} = 50$ , number of epochs for optimizing task loss $N_{E_{t a s k}} = 50$ , on analogy problems, and for $N_{E_{D P P}} = 500$ , $N_{E_{t a s k}} = 500,$ on arithmetic problems with a batch size of 256, using the ADAM optimizer (Kingma and Ba, 2014), and a learning rate of $1 e^{- 3}$ . For networks using a transformer in the inference module, we trained with a batch size of 128 on analogy with a learning rate of $5 e^{- 4}$ , and on arithmetic problems with a learning rate of $5 e^{- 5}$ . More details can be found in Appendix 1: ‘More experimental details’.

Comparison models

To evaluate how the grid cell code coupled with DPP-A compares with other architectures and approaches to generalization, and the extent to which each of these components contributed to the performance of the model, we compared it with several alternative models for performing the analogy and arithmetic tasks. First, we compared it with the TCN model (Webb et al., 2020) (see Related work), but modiﬁed so as to use the grid cell code as input. We passed the grid cell embeddings for each input through a shared feedforward encoder which consisted of two fully connected layers with 256 units per layer. ReLU nonlinearities were used in both the layers. The ﬁnal embedding was generated with a linear layer of 256 units. TCN was applied to these embeddings and then passed as a sequence for each candidate completion to the inference module. This implementation of TCN involved a learned encoder on top of the grid cell embeddings, so it is closely analogous to the original TCN.

Next, we compared our model to one that used variational dropout (Gal and Ghahramani, 2016), which is shown to be more eﬀective in generalization compared to naive dropout (Srivastava et al., 2014). We randomly sampled a dropout mask (50% dropout), zeroing out the contribution of some of the grid cell code in the input to the inference module. We then use that locked dropout mask for every time step in the sequence. We also compared DPP-A to a model that had an additional penalty (L1 regularization) proportional to the absolute sum of the attention gates $g$ along with the task-speciﬁc loss. We experimented with diﬀerent values of λ, which controlled the strength of the penalty relative to the cross-entropy loss. We report accuracy values for λ that achieved the best performance on the validation set. Accuracy values for various λs are provided in Appendix 1: ‘Eﬀect of L1 Regularization strength (λ)’. Dropout and L1 regularization were chosen as a proxy for DPP-A and hence we used the same input data for fair comparison. Finally, we compared to a model that used the complete grid cell code, that is no DPP-A.

Related work

A body of recent computational work has shown that representations similar to grid cells can be derived using standard analytical techniques for dimensionality reduction (Dordek et al., 2016; Stachenfeld et al., 2017), as well as from error-driven learning paradigms (Cueva and Wei, 2018; Banino et al., 2018; Whittington et al., 2020; Sorscher et al., 2023). Previous work has also shown that grid cells emerge in networks trained to generalize to novel location/object combinations, and support transitive inference (Whittington et al., 2020). Here, we show that grid cells enable strong OOD generalization when coupled with the appropriate attentional mechanism. Our proposed method is thus complementary to these previous approaches for obtaining grid cell representations from raw data.

In the ﬁeld of machine learning, DPPs have been used for supervised video summarization (Gong et al., 2014), diverse recommendations (Chen et al., 2018), selecting a subset of diverse neurons to prune a neural network without hurting performance (Mariet and Sra, 2015), and diverse minibatch attention for stochastic gradient descent (Zhang et al., 2017). Recently, Mariet et al., 2019 generated approximate DPP samples by proposing an inhibitive attention mechanism based on transformer networks as a proxy for capturing the dissimilarity between feature vectors, and Perez-Beltrachini and Lapata, 2021 used DPP-based attention with seq-to-seq architectures for diverse and relevant multi-document summarization. To our knowledge, however, DPPs have not previously been combined with the grid cell code that we employ here, and have not been used to enable OOD generalization.

Various approaches have been proposed to prevent deep learning systems from overﬁtting, and enable them to generalize. A commonly employed technique is weight decay (Krogh and Hertz, 1992). Srivastava et al., 2014 proposed dropout, a regularization technique which reduces overﬁtting by randomly zeroing units from the neural network during training. Recently, Webb et al., 2020 proposed TCN in which a normalization similar to batch normalization (Ioﬀe and Szegedy, 2015) was applied over the temporal dimension instead of the batch dimension. However, unlike these previous approaches, the method reported here achieves nearly perfect OOD generalization when operating over the appropriate representation, as we show in the results. Our proposed method also has the virtue of being based on a well understood, and biologically plausible, encoding scheme (grid cell code).

Results

Analogy

We ﬁrst present results on the analogy task for two types of testing data, translation and scaling using two types of inference module, LSTM and transformer. We trained three networks for each method and report mean accuracy along with standard error of the mean. Figure 3 shows the results for the analogy task using an LSTM in the inference module. The left panel shows results for the translation regime, and the right panel shows results for the scaling regime. Both plots show accuracy on the training and validation sets, and on a series of nine (increasingly distant) OOD generalization test regions. DPP-A (shown in blue) achieves nearly perfect accuracy on all of the test regions, considerably outperforming the other models.

Figure 3

Download asset Open asset

Results on analogy on each region for translation and scaling using LSTM in the inference module.

For the case of translation, using all the grid cell code with no DPP-A (shown in purple) led to the worst OOD generalization performance, overﬁtting on the training set. Locked dropout (denoted by green) and L1 regularization (denoted by red) reduced overﬁtting and demonstrated better OOD generalization performance than no DPP-A but still performed considerably worse than DPP-A. For the case of scaling, locked dropout and L1 regularization performed slightly better than TCN, achieving marginally higher test accuracy, but DPP-A still substantially outperformed all other models, with a nearly 70% improvement in average test accuracy.

To test the generality of the grid cell code and DPP-A across network architectures, we also tested a transformer (Vaswani et al., 2017) in place of the LSTM in the inference module. Previous work has suggested that transformers are particularly useful for extracting structure in sequential data and have been used for OOD generalization (Saxton et al., 2019). Figure 4 shows the results for the analogy task using a transformer in the inference module. With no explicit attention (no DPP-A) over the grid cell code (shown in orange), the transformer did poorly on the analogies on the test regions. This suggests that simply using a more sophisticated architecture with standard forms of attention is not suﬃcient to enable OOD generalization based on the grid cell code. With DPP-A (shown in blue), the OOD generalization performance of the transformer is nearly perfect for both translation and scaling. These results also demonstrate that grid cell code coupled with DPP-A can be exploited for OOD generalization eﬀectively by diﬀerent architectures.

Figure 4

Download asset Open asset

Results on analogy on each region for translation and scaling using the transformer in the inference module.

Arithmetic

We next present results on the arithmetic task for two types of operations, addition and multiplication using two types of inference module, LSTM and transformer. We trained three networks for each method and report mean accuracy along with standard error of the mean.

Figure 5 shows the results for arithmetic problems using an LSTM in the inference module. The left panel shows results for addition problems, and the right panel shows results for multiplication problems. DPP-A achieves higher accuracy for addition than multiplication on the test regions. However, in both cases DPP-A signiﬁcantly outperforms the other models, achieving nearly perfect OOD generalization for addition, and 65% accuracy for multiplication, compared with under 20% accuracy for all the other models. We found that grid cell embeddings obtained after the ﬁrst step in Algorithm 1 are not able to fully preserve the relational structure for multiplication problems on the test regions (more details in Appendix 1: ‘Why is OOD generalization performance worse for the multiplication task?’), but still it aﬀords superior capacity for OOD generalization than any of the other models. Thus, while it does not match the generalizability of a genuine algorithmic (i.e., symbolic) arithmetic function, it may be suﬃcient for some tasks (e.g., approximate multiplication ability in young children; Qu et al., 2021).

Figure 5

Download asset Open asset

Results on arithmetic on each region using LSTM in the inference module.

Figure 6 shows the results for arithmetic problems using a transformer in the inference module. With no DPP-A over the grid cell code the transformer did poorly on addition and multiplication on the test regions, achieving around 20–30% accuracy. With DPP-A, the OOD generalization performance of the transformer shows a pattern similar to that for LSTM: it is nearly perfect for addition and, though not as good on multiplication, nevertheless shows approximately 40% better performance than the transformer multiplication.

Figure 6

Download asset Open asset

Results on arithmetic on each region using the transformer in the inference module.

Ablation study

To determine the individual importance of the grid cell code and the DPP-A attention objective, we carried out several ablation studies. First, to evaluate the importance of grid cell code, we analyzed the eﬀect of DPP-A with other embeddings, using either one-hot or smoothed one-hot embeddings (similar to place cell coding) with standard deviations of 1, 10, and 100, each passed through a learned feedforward encoder, which consisted of two fully connected layers with 1024 units per layer, and ReLU nonlinearities. The ﬁnal embedding was generated with a fully connected layer with 1024 units and sigmoid nonlinearity. Since these embeddings do not have a frequency component, the training procedure with DPP-A consisted of only one step: minimizing the loss function $L = L_{t a s k} - λ * \hat{F} (g, V)$ . We tried diﬀerent values of λ (0.001, 0.01, 0.1, 1, 10, 100, 1000, 10,000). For each type of embedding (one-hots and smoothed one-hots with each value of standard deviation), we trained three networks and report for the model that achieved best performance on the validation set. Note that, given the much higher dimensionality and therefore memory demands of embeddings based on one-hots and smoothed one-hots, we had to limit the evaluation to a subset of the total space, resulting in evaluation on only two test regions (i.e., $K \in [1, 3]$ ).

Figure 7 shows the results for the analogy task (results for the arithmetic task are in Appendix 1: ‘Ablation study on arithmetic task’, Appendix 1—figure 3) using an LSTM in the inference module. The average accuracy on the test regions for translation and scaling using smoothed one-hots passed through an encoder (shown in green) is nearly 30% better than simple one-hot embeddings passed through an encoder (shown in orange), but both still achieve signiﬁcantly lower test accuracy than grid cell code which support perfect OOD generalization.

Figure 7

Download asset Open asset

Results on analogy on each region using determinantal point process attention (DPP-A), an LSTM in the inference module, and different embeddings (grid cell code, one-hots, and smoothed one-hots passed through a learned encoder) for translation (left) and scaling (right).

Each point is mean accuracy over three networks, and bars show standard error of the mean.

With respect to the importance of the DPP-A, we note that the simulations reported earlier show that replacing the DPP-A mechanism with either other forms of regularization (dropout and L1 regularization; see Comparison models) or a transformer (Figure 4 in Analogy for analogy and Figure 6 in Arithmetic for arithmetic tasks) failed to achieve the same level of OOD generalization as the network that used DPP-A. The results using a transformer are particularly instructive, as that incorporates a powerful mechanism for learned attention, but, even when provided with grid cell embeddings, failed to produce results comparable to DPP-A, suggesting that the latter provides a simple but powerful form of attention objective, at least when used in conjunction with grid cell embeddings.

Finally, for completeness, we also carried out a set of simulations that examined the performance of networks with various embeddings (grid cell code, and one-hots or smoothed one-hots with or without a learned feedforward encoder), but no attention or regularization (i.e., neither DPP-A, transformer, nor TCN, L1 Regularization, or Dropout). Figure 8 shows the results for the diﬀerent embeddings on the analogy task (results for the arithmetic task are in Appendix 1: ‘Ablation study on arithmetic task’, Appendix 1—figure 4). For translation (left), the average accuracy over the test regions using grid cell code (shown in blue) is nearly 25% more compared to other embeddings, which all yield performance near chance (∼15%). For scaling (right), although other embeddings achieve higher performance than chance (except smoothed one-hots), they still achieve lower test accuracy than grid cell code. More ablation studies can be found in Appendix 1: ‘Ablation study on choice of frequency’, ‘Ablation on DPP-A’, and Figure A5.

Figure 8

Download asset Open asset

Results on analogy on each region using different embeddings (grid cell code, and one-hots or smoothed one-hots with and without an encoder) and an LSTM in the inference module, but without determinantal point process attention (DPP-A), temporal context normalization (TCN), L1 Regularization, or Dropout for translation (left) and scaling (right).

Discussion

We have identiﬁed how particular properties of processing observed in the brain can be used to achieve strong OOD generalization, and introduced a two-component algorithm to promote OOD generalization in deep neural networks. The ﬁrst component is a structured representation of the training data, modeled closely on known properties of grid cells – a population of cells that collectively represent abstract position using a periodic code. However, despite their intrinsic structure, we ﬁnd that grid cell code and standard error-driven learning alone are insuﬃcient to promote OOD generalization, and standard approaches for preventing overﬁtting oﬀer only modest gains. This is addressed by the second component, using DPP-A to implement attentional regularization over the grid cell code. DPP-A allows only a relevant and diverse subset of cells to inﬂuence downstream computation in the inference module using the statistics of the training data. For proof of concept, we started with two challenging cognitive tasks (analogy and arithmetic), and showed that the combination of grid cell code and DPP-A promotes OOD generalization across both translation and scale when incorporated into common architectures (LSTM and transformer). It is worth noting that the DPP-A attentional mechanism is diﬀerent from the attentional mechanism in the transformer module, and both are needed for strong OOD generalization in the case of transformers. Use of the standard self-attention mechanism in transformers over the inputs (i.e., A, B, C, and D for the analogy task) in place of DPP-A would lead to weightings of grid cell embeddings over all frequencies and phases. The objective function for the DPP-A represents an inductive bias, that selectively assigns the greatest weight to all grid cell embeddings (i.e., for all phases) of the frequency for which the determinant of the covariance matrix is greatest computed over the training space. The transformer inference module then attends over the inputs with the selected grid cell embeddings based on the DPP-A objective.

The current approach has some limitations and presents interesting directions for future work. First, we derive the grid cell code from known properties of neural systems, rather than obtaining the code directly from real-world data. Here, we are encouraged by the body of work providing evidence for grid cell code in the hidden layers of neural networks in a variety of task contexts and architectures (Wei et al., 2015; Cueva and Wei, 2018; Banino et al., 2018; Whittington et al., 2020). This suggests reason for optimism that DPP-A may promote strong generalization in cases where grid cell code naturally emerge: for example, navigation tasks (Banino et al., 2018) and reasoning by transitive inference (Whittington et al., 2020). Integrating our approach with structured representations acquired from high-dimensional, naturalistic datasets remains a critical next step which would have signiﬁcant potential for broader future practical applications. So too does application to more complex transformations beyond translation and scale, such as rotation, and complex forms of representations, and analogical reasoning tasks (Holyoak, 2012; Webb et al., 2023; Lu et al., 2022). Second, it is not clear how DPP-A might be implemented in a neural network. In that regard, Bozkurt et al., 2022 recently proposed a biologically plausible neural network algorithm using a weighted similarity matrix approach to implement a determinant maximization criterion, which is the core idea underlying the objective function we use for DPP-A (Equation 6), suggesting that the DPP-A mechanism we describe may also be biologically plausible. This could be tested experimentally by exposing individuals (e.g., rodents or humans) to a task that requires consistent exposure to a subregion, and evaluating the distribution of activity over the grid cells. Our model predicts that high-frequency grid cells should increase their ﬁring rate more than low-frequency cells, since the high-frequency grid cells maximize the determinant of the covariance matrix of the grid cell embeddings. It is also worth noting that Frankland and Cohen, 2020 have suggested that the use of DPPs may also help explain a mutual exclusivity bias observed in human word learning and reasoning. While this is not direct evidence of biological plausibility, it is consistent with the idea that the human brain selects representations for processing that maximize the volume of the representational space, which can be achieved by maximizing the DPP-A objective function deﬁned in Equation 6. Third, we compared grid cell code to only one-hots and place cell code. Future work could compare to a broader range of potential biological coding schemes for the overall space, for example boundary vector cell coding (Barry et al., 2006), band cell coding, or egocentric boundary cell coding (Hinman et al., 2019).

Finally, we focus on analogies in linear spaces which limits the generality of our approach in nonlinear spaces. In that regard, there are at least two potential directions that could be pursued. One is to directly encode nonlinear structures (such as trees and rings) with grid cells, to which DPP-A could be applied as described in our model. The Tolman-Eichenbaum Machine (TEM) model (Whittington et al., 2020) suggests that grid cells in the medial entorhinal may form a basis set that captures structural knowledge for such nonlinear spaces, such as social hierarchies and transitive inference when formalized as a connected graph. Another would be to use eigen decomposition of the successor representation (Dayan, 1993), a learnable predictive representation of possible future states that has been shown by Stachenfeld et al., 2017 to provide an abstract structured representation of a space that is analogous to the grid cell code. This general-purpose mechanism could be applied to represent analogies in nonlinear spaces (Frankland et al., 2019), for which there may not be a clear factorization in terms of grid cells (i.e., distinct frequencies and multiple phases within each frequency). Since the DPP-A mechanism, as we have described it, requires representations to be factored in this way it would need to be modiﬁed for such purpose. Either of these approaches, if successful, would allow our model to be extended to domains containing nonlinear forms of structure. To the extent that diﬀerent coding schemes (i.e., basis sets) are needed for diﬀerent forms of structure, the question of how these are identiﬁed and engaged for use in a given setting is clearly an important one, that is not addressed by the current work. We imagine that this is likely subserved by monitoring and selection mechanisms proposed to underlie the capacity for selective attention and cognitive control (Shenhav et al., 2013), though the speciﬁc computational mechanisms that underlie this function remain an important direction for future research.

Appendix 1

Code and data availability

The code and the data can be downloaded from https://github.com/Shanka123/DPP-Attention_Grid-Cell-Code (copy archived at Mondal, 2024).

More experimental details

The size of the training region was 100. For the analogy task, we used 653,216 training samples, 163,304 validation samples, and 20,000 testing samples for each of the nine regions. For the arithmetic task, we used 80,000 training samples, 20,000 validation samples, and 20,000 testing samples for each of the nine regions with equal number of addition and multiplication problems. We used the PyTorch library (Paszke et al., 2017) for all experiments. For each network, the training epoch that achieved the best validation accuracy was used to report performance accuracy for the training stimulus sets, validation sets (held out stimuli from the training range), and OOD generalization test sets (from regions beyond the range of the training data).

Why is OOD generalization performance worse for the multiplication task?

In an eﬀort to understand why DPP-A achieved around 65% average test accuracy on multiplication compared to nearly perfect accuracy for addition and analogy tasks, we analyzed the distribution of the grid cell embeddings for the frequency which had the maximum within-frequency determinant at the end of the ﬁrst step in Algorithm 1. More speciﬁcally for $A$ , $B$ , and the correct answer $C$ , we analyzed the distribution of each grid cell for the training and the nine test regions. Note that since the total number of grid cells was 900 and there were nine frequencies, the dimension of the grid cell embeddings corresponding to $f_{m a x_{D P P}}$ grid cell frequency was 100. To quantify the similarity between training and the test distributions, we computed cosine distance (1 − cosine similarity), and averaged it over the 100 dimensions and nine test regions. We found that the average cosine distance is 5× greater for multiplication than addition problem (0.0002 for addition: 0.001 for multiplication). In this respect, grid cell code does not perfectly preserve the relational structure of the multiplication problem, which we would expect to limit DPP-A’s OOD generalization ability in that task domain.

Ablation study on choice of frequency

Appendix 1—figure 1

Download asset Open asset

Results on analogy on each region using LSTM in the inference module for choosing top $K$ frequencies with ${\hat{F}}_{f}$ in Algorithm 1.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

Baseline using dynamic attention across frequencies

Appendix 1—figure 2

Download asset Open asset

Ablation study on arithmetic task

Appendix 1—figure 3

Download asset Open asset

Results on arithmetic with different embeddings (with determinantal point process attention [DPP-A]) using LSTM in the inference module.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

Appendix 1—figure 4

Download asset Open asset

Results on arithmetic with different embeddings (without determinantal point process attention [DPP-A], temporal context normalization [TCN], L1 Regularization, or Dropout) using LSTM in the inference module.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

Appendix 1—figure 5

Download asset Open asset

Results on arithmetic for increasing number of grid cell frequencies $N_{f}$ on each region using LSTM in the inference module.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

Regression formulation

Appendix 1—figure 6

Download asset Open asset

Results for regression on analogy using LSTM in the inference module.

Results show mean squared error on each region averaged over three trained networks along with errorbar (standard error of the mean).

Appendix 1—figure 7

Download asset Open asset

Results for regression on arithmetic on each region using LSTM in the inference module.

Results show mean squared error on each region averaged over three trained networks along with errorbar (standard error of the mean).

We also tried formulating the analogy and arithmetic tasks as regression instead of classiﬁcation via a scoring mechanism. For DPP-A, the inference module was trained to generate the grid cell embeddings belonging to the $f_{m a x_{D P P}}$ frequency, which had the maximum within-frequency determinant at the end of the ﬁrst step in Algorithm 1 for the correct completion, given as input the $f_{m a x_{D P P}}$ frequency grid cell embeddings for $A$ , $B$ , $C$ for the analogy task and $A$ , $B$ for the arithmetic task. A linear layer with 100 units and sigmoid activation was used to generate the output of the inference module and was trained to minimize the mean squared error with the $f_{m a x_{D P P}}$ frequency grid cell embeddings of the correct completion. We compared DPP-A with a version that didn’t use the attentional objective (no DPP-A), where the inference module was trained to generate the grid cell embeddings for all the frequencies, but was evaluated on only the $f_{m a x_{D P P}}$ frequency grid cell embeddings for fair comparison with the DPP-A version. Appendix 1—figure 6 shows the results for the analogy task using an LSTM in the inference module. For both the translation (left) and scaling (right) regimes, DPP-A achieves nearly zero mean squared error on all the test regions, considerably outperforming the no DPP-A which achieves a much higher error. Appendix 1—figure 7 shows the results for arithmetic problems using an LSTM in the inference module. For addition problems, shown on the left, DPP-A achieves nearly zero mean squared error on the test regions. For multiplication problems, shown on the right, DPP-A achieves a lower mean squared error on the test regions, 0.11, compared to no DPP-A which achieves around 0.17.

Eﬀect of L1 regularization strength (λ)

Appendix 1—figure 8

Download asset Open asset

Results on analogy for L1 regularization for various λs for translation and scaling using LSTM in the inference module.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

Appendix 1—figure 9

Download asset Open asset

Results on arithmetic for L1 regularization for various λs using LSTM in the inference module.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

Ablation on DPP-A

Appendix 1—figure 10

Download asset Open asset

Results on analogy for one step determinantal point process attention (DPP-A) over the complete grid cell code for various λs for translation and scaling using LSTM in the inference module.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

Appendix 1—figure 11

Download asset Open asset

Results on analogy for one step determinantal point process attention (DPP-A) within frequencies for various λs for translation and scaling using LSTM in the inference module.

Results show mean accuracy on each region averaged over three trained networks along with errorbar (standard error of the mean).

The proposed DPP-A method (Algorithm 1) consists of two steps with $L_{D P P}$ in the ﬁrst step and $L_{t a s k}$ in the second step. We considered two ablation experiments which consists of a single step. In one case, we maximized the objective function, $\hat{F} (g, V) = \log det (d i a g (σ (g)) (V - I) + I)$ , over the grid cell embeddings of all frequencies and phases (instead of summing $\hat{F}$ corresponding to the grid cell embeddings from each frequency independently as done in the ﬁrst step of Algorithm 1), using stochastic gradient ascent, along with minimizing $L_{t a s k}$ , which would use all the attended grid cell embeddings (instead of using $f_{m a x_{D P P}}$ frequency grid cell embeddings as done in the second step of Algorithm 1). So total loss, $L = L_{t a s k} - λ * \hat{F} (g, V)$ . We refer to this ablation experiment as one step DPP-A over the complete grid cell code. The results on the analogy for this ablation experiment are shown in Appendix 1—figure 10. We see that the accuracy on test analogies for translation for various λs is around 30–60%, and for scaling around 20–40%, which is much lower than the nearly perfect accuracy achieved by the proposed DPP-A method. In the other case, we maximized the objective function $\hat{F} (g, V, N_{f}, N_{p}) = \sum_{f = 1}^{N_{f}} \log det (d i a g (σ (g_{f})) (V_{f} - I) + I)$ , using stochastic gradient ascent, which is same as $L_{D P P}$ in the ﬁrst step of Algorithm 1, along with minimizing $L_{t a s k}$ , which would use all the attended grid cell embeddings. So total loss, $L = L_{t a s k} - λ * \hat{F} (g, V, N_{f}, N_{p})$ . We refer to this ablation experiment as one step DPP-A within frequencies. As shown in Appendix 1—figure 11, the accuracy on test analogies for both translation and scaling for various λs is in a similar range to one step DPP-A over the complete grid cell code, and is much lower than the nearly perfect accuracy achieved by the proposed DPP-A method.

DPP-A attentional modulation

Appendix 1—figure 12

Download asset Open asset

Approximate maximum log determinant of the covariance matrix over the grid cell embeddings (y-axis) for each frequency (x-axis), obtained after maximizing Equation 6.

The gates within each frequency were optimized (independent of the task inputs), according to Equation 6, to compute the approximate maximum within-frequency log determinant of the covariance matrix over the grid cell embeddings individually for each frequency. We then used the grid cell embeddings belonging to the frequency that had the maximum within-frequency log determinant for training the inference module, which always happened to be grid cells within the top three frequencies. Appendix 1—figure 12 shows the approximate maximum log determinant (on the y-axis) for the diﬀerent frequencies (on the x-axis). The intuition behind why DPP-A identiﬁed grid cell embeddings corresponding to the highest spatial frequencies, and why this produced the best OOD generalization (i.e., extrapolation on our analogy tasks), is because those grid cell embeddings exhibited greater variance over the training data than the lower-frequency embeddings, while at the same time the correlations among those grid cell embeddings were lower than the correlations among the lower-frequency grid cell embeddings. The determinant of the covariance matrix of the grid cell embeddings is maximized when the variances of the grid cell embeddings are high (they are ‘expressive’) and the correlation among the grid cell embeddings is low (they ‘cover the representational space’). As a result, the higher-frequency grid cell embeddings more eﬃciently covered the representational space of the training data, allowing them to eﬃciently capture the same relational structure across training and test distributions which is required for OOD generalization.

To demonstrate that the higher grid cell frequencies more eﬃciently cover the representational space, we analyzed in Appendix 1—figure 13, the results after the summation of the multiplication of the grid cell embeddings over the 2D space of 1000 × 1000 locations, with their corresponding gates for three representative frequencies (left, middle, and right panels showing results for the lowest, middle, and highest grid cell frequencies, respectively, of the nine used in the model), obtained after maximizing Equation 6 for each grid cell frequency. The color code indicates the responsiveness of the grid cells to diﬀerent X and Y locations in the input space (lighter color corresponding to greater responsiveness). Note that the dark blue area (denoting regions of least responsiveness to any grid cell) is greatest for the lowest frequency and nearly zero for the highest frequency, illustrating that grid cell embeddings belonging to the highest frequency more eﬃciently cover the representational space which allows them to capture the same relational structure across training and test distributions as required for OOD generalization.

Appendix 1—figure 13

Download asset Open asset

Each panel shows the results after summation of the multiplication of the grid cell embeddings over the 2D space of 1000 × 1000 locations, with their corresponding gates for a particular frequency, obtained after maximizing Equation 6 for each grid cell frequency.

The left, middle, and right panels show results for the lowest, middle, and highest grid cell frequencies, respectively, of the nine used in the model. Lighter color in each panel corresponds to greater responsiveness of grid cells at that particular location in the 2D space.

Data availability

All data generated during this study are uploaded as source data files. Modeling code is uploaded as Source code 1.

References

(2017) Mapping of a non-spatial dimension by the hippocampal-entorhinal circuit
Nature 543:719–722.

https://doi.org/10.1038/nature21692
- PubMed
- Google Scholar
Preprint
1. Ba JL
2. Kiros JR
3. Hinton GE
(2016) Layer Normalization
arXiv.

https://arxiv.org/abs/1607.06450
- Google Scholar
1. Banino A
2. Barry C
3. Uria B
4. Blundell C
5. Lillicrap T
6. Mirowski P
7. Pritzel A
8. Chadwick MJ
9. Degris T
10. Modayil J
11. Wayne G
12. Soyer H
13. Viola F
14. Zhang B
15. Goroshin R
16. Rabinowitz N
17. Pascanu R
18. Beattie C
19. Petersen S
20. Sadik A
21. Gaffney S
22. King H
23. Kavukcuoglu K
24. Hassabis D
25. Hadsell R
26. Kumaran D
(2018) Vector-based navigation using grid-like representations in artificial agents
Nature 557:429–433.

https://doi.org/10.1038/s41586-018-0102-6
- PubMed
- Google Scholar
1. Bao X
2. Gjorgieva E
3. Shanahan LK
4. Howard JD
5. Kahnt T
6. Gottfried JA
(2019) Grid-like neural representations support olfactory navigation of a two-dimensional odor space
Neuron 102:1066–1075.

https://doi.org/10.1016/j.neuron.2019.03.034
- PubMed
- Google Scholar
Conference
1. Barrett D
2. Hill F
3. Santoro A
4. Morcos A
5. Lillicrap T
(2018)
Measuring abstract reasoning in neural networks

Proceedings of the 35th International Conference on Machine Learning. pp. 511–520.
- Google Scholar
1. Barry C
2. Lever C
3. Hayman R
4. Hartley T
5. Burton S
6. O’Keefe J
7. Jeffery K
8. Burgess N
(2006) The boundary vector cell model of place cell firing and spatial memory
Reviews in the Neurosciences 17:71–97.

https://doi.org/10.1515/revneuro.2006.17.1-2.71
- PubMed
- Google Scholar
1. Barry C
2. Hayman R
3. Burgess N
4. Jeffery KJ
(2007) Experience-dependent rescaling of entorhinal grids
Nature Neuroscience 10:682–684.

https://doi.org/10.1038/nn1905
- PubMed
- Google Scholar
Software
1. Bicanski A
(2019) VisualGridsRecognitionMem, version 9d56e64
GitHub.

https://github.com/bicanski/VisualGridsRecognitionMem
1. Bicanski A
2. Burgess N
(2019) A computational model of visual recognition memory via grid cells
Current Biology 29:979–990.

https://doi.org/10.1016/j.cub.2019.01.077
- PubMed
- Google Scholar
Conference
(2022)
Biologically-plausible determinant maximization neural networks for blind separation of correlated sources

Advances in Neural Information Processing Systems. pp. 13704–13717.
- Google Scholar
(2011) Reduction of theta rhythm dissociates grid cell spatial periodicity from directional tuning
Science 332:595–599.

https://doi.org/10.1126/science.1201652
- PubMed
- Google Scholar
1. Braver TS
2. Cohen JD
(2000)
On the control of control: The role of dopamine in regulating prefrontal function and working memory

Attention and Performance 18:712–737.
- Google Scholar
1. Bush D
2. Barry C
3. Manson D
4. Burgess N
(2015) Using grid cells for navigation
Neuron 87:507–520.

https://doi.org/10.1016/j.neuron.2015.07.006
- PubMed
- Google Scholar
Preprint
(2023) High-Capacity FLexible Hippocampal Associative and Episodic Memory Enabled by Prestructured”spatial”representations
bioRxiv.

https://doi.org/10.1101/2023.11.28.568960
- Google Scholar
Conference
1. Chen L
2. Zhang G
3. Zhou H
(2018)
Fast greedy map inference for determinantal point process to improve recommendation diversity

Proceedings of the 32nd International Conference on Neural Information Processing Systems. pp. 5627–5638.
- Google Scholar
(2001) Rostrolateral prefrontal cortex involvement in relational integration during reasoning
NeuroImage 14:1136–1149.

https://doi.org/10.1006/nimg.2001.0922
- PubMed
- Google Scholar
(2016) Organizing conceptual knowledge in humans with a gridlike code
Science 352:1464–1468.

https://doi.org/10.1126/science.aaf0941
- PubMed
- Google Scholar
Preprint
1. Cueva CJ
2. Wei XX
(2018) Emergence of Grid-like Representations by Training Recurrent Neural Networks to Perform Spatial Localization
arXiv.

https://arxiv.org/abs/1803.07770
- Google Scholar
1. Dayan P
(1993) Improving generalization for temporal difference learning: The successor representation
Neural Computation 5:613–624.

https://doi.org/10.1162/neco.1993.5.4.613
- Google Scholar
Preprint
1. Devlin J
2. Chang MW
3. Lee K
4. Toutanova K
(2018) Bert: Pre-Training of Deep Bidirectional Transformers for Language Understanding
arXiv.

https://arxiv.org/abs/1810.04805
- Google Scholar
(2010) Evidence for grid cells in a human memory network
Nature 463:657–661.

https://doi.org/10.1038/nature08704
- PubMed
- Google Scholar
1. Dordek Y
2. Soudry D
3. Meir R
4. Derdikman D
(2016) Extracting grid cell characteristics from place cell inputs using non-negative principal component analysis
eLife 5:e10094.

https://doi.org/10.7554/eLife.10094
- PubMed
- Google Scholar
1. Erdem UM
2. Hasselmo ME
(2014) A biologically inspired hierarchical goal directed navigation model
Journal of Physiology, Paris 108:28–37.

https://doi.org/10.1016/j.jphysparis.2013.07.002
- PubMed
- Google Scholar
(2001) Interactions between frontal cortex and basal ganglia in working memory: A computational model
Cognitive, Affective, & Behavioral Neuroscience 1:137–160.

https://doi.org/10.3758/CABN.1.2.137
- Google Scholar
Conference
(2019)
Extracting and utilizing abstract, structured representations for analogy

Proceedings of the 41st Annual Meeting of the Cognitive Science Society: Creativity + Cognition + Computation, CogSci 2019. pp. 1766–1772.
- Google Scholar
Conference
1. Frankland S
2. Cohen J
(2020)
Determinantal point processes for memory and structured inference

42nd Annual Meeting of the Cognitive Science Society: Developing a Mind: Learning in Humans, Animals, and Machines, CogSci 2020. pp. 3302–3308.
- Google Scholar
Conference
1. Gal Y
2. Ghahramani Z
(2016)
A theoretically grounded application of dropout in recurrent neural networks

Advances in Neural Information Processing Systems. pp. 1019–1027.
- Google Scholar
Conference
(2012)
Near-optimal map inference for determinantal point processes

Advances in Neural Information Processing Systems. pp. 2744–2752.
- Google Scholar
1. Giocomo LM
2. Hussaini SA
3. Zheng F
4. Kandel ER
5. Moser M-B
6. Moser EI
(2011) Grid cells use HCN1 channels for spatial scaling
Cell 147:1159–1170.

https://doi.org/10.1016/j.cell.2011.08.051
- PubMed
- Google Scholar
Conference
1. Gong B
2. Chao W-L
3. Grauman K
4. Sha F
(2014)
Diverse sequential subset selection for supervised video summarization

Advances in Neural Information Processing Systems. pp. 2069–2077.
- Google Scholar
1. Hafting T
2. Fyhn M
3. Molden S
4. Moser M-B
5. Moser EI
(2005) Microstructure of a spatial map in the entorhinal cortex
Nature 436:801–806.

https://doi.org/10.1038/nature03721
- PubMed
- Google Scholar
Software
1. Han I
(2017) Faster greedy MAP inference for determinantal point processes, version db7a28c
GitHub.

https://github.com/insuhan/fastdppmap
Conference
1. He K
2. Zhang X
3. Ren S
4. Sun J
(2016) Deep residual learning for image recognition
Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. pp. 770–778.

https://doi.org/10.1109/CVPR.2016.90
- Google Scholar
Conference
1. He K
2. Gkioxari G
3. Dollár P
4. Girshick R
(2017) Mask R-CNN
Proceedings of the IEEE International Conference on Computer Vision. pp. 2961–2969.

https://doi.org/10.1109/ICCV.2017.322
- Google Scholar
Preprint
1. Hill F
2. Santoro A
3. Barrett DG
4. Morcos AS
5. Lillicrap T
(2019) Learning to Make Analogies by Contrasting Abstract Relational Structure
arXiv.

https://arxiv.org/abs/1902.00120
- Google Scholar
(2019) Neuronal representation of environmental boundaries in egocentric coordinates
Nature Communications 10:2772.

https://doi.org/10.1038/s41467-019-10722-y
- PubMed
- Google Scholar
Conference
1. Hochreiter S
2. Schmidhuber J
(1997)
Lstm can solve hard long time lag problems

Advances in Neural Information Processing Systems. pp. 473–479.
- Google Scholar
Book
1. Holyoak KJ
(2012) Analogy and relational reasoning
In: Holyoak KJ, Morrison RG, editors. The Oxford Handbook of Thinking and Reasoning. Oxford University Press. pp. 234–259.

https://doi.org/10.1093/oxfordhb/9780199734689.013.0013
- Google Scholar
(2005) The temporal context model in spatial navigation and relational learning: toward a common explanation of medial temporal lobe function across domains
Psychological Review 112:75–116.

https://doi.org/10.1037/0033-295X.112.1.75
- PubMed
- Google Scholar
Conference
1. Ioﬀe S
2. Szegedy C
(2015)
Batch normalization: Accelerating deep network training by reducing internal covariate shift

International Conference on Machine Learning. pp. 448–456.
- Google Scholar
Website
1. Kazemnejad A
(2019) Transformer architecture: The positional encoding
Accessed September 20, 2019.

https://kazemnejad.com/blog/transformer_architecture_positional_encoding/
Preprint
1. Kingma DP
2. Ba J
(2014) Adam: A Method for Stochastic Optimization
arXiv.

https://arxiv.org/abs/1412.6980
- Google Scholar
(2012) A neurocomputational system for relational reasoning
Trends in Cognitive Sciences 16:373–381.

https://doi.org/10.1016/j.tics.2012.06.002
- Google Scholar
1. Ko C-W
2. Lee J
3. Queyranne M
(1995) An exact algorithm for maximum entropy sampling
Operations Research 43:684–691.

https://doi.org/10.1287/opre.43.4.684
- Google Scholar
Conference
1. Krogh A
2. Hertz JA
(1992)
A simple weight decay can improve generalization

Advances in Neural Information Processing Systems. pp. 950–957.
- Google Scholar
Preprint
1. Kulesza A
2. Taskar B
(2012) Determinantal Point Processes for Machine Learning
arXiv.

https://arxiv.org/abs/1207.6083
- Google Scholar
Conference
1. Lake B
2. Baroni M
(2018)
Generalization without systematicity: On the compositional skills of sequence-to-sequence recurrent networks

International Conference on Machine Learning. pp. 2873–2882.
- Google Scholar
1. Lu H
2. Ichien N
3. Holyoak KJ
(2022) Probabilistic analogical mapping with semantic relation networks
Psychological Review 129:1078–1103.

https://doi.org/10.1037/rev0000358
- PubMed
- Google Scholar
1. Macchi O
(1975) The coincidence approach to stochastic point processes
Advances in Applied Probability 7:83–122.

https://doi.org/10.2307/1425855
- Google Scholar
Preprint
1. Mariet Z
2. Sra S
(2015) Diversity Networks: Neural Network Compression Using Determinantal Point Processes
arXiv.

https://arxiv.org/abs/1511.05077
- Google Scholar
Preprint
(2019) Dppnet: Approximating Determinantal Point Processes with Deep Networks
arXiv.

https://arxiv.org/abs/1901.02051
- Google Scholar
(2012) Optimal population codes for space: grid cells outperform place cells
Neural Computation 24:2280–2317.

https://doi.org/10.1162/NECO_a_00319
- PubMed
- Google Scholar
(1995) Why there are complementary learning systems in the hippocampus and neocortex: insights from the successes and failures of connectionist models of learning and memory
Psychological Review 102:419–457.

https://doi.org/10.1037/0033-295X.102.3.419
- PubMed
- Google Scholar
(2022) Compositional sequence generation in the entorhinal-hippocampal system
Entropy 24:1791.

https://doi.org/10.3390/e24121791
- PubMed
- Google Scholar
Preprint
1. Mikolov T
2. Chen K
3. Corrado G
4. Dean J
(2013) Eﬃcient Estimation of Word Representations in Vector Space
arXiv.

https://arxiv.org/abs/1301.3781
- Google Scholar
Software
1. Mondal SS
(2024) DPP-Attention_Grid-cell-code, version swh:1:rev:f7c2d97feddd5427992a25b4bedf4f7af453c845
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:254a2d2ddc55f91afd102943f115ab68c126f5ae;origin=https://github.com/Shanka123/DPP-Attention_Grid-Cell-Code;visit=swh:1:snp:061922fece7d2d628477865513968f27ea71c634;anchor=swh:1:rev:f7c2d97feddd5427992a25b4bedf4f7af453c845
(2015) Place cells, grid cells, and memory
Cold Spring Harbor Perspectives in Biology 7:a021808.

https://doi.org/10.1101/cshperspect.a021808
- PubMed
- Google Scholar
Conference
1. Paszke A
2. Gross S
3. Chintala S
4. Chanan G
5. Yang E
6. DeVito Z
7. Lin Z
8. Desmaison A
9. Antiga L
10. Lerer A
(2017)
Automatic differentiation in PyTorch

Advances in Neural Information Processing Systems.
- Google Scholar
1. Perez-Beltrachini L
2. Lapata M
(2021) Multi-document summarization with determinantal point process attention
Journal of Artificial Intelligence Research 71:371–399.

https://doi.org/10.1613/jair.1.12522
- Google Scholar
(2021) Approximate multiplication in young children prior to multiplication instruction
Journal of Experimental Child Psychology 207:105116.

https://doi.org/10.1016/j.jecp.2021.105116
- PubMed
- Google Scholar
1. Rumelhart DE
2. Abrahamson AA
(1973) A model for analogical reasoning
Cognitive Psychology 5:1–28.

https://doi.org/10.1016/0010-0285(73)90023-6
- Google Scholar
Preprint
(2019) Analysing Mathematical Reasoning Abilities of Neural Models
arXiv.

https://arxiv.org/abs/1904.01557
- Google Scholar
(2013) The expected value of control: an integrative theory of anterior cingulate cortex function
Neuron 79:217–240.

https://doi.org/10.1016/j.neuron.2013.07.007
- PubMed
- Google Scholar
1. Silver D
2. Schrittwieser J
3. Simonyan K
4. Antonoglou I
5. Huang A
6. Guez A
7. Hubert T
8. Baker L
9. Lai M
10. Bolton A
11. Chen Y
12. Lillicrap T
13. Hui F
14. Sifre L
15. van den Driessche G
16. Graepel T
17. Hassabis D
(2017) Mastering the game of Go without human knowledge
Nature 550:354–359.

https://doi.org/10.1038/nature24270
- PubMed
- Google Scholar
1. Sorscher B
2. Mel GC
3. Ocko SA
4. Giocomo LM
5. Ganguli S
(2023) A unified theory for the computational and mechanistic origins of grid cells
Neuron 111:121–137.

https://doi.org/10.1016/j.neuron.2022.10.003
- PubMed
- Google Scholar
1. Sreenivasan S
2. Fiete I
(2011) Grid cells generate an analog error-correcting code for singularly precise neural computation
Nature Neuroscience 14:1330–1337.

https://doi.org/10.1038/nn.2901
- PubMed
- Google Scholar
(2014)
Dropout: a simple way to prevent neural networks from overﬁtting

The Journal of Machine Learning Research 15:1929–1958.
- Google Scholar
(2017) The hippocampus as a predictive map
Nature Neuroscience 20:1643–1653.

https://doi.org/10.1038/nn.4650
- PubMed
- Google Scholar
1. Stensola H
2. Stensola T
3. Solstad T
4. Frøland K
5. Moser M-B
6. Moser EI
(2012) The entorhinal grid map is discretized
Nature 492:72–78.

https://doi.org/10.1038/nature11649
- PubMed
- Google Scholar
(2020) Structure learning and the posterior parietal cortex
Progress in Neurobiology 184:101717.

https://doi.org/10.1016/j.pneurobio.2019.101717
- PubMed
- Google Scholar
Conference
1. Vaswani A
2. Shazeer N
3. Parmar N
4. Uszkoreit J
5. Jones L
6. Gomez AN
7. Kaiser Ł
8. Polosukhin I
(2017)
Attention is all you need

Advances in Neural Information Processing Systems. pp. 5998–6008.
- Google Scholar
(1999) A system for relational reasoning in human prefrontal cortex
Psychological Science 10:119–125.

https://doi.org/10.1111/1467-9280.00118
- Google Scholar
Conference
1. Webb T
2. Dulberg Z
3. Frankland S
4. Petrov A
5. O’Reilly R
6. Cohen J
(2020)
Learning representations that support extrapolation

International Conference on Machine Learning. pp. 10136–10146.
- Google Scholar
1. Webb T
2. Fu S
3. Bihl T
4. Holyoak KJ
5. Lu H
(2023) Zero-shot visual reasoning through probabilistic analogical mapping
Nature Communications 14:5144.

https://doi.org/10.1038/s41467-023-40804-x
- PubMed
- Google Scholar
(2015) A principle of economy predicts the functional architecture of grid cells
eLife 4:e08362.

https://doi.org/10.7554/eLife.08362
- PubMed
- Google Scholar
1. Whittington JC
2. Muller TH
3. Mark S
4. Chen G
5. Barry C
6. Burgess N
7. Behrens TE
(2020) The tolman-eichenbaum machine: Unifying space and relational memory through generalization in the hippocampal formation
Cell 183:1249–1263.

https://doi.org/10.1016/j.cell.2020.10.024
- PubMed
- Google Scholar
Preprint
1. Wu Y
2. Schuster M
3. Chen Z
4. Le QV
5. Norouzi M
6. Macherey W
7. Krikun M
8. Cao Y
9. Gao Q
10. Macherey K
(2016) Google’s Neural Machine Translation System: Bridging the Gap between Human and Machine Translation
arXiv.

https://arxiv.org/abs/1609.08144
- Google Scholar
Preprint
(2017) Determinantal Point Processes for Mini-Batch Diversiﬁcation
arXiv.

https://arxiv.org/abs/1705.00607
- Google Scholar

Article and author information

Author details

Shanka Subhra Mondal

Department of Electrical and Computer Engineering, Princeton University, Princeton, United States

Contribution
Data curation, Formal analysis, Visualization, Methodology, Writing - original draft

For correspondence
smondal@princeton.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-1252-0129
Steven Frankland

Princeton Neuroscience Institute, Princeton University, Princeton, United States

Contribution
Conceptualization, Supervision, Writing – review and editing

Competing interests
No competing interests declared
Taylor W Webb

Department of Psychology, University of California, Los Angeles, Los Angeles, United States

Contribution
Conceptualization, Supervision, Writing – review and editing

Competing interests
No competing interests declared
Jonathan D Cohen

Princeton Neuroscience Institute, Princeton University, Princeton, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Project administration, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-2316-0763

Funding

Office of Naval Research

Shanka Subhra Mondal
Steven Frankland
Taylor W Webb
Jonathan D Cohen

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Version history

Preprint posted: May 28, 2023
Sent for peer review: June 17, 2023
Reviewed Preprint version 1: August 23, 2023
Reviewed Preprint version 2: March 21, 2024
Version of Record published: August 1, 2024

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.89911. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

722

views
90

downloads
0

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Article PDF

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Shanka Subhra Mondal
Steven Frankland
Taylor W Webb
Jonathan D Cohen

(2024)

Determinantal point process attention over grid cell code supports out of distribution generalization

eLife 12:RP89911.

https://doi.org/10.7554/eLife.89911.3

Share this article

Cite this article

Schematic of the overall framework.

Generation of test analogies from training analogies (region marked in blue) by: (a) translating both dimension values of A,B,C,D by the same amount; and (b) scaling both dimension values of A,B,C,D by the same amount.

Figure 2—source data 1

Results on analogy on each region for translation and scaling using LSTM in the inference module.

Results on analogy on each region for translation and scaling using the transformer in the inference module.

Results on arithmetic on each region using LSTM in the inference module.

Results on arithmetic on each region using the transformer in the inference module.

Results on analogy on each region using determinantal point process attention (DPP-A), an LSTM in the inference module, and different embeddings (grid cell code, one-hots, and smoothed one-hots passed through a learned encoder) for translation (left) and scaling (right).

Results on analogy on each region using LSTM in the inference module for choosing top K frequencies with F^f in Algorithm 1.

Results on analogy on each region for translation and scaling using the transformer in the inference module.

Results on arithmetic with different embeddings (with determinantal point process attention [DPP-A]) using LSTM in the inference module.

Results on arithmetic with different embeddings (without determinantal point process attention [DPP-A], temporal context normalization [TCN], L1 Regularization, or Dropout) using LSTM in the inference module.

Results on arithmetic for increasing number of grid cell frequencies Nf on each region using LSTM in the inference module.

Results for regression on analogy using LSTM in the inference module.

Results for regression on arithmetic on each region using LSTM in the inference module.

Results on analogy for L1 regularization for various λs for translation and scaling using LSTM in the inference module.

Results on arithmetic for L1 regularization for various λs using LSTM in the inference module.

Results on analogy for one step determinantal point process attention (DPP-A) over the complete grid cell code for various λs for translation and scaling using LSTM in the inference module.

Results on analogy for one step determinantal point process attention (DPP-A) within frequencies for various λs for translation and scaling using LSTM in the inference module.

Approximate maximum log determinant of the covariance matrix over the grid cell embeddings (y-axis) for each frequency (x-axis), obtained after maximizing Equation 6.

Each panel shows the results after summation of the multiplication of the grid cell embeddings over the 2D space of 1000 × 1000 locations, with their corresponding gates for a particular frequency, obtained after maximizing Equation 6 for each grid cell frequency.

Author details

Shanka Subhra Mondal

Contribution

For correspondence

Competing interests

Steven Frankland

Contribution

Competing interests

Taylor W Webb

Contribution

Competing interests

Jonathan D Cohen

Contribution

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

Generation of test analogies from training analogies (region marked in blue) by: (a) translating both dimension values of $A, B, C, D$ by the same amount; and (b) scaling both dimension values of $A, B, C, D$ by the same amount.

Results on analogy on each region using LSTM in the inference module for choosing top $K$ frequencies with ${\hat{F}}_{f}$ in Algorithm 1.

Results on arithmetic for increasing number of grid cell frequencies $N_{f}$ on each region using LSTM in the inference module.