BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Peter J. Thomas

Peter J. Thomas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Information Theorymath.ITNeurons and Cognitionmath.DSMolecular NetworksQuantitative MethodsBiological PhysicsCell Behaviorcond-mat.stat-mecheess.SPmath.PRnlin.CDOther Quantitative Biologyphysics.app-phphysics.ins-det

Catalog footprint

What is connected

15works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

15 published item(s)

preprint2023arXiv

Explicitly Solvable Continuous-time Inference for Partially Observed Markov Processes

Many natural and engineered systems can be modeled as discrete state Markov processes. Often, only a subset of states are directly observable. Inferring the conditional probability that a system occupies a particular hidden state, given the partial observation, is a problem with broad application. In this paper, we introduce a continuous-time formulation of the sum-product algorithm, which is a well-known discrete-time method for finding the hidden states' conditional probabilities, given a set of finite, discrete-time observations. From our new formulation, we can explicitly solve for the conditional probability of occupying any state, given the transition rates and observations within a finite time window. We apply our algorithm to a realistic model of the cystic fibrosis transmembrane conductance regulator (CFTR) protein for exact inference of the conditional occupancy probability, given a finite time series of partial observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Variational and phase response analysis for limit cycles with hard boundaries, with applications to neuromechanical control problems

Motor systems show an overall robustness, but because they are highly nonlinear, understanding how they achieve robustness is difficult. In many rhythmic systems, robustness against perturbations involves response of both the shape and the timing of the trajectory. This makes the study of robustness even more challenging. To understand how a motor system produces robust behaviors in a variable environment, we consider a neuromechanical model of motor patterns in the feeding apparatus of the marine mollusk \textit{Aplysia californica} \citep{shaw2015,lyttle2017}. We established in \citep{WGCT2021} the tools for studying combined shape and timing responses of limit cycle systems under sustained perturbations and here apply them to study robustness of the neuromechanical model against increased mechanical load during swallowing. Interestingly, we discover that nonlinear biomechanical properties confer resilience by immediately increasing resistance to applied loads. In contrast, the effect of changed sensory feedback signal is significantly delayed by the firing rates' hard boundary properties. Our analysis suggests that sensory feedback contributes to robustness in swallowing primarily by shifting the timing of neural activation involved in the power stroke of the motor cycle (retraction). This effect enables the system to generate stronger retractor muscle forces to compensate for the increased load, and hence achieve strong robustness. The approaches that we are applying to understanding a neuromechanical model in \textit{Aplysia}, and the results that we have obtained, are likely to provide insights into the function of other motor systems that encounter changing mechanical loads and hard boundaries, both due to mechanical and neuronal firing properties.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An analytical approach to the Mean--Return-Time~Phase of isotropic stochastic oscillators

One notion of phase for stochastic oscillators is based on the mean return-time (MRT): a set of points represents a certain phase if the mean time to return from any point in this set to this set after one rotation is equal to the mean rotation period of the oscillator (irrespective of the starting point). For this so far only algorithmically defined phase, we derive here analytical expressions for the important class of isotropic stochastic oscillators. This allows us to evaluate cases from the literature explicitly and to study the behavior of the MRT phase in the limits of strong noise. We also use the same formalism to show that lines of constant return time variance (instead of constant mean return time) can be defined, and that they in general differ from the MRT-isochrons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Trapped air metamaterials for ultrasonic sub-wavelength imaging in water

Acoustic metamaterials constructed from conventional base materials can exhibit exotic phenomena such as negative refractive index, extraordinary transmission-absorption and sub-wavelength imaging. These are typically achieved by combining geometrical and resonance effects. Holey-structured acoustic metamaterials have already shown potential for sub-wavelength acoustic imaging in air, where it is relatively simple to achieve the required difference in acoustic impedance needed to create resonances within the holes. However, the use of polymers for metamaterial operation in water is difficult, due to the much lower difference in acoustic impedance between water and the many polymers used, for example, in additive manufacturing. Hence, metals are commonly used. Here we show that this can be overcome by trapping air layers so that they surround each water-filled channel, making sub-wavelength imaging in water possible. The operation of such a Trapped air design is confirmed at 200-300 kHz ultrasonic frequencies via both finite element modelling and experimental measurements in a water tank. It is also shown quantitatively that the trapped air design outperforms its metal counterpart. The results indicate a way forward for exploiting additive-manufacturing for realising acoustic metamaterials in water at ultrasonic frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Capacity of a Simple Intercellular Signal Transduction Channel

We model biochemical signal transduction, based on a ligand-receptor binding mechanism, as a discrete-time finite-state Markov channel, which we call the BIND channel. We show how to obtain the capacity of this channel, for the case of binary output, binary channel state, and arbitrary finite input alphabets. We show that the capacity-achieving input distribution is IID. Further, we show that feedback does not increase the capacity of this channel. We show how the capacity of the discrete-time channel approaches the capacity of Kabanov's Poisson channel, in the limit of short time steps and rapid ligand release.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Entropy of the Sum of Two Independent, Non-Identically-Distributed Exponential Random Variables

In this letter, we give a concise, closed-form expression for the differential entropy of the sum of two independent, non-identically-distributed exponential random variables. The derivation is straightforward, but such a concise entropy has not been previously given in the literature. The usefulness of the expression is demonstrated with examples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite-State Channel Models for Signal Transduction in Neural Systems

Information theory provides powerful tools for understanding communication systems. This analysis can be applied to intercellular signal transduction, which is a means of chemical communication among cells and microbes. We discuss how to apply information-theoretic analysis to ligand-receptor systems, which form the signal carrier and receiver in intercellular signal transduction channels. We also discuss the applications of these results to neuroscience.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shannon Capacity of Signal Transduction for Multiple Independent Receptors

Cyclic adenosine monophosphate (cAMP) is considered a model system for signal transduction, the mechanism by which cells exchange chemical messages. Our previous work calculated the Shannon capacity of a single cAMP receptor; however, a typical cell may have thousands of receptors operating in parallel. In this paper, we calculate the capacity of a cAMP signal transduction system with an arbitrary number of independent, indistinguishable receptors. By leveraging prior results on feedback capacity for a single receptor, we show (somewhat unexpectedly) that the capacity is achieved by an IID input distribution, and that the capacity for n receptors is n times the capacity for a single receptor.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymptotic Phase for Stochastic Oscillators

Oscillations and noise are ubiquitous in physical and biological systems. When oscillations arise from a deterministic limit cycle, entrainment and synchronization may be analyzed in terms of the asymptotic phase function. In the presence of noise, the asymptotic phase is no longer well defined. We introduce a new definition of asymptotic phase in terms of the slowest decaying modes of the Kolmogorov backward operator. Our stochastic asymptotic phase is well defined for noisy oscillators, even when the oscillations are noise dependent. It reduces to the classical asymptotic phase in the limit of vanishing noise. The phase can be obtained either by solving an eigenvalue problem, or by empirical observation of an oscillating density's approach to its steady state.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stochastic Representations of Ion Channel Kinetics and Exact Stochastic Simulation of Neuronal Dynamics

In this paper we provide two representations for stochastic ion channel kinetics, and compare the performance of exact simulation with a commonly used numerical approximation strategy. The first representation we present is a random time change representation, popularized by Thomas Kurtz, with the second being analogous to a "Gillespie" representation. Exact stochastic algorithms are provided for the different representations, which are preferable to either (a) fixed time step or (b) piecewise constant propensity algorithms, which still appear in the literature. As examples, we provide versions of the exact algorithms for the Morris-Lecar conductance based model, and detail the error induced, both in a weak and a strong sense, by the use of approximate algorithms on this model. We include ready-to-use implementations of the random time change algorithm in both XPP and Matlab. Finally, through the consideration of parametric sensitivity analysis, we show how the representations presented here are useful in the development of further computational methods. The general representations and simulation strategies provided here are known in other parts of the sciences, but less so in the present setting.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Lower Bound for the First Passage Time Density of the Suprathreshold Ornstein-Uhlenbeck Process

We prove that the first passage time density $ρ(t)$ for an Ornstein-Uhlenbeck process $X(t)$ obeying $dX=-βX dt + σdW$ to reach a fixed threshold $θ$ from a suprathreshold initial condition $x_0>θ>0$ has a lower bound of the form $ρ(t)>k \exp\left[-p e^{6βt}\right]$ for positive constants $k$ and $p$ for times $t$ exceeding some positive value $u$. We obtain explicit expressions for $k, p$ and $u$ in terms of $β$, $σ$, $x_0$ and $θ$, and discuss application of the results to the synchronization of periodically forced stochastic leaky integrate-and-fire model neurons.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase Resetting in an Asymptotically Phaseless System: On the Phase Response of Limit Cycles Verging on a Heteroclinic Orbit

Rhythmic behaviors in neural systems often combine features of limit cycle dynamics (stability and periodicity) with features of near heteroclinic or near homoclinic cycle dynamics (extended dwell times in localized regions of phase space). Proximity of a limit cycle to one or more saddle equilibria can have a profound effect on the timing of trajectory components and response to both fast and slow perturbations, providing a possible mechanism for adaptive control of rhythmic motions. Reyn showed that for a planar dynamical system with a stable heteroclinic cycle (or separatrix polygon), small perturbations satisfying a net inflow condition will generically give rise to a stable limit cycle (Reyn, 1980; Guckenheimer and Holmes, 1983). Here we consider the asymptotic behavior of the infinitesimal phase response curve (iPRC) for examples of two systems satisfying Reyn's inflow criterion, (i) a smooth system with a chain of four hyperbolic saddle points and (ii) a piecewise linear system corresponding to local linearization of the smooth system about its saddle points. For system (ii), we obtain exact expressions for the limit cycle and the iPRC as a function of a parameter, mu>0, representing the distance from a heteroclinic bifurcation point. In the limit, as mu approaches zero, we find that perturbations parallel to the unstable eigenvector direction in a piecewise linear region lead to divergent phase response, as previously observed (Brown, Moehlis and Holmes (2004), Neural Computation). In contrast to previous work, we find that perturbations parallel to the stable eigenvector direction can lead to either divergent or convergent phase response, depending on the phase at which the perturbation occurs. In the smooth system (i), we show numerical evidence of qualitatively similar phase specific sensitivity to perturbation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Information Maximization Fails to Maximize Expected Utility in a Simple Foraging Model

Information theory has explained the organization of many biological phenomena, from the physiology of sensory receptive fields to the variability of certain DNA sequence ensembles. Some scholars have proposed that information should provide the central explanatory principle in biology, in the sense that any behavioral strategy that is optimal for an organism's survival must necessarily involve efficient information processing. We challenge this view by providing a counterexample. We present an analytically tractable model for a particular instance of a perception-action loop: a creature searching for a food source confined to a one-dimensional ring world. The model incorporates the statistical structure of the creature's world, the effects of the creature's actions on that structure, and the creature's strategic decision process. The model takes the form of a Markov process on an infinite dimensional state space. To analyze it we construct an exact coarse graining that reduces the model to a Markov process on a finite number of "information states". This technique allows us to make quantitative comparisons between the performance of an information-theoretically optimal strategy with other candidate strategies on a food gathering task. We find that: 1. Information optimal search does not necessarily optimize utility (expected food gain). 2. The rank ordering of search strategies by information performance does not predict their ordering by expected food obtained. 3. The relative advantage of different strategies depends on the statistical structure of the environment, in particular the variability of motion of the source. We conclude that there is no simple relationship between information and utility. Behavioral optimality does not imply information efficiency, nor is there a simple tradeoff between gaining information about a food source versus obtaining the food itself.

0 citations0 reviewsNo code linkedOpen access