RESUMEN

Periodic travelling waves (PTW) are a common solution type of partial differential equations. Such models exhibit multistability of PTWs, typically visualised through the Busse balloon, and parameter changes typically lead to a cascade of wavelength changes through the Busse balloon. In the past, the stability boundaries of the Busse balloon have been used to predict such wavelength changes. Here, motivated by anecdotal evidence from previous work, we provide compelling evidence that the Busse balloon provides insufficient information to predict wavelength changes due to a delayed loss of stability phenomenon. Using two different reaction-advection-diffusion systems, we relate the delay that occurs between the crossing of a stability boundary in the Busse balloon and the occurrence of a wavelength change to features of the essential spectrum of the destabilised PTW. This leads to a predictive framework that can estimate the order of magnitude of such a time delay, which provides a novel "early warning sign" for pattern destabilization. We illustrate the implementation of the predictive framework to predict under what conditions a wavelength change of a PTW occurs.

RESUMEN

This paper concerns the invasion dynamics of the lattice pioneer-climax competition model with parameter regions in which the system is non-monotone. We estimate the spreading speeds and establish appropriate conditions under which the spreading speeds are linearly selected. Moreover, the existence of travelling waves is determined by constructing suitable upper and lower solutions. It shows that the spreading speed coincides with the minimum wave speed of travelling waves if the diffusion rate of the invasive species is larger or equal to that of the native species. Our results are new to estimate the spreading speed of non-monotone lattice pioneer-climax systems, and the techniques developed in this work can be used to study the invasion dynamics of the pioneer-climax system with interaction delays, which could extend the results in the literature. The analysis replies on the construction of auxiliary systems, upper and lower solutions, and the monotone dynamical system approach.

Asunto(s)

RESUMEN

The emergence and maintenance of tree species diversity in tropical forests is commonly attributed to the Janzen-Connell (JC) hypothesis, which states that growth of seedlings is suppressed in the proximity of conspecific adult trees. As a result, a JC distribution due to a density-dependent negative feedback emerges in the form of a (transient) pattern where conspecific seedling density is highest at intermediate distances away from parent trees. Several studies suggest that the required density-dependent feedbacks behind this pattern could result from interactions between trees and soil-borne pathogens. However, negative plant-soil feedback may involve additional mechanisms, including the accumulation of autotoxic compounds generated through tree litter decomposition. An essential task therefore consists in constructing mathematical models incorporating both effects showing the ability to support the emergence of JC distributions. In this work, we develop and analyse a novel reaction-diffusion-ODE model, describing the interactions within tropical tree species across different life stages (seeds, seedlings, and adults) as driven by negative plant-soil feedback. In particular, we show that under strong negative plant-soil feedback travelling wave solutions exist, creating transient distributions of adult trees and seedlings that are in agreement with the Janzen-Connell hypothesis. Moreover, we show that these travelling wave solutions are pulled fronts and a robust feature as they occur over a broad parameter range. Finally, we calculate their linear spreading speed and show its (in)dependence on relevant nondimensional parameters.

Asunto(s)

RESUMEN

With over 2 million people in the UK suffering from chronic wounds, understanding the biochemistry and pharmacology that underpins these wounds and wound healing is of high importance. Chronic wounds are characterised by high levels of matrix metalloproteinases (MMPs), which are necessary for the modification of healthy tissue in the healing process. Overexposure of MMPs, however, adversely affects healing of the wound by causing further destruction of the surrounding extracellular matrix. In this work, we propose a mathematical model that focuses on the interaction of MMPs with dermal cells using a system of partial differential equations. Using biologically realistic parameter values, this model gives rise to travelling waves corresponding to a front of healthy cells invading a wound. From the arising travelling wave analysis, we observe that deregulated apoptosis results in the emergence of chronic wounds, characterised by elevated MMP concentrations. We also observe hysteresis effects when both the apoptotic rate and MMP production rate are varied, providing further insight into the management (and potential reversal) of chronic wounds.

Asunto(s)

RESUMEN

We derive rigorous estimates on the speed of invasion of an advantageous trait in a spatially advancing population in the context of a system of one-dimensional F-KPP equations. The model was introduced and studied heuristically and numerically in a paper by Venegas-Ortiz et al. (Genetics 196:497-507, 2014). In that paper, it was noted that the speed of invasion by the mutant trait is faster when the resident population is expanding in space compared to the speed when the resident population is already present everywhere. We use the Feynman-Kac representation to provide rigorous estimates that confirm these predictions.

Asunto(s)

RESUMEN

The main goal of the presented study is to introduce a model of a pairwise invasion interaction with a nonlinear diffusion and advection. The new equation is based on the further general works introduced by Bramson (1988) to describe the invasive-invaded dynamics. This type of model is made particular with a density dependent diffusion along with an advection term. The new resulting model is then analyzed to explore the regularity, existence and uniqueness of solutions. It is well known that density dependent diffusion operators induce a propagating front with finite speed for compactly supported functions. Based on this, we introduce an analytical approach to determine the evolution of such a propagating front in the invasion dynamics. Afterward, we study the problem with travelling wave profiles and a numerical assessment. As a main finding to remark: When both species propagate with significantly different travelling wave speeds, the interaction becomes unstable, while when the species propagate with similar low speeds, the interaction stabilizes.

RESUMEN

We introduce a broad class of mechanistic spatial models to describe how spatially heterogeneous populations live, die, and reproduce. Individuals are represented by points of a point measure, whose birth and death rates can depend both on spatial position and local population density, defined at a location to be the convolution of the point measure with a suitable non-negative integrable kernel centred on that location. We pass to three different scaling limits: an interacting superprocess, a nonlocal partial differential equation (PDE), and a classical PDE. The classical PDE is obtained both by a two-step convergence argument, in which we first scale time and population size and pass to the nonlocal PDE, and then scale the kernel that determines local population density; and in the important special case in which the limit is a reaction-diffusion equation, directly by simultaneously scaling the kernel width, timescale and population size in our individual based model. A novelty of our model is that we explicitly model a juvenile phase. The number of juveniles produced by an individual depends on local population density at the location of the parent; these juvenile offspring are thrown off in a (possibly heterogeneous, anisotropic) Gaussian distribution around the location of the parent; they then reach (instant) maturity with a probability that can depend on the local population density at the location at which they land. Although we only record mature individuals, a trace of this two-step description remains in our population models, resulting in novel limits in which the spatial dynamics are governed by a nonlinear diffusion. Using a lookdown representation, we are able to retain information about genealogies relating individuals in our population and, in the case of deterministic limiting models, we use this to deduce the backwards in time motion of the ancestral lineage of an individual sampled from the population. We observe that knowing the history of the population density is not enough to determine the motion of ancestral lineages in our model. We also investigate (and contrast) the behaviour of lineages for three different deterministic models of a population expanding its range as a travelling wave: the Fisher-KPP equation, the Allen-Cahn equation, and a porous medium equation with logistic growth.

RESUMEN

Growing evidence suggests that travelling waves are functionally relevant for cognitive operations in the brain. Several electroencephalography (EEG) studies report on a perceptual alpha-echo, representing the brain response to a random visual flicker, propagating as a travelling wave across the cortical surface. In this study, we ask if the propagating activity of the alpha-echo is best explained by a set of discrete sources mixing at the sensor level rather than a cortical travelling wave. To this end, we presented participants with gratings modulated by random noise and simultaneously acquired the ongoing MEG. The perceptual alpha-echo was estimated using the temporal response function linking the visual input to the brain response. At the group level, we observed a spatial decay of the amplitude of the alpha-echo with respect to the sensor where the alpha-echo was the largest. Importantly, the propagation latencies consistently increased with the distance. Interestingly, the propagation of the alpha-echoes was predominantly centro-lateral, while EEG studies reported mainly posterior-frontal propagation. Moreover, the propagation speed of the alpha-echoes derived from the MEG data was around 10 m/s, which is higher compared to the 2 m/s reported in EEG studies. Using source modelling, we found an early component in the primary visual cortex and a phase-lagged late component in the parietal cortex, which may underlie the travelling alpha-echoes at the sensor level. We then simulated the alpha-echoes using realistic EEG and MEG forward models by placing two sources in the parietal and occipital cortices in accordance with our empirical findings. The two-source model could account for both the direction and speed of the observed alpha-echoes in the EEG and MEG data. Our results demonstrate that the propagation of the perceptual echoes observed in EEG and MEG data can be explained by two sources mixing at the scalp level equally well as by a cortical travelling wave. Importantly, these findings should not be directly extrapolated to intracortical recordings, where travelling waves gradually propagate at a sub-millimetre scale.

Asunto(s)

RESUMEN

We present an experimental and two-phase computational study of convection in a liquid bridge ([Formula: see text]) that develops under the action of a parallel gas flow. The study focuses on tracking the evolution of hydrothermal waves by increasing the applied temperature difference [Formula: see text] and the temperature of gas moving at the velocity [Formula: see text]. Our experiments revealed certain regularity in the change of oscillatory states with an increase in the control parameters. Above the instability threshold, the nonlinear dynamics passes through three oscillatory regimes, which are repeated in a somewhat similar way at higher values of the control parameters. They are periodic, quasi-periodic with two or three frequencies and multi-frequency state when the Fourier spectrum is filled with clusters of duplex, triplex or higher numbers of frequencies. Three-dimensional numerical simulation, complemented by a deep spectral analysis, sheds light on the evolution of the flow pattern observed in experiments. The developed methodology identified conditions for the existence of a multi-frequency regime such as the presence of a weak low-frequency mode that can modulate strong high-frequency modes, the existence of strong azimuthal modes with different wavenumbers and the [Formula: see text] mode, and the structured combination of peaks in the Fourier spectrum. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.

RESUMEN

While seizure activity may be electrographically widespread, increasing evidence has suggested that ictal discharges may in fact represent travelling waves propagated from a focal seizure source. Interictal epileptiform discharges (IEDs) are an electrographic manifestation of excessive hypersynchronization of cortical activity that occur between seizures and are considered a marker of potentially epileptogenic tissue. The precise relationship between brain regions demonstrating IEDs and those involved in seizure onset, however, remains poorly understood. Here, we hypothesize that IEDs likewise reflect the receipt of travelling waves propagated from the same regions which give rise to seizures. Forty patients from our institution who underwent invasive monitoring for epilepsy, proceeded to surgery and had at least one year of follow-up were included in our study. Interictal epileptiform discharges were detected using custom software, validated by a clinical epileptologist. We show that IEDs reach electrodes in sequences with a consistent temporal ordering, and this ordering matches the timing of receipt of ictal discharges, suggesting that both types of discharges spread as travelling waves. We use a novel approach for localization of ictal discharges, in which time differences of discharge receipt at nearby electrodes are used to compute source location; similar algorithms have been used in acoustics and geophysics. We find that interictal discharges co-localize with ictal discharges. Moreover, interictal discharges tend to localize to the resection territory in patients with good surgical outcome and outside of the resection territory in patients with poor outcome. The seizure source may originate at, and also travel to, spatially distinct IED foci. Our data provide evidence that interictal discharges may represent travelling waves of pathological activity that are similar to their ictal counterparts, and that both ictal and interictal discharges emerge from common epileptogenic brain regions. Our findings have important clinical implications, as they suggest that seizure source localizations may be derived from interictal discharges, which are much more frequent than seizures.

Asunto(s)

RESUMEN

Spatiotemporal patterns of vegetation are a ubiquitous feature of semi-arid ecosystems. On sloped terrain, vegetation patterns occur as stripes perpendicular to the contours. Field studies report contrasting long-term dynamics between different observation sites; some observe slow uphill migration of vegetation bands while some report stationary patterns. In this paper, we show that long-range seed dispersal provides a mechanism that enables the occurrence of both migrating and stationary patterns. We utilise a nonlocal PDE model in which seed dispersal is accounted for by a convolution term. The model represents vegetation patterns as periodic travelling waves and numerical continuation shows that both migrating and almost stationary patterns are stable if seed dispersal distances are sufficiently large. We use a perturbation theory approach to obtain analytical confirmation of the existence of almost stationary patterned solutions and provide a biological interpretation of the phenomenon.

Asunto(s)

RESUMEN

Cyclic predator-prey systems are often observed in nature. In a spatial setting, these can manifest as periodic traveling waves (PTW). Environmental change and direct human activity are known to, among other effects, increase the heterogeneity of the physical environment, which prey and predator inhabit. Aiming to understand the effects of heterogeneity on predator-prey PTWs, we consider a one-dimensional infinite landscape Rosenzweig-MacArthur reaction-diffusion model, with alternating patch types, and study the PTWs in this system. Applying the method of homogenisation, we show how heterogeneity can affect the stability of PTW solutions. We illustrate how the effects of heterogeneity can be understood and interpreted using Turchin's concept of residence index (encapsuling diffusion rate and patch preference). In particular, our results show that prey heterogeneity acts to modulate the effects of predator heterogeneity, by this we mean that as prey increasingly spend more time in one patch type over another the stability of the PTWs becomes more sensitive to heterogeneity in predator movement and behaviour.

Asunto(s)

RESUMEN

We analyze an 'up-the-gradient' model for the formation of transport channels of the phytohormone auxin, through auxin-mediated polarization of the PIN1 auxin transporter. We show that this model admits a family of travelling wave solutions that is parameterized by the height of the auxin-pulse. We uncover scaling relations for the speed and width of these waves and verify these rigorous results with numerical computations. In addition, we provide explicit expressions for the leading-order wave profiles, which allows the influence of the biological parameters in the problem to be readily identified. Our proofs are based on a generalization of the scaling principle developed by Friesecke and Pego to construct pulse solutions to the classic Fermi-Pasta-Ulam-Tsingou model, which describes a one-dimensional chain of coupled nonlinear springs.

Asunto(s)

RESUMEN

Spatiotemporal patterns are common in biological systems. For electrically coupled cells, previous studies of pattern formation have mainly used applied current as the primary bifurcation parameter. The purpose of this paper is to show that applied current is not needed to generate spatiotemporal patterns for smooth muscle cells. The patterns can be generated solely by external mechanical stimulation (transmural pressure). To do this we study a reaction-diffusion system involving 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 that the spatially extended model exhibits spatiotemporal chaos. We also use travelling wave coordinates to analyse travelling waves.

Asunto(s)

RESUMEN

Nonlocal reaction-diffusion equations describe various biological and biomedical applications. Their mathematical properties are essentially different in comparison with the local equations, and this difference can lead to important biological implications. This review will present the state of the art in the investigation of nonlocal reaction-diffusion models in biomedical applications. We will consider various models arising in mathematical immunology, neuroscience, cancer modelling, and we will discuss their mathematical properties, nonlinear dynamics, resulting spatiotemporal patterns and biological significance.

Asunto(s)

RESUMEN

Background: Kinetic Inductance Travelling Wave Parametric Amplifiers (KITWPAs) are a variant of superconducting amplifier that can potentially achieve high gain with quantum-limited noise performance over broad bandwidth, which is important for many ultra-sensitive experiments. In this paper, we present a novel modelling technique that can better capture the electromagnetic behaviour of a KITWPA without the translation symmetry assumption, allowing us to flexibly explore the use of more complex transmission line structures and better predict their performance. Methods: In order to design a KITWPA with optimal performance, we investigate the use of different superconducting thin film materials, and compare their pros and cons in forming a high-gain low-loss medium feasible for amplification. We establish that if the film thickness can be controlled precisely, the material used has less impact on the performance of the device, as long as it is topologically defect-free and operating within its superconducting regime. With this insight, we propose the use of Titanium Nitride (TiN) film for our KITWPA as its critical temperature can be easily altered to suit our applications. We further investigate the topological effect of different commonly used superconducting transmission line structures with the TiN film, including the effect of various non-conducting materials required to form the amplifier. Results: Both of these comprehensive studies led us to two configurations of the KITWPA: 1) A low-loss 100 nm thick TiN coplanar waveguide amplifier, and 2) A compact 50 nm TiN inverted microstrip amplifier. We utilise the novel modelling technique described in the first part of the paper to explore and investigate the optimal design and operational setup required to achieve high gain with the broadest bandwidth for both KITWPAs, including the effect of loss. Conclusions: Finally, we conclude the paper with the actual layout and the predicted gain-bandwidth product of our KITWPAs.

RESUMEN

Glycolysis is the central metabolic pathway of almost every cell and organism. Under appropriate conditions, glycolytic oscillations may occur in individual cells as well as in entire cell populations or tissues. In many biological systems, glycolytic oscillations drive coherent oscillations of other metabolites, for instance in cardiomyocytes near anorexia, or in pancreas where they lead to a pulsatile release of insulin. Oscillations at the population or tissue level require the cells to synchronize their metabolism. We review the progress achieved in studying a model organism for glycolytic oscillations, namely yeast. Oscillations may occur on the level of individual cells as well as on the level of the cell population. In yeast, the cell-to-cell interaction is realized by diffusion-mediated intercellular communication via a messenger molecule. The present mini-review focuses on the synchronisation of glycolytic oscillations in yeast. Synchronisation is a quorum-sensing phenomenon because the collective oscillatory behaviour of a yeast cell population ceases when the cell density falls below a threshold. We review the question, under which conditions individual cells in a sparse population continue or cease to oscillate. Furthermore, we provide an overview of the pathway leading to the onset of synchronized oscillations. We also address the effects of spatial inhomogeneities (e.g., the formation of spatial clusters) on the collective dynamics, and also review the emergence of travelling waves of glycolytic activity. Finally, we briefly review the approaches used in numerical modelling of synchronized cell populations.

Asunto(s)

RESUMEN

Stimuli-evoked and spontaneous brain activity propagates across the cortex in diverse spatiotemporal patterns. Despite extensive studies, the relationship between spontaneous and evoked activity is poorly understood. We investigate this relationship by comparing the amplitude, speed, direction, and complexity of propagation trajectories of spontaneous and evoked activity elicited with visual, auditory, and tactile stimuli using mesoscale wide-field imaging in mice. For both spontaneous and evoked activity, the speed and direction of propagation is modulated by the amplitude. However, spontaneous activity has a higher complexity of the propagation trajectories. For low stimulus strengths, evoked activity amplitude and speed is similar to that of spontaneous activity but becomes dissimilar at higher stimulus strengths. These findings are consistent with observations that primary sensory areas receive widespread inputs from other cortical regions, and during rest, the cortex tends to reactivate traces of complex multisensory experiences that might have occurred in exhibition of different behaviors.

Asunto(s)

RESUMEN

Since Noble adapted in 1962 the model of Hodgkin and Huxley to fit Purkinje fibres, the refinement of models for cardiomyocytes has continued. Most of these models are high-dimensional systems of coupled equations so that the possible mathematical analysis is quite limited, even numerically. This has inspired the development of reduced, phenomenological models that preserve qualitatively the main feature of cardiomyocyte's dynamics. In this paper, we present a systematic comparison of the dynamics between two notable low-dimensional models, the FitzHugh-Nagumo model (FitzHugh in Bull Math Biophys 17:257-269, 1955, J Gen Physiol 43:867-896, 1960, Biophys J 1:445-466, 1961) as a prototype of excitable behaviour and a polynomial version of the Karma model (Karma in Phys Rev Lett 71(7):16, 1993, Chaos 4:461, 1994) which is specifically developed to fit cardiomyocyte's behaviour well. We start by introducing the models and considering their pure ODE versions. We analyse the ODEs employing the main ideas and steps used in the setting of geometric singular perturbation theory. Next, we turn to the spatially extended models, where we focus on travelling wave solutions in 1D. Finally, we perform numerical simulations of the 1D PDE Karma model varying model parameters in order to systematically investigate the impact on wave propagation velocity and shape. In summary, our study provides a reference regarding key similarities as well as key differences of the two models.

Asunto(s)

RESUMEN

Alpha rhythms (∼10Hz) in the human brain are classically associated with idling activities, being predominantly observed during quiet restfulness with closed eyes. However, recent studies demonstrated that alpha (∼10Hz) rhythms can directly relate to visual stimulation, resulting in oscillations, which can last for as long as one second. This alpha reverberation, dubbed perceptual echoes (PE), suggests that the visual system actively samples and processes visual information within the alpha-band frequency. Although PE have been linked to various visual functions, their underlying mechanisms and functional role are not completely understood. In this study, we investigated the relationship between conscious perception and the generation and the amplitude of PE. Specifically, we displayed two coloured Gabor patches with different orientations on opposite sides of the screen, and using a set of dichoptic mirrors, we induced a binocular rivalry between the two stimuli. We asked participants to continuously report which one of two Gabor patches they consciously perceived, while recording their EEG signals. Importantly, the luminance of each patch fluctuated randomly over time, generating random sequences from which we estimated two impulse-response functions (IRFs) reflecting the PE generated by the perceived (dominant) and non-perceived (suppressed) stimulus, respectively. We found that the alpha power of the PE generated by the consciously perceived stimulus was comparable with that of the PE generated during monocular vision (control condition) and higher than the PE induced by the suppressed stimulus. Moreover, confirming previous findings, we found that all PEs propagated as a travelling wave from posterior to frontal brain regions, irrespective of conscious perception. All in all our results demonstrate a correlation between conscious perception and PE, suggesting that the synchronization of neural activity plays an important role in visual sampling and conscious perception.