Source author record

Zachary P. Kilpatrick

Zachary P. Kilpatrick appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2024arXiv

Fast decisions reflect biases, slow decisions do not

Decisions are often made by heterogeneous groups of individuals, each with distinct initial biases and access to information of different quality. We show that in large groups of independent agents who accumulate evidence the first to decide are those with the strongest initial biases. Their decisions align with their initial bias, regardless of the underlying truth. In contrast, agents who decide last make decisions as if they were initially unbiased, and hence make better choices. We obtain asymptotic expressions in the large population limit that quantify how agents' initial inclinations shape early decisions. Our analysis shows how bias, information quality, and decision order interact in non-trivial ways to determine the reliability of decisions in a group.

preprint2016arXiv

Persistent search in confined domains: a velocity-jump process model

We analyze velocity-jump process models of persistent search for a single target on a bounded domain. The searcher proceeds along ballistic trajectories and is absorbed upon collision with the target boundary. When reaching the domain boundary, the searcher chooses a random direction for its new trajectory. For circular domains and targets, we can approximate the mean first passage time (MFPT) using a Markov chain approximation of the search process. Our analysis and numerical simulations reveal that the time to find the target decreases for targets closer to the domain boundary. When there is a small probability of direction-switching within the domain, we find the time to find the target decreases slightly with the turning probability. We also extend our exit time analysis to the case of partitioned domains, where there is a single target within one of multiple disjoint subdomains. Given an average time of transition between domains $\langle T \rangle$, we find that the optimal rate of transition that minimizes the time to find the target obeys $β_{\rm min} \propto 1/ \sqrt{\langle T \rangle}$.

preprint2015arXiv

Entrainment in up and down states of neural populations: non-smooth and stochastic models

We study the impact of noise on a neural population rate model of up and down states. Up and down states are typically observed in neuronal networks as a slow oscillation, where the population switches between high and low firing rates (Sanchez-Vivez and McCormick, 2000). A neural population model with spike rate adaptation is used to model such slow oscillations, and the timescale of adaptation determines the oscillation period. Furthermore, the period depends non-monotonically on the background tonic input driving the population, having long periods for very weak and very strong stimuli. Using both linearization and fast-slow timescale separation methods, we can compute the phase sensitivity function of the slow oscillation. We find that the phase response is most strongly impacted by perturbations to the adaptation variable. Phase sensitivity functions can then be utilized to quantify the impact of noise on oscillating populations. Noise alters the period of oscillations by speeding up the rate of transition between the up and down states. When common noise is presented to two distinct populations, their transitions will eventually become entrained to one another through stochastic synchrony.

preprint2015arXiv

Ghosts of bump attractors in stochastic neural fields: Bottlenecks and extinction

We study the effects of additive noise on stationary bump solutions to spatially extended neural fields near a saddle-node bifurcation. The integral terms of these evolution equations have a weight kernel describing synaptic interactions between neurons at different locations of the network. Excited regions of the neural field correspond to parts of the domain whose fraction of active neurons exceeds a sharp threshold of a firing rate nonlinearity. For sufficiently low firing threshold, a stable bump coexists with an unstable bump and a homogeneous quiescent state. As the threshold is increased, the stable and unstable branch of bump solutions annihilate in a saddle node bifurcation. Near this criticality, we derive a quadratic amplitude equation that describes the slow evolution of the even mode (bump contractions) as it depends on the distance from the bifurcation. Beyond the bifurcation, bumps eventually become extinct, and the time it takes for this to occur increases for systems nearer the bifurcation. When noise is incorporated, a stochastic amplitude equation for the even mode can be derived, which can be analyzed to reveal bump extinction time both below and above the saddle-node.

preprint2015arXiv

Networks that learn the precise timing of event sequences

Neuronal circuits can learn and replay firing patterns evoked by sequences of sensory stimuli. After training, a brief cue can trigger a spatiotemporal pattern of neural activity similar to that evoked by a learned stimulus sequence. Network models show that such sequence learning can occur through the shaping of feedforward excitatory connectivity via long term plasticity. Previous models describe how event order can be learned, but they typically do not explain how precise timing can be recalled. We propose a mechanism for learning both the order and precise timing of event sequences. In our recurrent network model, long term plasticity leads to the learning of the sequence, while short term facilitation enables temporally precise replay of events. Learned synaptic weights between populations determine the time necessary for one population to activate another. Long term plasticity adjusts these weights so that the trained event times are matched during playback. While we chose short term facilitation as a time-tracking process, we also demonstrate that other mechanisms, such as spike rate adaptation, can fulfill this role. We also analyze the impact of trial-to-trial variability, showing how observational errors as well as neuronal noise result in variability in learned event times. The dynamics of the playback process determine how stochasticity is inherited in learned sequence timings. Future experiments that characterize such variability can therefore shed light on the neural mechanisms of sequence learning.

preprint2015arXiv

Nonlinear Langevin equations for wandering patterns in stochastic neural fields

We analyze the effects of additive, spatially extended noise on spatiotemporal patterns in continuum neural fields. Our main focus is how fluctuations impact patterns when they are weakly coupled to an external stimulus or another equivalent pattern. Showing the generality of our approach, we study both propagating fronts and stationary bumps. Using a separation of time scales, we represent the effects of noise in terms of a phase-shift of a pattern from its uniformly translating position at long time scales, and fluctuations in the pattern profile around its instantaneous position at short time scales. In the case of a stimulus-locked front, we show that the phase-shift satisfies a nonlinear Langevin equation (SDE) whose deterministic part has a unique stable fixed point. Using a linear-noise approximation, we thus establish that wandering of the front about the stimulus-locked state is given by an Ornstein-Uhlenbeck (OU) process. Analogous results hold for the relative phase-shift between a pair of mutually coupled fronts, provided that the coupling is excitatory. On the other hand, if the mutual coupling is given by a Mexican hat function (difference of exponentials), then the linear-noise approximation breaks down due to the co-existence of stable and unstable phase-locked states in the deterministic limit. Similarly, the stochastic motion of mutually coupled bumps can be described by a system of nonlinearly coupled SDEs, which can be linearized to yield a multivariate OU process. As in the case of fronts, large deviations can cause bumps to temporarily decouple, leading to a phase-slip in the bump positions.

preprint2015arXiv

Stochastic models of evidence accumulation in changing environments

Organisms and ecological groups accumulate evidence to make decisions. Classic experiments and theoretical studies have explored this process when the correct choice is fixed during each trial. However, we live in a constantly changing world. What effect does such impermanence have on classical results about decision making? To address this question we use sequential analysis to derive a tractable model of evidence accumulation when the correct option changes in time. Our analysis shows that ideal observers discount prior evidence at a rate determined by the volatility of the environment, and the dynamics of evidence accumulation is governed by the information gained over an average environmental epoch. A plausible neural implementation of an optimal observer in a changing environment shows that, in contrast to previous models, neural populations representing alternate choices are coupled through excitation. Our work builds a bridge between statistical decision making in volatile environments and stochastic nonlinear dynamics.

preprint2015arXiv

Stochastic motion of bumps in planar neural fields

We analyze the effects of spatiotemporal noise on stationary pulse solutions (bumps) in neural field equations on planar domains. Neural fields are integrodifferential equations whose integral kernel describes the strength and polarity of synaptic interactions between neurons at different spatial locations of the network. Fluctuations in neural activity are incorporated by modeling the system as a Langevin equation evolving on a planar domain. Noise causes bumps to wander about the domain in a purely diffusive way. Utilizing a small noise expansion along with a solvability condition, we can derive an effective stochastic equation describing the bump dynamics as two-dimensional Brownian motion. The diffusion coefficient can then be computed explicitly. We also show that weak external inputs can pin the bump so it no longer wanders diffusively. Inputs reshape the effective potential that guides the dynamics of the bump position, so it tends to lie near attractors which can be single points or contours in the plane. Perturbative analysis shows the bump position evolves as a multivariate Ornstein-Uhlenbeck process whose relaxation constants are determined by the shape of the input. Our analytical approximations all compare well to statistics of bump motion in numerical simulations.

preprint2015arXiv

Stochastic synchronization of neural activity waves

We demonstrate that waves in distinct layers of a neuronal network can become phase-locked by common spatiotemporal noise. This phenomenon is studied for stationary bumps, traveling waves, and breathers. A weak noise expansion is used to derive an effective equation for the position of the wave in each layer, yielding a stochastic differential equation with multiplicative noise. Stability of the synchronous state is characterized by a Lyapunov exponent, which we can compute analytically from the reduced system. Our results extend previous work on limit-cycle oscillators, showing common noise can synchronize waves in a broad class of models.

preprint2014arXiv

Delay stabilizes stochastic motion of bumps in layered neural fields

We study the effects of propagation delays on the stochastic dynamics of bumps in neural fields with multiple layers. In the absence of noise, each layer supports a stationary bump. Using linear stability analysis, we show that delayed coupling between layers causes translating perturbations of the bumps to decay in the noise-free system. Adding noise to the system causes bumps to wander as a random walk. However, coupling between layers can reduce the variability of this stochastic motion by canceling noise that perturbs bumps in opposite directions. Delays in interlaminar coupling can further reduce variability, since they couple bump positions to states from the past. We demonstrate these relationships by deriving an asymptotic approximation for the effective motion of bumps. This yields a stochastic delay-differential equation where each delayed term arises from an interlaminar coupling. The impact of delays is well approximated by using a small delay expansion, which allows us to compute the effective diffusion in bumps' positions, accurately matching results from numerical simulations.

preprint2014arXiv

Molecular mechanisms that regulate the coupled period of the mammalian circadian clock

In mammals, most cells in the brain and peripheral tissues generate circadian (~24hr) rhythms autonomously. These self-sustained rhythms are coordinated and entrained by a master circadian clock in the suprachiasmatic nucleus (SCN). Within the SCN, the individual rhythms of each neuron are synchronized through intercellular signaling. One important feature of SCN is that the synchronized period is close to the cell population mean of intrinsic periods. In this way, the synchronized period of the SCN stays close to the periods of cells in peripheral tissues. This is important for SCN to entrain cells throughout the body. However, the mechanism that drives the period of the coupled SCN cells to the population mean is not known. We use mathematical modeling and analysis to show that the mechanism of transcription repression plays a pivotal role in regulating the coupled period. Specifically, we use phase response curve analysis to show that the coupled period within the SCN stays near the population mean if transcriptional repression occurs via protein sequestration. In contrast, the coupled period is far from the mean if repression occurs through highly nonlinear Hill-type regulation (e.g. oligomer- or phosphorylation-based repression). Furthermore, we find that the timescale of intercellular coupling needs to be fast compared to that of intracellular feedback to maintain the mean period. These findings reveal the important relationship between the intracellular transcriptional feedback loop and intercellular coupling. This relationship explains why transcriptional repression appears to occur via protein sequestration in multicellular organisms, mammals and Drosophila, in contrast with the phosphorylation-based repression in unicellular organisms. That is, transition to protein sequestration is essential for synchronizing multiple cells with a period close to the population mean (~24hr).

preprint2013arXiv

Coupling layers regularizes wave propagation in laminar stochastic neural fields

We study the effects of coupling between layers of stochastic neural field models with laminar structure. In particular, we focus on how the propagation of waves of neural activity in each layer is affected by the coupling. Synaptic connectivities within and between each layer are determined by integral kernels of an integrodifferential equation describing the temporal evolution of neural activity. Excitatory neural fields, with purely positive connectivities, support traveling fronts in each layer, whose speeds are increased when coupling between layers is considered. Studying the effects of noise, we find coupling also serves to reduce the variance in the position of traveling fronts, as long as the noise sources to each layer are not completely correlated. Neural fields with asymmetric connectivity support traveling pulses whose speeds are decreased by interlaminar coupling. Again, coupling reduces the variance in traveling pulse position, when noise is considered that is not totally correlated between layers. To derive our stochastic results, we employ a small-noise expansion, also assuming interlaminar connectivity scales similarly. Our asymptotic results agree reasonably with accompanying numerical simulations.

preprint2013arXiv

Optimizing working memory with heterogeneity of recurrent cortical excitation

A neural correlate of parametric working memory is a stimulus specific rise in neuron firing rate that persists long after the stimulus is removed. Network models with local excitation and broad inhibition support persistent neural activity, linking network architecture and parametric working memory. Cortical neurons receive noisy input fluctuations which causes persistent activity to diffusively wander about the network, degrading memory over time. We explore how cortical architecture that supports parametric working memory affects the diffusion of persistent neural activity. Studying both a spiking network and a simplified potential well model, we show that spatially heterogeneous excitatory coupling stabilizes a discrete number of persistent states, reducing the diffusion of persistent activity over the network. However, heterogeneous coupling also coarse-grains the stimulus representation space, limiting the capacity of parametric working memory. The storage errors due to coarse-graining and diffusion tradeoff so that information transfer between the initial and recalled stimulus is optimized at a fixed network heterogeneity. For sufficiently long delay times, the optimal number of attractors is less than the number of possible stimuli, suggesting that memory networks can under-represent stimulus space to optimize performance. Our results clearly demonstrate the effects of network architecture and stochastic fluctuations on parametric memory storage.

preprint2012arXiv

Wandering bumps in stochastic neural fields

We study the effects of noise on stationary pulse solutions (bumps) in spatially extended neural fields. The dynamics of a neural field is described by an integrodifferential equation whose integral term characterizes synaptic interactions between neurons in different spatial locations of the network. Translationally symmetric neural fields support a continuum of stationary bump solutions, which may be centered at any spatial location. Random fluctuations are introduced by modeling the system as a spatially extended Langevin equation whose noise term we take to be multiplicative or additive. For nonzero noise, these bumps are shown to wander about the domain in a purely diffusive way. We can approximate the effective diffusion coefficient using a small noise expansion. Upon breaking the (continuous) translation symmetry of the system using a spatially heterogeneous inputs or synapses, bumps in the stochastic neural field can become temporarily pinned to a finite number of locations in the network. In the case of spatially heterogeneous synaptic weights, as the modulation frequency of this heterogeneity increases, the effective diffusion of bumps in the network approaches that of the network with spatially homogeneous weights.