2026-07-14 | | Total: 54
We consider a two-species reaction-diffusion system in one space dimension that is derived from an epidemiological model in a spatially periodic environment with two types of pathogens: the wild type and the mutant. The system is of a hybrid nature, partly cooperative and partly competitive, but neither of these entirely. As a result, the comparison principle does not hold. In the previous work, we studied the propagation properties of the solutions to the Cauchy problem for this system and showed, among other things, that the spreading speeds of the fronts to the right and to the left directions, denoted by $ c^*_R$ and $ c^*_L$, can be characterized by using certain principal eigenvalues, and studied the homogenization limit as the spatial period $L$ tends to $0$, and also discussed the long-time behavior of solutions behind the fronts. In the present paper we prove the existence of pulsating traveling waves in the right direction (resp. left direction) with speed $c$ for any $c\geq c^*_R$ (resp. $c\geq c^*_L$), where $c^*_R$ and $c^*_L$ denote the aforementioned spreading speeds in the right and left directions. We also prove that the leading edge of any traveling wave has the exponential decay rate that is anticipated from formal linear analysis, thus extending part of the results of Hamel 2008 to systems of equations. Finally, we present an example in which the two speeds $c^*_R$ and $c^*_L$ are different. This is done by considering a multi-scale singular limit problem. This result highlights a marked difference between our system and scalar KPP type equations.
In earlier work, we posed a stability question for Serrin's overdetermined problem under Dirichlet perturbations and proved that the answer is negative in dimensions $n\ge3$. Here we resolve the question in the planar convex class and obtain a sharp quantitative theory without any a priori geometric nondegeneracy. Let $u_Ω$ solve \[ -Δu_Ω=1\ \text{in }Ω,\qquad \partial_νu_Ω=-\frac{|Ω|}{P(Ω)}\ \text{on }\partialΩ, \qquad \int_{\partialΩ}u_Ω\,dσ=0, \] and set $O(Ω):=\text{osc}_{\partial Ω}u_Ω$. We construct fixed-area annuli with $O(Ω_k)\to0$ that remain far from every disk, showing that convexity is essential in dimension two. By contrast, if $Ω_k\subset\mathbb R^2$ are convex, $|Ω_k|=π$, and $O(Ω_k)\to0$, then, up to translations, $Ω_k$ converges in Hausdorff distance to the unit disk. Moreover, \[ R_Ω-r_Ω+\inf_{z\in\mathbb R^2}d_H(Ω,B_1(z)) \le C\,O(Ω) \] for all planar convex $Ω$ with $|Ω|=π$ and sufficiently small $O(Ω)$, and the linear order is optimal. The proof combines a new mechanism excluding long-thin degeneration, the rough-domain Serrin rigidity theorem of Figalli--Zhang, new tangential-gradient and linear boundary-growth estimates, a boundary $P$-function estimate, and the reverse-Serrin identity of Magnanini--Molinarolo--Poggesi. We also study the weaker deficit \[ A(Ω):=\frac1{P(Ω)}\int_{\partialΩ}u_Ω,dσ-\min_{\partialΩ}u_Ω. \] In the planar convex class, $A(Ω_k)\to0$ still forces convergence to a disk, and \[ R_Ω-r_Ω+\inf_z d_H(Ω,B_1(z)) \le C A(Ω)^{2/3} \] for $|Ω|=π$ and sufficiently small $A(Ω)$.
We prove global well-posedness and asymptotic stability for the three-dimensional Peskin problem, which models a closed, elastic membrane immersed in an incompressible Stokes fluid. We work with initial data in the optimal regularity space $W^{1,\infty}(\mathbb{S}^2)$, which may contain infinitely many corners. These initial configurations are instantly desingularized by the flow's parabolic smoothing effect, becoming smooth for all $t > 0$. Then we establish that the solutions converge exponentially in the $C^1$ topology to a translated and dilated conformal sphere. The stability is achieved by combining our nonlinear estimates with an exact structural decoupling of the 10-dimensional manifold of conformal steady states, demonstrating that the infinite-dimensional dissipative perturbation is strictly controlled. The core of our analysis is a functional framework on the sphere $\mathbb{S}^2$ that uses spectral Littlewood-Paley projections to control the highly singular multilinear operators arising from the fluid nonlinearity
This paper investigates the dependence of the control cost for a wave equation with respect to perturbation by a time-independent potential $\lmbd V$ scaled by a large parameter $\lmbd$ on a compact Riemannian manifold. We introduce the geometric control condition~\eqref{GCC+}, a variant of the geometric control condition of Bardos--Lebeau--Rauch--Taylor, tailored to accommodate the influence of the potential $V$. We show that~\eqref{GCC+} is necessary and sufficient for the existence of a uniform \emph{observability cost} with respect to the large parameter $\lmbd$. We provide geometric examples satisfying~\eqref{GCC+} and estimate the blow-up rate of the \emph{observability cost} in situations where it fails. The proofs rely on semiclassical and second microlocal defect measures.
We study a densely connected system of excitatory integrate-and-fire neurons which are subject to common noise. The system incorporates a gradual transmission of action potentials and captures the re- and hyperpolarization phases through a random refractory period followed by a reset to a randomized level below the rest potential. As the number of neurons tends to infinity, we show that there is weak convergence to a unique membrane potential density governed by a stochastic Fokker--Planck equation with a well-defined spike transmission rate. The latter is driven by the mean cumulative spike count, which is non-differentiable but shown to satisfy a generalized flux condition. We obtain the uniqueness of the Fokker--Planck equation from energy estimates in the dual of the first Sobolev space. Finally, we give a conditional McKean--Vlasov representation of the membrane potential density as the law of a representative neuron given the common noise.
We construct homogeneous solutions with non-Lipschitz Hessian for finite, constant-coefficient Bellman equations. First, for every $σ\in(0,1)$, we find two uniformly elliptic matrices $A_1,A_2\in\mathcal{S}^4$ and a nonzero $(2+σ)$-homogeneous solution $u$ of \[\max\bigl\{{\rm tr}\,(A_1D^2u),{\rm tr}\,(A_2D^2u)\bigr\}=0 \qquad\text{in }\mathbb{R}^4.\] Second, in $\mathbb{R}^2$ we construct three matrices satisfying ${\rm Id}_2\leq A_j\leq3{\rm Id}_2$ for which the corresponding Bellman equation admits a homogeneous solution with a non-Lipschitz Hessian. In particular, solutions to convex fully nonlinear uniformly elliptic equations are not in $C^{2,1}$, and not even in $C^{2, 1-\varepsilon}$ for $\varepsilon > 0$ small.
In this work, we study the optimization problem for the first eigenvalue of mixed local nonlocal operators plus a potential $V$. We begin by investigating the existence and fundamental properties of the eigenvalues, with special emphasis on the first eigenvalue. Finally, we discuss the dependence of the first eigenvalue on the potential function and establish the existence of optimal potentials within certain admissible classes. In particular, we show the existence of a unique maximizer and a minimizer of the first eigenvalue on any bounded, closed, and convex subset of $L^q(Ω)$. Moreover, these results enable us to characterize the maximizers and minimizers in the closed unit ball of $L^q(Ω)$, as well as in the class of rearrangements of any $V\in L^q(Ω)$.
We introduce graph-compatible powers determined by the coupling graph for multicomponent weakly coupled elliptic systems and establish concavity of the associated transformed components through a concave-envelope method and a weighted viscosity comparison principle. We further establish a componentwise constant-rank theorem for the transformed Hessians, which yields strict power concavity under additional assumptions.
This paper discusses about an inverse scattering problem to the nonlocal nonlinear Schrodinger equation. We extend the previous result treated by Y. Zhao and E. Fan of inverse scattering result at the point with slowly decaying data. More precisely, we show the global existence of solution in $H^{1,s}$ for $s>1/2$ by assuming small data in $L^1$. The proof relies on the inverse scattering method based on the ZS-AKNS (Zakharov-Shabat, Ablowitz-Kaup-Newell-Segur) eigenvalue problem on the spectral theory of ordinary differential equations upon the framework of slowly decaying data.
We consider nonlinear diffusion equations on an interval where the diffusion coefficient in a small region near the center is scaled such that it approximates a transmission condition at a membrane. While the limiting behavior of the solutions is well-understood, we study the convergence in the sense of the energy-dissipation principle (EDP) of associated gradient structures given in terms of a free energy and a dissipation potential. EDP-convergence provides a uniquely specified limiting gradient structure that reformulates the transmission condition in terms of an effective kinetic relation for the membrane, which relates the jump of the chemical potential and the flux through the membrane. We show how \AAA properties of the chosen free energy and the mobility of small-scale diffusion migrate to the effective kinetic relation. A surprising result is that starting from the linear Onsager relation of Otto's gradient structure for the linear diffusion equation, one obtains an exponentially growing kinetic relation, the so-called Marcelin-De Donder kinetic.
We prove the Harnack inequality (with tails) and local Hölder regularity for the fractional logarithmic $p$-Laplace operator, which is derived by differentiating the fractional $p$-Laplace operator with respect to its order. To be more precise, for a suitable function $u,$ the operator reads as the first order derivative \begin{align*} (-Δ_p)^{s+\log} u:= \frac{\rm d}{{\rm d}t}(-Δ_p)^t u \Big|_{t=s} \end{align*} at any arbitrary order $s\in (0, 1).$ The kernel of this operator involves a logarithmic factor that changes sign at large scales and, near the diagonal, is more singular than the kernel of the fractional $p$-Laplacian. To achieve our regularity estimates, we adopt the classical De Giorgi-Nash-Moser techniques in this setting. We also construct an example showing that the Harnack inequality fails without tail terms. Our results are new even in the linear setup $p=2$.
We prove the minimisers of the Lieb-Oxford bound (if any) for a fixed (and finite) number of particles $N \geq 1$ are necessarily compactly supported, extending a result proved by E. H. Lieb and S. Oxford in the one-particle case.
We study Nakao's problem in two and three space dimensions, a weakly coupled system consisting of a semilinear damped wave equation and a semilinear undamped wave equation. We establish the global in-time existence and uniqueness of small data mild Sobolev solutions in new admissible ranges, together with time-dependent estimates matching those for the corresponding linearized problems. The proof combines diffusion-type $L^m-L^r$ estimates for the damped component with the Poisson and Kirchhoff formulas for the undamped component. Dimension-dependent solution spaces are introduced to incorporate the weaker wave decay and logarithmic $L^2$ growth in two space dimensions. In three space dimensions, a two-level fixed point argument avoids an artificial restriction on $p$ by establishing the self-map property in a strong space and the contraction in a weaker metric. Furthermore, the undamped component scatters to a free wave in $H^2\times H^1$.
This study presents a three-phase suspension lava flow model with a Herschel-Bulkley rheology. The suspension contains crystals and gas bubbles, and two closures for the evolution of the crystal volume fraction are considered. The first closure minimizes the complexity of the system by treating crystal fraction as a transported quantity subject to relaxation towards an equilibrium state. This closure avoids a parametrization of the many heat transfer mechanisms a lava flow is subjected to. The other closure is on the lava temperature considering four heat transfer mechanisms and a prescribed temperature--crystallinity relationship. We deduce from this system and solve numerically a one-dimensional depth-averaged model. A comparison with a pre-existing model based on real lava flow data suggests that the prediction of flow parameters done with the multi-parametric evolution of temperature yields more accurate results than those obtained with the other closure. The transient nature of our model correctly predicts that confined lava traveling down an irregular steep slope yields a series of cascading fill-then-breakout lumps that causes the overall flow to be pulsatory. These pulses dominate the local dynamics and preclude a strict steady state to be reached. Theoretically, taking gas bubbles into account is best done with a general rheological relationship valid at any capillary number. In the conditions explored herein, bubbles modulate viscosity within a factor 2 with a shear thinning behavior, decelerating slow flows and accelerating fast flows. When a simplified rheology treating bubbles as hard spheres was used, the only dynamic parameter affected was bulk viscosity.
Let $H=-Δ+V(|x|)$ be a nonnegative radial Schrödinger operator on $\mathbb{R}^d$, $d\ge 2$, whose positive harmonic function satisfies $U(r)\simeq r^{-σ_0}$ for $0<r\le 1$ and $U(r)\simeq r^{-σ_\infty}$ for $r\ge 1$, with $-d/2<σ_0,σ_\infty<d/2$. Assuming the two-sided ground-state heat-kernel estimate of Ishige, Kabeya, and Ouhabaz, we determine the maximal open range in which the kernel of $H^{-s/2}$ admits the clean two-sided estimate $K_s^H(x,y)\simeq |x-y|^{s-d}U(|x|)U(|y|)/[U(|x|+|x-y|)U(|y|+|x-y|)]$, namely $0<s<\min\{d,d-2σ_0,d-2σ_\infty\}$. In this range we give a complete necessary-and-sufficient classification of the broken-power estimate $\|w_{-β_0,-β_\infty}H^{-s/2}f\|_{L^{q,v}}\lesssim \|w_{α_0,α_\infty}f\|_{L^{p,u}}$ for $1<p,q<\infty$ and $1\le u,v\le\infty$. The result covers signed ground-state exponents, the full range $q<p$, all one-sided weight equalities, both scale equalities, and simultaneous endpoint corners. Scale equality is governed by $u\le v$, including when $q<p$; an input or output power endpoint requires respectively $u=1$ or $v=\infty$; and at a same-side power/scale corner the only admissible pair is $(u,v)=(1,\infty)$. The proof combines clean-kernel analysis, local Lorentz-Hardy-Littlewood-Sobolev estimates, rank-one endpoint arguments, geometric annular sequence spaces, a triangular matrix theorem, and a nine-block decomposition.
In this paper, we investigate the existence, multiplicity, and semiclassical concentration of normalized solutions to a critical Schrödinger--Poisson system with a singular Hardy potential in \(\mathbb{R}^3\). More precisely, we consider \[ \begin{cases} -\varepsilon^2Δu+ \left(V(x)-\dfrac{κ\varepsilon^2}{|x|^2}\right)u -φ|u|^3u =λu+μ|u|^{q-2}u+|u|^4u, & \text{in } \mathbb{R}^3, \\[1mm] -\varepsilon^2Δφ=|u|^5, & \text{in } \mathbb{R}^3, \end{cases} \] under the prescribed mass constraint \[ \int_{\mathbb{R}^3}|u|^2\,dx=a^2\varepsilon^3, \] where \(a,μ>0\), \(q\in(2,10/3)\), \(\varepsilon>0\) is a small semiclassical parameter, and \(0<κ<1/4\). The parameter \(λ\in\mathbb{R}\) appears as a Lagrange multiplier associated with the mass constraint, while \(V:\mathbb{R}^3\to(0,+\infty)\) is a continuous electric potential whose minimum set is assumed to be nonempty and compact. The main difficulty stems from the simultaneous presence of the inverse-square Hardy singularity, the mass constraint, the critical local nonlinearity, and the nonlocal Poisson interaction. By combining the Hardy inequality, constrained variational methods, suitable truncation arguments, and concentration-compactness techniques, we first establish the existence of a normalized ground state for sufficiently small mass and sufficiently small \(\varepsilon\). We then employ Ljusternik--Schnirelmann category theory to obtain multiple normalized solutions whose number is related to the topology of the minimum set of \(V\). Finally, we show that the corresponding semiclassical states concentrate near the global minimum set of the electric potential as \(\varepsilon\to0\).
In this paper, we introduce and analyse a model in which a shared, age-structured resource is exploited by two competing populations through age-dependent interaction kernels. It provides a general framework for studying resource-mediated competition in epidemiology, immunology and ecology. We use integrated semigroups theory to establish the well-posedness of the model and its fundamental properties, including positivity and the existence of a global attractor. Then, by means of Lyapunov functionals, persistence results and generalised LaSalle principle, we fully characterize the asymptotic dynamics in terms of the competitors's basic reproduction numbers and invasive fitnesses. Our results demonstrate that in addition to the classical competitive exclusion, the system can exhibit globally asymptotically stable coexistence, an outcome allowed by a mechanism of age-mediated niche partitioning.
We prove existence and uniqueness in BV for a Vlasov-Poisson boundary value problem which takes into account a nonlinear boundary condition at the boundary on the electric potential. It models the dynamical interaction between a plasma and an electrically isolated metallic wall. The nonlinear boundary condition is the so-called floating potential equation of plasma physics.
We consider the inhomogeneous nonlinear heat equation \[ \partial_t u-Δu=|u|^{p-1}u+f(x),\qquad x\in\mathbb{R}^3,\quad p>5, \] where \(f\in L^\infty\cap C^{0,1}(\mathbb{R}^3)\). For every sufficiently large integer \(n\), we construct a codimension-\(n\) Lipschitz manifold of non-radial initial data whose corresponding solutions blow up in finite time and whose rescaled profiles converge to the prescribed self-similar profile \(Φ_n\) of the homogeneous equation. The main novelty is to show that the finite-codimensional stability mechanism for self-similar blow-up, developed in the work of Collot, Raphaël and Szeftel [Mem. Amer. Math. Soc. (2019)] for the homogeneous equation, is robust under the addition of a bounded, Lipschitz spatially inhomogeneous source term. In contrast with the homogeneous problem, the equation considered here has no exact scaling invariance, which is a key ingredient in many previous constructions. We expect that the framework developed in this work may also be useful for related problems in which exact scaling invariance is broken.
In this paper, we consider the famous Brezis-Nirenberg equation \begin{eqnarray*} \left\{ \aligned &-Δu=λu+|u|^{\frac{4}{N-2}}u,\quad&\mbox{in}\,\, Ω,\\ &u=0,\quad&\mbox{on}\,\, \partialΩ, \endaligned \right. \end{eqnarray*} where $N\geq3$ is the dimension, $Ω\subset\mathbb{R}^N$ is a bounded domain with smooth boundary $\partialΩ$ and $λ>0$ is a parameter. By developing a refined blow-up analysis based on the inverse reduction argument developed in \cite{WW2019,WW2019-2}, we classify, for the fist time, the Struwe decomposition of the Brezis-Nirenberg equation in the one-bubble case as the parameter $λ$ varies for $N\geq4$. As applications, we prove that the $4d$ Brezis-Nirenberg equation has a nontrivial solution (least energy solution) for $λ\inσ(-Δ)$ in general bounded domains, where $σ(-Δ)$ is the spectrum of $-Δ$ in $H^1_0(Ω)$. Our result completes the existence theory of the Brezis-Nirenberg equation for $N\geq4$ in \cite{AP2025,CFP1985,CFS,CSS1986,CW2005,CSZ2012,SWW2009,TYZ2022} since 1984.
This paper establishes the direct scattering transform for the Calogero-Moser derivative nonlinear Schrödinger (CMDNLS) equation on the real line. For potentials in a weighted Hardy-Sobolev space, we prove the existence and uniqueness of the Jost functions associated with the Lax operator, and construct the corresponding scattering coefficients including the transmission coefficient and the reflection coefficient. Under the CMDNLS flow, we determine the time evolution of these Jost functions and scattering data, revealing a simple dynamics: the transmission coefficient is invariant, and the reflection coefficient acquires a purely rotational phase. Furthermore, we establish the invariance of a suitable weighted Sobolev space under the flow, which guarantees the persistence of the required regularity for the solutions.
We investigate the existence of normalized solutions to the gradient-type Schrödinger system \begin{equation*} \begin{cases} -Δu+ V_1(x)u+λu= uv^2 & \text{ in } Ω,\\ -Δv+ V_2(x)v+λv= u^2v & \text{ in } Ω%\frac{\partial u}{\partial ν}=\frac{\partial v}{\partial ν}=0 \, & \text{ on } \partial Ω\end{cases} \end{equation*} subject to the mass constraint $\int_Ω\left(|u|^2+|v|^2 \right)dx=a>0$ and Neumann boundary conditions, where $Ω\subset \mathbb{R}^3$ is a smooth bounded domain, each $V_i$ is continuous, and $λ$ is a Lagrange multiplier. Applying a minimax principle that incorporates Morse index information, we establish the existence of nontrivial normalized solutions of mountain pass type. The proof is based on a refined blow-up analysis adapted to such gradient-type systems, together with new Liouville-type theorems for finite Morse index solutions of the associated limit systems in $\mathbb{R}^3$ and $\mathbb{R}^3_+$.
This paper provides the study of the homogenization of the 3D non-homogeneous incompressible heat-conducting fluid in perforated domains with holes of subcritical size, where the viscosity and the heat conductivity coefficient are assumed to depend on the temperature. The diameter of the holes is of order $\varepsilon^α \ (α>3)$, where $\varepsilon > 0$ is a small parameter that measures the mutual distance between the holes. We prove that as $\varepsilon\to 0$, the limit behavior of velocity, density and temperature is governed by the original system in homogeneous domain without holes.
We prove the existence of global weak entropy solutions to the $p$-system with the piecewise affine flux function. The construction is based on a wave-front tracking scheme for which we identify a necessary and sufficient condition for the occurrence of infinitely many fronts. For smooth fluxes, the theory is well established. However, the lower Lipschitz regularity of the flux requires a new method to obtain the interaction estimates. We also prove the stability of the solution in the sense of \cite{bia-col-02}, showing that solutions are Lipschitz continuous with respect to the $L^{\infty}$ norm of the difference of the flux derivatives. We obtain the convergence rate of the solutions from the piecewise affine flux to smooth flux by using the stability estimate.
For nonlinear Fokker-Planck equations with mobility, the Wasserstein gradient flow structure is described by the generalized relative entropy as the energy functional and the modified Wasserstein metric $W_h$ as the associated metric structure. This work investigates the nonlinear effects induced by mobility and establishes the corresponding inequalities. For nonlinear diffusion, we establish a logarithmic Sobolev inequality, which yields the convergence rate of the free energy functional and the Talagrand inequality. By further exploiting the relationship between the weighted homogeneous Sobolev norm and the $W_h$ metric, we derive an HWI inequality relating the relative entropy, the $W_h$ metric, and the Fisher information. In the case of mobility dependent drift and linear diffusion, the convergence rate in the $W_h$ metric is also obtained by applying the Girsanov theorem.