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

Matthew J. Simpson

Matthew J. Simpson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Biological Physicsmath.DSQuantitative MethodsComputationmath.APMethodologyMolecular Networksnlin.PSphysics.comp-phPopulations and Evolution

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

9 published item(s)

preprint2022arXiv

Approximate analytical solution for transient heat and mass transfer across an irregular interface

Motivated by practical applications in heat conduction and contaminant transport, we consider heat and mass diffusion across a perturbed interface separating two finite regions of distinct diffusivity. Under the assumption of continuity of the solution and diffusive flux at the interface, we use perturbation theory to develop an asymptotic expansion of the solution valid for small perturbations. Each term in the asymptotic expansion satisfies an initial-boundary value problem on the unperturbed domain subject to interface conditions depending on the previously determined terms in the asymptotic expansion. Demonstration of the perturbation solution is carried out for a specific, practically-relevant set of initial and boundary conditions with semi-analytical solutions of the initial-boundary value problems developed using standard Laplace transform and eigenfunction expansion techniques. Results for several choices of the perturbed interface confirm the perturbation solution is in good agreement with a standard numerical solution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient inference and identifiability analysis for differential equation models with random parameters

Heterogeneity is a dominant factor in the behaviour of many biological processes. Despite this, it is common for mathematical and statistical analyses to ignore biological heterogeneity as a source of variability in experimental data. Therefore, methods for exploring the identifiability of models that explicitly incorporate heterogeneity through variability in model parameters are relatively underdeveloped. We develop a new likelihood-based framework, based on moment matching, for inference and identifiability analysis of differential equation models that capture biological heterogeneity through parameters that vary according to probability distributions. As our novel method is based on an approximate likelihood function, it is highly flexible; we demonstrate identifiability analysis using both a frequentist approach based on profile likelihood, and a Bayesian approach based on Markov-chain Monte Carlo. Through three case studies, we demonstrate our method by providing a didactic guide to inference and identifiability analysis of hyperparameters that relate to the statistical moments of model parameters from independent observed data. Our approach has a computational cost comparable to analysis of models that neglect heterogeneity, a significant improvement over many existing alternatives. We demonstrate how analysis of random parameter models can aid better understanding of the sources of heterogeneity from biological data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Interpreting how nonlinear diffusion affects the fate of bistable populations using a discrete modelling framework

Understanding whether a population will survive and flourish or become extinct is a central question in population biology. One way of exploring this question is to study population dynamics using reaction-diffusion equations, where migration is usually represented as a linear diffusion term, and birth-death is represented with a bistable source term. While linear diffusion is most commonly employed to study migration, there are several limitations of this approach, such as the inability of linear diffusion-based models to predict a well-defined population front. One way to overcome this is to generalise the constant diffusivity, $D$, to a nonlinear diffusivity function $D(C)$, where $C>0$ is the density. While it has been formally established that the choice of $D(C)$ affects long-term survival or extinction of a bistable population, working solely in a classical continuum framework makes it difficult to understand precisely how the choice of $D(C)$ affects survival or extinction. Here, we address this question by working with a simple discrete simulation model that is easy to interpret. The continuum limit of the discrete model is a nonlinear reaction-diffusion equation, where the flux involves a nonlinear diffusion term and the source term is given by the strong Allee effect bistable model. We study population extinction/survival using this very intuitive discrete framework together with numerical solutions of the reaction-diffusion continuum limit equation. This approach provides clear insight into how the choice of $D(C)$ either encourages or suppresses population extinction relative to the classical linear diffusion model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes

Models of stochastic processes are widely used in almost all fields of science. Theory validation, parameter estimation, and prediction all require model calibration and statistical inference using data. However, data are almost always incomplete observations of reality. This leads to a great challenge for statistical inference because the likelihood function will be intractable for almost all partially observed stochastic processes. This renders many statistical methods, especially within a Bayesian framework, impossible to implement. Therefore, computationally expensive likelihood-free approaches are applied that replace likelihood evaluations with realisations of the model and observation process. For accurate inference, however, likelihood-free techniques may require millions of expensive stochastic simulations. To address this challenge, we develop a new method based on recent advances in multilevel and multifidelity. Our approach combines the multilevel Monte Carlo telescoping summation, applied to a sequence of approximate Bayesian posterior targets, with a multifidelity rejection sampler to minimise the number of computationally expensive exact simulations required for accurate inference. We present the derivation of our new algorithm for likelihood-free Bayesian inference, discuss practical implementation details, and demonstrate substantial performance improvements. Using examples from systems biology, we demonstrate improvements of more than two orders of magnitude over standard rejection sampling techniques. Our approach is generally applicable to accelerate other sampling schemes, such as sequential Monte Carlo, to enable feasible Bayesian analysis for realistic practical applications in physics, chemistry, biology, epidemiology, ecology and economics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Effect of Geometry on Survival and Extinction in a Moving-Boundary Problem Motivated by the Fisher-KPP Equation

The Fisher-Stefan model involves solving the Fisher-KPP equation on a domain whose boundary evolves according to a Stefan-like condition. The Fisher-Stefan model alleviates two practical limitations of the standard Fisher-KPP model when applied to biological invasion. First, unlike the Fisher-KPP equation, solutions to the Fisher-Stefan model have compact support, enabling one to define the interface between occupied and unoccupied regions unambiguously. Second, the Fisher-Stefan model admits solutions for which the population becomes extinct, which is not possible in the Fisher-KPP equation. Previous research showed that population survival or extinction in the Fisher-Stefan model depends on a critical length in one-dimensional Cartesian or radially-symmetric geometry. However, the survival and extinction behaviour for general two-dimensional regions remains unexplored. We combine analysis and level-set numerical simulations of the Fisher-Stefan model to investigate the survival-extinction conditions for rectangular-shaped initial conditions. We show that it is insufficient to generalise the critical length conditions to critical area in two-dimensions. Instead, knowledge of the region geometry is required to determine whether a population will survive or become extinct.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Biologically-informed neural networks guide mechanistic modeling from sparse experimental data

Biologically-informed neural networks (BINNs), an extension of physics-informed neural networks [1], are introduced and used to discover the underlying dynamics of biological systems from sparse experimental data. In the present work, BINNs are trained in a supervised learning framework to approximate in vitro cell biology assay experiments while respecting a generalized form of the governing reaction-diffusion partial differential equation (PDE). By allowing the diffusion and reaction terms to be multilayer perceptrons (MLPs), the nonlinear forms of these terms can be learned while simultaneously converging to the solution of the governing PDE. Further, the trained MLPs are used to guide the selection of biologically interpretable mechanistic forms of the PDE terms which provides new insights into the biological and physical mechanisms that govern the dynamics of the observed system. The method is evaluated on sparse real-world data from wound healing assays with varying initial cell densities [2].

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Diffusion in heterogeneous discs and spheres: new closed-form expressions for exit times and homogenization formulae

Mathematical models of diffusive transport underpin our understanding of chemical, biochemical and biological transport phenomena. Analysis of such models often focusses on relatively simple geometries and deals with diffusion through highly idealised homogeneous media. In contrast, practical applications of diffusive transport theory inevitably involve dealing with more complicated geometries as well as dealing with heterogeneous media. One of the most fundamental properties of diffusive transport is the concept of mean particle lifetime or mean exit time, which are particular applications of the concept of first passage time, and provide the mean time required for a diffusing particle to reach an absorbing boundary. Most formal analysis of mean particle lifetime applies to relatively simple geometries, often with homogeneous (spatially-invariant) material properties. In this work, we present a general framework that provides exact mathematical insight into the mean particle lifetime, and higher moments of particle lifetime, for point particles diffusing in heterogeneous discs and spheres with radial symmetry. Our analysis applies to geometries with an arbitrary number and arrangement of distinct layers, where transport in each layer is characterised by a distinct diffusivity. We obtain exact closed-form expressions for the mean particle lifetime for a diffusing particle released at an arbitrary location and we generalise these results to give exact, closed-form expressions for any higher-order moment of particle lifetime for a range of different boundary conditions. Finally, using these results we construct new homogenization formulae that provide an accurate simplified description of diffusion through heterogeneous discs and spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

New travelling wave solutions of the Porous-Fisher model with a moving boundary

We examine travelling wave solutions of the Porous-Fisher model, $\partial_t u(x,t)= u(x,t)\left[1-u(x,t)\right] + \partial_x \left[u(x,t) \partial_x u(x,t)\right]$, with a Stefan-like condition at the moving front, $x=L(t)$. Travelling wave solutions of this model have several novel characteristics. These travelling wave solutions: (i) move with a speed that is slower than the more standard Porous-Fisher model, $c<1/\sqrt{2}$; (ii) never lead to population extinction; (iii) have compact support and a well-defined moving front, and (iv) the travelling wave profiles have an infinite slope at the moving front. Using asymptotic analysis in two distinct parameter regimes, $c \to 0^+$ and $c \to 1/\sqrt{2}\,^-$, we obtain closed-form mathematical expressions for the travelling wave shape and speed. These approximations compare well with numerical solutions of the full problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Travelling wave solutions in a negative nonlinear diffusion-reaction model

We use a geometric approach to prove the existence of smooth travelling wave solutions of a nonlinear diffusion-reaction equation with logistic kinetics and a convex nonlinear diffusivity function which changes sign twice in our domain of interest. We determine the minimum wave speed, c*, and investigate its relation to the spectral stability of the travelling wave solutions.

0 citations0 reviewsNo code linkedOpen access