Researcher profile

Matthew J Simpson

Matthew J Simpson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
11topics
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

7 published item(s)

preprint2024arXiv

Front stability of infinitely steep travelling waves in population biology

Reaction-diffusion models are often used to describe biological invasion, where populations of individuals that undergo random motility and proliferation lead to moving fronts. Many models of biological invasion are extensions of the Fisher-KPP model that describes the evolution of a 1D population density as a result of linear diffusion and logistic growth. In 2020 Fadai introduced a new model of biological invasion that was formulated as a moving boundary problem with a nonlinear degenerate diffusive flux. Fadai's model leads to travelling wave solutions with infinitely steep, well-defined fronts at the moving boundary, and the model has the mathematical advantage of being analytically tractable in certain parameter limits. We aim to provide general insight by first presenting two key extensions by considering: (i) generalised nonlinear degenerate diffusion with flux; and, (ii) solutions describing both biological invasion, and biological recession. We establish the existence of travelling wave solutions for these two extensions, and then consider stability of the travelling wave solutions by introducing a lateral perturbation of the travelling wavefront. Full 2D time-dependent level-set numerical solutions indicate that invasive travelling waves are stable to small lateral perturbations, whereas receding travelling waves are unstable. These preliminary numerical observations are corroborated through a linear stability analysis that gives more formal insight into short time growth/decay of wavefront perturbation amplitude.

preprint2022arXiv

Exact time-dependent solutions of a Fisher-KPP-like equation obtained with nonclassical symmetry analysis

We consider a family of exact solutions to a nonlinear reaction-diffusion model, constructed using nonclassical symmetry analysis. In a particular limit, the mathematical model approaches the well-known Fisher-KPP model, which means that it is related to various applications including cancer progression, wound healing and ecological invasion. The exact solution is mathematically interesting since exact solutions of the Fisher-KPP model are rare, and often restricted to long-time travelling wave solutions for special values of the travelling wave speed.

preprint2022arXiv

Geometric analysis enables biological insight from complex non-identifiable models using simple surrogates

An enduring challenge in computational biology is to balance data quality and quantity with model complexity. Tools such as identifiability analysis and information criterion have been developed to harmonise this juxtaposition, yet cannot always resolve the mismatch between available data and the granularity required in mathematical models to answer important biological questions. Often, it is only simple phenomenological models, such as the logistic and Gompertz growth models, that are identifiable from standard experimental measurements. To draw insights from the complex, non-identifiable models that incorporate key biological mechanisms of interest, we study the geometry of a map in parameter space from the complex model to a simple, identifiable, surrogate model. By studying how non-identifiable parameters in the complex model quantitatively relate to identifiable parameters in surrogate, we introduce and exploit a layer of interpretation between the set of non-identifiable parameters and the goodness-of-fit metric or likelihood studied in typical identifiability analysis. We demonstrate our approach by analysing a hierarchy of mathematical models for multicellular tumour spheroid growth. Typical data from tumour spheroid experiments are limited and noisy, and corresponding mathematical models are very often made arbitrarily complex. Our geometric approach is able to predict non-identifiabilities, subset non-identifiable parameter spaces into identifiable parameter combinations that relate to individual data features, and overall provide additional biological insight from complex non-identifiable models.

preprint2022arXiv

Non-vanishing sharp-fronted travelling wave solutions of the Fisher-Kolmogorov model

The Fisher-KPP model, and generalisations thereof, is a simple reaction-diffusion models of biological invasion that assumes individuals in the population undergo linear diffusion with diffusivity $D$, and logistic proliferation with rate $λ$. Biologically-relevant initial conditions lead to long-time travelling wave solutions that move with speed $c=2\sqrt{λD}$. Despite these attractive features, there are several biological limitations of travelling wave solutions of the Fisher-KPP model. First, these travelling wave solutions do not predict a well-defined invasion front. Second, biologically-relevant initial conditions lead to travelling waves that move with speed $c=2\sqrt{λD} > 0$. This means that, for biologically-relevant initial data, the Fisher-KPP model can not be used to study invasion with $c \ne 2\sqrt{λD}$, or retreating travelling waves with $c < 0$. Here, we reformulate the Fisher-KPP model as a moving boundary problem on $x < s(t)$, and we show that this reformulated model alleviates the key limitations of the Fisher-KPP model. Travelling wave solutions of the moving boundary problem predict a well-defined front, and can propagate with any wave speed, $-\infty < c < \infty$. Here, we establish these results using a combination of high-accuracy numerical simulations of the time-dependent partial differential equation, phase plane analysis and perturbation methods. All software required to replicate this work is available on GitHub.

preprint2022arXiv

Parameter estimation and uncertainty quantification using information geometry

In this work we: (1) review likelihood-based inference for parameter estimation and the construction of confidence regions; and, (2) explore the use of techniques from information geometry, including geodesic curves and Riemann scalar curvature, to supplement typical techniques for uncertainty quantification such as Bayesian methods, profile likelihood, asymptotic analysis and bootstrapping. These techniques from information geometry provide data-independent insights into uncertainty and identifiability, and can be used to inform data collection decisions. All code used in this work to implement the inference and information geometry techniques is available on GitHub.

preprint2020arXiv

A sharp-front moving boundary model for malignant invasion

We analyse a novel mathematical model of malignant invasion which takes the form of a two-phase moving boundary problem describing the invasion of a population of malignant cells into a population of background tissue, such as skin. Cells in both populations undergo diffusive migration and logistic proliferation. The interface between the two populations moves according to a two-phase Stefan condition. Unlike many reaction-diffusion models of malignant invasion, the moving boundary model explicitly describes the motion of the sharp front between the cancer and surrounding tissues without needing to introduce degenerate nonlinear diffusion. Numerical simulations suggest the model gives rise to very interesting travelling wave solutions that move with speed $c$, and the model supports both malignant invasion and malignant retreat, where the travelling wave can move in either the positive or negative $x$-directions. Unlike the well-studied Fisher-Kolmogorov and Porous-Fisher models where travelling waves move with a minimum wave speed $c \ge c^* > 0$, the moving boundary model leads to travelling wave solutions with $|c| < c^{**}$. We interpret these travelling wave solutions in the phase plane and show that they are associated with several features of the classical Fisher-Kolmogorov phase plane that are often disregarded as being nonphysical. We show, numerically, that the phase plane analysis compares well with long time solutions from the full partial differential equation model as well as providing accurate perturbation approximations for the shape of the travelling waves.

preprint2020arXiv

Critical length for the spreading-vanishing dichotomy in higher dimensions

We consider an extension of the classical Fisher-Kolmogorov equation, called the \textit{Fisher-Stefan} model, which is a moving boundary problem on $0 < x < L(t)$. A key property of the Fisher-Stefan model is the \textit{spreading-vanishing dichotomy}, where solutions with $L(t) > L_{\textrm{c}}$ will eventually spread as $t \to \infty$, whereas solutions where $L(t) \ngtr L_{\textrm{c}}$ will vanish as $t \to \infty$. In one dimension is it well-known that the critical length is $L_{\textrm{c}} = π/2$. In this work we re-formulate the Fisher-Stefan model in higher dimensions and calculate $L_{\textrm{c}}$ as a function of spatial dimensions in a radially symmetric coordinate system. Our results show how $L_{\textrm{c}}$ depends upon the dimension of the problem and numerical solutions of the governing partial differential equation are consistent with our calculations.