Showing 1–50 of 201 results for author: Simpson, M

Candidate Distant Trans-Neptunian Objects Detected by the New Horizons Subaru TNO Survey

Authors: Wesley C. Fraser, Simon B. Porter, Lowell Peltier, JJ Kavelaars, Anne J. Verbiscer, Marc W. Buie, S. Alan Stern, John R. Spencer, Susan D. Benecchi, Tsuyoshi Terai, Takashi Ito, Fumi Yoshida, David W. Gerdes, Kevin J. Napier, Hsing Wen Lin, Stephen D. J. Gwyn, Hayden Smotherman, Sebastien Fabbro, Kelsi N. Singer, Amanda M. Alexander, Ko Arimatsu, Maria E. Banks, Veronica J. Bray, Mohamed Ramy El-Maarry, Chelsea L. Ferrell , et al. (19 additional authors not shown)

Abstract: We report the detection of 239 trans-Neptunian Objects discovered through the on-going New Horizons survey for distant minor bodies being performed with the Hyper Suprime-Cam mosaic imager on the Subaru Telescope. These objects were discovered in images acquired with either the r2 or the recently commissioned EB-gri filter using shift and stack routines. Due to the extremely high stellar density o… ▽ More We report the detection of 239 trans-Neptunian Objects discovered through the on-going New Horizons survey for distant minor bodies being performed with the Hyper Suprime-Cam mosaic imager on the Subaru Telescope. These objects were discovered in images acquired with either the r2 or the recently commissioned EB-gri filter using shift and stack routines. Due to the extremely high stellar density of the search region down stream of the spacecraft, new machine learning techniques had to be developed to manage the extremely high false positive rate of bogus candidates produced from the shift and stack routines. We report discoveries as faint as r2$\sim26.5$. We highlight an overabundance of objects found at heliocentric distances $R\gtrsim70$~au compared to expectations from modelling of the known outer Solar System. If confirmed, these objects betray the presence of a heretofore unrecognized abundance of distant objects that can help explain a number of other observations that otherwise remain at odds with the known Kuiper Belt, including detections of serendipitous stellar occultations, and recent results from the Student Dust Counter on-board the New Horizons spacecraft. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: Accepted for publication in the Planetary Science Journal, 28 pages, 7 figures, 3 tables

Mathematical modelling and uncertainty quantification for analysis of biphasic coral reef recovery patterns

Authors: David J. Warne, Kerryn Crossman, Grace E. M. Heron, Jesse A. Sharp, Wang Jin, Paul Pao-Yen Wu, Matthew J. Simpson, Kerrie Mengersen, Juan-Carlos Ortiz

Abstract: Coral reefs are increasingly subjected to major disturbances threatening the health of marine ecosystems. Substantial research underway to develop intervention strategies that assist reefs in recovery from, and resistance to, inevitable future climate and weather extremes. To assess potential benefits of interventions, mechanistic understanding of coral reef recovery and resistance patterns is ess… ▽ More Coral reefs are increasingly subjected to major disturbances threatening the health of marine ecosystems. Substantial research underway to develop intervention strategies that assist reefs in recovery from, and resistance to, inevitable future climate and weather extremes. To assess potential benefits of interventions, mechanistic understanding of coral reef recovery and resistance patterns is essential. Recent evidence suggests that more than half of the reefs surveyed across the Great Barrier Reef (GBR) exhibit deviations from standard recovery modelling assumptions when the initial coral cover is low ($\leq 10$\%). New modelling is necessary to account for these observed patterns to better inform management strategies. We consider a new model for reef recovery at the coral cover scale that accounts for biphasic recovery patterns. The model is based on a multispecies Richards' growth model that includes a change point in the recovery patterns. Bayesian inference is applied for uncertainty quantification of key parameters for assessing reef health and recovery patterns. This analysis is applied to benthic survey data from the Australian Institute of Marine Sciences (AIMS). We demonstrate agreement between model predictions and data across every recorded recovery trajectory with at least two years of observations following disturbance events occurring between 1992--2020. This new approach will enable new insights into the biological, ecological and environmental factors that contribute to the duration and severity of biphasic coral recovery patterns across the GBR. These new insights will help to inform managements and monitoring practice to mitigate the impacts of climate change on coral reefs. △ Less

Submitted 27 June, 2024; originally announced June 2024.

MSC Class: 62P12 (Primary)

Mechanical cell interactions on curved interfaces

Authors: Pascal R. Buenzli, Shahak Kuba, Ryan J. Murphy, Matthew J. Simpson

Abstract: We propose a simple mathematical model to describe the mechanical relaxation of cells within a curved epithelial tissue layer represented by an arbitrary curve in two-dimensional space. The model represents the mechanics of the cell body either by straight springs between points of the curve, or by curved springs whose shape follows the curve. To understand the collective behaviour of these discre… ▽ More We propose a simple mathematical model to describe the mechanical relaxation of cells within a curved epithelial tissue layer represented by an arbitrary curve in two-dimensional space. The model represents the mechanics of the cell body either by straight springs between points of the curve, or by curved springs whose shape follows the curve. To understand the collective behaviour of these discrete models of cells at the broader tissue scale, we devise an appropriate continuum limit in which the number of cells is constant but the number of springs tends to infinity. The continuum limit shows that (i)~the straight spring model and the curved spring model converge to the same dynamics; and (ii)~the density of cells becomes governed by a diffusion equation in arc length space with second-order accuracy, where diffusion may be linear or nonlinear depending on the choice of the spring restoring force law. Our derivation of the continuum limit justifies that to reach consistent dynamics as the number of springs increases, the spring restoring force laws must be rescaled appropriately. Despite mechanical relaxation occurring within a curved tissue layer, we find that the curvature of the tissue does not affect tangential stress nor the mechanics-induced redistribution of cells within the layer in the continuum limit. However, the cell's normal stress does depend on curvature due to surface tension induced by the tangential forces. By characterising the full stress state of a cell, these models provide a basis to represent further mechanobiological processes. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 34 pages, 13 figures; v2: minor edits

Likelihood-based inference, identifiability and prediction using count data from lattice-based random walk models

Authors: Yihan Liu, David J Warne, Matthew J Simpson

Abstract: In vitro cell biology experiments are routinely used to characterize cell migration properties under various experimental conditions. These experiments can be interpreted using lattice-based random walk models to provide insight into underlying biological mechanisms, and continuum limit partial differential equation (PDE) descriptions of the stochastic models can be used to efficiently explore mod… ▽ More In vitro cell biology experiments are routinely used to characterize cell migration properties under various experimental conditions. These experiments can be interpreted using lattice-based random walk models to provide insight into underlying biological mechanisms, and continuum limit partial differential equation (PDE) descriptions of the stochastic models can be used to efficiently explore model properties instead of relying on repeated stochastic simulations. Working with efficient PDE models is of high interest for parameter estimation algorithms that typically require a large number of forward model simulations. Quantitative data from cell biology experiments usually involves non-negative cell counts in different regions of the experimental images, and it is not obvious how to relate finite, noisy count data to the solutions of continuous PDE models that correspond to noise-free density profiles. In this work we illustrate how to develop and implement likelihood-based methods for parameter estimation, parameter identifiability and model prediction for lattice-based models describing collective migration with an arbitrary number of interacting subpopulations. We implement a standard additive Gaussian measurement error model as well as a new physically-motivated multinomial measurement error model that relates noisy count data with the solution of continuous PDE models. Both measurement error models lead to similar outcomes for parameter estimation and parameter identifiability, whereas the standard additive Gaussian measurement error model leads to non-physical prediction outcomes. In contrast, the new multinomial measurement error model involves a lower computational overhead for parameter estimation and identifiability analysis, as well as leading to physically meaningful model predictions. △ Less

Submitted 23 June, 2024; originally announced June 2024.

Comments: 34 pages, 7 figures

MSC Class: 92B99; 82M99; 62M99

History-enhanced ICT For Sustainability education: Learning together with Business Computing students

Authors: Ian Brooks, Laura Harrison, Mark Reeves, Martin Simpson, Rose Wallis

Abstract: This research explores the use of History to enhance education in the field of ICT For Sustainability ICT4S in response to a challenge from the ICT4S 2023 conference. No previous studies were found in ICT4S but the literature on History and Education for Sustainable Development is reviewed. An ICT4S lecturer collaborated with History lecturers to add an historic parallel to each weeks teaching on… ▽ More This research explores the use of History to enhance education in the field of ICT For Sustainability ICT4S in response to a challenge from the ICT4S 2023 conference. No previous studies were found in ICT4S but the literature on History and Education for Sustainable Development is reviewed. An ICT4S lecturer collaborated with History lecturers to add an historic parallel to each weeks teaching on a Sustainable Business and Computing unit for final year undergraduate BSc Business Computing students. A list of the topics and rationale is provided. Student perceptions were surveyed before and after the teaching and semi-structured interviews carried out. A majority of students saw relevance to their degree and career. There was an increase in the proportion of students with interest in History. The paper explores the lessons learned from the interdisciplinary collaboration, including topic choice, format and perceived value. The project has enhanced the way we approach our subjects as computing and history educators. We believe this is the first empirical, survey-based study of the use of history to enhance ICT4S education. The team will extend the research to a larger unit covering a wider range of computing degrees. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Parameter identifiability, parameter estimation and model prediction for differential equation models

Authors: Matthew J Simpson, Ruth E Baker

Abstract: Interpreting data with mathematical models is an important aspect of real-world applied mathematical modeling. Very often we are interested to understand the extent to which a particular data set informs and constrains model parameters. This question is closely related to the concept of parameter identifiability, and in this article we present a series of computational exercises to introduce tools… ▽ More Interpreting data with mathematical models is an important aspect of real-world applied mathematical modeling. Very often we are interested to understand the extent to which a particular data set informs and constrains model parameters. This question is closely related to the concept of parameter identifiability, and in this article we present a series of computational exercises to introduce tools that can be used to assess parameter identifiability, estimate parameters and generate model predictions. Taking a likelihood-based approach, we show that very similar ideas and algorithms can be used to deal with a range of different mathematical modelling frameworks. The exercises and results presented in this article are supported by a suite of open access codes that can be accessed on GitHub. △ Less

Submitted 28 May, 2024; v1 submitted 13 May, 2024; originally announced May 2024.

Comments: 22 pages, 6 Figures

MSC Class: 97M10; 97M50; 00A71

A generalised sigmoid population growth model with energy dependence: application to quantify the tipping point for Antarctic shallow seabed algae

Authors: Elise Mills, Graeme F. Clark, Matthew J. Simpson, Mark Baird, Matthew P. Adams

Abstract: Sigmoid growth models are often used to study population dynamics. The size of a population at equilibrium commonly depends explicitly on the availability of resources, such as an energy or nutrient source, which is not explicit in standard sigmoid growth models. A simple generalised extension of sigmoid growth models is introduced that can explicitly account for this resource-dependence, demonstr… ▽ More Sigmoid growth models are often used to study population dynamics. The size of a population at equilibrium commonly depends explicitly on the availability of resources, such as an energy or nutrient source, which is not explicit in standard sigmoid growth models. A simple generalised extension of sigmoid growth models is introduced that can explicitly account for this resource-dependence, demonstrated by three examples of this family of models of increasing mathematical complexity. Each model is calibrated and compared to observed data for algae under sea-ice in Antarctic coastal waters. It was found that through careful construction, models satisfying the proposed framework can estimate key properties of a sea-ice break-out controlled tipping point for the algae, which cannot be estimated using standard sigmoid growth models. The proposed broader family of energy-dependent sigmoid growth models likely has usage in many population growth contexts where resources limit population size. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Structured methods for parameter inference and uncertainty quantification for mechanistic models in the life sciences

Authors: Michael J. Plank, Matthew J. Simpson

Abstract: Parameter inference and uncertainty quantification are important steps when relating mathematical models to real-world observations, and when estimating uncertainty in model predictions. However, methods for doing this can be computationally expensive, particularly when the number of unknown model parameters is large. The aim of this study is to develop and test an efficient profile likelihood-bas… ▽ More Parameter inference and uncertainty quantification are important steps when relating mathematical models to real-world observations, and when estimating uncertainty in model predictions. However, methods for doing this can be computationally expensive, particularly when the number of unknown model parameters is large. The aim of this study is to develop and test an efficient profile likelihood-based method, which takes advantage of the structure of the mathematical model being used. We do this by identifying specific parameters that affect model output in a known way, such as a linear scaling. We illustrate the method by applying it to three caricature models from different areas of the life sciences: (i) a predator-prey model from ecology; (ii) a compartment-based epidemic model from health sciences; and, (iii) an advection-diffusion-reaction model describing transport of dissolved solutes from environmental science. We show that the new method produces results of comparable accuracy to existing profile likelihood methods, but with substantially fewer evaluations of the forward model. We conclude that our method could provide a much more efficient approach to parameter inference for models where a structured approach is feasible. Code to apply the new method to user-supplied models and data is provided via a publicly accessible repository. △ Less

Submitted 8 July, 2024; v1 submitted 3 March, 2024; originally announced March 2024.

Fisher-KPP-type models of biological invasion: Open source computational tools, key concepts and analysis

Authors: Matthew J Simpson, Scott W McCue

Abstract: This review provides open-access computational tools that support a range of mathematical approaches to analyse three related scalar reaction-diffusion models used to study biological invasion. Starting with the classic Fisher-Kolmogorov (Fisher-KPP) model, we illustrate how computational methods can be used to explore time-dependent partial differential equation (PDE) solutions in parallel with p… ▽ More This review provides open-access computational tools that support a range of mathematical approaches to analyse three related scalar reaction-diffusion models used to study biological invasion. Starting with the classic Fisher-Kolmogorov (Fisher-KPP) model, we illustrate how computational methods can be used to explore time-dependent partial differential equation (PDE) solutions in parallel with phase plane and regular perturbation techniques to explore invading travelling wave solutions moving with dimensionless speed $c \ge 2$. To overcome the lack of a well-defined sharp front in solutions of the Fisher-KPP model, we also review two alternative modeling approaches. The first is the Porous-Fisher model where the linear diffusion term is replaced with a degenerate nonlinear diffusion term. Using phase plane and regular perturbation methods, we explore the distinction between sharp- and smooth-fronted invading travelling waves that move with dimensionless speed $c \ge 1/\sqrt{2}$. The second alternative approach is to reformulate the Fisher-KPP model as a moving boundary problem on $0 < x < L(t)$, leading to the Fisher-Stefan model with sharp-fronted travelling wave solutions arising from a PDE model with a linear diffusion term. Time-dependent PDE solutions and phase plane methods show that travelling wave solutions of the Fisher-Stefan model can describe both biological invasion $(c > 0)$ and biological recession $(c < 0)$. Open source Julia code to replicate all computational results in this review is available on GitHub; we encourage researchers to use this code directly or to adapt the code as required for more complicated models. △ Less

Submitted 24 April, 2024; v1 submitted 3 March, 2024; originally announced March 2024.

Comments: 55 Pages, 14 Figures

MSC Class: 35Q92

Front stability of infinitely steep travelling waves in population biology

Authors: Matthew J Simpson, Nizhum Rahman, Alexander KY Tam

Abstract: 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… ▽ More 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. △ Less

Submitted 5 January, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 8 Figures, 25 pages

MSC Class: 35C07; 35R37; 92D25

Discrete and continuous mathematical models of sharp-fronted collective cell migration and invasion

Authors: Matthew J Simpson, Keeley M Murphy, Scott W McCue, Pascal R Buenzli

Abstract: Mathematical models describing the spatial spreading and invasion of populations of biological cells are often developed in a continuum modelling framework using reaction-diffusion equations. While continuum models based on linear diffusion are routinely employed and known to capture key experimental observations, linear diffusion fails to predict well-defined sharp fronts that are often observed… ▽ More Mathematical models describing the spatial spreading and invasion of populations of biological cells are often developed in a continuum modelling framework using reaction-diffusion equations. While continuum models based on linear diffusion are routinely employed and known to capture key experimental observations, linear diffusion fails to predict well-defined sharp fronts that are often observed experimentally. This observation has motivated the use of nonlinear degenerate diffusion, however these nonlinear models and the associated parameters lack a clear biological motivation and interpretation. Here we take a different approach by developing a stochastic discrete lattice-based model incorporating biologically-inspired mechanisms and then deriving the reaction-diffusion continuum limit. Inspired by experimental observations, agents in the simulation deposit extracellular material, that we call a substrate, locally onto the lattice, and the motility of agents is taken to be proportional to the substrate density. Discrete simulations that mimic a two--dimensional circular barrier assay illustrate how the discrete model supports both smooth and sharp-fronted density profiles depending on the rate of substrate deposition. Coarse-graining the discrete model leads to a novel partial differential equation (PDE) model whose solution accurately approximates averaged data from the discrete model. The new discrete model and PDE approximation provides a simple, biologically motivated framework for modelling the spreading, growth and invasion of cell populations with well-defined sharp fronts △ Less

Submitted 19 January, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 47 Pages, 8 Figures

MSC Class: 92-10

The DECam Ecliptic Exploration Project (DEEP): V. The Absolute Magnitude Distribution of the Cold Classical Kuiper Belt

Authors: Kevin J. Napier, Hsing-Wen Lin, David W. Gerdes, Fred C. Adams, Anna M. Simpson, Matthew W. Porter, Katherine G. Weber, Larissa Markwardt, Gabriel Gowman, Hayden Smotherman, Pedro H. Bernardinelli, Mario Jurić, Andrew J. Connolly, J. Bryce Kalmbach, Stephen K. N. Portillo, David E. Trilling, Ryder Strauss, William J. Oldroyd, Chadwick A. Trujillo, Colin Orion Chandler, Matthew J. Holman, Hilke E. Schlichting, Andrew McNeill, the DEEP Collaboration

Abstract: The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth,… ▽ More The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth, enabling quantitative leaps forward in our understanding of the Kuiper Belt populations. This work reports results from an analysis of twenty 3 sq.\ deg.\ DECam fields along the invariable plane. We characterize the efficiency and false-positive rates for our moving-object detection pipeline, and use this information to construct a Bayesian signal probability for each detected source. This procedure allows us to treat all of our Kuiper Belt Object (KBO) detections statistically, simultaneously accounting for efficiency and false positives. We detect approximately 2300 candidate sources with KBO-like motion at S/N $>6.5$. We use a subset of these objects to compute the luminosity function of the Kuiper Belt as a whole, as well as the Cold Classical (CC) population. We also investigate the absolute magnitude ($H$) distribution of the CCs, and find consistency with both an exponentially tapered power-law, which is predicted by streaming instability models of planetesimal formation, and a rolling power law. Finally, we provide an updated mass estimate for the Cold Classical Kuiper Belt of $M_{CC}(H_r < 12) = 0.0017^{+0.0010}_{-0.0004} M_{\oplus}$, assuming albedo $p = 0.15$ and density $ρ= 1$ g cm$^{-3}$. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: Accepted by PSJ

Pushing coarse-grained models beyond the continuum limit using equation learning

Authors: Daniel J. VandenHeuvel, Pascal R. Buenzli, Matthew J. Simpson

Abstract: Mathematical modelling of biological population dynamics often involves proposing high fidelity discrete agent-based models that capture stochasticity and individual-level processes. These models are often considered in conjunction with an approximate coarse-grained differential equation that captures population-level features only. These coarse-grained models are only accurate in certain asymptot… ▽ More Mathematical modelling of biological population dynamics often involves proposing high fidelity discrete agent-based models that capture stochasticity and individual-level processes. These models are often considered in conjunction with an approximate coarse-grained differential equation that captures population-level features only. These coarse-grained models are only accurate in certain asymptotic parameter regimes, such as enforcing that the time scale of individual motility far exceeds the time scale of birth/death processes. When these coarse-grained models are accurate, the discrete model still abides by conservation laws at the microscopic level, which implies that there is some macroscopic conservation law that can describe the macroscopic dynamics. In this work, we introduce an equation learning framework to find accurate coarse-grained models when standard continuum limit approaches are inaccurate. We demonstrate our approach using a discrete mechanical model of epithelial tissues, considering a series of four case studies that consider problems with and without free boundaries, and with and without proliferation, illustrating how we can learn macroscopic equations describing mechanical relaxation, cell proliferation, and the equation governing the dynamics of the free boundary of the tissue. While our presentation focuses on this biological application, our approach is more broadly applicable across a range of scenarios where discrete models are approximated by approximate continuum-limit descriptions. All code and data to reproduce this work are available at https://github.com/DanielVandH/StepwiseEQL.jl. △ Less

Submitted 18 December, 2023; v1 submitted 21 August, 2023; originally announced August 2023.

Comments: Fixed some incorrect references to a figure

MSC Class: 35R30 (Primary); 92D25 (Secondary)

Implementing measurement error models with mechanistic mathematical models in a likelihood-based framework for estimation, identifiability analysis, and prediction in the life sciences

Authors: Ryan J. Murphy, Oliver J. Maclaren, Matthew J. Simpson

Abstract: Throughout the life sciences we routinely seek to interpret measurements and observations using parameterised mechanistic mathematical models. A fundamental and often overlooked choice in this approach involves relating the solution of a mathematical model with noisy and incomplete measurement data. This is often achieved by assuming that the data are noisy measurements of the solution of a determ… ▽ More Throughout the life sciences we routinely seek to interpret measurements and observations using parameterised mechanistic mathematical models. A fundamental and often overlooked choice in this approach involves relating the solution of a mathematical model with noisy and incomplete measurement data. This is often achieved by assuming that the data are noisy measurements of the solution of a deterministic mathematical model, and that measurement errors are additive and normally distributed. While this assumption of additive Gaussian noise is extremely common and simple to implement and interpret, it is often unjustified and can lead to poor parameter estimates and non-physical predictions. One way to overcome this challenge is to implement a different measurement error model. In this review, we demonstrate how to implement a range of measurement error models in a likelihood-based framework for estimation, identifiability analysis, and prediction, called Profile-Wise Analysis. This frequentist approach to uncertainty quantification for mechanistic models leverages the profile likelihood for targeting parameters and understanding their influence on predictions. Case studies, motivated by simple caricature models routinely used in systems biology and mathematical biology literature, illustrate how the same ideas apply to different types of mathematical models. Open-source Julia code to reproduce results is available on GitHub. △ Less

Submitted 8 November, 2023; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: Main: 43 pages, 12 figures. Supplementary: 24 pages, 12 figures

MSC Class: 62F99; 92B99; 92C42

Survival, extinction, and interface stability in a two--phase moving boundary model of biological invasion

Authors: Matthew J Simpson, Nizhum Rahman, Scott W McCue, Alexander KY Tam

Abstract: We consider a moving boundary mathematical model of biological invasion. The model describes the spatiotemporal evolution of two adjacent populations: each population undergoes linear diffusion and logistic growth, and the boundary between the two populations evolves according to a two--phase Stefan condition. This mathematical model describes situations where one population invades into regions o… ▽ More We consider a moving boundary mathematical model of biological invasion. The model describes the spatiotemporal evolution of two adjacent populations: each population undergoes linear diffusion and logistic growth, and the boundary between the two populations evolves according to a two--phase Stefan condition. This mathematical model describes situations where one population invades into regions occupied by the other population, such as the spreading of a malignant tumour into surrounding tissues. Full time--dependent numerical solutions are obtained using a level--set numerical method. We use these numerical solutions to explore several properties of the model including: (i) survival and extinction of one population initially surrounded by the other; and (ii) linear stability of the moving front boundary in the context of a travelling wave solution subjected to transverse perturbations. Overall, we show that many features of the well--studied one--phase single population analogue of this model can be very different in the more realistic two--phase setting. These results are important because realistic examples of biological invasion involve interactions between multiple populations and so great care should be taken when extrapolating predictions from a one--phase single population model to cases for which multiple populations are present. Open source Julia--based software is available on GitHub to replicate all results in this study. △ Less

Submitted 1 September, 2023; v1 submitted 27 June, 2023; originally announced June 2023.

Comments: 36 pages. 10 figures

MSC Class: 35C07; 35R37; 92D25

Exact sharp-fronted solutions for nonlinear diffusion on evolving domains

Authors: Stuart T. Johnston, Matthew J. Simpson

Abstract: Models of diffusive processes that occur on evolving domains are frequently employed to describe biological and physical phenomena, such as diffusion within expanding tissues or substrates. Previous investigations into these models either report numerical solutions or require an assumption of linear diffusion to determine exact solutions. Unfortunately, numerical solutions do not reveal the relati… ▽ More Models of diffusive processes that occur on evolving domains are frequently employed to describe biological and physical phenomena, such as diffusion within expanding tissues or substrates. Previous investigations into these models either report numerical solutions or require an assumption of linear diffusion to determine exact solutions. Unfortunately, numerical solutions do not reveal the relationship between the model parameters and the solution features. Additionally, experimental observations typically report the presence of sharp fronts, which are not captured by linear diffusion. Here we address both limitations by presenting exact sharp-fronted solutions to a model of degenerate nonlinear diffusion on a growing domain. We obtain the solution by identifying a series of transformations that converts the model of a nonlinear diffusive process on an evolving domain to a nonlinear diffusion equation on a fixed domain, which admits known exact solutions for certain choices of diffusivity functions. We determine expressions for critical time scales and domain growth rates such that the diffusive population never reaches the domain boundaries and hence the solution remains valid. △ Less

Submitted 5 October, 2023; v1 submitted 12 June, 2023; originally announced June 2023.

Generalised likelihood profiles for models with intractable likelihoods

Authors: David J. Warne, Oliver J. Maclaren, Elliot J. Carr, Matthew J. Simpson, Christopher Drovandi

Abstract: Likelihood profiling is an efficient and powerful frequentist approach for parameter estimation, uncertainty quantification and practical identifiablity analysis. Unfortunately, these methods cannot be easily applied for stochastic models without a tractable likelihood function. Such models are typical in many fields of science, rendering these classical approaches impractical in these settings. T… ▽ More Likelihood profiling is an efficient and powerful frequentist approach for parameter estimation, uncertainty quantification and practical identifiablity analysis. Unfortunately, these methods cannot be easily applied for stochastic models without a tractable likelihood function. Such models are typical in many fields of science, rendering these classical approaches impractical in these settings. To address this limitation, we develop a new approach to generalising the methods of likelihood profiling for situations when the likelihood cannot be evaluated but stochastic simulations of the assumed data generating process are possible. Our approach is based upon recasting developments from generalised Bayesian inference into a frequentist setting. We derive a method for constructing generalised likelihood profiles and calibrating these profiles to achieve desired frequentist coverage for a given coverage level. We demonstrate the performance of our method on realistic examples from the literature and highlight the capability of our approach for the purpose of practical identifability analysis for models with intractable likelihoods. △ Less

Submitted 19 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

MSC Class: 62M20 (Primary) 62-08; 62F40 (Secondary)

Simplified models of diffusion in radially-symmetric geometries

Authors: Luke P. Filippini, Matthew J. Simpson, Elliot J. Carr

Abstract: We consider diffusion-controlled release of particles from $d$-dimensional radially-symmetric geometries. A quantity commonly used to characterise such diffusive processes is the proportion of particles remaining within the geometry over time, denoted as $P(t)$. The stochastic approach for computing $P(t)$ is time-consuming and lacks analytical insight into key parameters while the continuum appro… ▽ More We consider diffusion-controlled release of particles from $d$-dimensional radially-symmetric geometries. A quantity commonly used to characterise such diffusive processes is the proportion of particles remaining within the geometry over time, denoted as $P(t)$. The stochastic approach for computing $P(t)$ is time-consuming and lacks analytical insight into key parameters while the continuum approach yields complicated expressions for $P(t)$ that obscure the influence of key parameters and complicate the process of fitting experimental release data. In this work, to address these issues, we develop several simple surrogate models to approximate $P(t)$ by matching moments with the continuum analogue of the stochastic diffusion model. Surrogate models are developed for homogeneous slab, circular, annular, spherical and spherical shell geometries with a constant particle movement probability and heterogeneous slab, circular, annular and spherical geometries, comprised of two concentric layers with different particle movement probabilities. Each model is easy to evaluate, agrees well with both stochastic and continuum calculations of $P(t)$ and provides analytical insight into the key parameters of the diffusive transport system: dimension, diffusivity, geometry and boundary conditions. △ Less

Submitted 18 September, 2023; v1 submitted 28 April, 2023; originally announced April 2023.

Comments: 22 pages, 3 figures, accepted version

Journal ref: Physica A: Statistical Mechanics and its Applications, 626 (2023) 129067

Exact solutions for diffusive transport on heterogeneous growing domains

Authors: Stuart T. Johnston, Matthew J. Simpson

Abstract: From the smallest biological systems to the largest cosmological structures, spatial domains undergo expansion and contraction. Within these growing domains, diffusive transport is a common phenomenon. Mathematical models have been widely employed to investigate diffusive processes on growing domains. However, a standard assumption is that the domain growth is spatially uniform. There are many rel… ▽ More From the smallest biological systems to the largest cosmological structures, spatial domains undergo expansion and contraction. Within these growing domains, diffusive transport is a common phenomenon. Mathematical models have been widely employed to investigate diffusive processes on growing domains. However, a standard assumption is that the domain growth is spatially uniform. There are many relevant examples where this is not the case, such as the colonisation of growing gut tissue by neural crest cells. As such, it is not straightforward to disentangle the individual roles of heterogeneous growth and diffusive transport. Here we present exact solutions to models of diffusive transport on domains undergoing spatially non-uniform growth. The exact solutions are obtained via a combination of transformation, convolution and superposition techniques. We verify the accuracy of these solutions via comparison with simulations of a corresponding lattice-based random walk. We explore various domain growth functions, including linear growth, exponential growth and contraction, and oscillatory growth. Provided the domain size remains positive, we find that the derived solutions are valid. The exact solutions reveal the relationship between model parameters, such as the diffusivity and the type and rate of domain growth, and key statistics, such as the survival and splitting probabilities. △ Less

Submitted 26 June, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Tunneling in a very slow ion-molecule reaction

Authors: Robert Wild, Markus Nötzold, Malcolm Simpson, Thuy Dung Tran, Roland Wester

Abstract: Quantum tunneling reactions play a significant role in chemistry when classical pathways are energetically forbidden, be it in gas phase reactions, surface diffusion, or liquid phase chemistry. In general, such tunneling reactions are challenging to calculate theoretically, given the high dimensionality of the quantum dynamics, and also very difficult to identify experimentally. Hydrogenic systems… ▽ More Quantum tunneling reactions play a significant role in chemistry when classical pathways are energetically forbidden, be it in gas phase reactions, surface diffusion, or liquid phase chemistry. In general, such tunneling reactions are challenging to calculate theoretically, given the high dimensionality of the quantum dynamics, and also very difficult to identify experimentally. Hydrogenic systems, however, allow for accurate first-principles calculations. In this way the rate of the gas phase proton transfer tunneling reaction of hydrogen molecules with deuterium anions, H_2 + D^- --> H^- + HD, has been calculated, but has so far lacked experimental verification. Here we present high-sensitivity measurements of the reaction rate carried out in a cryogenic 22-pole ion trap. We observe an extremely low rate constant of (5.2 +- 1.6) x 10^(-20) cm^3/s. This measured value agrees with quantum tunneling calculations, serving as a benchmark for molecular theory and advancing the understanding of fundamental collision processes. A deviation of the reaction rate from linear scaling, which is observed at high H_2 densities, can be traced back to previously unobserved heating dynamics in radiofrequency ion traps. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Journal ref: Nature 615, 425 (2023)

An ever-present $Gaia$ snail shell triggered by a dark matter wake

Authors: Robert J. J. Grand, Rüdiger Pakmor, Francesca Fragkoudi, Facundo A. Gómez, Wilma Trick, Christine M. Simpson, Freeke van de Voort, Rebekka Bieri

Abstract: We utilize a novel numerical technique to model star formation in cosmological simulations of galaxy formation - called Superstars - to simulate a Milky Way-like galaxy with $\gtrsim10^8$ star particles to study the formation and evolution of out-of-equilibrium stellar disc structures in a full cosmological setting. In the plane defined by the coordinate and velocity perpendicular to the mid-plane… ▽ More We utilize a novel numerical technique to model star formation in cosmological simulations of galaxy formation - called Superstars - to simulate a Milky Way-like galaxy with $\gtrsim10^8$ star particles to study the formation and evolution of out-of-equilibrium stellar disc structures in a full cosmological setting. In the plane defined by the coordinate and velocity perpendicular to the mid-plane (vertical phase space, $\{Z,V_Z\}$), stars in Solar-like volumes at late times exhibit clear spirals qualitatively similar in shape and amplitude to the $Gaia$ ``Snail shell'' phase spiral. We show that the phase spiral forms at a look back time of $\sim 6$ Gyr during the pericentric passage of a $\sim10^{10}$ $\rm M_{\odot}$ satellite on a polar orbit. This satellite stimulates the formation of a resonant wake in the dark matter halo while losing mass at a rate of $\sim0.5$-$1$ dex per orbit loop. The peak magnitude of the wake-induced gravitational torque at the Solar radius is $\sim 8$ times that from the satellite, and triggers the formation of a disc warp that wraps up into a vertical phase spiral over time. As the wake decays, the phase spiral propagates several Gigayears to present-day and can be described as ``ever-present'' once stable disc evolution is established. These results suggest an alternative scenario to explain the $Gaia$ phase spiral which does not rely on a perturbation from bar buckling or a recent direct hit from a satellite. △ Less

Submitted 28 June, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: accepted for pub in mnras. 11 + 5 pages. For a high-cadence animated version of figure 1, see https://wwwmpa.mpa-garching.mpg.de/auriga/movies/multi_halo_6_sf64.mp4

Colour and infall time distributions of satellite galaxies in simulated Milky-Way analogs

Authors: Yue Pan, Christine M. Simpson, Andrey Kravtsov, Facundo A. Gómez, Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Viraj Manwadkar, Clarke J. Esmerian

Abstract: We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 M_\odot \lesssim M_\star \lesssim 10^{11} M_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulati… ▽ More We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 M_\odot \lesssim M_\star \lesssim 10^{11} M_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulations, which was originally designed to model Milky Way-like host galaxies, resemble the populations in the Exploration of Local VolumE Satellites (ELVES) Survey and the Satellites Around Galactic Analogs (SAGA) survey in their luminosity function in the luminosity range $-12 \lesssim M_V \lesssim -15$ and resemble ELVES in their quenched fraction and colour--magnitude distribution in the luminosity range $-12 \lesssim M_g \lesssim -15$. We find that satellites transition from blue colours to red colours at the luminosity range $-15 \lesssim M_g \lesssim -12$ in both the simulations and observations and we show that this shift is driven by environmental effects in the simulations. We demonstrate also that the colour distribution in both simulations and observations can be decomposed into two statistically distinct populations based on their morphological type or star-forming status that are statistically distinct. In the simulations, these two populations also have statistically distinct infall time distributions. The comparison presented here seems to indicate that the tension between the quenched fraction in SAGA and simulations is resolved by the improved target selection of ELVES, but there are still tensions in understanding the colours of faint galaxies, of which ELVES appears to have a significant population of faint blue satellites not recovered in Auriga. △ Less

Submitted 30 January, 2023; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 15 pages, 10 figures. Published in MNRAS on December 15, 2022

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

Authors: Alexander P Browning, Matthew J Simpson

Abstract: 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 onl… ▽ More 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. △ Less

Submitted 3 August, 2022; originally announced August 2022.

MSC Class: 92B99; 62F99

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

Authors: Alexander P. Browning, Christopher Drovandi, Ian W. Turner, Adrianne L. Jenner, Matthew J. Simpson

Abstract: 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 under… ▽ More 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. △ Less

Submitted 27 October, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: Minor changes to text. Additional results in supplementary material. Additional statistics regarding results given in main and supplementary material

MSC Class: 92B99; 62F99

Machine Learning for Galactic Archaeology: A chemistry-based neural network method for identification of accreted disc stars

Authors: Thorold Tronrud, Patricia B. Tissera, Facundo A. Gómez, Robert J. J. Grand, Ruediger Pakmor, Federico Marinacci, Christine M. Simpson

Abstract: We develop a method ('Galactic Archaeology Neural Network', GANN) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy's own local environment defined by the stellar halo and the surviving satellites. We demonstrate that th… ▽ More We develop a method ('Galactic Archaeology Neural Network', GANN) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy's own local environment defined by the stellar halo and the surviving satellites. We demonstrate that this approach allows the detection of accreted stars that are spatially mixed into the disc. Two performance measures are defined - recovery fraction of accreted stars, and the probability that a star with a positive (accreted) classification is a true-positive result, P(TP). As the NNM output is akin to an assigned probability, we are able to determine positivity based on flexible threshold values that can be adjusted easily to refine the selection of presumed-accreted stars. We find that GANN identifies accreted disc stars within simulated galaxies, with high recovery fraction and/or high P(TP). We also find that stars in Gaia-Enceladus-Sausage (GES) mass systems are over 50% recovered by our NNMs in the majority (18/24) of cases. Additionally, nearly every individual source of accreted stars is detected at 10% or more of its peak stellar mass in the disc. We also demonstrate that a conglomerated NNM, trained on the halo and satellite stars from all of the Auriga galaxies provides the most consistent results, and could prove to be an intriguing future approach as our observational capabilities expand. △ Less

Submitted 19 July, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 19 pages, 12 figures

Pattern formation and front stability for a moving-boundary model of biological invasion and recession

Authors: Alexander K. Y. Tam, Matthew J. Simpson

Abstract: We investigate pattern formation in a two-dimensional (2D) Fisher--Stefan model, which involves solving the Fisher--KPP equation on a compactly-supported region with a moving boundary. By combining the Fisher--KPP and classical Stefan theory, the Fisher--Stefan model alleviates two limitations of the Fisher--KPP equation for biological populations. In this work, we investigate whether the 2D Fishe… ▽ More We investigate pattern formation in a two-dimensional (2D) Fisher--Stefan model, which involves solving the Fisher--KPP equation on a compactly-supported region with a moving boundary. By combining the Fisher--KPP and classical Stefan theory, the Fisher--Stefan model alleviates two limitations of the Fisher--KPP equation for biological populations. In this work, we investigate whether the 2D Fisher--Stefan model predicts pattern formation, by analysing the linear stability of planar travelling wave solutions to sinusoidal transverse perturbations. Planar fronts of the Fisher--KPP equation are linearly stable. Similarly, we demonstrate that invading planar fronts ($c > 0$) of the Fisher--Stefan model are linearly stable to perturbations of all wave numbers. However, our analysis demonstrates that receding planar fronts ($c < 0$) of the Fisher--Stefan model are linearly unstable for all wave numbers. This is analogous to unstable solutions for planar solidification in the classical Stefan problem. Introducing a surface tension regularisation stabilises receding fronts for short-wavelength perturbations, giving rise to a range of unstable modes and a most unstable wave number. We supplement linear stability analysis with level-set numerical solutions that corroborate theoretical results. Overall, front instability in the Fisher--Stefan model suggests a new mechanism for pattern formation in receding biological populations. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: 35 pages, 8 figures

MSC Class: 35C07; 35R37; 92D25

How is model-related uncertainty quantified and reported in different disciplines?

Authors: Emily G. Simmonds, Kwaku Peprah Adjei, Christoffer Wold Andersen, Janne Cathrin Hetle Aspheim, Claudia Battistin, Nicola Bulso, Hannah Christensen, Benjamin Cretois, Ryan Cubero, Ivan A. Davidovich, Lisa Dickel, Benjamin Dunn, Etienne Dunn-Sigouin, Karin Dyrstad, Sigurd Einum, Donata Giglio, Haakon Gjerlow, Amelie Godefroidt, Ricardo Gonzalez-Gil, Soledad Gonzalo Cogno, Fabian Grosse, Paul Halloran, Mari F. Jensen, John James Kennedy, Peter Egge Langsaether , et al. (18 additional authors not shown)

Abstract: How do we know how much we know? Quantifying uncertainty associated with our modelling work is the only way we can answer how much we know about any phenomenon. With quantitative science now highly influential in the public sphere and the results from models translating into action, we must support our conclusions with sufficient rigour to produce useful, reproducible results. Incomplete considera… ▽ More How do we know how much we know? Quantifying uncertainty associated with our modelling work is the only way we can answer how much we know about any phenomenon. With quantitative science now highly influential in the public sphere and the results from models translating into action, we must support our conclusions with sufficient rigour to produce useful, reproducible results. Incomplete consideration of model-based uncertainties can lead to false conclusions with real world impacts. Despite these potentially damaging consequences, uncertainty consideration is incomplete both within and across scientific fields. We take a unique interdisciplinary approach and conduct a systematic audit of model-related uncertainty quantification from seven scientific fields, spanning the biological, physical, and social sciences. Our results show no single field is achieving complete consideration of model uncertainties, but together we can fill the gaps. We propose opportunities to improve the quantification of uncertainty through use of a source framework for uncertainty consideration, model type specific guidelines, improved presentation, and shared best practice. We also identify shared outstanding challenges (uncertainty in input data, balancing trade-offs, error propagation, and defining how much uncertainty is required). Finally, we make nine concrete recommendations for current practice (following good practice guidelines and an uncertainty checklist, presenting uncertainty numerically, and propagating model-related uncertainty into conclusions), future research priorities (uncertainty in input data, quantifying uncertainty in complex models, and the importance of missing uncertainty in different contexts), and general research standards across the sciences (transparency about study limitations and dedicated uncertainty sections of manuscripts). △ Less

Submitted 1 July, 2022; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 40 Pages (including supporting information), 3 Figures, 2 Boxes, 1 Table

Misinformation Mitigation under Differential Propagation Rates and Temporal Penalties

Authors: Michael Simpson, Farnoosh Hashemi, Laks V. S. Lakshmanan

Abstract: We propose an information propagation model that captures important temporal aspects that have been well observed in the dynamics of fake news diffusion, in contrast with the diffusion of truth. The model accounts for differential propagation rates of truth and misinformation and for user reaction times. We study a time-sensitive variant of the \textit{misinformation mitigation} problem, where… ▽ More We propose an information propagation model that captures important temporal aspects that have been well observed in the dynamics of fake news diffusion, in contrast with the diffusion of truth. The model accounts for differential propagation rates of truth and misinformation and for user reaction times. We study a time-sensitive variant of the \textit{misinformation mitigation} problem, where $k$ seeds are to be selected to activate a truth campaign so as to minimize the number of users that adopt misinformation propagating through a social network. We show that the resulting objective is non-submodular and employ a sandwiching technique by defining submodular upper and lower bounding functions, providing data-dependent guarantees. In order to enable the use of a reverse sampling framework, we introduce a weighted version of reverse reachability sets that captures the associated differential propagation rates and establish a key equivalence between weighted set coverage probabilities and mitigation with respect to the sandwiching functions. Further, we propose an offline reverse sampling framework that provides $(1 - 1/e - ε)$-approximate solutions to our bounding functions and introduce an importance sampling technique to reduce the sample complexity of our solution. Finally, we show how our framework can provide an anytime solution to the problem. Experiments over five datasets show that our approach outperforms previous approaches and is robust to uncertainty in the model parameters. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 19 pages, VLDB 2022

How Cosmic Rays Mediate the Evolution of the Interstellar Medium

Authors: Christine M. Simpson, Rüdiger Pakmor, Christoph Pfrommer, Simon C. O. Glover, Rowan Smith

Abstract: We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or e… ▽ More We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or explosions in regions that have been cleared by prior SN, stellar winds, or radiation. We explore various mixed models by adjusting these fractions relative to each other. We find that in the simple system of a periodic stratified gas layer the ISM structure will evolve to one of two solutions: a "peak driving" state where warm gas is volume filling or a "thermal runaway" state where hot gas is volume filling. CR pressure and transport are important factors that strongly influence the solution state the ISM reaches and have the ability to flip the ISM between solutions. Observable signatures such as gamma ray emission and HI gas are explored. We find that gamma ray luminosity from pion decay is largely consistent with observations for a range of model parameters. The thickness of the HI gas layer may be too compact, however, this may be due to a large cold neutral fraction of midplane gas. The volume fraction of hot gas evolves to stable states in both solutions, but neither settles to a Milky Way-like configuration, suggesting that additional physics which is omitted here (e.g. a cosmological circum-galactic medium, radiation transport, or spectrally resolved and spatially varying CR transport) may be required. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 27 pages, 21 Figures, submitted to MNRAS, comments welcome

Formation and fate of low metallicity stars in IllustrisTNG50

Authors: Ruediger Pakmor, Christine M. Simpson, Freeke van de Voort, Lars Hernquist, Lieke van Son, Martyna Chruślińska, Rebekka Bieri, Selma E. de Mink, Volker Springel

Abstract: Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in… ▽ More Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in particular to early enrichment and to later accretion and mixing of lower metallicity gas. Because low metallicity stars are difficult to observe directly, cosmological simulations are crucial for understanding their formation. Here we quantify the rates and locations of low metallicity star formation using the high-resolution TNG50 magnetohydrodynamical cosmological simulation, and we examine where low metallicity stars end up at $z=0$. We find that $20\%$ of stars with $Z_*<0.1\,\mathrm{Z_\odot}$ form after $z=2$, and that such stars are still forming in galaxies of all masses at $z=0$ today. Moreover, most low-metallicity stars at $z=0$ reside in massive galaxies. We analyse the radial distribution of low metallicity star formation, and discuss the curious case of seven galaxies in TNG50 that form stars from primordial gas even at $z=0$. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 15 pages, 11 figures, accepted by MNRAS, comments welcome

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

Authors: Scott W McCue, Bronwyn H Bradshaw-Hajek, Matthew J Simpson

Abstract: 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… ▽ More 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. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 9 pages, 2 figures

MSC Class: 35-XX

An ALMA search for high albedo objects among the mid-sized Jupiter Trojan population

Authors: Anna M. Simpson, Michael E. Brown, Madeline J. Schemel, Bryan J. Butler

Abstract: We use ALMA measurements of 870 $μ$m thermal emission from a sample of mid-sized (15-40 km diameter) Jupiter Trojan asteroids to search for high albedo objects in this population. We calculate the diameters and albedos of each object using a thermal model which also incorporates {contemporaneous} Zwicky Transient Facility photometry to accurately measure the absolute magnitude at the time of the A… ▽ More We use ALMA measurements of 870 $μ$m thermal emission from a sample of mid-sized (15-40 km diameter) Jupiter Trojan asteroids to search for high albedo objects in this population. We calculate the diameters and albedos of each object using a thermal model which also incorporates {contemporaneous} Zwicky Transient Facility photometry to accurately measure the absolute magnitude at the time of the ALMA observation. We find that while many albedos are lower than reported from WISE, several small Trojans have high albedos independently measured both from ALMA and from WISE. The number of these high albedo objects is approximately consistent with expectations of the number of objects that recently have undergone large-scale impacts, suggesting that the interiors of freshly-crated Jupiter Trojans could contain high albedo materials such as ices. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: AJ, in press

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

Authors: Alexander K. Y. Tam, Matthew J. Simpson

Abstract: 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 interf… ▽ More 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. △ Less

Submitted 16 January, 2022; originally announced January 2022.

Comments: 20 pages, 5 figures

MSC Class: 35R37; 92D25

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

Authors: Yifei Li, Pascal R. Buenzli, Matthew J. Simpson

Abstract: 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 employ… ▽ More 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. △ Less

Submitted 6 January, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 40 pages, 11 figures, 1 supplementary material document

MSC Class: 92B05

Curvature dependences of wave propagation in reaction-diffusion models

Authors: Pascal R. Buenzli, Matthew J. Simpson

Abstract: Reaction-diffusion waves in multiple spatial dimensions advance at a rate that strongly depends on the curvature of the wave fronts. These waves have important applications in many physical, ecological, and biological systems. In this work, we analyse curvature dependences of travelling fronts in a single reaction-diffusion equation with general reaction term. We derive an exact, non-perturbative… ▽ More Reaction-diffusion waves in multiple spatial dimensions advance at a rate that strongly depends on the curvature of the wave fronts. These waves have important applications in many physical, ecological, and biological systems. In this work, we analyse curvature dependences of travelling fronts in a single reaction-diffusion equation with general reaction term. We derive an exact, non-perturbative curvature dependence of the speed of travelling fronts that arises from transverse diffusion occurring parallel to the wave front. Inward-propagating waves are characterised by three phases: an establishment phase dominated by initial and boundary conditions, a travelling-wave-like phase in which normal velocity matches standard results from singular perturbation theory, and a dip-filling phase where the collision and interaction of fronts create additional curvature dependences to their progression rate. We analyse these behaviours and additional curvature dependences using a combination of asymptotic analyses and numerical simulations. △ Less

Submitted 29 October, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 29 pages, 14 figures; v3: minor changes; additional discussion of establishment and travelling phases; new figure (Fig 9) showing simulations on a larger circular pore

Parameter estimation and uncertainty quantification using information geometry

Authors: Jesse A Sharp, Alexander P Browning, Kevin Burrage, Matthew J Simpson

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 59 pages (exc. references), 15 figures. Review

MSC Class: 62-08

A continuum mathematical model of substrate-mediated tissue growth

Authors: Maud El-Hachem, Scott W McCue, Matthew J Simpson

Abstract: We consider a continuum mathematical model of biological tissue formation inspired by recent experiments describing thin tissue growth in 3D-printed bioscaffolds. The continuum model involves a partial differential equation describing the density of tissue, $\hat{u}(\hat{\mathbf{x}},\hat{t})$, that is coupled to the concentration of an immobile extracellular substrate,… ▽ More We consider a continuum mathematical model of biological tissue formation inspired by recent experiments describing thin tissue growth in 3D-printed bioscaffolds. The continuum model involves a partial differential equation describing the density of tissue, $\hat{u}(\hat{\mathbf{x}},\hat{t})$, that is coupled to the concentration of an immobile extracellular substrate, $\hat{s}(\hat{\mathbf{x}},\hat{t})$. Cell migration is modelled with a nonlinear diffusion term, where the diffusive flux is proportional to $\hat{s}$, while a logistic growth term models cell proliferation. The extracellular substrate $\hat{s}$ is produced by cells, and undergoes linear decay. Preliminary numerical simulations show that this mathematical model, which we call the \textit{substrate model}, is able to recapitulate key features of recent tissue growth experiments, including the formation of sharp fronts. To provide a deeper understanding of the model we then analyse travelling wave solutions of the substrate model, showing that the model supports both sharp-fronted travelling wave solutions that move with a minimum wave speed, $c = c_{\rm{min}}$, as well as smooth-fronted travelling wave solutions that move with a faster travelling wave speed, $c > c_{\rm{min}}$. We provide a geometric interpretation that explains the difference between smooth- and sharp-fronted travelling wave solutions that is based on a slow manifold reduction of the desingularised three-dimensional phase space. In addition to exploring the nature of the smooth- and sharp-fronted travelling waves, we also develop and test a series of useful approximations that describe the shape of the travelling wave solutions in various limits. These approximations apply to both the sharp-fronted travelling wave solutions, and the smooth-fronted travelling wave solutions. Software to implement all calculations is available on GitHub. △ Less

Submitted 22 November, 2021; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: 46 pages, 10 figures, 1 supplementary material document

MSC Class: 92Bxx

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

Authors: David J. Warne, Thomas P. Prescott, Ruth E. Baker, Matthew J. Simpson

Abstract: 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 par… ▽ More 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. △ Less

Submitted 1 June, 2022; v1 submitted 26 October, 2021; originally announced October 2021.

MSC Class: 65C05

Mean exit time in irregularly-shaped annular and composite disc domains

Authors: Elliot J. Carr, Daniel J. VandenHeuvel, Joshua M. Wilson, Matthew J. Simpson

Abstract: Calculating the mean exit time (MET) for models of diffusion is a classical problem in statistical physics, with various applications in biophysics, economics and heat and mass transfer. While many exact results for MET are known for diffusion in simple geometries involving homogeneous materials, calculating MET for diffusion in realistic geometries involving heterogeneous materials is typically l… ▽ More Calculating the mean exit time (MET) for models of diffusion is a classical problem in statistical physics, with various applications in biophysics, economics and heat and mass transfer. While many exact results for MET are known for diffusion in simple geometries involving homogeneous materials, calculating MET for diffusion in realistic geometries involving heterogeneous materials is typically limited to repeated stochastic simulations or numerical solutions of the associated boundary value problem (BVP). In this work we derive exact solutions for the MET in irregular annular domains, including some applications where diffusion occurs in heterogenous media. These solutions are obtained by taking the exact results for MET in an annulus, and then constructing various perturbation solutions to account for the irregular geometries involved. These solutions, with a range of boundary conditions, are implemented symbolically and compare very well with averaged data from repeated stochastic simulations and with numerical solutions of the associated BVP. Software to implement the exact solutions is available at https://github.com/ProfMJSimpson/Exit_time. △ Less

Submitted 2 March, 2022; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 18 pages, 6 figures, accepted version

Journal ref: Journal of Physics A: Mathematical and Theoretical, 55 (2022) 105002

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

Authors: Maud El-Hachem, Scott W McCue, Matthew J Simpson

Abstract: 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,… ▽ More 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. △ Less

Submitted 23 January, 2022; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 41 pages, 16 figures

MSC Class: 92Bxx

Understanding the chemical mechanism behind photo-induced enhanced Raman spectroscopy

Authors: Junzhi Ye, Rakesh Arul, Michel K. Nieuwoudt, Junzhe Dong, Wei Gao, M. Cather Simpson

Abstract: Photo-Induced Enhanced Raman Spectroscopy (PIERS) is a new surface enhanced Raman spectroscopy (SERS) modality with an order-of-magnitude Raman signal enhancement of adsorbed analytes over that of typical SERS substrates. Despite the impressive PIERS enhancement factors and explosion in recent demonstrations of its utility, the detailed enhancement mechanism remains undetermined. Using a range of… ▽ More Photo-Induced Enhanced Raman Spectroscopy (PIERS) is a new surface enhanced Raman spectroscopy (SERS) modality with an order-of-magnitude Raman signal enhancement of adsorbed analytes over that of typical SERS substrates. Despite the impressive PIERS enhancement factors and explosion in recent demonstrations of its utility, the detailed enhancement mechanism remains undetermined. Using a range of optical and X-ray spectroscopies, supported by density functional theory calculations, we elucidate the chemical and atomic-scale mechanism behind the PIERS enhancement. Stable PIERS substrates with enhancement factors of 10^6 were fabricated using self-organized hexagonal arrays of TiO2 nanotubes that were defect-engineered via annealing in inert atmospheres, and silver nanoparticles were deposited by magnetron sputtering and subsequent thermal dewetting. We identified the key source of the enhancement of PIERS vs. SERS in these structures as an increase in the Raman polarizability of the adsorbed probe molecule upon photo-induced charge transfer. A balance between crystallinity, which enhances charge transfer due to higher electron mobility in anatase-rutile heterostructures but decreases visible light absorption, and oxygen vacancy defects, which increase visible light absorption and photo-induced electron transfers, was critical to achieve high PIERS enhancements. △ Less

Submitted 30 May, 2021; originally announced May 2021.

Simulating radio synchrotron emission in star-forming galaxies: small-scale magnetic dynamo and the origin of the far infrared-radio correlation

Authors: Christoph Pfrommer, Maria Werhahn, Rüdiger Pakmor, Philipp Girichidis, Christine M. Simpson

Abstract: In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). We examine magneto-hydrodynamic galaxy simulations with cosmic rays (CRs), accounting for their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shoc… ▽ More In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). We examine magneto-hydrodynamic galaxy simulations with cosmic rays (CRs), accounting for their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shock, which injects turbulence and drives a small-scale magnetic dynamo. As the shock propagates outwards and the associated turbulence decays, the large velocity shear between the supersonically rotating cool disc with respect to the (partially) pressure-supported hot circumgalactic medium excites Kelvin-Helmholtz surface and body modes. Those inject turbulence and drive multiple small-scale dynamos, which exponentially amplify magnetic fields. They grow in scale to reach equipartition with thermal and CR energies in Milky Way-mass galaxies. In small galaxies, the magnetic energy saturates at the turbulent energy while it fails to reach equipartition with thermal and CR energies. We solve for steady-state spectra of CR protons, secondary electrons/positrons from hadronic CR-proton interactions with the interstellar medium, and primary shock-accelerated electrons at supernovae. The radio-synchrotron emission is dominated by primary electrons, irradiates the magnetised disc, bulge, and bubble-shaped magnetically-loaded outflows of our simulated Milky Way-mass galaxy. Our star-forming and star-bursting galaxies with saturated magnetic fields match the global FIR-radio correlation (FRC) across four orders of magnitude. Its intrinsic scatter arises due to (i) different magnetic saturation levels that result from different seed magnetic fields, (ii) different radio synchrotron luminosities for different specific SFRs at fixed SFR and (iii) a varying radio intensity with galactic inclination. (abridged) △ Less

Submitted 5 July, 2022; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 37 pages, 26 figures, accepted for publication in MNRAS

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

Authors: Elliot J. Carr, Dylan J. Oliver, Matthew J. Simpson

Abstract: 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. E… ▽ More 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. △ Less

Submitted 10 January, 2022; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: 14 pages, 3 figures, accepted version of paper published in Communications in Nonlinear Science and Numerical Simulation

Journal ref: Communications in Nonlinear Science and Numerical Simulation, 108 (2022) 106211

Satellites Around Milky Way Analogs: Tension in the Number and Fraction of Quiescent Satellites Seen in Observations Versus Simulations

Authors: Ananthan Karunakaran, Kristine Spekkens, Kyle A. Oman, Christine M. Simpson, Azadeh Fattahi, David J. Sand, Paul Bennet, Denija Crnojević, Carlos S. Frenk, Facundo A. Gómez, Robert J. J. Grand, Michael G. Jones, Federico Marinacci, Burçin Mutlu-Pakdil, Julio F. Navarro, Dennis Zaritsky

Abstract: We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APO… ▽ More We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APOSTLE and Auriga. We compare our detection rates from the NUV and FUV bands to the SAGA-II H$α$ detections and find that they are broadly consistent with over $85\%$ of observed satellites detected in all three tracers. We apply the same spatial selection criteria used around SAGA-II hosts to select satellites around the MW-like hosts in APOSTLE and Auriga. We find very good overall agreement in the derived SFRs for the star-forming satellites as well as the number of star-forming satellites per host in observed and simulated samples. However, the number and fraction of quenched satellites in the SAGA-II sample are significantly lower than those in APOSTLE and Auriga below a stellar mass of $M_*\sim10^{8}\,M_{\odot}$, even when the SAGA-II incompleteness and interloper corrections are included. This discrepancy is robust with respect to the resolution of the simulations and persists when alternative star-formation tracers are employed. We posit that this disagreement is not readily explained by vagaries in the observed or simulated samples considered here, suggesting a genuine discrepancy that may inform the physics of satellite populations around MW analogs. △ Less

Submitted 19 July, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: Accepted for publication in ApJ Letters after minor changes to text

Linking the brightest stellar streams with the accretion history of Milky Way-like galaxies

Authors: Alex Vera-Casanova, Facundo A. Gómez, Antonela Monachesi, Ignacio Gargiulo, Diego Pallero, Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Carlos S. Frenk, Gustavo Morales

Abstract: According to the current galaxy formation paradigm, mergers and interactions play an important role in shaping present-day galaxies. The remnants of this merger activity can be used to constrain galaxy formation models. In this work we use a sample of thirty hydrodynamical simulations of Milky Way-mass halos, from the AURIGA project, to generate surface brightness maps and search for the brightest… ▽ More According to the current galaxy formation paradigm, mergers and interactions play an important role in shaping present-day galaxies. The remnants of this merger activity can be used to constrain galaxy formation models. In this work we use a sample of thirty hydrodynamical simulations of Milky Way-mass halos, from the AURIGA project, to generate surface brightness maps and search for the brightest stream in each halo as a function of varying limiting magnitude. We find that none of the models shows signatures of stellar streams at $μ_{r}^{lim} \leq 25$ mag arcsec$^{-2}$. The stream detection increases significantly between 27 and 28 mag arcsec$^{-2}$. Nevertheless, even at 30 mag arcsec$^{-2}$, 13 percent of our models show no detectable streams. We study the properties of the brightest streams progenitors (BSPs). We find that BSPs are accreted within a broad range of infall times, from 1.6 to 10 Gyr ago, with only 25 percent accreted within the last 5 Gyrs; thus most BSPs correspond to relatively early accretion events. We also find that 37 percent of the BSPs survive to the present day. The median infall times for surviving and disrupted BSPs are 5.6 and 6.7 Gyr, respectively. We find a clear relation between infall time and infall mass of the BSPs, such that more massive progenitors tend to be accreted at later times. However, we find that the BSPs are not, in most cases, the dominant contributor to the accreted stellar halo of each galaxy. △ Less

Submitted 9 June, 2022; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS, 15 pages, 11 figures

Travelling wave analysis of cellular invasion into surrounding tissues

Authors: Maud El-Hachem, Scott W McCue, Matthew J Simpson

Abstract: Single-species reaction-diffusion equations, such as the Fisher-KPP and Porous-Fisher equations, support travelling wave solutions that are often interpreted as simple mathematical models of biological invasion. Such travelling wave solutions are thought to play a role in various applications including development, wound healing and malignant invasion. One criticism of these single-species equatio… ▽ More Single-species reaction-diffusion equations, such as the Fisher-KPP and Porous-Fisher equations, support travelling wave solutions that are often interpreted as simple mathematical models of biological invasion. Such travelling wave solutions are thought to play a role in various applications including development, wound healing and malignant invasion. One criticism of these single-species equations is that they do not explicitly describe interactions between the invading population and the surrounding environment. In this work we study a reaction-diffusion equation that describes malignant invasion which has been used to interpret experimental measurements describing the invasion of malignant melanoma cells into surrounding human skin tissues. This model explicitly describes how the population of cancer cells degrade the surrounding tissues, thereby creating free space into which the cancer cells migrate and proliferate to form an invasion wave of malignant tissue that is coupled to a retreating wave of skin tissue. We analyse travelling wave solutions of this model using a combination of numerical simulation, phase plane analysis and perturbation techniques. Our analysis shows that the travelling wave solutions involve a range of very interesting properties that resemble certain well-established features of both the Fisher-KPP and Porous-Fisher equations, as well as a range of novel properties that can be thought of as extensions of these well-studied single-species equations. Of particular interest is that travelling wave solutions of the invasion model are very well approximated by trajectories in the Fisher-KPP phase plane that are normally disregarded. This observation establishes a previously unnoticed link between coupled multi-species reaction diffusion models of invasion and a different class of models of invasion that involve moving boundary problems. △ Less

Submitted 19 July, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: 30 pages, 8 figures

MSC Class: 92Bxx

Determining the full satellite population of a Milky Way-mass halo in a highly resolved cosmological hydrodynamic simulation

Authors: Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Ashley J. Kelly, Facundo A. Gómez, Adrian Jenkins, Volker Springel, Carlos S. Frenk, Simon D. M. White

Abstract: We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magneto-hydrodynamic cosmological zoom-in simulation (baryonic mass resolution $m_b =$ 800 $\rm M_{\odot}$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are: i) robust to stochastic vari… ▽ More We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magneto-hydrodynamic cosmological zoom-in simulation (baryonic mass resolution $m_b =$ 800 $\rm M_{\odot}$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are: i) robust to stochastic variations associated with the so-called ``Butterfly Effect''; and ii) well converged over 3.5 orders of magnitude in mass resolution. We find that there are approximately five times as many satellite galaxies at this high resolution compared to a standard ($m_b\sim 10^{4-5}\, \rm M_{\odot}$) resolution simulation of the same system. This is primarily because 2/3rds of the high resolution satellites do not form at standard resolution. A smaller fraction (1/6th) of the satellites present at high resolution form and disrupt at standard resolution; these objects are preferentially low-mass satellites on intermediate- to low-eccentricity orbits with impact parameters $\lesssim 30$ kpc. As a result, the radial distribution of satellites becomes substantially more centrally concentrated at higher resolution, in better agreement with recent observations of satellites around Milky Way-mass haloes. Finally, we show that our galaxy formation model successfully forms ultra-faint galaxies and reproduces the stellar velocity dispersion, half-light radii, and $V$-band luminosities of observed Milky Way and Local Group dwarf galaxies across 6 orders of magnitude in luminosity ($10^3$-$10^{9}$ $\rm L_{\odot}$). △ Less

Submitted 3 September, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: accepted for publication in MNRAS. Main changes include new figure 4 and figure 11

Revisiting the Color-Color Selection: Submillimeter and AGN Properties of NUV-r-J Selected Quiescent Galaxies

Authors: Yu-Hsuan Hwang, Wei-Hao Wang, Yu-Yen Chang, Chen-Fatt Lim, Chian-Chou Chen, Zhen-Kai Gao, James S. Dunlop, Yu Gao, Luis C. Ho, Ho Seong Hwang, Maciej Koprowski, Michał J. Michałowski, Ying-jie Peng, Hyunjin Shim, James M. Simpson, Yoshiki Toba

Abstract: We examine the robustness of the color-color selection of quiescent galaxies (QGs) against contamination of dusty star-forming galaxies using the latest submillimeter data. We selected 18,304 QG candidates out to $z\sim$ 3 using the commonly adopted $NUV-r-J$ selection based on the high-quality multi-wavelength COSMOS2015 catalog. Using extremely deep 450 and 850 $μ$m catalogs from the latest JCMT… ▽ More We examine the robustness of the color-color selection of quiescent galaxies (QGs) against contamination of dusty star-forming galaxies using the latest submillimeter data. We selected 18,304 QG candidates out to $z\sim$ 3 using the commonly adopted $NUV-r-J$ selection based on the high-quality multi-wavelength COSMOS2015 catalog. Using extremely deep 450 and 850 $μ$m catalogs from the latest JCMT SCUBA-2 Large Programs, S2COSMOS, and STUDIES, as well as ALMA submillimeter, VLA 3 GHz, and $Spitzer$ MIPS 24 $μ$m catalogs, we identified luminous dusty star-forming galaxies among the QG candidates. We also conducted stacking analyses in the SCUBA-2 450 and 850 $μ$m images to look for less-luminous dusty galaxies among the QG candidates. By cross-matching to the 24 $μ$m and 3 GHz data, we were able to identify a sub-group of "IR-radio-bright" QGs who possess a strong 450 and 850 $μ$m stacking signal. The potential contamination of these luminous and less-luminous dusty galaxies accounts for approximately 10% of the color-selected QG candidates. In addition, there exists a spatial correlation between the luminous star-forming galaxies and the QGs at a $\lesssim60$ kpc scale. Finally, we found a high QG fraction among radio AGNs at $z<$ 1.5. Our data show a strong correlation between QGs and radio AGNs, which may suggest a connection between the quenching process and the radio-mode AGN feedback. △ Less

Submitted 29 March, 2021; v1 submitted 26 March, 2021; originally announced March 2021.

Comments: This paper is accepted for publication on ApJ

Travelling waves, blow-up and extinction in the Fisher-Stefan model

Authors: Scott W. McCue, Maud El-Hachem, Matthew J. Simpson

Abstract: While there is a long history of employing moving boundary problems in physics, in particular via Stefan problems for heat conduction accompanied by a change of phase, more recently such approaches have been adapted to study biological invasion. For example, when a logistic growth term is added to the governing partial differential equation in a Stefan problem, one arrives at the Fisher-Stefan mod… ▽ More While there is a long history of employing moving boundary problems in physics, in particular via Stefan problems for heat conduction accompanied by a change of phase, more recently such approaches have been adapted to study biological invasion. For example, when a logistic growth term is added to the governing partial differential equation in a Stefan problem, one arrives at the Fisher-Stefan model, a generalisation of the well-known Fisher-KPP model, characterised by a leakage coefficient $κ$ which relates the speed of the moving boundary to the flux of population there. This Fisher-Stefan model overcomes one of the well-known limitations of the Fisher-KPP model, since time-dependent solutions of the Fisher-Stefan model involve a well-defined front with compact support which is more natural in terms of mathematical modelling. Almost all of the existing analysis of the standard Fisher-Stefan model involves setting $κ> 0$, which can lead to either invading travelling wave solutions or complete extinction of the population. Here, we demonstrate how setting $κ< 0$ leads to retreating travelling waves and an interesting transition to finite-time blow-up. For certain initial conditions, population extinction is also observed. Our approach involves studying time-dependent solutions of the governing equations, phase plane and asymptotic analysis, leading to new insight into the possibilities of travelling waves, blow-up and extinction for this moving boundary problem. Matlab software used to generate the results in this work are available on Github. △ Less

Submitted 22 September, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Bayesian uncertainty quantification for data-driven equation learning

Authors: Simon Martina-Perez, Matthew J. Simpson, Ruth E. Baker

Abstract: Equation learning aims to infer differential equation models from data. While a number of studies have shown that differential equation models can be successfully identified when the data are sufficiently detailed and corrupted with relatively small amounts of noise, the relationship between observation noise and uncertainty in the learned differential equation models remains unexplored. We demons… ▽ More Equation learning aims to infer differential equation models from data. While a number of studies have shown that differential equation models can be successfully identified when the data are sufficiently detailed and corrupted with relatively small amounts of noise, the relationship between observation noise and uncertainty in the learned differential equation models remains unexplored. We demonstrate that for noisy data sets there exists great variation in both the structure of the learned differential equation models as well as the parameter values. We explore how to combine data sets to quantify uncertainty in the learned models, and at the same time draw mechanistic conclusions about the target differential equations. We generate noisy data using a stochastic agent-based model and combine equation learning methods with approximate Bayesian computation (ABC) to show that the correct differential equation model can be successfully learned from data, while a quantification of uncertainty is given by a posterior distribution in parameter space. △ Less

Submitted 29 September, 2021; v1 submitted 23 February, 2021; originally announced February 2021.