2025-06-16 | | Total: 23
We study minimizers $\boldsymbol{m}\colon \mathbb R^2\to\mathbb S^2$ of the energy functional \begin{align*} E_\sigma(\boldsymbol{m}) = \int_{\mathbb R^2} \bigg(\frac 12 |\nabla\boldsymbol{m}|^2 +\sigma^2 \boldsymbol{ m} \cdot \nabla \times\boldsymbol{m} +\sigma^2 m_3^2 \bigg)\, dx\,, \end{align*} for $0<\sigma\ll 1$, with prescribed topological degree \begin{align*} Q(\boldsymbol{m})=\frac{1}{4\pi} \int_{\mathbb R^2}\boldsymbol{m} \cdot \partial_1 \boldsymbol{m}\times\partial_2\boldsymbol{m}\, dx =\pm 1\,. \end{align*} This model arises in thin ferromagnetic films with Dzyaloshinskii-Moriya interaction and easy-plane anisotropy, where these minimizers represent bimeron configurations. We prove their existence, and describe them precisely as perturbations of specific Möbius maps: we establish in particular that they are localized at scale of order $1/|\ln(\sigma^2)|$. The proof follows a strategy introduced by Bernand-Mantel, Muratov and Simon (Arch. Ration. Mech. Anal., 2021) for a similar model with easy-axis anisotropy, but requires several adaptations to deal with the less coercive easy-plane anisotropy and different symmetry properties.
The natural $H^1(\Omega)$ energy norm for Helmholtz problems is weighted with the wavenumber modulus $\sigma$ and induces natural weighted norms on the trace spaces $H^{\pm1/2}(\Gamma)$ by minimial extension to $\Omega\subset\mathbb R^n$. This paper presents a rigorous analysis for these trace norms with an explicit characterisation by weighted Sobolev-Slobodeckij norms and scaling estimates, highlighting their dependence on the geometry of the extension set $\Omega\subset\mathbb R^n$ and the weight $\sigma$. The analysis identifies conditions under which these trace norms are intrinsic to the isolated boundary component $\Gamma\subset\partial\Omega$ and provides $\sigma$-explicit estimates for trace inequalities in weighted spaces. In these natural wavenumber-weighted norms, the boundary integral operators allow improved continuity estimates that do \emph{not} deterioriate as $\sigma\to 0$.
We study a particular class of mean field games whose solutions can be formally connected to a scalar transport equation on the Wasserstein space of measures. For this class, we construct some interesting explicit examples of non-uniqueness of Nash equilibria. We then address the selection problem of finding rational criteria by which to choose one equilibrium over others. We show that when the theory of entropy solutions is used, we can obtain explicit error estimates for the ``vanishing noise limit,'' where the error is measured in a certain norm that measures the distance between two functions over the set of empirical measures.
Mean Field Game (MFG) models implicitly assume "rational expectations", meaning that the heterogeneous agents being modeled correctly know all relevant transition probabilities for the complex system they inhabit. When there is common noise, this assumption results in the "Master equation" (a.k.a. "Monster equation"), a Hamilton-Jacobi-Bellman equation in which the infinite-dimensional density of agents is a state variable. The rational expectations assumption and the implication that agents solve Master equations is unrealistic in many applications. We show how to instead formulate MFGs with non-rational expectations. Departing from rational expectations is particularly relevant in "MFGs with a low-dimensional coupling", i.e. MFGs in which agents' running reward function depends on the density only through low-dimensional functionals of this density. This happens, for example, in most macroeconomics MFGs in which these low-dimensional functionals have the interpretation of "equilibrium prices." In MFGs with a low-dimensional coupling, departing from rational expectations allows for completely sidestepping the Master equation and for instead solving much simpler finite-dimensional HJB equations. We introduce an adaptive learning model as a particular example of non-rational expectations and discuss its properties.
This paper addresses elliptic problems involving the superposition of nonlinear fractional operators with the critical Sobolev exponent in sublinear regimes. We establish the existence of infinitely many nontrivial weak solutions using a variational framework that combines a truncation argument with the notion of genus. A central part of our analysis is the verification of the Palais--Smale $(\mathrm{PS})_c$ condition for every $q \in (1, p_{s_\sharp}^*)$, despite the challenges posed by the lack of compactness due to the critical exponent. On the one hand, our approach extends and generalizes earlier results by García Azorero and Peral Alonso \emph{[Trans. Amer. Math. Soc., 1991]} and by Da Silva, Fiscella, and Viloria \emph{[J. Differential Equations, 2024]}; on the other hand, the results obtained here complement the study of the Brezis--Nirenberg-type problem by Dipierro, Perera, Sportelli, and Valdinoci \emph{[Commun. Contemp. Math., 2024]} for $q = p$ and by Aikyn, Ghosh, Kumar, and Ruzhansky \emph{[arXiv:2504.05105, 2025]} for $p < q < p_{s_\sharp}^*$. Notably, our results are new even in the classical case $p = 2$, highlighting the broader applicability of the methods developed here.
We consider the direct and inverse scattering problem for a penetrable, isotropic obstacle with a second-order Robin boundary condition which asymptotically models the delamination of boundary of the scatterer. We develop a direct sampling method to solve the inverse shape problem by numerically recovering the scatterer. Here we assume that the corresponding Cauchy data is measured on the boundary of a region which fully contains the scatterer. Similar methods have been applied to other inverse shape problems, but they have not been studied for a penetrable, isotropic scatterer with a second-order Robin boundary condition. We also initiate the study of the corresponding transmission eigenvalue problem, which is derived from assuming zero Cauchy data is measured on the boundary of the region which fully contains the scatterer. We prove that the transmission eigenvalues for this problem are at most a discrete set. Numerical examples will be presented for the inverse shape problem and transmission eigenvalue problem in two dimensions for circular and non-circular scatterers.
In this paper, we establish a new global Hessian matrix estimate for heat-type equations on Riemannian manifolds using a Bismut-type Hessian formula. Our results feature fully explicit coefficients as well as delay / growth rate functions. These estimates yield two key applications: a novel backward weak Harnack inequality and a precise pointwise Hessian estimate for eigenfunctions.
This paper investigates the connection between the chemotaxis--Navier--Stokes system with porous medium type nonlinear diffusion and the Hele--Shaw problem in $\mathbb{R}^d$ ($d\geq2$). First, we prove the global-in-time existence of weak solutions for the Cauchy problem of the chemotaxis-Navier-Stokes system with the general initial data, uniformly in the diffusion range $m\in [3,\infty)$. Then, we rigorously justify the Hele--Shaw limit for this system as $m\rightarrow\infty$, showing the convergence to a free boundary problem of Hele--Shaw type, where the bacterium (cell) diffusion is governed by the stiff pressure law. Moreover, the complementarity relation characterizing the limiting bacterium (cell) pressure via a degenerate elliptic equation is verified by a novel application of the Hele--Shaw framework.
In this paper, we investigate the regularity of weak solutions $u\colon\Omega\to\mathbb{R}$ to elliptic equations of the type \begin{equation*} \mathrm{div}\, \nabla \mathcal{F}(x,Du) = f\qquad\text{in $\Omega$}, \end{equation*} whose ellipticity degenerates in a fixed bounded and convex set $E\subset\mathbb{R}^n$ with $0\in \mathrm{Int}\, E$. Here, $\Omega\subset\mathbb{R}^n$ denotes a bounded domain, and $\mathcal{F} \colon \Omega\times\mathbb{R}^n \to\mathbb{R}_{\geq 0}$ is a function with the properties: for any $x\in\Omega$, the mapping $\xi\mapsto \mathcal{F}(x,\xi)$ is regular outside $E$ and vanishes entirely within this set. Additionally, we assume $f\in L^{n+\sigma}(\Omega)$ for some $\sigma > 0$, representing an arbitrary datum. Our main result establishes the regularity \begin{equation*} \mathcal{K}(Du)\in C^0(\Omega) \end{equation*} for any continuous function $\mathcal{K}\in C^0(\mathbb{R}^n)$ vanishing on $E$.
Let $n\ge 3$, $0<m<\frac{n-2}{n}$, $\alpha=\frac{2\beta-1}{1-m}$ and $\frac{2}{1-m}<\frac{\alpha}{\beta}<\frac{n-2}{m}$. We give a new direct proof using fixed point method on the existence of singular radially symmetric forward self-similar solution of the form $V(x,t)=t^{-\alpha} f(t^{-\beta}x)$ $\forall x\in\mathbb{R}^n\setminus\{0\}$, $t>0$, for the fast diffusion equation $u_t=\Delta (u^m/m)$ in $(\mathbb{R}^n\setminus\{0\})\times (0,\infty)$, where $f$ satisfies \begin{equation*} \Delta (f^m/m) + \alpha f + \beta x \cdot \nabla f =0 \quad \text{in} \; \mathbb{R}^n\setminus\{0\} \end{equation*} with $\lim_{|x| \to 0} |x|^{ \frac{\alpha}{\beta}}f(x)=A$ and $\lim_{|x| \to \infty}f(x) = D_A$ for some constants $A>0$, $D_A > 0$. We also obtain an asymptotic expansion of such singular radially symmetric solution $f$ near the origin. We will also prove the asymptotic large time behaviour of the singular solutions of the fast diffusion equation $u_t= \Delta (u^m/m)$ in $(\mathbb{R}^n\setminus\{0\})\times (0,\infty)$, $u(x,0)=u_0(x)$ in $\mathbb{R}^n\setminus\{0\}$, satisfying the condition $A_1|x|^{-\gamma}\leq u_0(x)\leq A_2|x|^{-\gamma}$ in $\mathbb{R}^n\setminus\{0\}$, for some constants $A_2>A_1>0$ and $n\le\gamma<\frac{n-2}{m}$.
In this paper we prove existence and uniqueness of energy solutionns for singular problems with absorption driven by local-nonlocal operators. Moreover, we establish a comparison principle à la Talenti, leading to a gain of summability result for the solutions of these problems.
We prove that instantaneous gelation (i.e., instantaneous loss of mass) occurs for coagulation equations with sum-type kernels of homogeneity greater than one which vanish on the diagonal. Our proof includes solutions that are Radon measures if we exclude the case of initial data that are a single Dirac delta.
We analyse the scattering operator associated with the defocusing nonlinear Schrödinger equation which captures the evolution of solutions over an infinite time-interval under the nonlinear flow of this equation. The asymptotic nature of the scattering operator (involving unbounded time) makes its computation particularly challenging. We overcome this by exploiting the space-time compactification provided by the lens transform, marking the first use of this technique in numerical simulations. This results in a highly efficient and reliable methodology for computing the scattering operator in various regimes. In developing this approach we introduce and prove several new identities and theoretical properties of the scattering operator. We support our construction with several numerical experiments which we show to agree with known analytical properties of the scattering operator, and also address the case of long-range scattering for the one-dimensional cubic Schrödinger equation. Our simulations permit us to further explore regimes beyond current analytical understanding, and lead us to formulate new conjectures concerning fixed and rotating points of the operator, as well as its existence in the long-range setting for both defocusing and focusing cases.
The work is devoted to the study of the inverse problem of determining the right-hand side of a nonlinear subdiffusion equation with a Caputo derivative with respect to time. Nonlinearity of the equation means that the right-hand side of the equation depends nonlinearly on the solution of the equation. The inverse problem consists of reconstructing the coefficient of the right-hand side, which depends on both time and spatial variables, under a measurement in an integral form. Similar inverse problems were previously studied in the case when the right-hand side depends only on time or on a spatial variable. A weak solution is sought by the Galerkin method. A priori estimates are proved, and with their help, the existence and uniqueness of a solution to the inverse problem under consideration are established. It is noteworthy that the results obtained are new for diffusion equations as well.
In this paper, we consider the initial boundary value problem of a doubly nonlinear parabolic equation with nonlinear perturbation. We impose the homogeneous Dirichlet condition on this problem. We aim to reduce the growth condition of the nonlinear term and the largeness of exponent as possible so that we can apply our result to both singular and degenerate type parabolic equations. We use in our proof the $L ^{\infty }$-estimate of the time-discrete equation derived in the previous work and apply the so-called $L ^{\infty }$-energy method to the time-discrete problem. We also discuss the uniqueness of solution by using the previous result.
In this work, we concern ourselves with the evolution of a droplet of an ideal fluid in two dimensions, which has nontrivial bulk vorticity and is only subject to the effects of surface tension. We construct initial data with initial domain being a disk and initial velocity being arbitrarily small, such that the vorticity Hessian grows exponentially infinitely in time.
In this paper, we establish a sharp stability inequality on the Heisenberg group for functions that are close to the sum of m weakly interacting Jerison-Lee bubbles. As a consequence, we obtain a sharp quantitative stability of global compactness for non-negative functions on Heisenberg group.
We study the propagation of plasmons on graphene. The problem is considered in two dimensions with a transverse magnetic (TM) electromagnetic field. The graphene material is assumed to be flat and is modeled as a conductive sheet. This leads to a jump condition for the magnetic field on the sheet where it is related to the current density on the sheet. The current density itself satisfies Drude's law. The model then consists of Maxwell's equations coupled to current density on the sheet. To make the problem more computationally approachable, we develop a time-dependent integro-differential equation for the current density. This effectively reduces the problem to one space dimension. A finite difference method is proposed to solve the resulting equation. Numerical examples are given to illustrate previously unreported behavior of the system.
We study a very general class of first-order linear hyperbolic systems that both become weakly hyperbolic and contain lower-order coefficients that blow up at a single time $t = 0$. In "critical" weakly hyperbolic settings, it is well-known that solutions lose a finite amount of regularity at the degenerate time $t = 0$. In this paper, we both improve upon the results in the weakly hyperbolic setting, and we extend this analysis to systems containing critically singular coefficients, which may also exhibit significantly modified asymptotics at $t = 0$. In particular, we give precise quantifications for (1) the asymptotics of solutions as $t$ approaches $0$; (2) the scattering problem of solving the system with asymptotic data at $t = 0$; and (3) the loss of regularity due to the degeneracies at $t = 0$. Finally, we discuss a variety of applications for these results, including to weakly hyperbolic and singular wave equations, equations of higher order, and equations arising from relativity and cosmology, e.g. at big bang singularities.
In this paper, we consider a space-time fractional partial differential equation with a reactive term. We describe the speed of invasion of its fundamental solution, extending recent results in this topic, which had been proved for the one dimensional spatial setting and some fractional parameters involved in the equation. The key tool to achieve these results in a wide range of cases (any spatial dimension, any time and spatial fractional differential parameters, and any polynomial speed of invasion) is the subordination.
Fractional differential equations are powerful mathematical descriptors for intricate physical phenomena in a compact form. However, compared to integer ordinary or partial differential equations, solving fractional differential equations can be challenging considering the intricate details involved in their numerical solutions. Robust data-driven solutions hence can be of great interest for solving fractional differential equations. In the recent years, fractional physics-informed neural network has appeared as a platform for solving fractional differential equations and till now, efforts have been made to improve its performance. In this work, we present a fully detailed interrogation of fractional physics-informed neural networks with different foundations to solve different categories of fractional differential equations: fractional ordinary differntial equation, as well as two and three dimensional fractional partial differential equations. These equations are solved employing two numerical methods based on the Caputo formalism. We show that these platforms are generally able to accurately solve the equations with minor discrepancies at initial times. Nonetheless, since in Caputo formalism, the value of a fractional derivative at each point requires the function's value in all of its previous history, it is computationally burdensome. Here, we discuss strategies to improve accuracy of fractional physics-informed neural networks solutions without imposing heavy computational costs.
The paper introduces a Poisson-type problem on a mixed-dimensional structure combining a Euclidean domain and a lower-dimensional self-similar component touching a compact surface (interface). The lower-dimensional piece is a so-called infinite metric tree (one-dimensional branching structure), and the key ingredient of the study is a rigorous definition of the gluing conditions between the two components. These constructions are based on the recent concept of embedded trace maps and some abstract machineries derived from a suitable Green-type formula. The problem is then reduced to the study of Fredholm properties of a linear combination of Dirichlet-to-Neumann maps for the tree and the Euclidean domain, which yields desired existence and uniqueness results. One also shows that finite sections of tree can be used for an efficient approximation of the solutions.
Constant Mean Curvature (CMC) 1-immersions of surfaces into hyperbolic 3-manifolds are natural and yet rather curious objects in hyperbolic geometry with interesting applications. Firstly, Bryant revealed surprising relations between (CMC) $1$-immersions of surfaces into $\mathbb H^3$ (Bryant surfaces) and (cousins) minimal immersions into $\mathbb E^3.$ In addition, the interest to (CMC) immersions of a surface $S$ (closed, orientable, with genus $\mathfrak{g} \geq2$) into hyperbolic 3-manifolds was motivated by Uhlenbeck in connection to irreducible representations of the fundamental group $\pi_{1}(S)$ into $PSL(2,\mathbb{C}).$ However a (CMC) 1-immersed compact surface is likely to develop singularities (punctures at finitely many points), and indeed in our analysis the prescribed value 1 of the mean curvature enters as a "critical" parameter. In fact, Huang-Lucia-Tarantello showed that (CMC) $c$-immersions of $S$ into hyperbolic 3-manifolds exist for $|c | <1$ and are parametrized by elements of the tangent bundle of the Teichmueller space of $S.$ More importantly, (CMC) $1$-immersions are attained only as "limits" for $|c| \to 1^-$ . In general the passage to the limit can be prevented by possible blow-up phenomena captured in terms of the Kodaira map and its suitable extension respectively for genus $\mathfrak{g}=2$ and $\mathfrak{g}=3.$ Here we handle the case of surfaces of any genus. In Theorem , we are able to encompass the blow up situation in terms of an appropriate "orthogonality" condition. Subsequently, we can provide the existence and uniqueness of (CMC) 1-immersions under an appropriate "generic" condition, see Theorem 2.