We study the global well-posedness and uniform boundedness of a two-dimensional reaction–advection–diffusion system with nonlinear advection. This strongly coupled system of nonlinear partial differential equations represents the continuum of a 2D lattice model designed to describe residential burglary, where each location is characterised by a tractability value that varies in both space and time. We show that the model with sublinear advection enhancement is globally well-posed, with a unique solution that is classical and uniformly bounded in time. Our results provide valuable insights into the development of urban crime models with nonlinear advection enhancements, making them suitable for broader applications, including nonlocal or heterogeneous near-repeat victimisation effects.

We study the community detection problem on a Gaussian mixture model, in which vertices are divided into $k\geq 2$ distinct communities. The major difference in our model is that the intensities for Gaussian perturbations are different for different entries in the observation matrix, and we do not assume that every community has the same number of vertices. We explicitly find the necessary and sufficient conditions for the exact recovery of the maximum likelihood estimation, which can give a sharp phase transition for the exact recovery even though the Gaussian perturbations are not identically distributed; see Section 7. Applications include the community detection on hypergraphs.

In this work, we study early warning signs for stochastic partial differential equations (SPDEs), where the linearisation around a steady state is characterised by continuous spectrum. The studied warning sign takes the form of qualitative changes in the variance as a deterministic bifurcation threshold is approached via parameter variation. Specifically, we focus on the scaling law of the variance near the transition. Since we are dealing here, in contrast to previous studies, with the case of continuous spectrum and quantitative scaling laws, it is natural to start with linearisations of the drift operator that are multiplication operators defined by analytic functions. For a one-dimensional spatial domain, we obtain precise rates of divergence. In the case of the two- and three-dimensional domains, an upper bound to the rate of the early warning sign is proven. These results are cross-validated by numerical simulations. Our theory can be generically useful for several applications, where stochastic and spatial aspects are important in combination with continuous spectrum bifurcations.

Motivated by the impact of worsening climate conditions on vegetation patches, we study dynamic instabilities in an idealised Ginzburg–Landau model. Our main results predict time instances of sudden drops in wavenumber and the resulting target states. The changes in wavenumber correspond to the annihilation of individual vegetation patches when resources are scarce and cannot support the original number of patches. Drops happen well after the primary pattern has destabilised at the Eckhaus boundary and key to distinguishing between the disappearance of 1,2 or more patches during the drop are complex spatio-temporal resonances in the linearisation at the unstable pattern. We support our results with numerical simulations and expect our results to be conceptually applicable universally near the Eckhaus boundary, in particular in more realistic models.

A wake model is pursued for potential flow past a submerged, finite-length plate that is perpendicular to a uniform, horizontal stream bounded above by a free surface. The effects of gravity are included along the free surface. The approach is to adopt an open-wake model such that the wake boundaries become parallel to the undisturbed stream at some (unknown) point downstream. Boundary integral equations are formed and then discretised along the wake boundaries and free surface in order to obtain a solution numerically. In terms of the dependency of the solution on various parameters, the problem will be formulated in two ways. First, for a given Froude number and ratio of the length of the vertical plate to the draft (the depth of the bottom of the vertical plate relative to the undisturbed free surface), the effect of the wake underpressure coefficient on the size of the wake will be considered. Then, the problem will be discussed where we instead (more naturally) fix the Froude number, draft and length of the vertical, submerged plate. The dependencies of the solution on these parameters will be analysed regarding the effects on several factors, including the size of the wake, the relative lengths of the upper and lower wake boundaries, and the resulting wake underpressure coefficient.

We investigate a recent model proposed in the literature elucidating patterns driven by chemotaxis, similar to viscous fingering phenomena. Notably, this model incorporates a singular advection term arising from a modified formulation of Darcy’s law. It is noteworthy that this type of advection can also be well interpreted as a description of a radial fluid flow source surrounding an aggregation of cells. For the two-dimensional scenario, we establish a precise threshold delineating between blow-up and global solution existence. This threshold is contingent upon the pressure magnitude and the initial total mass of the aggregating cells.

The connection between Residual Neural Networks (ResNets) and continuous-time control systems (known as NeurODEs) has led to a mathematical analysis of neural networks, which has provided interesting results of both theoretical and practical significance. However, by construction, NeurODEs have been limited to describing constant-width layers, making them unsuitable for modelling deep learning architectures with layers of variable width. In this paper, we propose a continuous-time Autoencoder, which we call AutoencODE, based on a modification of the controlled field that drives the dynamics. This adaptation enables the extension of the mean-field control framework originally devised for conventional NeurODEs. In this setting, we tackle the case of low Tikhonov regularisation, resulting in potentially non-convex cost landscapes. While the global results obtained for high Tikhonov regularisation may not hold globally, we show that many of them can be recovered in regions where the loss function is locally convex. Inspired by our theoretical findings, we develop a training method tailored to this specific type of Autoencoders with residual connections, and we validate our approach through numerical experiments conducted on various examples.

In this paper, we consider a delayed discrete single population patch model in advective environments. The individuals are subject to both random and directed movements, and there is a net loss of individuals at the downstream end due to the flow into a lake. Choosing time delay as a bifurcation parameter, we show the existence of Hopf bifurcations for the model. In homogeneous non-advective environments, it is well known that the first Hopf bifurcation value is independent of the dispersal rate. In contrast, for homogeneous advective environments, the first Hopf bifurcation value depends on the dispersal rate. Moreover, we show that the first Hopf bifurcation value in advective environments is larger than that in non-advective environments if the dispersal rate is large or small, which suggests that directed movements of the individuals inhibit the occurrence of Hopf bifurcations.

We study the poor-biased model for money exchange introduced in Cao & Motsch ((2023) Kinet. Relat. Models 16(5), 764–794.): agents are being randomly picked at a rate proportional to their current wealth, and then the selected agent gives a dollar to another agent picked uniformly at random. Simulations of a stochastic system of finitely many agents as well as a rigorous analysis carried out in Cao & Motsch ((2023) Kinet. Relat. Models 16(5), 764–794.), Lanchier ((2017) J. Stat. Phys. 167(1), 160–172.) suggest that, when both the number of agents and time become large enough, the distribution of money among the agents converges to a Poisson distribution. In this manuscript, we establish a uniform-in-time propagation of chaos result as the number of agents goes to infinity, which justifies the validity of the mean-field deterministic infinite system of ordinary differential equations as an approximation of the underlying stochastic agent-based dynamics.

We study an epidemic patch model that describes the disease spread in population with variable latency due to the differences in immunologic tolerance between individuals. We focus on whether the disease can spread in space that leads to the emergence of epidemic wave, that is the travelling wave solution with constant speed. We first establish some properties of the linearized wave profile equations, which are helpful in obtaining the priori estimates of travelling waves and wave speeds. Then, applying the truncation method and limiting arguments, we can obtain threshold propagation dynamics of the epidemic model. Our result gives a complete characterization of the existence, nonexistence and minimal wave speed of travelling waves. To the best of our knowledge, this is the first time to characterize the propagation dynamics of epidemic patch model with variable latency, which contributes to the understanding of the transmission phenomenon of disease.

This article offers an advanced and novel investigation into the intricate propagation dynamics of the Belousov–Zhabotinsky system with non-local delayed interaction, which exhibits dynamical transition structure from bistable to monostable. We first solved the enduring open problem concerning the existence, uniqueness and the speed sign of the bistable travelling waves. In the monostable case, we developed and derived new results for the minimal wave speed selection, which, as an application, further improved the existing investigations on pushed and pulled wavefronts. Our results can provide new estimate to the minimal speed as well as to the determinacy of the transition parameters. Moreover, these results can be directly applied to standard localised models and delayed reaction diffusion models by choosing appropriate kernel functions.

This paper investigates the separation property in binary phase-segregation processes modelled by Cahn-Hilliard type equations with constant mobility, singular entropy densities and different particle interactions. Under general assumptions on the entropy potential, we prove the strict separation property in both two and three-space dimensions. Namely, in 2D, we notably extend the minimal assumptions on the potential adopted so far in the literature, by only requiring a mild growth condition of its first derivative near the singular points $\pm 1$, without any pointwise additional assumption on its second derivative. For all cases, we provide a compact proof using De Giorgi’s iterations. In 3D, we also extend the validity of the asymptotic strict separation property to the case of fractional Cahn-Hilliard equation, as well as show the validity of the separation when the initial datum is close to an ‘energy minimizer’. Our framework offers insights into statistical factors like particle interactions, entropy choices and correlations governing separation, with broad applicability.

The purpose of the paper is twofold. First, we show that partial-data transmission eigenfunctions associated with a conductive boundary condition vanish locally around a polyhedral or conic corner in $\mathbb{R}^n$, $n=2,3$. Second, we apply the spectral property to the geometrical inverse scattering problem of determining the shape as well as its boundary impedance parameter of a conductive scatterer, independent of its medium content, by a single far-field measurement. We establish several new unique recovery results. The results extend the relevant ones in [26] in two directions: first, we consider a more general geometric setup where both polyhedral and conic corners are investigated, whereas in [26] only polyhedral corners are concerned; second, we significantly relax the regularity assumptions in [26] which is particularly useful for the geometrical inverse problem mentioned above. We develop novel technical strategies to achieve these new results.

We formulate haptotaxis models of cancer invasion wherein the infiltrating cancer cells can occupy a spectrum of states in phenotype space, ranging from ‘fully mesenchymal’ to ‘fully epithelial’. The more mesenchymal cells are those that display stronger haptotaxis responses and have greater capacity to modify the extracellular matrix (ECM) through enhanced secretion of matrix-degrading enzymes (MDEs). However, as a trade-off, they have lower proliferative capacity than the more epithelial cells. The framework is multiscale in that we start with an individual-based model that tracks the dynamics of single cells, which is based on a branching random walk over a lattice representing both physical and phenotype space. We formally derive the corresponding continuum model, which takes the form of a coupled system comprising a partial integro-differential equation for the local cell population density function, a partial differential equation for the MDE concentration and an infinite-dimensional ordinary differential equation for the ECM density. Despite the intricacy of the model, we show, through formal asymptotic techniques, that for certain parameter regimes it is possible to carry out a detailed travelling wave analysis and obtain invading fronts with spatial structuring of phenotypes. Precisely, the most mesenchymal cells dominate the leading edge of the invasion wave and the most epithelial (and most proliferative) dominate the rear, representing a bulk tumour population. As such, the model recapitulates similar observations into a front to back structuring of invasion waves into leader-type and follower-type cells, witnessed in an increasing number of experimental studies over recent years.

We study the Cauchy problem on the real line for the nonlocal Fisher-KPP equation in one spatial dimension,

\begin{equation*} u_t = D u_{xx} + u(1-\phi *u), \end{equation*}

$\phi *u$

$\phi (y) \equiv H\left (\frac{1}{4}-y^2\right )$

$u=1$

$D$

$\Delta _1 \approx 0.00297$

$D \ll 1$

$O(1)$

$u=O(1)$

$u$

$D \to 0^+$

$D \geq \Delta _1$

$2 \sqrt{D}$

$0 \lt D \lt \Delta _1$

$u=0$

$D$

Opinion dynamics is an important and very active area of research that delves into the complex processes through which individuals form and modify their opinions within a social context. The ability to comprehend and unravel the mechanisms that drive opinion formation is of great significance for predicting a wide range of social phenomena such as political polarisation, the diffusion of misinformation, the formation of public consensus and the emergence of collective behaviours. In this paper, we aim to contribute to that field by introducing a novel mathematical model that specifically accounts for the influence of social media networks on opinion dynamics. With the rise of platforms such as Twitter, Facebook, and Instagram and many others, social networks have become significant arenas where opinions are shared, discussed and potentially altered. To this aim after an analytical construction of our new model and through incorporation of real-life data from Twitter, we calibrate the model parameters to accurately reflect the dynamics that unfold in social media, showing in particular the role played by the so-called influencers in driving individual opinions towards predetermined directions.

In a smoothly bounded domain $\Omega \subset \mathbb{R}^n$, $n\ge 1$, this manuscript considers the homogeneous Neumann boundary problem for the chemotaxis system

\begin{eqnarray*} \left \{ \begin{array}{l} u_t = \Delta u - \nabla \cdot (u\nabla v), \\[5pt] v_t = \Delta v + u - \alpha uv, \end{array} \right . \end{eqnarray*}

$\alpha \gt 0$

It is shown that there exists $\delta _\star =\delta _\star (n)\gt 0$ such that for any given $\alpha \ge \frac{1}{\delta _\star }$ and for any suitably regular initial data satisfying $v(\cdot, 0)\le \delta _\star$, this problem admits a unique classical solution that stabilizes to the constant equilibrium $(\frac{1}{|\Omega |}\int _\Omega u(\cdot, 0), \, \frac{1}{\alpha })$ in the large time limit.

Contemporary epidemiological models often involve spatial variation, providing an avenue to investigate the averaged dynamics of individual movements. In this work, we extend a recent model by Vaziry, Kolokolnikov, and Kevrekidis [Royal Society Open Science 9 (10), 2022] that included, in both infected and susceptible population dynamics equations, a cross-diffusion term with the second spatial derivative of the infected population density. Diffusion terms of this type occur, for example, in the Keller–Siegel chemotaxis model. The presented model corresponds to local orderly commute of susceptible and infected individuals and is shown to arise in two dimensions as a limit of a discrete process. The present contribution identifies and studies specific features of the new model’s dynamics, including various types of infection waves and buffer zones protected from the infection. The model with vital dynamics additionally exhibits complex spatio-temporal behaviour that involves the generation of quasiperiodic infection waves and emergence of transient strongly heterogeneous patterns.

$c$-cyclical monotonicity is the most important optimality condition for an optimal transport plan. While the proof of necessity is relatively easy, the proof of sufficiency is often more difficult or even elusive. We present here a new approach, and we show how known results are derived in this new framework and how this approach allows to prove sufficiency in situations previously not treatable.