Showing 51–100 of 428 results for author: Simpson, J

Anthropogenic Mixing of Seasonally Stratified Shelf Seas by Offshore Wind Farm Infrastructure

Authors: Robert Dorrell, Charlie Lloyd, Ben Lincoln, Tom Rippeth, John Taylor, Colm-cille Caulfield, Jonathan Sharples, Jeff Polton, Brian Scannell, Deborah Greaves, Rob Hall, John Simpson

Abstract: The offshore wind energy sector has rapidly expanded over the past two decades, providing a renewable energy solution for coastal nations. Sector development has been led in Europe, but is growing globally. Most developments to date have been in well-mixed, i.e. unstratified, shallow-waters near to shore. Sector growth is, for the first time, pushing developments to deep water, into a brand new en… ▽ More The offshore wind energy sector has rapidly expanded over the past two decades, providing a renewable energy solution for coastal nations. Sector development has been led in Europe, but is growing globally. Most developments to date have been in well-mixed, i.e. unstratified, shallow-waters near to shore. Sector growth is, for the first time, pushing developments to deep water, into a brand new environment: seasonally stratified shelf seas. Seasonally stratified shelf seas, where water density varies with depth, have a disproportionately key role in primary production, marine ecosystem and biochemically cycles. Infrastructure will directly mix stratified shelf seas. The magnitude of this mixing, additional to natural background processes, has yet to be fully quantified. If large enough it may erode shelf sea stratification. Therefore, offshore wind growth may destabilize and fundamentally change shelf sea systems. However, enhanced mixing may also positively impact some marine ecosystems. This paper sets the scene for sector development into this new environment, reviews the potential physical and environmental benefits and impacts of large scale industrialization of seasonally stratified shelf seas and identifies areas where research is required to best utilise, manage and mitigate environmental change. △ Less

Submitted 18 December, 2021; originally announced December 2021.

Comments: 31 pages, 14 figures

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

Handling missing data when estimating causal effects with Targeted Maximum Likelihood Estimation

Authors: S. Ghazaleh Dashti, Katherine J. Lee, Julie A. Simpson, Ian R. White, John B. Carlin, Margarita Moreno-Betancur

Abstract: Targeted Maximum Likelihood Estimation (TMLE) is increasingly used for doubly robust causal inference, but how missing data should be handled when using TMLE with data-adaptive approaches is unclear. Based on the Victorian Adolescent Health Cohort Study, we conducted a simulation study to evaluate eight missing data methods in this context: complete-case analysis, extended TMLE incorporating outco… ▽ More Targeted Maximum Likelihood Estimation (TMLE) is increasingly used for doubly robust causal inference, but how missing data should be handled when using TMLE with data-adaptive approaches is unclear. Based on the Victorian Adolescent Health Cohort Study, we conducted a simulation study to evaluate eight missing data methods in this context: complete-case analysis, extended TMLE incorporating outcome-missingness model, missing covariate missing indicator method, five multiple imputation (MI) approaches using parametric or machine-learning models. Six scenarios were considered, varying in exposure/outcome generation models (presence of confounder-confounder interactions) and missingness mechanisms (whether outcome influenced missingness in other variables and presence of interaction/non-linear terms in missingness models). Complete-case analysis and extended TMLE had small biases when outcome did not influence missingness in other variables. Parametric MI without interactions had large bias when exposure/outcome generation models included interactions. Parametric MI including interactions performed best in bias and variance reduction across all settings, except when missingness models included a non-linear term. When choosing a method to handle missing data in the context of TMLE, researchers must consider the missingness mechanism and, for MI, compatibility with the analysis method. In many settings, a parametric MI approach that incorporates interactions and non-linearities is expected to perform well. △ Less

Submitted 3 May, 2024; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 31 pages, 2 tables, 5 figures, 9 supplementary tables

Journal ref: Am J Epidemiol. 2024 Feb 22:kwae012. Epub ahead of print. PMID: 38400653

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

The GALAH Survey: Improving our understanding of confirmed and candidate planetary systems with large stellar surveys

Authors: Jake T. Clark, Duncan J. Wright, Robert A. Wittenmyer, Jonathan Horner, Natalie R. Hinkel, Mathieu Clerté, Brad D. Carter, Sven Buder, Michael R. Hayden, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaž Zwitter , et al. (2 additional authors not shown)

Abstract: Pioneering photometric, astrometric, and spectroscopic surveys are helping exoplanetary scientists better constrain the fundamental properties of stars within our galaxy, and the planets these stars host. In this study, we use the third data release from the stellar spectroscopic GALAH Survey, coupled with astrometric data of eDR3 from the \textit{Gaia} satellite, and other data from NASA's Exopla… ▽ More Pioneering photometric, astrometric, and spectroscopic surveys are helping exoplanetary scientists better constrain the fundamental properties of stars within our galaxy, and the planets these stars host. In this study, we use the third data release from the stellar spectroscopic GALAH Survey, coupled with astrometric data of eDR3 from the \textit{Gaia} satellite, and other data from NASA's Exoplanet Archive, to refine our understanding of 279 confirmed and candidate exoplanet host stars and their exoplanets. This homogenously analysed data set comprises 105 confirmed exoplanets, along with 146 K2 candidates, 95 TESS Objects of Interest (TOIs) and 52 Community TOIs (CTOIs). Our analysis significantly shifts several previously (unknown) planet parameters while decreasing the uncertainties for others; Our radius estimates suggest that 35 planet candidates are more likely brown dwarfs or stellar companions due to their new radius values. We are able to refine the radii and masses of WASP-47 e, K2-106 b, and CoRoT-7 b to their most precise values yet, to less than 2.3\% and 8.5\% respectively. We also use stellar rotational values from GALAH to show that most planet candidates will have mass measurements that will be tough to obtain with current ground-based spectrographs. With GALAH's chemical abundances, we show through chemo-kinematics that there are five planet-hosts that are associated with the galaxy's thick disc, including NGTS-4, K2-183 and K2-337. Finally, we show there is no statistical difference between the chemical properties of hot Neptune and hot rocky exoplanet hosts, with the possibility that short-period rocky worlds might be the remnant cores of hotter, gaseous worlds. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 21 pages with 11 figures and 5 tables. Accepted into MNRAS with appendix tables available within this arXiv submission

Narrow resonances in the continuum of the unbound nucleus $^{15}$F

Authors: V. Girard-Alcindor, A. Mercenne, I. Stefan, F. de Oliveira Santos, N. Michel, M. Płoszajczak, M. Assié, A. Lemasson, E. Clément, F. Flavigny, A. Matta, D. Ramos, M. Rejmund, J. Dudouet, D. Ackermann, P. Adsley, M. Assunção, B. Bastin, D. Beaumel, G. Benzoni, R. Borcea, A. J. Boston, L. Cáceres, B. Cederwall, I. Celikovic , et al. (78 additional authors not shown)

Abstract: The structure of the unbound $^{15}$F nucleus is investigated using the inverse kinematics resonant scattering of a radioactive $^{14}$O beam impinging on a CH$_2$ target. The analysis of $^{1}$H($^{14}$O,p)$^{14}$O and $^{1}$H($^{14}$O,2p)$^{13}$N reactions allowed the confirmation of the previously observed narrow $1/2^{-}$ resonance, near the two-proton decay threshold, and the identification o… ▽ More The structure of the unbound $^{15}$F nucleus is investigated using the inverse kinematics resonant scattering of a radioactive $^{14}$O beam impinging on a CH$_2$ target. The analysis of $^{1}$H($^{14}$O,p)$^{14}$O and $^{1}$H($^{14}$O,2p)$^{13}$N reactions allowed the confirmation of the previously observed narrow $1/2^{-}$ resonance, near the two-proton decay threshold, and the identification of two new narrow 5/2$^{-}$ and 3/2$^{-}$ resonances. The newly observed levels decay by 1p emission to the ground of $^{14}$O, and by sequential 2p emission to the ground state (g.s.) of $^{13}$N via the $1^-$ resonance of $^{14}$O. Gamow shell model (GSM) analysis of the experimental data suggests that the wave functions of the 5/2$^{-}$ and 3/2$^{-}$ resonances may be collectivized by the continuum coupling to nearby 2p- and 1p- decay channels. The observed excitation function $^{1}$H($^{14}$O,p)$^{14}$O and resonance spectrum in $^{15}$F are well reproduced in the unified framework of the GSM. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 7 pages, 4 figures

Robust Devaney chaos in the two-dimensional border-collision normal form

Authors: Indranil Ghosh, David J. W. Simpson

Abstract: The collection of all non-degenerate, continuous, two-piece, piecewise-linear maps on $\mathbb{R}^2$ can be reduced to a four-parameter family known as the two-dimensional border-collision normal form. We prove that throughout an open region of parameter space this family has an attractor satisfying Devaney's definition of chaos. This strengthens existing results on the robustness of chaos in piec… ▽ More The collection of all non-degenerate, continuous, two-piece, piecewise-linear maps on $\mathbb{R}^2$ can be reduced to a four-parameter family known as the two-dimensional border-collision normal form. We prove that throughout an open region of parameter space this family has an attractor satisfying Devaney's definition of chaos. This strengthens existing results on the robustness of chaos in piecewise-linear maps. We further show that the stable manifold of a saddle fixed point, despite being a one-dimensional object, densely fills an open region containing the attractor. Finally we identify a heteroclinic bifurcation, not described previously, at which the attractor undergoes a crisis and may be destroyed. △ Less

Submitted 24 November, 2021; originally announced November 2021.

MSC Class: 37G35; 39A28

Inclusion of higher-order terms in the border-collision normal form: persistence of chaos and applications to power converters

Authors: David J. W. Simpson, Paul A. Glendinning

Abstract: The dynamics near a border-collision bifurcation are approximated to leading order by a continuous, piecewise-linear map. The purpose of this paper is to consider the higher-order terms that are neglected when forming this approximation. For two-dimensional maps we establish conditions under which a chaotic attractor created in a border-collision bifurcation persists for an open interval of parame… ▽ More The dynamics near a border-collision bifurcation are approximated to leading order by a continuous, piecewise-linear map. The purpose of this paper is to consider the higher-order terms that are neglected when forming this approximation. For two-dimensional maps we establish conditions under which a chaotic attractor created in a border-collision bifurcation persists for an open interval of parameters beyond the bifurcation. We apply the results to a prototypical power converter model to prove the model exhibits robust chaos. △ Less

Submitted 23 November, 2021; originally announced November 2021.

MSC Class: 39A28; 37G35

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

Pattern Formation in a Spatially-Extended Model of Pacemaker Dynamics in Smooth Muscle Cells

Authors: H. O. Fatoyinbo, R. G. Brown, D. J. W. Simpson, B. van Brunt

Abstract: Spatiotemporal patterns are common in biological systems. For electrically-coupled cells previous studies of pattern formation have mainly used external forcing as the main bifurcation parameter. The purpose of this paper is to show that spatiotemporal patterns in electrically-coupled smooth muscle cells occur even in the absence of forcing. We study a reaction-diffusion system with the Morris-Lec… ▽ More Spatiotemporal patterns are common in biological systems. For electrically-coupled cells previous studies of pattern formation have mainly used external forcing as the main bifurcation parameter. The purpose of this paper is to show that spatiotemporal patterns in electrically-coupled smooth muscle cells occur even in the absence of forcing. We study a reaction-diffusion system with the Morris-Lecar equations and observe a wide range of spatiotemporal patterns for different values of the model parameters. Some aspects of these patterns are explained via a bifurcation analysis of the system without coupling -- in particular Type I and Type II excitability both occur. We show the patterns are not due to a Turing instability and use travelling wave coordinates to analyse travelling waves. △ Less

Submitted 29 October, 2021; originally announced November 2021.

MSC Class: 39A28; 37N25; 74H60; 7M20; 37C75; 92C30

Counting Subwords in Circular Words and Their Parikh Matrices

Authors: Ghajendran Poovanandran, Jamie Simpson, Wen Chean Teh

Abstract: The word inference problem is to determine languages such that the information on the number of occurrences of those subwords in the language can uniquely identify a word. A considerable amount of work has been done on this problem, but the same cannot be said for circular words despite growing interests on the latter due to their applications -- for example, in splicing systems. Meanwhile, Parikh… ▽ More The word inference problem is to determine languages such that the information on the number of occurrences of those subwords in the language can uniquely identify a word. A considerable amount of work has been done on this problem, but the same cannot be said for circular words despite growing interests on the latter due to their applications -- for example, in splicing systems. Meanwhile, Parikh matrices are useful tools and well established in the study of subword occurrences. In this work, we propose two ways of counting subword occurrences in circular words. We then extend the idea of Parikh matrices to the context of circular words and investigate this extension. Motivated by the word inference problem, we study ambiguity in the identification of a circular word by its Parikh matrix. Accordingly, two rewriting rules are developed to generate ternary circular words which share the same Parikh matrix. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: Preprint, Submitted for journal consideration

MSC Class: 68R15; 05A05

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

$S^5$: The Orbital and Chemical Properties of One Dozen Stellar Streams

Authors: Ting S. Li, Alexander P. Ji, Andrew B. Pace, Denis Erkal, Sergey E. Koposov, Nora Shipp, Gary S. Da Costa, Lara R. Cullinane, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Jeffrey D. Simpson, Daniel B. Zucker, Peter S. Ferguson, Sarah L. Martell, Joss Bland-Hawthorn, Eduardo Balbinot, Kiyan Tavangar, Alex Drlica-Wagner, Gayandhi M. De Silva, Joshua D. Simon, S5 Collaboration

Abstract: We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full… ▽ More We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ${\sim}10-50$ kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass-metallicity relation, the luminosities of the progenitors of the DG streams range between Carina and Ursa Major I ($-9.5\lesssim M_V\lesssim-5.5$). Four of the six GC streams have mean metallicities of [Fe/H]$< -2$, more metal-poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the 300S and Jet GC streams are the only streams on retrograde orbits in our dozen stream sample. Finally, we compare the orbital properties of the streams with known DGs and GCs in the MW, finding several possible associations. Some streams appear to have been accreted with the recently discovered Gaia-Enceladus-Sausage system, and others suggest that GCs were formed in and accreted together with the progenitors of DG streams whose stellar masses are similar to Draco to Carina ($\sim10^5-10^6M_\odot$). △ Less

Submitted 2 January, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 24 pages, 8 Figures, 2 Tables. Accepted for publication in ApJ

Renormalisation of the two-dimensional border-collision normal form

Authors: Indranil Ghosh, David J. W. Simpson

Abstract: We study the two-dimensional border-collision normal form (a four-parameter family of continuous, piecewise-linear maps on $\mathbb{R}^2$) in the robust chaos parameter region of [S. Banerjee, J.A. Yorke, C. Grebogi, Robust Chaos, Phys. Rev. Lett. 80(14):3049--3052, 1998]. We use renormalisation to partition this region by the number of connected components of a chaotic Milnor attractor. This reve… ▽ More We study the two-dimensional border-collision normal form (a four-parameter family of continuous, piecewise-linear maps on $\mathbb{R}^2$) in the robust chaos parameter region of [S. Banerjee, J.A. Yorke, C. Grebogi, Robust Chaos, Phys. Rev. Lett. 80(14):3049--3052, 1998]. We use renormalisation to partition this region by the number of connected components of a chaotic Milnor attractor. This reveals previously undescribed bifurcation structure in a succinct way. △ Less

Submitted 19 September, 2021; originally announced September 2021.

MSC Class: 37G35; 39A28

Agile, Antifragile, Artificial-Intelligence-Enabled, Command and Control

Authors: Jacob Simpson, Rudolph Oosthuizen, Sondoss El Sawah, Hussein Abbass

Abstract: Artificial Intelligence (AI) is rapidly becoming integrated into military Command and Control (C2) systems as a strategic priority for many defence forces. The successful implementation of AI is promising to herald a significant leap in C2 agility through automation. However, realistic expectations need to be set on what AI can achieve in the foreseeable future. This paper will argue that AI could… ▽ More Artificial Intelligence (AI) is rapidly becoming integrated into military Command and Control (C2) systems as a strategic priority for many defence forces. The successful implementation of AI is promising to herald a significant leap in C2 agility through automation. However, realistic expectations need to be set on what AI can achieve in the foreseeable future. This paper will argue that AI could lead to a fragility trap, whereby the delegation of C2 functions to an AI could increase the fragility of C2, resulting in catastrophic strategic failures. This calls for a new framework for AI in C2 to avoid this trap. We will argue that antifragility along with agility should form the core design principles for AI-enabled C2 systems. This duality is termed Agile, Antifragile, AI-Enabled Command and Control (A3IC2). An A3IC2 system continuously improves its capacity to perform in the face of shocks and surprises through overcompensation from feedback during the C2 decision-making cycle. An A3IC2 system will not only be able to survive within a complex operational environment, it will also thrive, benefiting from the inevitable shocks and volatility of war. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 12 pages, 7 figures, included in the 26th International Command and Control Research and Technology Symposium (ICCRTS)

Evidence for enhanced neutron-proton correlations from the level structure of the $N=Z+1$ nucleus $^{87}_{43}$Tc$^{\ }_{44}$

Authors: X. Liu, B. Cederwall, C. Qi, R. A. Wyss, Ö. Aktas, A. Ertoprak, W. Zhang, E. Clément, G. de France, D. Ralet, A. Gadea, A. Goasduff, G. Jaworski, I. Kuti, B. M. Nyakó, J. Nyberg, M. Palacz, R. Wadsworth, J. J. Valiente-Dobón, H. Al-Azri, A. Ataç Nyberg, T. Bäck, G. de Angelis, M. Doncel, J. Dudouet , et al. (47 additional authors not shown)

Abstract: The low-lying excited states in the neutron-deficient $N=Z+1$ nucleus $^{87}_{43}$Tc$^{\ }_{44}$ have been studied via the fusion-evaporation reaction $^{54}$Fe($^{36}$Ar, $2n1p$)$^{87}$Tc at the Grand Accélérateur National d'Ions Lourds (GANIL), France. The AGATA spectrometer was used in conjunction with the auxiliary NEDA, Neutron Wall, and DIAMANT detector arrays to measure coincident prompt… ▽ More The low-lying excited states in the neutron-deficient $N=Z+1$ nucleus $^{87}_{43}$Tc$^{\ }_{44}$ have been studied via the fusion-evaporation reaction $^{54}$Fe($^{36}$Ar, $2n1p$)$^{87}$Tc at the Grand Accélérateur National d'Ions Lourds (GANIL), France. The AGATA spectrometer was used in conjunction with the auxiliary NEDA, Neutron Wall, and DIAMANT detector arrays to measure coincident prompt $γがんま$-rays, neutrons, and charged particles emitted in the reaction. A level scheme of $^{87}$Tc from the (9/2$^{+}_{g.s.}$) state to the (33/2$^{+}_{1}$) state was established based on 6 mutually coincident $γがんま$-ray transitions. The constructed level structure exhibits a rotational behavior with a sharp backbending at $\hbarωおめが\approx 0.50$ MeV. A decrease in alignment frequency and increase in alignment sharpness in the odd-mass isotonic chains around $N=44$ is proposed as an effect of the enhanced isoscalar neutron-proton interactions in odd-mass nuclei when approaching the $N=Z$ line. △ Less

Submitted 13 September, 2021; originally announced September 2021.

The GALAH Survey: Chemical tagging and chrono-chemodynamics of accreted halo stars with GALAH+ DR3 and $Gaia$ eDR3

Authors: Sven Buder, Karin Lind, Melissa K. Ness, Diane K. Feuillet, Danny Horta, Stephanie Monty, Tobias Buck, Thomas Nordlander, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael R. Hayden, Janez Kos, Sarah L. Martell, Geraint F. Lewis, Jane Lin, Katharine. J. Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaz Zwitter, Ioana Ciuca , et al. (5 additional authors not shown)

Abstract: Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early "building block" given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we invest… ▽ More Since the advent of $Gaia$ astrometry, it is possible to identify massive accreted systems within the Galaxy through their unique dynamical signatures. One such system, $Gaia$-Sausage-Enceladus (GSE), appears to be an early "building block" given its virial mass $> 10^{10}\,\mathrm{M_\odot}$ at infall ($z\sim1-3$). In order to separate the progenitor population from the background stars, we investigate its chemical properties with up to 30 element abundances from the GALAH+ Survey Data Release 3 (DR3). To inform our choice of elements for purely chemically selecting accreted stars, we analyse 4164 stars with low-$αあるふぁ$ abundances and halo kinematics. These are most different to the Milky Way stars for abundances of Mg, Si, Na, Al, Mn, Fe, Ni, and Cu. Based on the significance of abundance differences and detection rates, we apply Gaussian mixture models to various element abundance combinations. We find the most populated and least contaminated component, which we confirm to represent GSE, contains 1049 stars selected via [Na/Fe] vs. [Mg/Mn] in GALAH+ DR3. We provide tables of our selections and report the chrono-chemodynamical properties (age, chemistry, and dynamics). Through a previously reported clean dynamical selection of GSE stars, including $30 < \sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} < 55$, we can characterise an unprecedented 24 abundances of this structure with GALAH+ DR3. Our chemical selection allows us to prevent circular reasoning and characterise the dynamical properties of the GSE, for example mean $\sqrt{J_R~/~\mathrm{kpc\,km\,s^{-1}}} = 26_{-14}^{+9}$. We find only $(29\pm1)\%$ of the GSE stars within the clean dynamical selection region. Our methodology will improve future studies of accreted structures and their importance for the formation of the Milky Way. △ Less

Submitted 5 January, 2022; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: 29 pages, 18 figures, 6 tables; publicshed in MNRAS; all code to recreate results (based on public data described in arXiv:2011.02505) and figures available at https://github.com/svenbuder/Accreted-stars-in-GALAH-DR3

Journal ref: MNRAS, 510, 2407 (2022)

Signature of a massive rotating metal-poor star imprinted in the Phoenix stellar stream

Authors: Andrew R. Casey, Alexander P. Ji, Terese T. Hansen, Ting S. Li, Sergey E. Koposov, Gary S. Da Costa, Joss Bland-Hawthorn, Lara Cullinane, Denis Erkal, Geraint F. Lewis, Kyler Kuehn, Dougal Mackey, Sarah L. Martell, Andrew B. Pace, Jeffrey D. Simpson, Daniel B. Zucker

Abstract: The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low mass globular cluster. In this work we use Magellan/MIKE high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($σしぐま(\mathrm{[Fe~II/H]}) = 0.04^{+0.11}_{-0.03}$) and… ▽ More The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low mass globular cluster. In this work we use Magellan/MIKE high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ($σしぐま(\mathrm{[Fe~II/H]}) = 0.04^{+0.11}_{-0.03}$) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr II/Fe] spanning $\sim$1 dex, while [Ba II/Fe] shows nearly no intrinsic spread ($σしぐま(\mathrm{[Ba~II/H]}) = {0.03}^{+0.10}_{-0.02}$). This abundance signature is best interpreted as slow neutron capture element production from a massive fast-rotating metal-poor star ($15-20 \mathrm{M}_\odot$, $v_\mathrm{ini}/v_\mathrm{crit} = 0.4$, $[\mathrm{Fe/H}] = -3.8$). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Neutron capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant ($A(\mathrm{Li}) = 3.1 \pm 0.1$). At $[\mathrm{Fe/H}] = -2.93$ it is among the most metal-poor lithium-rich giants known. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: Accepted to ApJ 2021-07-09

Unfolding globally resonant homoclinic tangencies

Authors: Sishu Shankar Muni, Robert I. McLachlan, David J. W. Simpson

Abstract: Global resonance is a mechanism by which a homoclinic tangency of a smooth map can have infinitely many asymptotically stable, single-round periodic solutions. To understand the bifurcation structure one would expect to see near such a tangency, in this paper we study one-parameter perturbations of typical globally resonant homoclinic tangencies. We assume the tangencies are formed by the stable a… ▽ More Global resonance is a mechanism by which a homoclinic tangency of a smooth map can have infinitely many asymptotically stable, single-round periodic solutions. To understand the bifurcation structure one would expect to see near such a tangency, in this paper we study one-parameter perturbations of typical globally resonant homoclinic tangencies. We assume the tangencies are formed by the stable and unstable manifolds of saddle fixed points of two-dimensional maps. We show the perturbations display two infinite sequences of bifurcations, one saddle-node the other period-doubling, between which single-round periodic solutions are asymptotically stable. Generically these scale like $|λらむだ|^{2 k}$, as $k \to \infty$, where $-1 < λらむだ< 1$ is the stable eigenvalue associated with the fixed point. If the perturbation is taken tangent to the surface of codimension-one homoclinic tangencies, they instead scale like $\frac{|λらむだ|^k}{k}$. We also show slower scaling laws are possible if the perturbation admits further degeneracies. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: 19 pages, 7 figures

Chaos in the border-collision normal form: A computer-assisted proof using induced maps and invariant expanding cones

Authors: P. A. Glendinning, D. J. W. Simpson

Abstract: In some maps the existence of an attractor with a positive Lyapunov exponent can be proved by constructing a trapping region in phase space and an invariant expanding cone in tangent space. If this approach fails it may be possible to adapt the strategy by considering an induced map (a first return map for a well-chosen subset of phase space). In this paper we show that such a construction can be… ▽ More In some maps the existence of an attractor with a positive Lyapunov exponent can be proved by constructing a trapping region in phase space and an invariant expanding cone in tangent space. If this approach fails it may be possible to adapt the strategy by considering an induced map (a first return map for a well-chosen subset of phase space). In this paper we show that such a construction can be applied to the two-dimensional border-collision normal form (a continuous piecewise-linear map) if a certain set of conditions are satisfied and develop an algorithm for checking these conditions. The algorithm requires relatively few computations, so it is a more efficient method than, for example, estimating the Lyapunov exponent from a single orbit in terms of speed, numerical accuracy, and rigor. The algorithm is used to prove the existence of an attractor with a positive Lyapunov exponent numerically in an area of parameter space where the map has strong rotational characteristics and the consideration of an induced map is critical for the proof of robust chaos. △ Less

Submitted 12 August, 2021; originally announced August 2021.

MSC Class: 37G35; 39A28

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

Can non-specialists provide high quality gold standard labels in challenging modalities?

Authors: Samuel Budd, Thomas Day, John Simpson, Karen Lloyd, Jacqueline Matthew, Emily Skelton, Reza Razavi, Bernhard Kainz

Abstract: Probably yes. -- Supervised Deep Learning dominates performance scores for many computer vision tasks and defines the state-of-the-art. However, medical image analysis lags behind natural image applications. One of the many reasons is the lack of well annotated medical image data available to researchers. One of the first things researchers are told is that we require significant expertise to reli… ▽ More Probably yes. -- Supervised Deep Learning dominates performance scores for many computer vision tasks and defines the state-of-the-art. However, medical image analysis lags behind natural image applications. One of the many reasons is the lack of well annotated medical image data available to researchers. One of the first things researchers are told is that we require significant expertise to reliably and accurately interpret and label such data. We see significant inter- and intra-observer variability between expert annotations of medical images. Still, it is a widely held assumption that novice annotators are unable to provide useful annotations for use by clinical Deep Learning models. In this work we challenge this assumption and examine the implications of using a minimally trained novice labelling workforce to acquire annotations for a complex medical image dataset. We study the time and cost implications of using novice annotators, the raw performance of novice annotators compared to gold-standard expert annotators, and the downstream effects on a trained Deep Learning segmentation model's performance for detecting a specific congenital heart disease (hypoplastic left heart syndrome) in fetal ultrasound imaging. △ Less

Submitted 30 July, 2021; originally announced July 2021.

Comments: Accepted at the FAIR workshop in conjunction with MICCAI'21

Measuring the Mass of the Large Magellanic Cloud with Stellar Streams Observed by ${S}^5$

Authors: Nora Shipp, Denis Erkal, Alex Drlica-Wagner, Ting S. Li, Andrew B. Pace, Sergey E. Koposov, Lara R. Cullinane, Gary S. Da Costa, Alexander P. Ji, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Jeffrey D. Simpson, Zhen Wan, Daniel B. Zucker, Joss Bland-Hawthorn, Peter S. Ferguson, Sophia Lilleengen, S5 Collaboration

Abstract: Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (${S}^5$), Gaia EDR3, and the Dark Energy Survey (DES), to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble o… ▽ More Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (${S}^5$), Gaia EDR3, and the Dark Energy Survey (DES), to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from 14 to $19 \times 10^{10}\ \mathrm{M}_{\odot}$, probed over a range of closest approach times and distances. With the most constraining stream (Orphan-Chenab), we measure an LMC mass of $18.8^{+ 3.5}_{- 4.0} \times 10^{10}\ \mathrm{M}_{\odot}$, probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC's influence on structures in the Milky Way. Using this sample of streams, we find that the LMC's effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on the impulse approximation and we show that the LMC's effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: 23 pages, 10 figures, submitted to ApJ

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

Detecting Hypo-plastic Left Heart Syndrome in Fetal Ultrasound via Disease-specific Atlas Maps

Authors: Samuel Budd, Matthew Sinclair, Thomas Day, Athanasios Vlontzos, Jeremy Tan, Tianrui Liu, Jaqueline Matthew, Emily Skelton, John Simpson, Reza Razavi, Ben Glocker, Daniel Rueckert, Emma C. Robinson, Bernhard Kainz

Abstract: Fetal ultrasound screening during pregnancy plays a vital role in the early detection of fetal malformations which have potential long-term health impacts. The level of skill required to diagnose such malformations from live ultrasound during examination is high and resources for screening are often limited. We present an interpretable, atlas-learning segmentation method for automatic diagnosis of… ▽ More Fetal ultrasound screening during pregnancy plays a vital role in the early detection of fetal malformations which have potential long-term health impacts. The level of skill required to diagnose such malformations from live ultrasound during examination is high and resources for screening are often limited. We present an interpretable, atlas-learning segmentation method for automatic diagnosis of Hypo-plastic Left Heart Syndrome (HLHS) from a single `4 Chamber Heart' view image. We propose to extend the recently introduced Image-and-Spatial Transformer Networks (Atlas-ISTN) into a framework that enables sensitising atlas generation to disease. In this framework we can jointly learn image segmentation, registration, atlas construction and disease prediction while providing a maximum level of clinical interpretability compared to direct image classification methods. As a result our segmentation allows diagnoses competitive with expert-derived manual diagnosis and yields an AUC-ROC of 0.978 (1043 cases for training, 260 for validation and 325 for testing). △ Less

Submitted 6 July, 2021; originally announced July 2021.

Comments: MICCAI'21 Main Conference

Detecting Outliers with Poisson Image Interpolation

Authors: Jeremy Tan, Benjamin Hou, Thomas Day, John Simpson, Daniel Rueckert, Bernhard Kainz

Abstract: Supervised learning of every possible pathology is unrealistic for many primary care applications like health screening. Image anomaly detection methods that learn normal appearance from only healthy data have shown promising results recently. We propose an alternative to image reconstruction-based and image embedding-based methods and propose a new self-supervised method to tackle pathological an… ▽ More Supervised learning of every possible pathology is unrealistic for many primary care applications like health screening. Image anomaly detection methods that learn normal appearance from only healthy data have shown promising results recently. We propose an alternative to image reconstruction-based and image embedding-based methods and propose a new self-supervised method to tackle pathological anomaly detection. Our approach originates in the foreign patch interpolation (FPI) strategy that has shown superior performance on brain MRI and abdominal CT data. We propose to use a better patch interpolation strategy, Poisson image interpolation (PII), which makes our method suitable for applications in challenging data regimes. PII outperforms state-of-the-art methods by a good margin when tested on surrogate tasks like identifying common lung anomalies in chest X-rays or hypo-plastic left heart syndrome in prenatal, fetal cardiac ultrasound images. Code available at https://github.com/jemtan/PII. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Comments: MICCAI 2021

Kinematics of Antlia 2 and Crater 2 from The Southern Stellar Stream Spectroscopic Survey (S5)

Authors: Alexander P. Ji, Sergey E. Koposov, Ting S. Li, Denis Erkal, Andrew B. Pace, Joshua D. Simon, Vasily Belokurov, Lara R. Cullinane, Gary S. Da Costa, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Nora Shipp, Jeffrey D. Simpson, Daniel B. Zucker, Terese T. Hansen, Joss Bland-Hawthorn, S5 Collaboration

Abstract: We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia… ▽ More We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia EDR3 proper motions, detects a clear velocity gradient in Antlia 2 and a tentative velocity gradient in Crater 2. The velocity gradient magnitudes and directions are consistent with particle stream simulations of tidal disruption. Furthermore, the orbit and kinematics of Antlia 2 require a model that includes the reflex motion of the Milky Way induced by the Large Magellanic Cloud. We also find that Antlia 2's metallicity was previously overestimated, so it lies on the empirical luminosity-metallicity relation and is likely only now experiencing substantial stellar mass loss. Current dynamical models of Antlia 2 require it to have lost over 90% of its stars to tides, in tension with the low stellar mass loss implied by the updated metallicity. Overall, the new kinematic measurements support a tidal disruption scenario for the origin of these large and extended dwarf spheroidal galaxies. △ Less

Submitted 21 September, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: 23 pages (+appendix), 11 figures, accepted to ApJ. Updated to accepted version (1 new figure, minor changes to text) Video of Antlia 2 tidal disruption here: https://www.youtube.com/watch?v=wvt-Q4kRq90

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

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

The GALAH Survey: No chemical evidence of an extragalactic origin for the Nyx stream

Authors: Daniel B. Zucker, Jeffrey D. Simpson, Sarah L. Martell, Geraint F. Lewis, Andrew R. Casey, Yuan-Sen Ting, Jonathan Horner, Thomas Nordlander, Rosemary F. G. Wyse, Tomaz Zwitter, Joss Bland-Hawthorn, Sven Buder, Martin Asplund, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael R. Hayden, Janez Kos, Jane Lin, Karin Lind, Katharine J. Schlesinger, Sanjib Sharma, Dennis Stello

Abstract: The results from the ESA Gaia astrometric mission and deep photometric surveys have revolutionized our knowledge of the Milky Way. There are many ongoing efforts to search these data for stellar substructure to find evidence of individual accretion events that built up the Milky Way and its halo. One of these newly identified features, called Nyx, was announced as an accreted stellar stream travel… ▽ More The results from the ESA Gaia astrometric mission and deep photometric surveys have revolutionized our knowledge of the Milky Way. There are many ongoing efforts to search these data for stellar substructure to find evidence of individual accretion events that built up the Milky Way and its halo. One of these newly identified features, called Nyx, was announced as an accreted stellar stream traveling in the plane of the disk. Using a combination of elemental abundances and stellar parameters from the GALAH and APOGEE surveys, we find that the abundances of the highest likelihood Nyx members are entirely consistent with membership of the thick disk, and inconsistent with a dwarf galaxy origin. We conclude that the postulated Nyx stream is most probably a high-velocity component of the Milky Way's thick disk. With the growing availability of large data sets including kinematics, stellar parameters, and detailed abundances, the probability of detecting chance associations increases, and hence new searches for substructure require confirmation across as many data dimensions as possible. △ Less

Submitted 17 April, 2021; originally announced April 2021.

Comments: Accepted for publication by ApJ Letters

Fluorescence-Detected Mid-Infrared Photothermal Microscopy

Authors: Minghe Li, Aleksandr Razumtcev, Ruochen Yang, Youlin Liu, Jiayue Rong, Andreas C. Geiger, Romain Blanchard, Cristian Pfluegl, Lynne S. Taylor, Garth J. Simpson

Abstract: We demonstrate instrumentation and methods to enable fluorescence-detected photothermal infrared (F-PTIR) microscopy, then demonstrate the utility of F-PTIR to characterize the composition within phase-separated domains of model amorphous solid dispersions (ASDs) induced by water sorption. In F-PTIR, temperature-dependent changes in fluorescence quantum efficiency are shown to sensitively report o… ▽ More We demonstrate instrumentation and methods to enable fluorescence-detected photothermal infrared (F-PTIR) microscopy, then demonstrate the utility of F-PTIR to characterize the composition within phase-separated domains of model amorphous solid dispersions (ASDs) induced by water sorption. In F-PTIR, temperature-dependent changes in fluorescence quantum efficiency are shown to sensitively report on highly localized absorption of mid-infrared radiation. The spatial resolution with which infrared spectroscopy can be performed is dictated by fluorescence microscopy, rather than the infrared wavelength. Following proof of concept F-PTIR demonstration on model systems of polyethylene glycol (PEG) and silica gel, F-PTIR enabled the characterization of chemical composition within inhomogeneous ritonavir / polyvinylpyrrolidone-vinyl acetate (PVPVA) amorphous dispersions. Phase separation is implicated in the observation of critical behaviors in ASD dissolution kinetics, with the results of F-PTIR supporting the formation of phase-separated drug-rich domains upon water absorption in spin-cast films. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 5 pages, 4 figures

The GALAH Survey and Symbiotic Stars. I. Discovery and follow-up of 33 candidate accreting-only systems

Authors: U. Munari, G. Traven, N. Masetti, P. Valisa, G. -L. Righetti, F. -J. Hambsch, A. Frigo, K. Cotar, G. M. De Silva, K. C. Freeman, G. F. Lewis, S. L. Martell, S. Sharma, J. D. Simpson, Y. -S. Ting, R. A. Wittenmyer, D. B. Zucker

Abstract: We have identified a first group of 33 new candidates for symbiotic stars (SySt) of the accreting-only variety among the 600,255 stars so far observed by the GALAH high-resolution spectroscopic survey of the Southern Hemisphere, more than doubling the number of those previously known. GALAH aims to high latitudes and this offers the possibility to sound the Galaxy for new SySt away from the usual… ▽ More We have identified a first group of 33 new candidates for symbiotic stars (SySt) of the accreting-only variety among the 600,255 stars so far observed by the GALAH high-resolution spectroscopic survey of the Southern Hemisphere, more than doubling the number of those previously known. GALAH aims to high latitudes and this offers the possibility to sound the Galaxy for new SySt away from the usual Plane and Bulge hunting regions. In this paper we focus on SySt of the M spectral type, showing an Halpha emission with a peak in excess of 0.5 above the adjacent continuum level, and not affected by coherent radial pulsations. These constraints will be relaxed in future studies. The 33 new candidate SySt were subjected to a vast array of follow-up confirmatory observations (X-ray/UV observations with the Swift satellite, search for optical flickering, presence of a near-UV upturn in ground-based photometric and spectroscopic data, radial velocity changes suggestive of orbital motion, variability of the emission line profiles). According to Gaia eDR3 parallaxes, the new SySt are located at the tip of the Giant Branch, sharing the same distribution in M(Ks) of the well established SySt. The accretion luminosities of the new SySt are in the range 1-10 Lsun, corresponding to mass-accretion rates of 0.1-1x10(-9) Msun/yr for WDs of 1 Msun. The M giant of one of the new SySt presents a large Lithium over-abundance. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: MNRAS, revised version

Leveraging pre-trained representations to improve access to untranscribed speech from endangered languages

Authors: Nay San, Martijn Bartelds, Mitchell Browne, Lily Clifford, Fiona Gibson, John Mansfield, David Nash, Jane Simpson, Myfany Turpin, Maria Vollmer, Sasha Wilmoth, Dan Jurafsky

Abstract: Pre-trained speech representations like wav2vec 2.0 are a powerful tool for automatic speech recognition (ASR). Yet many endangered languages lack sufficient data for pre-training such models, or are predominantly oral vernaculars without a standardised writing system, precluding fine-tuning. Query-by-example spoken term detection (QbE-STD) offers an alternative for iteratively indexing untranscri… ▽ More Pre-trained speech representations like wav2vec 2.0 are a powerful tool for automatic speech recognition (ASR). Yet many endangered languages lack sufficient data for pre-training such models, or are predominantly oral vernaculars without a standardised writing system, precluding fine-tuning. Query-by-example spoken term detection (QbE-STD) offers an alternative for iteratively indexing untranscribed speech corpora by locating spoken query terms. Using data from 7 Australian Aboriginal languages and a regional variety of Dutch, all of which are endangered or vulnerable, we show that QbE-STD can be improved by leveraging representations developed for ASR (wav2vec 2.0: the English monolingual model and XLSR53 multilingual model). Surprisingly, the English model outperformed the multilingual model on 4 Australian language datasets, raising questions around how to optimally leverage self-supervised speech representations for QbE-STD. Nevertheless, we find that wav2vec 2.0 representations (either English or XLSR53) offer large improvements (56-86% relative) over state-of-the-art approaches on our endangered language datasets. △ Less

Submitted 13 September, 2021; v1 submitted 26 March, 2021; originally announced March 2021.

Comments: Accepted at ASRU 2021

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

Cosmic evolution of the H2 mass density and the epoch of molecular gas

Authors: T. K. Garratt, K. E. K. Coppin, J. E. Geach, O. Almaini, W. G. Hartley, D. T. Maltby, C. J. Simpson, A. Wilkinson, C. J. Conselice, M. Franco, R. J. Ivison, M. P. Koprowski, C. C. Lovell, A. Pope, D. Scott, P. van der Werf

Abstract: We present new empirical constraints on the evolution of $ρろー_{\rm H_2}$, the cosmological mass density of molecular hydrogen, back to $z\approx2.5$. We employ a statistical approach measuring the average observed $850μみゅー{\rm m}$ flux density of near-infrared selected galaxies as a function of redshift. The redshift range considered corresponds to a span where the $850μみゅー{\rm m}$ band probes the Rayleig… ▽ More We present new empirical constraints on the evolution of $ρろー_{\rm H_2}$, the cosmological mass density of molecular hydrogen, back to $z\approx2.5$. We employ a statistical approach measuring the average observed $850μみゅー{\rm m}$ flux density of near-infrared selected galaxies as a function of redshift. The redshift range considered corresponds to a span where the $850μみゅー{\rm m}$ band probes the Rayleigh-Jeans tail of thermal dust emission in the rest-frame, and can therefore be used as an estimate of the mass of the interstellar medium (ISM). Our sample comprises of ${\approx}150,000$ galaxies in the UKIDSS-UDS field with near-infrared magnitudes $K_{\rm AB}\leq25$ mag and photometric redshifts with corresponding probability distribution functions derived from deep 12-band photometry. With a sample approximately 2 orders of magnitude larger than in previous works we significantly reduce statistical uncertainties on $ρろー_{\rm H_2}$ to $z\approx2.5$. Our measurements are in broad agreement with recent direct estimates from blank field molecular gas surveys, finding that the epoch of molecular gas coincides with the peak epoch of star formation with $ρろー_{\rm H_2}\approx2\times10^7\,{\rm M_\odot}\,{\rm Mpc^{-3}}$ at $z\approx2$. We demonstrate that $ρろー_{\rm H_2}$ can be broadly modelled by inverting the star-formation rate density with a fixed or weakly evolving star-formation efficiency. This 'constant efficiency' model shows a similar evolution to our statistically derived $ρろー_{\rm H_2}$, indicating that the dominant factor driving the peak star formation history at $z\approx2$ is a larger supply of molecular gas in galaxies rather than a significant evolution of the star-formation rate efficiency within individual galaxies. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: Accepted for publication in 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.

Hunting for Planetary Nebulae toward the Galactic Center

Authors: Jihye Hong, Janet P. Simpson, Deokkeun An, Angela S. Cotera, Solange V. Ramírez

Abstract: We present near-infrared (IR) spectra of two planetary nebula (PN) candidates in close lines of sight toward the Galactic center (GC) using the Gemini Near-Infrared Spectrograph (GNIRS) at Gemini North. High-resolution images from radio continuum and narrow-band IR observations reveal ringlike or barrel-shaped morphologies of these objects, and their mid-IR spectra from the Spitzer Space Telescope… ▽ More We present near-infrared (IR) spectra of two planetary nebula (PN) candidates in close lines of sight toward the Galactic center (GC) using the Gemini Near-Infrared Spectrograph (GNIRS) at Gemini North. High-resolution images from radio continuum and narrow-band IR observations reveal ringlike or barrel-shaped morphologies of these objects, and their mid-IR spectra from the Spitzer Space Telescope exhibit rich emission lines from highly-excited species such as [S IV], [Ne III], [Ne V], and [O IV]. We also derive elemental abundances using the Cloudy synthetic models, and find an excess amount of the $s$-process element Krypton in both targets, which supports their nature as PN. We estimate foreground extinction toward each object using near-IR hydrogen recombination lines, and find significant visual extinctions ($A_V > 20$). The distances inferred from the size versus surface brightness relation of other PNe are $9.0\pm1.6$ kpc and $7.6\pm1.6$ kpc for SSTGC 580183 and SSTGC 588220, respectively. These observed properties along with abundance patterns and their close proximity to Sgr A$^*$ (projected distances $<20$ pc) make it highly probable that these objects are the first confirmed PN objects in the nuclear stellar disk. The apparent scarcity of such objects resembles the extremely low rate of PN formation in old stellar systems, but is in line with the current rate of the sustained star formation activity in the Central Molecular Zone. △ Less

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

Comments: 20 pages, 9 figures, accepted for publication in the Astronomical Journal

Sagittarius B1 -- A Patchwork of H II Regions and PhotoDissociation Regions

Authors: Janet P. Simpson, Sean W. J. Colgan, Angela S. Cotera, Michael J. Kaufman, Susan R. Stolovy

Abstract: Sgr B1 is a luminous H II region in the Galactic Center immediately next to the massive star-forming giant molecular cloud Sgr B2 and apparently connected to it from their similar radial velocities. In 2018 we showed from SOFIA FIFI-LS observations of the [O III] 52 and 88 micron lines that there is no central exciting star cluster and that the ionizing stars must be widely spread throughout the r… ▽ More Sgr B1 is a luminous H II region in the Galactic Center immediately next to the massive star-forming giant molecular cloud Sgr B2 and apparently connected to it from their similar radial velocities. In 2018 we showed from SOFIA FIFI-LS observations of the [O III] 52 and 88 micron lines that there is no central exciting star cluster and that the ionizing stars must be widely spread throughout the region. Here we present SOFIA FIFI-LS observations of the [O I] 146 and [C II] 158 micron lines formed in the surrounding photodissociation regions (PDRs). We find that these lines correlate neither with each other nor with the [O III] lines although together they correlate better with the 70 micron Herschel PACS images from Hi-GAL. We infer from this that Sgr B1 consists of a number of smaller H II regions plus their associated PDRs, some seen face-on and the others seen more or less edge-on. We used the PDR Toolbox to estimate densities and the far-ultraviolet intensities exciting the PDRs. Using models computed with Cloudy, we demonstrate possible appearances of edge-on PDRs and show that the density difference between the PDR densities and the electron densities estimated from the [O III] line ratios is incompatible with pressure equilibrium unless there is a substantial pressure contribution from either turbulence or magnetic field or both. We likewise conclude that the hot stars exciting Sgr B1 are widely spaced throughout the region at substantial distances from the gas with no evidence of current massive star formation. △ Less

Submitted 13 February, 2021; originally announced February 2021.

Comments: 18 pages, 8 figures. Accepted by The Astrophysical Journal

On the stability of boundary equilibria in Filippov systems

Authors: David J. W. Simpson

Abstract: The leading-order approximation to a Filippov system $f$ about a generic boundary equilibrium $x^*$ is a system $F$ that is affine one side of the boundary and constant on the other side. We prove $x^*$ is exponentially stable for $f$ if and only if it is exponentially stable for $F$ when the constant component of $F$ is not tangent to the boundary. We then show exponential stability and asymptoti… ▽ More The leading-order approximation to a Filippov system $f$ about a generic boundary equilibrium $x^*$ is a system $F$ that is affine one side of the boundary and constant on the other side. We prove $x^*$ is exponentially stable for $f$ if and only if it is exponentially stable for $F$ when the constant component of $F$ is not tangent to the boundary. We then show exponential stability and asymptotic stability are in fact equivalent for $F$. We also show exponential stability is preserved under small perturbations to the pieces of $F$. Such results are well known for homogeneous systems. To prove the results here additional techniques are required because the two components of $F$ have different degrees of homogeneity. The primary function of the results is to reduce the problem of the stability of $x^*$ from the general Filippov system $f$ to the simpler system $F$. Yet in general this problem remains difficult. We provide a four-dimensional example of $F$ for which orbits appear to converge to $x^*$ in a chaotic fashion. By utilising the presence of both homogeneity and sliding motion the dynamics of $F$ can in this case be reduced to the combination of a one-dimensional return map and a scalar function. △ Less

Submitted 11 January, 2021; originally announced January 2021.

MSC Class: 34A36; 37G10

Mean exit time for diffusion on irregular domains

Authors: Matthew J Simpson, Daniel J Vandenheuvel, Joshua M Wilson, Scott W McCue, Elliot J Carr

Abstract: Many problems in physics, biology, and economics depend upon the duration of time required for a diffusing particle to cross a boundary. As such, calculations of the distribution of first passage time, and in particular the mean first passage time, is an active area of research relevant to many disciplines. Exact results for the mean first passage time for diffusion on simple geometries, such as l… ▽ More Many problems in physics, biology, and economics depend upon the duration of time required for a diffusing particle to cross a boundary. As such, calculations of the distribution of first passage time, and in particular the mean first passage time, is an active area of research relevant to many disciplines. Exact results for the mean first passage time for diffusion on simple geometries, such as lines, discs and spheres, are well--known. In contrast, computational methods are often used to study the first passage time for diffusion on more realistic geometries where closed--form solutions of the governing elliptic boundary value problem are not available. Here, we develop a perturbation solution to calculate the mean first passage time on irregular domains formed by perturbing the boundary of a disc or an ellipse. Classical perturbation expansion solutions are then constructed using the exact solutions available on a disc and an ellipse. We apply the perturbation solutions to compute the mean first exit time on two naturally--occurring irregular domains: a map of Tasmania, an island state of Australia, and a map of Taiwan. Comparing the perturbation solutions with numerical solutions of the elliptic boundary value problem on these irregular domains confirms that we obtain a very accurate solution with a few terms in the series only. Matlab software to implement all calculations is available on GitHub. △ Less

Submitted 5 February, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: 31pages, 12 figures

MSC Class: 92Bxx; 92-08

Extinction of bistable populations is affected by the shape of their initial spatial distribution

Authors: Yifei Li, Stuart T. Johnston, Pascal R. Buenzli, Peter van Heijster, Matthew J. Simpson

Abstract: The question of whether biological populations survive or are eventually driven to extinction has long been examined using mathematical models. In this work we study population survival or extinction using a stochastic, discrete lattice-based random walk model where individuals undergo movement, birth and death events. The discrete model is defined on a two-dimensional hexagonal lattice with perio… ▽ More The question of whether biological populations survive or are eventually driven to extinction has long been examined using mathematical models. In this work we study population survival or extinction using a stochastic, discrete lattice-based random walk model where individuals undergo movement, birth and death events. The discrete model is defined on a two-dimensional hexagonal lattice with periodic boundary conditions. A key feature of the discrete model is that crowding effects are introduced by specifying two different crowding functions that govern how local agent density influences movement events and birth/death events. The continuum limit description of the discrete model is a nonlinear reaction-diffusion equation, and we focus on crowding functions that lead to linear diffusion and a bistable source term that is often associated with the strong Allee effect. Using both the discrete and continuum modelling tools we explore the complicated relationship between the long-term survival or extinction of the population and the initial spatial arrangement of the population. In particular, we study different spatial arrangements of initial distributions: (i) a well-mixed initial distribution where the initial density is independent of position in the domain; (ii) a vertical strip initial distribution where the initial density is independent of vertical position in the domain; and, (iii) several forms of two-dimensional initial distributions where the initial population is distributed in regions with different shapes. Our results indicate that the shape of the initial spatial distribution of the population affects extinction of bistable populations. All software required to solve the discrete and continuum models used in this work are available on GitHub. △ Less

Submitted 14 September, 2021; v1 submitted 5 January, 2021; originally announced January 2021.

Comments: 14 figures

The GALAH+ Survey: A New Library of Observed Stellar Spectra Improves Radial Velocities and Hints at Motions within M67

Authors: Tomaž Zwitter, Janez Kos, Sven Buder, Klemen Čotar, Martin Asplund, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Kenneth C. Freeman, Michael R. Hayden, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Kevin L. Beeson, Richard de Grijs, Thomas Nordlander, Yuan-Sen Ting, Gregor Traven , et al. (3 additional authors not shown)

Abstract: GALAH+ is a magnitude-limited survey of high resolution stellar spectra obtained by the HERMES spectrograph at the Australian Astronomical Observatory. Its third data release provides reduced spectra with new derivations of stellar parameters and abundances of 30 chemical elements for 584,015 dwarfs and giants, 88% of them in the Gaia magnitude range 11 < G < 14. Here we use these improved values… ▽ More GALAH+ is a magnitude-limited survey of high resolution stellar spectra obtained by the HERMES spectrograph at the Australian Astronomical Observatory. Its third data release provides reduced spectra with new derivations of stellar parameters and abundances of 30 chemical elements for 584,015 dwarfs and giants, 88% of them in the Gaia magnitude range 11 < G < 14. Here we use these improved values of stellar parameters to build a library of observed spectra which is useful to study variations of individual spectral lines with stellar parameters. This and other improvements are used to derive radial velocities with uncertainties which are generally within 0.1 km/s or ~25% smaller than in the previous release. Median differences in radial velocities measured here and by the Gaia DR2 or APOGEE DR16 surveys are smaller than 30 m/s, a larger offset is present only for Gaia measurements of giant stars. We identify 4483 stars with intrinsically variable velocities and 225 stars for which the velocity stays constant over >=3 visits spanning more than a year. The combination of radial velocities from GALAH+ with distances and sky plane motions from Gaia enables studies of dynamics within streams and clusters. For example, we estimate that the open cluster M67 has a total mass of ~3300 Msun and its outer parts seem to be expanding, though astrometry with a larger time-span than currently available from Gaia eDR3 is needed to judge if the latter result is real. △ Less

Submitted 14 September, 2021; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 15 pages, 4 tables, 17 figures + online data files, accepted for publication in MNRAS

Journal ref: 2021, MNRAS, 508, 4202-4215

Learning normal appearance for fetal anomaly screening: Application to the unsupervised detection of Hypoplastic Left Heart Syndrome

Authors: Elisa Chotzoglou, Thomas Day, Jeremy Tan, Jacqueline Matthew, David Lloyd, Reza Razavi, John Simpson, Bernhard Kainz

Abstract: Congenital heart disease is considered as one the most common groups of congenital malformations which affects $6-11$ per $1000$ newborns. In this work, an automated framework for detection of cardiac anomalies during ultrasound screening is proposed and evaluated on the example of Hypoplastic Left Heart Syndrome (HLHS), a sub-category of congenital heart disease. We propose an unsupervised approa… ▽ More Congenital heart disease is considered as one the most common groups of congenital malformations which affects $6-11$ per $1000$ newborns. In this work, an automated framework for detection of cardiac anomalies during ultrasound screening is proposed and evaluated on the example of Hypoplastic Left Heart Syndrome (HLHS), a sub-category of congenital heart disease. We propose an unsupervised approach that learns healthy anatomy exclusively from clinically confirmed normal control patients. We evaluate a number of known anomaly detection frameworks together with a model architecture based on the $αあるふぁ$-GAN network and find evidence that the proposed model performs significantly better than the state-of-the-art in image-based anomaly detection, yielding average $0.81$ AUC \emph{and} a better robustness towards initialisation compared to previous works. △ Less

Submitted 9 September, 2021; v1 submitted 15 November, 2020; originally announced December 2020.

The GALAH Survey: Dependence of elemental abundances on age and metallicity for stars in the Galactic disc

Authors: Sanjib Sharma, Michael R. Hayden, Joss Bland-Hawthorn, Dennis Stello, Sven Buder, Joel C. Zinn, Lorenzo Spina, Thomas Kallinger, Martin Asplund, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Klemen Cotar, Boquan Chen, Prajwal R. Kafle, Shourya Khanna , et al. (2 additional authors not shown)

Abstract: Using data from the GALAH survey, we explore the dependence of elemental abundances on stellar age and metallicity among Galactic disc stars. We find that the abundance of most elements can be predicted from age and [Fe/H] with an intrinsic scatter of about 0.03 dex. We discuss the possible causes for the existence of the abundance-age-metallicity relations. Using a stochastic chemical enrichment… ▽ More Using data from the GALAH survey, we explore the dependence of elemental abundances on stellar age and metallicity among Galactic disc stars. We find that the abundance of most elements can be predicted from age and [Fe/H] with an intrinsic scatter of about 0.03 dex. We discuss the possible causes for the existence of the abundance-age-metallicity relations. Using a stochastic chemical enrichment scheme based on the size of Supernovae remnants, we show the intrinsic scatter is expected to be small, about 0.05 dex or even smaller if there is additional mixing in the ISM. Elemental abundances show trends with both age and metallicity and the relationship is well described by a simple model in which the dependence of abundance ([X/Fe]) on age and [Fe/H] are additively separable. Elements can be grouped based on the direction of their abundance gradient in the (age,[Fe/H]) plane and different groups can be roughly associated with three distinct nucleosynthetic production sites, the exploding massive stars, the exploding white dwarfs and the AGB stars. However, the abundances of some elements, like Co, La, and Li, show large scatter for a given age and metallicity, suggesting processes other than simple Galactic chemical evolution are at play. We also compare the abundance trends of main-sequence turn-off stars against that of giants, whose ages were estimated using asteroseismic information from the K2 mission. For most elements, the trends of main-sequence turn-off stars are similar to that of giants. The existence of abundance relations implies that we can estimate the age and birth radius of disc stars, which is important for studying the dynamic and chemical evolution of the Galaxy. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 21 pages, 11 Figures

The GALAH Survey: Chemical Clocks

Authors: Michael R. Hayden, Sanjib Sharma, Joss Bland-Hawthorn, Lorenzo Spina, Sven Buder, Martin Asplund, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Boquan Chen, Klemen Cotar, Diane Feuillet, Jonti Horner, Meridith Joyce, Thomas Nordlander , et al. (5 additional authors not shown)

Abstract: Previous studies have found that the elemental abundances of a star correlate directly with its age and metallicity. Using this knowledge, we derive ages for a sample of 250,000 stars taken from GALAH DR3 using only their overall metallicity and chemical abundances. Stellar ages are estimated via the machine learning algorithm $XGBoost$, using main sequence turnoff stars with precise ages as our i… ▽ More Previous studies have found that the elemental abundances of a star correlate directly with its age and metallicity. Using this knowledge, we derive ages for a sample of 250,000 stars taken from GALAH DR3 using only their overall metallicity and chemical abundances. Stellar ages are estimated via the machine learning algorithm $XGBoost$, using main sequence turnoff stars with precise ages as our input training set. We find that the stellar ages for the bulk of the GALAH DR3 sample are accurate to 1-2 Gyr using this method. With these ages, we replicate many recent results on the age-kinematic trends of the nearby disk, including the age-velocity dispersion relationship of the solar neighborhood and the larger global velocity dispersion relations of the disk found using $Gaia$ and GALAH. The fact that chemical abundances alone can be used to determine a reliable age for a star have profound implications for the future study of the Galaxy as well as upcoming spectroscopic surveys. These results show that the chemical abundance variation at a given birth radius is quite small, and imply that strong chemical tagging of stars directly to birth clusters may prove difficult with our current elemental abundance precision. Our results highlight the need of spectroscopic surveys to deliver precision abundances for as many nucleosynthetic production sites as possible in order to estimate reliable ages for stars directly from their chemical abundances. Applying the methods outlined in this paper opens a new door into studies of the kinematic structure and evolution of the disk, as ages may potentially be estimated for a large fraction of stars in existing spectroscopic surveys. This would yield a sample of millions of stars with reliable age determinations, and allow precise constraints to be put on various kinematic processes in the disk, such as the efficiency and timescales of radial migration. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 13 pages, 15 figures, submitted to MNRAS

Magnon-phonon hybridization in quasi-2D antiferromagnet MnPSe$_3$

Authors: Thuc T. Mai, K. F. Garrity, A. McCreary, J. Argo, J. R. Simpson, V. Doan-Nguyen, R. Valdes Aguilar, A. R. Hight Walker

Abstract: Magnetic excitations in van der Waals (vdW) materials, especially in the two-dimensional (2D) limit, are an exciting research topic from both the fundamental and applied perspectives. Using temperature-dependent, magneto-Raman spectroscopy, we identify the hybridization of two-magnon excitations with two separate E$_\mathrm{g}$ phonons in MnPSe$_3$, a magnetic vdW material that could potentially h… ▽ More Magnetic excitations in van der Waals (vdW) materials, especially in the two-dimensional (2D) limit, are an exciting research topic from both the fundamental and applied perspectives. Using temperature-dependent, magneto-Raman spectroscopy, we identify the hybridization of two-magnon excitations with two separate E$_\mathrm{g}$ phonons in MnPSe$_3$, a magnetic vdW material that could potentially host 2D antiferromagnetism. Results from first principles calculations of the phonon and magnon spectra further support our identification. The Raman spectra's rich temperature dependence through the magnetic transition displays an avoided-crossing behavior in the phonons' frequency and a concurrent decrease in their lifetimes. We construct a model based on the interaction between a discrete level and a continuum that reproduces these observations. The strong magnon-phonon hybridization reported here highlights the need to understand its effects on spin transport experiments in magnetic vdW materials. △ Less

Submitted 9 December, 2020; v1 submitted 25 November, 2020; originally announced November 2020.

FIR-luminous [CII] emitters in the ALMA-SCUBA-2 COSMOS survey (AS2COSMOS): The nature of submillimeter galaxies in a 10 comoving Mpc-scale structure at z~4.6

Authors: Ikki Mitsuhashi, Yuichi Matsuda, Ian Smail, Natsuki Hayatsu, James Simpson, Mark Swinbank, Hideki Umahata, Ugne Dudzevičiūtė, Jack Birkin, Soh Ikarashi, Chian-Chou Chen, Ken-ichi Tadaki, Hidenobu Yajima, Yuichi Harikane, Hanae Inami, Scott Chapman, Bunyo Hatsukade, Daisuke Iono, Andrew Bunker, Yiping Ao, Tomoki Saito, Junko Ueda, Seiichi Sakamoto

Abstract: We report the discovery of a 10 comoving Mpc-scale structure traced by massive submillimeter galaxies (SMGs) at z~4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources ($S_{850μみゅー{\rm m}}\geq$ 6.2 mJy) across 1.6 degrees$^2$ in the COSMOS field. We identify four [CII] emitting SMGs and two probable [CII] emitting SMG ca… ▽ More We report the discovery of a 10 comoving Mpc-scale structure traced by massive submillimeter galaxies (SMGs) at z~4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources ($S_{850μみゅー{\rm m}}\geq$ 6.2 mJy) across 1.6 degrees$^2$ in the COSMOS field. We identify four [CII] emitting SMGs and two probable [CII] emitting SMG candidates at z=4.60-4.64 with velocity-integrated signal-to-noise ratio of SNR>8. Four of the six emitters are near-infrared blank SMGs. After excluding one SMG whose emission line is falling at the edge of the spectral window, all galaxies show clear velocity gradients along the major axes that are consistent with rotating gas disks. The estimated rotation velocities of the disks are 330-550 km s$^{-1}$ and the inferred host dark-matter halo masses are ~2-8 $\times$ 10$^{12}$M$_{\odot}$. From their estimated halo masses and [CII] luminosity function, we suggest that these galaxies have a high (50-100%) duty cycle and high (~0.1) baryon conversion efficiency (SFR relative to baryon accretion rate), and that they contribute $\simeq$2% to the total star-formation rate density at z=4.6. These SMGs are concentrated within just 0.3% of the full survey volume, suggesting they are strongly clustered. The extent of this structure and the individual halo masses suggest that these SMGs will likely evolve into members of a ~10$^{15}$M$_{\odot}$ cluster at z=0. This survey reveals synchronized dusty starburst in massive halos at z>4, which could be driven by mergers or fed by smooth gas accretion. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Learning differential equation models from stochastic agent-based model simulations

Authors: John T. Nardini, Ruth E. Baker, Matthew J. Simpson, Kevin B. Flores

Abstract: Agent-based models provide a flexible framework that is frequently used for modelling many biological systems, including cell migration, molecular dynamics, ecology, and epidemiology. Analysis of the model dynamics can be challenging due to their inherent stochasticity and heavy computational requirements. Common approaches to the analysis of agent-based models include extensive Monte Carlo simula… ▽ More Agent-based models provide a flexible framework that is frequently used for modelling many biological systems, including cell migration, molecular dynamics, ecology, and epidemiology. Analysis of the model dynamics can be challenging due to their inherent stochasticity and heavy computational requirements. Common approaches to the analysis of agent-based models include extensive Monte Carlo simulation of the model or the derivation of coarse-grained differential equation models to predict the expected or averaged output from the agent-based model. Both of these approaches have limitations, however, as extensive computation of complex agent-based models may be infeasible, and coarse-grained differential equation models can fail to accurately describe model dynamics in certain parameter regimes. We propose that methods from the equation learning field provide a promising, novel, and unifying approach for agent-based model analysis. Equation learning is a recent field of research from data science that aims to infer differential equation models directly from data. We use this tutorial to review how methods from equation learning can be used to learn differential equation models from agent-based model simulations. We demonstrate that this framework is easy to use, requires few model simulations, and accurately predicts model dynamics in parameter regions where coarse-grained differential equation models fail to do so. We highlight these advantages through several case studies involving two agent-based models that are broadly applicable to biological phenomena: a birth-death-migration model commonly used to explore cell biology experiments and a susceptible-infected-recovered model of infectious disease spread. △ Less

Submitted 27 April, 2021; v1 submitted 16 November, 2020; originally announced November 2020.

Journal ref: Journal of the Royal Society Interface 18 (176) 2021