BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Paul C. Bressloff

Paul C. Bressloff contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechQuantitative MethodsNeurons and CognitionCell Behaviorcond-mat.dis-nnmath.APnlin.AOnlin.PSphysics.chem-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
13works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2023arXiv

Renewal equations for single-particle diffusion in multi-layered media

In this paper we develop a probabilistic model of single-particle diffusion in 1D multi-layered media by constructing a multi-layered version of so-called snapping out Brownian motion (BM). The latter sews together successive rounds of reflected BM, each of which is restricted to a single layer. Each round of reflected BM is killed when the local time at one end of the layer exceeds an independent, exponentially distributed random variable. (The local time specifies the amount of time a reflected Brownian particle spends in a neighborhood of a boundary.) The particle then immediately resumes reflected BM in the same layer or the layer on the other side of the boundary with equal probability, and the process is iterated We proceed by constructing a last renewal equation for multi-layered snapping out BM that relates the full probability density to the probability densities of partially reflected BM in each layer. We then show how transfer matrices can be used to solve the Laplace transformed renewal equation, and prove that the renewal equation and corresponding multi-layer diffusion equation are equivalent. We illustrate the theory by analyzing the first passage time (FPT) problem for escape at the exterior boundaries of the domain. Finally, we use the renewal approach to incorporate a generalization of snapping out BM based on the encounter-based method for surface absorption; each round of reflected BM is now killed according to a non-exponential distribution for each local time threshold. This is achieved by considering a corresponding first renewal equation that relates the full probability density to the FPT densities for killing each round of reflected BM. We show that for certain configurations, non-exponential killing leads to an effective time-dependent permeability that is normalizable but heavy-tailed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Diffusion-mediated surface reactions and stochastic resetting

In this paper, we investigate the effects of stochastic resetting on diffusion in $\R^d\backslash \calU$, where $\calU$ is a bounded obstacle with a partially absorbing surface $\partial \calU$. We begin by considering a Robin boundary condition with a constant reactivity $κ_0$, and show how previous results are recovered in the limits $κ_0\rightarrow 0,\infty$. We then generalize the Robin boundary condition to a more general probabilistic model of diffusion-mediated surface reactions using an encounter-based approach. The latter considers the joint probability density or propagator $P(\x,\ell,t|\x_0)$ for the pair $(\X_t,\ell_t)$ in the case of a perfectly reflecting surface, where $\X_t$ and $\ell_t$ denote the particle position and local time, respectively. The local time determines the amount of time that a Brownian particle spends in a neighborhood of the boundary. The effects of surface reactions are then incorporated via an appropriate stopping condition for the boundary local time. We construct the boundary value problem (BVP) satisfied by the propagator in the presence of resetting, and use this to derive implicit equations for the marginal density of particle position and the survival probability. We highlight the fact that these equations are difficult to solve in the case of non-constant reactivities, since resetting is not governed by a renewal process. We then consider a simpler problem in which both the position and local time are reset. In this case, the survival probability with resetting can be expressed in terms of the survival probability without resetting, which allows us to explore the dependence of the MFPT on the resetting rate $r$ and the type of surface reactions. The theory is illustrated using the example of a spherically symmetric surface.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Diffusion-mediated surface reactions, Brownian functionals and the Feynman-Kac formula

Many processes in cell biology involve diffusion in a domain $Ω$ that contains a target $\calU$ whose boundary $\partial \calU$ is a chemically reactive surface. Such a target could represent a single reactive molecule, an intracellular compartment or a whole cell. Recently, a probabilistic framework for studying diffusion-mediated surface reactions has been developed that considers the joint probability density or propagator for the particle position and the so-called boundary local time. The latter characterizes the amount of time that a Brownian particle spends in the neighborhood of a point on a totally reflecting boundary. The effects of surface reactions are then incorporated via an appropriate stopping condition for the boundary local time. In this paper we generalize the theory of diffusion-mediated surface reactions to cases where the whole interior target domain $\calU$ acts as a partial absorber rather than the target boundary $\partial \calU$. Now the particle can freely enter and exit $\calU$, and is only able to react (be absorbed) within $\calU$. The appropriate Brownian functional is then the occupation time (accumulated time that the particle spends within $\calU$) rather than the boundary local time. We show that both cases can be considered within a unified framework by using a Feynman-Kac formula to derive a boundary value problem (BVP) for the propagator of the corresponding Brownian functional, and introducing an associated stopping condition. We illustrate the theory by calculating the mean first passage time (MFPT) for a spherical target $\calU$ located at the center of a spherical domain $Ω$. This is achieved by solving the propagator BVP directly, rather than using spectral methods. We find that if the first moment of the stopping time density is infinite, then the MFPT is also infinite, that is, the spherical target is not sufficiently absorbing.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local accumulation time for diffusion in cells with gap junction coupling

In this paper we analyze the relaxation to steady-state of intracellular diffusion in a pair of cells with gap-junction coupling. Gap junctions are prevalent in most animal organs and tissues, providing a direct diffusion pathway for both electrical and chemical communication between cells. Most analytical models of gap junctions focus on the steady-state diffusive flux and the associated effective diffusivity. Here we investigate the relaxation to steady state in terms of the so-called local accumulation time. The latter is commonly used to estimate the time to form a protein concentration gradient during morphogenesis. The basic idea is to treat the fractional deviation from the steady-state concentration as a cumulative distribution for the local accumulation time. One of the useful features of the local accumulation time is that it takes into account the fact that different spatial regions can relax at different rates. We consider both static and dynamic gap junction models. The former treats the gap junction as a resistive channel with effective permeability $μ$, whereas the latter represents the gap junction as a stochastic gate that randomly switches between an open and closed state. The local accumulation time is calculated by solving the diffusion equation in Laplace space and then taking the small-$s$ limit. We show that the accumulation time is a monotonically increasing function of spatial position, with a jump discontinuity at the gap junction. This discontinuity vanishes in the limit $μ\rightarrow \infty$ for a static junction and $β\rightarrow 0$ for a stochastically-gated junction, where $β$ is the rate at which the gate closes. Finally, our results are generalized to the case of a linear array of cells with nearest neighbor gap junction coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local accumulation times in a diffusion-trapping model of synaptic receptor dynamics

The lateral diffusion and trapping of neurotransmitter receptors within the postsynaptic membrane of a neuron plays a key role in determining synaptic strength and plasticity. Trapping is mediated by the reversible binding of receptors to scaffolding proteins (slots) within a synapse. In this paper we introduce a new method for analyzing the transient dynamics of synapses in a diffusion-trapping model of receptor trafficking. Given a population of spatially distributed synapses, each of which has a fixed number of slots, we calculate the rate of relaxation to the steady-state distribution of bound slots (synaptic weights) in terms of a set of local accumulation times. Assuming that the rates of exocytosis and endocytosis are sufficiently slow, we show that the steady-state synaptic weights are independent of each other (purely local). On the other hand, the local accumulation time of a given synapse depends on the number of slots and the spatial location of all the synapses, indicating a form of transient heterosynaptic plasticity. This suggests that local accumulation time measurements could provide useful information regarding the distribution of synaptic weights within a dendrite.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The narrow capture problem: an encounter-based approach to partially reactive targets

A general topic of current interest is the analysis of diffusion problems in singularly perturbed domains with small interior targets or traps (the narrow capture problem). One major application is to intracellular diffusion, where the targets typically represent some form of reactive biochemical substrate. Most studies of the narrow capture problem treat the target boundaries as totally absorbing. In this paper, we analyze the three-dimensional narrow capture problem in the more realistic case of partially reactive target boundaries. We begin by considering classical Robin boundary conditions. Matching inner and outer solutions of the single-particle probability density, we derive an asymptotic expansion of the Laplace transformed flux into each reactive surface in powers of $ε$, where $ερ$ is a given target size. In turn, the fluxes determine the splitting probabilities for target absorption. We then extend our analysis to more general types of reactive targets by combining matched asymptotic analysis with an encounter-based formulation of diffusion-mediated surface reactions. That is, we derive an asymptotic expansion of the joint probability density for particle position and the boundary local time. The effects of surface reactions are then incorporated via an appropriate stopping condition for the boundary local time. Finally, we illustrate the theory by exploring how the leading-order contributions to the splitting probabilities depend on the choice of surface reactions. In particular, we show that there is an effective renormalization of the target radius of the form $ρ\rightarrow ρ-\widetildeΨ(1/ρ)$, where $\widetildeΨ$ is the Laplace transform of the stopping local time distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Threshold surface reactions and local time resetting

In this paper we consider a threshold surface absorption mechanism for a particle diffusing in a domain containing a single target $\calU $. The target boundary $\partial \calU$ is taken to be a reactive surface that modifies an internal state $Z_t$ of the particle when in contact with the surface at time $t$, with $Z_0=h$. The state $Z_t$ is taken to be a monotonically decreasing function of the so-called boundary local time, and absorption occurs as soon as $Z_t$ reaches zero. (The boundary local time is a Brownian functional that determines the amount of time that the particle spends in a neighborhood of $\partial \calU$.) We first show how to analyze threshold surface absorption in terms of the joint probability density or generalized propagator $P_0(\x,\ell,t|\x_0)$ for the pair $(\X_t,\ell_t)$ in the case of a perfectly reflecting surface, where $\X_t$ and $\ell_t$ denote the particle position and local time at time $t$, respectively, and $\x_0$ is the initial position. We then introduce a generalized stochastic resetting protocol in which both the position $\X_t$ and the internal state $Z_t$ are reset to their initial values, $\X_t\rightarrow \x_0$ and $Z_t\rightarrow h$, at a Poisson rate $r$. The latter is mathematically equivalent to resetting the boundary local time, $\ell_t\rightarrow 0$. Since resetting is governed by a renewal process, the survival probability with resetting can be expressed in terms of the survival probability without resetting, which means that the statistics of absorption can be determined by calculating the Laplace transform of $P_0(\x,\ell,t|\x_0)$ with respect to $t$. We contrast this with the case where only particle position is reset, which is not governed by a renewal process. We illustrate the theory using the simple examples of diffusion on the half-line and a spherical target in a spherical domain.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastic Turing pattern formation in a model with active and passive transport

We investigate Turing pattern formation in a stochastic and spatially discretized version of a reaction diffusion advection (RDA) equation, which was previously introduced to model synaptogenesis in \textit{C. elegans}. The model describes the interactions between a passively diffusing molecular species and an advecting species that switches between anterograde and retrograde motor-driven transport (bidirectional transport). Within the context of synaptogenesis, the diffusing molecules can be identified with the protein kinase CaMKII and the advecting molecules as glutamate receptors. The stochastic dynamics evolves according to an RDA master equation, in which advection and diffusion are both modeled as hopping reactions along a one-dimensional array of chemical compartments. Carrying out a linear noise approximation of the RDA master equation leads to an effective Langevin equation, whose power spectrum provides a means of extending the definition of a Turing instability to stochastic systems, namely, in terms of the existence of a peak in the power spectrum at a non-zero spatial frequency. We thus show how noise can significantly extend the range over which spontaneous patterns occur, which is consistent with previous studies of RD systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Coupling sample paths to the partial thermodynamic limit in stochastic chemical reaction networks

Many biochemical systems appearing in applications have a multiscale structure so that they converge to piecewise deterministic Markov processes in a thermodynamic limit. The statistics of the piecewise deterministic process can be obtained much more efficiently than those of the exact process. We explore the possibility of coupling sample paths of the exact model to the piecewise deterministic process in order to reduce the variance of their difference. We then apply this coupling to reduce the computational complexity of a Monte Carlo estimator. In addition to rigorous results concerning the asymptotic computational complexity of the Monte Carlo estimator, numerical simulations are performed on some simple biological models confirming that computational gains are made.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Breakdown of fast-slow analysis in an excitable system with channel noise

We consider a stochastic version of an excitable system based on the Morris-Lecar model of a neuron, in which the noise originates from stochastic Sodium and Potassium ion channels opening and closing. One can analyze neural excitability in the deterministic model by using a separation of time scales involving a fast voltage variable and a slow recovery variable, which represents the fraction of open Potassium channels. In the stochastic setting, spontaneous excitation is initiated by ion channel noise. If the recovery variable is constant during initiation, the spontaneous activity rate can be calculated using Kramer's rate theory. The validity of this assumption in the stochastic model is examined using a systematic perturbation analysis. We find that in most physically relevant cases, this assumption breaks down, requiring an alternative to Kramers theory for excitable systems with one deterministic fixed point. We also show that an exit time problem can be formulated in an excitable system by considering maximum likelihood trajectories of the stochastic process.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Metastability in a stochastic neural network modeled as a velocity jump Markov process

One of the major challenges in neuroscience is to determine how noise that is present at the molecular and cellular levels affects dynamics and information processing at the macroscopic level of synaptically coupled neuronal populations. Often noise is incorprated into deterministic network models using extrinsic noise sources. An alternative approach is to assume that noise arises intrinsically as a collective population effect, which has led to a master equation formulation of stochastic neural networks. In this paper we extend the master equation formulation by introducing a stochastic model of neural population dynamics in the form of a velocity jump Markov process. The latter has the advantage of keeping track of synaptic processing as well as spiking activity, and reduces to the neural master equation in a particular limit. The population synaptic variables evolve according to piecewise deterministic dynamics, which depends on population spiking activity. The latter is characterised by a set of discrete stochastic variables evolving according to a jump Markov process, with transition rates that depend on the synaptic variables. We consider the particular problem of rare transitions between metastable states of a network operating in a bistable regime in the deterministic limit. Assuming that the synaptic dynamics is much slower than the transitions between discrete spiking states, we use a WKB approximation and singular perturbation theory to determine the mean first passage time to cross the separatrix between the two metastable states. Such an analysis can also be applied to other velocity jump Markov processes, including stochastic voltage-gated ion channels and stochastic gene networks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hebbian learning of recurrent connections: a geometrical perspective

We show how a Hopfield network with modifiable recurrent connections undergoing slow Hebbian learning can extract the underlying geometry of an input space. First, we use a slow/fast analysis to derive an averaged system whose dynamics derives from an energy function and therefore always converges to equilibrium points. The equilibria reflect the correlation structure of the inputs, a global object extracted through local recurrent interactions only. Second, we use numerical methods to illustrate how learning extracts the hidden geometrical structure of the inputs. Indeed, multidimensional scaling methods make it possible to project the final connectivity matrix on to a distance matrix in a high-dimensional space, with the neurons labelled by spatial position within this space. The resulting network structure turns out to be roughly convolutional. The residual of the projection defines the non-convolutional part of the connectivity which is minimized in the process. Finally, we show how restricting the dimension of the space where the neurons live gives rise to patterns similar to cortical maps. We motivate this using an energy efficiency argument based on wire length minimization. Finally, we show how this approach leads to the emergence of ocular dominance or orientation columns in primary visual cortex. In addition, we establish that the non-convolutional (or long-range) connectivity is patchy, and is co-aligned in the case of orientation learning.

0 citations0 reviewsNo code linkedOpen access